SciELO - Scientific Electronic Library Online

 
vol.40 número2Inhibitive Action of Aniline on Zinc Corrosion in a H2SO4 Solution: Electrochemical Study índice de autoresíndice de assuntosPesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

  • Não possue artigos similaresSimilares em SciELO

Compartilhar


Portugaliae Electrochimica Acta

versão impressa ISSN 0872-1904

Port. Electrochim. Acta vol.40 no.2 Coimbra abr. 2022  Epub 30-Abr-2022

https://doi.org/10.4152/pea.2022400206 

Research Article

Multi-component Theoretical Assessment of the Corrosion Inhibitory Performance of some Macrocyclic Polyethers Containing 1,3,4 Thiadiazole on the Iron Surface

A. Mahsoune1 

K. Sadik1 

S. Byadi1 

A. Aboulmouhajir1  2 

1Team of Molecular Modelling and Spectroscopy, Sciences Faculty, Chouaib Doukkali University, Morocco

2Organic Synthesis, Extraction and Valorization Laboratory, Team of Extraction, Spectroscopy and Valorization, Sciences Faculty of Ain Chock, Hassan II University, Morocco


Abstract

The adsorption capability of the series of some Macrocyclic Polyether compounds containing 1, 3, 4-thiadiazole entity n-MCTH (n = 1-5), especially 3- MCTH, 4-MCTH and 5 MCTH, and their protonated forms, was studied in the light of DFT quantum modeling and Monte Carlo dynamics calculations. Sensitivity to corrosion has been quantified using the degree of planarity, global and local electronic proprieties, as well as the inhibitor strength of interaction, in neutral and protonated forms, with the (111) iron surface in the metallic complex. The results of both approaches showed the supremacy of the interactions of neutral and proton variants of the 5-MCTH-Fe complexes, compared to their homologues of 3-MCTH-Fe and 4-MCTH-Fe, due to the significant involvement of aryl rings, in addition to the thiadiazole ring, in the process of electron donation and acceptance.

Keywords: DFT calculations; Monte Carlo dynamics calculations; global electronic proprieties; local electronic proprieties; complexation strength

Introduction

The use of inhibitors is one of the best options for metals protection against corrosion (1, 2). However, the process mechanistic is very complex, hence, the usefulness of computational tools to elucidate the interaction between the adsorbed inhibitor molecules and the corroding metal surface.

Most corrosion inhibitors, which slow down the corrosion rate of a metal through the mechanism of adsorption onto the metal surface , (3) are organic compounds containing hetero atoms (such as N, O, S and P), as well as the π orbital character of donating electrons into their carbon chain systems.

Heterocyclic systems rich in electrons, such as phenyls or thiadiazoles, have constantly attracted the attention of researchers in the corrosion field. (4, 5) The series of 2,5-bis (n-pyridyl) -1,3,4-thiadiazoles (n-PTH: n = 2-4), 2,5-bis (phenyl) -1,3,4-thiadiazole (DPTH), 2,5-bis (4-methoxyphenyl) -1,3 molecules, 4-thiadiazole (4-MTH), 2,5-bis (4-dimethylaminophenyl) -1,3,4-thiadiazole (4-DATH), 2,5-bis (4-methylphenyl) -1,3,4- thiadiazole (4-MPTH), 2,5-bis (4-nitrophenyl) -1,3,4-thiadiazole (4-NPTH) and 2,5-bis (4-chlorophenyl) -1,3,4-thiadiazole (4-CPTH), are known for their corrosion inhibiting properties (6). The Electrochemical Impedance Spectroscopy (EIS) measurements of n-PTH (n = 2 - 4) corrosion inhibition on mild steel, in different acidic media (1 M HCl, 0.5 M H2SO4 and 1 M HClO4), have detected a noticeable inhibition, which was assigned, through primary quantum calculations, to the electron density of heterocyclic rings N and thiadiazole ring S atoms (7) In the same way, the inhibition performance, on mild steel in 1 M HCl solutions, of DPTH, 4- MTH, 4- DATH, 4-MPTH, 4-NPTH and 4-CPTH, has been theoretically investigated by using a limited number of quantum molecular descriptors, such as dipole moment (μ), highest occupied (EHOMO) and lowest unoccupied (ELUMO) molecular orbitals and the differences between them (HOMO-LUMO gap) . (8) Similarly, the macrocyclic polyethers, such the series of macrocyclic polyether compounds (n-MCTH, n: 1-5), which are shown in Fig. 1, synthesized by Bentiss et al. (9), also exhibit an electronic richness, since they contain one thiadiazole ring, polyether macrocyclic oxygen and two phenyls rings. Their corrosion inhibitory effect on carbon steel in 1 M HCl has been reported earlier by the Weight Loss Method (WLM), EIS and potentiodynamic polarization (10, 11, 12). It was shown that the macrocyclic polyethers inhibition efficiency increased with their concentration in the order of 5-MCTH > 4-MCTH > 3-MCTH > 2-MCTH > 1-MCTH. The same authors tried to predict the experimental inhibition efficiency classification via B3LYP/6-31G (d, p) DFT level. Based only on a limited number of global descriptors (EHOMO, ELUMO, ΔE, μ, molecular area and CCCS dihedral angle), the prediction remains very incomplete towards a deep understanding of the anti-corrosive adsorption. That is why we have largely explored the ab initio and molecular dynamics interaction of the iron surface with n-MCTH, to elucidate the impact of the thiadiazole nucleus enriched by the aromatic rings electrons, and especially that of the oxygen atoms of the macrocycle polyether part.

Figure 1 Series of macrocyclic polyether compounds structures (n-MCTH, n: 1-5). 

Recently (13), we have tried to elucidate the mechanism of two macrocyclic polyether compounds, 1-MCTH and 2-MCTH, and their corrosion inhibition protonated forms, as well as the intrinsic reasons for their inhibition efficiency classification, by using DFT and Dynamic Monte Carlo calculations. Therefore, as a logical continuation of this work, we present ab initio and molecular dynamics investigations of iron surface interactions with the remaining molecules of the n-MCTH (n = 1-5) series, which are 3- MCTH, 4-MCTH and 5 MCTH.

Calculation details

Reactivity quantum global and local descriptors

All geometry optimizations and quantum chemical calculations were performed using B3LYP/6-31G (d,p) level (14), by means of the Gaussian 2009 (15) and Gaussview 05 package. To model the solvation effect, the PCM (Polarizable Continuum Model), considering the solvent as a continuum and the solute in its cavities, was adopted . (16)

As well described in our earlier work (13), quantum electronic parameters investigation concerns:

  • The molecular descriptors such as the energy of highest occupied molecular orbital (EHOMO), the energy of the lowest unoccupied molecular orbital (ELUMO), the energy gap (∆E) (17), the global hardness (η) (18), the absolute electronegativity (χ), the overall softness (S), the electrophile index (ω), the fraction of electrons transferred (ΔN) (19) and the total exchange energy (ΔE) back-donation.

  • The local descriptors, such as molecular electrostatic potential (MEP), NBO charges, NBO energy of hyper-conjugative interactions (17) and Fukui indexes (fk-, fk+ and fk0) of the nucleophilic, electrophilic and radical attack, respectively (20) (19), knowing that fk- higher value gives rise to a nucleophilic site, and that fk+ higher value gives rise to an electrophilic site. The formulations of these electronic characteristics were mentioned elsewhere (21, 22).

Dynamic Monte Carlo simulation

Likewise for n-MCTH (n:1-2) (13), the Metropolis Dynamic Monte Carlo simulation, using the adsorption locator and Forcite codes implemented in the Material Studio 7.0 software (23), was carried out to locate the stronger adsorption between the inhibitor molecules and the iron surface, i.e. the higher negative interaction energy value. Although Fe (110) plane is a densely-packed surface, the crystallographic Fe (111) surface was chosen because of the lower value of its formation energy (24). The structures of all components of the inhibitor-iron surface system were optimized with COMPASS II force field, and the simulation box (17.89×17.89×38.34Å) was used with periodic boundary conditions, in order to simulate a representative part of an interface. The Fe (111) plane was next enlarged to a (9 × 9) super cell. After that, a vacuum slab with 3.0 nm thickness was built above the Fe (111) plane. To mimic the real experimental corrosion environment, 100 molecules of water were added to the simulation box.

Results and discussion

In addition to the solvent effect, the structural properties, mainly, the level of planarity of n-MCTH (n = 3-5) molecules, and their competitive protonated forms, on the one hand, and their electronic properties, on the other hand, have an impact on their inhibition efficiency.

Regio-selectivity of neutral species isolated and in an aqueous solution

The optimized geometries of 3-MCTH, 4-MCTH and 5-MCTH are illustrated in Fig. 2.

Figure 2 Schematic representation of 2D structures and 3-MCTH, 4-MCTH and 5-MCTH optimized structures, at B3LYP / 6-31G **. 

Selected geometrical parameters

As confirmed by a structural comparison of 1-MCTH and 2-MCTH, the increase in the number of oxygen atoms in 3-MCTH, 4-MCTH and 5-MCTH makes the ring containing oxygen less constrained. Indeed, the distance between the two aryls increases (C24-C25 distance increases up to 0.15 A° from 3-MCTH to 4-MCTH and 5-MCTH), and the dihedral angle between the two aryls described by C20-C3-C11-C22 decreases, until making them quasi-coplanar for 5-MCTH (Tables 1 and 2). It is the same with the thiadiazole ring that also becomes coplanar with aryls, as C24-C3-C6-S7 and C25-C11-C8-S7 have reached almost zero. However, the aqueous solvent acts against the planarity of these molecules, as C20-C3-C11-C22 and C25-C11-C8-S7 increase up to 20° for 5-MCTH, e.g.

Table 1. Bond distance (Å) for the neutral inhibitors 3-MCTH, 4-MCTH and 5-MCTH optimized structures and their double protonated forms (N9H+N10H+), in vacuum (G) and in an aqueous solution (A), at B3LYP/6-31G**. 

Table 2 Torsional angle (º) for the 3-MCTH, 4-MCTH and 5-MCTH neutral inhibitors optimized structures and their double protonated forms (N9H+N10H+), in vacuum (G) and in an aqueous solution (A), at B3LYP/6-31G**. 

Taking into account the structural aspect, 5-MCTH can be considered the most appropriate corrosion inhibitor, because it has a very planar first molecular part, and a second part with a large recovering surface. In comparison with the available experimental data, 4-MCTH X-ray crystallographic data show that its structure belongs to the C 2/c space group of the monoclinic system, and that it is practically plane. Moreover, the phenyl rings are shifted relative to the thiadiazole ring by an angle equal to 4.93°. In addition, a weak hydrogen bond connects the polyether macrocyclic oxygen and phenyls, providing two-dimensional stability of the crystalline system (10).

Molecular quantum chemical parameters

The results presented in Fig. 3 show that HOMO is delocalized throughout all aryl and thiadiazole rings, except the S atom of thiadiazole, and C16, C18, C20 and C22 of aryl rings, which are excluded from the electron donation regions. LUMO is spread out on all the aryl and thiadiazole rings, providing them high tendency to accept electrons from nucleophilic species. However, the formation of the LUMO orbital does not involve the Pz orbitals participation. Moreover, oxygen atoms are less concerned by the electrons transfer effect.

Figure 3 HOMO and LUMO 3-MCTH, 4-MCTH and 5-MCTH neutral inhibitors and their double protonated forms (N9H+N10H+), in vacuum, at B3LYP/6-31G**. 

Molecular quantum parameters allow us to have more information about the global reactivity of inhibitors and their potential to interact with the metal. From Table 3, we notice that 5-MCTH has the highest EHOMO, indicating its ability to easily give electrons to electrophilic sites, and its higher reactivity with the lower ∆E gap, as ELUMO remains practically constant. The molecular reactivity could also be expressed by the global softness values and the fraction of electrons (∆N) transferred from the inhibitor to iron, knowing that the latter would preferentially interact with inhibitors having high S and low ƞ (ƞ= ∆E /2) values.

Table 3 Some quantum chemical parameters for the 3-MCTH, 4-MCTH and 5-MCTH neutral inhibitors, in vacuum (G) and in an aqueous solution (A), at B3LYP/6-31G**. 

With higher S and ∆N, 5-MCTH is more reactive than the other inhibitors, which is confirmed by its lower electronegativity and larger dipole moment that inform about the molecule polarity (25). This decrease in electronegativity led to a reduction in the global electrophilicity index, and, therefore, to an increase in the nucleophylic character, as the nucleophile is characterized by a lower ω value.

All energy values are in ev, the dipole moment (μ) is in Debye, S is the global softness in ev-1, and EEIS and EWL are the average experimental percentage inhibition efficiency (%), obtained, respectively, from electrochemical impedance spectroscopy and weight loss method, as given in (10,12).

As expected, the solvent addition substantially increases the molecules dipole moment, as their polarity changes seriously. It also reduces significantly the electrophilicity index (ω), while only a slight fluctuation has reached the other parameters, without any change in the reactivity classification order.

Combining these results with those of the frontier orbital picture analysis for molecules isolated and in the aqueous solution, we can conclude that 5-MCTH is more reactive in both gaseous and aqueous phases.

Table 4 Condensed Fukui functions on the 3-MCTH, 4-MCTH and 5-MCTH neutral inhibitors atoms, in vacuum (G) and in an aqueous solution (A), at B3LYP/6-31G**. 

Local reactivity and population analyses

Local reactivity based on FUKUI indices

It can be seen from Table 4, for the three molecules, that the atoms which are the most susceptible sites for the nucleophilic or electrophilic attacks belong to the thiadiazole and aryls rings, as they present fk+ and fk- high values, such as C1, C3, C6, C8, C11, N9, N10, C11 and C14.

Molecular Electrostatic Potential

The MEP surface picture of the three molecules, given in Fig. 4, reveals that the region of high negative charges (sites rich in electrons, represented by the red colour) is seen around the N9 and N10 nitrogen atoms that are susceptible to an electrophilic attack.

Figure 4 Molecular electrostatic potential contour map (MPE) for the 3-MCTH, 4-MCTH and 5-MCTH neutral inhibitors optimized structures and their double protonated N9H+N10H+ optimized forms, in vacuum (G), at B3LYP/6-31G**. (Red: strong negative electrostatic potential (EP); yellow: moderately negative EP; blue: strong positive EP; green: moderately positive EP). 

The protonation process of these two sites confirms this, since it weakens their electronic abundance, making them strongly positive, with high electrostatic potential energy value, represented by the blue colour. Moreover, the MEP surface also shows a substantial distribution of negative charge in the central oxygen atom, compared to the peripheral oxygen atom for 4-MCTH, as found for 2-MCTH in previous work, i.e. when the oxygen atoms number is odd(13) .

It should be noted, though, that C8 fk+ and fk- became more important for 4-MCTH and 5-MCTH, compared to those of 3-MCTH (about 0.15, instead of 0.04), and that C14 fk+ and fk- decreased remarkably from 0.05 (3-MCTH) to 0.02 (4-MCTH and 5-MCTH), while their values increased from 0.02 to 0.05, for C12 and C16. On the other hand, the sulphur site has important fk+ values, of the order of 0.18, for 4-MCTH and 5-MCTH, and of 0.17, for 3-MCTH, which implies that sulphur is more susceptible to receive electrons than to donate them to the metal surface. Concerning the solvent effect on these indices, we noted that water, as solvent, had no significant effect on the local reactivity.

As it can be seen, the results of the local reactivity based on Fukui indices corroborate those obtained by the HOMO and LUMO analysis.

Figure 5 NBO charges of the optimized neutral inhibitors 3-MCTH, 4-MCTH and 5-MCTH in vacuum, at B3LYP/6-31G**. 

Natural population atomic charges and natural bond orbital analyses

The distribution of the molecular charge in terms of NPA (Natural Charge Population Analysis) is shown in Fig. 5. The partial atomic charges of nitrogen take a value of the order of -0.3, while those of oxygen are comparable to -0.5, for the three inhibitors .

With regard to the NBO analyzes, Table 5 shows that the strongest intramolecular charge transfer concerns the πCC - π*CC transition within the same phenyl radical, with a stabilization energy of approximately 24 kcal/mol, which concerns πC12C22- π*C14C18, for 3-MCTH and 5-MCTH, and πC3C20 - π*C16C24, for 4-MCTH.

Table 5 Second Order Perturbation theory analysis of Fock matrix in NBO basis for 3-MCTH, 4-MTCH and 5- MCTH, at B3LYP/6-31G**. (E(2) means energy of hyper conjugative interactions). 

Regioselectivity of isolated protonated species and in an aqueous solution

By comparing the proton affinity (PA) at the different possible sites for 1-MCTH and 2-MCTH, we have shown that the protonation and double protonation on N9 and N10 nitrogen atoms are preferred to those on O26 and O27 oxygen atoms. Moreover, the aqueous solvation had a stabilizing effect on the protonation process, and preserved the same order of protonation preference. For these reasons, we have considered, for these three molecules, only the double protonation on nitrogen atoms.

The 3-MCTHN9H+N10H+, 4-MCTHN9H+N10H+ and 5-MCTHN9H+N10H+ optimized geometries are shown in Fig. 6. The evaluations of the proton affinity (26) show that the protonation process is exothermic, meaning that all inhibitors have a tendency for protonation. It should be also noted that the proton affinity (PA) became more important for 5-MCTHN9H+N10H+ and 4-MCTHN9H+N10H+.

Figure 6 Double protonated N-N species of 3-MCTHN9H+N10H+, 4- MCTHN9H+N10H+ and 5- MCTHN9H+N10H+, in vacuum (G) and in an aqueous solution (A), at B3LYP/6-31G**, and the corresponding proton affinity values (kcal/mol). 

Selected geometrical parameters

The geometric parameters selected for the double protonated species are also reported in Tables 1 and 2. The double protonation on N9-N10 resulting in n-MCTHN9H + N10H + (n = 3 - 5) did not show any significant variation in binding lengths, including those sensitive to double protonation in 1-MCTH and 2 -MCTH, such as C3-C6 and C11-C8, separating the aromatic rings from the thiadiazole ring, as well as C6-N9 and C8-N10, within the thiadiazole ring, and the interaomic distance between the two C24-C25 phenyls. In the same context, the double protonation on the two nitrogen atoms had a variation not exceeding 5° on the torsion angles between the two aromatic rings, expressed by C20-C3-C11-C22 or C16-C24-C25-C18, and the twist angle between the thiadiazole ring and the aromatic rings, expressed by C24-C3-C6-S7 or C25-C11-C8-S7.

As a result, and contrary to 1-MCTH and 2-MCTH, n-MCTHN9H + N10H (n = 3 - 5) does not appear to have a degree of planarity very different from its neutral counterparts (Fig. 7). The same observation can be made for the solvation of protonated double nitrogen species, since we have not found any significant effect on the degree of planarity, although it acts on the bond lengths.

Figure 7 Planarity influencing factor in 3-MCTH, 4-MCTH and 5-MCTH neutral inhibitors and in their double protonated forms, at B3LYP/6-31G**. The region between the dotted lines represents the metal surface. 

Quantum chemical parameters of protonated forms isolated and in aqueous solutions

Table 6 reports quantum chemical parameters values of n-MCTHN9H + N10H+ (n = 3 - 5) optimized structures, in gaseous (G) and in aqueous (A) solutions. We note that EHOMO and ELUMO are lower than those of the neutral forms, suggesting that protonation decreases the inhibitor tendency to donate electrons, and increases its tendency to accept electrons, which is confirmed by ΔN negative values, indicating that the electronic exchange is more favorable from the metal towards the protonated inhibitor. With a lower ΔE value and, thus, a higher global softness, S, the protonated species are more reactive than the neutral forms. The electrophilicity character also increases considerably with electronegativity.

With regard to the solvent effect on the protonated species molecular properties, the reactivity expressed by the protonated species ΔE or S is lower than that of the protonated species without solvent. Additionally, the electrophilicity character also decreases drastically, as well as electronegativity. It should also be noted that the ΔN values, which are negative for the isolated protonated forms, become positive under the solvent effect, reflecting the possibility of inversion on the exchange direction, i.e., from the metal to the inhibitor. In addition, the protonation reduces the HOMO orbital space, essentially on the aryl and thiadiazole ring side, thus reducing the cation ability to give electrons, while the same zone remains enriched by electrons in LUMO (Fig. 3).

We can conclude, in the light of the quantum parameters obtained for the n-MCTH (n = 3 - 5) protonated forms, that these are more reactive, in terms of electron acceptance, than the neutral forms; however, the solvent moderates their reactivity. This can be confirmed by the HOMO and LUMO orbitals shape.

Table 6 Some quantum chemical parameters for the 3-MCTH, 4-MCTH and 5-MCTH optimized protonated inhibitors, in vacuum (G) and in an aqueous solution (A), at B3LYP/6-31G**. 

All energy values are in ev, the dipole moment (μ) is in Debye, S is the global softness in ev-1, and EEIS and EWL are the average experimental percentage inhibition efficiency (%), obtained, respectively, from electrochemical impedance spectroscopy and weight loss method, as given in (10, 12).

Table 7 B3LYP/6-31G** interaction energies of the 3-MCTH, 4- MCTH and 5-MCTH neutral inhibitors optimized structures, in vacuum (G) and in an aqueous solution (A), with the iron surface. 

Metal-inhibitor interaction mechanism in vacuum and with the solvent effect

After having studied the competitiveness between active sites within the same inhibitor, the competitiveness between neutral and protonated forms, and the solvent effect, we were able to deepen our understanding of the inhibition mechanism through the DFT and dynamic Monte Carlo (MC) study of complexation between corrosion inhibitor and iron surface (27).

Table 8 B3LYP/6-31G** interaction energies of the neutral inhibitors 3-MCTH, 4- MCTH double nitrogen protonated forms, in vacuum (G) and in an aqueous solution (A), with the iron surface. 

Figure 8 B3LYP/6-31G** interaction energies of the 3-MCTH, 4-MCTH and 5-MCTH neutral inhibitors optimized structures and the optimized double nitrogen protonated forms, in vacuum (G) and in an aqueous solution (A). 

Complexation DFT calculation

In this part, we evaluated the interaction energy of neutral forms, and of their double nitrogen protonated forms, with the iron metal reduced to a single Fe atom. From the results collected in Tables 7 and 8 and illustrated as a histogram in Fig. 8, we noted, for each molecule, the importance of the interaction involving two nitrogen atoms, compared to that involving one of them, and the importance of the interaction involving a neutral form compared to that of a protonated form. Moreover, we noted the supremacy of the interaction of 5-MCTH neutral and protonated variants over that of their 3-MCTH and 4-MCTH counterparts, which in turn are superior to the interaction of 2-MCTH and 1-MCTH (13).

The comparison of the molecular quantum descriptors of the formed complexes is an additional confirmation (Table 9). Indeed, the complexes ΔE and hardness (ƞ) increase, and their electronegativity (χ) and electrophilicity (Ω) decrease, from 3-MCTH / (N9 N10) -Fe to 5-MCTH / (N9 N10) -Fe, confirming the stability of the latter. Table 8 results confirmed also the strength of the complexation due to neutral forms, compared to that due to protonated ones.

We also found that the lengths of the N-Fe bonds decreased with the increase in interaction force, particularly when the Fe atom was simultaneously bonded to N9 and N10. Indeed, Fig. 9 shows this decrease, ranging from 1.953A° to 1.897A° (taken on the average of the two NH lengths), in 3-MCTH and 5-MCTH, respectively, for the neutral forms.

Fig. 9 also shows HOMO, LUMO and total electron density for Fe inhibitor complexes. It confirms the stabilizing effect of the Fe-(NN) interaction in the complex, compared to the Fe-N interaction. However, the supremacy of the interaction of 5-MCTH, compared to the two others, remains hardly perceptible in a visual way, if we exclude the slight reduction of the molecular orbital HOMO or LUMO, from 3-MCTH to 5-MCTH.

Table 9 Some optimized quantum chemical parameters for the Fe-inhibitor complexes and Fe-protonated inhibitor complexes, in vacuum (G). 

Dynamic Monte Carlo study of complexation

We point out first that the Dynamic Monte Carlo was performed to look for the lowest configuration adsorption energy of the inhibitor molecules interactions with the iron surface. It provides some vital parameters, such as total surface energy (Etotal), which is the summation of all the adsorbate components energies (28)(29). It also furnishes the inhibitor adsorption energy (Eads) and dEad/dNi energy required when one mole of it is adsorbed onto the metallic surface (30). This informs about the interaction strength between the inhibitor and the metallic surface, knowing that an inhibitor with more negative Eads or dEad/dNi values acts as better adsorbate onto the Fe(111) surface (31).

Figure 9 HOMO and LUMO optimized structures and total charge density (T. D.) of the Fe-inhibitor complexes in vacuum . 

It is quite clear from Table 10 that the total energies of the neutral and protonated forms in vacuum and aqueous phases increased in the order: 5-MCTH <4-MCTH < 3-MCTH. This ordering indicates that 5-MCTH is the most stable among the three inhibitors. Furthermore, we notice that all the Eads values are negative, which suggests that the neutral and protonated forms can strongly adsorb onto the Fe (111) surface. It should be also noted that 5-MCTH-Fe and its protonated 5-MCTHN9H+N10H+-Fe counterpart have the lowest adsorption energy value in both vacuum (G) and aqueous phases, which confirms the supremacy of their interactions over their 3-MCTH-Fe and 4-MCTH-Fe homologues complexes.

Table 10 Evaluation by Monte Carlo simulation, in vacuum and aqueous phases, of the descriptors of the adsorption configurations onto Fe (111), for 3-MCTH, 4-MCTH and 5-MCTH and their doubly protonated forms at the nitrogen level (all values in kcal / mol). 

The top and side views of the equilibrium adsorption configurations of 3-MCTH, 4-MCTH and 5-MCTH inhibitors and their protonated forms, onto the Fe (111) surface, in both vacuum (G) and in the aqueous phase (A) (100 H2O), are shown in Fig. 10. Here, we notice easily that the maximum surface recovery is achieved by the parallel adsorption of all inhibitor molecules or protonated forms onto the Fe (111) surface, in both vacuum (G) and the aqueous phases (A). That is why the degree of planarity is a key element in the adsorption efficiency.

From Table 11, collecting the distance between nitrogen atoms and the nearest Fe atom of the Fe (111) surface, it can be emphasized that 5-MCTH, in both neutral and protonated forms, has the lowest distance between nitrogen and the nearest Fe-N surface.

Table 11 Distance between nitrogen atoms and the nearest Fe atom of the Fe (111) surface. 

*When N9 and N10 are not equivalent, the distance considered is that separating the nearest Fe atom from the middle of the NN bond.

Figure 10 Top and side views of the stable adsorption configurations of 3-MCTH, 4-MCTH and 5-MCTH inhibitors and their protonated forms onto the Fe(111) iron surface, in vacuum (G) and in the aqueous phase (A) (100 H2O). 

Conclusion

By combining local and global electronic properties for molecules isolated and in an aqueous solution, we can conclude that 5-MCTH is more reactive in both vacuum and aqueous phase, due to the significant involvement of aryl and thiadiazole rings, in the process of electron donation and acceptance. Moreover, based on the quantum parameters evaluation, the protonated forms are more reactive, in terms of electron acceptance (with a notable superiority of the 5-MCTHN9H +N10H +), than the neutral forms; however, the solvent moderates their reactivity.

In deepening our understanding of the Fe (111) surface complexation mechanism, the results of both DFT and molecular dynamics Monte Carlo approaches showed the importance of the interaction involving a neutral form compared to that of a protonated form, and the supremacy of the interaction of 5-MCTH neutral and protonated variants, compared to their 3-MCTH-Fe and 4-MCTH-Fe complexes homologues.

Authors’ contribution

A. Mahsoune: collected the data; contributed with data or analysis tools; wrote the paper. K. Sadik: contributed with data or analysis tools; performed the analysis; wrote the paper. S. Byadi: collected the data. A. Aboulmouhajir: conceived and designed the analysis; contributed with data or analysis tools; performed the analysis; wrote the paper.

References

1. Liao LL, Mo S, Luo HQ, et al. Corrosion protection for mild steel by extract from the waste of lychee fruit in HCl solution : experimental and theoretical studies. J Colloid Interf Sci. 2018;520:41-49. Doi: https://doi.org/10.1016/j.jcis.2018.02.071 [ Links ]

2. Lgaz H, Chung IM, Salghi R, et al. On the understanding of the adsorption of Fenugreek gum on mild steel in an acidic medium: Insights from experimental and computational studies. Appl Surf Sci. 2019;463:647-58. Doi: https://doi.org/10.1016/j.apsusc.2018.09.001 [ Links ]

3. Hamadi L, Mansouri S, Oulmi K, et al. The use of amino acids as corrosion inhibitors for metals: A review. Egyp J Petr. 2018;27:1157-1165. Doi: https://doi.org/10.1016/j.ejpe.2018.04.004 [ Links ]

4. Quraishi MA, Chauhan DS, Saji VS. Heterocyclic corrosion inhibitors. In: Heterocyclic Organic Corrosion Inhibitors. Elsevier; 2020. P. 1-19. Doi: https://doi.org/10.1016/b978-0-12-818558-2.00001-1. [ Links ]

5. Zhang QH, Hou BS, Zhang GA. Inhibitive and adsorption behavior of thiadiazole derivatives on carbon steel corrosion in CO2-saturated oilfield produced water: Effect of substituent group on efficiency. J Colloid Interf Sci. 2020 Jul 15;572:91-106. Doi: https://doi.org/10.1016/j.jcis.2020.03.065 [ Links ]

6. Bentiss F, Lebrini M, Lagrenée M, et al. The influence of some new 2,5-disubstituted 1,3,4-thiadiazoles on the corrosion behaviour of mild steel in 1 M HCl solution: AC impedance study and theoretical approach. Port Electrochim Acta. 2007;52(24):6865-72. Doi: https://doi.org/10.1016/j.electacta.2007.04.111 [ Links ]

7. Bentiss F, Mernari B, Traisnel M, et al. On the relationship between corrosion inhibiting effect and molecular structure of 2 , 5-bis ( n -pyridyl ) -1 , 3 , 4-thiadiazole derivatives in acidic media : Ac impedance and DFT studies. Corr Sci. 2011;53:487-95. Doi: https://doi.org/10.1016/j.corsci.2010.09.063 [ Links ]

8. Bentiss F, Lebrini M, Lagrenée M, et al. The influence of some new 2 , 5-disubstituted 1 , 3 , 4-thiadiazoles on the corrosion behaviour of mild steel in 1 M HCl solution : AC impedance study and theoretical approach. Port Electrochim Acta. 2007;52:6865-72. Doi: https://doi.org/10.1016/j.electacta.2007.04.111 [ Links ]

9. Bentiss F, Mernari B, Traisnel M, et al. On the relationship between corrosion inhibiting effect and molecular structure of 2,5-bis(n-pyridyl)-1,3,4-thiadiazole derivatives in acidic media: Ac impedance and DFT studies. Corr Sci. 2011 Jan 1;53(1):487-95. Doi: https://doi.org/10.1016/j.corsci.2010.09.063 [ Links ]

10. Lebrini M. Synthèses et études physicochimiques de nouveaux thiadiazoles inhibiteurs de corrosion de l'acier en milieu acide. L'Université des Sciences et Tehnologies de Lille. 2005. [ Links ]

11. Bentiss F, Lebrini M, Vezin H, et al. Enhanced corrosion resistance of carbon steel in normal sulfuric acid medium by some macrocyclic polyether compounds containing a 1,3,4-thiadiazole moiety: AC impedance and computational studies. Corr Sci. 2009;51(9):2165-73. Doi: http://dx.doi.org/10.1016/j.corsci.2009.05.049 [ Links ]

12. Bentiss F, Lebrini M, Vezin H, et al. Experimental and theoretical study of 3-pyridyl-substituted 1,2,4-thiadiazole and 1,3,4-thiadiazole as corrosion inhibitors of mild steel in acidic media. Mat Chem and Phys. 2004;87(1):18-23. Doi: https://doi.org/10.1016/j.matchemphys.2004.05.040 [ Links ]

13. Mahsoune A, Sadik K, Belghiti ME, et al. Toward a Theoretical Understanding of the Corrosion Inhibitive Performance on Iron Surface by Some Macrocyclic Polyether Compounds Containing 1 , 3 , 4-thiadiazole Entity. Int J Electrochem Sci. 2018;13:8396-427. Doi: https://doi.org/10.20964/2018.09.13 [ Links ]

14. Hertwig RH, Koch W. On the parameterization of the local correlation functional. What is Becke-3-LYP? Chem Physics Lett. 1997;268(5-6):345-51. Doi: https://doi.org/10.1016/S0009-2614(97)00207-8 [ Links ]

15. Frisch MJ, Trucks GW, Schlegel HB, et al. Gaussian 09 Citation. Gaussian.com (2009). http://gaussian.com/g09citation. [ Links ]

16. Cammi R. The Quantum Chemical Study of Chemical Reactions at Extreme High Pressure by Means of the Extreme-Pressure Polarizable Continuum Model (Internet(. 1st ed. Vol. 13. Annual Reports in Computational Chemistry. Elsevier B.V. 2017. P 117-135. Doi: http://dx.doi.org/10.1016/bs.arcc.2017.06.001 [ Links ]

17. Hachim ME, Sadik K, Byadi S, et al. Ab initio study on the six lowest energy conformers of iso-octane : conformational stability , barriers to internal rotation , natural bond orbital and first-order hyperpolarizability analyses , UV and NMR predictions, spectral temperature sensitivity. J Molec Mod. 2019;25:254. Doi: https://doi.org/10.1007/s00894-019-4105-5 [ Links ]

18. Sanaei Z, Bahlakeh G, Ramezanzadeh B, et al. Application of green molecules from Chicory aqueous extract for steel corrosion mitigation against chloride ions attack ; the experimental examinations and electronic / atomic level computational studies. J Molec Liq. 2019;290:111-176. Doi: https://doi.org/10.1016/j.molliq.2019.111176 [ Links ]

19. Ramezanzadeh M, Bahlakeh G, Sanaei Z, et al. Corrosion inhibition of mild steel in 1 M HCl solution by ethanolic extract of eco-friendly Mangifera indica (mango) leaves: Electrochemical, molecular dynamics, Monte Carlo and ab initio study. Appl Surf Sci. 2019;463:1058-77. Doi: https://doi.org/10.1016/j.apsusc.2018.09.029 [ Links ]

20. Singh A, Ansari KR, Haque J, et al. Effect of electron donating functional groups on corrosion inhibition of mild steel in hydrochloric acid: Experimental and quantum chemical study. J Taiwan Instit Chem Engin. 2018;82:233-51. Doi: https://doi.org/10.1016/j.jtice.2017.09.021 [ Links ]

21. Dehghani A, Bahlakeh G, Ramezanzadeh B, et al. A combined experimental and theoretical study of green corrosion inhibition of mild steel in HCl solution by aqueous Citrullus lanatus fruit (CLF) extract. J Molec Liq. 2019; 279:603-624. Doi: https://doi.org/10.1016/j.molliq.2019.02.010 [ Links ]

22. Ramezanzadeh M, Bahlakeh G, Ramezanzadeh B. Elucidating detailed experimental and fundamental understandings concerning the green organic-inorganic corrosion inhibiting molecules onto steel in chloride solution. J Molec Liq. 2019;290:111-212. Doi: https://doi.org/10.1016/j.molliq.2019.111212 [ Links ]

23. S.D. Accelrys Software Inc. SD Accelrys Software Inc. 2009. [ Links ]

24. Shahraki M, Dehdab M, Elmi S. Theoretical studies on the corrosion inhibition performance of three amine derivatives on carbon steel : Molecular dynamics simulation and density functional theory approaches. J Taiwan Instit Chem Engin. 2016;62:313-321. Doi: http://dx.doi.org/10.1016/j.jtice.2016.02.010 [ Links ]

25. Behzadi H, Roonasi P, Momeni MJ, et al. A DFT study of pyrazine derivatives and their Fe complexes in corrosion inhibition process. J Mol Struc. 2015;1086:64-72. Doi: http://dx.doi.org/10.1016/j.molstruc.2015.01.008 [ Links ]

26. El Ibrahimi B, Soumoue A, Jmiai A, et al. Computational study of some triazole derivatives (un- and protonated forms) and their copper complexes in corrosion inhibition process. J Mol Struc. 2016;1125:93-102. Doi: http://dx.doi.org/10.1016/j.molstruc.2016.06.057 [ Links ]

27. Verma C, Lgaz H, Verma DK, et al. Molecular dynamics and Monte Carlo simulations as powerful tools for study of interfacial adsorption behavior of corrosion inhibitors in aqueous phase: A review. J Mol Liq. 2018,260:99-120. Doi: https://doi.org/10.1016/j.molliq.2018.03.045 [ Links ]

28. Madkour LH, Kaya S, Guo L, et al. Quantum chemical calculations, molecular dynamic (MD) simulations and experimental studies of using some azo dyes as corrosion inhibitors for iron. Part 2: Bis-azo dye derivatives. J Mol Struc. 2018;1163:397-417. Doi: https://doi.org/10.1016/j.molstruc.2018.03.013. [ Links ]

29. Guo L. Toward understanding the anticorrosive mechanism of some thiourea derivatives for carbon steel corrosion : A combined DFT and molecular dynamics investigation. J Colloid Interf Sci. 2017;506:478-485. Doi: https://doi.org/10.1016/j.jcis.2017.07.082 [ Links ]

30. Madkour LH, Kaya S, Obot IB. Computational, Monte Carlo simulation and experimental studies of some arylazotriazoles (AATR) and their copper complexes in corrosion inhibition process. J Mol Liq. 2018;260:351-74. Doi: https://doi.org/10.1016/j.molliq.2018.01.055 [ Links ]

31. Belghiti ME, Echihi S, Mahsoune A, et al. Piperine derivatives as green corrosion inhibitors on iron surface; DFT, Monte Carlo dynamics study and complexation modes. J Mol Liq. 2018;261(2017):62-75. Doi: https://doi.org/10.1016/j.molliq.2018.03.127 [ Links ]

Received: November 24, 2019; Accepted: July 15, 2021

Corresponding author: aboulmouhajir@gmail.com

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License