Advanced
Experimental and Theoretical Studies on Corrosion Inhibition Performance of Phenanthroline for Cast Iron in Acid Solution
Experimental and Theoretical Studies on Corrosion Inhibition Performance of Phenanthroline for Cast Iron in Acid Solution
Journal of Electrochemical Science and Technology. 2018. Dec, 9(4): 260-275
Copyright © 2018, The Korean Electrochemical Society
  • Received : May 12, 2018
  • Accepted : July 16, 2018
  • Published : December 31, 2018
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
About the Authors
B. Idir
Research center in industrial technologies (CRTI), P.O. Box 64, Cheraga 16014 Algiers, Algeria
idirbrahim2016@gmail.com
F. Kellou-Kerkouche
USTHB, Laboratoire d’Electrochimie-Corrosion, Métallurgie et Chimie Minérale Faculté de Chimie BP32 El-Alia Bab-Ezzouar Alger-Algérie

Abstract
The corrosion inhibition of cast iron in 1 M HCl by Phenanthroline (Phen) was investigated using potentiodynamic polarization (PDP) curves, electrochemical impedance spectroscopy (EIS), surface analysis and theoretical calculations. It is found that Phen exhibits high inhibition activity towards the corrosive action of HCl and its adsorption obeys the Langmuir adsorption isotherm model. The results showed that inhibition efficiency increases with Phen concentration up to a maximum value of 96% at 1.4 mM, and decreases slightly with the increase in temperature. The free adsorption energy value indicates that Phen adsorbs on cast iron surface in 1 M HCl via a simultaneous physisorption and chemisorption mechanism. Scanning electron microscopy (SEM) micrographs, atomic force microscopy (AFM) and FTIR analysis confirmed the formation of a protective film on cast iron surface, resulting in the improvement of its corrosion resistance in the studied aggressive solution. Quantum chemical calculations at the DFT level were achieved to correlate electronic structure parameters of Phen molecules with their adsorption mode.
Keywords
1. Introduction
Ductile cast iron is specifically useful in pipelines, agriculture, some vehicle components and wind turbine parts [1 - 3] . It has been recognized as the most appropriate material for a wide variety of industrial applications due to its good machinability, high castability, and other various advantages [4 - 6] . With its unique microstructure, this material provides manufacturers an exceptional combination of toughness and reliability at a considerably lower cost [7] .
Acidic solutions are widely employed in the industry for acid pickling and acid cleaning; they are also added to oil and gas production. It should be emphasized that hydrochloric acid (HCl) is the most used for this purpose [8 , 9] . Unfortunately, because of the presence of Cl ions this acid is highly corrosive. The use of organic inhibitors is the most practical and effective method for the protection of ferrous materials [10 - 17] . Generally, molecules containing nitrogen, sulfur, oxygen and aromatic rings are preferred [18 - 20] . The efficiency of these molecules depends mainly on their ability to adhere to the metal surface [21] .
Phenanthroline is a classic chelating bidentate ligand for transition metal ions, which has played an important role in the development of coordination chemistry; the presence of nitrogen and multiple bonds with high electron density in its flat molecular structure facilitates its ability to inhibit corrosion of most metals. Banerjee [22] reported the corrosion inhibition of Phenanthroline on mild steel in aerated 0.5 M H 2 SO 4 , whereas its adsorption on the aluminum surface was investigated by Xianghong Li in 1 M HCl solution [23] . Obi-egbedi et al have tested a synthetized Phenanthroline derivative as a corrosion inhibitor for mild steel in 0.5 M H 2 SO 4 solution [24] . Furthermore, Xia Liu et al. [25] have studied the effect of Phenanthroline and its cobalt complex on the corrosion of mild steel in 0.1 M H 2 SO 4 solution. However, an overview of the literature revealed that few studies were devoted to the corrosion inhibition of various categories of cast iron in acidic medium. One can cite the influence of sodium dodecyl benzene sulfonate on the corrosion behavior of an elaborated cast iron alloy in sulfuric acid [26] . Recently, many studies [27 - 30] have reported that some expired drugs, carbohydrates, Schiff base derivatives and some plant extracts are effective corrosion inhibitors of cast iron corrosion in various environments.
In this context, the current work describes the inhibitive performance of Phenanthroline on cast iron corrosion in aerated hydrochloric acid, using electrochemical measurements. The inhibited metallic surface was examined by scanning electron microscopy, atomic force microscopy, and FTIR spectroscopy. The donor-acceptor relationship between inhibitor molecules and cast iron surface was achieved by the calculated theoretical parameters of the neutral Phen molecule and its protonated specie, in aqueous phase, using quantum chemical calculations.
2. Experimental
- 2.1 Material
The chemical composition (in wt.%) of the commercial studied cast iron is: 3.70% C, 1.98% Si, 0.25% Mn, 0.019% P, 0.017% S, 0.093% Cr, 0.01% Al, 0.021% Ni, 0.007% Cu, 0.059% Mg, 0.01% Ti and balance Fe. A round rod with an exposed area of 1.76 cm 2 was assembled as the working electrode. Before introducing it in the electrochemical cell, the exposed face of the electrode was ground mechanically with silicon carbide papers of graded grit ranging from 600 to 1200 grit, then rinsed with bidistilled water and degreased with acetone.
- 2.2 Corrosive and inhibitive solutions
The corrosive 1 M HCl solution was prepared using analytical reagent HCl in double distilled water. Phenanthroline ( Fig. 1 ) was added to the acid solution at different concentrations. Phenanthroline and its derivatives have been employed as a ligand in coordination chemistry; forming stable complexes with most metal ions [31] .
PPT Slide
Lager Image
Chemical molecular structure of Phenanthroline
- 2.3 Electrochemical methods
Electrochemical experiments were performed in a conventional three-electrode cell. The cast iron sample was used as the working electrode (WE); a saturated Ag/AgCl electrode and a platinum foil (surface area of 1.0 cm 2 ) were used as reference and counter electrode, respectively. Electrochemical impedance spectroscopy measurements were achieved at the open circuit potential for the frequency range of 100 KHz to 10 mHz, using a peak-to-peak voltage excitation of 10 mV. Potentiodynamic polarization curves were plotted from −300 to +300 mV versus open circuit potential ( E cor ) at the scan rate of 1 mV/s.
All electrochemical measurements were carried out using a Potentiostat-Galvanostat (Princeton Applied Research Versastat 3). Corr-View and Z -view software’s were used for Tafel extrapolations and the fit of EIS diagrams, respectively. Before performing EIS diagrams and potentiodynamic polarization curves, the sample has been immersed in 1 M HCl solution with and without addition of inhibitor for 30 minutes until a steady state open circuit potential ( E cor ) is reached.
- 2.4 Quantum chemical calculations
Theoretical calculations were performed using density functional theory (DFT) in Spartan 08 program package, at B3LYP/6-31 G (d, p) level of theory. Optimization of the molecules structure was accomplished in a water phase [32] . It is believed that this level optimizes accurately and provides right electronic properties for a wide range of organic molecules. Some parameters such as the energy of the highest occupied molecular orbital ( E HOMO ), the energy of the lowest unoccupied molecular orbital ( E LUMO ) and dipole moments ( μ ) were obtained. Absolute electronegativity ( χ ), hardness ( η ), global electrophilicity ( ω ) and global nucleophilicity ( ε ) were calculated from the following equations:
PPT Slide
Lager Image
PPT Slide
Lager Image
Where ionization potential I = − E HOMO and electron affinity A = − E LUMO .
Global electrophilicity index ( ω ) is as follow [33] :
PPT Slide
Lager Image
Nucleophilicity ( ε ) is the inverse of electrophilicity (1/ ω ).
To get more insight about ability of inhibitor molecules to interact with iron atoms, and to predict the most favorable adsorbed inhibitor, Pearson [34] has shown that the fraction of electrons transferred (Δ N ) between the inhibitor molecules and the iron atom could be calculated by the following equation:
PPT Slide
Lager Image
Where χ Fe , χ inh , η Fe , η inh are electronegativities and hardness values of iron and inhibitor molecule, respectively. Recently, Kokalj [35] has shown that the use of work function ( Φ ) of the iron surface instead of χ Fe is a more appropriate value of iron electronegativity. In this study, we have considered Fe(110) because it is the most stable plane and it presents a dense surface package, hence derived theoretical values of iron (110) have been considered: Φ Fe (110) = 4.82 eV and η Fe = 0 eV [35] .
The Fukui indices are calculated by taking the finite difference approximations from Mulliken population’s analysis (MPA) of atoms for Phen and Phen-H + as follows:
PPT Slide
Lager Image
PPT Slide
Lager Image
Where qk ( N ), qk ( N + 1) and qk ( N − 1) are the atomic charges of the systems with N , N +1, and N −1 electrons, respectively.
- 2.5 Surface characterization
The cast iron specimens were in first polished with a series of emery paper to obtain a surface mirror, rinsed with double distilled water and thereafter degreased with acetone. The cleaned specimens were immersed in 1 M HCl solution with and without the optimum concentration of Phen (1.4 mM). After 04 hours of exposure to the tested solutions, the cast iron specimens were taken out, cleaned with double distilled water and finally dried with a cold air blower. Surface characterization was carried out using a ZEISS EVO MA 25 scanning electron microscope at an accelerating voltage of 20 kV.
The topology measurements (on the samples immersed for 03 hours in the above testing solutions) were performed using a Bruker Dimension Icon AFM atomic force microscope. Before starting the tests, the samples were cleaned and dried thoroughly with a cold air blower to eliminate all corrosion products.
Infrared spectra of the pure Phen and the scraped layer collected from specimen immersed during 6 hours in HCl solution containing 1.4 mM of inhibitor (Phen) were recorded via FTIR spectrometer (Agilent Cary 630) measured in the range of 4000-600 cm −1 .
3. Results and Discussion
- 3.1 Electrochemical measurements
- 3.1.1 Polarization studies at 25°C
The inhibitive properties of Phenanthroline were analyzed by potentiodynamic polarization curves and electrochemical impedance spectroscopy diagrams. The polarization curves of cast iron in 1 M HCl solution in the absence and presence of various concentrations of Phen are illustrated in Fig. 2 . Different electrochemical kinetic parameters were determined by extrapolating the linear Tafel regions to the corrosion potential ( Table 1 ): the corrosion current density ( i cor ), the corrosion potential ( E cor ), cathodic Tafel slope ( b c ), anodic Tafel slope ( b a ), as well as the inhibition efficiency ( IE %) which is defined as follows:
PPT Slide
Lager Image
Where: i cor and i cor(inh) are the uninhibited and inhibited corrosion current densities, respectively.
As can be seen in Fig. 2 , anodic and cathodic current densities decrease with increasing inhibitor concentration. This is an indication that the organic compound under study reduces the hydrogen evolution reaction as well as the anodic cast iron dissolution. The parallel Tafel lines obtained from cathodic branches indicate that hydrogen evolution reaction is activation controlled and the addition of Phen does not modify the mechanism of this process, however, the proton discharge process is delayed. The anodic curves of cast iron in 1 M HCl in the presence of organic compound show that the tested molecule has markedly hindered the anodic reaction in the potential range from E cor to −200 mV, however, at the highest anodic potentials the curves ideally overlapped; this phenomenon can be attributed to the desorption of the inhibitive molecules from iron surface [36] .
PPT Slide
Lager Image
Potentiodynamic polarization curves of cast iron in 1 M HCl containing Phen at different concentrations: 0.36, 0.76, 1 and 1.40 mM at 25°C with a scan rate of 1 mV/s.
The data in Table 1 showed that corrosion current densities ( i cor ) are significantly decreased with increasing Phen concentration, which means that the iron sample has become more resistant to corrosion. The values of the anodic and cathodic Tafel slopes, b a and b c , respectively, have nearly the same order of magnitude in the free and inhibitive solutions, which indicates that the corrosion mechanism of the cast iron is unaffected by the addition of Phen. On the other hand, the corrosion potential ( E cor ) has not undergone any significant variation, meaning that the organic inhibitor is adsorbed on the anodic and cathodic sites of the metal surface. So, it can be recognized that Phen acts as a mixed-type inhibitor [37] . The value of inhibition efficiency ( IE %) increases up to 94% with increasing the concentration of Phen at the optimum value of 1.4 mM.
Tafel and EIS parameters of cast iron in 1 M HCl solution with and without Phen at 25°C
PPT Slide
Lager Image
Tafel and EIS parameters of cast iron in 1 M HCl solution with and without Phen at 25°C
- 3.1.2 EIS studies at 25°C
The impedance diagrams ( Fig. 3a ) represented in a complex plane, describe the behavior of cast iron in a solution of HCl without and with inhibitor at different concentrations. It can be seen that all EIS diagrams display a single depressed capacitive semicircle. These imperfect semicircles may be attributed to the surface roughness, the frequency dispersion, as well as the heterogeneity of the cast iron surface [38] . Nyquist plots of the blank and inhibited solutions show similar behavior, indicating that the cast iron corrosion reaction is mainly controlled by charge transfer, and the addition of the inhibitor does not modify the mechanism of the process [39] . The semicircles diameters in the presence of Phen are larger than that observed in its absence; this phenomenon could be attributed to the adsorption of the investigated inhibitor on the cast iron surface by blocking its active sites, leading to the improvement of cast iron protection against corrosion [40] .
The impedance plots are fitted accurately with a good chi-squared coefficient ( X 2 ) using the Randles equivalent circuit ( R s ( R ct / CPE ) having a single time constant, as shown in the insert of Fig. 3b . In this circuit, R s is the solution resistance, R ct is the charge transfer resistance related to the OCP corrosion reaction, while the CPE represents a constant phase element related to the non-ideal capacity ( C dl ). A similar circuit was proposed in several works [41 - 44] . The constant phase element ( CPE ) is introduced to simulate the non-ideal capacitive behavior of cast iron/HCl interface, and its impedance is defined as follows [43] :
PPT Slide
Lager Image
Where Y 0 is a CPE constant, j is the imaginary number, ω is the angular frequency ( ω = 2π f , f represents the AC frequency in Hz), n is a phase shift, which is related to the system homogeneity, when the CPE display a pure capacitance, n = 1. The electrical double layer capacitance ( C dl ) was determined by the following equation [45] :
PPT Slide
Lager Image
The EIS parameters of the cast iron immersed in uninhibited and inhibited solutions are summarized in Table 1 . The results obtained are consistent with those deduced from the Tafel method. Indeed, the charge transfer resistance increases with the increase in Phen concentration, while the capacitance decreases, indicating that this inhibitor has reduced the aggressiveness of the acid solution.
PPT Slide
Lager Image
a) Nyquist plots and b) Bode plots of cast iron corrosion in 1 M HCl without and with different concentrations of Phen: 0.36, 0.76, 1 and 1.40 mM at 25°C.
The inhibition efficiency ( IE %) was calculated, using the charge transfer resistance values ( R ct ), by the following equation:
PPT Slide
Lager Image
Where R ct(inh) and R ct represent the charge transfer resistance derived from the inhibited and uninhibited solutions, respectively.
The highest value of R ct bound to lowest C dl value corresponds to the highest inhibition efficiency recorded at the optimum concentration of Phen (1.4 mM); which means the effective adsorption of Phen which causes the slowing of the cast iron corrosion process.
According to the Helmholtz model, the decrease of the C dl value (Eq.11) indicates an increase in the double layer thickness ( d ), which can be attributed to the development of a compact protective film on the metal surface by the inhibitor adsorption [46] .
PPT Slide
Lager Image
Where d is the thickness of the film, ε ° is the permittivity of the air, ε is the medium dielectric constant, and S represents the surface area of the working electrode.
Bode and phase angle plots are given in Fig. 3b , they can provide complementary information on the complicated corrosion process. The impedance modulus, at low frequencies, increases with the increase of the inhibitor concentration, indicating that the adsorption of the inhibitive molecules improves the corrosion resistance of the studied cast iron in the corrosive HCl solution [47] . It is also evident, that the mid-frequency phase angle values have evolved in the negative direction with increasing inhibitor concentration, indicating the improvement of Phen adsorption on the metal surface. In addition, there is only one peak in the phase angle plots that confirms the existence of a single time constant at the metal/solution interface.
The results obtained from the electrochemical measurements are in reasonable agreement and confirm once more that Phenanthroline shows good inhibition performance (above 95%) for cast iron in hydrochloric acid solution. A comparison of corrosion inhibition efficiency of phenanthroline with various phenanthroline derivatives reported as corrosion inhibitors is presented in Table 2 [22 - 25 , 48 - 52] .
Comparison of corrosion inhibition efficiency of phenanthroline with other phenanthroline derivatives reported as corrosion inhibitors in the literature.
PPT Slide
Lager Image
Comparison of corrosion inhibition efficiency of phenanthroline with other phenanthroline derivatives reported as corrosion inhibitors in the literature.
- 3.2 Adsorption isotherm
The electrochemical measurements evidenced the Phen adsorption on the cast iron surface and therefore the inhibition of general corrosion of the metal surface, by forming a barrier for the hydrated corrosive ions. Consequently, it would be very useful to describe the adsorption process by an appropriate adsorption isotherm. The relative coverage area ( θ ) is proportional to the inhibition efficiency ( IE %), which is determined directly as follow [47] :
PPT Slide
Lager Image
Adsorption isotherms are usually used to describe the adsorption process of inhibitors. Their determination can provide important clues about the nature of metal-inhibitor interactions [53] . The surface coverage ( θ ) values obtained from the Tafel and EIS methods, corresponding to different Phen concentrations at 25°C, were used to deduce the best isotherm. Several attempts have been made to identify among all the adsorption isotherms; the best fit was obtained with the Langmuir adsorption isotherm modeled by the following equation [54] :
PPT Slide
Lager Image
The plot of C inh / θ vs C inh is characterized by a straight line ( Fig. 4 ).
PPT Slide
Lager Image
Adsorption isotherms of inhibitor on cast iron surface in 1 M HCl.
The slope of the plot was close to 1, indicating that Phen molecules interact significantly with cast iron surface to form an inhibiting film that corresponds to a single layer [55] . The hypothesis of this isotherm implies that all the adsorption sites are equivalent and that the adsorption capacity of the molecules on a given site is independent for the occupation of neighboring sites [56] .
In addition, with using this model (Langmuir isotherm), the Gibbs free energy (∆ G o ads ), which can be used to describe the stability of the adsorption bond between the compound and the metal, was determined using K ads in the following equation:
PPT Slide
Lager Image
It is generally accepted that the values of Δ Go ads revolving around −20 kJmol −1 , indicate that the adsorption is considered as a physisorption, the inhibition occurs via electrostatic interactions between the charged molecules and the charged metal [57] , while the values around −40 kJ.mol −1 or smaller (more negative) are seen as chemisorption, which is due to the charge sharing or electrons transfer from the inhibitor molecules to the metal surface [58] . The values of Δ Go ads deduced from the potentiodynamic polarization curves and the EIS measurement are −35.53 kJ mol −1 and −36.38 kJ mol −1 , respectively. A large negative value of Δ Go ads indicates the greater adsorption capacity of Phen molecules on cast iron surface. In addition, calculated ∆ Go ads values exceed electrostatic interactions range and are close to the chemical interactions region. This confirms the mixed mode of adsorption where physisorption and chemisorption contribute to the whole adsorption process [59] .
- 3.3 Temperature effect and thermodynamic parameters
Temperature has a great influence on the corrosion process and generally, it increases the corrosion rate [60] . For this purpose, EIS measurements were performed in the temperature range of 25-55°C for the cast iron sample immersed in 1 M HCl solution, without and with the optimal Phen concentration ( Fig. 5a - b ).
Examination of Nyquist plots in Fig. 5a - b reveals that the rising in temperature does not change the profile of diagrams, which means that the corrosion process is charge transfer controlled [61] . The non-ideal circular shape of these diagrams could be attributed to the frequency dispersion and heterogeneity of the cast iron surface resulting from the interfacial phenomena [61 - 63] . The size of these imperfect semi-circles decreases with increasing temperature; this is probably attributed to the decrease of cast iron corrosion resistance resulting from the partial desorption of the inhibitive molecules [63] . The equivalent electrical circuit used to fit the EIS data is shown in the inset of Fig. 5a - b . The deduced results are summarized in Table 3 .
PPT Slide
Lager Image
Nyquist plots of cast iron in 1 M HCl (a) in absence and (b) in presence of 1.4 Mm of Phen at various temperatures (between 25 and 55°C).
It can be seen from table 3 that the temperature increase reduced the value of the charge transfer resistance R ct in the presence or absence of inhibitor and simultaneously the double layer capacitance value increased. However, it should be noted that the inhibition efficiency remains almost constant, which provides evidence for the existence of chemisorption because this adsorption type has a great adsorption heat and it is usually irreversible [64] . It follows from the above that Phen is an effective corrosion inhibitor for cast iron immersed in the acidic solution since even at 55°C its inhibition efficiency is about 95%.
Temperature effect on the EIS parameters for the cast iron immersed in 1 M HCl with and without the optimal concentration of Phen
PPT Slide
Lager Image
Temperature effect on the EIS parameters for the cast iron immersed in 1 M HCl with and without the optimal concentration of Phen
It would be interesting to investigate the corrosion kinetics of cast iron in HCl solution in the presence and absence of inhibitor, using the Arrhenius equations as follows:
PPT Slide
Lager Image
PPT Slide
Lager Image
Where A (A cm −2 ) represents a pre-exponential factor, E a (J mol −1 ) is the apparent activation energy, R (8.314 J mol −1 K −1 ) is the gas constant and T (K) is the absolute temperature, N is the Avogadro’s number, h is the plank’s constant,
PPT Slide
Lager Image
and
PPT Slide
Lager Image
are the enthalpy and entropy of activation, and i cor (A cm −2 ) represents the corrosion current densities computed from the charge transfer resistance R ct at various temperatures in the presence and absence of Phen using the following equation:
PPT Slide
Lager Image
Herein; R ct represents the charge transfer resistance, z is the valence of iron ( z =2), F is the Faraday constant (F = 96485 C), R and T retain the earlier meaning in the previous equation.
The plots ln ( i cor ) versus
PPT Slide
Lager Image
( Fig. 6a ) corresponding to the blank and inhibitive solutions are linear throughout the temperature range, which indicates that the corrosion and corrosion inhibition processes follow the kinetics Arrhenius. The apparent activation energy values deduced from the slopes of the two straight lines are: 30.50 kJ mol −1 and 40.78 kJ mol −1 in the blank and inhibitive solution respectively. It has been reported [65 - 67] that inhibitors whose inhibition efficiency decreases with increasing temperature have a higher apparent activation energy value than that of the blank solution, due to physical adsorption of inhibitors on the metal surface. In the present study, the increase of the activation energy value in the presence of Phen indicates that the adsorption of Phen molecules occurs primarily through physical interactions.
The representative plots of ln (
PPT Slide
Lager Image
) versus
PPT Slide
Lager Image
corresponding to the acidic and inhibitive solutions are linears ( Fig. 6b ), with a great coefficient of correlation. The slope
PPT Slide
Lager Image
gives
PPT Slide
Lager Image
and the intercept ln
PPT Slide
Lager Image
gives
PPT Slide
Lager Image
. The calculated values are gathered in Table 4 .
PPT Slide
Lager Image
a) Arrhenius plots and b) transition state plots for cast iron in 1 M HCl (blank) and in the solution containing 1.4 mM of Phen.
Activation parameters for the corrosion of cast iron in 1 M HCl in the presence and absence of Phen
PPT Slide
Lager Image
Activation parameters for the corrosion of cast iron in 1 M HCl in the presence and absence of Phen
The positive values of
PPT Slide
Lager Image
reflected the endothermic nature of the cast iron dissolution process. The value of
PPT Slide
Lager Image
for the uninhibited solution is lower than that for the inhibited solution. This phenomenon means a decrease in ramdomness has occurred going from reactants to the activated complex. This could be the result of the adsorption of the organic inhibitor on the electrode surface. The great negative value of entropy in the inhibited solution means that the disorder on cast iron surface decreased in the presence of Phen.
- 3.4 Theoretical study
- 3.4.1 Molecular orbital theory and quantum molecular descriptors
Theoretical calculations can provide detailed information for understanding the reactive behavior of corrosion inhibiting molecules. Density Functional Theory (DFT) has become a very useful tool for analyzing experimental data to elucidate the corrosion inhibition mechanism [68] . The optimized molecular structure and the frontier orbitals density distributions (HOMO and LUMO) of neutral Phen and Phen-H + molecules in aqueous phase are shown in Fig. 7 . It can be seen that the two chemical structures are planar which favors their efficient adsorption leading to high coverage of cast iron surface. The highest occupied molecular orbital (HOMO) is related to the ability of a molecule to donate electrons to empty metal orbitals, while the lowest unoccupied molecular orbital (LUMO) represents the capacity of the molecule to accept electrons from full metal orbitals ( Fig. 7 ). It can be seen that the electron densities HOMO and LUMO of Phen and Phen-H + are distributed in the whole area of the molecules, thus, the most likely adsorption mode is in parallel on cast iron surface [69] .
PPT Slide
Lager Image
The optimized chemical structure, topology of HOMO and LUMO orbitals and Molecular Electrostatic Potential (MEP) for neutral and protonated Phen.
Inspection of determined quantum parameters, listed in Table 5 , reveals that after Phen protonation, E HOMO becomes more negative, suggesting that protonation reduces the ability to donate electrons to unoccupied metal orbitals [70] . It is well known that the higher the HOMO energy of the inhibitor, the higher its ability to donate electrons and the higher its corrosion inhibition efficiency. In addition, the E LUMO value of Phen-H + shifts to the negative value compared to the E LUMO value of Phen neutral, suggesting that the protonated molecule has more ability to accept electrons from iron [71] . On the other hand, the energy gap (Δ E ) of the protonated molecule is slightly lower than that of the neutral molecule; this suggests that it has more ability to react than the neutral molecule. Recently several authors introduced the electrophilicity index ( ω ), which describes the electrophilic capacity of a molecule [33] . In the current study, Phen-H + has the highest value of ω , which implies the protonated form’s greater ability to accept electrons from the iron atoms. According to the literature [71 , 72] , a molecule with a low electrophilicity index is a good corrosion inhibitor, for this reason, it can be stated that the neutral inhibitor is more effective against the corrosion of iron in acid solution.
Quantum chemical parameters of Phen and Phen-H+
PPT Slide
Lager Image
Quantum chemical parameters of Phen and Phen-H+
The organic inhibitor adsorption is usually considered as a substitution process between the organic molecules in the aqueous solution Org (sol) and the water molecules H 2 O (ads) adsorbed on the electrode surface [73] :
Org(sol) + xH2O(ads) → Org(ads) + xH2O(sol)
Where: Org (sol) and Org (ads) are the organic molecules in the solution and adsorbed on the metal surface, respectively. H 2 O (ads) is the water molecules adsorbed on the metal surface, and x is the size ratio representing the number of water molecules replaced by an adsorbed organic molecule.
The dipole moment ( μ ) values of Phen and PhenH+ shown in Table 4 are higher than that of water ( μ H 2 O =1.84 Debye), which implies that the displacement of the water molecules from the metal surface is favorable [74] .
According to Lukovits et al. [75] , the Δ N value indicates the donating ability of electrons from molecule to metal if Δ N > 0 and vice versa if Δ N < 0. They also suggest that if Δ N < 3.6, the inhibition efficiency increases with increasing electron donating ability. In the present study, the results in Table 5 show that Phen and its protonated form are electrons donors, which means the strong and stable adsorption of these species on the cast iron surface in hydrochloric solution. Therefore, all of these factors indicate that the inhibitive molecules could involve both electrostatic interaction and chemical bond, resulting in higher inhibition efficiency.
- 3.4.2 Local reactivity properties
Molecular electrostatic potential (MEP) illustrated in Fig.7 is a visual method that allows us to distinguish the location of the electron density. In these maps, different colors are observed, where the red represents the zone with the negative potential of MEP, associated with reactive electrophilic sites, the blue color is adapted to the zone with the positive potential and represents the suitable center of the nucleophilic attacks [76] . The red colored region in the neutral molecule is concentrated mainly on the hetero-atoms, while the blue region is located around the phenanthrene. However, in the case of Phen-H + , the positive charge is obviously concentrated on the protonated nitrogen and concerning the negative charge, it is localized on the non-protonated nitrogen atom and also on the adjacent carbon atoms. This presentation separates electron density and electrostatic potential surfaces with clarity and compactness. The main disadvantage is that it provides information only on the contact surface and does not reveal how far electron-rich and electron-poor areas extend beyond the surface.
In addition to the MEP surfaces, the Fukui indices are more accurate indicators of the molecule parts disposed to electrophilic or nucleophilic attacks, they describe changes in atomic electron density as a result of adding or removing a charge.
The Fukui indices
PPT Slide
Lager Image
and
PPT Slide
Lager Image
shown in Table 6 , were calculated to predict the most likely atomic sites for nucleophilic and electrophilic attacks respectively. They can also explain the charge transfer as they indicate the reactive sites in the molecule. Atoms that have a high absolute
PPT Slide
Lager Image
are more susceptible to nucleophilic attacks, while higher absolute
PPT Slide
Lager Image
are more prone to attack by species that have electrons deficiencies.
The values of Fukui indices for Phen and Phen-H+considering Mulliken population analysis (MPA)
PPT Slide
Lager Image
The values of Fukui indices for Phen and Phen-H+ considering Mulliken population analysis (MPA)
For nucleophilic attack, the most reactive sites are N10 and N14 for the neutral Phenanthroline and C11 and C13 for the protonated Phenanthroline, which can accept electrons from cast iron surface to form back donating bond. On the other hand, the highest
PPT Slide
Lager Image
coincide on C5 and C6 for both the neutral and protonated Phen, which may denote electrons transfer to the metal to form a coordinate bond. According to these findings, it can be suspected a likely donoracceptor interaction engaged by inhibitor molecules and iron surface. Overall, the chemical reactivity of Phen and Phen-H + can be regarded as π-excessive, thus the adsorption could be achieved simply by the geometric blocking effect. The theoretical study demonstrates the ability of Phen molecules to offer electrons to unoccupied d - orbital of iron, and also to accept electrons from the iron surface, which is consistent with the obtained adsorption isotherm.
- 3.5 Surface analysis
Scanning electron microscope (SEM) was used in order to demonstrate that the corrosion inhibition is mainly due to the formation of an adsorptive layer on the metal surface. Fig. 8a shows the surface of the freshly polished cast iron alloy, which consists of iron matrix embedded graphite nodules; this surface is very smooth, except some defects which had arisen during polishing.
After immersion in the uninhibited solution ( Fig. 8b ), the cast iron surface is strongly damaged due to the uniform dissolution of the matrix. In addition, graphite nodules peeled off from the surface caused by galvanic corrosion leading to wide and deep holes, as shown in the magnified image ( Fig. 8b ’). However, in the presence of Phen at its optimal concentration ( Fig. 8c ) the cast iron surface is smooth, intact and well improved. These observations confirm that Phenanthroline forms a protective layer by adsorption, which blocks the matrix dissolution, leading to the lowering of cast iron corrosion rate.
PPT Slide
Lager Image
SEM micrographs of cast iron samples a) As polished, b) immersed in 1 M HCl during four hours without inhibitor, b’) Magnified image of cast iron sample immersed in the uninhibited solution and c) in the presence of 1.4 mM of Phen.
The 3D AFM image of the freshly polished sample ( Fig. 9a ) shows a regular and smooth surface, his average roughness (Ra) is 17 nm. After immersing the sample in the blank solution ( Fig. 9b ), the average roughness increases to 275 nm , as a result of severe damage of the surface. However, in the presence of Phen ( Fig. 9c ) the average roughness decreases to 97 nm, due to the formation of a protective film by adsorption of Phen molecules on the metal surface.
PPT Slide
Lager Image
3D AFM micrographs of cast iron samples a) as polished, b) immersed in 1 M HCl during three hours in the absence of inhibitor and c) in the presence of 1.4 mM of Phen.
The FTIR spectra for pure Phen and scrapped adsorbed layer on the cast iron ( Fig. 10 ) have been recorded to determine if Phen effectively adsorbs on the metal surface. The FTIR spectrum of pure phen shows a broad water band at 3365 cm −1 attributed to O-H stretching. The absorption band at 3050 cm −1 is related to the aromatic C-H stretching mode, and the strong absorption bands around 1585 cm −1 , 1558 cm −1 are assigned to C=N stretching vibrations. The peaks at 1502 and 1418 cm −1 are associated with the carbocyclic ring vibrations (C=C), the presence of band at 1134 cm −1 corresponds probably to the in-plane hydrogen deformation modes or ring vibrations, and other peaks at 850 and 736 cm −1 correspond to the hydrogen (C-H) out of plane bending on the benzene and the heterocyclic rings. A full vibrational assignment of phenanthroline molecule has been previously published [77] . The most bands observed in the FTIR spectra of the adsorbed layer on cast iron surface closely resemble those appearing in the pure Phen spectra. A broad band around 3300 cm −1 attributed to O-H stretching, which further indicates the formation of FeOOH on the adsorbed layer. the aromatic C-H stretching band has vanished from 3050 to 3029 cm −1 . The shifted to higher wavenumbers of the C=N and C=C bands confirm the involvements of these centers in the adsorption. in addition, the absorption bands at 850 and 738 cm −1 assigned to the out of plane C-H vibrations are shifted to 835 and 715 cm −1 , respectively. Reference has been made to the fact that these bands (out of plane vibrations) move to lower wavenumbers on coordination [78] . All the informations indicate clearly towards adsorption of Phen molecules on the cast iron surface, similar conclusions have been made by many other researchers [79 - 80] .
PPT Slide
Lager Image
(a, a’) FTIR spectra of pure Phen and adsorbed layer formed on cast iron surface after 6 hours’ immersion in 1 M HCl solution with 1.4 mM of Phen.
- 3.6 Inhibition mechanism
The adsorption process of organic molecules on a metal surface depends on their chemical structure, the distribution charge in the molecule, the aggressive medium and the topography of the metal surface. Therefore, on the basis of the experimental results, FTIR in particular which indicates the presence of most functional groups of Phen on the adsorbed layer, a corrosion inhibition mechanism could be proposed:
In HCl solution, Phenanthroline may exist in the protonated form in equilibrium with its molecular form, thus, two modes of adsorption could be envisaged. According to the literature, generally ferrous alloys carry an excess positive charge in acid solution, thus, the acidic anions (Cl ) specifically adsorbed on the metal surface create an excess negative charge which favors the adsorption of Phen-H + , by involving electrostatic forces.
In addition to the physical adsorption, the neutral Phen may be adsorbed on the metal surface through chemical adsorption by sharing electrons on the basis of donor-acceptor between the π-electrons of the phenanthroline ring and the vacant d-orbital of iron surface atoms. The molecular structure of Phen-H + remains almost unchanged with respect to its neutral form, it can be stated that when Phen-H + is adsorbed on the metal surface, the coordination bond can also be formed by partial transference of electrons from the polar Nitrogen atom to the metal surface.
4. Conclusions
The corrosion inhibition effect of Phenanthroline on the corrosion of ductile iron in 1 M HCl was investigated using electrochemical and theoretical calculations. On the basis of the results of these studies, the following conclusions can be drawn:
· Analysis of the polarization curves showed that Phenanthroline acts as a mixed-type inhibitor; its inhibition efficiency increases with the concentration up to a limit value obtained at the optimum concentration.
· The adsorption of Phen molecules on the metal surface obeys the Langmuir adsorption isotherm model.
· The high negative adsorption free energy value and the inhibition efficiency maintained at high temperatures indicate that Phen adsorption on the cast iron surface could be of mixed type where both physisorption and chemisorption contribute to the whole process.
· SEM and AFM micrographs have shown that the cast iron surface damage is greatly reduced in the presence of Phenanthroline which forms a protective layer. The Phen adsorption on the cast iron surface was confirmed by FTIR analysis.
· The Fukui indices analysis of inhibitor molecules has successfully described the nucleophilic and electrophilic sites responsible for the donation and acceptance of electrons.
· In summary, it is concluded that the results obtained from the analysis of chemical quantum data corroborate those achieved by electrochemical measurements.
References
Olofsson J. , Svensson I. 2013 Mater Des 43 264 - 271    DOI : 10.1016/j.matdes.2012.07.006
Górny M. , Tyrała E. 2013 J Mater Eng Perform 22 (1) 300 - 305    DOI : 10.1007/s11665-012-0233-0
Nilsson K. F. , Blagoeva D. , Moretto P. 2006 Eng Fract Mech 73 (9) 1133 - 1157    DOI : 10.1016/j.engfracmech.2005.12.005
Geier G.F. , Bauer W. , McKay B.J. , Schumacher P. 2005 Mater Sci Eng A 413 339 - 345
Kim S. , Cockcroft S.L. , Omran A.M. 2009 J Alloys Compd 476 (1-2) 728 - 732    DOI : 10.1016/j.jallcom.2008.09.082
Wessén M. , Svensson I.L. , Aagaard R. 2003 Int J Cast Met Res 16 (1-3) 119 - 124    DOI : 10.1080/13640461.2003.11819569
Putatunda S.K. 2010 Mater Manuf Process 25 (8) 749 - 757    DOI : 10.1080/10426910903367394
El-Askalany A.H. , Mostafa S.I. , Shalabi K. , Eid A.M. , Shaaban S. 2016 J Mol Liq 223 497 - 508    DOI : 10.1016/j.molliq.2016.08.088
Musa A.Y. , Mohamad A.B. , Kadhum A.A.H. , Takriff M.S. , Tien L.T. 2011 Corros Sci 53 (11) 3672 - 3677    DOI : 10.1016/j.corsci.2011.07.010
Noor E. A. , Al-Moubaraki A. H. 2008 Mater Chem Phys 110 145 - 154    DOI : 10.1016/j.matchemphys.2008.01.028
Slimane M. , Kellou F. , Kellou-Tairi S. 2015 Res Chem Intermed 41 (11) 8571 - 8590    DOI : 10.1007/s11164-014-1912-2
Mello L. D. , Gonçalves R.S. 2001 Corros Sci 43 (3) 457 - 470
Abdeli M. , Ahmadi N. P. , Khosroshahi R. A. 2010 J Solid State Electrochem 14 (7) 1317 - 1324    DOI : 10.1007/s10008-009-0925-z
Abdallah M. , Al-Tass H. M. , Jahdaly B.A.A.L. , Fouda A. S. 2016 J Mol Liq 216 590 - 597    DOI : 10.1016/j.molliq.2016.01.077
Morales-Gil P. , Walczak M.S. , Camargo C. R. , Cottis R. A. , Romero J.M. , Lindsay R. 2015 Corros Sci 101 47 - 55    DOI : 10.1016/j.corsci.2015.08.032
Yan Y. , Li W. , Cai L. , Hou B. 2008 Electrochem Acta 53 (20) 5953 - 5960    DOI : 10.1016/j.electacta.2008.03.065
Finsgar M. , Jackson J. 2014 Corros Sci 86 17 - 41    DOI : 10.1016/j.corsci.2014.04.044
Shabani-Nooshabadi M. , Behpour M. , Razavi F. S. , Hamadanian M. , Nejadshafiee V. 2015 RSC Adv 5 (30) 23357 - 23366    DOI : 10.1039/C5RA00561B
Bai L. , Feng L. J. , Wang H. Y. , Lu Y. B. , Lei X. W. , Bai F. L. 2015 RSC Adv 5 (6) 4716 - 4726    DOI : 10.1039/C4RA12286K
Abboud Y. , Abourriche A. , Saffaj T. , Berrada M. , Charrouf M. , Bennamara A. , Al Himidi N. , Hannache H. 2007 Mater Chem Phys 105 (1) 1 - 5    DOI : 10.1016/j.matchemphys.2007.03.037
Quraishi M. A. , Rawat J. , Ajmal M. 2000 J Appl Electrochem 30 745 - 751    DOI : 10.1023/A:1004099412974
Banerjee S. N. , Misra S. 1989 Corrosion 45 (9) 780 - 783    DOI : 10.5006/1.3585034
Li X. , Deng S. , Xie X. 2014 J Taiwan Inst Chem Eng 45 (4) 1865 - 1875    DOI : 10.1016/j.jtice.2013.10.007
Obi-Egbedi N. O. , Obot I. B. , Eseola A. O. 2014 Arab J Chem 7 197 - 207    DOI : 10.1016/j.arabjc.2010.10.025
Liu X. , Okafor P. C. , Jiang B. , Hu H. , Zheng Y. 2015 J Mater Eng Perform 24 (9) 3599 - 3606    DOI : 10.1007/s11665-015-1608-9
Kellou-Kerkouche F. , Benchettara A. , Amara S. 2008 Mater Chem Phys 110 (1) 26 - 33    DOI : 10.1016/j.matchemphys.2008.01.005
Rajeswari V. , Viswanathamurthi P. , Kesavan D. 2017 Res Chem Intermed 43 (7) 3893 - 3913    DOI : 10.1007/s11164-016-2852-9
Rajeswari V. , Kesavan D. , Gopiraman M. , Viswanathamurthi P. 2013 Carbohydr Polym 95 288 - 294    DOI : 10.1016/j.carbpol.2013.02.069
Rajeswari V. , Kesavan D. , Gopiraman M. , Viswanathamurthi P. 2013 J Surfact Deterg 16 (4) 571 - 580    DOI : 10.1007/s11743-013-1439-3
Rajeswari V. , Kesavan D. , Gopiraman M. , Viswanathamurthi P. , Poonkuzhali K. , Palvannan T. 2014 Appl Surf Sci 314 537 - 545    DOI : 10.1016/j.apsusc.2014.07.017
Chikira M. , Tomizawa Y. , Fukita D. , Sugizaki T. , Sugawara N. , Yamazaki T. , Sasano A. , Shindo H. , Palaniandavar M. , Antholine W. 2002 J Inorg Biochem 89 163 - 173    DOI : 10.1016/S0162-0134(02)00378-1
Li X. , Deng S. , Lin T. , Xie X. , Du G. 2017 Corros Sci 118 202 - 216    DOI : 10.1016/j.corsci.2017.02.011
Parr R. G. , Szentpaly L. , Liu S. 1999 J Am Chem Soc 121 (9) 1922 - 1924    DOI : 10.1021/ja983494x
Pearson R. G. 1988 Inorg. Chem. 27 (4) 734 - 740    DOI : 10.1021/ic00277a030
Kokalj A. 2001 Electrochim Acta 56 745 - 755
Qu Q. , Hao Z. , Li L. , Bai W. , Liu Y. , Ding Z. 2009 Corros Sci 51 (3) 569 - 574    DOI : 10.1016/j.corsci.2008.12.010
Ostovari A. , Hoseinieh S. M. , Peikari M. , Shadizadeh S. R. , Hashemi S. J. 2009 Corros Sci 51 (9) 1935 - 1949    DOI : 10.1016/j.corsci.2009.05.024
Singh A. K. , Quraishi M. A. 2010 Corros Sci 52 (4) 1373 - 1385    DOI : 10.1016/j.corsci.2010.01.007
Oguzie E. E. , Li Y. , Wang F. H. 2007 J Colloid Interf Sci 310 (1) 90 - 98    DOI : 10.1016/j.jcis.2007.01.038
Bommersbach P. , Dumont C. A. , Millet J. P. , Normand B. 2006 Electrochim Acta 51 (19) 4011 - 4018    DOI : 10.1016/j.electacta.2005.11.020
Abdel-Gaber A. M. , Abd-El-Nabey B. A. , Sidahmed I. M. , El-Zayady A. M. , Saadawy M. 2006 Corros Sci 48 (9) 2765 - 2779    DOI : 10.1016/j.corsci.2005.09.017
Li X. , Deng S. , Fu H. 2012 Corros Sci 62 163 - 175    DOI : 10.1016/j.corsci.2012.05.008
Xu B. , Yang W. , Liu Y. , Yin X. , Gong W. , Chen Y. 2014 Corros Sci 78 260 - 268    DOI : 10.1016/j.corsci.2013.10.007
Khan G. , Basirun W. J. , Kazi S. N. , Ahmed P. , Magaji L. , Ahmed S. M. , Khan G. M. , Rehman M. A. 2017 J Colloid Interf Sci 502 134 - 145    DOI : 10.1016/j.jcis.2017.04.061
Hsu C H. , Mansfeld F. 2001 Corrosion 57 (9) 747 - 748    DOI : 10.5006/1.3280607
Khaled K. F. 2008 Electrochim Acta 53 (9) 3484 - 3492    DOI : 10.1016/j.electacta.2007.12.030
Abd- El Rehim S. S. , Ibrahim M. A. M. , Khalid K.F. 2001 Mater Chem Phys 70 (3) 268 - 273    DOI : 10.1016/S0254-0584(00)00462-4
Li X. , Tang L. , Li L. , Mu G. , Liu G. 2006 Corros Sci 48 (2) 308 - 321    DOI : 10.1016/j.corsci.2004.11.029
Mu G.N , Li X , Li F 2004 Mater Chem Phys 86 (1) 59 - 68    DOI : 10.1016/j.matchemphys.2004.01.041
Obot I. B. , Obi-Egbedi N.O , Ebenso E.E , Afolabi A.S. , Oguzie E.E. 2013 Res. Chem. Intermed 39 (5) 1927 - 1948    DOI : 10.1007/s11164-012-0726-3
Obot I.B. , Obi-Egbedi N.O. , Eseola AO. 2011 Ind. Eng. Chem. Res 50 (4) 2098 - 2110    DOI : 10.1021/ie102034c
Lei X , Wang H , Feng Y , Zhang J , Sun X , Lai S , Wang Z , Kang S 2015 RSC. Adv 5 (120) 99084 - 99094    DOI : 10.1039/C5RA15002G
Ashassi-Sorkhabi H. , Seifzadeh D. 2008 J. Appl. Electrochem 38 (11) 1545 - 1552    DOI : 10.1007/s10800-008-9602-7
Yadav M. , Kumar S. , Gope L. 2014 J. Adhes. Sci. Technol 28 (11) 1072 - 1089    DOI : 10.1080/01694243.2014.884754
Badiea A. M. , Mohana K. N. 2009 Corros. Sci 51 (9) 2231 - 2241    DOI : 10.1016/j.corsci.2009.06.011
Solomon M. M. , Umoren S. A. , Udosoro I. I. , Udoh A.P. 2010 Corros. Sci 52 (4) 1317 - 1325    DOI : 10.1016/j.corsci.2009.11.041
El-Awady A. A. , Abd-El-Nabey B. A. , Aziz S. G. 1992 J. Electrochem. Soc 139 (8) 2149 - 2154    DOI : 10.1149/1.2221193
Obot I. B. , Ebenso E. E. , Obi-Egbedi N. O. , Afolabi A. S. , Gasem Z. M. 2012 Res. Chem. Intermed 38 (8) 1761 - 1779    DOI : 10.1007/s11164-012-0501-5
Abdeli M. , Ahmadi N. P. , Khosroshahi R. A. 2011 J Solid State Electrochem 15 (9) 1867 - 1873    DOI : 10.1007/s10008-010-1162-1
Nesic S. 2007 Corros. Sci 49 (12) 4308 - 4338    DOI : 10.1016/j.corsci.2007.06.006
Goulart C. M. , Esteves-Souza A. , Martinez-Huitle C.A. , Ferreira C. J. , Medeiros M. A. , Echevarria A. 2013 Corros. Sci 67 281 - 291    DOI : 10.1016/j.corsci.2012.10.029
Lebrini M. , Lagrenée M. , Vezin H. , Traisnel M. , Bentiss F. 2007 Corros Sci 49 (5) 2254 - 2269    DOI : 10.1016/j.corsci.2006.10.029
Behpour M. , Ghoreishi S. M. , Soltani N. , Salavati-Niasari M. , Hamadanian M. , Gandomi A. 2008 Corros. Sci 50 (8) 2172 - 2181    DOI : 10.1016/j.corsci.2008.06.020
Li X. , Deng S. , Fu H. 2012 Corros. Sci 55 280 - 288    DOI : 10.1016/j.corsci.2011.10.025
Szauer T. , Brand A. 1981 Electrochim. Acta 26 (9) 1253 - 1256    DOI : 10.1016/0013-4686(81)85107-9
Hamani H. , Douadi T. , Daoud D. , Al-Noaimi M. , Rikkouh R. , Chafaa S. 2017 J. electroanal. Chem 801 425 - 438    DOI : 10.1016/j.jelechem.2017.08.031
Riggs OL , Hurd RM 1967 Corrosion 23 (8) 252 - 260    DOI : 10.5006/0010-9312-23.8.252
Gece G. 2008 Corros Sci 50 (11) 2981 - 2992    DOI : 10.1016/j.corsci.2008.08.043
Martinez S. , Stagljar I. 2003 J Mol Struc Theochem 640 167 - 174    DOI : 10.1016/j.theochem.2003.08.126
Lukovits I. , Kalman E. , Zucchi F. 2001 Corrosion 57 (1) 3 - 8    DOI : 10.5006/1.3290328
Al-Sabagh A. M. , Nasser N. M. , Farag A. A. , Migahed M. A. , Eissa A. M. F , Mahmoud T. 2013 Egypt. J. Pet 22 (1) 101 - 116    DOI : 10.1016/j.ejpe.2012.09.004
El-Din M. R. N. , Khamis E. A. 2014 J. Surfact. Deterg 17 (4) 795 - 805    DOI : 10.1007/s11743-014-1565-6
Obot I. B. , Kaya S. , Kaya C. , Tuzun B. 2016 Res. Chem. Intermed 42 (5) 4963 - 4983    DOI : 10.1007/s11164-015-2339-0
Hegazy M. A. 2015 J. Mol. Liq 208 227 - 236    DOI : 10.1016/j.molliq.2015.04.042
Abd El-Lateef H. M. , Tantawy A. H. 2016 RSC. Adv 6 (11) 8681 - 8700    DOI : 10.1039/C5RA21626E
Madkour L. , Kaya S. , Kaya C. , Guo L. 2016 J Taiwan Inst Chem E 68 461 - 480    DOI : 10.1016/j.jtice.2016.09.015
Thornton D.A. , Watkins G.M. 1991 Spectrochim. Acta A 47 (8) 1085 - 1096    DOI : 10.1016/0584-8539(91)80039-L
Schilt A.A. , Taylor R.C 1959 J. Inorg. Nucl. Chem. 9 (3-4) 211 - 221    DOI : 10.1016/0022-1902(59)80224-4
Vengatesh G. , Karthik G. , Sundaravadivelu M. 2017 Egypt. J. Petrol 26 (3) 705 - 719    DOI : 10.1016/j.ejpe.2016.10.011
Li L. , Zhang X. , Lei J. , He J. , Zhang S. , Pan F 2012 Corros Sci 63 82 - 90    DOI : 10.1016/j.corsci.2012.05.026