Corrosion Study of Mild Steel in Aqueous Sulfuric Acid Solution Using 4-Methyl-4H-1,2,4-Triazole-3-Thiol and 2-Mercaptonicotinic Acid—An Experimental and Theoretical Study

The corrosion behavior of mild steel in 0.1 M aqueous sulfuric acid medium has been studied using weight loss, potentiodynamic polarization measurements, quantum chemical calculations, and molecular dynamic simulations in the presence and absence of 4-methyl-4H-1,2,4-triazole-3-thiol and 2-mercaptonicotinic acid. Potentiodynamic measurements indicate that these compounds mostly act as mixed inhibitors due to their adsorption on the mild steel surface. The goal of the study was to use theoretical calculations to better understand the inhibition. Monte Carlo simulation was used to calculate the adsorption behavior of the studied molecules onto Fe (1 1 1) and Fe₂O₃ (1 1 1) surface. The molecules were also studied with the density functional theory (DFT), using the B3LYP functional in order to determine the relationship between the molecular structure and the corrosion inhibition behavior. More accurate adsorption energies between the studied molecules and iron or iron oxide were calculated by using DFT with periodic boundary conditions. The calculated theoretical parameters gave important assistance into the understanding the corrosion inhibition mechanism expressed by the molecules and are in full agreement with the experimental results.

Introduction

Mild steel is a valuable construction material used in myriad diverse industries mostly for to its outstanding mechanical properties and its almost insignificant cost compared to other materials (Al-Amiery et al., 2014; Su et al., 2016). However, mild steel as many other metals is prone to corrosion, hence their surface must be protected from this undesired process (Migahed et al., 2004). The protection of metals, apart from the use of classic inhibitors (Sanyal, 1981; Selvi et al., 2003; Fouda and Ellithy, 2009; Obi-Egbedi et al., 2011b; Finšgar and Jackson, 2014; Berisha et al., 2015b; Mohsenifar et al., 2016), might be achieved by chemical or electrochemical surface modification such as, SAM's (Self-Assembled Monolayers) formed from silanes (Palanivel et al., 2005; Van Ooij et al., 2005), phosphonic acids (Abohalkuma and Telegdi, 2015; Kosian et al., 2016), sodium oleate (Shubha et al., 2013), or electrochemical reduction of aryldiazonium salts on metals (Chaussé et al., 2002; Berisha et al., 2011, 2015a). Nearly all of the important acid inhibitors are organic molecules containing heteroatoms like nitrogen, oxygen, phosphorous, sulfur, etc. It has been claimed that the corrosion protection performance exhibited by these molecules decreases following the order: O > N > S > P (Berisha et al., 2015b; Mohsenifar et al., 2016). As triazole derivatives represent a widely known class of molecules that exhibit very good corrosion inhibition properties toward mild steel we selected one of such molecule for our study (Bentiss et al., 1999; Selvi et al., 2003; Döner et al., 2011). Density functional theory has become a convenient method to decipher experimental results, permitting to obtain reliable structural parameters for molecules (Geerlings et al., 2003; Gece, 2008; Berisha et al., 2015b). In corrosion studies, this method makes it possible to accurately predict the inhibition efficiency of organic corrosion inhibitors on the basis of electronic and molecular properties as well as reactivity indexes (Obot et al., 2015). In this study, two different mercapto compounds: (a) 2-mercaptonicotinic acid (2MA) and (b) 4-methyl-4H-1,2,4-triazole-3-thiol (4MT) were used as corrosion inhibitors of mild steel sulfuric acid solution (c = 0.1 M). The adsorption mechanism and inhibition performance of 4MT and 2MA molecules (in the neutral and protonated forms) acid were examined as corrosion inhibitors by means of DFT at the B3LYP/6-31G (d,p) basis set level. Moreover, molecular dynamics simulations where used to calculate the adsorption geometries of the adsorbate and DFT with a plane wave basis set calculations to evaluate more precisely adsorption energies and the interaction with the iron or iron oxide surface.

GRAPHICAL ABSTRACT

Graphical Abstract. Adsorption of 2MA and 4MT acts as proton barrier that decreases the corrosion rate of mild steel in the acidic aqueous media.

Experimental

For the electrochemical measurements, the electrode was prepared by embedding a mild steel wire (d = 2 mm, l = 10 mm) inside a Teflon^® (d = 1 cm, l = 6 cm) tube with epoxy resin. Prior to its use, the electrode was polished on silicon carbide abrasive paper (medium grain diameter 6.5–15.3 microns), then on a (DP-Nap) cloth with an aluminum oxide (0.3 micron particle size) suspension, then the electrode was washed and sonicated in water. The chemical composition of the electrode and the coupons was as follows: iron 99.5494%, carbon 0.1252%, phosphorous 0.0316%, manganese 0.1831%, silicon 0.0561%, chromium 0.0124%, sulfur 0.0282%, molybdenum 0.0125%, and nickel 0.0015%.

Electrochemical Test

Electrochemical studies: A PalmSens3 potentiostat was used along with a three-electrode cell at 298K. A graphite rod (d = 3 mm, l = 4 cm) served as an auxiliary electrode and the saturated calomel electrode (SCE) as a reference electrode.

The potentiodynamic polarization curves were obtained by scanning the electrode potential at least 250 mV vs. E_OCP with a sweep rate of 1 mVs⁻¹. The measurements were conducted under atmospheric conditions. To check the reproducibility, every experiment was repeated three times.

Weight Loss Measurements

The weight loss tests (repeated three times) were performed at 298K using 100 ml of the aerated corrosion solution (0.1 M H₂SO₄) in the absence or presence of the studied inhibitors (Obi-Egbedi et al., 2011a). Prior to immersion of mild steel coupons [Size (W × L × D): 12.7 mm × 50.8 mm × 6 mm], the mirror like coupons (abraded with emery papers of various grade sizes: 400, 600, 1,000, and 1,500) were rinsed with double distilled water, cleaned in a sonicating acetone bath for 15 min, followed by sonication in ethanol. The weight difference (weighed using Scaltec Analytical Balance model SBC 31) between the mild steel coupons weight at 6 h time and the initial weight of the coupons was taken as the weight loss which was used to calculate the corrosion rate given by:

$\begin{array}{ll}\rho =\frac{\Delta W}{\Delta t}& (1)\end{array}$

where:

ρ is the corrosion rate, Δw is the weight loss, A is the area of the coupon and t is the corrosion time.

Computational Details

Molecular Dynamics Simulation

Adsorption Locator module in Materials studio 7.0 has been used to build 4MT, 2MA molecules, Fe(1 1 1) (Khaled et al., 2012) and Fe₂O₃ (1 1 1) (Bowker et al., 2012) surface (Akkermans et al., 2013). The molecular dynamic simulations of the interaction between the studied inhibitors and the two surfaces were carried out in the simulation box [Fe(111)–16.21Å × 16.21Å×7.46Å; Fe₂O₃(111)–14.39Å × 13.46Å × 5.88Å) using periodic boundary conditions with a of 20 Å vacuum along the C-axis. The solvent (water) effect was simulated by loading 50 water molecules (geometrically optimized using COMPASS forcefield) onto the simulation box together with the studied molecules (using the same optimization algorithm). For the charged structures (protonated forms of triazole and pyridine rings; Supplementary Figure 1), a positive charge was applied on the protonated nitrogens in both of the structures (using QEq method).

The Metropolis Monte Carlo method was used to determine the adsorption configurations (COMPASS, force field) of the interaction between the 4MT, 2MA, and the Fe (1 1 1) or Fe₂O₃ (1 1 1) surface (Khaled et al., 2012; Yesudass et al., 2016).

Quantum Chemical Calculations

The quantum chemical calculations (Obi-Egbedi et al., 2011b; Obot et al., 2015) were performed by using DFT: B3LYP/6-31G* as implemented in the Spartan 16 software (Wavefunction Inc., Irvine, CA). All energy minima were characterized by performing a vibrational analysis to ensure the lack of imaginary frequencies (Berisha et al., 2017). The DFT: B3LYP/6-31G* permits calculating, in a reproducible manner, the molecular parameters (electronegativity, global hardness, and softness, electron affinity and ionization potential, etc.; Geerlings et al., 2003; Gece, 2008; Berisha et al., 2015b; Obot et al., 2015).

DFT molecular parameters, their meaning and the equations used for their calculations are presented in the Table 1.

TABLE 1

Table 1. Molecular DFT parameters definitions.

DFT Calculations with a Plane Wave Basis Set

Geometry optimization with CASTEP (Clark et al., 2005) is done through generalized gradient approximations (GGE) using PBE functional (Hammer et al., 1999). Self-consistent iteration method (SCF; Payne et al., 1992) was used for geometry optimization (1,000 iteration steps, using energy convergence of 2.0e-5 eV/atom). A Fe cluster of 3.25 Å (in a vacuum box of 10 ų) and Fe₂O₃ cluster of 3.5 Å in diameter (in a vacuum box of 15 ų) were used for the calculations. The same energy optimization parameters and vacuum box sizes were used also to compute energies of the adsorbate molecules and the DOS (Density of States) plots.

The adsorption energy (E_ads.) (Gao et al., 2015) was calculated as:

$\begin{array}{ll}E(ads)=-(Etotal+Ecluster-Eadsorbate)& (2)\end{array}$

where:

E_total is the energy of the Fe or Fe₂O₃ cluster/ adsorbate, E_adsorbate is the energy of the isolated adsorbate molecules (2MA or 4MA) and E_cluster is the energy of the isolated, fully relaxed Fe or Fe₂O₃ clusters.

Results and Discussion

In Figures 1, 2, the anodic and cathodic polarization curves of a mild steel electrode in H₂SO₄ solution (c = 0.1 M) are presented in the absence and presence of 1 × 10⁻³ M, 1 × 10⁻⁴ M of 2MA and 4MT, respectively, at 298 K. The IE (%) was calculated using the following Equation (9):

$\begin{array}{ll}IE(\%)=\frac{\left[i_{absenceofinhibitor}^{corr.}-i_{presenceofinhibitor}^{corr.}\right]}{\left[i_{absenceofinhibitor}^{corr.}\right]}100& (3)\end{array}$

The electrochemical parameters: corrosion potential (E_corr) and corrosion current density (i_corr), were determined from the intersection of anodic and cathodic Tafel slopes and are presented in Table 2.

FIGURE 1

Figure 1. Tafel plot of the mild steel electrode measured in H₂SO₄ solution (c = 0.1 M): in the absence and in the presence of 1 × 10⁻³ M and 1 × 10⁻⁴ M of 2-mercaptonicotinic acid (2MA).

FIGURE 2

Figure 2. The Tafel plot of the mild steel electrode measured in H₂SO₄ solution (c = 0.1 M): in the absence and in the presence of 1 × 10⁻³ M and 1 × 10⁻⁴ M of 4-methyl-4H-1,2,4-triazole-3-thiol (4MT).

TABLE 2

Table 2. The inhibition efficiency and other electrochemical parameters of the 2MA and 4MT inhibitors at 1 × 10⁻³ M and 1 × 10⁻⁴ M toward the mild steel in 0.1 M H₂SO₄corrosion media.

From the Tafel plot in Figures 1, 2, it is evident that the adsorption of the 2MA and 4MT, respectively, molecules onto the mild steel surface strongly lowers the corrosion current of the mild steel in this aggressive media, reflecting a high value of the corrosion inhibition efficiency up to 70.13 ± 1.3% for the 2MA molecule and 63.17 ± 1.4% 4MT. This is in good agreement with the corrosion inhibition efficiency measured through the weight loss measurements (69.85 ± 2.2% for 2MA and 64.52 ± 2.6% for 4MT).

The Resistance of the polarization (Rp) from Tafel extrapolation method was calculated by employing the Stern–Geary Equation (Larabi et al., 2004) Equation (10):

$\begin{array}{ll}{R}_p=\frac{{\beta }_a{\beta }_c}{2.303({\beta }_a{\beta }_c)}\frac{1}{i_{corr.}}& (4)\end{array}$

Table 2 shows that upon increasing the inhibitor concentration, a larger polarization resistance is obtained with both compounds, this reflects the adsorption of the inhibitor on the metal surface which passivates efficiently the active sites and inhibits the corrosion (Mansfeld, 1981; Mohsenifar et al., 2016).

The values of cathodic and anodic Tafel slopes (bc, ba) change when 2MA or 4MT are present in the solution. The Tafel slope differences imply that both of molecules affect the kinetics of the hydrogen evolution reaction (Quartarone et al., 2006). This corresponds to a higher energy barrier for proton discharge, resulting in less gas evolution (Selvi et al., 2003). The studied molecules do not change significantly the corrosion potential, indicating a mixed inhibitor (Berisha et al., 2015b).

From the data presented in Table 3 (calculated from Figure 3) it is obvious that the use of the mixed inhibitors, in fact, does not have an important impact on the IE.

TABLE 3

Table 3. The inhibition efficiency of the mixed inhibitors (2MA/4MT) at different molar ratios from 0 to 1 (each 1 × 10⁻³ M) toward the mild steel in 0.1 M H₂SO₄ corrosion media.

FIGURE 3

Figure 3. The Tafel plot of the mild steel electrode measured in H₂SO₄ (c = 0.1 M) solution: with and without different molar ratio of x(2MA/(4MT) = 1, 0.66, 0.5, 0.33, and x(4MT/(2MA) = 1 (c = 1 × 10⁻³ M for both of the molecules).

Monte Carlo Results

The most stable adsorption configuration of 2MA and 4MT molecules onto Fe (111) and Fe₂O₃ (111) obtained by simulated annealing using the Adsorption Locator module are presented in Figure 4. Both of the studied molecules prefer planar adsorption onto the studied surfaces in each case with the hydrogen atom of the thiol group pointing toward the surface (in the case of 4MT and 2 MA) and the hydrogen atom (2MA) pointing out of the surface.

FIGURE 4

Figure 4. Most likely configuration for the adsorption of 4MT and 2MA on Fe (111) and Fe₂O₃ (111) substrate attained by the adsorption locator module.

The results of the Monte Carlo simulation are presented in Tables 4, 5. In this study, the substrate energy (iron surface) is taken as zero. In addition, adsorption energy in kcal/mol, represents the energy released upon the adsorption of the relaxed adsorbate molecules onto the substrate.

TABLE 4

Table 4. The results calculated by the Monte Carlo simulation of 4MT and 2MA conformations on iron (111) surface.

TABLE 5

Table 5. The results calculated by the Monte Carlo simulation of 4MT and 2MA conformations on iron oxide (111) surface.

The adsorption energy is defined as the sum of the rigid adsorption energy and the deformation energy for the adsorbate components. The rigid adsorption energy reports the energy, in kcal/mol, released (or required) when the unrelaxed adsorbate components (i.e., before the geometry optimization step) are adsorbed on the substrate.

Figure 6 presents the adsorption energy distribution of the 4MT and 2MA molecules on Fe (1 1 1), respectively, Fe₂O₃ (1 1 1). The adsorption energy of 4MT is −68.01 kcal/mol onto Fe (1 1 1) and −64.38 kcal/mol onto Fe₂O₃ (1 1 1) surface, whereas the adsorption energy for 2MA reaches −98.92 kcal/mol in the case Fe (1 1 1) and −109.46 kcal/mol for the Fe₂O₃ (1 1 1) surface. These relatively large values for both molecules reflect their strong adsorption on the studied surfaces. The adsorption energies through Monte Carlo simulations (Table 6) were calculated also with neutral and protonated 4MT and 2MA molecules in the absence and presence of water (the adsorption energy distribution of the adsorbate and the appropriate configuration are found in the Supplementary Figures 2–5).

TABLE 6

Table 6. Adsorption energies calculated by the Monte Carlo simulation of 4MT and 2MA (in neutral and protonated form) on Fe (1 1 1) and Fe₂O₃ (1 1 1) surface, in the presence and absence of water.

The adsorption geometries for 2MA and 4MT are in all cases in the planar configurations. From the calculated results, it is evident that the most important energies, are obtained for the protonated forms of 2MA and 4MT in the presence of water. In this case for the Fe (111) surface the adsorption energy of 2MA-H⁺ has a relatively large value −254.99 kcal/mol, to be compared with that 4MT-H⁺ −159.15 kcal/mol. These large values indicate a strong interaction of these molecules with the studied surfaces.

Quantum Chemical Calculations

The shape and symmetry of the HOMO and the LUMO are the central factors for the estimation of the reactivity of a compound (Gece, 2008; Obot et al., 2015). The analysis of the HOMO highlights the areas of the molecule that can donate electrons to electrophilic species while the analysis of the LUMO predicts the regions of the molecule with high affinity to accept electrons from nucleophilic species. In 2MA (Figure 5), the HOMO orbital has the highest electron density on the N1-C6, C4-C5 single bonds, and C2 = C6 double bond, this points out that these are the areas of the molecule with the highest affinity to donate electrons; the LUMO is localized on C3, C5, and C6 atoms. In the 4MT molecule, the HOMO orbital is extended over a large region between the C2-N3, N1-N2 single bonds, C2 = N1 double bond and S atom; the LUMO is strongly localized at the −SH group. Ionization potential is directly related to the energy of HOMO orbital and the electron affinity to that of LUMO orbital. Our results show that the HOMO and the LUMO are delocalized throughout molecule. In the investigated compounds (Table 7), the E_HOMO is larger in 4MT, consequently, 4MT would have the highest affinity to adsorb onto the metal surface and to offer electrons (the lone pair electrons on S and N atoms) to the unoccupied iron d-orbital (Fouda and Ellithy, 2009).

FIGURE 5

Figure 5. Optimized structures: A1; HOMO and LUMO; A2, Electrostatic potential map and A3, Mulliken charge for 2MA and 4MT molecules. (B3LYP/6-31G* results in vacuo).

TABLE 7

Table 7. The molecular properties of the investigated compounds using: DFT, B3LYP/6-31G (d,p). Results in vacuo.

These results are in good agreement with experimentally observed inhibition efficiencies. A negative value for E_HOMO suggests a physical adsorption of the organic molecule (Yurt et al., 2006). The energy of the LUMO shows the tendency of the molecule to accept electrons from an electron-rich species. The molecule with the lowest energy of the LUMO has the highest tendency to accept electrons. This is in perfect agreement with the calculated values E_LUMO 4MT < E_LUMO 2MA implying that 4MT would preferentially accept electrons from the metal surface more than 2MA.

The energy difference between the E_HOMO and the E_LUMO, ΔE, expresses the reactivity of the molecule by comparison with other compounds. The molecules with the lowest ΔE value express the highest affinity reactivity and would favorably interact with the surface, this is not in agreement with what we found ΔE (4MT) > ΔE (2MA).

The global hardness and softness parameters are commonly described in terms of hard-soft-acid-base (HSAB) theory; by this theory, soft acids act preferentially with soft bases and hard ones with hard bases. Metals are normally considered to be soft acids (Obi-Egbedi et al., 2011b); therefore, they would favorably interact with the inhibitors that display low η and high σ values.

The σ value for the investigated compounds is higher for 4MT > 2MA which fully agrees with the experimentally determined inhibition efficiency of the inhibitors. The dipole moment tells about the polarity of the molecule. It has been established that an increased inhibition efficiency is related to an increased dipole moment (Obi-Egbedi et al., 2011b; Berisha et al., 2015b). The dipole moment for the studied compounds is the higher for 4MT, which fits with our experimental results. Certain molecular properties do not only specify the reactivity of molecules but they also designate the position selectivity in the molecules, i.e., the areas on which particular types of reactions are more likely to take place. One of the parameters in this group is the partial atomic charge on the atoms. Frequently the interaction between the metal and the inhibitor is occurs favorably on the atom bearing the utmost negative charge. The Mulliken atomic charges for 4MT and 2MA are shown in Figure 5 and the highest values of the negative charge on the molecules are to be found on the heteroatoms, suggesting that these centers have the maximum electron density and would preferentially interact with the metal surface. In order to obtain more realistic results, DFT calculations were also performed using water as a solvent with 2MA and 4MT in neutral and protonated state. The results are very similar and lead to the same conclusions as above (Supplementary Table 1) indicating that solvation of the surface and molecules is not determining in the corrosion process.

DFT Calculations with a Plane Wave Basis Set

As Monte Carlo calculations are based on the use of the force fields and as there is no electron correlation (Kalos and Whitlock, 2008; Akkermans et al., 2013; Schattke and Muiño, 2013), more accurate calculations are used to estimate a correct adsorption energy and geometries. This is performed through DFT periodic calculations with a plane wave basis set.

The planar geometry is preferred for the adsorption of the 2MA molecule onto the Fe surface (Figure 6). The S atom of the thiol group of the 2MA molecule is positioned at 2.33 Å from the Fe nanocluster, whereas the N of the pyridine ring is at 1.95 Å. The hydrogen atom of the thiol group is oriented away from the surface, indicating that the adsorption is taking place through the S atom, pyridine ring and the carboxyl group (positioned almost flat onto the cluster surface).

FIGURE 6

Figure 6. The optimized geometries with adsorption distances of 2MA molecule (A1/B1 neutral and A2/B2 protonated form) onto: (A) Fe and (B) Fe₂O₃ cluster.

The computed (by periodic DFT calculations) adsorption energy of −122.57 kcal/mole (Figure 7) supports this strong interaction of this molecule with the iron cluster surface, this also evidenced by analyzing the DOS plot (Figure 8A).

FIGURE 7

Figure 7. Adsorption energies for 2MA and 4MT (in neutral and protonated forms) onto Fe and Fe₂O₃ clusters.

FIGURE 8

Figure 8. The DOS plot for: (A) 2MA, Fe and 2MA/Fe and (B). (A) 2MA-H+, Fe and 2MA-H+/Fe cluster system.

The Fermi level of the Fe cluster after the interaction with the 2MA molecule shifts upwards, with a peak broadening due to the overlap of the valence band orbitals of the Fe cluster with those of the 2MA molecule (Cho et al., 2008; Raouafi et al., 2016).

In the protonated form, the distances of the 2MA molecules onto the cluster are slightly increased, the S atom is positioned at 2.34 Å, the N of the pyridine ring 2.0 Å and the adsorption energy in this case is increased up to −188.30 kcal/mol. The pyridine ring interacts with iron surface through the nitrogen atom via coordinative Fe-N bonding (Gu et al., 2001). This interaction is evidenced in the DOS plot (Figure 8B) by a more pronounced peak broadening compared with the neutral form.

The adsorption of 2MA onto the Fe₂O₃ surface is planar as with the 4MA. The S atom distance from the surface of this cluster is at 3.71 Å, whereas the N of the pyridine ring is at 2.14 Å. These distances are larger than in the case of Fe (1.375 Å for S atom and 0.162 Å for pyridine ring nitrogen). The hydrogen atom of the thiol group is parallel to the surface, as well as the pyridine ring and the carboxyl group. The relatively insignificant value of adsorption energy −6.1 kcal/mole (2MA) and −25.6 kcal/mol (2MA-H⁺) points out a weaker interaction at this surface as proved also via evaluating the DOS plot (Guo et al., 2017b) and the increase of the adsorption distance of the molecule compared to the Fe. The upward shift of the Fermi level after the interaction with the 2MA or 2MA-H⁺ is less apparent (Figures 9A,B). Minor peak broadening occurs as a result of the slight overlap of the valence band orbitals of the Fe₂O₃ cluster with the valence orbitals of 2MA molecule (Jarvis et al., 2015).

FIGURE 9

Figure 9. The DOS plot for: (A) 2MA, Fe₂O₃ and 2MA/Fe₂O₃ and (B) (A) 2MA-H+, Fe₂O₃ and 2MA-H+/ Fe₂O₃ cluster system.

The 4MT molecule adsorbs flat onto the Fe surface (Figure 10). The S atom distance of the thiol group from the Fe surface is 2.43 Å, the N of the triazole ring next to the thiol group is at 1.99 Å. Contrary to the 2MA adsorption, the hydrogen atom of the thiol group is facing the surface, indicating that the adsorption is mainly due to the nitrogen atoms of the azole ring. The calculated adsorption energy is lower by about −19.18 kcal/mole compared to the 2MA. This surface/adsorbate interaction is supported by analyzing the upward shift of the Fermi level (with peak broadening) of the 4MT or 4MT-H+ /Fe surface (Figures 11, 12) (Jarvis et al., 2015; Raouafi et al., 2016).

FIGURE 10

Figure 10. The optimized geometries with adsorption distances of 4MT molecule (A1/B1 neutral and A2/B2 protonated form) onto: (A) Fe and (B) Fe₂O₃ cluster.

FIGURE 11

Figure 11. The DOS plot for: (A) 4MT, Fe and 4MT/Fe and (B) (A) 4MT -H+, Fe and 4MT-H+/Fe cluster system.

FIGURE 12

Figure 12. The DOS plot for: (A) 4MT, Fe₂O₃ and 4MT/ Fe₂O₃ and (B) (A) 4MT -H+, Fe₂O₃ and 4MT-H+/ Fe₂O₃ cluster system.

In the protonated form, the distances of the 4MT molecule onto the cluster are slightly decreased, the S atom is positioned at 2.25 Å, the N of the ring is at 1.98 Å. The hydrogen atom of the thiol group is facing out the surface, exposing the S atom to the surface.

The flat geometry is the favored one for both 4MT and 4MT-H⁺ (Figure 10) onto the Fe₂O₃ cluster. The S atom of the thiol group of the 4MT molecule is at 3.47 Å (for 4MT), respectively, for 3.67 Å from the clusters surface. The N atom (next to the thiol group) of the triazole ring is at 3.30 Å for 4MT and 3.14 Å for 4MT-H⁺. The adsorption energies are greater for the protonated than the neutral form of 4MT.

The Inhibition Mechanism

In general the corrosion inhibition mechanism in acid medium is the adsorption of inhibitor molecule onto the metals surface by one of the four types of adsorption: (1) electrostatic attraction between charged molecules and the charged metal, (2) interaction of unshared electron pairs in the molecule with the metal, (3) interaction of π-electrons with the metal, and (4) a combination of the above (Bhajiwala and Vashi, 2001; Singh et al., 2011). For the physical adsorption to take place in acidic medium, the presence of both a metal surface (with vacant low-energy electron orbitals) and charged species in the solution (a molecule having reasonably loosely bound electrons or heteroatom with lone pair electrons) are necessary.

Coordinate covalent bond formation between electron pairs of unprotonated S atom and N atoms of aromatic rings with metal surface can take place for both of the molecules during their adsorption (Gu et al., 2001). This leads to a protective film on mild steel surface that reduces the corrosion rate in good agreement with the experimental results. Through the above calculation we show that the main interactions take place through the heteroatoms. They show that that these molecules lie flat on the surface (Guo et al., 2017a). Finally an organic film forms on the surface with a interaction energy in the range of −122.57 kcal/mol for 2MA-H⁺, respectively, −103.39 for 4MT-H⁺ that prevents the corroding species from reaching the metallic surface (Singh and Quraishi, 2010; Karthik and Sundaravadivelu, 2013).

Conclusion

The inhibition efficiency of 4MT and 2MA has been explored experimentally and by quantum calculations. Potentiodynamic measurements reveals a strong inhibition behavior of these molecules toward mild steel corrosion. Different calculations were performed to rationalize these results: theoretical parameters [E_HOMO, E_LUMO, ΔE (HOMO-LUMO), global hardness (η), global softness (σ), and the dipole moment (μ)], Monte Carlo simulations and DFT quantum chemical calculations using periodic systems.

The calculations permit to shed some light on the fundamental reasons for this reduced corrosion rates. The active centers which would interact with the metal surface are determined as well as the adsorption geometries and energies. Experimental and theoretical results show that the inhibition efficiency is higher for 2MA due to the presence of the pyridine nitrogen, carboxyl, and the sulfur atom in the thiol group interacting in a planar geometry with the surface. These results show the interest of theoretical calculation in the interpretation of corrosion measurements.

Author Contributions

AB proposed the idea and run the theoretical calculations. VM realized the experimental results. The article was written jointly by both of the authors.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are grateful to the OeAD for their financial support. The Frontiers Waivers support is kindly appreciated by the authors. The authors gratefully acknowledge also the support from Ministry of Education, Science and Technology of Kosovo (Nr.2-5069) for providing them with the computing resources.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fchem.2017.00061/full#supplementary-material

References

Abohalkuma, T., and Telegdi, J. (2015). Corrosion protection of carbon steel by special phosphonic acid nano-layers. Mater. Corros. 66, 1382–1390. doi: 10.1002/maco.201508304

CrossRef Full Text | Google Scholar

Akkermans, R. L., Spenley, N. A., and Robertson, S. H. (2013). Monte Carlo methods in materials studio. Mol. Simul. 39, 1153–1164. doi: 10.1080/08927022.2013.843775

CrossRef Full Text | Google Scholar

Al-Amiery, A. A., Kadhum, A. A. H., Kadihum, A., Mohamad, A. B., How, C. K., and Junaedi, S. (2014). Inhibition of mild steel corrosion in sulfuric acid solution by new Schiff base. Materials 7, 787–804. doi: 10.3390/ma7020787

PubMed Abstract | CrossRef Full Text | Google Scholar

Bentiss, F., Lagrenee, M., Traisnel, M., and Hornez, J. (1999). The corrosion inhibition of mild steel in acidic media by a new triazole derivative. Corros. Sci. 41, 789–803. doi: 10.1016/S0010-938X(98)00153-X

CrossRef Full Text | Google Scholar

Berisha, A., Combellas, C., Kanoufi, F., Decorse, P., Oturan, N., Médard, J., et al. (2017). Some theoretical and experimental insights on the mechanistic routes leading to the spontaneous grafting of gold surfaces by diazonium salts. Langmuir. doi: 10.1021/acs.langmuir.7b01371

PubMed Abstract | CrossRef Full Text | Google Scholar

Berisha, A., Combellas, C., Kanoufi, F., Pinson, J., and Podvorica, F. I. (2011). Physisorption vs grafting of aryldiazonium salts onto iron: a corrosion study. Electrochim. Acta 56, 10762–10766. doi: 10.1016/j.electacta.2011.01.049

CrossRef Full Text | Google Scholar

Berisha, A., Chehimi, M. M., Pinson, J., and Podvorica, F. I. (2015a). “Electrode surface modification using diazonium salts”, in Electroanalytical Chemistry, A Series of Advances, Vol. 26, eds A. J. Bard and C. G. Zoski (Boca Raton, FL: CRC Press), 115–224. doi: 10.1201/b19196-4

CrossRef Full Text

Berisha, A., Podvorica, F., Mehmeti, V., Syla, F., and Vataj, D. (2015b). Theoretical and experimental studies of the corrosion behavior of some thiazole derivatives toward mild steel in sulfuric acid media. Macedonian J. Chem. Chem. Eng. 34, 287–294. doi: 10.20450/mjcce.2015.576

CrossRef Full Text | Google Scholar

Bhajiwala, H., and Vashi, R. (2001). Ethanolamine, diethanolamine and triethanolamine as corrosion inhibitors for zinc in binary acid mixture [HNO3+ H3PO4]. Bull. Electrochem. 17, 441–448.

Google Scholar

Bowker, M., Hutchings, G., Davies, P. R., Edwards, D., Davies, R., Shaikhutdinov, S., et al. (2012). Surface structure of γ-Fe2O3(111). Surf. Sci. 606, 1594–1599. doi: 10.1016/j.susc.2012.06.010

CrossRef Full Text | Google Scholar

Chaussé, A., Chehimi, M. M., Karsi, N., Pinson, J., Podvorica, F., and Vautrin-Ul, C. (2002). The electrochemical reduction of diazonium salts on iron electrodes. The formation of covalently bonded organic layers and their effect on corrosion. Chem. Mater. 14, 392–400. doi: 10.1021/cm011212d

CrossRef Full Text | Google Scholar

Cho, Y., Kim, C., Moon, H., Choi, Y., Park, S., Lee, C.-K., et al. (2008). Electronic structure tailoring and selective adsorption mechanism of metal-coated nanotubes. Nano Lett. 8, 81–86. doi: 10.1021/nl0720051

PubMed Abstract | CrossRef Full Text | Google Scholar

Clark, S. J., Segall, M. D., Pickard, C. J., Hasnip, P. J., Probert, M. I., Refson, K., et al. (2005). First principles methods using CASTEP. Z. für Kristallographie Crystalline Mater. 220, 567–570. doi: 10.1524/zkri.220.5.567.65075

CrossRef Full Text | Google Scholar

Döner, A., Solmaz, R., Özcan, M., and Kardaş, G. (2011). Experimental and theoretical studies of thiazoles as corrosion inhibitors for mild steel in sulphuric acid solution. Corros. Sci. 53, 2902–2913. doi: 10.1016/j.corsci.2011.05.027

CrossRef Full Text | Google Scholar

Finšgar, M., and Jackson, J. (2014). Application of corrosion inhibitors for steels in acidic media for the oil and gas industry: a review. Corros. Sci. 86, 17–41. doi: 10.1016/j.corsci.2014.04.044

CrossRef Full Text | Google Scholar

Fouda, A., and Ellithy, A. (2009). Inhibition effect of 4-phenylthiazole derivatives on corrosion of 304L stainless steel in HCl solution. Corros. Sci. 51, 868–875. doi: 10.1016/j.corsci.2009.01.011

CrossRef Full Text | Google Scholar

Gao, W., Cui, T. T., Zhu, Y. F., Wen, Z., Zhao, M., Li, J. C., et al. (2015). Design principles of inert substrates for exploiting gold clusters' intrinsic catalytic reactivity. Sci. Rep. 5:15095. doi: 10.1038/srep15095

PubMed Abstract | CrossRef Full Text | Google Scholar

Gece, G. (2008). The use of quantum chemical methods in corrosion inhibitor studies. Corros. Sci. 50, 2981–2992. doi: 10.1016/j.corsci.2008.08.043

CrossRef Full Text | Google Scholar

Geerlings, P., De Proft, F., and Langenaeker, W. (2003). Conceptual density functional theory. Chem. Rev. 103, 1793–1874. doi: 10.1021/cr990029p

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, R., Cao, P., Yao, J., Ren, B., Xie, Y., Mao, B., et al. (2001). Surface Raman spectroscopic studies on the adsorption of pyridine at bare iron electrodes. J. Electroanal. Chem. 505, 95–99. doi: 10.1016/S0022-0728(01)00478-8

CrossRef Full Text | Google Scholar

Guo, L., Obot, I. B., Zheng, X., Shen, X., Qiang, Y., Kaya, S., et al. (2017a). Theoretical insight into an empirical rule about organic corrosion inhibitors containing nitrogen, oxygen, and sulfur atoms. Appl. Surf. Sci. 406, 301–306. doi: 10.1016/j.apsusc.2017.02.134

CrossRef Full Text | Google Scholar

Guo, L., Qi, C., Zheng, X., Zhang, R., Shen, X., and Kaya, S. (2017b). Toward understanding the adsorption mechanism of large size organic corrosion inhibitors on an Fe (110) surface using the DFTB method. RSC Adv. 7, 29042–29050. doi: 10.1039/C7RA04120A

CrossRef Full Text | Google Scholar

Hammer, B., Hansen, L. B., and Nørskov, J. K. (1999). Improved adsorption energetics within density-functional theory using revised Perdew-Burke-Ernzerhof functionals. Phys. Rev. B 59:7413. doi: 10.1103/PhysRevB.59.7413

CrossRef Full Text | Google Scholar

Jarvis, S. P., Taylor, S., Baran, J. D., Thompson, D., Saywell, A., Mangham, B., et al. (2015). Physisorption controls the conformation and density of states of an adsorbed porphyrin. J. Phys. Chem. C 119, 27982–27994. doi: 10.1021/acs.jpcc.5b08350

CrossRef Full Text | Google Scholar

Kalos, M. H., and Whitlock, P. A. (2008). Monte Carlo Methods. John Wiley & Sons.

PubMed Abstract | Google Scholar

Karthik, G., and Sundaravadivelu, M. (2013). Inhibition of mild steel corrosion in sulphuric acid using esomeprazole and the effect of iodide ion addition. ISRN Electrochemistry 2013, 10. doi: 10.1155/2013/403542

CrossRef Full Text | Google Scholar

Khaled, K., Abdelshafi, N., El-Maghraby, A., Aouniti, A., Al-Mobarak, N., and Hammouti, B. (2012). Alanine as corrosion inhibitor for iron in acid medium: a molecular level study. Int. J. Electrochem. Sci. 7, 12706–12719.

Google Scholar

Kosian, M., Smulders, M. M., and Zuilhof, H. (2016). Structure and long-term stability of alkylphosphonic acid monolayers on SS316L stainless steel. Langmuir 32, 1047–1057. doi: 10.1021/acs.langmuir.5b04217

PubMed Abstract | CrossRef Full Text | Google Scholar

Larabi, L., Harek, Y., Traisnel, M., and Mansri, A. (2004). Synergistic influence of poly (4-vinylpyridine) and potassium iodide on inhibition of corrosion of mild steel in 1M HCl. J. Appl. Electrochem. 34, 833–839. doi: 10.1023/B:JACH.0000035609.09564.e6

CrossRef Full Text | Google Scholar

Mansfeld, F. (1981). Recording and analysis of AC impedance data for corrosion studies. Corrosion 37, 301–307. doi: 10.5006/1.3621688

CrossRef Full Text | Google Scholar

Migahed, M., Azzam, E., and Al-Sabagh, A. (2004). Corrosion inhibition of mild steel in 1 M sulfuric acid solution using anionic surfactant. Mater. Chem. Phys. 85, 273–279. doi: 10.1016/j.matchemphys.2003.12.027

CrossRef Full Text | Google Scholar

Mohsenifar, F., Jafari, H., and Sayin, K. (2016). Investigation of thermodynamic parameters for steel corrosion in acidic solution in the presence of N, N′-Bis (phloroacetophenone)-1, 2 propanediamine. J. Bio Tribo Corros. 2, 1–13. doi: 10.1007/s40735-015-0031-y

CrossRef Full Text | Google Scholar

Obi-Egbedi, N., Essien, K., Obot, I., and Ebenso, E. (2011a). 1, 2-Diaminoanthraquinone as corrosion inhibitor for mild steel in hydrochloric acid: weight loss and quantum chemical study. Int. J. Electrochem. Sci. 6, 913–930.

Google Scholar

Obi-Egbedi, N., Obot, I., and El-Khaiary, M. I. (2011b). Quantum chemical investigation and statistical analysis of the relationship between corrosion inhibition efficiency and molecular structure of xanthene and its derivatives on mild steel in sulphuric acid. J. Mol. Struct. 1002, 86–96. doi: 10.1016/j.molstruc.2011.07.003

CrossRef Full Text | Google Scholar

Obot, I., Macdonald, D., and Gasem, Z. (2015). Density functional theory (DFT) as a powerful tool for designing new organic corrosion inhibitors. Part 1: an overview. Corros. Sci. 99, 1–30. doi: 10.1016/j.corsci.2015.01.037

CrossRef Full Text | Google Scholar

Palanivel, V., Huang, Y., and Van Ooij, W. J. (2005). Effects of addition of corrosion inhibitors to silane films on the performance of AA2024-T3 in a 0.5 M NaCl solution. Prog. Org. Coatings 53, 153–168. doi: 10.1016/j.porgcoat.2003.07.008

CrossRef Full Text | Google Scholar

Parr, R. G., Szentpaly, L. V., and Liu, S. (1999). Electrophilicity index. J. Am. Chem. Soc. 121, 1922–1924. doi: 10.1021/ja983494x

CrossRef Full Text | Google Scholar

Payne, M. C., Teter, M. P., Allan, D. C., Arias, T., and Joannopoulos, J. (1992). Iterative minimization techniques for ab initio total-energy calculations: molecular dynamics and conjugate gradients. Rev. Mod. Phys. 64:1045. doi: 10.1103/RevModPhys.64.1045

CrossRef Full Text | Google Scholar

Quartarone, G., Bonaldo, L., and Tortato, C. (2006). Inhibitive action of indole-5-carboxylic acid towards corrosion of mild steel in deaerated 0.5 M sulfuric acid solutions. Appl. Surf. Sci. 252, 8251–8257. doi: 10.1016/j.apsusc.2005.10.051

CrossRef Full Text | Google Scholar

Raouafi, F., Seydou, M., Lassoued, K., Taleb, A., and Diawara, B. (2016). First-principles investigation of methanethiol adsorption and dissociation mechanisms on the high-Miller-index vicinal surface Cu (4 1 0). J. Phys. Condens. Matter 28:175001. doi: 10.1088/0953-8984/28/17/175001

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanyal, B. (1981). Organic compounds as corrosion inhibitors in different environments—a review. Prog. Org. Coatings 9, 165–236. doi: 10.1016/0033-0655(81)80009-X

CrossRef Full Text | Google Scholar

Schattke, W., and Muiño, R. D. (2013). Quantum Monte-Carlo Programming: for Atoms, Molecules, Clusters, and Solids. John Wiley & Sons.

Google Scholar

Selvi, S. T., Raman, V., and Rajendran, N. (2003). Corrosion inhibition of mild steel by benzotriazole derivatives in acidic medium. J. Appl. Electrochem. 33, 1175–1182. doi: 10.1023/B:JACH.0000003852.38068.3f

CrossRef Full Text | Google Scholar

Shubha, H., Venkatesha, T., Vathsala, K., Pavitra, M., and Punith Kumar, M. (2013). Preparation of self assembled sodium oleate monolayer on mild steel and its corrosion inhibition behavior in saline water. ACS Appl. Mater. Interfaces 5, 10738–10744. doi: 10.1021/am4028857

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, A. K., and Quraishi, M. (2010). Effect of cefazolin on the corrosion of mild steel in HCl solution. Corros. Sci. 52, 152–160. doi: 10.1016/j.corsci.2009.08.050

CrossRef Full Text | Google Scholar

Singh, A. K., Shukla, S. K., and Quraishi, M. (2011). Corrosion behaviour of mild steel in sulphuric acid solution in presence of ceftazidime. Int. J. Electrochem. Sci. 6, 5802–55814.

Google Scholar

Su, X., Lai, C., Peng, L., Zhu, H., Zhou, L., Zhang, L., et al. (2016). A dialkyldithiophosphate derivative as mild steel corrosion inhibitor in sulfuric acid solution. Int. J. Electrochem. Sci. 11, 6412–6412. doi: 10.20964/2016.07.101

CrossRef Full Text | Google Scholar

Van Ooij, W., Zhu, D., Stacy, M., Seth, A., Mugada, T., Gandhi, J., et al. (2005). Corrosion protection properties of organofunctional silanes—an overview. Tsinghua Sci. Technol. 10, 639–664. doi: 10.1016/S1007-0214(05)70134-6

CrossRef Full Text | Google Scholar

Yesudass, S., Olasunkanmi, L. O., Bahadur, I., Kabanda, M. M., Obot, I., and Ebenso, E. E. (2016). Experimental and theoretical studies on some selected ionic liquids with different cations/anions as corrosion inhibitors for mild steel in acidic medium. J. Taiwan Inst. Chem. Eng. 64, 252–268. doi: 10.1016/j.jtice.2016.04.006

CrossRef Full Text | Google Scholar

Yurt, A., Ulutas, S., and Dal, H. (2006). Electrochemical and theoretical investigation on the corrosion of aluminium in acidic solution containing some Schiff bases. Appl. Surf. Sci. 253, 919–925. doi: 10.1016/j.apsusc.2006.01.026

CrossRef Full Text | Google Scholar