doi: 10.4152/pea.201801035

Experimental and Theoretical Studies of the Corrosion Inhibition of 4-amino-2-(4-chlorophenyl)-8-(2, 3-dimethoxyphenyl)-6-oxo-2, 6-dihydropyrimido [2, 1-b][1, 3] thiazine-3,7-dicarbonitrile on Carbon Steel in a 1.0 M HCl Solution

H. Serrar^a,* , M. Larouj^b,c , H.L. Gaz^b,c , Z. Benzekri^a , A. Zarguil^a , H. Essebaai^d , S. Boukhris^a , H. Oudda^b , R. Salghi^c , A. Hassikou^a and A. Souizi^a

^a Laboratory of Organic, Organometallic and Theoretical Chemistry, Department of Chemistry, Faculty of Sciences, University Ibn Tofail, 14000 Kenitra, Morocco

^b Laboratory of Separation Processes, Central Post, Department of Chemistry, Faculty of Sciences, University Ibn Tofail, 14000 Kenitra, Morocco

^c Equipe de Genie de l'Environnement et Biotechnologie, ENSA, Universite Ibn Zohr, BP1136 Agadir, Morocco

^d Laboratory of Organic Synthesis and Extraction Processes, Department of Chemistry, Faculty of Sciences, University Ibn Tofail, 14000 Kenitra, Morocco

The corrosion inhibition of carbon steel in a 1.0 M HCl solution, using 4-amino-2-(4 chlorophenyl)-8-(2,3-dimethoxyphenyl)-6-oxo-2,6-dihydropyrimido[2,1b][1,3]thiazine- 3,7-dicarbonitrile (ACMPT) was investigated by weight loss, potentiodynamic polarization, electrochemical impedance spectroscopy (EIS) and quantum chemical calculations. Polarization curves indicate that the studied compound was acting as a mixed inhibitor with predominant cathodic effectiveness. The inhibition efficiency decreased with an increased temperature, and the thermodynamic and activation parameters obtained from this study were discussed. The adsorption behavior of ACMT follows Langmuir's isotherm. In addition, Density Function Theory (DFT) calculations were performed on the studied molecule. The theoretical parameters obtained from this method are in good agreement with the experimental results.

Keywords: pyrimidothiazine, carbon steel, corrosion inhibition, theoretical studies, DFT.

Introduction

Mild steel is widely used in different areas, such as in chemical processing, marine application, petroleum production and refining, and construction [1-4]. Furthermore, hydrochloric acid solutions are widely used for the pickling, cleaning, desaling and etching of mild steel [5-7]. In the other hand, hydrochloric acid is the most difficult to handle of the common acids, from the standpoint of corrosion and construction materials. It is required an extreme care to choose the right material to handle this acid by itself, even in relatively dilute concentrations, or in process solutions containing it in appreciable amounts. This acid is well-known for its corrosive capacity towards the most common metals and alloys [8]. One of the most used methods for the protection against corrosion, especially in acidic solutions, is the use of inhibitors [9]. The development of organic inhibitors containing nitrogen, sulphur and oxygen atoms is of growing interest in the field of corrosion and industrial chemistries, as corrosion poses a serious problem to the iron industry [10, 11]. In the other hand, quantum chemical calculations became widely reputed and used to study reaction mechanisms, and to interpret the experimental results, as well as to resolve chemical ambiguities [12]. The choice of pyrimidothiazine as an inhibitor is based on the good results obtained with this type of heterocyclic compounds in different media [13-16]. In the present work, 4-amino-2-(4-chlorophenyl)-8-(2,3- dimethoxyphenyl)-6-oxo-2,6-dihydropyrimido[2,1-b] [1,3] thiazine-3,7- dicarbonitrile (ACMPT) was investigated as a corrosion inhibitor for carbon steel in 1.0 M HCl, using detailed gravimetric measurements, electrochemical impedance spectroscopy and quantum chemical calculations. The molecular structure of this pyrimidothiazine derivative is shown in Fig. 1.

 

 

Experimental methods

]]> Materials and methods

Materials

The metal used in this study is carbon steel (CS) (Euronorm: C35E carbon steel and US specification: SAE 1035); Table 1 shows its chemical composition.

 

Synthesis

Pyrimidin-4-one (1 mmol) was added to a solution of 2-(2, 3-dimethoxybenzy- lidene) malononitrile (1 mmol) in acetonitrile (20 mL). The mixture was then refluxed for 30 h. The solvent was removed under reduced pressure, and the crude product was treated with a mixture of ether/petroleum ether and the 4- amino-2-(4-chlorophenyl)-8-(2,3-dimethoxyphenyl)-6-oxo-2,6- dihydropyrimido[2,1-b][1,3]thiazine-3,7-dicarbonitrile (ACMPT) precipitate, and it was recrystallized from EtOH. The synthetic approaches adopted to obtain the target compounds (ACMOPT) are shown in Scheme 1.

 

Solution

]]> The aggressive solution of 1.0 M HCl was prepared by the dilution of analytical grade 37% HCl with distilled water. The concentration range of used 4-amino-2- (4-chlorophenyl)-8-(2, 3-dimethoxyphenyl)-6-oxo-2, 6-dihydropyrimido[2, 1b][1, 3] thiazine-3,7-dicarbonitrile (ACMPT) was 10^-6 to 10^-3 M.

Weight loss measurements

Gravimetric measurements were carried out at a definite time interval of 4 h, at room temperature, using an analytical balance (precision ± 0.1 mg). The used carbon steel specimens had a rectangular form (length = 1.6 cm, width = 1.6 cm, thickness = 0.07 cm). Gravimetric experiments were carried out in a double glass cell equipped with a thermostated cooling condenser containing 80 mL of the non-de-aerated test solution. After the immersion period, the steel specimens were withdrawn, carefully rinsed with bidistilled water, ultrasonic cleaned in acetone, dried at room temperature and then weighed. Triplicate experiments were performed in each case, and the mean value of the weight loss was calculated.

Electrochemical measurements

Electrochemical impedance spectroscopy

The electrochemical measurements were carried out with a Voltalab (Tacussel- Radiometer PGZ 100) potentiostat, and controlled by Tacussel corrosion analysis software model (Voltamaster 4) under static conditions. The used corrosion cell had three electrodes. The reference electrode was a saturated calomel electrode (SCE). A platinum electrode was used as auxiliary electrode, with a surface area of 1 cm². The working electrode was carbon steel.

All potentials given in this study were referred to this reference electrode. The working electrode was immersed in a test solution for 30 min, to establish a steady state open circuit potential (Eocp). After measuring the Eocp, the electrochemical measurements were performed. All electrochemical tests have been performed in aerated solutions at 303 K. The ESI experiments were conducted in the frequency range with a high limit of 100 KHz, and a different low limit of 0.1 Hz at open circuit potential, with 10 points per decade, at the rest potential after 30 min of immersion in the acid, by applying 10 mV ac voltage peak-to-peak. Nyquist plots were made from these experiments. The best semicircle could be fit through the data points in the Nyquist plot, using a non- linear least square fit, so as to give the intersections with the X-axis. The inhibition efficiency was calculated from the charge transfer resistance values using the following equation [17]:

 

 

where R⁰_ct and Rⁱ_ct are the charge transfer resistance, respectively in the absence and presence of inhibitor.

]]> Potentiodynamic polarization

The electrochemical behavior of a carbon steel sample in inhibited and uninhibited solutions was studied by recording anodic and cathodic potentiodynamic curves. Measurements were performed in the 1.0 M HCl solution containing different concentrations of the tested inhibitor, by automatically changing the electrode potential from 800 to -100 mV versus the corrosion potential at a scan rate of 2 mV.s^-1. The linear Tafel segments of anodic and cathodic curves were extrapolated to the corrosion potential to obtain corrosion current densities (Icorr). From the obtained polarization curves, the corrosion current (Icorr) was calculated by the curve fitting using the following equation:

 

 

The inhibition efficiency was evaluated from the measured Icorr values using the following relationship [18]:

 

 

where I⁰_corr and Iⁱ_corr are the corrosion current density, respectively in the absence and presence of the inhibitor.

Computational details

Density Functional Theory (DFT) has been recently used [18-22] to describe the interaction between the inhibitor molecule and the surface, as well as the properties of these inhibitors concerning their reactivity. The molecular hand gap was computed as the first vertical electronic excitation energy from the ground state using the time dependent density functional theory (TD-DFT) approach, as implemented in Gaussian 03 [23]. For this purpose, some molecular descriptors, such as HOMO and LUMO energy values, frontier orbital energy gap, molecular dipole moment electronegativity (χ), global hardness (η), softness (S) and the fraction of transferred electrons (ΔN) were calculated using DFT method, and have been used to understand the inhibitors' properties, to help in the explanation of the experimental data obtained for the corrosion process. According to Koopman's theorem [24], the ionization potential (IE) and electron affinity (EA) of the inhibitors are calculated using the following equations:

]]>  

 

 

Thus, the values of electronegativity (χ) and chemical hardness (η) that, according to Pearson, are operational and approximate definitions, can be evaluated using the following relation [25]:

 

 

 

The number of transferred electrons (ΔN) was also calculated depending on the quantum chemical method [26, 27].

 

 

]]> where χFe and χinh donate the absolute electronegativity of iron and of the inhibitor molecule, and ηFe and ηinhi donate the absolute hardness of iron and of the inhibitor molecule. In this study, we use the theoretical values of χFe for the absolute hardness of iron and of the inhibitor molecule, and we use the theoretical value of χFe = 7.0 eV and ηFe = 0 to calculate the number of transferred electrons.

 

Results and discussion

Polarization curves

In the present paper, a potentiodynamic polarization study has been made in order to understand the role of the inhibitor in biasing anodic and cathodic reactions. Fig. 2 shows the anodic and cathodic polarization plots of carbon steel in a 1.0 M HCl solution, in the absence and presence of different concentrations of the ACMPT inhibitor at 298 K, while Table 2 shows electrochemical corrosion parameters, such as corrosion potential (Ecorr), corrosion current density (Icorr) and inhibition efficiency (ηIE%).

 

 

It is noted from Fig. 2 that the polarization curves in the 1.0 M HCl solution, containing different concentrations of the studied pyrimidothiazine, are nearly the same. In addition, cathodic and anodic current densities of the polarization curves are both reduced [28, 29].

The above results indicate that retardation of the samples corrosion processes at cathodic and anodic reaction sites occurred, and this was due to the adsorption of these sites by the inhibitor molecules [30-33].

In addition to what was previously described, an inhibitor can generally be classified as from the cathodic or anodic type, if the shift of corrosion potential in the presence of the inhibitor is more than 85 mV, with respect to the corrosion potential of the blank solution.

In the present study, the maximum displacement was 30 mV, suggesting that ACMPT can be arranged as a mixed type inhibitor, with predominant cathodic effectiveness [36, 37].

Also, the examination of Table 2 concludes that the corrosion current density (Icorr) values of carbon steel, in the presence of the ACMPT inhibitor, are lower than those found without inhibitor.

It is noted that the inhibitory efficiency (IE%) increased with the increase in the inhibitor's concentration, and reached a maximum value of 92.35% for a concentration of 10^-3 M of the studied inhibitor.

Electrochemical impedance spectroscopy measurements

The corrosion behavior of carbon steel in 1.0 M HCl, in the absence and presence of various concentrations of ACMPT, was also investigated by EIS technique. The resultant Nyquist plots are shown in Fig. 3.

 

The values of inhibition efficiency ((ηZ%) were calculated by the equation (1).

 

Nyquist plots contain a semicircle with the center under real axis. The size of the semicircle increases with the inhibitor's concentration, indicating the charge transfer process as the main controlling factor for the corrosion inhibition of carbon steel. As we could see from the plots, the impedance of the inhibited solution has increased with the increase in the inhibitor's concentration.

The experimental results of EIS measurements for the corrosion of carbon steel in 1.0 M HCl, in the absence and presence of inhibitor, are given in Table 3.

 

As it can be observed from Table 3, the charge transfer resistance (Rct) value increased with the increase in the inhibitor's concentration. Inversely, values of the interface's capacitance (Cdl) started decreasing, with the increase in the inhibitor's concentration, which is most probably due to the decrease in local dielectric constant and/or to the increase in thickness of the electrical double layer. This suggests that the inhibitor acts via adsorption at the metal/solution interface [38], and that the decrease in the Cdl values is caused by the gradual replacement of water molecules by the adsorption of the inhibitor molecules onto the electrode surface, which decreases the extent of metal dissolution [39]. A keen observation of these plots has revealed that capacitive loops are not perfect semicircles, which is possibly due to the frequency dispersion, roughness, and inhomogeneity of the metal surface, and to impurities, grain boundaries, and distribution of the surface's active sites. The existence of a single semicircle can be attributed to the single charge transfer reaction during the metal dissolution process [40].

The equivalent circuit that describes the present metal/electrolyte corroding system is the simple Randles model shown in Fig. 4, where Rs, Rct and Cdl are, respectively, the resistance in solution, the charge transfer resistor and the double layer capacitance of the interface.

]]>  

Fig. 5 indicates the variation of ACMPT's corrosion rate as a function of its concentration in 1.0 M HCl at 298 K, from weight loss measurements.

 

This figure reveals that the rate of mild steel corrosion in 1.0 M HCl decreases with an increase in the inhibitor's concentration at 298 K.

The increase in inhibition efficiency with an increase in the concentration of the studied compound can be explained by the increased adsorption of ACMPT on the metal surface (Table 4).

 

In the inhibition process, the first step in the inhibition mechanism is the adsorption at the metal/aggressive solution interface [41]. Given the dependence of inhibition efficiency on the concentration, as represented in Fig. 6, the inhibitor decreases the active center for steel dissolution. The adsorption process is made possible by the presence of heteroatoms such as N, S and O, which are regarded as active adsorption centers.

]]>  

The adsorption of ACMPT on the surface of mild steel can be achieved by the interaction between the pair of single electrons of the heteroatoms or between the electron-rich π electron systems and the metal surface. This, as earlier reported by Umoren and Ebenso [42], may be facilitated by the presence of the vacant d- orbital of iron that constitutes steel, as observed in d group metals or in the transition element. In addition to the molecular form, ACMPT can be present in protonated species in an acidic solution. This protonated form facilitates the adsorption of the compound onto the metal surface by electrostatic interaction between the organic molecules and the metal surface [43].

Effect of temperature

The influence of temperature is also studied by potentiodynamic polarization. We have carried out a study in the range of 298-328 K, using stationary electrochemical measurements in potentiostatic mode.

Figs. 7 and 8 show the polarization curves of steel in a 1.0 M HCl medium, respectively in the absence and presence of 10^-3 M of ACMPT, in the studied temperature range.

 

 

The results deduced from the polarization curves are given in Table 5.

]]>  

As it can be seen, raising the temperature increases both anodic and cathodic reactions of the carbon steel electrode, both in the absence and presence of the organic compound.

The inhibition efficiency of ACMPT reaches the maximum values (92.35%) at 298 K, upon an increased temperature; ηIE% is slightly decreased, and the pyrimidothiazine derivative becomes less effective at 328 K. This behavior again shows the physical nature of ACMPT's adsorption in a HCl solution.

The activation energy, Ea, was calculated from an Arrhenius-type plot (Eq. 10) [44, 45].

 

where E is the apparent activation corrosion energy, T is the absolute temperature, k is the Arrhenius pre-exponential constant and R is the universal gas constant.

The enthalpy of activation, ΔHa, and the entropy of activation, ΔSa, were calculated from equation 11) [44, 45].

 

]]> where h, N, ΔSa and ΔHa are, respectively, the Planck constant, the Avogadro number, the activation entropy and the activation entropy.

Plots of ln (Icorr) vs. 1000/T and ln (Icorr/T) vs. 1000/T gave straight lines, respectively with slopes of -Ea/R and -ΔHa/R. The intercepts were A and [ln (R/Nh)+ (ΔSa/R)], respectively for the Arrhenius and transition state equations. Figs. 9 and 10 represent the data plots of Ln (Icorr) vs. 1000/T and Ln (Icorr/T) vs. 1000/T, in the absence and presence of 10^-3 M ACMPT, as a representative example.

 

 

The calculated values from both methods of the activation energy, Ea, the enthalpy of activation, ΔHa, and the entropy of activation, ΔSa, are shown in Table 6.

 

From this table, it could be clearly seen that the values of Ea and of ΔHa, obtained in the presence of ACMPT, are higher than those obtained in the solution without inhibitor. The positive sign of ΔHa reflects the endothermic nature of carbon steel dissolution process, suggesting that the dissolution of carbon steel is slower in the inhibitor's presence [46]. It is also clear that Ea, in the presence of 10^-3 M of ACMPT, is higher than that obtained in the blank. This observation has been reported to be indicative of physical adsorption mechanism, while lower values of Ea suggest a chemisorption mechanism [47, 48].

The presence of the test compound results in an increase in the ΔSa values. This result indicates that corrosion changes from a more ordered system into a disordered system, with the increase in the inhibitory nature of the studied compounds [49]. In the presence of the inhibitor, the value of ΔSa increases, which is generally interpreted as an increase in disorder, as the reactants are converted to the activated complexes [50].

]]> Adsorption isotherm and standard adsorption of free energy

The adsorption isotherm that describes the adsorptive behavior of organic inhibitors is important, in order to know the mechanism of corrosion inhibition. Basic information on the interaction between the inhibitor molecules and the metal surface can be provided by adsorption isotherms. Several adsorption isotherms were attempted to fit the degree of surface coverage values (θ), including Frumkin, Temkin, Freundlich and Langmuir isotherms. The θ values for various concentrations of inhibitors in acidic media have been evaluated from the polarization measurements. The best fit was obtained in the case of Langmuir isotherm, which assumes that the solid surface contains a fixed number of adsorption sites, and that each site holds one adsorbed species [51]. The plot of Cinh/θ vs. Cinh (Fig. 11) yields a straight line with a correlation coefficient of 0.99998, providing that the adsorption of ACMPT from the 1.0 M HCl solution on the carbon steel surface obeys Langmuir adsorption isotherm.

 

This isotherm can be represented as:

 

 

where Cinh is the molar concentration of the inhibitor and Kads is the equilibrium constant or the adsorption-desorption process.

The value of Kads was found to be 654900.3 M^-1 (Table 7).

]]>  

The relatively high value of the adsorption equilibrium constant reflects the high adsorption ability of ACMPT on the carbon steel surface [52].

Kads is related to the standard free energy of adsorption, ΔGads , by the following equation:

 

 

where R is the gas constant, T is the absolute temperature of the experiment, and the constant value of 55.5 is the concentration of water in the solution in mol L-1. The ΔG⁰ads was calculated as -43.13 KJ mol^-1. The negative value of ΔG⁰ads indicates the spontaneity of the adsorption process and the stability of the adsorbed layer on the carbon steel surface.

It is well known that the values of ΔG⁰ads of the order of -20 KJ mol^-1 or lower indicate a physisorption; those of the order of -40 kJ mol^-1 or higher involve a charge sharing or transfer from the inhibitor molecules to the metal surface, to form a coordinate type of bond (chemisorption) [53-54]. On the other hand, the adsorption phenomenon of an organic molecule is not just considered as a purely physical or chemical adsorption phenomenon [55, 56]. A wide spectrum of conditions, from the dominance of chemisorption to electrostatic effects, arises from other adsorptions experimental data [57]. The value of -43.13 kJ mol^-1 may suggest chemisorption mode.

Computational studies

Computational methods have a strong impact towards the design and development of organic corrosion inhibitors.

Nowadays, density function theory (DFT) has been used to analyze the characteristics of the inhibitor/surface mechanism, and to describe the structural nature of the inhibitor on the corrosion process. Furthermore, DFT is considered to be a very useful technique to probe the inhibitor/surface interaction, as well as to analyze the experimental data [58]. Thus, in our present investigation, DFT method was employed to give some insight into the inhibition action of ACMPT molecule on the carbon steel surface. Quantum chemical parameters, such as EHOMO, ELUMO, energy gap, ΔE (ELUMO - EHOMO) and dipole moment (μ), were obtained, for the neutral molecule to predict their activity towards the metal surface. These quantum chemical parameters were generated after geometric optimization with respect to all nuclear coordinates.

]]>

 

Frontier orbital density distribution is useful in predicting adsorption centers of the ACMPT molecule responsible for the interaction with metal surface atoms. As it can be seen from this figure, HOMO is distributed over the entire molecule, but the LUMO density is mainly localized on the dimethoxyphenyl group. The presence of these adsorption centers can cause a flat orientation of ACMPT molecules on the steel surface; thus, a high degree of surface coverage and inhibition efficiency is expected for ACMPT, from the theoretical point of view. These results suggest that O and S atoms and two N atoms are probably the most favorable reactive sites for the adsorption of ACMPT onto the metal surface.

The calculated molecular parameters are listed in Table 8.

 

EHOMO is often associated with the electron-donating ability of a molecule, and its high value (-4.9593) is likely to indicate a tendency to donate electrons to appropriate low-energy acceptor states. Increasing EHOMO values facilitates adsorption (and therefore inhibition), by influencing the transport process through the adsorbed layer. ELUMO indicates the ability of the molecule to accept electrons; hence, these are the acceptor states. The lower the value (-2.8882) of ELUMO, the more probable it is that the molecule would accept electrons [59]. As for the values of ΔE (ELUMO - EHOMO), lower values (2.0711) of the energy difference, ΔE, will cause a higher inhibition efficiency, because the energy to release an electron from the last occupied orbital will be low [60]. For the dipole moment (μ), a higher value (6.6028) will favor a strong interaction of the inhibitor molecules with the metal surface.

In addition, the number of transferred electrons (ΔN) was also calculated based on the quantum chemical method [61]:

 

]]>  

where χFe and χinh denote, respectively, the absolute electronegativity of iron and the inhibitor molecule; and ηFe and ηinh denote, respectively, the absolute hardness of iron and the inhibitor molecule. These quantities are related to the electron affinity (A) and ionization potential (I).

Absolute electro negativity, (χ), and absolute hardness, (η), of the inhibitor molecules are given by [62]:

 

 

IE and EA are related [63] to EHOMO and ELUMO, as follows:

 

 

 

For an iron atom, a theoretical χ, with the value of 7 eV mol^-1, and η, with the value of 0 eV mol^-1, were used [64] to calculate the number of electrons transferred (ΔN) from the inhibitor to the iron atom. The number of transferred electrons strongly depends on what the actual quantum chemical method has employed for computation. Furthermore, the expression ''number of transferred electrons'' is equal to the wording ''electron-donating ability'', which does not imply that the figures of ΔN actually indicate the number of electrons leaving the donor and entering the acceptor molecule.

]]> The value of electron-donating ability (ΔN) was calculated, and it is given in

Conclusions

From the results of this paper, several conclusions could be listed:

• ACMPT was found to be a good inhibitor for carbon steel in 1.0 M HCl. The inhibition efficiency increases with the inhibitor concentration.

• ACMPT acts as a mixed type inhibitor, with predominant cathodic effectiveness.

• The temperature affects the inhibitory efficiency; this is clearly demonstrated by the decrease in the studied inhibitor effect, with the increase in temperature.

• The adsorption of ACMPT at the carbon steel surface has obeyed the Langmuir's adsorption isotherm.

• The inhibition efficiencies obtained from EIS, Tafel polarization and weight loss methods are in reasonable agreement with each other.

• The calculated parameters obtained from the quantum chemical calculations were correlated with experimental results.

References

1. Saravanamoorthy S, Velmathi S. Prog Org Coat. 2013;76:1527.         [ Links ]

2. Nam DN, Bui VQ, Mathesh M, et al. Corros Sci. 2013;76:257.         [ Links ]

3. Nam DN, Somers A, Mathesh M, et al. Corros Sci. 80;2014:128.

4. Bobina M, Kellenberger A, Millet JP, et al. Corros. Sci. 2013; 69:389.         [ Links ]

5. Bilgic S, Yilmaz H. Mater Chem Phys. 2007;79:5.         [ Links ]

6. Ekpe JU, Okafor CP, Ebenso EE, et al. Bull Electrochem. 2001;17:135.         [ Links ]

7. Odoemelam AS, Eddy ON. J Surf Sci Technol. 2008;24:65.         [ Links ]

8. Noor AE, Al-Moubaraki AH. Int J Electrochem Sci. 2008;3:806.         [ Links ]

9. Anejjar A, Salghi R, Zarrouk A, et al. Res Chem Intermediat. 2013;41:913.         [ Links ]

10. Awad KM, Metwally SM., Soliman AS, et al. J Ind Eng Chem. 2014;20:796.         [ Links ]

11. Soliman AS, Metwally SM, Selim SR, et al. J Ind Eng Chem. 2014;20:4311.         [ Links ]

12. Awe EF, Idris OS, Abdulwahab M, et al. Cogent Chem. 2015;1.         [ Links ]

13. Belayachi M, Serrar H, Zarrok H, et al. Int J Electrochem Sci. 2015;10:3010.         [ Links ]

14. Belayachi M, Serrar H, Assyry AE, et al. Int J Electrochem Sci. 2015;10:3038.         [ Links ]

15. Larouj M, Elaoufir Y, Serrar H, et al. Pharm Lett. 2014;6:324.         [ Links ]

16. Larouj M, Belkhaouda M, Lgaz H, et al. Pharm Lett. 2016;8:114.         [ Links ]

17. Zarrok H, Zarrouk A, Hammouti B, et al. Corros. Sci. 2012;64:243.         [ Links ]

18. Obot BI, Obi-Egbedi ON. Curr Appl Phys. 2011;11:382.         [ Links ]

19. Ma H, Chen S, Liu Z, et al. J Mol Struct THEOCHEM. 2006:774;19.         [ Links ]

20. Henriquez-Roman HJ, Padilla-Campos L, Paez AM, et al. J Mol Struct THEOCHEM. 2005;757:1.         [ Links ]

21. Rodriguez-Valdez ML, Martinez-Villafañe A, Glossman-Mitnik D. J Mol Struct THEOCHEM. 2005;713:65.

22. Feng Y, Chen S, Guo W, et al. J Electroanal Chem. 2007;602;115.         [ Links ]

23. Frisch JM, Trucks WG, Schlegel BH, et al. Gaussian 03, Revision E.01, Gaussian, Inc, Wallingford CT; 2004.         [ Links ]

24. Dewar SJM, Thiel W. J Am Chem Soc. 1977;99:4899.         [ Links ]

25. Pearson GR. Inorg Chem. 1988;27:734.         [ Links ]

26. Sastri SV, Perumareddi RJ. Corrosion. 1997;53:617.         [ Links ]

27. Lukovits I, Kalman E, Zucchi F. Corrosion. 2001;57:3.         [ Links ]

28. Fu J, Zang H, Wang Y, et al. Ind Eng Chem Res. 2012;51:6377.         [ Links ]

29. Qiang Y, Zhang S, Xu S, et al. J Colloid Interface Sci. 2016;472:52.         [ Links ]

30. Abd El-Lateef MH, Abbasov MV, Aliyeva IL, et al. Mater Chem Phys. 2013;142:502.         [ Links ]

31. Abd El-Lateef MH, Abbasov MV, Aliyeva IL, et al. J Korean Chem Soc. 2013;57(1):25.         [ Links ]

32. Abd El-Lateef MH. Corros Sci. 2015;92:104.         [ Links ]

33. Abd El-Lateef MH, Ismael M, Mohamed I. Corros Rev. 2015;33:77.         [ Links ]

34. Marusic K, Curkovic OH, Takenouti H. Electrochim Acta. 2011;56:7491.         [ Links ]

35. Khaled KF, El-Maghraby A. Arabian J Chem. 2014;7.3:319.         [ Links ]

36. El Azzouzi M, Aouniti A, Tighadouin S, et al. J Mol Liq. 2016;221:633.         [ Links ]

37. Elmsellem H, Aouniti A, Khoutoul M, et al. J Chem Pharm Res. 2014;6:1216.         [ Links ]

38. Khaled FK. Electrochimica Acta. 2003;48:2493.         [ Links ]

39. Aljourani J, Raeissi K, Golozar AM. Corros Sci. 2009;51:1836.         [ Links ]

40. Chaitra KT, Mohana SNK, Tandon CH. J Mol Liq. 2015;211:1026.         [ Links ]

41. Niu L, Zhang H, Wei F, et al. Appl Surf Sci. 2005;252:1634.         [ Links ]

42. Umoren AS, Ebenso E. Mater Chem Phys. 2007;106:387.         [ Links ]

43. Feng Y, Siow SK, Teo KW, et al. Corros Sci. 1999;41:829.         [ Links ]

44. Wenner M, Branan N, Stern V, et al. Sci Am Mind. 2008;19:82.         [ Links ]

45. Li P, Lin YJ, Tan LK, et al. Electrochimica Acta. 1997;42:605.         [ Links ]

46. Mu NG, Li X, Li F. Mater Chem Phys. 2004;86:59.         [ Links ]

47. Lebrini M, Robert F, Roos C. Int J Electrochem Sci. 2010;5:1698.         [ Links ]

48. Moretti G, Quartarone G, Tassan A, et al. Electrochimica Acta. 1996;41:1971.         [ Links ]

49. Growcock BF, Lopp VR. Corros Sci.1988;28:397.         [ Links ]

50. Zarrouk A, El Ouali I, Bouachrine M, et al. Res Chem Intermediat. 2012;39:1125.         [ Links ]

51. Ali AS, Saeed MT, Rahman US. Corros Sci. 2003;45:253.         [ Links ]

52. Migahed AM. Mater Chem Phys. 2005;93:48.         [ Links ]

53. Ahamad I, Prasad R, Quraishi AM. Corros Sci. 2010;52:933.         [ Links ]

54. Noor AE, Al-Moubaraki HA. Mater Chem Phys. 2008;110:145.         [ Links ]

55. Ozcan M, Solmaz R, Kardas G, et al. Colloids Surf A. Physicochem Eng Asp. 2008;325:57.         [ Links ]

56. Solmaz R. Corros Sci. 2010;52:3321.         [ Links ]

57. Doner A, Solmaz R, ozcan M, et al. Corros Sci. 2011;53:2902.         [ Links ]

58. Moretti G, Guidi F, Grion G. Corros Sci. 2004;46:387.         [ Links ]

59. Guadalupe JH, Garcia-Ochoa E, Maldonado-Rivas JP, et al. J Electroanal Chem. 2011;655;164.         [ Links ]

60. Gece G. Corros Sci. 2008;50:2981.         [ Links ]

61. Issa MR, Awad KM, Atlam MF. Appl Surf Sci. 2008;255:2433.         [ Links ]

62. Sastri SV, Perumareddi RJ. Corrosion. 1997;53:617.         [ Links ]

63. Pearson GR. Inorg Chem. 1988;27:734.         [ Links ]

64. Bentiss F, Traisnel M, Lagrenee M. Corros Sci. 2000;42:127.         [ Links ]

65. Bentiss F, Traisnel M, Lagrenee M. J Appl Electrochem. 2001;31:41.         [ Links ]

66. Lukovits I, Kalman E, Zucchi F. Corrosion. 2001;57:3-8.         [ Links ]

^*Corresponding author. E-mail address: serrar.doc2015@gmail.com

Received January 5, 2017; accepted May 25, 2017

www.peacta.org