Next Article in Journal
Polymer Supercritical CO2 Foaming under Peculiar Conditions: Laser and Ultrasound Implementation
Next Article in Special Issue
In Vitro Optical and Physical Stability of Resin Composite Materials with Different Filler Characteristics
Previous Article in Journal
Effect of UV-C Radiation on 3D Printed ABS-PC Polymers
Previous Article in Special Issue
Influence of Surface Discharge on Resin Degradation in Decay-like Fracture of Composite Insulators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Combined Experimental and Computational (DFT, RDF, MC and MD) Investigation of Epoxy Resin as a Potential Corrosion Inhibitor for Mild Steel in a 0.5 M H2SO4 Environment

1
Laboratory of Organic Chemistry, Bioorganic and Environment, Chemistry Department, Faculty of Sciences, Chouaib Doukkali University, BP 20, El Jadida 24000, Morocco
2
Laboratory of Organic Chemistry, Catalysis and Environment, Department of Chemistry, Faculty of Sciences, Ibn Tofail University, BP 242, Kenitra 14000, Morocco
3
Euromed Polytechnic School, Euromed Research Center, Euromed University of Fes, Fès 30000, Morocco
4
Laboratory of Advanced Materials and Process Engineering, Faculty of Sciences, Ibn Tofail University, BP 242, Kenitra 14000, Morocco
5
Chemistry Department, Faculty of Science, Al Azhar University-Gaza, Gaza BP 1277, Palestine
6
Chemistry Department, Faculty of Science, King Abdulaziz University, BP 42805, Jeddah 21589, Saudi Arabia
7
Department of Chemistry, Faculty of Natural and Mathematics Science, University of Prishtina, 10000 Prishtina, Kosovo
*
Author to whom correspondence should be addressed.
Polymers 2023, 15(8), 1967; https://doi.org/10.3390/polym15081967
Submission received: 13 March 2023 / Revised: 14 April 2023 / Accepted: 17 April 2023 / Published: 21 April 2023
(This article belongs to the Special Issue Resin-Based Polymers and Composite Materials)

Abstract

:
In this work, a tetrafunctional epoxy resin entitled 2,3,4,5-tetraglycidyloxy pentanal (TGP) was tested and investigated as a potential corrosion inhibitor for mild steel (MS) in 0.5 M H2SO4 solution. The corrosion inhibition process for mild steel was employed alongside various techniques, such as potentiodynamic polarization (PDP), electrochemical impedance spectroscopy (EIS), temperature effect (TE), scanning electron microscopy (SEM), energy-dispersive X-ray spectroscopy (EDS) and theoretical approaches (DFT, MC, RDF and MD). Further, the corrosion efficacies obtained at the optimum concentration (10−3 M of the TGP) were 85.5% (EIS) and 88.6% (PDP), respectively. The PDP results indicated that the TGP tetrafunctional epoxy resin acted the same as an anodic inhibitor type in 0.5 M H2SO4 solution. SEM and EDS analyses found that the protective layer formed on the MS electrode surface in the presence of TGP could prevent the attack of the sulfur ions. The DFT calculation provided more information regarding the reactivity, geometric properties and the active centers of the corrosion inhibitory efficiency of the tested epoxy resin. RDF, MC and MD simulations showed that the investigated inhibitory resin have a maximum inhibition efficiency in 0.5 M H2SO4 solution.

1. Introduction

Metallic steel is widely employed and investigated in various industrial fields of application, including aeronautics, automotives, railways, metal construction, etc., due to its advantageous availability and excellent mechanical resistance [1,2]. In addition, several methods have been investigated and developed in order to control and prevent the corrosion process in destructive environments, including sacrificial anode protection, galvanization, and the application of polymer composite coatings and/or the use of anticorrosion protection [3,4,5]. One of the most cost-effective and practical approaches is the implementation of organic inhibitory resins. The incorporation of polyepoxide inhibitory resins that possess heterocyclic compounds, epoxy systems, heteroatoms (phosphorus, sulfur, nitrogen and oxygen) and electronegative polar groups is the method used most in order to slow and stop the corrosion of metallic samples in aggressive environments [6,7]. Indeed, they act through an adsorption mechanism on the metal substrate area by forming a protective layer barrier against corrosion in aggressive environments. Among these various compounds, epoxide resins, which result from the reaction between glycidol and halogenated compounds and/or the condensation of epichlorohydrin with compounds, possess mobile protons due to the presence of heterocyclic compounds, aromatic rings and epoxy systems in their structure [8]. Several epoxide resins have been reported and investigated for their anticorrosive protection of metallic samples in aggressive environments [9]. Furthermore, several publications on the evaluation and performance of epoxide resins and their corrosion inhibition activity have investigated the fact that these macromolecules are mostly employed owing to their important inhibitory efficiency [10].
Recently, and based on the great advances made in supercomputers, quantum chemical methods, such as density functional theory (DFT), molecular dynamic simulation and Monte Carlo methods, have become very useful tools in the determination of the molecular structure, electronic structure and molecular reactivity of chemical compounds. Therefore it has become common practice to perform quantum chemistry calculations in order to predict the molecular corrosion inhibition performance and efficiency of different molecules and to understand their inhibitory mechanisms.
In this potential study, an inhibitory resin entitled 2,3,4,5-tetraglycidyloxy pentanal was employed; it performed excellently in terms of providing anticorrosion protection to steel substrate electrodes in 0.5 M H2SO4, including when analyzed using different characterization techniques, such as potentiodynamic polarization, electrochemical impedance spectroscopy, temperature effect, scanning electron microscopy and X-ray energy-dispersive spectroscopy [11]. To know and understand the interaction mechanism between the investigated inhibitory resin and the surface of the metallic substrates, several global chemical quantum descriptors, such as the electronegativity, the absolute hardness, the electrophilicity index, the softness, the fraction of electrons transferred, the condensed nucleophiles and electrophiles, the Fukui functions, the local softness, the local electrophilicity, the dual Fukui functions, the dual softness and the dual philicity, were used and discussed further. In addition, the DFT method is an important and effective approach to describing the structural and geometric behaviors of a huge number of inhibitors and their role in the protection of several metals [12]. Additionally, the RDF, MC and MD approaches were performed and investigated in order to evaluate the adsorption energies and the adsorption configuration of the inhibitory resin on the Fe(110) area in the simulated aggressive environment [13].

2. Experimental Part

2.1. Inhibitor Used Material and Characterization Techniques

The inhibitory resin, namely 2,3,4,5-tetraglycidyloxy pentanal, was synthetized, identified and confirmed by Hsissou et al., according to the procedures illustrated several of the previously published studies in this area [6,14]. In addition, the chemical composition of the mild steel (MS) samples used in the experimental section was C (3.66 ± 0.03%), O (22.62 ± 0.06%), S (5.89 ± 0.03%), Cl (0.24 ± 0.01%) and Fe (67.60 ± 0.14%). Moreover, the MS samples with a 1 cm × 1 cm surface were abraded with various grades of sandpaper from 180 to 2000, rinsed with bi-distilled water, washed using acetone, and then dried at room temperature. Then, to protect the MS surface, a TGP concentration range of 10−3 to 10−6 M in 0.5 M H2SO4 solution was maintained in order to obtain the potentiodynamic polarization curves (PDP) and electrochemical impedance spectroscopy (EIS) measurements. In addition, the electrochemical measurements were investigated using a Potentiostat/Galvanostat/ZRA (Gamry apparatus). Furthermore, the electrochemical cell was composed of the mild steel, which was used as the work electrode, the platinum, which was used as the counter electrode, and the saturated calomel electrode, which was used as the reference electrode. The open circuit potential (OCP) was maintained for 30 min in each assay after reaching the equilibrium state. Then, the PDP plots were measured at 0.3 mV/s when the potential range was between −900 and −100 mV. The EIS measurements were carried out when the OCP was in the frequency range of 100 KHz to 0.1 Hz with an AC voltage of 10 mV. In addition, the surfaces of the MS samples immersed in 0.5 M H2SO4 solution were uninhibited and inhibited with 10−3 M of TGP; this was performed using a scanning electron microscope and energy-dispersive X-ray spectroscopy. In addition, the SEM technique was investigated for its potential ability to create images of the uninhibited and inhibited mild steel surface in the presence of TGP at 10−3 M, using a JEOL-JSM-5500-type microscope. Moreover, EDS analysis was used in order to obtain the elemental chemical composition. Furthermore, the measurement conditions were as follows: the magnification was x2000, the process time was T4, the acceleration voltage was 20.00 kV, the measurement was first, the real time was 38.74 s, the live time was 30.00 s, the count rate was 29,643.00 CPS and the dead time was 22.00%. Additionally, the corrosion efficiencies for the potentiodynamic polarization curves and for electrochemical impedance spectroscopy were measured according to Equations (1) and (2), as follows:
% η PDP = i corr 0   -   i corr i corr 0 × 100
% η EIS = R ct   -   R ct 0 R ct × 100  
where R ct 0 and R ct are the charge transfer resistance for the uninhibited and inhibited MS with different concentrations of TGP, respectively. In addition, i corr 0 and i corr are the corrosion current densities for the uninhibited and inhibited MS in the presence of TGP at different concentrations, respectively.

2.2. DFT Details

In this work, density functional theory (DFT) was used in order to analyze the electronic behaviors of the investigated tetrafunctionnal epoxy resin (TGP) and its protonated form (Scheme 1). Optimization of the proposed geometries of all species was performed by applying Becke’s three-parameter Lee–Yang–Par correlation functional (the B3LYP hybrid functional) [15] alongside the 6-31+G(d,p) basis set. Since the electrochemical corrosion takes place in the aqueous phase and because this study aimed to model the solvent effect, all the calculations were performed in an aqueous environment using the polarizable continuum model (PCM) of solvation [16]. The PCM is a method commonly used to model solvation effects in computational chemistry. It is necessary to consider each solvent molecule as a separate molecule. Two types of PCM have been popularly used: the dielectric PCM (D-PCM), in which the continuum is polarizable, and the conductor-like PCM (C-PCM), in which the continuum is like a conductor. Because corrosion occurs in an aqueous environment, the calculation is more interesting and correlated with the experimental results when the solvent effects are considered.
The minimum value of the potential energy areas of these samples was found using the absence of imaginary frequencies. The molecular modeling studies were carried out by using Gaussian 09 [17]. Gaussview [18] programs were used to build up the input files, to visualize the highest occupied and lowest occupied molecular orbitals (HOMO and LUMO), and to visualize the electrostatic potential maps (ESP). The energies of HOMO (I = −EHOMO) and the LUMO (A = −ELUMO) were investigated in order to illustrate the energy gap (∆Egap) of the neutral and protonated species (∆Egap = ELUMOEHOMO). Several quantum global reactivity descriptors of the investigated inhibitor, including electronegativity (χ), absolute hardness (η), the electrophilicity index (ω), softness (S) and the fraction of electrons transferred (∆N), were evaluated in terms of the vertical ionization potential (Iv) and vertical electron affinity (Av); these were calculated as functions of Equation (3), as follows [19,20,21]:
I v E N 1 E N   and   A v E N E N + 1 μ o = χ = I v + A v 2 ,   η = I v A v 2   , S = 1 η , ω   = μ o 2 2 η   , Δ N = Φ F e χ i n h 2 η F e + η i n h Δ E b d = η i n h 4
In the above set of equations, E(N), E(N − 1) and E(N + 1), and ΦFe, η F e and ΔEb−d, are the total energy of the optimized structures, the cationic and anionic energies at the optimized neutral geometries, the work function (4.82 eV), the hardness of Fe(110) (0.0 eV) and the electronic back-donation energy, respectively [22]. If ΔEb−d < 0, the mechanism is energetically favored, otherwise the process is energetically unfavored [23].
Fukui functions can give a useful indication of the local reactivity centers regarding nucleophilic and electrophilic attacks when using the investigated inhibitory resin [24]. To compute the Fukui indices, single-point energy calculations for the neutral, cation and anion radicals that comprise the geometrical structure of the neutral species were performed in order to generate the required wave functions. The condensed nucleophilic ( f k + ) and electrophilic ( f k ) Fukui functions, the local softness s k ± , the local electrophilicity ( ω k ± ), the dual Fukui functions Δ f r , the dual softness and dual philicity, which can give a numerical value for the chemical reactivity of all sites in the molecule, were evaluated using the following expressions (Equation (4)) [25]:
f k + = q k N + 1 q k N f k = q k N q k N 1 s k ± = S f k ± r ,   ω k ± = ω f k ± r f k 2 = Δ f r = f k + f k Δ s k = s k + s k Δ ω k = ω k + ω k

2.3. Monte Carlo and Molecular Dynamics Computational

Using the Material Studio package’s Forcite module, the ability of neutral and protonated inhibitory resins to adhere to the area of iron was studied via MC and MD simulations. This was performed in a simulation box (2.4823 nm × 2.4823 nm × 1.824 nm + 3.5 nm vacuum layer) with regular conditions at the edges. For the system being studied, the Fe(110) area with a slab of 15 was chosen; this is because it is the most stable of the three common iron areas [26]. Details about the MS simulations were described previously [27,28]. The slab model in the simulation contained the following: 1 inhibitory resin (in neutral or protonated state), 900 water molecules, and 10 hydroniums plus 5 sulfate ions. In short, these simulations used the COMPASSIII force field [29], according to the following parameters: the NVT ensemble, 295 K, a time step of 1 fs, and a simulation time of 0.8 ns [30].

3. Results and Discussion Part

3.1. Electrochemical Investigation

In order to determine the interaction between the tetrafunctionnal inhibitory resin and the surface of the MS samples uninhibited and inhibited with TGP at different concentrations in 0.5 M H2SO4 solution at 298 K, PD was used and investigated, as illustrated in Figure 1 and Table 1. The data listed in Table 1 demonstrate that the icorr decreases to lower values that ranging between 1850 and 211 µA cm−2 with as the TGP concentration increases (from 10−6 to 10−3 M) [6]. In addition, this diminishes due to the development of a covalent bond between the methoxy group, the ketone function, the oxiran systems and the free electrons of oxygen heteroatoms that possess the vacant iron orbitals of the MS surface [7]. These data demonstrate that the TGP inhibitor, at a low concentration, was able to form a protective layer on the surface of the metallic samples by diminishing the number of active sites. These data also demonstrate that it can be employed as an excellent inhibitory resin for the mild steel electrode in 0.5 M H2SO4 solution. Additionally, all the corrosion potential displacements (ΔEcorr) for MS uninhibited and inhibited with different concentrations of TGP (97 to 112 mV) are higher than ±85 mV, indicating that the TGP molecule acted as an anodic inhibitor type [31]. In addition, after carefully examining the results illustrated in Table 1, it is evident that the cathodic and anodic slopes are modified after the addition of the TGP to the aggressive environment, compared to the blank [32]. Further, from the data shown in Table 1, it can be seen that, as the concentration of TGP increases, the inhibitory efficacy increases and reaches a high value of 88.6% at a high concentration, resulting in the TGP potentially being adsorbed on the MS electrode surface.
In order to determine and understand the corrosion inhibition of the metallic surface, TGP was employed and investigated. Figure 2 illustrates the Nyquist plots for the anticorrosive protective efficiency of TGP in the 0.5 M H2SO4 environment, uninhibited and inhibited with TGP at different concentrations. Further, a careful examination of the Nyquist plots reveals that the corrosion processes are the same for both the uninhibited and inhibited MS surfaces. According to the data displayed in Figure 2, the charge transfer resistance and the diameter of the semi-circles increase after the addition of different concentrations of TGP in 0.5 M H2SO4 solution, which indicates that the TGP can be adsorbed into the metallic substrates by forming a protective layer and by preventing the dissolution of the metal [33]. According to the data illustrated in Figure 2, the Bode plots increase as the concentration of TGP increases, and also decreases according to the frequency. This result implies that the TGP contributes to the formation of a protective layer on the MS surface. In addition, electrochemical data were simulated by using the equivalent circuit (EC), including the electrolyte resistance (Re), the constant-phase element (CPE(Cdl and ndl)), the charge transfer resistance (Rct) and the inductance element (RL and L), as illustrated in Figure 2. A careful examination of the data grouped in Table 2 reveals that the charge transfer resistance values increase as the TGP concentration increases from 11.7 to 81.37 Ω cm2, and also that the inductive element values increase from 2.17 to 22.35 Ω cm2 for the MS surface uninhibited and inhibited with 10−3 M of TGP; this reflects the highest blockage of sulfur ions by using TGP at the MS/0.5 M H2SO4 interface [34]. These data suggest that the inhibitory resin is able to protect mild steel substrates to an excellent standard by preventing the ion’s diffusion. From the data listed in Table 2, it is evident that a strong agreement exists between the Q values and the anticorrosive inhibitory efficacy. For the mild steel electrode, the Q values decrease as the TGP concentration increases, but only when compared to the blank; this is due to an increase in the thickness of the protective layer and/or a decrease in the dielectric constant. This protective property is predominantly shown in the differences between the surface states of the metallic substrates in terms of the structure and roughness of the samples after anticorrosive protection. For the inhibitory resin investigated, a more pronounced effectiveness is achieved with 10−3 M of TGP, which is in excellent correlation with the data obtained using the potentiodynamic polarization technique. For all concentrations, the ndl values for the MS surfaces uninhibited and inhibited with inhibitory resin in 0.5 M H2SO4 solution are less than 1, resulting in the pseudo-capacitive system of the CPE. These data suggest that the inhibitory resin offers excellent anticorrosion protection compared to other epoxy resins that have been investigated for their potential anticorrosion protection [35].

3.2. SEM/EDS Investigation

In order to display and evaluate the adsorption phenomena between the surface of metallic substrates and the TGP employed, SEM and EDS characterization analyses were used and discussed further. Figure 3 displays the SEM morphology of the uninhibited and inhibited MS surfaces immersed in 0.5 M H2SO4, with inhibitory resin at a higher concentration (10−3 M). In addition, the surface of the MS electrode shows cracks, holes and a severe surface. However, the MS samples display the protective layer, and are smooth with fewer cracks and holes with 10−3 M of TGP, suggesting that the TGP is a good inhibitor and has a more significant protective effect [6]. Additionally, the EDS characterization data of the uninhibited and inhibited surfaces of the MS samples with 10−3 M of TGP in 0.5 M H2SO4 are shown in Figure 4 and Table 3. Considering the data illustrated in Figure 4 and after a careful analysis of Table 3, it can be remarked that the percentage of Fe present on the MS electrode surface increases from 67.60 to 84.24%, and that the elemental peak of Sulfur (S) decreases from 5.89 to 0.97% after the addition of the inhibitory resin TGP. In addition, the elemental peaks in the content of oxygen and carbon increase with the addition of the 10−3 M of TGP (from 9.51 to 34.66% (for oxygen) and from 5.04 to 6.89% (for carbon)) due to the heteroatoms of oxygen and carbon that are present in the inhibitory resin studied; this suggests that the TGP reduces the density of the active sites on the surface of the metallic substrates by forming a protective layer, thus delaying the corrosion inhibition.

3.3. Adsorption Study

In order to determine and evaluate the adsorption mechanism mode that functions between the inhibitory resin and the MS surface in the 0.5 M H2SO4 solution, the experimental results obtained using the PDP technique were evaluated by applying the various adsorption isotherms equations and ascertaining the best fit, which in this case, was the Langmuir adsorption mechanism mode (Figure S1). According to the results found using the PDP plots, it appeared that the TGP was adsorbed on the surface of the MS samples via the formation of a protective layer, which protects the metallic substrates area in an aggressive environment. In addition, Table 4 illustrates the various parameters determined using the PDP technique applied for the studied inhibitory resin. According to data obtained using the PDP technique, the coefficient value (R2) is 1, which is in good agreement with the Langmuir adsorption mechanism. In addition, this information indicates that the inhibitory resin TGP could be adsorbed on the MS surface, due to the fact that the best mode could be used in order to develop a certain experience. Moreover, the Gibbs energy variation (ΔGads) values for adsorption offer more data concerning the adsorption mechanism performed by the inhibitory resin and the MS surface. Additionally, a ΔGads value that is higher than −40 KJ/mol indicates chemical adsorption. However, a ΔGads value that is less than −20 kJ/mol suggests the physical adsorption mode [6]. Furthermore, the Langmuir adsorption isotherm mode and the Gibbs energy variation for adsorption can be determined by using Equations (5) and (6), as follows [13]:
C inh θ = 1 K ads + C inh
Δ G ads = R × T × Ln   ( 55 . 5 × K ads )
where Kads, θ, Cinh, T and R are the adsorption equilibrium constant, the surface coverage degree, the inhibitory resin concentration, the employed temperature and the perfect gas constant, respectively.
A careful analysis of Table 4 reveals that the ΔGads value is equal to −44.74 KJ mol−1, which suggests that the TGP inhibitory resin could be adsorbed on the MS area through the chemical adsorption isotherm model. Therefore, the chemical adsorption of the inhibitory resin studied could be due to the nature the bond between the Fe ions and the epoxy systems, and the methoxy function and the free electron of the oxygen heteroatoms present in the TGP inhibitory resin employed.

3.4. Temperature Effect and Kinetics Parameters

In order to determine and evaluate the adsorption mechanism, the temperature effect of the system was one of the essential factors used to study the corrosion inhibition of an uninhibited and inhibited steel electrode in a 0.5 M H2SO4 medium with an optimal concentration of TGP at different temperatures. In addition, in order to study the temperature effect on the inhibitory resin employed, potentiodynamic polarization analysis was carried out at varying temperatures, as follows: 298, 308, 318 and 328 K. Figure 5 illustrates the temperature effect of the PDP plots obtained in an aggressive environment for MS uninhibited and inhibited with 10−3 M of inhibitory resin and tested at different temperatures. The different electrochemical parameters obtained using the PDP curves for different temperatures are displayed in Table 5. A careful examination of the data illustrated in Table 5 reveals the difference in the corrosion potential between the MS uninhibited and inhibited with 10−3 M at a TGP concentration higher than 85 mV/SCE, for all the temperatures investigated. Further, it is evident that the cathodic slops within 10−3 M of TGP at different temperatures were higher compared to those of the acidic environment only; these data indicate that the TGP studied acted as an excellent anodic-type inhibitor. According to the data listed in Table 5, it is evident an increase in the temperature led to an increase in the corrosion current densities and to a decrease in the inhibitor efficiencies. This result confirms that the mild steel electrode can be dissolved by increasing the temperature. Additionally, these decreases could be explained by the phenomenon that sees the adsorption of the TGP on the surface of the metallic substrate.
In order to determine the activation energies of the corrosion process of the metallic electrode surface when it is plunged in an aggressive environment (0.5 M H2SO4) at 10−3 M of employed inhibitory resin, the Arrhenius formula between the icorr values and the temperature applied was used and investigated according to Equation (7). Additionally, an alternative formula for the Arrhenius relationship was investigated in order to calculate the activation standard enthalpy (∆Ha) and the activation standard entropy (∆Sa), according to Equation (8), as follows [36]:
i corr = Aexp E a R T
i corr = R T h N exp Δ S a T exp Δ H a R T
where A, h, N, and R denote the Arrhenius constant (pre-exponential factor), the Planck number, the Avogadro number and the perfect gas constant, respectively.
Furthermore, the activation energy values were measured using the slopes of the lines. In addition, Figure 6 illustrates the graphical representations of ln(icorr) as a function of 1000/T for MS uninhibited and inhibited with the investigated inhibitory resin at different concentrations. Table 5 illustrates the activation energy values calculated using the straight lines of the Arrhenius equation for MS uninhibited and inhibited with TGP. According to the results reported in several previously published studies, when the activation energy of the studied inhibitory resin is greater than that of the 0.5 M H2SO4, only ( E inh a > Ea) can be adsorbed on the electrode area of the metallic sample via electrostatic bonds (the physisorption mode). However, the activation energy of the inhibitory resin is lower than that of the 0.5 M H2SO4 alone ( E inh a < 𝐸a), which could be adsorbed on the metallic sample area by the covalent bonds (the chemisorption mode). A careful analysis of the data illustrated in Table 5 reveals that the activation energy value of the TGP is higher than that of the 0.5 M H2SO4 alone. These data indicate that the tested inhibitory resin could be adsorbed via the physico-chemical adsorption mode on the surface of the steel electrode. In addition, Figure 6 illustrates the variations in the ln(icorr/T) as a function of the 1000/T in an aggressive environment for MS uninhibited and inhibited with 10−3 M of the investigated inhibitory resin. In addition, the slopes of the straight lines were investigated in order to obtain the values of the activation standard enthalpy; meanwhile, the intercept allows us to determine the activation standard entropy. Furthermore, according to the results illustrated in Table 5, the positive value of ∆Ha reflects the endothermic nature of the surface of the steel electrode and the dissolution of the corrosion inhibition process. In addition, the negative signs of ∆Sa in the aggressive environment for uninhibited MS mean that the activated TGP molecule was in a lower order state compared to the initial step. However, in the presence of 10−3 M of TGP inhibitory resin, the activation standard entropy (∆Sa) is positive; this result indicates that the activated molecule was in a higher order state compared to the initial step. Furthermore, according to the data listed in Table 5, the activation standard entropy value with 10−3 M of TGP is higher than that of the uninhibited MS (0.5 M H2SO4 alone); thus, there is an increase in the disorder during the adsorption process.

3.5. DFT Results

In order to confirm the data obtained during the experimental part, the DFT descriptors, namely EHOMO, ELUMO, ∆Egap, χ, μ, S, ω, ΔN and ∆Eb−d, were calculated so as to study the inhibitory efficacy of the employed TGP inhibitory resin used to protect the mild steel electrode area from corrosion in the 0.5 M H2SO4 solution. For the sake of comparison, the neutral and protonated forms of the investigated resin were used in an aqueous environment. Figure 7 illustrates the optimized structures, the HOMO and LUMO energy molecular orbitals, electrostatic potential maps (ESP) and the energy gap profiles of the neutral epoxy resin (left) and its protonated form (right). In addition, the data obtained are displayed in Table 6.
The distribution of the HOMO surface of the neutral/protonated forms on all of the parts of the molecules, except the oxirane rings, implies that the carbonyl oxygen atoms and the oxygen ether group atoms are able to donate their electrons to the vacant d-orbital of the metal. However, the LUMO surface is mostly distributed on the carbonyl function and its neighboring carbon atoms, implying that this part of the epoxide resin has the ability to accept electrons from the occupied d-orbital of the metallic substrates area. These data are also supported by the electrostatic potential map.
As is known, EHOMO is the most commonly used quantum chemical descriptor, and is associated with the tendency to donate electrons to the vacant d-orbitals of the metallic area. The high values of the EHOMO suggest the high propensity of electrons to transfer to the lowest virtual d-orbital of the acceptor [37]. On the other hand, ELUMO indicates that inhibitors can accept electrons from the vacant d-orbital of the metallic area. It has been outlined that the lower that the ELUMO of the molecule is, the greater its tendency to accept electrons from the donor. Therefore, the adsorptive mechanism can be supported by the donation and acceptance competence to/from the inhibitor/metal and can also help to protect the metallic area from corrosion. The chemical reactivity of the inhibitors also depends on the energy gap between the EHOMO and ELUMO [38]. The smaller the value of the energy gap, the more reactive the inhibitory resin is. Additionally, small hardness values and higher softness values also have a critical role in increasing the anticorrosive properties of the inhibitory resins. The energy gap of the investigated compounds is small enough compared to the other compounds, which makes it reactive and able to protect the metallic substrates from corrosion. It was also illustrated that when η > 0 or ΔEb−d < 0, back-donation from the inhibitory resin to the metallic area is energetically favored [23]. The anticorrosive properties of the inhibitors are also affected by the higher values of the absolute electronegativity, which show that inhibitory resin cannot share electrons with others easily. The results in Table 6 show that the employed inhibitory resin has a comparable electronegativity value (4.264/4.228 eV (neutral/protonated), which enhances its inhibition efficiency. The positive values of the ΔN(0.108/0.113 (neutral/protonated)), which are lower than 3.6, suggest that the transfer of electron takes place between the inhibitory resin and the metallic surface [39]. The dipole moments of the synthesized resin in their neutral and protonated forms of 6.3 and 13.6 Debye, respectively, are higher than those of water of 1.8 Debye. This result displays the tendency of the inhibitory resin to expel water from the area of the metallic substrates and increase its inhibition efficiency. Furthermore, it was found that high values in the dipole moment shows a high value of polarizability, which is related to the inhibition efficacy of the inhibitory resin. As is known, if the resin easily polarizes, the tested resin can donate to or share electrons with the metallic substrates area.
Since we are dealing with the donation/acceptance of electrons to/from metal, it is helpful to estimate two important parameters that explain the metal–inhibitor interaction (ΔE1 and ΔE2) [38]. These descriptors are able to determine whether the electrons are transferred from the epoxide resin to the mild steel electrode area or vice versa. ΔE1 ( E L U M O i n h E H O M O   F e ) denotes the mechanism through which the electron is donated from the occupied d-orbitals of the Fe metallic samples area to the LUMO of the inhibitory resin. Meanwhile, ∆E2 ( E L U M O F e E H O M O i n h ) denotes the mechanism through which the electrons are donated from the HOMO of the inhibitory resin to the unoccupied d-orbitals of the Fe metallic substrates area [18]. In this study, the EHOMO and ELUMO of iron were found, using the 88th Edition of the CRC Handbook of Chemistry and Physics, to be −7.9024 and −0.151 eV, respectively [40]. For the neutral/protonated form, the data displayed in Table 6 illustrates that the ΔE1 (6.368/6.441 eV) values are slightly lower than those of ΔE2 (6.861/6.843 eV). These results reveal that the electron’s transfer from the metal ion to the inhibitor is energetically favored.
The 3D plots of the iso-area of the nucleophilic and electrophilic Fukui functions, and the dual Fukui descriptors of the employed inhibitory epoxide resin in its neutral and protonated forms in an aqueous environment with an iso value = 0.02 au and a density = 0.005 au, are displayed in Figure 8. A careful examination of the data illustrated in Figure 8 (in both forms) reveals that the electrophilic attack sites are more important compared to the nucleophilic attack sites, which might be explained by the presence of oxygen atoms. Figure 9 shows the condensed nucleophilic ( f k + ) and electrophilic ( f k ) Fukui functions, the local softness ( S k ± ), local electrophilicity ( ω k ± ) and their dual local descriptors ( f k 2 , Δsk and Δ w k ). All sets of numerical values are given in the supporting information (Tables S1 and S2) (note that the numbering system of the atoms is shown in Scheme 1). For both forms, the results reveal that the reactive nucleophilic attacks are produced by the carbonyl group (C1, O18, C2 and C19), whereas the more reactive sites for the electrophilic attacks are the oxygen heteroatoms (O21, O18 and O21 and O20) (Table 7).

3.6. Monte Carlo and Molecular Dynamic Approaches

Starting with the Fe(110) area makes it easy to compute the system’s adsorption energies. Applying Equation (9) to a molecule yields its adsorption energy, or Eads [41]:
E a d s o r p t i o n = E Fe 110 | | i n h i b i t o r   Fe 110 + E i n h i b i t o r  
where E Fe 110 | | i n h i b i t o r is the total energy of the simulated system, EFe, and Einhibitor are the total energy of the Fe(110) area and the corresponding free inhibitory resins, respectively.
After the MC approaches were successful (Figure 10), a considerable amount of research was conducted on the adsorption geometry of the inhibitory resin to make sure that the results were valid. This was performed to make sure that the data were correct. The MC simulation’s ability to reach a state of equilibrium can be judged by comparing the energy values at the steady state to the energy values at the start of the simulation. The simulation had reached the point at which the system could operate with the least amount of energy needed to keep going. The actual arrangement of the adsorbent inhibitory resin is displayed in Figure 10. This picture is based on a simulated plane of the type Fe-110. During the MD technique, the inhibitor decorates the area of Fe(110) in a way that is almost parallel to how the inhibitory resin is oriented on the side that maximizes contact to its oxygen atoms. We think that the backbone of an inhibitory resin that is willing to stick to the surface atoms of the Fe(110) plane causes this adsorption pattern, as shown in Figure 10. This theory is supported by the evidence that has been gathered. The tendency of resins to bring their heteroatoms (O) and electron rings closer to the area makes it possible for those inhibitory resins to adsorb, and the fact that they can do so is what gives those epoxide resins their ability to adsorb [42].
The adsorption of inhibitory resin onto the metallic substrates area offers data on the formation of the substantial Eads (Figure 11), which are visible on the metallic samples area. The inhibitory resin has an important adsorption interaction with the metallic area, owing to the extraordinarily high adsorption energies that it possesses. This activity, which is evidenced in the creation of a protective thin layer on the metallic substrates area, makes the electrode area of the metallic samples resistant to corrosion, which protects the metal’s original tate [26]. MD is widely considered to be a more specific illustration of the adsorption dynamics. In addition, after running the NVT simulation for several hundreds of picoseconds, it is abundantly obvious that the inhibitory resins depicted in Figure 11 adopt a configuration that is somewhat flat on the sides of the oxiran systems onto the metallic substrates area. Because of this, the inhibitory resins are powerfully absorbed onto the iron area. In addition, the radial distribution function (RDF) approach of the MD trajectory carried out during corrosion simulations is a simple technique that may be used to investigate the anticorrosion inhibitor that is adsorbed on the metallic substrates area (Figure 12). These results were obtained by running the simulations. This approach is uncomplicated and does not entail any tricky steps [43]. In MD simulations, the RDF (g) is a great way to assess the interaction between molecules and surfaces of iron [30]. The peak for chemisorption is between 1 and 3.5 Å, while the peak for physisorption is higher than 3.5 Å [28]. Figure 12 shows the radial distribution function of the inhibitor. This shows that the distance between the Fe (surface) and the inhibitor (2.81 and 2.95 Å for neutral and protonated inhibitor, respectively) is less than 3.5 Å. The negative energy value and the RDF peaks, indicates that the inhibitory resin has a significant effect on the metallic substrates area. A careful analysis of the data illustrated in Figure 12 reveals that the inhibitory resin mostly interacts with the iron area via oxygen atoms.

4. Conclusions

The corrosion inhibition of the employed inhibitory resin has been investigated via electrochemical analyses, morphology and computational approaches. PDP and EIS techniques illustrated that the inhibitory resin possessed high inhibitory efficacies at its optimum concentrations, such as at 88.6% (PDP) and 85.6% (EIS). SEM/EDS analyses revealed the potential ability of inhibitory resin to protect a mild steel surface in 0.5 M H2SO4. Data obtained using the PDP technique suggest that the employed epoxide resin can be adsorbed strongly in the metallic substrates area, as an anodic-type inhibitor. Additionally, the DFT calculation gives further data about the reactivity and energetic properties of the interaction between the inhibitory resin and the iron area. According to the data obtained using the MC and MD approaches, the inhibitory resin has a somewhat flat-adsorbed configuration on the side walls of the resin, which maximizes its contact with the oxygen heteroatoms. The concept that the inhibitory resin contains important negative adsorption energies is given further credence by the evidence that was gathered through experimental testing.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/polym15081967/s1, Figure S1: Langmuir adsorption isotherm model of TGP for MS in 0.5 M H 2SO4 at 298 K; Table S1: Hirshfeld charges, condensed Fukui functions and condensed dual descriptors for the neutral form of the TGP inhibitor molecules as obtained at the B3LYP/6-31G(d,p)/H2O/PCM model of theory. Units used below are “e” (elementary charge); Table S2: Hirshfeld charges, condensed Fukui functions and condensed dual descriptors for the protonated form of the TGP inhibitor molecules as obtained at the B3LYP/6-31G(d,p)/H2O/PCM model of theory. Units used below are “e” (elementary charge).

Author Contributions

Investigation, writing original draft preparation, supervision, review and editing, R.H.; formal analysis, K.D., A.E.M., A.H., N.D. and M.G.; Software, performance of experiment, writing, review and editing of the DFT part, Z.S. and N.W.; Software, performance of experiment, writing, review and editing of the MC and MD part, A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

Nuha Wazzan gratefully acknowledges King Abdulaziz University’s High-Performance Computing Centre (Aziz Supercomputer) for assisting with the calculations performed in this study. Avni Berisha gratefully acknowledges the support of the Ministry of Education, Science and Technology of Kosovo (Nr.2-5069) for providing the computing resources. Zaki Safi gratefully acknowledges the computational center at the Universidad Autonoma de Madrid, Spain, for assisting with the calculations performed in this study. Savaş Kaya from Cumhuriyet University, Faculty of Science, Department of Chemistry, Sivas, Turkey, is also gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lazrak, J.; Ech-chihbi, E.; Salim, R.; Saffaj, T.; Rais, Z.; Taleb, M. Insight into the corrosion inhibition mechanism and adsorption behavior of aldehyde derivatives for mild steel in 1.0M HCl and 0.5M H2SO4. Colloids Surf. A Physicochem. Eng. Asp. 2023, 664, 131148. [Google Scholar] [CrossRef]
  2. Loto, R.T.; Solomon, M.M. Application of ginger and grapefruit essential oil extracts on the corrosion inhibition of mild steel in dilute 0.5 M H2SO4 electrolyte. Sci. Afr. 2023, 19, e01489. [Google Scholar] [CrossRef]
  3. Ganjoo, R.; Sharma, S.; Verma, C.; Quraishi, M.A.; Kumar, A. Heteropolysaccharides in sustainable corrosion inhibition: 4E (Energy, Economy, Ecology, and Effectivity) dimensions. Int. J. Biol. Macromol. 2023, 235, 123571. [Google Scholar] [CrossRef] [PubMed]
  4. Devi, G.N.; Saranya, J.; Manjubaashini, N.; Thangadurai, T.D.; Roopan, S.M.; Chitra, S. Polyamidoaminoepichlorohydrin resin a novel synthetic anti-corrosive water soluble polymer for mild steel. Prog. Org. Coat. 2017, 109, 117–125. [Google Scholar] [CrossRef]
  5. Tamilselvi, B.; Bhuvaneshwari, D.S.; Padmavathy, S.; Raja, P.B. Corrosion inhibition of Pichia sp. biofilm against mild steel corrosion in 1 M H2SO4. J. Mol. Liq. 2022, 359, 119359. [Google Scholar] [CrossRef]
  6. Hsissou, R.; Lachhab, R.; El Magri, A.; Echihi, S.; Vanaei, H.R.; Galai, M.; Ebn Touhami, M.; Rafik, M. Synthesis Characterization and Highly Protective Efficiency of Tetraglycidyloxy Pentanal Epoxy Prepolymer as a Potential Corrosion Inhibitor for Mild Steel in 1 M HCl Medium. Polymers 2022, 14, 3100. [Google Scholar] [CrossRef]
  7. Hsissou, R.; Benassaoui, H.; Benhiba, F.; Hajjaji, N.; El Harfi, A. Application of a new trifunctional epoxy prepolymer, triglycidyl ethylene ether of bisphenol A, in the coating of E24 steel in 3.5% NaCl. J. Chem. Technol. Metall. 2017, 52, 431–438. [Google Scholar]
  8. Hsissou, R.; Benhiba, F.; El Aboubi, M.; Abbout, S.; Benzekri, Z.; Safi, Z.; Rafik, M.; Bahaj, H.; Kaba, M.; Galai, M.; et al. Synthesis and performance of two ecofriendly epoxy resins as a highly efficient corrosion inhibition for carbon steel in 1 M HCl solution: DFT, RDF, FFV and MD approaches. Chem. Phys. Lett. 2022, 806, 139995. [Google Scholar] [CrossRef]
  9. Dohare, P.; Quraishi, M.; Obot, I. A combined electrochemical and theoretical study of pyridine-based Schiff bases as novel corrosion inhibitors for mild steel in hydrochloric acid medium. J. Chem. Sci. 2018, 130, 8. [Google Scholar] [CrossRef]
  10. Rbaa, M.; Galai, M.; Ouakki, M.; Hsissou, R.; Berisha, A.; Kaya, S.; Berdimurodov, E.; Lakhrissi, B.; Zarrouk, A. Synthesis of new halogenated compounds based on 8-hydroxyquinoline derivatives for the inhibition of acid corrosion: Theoretical and experimental investigations. Mater. Today Commun. 2022, 33, 104654. [Google Scholar] [CrossRef]
  11. Hsissou, R.; Azogagh, M.; Benhiba, F.; Echihi, S.; Galai, M.; Shaim, A.; Bahaj, H.; Briche, S.; Kaya, S.; Serdaroğlu, G.; et al. Insight of development of two cured epoxy polymer composite coatings as highly protective efficiency for carbon steel in sodium chloride solution: DFT, RDF, FFV and MD approaches. J. Mol. Liq. 2022, 360, 119406. [Google Scholar] [CrossRef]
  12. Guo, L.; Safi, Z.S.; Kaya, S.; Shi, W.; Tüzün, B.; Altunay, N.; Kaya, C. Anticorrosive effects of some thiophene derivatives against the corrosion of iron: A computational study. Front. Chem. 2018, 6, 155. [Google Scholar] [CrossRef]
  13. Molhi, A.; Hsissou, R.; Damej, M.; Berisha, A.; Thaçi, V.; Belafhaili, A.; Benmessaoud, M.; Labjar, N.; El Hajjaji, S. Contribution to the corrosion inhibition of C38 steel in 1 M hydrochloric acid medium by a new epoxy resin PGEPPP. Int. J. Corros. Scale Inhib. 2021, 10, 399–418. [Google Scholar]
  14. Hsissou, R.; Dagdag, O.; Berradi, M.; El Bouchti, M.; Assouag, M.; El Bachiri, A.; Elharfi, A. Investigation of structure and rheological behavior of a new epoxy polymer pentaglycidyl ether pentabisphenol A of phosphorus and of its composite with natural phosphate. SN Appl. Sci. 2019, 1, 869. [Google Scholar] [CrossRef]
  15. Becke, A.D. A new mixing of Hartree–Fock and local density-functional theories. J. Chem. Phys. 1993, 98, 1372–1377. [Google Scholar] [CrossRef]
  16. Tomasi, J.; Mennucci, B.; Cammi, R. Quantum mechanical continuum solvation models. Chem. Rev. 2005, 105, 2999–3094. [Google Scholar] [CrossRef]
  17. Frisch, M.; Trucks, G.; Schlegel, H.; Scuseria, G.; Robb, M.; Cheeseman, J.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. Gaussian 09, revision A. 02, Gaussian, Inc., Wallingford. Phys. Rev. B Condens. Matter Mater. Phys. 1988, 37, 785. [Google Scholar]
  18. About, S.; Hsissou, R.; Erramli, H.; Chabebe, D.; Salim, R.; Kaya, S.; Hajjaji, N. Gravimetric, electrochemical and theoretical study, and surface analysis of novel epoxy resin as corrosion inhibitor of carbon steel in 0.5 M H2SO4 solution. J. Mol. Struct. 2021, 1245, 131014. [Google Scholar] [CrossRef]
  19. Özlem, U.; Gümüş, M.; Yusuf, S.; İrfan, K.; Atıf, K. Utilization of pyrazole-perimidine hybrids bearing different substituents as corrosion inhibitors for 304 stainless steel in acidic media. J. Mol. Struct. 2022, 1262, 133025. [Google Scholar]
  20. Bondarev, N.V.; Katin, K.P.; Merinov, V.B.; Kochaev, A.I.; Kaya, S.; Maslov, M.M. Probing of Neural Networks as a Bridge from Ab Initio Relevant Characteristics to Differential Scanning Calorimetry Measurements of High-Energy Compounds. Phys. Status Solidi (RRL)–Rapid Res. Lett. 2022, 16, 2100191. [Google Scholar] [CrossRef]
  21. Safi, Z.S.; Wazzan, N.; Aqel, H. Calculation of vertical and adiabatic ionization potentials for some benzaldehydes using hybrid DFT, multilevel G3B3 and MP2 methods. Chem. Phys. Lett. 2022, 791, 139349. [Google Scholar] [CrossRef]
  22. Molhi, A.; Hsissou, R.; Damej, M.; Berisha, A.; Bamaarouf, M.; Seydou, M.; Benmessaoud, M.; El Hajjaji, S. Performance of two epoxy compounds against corrosion of C38 steel in 1 M HCl: Electrochemical, thermodynamic and theoretical assessment. Int. J. Corros. Scale Inhib. 2021, 10, 812–837. [Google Scholar]
  23. Gómez, B.; Likhanova, N.V.; Domínguez-Aguilar, M.A.; Martínez-Palou, R.; Vela, A.; Gázquez, J.L. Quantum chemical study of the inhibitive properties of 2-pyridyl-azoles. J. Phys. Chem. B 2006, 110, 8928–8934. [Google Scholar] [CrossRef] [PubMed]
  24. Chandra, A.; Nguyen, M. Fukui Function and Local Softness; Taylor and Francis: New York, NY, USA, 2008; pp. 163–178. [Google Scholar]
  25. Erazua, E.; Adeleke, B. A Computational Study of Quinoline Derivatives as Corrosion Inhibitors for Mild Steel in Acidic Medium. J. Appl. Sci. Environ. Manag. 2019, 23, 1819–1824. [Google Scholar] [CrossRef]
  26. Mehmeti, V.V.; Berisha, A.R. 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. Front. Chem. 2017, 5, 61. [Google Scholar] [CrossRef]
  27. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef]
  28. Ouass, A.; Galai, M.; Ouakki, M.; Ech-Chihbi, E.; Kadiri, L.; Hsissou, R.; Essaadaoui, Y.; Berisha, A.; Cherkaoui, M.; Lebkiri, A. Poly (sodium acrylate) and Poly (acrylic acid sodium) as an eco-friendly corrosion inhibitor of mild steel in normal hydrochloric acid: Experimental, spectroscopic and theoretical approach. J. Appl. Electrochem. 2021, 51, 1009–1032. [Google Scholar] [CrossRef]
  29. Akkermans, R.L.; Spenley, N.A.; Robertson, S.H. COMPASS III: Automated fitting workflows and extension to ionic liquids. Mol. Simul. 2021, 47, 540–551. [Google Scholar] [CrossRef]
  30. Rahimi, A.; Farhadian, A.; Berisha, A.; Shaabani, A.; Varfolomeev, M.A.; Mehmeti, V.; Zhong, X.; Yousefzadeh, S.; Djimasbe, R. Novel sucrose derivative as a thermally stable inhibitor for mild steel corrosion in 15% HCl medium: An experimental and computational study. Chem. Eng. J. 2022, 446, 136938. [Google Scholar] [CrossRef]
  31. El-Aouni, N.; Hsissou, R.; Safi, Z.; Abbout, S.; Benhiba, F.; El Azzaoui, J.; Haldhar, R.; Wazzan, N.; Guo, L.; Erramli, H. Performance of two new epoxy resins as potential corrosion inhibitors for carbon steel in 1MHCl medium: Combining experimental and computational approaches. Colloids Surf. A Physicochem. Eng. Asp. 2021, 626, 127066. [Google Scholar] [CrossRef]
  32. Echihi, S.; Hsissou, R.; Benzbiria, N.; Afrokh, M.; Boudalia, M.; Bellaouchou, A.; Guenbour, A.; Azzi, M.; Tabyaoui, M. Performance of methanolic extract of artemisia herba alba as a potential green inhibitor on corrosion behavior of mild steel in hydrochloric acid solution. Biointerface Res. Appl. Chem 2021, 11, 14751–14763. [Google Scholar]
  33. Damej, M.; Hsissou, R.; Berisha, A.; Azgaou, K.; Sadiku, M.; Benmessaoud, M.; Labjar, N. New epoxy resin as a corrosion inhibitor for the protection of carbon steel C38 in 1M HCl. experimental and theoretical studies (DFT, MC, and MD). J. Mol. Struct. 2022, 1254, 132425. [Google Scholar] [CrossRef]
  34. Hsissou, R.; Benzidia, B.; Hajjaji, N.; Elharfi, A. Elaboration and electrochemical studies of the coating behavior of a new pentafunctional epoxy polymer: Pentaglycidyl ether pentabisphenol A phosphorus on E24 carbon Steel in 3.5% NaCl. J. Chem. Technol. Metall. 2018, 53, 898–905. [Google Scholar]
  35. Abd-El-Nabey, B.; Abdel-Gaber, A.; Elewady, G.; Sadeek, M.; Abd-El-Rhman, H. Inhibitive action of benzaldehyde thiosemicarbazones on the corrosion of mild steel in 3 M H3PO4. Int. J. Electrochem. Sci. 2012, 7, 11718–11733. [Google Scholar]
  36. Benhiba, F.; Sebbar, N.; Bourazmi, H.; Belghiti, M.; Hsissou, R.; Hökelek, T.; Bellaouchou, A.; Guenbour, A.; Warad, I.; Oudda, H. Corrosion inhibition performance of 4-(prop-2-ynyl)-[1, 4]-benzothiazin-3-one against mild steel in 1 M HCl solution: Experimental and theoretical studies. Int. J. Hydrog. Energy 2021, 46, 25800–25818. [Google Scholar] [CrossRef]
  37. Abboud, Y.; Abourriche, A.; Saffaj, T.; Berrada, M.; Charrouf, M.; Bennamara, A.; Al Himidi, N.; Hannache, H. 2, 3-Quinoxalinedione as a novel corrosion inhibitor for mild steel in 1 M HCl. Mater. Chem. Phys. 2007, 105, 1–5. [Google Scholar] [CrossRef]
  38. Kadiri, L.; Ouass, A.; Hsissou, R.; Safi, Z.; Wazzan, N.; Essaadaoui, Y.; Lebkiri, I.; El Khattabi, O.; Lebkiri, A. Adsorption properties of coriander seeds: Spectroscopic kinetic thermodynamic and computational approaches. J. Mol. Liq. 2021, 343, 116971. [Google Scholar] [CrossRef]
  39. Lukovits, I.; Kalman, E.; Zucchi, F. Corrosion inhibitors—Correlation between electronic structure and efficiency. Corrosion 2001, 57, 3–8. [Google Scholar] [CrossRef]
  40. Benhiba, F.; Hsissou, R.; Benzikri, Z.; Echihi, S.; El Blilk, J.; Boukhris, S.; Bellaouchou, A.; Guenbour, A.; Oudda, H.; Warad, I.; et al. DFT/electronic scale, MD simulation and evaluation of 6-methyl-2-(p-tolyl)-1,4-dihydroquinoxaline as a potential corrosion inhibition. J. Mol. Liq. 2021, 335, 116539. [Google Scholar] [CrossRef]
  41. Guo, L.; Zhang, S.; Li, W.; Hu, G.; Li, X. Experimental and computational studies of two antibacterial drugs as corrosion inhibitors for mild steel in acid media. Mater. Corros. 2014, 65, 935–942. [Google Scholar] [CrossRef]
  42. Berisha, A.; Podvorica, F.; Mehmeti, V.; Syla, F.; Vataj, D. Theoretical and experimental studies of the corrosion behavior of some thiazole derivatives toward mild steel in sulfuric acid media. Maced. J. Chem. Chem. Eng. 2015, 34, 287–294. [Google Scholar] [CrossRef]
  43. Yu, L.; Zhang, J.; Qiao, G.; Yan, Y.; Ti, Y.; Zhang, Y. Effect of alkyl chain length on inhibition performance of imidazoline derivatives investigated by molecular dynamics simulation. Mater. Corros. 2013, 64, 225–230. [Google Scholar] [CrossRef]
Scheme 1. Chemical structures of the TGP (neutral (left) and protonated (right)). The numbers do not follow the IUPAC system.
Scheme 1. Chemical structures of the TGP (neutral (left) and protonated (right)). The numbers do not follow the IUPAC system.
Polymers 15 01967 sch001
Figure 1. PDP plots for MS in an 0.5 M H2SO4 environment, uninhibited and inhibited with inhibitory resin at different concentrations.
Figure 1. PDP plots for MS in an 0.5 M H2SO4 environment, uninhibited and inhibited with inhibitory resin at different concentrations.
Polymers 15 01967 g001
Figure 2. Nyquist and Bode plots and EC for MS immersed in 0.5 M H2SO4, uninhibited and inhibited by inhibitory resin at different concentrations.
Figure 2. Nyquist and Bode plots and EC for MS immersed in 0.5 M H2SO4, uninhibited and inhibited by inhibitory resin at different concentrations.
Polymers 15 01967 g002
Figure 3. SEM images for uninhibited and inhibited MS substrates immersed in 0.5 M H2SO4 with 10−3 M of inhibitory resin.
Figure 3. SEM images for uninhibited and inhibited MS substrates immersed in 0.5 M H2SO4 with 10−3 M of inhibitory resin.
Polymers 15 01967 g003
Figure 4. EDS analysis for the uninhibited and inhibited MS electrode area after immersion in 0.5 M H2SO4 with 10−3 M of inhibitory resin.
Figure 4. EDS analysis for the uninhibited and inhibited MS electrode area after immersion in 0.5 M H2SO4 with 10−3 M of inhibitory resin.
Polymers 15 01967 g004
Figure 5. PDP plots of MS uninhibited and inhibited in 0.5 M H2SO4with inhibitory TGP resin (10−3 M) at different temperatures.
Figure 5. PDP plots of MS uninhibited and inhibited in 0.5 M H2SO4with inhibitory TGP resin (10−3 M) at different temperatures.
Polymers 15 01967 g005
Figure 6. Ln(icorr) as a function of 1000/T (left) and Ln(icorr/T) as a function of 1000/T (right) curves for MS uninhibited and inhibited in 0.5 M H2SO4 with 10−3 M of TGP.
Figure 6. Ln(icorr) as a function of 1000/T (left) and Ln(icorr/T) as a function of 1000/T (right) curves for MS uninhibited and inhibited in 0.5 M H2SO4 with 10−3 M of TGP.
Polymers 15 01967 g006
Figure 7. Optimized structures, 3D-plot electrostatic potential map and 3D-Iso-surface densities of HOMOs and LUMOs, and the energies and energy gaps (in eV) (Isovalue: 0.02 au) of the investigated epoxy resin in its neutral form (left) and protonated form (right) in an aqueous solution using a B3LYP/6-31+G(d,p) level of theory with a PCM solvation model.
Figure 7. Optimized structures, 3D-plot electrostatic potential map and 3D-Iso-surface densities of HOMOs and LUMOs, and the energies and energy gaps (in eV) (Isovalue: 0.02 au) of the investigated epoxy resin in its neutral form (left) and protonated form (right) in an aqueous solution using a B3LYP/6-31+G(d,p) level of theory with a PCM solvation model.
Polymers 15 01967 g007
Figure 8. The 3D plots of the iso-surfaces of the electrophilic ( f r (left), nucleophilic ( f + r ) (middle) and dual Fukui descriptor (right) of the TGP inhibitor, as found using the B3LYP/6-31+G(d,p)/H2O/PCM model of theory. Blue and negative lobes correspond to positive and negative values, respectively.
Figure 8. The 3D plots of the iso-surfaces of the electrophilic ( f r (left), nucleophilic ( f + r ) (middle) and dual Fukui descriptor (right) of the TGP inhibitor, as found using the B3LYP/6-31+G(d,p)/H2O/PCM model of theory. Blue and negative lobes correspond to positive and negative values, respectively.
Polymers 15 01967 g008
Figure 9. Condensed Fukui indices and the dual local descriptors (all values are in Hartree) of the TGP inhibitory resin. For the condensed Fukui functions: Open red and orange squares correspond to f k ,   filled blue circles correspond to f k + . For dual descriptors: Filled green squares correspond to f k 2 , filled red squares correspond to Δsk and blue x marks correspond to Δ w k .
Figure 9. Condensed Fukui indices and the dual local descriptors (all values are in Hartree) of the TGP inhibitory resin. For the condensed Fukui functions: Open red and orange squares correspond to f k ,   filled blue circles correspond to f k + . For dual descriptors: Filled green squares correspond to f k 2 , filled red squares correspond to Δsk and blue x marks correspond to Δ w k .
Polymers 15 01967 g009
Figure 10. (Left): MC and (Right): MD poses the lowest adsorption configurations for inhibitory resin adsorption in the simulated aggressive environment on the Fe(110) area.
Figure 10. (Left): MC and (Right): MD poses the lowest adsorption configurations for inhibitory resin adsorption in the simulated aggressive environment on the Fe(110) area.
Polymers 15 01967 g010
Figure 11. Distribution of the adsorption energy of the inhibitory resin employed in the simulated aggressive environment obtained using MC approaches.
Figure 11. Distribution of the adsorption energy of the inhibitory resin employed in the simulated aggressive environment obtained using MC approaches.
Polymers 15 01967 g011
Figure 12. RDF oxygen atoms for the inhibitor (neutral and protonated) on the Fe(110) area accomplished using MD trajectory approach.
Figure 12. RDF oxygen atoms for the inhibitor (neutral and protonated) on the Fe(110) area accomplished using MD trajectory approach.
Polymers 15 01967 g012
Table 1. PDP parameters for MS in 0.5 M H2SO4, with and without inhibitory resin at different concentrations.
Table 1. PDP parameters for MS in 0.5 M H2SO4, with and without inhibitory resin at different concentrations.
Conc.
(M)
Ecorr
(mV/Ag/AgCl)
icorr
(µA cm−2)
βc
(mV dec−1)
βa
(mV dec−1)
ƞPDP
(%)
Uninhibited0.5451185099121-
Inhibited with TGP10−65604721479074.4
10−55634691468674.6
10−45482221557988.0
10−35562111487188.6
Table 2. EIS parameters for MS immersed in 0.5 M H2SO4 medium, uninhibited and inhibited by different concentrations of TGP.
Table 2. EIS parameters for MS immersed in 0.5 M H2SO4 medium, uninhibited and inhibited by different concentrations of TGP.
Conc.
(M)
Re
(Ω cm2)
Rct
(Ω cm2)
RL
(Ω cm2)
L
(H)
Q
(µF sn−1)
ndlθƞEIS (%)
Uninhibited0.51.611.72.170.484300.830--
Inhibited with TGP10−62.239.06.450.601950.8110.7070.0
10−52.140.977.530.711750.8280.71471.4
10−42.278.2717.792.551480.7620.85085.0
10−32.181.3722.353.471400.7940.85685.6
Table 3. EDS analysis of MS surface after immersion in 0.5 M H2SO4 at each point with and without 10−3 M of TGP.
Table 3. EDS analysis of MS surface after immersion in 0.5 M H2SO4 at each point with and without 10−3 M of TGP.
Each PointC (%)O (%)S (%)Cl (%)Fe (%)
0.5 M H2SO4Spc_0013.4722.923.430.2469.93
Spc_0024.4734.057.710.2153.56
Spc_0033.6622.625.890.2467.60
10−3 M of TGPSpc_0015.1124.635.820.3064.14
Spc_0025.5917.662.280.3974.07
Spc_0035.049.510.970.2284.26
Table 4. Adsorption parameters of the TGP for MS in 0.5 M H2SO4 at 298 K.
Table 4. Adsorption parameters of the TGP for MS in 0.5 M H2SO4 at 298 K.
Kads (L/mol)R²ΔGads (KJ/mol)
TGP1,265,760.291−44.74
Table 5. Electrochemical and kinetic parameters of MS uninhibited and inhibited in 0.5M H2SO4 with inhibitory TGP resin (10−3 M) at different temperatures.
Table 5. Electrochemical and kinetic parameters of MS uninhibited and inhibited in 0.5M H2SO4 with inhibitory TGP resin (10−3 M) at different temperatures.
T
K
Ecorr
mVAg/AgCl
icorr
µA cm−2
βc
mV dec−1
η PC % Ea
kJ mol−1
ΔHa
kJ mol−1
−ΔSa
J mol−1 K−1
MS substrates/0.5 M H2SO4
298451185099-1512.5140.5
308453225092-
318449248096-
3284423340102-
TGP/MS substrates/0.5 M H2SO4
29855621114888.628.3225.8113.6
30854032513785.5
31854943513282.4
32854561415281.6
Table 6. Various DFT parameters for the epoxy resin (neutral and protonated forms).
Table 6. Various DFT parameters for the epoxy resin (neutral and protonated forms).
NeutralProtonated
Total energy, E(N) (Hartree)−1340.3180−1646.2206
Dipole moment, μ, (Dybe)6.413.6
EHOMO (eV)−7.012−6.994
ELUMO (eV)−1.517−1.461
ΔEgap (eV)5.4955.533
Iv (eV)6.9496.934
Av (eV)1.5171.465
χ (eV)4.2334.200
η (eV)2.7162.735
S (eV−1)0.3680.366
ω (eV)3.2993.225
ΔN0.1080.113
ΔEb−d (eV)−0.679−0.684
ΔE1 (eV)6.3866.441
ΔE2 (eV)6.8616.843
Table 7. Fukui functions and the dual local descriptors.
Table 7. Fukui functions and the dual local descriptors.
Atoms f k + f k f k 2 s k + s k Δs w k + w k + Δ w
C10.0420.0270.0150.4110.2630.1493.6742.3471.327
C100.2960.0490.2472.9080.4792.42925.9834.28121.703
O110.2490.1000.1492.4440.9791.46621.8438.74613.097
O130.0390.052−0.0130.3830.512−0.1303.4194.580−1.160
O140.0040.083−0.0800.0340.816−0.7820.3087.296−6.988
O210.0040.067−0.0630.0410.656−0.6150.3695.863−5.494
O250.0020.041−0.0390.0240.405−0.3820.2113.621−3.411
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hsissou, R.; Dahmani, K.; El Magri, A.; Hmada, A.; Safi, Z.; Dkhireche, N.; Galai, M.; Wazzan, N.; Berisha, A. A Combined Experimental and Computational (DFT, RDF, MC and MD) Investigation of Epoxy Resin as a Potential Corrosion Inhibitor for Mild Steel in a 0.5 M H2SO4 Environment. Polymers 2023, 15, 1967. https://doi.org/10.3390/polym15081967

AMA Style

Hsissou R, Dahmani K, El Magri A, Hmada A, Safi Z, Dkhireche N, Galai M, Wazzan N, Berisha A. A Combined Experimental and Computational (DFT, RDF, MC and MD) Investigation of Epoxy Resin as a Potential Corrosion Inhibitor for Mild Steel in a 0.5 M H2SO4 Environment. Polymers. 2023; 15(8):1967. https://doi.org/10.3390/polym15081967

Chicago/Turabian Style

Hsissou, Rachid, Khadija Dahmani, Anouar El Magri, Abdelfettah Hmada, Zaki Safi, Nadia Dkhireche, Mouhsine Galai, Nuha Wazzan, and Avni Berisha. 2023. "A Combined Experimental and Computational (DFT, RDF, MC and MD) Investigation of Epoxy Resin as a Potential Corrosion Inhibitor for Mild Steel in a 0.5 M H2SO4 Environment" Polymers 15, no. 8: 1967. https://doi.org/10.3390/polym15081967

APA Style

Hsissou, R., Dahmani, K., El Magri, A., Hmada, A., Safi, Z., Dkhireche, N., Galai, M., Wazzan, N., & Berisha, A. (2023). A Combined Experimental and Computational (DFT, RDF, MC and MD) Investigation of Epoxy Resin as a Potential Corrosion Inhibitor for Mild Steel in a 0.5 M H2SO4 Environment. Polymers, 15(8), 1967. https://doi.org/10.3390/polym15081967

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop