ISSN : 0970 - 020X, ONLINE ISSN : 2231-5039
     FacebookTwitterLinkedinMendeley

Quantum chemical studies on the corrosion inhibition of some hector bases on mild steel in acidic medium

Majid Khodaei-Tehrani1* and Ali Niazi

1*Department of Chemistry, Arak branch, Islamic Azad University, Arak, Iran

DOI : http://dx.doi.org/10.13005/ojc/310150

Article Publishing History
Article Received on :
Article Accepted on :
Article Published : 27 Feb 2015
Article Metrics
ABSTRACT:

The density functional theory (DFT) at the B3LYP/6-31G (d,p) basis set level method were performed on three hector bases used as corrosion inhibitors; namely, 3-anilino-5-imino-4-phenyl-1, 2,4-thiadiazoline (AIPT), 3-anilino-5-imino-4-tolyl-1, 2,4-thiadiazoline (AITT), and 4-(4-chlorophenyl)-5-imino-N-phenyl-4,5-dihydro-1,2,4-thiadiazol-3-amine (AICT). They were used as corrosion inhibitors for mild steel in acidic medium in order to determine the relationship between molecular structure and their corresponding inhibition efficiency (%IE). The results of the quantum chemical calculations and experimental %IE were subjected to correlation analysis. This indicates that their inhibition effects are closely related to the highest occupied molecular orbital energy (EHOMO), the lowest unoccupied molecular orbital energy (ELUMO), the energy gap (ΔE), the hardness (η), the softness (σ), the electronegativity (χ), and the fraction of electrons transferred from the inhibitor molecule to the metal surface (ΔN). In addition, the local reactivity has been analyzed through the Fukui function. Two QSAR equations were developed and used to predict the corrosion inhibition efficiency for hector bases.

KEYWORDS:

hector bases; Inhibition of corrosion; Quantum chemical parameters; Fukui function

Download this article as: 

Copy the following to cite this article:

Khodaei-Tehrani M, Niazi A. Quantum chemical studies on the corrosion inhibition of some hector bases on mild steel in acidic medium. Orient J Chem 2015;31(1).


Copy the following to cite this URL:

Khodaei-Tehrani M, Niazi A. Quantum chemical studies on the corrosion inhibition of some hector bases on mild steel in acidic medium. Orient J Chem 2015;31(1). Available from: http://www.orientjchem.org/?p=7477


Introduction

The use of corrosion inhibitor is one of the most effective and highly cost-effective measures to mitigate corrosion of mild steel in acidic media [1-4]. Most of acid corrosion inhibitors are organic compounds containing electronegative atoms (such as sulfur, phosphorus, Nitrogen, and oxygen, etc.), the unsaturated bonds (such as, double bonds or triple bonds, etc.) and the plane conjugated systems including all kinds of aromatic cycles[5-11].It has been reported that the best efficiencies of N-containing inhibitors have been shown in hydrochloric acid [12]. The adsorption of the examined thiadiazole molecules can occur as a result of the formation of links between the d-orbital of iron atoms, which involve the displacement of water molecules from the metal surface by the lone sp2 electron pairs present on the N atoms of both the heterocyclic rings and S atom of thiadiazole ring. Quantum chemical calculations have been widely used to study the reaction mechanism of corrosion inhibition or even predict with confidence some experimentally unknown properties [13-16].

These studies have been successfully performed to link the corrosion inhibition efficiency and have provided specific information about charge distribution, molecular structure, the frontier molecular orbital energy (EHOMO, ELUMO, etc.) and adsorption processes of corrosion inhibitors.

The objective of this work is to investigate the dependence of inhibition efficiency of three hector bases used as corrosion inhibitors namely 3-anilino-5-imino-4-phenyl-1, 2,4-thiadiazoline (AIPT), 3-anilino-5-imino-4-tolyl-1, 2,4-thiadiazoline (AITT), and 4-(4-chlorophenyl)-5-imino-N-phenyl-4,5-dihydro-1,2,4-thiadiazol-3-amine (AICT)  and used as corrosion inhibitors towards mild steel in hydrochloric acid with their calculated quantum chemical parameters in the neutral forms using statistical tool .These inhibitors have been reported by M.A.Quraishi et al.[17].The study has been conducted by polarization measurements, impedance techniques and weight loss measurements. Using both electrochemical polarization and weight loss techniques, and found that all hector bases reduce the corrosion of steel.

 Table 1: Names, molecular structure, the abbreviation and atom numbering of the inhibitors investigated. Table1: Names, molecular structure, the abbreviation and atom numbering of the inhibitors investigated. 
Click here to View table

 

Computational Method

The structures optimized by using: DFT (density functional theory) with The Becke’s three-parameter hybrid functional [18] was combined with the Lee, Yang, and Parr (LYP) correlation functional [19] and denoted as B3 LYP theory and was employed in the DFT calculations using 6-31 G (d,p) basis set. In all cases, the quantum chemical calculations were performed with complete geometry optimizations are implemented by means of Gaussian 03 (Review B.05) program software package [20].

According to Koopman’s theorem [21], Frontier molecular orbitals (HOMO and LUMO) could be used to predict the adsorption centers of the inhibitors. For the simplest transfer of electrons, adsorption should occur at the part of the molecule where the softness, σ, which is a local property, has the highest value. The EHOMO and ELUMO of the inhibitor molecule are related to the electron affinity, A, and the ionization potential, I, respectively, are given by: I = EHOMO and A = – ELUMO .

Absolute electronegativity, χ, and absolute hardness, η, of the Inhibitor molecule were calculated by using the values of I and A and given by: χ = I+A/2 and η= I-A/2. And the global softness, σ, is the inverse of the global hardness: σ = 1/η

Thus the fraction of electrons transferred from the inhibitor to metallic surface, ΔN, as follow [22]:

formula1

Results snd Discussion

Quantum chemical calculations were carried out to research the effect of Geometric and electronic structural parameters on the inhibition efficiency of inhibitors and powerful tools for studying their adsorption mechanisms on the metal surface [23]. The optimized molecular structures with minimum energies obtained from the calculations are given in

 Fig. 1: The optimized molecular structures, HOMO and LUMO of the neutral inhibitor molecules using DFT/B3LYP/6-31G (d p). Figure1: The optimized molecular structures, HOMO and LUMO of the neutral inhibitor molecules using DFT/B3LYP/6-31G (d p). 

Click here to View figure

 

Quantum chemical parameters obtained from the calculations which are responsible for the inhibition efficiency of inhibitors, such as the energies of highest occupied molecular orbital (EHOMO), energy of lowest unoccupied molecular orbital (ELUMO), the separation energy (ELUMO-EHOMO), ΔE, representing the function of reactivity, the net charge on the functional group, dipole moment, , global hardness (η), softness ,(σ), and the fraction of electrons transferred from the inhibitor to iron surface(ΔN) for the neutral inhibitors are collected in Table 2 .

Table 2: The calculated quantum chemical parameters for the neutral inhibitors. Table2: The calculated quantum chemical parameters for the neutral inhibitors. 

Click here to View table

 

The hector bases may adsorb on the metal surface in the form of neutral molecules involving the displacement of water molecules from the metal surface and sharing of electrons between the nitrogen and sulfur atoms in the hector base molecule and the metal surface.

Neutral Inhibitor

Frontier molecular orbitals (HOMO and LUMO) may be used to predict the adsorption centers and involved in the activity properties of the inhibitor molecule. The higher the energetic level (HOMO) of the inhibitor, the greater the trend of offering electrons to unoccupied orbital of the iron metal, and the higher the corrosion inhibition efficiency for iron in HCl acid solutions; if the energy of HOMO level was determined for the inhibitor properties, the ranking of the compounds should be:

AICT >AITT >AIPT, EHOMO

This agrees well with the experimental observations. Another parameter of the molecular structure is energetic level of (LUMO), in addition, the lower the LUMO energy, the easier the acceptance of electrons from metal surface, the calculations also shows that AIPT has the lowest ELUMO, so it has the greatest ability to interact with the metal surface. If the energy of LUMO level was determined for the inhibitor properties, the ranking of the compounds should be:

AICT >AITT >AIPT; ELUMO

The trend across structures in the ELUMO of the compounds is not in good agreement with the trend in the inhibition efficiencies of the inhibitors.

The separation energy, (ΔE= ELUMO-EHOMO) is another important factor as a function of reactivity of the inhibitor molecule towards the adsorption on metallic surface. As ΔE decreases, the reactivity of the molecule increases leading to increase the inhibition efficiency of the molecule i.e. a molecule with a low energy gap is more polarizable and low kinetic stability and is generally associated with the high chemical activity and is named soft molecule[24]. The results obtained from quantum chemical calculation are listed in Table 2. The calculations indicate that AICT has the lowest value which means the highest reactivity among the other inhibitors. The order of reactivity in this case will be: AIPT >AITT >AICT; (ΔE)

The number of electrons transferred (ΔN) was also calculated depending on the quantum chemical method as in Eq. (1).Values of ΔN showed inhibition effect resulted from electrons donation .Using a theoretical χ value of 7 eV/mol according to Pearson’s electronegativity scale [25] and η value of 0eV/mol for iron atom [21], According to Lukovits’s study, [26], the inhibition efficiency increased with increasing electron-donating ability at the metal surface. In this study, the three selected organic compounds will be expected to inhibit the iron corrosion through donations of electrons to the iron surface that will be the electron acceptor as indicated in Table 2.

The order of ΔN is presented as follows; AICT >AITT >AIPT, and the order is the same as in the values of the experimental corrosion inhibition efficiency (%IE).

The dipole moment (μ) is an index that can also be used for the prediction of the direction of a corrosion inhibition process and provides information on the polarity of the molecule and it is also a good reactivity indicator. The dipole moment is a good reactivity indicator and also is related to the distribution of electrons in a molecule [27]. the literature is variable over the use of μ as a predictor for the direction of a corrosion inhibition reaction, it is agreed that the adsorption of polar compounds possessing high dipole moments on the metal surface should lead to better inhibition efficiency [28, 29] and should lead to worse inhibition efficiency but on the other hand, a survey of literature reveals that the dipole moment does not correlate well with the corrosion inhibition efficiencies of the inhibitors [30]. In the current work, the trend across structures in the dipole moment of the compounds is not in good agreement with the trend in the inhibition efficiencies of the inhibitors.

Pearson presented the quantities of electronic hardness (η) and softness (σ) in his hard–soft–acid–base principle (HSAB) [31] The species are classified as soft (hard) if their valence electrons are easy (hard) to polarize or to remove and a hard molecule has a large energy gap, and a soft molecule has a small energy gap. According to this theory, soft acids interact preferentially with soft bases and hard acids interact preferentially with hard bases. Metals are generally considered to be soft acid; consequently they would like interact with inhibitors that have high σ values and low η values [32, 33]. The σ values for the investigated compounds follow the trend AICT >AITT >AIPT, which is in good agreement with experimentally determined inhibition efficiency of the inhibitors.

The Table 2 shows the order of electronegativity as AIPT> AITT> AICT. Hence an increase in the difference of electronegativity between the metal and the inhibitor is observed in the order AICT >AITT >AIPT. According to Sander son’s electronegativity equalization principle [34], the electron flow happens from the molecule with the low electronegativity towards that with a higher value, until the chemical potentials are the same. AICT with a high electronegativity and low difference of electronegativity quickly reaches equalization and hence low reactivity is expected which in turn indicates low inhibition efficiency.

formula2

Local Reactivity

Fukui functions proposed by Parr and Yang [35], as measurement of local reactivity of the molecules indicative of the reactive regions and the nucleophilic and electrophilic behavior of the molecule. The Fukui function fk is defined as the first derivative of the electronic density p(r) with respect to the number of electrons N in a constant external potential v(r)

formula3,4

Where q(N), qk (N+1) and are determined as the atomic charges of the neutral, anionic and cationic species, respectively.

The fk+ measures the changes in density when the molecule gains electrons, and it corresponds to reactivity with respect to nucleophilic attack. The site for nucleophilic attack will be the place where the value of fkis a maximum. On the other hand, the site for electrophilic attack is controlled by the value of fk– and the highest value of fk– corresponds to a possible site for an electrophilic attack.

As can be seen from Table 3, the nitrogen atoms contain the surplus of negative charge. Among the nitrogen atoms, the highest negative charge is located on the N7 in AIPT and AITT. It is confirmed that the more negative the atomic charges of the adsorbed center, the more easily the atom donates its electrons to the unoccupied orbital of metal. Therefore, these atoms contain the surplus of negative charge should be the active adsorption sites. On the other hand, the sulfur atom contains the surplus of positive charge. The highest positive charge among the sulfur atoms related to sulfur atom of AICT. Thus, the N2, N4, N6 and N7 may be the most active adsorption sites. It is thought that the more positive the atomic charge leads to more easily the atom acceptance of metal electrons to the unoccupied orbital of sulfur atom. Further analysis of fukui functions (Table 3) show that the N1 and N3 are the most susceptible sites for the electrophilic attacks as they present the highest values fk– The highest value of fkand the most susceptible sites for the electrophilic attacks are S1 and N2 in AIPT and AITT and S1 in AICT.

 Table 3: The condensed Fukui functions on the atoms of the studied compounds. Table3: The condensed Fukui functions on the atoms of the studied compounds.  

Click here to View table

 

QSAR Study

In this investigation, quantitative structure and activity relationship (QSAR) has also been used to correlation between quantum chemical parameters and the observed inhibition efficiency of the studied inhibitors. The experimental determination the adsorption of the studied inhibitors on the mild steel surface follows Langmuir’s adsorption isotherm.

Formula5

The non-linear model has been employed:

formula6

Where IEcal% is the inhibition efficiency, A and B are the regression coefficients determined by regression analysis, Xj is a quantum chemical index characteristic for the molecule (j) and Cinh,i denotes the concentration of the inhibitor in experiment i.

Therefore, it is practical to use the linear and the non-linear mathematical models proposed by Lukovit,s et al. [38], in order to correlate the observed inhibition efficiency to the calculated quantum chemical parameters. The linear model approximates inhibition efficiency (IEcal %) has the form:

The best linear and non-linear equations resulting were using DFT/B3LYP/6-31 G (d, p) for from the combination of the EHOMO, ELUMO, ΔE, ΔN and μ parameters have these forms;

formula7

 

Formula8

The correlation coefficient values (R2) were obtained for the linear, resistance and non-linear models were 0.926 and 0.634. The best equation resulting from these models was obtained from linear model.

Conclusion

The correlation between the quantum chemical parameters and inhibition efficiency of some hector bases was investigated using DFT/B3LYP calculations. The inhibition efficiency of the inhibitor is closely related to the quantum chemical parameters, EHOMO, ƞ, σ, χ, and ΔN for the neutral inhibitors and no significant relationship was found with parameters, μ and ELUMO. It was shown by fukui functions that the N1 and N3 atoms are the most susceptible sites for the electrophilic attacks and the most susceptible sites for the electrophilic attacks are S1 and N2 atoms in AIPT and AITT molecules and S1 atom in AICT molecule. The inhibition efficiency of the inhibitors is closely related to the quantum chemical parameters and the linear model presented the best correlation with experimental data. The prediction of corrosion efficiencies of these compounds elegantly matched to the experimental measurements.

References

  1. Bentiss, F.; Lebrini, M.; Lagrenee, M. Corros. Sci. 2005, 47, 2915–2931
  2. Aljourani, J; Raeissi, K.; Golozar, M.A. Corros. Sci. 2009, 51, 1836–1843
  3. Obot, I.B.; Obi-Egbedi, N.O.; Umoren, S.A. Corros. Sci. 2009, 51, 276–282
  4. Zheludkevich, M.L.; Yasakau, K.A.; Poznyak, S.K.; Ferreira, M.G. Corros. Sci. 2005, 47, 3368–3383
  5. Wang, H.; Liu, R.; Xin, J. Corros. Sci. 2004, 46, 2455–2466
  6. Khaled, K.F.; Babic-Samardzija, K.; Hackerman, N. J. Appl. Electrochem. 2004, 34, 697–704
  7. Hosseini, M.; Mertens, S.; Ghorbani, M.; Arshadi, M. Mater. Chem. Phys. 2003, 78, 800-808
  8. Khaled, K.F.; Babic-Samradzija, K.; Hackerman, N. Electrochim. Acta 2005, 50, 2515–2520
  9. Abd El-Rehim, S.S.; Ibrahim, M.A.M.; Khaled, K.F. J. Appl. Electrochem 1999, 29, 593–599
  10. Bentiss, F.; Lagrenee, M.; Traisnel, M.; Hornez, J.C. Corros. Sci. 1999, 41, 789–803
  11. Sherif, E; Park, S.M. J Electrochem Soc.2005, 152, 428–433
  12. Schmitt, G. Br. Corros. J. 1984, 19, 165-176
  13. Awad, M. K.; Mustafa, M. R.; Abo Elnga, M. M. J. Mol. Struct. (Theochem) 2010, 959, 66-74
  14. Kosari, A.; Momeni, M.; Parvizi, R.; Zakeri, M.; Moayed, M.H.; Davoodi, A.; Eshghi, H. Corros. Sci. 2011, 53, 3058–3067.
  15. Khaled, K.F.; Amin, M.; Corros. Sci. 2009, 51, 1964–1975
  16. Growcock, F.B. Corros. Sci. 1989 45, 1003–1007
  17. Quraishi, M.A.; Sardar, R.; J.  Appl. Electrochem. 2003, 33, 1163–1168.
  18. Becke, A.D. J. Chem. Phys. 1993, 98, 5648 – 5652
  19. Lee, C.; Yang, W.; Parr, R.G. Phys. Rev. B: Condens. Matter 1988, 37, 785-789
  20. Frisch, M.J. et al. Gaussian 03, Revision B.05, Gaussian, Inc., Wallingford, CT, 2004.
  21. Sastri, V.S.; Perumareddi, J.R. Corrosion, 1996, 671, 617-622
  22. Pearson, R.G. Inorg. Chem. 1988, 27, 734–740.
  23. Fang, J; Li, J. J. Mol. Struct. (Theochem) 2002, 593, 179–185
  24. Fleming, Frontier Orbitals and Organic Chemical Reactions, John Wiley and Sons, New York, 1976.
  25. Pearson, R.G. Proc. Nati. Acad. Sci. 1986, 83, 8440-8441
  26. Lukovits, I.; Kalman, E.; Zucchi, F. 2001, 57, 3–8
  27. Issa, R.M.; Awad, M.K.; Atlam, F. M. Appl. Surf. Sci. 2008, 255, 2433-2441
  28. Ju, H.; Kai, Z.P.; Li, Y. Corros. Sci. 2008, 50, 865-871
  29. Popova, A.; Christov, M.; Raicheva, S.; Sokolova, E. Corros. Sci. 2004, 46, 1333–1350
  30. Ebenso, E.E. Bull. Electrochem. 2003, 19, 209–216
  31. Koch, E, Propellants Explos. Pyrotech. 2005, 30, 209-215


Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.