Combining Experimental and Quantum Chemical Study of 2-(5-Nitro-1,3-Dihydro Benzimidazol-2-Ylidene)-3-Oxo-3-(2-Oxo-2H-Chromen-3-yl) Propanenitrile as Copper Corrosion Inhibitor in Nitric Acid Solution ()
1. Introduction
Inhibitors corrosion search has become a sensitive issue [1] [2]. Although corrosion degrades certain materials and metal equipment, it is necessary to find inhibitors capable of reducing this phenomenon without harming the community and creating damage to environment. In fact, some inorganic substances including phosphates, chromates, dichromates, silicates, bromates, arsenates, tungstates, molybdates, chlorides and their derivatives were used for a long time to combat metals dissolution. Unfortunately these compounds [3] are mostly carcinogenic and pollute environment.
The current trend is to search for non-toxic organic molecules for metal inhibition corrosion, in accordance with the standards governing environmental protection [4] [5] [6] [7]. Thus this constraint has led several researchers to be interested in pharmaceutical and plant extracts molecules to prevent metals corrosion [8] [9] [10] [11]. Copper, which is a metal massively used because of its excellent mechanical and chemical properties, is not spared to corrosion [12] [13]. During its use, copper is covered with impurities or corrosion products. The removal of these corrosion products requires acidic solutions because they are often difficult to clean. These cleaning operations require some acids such as nitric acid. These acids often dissolve copper equipment, thus affecting the users on a safety and economic level.
Based on all this information, it is advisable to consider copper protection against corrosion by the use of inhibitors based on drugs or plant extracts which are generally very low toxicity and stable at high temperatures [14] [15] [16]. Some drugs such as antibacterial drugs [17], diclofenac sodium drug [18], paracetamol [8], melatonin drug [19], azithromycin [20], and vitamins [4] have been used to inhibit corrosion of metals.
For this reason, 2-(5-nitro-1,3-dihydrobenzimidazol-2-ylidene)-3-oxo-3-(2-oxo-2H-chromen-3-yl) propanenitrile was chosen for this work in order to evaluate its influence in corrosive environment to mitigate copper dissolution. This compound is derived from benzimidazolyl 2-cyanochalcones, it is a drug with the ability to treat infections caused by microscopic fungi and yeasts [21]. This molecule which has a therapeutic function could meet the standards of non-toxicity and environmental regulations based on eco-friendly and biodegradable. In addition, this molecule has in its molecular structure oxygen (O) and nitrogen (N) atoms which according to the literature [6] [9] can facilitate electronic transactions between inhibitor and copper.
It turns out that despite the high cost of experimental methods, these different electronic transactions between metal and inhibitor remain unexplained. However, density functional theory, which is a low cost computational method, is currently being used more and more to explain the inhibition mechanism [22] [23]. This theory relates that molecule inhibition performance depends on the protective layer quality which resulting from the chemical bond between metal and inhibitor molecules [24] [25].
In this work, it seems appropriate to study the inhibition power of this synthesized organic molecule which is 2-(5-nitro-1,3-dihydrobenzimidazol-2-ylidene)-3-oxo-3-(2-oxo-2H -chromen-3-yl) propanenitrile on copper corrosion in 1 M HNO3 solution by using combined experimental and theoretical approaches.
2. Experimental
2.1. Inhibitor Synthesis Process
The inhibitor used in this work was 2-(5-nitro-1,3-dihydrobenzimidazol-2-ylidene)-3-oxo-3-(2-oxo-2H-chromen-3-yl) propanenitrile (NDCP). Its preparation consists in first reacting ortho-phenylenediamines with ethyl cyanoacetate used as reagent and solvent. The product obtained is brought to reflux in toluene with 2-oxo-2H-chromene-3-carbonyl chloride presence. After precipitation, the precipitate obtained is filtered and recrystallized in water [26]. This synthesis process is illustrated by the following reaction scheme (Scheme 1):
Physico-chemical properties: NMR 1H, 13C (DMSO-d6, δ ppm):
• 1H NMR: (300 MHz, DMSO-d6, δ ppm): 12.98 (2H, br s, NH); 8.57 (1H, s, CO-C=CH)
• 13C NMR: (75 MHz, DMSO-d6, δ ppm): 187.0 (C=O ketone); 115.8 (C≡N); 63.5 (C=C-CN)
Figure 1 indicates synthesized molecule appearance.
Scheme 1. Synthesis procedure of 2-(5-nitro-1,3-dihydrobenzimidazol-2-ylidene)-3-oxo-3-(2-oxo-2H-chromen-3-yl) propanenitrile (NDCP).
Figure 1. Images of a synthesized sample of NDCP.
2.2. Chemicals
The essential chemicals used for copper samples treatment are:
• HNO3 from Merck with purity: 70% nitric. Using this commercial solution, a 1 M concentration serving as a blank for various tests was prepared.
• Acetone from Sigma Aldrich with purity: 99.5% was used to degrease copper coupons
2.3. Copper Pretreatment
Copper coupons were in form of rod measuring 10 mm in length and 2.2 mm of diameter; they were cut in commercial copper of purity 99.5%. To perform copper coupons pretreatment, these coupons were successively polished with fine quality emery papers with grain ranging from 150 to 600, cleaned in acetone, washed with double distilled water, then dried in a proofer (ASTEL) and weighed (m1).
2.4. Gravimetric Measurements
Initially pre-treated copper coupons were completely immersed in an Erlenmeyer flask containing 50 mL of 1 M HNO3 without and with different concentrations 2-(5-nitro-1,3-dihydrobenzimidazol-2-ylidene)-3-oxo-3-(2-oxo-2H-chromen-3-yl) propanenitrile. After one hour immersion, the coupons retrieved, washed again with double distilled water, dried in a proofer before re-weighing (m2). The tests were repeated at temperatures ranging from 298 K to 323 K. The temperature was controlled by a thermostat water bath (Frigitherm) and an analytical balance from KERN & SOHN GmbH (precision: ±0.1 mg) was used to perform the weighing. The mass loss (∆m) was the difference of initial mass of pre-treated coupons and the mass immersed in the blank solution without or with NDCP after treatment. This mass loss is the average value from the repetition of each test. From the mass loss results, the parameters such as corrosion rate (W), surface coverage (θ) and inhibition efficiency (IE) were calculated. The expressions for the calculation of these parameters are outlined as follows:
(1)
(2)
(3)
Δm: mass loss (g) m1 and m2 are respectively, the mass (g) before and after immersion in the solution test; t: immersion time (h); S: total surface of coupon (cm2); w0 and w; are respectively copper corrosion rates in the absence and presence of NDCP.
2.5. Quantum Chemical Computations
To evaluate molecule adsorption properties, DFT (Density functional theory) method were used in this work. The choice of this method is based on electronic correlation consideration where the wave function is replaced by electronic density. The exchange and correlation contributions are commonly treated separately and then combined to give a complete functional. Functionals that permits to link the electron density to energy are determined to access the calculations. This method reaches similar results to other methods (ab initio, Hartree-Fock, post Hartree-Fock.) with less computation time. All calculations were performed by Gaussian 09 W [27] software using B3LYP functional with three-parameter Becke (B3) Lee-Yang-Parr [28] [29], whose analytical form is given by the following relation:
(4)
where
.
To obtain a good accuracy and precision in these theoretical computations results and geometry optimization, basis 6 - 31 G (d, p) was used. The optimized minimum energy geometrical configuration using Gauss-View 5 software and NDCP molecular are given by Figure 2.
NDCP molecular structure
NDCP optimized structure
Figure 2. NDCP molecular and optimized structures.
Quantum chemical parameters involved in NDCP inhibition mechanism that were determined by this technique are: highest occupied molecular orbital energy (EHOMO), lowest unoccupied molecular orbital energy (ELUMO), energy gap (ΔE = ELUMO − EHOMO), dipole moment (μ), total energy (ET) electron affinity (A = −ELUMO) [30], ionization energy (I = −EHOMO) [30]. Some parameters such electronegativity (χ), hardness (η) softness (σ) affinity and electrophylicity index (ω) are related to ionization energy (I) and electron affinity (A). They are expressed from the following relationships [31] [32]:
(5)
(6)
(7)
(8)
Fraction of electron transferred (ΔN) from the inhibitor to copper surface was computed according to the following relation [33]:
(9)
where
[34] and
[35].
NDCP local selectivity parameters such as Fukui functions and dual descriptor responsible for nucleophilic and electrophilic attacks was also determined [36] [37] [38]. The following equations were used to determine these parameters:
(10)
(11)
(12)
where
and
are respectively nucleophilic and electrophilic Fukui functions,
,
and
are the electronic population of atomk in
,N and
electrons systems .
3. Results and Discussion
3.1. Gravimetric Results Interpretation
In order to have significant information about NDCP effect in copper corrosion, gravimetric tests were carried out. The results of these tests have permitted to plot the evolution of corrosion rate (W) and the inhibition efficiency (IE) as a function of NDCP concentration and corrosive solution temperature. These plots are illustrated in Figure 3 and Figure 4.
Figure 3. Corrosion rate versus NDPC concentration for various temperatures.
Figure 4. NDPC Inhibition efficiency (IE) versus concentrations for different temperatures.
It can be seen in Figure 3 that corrosion rate increases with increasing corrosive medium temperature and NDCP absence, while this rate decreases considerably in NDCP presence. Although the increase in temperature favors copper oxidation, it is clear that for high temperatures, NDCP presence extremely reduces copper dissolution. This dissolution is characterized by Cu2+ ions formation in reaction medium.
Figure 4 analysis shows that inhibition efficiency increases both increasing concentration and temperature. These observations relate that when temperature evolves, Cu2+ ions are massively produced. So, NDCP presence in nitric solution promotes a protective layer formation with these ions. This protective whose thickness depends on NDCP concentration and which forms on copper surface, preserves it from the aggressive environment, thus reducing its dissolution. However, increasing temperature and concentration facilitates NDCP adsorption on copper surface.
Based on the above results, it was found that at low concentration and high temperatures copper dissolution is completely reduced by NDCP. Moreover, for at T = 323 K and with Cinh = 0.2 mM, NDCP inhibition efficiency is 92.20%. Similar results were obtained by previous studies [8] [39] [40] [41].
3.2. Adsorption Isotherm and Thermodynamic Adsorption Parameters
Gravimetric results attest that copper corrosion inhibition comes from NDCP adsorption on said metal surface. It is then necessary to determine the isotherm model and thermodynamic adsorption parameters which permit to judge NDCP absorptive behavior.
To find the isotherm model suitable for NDCP behavior, several attempts have been made at some ordinary isotherms which are El-Awady, Freundlich, Temkin, Langmuir and Frumkin. The fits obtained with different graphs indicate that experimental data from gravimetric measurements are in better agreement with Langmuir isotherm. Indeed, regression coefficient (R2) of straight lines obtained during the study of this model are closer to unit. Langmuir adsorption isotherm is provided by following corresponding equation [42]:
(13)
Cinh: inhibitor concentration; Kads adsorption equilibrium constant; θ: surface coverage rate.
This model representation is given by the plot of Cinh/θ versus Cinh (Figure 5).
Moreover, the standard free energy of adsorption which is a thermodynamic parameter can be calculated according to the following relation [43]:
(14)
Figure 5. Langmuir adsorption isotherm of NDCP.
where R is gas constant, T is absolute temperature, Kadsadsorption equilibrium constant and 55.5 is water concentration ( mol/L) in solution.
Standard adsorption enthalpy (
) and entropy (
) are determined by applying the following equation [44]:
(15)
Figure 6 illustrates the plot of
versus temperature.
and
values were determined respectively from the intercept and the slope of the straight line obtained. All the quantities related to Langmuir isotherm and thermodynamic adsorption parameters are listed in Table 1.
Figure 6. Evolution of
as a function of temperature.
Table 1. NDCP adsorption parameters obtained by gravimetric estimations.
A review of Table 1 relates that the slopes of straight lines obtained are approximately equal to unity. In this case, NDCP adsorbs on copper according to Langmuir isotherm. This isotherm mentions that there is a fixed number of sites on copper surface and that each of these sites can adsorb only one particle, it is thus that a monolayer of NDCP is adsorbed on metal surface. Moreover, the interactions between adsorbed NDCP particles and metal surface are negligible and adsorption energy is constant.
As can be seen in Table 1, the negative
values suggest that NDCP adsorption onto copper is spontaneous and the layer from this adsorption is stable [45]. The results relating to
attest that the values of this quantity are between −20 kJ/mol and −40 kJ/mol, thus translating the existence of two adsorption modes: chemisorption and physisorption [45]. These adsorption modes indicate that covalent bonds are therefore established between the molecule and the metal surface as well as electrostatic interactions between the inhibitor and copper surface. As for the values of standard adsorption enthalpy (
) and entropy (
), they are positive, this translates respectively the endothermic character of NDCP adsorption and increasing disorder during NDCP molecules adsorption on copper surface [17]. This disorder is caused by water molecules desorption.
3.3. Corrosion Kinetic Analysis
Gravimetric tests proved that NDCP inhibition efficiency depends on temperature variations. These results confirm that the temperature is an element allowing to control copper dissolution process in the presence and absence of NDCP. Thus kinetic parameters of this process were determined from Arrhenius relationship and the transition state described as follows [46]:
(16)
(17)
W: corrosion rate, Ea: activation energy of metal dissolution reaction, R: universal gas constant, T: temperature, A: Arrhenius preexponential constant.
: activation entropy,
: activation enthalpy, R universal gas constant,
: Avogadro number, h: Planck’s constant.
Figure 7 shows the plot of
versus
. The slopes of the straight lines
obtained permit to determine activation energy (Ea) in absence and presence of NDCP.
and
are determined from the plot of
versus
, where the slopes
and intercepts
of the straight lines (Figure 8) lead to respective values of
and
. The values of the different activation quantities are recorded in Table 2.
Table 2. Activation parameters of copper dissolution in 1 M HNO3.
Inspection of the results obtained show that the activation energy (Ea) decreases as inhibitor concentration increases. It was observed that (Ea) value in blank solution is higher than the values of Ea in NDCP presence. The decrease of Ea in inhibited solutions attests that copper dissolution is strongly attenuated in NDCP presence [47]. Indeed in the corrosive solution copper dissolves rapidly by losing electrons, which electrons are replaced by those of the inhibitor from heteroatoms (N, O) and π bonds. These findings translate that NDCP is absorbed chemically on copper surface. Furthermore activation energy values of dissolution process in absence and NDCP presence are greater than 20 kJ∙mol−1, which means that the whole process is controlled by the surface reaction [47]. In this case, it appears that chemical and physical adsorption are present during NDCP inhibition process and this is done mainly by chemical adsorption. Activation enthalpy values are positive mentioning the endothermic character of copper dissolution process of in the absence and presence of NDCP [48]. These values decrease in NDCP presence, it reveals that copper dissolution becomes slow in studied molecule presence. It is also found that the difference:
is practically constant, translating that
and
vary in the same way. Consequently, this result obeys the thermodynamic relation between these two quantities:
[49]. Negative values of
which decrease in NDCP presence reveal that Cu-inhibitor complex formation considerably reduces the disorder created by water molecules desorption [48].
3.4. Adsorption Type Analysis
Adsorption and activation parameters have meant that both types of adsorption exist in NDCP inhibition. To have complete information to distinguish physisorption and chemisorption, Adejo-Ekwenchi isotherm has been used [50]. The relationship that characterizes this isotherm is reported below.
(18)
The representation associated with this model is exemplified in Figure 9. The parameters corresponding to this model are listed in Table 3.
Figure 9. Adejo-Ekwenchi adsorption model plots for NDCP for various temperatures.
Table 3. Adejo-Ekwenchi model parameters.
According to the results from Adejo-Ekwenchi model (Table 3), it appears that parameter b values decreases from 298 K to 308 K. This decrease relates that NDCP adsorbs on copper surface by electrostatic bonds in range of 298 K to 308 K [50]. Inhibition efficiencies values obtained by means of gravimetric tests in this range of temperatures are low. These results attest that the electrostatic interactions that occur therefore do not properly promote good corrosion inhibition of copper, because these electrostatic bonds are weak and are sensitive to temperature increase. In contrast, for temperatures range from 308 K to 323 K, parameter b increases, indicating that NDCP adsorbs onto metal surface through chemical bonds [50]. These bonds are formed by electron sharing between NDCP molecules and unsaturated copper d-orbitals, thus forming covalent bonds that are strong bonds. These bonds allow to maintain a good inhibition efficiency when the temperature increases, which justifies the high IE% values obtained at high temperature.
3.5. Quantum Chemical Interpretation
3.5.1. Global Reactivity and HOMO-LUMO Orbitals
Gravimetric estimations revealed that NDCP has a good inhibition performance. This performance is due to its absorptive properties. In order to deepen the inhibition mechanism characteristics, quantum chemical calculations were performed. Thus, the global parameters resulting from these calculations were requested to explain experimental results obtained and to better understand NDCP inhibition action. These parameters were determined with DFT at B3LYP/6-31G (d, p) basis set and are recorded in Table 4.
An organic compound reactivity depends on its ability to donate or receive electrons [51]. The parameters indicative of electron donation and acceptance are respectively HOMO and LUMO orbitals energies [52]. Indeed, according to
Table 4. NDCP global reactivity parameters calculated using B3LYP/6-31G(d, p).
the values obtained in literature [53] [54] it was found that NDCP EHOMO value is high, which results in NDCP having a greater ability to donate electrons to copper which has low energy empty molecular orbitals. In contrast, a lower ELUMO value as displayed by the one obtained with NDCP means that this molecule is likely to accept electrons from copper appropriate occupied orbitals. These donation and acceptance properties that NDCP possesses justify its strong adsorption to copper surface. An organic compound reactivity evolution is also controlled by energy gap (ΔE). In effect when ΔE decreases, molecule reactivity increases and when ΔE increases molecule reactivity decreases [55]. ΔE value obtained in this study reveals that the electronic exchanges molecule-copper are very favorable, which contributes to establish chemical bonds with copper. These bonds reinforce its adsorption capacity on the metal surface thus justifying the high efficiencies obtained experimentally. Figure 10 illustrates the molecular energy gap (molecular energy gap) which permits to judge the donation and electron acceptance capacities of NDCP.
Although some authors argue that a lower dipole moment (μ) facilitates inhibitor adsorption to metal surface [56], others report that a high dipole moment, ensures a good reactivity of the molecule which can increase the corrosion inhibition efficiency [17] [57]. In our study, the value of the dipole moment for NDCP (6.303 D) comparing to other values reported in previous works is high [58] [59], which justifies the better inhibition performance of the compound in tested solution.
Molecule electron donation and acceptance properties can also be discussed with ionization potential (I) and electron affinity (A) [60]. The low value of I and high value of A displayed in Table 4 confirms these properties for NDCP.
The electronegativity (χ) represents the electron attraction capacity of a molecule [61]. Electronegativity value obtained is lower than that of copper, this indicates NDCP electrons are strongly attracted to copper. This electron transfer is corroborated by the positive value of fraction of transferred electrons (∆N) and low hardness (η) [62]. Therefore, the electron loss that copper undergoes in 1 M HNO3 is systematically reduced by an electron donation from NDCP. This electron donation contributes to covalent bonds formation with copper. The high value of global softness (σ) relates that NDCP is a soft molecule therefore, it
Figure 10. Energy diagram of HOMO-LUMO orbitals for NDCP.
possesses a high chemical reactivity [63]. Which high reactivity would favor its strong adsorption on copper surface hence a good inhibition efficiency.
Electrophilicity index (ω) is a parameter which permits to estimate a compound capacity to accept or to give electrons [64]. In our work high value of ω attests that NDCP is a good electrophile and has a good ability to accept electrons from copper. Global reactivity parameters certified that electron transfer from NDCP to copper and vice versa are very favorable. Negative value total energy claims this assertion [65]. Analysis global reactivity descriptors are consistent with gravimetric estimations.
3.5.2. Local Reactivity
Global reactivity study indicated that NDCP adsorption on copper surface as estimated by gravimetric tests is related to its propensity to donate and accept electrons from appropriate occupied copper orbitals. Then, it is necessary to study local molecule reactivity. This reactivity consists in locating within molecule the atoms which are responsible for nucleophilic and electrophilic attacks [66]. Indeed NDCP molecular configuration indicates that this molecule has 38 atoms among which are 9 heteroatoms (N, O). The presence of these heteroatoms and aromatic rings can influence the behavior of certain atoms in the molecule. Therefore the centers of attack which permit to clarify the electronic properties will be elucidated by Fukui functions and dual descriptor. The values of these different parameters are listed in Table 5.
According to the local reactivity indices values reported in Table 5, C(26) and N(36) atoms are centers for electrophilic attacks associated with HOMO region and nucleophilic attacks associated with LUMO region, respectively. These centers are confirmed by
high value and the lowest value of
possessed by C(26) and the highest values
and
that N(36) has. The most favourable reactivity sites for adsorption are indicated by Figure 11.
Table 5. Mulliken atomic charges, Fukui function and dual descriptor.
Figure 11. Optimized geometry for NDCP obtained by DFT at B3LYP/6-31G(d, p).
In fact carbon C(26) surrounded by two heteroatoms O(37) and O(38) having solitary electron pairs will be filled with electrons. Thus this electron rich site is suitable to bind to metal surface by giving electrons to Cu2+ ions ([Ar]3d9). Whereas nitrogen N(36) which is only linked to C(27) will be depleted in electrons because of its electron pair delocalization. This site will therefore receive electrons from copper. Finally, the identified sites are responsible for chemical bonds formation.
3.5.3. Mechanism of Inhibition
Global reactivity study explained the experimental results, while NDCP reactivity sites were located through local parameters. Therefore, it is significant to clarify NDCP mechanism inhibition for copper corrosion in 1 M HNO3. In fact, copper oxidizes in 1 M HNO3 to Cu2+ while NDCP inhibitor molecules can easily be protonated in the said solution according to the equation below:
When negatively charged
ions attach to copper surface, electrostatic interactions are created between these ions and inhibitor protonated form. Furthermore, when NDCP neutral and protonated forms adsorb onto copper surface, coordination bonds are formed. These bonds arise from donor-acceptor interactions between π-electrons of aromatic rings, N(36) and C(26) atoms of inhibitor and copper vacant d-orbitals. This chemical adsorption process creates a protective film that will slow down copper oxidation. This protective film becomes thicker when temperature increases. Indeed, the increase of the temperature facilitates copper dissolution, thus a great Cu2+ ions number is created, metal electrons and those of NDCP become very mobile. This high mobility reinforces donor-acceptor interactions which increases protective film formation on copper surface, which explains the increase of inhibition efficiency at high temperature.
4. Conclusions
The following highlights were deduced from combining theoretical and experimental results obtained in this work:
• Gravimetric estimations mentioned that NDCP has good inhibition potentiality for copper corrosion in 1 M HNO3.
• NDCP corrosion inhibition ability generally results from a protective film formation on copper surface whose thickness increases with increasing temperature and inhibitor concentration.
• NDCP adsorption on copper surface obeys Langmuir model isotherm and adsorption parameters indicate that this process is spontaneous and endothermic.
•
negative values, the activation parameters and Adejo-Ekwenchi model confirmed that NDCP adsorption on copper surface occurred chemically and physically.
• The study of the corrosion kinetic parameters permitted to verify thermodynamic relation which is:
• Through DFT-based quantum chemical calculations at B3LYP/6-31-G(d, p) basis set, a correlation between the parameters related to electronic and molecular structures of NDCP and its ability to delay copper corrosion process was established
• Theoretical calculations provided a better understanding of NDCP inhibition action mechanism and are in full agreement with experimental results.
• Future work will focus on QSPR modeling (at 25˚C), which consists in establishing a mathematical relationship between quantum chemical parameters and inhibition efficiency and also perform surface studies before and after copper immersion in study solution.