Abstract

The density functional theory (DFT) was used to investigate the corrosion inhibition of three inhibitors on zinc. Quantum chemical parameters such as the energy of highest occupied molecular orbital (), the energy of lowest unoccupied molecular orbital (), energy gap , hardness , softness , electrophilicity index , the fraction of electrons transferred from inhibitor molecule to the metal surface, energy change when both processes occur, namely, charge transfer to the molecule and backdonation from the molecule (), natural charge , and Fukui functions have been calculated by using B3LYP/6-31+G(d) basis set. The relation between the inhibition efficiency and quantum chemical parameters has been discussed in order to elucidate the inhibition mechanism of the chloro-N-(4-methoxybenzylidene)aniline derivatives.

1. Introduction

Quantitative structure activity relationship (QSAR) has been a subject of intense interest in many disciplines of chemistry. In addition, the results of quantum chemical calculations are obtained without laboratory measurements, thus saving time and equipment and alleviating safety and disposal concerns [1, 2]. Density functional theory (DFT) has a potential application towards the design and development of organic corrosion inhibitors in corrosion field [3]. These facts have made quantum calculations to be very powerful tool for studying corrosion inhibition mechanism. Zinc is one of the most vital nonferrous metals having extensive use in metallic coating [4]. The use of zinc sheets as roofing materials and in galvanization of steel makes the study of its corrosion an important theme [5]. Zinc is often attacked by aggressive media, such as acids, bases, and salt solutions [68].

In acidic solutions, zinc has high solubility and dissolves with the formation of Zn2+ ions. Since the solubility of zinc decreases with increasing pH in acidic solutions, precipitation of Zn(OH)2 occurs when a certain pH value is reached. In alkaline solutions, with , the solubility of zinc increases with increasing pH, and in the high pH range, zinc oxide and hydroxides tend to dissolve with the formation of zincate ions. Zinc oxide is a common corrosion product in diluted solutions of salts such as NaCl or Na2SO4 with zinc hydroxide present in various amounts as a minor component. For scale removal and cleaning of zinc surface with acid solutions, the use of organic inhibitors is one of the most practical methods for protection against corrosion, especially in acidic media [913]. The aim of this work is to study the effect of quantum chemical parameters of 2-Chloro-N-(4-methoxybenzylidene)aniline (2-CNMBA), 3-chloro-N-(4-methoxybenzylidene)aniline (3-CNMBA), and 4-chloro-N-(4-methoxybenzylidene)aniline (4-CNMBA) on the inhibition efficiency of corrosion of zinc by using density functional theory (DFT) methods. For this purpose, the relation between the inhibition efficiencies of these compounds with quantum parameters has been investigated. The Fukui functions and local softness(s) values obtained were utilised to elucidate the reactive sites in the studied molecules.

2. Materials and Methods

2.1. Preparation of Schiff Bases

The Schiff bases were synthesized by condensation of 4-methoxybenzaldehyde with 2-chloroaniline, 3-chloroaniline, and 4-chloroaniline in the presence of ethyl alcohol as per the procedure described by Shah et al. [14]. 2-CNMBA (boiling point: 362.4°C, registry number is not available) is dark reddish liquid, 3-CNMBA (boiling point: 335.9°C, registry number: 24776-57-6) is brownish liquid, and 4-CNMBA (melting point: 335.9°C, registry number: 15485-22-0) is a white crystalline solid. These compounds are insoluble in water but soluble in ethanol, methanol, ethyl acetate, and THF. These compounds have previously been tested as inhibitors for Al-Mg alloy in HCl [15]. The molecular structures of the Schiff bases are shown in Figures 1(a), 1(b), and 1(c). The compounds were characterised through IR spectroscopy as shown in Figures 2(a), 2(b), and 2(c). 4-CNMBA has −C−H stretching of aromatic ring (3131 cm−1), >C=N– (1648 cm−1), –C–O–C– stretching in aromatic ring (1257 cm−1), –C–Cl stretching in aromatic ring (766 cm−1) and 3-CNMBA has −C−H stretching of aromatic ring (3164 cm−1), >C=N– (1634 cm−1), –C–O–C– stretching in aromatic ring (1240 cm−1) and –C–Cl stretching in aromatic ring (737 cm−1). 2-CNMBA has −C−H stretching aromatic ring (3120 cm−1), >C=N (1642 cm−1), –C–O–C– stretching in aromatic ring (1263 cm−1) and –C–Cl stretching in aromatic ring (744 cm−1).

2.2. Weight Loss Measurements

Rectangular specimens of electrolytic zinc of size (thickness = 0.139 cm) with a small hole of about 2 mm diameter just near one end of the specimen were used for the determination of the corrosion rate. The specimens were grinded using successively “0” to “0000” Oakey emery paper. The polishing was done using jewellers rouge, which gave mirror-like finish. The specimens were washed with ethanol and finally with distilled water and then dried at room temperature. The dried and weighed specimens were exposed to 230 mL of 0.5 M HCl solution containing controlled additions of various Schiff bases in the range of 0.0–0.5% (w/v). One specimen was suspended by a pyrex glass hook in each beaker of the test solution, which was open to the air at °C, to the same depth of about 1.5 cm below the surface of the liquid. The weight of the specimen before and after immersion was determined using mettler balance—M5 type.

The percentage inhibition efficiencies were calculated using the following equation: where is weight loss of zinc in uninhibited acid and is weight loss of zinc in inhibited acid.

2.3. Quantum Chemical Calculations

All the quantum chemical calculations have been carried out with Gaussian 09 programme package [16]. In our calculation we have used B3LYP, a hybrid functional of the DFT method, which consists of the Becke’s three parameters; exact exchange functional B3 combined with the nonlocal gradient corrected correlation functional of Lee-Yang-Par (LYP) has been used along with 6-31+G(d) basis set. In the process of geometry optimisation for the fully relaxed method, convergence of all the calculations has been confirmed by the absence of imaginary frequencies. The aim of our calculation is to calculate the following quantum chemical indices: the energy of highest occupied molecular orbital , the energy of lowest unoccupied molecular orbital , energy gap , hardness , softness , electrophilicity index , the fraction of electrons transferred from inhibitor molecule to the metal surface, and energy change when both processes occur, namely, charge transfer to the molecule and backdonation from the molecule , and correlate these with the experimental observations.

The number of electrons transferred from the inhibitor molecule to the metal surface was calculated according to Pearson [17, 18] as shown in (2) Here and denote the absolute electronegativity of zinc and the inhibitor molecules, respectively. and denote the absolute hardness of zinc and the inhibitor molecules, respectively. The fraction of electron transferred is calculated for zinc metal, and the experimental work function of polycrystalline zinc (8.126 eV) [19, 20] was employed for electronegativity, and global hardness of zinc was taken as zero since for bulk zinc. Ionisation potential and electron affinity are related to and as follows [21, 22]: The quantities absolute electronegativity and global hardness are related to electron affinity and ionization energy as follows [23]: Global softness is defined as [24] Recently, a new global chemical reactivity parameter has been introduced and known as electrophilicity index . It is defined as [19, 25] This is proposed as a measure of the electrophilic power of a molecule. When two molecules react, which one will act as an electrophile or nucleophile will depend on the fact that which one has a higher or lower electrophilicity index. The electrophilicity is a measure of reactivity that quantitatively classifies the global electrophilic nature of a molecule within a reactive scale. Partial atomic charges and atomic Fukui functions are good indicators of selectivity, that is, the region on the molecule on which certain type of reactions is likely to occur. The regions of a molecule where the Fukui function is large are chemically softer than the regions where the Fukui function is small. By involving the hard soft acid base (HSAB) principle in the local sense, one may establish the behaviour of the different sites with respect to hard or soft reagents [26]. The Fukui functions have been calculated by taking the finite difference approximations from natural population analysis of atoms in all the Schiff bases: Here “” is natural charge of atom in the molecule, that is, the electron density at a point “” in space around the molecule. The “” corresponds to the number of electrons in the neutral molecule. “” corresponds to an anion, with an electron added to the LUMO of the neutral molecule. “” corresponds to a cation, with an electron removed from the HOMO of the neutral molecule. All the calculations were performed at the ground state geometry.

Reactivity of similar atoms of different molecules can also be compared by calculating condensed softness indices from the Fukui functions and global softness by using the following expression [27]: According to the simple charge transfer model for donation and backdonation of charges, proposed by Gomez et al. [28, 29], if both processes occur (namely, charge transfer to the molecule and backdonation from the molecule) the energy change is directly proportional to the hardness of the molecule as shown above in  (10).

3. Results and Discussions

3.1. Weight Loss Measurements

The effect of addition of different chloro-N-(4-methoxybenzylidene)aniline derivatives at various concentrations on the zinc corrosion in 0.5 M HCl solutions is studied by weight loss measurements at 35°C ± 0.5°C after 30 min immersion. The values of inhibition efficiencies (%IE) obtained from (1) and corrosion rates obtained from weight loss measurements for the three chloro-N-(4-methoxybenzylidene)aniline derivatives at different concentrations in 0.5 M HCl are listed in Table 1. It is found that inhibition efficiency increases with increasing inhibitor concentration, while corrosion rate decreases with inhibitor concentration. The inhibition of zinc corrosion can be attributed to the adsorption of the inhibitors at the zinc/HCl solution interface. All Schiff bases are good inhibitors showing more than 99% inhibition efficiency at 0.5% (w/v) concentrations. Good performance of these compounds as corrosion inhibitors for zinc in 0.5 M HCl solutions may be due to the presence of the >C=N, OCH3 groups and benzene rings in their structures. From Table 1, it is clear that the order of the inhibition efficiency of these Schiff bases is as follows: 4-CNMBA ≥ 3-CNMBA > 2-CNMBA.

3.2. Computational Study

We have used DFT method in order to understand if any structural differences, induced by the position of chlorine on the N-(4-methoxybenzylidene) aniline, are related to the observed differences of the corrosion inhibition efficiency. DFT calculations and correlations were performed on the three inhibitors, that is, 2-CNMBA, 3-CNMBA, and 4-CNMBA. It is generally accepted that the values of indicate the electron donating ability of the molecule and the inhibition efficiency increases with increasing value. High values indicate that the molecule has a tendency to donate electrons to the unoccupied orbitals of the metal with low energy empty orbital [30]. An increase in values facilitates the adsorption and therefore corrosion inhibition. indicates the ability of the molecules to accept electrons; with lower values, a greater adsorption ability and better corrosion inhibition can be expected [31]. The optimized structures of all three Schiff bases, namely, (a) 4-CNMBA, (b) 3-CNMBA, and (c) 2-CNMBA with numbering schemes are depicted in Figure 3 which is shown below.

The total energy values given in Table 2 show the stability of the molecules which is 4-CNMBA > 3-CNMBA ≥ 2-CNMBA. This trend is the same as that experimentally.

The value of of 4-CNMBA is the highest which correlates with the experimentally determined inhibition efficiency and the values of favour 3-CNMBA. The separation energy, , is an important parameter and it is a function of reactivity of the inhibitor molecule towards the adsorption on metallic surface. As decreases, the reactivity of the molecule increases leading to increase in the inhibition efficiency of the molecule [31]. The calculation from Table 2 shows the decreasing trend for the : 2-CNMBA (4.2179 eV) > 3-CNMBA (4.1906 eV) > 4-CNMBA (4.1362 eV) which follows the same order of inhibition efficiency as obtained experimentally.

The calculations show an obvious correlation between the molecular area of the molecules and the inhibition efficiency. The inhibition efficiency increases as molecular area increases due to the increase of the contact area between molecule and surface. It is clear from Table 2 that 4-CNMBA has the highest molecular area that probably increases its adsorption on the metal surface and hence increases the inhibition efficiency. Global hardness and global softness are important properties to measure the molecular stability and reactivity. A hard molecule has a large energy gap and a soft molecule has a small energy gap. Soft molecules are more reactive than hard ones because they could easily offer electrons to an acceptor. In the corrosion system, the inhibitor acts as a Lewis base, while the metal acts as a Lewis acid. Bulk metals are soft acids and thus soft base inhibitors are most effective for acidic corrosion inhibition of these metals. It is also in accordance with HSAB principle. Normally, the molecule with least value of global hardness (hence highest value of global softness) is expected to have the highest inhibition efficiency [27]. It is clear from the calculation that 4-CNMBA has the lowest hardness and the highest softness which follow the same trend as obtained experimentally.

The electrophilicity index , which shows the ability of the inhibitor molecule to accept electrons, follows the trend: 3-CNMBA (3.766) > 4-CNMBA (3.7123) > 2-CNMBA (3.6152). 3-CNMBA exhibits the highest value of electrophilicity (Table 2) which confirms its highest capacity to accept electrons. There is strong correlation between this property and inhibition efficiency.

The number of electrons transferred was also calculated and tabulated in Table 2. The indicates the tendency of a molecule to donate electrons to the metal surface; therefore, the inhibition efficiency increases with increasing electron donating ability of these inhibitors to the metal surface [29] and it follows the order 3-CNMBA (0.991) ≤ 2-CNMBA (1.008) < 4-CNMBA (1.0172). The results indicate that value of 4-CNMBA is highest which correlates strongly with experimental inhibition efficiencies. Thus, highest fraction of electrons transferred is associated with the best inhibitor (4-CNMBA). Ionization energy is a fundamental descriptor of the chemical reactivity of atoms and molecules. High ionization energy indicates high stability and chemical inertness, while small ionization energy indicates high reactivity of the atoms and molecules. The low ionization energy of 4-CNMBA indicates its high inhibition efficiency.

The dipole moment provides information about the polarity of the whole molecule. High dipole moment means greater molecular polarity, which probably gives rise to high chemical reactivity [32]. It relates to the dipole-dipole interaction of inhibitor molecule and metal surface. Table 2 shows that 4-CNMBA has the highest dipole moment. it is due to electron releasing effect (+R effect) of –OCH3 and electron withdrawing effect of –Cl atom (−I effect), both of which are at the para position to the benzene rings attached to the > C=N group which will increase its molecular reactivity, leading to stronger adsorption onto the metal surface compared to the adsorption of 3-CNMBA and 2-CNMBA.

The unoccupied orbitals of zinc can accept electrons from inhibitor molecules to form a coordinate bond. Also the inhibitor molecules can accept electrons from zinc atom using its antibonding orbitals to form backdonating bond. These donation and backdonation processes strengthen the adsorption of these inhibitors onto the zinc surface. When and , the charge transferred to a molecule followed by backdonation from a molecule is energetically favourable. In Table 2, the calculated values for the inhibitors followed the order 4-CNMBA > 3-CNMBA > 2-CNMBA which indicates that backdonation is favoured for the 4-CNMBA, which is the best inhibitor.

As we know, frontier orbital theory is useful in predicting the adsorption centres of the inhibitors responsible for the interaction with surface metal atoms. Figures 4(a), 4(b), and 4(c) show the HOMO and LUMO orbital contributions for the neutral species of the studied inhibitor molecules 4-CNMBA, 3-CNMBA, and 2-CNMBA, respectively. For all the molecules, the HOMO densities were concentrated on both rings, Cl atom, >C=N group, and O-atom of methoxy group. This means that these are active sites of the molecules responsible for interaction with metal surface.

The LUMO molecular orbital was uniformly distributed over rings as well as on >C=N group, while contribution from Cl atom is negligible. Due to the presence of an iminic (>C=N–) group and rings, these molecules could function as good inhibitors.

Local charges are important in many chemical reactions and physicochemical properties of compound [26]. Tables 3, 4, and 5 show that most of the carbon atoms (C1, C2, C3, C4, C5, C12, C20, C21, C22, C23, and C25), oxygen atom (O11), nitrogen atom (N18), and chlorine atom (Cl29 in case of 3-CNMBA and 4-CNMBA) possess the excess of negative charge. Among them, highest negative charge is located on the nitrogen (N18) and oxygen atom (O11), in all the three studied inhibitors. It confirmed that electron donation ability of atoms to the unoccupied orbital of metal is directly related to negative atomic charges of adsorption sites. Therefore, these atoms should be the active adsorption sites.

The preferred site for nucleophilic attack is the atom in the molecule where Fukui function has the highest value as it is associated with the LUMO and measures the reactivity towards a donor reagent. The preferred site for electrophilic attack is the atom in the molecule where the value of Fukui function is the highest as it is associated with the HOMO and measures reactivity toward an acceptor reagent. For radical attack, the preferred site is one with highest value. An analysis of Fukui functions depicted in Tables 35 shows that the atoms O11 and N18 are preferred sites for electrophilic attack, as these sites have higher values of Fukui function . The atoms C6, C16, and O11 are most susceptible sites for nucleophilic attack since these sites have highest values of Fukui function . Chlorine Cl29 has very low values of Fukui indices. For the simplest transfer of electrons, adsorption could occur at the part of the molecule, where local softness has the highest value. A high value of indicates high nucleophilicity and a high value of indicates high electrophilicity. The results shown in Tables 35 indicate that the condensed local softness indices and follow the same trend as that of Fukui functions.

3.3. Quantitative Structure and Activity Relationship (QSAR) Consideration

Theoretical results discussed above show that there is no simple relation between quantum chemical parameters and inhibition efficiency of studied inhibitors. In order to correlate the quantum chemical parameters of the inhibitors and their experimental inhibition efficiencies, attempts were made to fit % IE to linear and nonlinear models by using XLSTAT [19, 33]. The following quantum chemical parameters were considered; these were the energy of the highest occupied molecular orbital (eV), energy of lowest unoccupied molecular orbital (eV), energy gap (eV), and the dipole moment (D). The nonlinear regression using polynomial equation of degree two yielded nearly as good results as found by experiments. The polynomial model yielded the following where is the concentration of inhibitor, is the coefficient of determination, and RMSE is defined as where IEpred is the predicted inhibition efficiency, is the experimental inhibition efficiency, and is the number of observations considered.

The experimental and calculated efficiencies are listed in Figure 5. The present model was successfully used to predict corrosion inhibition efficiencies of zinc in acidic medium for intermediate concentrations of inhibitors. It was clearly insufficient at very high concentrations where efficiencies do not increase with increase in concentrations of the inhibitors.

Inhibition efficiency is best described by dipole moment and remaining three variables , , and equally contribute in predicting inhibition efficiency.

4. Conclusions

4-CNMBA, 3-CNMBA, and 2-CNMBA are effective inhibitors of corrosion of zinc exposed to 0.5 M HCl solution. Their inhibition efficiencies increase in the order of 4-CNMBA > 3-CNMBA > 2-CNMBA. The inhibition efficiencies tend to increase by increasing the inhibitor concentration. Percent inhibition efficiencies were explained in terms of molecular parameters. Since the calculated quantum chemical properties can change by an addition of a functional group to the molecule, the correlations between inhibition efficiency and molecular parameters can be used for preselection of new inhibitors. Understanding the phenomena of adsorption is key point in corrosion problems. Computational studies help to find the most probable adsorption sites for a molecule. This information can help to gain further insight into corrosion system, such as the most probable site for inhibitor adsorption. A QSAR equation was developed and used to predict the corrosion inhibition efficiency for studied inhibitors. The experimental measurements match the prediction of corrosion efficiencies of these compounds. QSAR approach may be useful to predict and design the structure and molecule that is suitable to be an effective corrosion inhibitor.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are grateful to the Chemistry Department, School of Sciences, Gujarat University, Gujarat, for the laboratory facilities and to the Central University of Gujarat for providing computational support.