Corrosion behavior of copper in 1 M nitric acid containing either 2-mercaptobenzimidazole (MBI) or 2-thiomethylbenzimidazole (TMBI) was investigated experimentally and theoretically via weight loss method and quantum chemical approaches. It was found that the two compounds exhibit a very good performance as inhibitors for copper corrosion in the studied medium. Results show that the inhibition efficiencies increase with increasing temperature and increasing concentration of the inhibitors. It has been found that the studied compounds adsorb onto copper according to the modified Langmuir adsorption isotherm and the kinetic/thermodynamic isotherm of El-Awady. The thermodynamic adsorption parameters and kinetic corrosion parameters were determined and analyzed; on the bases of these parameters both physisorption and chemisorption were suggested for the studied compounds. Furthermore, the quantum chemical properties/descriptors most relevant to their potential action as corrosion inhibitors have been calculated. They include EHOMO, ELUMO, energy gap (ΔE), dipole moment (μ), hardness (η), softness (σ), the fractions of electrons transferred (ΔN), electrophilicity index (ω), and the total energy change (ΔET). The theoretical results were found to be consistent with the experimental data reported.
1. Introduction
Because of its excellent conductivity, good mechanical workability, and relatively low cost and reactivity, copper is one of the most important materials used widely in different industries, especially in central heating installations, car industry, energetic, oil refineries, sugar factories, marine environment, to name only a few of its various applications. Scale and corrosion products produced during the work of systems have some negative effects on their heat-exchange performance. Thus, acid washing is periodically carried out to descale and clean these systems. Acidic solutions are widely used in various industries for the cleaning of copper. The behavior of copper in acidic media is extensively investigated, and several ideas have been presented for the dissolution process [1, 2]. To avoid the base metal attack and to ensure the removal of corrosion products/scales alone, inhibitors are extensively used. Corrosion inhibitors are substances that protect metals against corrosion by decreasing the rate of corrosion processes. Many efficient inhibitors are heterocyclic organic compounds consisting of a π-system and/or O, N, P, or S heteroatoms [3–5]. Among them, azoles, benzotriazole, and its derivatives are the most often used organic inhibitors against copper corrosion. Unfortunately, it is proven that benzotriazole and its derivatives are highly toxic compounds; their use is limited therefore by the restrictions due to the protection of the environment [6]. Benzimidazole derivatives are well known because of their nontoxicity and good inhibitor effectiveness in various media [7, 8].
The efficiency of an organic compound as a successful inhibitor is mainly dependent on its ability to be adsorbed on the metal surface, which consists of the replacement of water molecule at a corroding interface. The adsorption of these compounds is influenced by the electronic structure of the inhibiting molecules [9, 10] and also by the steric factors, aromaticity, electron density at the donor atoms, and the presence of functional groups [11–14].
Adsorption can be described by two mean types of interactions [15].
Physisorption
It involves electrostatic forces between ionic charges or dipoles on the adsorbed species and the electric charge at the metal/solution interface. The heat of adsorption is low, and therefore this type of adsorption is stable only at relatively low temperatures.
Chemisorption
It involves charge sharing or charge transfer from the inhibitor molecule to the metal surface to form a coordinate type bond. In fact, electron transfer is typically for transition metals having vacant low-energy electron orbital. Chemisorption is typified by much stronger adsorption energy than physical adsorption. Such a bond is therefore more stable at higher temperature.
Thermodynamic adsorption parameters and kinetic corrosion parameters are a useful tool for clarifying the adsorption behaviour of an inhibitor. The reactivity of organic inhibitors is mainly interpreted by adsorption on the metal surface and depends on the molecular structure. In other words, the electronic structure of the organic compounds has a key influence on their corrosion inhibition efficiency. The molecular structure, including the electronic parameters, can be obtained by means of theoretical calculations by using the computational methodologies of quantum chemistry. Quantum chemical calculation has been used recently to explain the mechanism of corrosion inhibition [16–18] and proved to be a very powerful tool for studying the mechanism [19–21]. The survey of theoretical corrosion literature presented by Gece [19] demonstrates that quantum chemistry is a powerful tool to study the fundamental, molecular-level processes related to corrosion inhibition.The object of this paper is
to evaluate the thermodynamic adsorption and the activation parameters of the copper/inhibitor/1 M HNO3 system in the case of each of the selected compound. The choice of these compounds is based on molecular structure considerations; that is, these molecules are organic compounds having the same adsorption centers, but they only differ in the substituent CH3 on S in TMBI.
To correlate the effect of structure parameters of the two compounds through their quantum parameters (EHOMO,ELUMO,ΔE=EHOMO-EHOMO,μ,η,σ,ω) to explain their inhibition efficiency.
2. Materials and Methods2.1. Molecules Structures
Figure 1 gives the structures of 2-mercaptobenzimidazole (MBI) and 2-thiomethylbenzimidazole (TMBI). Their molecular weight is respectively M = 150.2 g/moL and M = 164, 04 g/moL. These organic compounds have been synthesized in the Laboratory according to Van Allan method [22]. Their molecular structures have been identified by RMN-1H and 13C spectroscopies and mass spectroscopy.
Optimized chemical structures of 2-thiomethylbenzimidazole (TMBI) and 2-mercaptobenzimidazole (MBI).
MBI
RMN 1H (DMSOd-6): 7.10–7.16 (2H,m), 7.22–7.34 (2H,m); 13C NMR (DMSOd-6): 114.08 (2C) 121.39 (2C), 128.44 (2C), 149.66 (C=N); SDM (E.I.): M+1 = 151. Analysis for C7H6N2S: C (55.97%), H (4.03%), N (18.65%), S (21.35%).
The weight loss method of monitoring corrosion rate is useful because of its simple application and reliability [23]. Several authors have reported on comparable agreement between weight-loss technique and other techniques of corrosion monitoring, including polarization measurement [24, 25], hydrogen evolution [26, 27], thermometric technique [28], and electrochemical impedance spectroscopy [29, 30].
The weight-loss measurements were performed with samples of copper in the form of rods measuring 10 mm in length and 2.2 mm in diameter that were cut from commercial pure copper (Cu 99.5%). The corrosive solution of 1 M HNO3 was prepared by dilution of analytical grade 65% HNO3 from MERCK with double-distilled water. The samples were polished successively with metallographic emery papers of increasing fineness of up to 600 grits and further with 5, 1, 0.5, and 0.3 μm alumina slurries (Buehler), washed thoroughly with double-distilled water, degreased, and dried with acetone. The samples were then kept in a desiccator, weighed and immersed in the corrosive medium (50 mL of 1 M HNO3) with or without the tested compounds. A water thermostat Selecta (Frigiterm) controlled to ±0.5°C maintained the temperatures ranging from 25°C to 65°C. After 1 hour, the specimens were removed, washed with double-distilled water, dried, kept in a desiccators, and then reweighed. All tests were made in aerated solutions and were run triplicate to guarantee the reliability of the results. The corrosion rate (W) was calculated using the following:
W=m1-m2St,
where m1 and m2 are, respectively, the weight (in g) before and after immersion in the test solution, S the total surface of the sample (in cm2), and t the immersion time (in h).
The inhibition efficiency IE(%) was then calculated using the following:
IE(%)=W0-WW0*100,
where W0 and W are, respectively, the corrosion rates of copper in the absence and presence of the tested compounds.
3. Results and Discussion3.1. Effect of the Concentration of the Inhibitors on Corrosion Rate of Copper
Figures 2(a) and 2(b) represent the corrosion rate as a function of inhibitors concentration at different temperatures. It is clear that, for all temperature, corrosion rate decreases with increasing concentration. In other words, the addition of the molecules to 1 M HNO3 solution retards copper corrosion, and the extend of retardation is concentration and temperature dependent. Similar observations have been seen in the literature [31]. However, we note that, for a given concentration, the corrosion rate increases with increasing temperature. This can be attributed to the fact that the rate of chemical reaction increases with increase in temperature.
Dependence of copper corrosion rate on the inhibitors concentration at different temperatures. (a): MBI (b): TMBI.
3.2. Effect of the Inhibitors Concentration on the Inhibiting Efficiency
Figures 3(a) and 3(b) show the evolution of the inhibiting efficiencies with inhibitors concentration at different temperatures. Inspections of these figures reveal that inhibition efficiency increases with increasing concentration of each inhibitor. We also note that inhibition efficiency increases with increasing temperature. Such behavior can be interpreted on the basis that the inhibitors exert their action by adsorbing themselves on the metal surface and that the increase in temperature leads to a shift of the equilibrium constant towards adsorption of the inhibitors molecules. For a certain inhibitor concentration, the inhibitive action of the studied compounds can be given in the following increasing order: MBI < TMBI. Figure 4 gives the example of C = 10−3 M. The ability of the molecule to adsorb on copper surface is probably dependent on the substituent on the sulfur atom. This result can be explained by orbital approach [32].
Variation of inhibition efficiency with inhibitors concentration for copper corrosion inhibition by (a): MBI (b): TMBI.
Evolution of the inhibiting efficiency of the studied compounds with the temperature for copper corrosion inhibition in 1 M HNO3.
3.3. Activation Parameters
Arrhenius suggested the famous equation which evaluates the temperature dependence of the rate constant as follows [33]:
logW=logA-Ea2.303RT.
Here, A is the frequency factor (in g·cm−2·h−1), Ea is the apparent activation energy (in J·mol−1), R is the perfect gas constant (R = 8.314 J·moL−1·K−1), and T is the absolute temperature (in K). Equation (3) predicts that the plot of log W versus 1/T should be a straight line. The slope of the line is (–Ea/2.303R), and the intercept of the line extrapolated (1/T = 0) gives logA.
On the other hand, activation thermodynamic parameters including change in enthalpy(ΔHa*), change in free energy (ΔGa*), and change in entropy (ΔSa*) of activation were calculated using the following [34]:
ΔHa*=Ea-RT,ΔGa*=2.303RT[log(kBTh)-logW],ΔSa*=(ΔHa*-ΔGa*T),
where KB is the Boltzmann constant (1.38·10−23J·K−1), h is the Planck constant (6.626·10−34 J·s), and R is the perfect gas constant. Figure 5 gives the representation of logW versus 1/T.
Arrhenius plots for copper corrosion in 1 M HNO3 (a): MBI (b): TMBI.
All the activation thermodynamic parameters (Ea, ΔHa*, ΔSa*, and ΔGa*) are estimated and listed in Table 1.
Corrosion kinetic parameters for copper in 1 M HNO3 in absence and presence of the investigated inhibitors.
Concentration (M)
Ea (kJ·mol−1)
ΔHa* (kJ·mol−1)
ΔSa*(J·mol−1·K−1)
ΔGa* (kJ·mol−1)
Blank
0
12.24
9.60
−233.3
83.85
MBI
5·10−5
6.49
3.85
−251.5
83.82
10−4
6.29
3.65
−254.1
84.45
5·10−4
6.14
3.49
−263.7
87.36
10−3
5.79
3.15
−269.2
88.77
TMBI
5·10−5
9.85
7.21
−241.0
83.86
10−4
9.48
6.84
−244.8
84.70
5·10−4
6.88
4.24
−261.3
87.32
10−3
4.92
2.27
−273.6
89.28
From Table 1 it can be seen that, in general, values of the apparent activation energy Ea for the inhibited solutions are lower than that for the uninhibited one, indicating a chemisorption process of adsorption [35]. The energetic barrier is lower, facilitating the formation of Cu2+ ions which interact with the studied inhibitors to form a protective film [36]. Indeed it is proven that copper is corroded in Cu2+ in HNO3 solution, and no oxide film is formed to protect the surface from corrosion [37, 38]. So in HNO3 solution the following reactions occur:
Cu2O+2H+⇄2Cu++H2OCuO+2H+⇄Cu2++H2O
Moreover, Cu+ ions undergo disproportionation according to the equation below:
2Cu+⇄Cu2++Cu
The average values of apparent activation energies for MBI and TMBI are, respectively, 6.18 and 7.78 kJ·moL−1. According to these values one can give the inhibitive character of the studied compounds in the following order: TMBI > MBI.
The positive sign of the activation enthalpy reflects the endothermic nature of copper dissolution process meaning that the dissolution of copper is difficult.
The negative values of ΔSa* in the inhibited and uninhibited systems impliy that activation complex in the rate-determining step represents an association rather than a dissociation step, meaning that a decrease disorder takes place on going from reactants to the activated complex [39].
3.4. Thermodynamic Adsorption Parameters
Since the corrosion inhibition process is based on the adsorption of the inhibitor molecules on the metal surface, it is essential to know the mode of adsorption and the adsorption isotherm that fits the experimental results. The values of surface coverage for various concentrations of the inhibitors under investigation must be used to explain the best adsorption isotherm. It is a widely held view by many authors that the adsorption of organic inhibitor molecules on metal surface is often a displacement reaction involving the removal of adsorbed water molecules from this surface [40, 41]
Orgsol+nH2Oads⟷Orgads+nH2Osol
The surface coverage (θ) was calculated as follows:
θ=IE(%)100.
Basic information on the interaction between the inhibitor and the copper surface can be provided by the adsorption isotherm. It is necessary to determine empirically which adsorption isotherm fits best to the surface coverage data in order to use the corrosion rate measurements to calculate the thermodynamic parameters pertaining inhibitor adsorption. In order to obtain the isotherm, linear relation between θ values and the inhibitor concentration C must be found. Attempts were made to fit the θ values to various isotherms including Langmuir, Temkin, Frumkin, Freundlich, and the kinetic/thermodynamic isotherm of El-Awady. The best fit was obtained with the Langmuir isotherm and that of El-Awady.
The Langmuir adsorption isotherm is described by the following:
Cθ=1K+C,
where C is the inhibitor concentration, θ is the fraction of the surface covered, and K the equilibrium constant of the adsorption process. Figure 6 shows the dependence of C/θ as a function of the concentration C for MBI and TMBI. The obtained plots are linear with correlation coefficients (R2) almost equal to unity. It is also found that all the values of the slopes are very close to unity (Table 2). The deviation of the slopes from unity is often interpreted as a sign that the adsorption species occupy more or less a typical adsorption site at the metal/solution interface [42]. The divergence of the slopes from unity is attributable to interactions between adsorbate species on the metal surface as well as changes in the adsorption heat with increasing surface coverage. The Langmuir adsorption isotherm cannot be applied rigorously. A modified Langmuir adsorption isotherm [43] can be applied to this phenomenon which is given by the corrected
Cθ=nKads+nC.
The equilibrium constant of adsorption is related to the standard free energy of adsorption ΔGads0by the following:
Kads=155.5⋅exp(-ΔGads0R.T),
where R is universal gas constant, T is the absolute temperature, and 55.5 is the molar concentration of water in the solution in mol·L−1. The adsorption equilibrium constant and the standard free energy of adsorption ΔGads0 (Table 2) were calculated for C = 10−3 M, using (10) and (11).
Calculated parameters from modified Langmuir isotherm and El-Awady adsorption isotherm at C = 10−3 M.
Inhibitor
TemperatureT(K)
Modified Langmuir isotherm
El-Awady isotherm
Kads
ΔGads°
Slope
R2
Kads
ΔGads°
1/y
R2
(105 M−1)
(kJ·mol−1)
(105 M−1)
(kJ·mol−1)
MBI
298
0.71
−31.93
0.968
0.999
0.63
−31.63
0.941
0.998
308
0.77
−33.19
0.992
0.999
0.71
−32.98
0.957
0.997
318
0.83
−34.48
1.005
0.999
0.85
−34.52
1.006
0.998
328
1.00
−36.06
1.004
0.999
0.91
−35.80
0.982
0.973
338
1.01
−37.16
1.016
0.999
1.16
−37.58
1.067
0.996
TMBI
298
0.80
−32.18
1.007
0.999
0.84
−32.33
1.029
0.984
308
0.96
−33.88
0.977
0.999
0.90
−33.58
0.949
0.994
318
1.11
−35.23
0.982
0.999
1.02
−35.00
0.964
0.990
328
1.15
−36.44
0.998
0.999
1.19
−36.55
1.009
0.995
338
1.25
−37.49
1.008
0.999
2.96
−40.20
1.051
0.996
Langmuir adsorption isotherm for the investigated inhibitors at different temperatures.
The experimental data were also well described by the El-Awady kinetic/thermodynamic adsorption isotherm. The characteristic equation of this model is given by
log(θ1-θ)=logK+ylogC,
where C is the concentration of the inhibitors, θ is the degree of surface coverage, Kads is the equilibrium constant of the adsorption process, and Kads=K1/y. In this model, the number of active sites y is included. Values of 1/y less than one imply multilayer adsorption, while 1/y greater than one suggests that a giving inhibitor molecule occupies more than one active site. Curve fitting of the data to the thermodynamic/kinetic model is shown in Figure 7. These plots give straight lines which clearly show that the data fit well to the isotherm. The values of Kads and ΔGads0 for C = 10−3 M are listed in Table 2. From the table, it can be seen that values of 1/y are very closed to the unity.
El-Awady adsorption isotherm for the investigated inhibitors at different temperatures.
The negative values of ΔGads0 indicate a spontaneous adsorption process of MBI and TMBI on copper [44]. It is generally accepted that, for the values of ΔGads0up to −20 kJ/moL, the type of adsorption was regarded as physisorption; the inhibition acts due to the electrostatic interactions between the charged molecules and the charged metallic surfaces, while the values around −40 kJ/moL or smaller were seen as chemisorption, which is due to the charge sharing or a transfer from the inhibitors molecules to the metal surface to form a covalent bond [45, 46]. The values of ΔGads0in our measurements range from −30 to −40 kJ·moL−1. it is suggested that the adsorption of these molecules involves both types of interactions, chemisorptions, and physisorption.
ΔGads0 is related to the enthalpy change and entropy change of adsorption process, ΔHads0 and ΔSads0, respectively by Gibbs equation
ΔGads0=ΔHads0-TΔSads0.
Using the obtained values of ΔGads0 from the modified Langmuir isotherm, we can plot ΔGads0 versus T (Figure 8). These plots give straight lines with slopes (−ΔSads0) and intercepts (ΔHads0). The obtained values of the lines parameters are listed in Table 3.
Gibbs plots parameters for copper corrosion in 1 M HNO3 for different concentrations.
Inhibitor C (M)
ΔHads0 (kJ·mol−1)
ΔSads0 (J·mol−1·K−1)
R2
MBI
5·10−5
15.93
159.2
0.999
10−4
12.47
147.6
0.996
5·10−4
9.26
138.3
0.996
10−3
7.82
133.3
0.997
TMBI
5·10−5
13.70
152.6
0.994
10−4
12.05
147.4
0.992
5·10−4
11.53
145.5
0.991
10−3
8.72
137.8
0.995
Gibbs plots for copper corrosion in 1 M HNO3 at different concentrations (a): MBI (b): TMBI.
The values of thermodynamic parameters for the adsorption of the tested molecules can provide valuable information about the mechanism of corrosion inhibition. The endothermic adsorption process (ΔHads0>0) is attributed unequivocally to chemisorption [47], while generally, an exothermic adsorption process (ΔHads0<0) may involve either physisorption or chemisorption or a mixture of both the processes. In the present case, the positive values of ΔHads0 indicate that the adsorption of the tested compounds is an endothermic process. The positive values of ΔSads0 indicate that the adsorption of these molecules is a process accompanied by an increase in entropy change, favoring the adsorption of the molecules onto the surface of the metal.
3.5. Mechanism of Inhibition
As far as the inhibition process is concerned, it is generally assumed that adsorption of an organic inhibitor at the metal/solution interface is the first step in the action mechanism of the organic compounds in aggressive acidic media. Four types of adsorption may take place during inhibition involving organic molecules at the metal/solution interface:
electrostatic attraction between charged molecules and the charged metal,
interaction of unshared electron pairs in the molecule with the metal,
interaction of π-electrons with the metal,
a combination of the above situations [48].
Concerning inhibitors, the inhibition efficiency depends on several factors, such as the number of adsorption sites and their charge density, molecular size, heat of hydrogenation, mode of interaction with the metal surface, and the formation metallic complexes [49].
In our case, the two molecules (MBI and TMBI) can be protonated in HNO3 solution. Thus, they become cations, existing in equilibrium with the corresponding molecular form [50–52]
On the other hand, in 1 M HNO3, copper is corroded in Cu2+, and the surface of copper is positively charged. NO3- ions are attracted by the charges on the copper surface which tends to be charged negatively.
The obtained results could be explained on the assumption that the negatively charged NO3- would attach to the positively charged surface. There may be a synergism between NO3- and the protonated inhibitors [MBIH]+ and [TMBIH]+ near the interface, and the concentrations of NO3- and that of the neutral forms (MBI and TMBI) and the protonated forms of the inhibitors ([MBIH]+ and [TMBIH]+) were probably much higher than those in the bulk solution; the protonated forms did attach electrostatically to the negative charges at the copper surface. When the neutral forms and the protonated forms of the inhibitors adsorb on the copper surface, coordinate bonds are formed by partial transference of electrons from the unprotonated N atoms, sulfur atoms, delocalized π electrons in the benzimidazole rings to the metal surface via vacant d orbitals of Cu2+ ions. So, in the process of adsorption, both physical and chemical adsorptions might take place. Similar mechanism has been proposed by Tang et al. [53]. The inhibitory mechanism on copper dissolution is related to the chelating effect of Cu2+ ions close to copper surface. Indeed, it is proven that copper is corroded to Cu2+ in HNO3 solution and no oxide film is formed to protect the surface from corrosion [37, 38]. It is also known that the organic compounds can be adsorbed by the interactions between the lone pairs of electrons of nitrogen, sulfur, or oxygen atoms with metal surface. These processes are facilitated by the presence of d vacant orbitals of low energy in the copper ions, as observed in transition group metals. Recently, it was found that the formation of donor-acceptor surface complexes between free electrons of an inhibitor and a vacant d orbital of a metal is responsible for the inhibition corrosion process [54]. The number of Cu2+ ions in proximity of copper increases with increasing temperature, leading to the formation of Cu2+-inhibitor complexes protective film (increase in surface coverage) which create a physical barrier between Cu surface and the electrolyte, retarding the dissolution of the metal [36]. The inhibition efficiency is then enhanced with increase in temperature.
In order to support experimental data, theoretical calculations were conducted in order to provide molecular-level understanding of the observed experimental behavior. Among quantum chemical methods for evaluation of corrosion inhibitors, density functional theory, DFT has shown significant promise [55] and appears to be adequate for pointing out the changes in electronic structure responsible for inhibitory action.
There is no doubt that the recent progress in DFT has provided a very useful tool for understanding molecular properties and for describing the behavior of atoms in molecules. DFT methods have become very popular in the last decade due to their accuracy that is similar to other methods in less time and with a smaller investment from the computational point of view. In agreement with the DFT, the energy of the fundamental state of polyelectronic systems can be expressed through the total electronic density, and in fact the use of the electronic density instead of the wave function for the calculation of the energy constitutes the fundamental base of DFT [56]. Several quantum chemical methods and molecular modeling techniques have been performed to correlate the inhibition efficiency of the inhibitors with their molecular properties [16, 41, 57–63]. The reactive ability of the inhibitor is closely linked to their frontier molecular orbitals (FMOs), including highest occupied molecular orbital, HOMO, and lowest unoccupied molecular orbital, LUMO, and the other parameters such as hardness and softness, and so forth. Density functional theory (DFT) has become an attractive theoretical method because it gives exact basic vital parameters for even huge complex molecules at low cost [54, 64, 65]. Furthermore, by using sophisticated computational tools, we can understand reactivity behavior of hard and soft acid-base (HSAB) theory that provides a systematic way for the analysis of the inhibitor/surface interaction [66]. Thus, DFT has become a main source of connecting some traditional empirical concepts with quantum mechanics. Therefore, DFT is a very powerful technique to probe the inhibitor/surface interaction and to analyze experimental data.
For an N-electrons system with total electronic energy (E) and an external potential V(r), chemical potential (μP), known as the negative of electronegativity (χ), has been defined as the first derivative of the total energy with respect to N at V(r) [67]
χ=-μP=-(∂E∂N)V(r).
Hardness (η) has been defined within the DFT theory as the second derivative of the total energy with respect to N at V(r) property which measures both the stability and reactivity of a molecule [68]
η=(∂2E∂N2)V(r)=(∂μP∂N)V(r).
In this formula, E is the electronic energy, N is the number of electrons, V(r) is the external potential due to the nuclei, and μP is the chemical potential.
The number of electrons transferred (ΔN) from the inhibitor molecule to the metal surface can be calculated using the following [69]:
ΔN=χCu-χinh2(ηCu+ηinh),
where χCu and χinh are the absolute electronegativity of copper and inhibitor molecule, respectively; ηCu and ηinh are the absolute hardness of copper and the inhibitor molecule, respectively.
The ionization potential (I) and the electron affinity (A) are related to EHOMO and ELUMO by the following [70]:
I=-EHOMO,A=-ELUMO.
These quantities are related to the absolute electronegativity and the absolute hardness by the equations below [70]
χ=I+A2,η=I-A2.
Global softness can also be defined as [71]
σ=1η=2I-A.
Actually, a new global chemical reactivity parameter has been introduced and is called electrophilicity index (ω). This parameter is defined as [72]
ω=μP22η=χ22η.
This parameter measures the electrophilic power of a molecule.
Let’s consider the situation corresponding to a molecule that is going to receive certain amount of charge at some center and is going to backdonate a certain amount of charge through the same center or another one. To describe the energy charge associated with these two processes, the hardness is fixed to the value of η=(μp+-μP-) in both situations. Thus, according to the simple charge transfer model for donation and backdonation of charges proposed recently [73], when a molecule receives a certain amount of charge ΔN+ associated with the energy changeΔE+ΔE+=μP+ΔN++12η(ΔN+)2.
When a molecule backdonates a certain amount of charge ΔN- associated with the energy changeΔE-, then
ΔE-=μP-ΔN-+12η(ΔN-)2.
If the total energy change is approximated by the sum of the contributions of (22) and (23), and assuming that the amount of charge backdonation is equal to the amount of charge received: ΔN-=-ΔN+, then
ΔET=ΔE++ΔE-=(μP+-μP-)ΔN++η(ΔN+)2.
The total energy change is minimum with respect toΔN+ and implies that ΔN+=-(μP++μP-)/2η and that
ΔET=-(μP+-μP-)24η=-η4.
3.6.2. Computational Details
The aim in this part of the work is to correlate the observed inhibition efficiency with the quantum chemical parameters of the investigated inhibitors. All the calculations were performed by resorting to density-functional theory (DFT) method using the Gaussian 03 W suite of programs. To establish correlation between experimental data and structural and electronic characteristics of the studied compounds, the geometry of the molecules (Figure 1) was optimized by the density functional theory (DFT) method level with B3LYP exchange correlation functional, using 6–31G (d, p) basis set. This basis set provided accurate geometry and electronic properties for a wide range of organic compounds [19, 74].
3.6.3. Correlation between Molecular Orbital Energy Level and Inhibition Efficiency
According to the frontier molecular orbital theory, the chemical reactivity is a function of interaction between HOMO and LUMO levels of the reacting species [62]. EHOMO is a quantum chemical parameter often associated with the electron donating ability of the molecule. High value of EHOMO indicates a tendency of the molecule to donate electrons to appropriate acceptor system of low empty molecular orbital energy [20]. The energy of lowest unoccupied molecular orbital, ELUMO, therefore, indicates the ability of the molecule to accept electrons [71]. So, the lower the value of ELUMO, the more probable the molecule would accept electrons. Thus, the binding ability of the inhibitor to the metal surface increases with increasing of the HOMO and decreasing of the LUMO energy values.
Other indicators are absolute electronegativity (χ), absolute hardness (η), and absolute softness (σ). The absolute electronegativity is the chemical property that describes the ability of a molecule to attract electron towards itself in a covalent bond, while the absolute hardness is measured by the energy gap between the lowest unoccupied and the highest occupied molecular orbitals.
Let’s consider two systems with different electronegativity as a metallic surface and an inhibitor molecule; the following mechanism will take place: the electron flow will happen from the system with the low electronegativity towards that of a higher value until the chemical potentials are the same.
In order to calculate the fraction of electrons transferred, a theoretical value of electronegativity of bulk copper was used χCu = 4.98 eV [75] and a global hardness assuming that for a metallic bulk I = A [76] because they are softer than the neutral metallic atoms. All the calculated quantum chemical parameters are listed in Table 4.
The calculated quantum chemical parameters and descriptors for MBI and TMBI obtained using DFT at B3LYP/6–31(d, p) basis set.
Parameter
MBI
TMBI
EHUMO (eV)
−6.033
−5.856
ELUMO (eV)
−0.639
−0.556
ΔE = ELUMO − EHUMO (eV)
5.394
5.300
μ(D)
3.8428
1.7862
I (eV)
6.033
5.856
A (eV)
0.639
0.556
X (eV)
3.257
3.084
η (eV)
2.783
2.729
ΔN (e)
0.305
0.334
σ (eV−1)
0.359
0.366
ω (eV)
2.063
1.939
ΔET (eV)
−0.696
−0.682
From Table 4, it can be clearly seen that the values of EHOMO are in descent order as TMBI > MBI. This is in agreement with the experimentally determined inhibition efficiencies (Figure 4). For ELUMO, we have the same order TMBI > MBI which follows the order of the inhibition efficiencies. These results show that TMBI has more tendencies to bond with copper by donating electrons.
The energy gap ΔE = ELUMO − EHOMO is an important parameter, and it is 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 in inhibition efficiency of the molecule [77]. The values of ΔE in Table 4 show the following relation: MBI > TMBI, which is in agreement with the inhibition efficiencies obtained for the inhibitors.
The bonding tendencies of the inhibitors towards copper can be discussed in terms of hard-soft-acid-base (HSAB) and the frontier controlled interaction concepts [78, 79]. General rule suggested by the principle of HSAB is that hard acids prefer to coordinate to hard bases, and soft acids prefer to coordinate to soft bases. On the other hand, metal atoms are known as soft acids. Hard molecules have high HOMO-LUMO gap, and soft molecule has a small HOMO-LUMO gap [80], and thus soft bases inhibitors are the most effective for metals [77]. In our case, TMBI, which has the lowest energy gap (ΔE) and the highest softness (σ), has the most inhibition efficiency [81]. We can observe from Table 4 that the values of the hardness (η) of the two molecules are in the following order: TMBI < MBI. This order is the reverse of that obtained for softness. So the inhibitor with the least value of global hardness (hence the highest value of global softness) is the best and vice versa. These results are consistent with the results obtained for experimental efficiencies.
In Table 4, we can also see the calculated values of the number of electrons transferred (ΔN). Values of this parameter show that the inhibition efficiency resulting from electron donation agrees with Lukovits’s study [69]. If ΔN<3.6, the inhibition efficiency increases by increasing electron-donating ability of the inhibitors to the metal surface. The calculated values of ΔN correlate strongly with experimental inhibition efficiencies. Thus, the highest fraction of electrons transferred is associated with the best inhibitor (TMBI), while the least fraction is associated with the inhibitor that has the least inhibition efficiency (MBI).
The electrophilicity index (ω) expresses the ability of the inhibitor molecules to accept electrons; in our case, the values are nearly the same for the studied inhibitors (Table 4); the two inhibitors have nearly the same capacity to accept electrons. The unoccupied d orbitals of Cu2+ ([Ar]3d9) can accept electrons from the inhibitors to form coordinate bond. On the other hand, the inhibitors molecules can accept electrons from Cu2+ ion with its antibonding orbitals to form backdonating bond. These donation and back donation processes strengthen the adsorption of the inhibitors molecules onto copper surface [16].
Signs of hardness (η>0) and that of the total energy change (ΔET<0) (Table 4) show that the charge transfer to the molecule followed by backdonation from the molecule is energetically favorable (charge transfer to the molecule and backdonation from the molecule). Similar observation has been reported in the literature [73]. However, we must note that (25) does not predict that a back donation process is going to occur; it only establishes if both processes occur, the energy change is directly proportional to the hardness of the molecule.
The dipole moment (μ) is another parameter of the electronic distribution in a molecule and is the measure of polarity of a polar covalent bond. According to Khalil [82], lower values of dipole moments (μ) will favour accumulation of the inhibitor in the surface layer and therefore higher inhibition efficiency. In our case TMBI has the lower value of dipole moment and MBI the highest one. So considering Khalil’s point of view, there is a good agreement between the experimental inhibition efficiency values and the obtained dipole moment values. According to Yurt and coworkers, in most cases, the adsorption power (IE (%) of the inhibitor) depends on certain physicochemical parameters (solubility, molecular size, etc.); it increases with lower dipole moments and with increasing nitrogen charge [46]. However, several authors state that the inhibition efficiency increases with increasing values of dipole moment [83, 84]. On the other hand, survey of the literature reveals that several irregularities appeared in case of correlation of dipole moment with inhibitor efficiency [85, 86]. So in general, there is no significant relationship between the dipole moment values and inhibition efficiencies [19].
3.6.4. Correlation between Atomic Charge and Inhibition Efficiency
Frontier orbital energy level indicates the tendency of bonding to the metal surface. Further study on the formation of chelating centers in an inhibitor requires the information of spatial distribution of electronic density of the molecule [87].
The structure of the molecule can affect the adsorption, influencing the electron density at the functional group. Generally, electrophile systems attack the molecule at negatives charged sites. It is proven that the electron density focused on heteroatoms (S, N, O, Se, P,…) which are the active centers [88]; these atoms have the strongest ability of bonding to the metal surface. We have listed in Table 5 the charges of polar atoms which are the centers of the studied compounds, which favor the adsorption of these molecules onto copper surface.
Atomic charges with hydrogen summed into heavy atoms obtained using DFT at B3LYP/6–31(d, p) basis set.
Compound
–N=
–NH–
–S–
MBI
−0.538551
−0.295407
−0.066943
TMBI
−0.651433
−0.294465
−0.257838
From Table 5, it can be seen that charges on nitrogen are nearly similar even if TMBI has the highest global partial charge and MBI the lowest one. We note that great differences exist between charges on S. The absolute values of these partial charges are in descent order as δS(TMBI)>δS(MBI). Inhibition efficiency and electronic charges on S atom have the same trend. The HOMO of TMBI and MBI is made up predominantly of 3px of S atom in each molecule which contribute, respectively, to the electronic charge density of their HOMO levels. S atom of TMBI has the highest contribution, while that of MBI has the lowest one. This can be associated with the dominative contribution of the S atom according to the frontier orbital approximation; the region of highest density is generally the site of electrophilic attack [88].
N and S heteroatoms have lonely electron pairs that are important for bonding unfilled 3d orbitals of copper ion Cu2+; they determine the adsorption of the molecules on the metal surface. TMBI has the highest inhibition efficiency, which result from the geometry change that led to HOMO energy increase and electron density distribution in the molecule. Based on the discussion above, it can be concluded that N and S atoms are the most possible sites of bonding metal surface by donating electrons to copper ions; the inhibition efficiency depends on the total charge (δT) on N and S of each molecule. On the other hand, the efficiency values of the two inhibitors can be explained by the fact that –CH3 attached to benzimidazole ring (TMBI) has the ability to donate more electrons to the ring compared to −H attached to benzimidazole ring (MBI).
4. Conclusion
From the present study, the main following conclusions can be drawn.
The two compounds (TMBI and MBI) are good inhibitors for copper corrosion in 1 M HNO3.
Their inhibition efficiencies increase with increasing concentration and increasing temperature.
The thermodynamic adsorption and activation parameters obtained support both physical and chemical adsorption mechanism.
Adsorption of the two inhibitors on copper surface follows modified Langmuir adsorption model and the thermodynamic/kinetic adsorption model of El-Awady.
The experimental results and the theoretical ones are in good agreement.
Acknowledgment
The authors acknowledge the financial support of the Laboratory of Physical Chemistry of University of Cocody-Abidjan (Côte d’Ivoire).
TrabanelliG.Inhibitors—an old remedy for a new challenge19914764104192-s2.0-0026170137JardyA.Lasalle-MolinA. L.KeddamM.TakenoutiH.Copper dissolution in acidic sulphate media studied by QCM and rrde under ac signal19923712219522012-s2.0-000126674610.1016/0013-4686(92)85111-WHosseiniS. M. A.AzimiA.The inhibition of mild steel corrosion in acidic medium by 1-methyl-3-pyridin-2-yl-thiourea20095147287322-s2.0-6284908774510.1016/j.corsci.2008.11.019OkaforP. C.ZhengY.Synergistic inhibition behaviour of methylbenzyl quaternary imidazoline derivative and iodide ions on mild steel in H2SO4 solutions20095148508592-s2.0-6284908906710.1016/j.corsci.2009.01.027FoudaA. S.EllithyA. S.Inhibition effect of 4-phenylthiazole derivatives on corrosion of 304 L stainless steel in HCl solution20095148688752-s2.0-6284908404510.1016/j.corsci.2009.01.011OtmačićH.Stupnišek-LisacE.Copper corrosion inhibitors in near neutral media20034889859912-s2.0-003742068510.1016/S0013-4686(02)00811-3GašparacR.MartinC. R.Stupnišek-LisacE.In situ studies of imidazole and its derivatives as copper corrosion inhibitors I. Activation energies and thermodynamics of adsorption200014725485512-s2.0-003821044510.1149/1.1393230ZhangD. Q.GaoL. X.ZhouG. D.Inhibition of copper corrosion in aerated hydrochloric acid solution by heterocyclic compounds containing a mercapto group20044612303130402-s2.0-464428466710.1016/j.corsci.2004.04.012BentissF.LagreneeM.TraisnelM.HornezJ. C.The corrosion inhibition of mild steel in acidic media by a new triazole derivative19994147898032-s2.0-003311772910.1016/S0010-938X(98)00153-XMcCaffertyE.PravdicV.ZettlemoyerA. C.Dielectric behaviour of adsorbed water films on the α-Fe2O3 surface197066172017312-s2.0-3704913174010.1039/TF9706601720BentissF.LebriniM.LagrenéeM.Thermodynamic characterization of metal dissolution and inhibitor adsorption processes in mild steel/2,5-bis(n-thienyl)-1,3,4-thiadiazoles/hydrochloric acid system20054712291529312-s2.0-2764444088310.1016/j.corsci.2005.05.034El-Taib HeakalF.HaruyamaS.Impedance studies of the inhibitive effect of benzotriazole on the corrosion of copper in sodium chloride medium19802078878982-s2.0-0019184456EkpeU. J.IbokU. J.ItaB. I.OffiongO. E.EbensoE. E.Inhibitory action of methyl and phenyl thiosemicarbazone derivatives on the corrosion of mild steel in hydrochloric acid199540287932-s2.0-0002841363Metikos-HukovicM.GrubacZ.Stupnisek-LisacE.Organic corrosion inhibitors for aluminum in perchloric acid19945021461512-s2.0-0028377583AdamsonA. W.1990New York, NY, USAJohn Wiley & SonsArslanT.KandemirliF.EbensoE. E.LoveI.AlemuH.Quantum chemical studies on the corrosion inhibition of some sulphonamides on mild steel in acidic medium200951135472-s2.0-5804921180410.1016/j.corsci.2008.10.016JamalizadehE.JafariA. H.HosseiniS. M. A.Semi-empirical and ab initio quantum chemical characterisation of pyridine derivatives as HCl inhibitors of aluminium surface20088701–323302-s2.0-5604910206910.1016/j.theochem.2008.08.033ZhangD. Q.CaiQ. R.GaoL. X.LeeK. Y.Effect of serine, threonine and glutamic acid on the corrosion of copper in aerated hydrochloric acid solution20085012361536212-s2.0-5664910800610.1016/j.corsci.2008.09.007GeceG.The use of quantum chemical methods in corrosion inhibitor studies20085011298129922-s2.0-5414909427710.1016/j.corsci.2008.08.043GeceG.BilgiçS.Quantum chemical study of some cyclic nitrogen compounds as corrosion inhibitors of steel in NaCl media2009518187618782-s2.0-6764979483410.1016/j.corsci.2009.04.003BouklahM.HammoutiB.LagrenéeM.BentissF.Thermodynamic properties of 2,5-bis(4-methoxyphenyl)-1,3,4-oxadiazole as a corrosion inhibitor for mild steel in normal sulfuric acid medium2006489283128422-s2.0-3374716791410.1016/j.corsci.2005.08.019Van AllanJ. A.DeagonB. D.19634New York, NY, USAJohn Wiley & SonsPopovaA.SokolovaE.RaichevaS.ChristovM.AC and DC study of the temperature effect on mild steel corrosion in acid media in the presence of benzimidazole derivatives200345133582-s2.0-003602647010.1016/S0010-938X(02)00072-0El-NaggarM. M.Corrosion inhibition of mild steel in acidic medium by some sulfa drugs compounds2007495222622362-s2.0-3384763072410.1016/j.corsci.2006.10.039LiY.ZhaoP.LiangQ.HouB.Berberine as a natural source inhibitor for mild steel in 1 M H2SO420052525124512532-s2.0-2794447822910.1016/j.apsusc.2005.02.094AbdallahM.Antibacterial drugs as corrosion inhibitors for corrosion of aluminium in hydrochloric solution2004468198119962-s2.0-294262554910.1016/j.corsci.2003.09.031UmorenS. A.EbensoE. E.The synergistic effect of polyacrylamide and iodide ions on the corrosion inhibition of mild steel in H2SO420071062-33873932-s2.0-3534900073210.1016/j.matchemphys.2007.06.018El-EtreA. Y.Inhibition of acid corrosion of aluminum using vanillin2001436103110392-s2.0-003537170110.1016/S0010-938X(00)00127-XLarabiL.BenaliO.MekellecheS. M.HarekY.2-Mercapto-1-methylimidazole as corrosion inhibitor for copper in hydrochloric acid20062533137113782-s2.0-3375121859210.1016/j.apsusc.2006.02.013LebriniM.BentissF.VezinH.LagrenéeM.The inhibition of mild steel corrosion in acidic solutions by 2,5-bis(4-pyridyl)-1,3,4-thiadiazole: structure-activity correlation2006485127912912-s2.0-3364600703910.1016/j.corsci.2005.05.001ObotI. B.Obi-EgbediN. O.UmorenS. A.Antifungal drugs as corrosion inhibitors for aluminium in 0.1 M HCl2009518186818752-s2.0-6764982339610.1016/j.corsci.2009.05.017LebriniM.LagrenéeM.VezinH.GengembreL.BentissF.Electrochemical and quantum chemical studies of new thiadiazole derivatives adsorption on mild steel in normal hydrochloric acid medium20054724855052-s2.0-954423072610.1016/j.corsci.2004.06.001PutilovaI. N.BalezinS. A.BarannikV. P.1960Elmsford, NY, USAPergamon PressElmorsiM. A.HassaneinA. M.Corrosion inhibition of copper by heterocyclic compounds19994112233723522-s2.0-003258655210.1016/S0010-938X(99)00061-XDehriI.ÖzcanM.The effect of temperature on the corrosion of mild steel in acidic media in the presence of some sulphur-containing organic compounds2006982-33163232-s2.0-3374750030410.1016/j.matchemphys.2005.09.020CicileoG. P.RosalesB. M.VarelaF. E.VilcheJ. R.Comparative study of organic inhibitors of copper corrosion1999417135913752-s2.0-003287259710.1016/S0010-938X(98)00190-5PourbaixM.1975Houston, Tex, USANACEJohnsonH. E.LejaJ.On the Potential/pH-Diagrams of Cu–NH3–H2O and Zn–NH3–H2O systems1965112663864110.1149/1.2423629Abd El-RehimS. S.RefaeyS. A. M.TahaF.SalehM. B.AhmedR. A.Corrosion inhibition of mild steel in acidic medium using 2-amino thiophenol and 2-cyanomethyl benzothiazole20013144294352-s2.0-003530364710.1023/A:1017592322277Ashassi-SorkhabiH.ShabaniB.AligholipourB.SeifzadehD.The effect of some Schiff bases on the corrosion of aluminum in hydrochloric acid solution200625212403940472-s2.0-3364619223110.1016/j.apsusc.2005.02.148EbensoE. E.IsabiryeD. A.EddyN. O.Adsorption and quantum chemical studies on the inhibition potentials of some thiosemicarbazides for the corrosion of mild steel in acidic medium2010116247324982-s2.0-7795481317110.3390/ijms11062473HosseiniM.MertensS. F. L.GhorbaniM.ArshadiM. R.Asymmetrical Schiff bases as inhibitors of mild steel corrosion in sulphuric acid media20037838008082-s2.0-003746949610.1016/S0254-0584(02)00390-5AliSk. A.SaeedM. T.RahmanS. U.The isoxazolidines: a new class of corrosion inhibitors of mild steel in acidic medium20034522532662-s2.0-003646262310.1016/S0010-938X(02)00099-9TalatiJ. D.GandhiD. K.N-heterocyclic compounds as corrosion inhibitors for aluminium-copper alloy in hydrochloric acid19832312131513322-s2.0-0020907278Szklarska-SmialowskaZ.MankowskiJ.Crevice corrosion of stainless steels in sodium chloride solution197818119539602-s2.0-0018056042YurtA.UlutasS.DalH.Electrochemical and theoretical investigation on the corrosion of aluminium in acidic solution containing some Schiff bases200625329199252-s2.0-3384540695310.1016/j.apsusc.2006.01.026DurnieW.De MarcoR.JeffersonA.KinsellaB.Development of a structure-activity relationship for oil field corrosion inhibitors19991465175117562-s2.0-003265513510.1149/1.1391837SchweinsbergD. P.GeorgeG. A.NanayakkaraA. K.SteinertD. A.The protective action of epoxy resins and curing agents-inhibitive effects on the aqueous acid corrosion of iron and steel198828133422-s2.0-0023844520FoudaA. S.MoussaM. N.TahaF. I.ElneanaaA. I.The role of some thiosemicarbazide derivatives in the corrosion inhibition of aluminium in hydrochloric acid19862697197262-s2.0-0000956539IvanovD.1964Sofia, BulgariaNauka i izkustvoHauptmanZ.GraefeU.RemaneH.1985Sofia, BulgariaNauka i izkustvoDyatkinaS. L.DamaskinB. B.Adsorption behavior of the different forms of dimethylbenzimidazole at the mercury electrode1986221283TangL.MuG.LiuG.The effect of neutral red on the corrosion inhibition of cold rolled steel in 1.0 M hydrochloric acid20034510225122622-s2.0-004206644910.1016/S0010-938X(03)00046-5ZhangD. Q.CaiQ. R.HeX. M.GaoL. X.KimG. S.Corrosion inhibition and adsorption behavior of methionine on copper in HCl and synergistic effect of zinc ions20091142-36126172-s2.0-5884914890110.1016/j.matchemphys.2008.10.007LópezN.IllasF.Ab initio modeling of the metal-support interface: the interaction of Ni, Pd, and Pt on MgO(100)19981028143014362-s2.0-0031998768AndrésJ.BeltranJ.2000Castellón de la Plana, EspañaUniversitat Jaume IAwadM. K.MahgoubF. M.El-iskandaraniM. M.Theoretical studies of the effect of structural parameters on the inhibition efficiencies of mercapto-1,2,4-triazoline derivatives20005311–31051172-s2.0-003470673410.1016/S0166-1280(00)00437-1AwadM. K.Semiempirical investigation of the inhibition efficiency of thiourea derivatives as corrosion inhibitors200456722192252-s2.0-244248199110.1016/j.jelechem.2003.12.028BereketG.ÖretirC.YurtA.Quantum mechanical calculations on some 4-methyl-5-substituted imidazole derivatives as acidic corrosion inhibitor for zinc20015711–31391452-s2.0-003595923210.1016/S0166-1280(01)00552-8GaoG.LiangC.Electrochemical and DFT studies of β-amino-alcohols as corrosion inhibitors for brass20075213455445592-s2.0-3384736949010.1016/j.electacta.2006.12.058JuH.KaiZ. P.LiY.Aminic nitrogen-bearing polydentate Schiff base compounds as corrosion inhibitors for iron in acidic media: a quantum chemical calculation20085038658712-s2.0-3984910702910.1016/j.corsci.2007.10.009MusaA. Y.KadhumA. A. H.MohamadA. B.RahomaA. A. B.MesmariH.Electrochemical and quantum chemical calculations on 4,4-dimethyloxazolidine-2-thione as inhibitor for mild steel corrosion in hydrochloric acid20109691–32332372-s2.0-7794948884510.1016/j.molstruc.2010.02.051ObotI. B.Obi-EgbediN. O.Theoretical study of benzimidazole and its derivatives and their potential activity as corrosion inhibitors20105226576602-s2.0-7264910094810.1016/j.corsci.2009.10.017ÖzcanM.AC impedance measurement of cystine adsorption at mild steel/sulfuric acid interface as corrosion inhibitor20081212165316612-s2.0-5254911794210.1007/s10008-008-0551-1ZorS.KandemirliF.BingulM.Inhibition effects of methionine and tyrosine on corrosion of iron in HCl solution: electrochemical, FTIR, and quantum-chemical study200945146532-s2.0-6024909241310.1134/S2070205109010079CruzJ.PandiyanT.García-OchoaE.A new inhibitor for mild carbon steel: electrochemical and DFT studies200558318162-s2.0-2384455721010.1016/j.jelechem.2005.02.026ParrR. G.DonnellyR. A.LevyM.PalkeW. E.Electronegativity: the density functional viewpoint1978688380138072-s2.0-36749117171ParrR. G.PearsonR. G.Absolute hardness: companion parameter to absolute electronegativity198310526751275162-s2.0-0347291894LukovitsI.KálmánE.ZucchiF.Corrosion inhibitors—correlation between electronic structure and efficiency2001571382-s2.0-13844295135PearsonR. G.Hard and soft acids and bases19638522353335392-s2.0-0345114077LesarA.MiloševI.Density functional study of the corrosion inhibition properties of 1,2,4-triazole and its amino derivatives20094834–61982032-s2.0-7164910636410.1016/j.cplett.2009.10.082ChattarajP. K.SarkarU.RoyD. R.Electrophilicity index20061066206520912-s2.0-3374575626010.1021/cr040109fGómezB.LikhanovaN. V.Domínguez-AguilarM. A.Martínez-PalouR.VelaA.GázquezJ. L.Quantum chemical study of the inhibitive properties of 2-pyridyl-azoles2006110188928893410.1021/jp057143yHenríquez-RománJ. H.Padilla-CamposL.PáezM. A.ZagalJ. H.RubioM. A.RangelC. M.CostamagnaJ.Cárdenas-JirónG.The influence of aniline and its derivatives on the corrosion behaviour of copper in acid solution: a theoretical approach20057571–31710.1016/j.theochem.2005.05.018MichaelsonH. B.The work function of the elements and its periodicity19774811472947332-s2.0-001755684610.1063/1.323539DewarM. J. S.ZoebischE. G.HealyE. F.StewartJ. J. P.Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model198510713390239092-s2.0-084234177110.1021/ja00299a024AwadM. K.MustafaM. R.ElngaM. M. AboComputational simulation of the molecular structure of some triazoles as inhibitors for the corrosion of metal surface20109591–366742-s2.0-7795726271110.1016/j.theochem.2010.08.008PearsonR. G.Absolute electronegativity and hardness: application to inorganic chemistry19882747347402-s2.0-33845279880KochE. C.Acid-base interactions in energetic materials—I. The hard and soft acids and bases (HSAB) principle—insights to reactivity and sensitivity of energetic materials20053015162-s2.0-2544451152210.1002/prep.200400080LiX.DengS.FuH.LiT.Adsorption and inhibition effect of 6-benzylaminopurine on cold rolled steel in 1.0 M HCl20095416408940982-s2.0-6524909262710.1016/j.electacta.2009.02.084HasanovR.SadikoǧluM.BilgiçS.Electrochemical and quantum chemical studies of some Schiff bases on the corrosion of steel in H2SO4 solution20072538391339212-s2.0-3384664164810.1016/j.apsusc.2006.08.025KhalilN.Quantum chemical approach of corrosion inhibition20034818263526402-s2.0-003873054610.1016/S0013-4686(03)00307-4LagrenéeM.MernariB.ChaibiN.TraisnelM.VezinH.BentissF.Investigation of the inhibitive effect of substituted oxadiazoles on the corrosion of mild steel in HCl medium20014359519622-s2.0-003534081210.1016/S0010-938X(00)00076-7QuraishiM. A.SardarR.Hector bases—a new class of heterocyclic corrosion inhibitors for mild steel in acid solutions20033312116311682-s2.0-034498391310.1023/B:JACH.0000003865.08986.fbKhaledK. F.Babić-SamardžijaK.HackermanN.Theoretical study of the structural effects of polymethylene amines on corrosion inhibition of iron in acid solutions20055012251525202-s2.0-1534434131610.1016/j.electacta.2004.10.079BereketG.HürE.ÖretirC.Quantum chemical studies on some imidazole derivatives as corrosion inhibitors for iron in acidic medium20025781–379882-s2.0-003707479810.1016/S0166-1280(01)00684-4FangJ.LiJ.Quantum chemistry study on the relationship between molecular structure and corrosion inhibition efficiency of amides20025931–31791852-s2.0-003718384510.1016/S0166-1280(02)00316-0JuH.KaiZ.-P.LiY.Aminic nitrogen-bearing polydentate Schiff base compounds as corrosion inhibitors for iron in acidic media: a quantum chemical calculation20085038658712-s2.0-3984910702910.1016/j.corsci.2007.10.009