Abstract

The corrosion inhibition of mild carbon steels in an acidic environment by using ampicillin (AMP) has been evaluated based on experimental and quantum chemistry techniques. The experimental results indicate that the inhibition efficiency goes up at higher AMP concentration. The highest inhibition efficiency reaches 84.9% for polarization measurement and 90.1% for electrochemical impedance spectroscopy with the inhibitor of concentration 100 ppm at 298 K. The surface characteristics (SEM) also reconfirm the steel corrosion inhibition ability of AMP. Some important chemical factors such as EHOMO (highest occupied molecular orbital energy), ELUMO (lowest unoccupied molecular orbital energy), ΔELH (energy gap), and Fukui functions were calculated based on the optimized configuration of AMP at the theoretical level of B3LYP/6-31+G(d,p). Moreover, Monte Carlo and molecular dynamics simulations were used to analyze the absorption behavior of inhibitor on the surface of Fe(110), which supplies the mechanism of inhibition corrosion processes. The obtained results showed that AMP is considered to be a potential corrosion inhibitor for mild steel in 1M HCl medium. Moreover, the protonated state of AMP plays an important role in the protection of Fe surface against the corrosive process.

1. Introduction

Corrosion of metals is a key problem in industrial processes that has fascinated a large number of studies [14]. Because mild steel has good mechanical strength and low cost, it is widely used for industrial applications [5]. Besides, the acidic medium is usually utilized in industrial processes such as metal surface cleaning. This is the cause of the steel corrosion to be enhanced [68]. To prevent unexpected metal dissolution, inhibitors are used against corrosion [3]. Most inhibitors are selected from organic species containing π-systems and heteroatoms such as N, O, S, and P in the aromatic rings or having a long carbon chain [9, 10].

The advantage of many synthetic organic compounds is that they are cheap and possess efficient anticorrosive activity, but their disadvantage is that they are harmful to humans and pollute the environment. To remedy this situation, green corrosion inhibitors originated from natural products and drug compounds are preferred because of their eco-friendly and nontoxic nature [11].

Antibiotics are used for the treatment of human and animal diseases. However, the antibiotic residues may be released into the environment through wastewater effluent, medical waste, etc., which are the causes of environmental pollution. Fortunately, these compounds have heteroatoms in their structures, as mentioned in the previous paragraph, which can get effectively absorbed on the metal surface and thus protect it from corrosion [2, 12, 13]. In this context, one of the fields of green chemistry is renewable raw materials to protect the metallic surface for reducing the environmental impact and the wastes.

Practically, the obtained results show that ampicillin (AMP) can be used as the metal protecting compound in acidic media [8, 13, 14]. In the acidic solutions, AMP molecule can coexist in both neutral and protonated forms because it contains heteroatoms such as oxygen, nitrogen, and sulfur. In this regard, rarely reports have considered the protonated form of AMP and made clear the mechanism of corrosion protection in the acidic environment at the molecular level. The structure of AMP is shown in Figure 1 with the commercial name of (2S,5R,6R)-6-([(2R)-2-amino-2-phenylacetyl]amino)-3,3-dimethyl-7-oxo-4-thia-1-azabicyclo[3.2.0] heptane-2-carboxylic acid.

Normally, experimental techniques have been used to determine corrosion protection activity. Furthermore, computational methods based on density functional theory (DFT), molecular dynamics (MD) simulation, and Monte Carlo (MC) simulation have also been used widely to clarify the anticorrosive mechanism in recent years [15]. Thus, this technique is a helpful tool to obtain essential information such as molecular reactivity, charge transfer, preferred adsorption configuration, inhibitor/metal interface interaction, and the relationship between structural characteristics and anticorrosive efficiency.

In the current study, a combination of the experimental and theoretical approaches to exploit the anticorrosive mechanism of AMP was carried out and presented in detail. In the first part, electrochemical measurements and scanning electron microscopy were utilized to evaluate the corrosion inhibitor effect of AMP at the various concentrations in 1M HCl. In the second part, three computational methods (DFT, MC, and MD) were used to provide vital insights into the inhibition process at the molecular level.

2. Experimental Details

2.1. Materials

The major ingredients of the mild steel CT3 (Dana Italia) include (in percentage mass) C = 0.2682, Mn = 0.5798, Si = 0.2653, S = 0.0211, Ni = 0.0685, Cr = 0.0946, Cu = 0.2340, Mo = 0.0096, V = 0.0041, P = 0.0226, and balance Fe. The specimens with cylindrical shape and exposing 1 cm2 surface area were prepared for electrochemical experiments.

The aggressive solution was prepared by diluting 37% HCl, d = 1.19 kg/L (purchased from Merck), with distilled water. The investigated concentrations of AMP (Sigma-Aldrich) are 20, 40, 60, 80, and 100 ppm, respectively.

2.2. Electrochemical Measurements

Potentiostat PGS-HH5 equipment (manufactured by Vietnam Academy of Science and Technology) with a three-electrode system including the reference electrode (Ag/AgCl), the counter electrode (stainless steel), and the working electrode (mild steel CT3) was used for electrochemical studies.

After the specimen was immersed in the 1M HCl solution for 60 minutes, potentiodynamic polarization investigations were performed with a scanning rate of 1 mV/s in the potential range of −600 mV to −300 mV. All electrochemical experiments were set up at 298 K and stabilized by the Stanhope-Seta, Model 22101-4, with ±0.1°C accuracy. The inhibition efficiency (%IE) values were obtained by using the following equation:where and Icorr are the corrosion current density in the absence and presence of inhibitor, respectively.

The parameters of electrochemical impedance spectroscopy (EIS) were performed at open-circuit potential with an alternating current amplitude of 10 mV using a frequency region of 10 mHz–100 kHz.

The EIS was recorded with Zahner Zennium (IM6) using the Thales 4.5 software package. Anticorrosion behavior (η%) of EIS experiments was expressed by the following equation:where the charge transfer resistance corresponding to without/with AMP is and Rct, respectively [16, 17].

2.3. Surface Analysis

Mild steel samples were prepared and then immersed in 1M HCl with the AMP concentrations of 0, 20, 40, 60, 80, and 100 ppm in 6 hours. Thereafter, the specimens were washed with distilled water, degreased with acetone, wiped with a clean towel paper, and finally air-dried. The surfaces of specimens were observed by using an EVO-ZEISS scanning electron microscope (SEM) at the Analysis Centre for Import Export Products, Danang Branch (Vietnam).

2.4. Computational Chemistry Approaches
2.4.1. Quantum Chemical Calculation

The B3LYP/6-31+G(d,p) method was utilized to calculate the corrosion descriptors for both neutral and protonated forms of AMP. The solvent effect on the enthalpies of this system was estimated by the IEF-PCM solvation method [18, 19]. Gaussian 09 software was used for all calculations [20].

Based on the above, the reactions are described as follows:

To determine the most favorable protonation site, protonation affinity (PA) and basicity (B) parameters were calculated via the following formula [21]:where H and are the enthalpies and Gibbs free energy of neutral form (AMP), protonated form (AMPH+), and proton (H+). H and of the proton at 298 K were taken from the literature [21, 22].

Adsorption centers of AMP are analyzed via the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). Energy gap (∆ELH) is estimated based on the following equation:

The local reactivity for nucleophilic (f+) and electrophilic (f) attacks was measured by Fukui functions [23]. These indexes were visualized by Multiwfn software [2426].

2.4.2. Monte Carlo and Molecular Dynamics Simulations

The simulation methods were performed by Material Studio 7.0. The Fe(110) surface with a thickness of 5 Å was relaxed by using the smart minimizer method, and then it was enlarged to a (10 × 10) supercell. First, MC simulations have been performed in the adsorption locator tool using COMPASS with an energy convergence tolerance of 2 × 10−5 kcal/mol and force convergence tolerance of 0.001 kcal/mol/Å. The proportion of water molecules to hydrogen chloride corresponding to 1M HCl solution was 500/9 [27]. The system (24.82 × 24.82 × 26.09 Å3) containing 491 H2O, 9 H3O+, 9 Cl, and 1 inhibitor molecule was constructed using the amorphous cell program. Finally, the corrosion system was built by placing the amorphous cell on Fe(110) supercell using the layer builder tool. Afterwards, MD simulations were executed in a simulation box (24.82 × 24.82 × 37.71 Å3). The COMPASS force field was applied for the simulation of all AMP species and the whole system. The parameters set up for the MD simulation such as the temperature of 298 K, the time step of the canonical ensemble (NVT) of 1.0 fs, and simulation time of 500 ps. The interaction energy (Einteraction) and binding energies (Ebinding) between the studied molecule and the Fe(110) surface were computed using the following equations:where Etotal, Esurface+solution, and Einhibitor denote the overall energy of the entire system, the energy of Fe(110) surface and solution without the inhibitor, and energy of inhibitor, respectively.

3. Results and Discussion

3.1. Experimental Results
3.1.1. Polarization Measurements

The results obtained from polarization experiments for mild steel with AMP’s various concentrations in HCl media are shown in Figure 2. It is clear that the polarization plots have tended to shift to lower corrosion current density regions when AMP concentration changes from 0 to 100 ppm. Depending on the Ecorr shift, there are three common types of inhibitors such as cathodic, anodic, and mixed type. According to the literature reviews [28, 29], the inhibitor belongs to a mixed class if Ecorr (inhibitor) is lower than 85 mV compared with Ecorr (blank). As can be seen in Figure 1, the shift of Ecorr (AMP) displacement is not remarkable; thus, AMP is regarded as a type of mixed inhibitor.

Electrochemical corrosion factors recorded from the polarization measurements include corrosion potential (Ecorr), corrosion current density (Icorr), and cathodic and anodic Tafel slopes (βc and βa) (displayed in Table 1). Furthermore, it can be observed from Table 1 that the values of βa and βc change significantly in the presence of the inhibitor. It indicates that AMP inhibits the anodic and cathodic reactions. It means that AMP interacts with the iron surface to form the film that can prevent the corrosion process occurring at the anodic electrode (mild steel) and make slowly the H+ ion reduction at the cathodic electrode.

The results of Figure 2 and Table 1 show that the inhibition efficiency (%IE) of AMP goes up when increasing concentration and reaches maximum behavior of 84.9% at 100 ppm.

3.1.2. Electrochemical Impedance Spectroscopy (EIS)

EIS is a helpful tool to provide more information about the anticorrosion activity of studied compounds on steel surfaces in an acidic environment. EIS of specimens with the absence/presence of AMP is depicted in Nyquist diagrams (Figure 3). These diagrams display that when increasing the AMP concentrations, the Nyquist diameters significantly increase as compared with the blank sample and there are nearly no changes in the semicircular shapes. It means that the anticorrosion behavior of AMP on the iron surface is remarkable in an acidic environment.

As can be seen in Figure 3, all curves show a typical Nyquist diagram with single semicircles shifted along the real impedance of the horizontal axis. This indicates that double-layer behavior and the charge transfer of the corrosion process primarily control the corrosion of low carbon steel.

For better visualization of the inflection points of the impedance module as well as the phase angle variation, Bode and phase angle plots are also analyzed and depicted in Figure 4. The value of the phase angle increases significantly as the inhibitor concentration increases. It is worth noting that the phase angle values of the blank and the inhibitor samples are always lower than −90°, which proves that the capacitor is nonideal.

To analyze the experimental results, an appropriate equivalent circuit model is required to correctly fit the impedance curves. The good fitted data are generally explained using the equivalent circuit shown in Figure 5, where Rs is the solution resistance, Rct is the charge transfer resistance, CPE is the constant phase element, L is the inductance, and RL is the inductive resistance. All relevant impedance parameters are calculated, and the results are given in Table 2. In order to observe better fit, the CPE element has been substituted for the capacitor. The impedance of the CPE can be given by the following equation [3032]:where Yo is the CPE constant, j is the imaginary unit, ω is the angular frequency (given by ω = 2πf), and n is the CPE exponent. If the electrode surface is homogeneous and plane, n is equal to 1 and the electrode surface can be treated as an ideal capacitance. The double-layer capacitance (Cdl) can be simulated via CPE using the following equation [31]:

The data obtained from Table 2 show that the resistance significantly changes with the increase of AMP concentration. Specifically, Rct value of the blank sample is 11.91 Ω·cm2 in HCl. In the presence of AMP with concentrations from 20 ppm to 100 ppm, Rct values change from 77.54 to 129.70 Ω·cm2. This loses control of electrical transfer, thereby reducing corrosion in hydrochloric acid. The possibility of anticorrosion increases gradually according to the AMP concentration.

When increasing AMP concentration, the Cdl value tends to decrease. Conforming to the Helmholtz model, the change of the Cdl value can make a valuable substitution of water molecules by the AMP molecules onto the iron surface, reducing the local dielectric constant or enhancing the double electrical layer thickness [32, 33]. The decrease of Cdl may be caused by the adsorption of AMP taking place on the surface of the steel sample [34].

In addition, when AMP is added into acidic solution, the n-value decreases, which indicates that the electrode surface becomes less homogeneous, relating to the existence of an adsorbent layer formed on the steel surface [35].

Besides, the corrosion inhibition behavior increased from 83.4% to 90.1% when the AMP concentration increased from 20 ppm to 100 ppm. It reveals that an adsorbent layer may form on the steel surface to protect the electrode corrosion [35]. Based on data in Figure 3 and Table 2, the results of the impedance spectroscopy are in the same trend with data of the polarization curves.

3.1.3. Surface Analysis

Surface analysis can be interpreted based on the iron surface images in an acidic environment without/with the presence of the AMP inhibitor. Figure 6 displays SEM images of the specimen with AMP’s optimal concentration before and after treating in 1M HCl solution. The surface of the specimen is very rough with several pits and cavities in 1M HCl because of acid attack (see Figure 6(b)). On the contrary, with the presence of AMP, the micrograph surface of iron becomes smoother and more homogeneous (see Figure 6(c)). This can be explained by an inhibiting film layer was formed by the AMP molecule adsorption on the iron surface, preventing the dissolution of mild steel electrode. This result also affirms the inhibition ability of AMP on mild steel corrosion.

3.2. Computational Study
3.2.1. Quantum Chemical Calculation

The structural characteristics of the inhibitors were studied by using computational chemistry approaches to further explain the pathway of inhibition. In acidic environment, heteroatoms (such as O, N, and S) in the corrosion inhibitor molecule can be easily protonated [21, 36]. Hence, to mimic the experimental conditions, both existed states of AMP were investigated. AMP has seven possible sites for protonation which are O19, O20, O24, N9, N12, N18, and S14 (Cartesian coordinates of all optimized structures are shown in Table S1). To find the most likely protonation position of AMP, proton affinity and basicity values were evaluated by following the B3LYP/6-31+G(d,p) method, and the obtained results are shown in Table 3. It is found show that the priorly probable position for protonation is N18 with PA = 227.4 kcal/mol and GB = 235.2 kcal/mol. For further explanation, optimized geometry structures, frontier molecular orbitals, and essential quantum chemical factors for both existed forms of AMP are discussed in detail.

The effect of the corrosion protection process depends on the adsorptive activity of an inhibitor on the metal surface, which was explained by the electron donor-acceptor ability between the organic compound and iron atom. HOMO electron distribution characterizes the electron donor ability, while LUMO electron distribution represents the acceptor ability. The optimized structures, HOMOs, and LUMOs of AMP inhibitor following neutral and protonated states are depicted in Figure 7.

Regarding the neutral form, the HOMO finds a major contribution from the benzene ring and the amino (-NH2) group, while the LUMO contributed around the ring containing sulfur atom and carboxylic (-COOH) group. For the protonated form, the electron density region allocated around two heterocycles containing nitrogen and sulfur atoms is denoted by HOMO, while the electron density of LUMO shows the main contribution from the aromatic ring and the amino group. It can be seen that the delocalized electron density of frontier molecular orbital of protonated AMP is reversed to that of neutral ones. Accordingly, these areas are almost active in the transfer of electrons between the metal surface and APM or vice versa. The frontier molecular orbital density distributions were useful in adsorptive orientation between the metal surface and AMP [37]. Moreover, the dominant solvent models [3840] for both neutral and cationic forms are also shown in Figure S1.

To indicate the strength of donor-acceptor interactions between the surface of the iron and APM, the energies of EHOMO, ELUMO, and ΔELH (energy gap) were calculated. These main parameters define the reactivity of an inhibitor molecule. If the EHOMO value is higher, the electron donating ability will be easier, whereas if the ELUMO value is lower, the electron accepting ability is more preferred. The obtained results in Figure 7 indicate that the EHOMO values of the neutral state (−6.636 eV and −6.794 eV) are higher than those of protonated state (−9.497 eV and −7.001 eV) in the gas phase and solvent, respectively. Hence, the neutral form tends to transfer electrons to the vacant 3d-orbitals of iron and forms an adsorptive layer on the surface of mild steel. However, the ELUMO values of the protonated types in the gas phase and water (−4.705 eV and −1.356 eV, respectively) are lower than those of the neutral types (−1.042 eV and −1.016 eV, respectively). It means that the protonated state is easy to gain the electrons from the d-orbital of the metal compared with the neutral state. In addition, the energy gap (ΔELH) also informs about the reactivity or stability of the inhibitor molecule. It is an essential descriptor that characterizes the binding ability of inhibitor molecules to the metal surface. With a lower ΔELH value, a molecule shows greater reactivity, leading to an increase of adsorption onto an iron surface [41, 42]. Figure 8 visualizes the density of states (DOS) spectra in the gas phase and water for both existed forms of AMP. The green lines are the occupied orbitals. The red lines are the unoccupied orbitals. Also, the gap between them is the band gap. Based on data from Figure 8, the ΔELH value of the protonated species is 0.802 eV and 0.132 eV lower than that of the neutral forms corresponding to the gas phase and water. Hence, the protonated forms of AMP are more expected to adsorb on the metal surface.

The nucleophilic/electrophilic sites (susceptible atoms and regions) in the molecular structures of AMP were localized by the Fukui function and . Figure 9 displays the electron density surfaces of and for two states of AMP. The areas of the molecules which have higher values of are more sensitive to be attacked by electron abundant groups. On the contrary, the positions of the molecules which have higher values of are interacted easily with electron poor groups. As can be shown in Figure 9, the site of the sulfur atom, the amino group, and the carboxyl group is most likely as electrophilic centers, which cannot interact with negatively charged steel, while the function distributes almost over phenyl ring that donates electrons to surface Fe atoms. Otherwise, the electron density surfaces for and of the protonated AMP are different from those of neutral species.

3.2.2. Molecular Simulations

MC and MD simulations were considered as useful theoretical tools to analyze the interaction between adsorbent and Fe(110) surface at a molecular level. The lowest energy configurations of neutral (AMP) and protonated (AMPH+) forms on Fe(110) obtained from the MC simulation are presented in Figure 10. It is clear that the molecular skeleton of neutral/protonated AMP is nearly parallel to the plane of the metal surface. The adsorption energy is more negative, and the studied system is more stable with strong interaction between the inhibitor and iron surface. As can be seen from Table 4, the obtained results indicate that the protonated form has a greater ability to interact with the metallic surface and adsorbs more strongly compared to the neutral form because the adsorption energy of AMPH+ (−196.15 kcal/mol) is more negative than that of AMP (−181.24 kcal/mol).

To imitate the real corrosive environment, MD simulations were carried out in a box which consisted of 491 H2O, 9 H3O+, and 9 Cl corresponding to 1M HCl solution to further examine the anticorrosion activity of AMP. From the most stable adsorption configurations of both neutral and protonated AMP on Fe(110)/HCl system (Figure 11), it can be seen that all the visualized configurations have a tendency to orient parallel to the iron surface in HCl solution. This is similar to the obtained results from MC simulations. It proves that the heteroatoms and π-electrons in the states of AMP (neutral and protonated) can interact strongly with vacant d-orbitals of Fe to form the stable coordination bonds. Then, the thin film layer can be formed covering the steel surface to protect the corrosion process. Furthermore, the interactions between the inhibitor and the steel surface were quantitatively studied by the interaction energy (EInteraction) and binding energy (EBinding), which are the key components to determine the binding ability of AMP onto the iron surface. According to data from Table 5, the EInteraction value of Fe(110)/AMPH+/HCl system is more negative than that of Fe(110)/AMP/HCl system. It indicates that the adsorption in the Fe(110)/AMPH+/HCl system is more stable and stronger leading to the higher corrosion inhibitor efficiency than Fe(110)/AMP/HCl system. Thus, the protonated form of AMP has significantly contributed to the protection of Fe surface against the corrosive process.

3.3. Proposed Mechanism of Adsorption and Inhibition

Based on the data obtained from the experimental and theoretical investigations, the suggested mechanism for the adsorption and inhibition of neutral/protonated AMP in 1M HCl solution is displayed in Figure 12. As per Figure 12, the adsorption process of AMP on the iron surface can occur by the electrostatic interaction, the interaction of lone electron pairs, π-electron, and vacant d-orbital of surface iron atoms. In the acid environment, AMP is protonated to become cationic species and the electrostatic interaction may occur with chloride ions (Cl) which gets absorbed in the positively charged iron surface (Fe2+). It can lead to forming a film layer of protective complex (FeClAMPH+). Moreover, the neutral form of AMP may exist in acidic solutions. Expectedly, the neutral AMP may get adsorbed on the mild steel surface by sharing electrons between heteroatoms (O, N, and S), π-electron, and vacant d-orbital of iron atoms.

4. Conclusions

In this context, the anticorrosion mechanism of AMP on mild steel in HCl 1M acid medium was explored using both experimental and theoretical approaches. Some results obtained from this paper can be highlighted as follows:(i)Experimental techniques (including electrochemical measurements and electrochemical impedance spectroscopy) demonstrate that AMP has potential inhibitive characteristics for the corrosion of mild steel in 1M HCl acid. The corrosion inhibition behavior of AMP increases following the concentration range of 0 ppm–100 ppm and reaches maximum efficiency at 100 ppm.(ii)SEM images also indicated that the iron surface was covered with the adsorbed protective molecules of AMP.(iii)To further explain the anticorrosion properties in the real environment, both neutral/protonated states (AMP/AMPH+) were studied in detail through quantum chemical calculations, Monte Carlo simulation, and molecular dynamics simulations. It can be noted that AMPH+ has significantly contributed to protecting the mild steel surfaces against the corrosive process.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by the Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 104.06-2020.17 (P.C.N).

Supplementary Materials

Table S1: optimized structure of neutral and protonated ampicillin in gas phase at the B3LYP/6-31+G(d,p) level. Figure S1: optimized structure, highest occupied molecular orbital (HOMO), and lowest unoccupied molecular orbital (LUMO) for both neutral (AMP) and protonated (AMPH+) forms under dominant solvent model. (Supplementary Materials)