# Diaryl Sulfide Derivatives as Potential Iron Corrosion Inhibitors: A Computational Study

Morad El-Hendawy, Asmaa Kamel, Mahmoud Mohamed, Rabah Boukherroub, Jacek Ryl, Mohammed Amin

► To cite this version:

Morad El-Hendawy, Asmaa Kamel, Mahmoud Mohamed, Rabah Boukherroub, Jacek Ryl, et al. Diaryl Sulfide Derivatives as Potential Iron Corrosion Inhibitors: A Computational Study. Molecules, 2021, 26 (20), pp.6312. 10.3390/molecules26206312. hal-03457664

HAL Id: hal-03457664

https://hal.science/hal-03457664

Submitted on 25 May 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Diaryl Sulfide Derivatives as Potential Iron Corrosion Inhibitors: A Computational Study

Morad M. El-Hendawy <fnref n="1" />,* ID, Asmaa M. Kamel <fnref n="1" />, Mahmoud M. A. Mohamed <fnref n="1" />, Rabah Boukherroub <fnref n="2" /> ID,
Jacek Ryl <fnref n="3" />,* ID and Mohammed A. Amin <fnref n="4" />,*

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** El-Hendawy, M.M.; Kamel, A.M.; Mohamed, M.M.A.; Boukherroub, R.; Ryl, J.; Amin, M.A. *Diaryl Sulfide Derivatives as Potential Iron Corrosion Inhibitors: A Computational Study.* *Molecules* **2021**, *26*, 6312. https://doi.org/10.3390/molecules26206312

Academic Editor: Bryan M. Wong

Received: 10 September 2021

Accepted: 14 October 2021

Published: 19 October 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** The present work aimed to assess six diaryl sulfide derivatives as potential corrosion inhibitors. These derivatives were compared with dapsone (4,4'-diaminodiphenyl sulfone), a common leprosy antibiotic that has been shown to resist the corrosion of mild steel in acidic media with a corrosion efficiency exceeding 90%. Since all the studied compounds possess a common molecular backbone (diphenyl sulfide), dapsone was taken as the reference compound to evaluate the efficiency of the remainder. In this respect, two structural factors were examined, namely, (i) the effect of replacement of the S-atom of diaryl sulfide by SO or SO₂ group, (ii) the effect of the introduction of an electron-withdrawing or an electron-donating group in the aryl moiety. Two computational chemical approaches were used to achieve the objectives: the density functional theory (DFT) and the Monto Carlo (MC) simulation. First, B3LYP/6-311+G(d,p) model chemistry was employed to calculate quantum chemical descriptors of the studied molecules and their geometric and electronic structures. Additionally, the mode of adsorption of the tested molecules was investigated using MC simulation. In general, the adsorption process was favorable for molecules with a lower dipole moment. Based on the adsorption energy results, five diaryl sulfide derivatives are expected to act as better corrosion inhibitors than dapsone.

**Keywords:** diaryl sulfides; DFT; MC simulation; corrosion inhibitor; dapsone

## 1. Introduction

Corrosion is a serious problem, as it threatens our life. Corrosives components can cause damage to not only metals but also to the human digestive and respiratory tracts, eyes, and skin. Using coating materials as inhibitors can lengthen the life and usability of metallic components, machinery, products, etc. Continuous advances in computer efficiency and computational chemistry software allowed examining different types of corrosion inhibitors. Density Functional Theory (DFT) and Monte Carlo (MC) simulation methods have recently been accepted as fast and powerful tools for estimating the relative corrosion inhibition activities of various molecules [1,2]. Since corrosion processes and their prevention by organic inhibitors is a very active research area, many researchers have focused their studies on this subject. Many researchers reported that some physicochemical and electronic properties of the organic inhibitor's functional groups, steric effects, electronic density of donor atoms, and orbital character of donating electrons affect the inhibition efficiency of metals [3,4].

Several organic and inorganic corrosion inhibitors are expensive, toxic, and harmful to the environment, restricting their use. Currently, researchers shifted their focus on expired

<footnote n="1">¹ Department of Chemistry, Faculty of Science, New Valley University, Kharga 72511, Egypt; asmaa.kamel@scinv.au.edu.eg (A.M.K.); mmhm802004@yahoo.com (M.M.A.M.)</footnote>

<footnote n="2">² University of Lille, CNRS, Centrale Lille, Université Polytechnique Hauts-de-France, UMR 8520-IEMN, F-59000 Lille, France; rabah.boukherroub@univ-lille.fr</footnote>

<footnote n="3">³ Institute of Nanotechnology and Materials Engineering, Gdansk University of Technology, Narutowicza 11/12, 80-233 Gdansk, Poland</footnote>

<footnote n="4">⁴ Department of Chemistry, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia</footnote>

<footnote>* Correspondence: morad.elhendawy@aucegypt.edu (M.M.E.-H.); jacek.ryl@pg.edu.pl (J.R.); mohamed@tu.edu.sa (M.A.A.)</footnote>
drugs as a safe alternative for protecting metals if they meet the following criteria [5]: (i) the chemical composition of the drug must have active centers such as N, O, and S atoms, (ii) it must be less hazardous and environmentally friendly, and (iii) it can be easily produced and purified. Several research studies reported different types of drugs (melatonin, cephapirin, tramadol, etc.) as corrosion inhibitors for various metals [6–10]. The findings indicated that these drugs form an insoluble complex on the metal surface, protecting it from corrosion.

Recently, Danaee et al. [6] studied experimentally and theoretically the inhibitory effect of an oxethazaine, 2,2'-(2-hydroxyethylimino)bis[N-(alphaalphadimethylphenethyl)-N-methylacetamide] drug on corrosion of mild steel in 1 M HCl. Theoretical results were found to be compatible with the experimental findings. Additionally, Sachin et al. demonstrated that L-Dopa, a drug used to treat Parkinson's disease, significantly minimized the corrosion of mild steel in both HCl and H₂SO₄. The findings indicate that inhibition results from the formation of a non-porous organic film on the metal surface, and the adsorption follows the Langmuir isotherm [7].

Dohare et al. [8] demonstrated the inhibitory effect of tramadol for protecting mild steel in acidic media by electrochemical impedance spectroscopy, potentiodynamic polarization, surface analysis (SEM and AFM), and DFT methods. The influence of cephapirin drugs on the corrosion of carbon steel was examined by weight loss and electrochemical methods [10]. EIS studies revealed that the inhibition process occurred through charge transfer. The quantum chemical calculations proved that cephapirin is a suitable corrosion inhibitor against carbon steel. Overall, research studies revealed that drugs as corrosion inhibitors are the best choice for protecting metals and lead to the management of expired drug waste.

Diaryl sulfide compounds and their sulfoxide and sulfone derivatives have long attracted the interest of researchers as effective corrosion inhibitors [11–16]. In 1980, Abelaal et al. demonstrated the anticorrosive inhibition of diphenyl sulfide (Ph₂S), diphenyl sulfoxide (Ph₂SO), and diphenyl sulfone (Ph₂SO₂) against tin and cadmium in acidic solutions [11]. The adsorption of inhibitors obeyed a Tomkin isotherm. The inhibiting effect decreased in the series: Ph₂S > Ph₂SO₂ > Ph₂SO. Earlier, the inhibiting efficiencies of some sulfoxides with aliphatic chains or aromatic rings were investigated in an acidic medium [12]. The anticorrosive activity followed the sequence: dibenzyl sulfoxide > di-n-butyl sulfoxide > di-p-tolyl sulfoxide > diphenyl sulfoxide > tetramethylene sulfoxide > dimethyl sulfoxide. The anticorrosive activity was influenced by the electron density of the sulfur atom and their rate of reduction to sulfide. Recently, James and Lalgudi reported in their patent that a metal-appended dichloro-diphenyl sulfone anti-corrosion additive could reduce the corrosion susceptibility of carbon steel and other metal substrates when added to topcoats or primer [13]. The influence of the addition of organic sulfur compounds (mercaptan, diphenyl sulfide, and diphenyl disulfide) was studied on the corrosion caused by elemental sulfur in synthetic naphtha [14]. The study indicated that diphenyl disulfide was the best corrosion inhibitor of elemental sulfur.

Although the mode of action of drugs within the mammalian body differs from their anticorrosive mechanism on a metal surface, a common factor uniting the two functions is the distinct electronic and geometric structures that facilitate chemical and physical attachment to a biological receptor or metal surface. For example, dapsone (4,4'-diamino diphenyl sulfone), the common antibiotic drug for treating leprosy, has proved its anticorrosive merit for mild steel in acidic media (1 M HCl and 0.5 H₂SO₄) [15]. The superiority of the anticorrosive activity of dapsone bearing the diphenyl sulfide backbone [15] motivated us to explore more compounds with the same molecular backbone. The finding revealed that inhibition occurs through drug adsorption on the metal surface without modifying the corrosion mechanism. Recently, Singh et al. reacted the expired drug, dapsone, with benzaldehyde and salicylaldehyde to produce dapsone-benzaldehyde and dapsone-salicylaldehyde, respectively [16]. The resulting products exhibited maximum corrosion inhibition efficiencies of 95.67% and 94.23% for mild steel at a concentration of 0.219 mM in an acidic medium. The addition of a tiny amount of KI further increased their efficiencies up to 99.03% and 97.98%, respectively. Additionally, diphenyl sulfides and their
sulfoxide and sulfone analogs proved their merit as excellent corrosion inhibitors [11–14]. The common structural factor among all these molecules is diphenyl sulfide, which reflects its superiority in anti-corrosion activity.

Sulfur-containing compounds, such as sulfides, sulfoxides, and sulfones, constitute a large portion of pharmaceutical drugs that possess various biological activities [17]. For example, Chan et al. assessed the antiviral activity against HIV-1 of a series of 2-amino-6-arylthiobenzonitriles and their sulfoxides and sulfones derivatives [18]. Based on the above, we chose six molecules from the long list of this series of 2-amino-6-arylthiobenzonitriles and their sulfoxide and sulfone derivatives, as shown in Figure 1, to examine their anticorrosive behavior on steel. We divided the studied molecules into two families to facilitate the discussion of results: X–CN and X–OMe, where X is S, SO, or SO₂ group. In this respect, the DFT method and MC simulation were applied in the absence and presence of iron metal, respectively. Several quantum chemical descriptors of the isolated inhibitors (Table 1) were calculated to study their global reactivity and selectivity toward other species. Furthermore, the adsorption descriptors and the mode of adsorption on the iron surface were calculated by MC simulation.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Chemical structures of investigated compounds with their IUPAC names as well as nicknames. The common molecular structure of the seven compounds is diphenyl sulfide.</figcaption></figure>

## 2. Computational Details

All DFT calculations were performed using Gaussian 16 suite [19]. In this context, the B3LYP/6-311+G(d,p) method was utilized for geometry optimization in the aqueous phase, followed by frequency calculations to confirm that the optimized geometry was not a saddle point. The polarizable continuum model (PCM) using the integral equation formalism variant (IEFPCM) was employed for modeling the aqueous phase [20]. Molecular properties related to the reactivity and selectivity of the compounds were estimated following Koopmans’s theorem [21] that depends on the energy of the frontier molecular orbitals. According to DFT-Koopmans’s theorem [21,22], the ionization potential (I) and electron affinity (A) can be approximated as the negative values of the E<sub>HOMO</sub> and E<sub>LUMO</sub>, respectively. The quantum chemical descriptors and the mathematical forms used for their estimation are collected in Table 1.

<table><caption>**Table 1.** The mathematical forms for estimating quantum chemical descriptors.</caption><thead><tr><th>Descriptor</th><th>Mathematical Form</th></tr></thead><tbody><tr><td>Ionization potential (I)</td><td>I = -E<sub>HOMO</sub></td></tr><tr><td>Electron affinity (A)</td><td>A = -E<sub>LUMO</sub></td></tr><tr><td>Energy gap (ΔE)</td><td>ΔE = E<sub>LUMO</sub> - E<sub>HOMO</sub></td></tr><tr><td>Electronegativity (χ) [23,24]</td><td></td></tr><tr><td>Chemical hardness (η) [23,24]</td><td></td></tr><tr><td>Global electrophilicity index (ω) [25]</td><td></td></tr><tr><td>The number of electrons transferred (ΔN) from the inhibitor to the iron surface. The work function of metal is the work function of iron [26]</td><td></td></tr><tr><td>The energy associated with a backing donation</td><td></td></tr></tbody></table>
Monte Carlo simulations [27] were conducted using the Adsorption Locator module [28] as implemented in BIOVIA Materials Studio 2017 [29] to locate the low energy adsorption sites of the potential corrosion inhibitors on the Fe surface. All components, such as clean iron surface, inhibitor, and the acidic medium, were optimized first using the Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies (COMPASS) force field [30]. Ewald and atom-based summation approaches were employed to calculate the electrostatic and van der Waals energies, respectively. The Fe(110) crystal is the most stable facet; thus, it was selected for the simulation [31]. An amount of 1 M HCl was commonly used in the experiment as an acidic medium; thus, 5 HCl/278 H₂O were used to represent its effect on adsorption mode [32]. To provide an appropriate surface for the interaction with the inhibitors, a 10 × 10 supercell was built with a 30 Å thick vacuum slab.

### 3. Results and Discussion

#### 3.1. *Quantum Chemical Study*

The geometry, electron density distribution of HOMO and LUMO, and ESP are presented in Figure 1. The findings indicate a bent molecular geometry of all studied molecules. The substitution of the S atom of diaryl sulfides by SO or SO₂ groups impacts the geometry as seen from the C–S bond length, <C–S–C bond angle, and <C–S–C–C dihedral angle (Table 2). For example, ongoing from S–CN to SO–CN and SO₂–CN, the bond length increases by 0.05 and 0.02 Å, and the bond angle changes by –7° and +1°, respectively. A tremendous and steady decrease in dihedral angles was observed ongoing from S–CN to SO–CN and SO₂–CN. Even small changes in the structural parameters can produce significant changes in the predicted energy and other molecular properties [33].

<table><caption>**Table 2.** Bond length, bond angle, and dihedral angle between sulfur atoms and sulfur-adjacent carbon atoms in the studied inhibitors.</caption><thead><tr><th></th><th>S–CN</th><th>SO–CN</th><th>SO<sub>2</sub>–CN</th><th>S–OMe</th><th>SO–OMe</th><th>SO<sub>2</sub>–OMe</th><th>Dapsone</th></tr></thead><tbody><tr><td>C<sub>M</sub>–S;</td><td>1.79;</td><td>1.84;</td><td>1.81;</td><td>1.79;</td><td>1.82;</td><td>1.79;</td><td>1.79;</td></tr><tr><td>C<sub>D</sub>–S (Å) *</td><td>1.79</td><td>1.84</td><td>1.82</td><td>1.79</td><td>1.84</td><td>1.82</td><td>1.79</td></tr><tr><td>&lt;C<sub>M</sub>–S–C<sub>D</sub> (°)</td><td>104</td><td>97</td><td>105</td><td>104</td><td>98</td><td>105</td><td>106</td></tr><tr><td>&lt;C<sub>M</sub>–C<sub>M</sub>–S–C<sub>D</sub> (°)</td><td>139</td><td>97</td><td>82</td><td>89</td><td>88</td><td>84</td><td>90</td></tr><tr><td>&lt;C<sub>D</sub>–C<sub>D</sub>–S–C<sub>M</sub> (°)</td><td>142</td><td>83</td><td>72</td><td>2</td><td>82</td><td>73</td><td>90</td></tr></tbody></table>

The frontier molecular orbital topologies of the studied inhibitors are illustrated in Figure 2. The HOMOs for all studied inhibitors are π-type molecular orbitals. The HOMO surface of the X–CN corrosion inhibitors is essentially localized over the disubstituted aryl moiety, indicating that it is electron-rich and may be involved in electron donation to the metal surface. Further, the LUMO plot shows a delocalization of the lowest vacant molecular orbitals over the entire skeleton, suggesting that the low-lying vacant orbitals of the inhibitors could be available for back-donation. A different situation is observed in the case of X–OMe inhibitors. For example, the HOMO surface of SO–OMe is mainly localized on the sulfoxide group and the mono-substituted aryl moiety, while the LUMO is shifted toward the other moiety. This behavior is a feature of the intramolecular charge transfer where the HOMO and LUMO are localized on the donor and acceptor moieties, respectively. Because of the symmetrical geometry of dapsone, the HOMO and LUMO electron density are delocalized over the entire molecular skeleton.

The molecular electrostatic potential (MEP) of a molecule represents its electrophilic and nucleophilic reactivities, according to the colored bar in Figure 2. For example, the color between red and green is susceptible to electrophilic attack; in contrast, green and blue are accessible to nucleophilic attack. At first glance, the most negative potential (red color) is located around the oxygen atoms of the SO and SO₂ groups, especially in the case

<footnote>* The symbols M and D in C<sub>M</sub> and C<sub>D</sub> refer to the carbon atom in the mono-substituted and disubstituted aryl moiety.</footnote>
of SO-OMe and SO₂-OMe, reflecting the role of the methoxy group at the para position to provide these atoms with excess electron density.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** The geometry, frontier molecular orbital electron density distribution, and electrostatic potential map of the investigated molecules. The colored bar of MEP indicates the positive and negative potentials.</figcaption></figure>

The global reactivity parameters of the studied inhibitors are gathered in Table 3. According to the FMO theory, the E<sub>HOMO</sub> and E<sub>LUMO</sub> are associated with the electron-donating and electron-accepting capacities of the molecule, respectively [34]. Therefore, a molecule with a higher E<sub>HOMO</sub> and a lower E<sub>LUMO</sub> is expected to adsorb strongly to the metallic surface. The results indicate that E<sub>HOMO</sub> and E<sub>LUMO</sub> steadily are reduced ongoing from the diaryl sulfides to diaryl sulfones. Thus, the electron-donating capacity decreases when the electron-accepting capabilities of the inhibitor increase in this direction. Similarly, the second structural factor of replacing X-CN with X-OMe gives the same trend. Therefore, the net inhibition score cannot be derived from both parameters.
<table><caption>**Table 3.** The DFT global reactivity descriptors of the studied inhibitors.</caption><thead><tr><th>Molecule</th><th>E<sub>HOMO</sub><br>(eV)</th><th>E<sub>LUMO</sub><br>(eV)</th><th>ΔE<br>(eV)</th><th>DM.<br>(D)</th><th>MV<br>(cm<sup>3</sup>/mol)</th><th>η<br>(eV)</th><th>χ<br>(eV)</th><th>ω<br>(eV)</th><th>ΔE<sub>b-d</sub><br>(eV)</th><th>ΔN<br>(e)</th></tr></thead><tbody><tr><td>S-CN</td><td>-6.28</td><td>-1.85</td><td>4.43</td><td>12.27</td><td>214</td><td>2.22</td><td>-4.07</td><td>3.73</td><td>-0.56</td><td>2.06</td></tr><tr><td>SO-CN</td><td>-6.39</td><td>-2.00</td><td>4.39</td><td>14.41</td><td>206</td><td>2.19</td><td>-4.19</td><td>4.01</td><td>-0.55</td><td>2.11</td></tr><tr><td>SO<sub>2</sub>-CN</td><td>-6.48</td><td>-2.29</td><td>4.19</td><td>12.89</td><td>195</td><td>2.10</td><td>-4.39</td><td>4.59</td><td>-0.53</td><td>2.25</td></tr><tr><td>S-OMe</td><td>-6.19</td><td>-1.75</td><td>4.44</td><td>6.09</td><td>186</td><td>2.22</td><td>-3.97</td><td>3.54</td><td>-0.56</td><td>2.04</td></tr><tr><td>SO-OMe</td><td>-6.35</td><td>-1.92</td><td>4.43</td><td>9.84</td><td>210</td><td>2.21</td><td>-4.13</td><td>3.86</td><td>-0.55</td><td>2.08</td></tr><tr><td>SO<sub>2</sub>-OMe</td><td>-6.43</td><td>-2.14</td><td>4.29</td><td>10.47</td><td>161</td><td>2.15</td><td>-4.28</td><td>4.28</td><td>-0.54</td><td>2.17</td></tr><tr><td>Dapsone</td><td>-5.93</td><td>-1.00</td><td>4.93</td><td>6.38</td><td>175</td><td>2.47</td><td>-3.47</td><td>2.44</td><td>-0.62</td><td>1.73</td></tr></tbody></table>

Since the hardness of a molecule is a function of the energy gap ΔE, we discuss the changes in the ΔE values. ΔE is a function of reactivity of the inhibitor molecule towards the adsorption on a metallic surface. With a decrease in ΔE, the reactivity of the molecule increases, which leads to improved inhibition efficiency of the molecule. According to Table 3, the value of ΔE steadily decreases when the S atom is replaced by SO and SO₂ either for X-CN or X-OMe molecules. Moreover, in general, X-CN molecules are more reactive and less rigid than X-OMe. On the other hand, the reactivity of dapsone is lower at least by 0.5 eV than those of the rest of the molecules due to the symmetrical structure.

The dipole moment of molecules can reflect their ability to protect the metal surface [35–37]. The dipole moment of X-CN molecules is higher than that of X-OMe molecules. It was reported that a molecule with a lower dipole moment favors its accumulation on the surface layer, leading to higher inhibition efficiency [38].

Molecular volume (MV) illustrates the possibility of the inhibitor covering the metal surface. A molecule that has a considerable MV value has the highest protection to the metal surface. There is no specific order when replacing the S atom with the SO or SO₂ group or when replacing the cyano group with the methoxy one, but the volume of the molecule depends on the degree of twisting in the molecule. S-CN is the highest in MV and SO₂-OMe is the lowest one.

ΔN is usually used to indicate the amount and direction of electron flow within molecular systems. Positive values of ΔN for all studied molecules imply the electron donation from the inhibitor to the metal surface atoms. The value of ΔN is comparable (~2 electrons) for all the tested compounds except for dapsone, which is 1.73 e. Furthermore, the back-donation can be evaluated by the value of ΔE<sub>b-d</sub>. The negative sign of ΔE<sub>b-d</sub> indicates that back-donation to the inhibitor is energetically favorable. All molecules have comparable values ~-0.55 eV. The donation and back-donation processes strengthen the adhesion of the inhibitor molecules to the iron surface.

## 3.2. Monte Carlo (MC) Simulations

The studied molecules contain two different aryl moieties: mono-substituted and disubstituted aryl moieties. Since the studied inhibitors have bent molecular geometry, both moieties cannot be adsorbed together on the surface. Figure 3 illustrates the most stable configurations for their adsorption on the Fe(110) surface. The findings indicate that the disubstituted aryl moiety of diaryl sulfides (S-CN and S-OMe) is horizontally loaded onto the iron surface, while the mono-substituted aryl one is almost perpendicular to it. This means that replacing the electron-withdrawing group (CN) with an electron-donating one (OMe) does not affect the adsorption mode. The preference for sticking the disubstituted moiety to the surface is attributed to its inclusion of more donor atoms than the mono-substituted one.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Side views of the most stable configurations of the adsorbed inhibitors on Fe(110) surface in acidic medium (1 M HCl).</figcaption></figure>

The mode of adsorption was changed after replacing the S atom of the diaryl sulfide with the SO group, as the mono-substituted aryl moiety with the SO moiety was loaded onto the surface instead of disubstituted aryl one. Additionally, it was observed that the oxygen atom of the SO fragment carries the highest negative charge (~-0.94). Thus, one can expect how such molecules adsorb on the metallic surface; first, the SO fragment attacks the surface, then the most coplanar moiety attaches to the surface. A similar situation can be expected in SO₂ molecules, because the oxygen atoms carry the highest negative atomic charges. For example, introducing an OMe group instead of the CN one in the SO-CN molecule does not affect adsorption, as the mono-substituted moiety and SO fragment adsorbed together on the iron surface. This is attributed to the smaller twisting angle between the mono-substituted aryl moiety and the SO fragment by 10° compared to the disubstituted counterpart for both molecules.

An unusual mode of adsorption is observed in the case of SO₂-based molecules, as it depends on the nature of the R-group, whether it is an electron withdrawer or an electron donor. For example, in the case of the SO₂-CN molecule, the dihedral angles between the mono and disubstituted aryl moieties with the SO₂ fragment are 180 and 145 degrees,
respectively. Thus, the more planar mono-substituted aryl ring has a higher tendency to adsorb on the surface, as indicated from the MC simulation (Figure 3). On the other hand, the two moieties of the SO₂–OMe molecule tilt by approximately 5 degrees with the SO₂ fragment, i.e., the exact extent of coplanarity. Therefore, the second factor (atomic charge) becomes the playmaker of the adsorption process. Accordingly, it is found that the disubstituted aryl ring adheres to the surface, because it possesses more than one donor atom ready to interact with the surface.

Based on the above, one can conclude that two factors influence the adsorption mode: the geometry of the inhibitor and its atomic charges. This finding could help the rational design of new corrosion inhibitors. In the case of dapsone, the adsorption mode is similar to the other sulfone molecules.

The MC adsorption parameters are listed in Table 4. The meaning of these parameters is discussed in our previous paper [37]. At first glance, all the adsorption energies of the studied systems are negative, indicating the automatic adsorption process. Additionally, the high negative adsorption energy suggests strong bonding between the adsorbates and the substrate [37]. The replacement of the S atom of the diaryl sulfide by SO or SO₂ group impacts the adsorption energy. For example, replacing the sulfur atom of the S–CN molecule with an SO₂ group increased the adsorption energy by 24 kcal/mol, while replacing it with an SO₂ group decreased it by 19 kcal/mol. Additionally, replacing the S atom of the S–OMe molecule with SO or SO₂ group decreased/increased the adsorption energy by 7 and 6 kcal/mol, respectively. These changes are attributed to the selection of the inhibitor to the mode of adsorption on the Fe surface, as discussed above. Interestingly, the adsorption energy of dapsone is higher than all other studied compounds except for SO–CN. Therefore, one can expect that the studied sulfides, sulfoxides (except SO–CN), and sulfones would have inhibition efficiency exceed 91% of dapsone [15].

<table><caption>**Table 4.** The output adsorption descriptors (kcal/mol) for inhibitor/Fe(110) systems in acidic medium.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Total Energy</th><th rowspan="2">Adsorption Energy</th><th rowspan="2">Rigid Adsorption Energy</th><th rowspan="2">Deformation Energy</th><th colspan="3">dE<sub>ad</sub>/dN<sub>i</sub></th></tr><tr><th>Inh</th><th>H<sub>2</sub>O</th><th>HCl</th></tr></thead><tbody><tr><td>S–CN</td><td>-3791</td><td>-3826</td><td>-4023</td><td>197</td><td>-148</td><td>-11</td><td>-9</td></tr><tr><td>SO–CN</td><td>-3760</td><td>-3802</td><td>-3993</td><td>191</td><td>-135</td><td>-12</td><td>-8</td></tr><tr><td>SO<sub>2</sub>–CN</td><td>-3799</td><td>-3845</td><td>-4040</td><td>195</td><td>-162</td><td>-14</td><td>-7</td></tr><tr><td>S–OMe</td><td>-3823</td><td>-3851</td><td>-4046</td><td>195</td><td>-146</td><td>-12</td><td>-9</td></tr><tr><td>SO–OMe</td><td>-3825</td><td>-3858</td><td>-4054</td><td>196</td><td>-125</td><td>-13</td><td>-9</td></tr><tr><td>SO<sub>2</sub>–OMe</td><td>-3808</td><td>-3845</td><td>-4040</td><td>196</td><td>-151</td><td>-12</td><td>-7</td></tr><tr><td>Dapsone</td><td>-3795</td><td>-3816</td><td>-4007</td><td>191</td><td>-137</td><td>-11</td><td>-7</td></tr></tbody></table>

On the other hand, the effect of insertion of the methoxy group instead of the cyano group is clearer, where the energy required for adsorption decreases in this direction. In such a case, the methoxy group increases the electron density on the phenyl ring and the sulfur atom, facilitating the interaction with the surface. Therefore, the inhibitors can be arranged according to the strength of their interaction with the surface as follows: SO–OMe > S–OMe > SO₂–OMe = SO₂–CN > S–CN > dapsone > SO–CN.

By comparing the average adsorption energies of the inhibitors (–3835 kcal/mol) with their average desorption energies (dE<sub>ad</sub>/dN<sub>i</sub>) (–143 kcal/mol), it was found that the difference between them is significant, which indicates that the adsorption process is irreversible and heavily favored.

## 4. Conclusions

The present work aimed to study the possibility of utilizing six non-nucleoside reverse transcriptase inhibitors of HIV-1 (diaryl sulfide derivatives) as anticorrosive materials for steel. Since dapsone has the same molecular backbone (diphenyl sulfide) as all studied molecules, and its inhibition efficiency has been tested experimentally, it was used here as a reference compound to evaluate the efficiency of the investigated compounds. In this
regard, the effect of two structural factors on the inhibition ability was studied: (i) the effect of replacing the sulfide group of diaryl sulfide with sulfoxide and sulfonyl ones, and (ii) the effect of R-substituent by replacing the cyano group with methoxy. Two computational approaches were employed to achieve this goal: DFT and MC simulation. Based on the findings, one can conclude the following points:

* The DFT geometry of the studied molecules is not flat and has a bent molecular shape, which leads to incomplete molecular adsorption on the surface.

* Two factors control the adsorption mode on the iron surface: the extent of coplanarity of aryl moiety with the X-group and the value of individual negative atomic charges.

* The dipole moment of the studied molecules correlates well with their adsorption ability. A molecule with a lower dipole moment has better adsorption on the iron surface.

* Although the X–CN molecules are more reactive and less rigid than X–OMe molecules, the latter have a stronger ability for adsorption because of the high electron-donating ability of the methoxy group.

* Based on the adsorption study, all the studied compounds, except for SO–CN, show higher inhibition efficiency than dapsone. Accordingly, we nominate these inhibitors as effective iron surface corrosion inhibitors.

**Author Contributions:** Conceptualization, M.M.E.-H.; methodology, M.M.E.-H., A.M.K.; software, M.M.E.-H., A.M.K.; validation, M.M.E.-H., A.M.K., M.A.A.; formal analysis, M.M.E.-H.; investigation, M.M.E.-H.; resources, M.M.E.-H., M.A.A.; data curation, M.M.E.-H.; writing—original draft preparation, M.M.E.-H.; writing—review and editing, M.M.E.-H., J.R., M.A.A.; visualization, M.M.E.-H., J.R., M.A.A.; supervision, M.A.A., M.M.A.M., R.B. and J.R.; project administration, M.A.A.; funding acquisition, M.A.A. All authors have read and agreed to the published version of the manuscript.

**Funding:** The authors are thankful to the Taif University Researchers Supporting Project number (TURSP-2020/03), Taif University, Taif, Saudi Arabia.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** All of the research data are available from the authors.

**Acknowledgments:** The computers used for quantum chemical calculations were allocated by the New Valley University, Egypt.

**Conflicts of Interest:** Authors declare no conflicts of interests.

## References

1. Kaya, S.; Banerjee, P.; Saha, S.K.; Tüzün, B.; Kaya, C. Theoretical evaluation of some benzotriazole and phospono derivatives as aluminum corrosion inhibitors: DFT and molecular dynamics simulation approaches. *RSC Adv.* **2016**, *6*, 74550–74559. [CrossRef]

2. Kabanda, M.M.; Obot, I.B.; Ebenso, E.E. Computational study of some amino acid derivatives as potential corrosion inhibitors for different metal surfaces and in different media. *Int. J. Electrochem. Sci.* **2013**, *8*, 10839–10850.

3. Kandemirli, S.G.; Yilmazer, M.I.; Saracoglu, M.; Kandemirli, F. Theoretical B3LYP study of contrast agent metrizoate. *Int. J. Chem. Technol.* **2020**, *4*, 60–70. [CrossRef]

4. Geerlings, P.; Proft, F.D.; Langenaeker, W. Conceptual Density functional theory. *Chem. Rev.* **2003**, *103*, 1793–1874. [CrossRef] [PubMed]

5. Ebenso, E.E.; Arslan, T.; Kandemirli, F.; Caner, N.; Love, I. Quantum chemical studies of some rhodanine azosulpha drugs as corrosion inhibitors for mild steel in acidic medium. *Int. J. Quantum Chem.* **2010**, *110*, 1003–1018. [CrossRef]

6. Danaee, I.; RameshKumar, S.; RashvandAvei, M.; Vijayan, M. Electrochemical and quantum chemical studies on corrosion inhibition performance of 2,2'-(2-hydroxyethylimino)bis[N-(alphaalpha-dimethylphenethyl)-N-methylacetamide] on mild steel corrosion in 1 M HCl Solution. *Mater. Res.* **2020**, *23*, e20180610. [CrossRef]

7. Sachin, H.P.; Moinuddin Khan, M.H.; Raghavendra, S.; Bhujangaiah, N.S. L-Dopa as corrosion inhibitor for mild steel in mineral acid medium. *Open Electrochem. J.* **2009**, *1*, 15–18. [CrossRef]

8. Dohare, P.; Chauhan, D.; Sorour, A.; Quraishi, M. DFT and experimental studies on the inhibition potentials of expired tramadol drug on mild steel corrosion in hydrochloric acid. *Mater. Discov.* **2017**, *9*, 30–41. [CrossRef]

9. Al-Fahemia, J.H.; Abdallaha, M.; Elshafie, Gad, A.M.; Jahdaly, B. Experimental and theoretical approach studies for melatonin drug as safely corrosion inhibitors for carbon steel using DFT. *J. Mol. Liq.* **2016**, *222*, 1157–1163. [CrossRef]
10. EL-Haddad, M.N.; Fouda, A.S.; Hassan, A.F. Data from chemical, electrochemical and quantum chemical studies for interaction between cephalpirin drug as an eco-friendly corrosion inhibitor and carbon steel surface in acidic medium. *Chem. Data Collect.* **2019**, *22*, 100251. [CrossRef]

11. Abdelaal, M.S.; Abdel-Wahab, A.A.; Assaf, F.H. Inhibition of corrosion by zinn und cadmium $H_2SO_4$ losses during verschiedene organische sulfone, sulfoxide und sulfide. *Metalloberflaeche* **1980**, *34*, 323–327.

12. Trabanelli, G.; Zucchi, F.; Gullini, G.; Carassiti, V. Correlation of the Structure and the Inhibitive Action of Some Sulphoxides. *Br. Corros. J.* **1969**, *4*, 212–215. [CrossRef]

13. James, J.P.; Lalgudi, R.S. Corrosion Inhibiting Coating Additive. U.S. Patent Application No. 16/347,247, 22 August 2019.

14. Monzo, J.; Garcia-Anton, J.; Guinon, J.L. Study of corrosion on copper strips by mixtures of mercaptans, sulphides and disulphides with elemental sulphur in the ASTM D-130 test by means of electron microscopy (SEM) and energy dispersive X-ray (EDX). *J. Anal. Chem.* **1992**, *343*, 593–596. [CrossRef]

15. Singh, A.; Singh, A.K.; Quraishi, M.A. Dapsone: A novel corrosion inhibitor for mild steel in acid media. *Open Corros. J.* **2010**, *2*, 43–51.

16. Singh, P.; Chauhan, D.S.; Chauhan, S.S.; Singh, G.; Quraishi, M.A. Chemically modified expired Dapsone drug as environmentally benign corrosion inhibitor for mild steel in sulphuric acid useful for industrial pickling process. *J. Mol. Liq.* **2019**, *286*, 110903. [CrossRef]

17. Zhang, R.; Ding, H.; Pu, X.; Qian, Z.; Xiao, Y. Recent advances in the synthesis of sulfides, sulfoxides and sulfones via C-S bond construction from non-halide substrates. *Catalysts* **2020**, *10*, 1339. [CrossRef]

18. Chan, J.H.; Hong, J.S.; Hunter, R.N., III; Orr, G.F.; Cowan, J.R.; Sherman, D.B.; Sparks, S.M.; Reitter, B.E.; Andrews, C.W., III; Hazen, R.J.; et al. 2-Amino-6-arylsulfonylbenzonitriles as non-nucleoside reverse transcriptase inhibitors of HIV-1. *J. Med. Chem.* **2001**, *44*, 1866–1882. [CrossRef]

19. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. (Eds.) *Gaussian 16, Revision C.01*; Gaussian, Inc.: Wallingford, CT, USA, 2016.

20. Tomasi, J.; Mennucci, B.; Cammi, R. Quantum mechanical continuum solvation models. *Chem. Rev.* **2005**, *105*, 2999–3093. [CrossRef]

21. Pearson, R.G. Hard and soft acids and bases—the evolution of a chemical concept, Coord. *Chem. Rev.* **1990**, *100*, 403–425. [CrossRef]

22. Dewar, M.J.; Thiel, W. Ground states of molecules. 38. The MNDO method. Approximations and parameters. *J. Am. Chem. Soc.* **1977**, *99*, 4899–4907. [CrossRef]

23. Pauling, L. *The Nature of the Chemical Bond*; Cornell University Press: New York, NY, USA, 1960.

24. Parr, R.G.; Pearson, R.G. Absolute hardness: Companion parameter to absolute electronegativity. *J. Am. Chem. Soc.* **1983**, *105*, 7512–7516. [CrossRef]

25. Chattaraj, P.K.; Sarkar, U.; Roy, D.R. Electrophilicity index. *Chem. Rev.* **2006**, *106*, 2065–2091. [CrossRef]

26. Kokalj, A. On the HSAB based estimate of charge transfer between adsorbates and metal surfaces. *Chem. Phys.* **2012**, *393*, 1–12. [CrossRef]

27. Rubinstein, R.Y.; Ridder, A.; Vaisman, R. *Fast Sequential Monte Carlo Methods for Counting and Optimization*; John Wiley & Sons: New York, NY, USA, 2014.

28. Frenkel, D.; Smit, B. *Understanding Molecular Simulation: From Algorithms to Applications*, 2nd ed.; Academic Press: San Diego, CA, USA, 2002.

29. Dassault Systèmes. *Dassault Systèmes BIOVIA, Materials Studio*, 17.1.0.48; Dassault Systèmes: San Diego, CA, USA, 2017.

30. Sun, H.; Ren, P.; Fried, J.R. The COMPASS force field: Parameterization and validation for phosphazenes. *Comput. Theor. Polym. Sci.* **1998**, *8*, 229–246. [CrossRef]

31. Guo, L.; Qi, C.; Zheng, X.; Zhang, R.; Shen, X.; Kaya, S. Toward understanding the adsorption mechanism of large size organic corrosion inhibitors on an Fe(110) surface using the DFTB method. *RSC Adv.* **2017**, *7*, 29042–29050. [CrossRef]

32. Fouda, A.S.; Ismail, M.A.; EL-ewady, G.Y.; Abousalem, A.S. Evaluation of 4-amidinophenyl-2,2'-bithiophene and its aza-analogue as novel corrosion inhibitors for CS in acidic media: Experimental and theoretical study. *J. Mol. Liq.* **2017**, *240*, 372–388. [CrossRef]

33. Foresman, J.B.; Frisch, E. *Exploring Chemistry with Electronic Structure Methods*, 3rd ed.; Gaussian, Inc.: Wallingford, CT, USA, 2015.

34. Pearson, R.G. Absolute electronegativity and hardness: Application to inorganic chemistry. *Inorg. Chem.* **1988**, *27*, 734–740. [CrossRef]

35. Lagrené, M.; Mernari, B.; Chaibi, N.; Traisnel, M.; Vezin, H.; Bentiss, F. Investigation of the inhibitive effect of substituted oxadiazoles on the corrosion of mild steel in HCl medium. *Corros. Sci.* **2001**, *43*, 951–962. [CrossRef]

36. Quraishi, M.A.; Sardar, R. Hector bases—A new class of heterocyclic corrosion inhibitors for mild steel in acid solutions. *J. Appl. Electrochem.* **2003**, *33*, 1163–1168. [CrossRef]

37. Ibrahim, M.M.; Saleh, D.I.; El-Hendawy, M.M.; Fallatah, A.M.; Mersal, G.A.M.; Boukherroub, R.; Wysocka, J.; Ryl, J.; Amin, M.A. Efficacious alkaline copper corrosion inhibition by a mixed ligand copper(II) complex of 2,2'-bipyridine and glycine: Electrochemical and theoretical studies. *ChemElectroChem* **2021**, *8*, 2052–2064. [CrossRef]

38. Khalil, N. Quantum chemical approach of corrosion inhibition. *Electrochim. Acta* **2003**, *48*, 2635–2640. [CrossRef]