DFT Calculations for Corrosion Inhibition of Ferrous Alloys by Pyrazolopyrimidine Derivatives ()
Keywords:Pyrazolopyrimidine Derivative; Corrosion Inhibitors; Adsorption; DFT Calculations; Fukui Function
1. Introduction
Corrosion of metals is a major industrial problem that has attracted many investigation and researches. The use of inhibitors is one of the most practical methods to protect metals against corrosion. Most efficient inhibitors are organic compounds containing electronegative functional groups and π-electrons in triple or conjugated double bonds. Researchers conclude that the adsorption on the metal surface depends mainly on the physicochemical properties of the inhibitor group, such as the functional group, molecular electronic structure, electronic density at the donor atom, π-orbital character and the molecular size [1-]">3]. A number of heterocyclic compounds containing nitrogen, oxygen, and sulphur either in the aromatic or long chain carbon system have been reported to be effective inhibitors [4,]">5] . The inhibition efficiency has been closely related to the inhibitor adsorption abilities and the molecular properties for different kinds of organic compounds [6,]">7] . Organic compounds, which can donate electrons to unoccupied d-orbital of metal surface to form coordinate covalent bonds and can also accept free electrons from the metal surface by using their anti-bonding orbital to form feedback bonds, constitute excellent corrosion inhibitors [8,]">9].
A lot of work on heterocyclic inhibitors has been studied experimentally [10-]"> 14]. On the other hand, many reported theoretical studies in order to correlate between structural and electronic parameters and the inhibition efficiency [6,8,15-1]">9 ]. The results from these researches have been used to interpret very well many experimental phenomena.
Density functional theory (DFT) [20,]">21] following the approach of Kohn and Sham has proven to be an important tool in modern quantum chemistry because of its ability to include some effects of electron correlation at a greatly reduced computational cost [22-]"> 24]. Recently, DFT methods have been used to analyze the characteristics of the inhibitor/surface mechanism and to describe the structural nature of the inhibitor in corrosion process [6,25-]">28].
The object of the present paper is to carry out a DFT calculations on the electronic parameters of two pyrazolo[1,5-c]pyrimidine derivatives and they are:
5-tolyl-2-phenylpyrazolo[1,5-c]pyrimidine-7(6H)thione (Tolyl)
5-tolyl-2-phenylpyrazolo[1,5-c]pyrimidine-7(6H)one (Inon)
used as inhibitors, and to determine a relationship between some quantum chemical parameters obtained from the structure of the compounds and the inhibition efficiencies of corrosion of ferrous alloys obtained experimentally by Mahgoub et al. [12]. The chemical structures of the studied inhibitors are shown in Figure 1.
Computational calculations were obtained by means of DFT/B3LYP/6-31 + G(d,p) level of theory. Parameters such as, , energy gap (), dipole moment (), ionization energy (), electron affinity, absolute electronegativity (), absolute hardness (), absolute softness (), the fraction of electron transferred (), and the total energy () were calculated. The local reactivity has been analyzed by means of the Fukui indices, since they indicate the reactive region, in the form of the nucleophilic and electrophilic behavior of each atom in the molecule.
2. DFT Calculations
The present calculations were performed using Gaussian
(a)(b)
Figure 1. The molecular structures of the investigated inhibitors: (a) 5-tolyl-2-phenylpyrazolo [1,5-c] pyrimidine- 7(6H) thione (Tolyl), (b) 5-tolyl-2-phenylpyrazolo [1,5-c] pyrimidine-7(6H)one (Inon).
09 program package [29]. Geometry optimizations were conducted by DFT using Becke’s three parameter exchange functional (B3) [30], and includes a mixture of HF with DFT exchange terms associated with the gradient corrected correlation functional of Lee, Yang, and Parr (LYP) [31] and the 6-31 + G(d,p) basis set.
In geometry optimizations every bond length, bond angle and dihedral angle was allowed to relax, free of constraints. The nature of the stationary points was confirmed by vibrational frequency analysis, to verify that only real frequencies values (i.e. no imaginary frequency) were obtained for all geometries. In order to calculate the Fukui functions, the electron populations for neutral Tolyl, Tolyl cation, Tolyl anion, neutral Inon, Inon cation, and Inon anion species were calculated (using: the keyword Pop = Full). For the ionic forms, the calculations were performed at the optimized geometries of the neutral forms (using: single point energy calculation; charge = ±1; and multiplicity = Doublet). This avoids the problem of trying to do a geometry optimization on a species that may not be a stationary point on the potential energy surface.
Frontier molecular orbitals; highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) was used to predict the adsorption centers of the inhibitor molecule. 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. According to Koopman’s theorem [32], the energies of the HOMO and the LUMO orbitals of the inhibitor molecule are related to the ionization potential, , and the electron affinity, , respectively, by the following relationships:
(1)
(2)
Absolute electronegativity, c, and absolute hardness, h, of the inhibitor molecule are given [33]:
(3)
(4)
The softness is the inverse of the hardness [33]:
(5)
Electronegativity, hardness, and softness have proved to be very useful quantities in chemical reactivity theory. When two systems, Fe and inhibitor, are brought together, electrons will flow from lower (inhibitor) to higher (Fe), until the chemical potentials become equal.
The number of transferred electrons () was also calculated by using the equation below [32]:
(6)
where and denote the absolute electronegativity of iron and inhibitor molecule, respectively, and denote the absolute hardness of iron and the inhibitor molecule, respectively.
In this study, we use the theoretical value of and for the computation of number of transferred electrons [32]. The absolute electrophilicity index is given by [34]:
(7)
According to the definition, this index measures the tendency of chemical species to accept electrons. More reactive nucleophile is characterized by lower value of, and; and conversely more reactive electrophile is characterized by a higher value of, and. This new reactivity index measures the stabilization in energy when the system acquires an additional electronic charge from the environment.
3. Results and Discussion
3.1. Experimental Findings [12]
F. M. Mahgoub et al. investigated experimentally the inhibition of corrosion of ferrous alloys using a group of pyrazolo [1,5-c] pyrimidine as corrosion inhibitors. Using electrochemical measurements and corrosion tests such as Tafel extrapolation and linear polarization resistance techniques they indicated clearly a decrease in the corrosion rate in the presence of Tolyl and Inon inhibitors. The inhibition efficiencies of the Tolyl compound are markedly higher than that of the Inon compound (Table 1).
3.2. Theoretical Findings
Figure 2 shows the optimized structures of Tolyl and Inon inhibitors along with atomic numbering. For Tolyl and Inon molecules, the computed quantum chemical properties such as, , , and dipole moment () are given in Table2
Frontier molecular orbital (FMO): The interaction between the inhibitor and the metal is through the donation of the electrons from the inhibitor occupied orbitals (mainly from the) to the d-orbital of the metal [27], and also through the acceptance of the electrons from the d-orbital of the metal to the unoccupied orbitals (mainly to the) of the inhibitor. Thus, measures the tendency of donating electron by a molecule [35]. Therefore, higher value of indicates better tendency of donating electron, and enhancing the
Table 1 . Corrosion rate and the protection efficiency (%P) of 0.01 M of pyrazolo [1,5-c] pyrimidine inhibitors from Tafel and linear polarization plots for carbon steel in stagnant cooling water (200 ppm) at 25˚C [12].
(a)(b)
Figure 2. Optimized geometries of (a) Tolyl and (b) Inon molecules along with atomic numbering calculated at B3LYP/6-31 + G(d, p) level of theory.
Table 2. Calculated quantum chemical parameters of Tolyl and Inon molecules calculated at B3LYP/6-31 + G(d,p) level of theory.
adsorption of the inhibitor on the metal surface and therefore better inhibition efficiency. On the other hand, indicates the ability of the molecule to accept electrons, therefore, lower value of indicates better ability to accept electrons, and also this will enhance the adsorption of the inhibitor on the metal surface and therefore better inhibition efficiency. The binding ability of the inhibitor to the metal surface increases with increasing of the and decreasing of the energy values. Frontier molecular orbital diagrams of Tolyl and Inon is shown in Figure 3. From Table 2, the highest value of −0.21463 eV of
Figure 3. Frontier molecular orbitals (FMO) of Tolyl and Inon molecules calculated at B3LYP/6-31 + G(d,p) level of theory.
Tolyl compared to that of Inon −0.21943 eV, also the lowest value of −0.06821 eV of Tolyl compared to that of Inon −0.06535eV indicates the better inhibition efficiency of Tolyl, a result which is in total a agreement with the experimental finding.
The energy gap () is an important parameter as a function of reactivity of the inhibitor molecule towards the adsorption on the metal surface. As decreases the reactivity of the molecule increases leading to better inhibition efficiency [28]. Table 2 shows that Tolyl inhibitor has the lowest energy gap 0.14642 eV compared to of Inon inhibitor 0.15408 eV, with a difference equals 0.0077 eV; this means that Tolyl molecule could have better performance as carrion inhibitor than Inon molecule, a conclusion which is in total agreement with the experimental finding.
Dipole moment: The dipole moment () is another important electronic parameter that result from non-uniform distribution of charges on the various atoms in the molecule. The high value of dipole moment probably increases the adsorption between the inhibitor and the metal surface [36]. The energy of the deformability increases with the increase in, making the molecule easier to adsorb at the ferrous surface. In addition, the volume of the inhibitor molecules also increases with the increase of, this increase the contact area between the molecule and the surface of iron and increasing the corrosion inhibition ability of the inhibitor. In our study, the value 7.2231 Debye of Tolyl indicates its better inhibition efficiency compared to 6.8992 Debye of Inon which totally agrees with the experimental findings (Table 2).
Other computed quantum chemical properties such as ionization energy (), electron affinity, absolute electronegativity (), absolute hardness (), absolute softness (), absolute electrophilicity (), the fraction of electron transferred (), and the total energy () for Tolyl and Inon inhibitors are given in Table3
Ionization energy: Ionization energy () is a fundamental descriptor of the chemical reactivity of a toms and molecules. High ionization energy indicates high stability and chemical inertness and vice versa [37]. The low ionization energy 0.21463 eV of Tolyl indicates its high inhibition efficiency compared to 0.21943 eV of Inon which totally agrees with the experimental findings (Table 3).
Electronegativity: Table 3 shows the order of electronegativity as Tolyl < Inon. Hence, an increase in the difference of electronegativity between the metal and the inhibitor is observed in the order Tolyl > Inon. According to Sanderson’s electronegativity equalization principal [38], Inon 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, also conclusion in total agreement with the experimental finding.
Hardness and softness: Absolute hardness and softness are important properties to measure the molecular stability and reactivity. It is apparent that the chemical hardness fundamentally signifies the resistance towards the deformation or polarization of the electron cloud of the atoms, ions, or molecules under small perturbation of chemical reaction. A hard molecule has a large energy gap (large value) and a soft molecule has a small energy gap (small value) [1]. In this study Tolyl with low hardness value eV compared to that of the Inon eV has a low energy gap. Normally, the inhibitor with the least value of absolute hardness (hence, the highest value of absolute softness) is expected to have the highest inhibition efficiency [39]. Tolyl with the softness value of 13.65934 eV
Table 3 . Calculated quantum chemical parameters of Tolyl and Inon molecules calculated at B3LYP/6-31 + G(d,p) level of theory.
has the highest inhibition efficiency compared to 12.98027 eV of Inon; this result is also in agreement with the experimental finding (Table 3)
Number of electrons transferred: The number of electrons transferred () was also calculated and tabulated in Table3 If, the inhibition efficiency increases by increasing the electron-donating ability of these inhibitors to donate electrons to the metal surface [40] and it increases in the following order: Tolyl > Inon. The results indicate that values correlate strongly with experimental inhibition efficiencies. Thus, the highest fraction of electron transferred is associated with the best inhibitor Tolyl (), while, the least fraction is associated with the inhibitor that has the least inhibition efficiency Inon ().
Total energy: The total energy of Tolyl inhibitor is, while for Inon inhibitor it is, hence, the difference in the total energies between the two inhibitors is This result indicated that Tolyl inhibitor is favorably adsorbed through the active centers of adsorption on the ferrous surface (Table 3).
Local selectivity: The local selectivity of a corrosion inhibitor is best analyzed using the Fukui function [41, 42]. The Fukui indices permit the distinction of each part of a molecule on the basis of its chemical behavior due to different substituent functional groups. The Fukui function can be formally defined as:
(8)
where indicates that the differentiation is carried out under constant external potentials, and the functional derivative must be taken at constant number of electron,. If it is assumed that the total energy, , is a function of and is an exact differential, then the Maxwell relations between derivatives can be applied to derive the following equation:
(9)
where is the electron density. Equation (9) is the most standard presentation of the Fukui function. The Fukui function is provoked by the fact that if an electron is transferred to an electron molecule, it will tend to distribute so as to minimize the energy of the resulting electron system. The resulting change in electron density is the nucleophilic () and electrophilic () Fukui functions, which can be calculated using the finite difference approximation as follows:
(10)
(11)
where, , and are the electronic densities of anionic, neutral, and cationic forms of the atom with, , and electrons.
Calculated values of, , , , and for Tolyl and Inon are presented in Tables 4 and 5. The measures the changes of density when the molecule gain electron/s and it corresponds to reactivity with respect to nuecleophilic attack, thus, the site for nucleophilic attack is the site where the value of is maximum. On the other hand, corresponds to reactivity with respect to electrophilic attack or when the molecule loss electron/s, thus, the site for electrophilic attack is the site where the value of is maximum. For the Tolyl inhibitor, Table 4, it can be deduced that the sites for nucleophilic attack is in the carbon atom (C13). However, the sites for electrophilic attack are in the carbon atoms (C26 and C29). For the Inon inhibitor, Table 5, the sites for nucleophilic attack is in the carbon atom (C3). However, the sites for electrophilic attack are in the carbon atoms (C13, C18, and C24).
The and orbitals of Tolyl and Inon, Figure 2, clearly reveal the information that governs the nucleophilic and electrophilic attacks on the studied inhibitors. The information obtained from the and LUMO orbitals are consistent with the findings obtained from the Fukui function. The FMO diagram of Tolyl and Inon indicates the lack of electron cloud in in C13 and C3 for Tolyl and Inon inhibitors, respectively, which is confirmed by the Fukui function. In case of of Tolyl the dense electron cloud around (C26 and C29) indicates the site of electrophilic attack. The same is the case around (C13 and C18, and C24) in Inon as confirmed by the Fukui function too.
The local softness, and, for an atom can be expressed as the product of the Fukui function, and
, and the global softness, , as follows:
(12)
(13)
The local softness contains information similar to those obtained from Fukui function plus additional information about the molecular softness, which is related to the global reactivity with respect to a reaction partner. Calculated values of and are presented in Tables 4 and 5. From the results obtained, the sites for electrophilic and nucleophilic attack in the studied molecules are not similar.
4. Conclusion
Using the DFT/B3LYP/6-31 + G(d,p) level of theory, the inhibition efficiency of two pyrazolo [1,5-c] pyrimidine
derivatives is investigated leading to the following conclusions:
1) Through DFT calculations, a correlation between parameters related to the electronic and molecular structures of these derivatives and their ability to inhibit the corrosion process could be established.
2) The inhibition efficiency of these derivatives obtained quantum chemically increases with the increase in, and decrease in and energy gap. Tolyl has the highest inhibition efficiency because it had the highest energy.
3) The parameters like hardness and softness confirm the inhibition efficiency in order of Tolyl > Inon.
4) The total energy of the best inhibitor, Tolyl, is lower than that of the Inon inhibitor, thus, the Tolyl is favorably adsorbed through the active center/s of adsorption on the ferrous surface.
5) A direct relationship between the inhibition efficiency and the dipole movement and the number of transferred electrons.
6) From the local reactivity indices, Fukui function shows the nucleophilic and electrophilic attacking sites in pyrazolo [1,5-c] pyrimidine derivatives.
Finally, this study displays a good correlation between the theoretical and experimental data which confirm the reliability of the quantum chemical methods to study the inhibition of corrosion of metal surfaces.
NOTES