Inhibition Performance of Some Sulfonylurea on Copper Corrosion in Nitric Acid Solution Evaluated Theoretically by DFT Calculations ()
1. Introduction
Corrosion is a chemical degradation of a material and the alteration of its properties by chemical reaction with the surrounding environment. It is a phenomenon that affects many structures by rendering them unusable for their intended application, particularly those made of metallic materials. Indeed, metallic materials, especially copper, are widely used in industry because of its remarkable physical, mechanical, anti-corrosion and biological properties [1]. During its use, copper is exposed to the action of aggressive media [2] [3], such as acid solutions which act as corrosive agents. These acid solutions, which very often contribute to copper degradation, are used in several areas [4] [5], namely: pickling or acid cleaning, oil well stimulation, removal of localized deposits and industrial synthesis processes. For reducing this phenomenon, corrosion inhibitors, especially organic molecules, are used. These first act by adsorption on the surface of metals before intervening in the corrosion reaction process to reduce the speed [6] [7]. Several works [8] - [13] have shown that organic compounds, corrosion inhibitors contain heteroatoms (S, P, O, N) and/or π bonds including all the aromatic rings in their molecular structure which can offer special active electrons or vacant orbitals capable of accepting or giving electrons. In order to preserve the environment, the solution of non-toxic organic inhibitors seems to be the most widely used method [14] [15] [16] in both industry and crafts. So, to inhibit the corrosion of metals, several researches have been directed towards natural products [17] [18] extracted from plant sources and pharmaceutical molecules [19] [20]. Since the adsorption of an organic compound is linked to its molecular structure, the use of quantum chemical methods is more than necessary to achieve the descriptive parameters of the molecules. Many studies in the literature [21] [22] [23] show the correlations between the organic molecules chemical reactivity and their descriptor parameters. Theoretical calculations [24] are today an indispensable tool for the control of molecular structures, the physico-chemical properties of matter and the understanding of metal-molecule interactions. This activity is favoured by the powerful computers and the quality of software (Gaussian, Gamess, Molpro, NWchem, etc.) which have enabled the development of numerical simulation and molecular modelling methods. These calculations [25] - [31] contribute to a better understanding of inhibition properties and reactive behavior of the molecules studied. So, they give access to the molecules parameters such as highest occupied molecular orbital energy (EHOMO), lowest unoccupied molecular orbital energy (ELUMO), energy gap (∆E), dipole moment (μ), electronegativity χ, hardness η, softness S, electrophylicity index ω . Local parameters (Fukui function
or
and dual descriptor (
or
) allowed to access the nucleophilic and electrophilic attack sites of the molecules in order to better understand the mechanism of corrosion inhibition by organic molecules.
The objective of this work is to study theoretically the inhibition properties of three antidiabetic drugs of the sufonylureas family which are chlorpropamide, tolazamide, and glipizide. The choice of these molecules is based on the existence of heteroatoms such as nitrogen, oxygen and sulfur; π bonds and aromatic rings which are important factors to be good inhibitors and they are not toxic at low concentration.
2. Materials and Methods
2.1. Studied Molecules
The molecular structure of each compound is given by Figure 1.
2.2. DFT Calculations
In this work, the quantum chemical calculations were performed with Gaussian 09 W software [32]. Density Functional Theory (DFT) is an application incorporated in the commercial software (Gaussian). This application has permitted to optimize the geometry of each molecule at B3LYP level, with 6-31G(d) basis set. Indeed B3LYP with Becke’s three parameter Lee-Yang Parr hybrid functional, provides good results which allow to describe with precision the behavior of organic compounds [33] [34]. These calculations provided to access to the global and local descriptor parameters, which will facilitate the understanding of metal-molecule interactions. The optimized structure of the studied compounds is given in Figure 2.
The determination of global reactivity parameters helps to explain the inhibition properties of organic molecules.
The electronegativity is related to the chemical potential by the following equation [35].
(1)
It is a global property of the molecular system. The chemical potential is equal to the slope of the total energy as a function of the number of electrons N at constant external potential ν(r).
According to Koopman’s theorem [36] the negative value of EHOMO is defined by the ionization potential (I) and is given by:
(2)
Figure 1. Chemical structure of the inhibitors investigated.
Chlorpropamide (CPA) Tolazamide(TZA) Glipizide (GP)
Figure 2. Optimized structure of CPA, TZA and GP with B3LYP/6-31G(d) method.
The lowest unoccupied molecular orbital energy ELUMO is another global parameter whose negative value is defined by the electronic affinity (A) according to the following expression:
(3)
Using finite difference approximation, μP chemical potential can be expressed as a function of the ionization potential I and the electron affinity A:
(4)
Hardness η and the softness S [37] [38] are obtained from the first derivative of chemical potential or the second derivative of the total energy:
(5)
These fundamental quantities in acids and bases theory developed by Pearson [39] have been used to interpret chemical reactions results. They can also be written as a function of ionization potential (I) and electron affinity (A). Softness reflects the capacity of an atom or molecule to retain an acquired charge [40] [41]. The hardness of an atom or molecule is the energy required for its dismutation; it is the resistance to charge transfer.
(6)
(7)
The electrophilicity index [42] is used to characterize the ability of a molecule to donate or accept electrons. It can then be expressed by:
(8)
The fraction of electrons transferred (ΔN) from the inhibitor molecule to the metal was calculated according to Pearson’s electronegativity relationship [43]:
(9)
where χCu and ηCu, χinh and ηinh denote the electronegativity and hardness of copper and the inhibitor molecule respectively. In this study, we use the theoretical value of χCu = 4.98 eV [44] and ηCu = 0, assuming that for a metallic charge I = A [45] because they are softer than the neutral metallic atoms.
The corrosion inhibition of a metal is influenced by the local reactivity based on Fukui functions. Their values are used to identify the atoms in the inhibitor that can undergo either electrophilic or nucleophilic attack.
Using finite difference approximation, the Fukui functions
and
are respectively equal to the electron densities of LUMO and HOMO orbitals [46] and are defined as follows:
Nucleophilic attack:
(10)
Electrophilic attack:
(11)
where
is the electron density at a point r in space around the molecule, N corresponds to the number of electrons in the neutral molecule, N + 1 corresponds to an anion with an electron added to LUMO of the neutral molecule and N − 1 corresponds to a cation with an electron removed from HOMO of the neutral molecule.
These functions are condensed on atom k where the electron density is replaced by an electronic population qk. In this case the previous expressions become [47]:
Nucleophilic Attack:
(12)
Electrophilic attack:
(13)
where
,
and
are the electronic population of atom k in (N + 1), N and (N − 1) of electrons systems.
A dual descriptor [31] has been introduced recently to determine the individual sites within the molecule with particular behavior. It is given by the following equation:
(14)
The condensed form of the dual descriptor is given by the following relation:
(15)
when
, the process is driven by a nucleophilic attack, and the atom k acts as an electrophile, for
the process is driven by an electrophilic attack and the atom k acts as a nucleophile.
3. Results and Discussion
3.1. Global Reactivity
The different formulas listed above have been used to determine the reactivity descriptors. DFT calculations permit to characterize the reactivity properties of chemical compounds in order to predict the correlation between their molecular structure and their possible behavior like corrosion inhibitors. The different values of these quantum chemical parameters of the molecules are recorded in Table 1.
The highest occupied molecular orbital energy (EHOMO) is a reactivity parameter of molecules associated with the ability to provide electrons. A high value of this parameter indicates a good tendency to donate electron to an appropriate acceptor with low empty energy orbital (metal) [48]. Indeed a metal has a high tendency to accept electrons from an electron donor (organic molecule) into its lowest unoccupied orbital, which could create an electron layer on the metal surface. In our case the high EHOMO values of the three molecules show their high ability to give electrons to copper and this justifies their good performance in inhibiting copper corrosion in nitric acid solution. The EHOMO values of the molecules studied increase in the following order: GP > TZA > CPA, which means that GP could have high inhibition efficiency.
Table 1. Quantum chemical parameters of the studied compounds calculated using B3LYP/6-31G(d).
The ability of a molecule to form bonds with a metal surface [49] also depends on the ELUMO (lowest unoccupied molecular orbital energy) value. Indeed the binding ability of the inhibitor to the metal surface increases with increasing of HOMO and decreasing of LUMO energy values. In general organic compounds do not only donate electron to the unoccupied d orbital of the metal ion such as Cu2+([Ar]3d9), but can also accept electron from the d orbital of the metal, leading to the formation of a feedback bond. According to the literature excellent corrosion inhibitors are usually those organic compounds which not only offer electrons to unoccupied orbital of the metal, but also accept free electrons from the metal [50]. In this work, the low ELUMO values of the studied molecule show that they have a tendency to accept electrons from the metal copper. In Table 1, GP has the lowest value of ELUMO, so it could be the best inhibitor. This transfer or sharing of electrons between the molecules and the unsaturated “d” orbitals of the metal surface allows the formation of covalent bonds (important bonds) reflecting the chemical adsorption.
The value of the energy gap ΔE = ELUMO − EHOMO is a parameter that gives information about the reactivity of a chemical species. Indeed a high value of ΔE means that the molecule is less reactive [51]. When ΔE value of a molecule is low, this favors the exchanges of electrons between this molecule and the metal because it is easier to remove an electron from HOMO orbital to LUMO. We observe that ΔE value of the studied compound are low, which justifies their good inhibition performance. Comparing the different values of ΔE in Table 1, GP has the lowest value therefore it could be the good inhibitor. This low value of ΔE for GP is justified by its large molecular structure with more heteroatoms (O, N and S) than the other two compounds. The value of ΔE for CPA is the highest, which could be justified by the presence of chlorine in this molecule. The HOMO-LUMO diagram in Figure 3 shows an approximation of the frontier molecular orbitals energy gap to be franked in order to remove an electron from HOMO toward LUMO for each molecule.
The dipole moment (μ) measures the polarity of a bond and is related to the distribution of electronic charges in the molecule and also reflects the ability of the molecule to adsorb to the metal surface. According to some authors [52] [53] the higher dipole moment of a molecule, the greater its ability to adsorb on the surface of a metal. However, many other authors [54] [55] state that low dipole moment values favour the adsorption process. In our case CPA and GP have larger dipole moments than TZA. Thus, taking into account the divergent views, there is no significant relationship between dipole moment and inhibition efficiency [56].
Ionization energy (I) and electron affinity (A), which are respectively associated with the HOMO and LUMO energies, are fundamental descriptors of the chemical reactivity of a molecule. A high ionization energy indicates that the molecule is stable and inert to any chemical reaction, while a low ionization energy indicates that the molecule is reactive [57]. The low ionization energy of
Figure 3. Energy diagram LUMO-HOMO of GP, TZA and CPA.
the three molecules indicates their high inhibition efficiency. The lowest value of GP confirms its greater inhibition efficiency.
The electronegativity (χ) of a molecule reflects its ability to attract electrons. The electronegativity values of the three molecules are lower than copper (4.98), which means that copper has the best capacity of attraction. This information indicates that there is a possible movement of electrons from each molecule towards copper.
The probable ability of a molecule to interact with a metal surface is translated by the global softness (S) and global hardness (η). A good inhibitor has a high softness value and a low hardness value [58]. The different molecules in our case have low hardness values and high softness values. So, these values show that the three (3) molecules studied are good inhibitors and GP which has the lowest hardness value (1.6424 eV) and the highest softness value (0.6089 (eV)−1) would have the best inhibition efficiency.
The Fraction of electron transferred (ΔN) of a molecule reflects its ability to donate electrons. According to Lukovits’s study, if ΔN < 3.6 [59] it means that the inhibition efficiency increases with the ability of molecules to donate electrons. In our study, the fractions of electrons transferred from molecules are below 3.6 and are in the order GP > TZA > CPA; this reflects that the good inhibition performance of these molecules is related to the transfer of electrons from these compound to the metal. GP has the best inhibition performance compared to the other two molecules.
The electrophilicity index (ω) measures the propensity of chemical species to accept electrons. A high value of ω [42] describes a good electrophile, while a low value of ω describes a good nucleophile. In our case the electrophilicity index values of the different molecules are high, which expresses that the molecule have a good capacity to accept electrons from copper.
The total energy is an indicator that provides information on metal-molecule interaction. The total energy of a system is the sum of the internal, potential and kinetic energy. Hohenberg and Kohn [60] have shown that the total energy of a system, including that of the many bodily effects of electrons (exchange and correlation) in the presence of a static external potential (atomic nuclei), is a unique functional of charge density. The exact electron density is the one that minimizes the functional energy. In our work, the total energy of the three molecules is less than zero (ET < 0) and η > 0, the charge transfer from each molecule to the metal is energetically favorable [61]. In this context, there is an interaction between the molecules and the metal surface.
3.2. Local Reactivity
Information on the local reactivity of the molecules is provided by the Fukui functions and the dual descriptor. Their values allow us to determine the probable electrophilic and nucleophilic attacks sites.
The atom having the highest value of
or
, represents the most likely centre for nucleophilic attack (electron receiving centre), this centre is associated with the LUMO orbital. The atom whose value of
is the highest or
, represents the most probable site of electrophilic attack (electron emission centre), this centre is associated with the HOMO orbital. The dual descriptor is defined in the interval
[31], which makes it easier to characterize the reactivity of a molecule. The different values are recorded in Tables 2-4.
According to Martínez-Araya’s study the dual descriptor is more accurate local reactivity descriptor than Fukui function [62]. Although the Fukui function has the ability to reveal nucleophilic and electrophilic sites in a molecule, the dual descriptor is able to unambiguously specify true sites for nucleophilic and electrophilic attacks; furthermore, the dual descriptor is less affected by the lack of relaxation terms than the Fukui function [62]. In this context we used the Fukui function and the dual descriptor to better explain the reactivity of the molecules.
In Table 2, the atom S (12) has the highest value of
and the positive value of
so, S (12) is the nucleophilic attack centre for CPA. On the other hand, O (19) has the highest value of
and does not have the negative value of
. Therefore the atom O (20) which has the negative value of
is the electrophilic attack centre for CPA. The reactivity parameters of TZA depicted in Table 3 implies that S (15) has the highest value of
and
so, it is the nucleophilic attack centre. The electrophilic attack centre
Table 2. Mulliken atomic charges, Fukui functions and dual descriptor by B3LYP/6-31G(d) of CPA.
Table 3. Mulliken atomic charges, Fukui functions and dual descriptor by B3LYP/6-31G(d) of TZA.
Table 4. Mulliken atomic charges, Fukui functions and dual descriptor by B3LYP/6-31G(d) of GP.
is controlled by O (16) because it has the highest value of
and the negative value of
. In that concerns GP listed in Table 4, the nucleophilic attack centre is S (27) because it has the largest value of
and the positive value of
while O (28) having the largest value of
and
is the electrophilic attack centre.
We note that for these three sulfonylureas molecules, the sulfur atom present in each compound is the most probable site for nucleophilic attack and the oxygen atom that shares a double bond with that sulfur atom is the most probable site for electrophilic attack. The different information permits to explain the copper corrosion inhibition mechanism in HNO3 by sulfonylureas.
4. Copper Corrosion Inhibition Mechanism
The inhibition of copper corrosion in nitric acid solution by sulfonylureas in particular chlorpropamide, tolazamide and glipizide, is favoured by a protective layer creation on the metal surface. This protective layer is due to the double transfer of electrons between the two entities: molecule copper and copper molecule; justified by ΔN and ω values respectively. These electron transfers justify the existence of the chemical adsorption. Moreover, in nitric acid solution some molecules of tolazamide and glipizide can be protonated because they contain heteroatoms:
(16)
Taking into the chlorine atom present in chloropropamide, it will be transformed in cationic form by releasing chloric ion and some of the cationic species could be protonated because they contain also heteroatoms:
(17)
(18)
There is an interaction between protonated species and NO− 3 ions adsorbed on the metal surface. In order to explain this phenomenon, a schematic mechanism has been proposed in Figure 4.
Figure 4. Schematic mechanism of copper corrosion inhibition by sulfonylureas.
5. Conclusion
The analysis of global descriptor parameters (ELUMO-EHOMO, η, S, ΔN, ω, …) of chloropropamide, tolazamide and glipizide has shown their inhibition properties. It generally appears that these molecules are able to give electrons to the metal, so they can inhibit copper corrosion in nitric acid solution. These descriptor parameters indicate that their inhibition efficiency increases in the following order: glipizide > tolazamide > chloropropamide. The determination of the local parameters (
,
and
of the studied molecules permitted to specify the centre of electrophilic and nucleophilic attack.
Acknowledgements
The authors gratefully acknowledged the support of Environmental Training and Research Unit of Daloa (Côte d’Ivoire) and the Laboratory of physical chemistry of Felix Houphouët Boigny university of Abidjan (Côte d’Ivoire).