Fang Luoa,
Pei Yaob,
Jiale Zhanga and
Entian Li*a
aSchool of Petroleum Engineering, Changzhou University, Wujin District, Changzhou 213164, Jiangsu, PR China. E-mail: let@cczu.edu.cn
bDepartment of Chemistry and Materials Engineering, Changzhou Vocational Institute of Engineering, Changzhou 213164, Jiangsu, PR China
First published on 1st August 2025
Using electrochemical measurements, the corrosion inhibition performance of pyridine bromide ionic liquid corrosion inhibitors on low carbon steel was evaluated under different concentrations and alkyl chain lengths of anions. The optimal process conditions were obtained: at 25 °C, in a 1 M HCl corrosion medium, the 1-dodecyl-3-methylpyridine bromide ionic liquid ([C12Py]Br) exhibited the best corrosion inhibition effect, with maximum inhibition rates of 94.1%. The study mechanism showed that the adsorption process was a mixed adsorption dominated by chemical adsorption. Surface analysis (SEM/EDS, XPS) and quantum chemistry results verified that [CnPy]Br could form a protective film on the surface of low carbon steel and inhibit its corrosion. In addition, molecular dynamics simulations were used to analyze the adsorption behavior of corrosion inhibitor molecules on the Fe(110) surface, and the formation sequence of bond and non-bond interactions in the Fe corrosion inhibitor molecular system on the Fe(110) surface was obtained, which was consistent with the experimentally determined corrosion inhibition efficiency. The radial distribution function (RDF) showed that the adsorption mode of the corrosion inhibitor on the metal surface was a mixed adsorption dominated by chemical adsorption and supplemented by physical adsorption.
The use of corrosion inhibitors in a corrosive environment effectively prevents severe damage caused by corrosion on metal surfaces. Corrosion inhibitors exhibit water solubility and maintain thermal stability even in highly aggressive acidic environments. These compounds can interact with metallic ions to create a metal-inhibitor complex that remains chemically, thermally, and thermodynamically stable, particularly in acidic environments.4 Next, the resulting complex is efficiently adsorbed onto the metal surface, creating a stable thin film that serves as a protective barrier, effectively isolating the metal from the corrosive environment. As a result, the protective film provides significant defense for the metallic surface, shielding it from corrosion processes effectively. The performance of corrosion inhibitors is attributed to their heteroatoms, such as nitrogen, oxygen, phosphorus, and sulfur, as well as their polar functional groups (–CN–, –NH2, –OH, –OCH3, –SH).
In the processes of adsorption and inhibition, the π-electrons of aromatic rings and lone p-electrons from heteroatoms are transferred into the empty d-orbitals of iron, forming strong covalent bonds.5 However, numerous organic compounds containing heteroatoms like phosphorus, oxygen, nitrogen, and sulfur have been identified as potential alternatives to toxic corrosion inhibitors. Consequently, a key area in contemporary corrosion inhibitor research is the development of innovative and eco-friendly solutions.6,7
A considerable number of researchers focus on utilizing environmentally friendly corrosion inhibitors, including natural extracts, pharmaceutical compounds, ionic liquids (ILs), and similar substances. In recent years, ionic liquids (ILs) have garnered significant interest owing to their distinctive characteristics, including minimal toxicity, high polarity, non-flammability, low vapor pressure, remarkable solubility, and exceptional thermal stability.1 In addition, their adjustable properties enable the creation of various “custom-designed” functional ionic liquids (ILs) tailored for specific applications and materials. For instance, the anti-corrosion properties of ionic liquids (ILs) can be tailored at the molecular level to safeguard materials with specific compositions, thereby enhancing the efficiency of inhibition.8 Among these compounds, imidazolium-based ionic liquids (ILs) with carbon chains of varying lengths have been suggested as effective candidates for corrosion inhibition. This is attributed to the aromatic ring, which serves as the bonding site between the metallic surface and the ILs.9,10 Conversely, the carbon chain tail contributes to the formation of a hydrophobic layer, which aids in protecting metallic surfaces by providing adequate coverage. For example, S. Shuncun et al.11 investigated how the length of the alkyl chain attached to the imidazolium ring in novel ionic liquids (ILs) influenced their corrosion inhibition performance in acidic solutions. They found that the inhibition efficiency improved as the carbon chain length of the alkyl group increased. They concluded that as the carbon chain length of the alkyl group increased, the inhibition efficiency also improved. In the study by Azeez et al.,12 three environmentally friendly ionic liquids (ILs)—1-methyl-3-propylimidazolium iodide (MPIMI), 1-butyl-3-methylimidazolium iodide (BMIMI), and 1-hexyl-3-methylimidazolium iodide (HMIMI)—were evaluated as inhibitors to reduce mild steel corrosion in a 1 M HCl solution.
Moreover, since pyridine-based ILs also have an aromatic ring, the corrosion inhibition effect of pyridine-based ILs corrosion inhibitors mild steel has been increasingly reported. In the study by Bousskri et al.,13 the protection of mild steel against corrosion in 1 M HCl was examined in the presence and absence of different concentrations of 1-(2-(4-chlorophenyl)-2-oxoethyl) pyridazinium bromide (CPEPB). The findings indicate that CPEPB is an effective inhibitor, with the inhibition efficiency reaching 91% at a concentration of 10−3 M. Despite the fact that existing studies have revealed the possible mechanisms of corrosion inhibition through various technical means, there is still a lack of in-depth analysis of the interaction between pyridine-based ILs corrosion inhibitors and metal surfaces, especially the impact of the molecular structure of the corrosion inhibitor on its performance.14,15
Due to the structural tunability and validated physicochemical stability of [CnPy] Br ionic liquids, we chose them. Pyridine based ionic liquids, such as [CnPy] Br, have excellent corrosion inhibition potential due to the strong adsorption promoted by π electrons through aromatic pyridine rings, and the hydrophobicity and surface coverage regulated by alkyl chains. Systematically changing the alkyl chain length (n = 4, 8, 12) can clearly understand the relationship between molecular structure and inhibitory performance. This design draws on the best case of pyridine based ionic liquids in the field of corrosion, such as Alrefaee et al.'s16 report that 1-dodecylpyridin-1-ium bromide (Pyr-C12-Br) has a corrosion inhibition rate of over 91% on low carbon steel in 1 M HCl. Prior studies confirm pyridinium ILs achieve >90% inhibition in acidic media (e.g., Bousskri et al. reported 91% for CPEPB in 1 M HCl17), yet none systematically investigated alkyl chain length effects on mild steel.
The synthesis of [CnPy] Br adopts a single step quaternization reaction (3-methylpyridine + bromoalkane) with a yield of >90%. This route conforms to the principles of green chemistry and avoids complex purification steps. Given its high efficiency, tunable structure, and low toxicity, [CnPy]Br holds promise for applications in industrial cooling water systems, oilfield acidizing processes, and metal processing, addressing the urgent need for sustainable corrosion inhibitors. Internationally, pyridinium ILs show promise in: oil/gas pipelines: Saudi Aramco field tests reduced corrosion rates by 94% in acidizing fluids;18 metal processing: Chinese steel mills adopted similar ILs for pickling baths, cutting waste treatment costs by 40%.19
In this study, mass loss tests and surface analysis techniques (scanning electron microscopy (SEM) and X-ray photoelectron spectroscopy (XPS)) were employed to assess the corrosion inhibition performance of [CnPy]Br ionic liquids on mild steel in 1 M HCl. The corrosion resistance properties and mechanisms of [CnPy]Br ionic liquids were further investigated using potentiodynamic polarization measurements (PDP) and electrochemical impedance spectroscopy (EIS). Lastly, the proposed corrosion inhibition mechanism was confirmed by adsorption analysis (adsorption isotherm model and quantum chemical calculations) and molecular dynamics simulations.
Chemical reagents | Providers | Purity |
---|---|---|
3-Methylpyridine | Macklin | ≥99% |
1-Bromobutane | Macklin | ≥99% |
1-Bromooctane | Macklin | ≥98% |
1-Bromododecane | Macklin | ≥98% |
Ethyl acetate | Sinopharm Chemical Reagents Co., Ltd | ≥99.5% |
Potassium bromide | Sinopharm Chemical Reagents Co., Ltd | ≥99% |
Acetone | Sinopharm Chemical Reagents Co., Ltd | ≥99.5% |
Hydrochloric acid | Jiangsu Yongfeng Chemical Reagent Co., Ltd | 37% |
Anhydrous ethanol | Jiangsu Yongfeng Chemical Reagent Co., Ltd | ≥99.5% |
The characterization of the synthesized products was conducted with a Thermo Scientific Nicolet iS50 FT-IR, as illustrated in Fig. 2. The peak observed at 3419 cm−1 is a water peak caused by moisture absorption. The peak at 3032 cm−1 is the result of the C–H stretching vibration of the pyridine aromatic ring. The peak at 2928 cm−1 is caused by C–H stretching vibrations of methyl groups.22,23 The peaks at 1633 cm−1 and 1492 cm−1 can be attributed to the stretching vibrations of CN and C
C in the pyridine aromatic ring, respectively. Asymmetric stretching vibration of methyl and methylene groups leads to the appearance of peak 1492 cm−1 (ref. 24). A peak observed at 1143 cm−1 indicates the in-plane deformation vibration of the C–H bond on the pyridine ring; an absorption peak around 1358 cm−1 confirms the presence of carbon–nitrogen single bonds. The increasing intensity of the hydroxyl peak at 3419 cm−1 with longer carbon chains is primarily attributed to the enhanced intermolecular hydrogen bonding. As the alkyl chain length increases, the hydrophobic effect promotes the aggregation of pyridinium cations, leading to a higher local concentration of hydroxyl groups from absorbed water molecules. This aggregation facilitates stronger hydrogen bonding networks, resulting in a more pronounced hydroxyl stretching vibration peak. Additionally, the increased molecular weight and surface area of longer-chain compounds may enhance water adsorption capacity, further contributing to the observed peak intensity increase. The presence of these spectral bands offer initial confirmation that the synthesis of compound [CnPy]Br has been successfully achieved.
C | Si | Mn | P | S | Fe |
---|---|---|---|---|---|
0.18 | 0.30 | 0.32 | 0.04 | 0.03 | Balance |
A 1 M HCl solution was made by mixing concentrated hydrochloric acid (37%) with purified water. Corrosion inhibitors were introduced into the 1 M HCl solution to create test solutions with inhibitor concentrations ranging from 0.25 mmol L−1 to 8 mmol L−1.
![]() | (1) |
![]() | (2) |
ηw = θ × 100% | (3) |
The EIS was conducted with an AC voltage amplitude of 10 mV and a frequency range of 100 kHz to 0.01 Hz. Eqn (4) was used to calculate the inhibition efficiency (ηE) using the charge transfer resistance found, Rct.
![]() | (4) |
PDP investigations were carried out by scanning the potential within a range of ±250 mV relative to the open circuit potential (OCP) at a scan rate of 0.5 mV s−1. The obtained corrosion current density Icorr was used to estimate the inhibition efficiency (ηP) by the formula (5).
![]() | (5) |
This study analyzed the optimized molecular structure and evaluated several theoretical parameters. These parameters encompassed the highest occupied molecular orbital energy (EHOMO), the lowest unoccupied molecular orbital energy (ELUMO), the energy gap (ΔEGAP), ionization potential (I), electron affinity (A), electronegativity (χ), hardness (η), softness (σ), and dipole moment (μ).28 The electronic parameters were calculated using the DFT-Koopmans theorem, where A = −ELUMO and I = −EHOMO. The formulas that followed were used to determine the remaining variables:
ΔE = ELUMO − EHOMO | (6) |
![]() | (7) |
![]() | (8) |
![]() | (9) |
![]() | (10) |
The Pearson method formula is employed to determine the quantity of electrons exchanged:
![]() | (11) |
The DFT-derived values for the φFe of Fe(100), Fe(110), and Fe(111) surfaces are 3.89, 4.82, and 3.78 eV, respectively. For this analysis, only the Fe(110) surface is considered, as it exhibits greater stability and a more densely packed structure.29
Einteraction = Etotal − (Esurface+solution + Einhibitor) | (12) |
Ebinding = −Einteraction | (13) |
Inhibitor | Conc. (mmol L−1) | CR (mg cm−2 h−1) | ηw (%) |
---|---|---|---|
Blank | — | 0.681 | — |
[C4Py]Br | 0.25 | 0.360 | 47.2 |
0.5 | 0.320 | 53.1 | |
1 | 0.266 | 60.9 | |
2 | 0.202 | 70.3 | |
4 | 0.198 | 70.9 | |
8 | 0.197 | 71.1 | |
[C8Py]Br | 0.25 | 0.269 | 60.5 |
0.5 | 0.233 | 65.8 | |
1 | 0.193 | 71.6 | |
2 | 0.146 | 78.5 | |
4 | 0.141 | 79.3 | |
8 | 0.133 | 80.4 | |
[C12Py]Br | 0.25 | 0.193 | 71.6 |
0.5 | 0.160 | 76.5 | |
1 | 0.133 | 80.4 | |
2 | 0.080 | 88.3 | |
4 | 0.060 | 91.2 | |
8 | 0.046 | 93.2 |
Based on these results, 2 mmol L−1 was selected as the optimal concentration for subsequent studies, balancing inhibition efficacy and practical applicability.
Temperature plays a crucial role in affecting the corrosion process of metals. As shown in Fig. 4 and Table 4, four different temperature ranges were investigated. A detailed analysis shows that as the temperature rises, there is a gradual reduction in inhibition efficiency. This phenomenon can be attributed to the fact that elevated temperatures increase the average kinetic energy of the reacting molecules, leading to the desorption of inhibitor molecules from the surface of mild steel.36,37
Inhibitor | Temp (°C) | Temp (K) | CR (mg cm−2 h−1) | ηw (%) |
---|---|---|---|---|
Blank | 25 | 298 | 0.681 | — |
35 | 308 | 0.735 | ||
45 | 318 | 0.793 | ||
55 | 328 | 0.868 | ||
[C4Py]Br | 25 | 298 | 0.202 | 70.7 |
35 | 308 | 0.213 | 68.7 | |
45 | 318 | 0.222 | 67.3 | |
55 | 328 | 0.233 | 65.8 | |
[C8Py]Br | 25 | 298 | 0.146 | 78.5 |
35 | 308 | 0.152 | 77.7 | |
45 | 318 | 0.162 | 76.1 | |
55 | 328 | 0.169 | 75.2 | |
[C12Py]Br | 25 | 298 | 0.082 | 87.9 |
35 | 308 | 0.093 | 86.3 | |
45 | 318 | 0.104 | 84.7 | |
55 | 328 | 0.116 | 83.0 |
![]() | ||
Fig. 5 Tafel curves for mild steel in 1 M HCl solution containing without and with different concentrations of [CnPy]Br ILs at 298 K. |
Inhibitor | Conc. (mmol L−1) | Ecorr (V) | Icorr (μA cm−2) | βa (V dec−1) | −βc (V dec−1) | ηp (%) |
---|---|---|---|---|---|---|
Blank | — | −0.455 | 373.8 | 13.694 | 7.441 | — |
[C4Py]Br | 0.25 | −0.465 | 111.8 | 15.206 | 7.024 | 70.1 |
0.5 | −0.467 | 74.26 | 14.989 | 6.793 | 80.1 | |
1 | −0.445 | 60.59 | 15.359 | 6.522 | 83.8 | |
2 | −0.456 | 39.15 | 9.218 | 4.333 | 89.5 | |
4 | −0.458 | 36.64 | 8.869 | 4.791 | 90.2 | |
8 | −0.453 | 35.84 | 9.402 | 4.728 | 90.4 | |
[C8Py]Br | 0.25 | −0.466 | 75.15 | 14.077 | 7.302 | 79.9 |
0.5 | −0.474 | 45.69 | 8.357 | 5.916 | 87.8 | |
1 | −0.460 | 43.4 | 8.814 | 4.519 | 88.4 | |
2 | −0.467 | 33.48 | 8.498 | 5.214 | 91.0 | |
4 | −0.475 | 29.43 | 8.003 | 5.888 | 92.1 | |
8 | −0.484 | 26.72 | 7.739 | 5.079 | 92.9 | |
[C12Py]Br | 0.25 | −0.468 | 51.49 | 8.504 | 5.827 | 86.2 |
0.5 | −0.454 | 38.72 | 8.273 | 4.641 | 89.6 | |
1 | −0.466 | 32.72 | 8.635 | 5.101 | 91.2 | |
2 | −0.465 | 28.5 | 8.476 | 5.861 | 92.4 | |
4 | −0.471 | 26.18 | 8.689 | 6.287 | 93.0 | |
8 | −0.485 | 22.14 | 7.854 | 6.138 | 94.1 |
The anodic Tafel slope (βa) and cathodic Tafel slope (−βc) are sensitive indicators of the inhibition degree of anodic dissolution and cathodic hydrogen evolution reactions, respectively. For [CnPy]Br inhibitors, we observed that: for anodic reactions, the βa value decreases with the increase of inhibitor concentration(e.g., for [C12Py]Br, βa decreased from 13.694 V dec−1 in the blank to 7.854 V dec−1 at 8 mmol L−1), indicating a moderate suppression of Fe oxidation (anodic dissolution). For cathodic reactions, the −βc values decreased slightly (e.g., for [C12Py]Br, −βc decreased from 7.441 V dec−1 in the blank to 6.138 V dec−1 at 8 mmol L−1), suggesting a stronger inhibition of the cathodic hydrogen evolution reaction (2H+ + 2e− → H2). This indicates that longer alkyl chains amplify this differential inhibition. [C12Py]Br showed a more pronounced decrease in −βc (18.8% reduction compared to the blank) than the change in βa (42.7% reduction), indicating that the longer hydrophobic alkyl chain preferentially hinders the diffusion of H+ to the metal surface (cathodic reaction) by forming a denser protective film, while still suppressing anodic dissolution through adsorption on active Fe sites. This differential inhibition confirms that [CnPy]Br acts as a mixed-type inhibitor with a slight preference for cathodic inhibition. The pyridinium ring adsorbs on Fe active sites to suppress anodic dissolution via coordinate bonds, while the alkyl chain forms a hydrophobic barrier that more effectively blocks H+ transport to the cathode, leading to stronger cathodic inhibition.
As seen in Fig. 5, in the presence of inhibitors, the curves of polarization exhibit a shift in the E_corr values, either to more anodic or cathodic orientations, relative to the blank solution. As the concentration of [CnPy]Br ILs increases, irregular variations in the Ecorr values are observed. This indicates that the inhibitors under investigation influence both the anodic and cathodic reactions.29
When a compound's change in Ecorr is less than 85 mV, it is typically regarded as a mixed-type inhibitor. On the other hand, if the shift in Ecorr surpasses 85 mV relative to the Ecorr of the uninhibited solution, the corrosion type will be categorized as either anodic or cathodic. In this research, the Ecorr values for the three ILs at varying concentrations are all below 85 mV, suggesting that these ILs fall under the mixed-type inhibitor category.38,39
Table 5 unequivocally demonstrates that Icorr is significantly reduced when [CnPy]Br inhibitors are present. The values of Icorr decreased from 373.8 μA cm−2 in the blank solution to 35.84 μA per cm2 ([C4Py]Br), 26.72 μA per cm2 ([C8Py]Br), and 22.14 μA per cm2 ([C12Py]Br) at a concentration of 8 mmol L−1, the corresponding ηp reached 90.4%, 92.9%, and 94.1%. Because of the increased adsorption of [CnPy]Br ILs molecules onto the active sites at the mild steel/solution interface as the amount of the inhibitor increases, there is a positive correlation between the concentration of added inhibitors and the ηp.36 Under identical concentration conditions, the inhibitory effectiveness escalates as the length of the cationic alkyl chain extends, corroborating the findings from weight loss assays. [CnPy]Br ILs efficiently inhibit both cathodic and anodic corrosion reactions, as seen by Table 5, which shows variations in the βa and −βc values with the introduction of inhibitors.
![]() | ||
Fig. 6 The Nyquist plots were recorded at 298 K for metal specimens in 1 M HCl with and without different doses of [CnPy]Br inhibitor. |
![]() | ||
Fig. 7 Bode graphs of metal samples at 298 K in 1 M HCl with and without varying [CnPy]Br inhibitor. |
The Nyquist plots (Fig. 6) distinctly show a flattened capacitive loop, suggesting that the corrosion mechanism is mainly controlled by the charge transfer at the interface between the mild steel surface and the solution. The imperfect semicircular shape of the capacitive loops is usually ascribed to frequency dispersion brought on by the mild steel's surface roughness and heterogeneity. The Nyquist diagrams display irregular semicircles, suggesting that a simple capacitance model, such as Cdl, is not suitable for direct calculations. Consequently, a constant phase element (CPEdl) is added to the circuit in its place.40 However, the similarity in the shapes of the curves for all tested inhibitors implies that the dissolution mechanism remains consistent, regardless of the presence or absence of [CnPy]Br ILs. The diameter of the capacitive loops is considerably increased by the inclusion of the synthesized ILs in comparison to the blank solution, as shown in Fig. 6, and the loop diameter increases as the inhibitor concentration increases. The diameter grows greater as the length of the cation chain increases at every concentration level.37 The increasing diameter of the capacitive loops indicates a rise in charge transfer resistance (Rct), which suggests that the inhibitor is effectively hindering the electron transfer between the metal surface and the corrosive medium. This behavior reflects the formation of a more stable and compact protective layer on the steel surface, thereby enhancing corrosion resistance. Such an increase in Rct is commonly associated with improved surface coverage and stronger adsorption of the inhibitor molecules.
As shown in Fig. 7, the phase angle observed for the blank medium was noticeably lower when compared to the values at different inhibitor concentrations. The incorporation of [CnPy]Br ILs resulted in elevated phase angle values, signifying an enhanced capacitive nature and diminished surface deterioration. This suggests that a protective layer has developed at the metal-solution contact.41,42 The |Z| value at low frequencies increases progressively with the rise in inhibitor concentration. This is explained by the improved resistance to polarization across the metal/solution interaction. Moreover, this leads to a reduction in the roughness of the metal surface. This further supports the idea that the [CnPy]Br ILs adsorb onto the surface. The impedance plot indicates that the system exhibits a single time constant. From the previous analysis, the experimental data can be fitted using a Randles equivalent circuit (Fig. 8), as illustrated in Fig. 7.30 Electrolyte resistance (Rs), charge transfer resistance (Rct), and CPEdl are the components of the circuit diagram, respectively. An accurate fit was obtained by analyzing the impedance using the CPEdl.43,44 The following formula is used to determine the capacitance of both layers at the metal/solution interface:
Cdl = T(1/n) × R(1−n)ns | (14) |
The following formula is used to determine the impedance related to the constant phase element (CPEdl):
![]() | (15) |
The parameters obtained through fitting with ZView software are presented in Table 6. As the dosage of corrosion inhibitors increases, the values of Cdl exhibit a consistent decrease. According to the Helmholtz model, this behavior can be explained by either a rise in the electric double layer's thickness or a decrease in local permittivity. This suggests that corrosion inhibitor molecules adsorb onto the metal surface, displacing the previously adsorbed water molecules or chloride ions.46 The results from the kinetic potential polarization method and the EIS measurements show a high degree of consistency in the fitting data. This demonstrates that mild steel can be effectively protected against corrosion by [CnPy]Br ILs in 1 M HCl.28
Inhibitor | Conc. (mmol L−1) | Rs (Ω cm2) | n | Rct (Ω cm2) | Cdl (μF cm−2) | χ2 × 10−3 | ηE (%) |
---|---|---|---|---|---|---|---|
Blank | — | 1.659 | 0.827 | 48.04 | 42.3 | 0.44 | — |
[C4Py]Br | 0.25 | 1.500 | 0.861 | 161.2 | 22.1 | 0.56 | 70.2 |
0.5 | 1.440 | 0.844 | 216.1 | 16.2 | 1.10 | 77.8 | |
1 | 1.308 | 0.851 | 327.2 | 15.2 | 1.34 | 85.3 | |
2 | 1.119 | 0.844 | 367 | 11.6 | 1.45 | 86.9 | |
4 | 0.969 | 0.847 | 431.6 | 9.98 | 0.76 | 88.9 | |
8 | 0.954 | 0.849 | 595.3 | 8.51 | 0.66 | 91.9 | |
[C8Py]Br | 0.25 | 1.484 | 0.934 | 267.3 | 44.7 | 0.36 | 82.0 |
0.5 | 1.235 | 0.811 | 380.2 | 25.1 | 0.80 | 87.4 | |
1 | 1.158 | 0.802 | 547.3 | 19.3 | 0.48 | 91.2 | |
2 | 1.078 | 0.816 | 627.8 | 16.8 | 0.33 | 92.4 | |
4 | 1.077 | 0.780 | 680.5 | 13.4 | 0.53 | 92.9 | |
8 | 0.942 | 0.838 | 757.3 | 7.6 | 0.32 | 93.7 | |
[C12Py]Br | 0.25 | 1.404 | 0.846 | 324 | 36.4 | 0.23 | 85.2 |
0.5 | 1.334 | 0.833 | 527.7 | 23.7 | 0.65 | 90.9 | |
1 | 1.277 | 0.808 | 637.3 | 15.9 | 0.89 | 92.5 | |
2 | 1.203 | 0.789 | 719.8 | 13.7 | 1.46 | 93.3 | |
4 | 1.183 | 0.790 | 761.9 | 12.8 | 0.93 | 93.7 | |
8 | 1.142 | 0.792 | 832.1 | 12.5 | 0.75 | 94.2 |
![]() | (16) |
ΔGads = −RT![]() | (17) |
Inhibitor | R2 | Slope | Kads | ΔGads (kJ mol−1) |
---|---|---|---|---|
[C4Py]Br | 0.99987 | 1.098 | 17![]() |
−34.213 |
[C8Py]Br | 0.99995 | 1.073 | 27![]() |
−35.257 |
[C12Py]Br | 0.99987 | 1.066 | 28![]() |
−35.348 |
The three ILs under investigation have a considerable potential to adsorb on the metal-electrolyte interfaces, as shown by the comparatively high Kads values in Table 7.21 Additionally, when the alkyl chains length on the cation increases, the Kads values increase as well. The powerful connection between corrosion inhibitor molecules and the metal surface, as well as the spontaneous nature of adsorption, are highlighted by the negative ΔGads values, which reveal a thermodynamically advantageous process. According to previous research, it is understood that when the ΔGads value reaches −20 kJ mol−1, the inhibitor's adsorption is classified as physisorption. Conversely, when the ΔGads value exceeds −40 kJ mol−1, the adsorption process is considered chemisorption. In this research, the ΔGads values of three studied ILs fall within the range of −34 to −36 kJ mol−1, substantiating that [CnPy]Br would adsorb on the metal surface by both physisorption and chemisorption. Thus, in a 1 M HCl solution, [CnPy]Br exhibits a mixed adsorption behavior on the metal surface.49
![]() | ||
Fig. 10 SEM images and corresponding EDS spectra are shown for (a and d) mild steel prior immersion in HCl, (b and e) after 24 h of immersion in HCl without the inhibitor, (c and f) with [C12Py]Br. |
![]() | ||
Fig. 11 XPS spectra of mild steel surface in presence of [C12Py]Br: (a) full scanned spectrum; (b) C 1s; (c) O 1s; (d) Fe 2p. |
The conventional peaks of Fe, O, and C are shown in Fig. 11(a), which is clearly caused by the process of adsorption of the [C12Py]Br inhibitor. The C 1s spectra are shown in Fig. 11(b), with three deconvoluted peaks at 284.8 eV, 286.4 eV, and 288.6 eV. The initial peak, observed around 284.8 eV, corresponds to the C–H and C–C bonds. The peak at 286.4 eV (C–N bonds) directly indicates interactions between the pyridinium ring of [C12Py]Br and the Fe surface. This bond forms via electron donation from the pyridinium N atom to Fe's empty d-orbitals, providing critical evidence for chemical adsorption—consistent with quantum chemical calculations showing high electron density on the pyridinium ring (Section 3.6). Though the N 1s signal is undetectable (due to low surface coverage), the C–N bond in C 1s spectra confirms the pyridinium ring's involvement in surface binding. Meanwhile, the second peak, appearing at 286.4 eV, is attributed to either C–O or C–N bonds. The final fitted peak, positioned near 288.6 eV, is characterized by a high binding energy and is likely associated with the CO group.52 As shown in Fig. 11(c), the peak of O 1s spectrum is also fitted into three peaks. The three fitted peaks are 530.1 eV, 531.9 eV, and 532.1 eV, respectively. The initial peak is assigned to O2−, theoretically linked to the bond between Fe3+ and Fe2O3. Meanwhile, the second peak corresponds to OH−, indicative of its presence in hydrous iron oxides such as FeOOH. The final peak is likely associated with the oxygen in adsorbed H2O molecules. For the Fe 2p spectra (Fig. 11(d)), the peaks at 711.2 eV and 713.8 eV are attributed to ferric compounds such as Fe2O3 and FeOOH, respectively.53 The peak at 716.2 eV associated with the satellite of Fe(III) and the peak at 725 eV is attributed to the presence of various iron species, including Fe3O4, α-F2O3, and FeOOH. Based on XPS analysis, it is evident that [C12Py]Br molecules undergo chemisorption onto the mild steel surface, a finding that aligns with the conclusions drawn from the thermodynamic study. Consequently, in acidic conditions, [CnPy]Br ILs demonstrate exceptional effectiveness as mild-steel inhibitors of corrosion.51
Based on frontier orbital theory, the HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) provide insights into potential adsorption sites for inhibitor molecules and help forecast the strength and stability of their interactions with metal surfaces.54,55 In general, the EHOMO indicates an inhibitor's ability to transfer electrons to the metal surface, while the ELUMO reflects the inhibitor's tendency to accept electrons from the metal. The energy gap (ΔE) was shown to correspond with the rate of corrosion. As ΔE increases, molecular stability improves, reducing the likelihood of adsorption. Conversely, reduced stability often improves the capacity of molecules to adhere to metal surfaces, effectively preventing corrosion. Table 8 shows the calculated quantum chemical descriptors. The results reveal that the energy gap follows this sequence: [C12Py]Br < [C8Py]Br < [C4Py]Br. This configuration indicates that [C12Py]Br is expected to show higher reactivity toward the metal surface, resulting in improved corrosion inhibition, which corresponds well with the experimentally observed trend in efficiency.56 Fig. 12 illustrates the optimized geometries, frontier molecular orbital distributions, and molecular electrostatic potential maps of the [CnPy]Br ionic liquids. The findings from Fig. 12 reveal that the HOMO orbitals are primarily localized near the anions, enabling these regions to donate electrons to the metal surface and establish coordination bonds. Specifically, the HOMO is concentrated on the pyridinium ring (especially the N atom and adjacent C atoms), confirming this ring as the key electron-donating site for coordination with Fe. The MEP map (Fig. 12) further shows negative electrostatic potential on Br−, which facilitates its pre-adsorption on the positively charged Fe surface via electrostatic interactions—explaining why Br− enhances [CnPy]+ adsorption (supported by EIS data in Table 6 showing reduced Cdl with increasing [CnPy]Br concentration). These results clarify that both the pyridinium ring (chemical adsorption) and Br− (electrostatic attraction) are indispensable functional groups. The LUMO orbitals are mainly located around the cationic ring structure, which may serve as acceptors for electrons on the metal surface, forming feedback bonds. These regions correspond to the active sites of the corrosion inhibitors being studied.
Inhibitor | EHOMO (eV) | ELUMO (eV) | ΔE (eV) | I (eV) | A (eV) | μ (D) | χ (eV) | η (eV) | σ (eV) | ΔN (eV) |
---|---|---|---|---|---|---|---|---|---|---|
[C4Py]Br | −6.275 | −2.356 | 3.919 | 6.275 | 2.356 | −4.315 | 4.315 | 1.959 | 0.510 | 0.129 |
[C8Py]Br | −6.270 | −2.353 | 3.916 | 6.270 | 2.354 | −4.312 | 4.312 | 1.952 | 0.519 | 0.130 |
[C12Py]Br | −6.268 | −2.351 | 3.915 | 6.268 | 2.353 | −4.311 | 4.311 | 1.941 | 0.537 | 0.130 |
The calculations show that [C12Py]Br exhibits the lowest hardness values compared to [C4Py]Br and [C8Py]Br. It was observed that the inhibitor exhibiting the lowest global hardness (and consequently the highest global softness) is anticipated to demonstrate the greatest inhibition efficiency. Our findings reveal that [C12Py]Br, characterized by the highest softness value (0.537 eV), achieves the most effective inhibition performance. In contrast, [C4Py]Br, with the lowest softness value (0.510 eV), exhibits relatively weaker inhibition efficiency, as shown in Table 8.
Electronegativity (χ) is a crucial parameter to consider when comparing corrosion inhibitors with similar molecular structures. As per Pearson's theory, when two systems interact, electrons will flow from an area of higher chemical potential to one of lower potential. This process continues until equilibrium is achieved in the chemical potential. The work function for Fe(110) (4.8 eV) is greater than the electronegativity of all three [CnPy]Br molecules. Consequently, the [CnPy]Br molecules primarily adsorb onto the Fe(110) surface, facilitated by electron transfer from a metal orbital to an available free orbital. The electronegativity values of the molecules examined exhibit the trend [C12Py]Br < [C8Py]Br < [C4Py]Br, which aligns with the sequence of inhibitory effectiveness observed in the electrochemical experiments.28
The stability and reactivity of corrosion inhibitor molecules are largely determined by the rate of electron transfer (ΔN). Theoretically, the corrosion inhibitor molecules have a stronger internal affinity when the computed value of ΔN is greater than 0. Consequently, these molecules' electrons are moved to the metal atoms, exhibiting improved inhibitory qualities. The corrosion inhibitor molecules have enough electron-donating capacity to successfully stop mild steel from corroding if ΔN is less than 3.6.57,58 As presented in Table 8, the calculated values of ΔN for the [CnPy]Br molecules investigated in this study are all above zero but do not exceed 3.6. This indicates that the Fe atom's unoccupied d orbitals can receive electron donations from the [CnPy]Br molecules.
The molecular electrostatic potential (ESP), which provides a fundamental representation of a molecule's electronic structure, has been effectively used to identify and label the electrophilic sites of corrosion inhibitors.28 While conductive interactions are more likely to occur in the blue (positive) region of the ESP, nucleophilic attacks can occur in the red (negative) region. Fig. 12 shows that the ESP map identifies the anionic group as the primary electrophilic active site. These regions can therefore engage in electron-transferring interactions with the working electrode's surface. Conversely, the areas with the highest positive electrostatic potential are found to be the nitrogen atoms in the [CnPy]Br molecules.
![]() | ||
Fig. 13 Side and top views of stabilized adsorption configurations of corrosion inhibitors on the Fe(110) surface. |
In this study, the Fe(110) surface was chosen for DFT modeling due to its low surface energy and high stability, making it a dominant orientation in bcc Fe. Although real mild steel surfaces are polycrystalline and often partially oxidized, the Fe(110) facet serves as a reasonable approximation to study the intrinsic adsorption and inhibition behavior at the atomic scale.
The MD results for Einteraction and Ebinding are shown in Table 9, where all three [CnPy]Br molecules have negative Einteraction values. This suggests that each of these molecules strongly interacts with the Fe(110) surface. The corrosion inhibitor molecule has a greater propensity to contact intimately with the metal surface atoms if its Einteraction value is small.6,45 The corrosion inhibitor [C12Py]Br shows significant stability on the Fe(110) surface in an acidic media environment with strong contacts, as seen by Table 9, which shows that [C4Py]Br > [C8Py]Br > [C12Py]Br, i.e., [C12Py]Br has greater inhibitory properties. One important metric for assessing how strongly corrosion inhibitor molecules connect with metal surfaces is Ebinding. Table 9 reveals that [C12Py]Br exhibits a high positive value, reflecting its significant adsorption capabilities and high efficacy in corrosion inhibition. The higher Ebinding of [C12Py]Br (390.473) compared to [C8Py]Br (343.673) and [C4Py]Br (315.581) (Table 9) is attributed to stronger van der Waals interactions between its longer alkyl chain and the Fe(110) surface. This is validated by the RDF (Fig. 15), which shows a peak at 3.0–4.0 Å for C–H⋯Fe interactions—indicating enhanced physical adsorption of longer alkyl chains. Thus, the alkyl chain length directly modulates the hydrophobic barrier effect, a mechanism that cannot be deduced without MD simulations. Fig. 14 illustrates how temperature changes, with smaller Ebinding values occurring at higher temperatures. This variation can be mainly attributed to the inhibitor compound's decreasing attraction to the metal surface, leading to a repulsive effect as temperature increases.
Inhibitor | Einteraction (kcal mol−1) | Ebinding (kcal mol−1) |
---|---|---|
[C12Py]Br | −390.473 | 390.473 |
[C8Py]Br | −343.673 | 343.673 |
[C4Py]Br | −315.581 | 315.581 |
![]() | ||
Fig. 14 Calculation of adsorption energy of [C12Py]Br ILs on Fe(110) surface in aqueous HCl solution at different temperatures. |
The radial distribution function (RDF) method is employed to identify chemical bonds, or both bond types, while focusing on the structural analysis of the results obtained from MD simulations.37 The peaks seen in Fig. 15 can be used to infer the interaction between the Fe(110) surface and the [CnPy]Br ILs. The initial peak, observed within the 1–3.5 Å range, suggests chemical adsorption. This indicates that the N atoms in the [CnPy]Br molecules form strong interactions with the Fe(110) surface through hydrogen bonding or chemical bonding.59 On the mild steel surface, [CnPy]Br also demonstrates physical adsorption, as evidenced by the peaks observed at a distance of more than 3.5 Å.
During the corrosion inhibition process, chloride ions (Cl−) initially adsorb onto the metal surface, creating negatively charged regions. These regions attract protonated cations through electrostatic interactions, a mechanism known as physical adsorption. This interaction reduces the active sites available for the cathodic hydrogen evolution reaction, thereby slowing the overall corrosion rate. Inhibitor molecules form coordinate bonds by donating lone pair electrons from their heteroatoms to the vacant d-orbitals of iron. This interaction establishes a protective layer that shields the metal from corrosive environments, with effectiveness influenced by the inhibitor's structure and environmental conditions. In addition, the π-electron clouds on the aromatic ring engage in donor–acceptor interactions (π backbonding) with the ionized Fe atoms on the surface. The combined interaction of donation and π backbonding promotes the creation of protective layers on the metal surface, effectively safeguarding it from corrosion. The large surface area of [C12Py]Br promotes the formation of an extensive protective layer, which effectively reduces metal interaction with corrosive environments. The corrosion inhibition mechanism of [CnPy]Br relies on synergism between three functional groups: (1) the pyridinium ring forms coordinate bonds with Fe via N → Fe electron donation (chemisorption), confirmed by XPS C–N bonds and DFT HOMO localization. (2) Br− pre-adsorbs on Fe, creating a negatively charged layer that attracts [CnPy]+ (electrostatic interaction), as supported by EIS Cdl trends and MEP maps. (3) Longer alkyl chains (e.g., C12) enhance physical adsorption via stronger C–H⋯Fe van der Waals interactions (MD RDF data), forming a denser hydrophobic barrier. This multi-functional synergy is novel and cannot be inferred without experimental (weight loss, XPS) and computational (DFT, MD) evidence. Its adsorption properties and molecular structure further enhance its role as a highly efficient corrosion inhibitor.41,60 When compared with [C4Py]Br and [C8Py]Br, this is probably the main reason why this molecule exhibits better inhibitory activity.
(2) Electrochemical testing revealed that the three added ILs corrosion inhibitors can effectively reduce the corrosion current density in acidic corrosive media, indicating their significant corrosion inhibition performance. In acidic environments, these ILs corrosion inhibitors act as mixed corrosion inhibitors; in acidic media, the corrosion inhibitor achieved a maximum inhibition efficiency of 94.1%. The adsorption behavior of all ILs corrosion inhibitors follows the Langmuir isotherm adsorption model. The Gibbs free energy analysis of corrosion inhibitor adsorption reveals that its adsorption mechanism is a combination of physical and chemical adsorption.
(3) The results of surface analysis technology show that the samples with added corrosion inhibitors exhibit significant corrosion inhibition effects, and trace element analysis confirms the formation of a protective corrosion inhibitor film on the surface of low carbon steel. XPS further confirmed the binding effect of physical adsorption and chemical adsorption by detecting the composition of the sample surface.
(4) The quantum chemistry calculation results indicate that the electrons of the three ILs corrosion inhibitors are mainly concentrated on nitrogen atoms, which helps to determine the adsorption sites on the surface of low carbon steel, and further supports the experimental results through the analysis of frontier orbital distribution and reaction activity. The molecular dynamics simulation results indicate that ILs can form stable adsorption on the surface of low carbon steel through interfacial donor acceptor interactions, and the inhibition sequence obtained from the binding energy results is consistent with the experimental results. Notably, our work clarifies the distinct roles of the pyridinium ring (chemical adsorption), Br− (electrostatic pre-adsorption), and alkyl chain (physical barrier) in [CnPy]Br—revealing a synergistic mechanism unreported in prior studies. Key findings (e.g., the critical role of C12 alkyl chain and C–N bond formation) depend on experimental data and calculations, as they cannot be intuitively deduced. The radial distribution function (RDF) calculation confirmed that the adsorption type is a mixed adsorption mainly composed of chemical adsorption and supplemented by physical adsorption.
This journal is © The Royal Society of Chemistry 2025 |