Abstract

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 𝐸HOMO, 𝐸LUMO, energy gap (Δ𝐸), dipole moment (πœ‡), hardness (πœ‚), softness (𝜎), the fractions of electrons transferred (Δ𝑁), electrophilicity index (Ο‰), and the total energy change (Δ𝐸𝑇). 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 (i)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 𝑆 in TMBI.(ii)To correlate the effect of structure parameters of the two compounds through their quantum parameters (𝐸HOMO,𝐸LUMO,Δ𝐸=𝐸HOMOβˆ’πΈHOMO,πœ‡,πœ‚,𝜎,πœ”) to explain their inhibition efficiency.

2. Materials and Methods

2.1. Molecules Structures

Figure 1 gives the structures of 2-mercaptobenzimidazole (MBI) and 2-thiomethylbenzimidazole (TMBI). Their molecular weight is respectively 𝑀 = 150.2 g/moL and 𝑀 = 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.

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%).

TMBI
RMN 1H (acetone): 2.88 (s, 3H, CH3); 7.26–7.65 (2 m, 4H, C6H4); RMN 13C: 13.969 (CH3S); 115–121.4 (C6H4); 151.4 (C=N); SDM (E.I.): M+1 = 165. Analysis for C8H8N2S: C(58.51%), H(4.91%), N(17.06%), S(19.53%).

2.2. Weight-Loss Measurements

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 (π‘Š) was calculated using the following: π‘šπ‘Š=1βˆ’π‘š2,𝑆𝑑(1) where π‘š1 and π‘š2 are, respectively, the weight (in g) before and after immersion in the test solution, 𝑆 the total surface of the sample (in cm2), and 𝑑 the immersion time (in h).

The inhibition efficiency IE(%) was then calculated using the following: IEπ‘Š(%)=0βˆ’π‘Šπ‘Š0βˆ—100,(2) where π‘Š0 and π‘Š are, respectively, the corrosion rates of copper in the absence and presence of the tested compounds.

3. Results and Discussion

3.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.

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 𝐢 = 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].

3.3. Activation Parameters

Arrhenius suggested the famous equation which evaluates the temperature dependence of the rate constant as follows [33]: 𝐸logπ‘Š=logπ΄βˆ’π‘Ž.2.303𝑅𝑇(3) Here, 𝐴 is the frequency factor (in gΒ·cmβˆ’2Β·hβˆ’1), πΈπ‘Ž is the apparent activation energy (in JΒ·molβˆ’1), 𝑅 is the perfect gas constant (𝑅 = 8.314 JΒ·moLβˆ’1Β·Kβˆ’1), and 𝑇 is the absolute temperature (in 𝐾). Equation (3) predicts that the plot of log π‘Š versus 1/𝑇 should be a straight line. The slope of the line is (β€“πΈπ‘Ž/2.303R), and the intercept of the line extrapolated (1/𝑇 = 0) gives log𝐴.

On the other hand, activation thermodynamic parameters including change in enthalpy(Ξ”π»βˆ—π‘Ž), change in free energy (Ξ”πΊβˆ—π‘Ž), and change in entropy (Ξ”π‘†βˆ—π‘Ž) of activation were calculated using the following [34]: Ξ”π»βˆ—π‘Ž=πΈπ‘Žβˆ’π‘…π‘‡,Ξ”πΊβˆ—π‘Žξ‚Έξ‚΅π‘˜=2.303𝑅𝑇logπ΅π‘‡β„Žξ‚Άξ‚Ή,βˆ’logπ‘ŠΞ”π‘†βˆ—π‘Ž=ξ‚΅Ξ”π»βˆ—π‘Žβˆ’Ξ”πΊβˆ—π‘Žπ‘‡ξ‚Ά,(4) where 𝐾𝐡 is the Boltzmann constant (1.38Β·10βˆ’23JΒ·Kβˆ’1), β„Ž is the Planck constant (6.626Β·10βˆ’34 JΒ·s), and 𝑅 is the perfect gas constant. Figure 5 gives the representation of logπ‘Š versus 1/𝑇.

All the activation thermodynamic parameters (πΈπ‘Ž, Ξ”π»βˆ—π‘Ž, Ξ”π‘†βˆ—π‘Ž, and Ξ”πΊβˆ—π‘Ž) are estimated and listed in Table 1.

From Table 1 it can be seen that, in general, values of the apparent activation energy πΈπ‘Ž 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(5) Moreover, Cu+ ions undergo disproportionation according to the equation below: 2Cu+⇄Cu2++Cu(6) 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 Ξ”π‘†βˆ—π‘Ž 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+𝑛H2Oads⟷Orgads+𝑛H2Osol(7) The surface coverage (πœƒ) was calculated as follows: πœƒ=IE(%).100(8) 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 𝐢 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: πΆπœƒ=1𝐾+𝐢,(9) where 𝐢 is the inhibitor concentration, πœƒ is the fraction of the surface covered, and 𝐾 the equilibrium constant of the adsorption process. Figure 6 shows the dependence of 𝐢/πœƒ as a function of the concentration 𝐢 for MBI and TMBI. The obtained plots are linear with correlation coefficients (𝑅2) 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 πΆπœƒ=𝑛𝐾ads+𝑛𝐢.(10) The equilibrium constant of adsorption is related to the standard free energy of adsorption Δ𝐺0adsby the following: 𝐾ads=1ξƒ©βˆ’55.5β‹…expΔ𝐺0adsξƒͺ,𝑅.𝑇(11) where 𝑅 is universal gas constant, 𝑇 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 Δ𝐺0ads (Table 2) were calculated for 𝐢 = 10βˆ’3 M, using (10) and (11).

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βˆ’πœƒ=log𝐾+𝑦log𝐢,(12) where 𝐢 is the concentration of the inhibitors, πœƒ is the degree of surface coverage, 𝐾ads is the equilibrium constant of the adsorption process, and 𝐾ads=𝐾1/𝑦. In this model, the number of active sites 𝑦 is included. Values of 1/𝑦 less than one imply multilayer adsorption, while 1/𝑦 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 𝐾ads and Δ𝐺0ads for 𝐢 = 10βˆ’3 M are listed in Table 2. From the table, it can be seen that values of 1/𝑦 are very closed to the unity.

The negative values of Δ𝐺0ads indicate a spontaneous adsorption process of MBI and TMBI on copper [44]. It is generally accepted that, for the values of Δ𝐺0adsup 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 Δ𝐺0adsin 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.

Δ𝐺0ads is related to the enthalpy change and entropy change of adsorption process, Δ𝐻0ads and Δ𝑆0ads, respectively by Gibbs equation Δ𝐺0ads=Δ𝐻0adsβˆ’π‘‡Ξ”π‘†0ads.(13) Using the obtained values of Δ𝐺0ads from the modified Langmuir isotherm, we can plot Δ𝐺0ads versus 𝑇 (Figure 8). These plots give straight lines with slopes (βˆ’Ξ”π‘†0ads) and intercepts (Δ𝐻0ads). The obtained values of the lines parameters are listed in Table 3.

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 (Δ𝐻0ads>0) is attributed unequivocally to chemisorption [47], while generally, an exothermic adsorption process (Δ𝐻0ads<0) may involve either physisorption or chemisorption or a mixture of both the processes. In the present case, the positive values of Δ𝐻0ads indicate that the adsorption of the tested compounds is an endothermic process. The positive values of Δ𝑆0ads 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:(i) electrostatic attraction between charged molecules and the charged metal,(ii) interaction of unshared electron pairs in the molecule with the metal,(iii) interaction of πœ‹-electrons with the metal,(iv)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] 623754.fig.009(14)

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.

3.6. Quantum Chemistry Studies
3.6.1. Theoretical Details

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 𝑁-electrons system with total electronic energy (𝐸) and an external potential 𝑉(π‘Ÿ), chemical potential (πœ‡π‘ƒ), known as the negative of electronegativity (πœ’), has been defined as the first derivative of the total energy with respect to 𝑁 at 𝑉(π‘Ÿ) [67] πœ’=βˆ’πœ‡π‘ƒξ‚€=βˆ’πœ•πΈξ‚πœ•π‘π‘‰(π‘Ÿ).(15) Hardness (πœ‚) has been defined within the DFT theory as the second derivative of the total energy with respect to 𝑁 at 𝑉(π‘Ÿ) property which measures both the stability and reactivity of a molecule [68] ξ‚΅πœ•πœ‚=2πΈπœ•π‘2𝑉(π‘Ÿ)=ξ‚΅πœ•πœ‡π‘ƒξ‚Άπœ•π‘π‘‰(π‘Ÿ).(16) In this formula, 𝐸 is the electronic energy, 𝑁 is the number of electrons, 𝑉(π‘Ÿ) is the external potential due to the nuclei, and πœ‡π‘ƒ is the chemical potential.

The number of electrons transferred (Δ𝑁) from the inhibitor molecule to the metal surface can be calculated using the following [69]: πœ’Ξ”π‘=Cuβˆ’πœ’inh2ξ€·πœ‚Cu+πœ‚inhξ€Έ,(17) 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 (𝐼) and the electron affinity (𝐴) are related to 𝐸HOMO and 𝐸LUMO by the following [70]: 𝐼=βˆ’πΈHOMO,𝐴=βˆ’πΈLUMO.(18) These quantities are related to the absolute electronegativity and the absolute hardness by the equations below [70] πœ’=𝐼+𝐴2,πœ‚=πΌβˆ’π΄2.(19) Global softness can also be defined as [71] 1𝜎=πœ‚=2.πΌβˆ’π΄(20) Actually, a new global chemical reactivity parameter has been introduced and is called electrophilicity index (πœ”). This parameter is defined as [72] πœ‡πœ”=2𝑃=πœ’2πœ‚2.2πœ‚(21) 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 πœ‚=(πœ‡π‘+βˆ’πœ‡π‘ƒβˆ’) 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 Δ𝑁+ associated with the energy changeΔ𝐸+Δ𝐸+=πœ‡π‘ƒ+Δ𝑁++12πœ‚(Δ𝑁+)2.(22) When a molecule backdonates a certain amount of charge Ξ”π‘βˆ’ associated with the energy changeΞ”πΈβˆ’, then Ξ”πΈβˆ’=πœ‡π‘ƒβˆ’Ξ”π‘βˆ’+12πœ‚(Ξ”π‘βˆ’)2.(23) 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: Ξ”π‘βˆ’=βˆ’Ξ”π‘+, then Δ𝐸𝑇=Δ𝐸++Ξ”πΈβˆ’=ξ€·πœ‡π‘ƒ+βˆ’πœ‡π‘ƒβˆ’ξ€ΈΞ”π‘++πœ‚(Δ𝑁+)2.(24) The total energy change is minimum with respect toΔ𝑁+ and implies that Δ𝑁+=βˆ’(πœ‡π‘ƒ++πœ‡π‘ƒβˆ’)/2πœ‚ and that Δ𝐸𝑇=βˆ’(πœ‡π‘ƒ+βˆ’πœ‡π‘ƒβˆ’)2πœ‚4πœ‚=βˆ’4.(25)

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]. 𝐸HOMO is a quantum chemical parameter often associated with the electron donating ability of the molecule. High value of 𝐸HOMO 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, 𝐸LUMO, therefore, indicates the ability of the molecule to accept electrons [71]. So, the lower the value of 𝐸LUMO, 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 𝐼 = 𝐴 [76] because they are softer than the neutral metallic atoms. All the calculated quantum chemical parameters are listed in Table 4.

From Table 4, it can be clearly seen that the values of  𝐸HOMO are in descent order as TMBI > MBI. This is in agreement with the experimentally determined inhibition efficiencies (Figure 4). For 𝐸LUMO, 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 Δ𝐸 = 𝐸LUMO βˆ’ 𝐸HOMO is an important parameter, and it is a function of reactivity of the inhibitor molecule towards the adsorption on metallic surface. As Δ𝐸 decreases, the reactivity of the molecule increases, leading to increase in inhibition efficiency of the molecule [77]. The values of Δ𝐸 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 (Δ𝐸) 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 (Δ𝑁). Values of this parameter show that the inhibition efficiency resulting from electron donation agrees with Lukovits’s study [69]. If Δ𝑁<3.6, the inhibition efficiency increases by increasing electron-donating ability of the inhibitors to the metal surface. The calculated values of  Δ𝑁 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 (Δ𝐸𝑇<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.

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 (𝛿𝑇) 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.(1)The two compounds (TMBI and MBI) are good inhibitors for copper corrosion in 1 M HNO3.(2)Their inhibition efficiencies increase with increasing concentration and increasing temperature.(3)The thermodynamic adsorption and activation parameters obtained support both physical and chemical adsorption mechanism.(4)Adsorption of the two inhibitors on copper surface follows modified Langmuir adsorption model and the thermodynamic/kinetic adsorption model of El-Awady.(5)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).