Translate this page into:
Thermodynamic, chemical and electrochemical investigations of 2-mercapto benzimidazole as corrosion inhibitor for mild steel in hydrochloric acid solutions
*Corresponding author at: Electrochemistry Research Laboratory, Chemistry Department, Faculty of Education, Ain Shams University, Roxy, Cairo, Egypt khaledrice2003@yahoo.com (K.F. Khaled)
-
Received: ,
Accepted: ,
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.
Available online 18 June 2010
Abstract
The inhibiting action of 2-mercapto benzimidazole (2MBI) on mild steel in 1.0 M hydrochloric acid has been investigated at 308 K using weight loss measurements and electrochemical techniques (impedance spectroscopy and potentiodynamic polarisation). Inhibition efficiency increases with 2MBI concentration to attain 98% at 10−3 M. Polarisation curves indicate that 2MBI acts as a mixed-type inhibitor. Inhibition efficiency values obtained from various methods were in good agreement. EIS measurements showed an increase of the transfer resistance with the inhibitor concentration. The temperature effect on the corrosion behaviour of steel in 1.0 M HCl without and with the 2MBI at various concentrations was studied in the temperature range from 308 to 353 K. Thermodynamic parameters such as heat of adsorption , entropy of adsorption and adsorption free energy have been calculated. Kinetic parameters for the corrosion reaction at different concentrations of 2MBI were determined. Adsorption of 2MBI on the mild steel surface in 1.0 M HCl follows the Langmuir isotherm model.
Keywords
Mild steel
2-Mercaptobenzimidazole
Acid inhibition
EIS
Adsorption
Thermodynamic parameters
1 Introduction
Hydrochloric acid is a strong inorganic acid that is used in many industrial processes. The most important areas of applications are acid pickling, acid descaling and oil well acidizing. During these processes, metal such as steel is subjected to serious acid corrosion and inhibitors are often needed to reduce corrosion rates in this media.
Most of the efficient inhibitors used in industry are organic compounds having multiple bonds in their molecules, which mainly contain nitrogen and sulphur. Nitrogen-containing compounds are more effectively in HCl (Bentiss et al., 2000; Abd El-Maksoud, 2003), whereas sulphur-containing compounds are preferred for H2SO4. These compounds can adsorb on metal surface and block the cathodic and/or anodic sites on the metal surface which lead to decrease the corrosion rate (Fisher, 1972; Riggs Jr., 1974). The mode of adsorption depends mainly on the chemical structure of the inhibitor, the chemical composition of the solution, the nature of the metal surface and the electrochemical potential of the metal-solution interface (Ita and Offiong, 1997a,b; Lukovits et al., 2001).
A survey of literature shows that nitrogen-containing organic compounds, such as amines, heterocyclic compounds (Ousslim et al., 2009; Sykes, 1984; Tang et al., 2003; Mernari et al., 1998; Bentiss et al., 2000; Walker, 1975; Xue et al., 1991; Wang et al., 2002; Kertit and Hammouti, 1996) and imidazole derivatives (Stupnicek-Lisac et al., 2002; Sivaraju et al., 2008; Khaled and Amin, 2009; Zhang et al., 2004; Christov and Popova, 2004; Khaled, 2003) offer good protection of metallic materials in acidic solutions. The encouraging results obtained with 2-mercapto benzimidazole (Wang et al., 2003; Popova et al., 2003; Trachli et al., 2002; Wang, 2001; Aljourani et al., 2009; Amar et al., 2007; Khaled and Fadl-Allah, 2009) have incited us to extend its use in the corrosion-inhibiting action on mild steel in HCl solution.
The aim of this work is to study the effect of 2-mercapto benzimidazole (2MBI) on the corrosion inhibition of mild steel in molar hydrochloric acid solution. The behaviour of steel in 1.0 M HCl with and without inhibitor is studied using gravimetric, potentiodynamic and EIS measurements. The thermodynamic parameters for the adsorption process and activation parameters for steel dissolution reactions are determined and discussed.
2 Experimental
2-Mercapto benzimidazole (2MBI) was analytical grade. The molecular structure is shown below:
Gravimetric measurements were carried out in a double walled glass cell equipped with a thermostat-cooling condenser. The solution volume was 100 ml. The steel specimens used had a rectangular form (2.5 cm × 2 cm × 0.05 cm). The immersion time for the weight loss was 6 h at 308 K and 1 h at the other temperatures. After the corrosion test, the specimens of steel were carefully washed in double-distilled water, dried and then weighed. The rinse removed loose segments of the film of the corroded samples. Triplicate experiments were performed in each case and the mean value of the weight loss is reported. Weight loss allowed us to calculate the mean corrosion rate as expressed in mg cm−2 h−1.
Electrochemical measurements were carried out in a conventional three-electrode electrolysis cylindrical Pyrex glass cell. The working electrode (WE) in the form of disc cut from steel has a geometric area of 1 cm2 and is embedded in polytetrafluoroethylene (PTFE). A saturated calomel electrode (SCE) and a disc platinum electrode were used as reference and auxiliary electrodes, respectively. The temperature was thermostatically controlled at 308 ± 1 K. The working electrode, WE was abraded with silicon carbide paper (grade P1200), degreased with acetone and rinsed with double-distilled water before use.
Running on a personal computer, the 352 Soft Corr™ III software communicates with EG&G Instruments potentiostat–galvanostat model 263 A at a scan rate of 0.5 mV s−1. Before recording the cathodic polarisation curves, the steel electrode is polarised at −800 mV for 10 min. For anodic curves, the potential of the electrode is swept from its corrosion potential after 30 min at free corrosion potential, to more positive values. The test solution is deaerated with pure nitrogen. Gas bubbling is maintained through the experiments.
The electrochemical impedance spectroscopy (EIS) measurements were carried out with the electrochemical system (Tacussel) which included a digital potentiostat model Voltalab PGZ 100 computer at Ecorr after immersion in solution without bubbling, the circular surface of steel exposing of 1 cm2 to the solution were used as working electrode. After the determination of steady-state current at a given potential, sine wave voltage (10 mV) peak to peak, at frequencies between 100 kHz and 10 mHz were superimposed on the rest potential. Computer programs automatically controlled the measurements performed at rest potentials after 30 min of exposure. The impedance diagrams are given in the Nyquist representation. Values of Rt and Cdl were obtained from Nyquist plots.
3 Results and discussion
3.1 Weight loss tests
Gravimetric measurements of steel were investigated in 1.0 M HCl in the absence and presence of various 2MBI concentrations at 6 h of immersion and at 308 K. The inhibition efficiency (α%) was calculated by the following relation:
Table 1 collects the corrosion rates and Fig. 1 shows the variation of inhibition efficiencies evaluated from weight loss measurements for different inhibitor concentrations in 1.0 M HCl. The corrosion rate decreases with 2MBI concentration and in turn the inhibition efficiency (α%) increases to attain 98%. From weight loss measurements, we can conclude that 2MBI is an excellent inhibitor.
Concentration (×103 M)
W (mg cm−2 h−1)
α%
0.000
1.151
–
0.001
1.113
03
0.005
0.943
18
0.010
0.661
43
0.050
0.374
67
0.100
0.203
82
0.500
0.041
96
1.000
0.026
98
Variation of inhibition efficiency (α%) with the concentration of 2MBI for mild steel in 1.0 M HCl.
3.2 Polarisation measurements
The polarisation curves of steel recorded in 1.0 M HCl in the absence and presence of 2MBI at different concentration at 308 K are presented in Fig. 2. The extrapolation method for the polarisation curves was applied and the data for corrosion potential (Ecorr), corrosion current density (Icorr), cathodic Tafel slopes (bc) and percentage inhibition efficiency (η%) are show in Table 2. The relation determines the inhibition efficiency (η%):
Potentiodynamic curves for mild steel in 1.0 M HCl in the presence of various 2MBI concentrations.
Concentration (×103 M)
Ecorr (mV)
bc (mV dec−1)
Icorr (μA cm−2)
η%
0.000
−418
188
508
–
0.001
−424
185
428
16
0.005
−420
178
286
44
0.010
−422
153
103
80
0.050
−429
145
65
87
0.100
−435
139
42
92
0.500
−434
137
23
96
1.000
−432
130
17
97
As it is show in Fig. 2 and Table 2, cathodic current–potential curves give rise to parallel Tafel lines indicating that the hydrogen evolution reaction is under activation controlled. The cathodic current density decreases with the concentration of 2MBI however, a slight effect is observed on the anodic portions. This result indicates that 2MBI is adsorbed on the metal surface on the cathodic sites and hence inhibition occurs.
These results demonstrate that the hydrogen reduction is inhibited and that the inhibition efficiency increases with inhibitor concentration to attain 97 at 10−3 M of 2MBI.
We remark that the inhibitor acts on the anodic portion and the anodic current density is reduced (Fig. 2). It seems also that the presence of the inhibitor change slightly the corrosion potential values in no definite direction. These results indicated that 2MBI acts as a mixed-type inhibitor.
3.3 Effect of 2MBI concentration on the corrosion rate of steel at various temperatures
The effect of temperature on the anti-corrosion effectiveness of 2MBI, was studied at various 2MBI concentrations in the temperature domain (313–353 K) after 1 h of immersion. The collected curves in Fig. 3 show the evolution of corrosion rate (W) with 2MBI concentration (C) at different temperatures.Variation of corrosion rate with the concentrations of 2MBI for mild steel in 1.0 M HCl at different temperatures.
Fig. 3 indicates that at a given 2MBI concentration the corrosion rate of steel increased with temperature. The increase is more pronounced at low concentrations. The results also indicate that for a given temperature, the corrosion rate of steel decreased with increasing inhibitor concentration. The values of inhibition efficiency obtained from the weight loss for different inhibitor concentrations and at various temperatures in 1.0 M HCl are given in Table 3 and Fig. 4. It is clear that inhibition efficiency increased with increase in inhibitor concentration. The maximum value of inhibition efficiency (α%) obtained for 10−3 M 2MBI is 98% at 308 K.
Temperature (K)
Concentration (×103 M)
W (mg cm−2 h−1)
α%
θ
313
0.000
1.431
–
–
0.050
0.667
53
0.53
0.100
0.367
74
0.74
0.500
0.111
92
0.92
1.000
0.046
97
0.97
323
0.000
3.134
–
–
0.050
2.035
35
0.35
0.100
1.295
59
0.59
0.500
0.392
87
0.87
1.000
0.208
93
0.93
333
0.000
6.618
–
–
0.050
5.402
18
0.18
0.100
3.202
52
0.52
0.500
0.995
85
0.85
1.000
0.521
92
0.92
343
0.000
11.37
–
–
0.050
10.801
05
0.05
0.100
5.57
51
0.51
0.500
1.672
85
0.85
1.000
0.989
91
0.91
353
0.000
20.45
–
–
0.050
20.02
02
0.02
0.100
14.34
30
0.3
0.500
3.84
81
0.81
1.000
2.513
88
0.88
Variation of inhibition efficiency (α%) with concentration of 2MBI for mild steel in 1.0 M HCl at different temperatures.
The inhibition efficiencies decrease slightly with increasing temperature indicating that higher temperature dissolution of steel predominates on adsorption of 2MBI at the metal surface.
3.4 Apparent activation energy (Ea) and pre-exponential factor (A)
In order to obtain more details on the corrosion process, activation kinetic parameters such energy (Ea), enthalpy
and entropy
evaluated from the effect of temperature using Arrhenius law (Eq. (3)) and the alternative formulation of Arrhenius equation (Eq. (4)) (Aljourani et al., 2009):
Fig. 5 shows Arrhenius plot of the corrosion rate for mild steel in 1.0 M HCl in the presence of different concentrations of 2MBI. Straight lines are obtained. The activation energies are estimated from the slopes (−Ea/R). Additionally The straight lines obtained from the transition state plot of the corrosion rate for mild steel in 1.0 M HCl in the presence of different concentrations of 2MBI (Fig. 6) permit the calculation of
and
, respectively, from slopes
and intercepts
. Ea,
and
are collected in Table 4.Arrhenius plot of the corrosion rate for mild steel in 1.0 M HCl in the presence of different concentrations of 2MBI.
Transition state plot of the corrosion rate for mild steel in 1.0 M HCl in the presence of different concentrations of 2MBI.
Concentration (×103 M)
Pre-exponential factor (kg m−2 s−1)
Linear regression coefficient
Ea (kJ mol−1)
(kJ mol−1)
(J mol−1 K−1)
Ea −
(kJ mol−1)
0.000
3.50 × 1010
0.998
59.39
56.65
−60.75
2.64
0.050
6.14 × 1013
0.997
80.96
78.23
1.41
2.73
0.010
1.10 × 1014
0.996
84.01
81.27
6.53
2.74
0.500
1.10 × 1014
0.997
87.49
84.75
6.26
2.74
1.000
1.33 × 1014
0.995
89.54
86.81
7.81
2.73
Previous studies (Aljourani et al., 2009; Elachouri et al., 1996; Ferreira et al., 2004) showed that compared with the activation energy in the absence of inhibitor, higher values for Ea were found in the presence of inhibitor. Other studies (Banerjee and Misra, 1989; Ferreira et al., 2004) showed that in the presence of 2MBI, the activation energy was lower than that in its absence.
It was clear that (see Table 4) the values of Ea in the presence of the 2MBI are higher than that in the uninhibited acid solution (59.39 kJ mol−1).
The decrease of inhibition efficiencies with increasing temperature and the increase of Ea in the presence of the inhibitor indicate the physical adsorption mechanism (Popova et al., 2003).
Fig. 7 indicates that the activation energy increases with the concentration of 2MBI. In the same way, the pre-exponential factor increases with the concentration of 2MBI (as shown in Table 4). The pre-exponential factor and activation energy varies in the same way. From Eq. (3), it can be seen that at a certain temperature, the value of the steel corrosion rate is jointly decided by the apparent activation energy and pre-exponential factor. As a whole, the steel corrosion rate basically decreases with an increase in concentration of 2MBI. As can be seen from Table 4, it was clear that both Ea and A increased in the presence of 2MBI (Mu et al., 2004).Variation of activation parameters Ea and
with 2MBI concentration for mild steel in 1.0 M HCl.
The variation of activation parameters Ea and
with 2MBI concentration is shown in Fig. 7. From the data obtained in Table 4, it seems that Ea increases with the inhibitor concentration. The positive sign of the enthalpy (
) reflects the endothermic nature of the steel dissolution process (Table 4). We remark that Ea and
values vary in the same way with the inhibitor concentration. This result permit to verify the known thermodynamic relation (Gomma and Wahdan, 1995) between Ea and
as shown also in Table 4:
The positive values of entropies ( ) in the presence of inhibitor imply that the activated complex in the rate determining step represents an association rather than a dissociation step, meaning that an increase in disordering takes place on going from reactants to the activated complex (Elachouri et al., 1996).
3.5 Adsorption isotherm and thermodynamic parameters
Assuming the corrosion inhibition was caused by the adsorption of 2MBI, and the values of surface coverage (θ) for different concentrations of 2MBI in 1 M HCl were evaluated from weight loss measurements using the Sekine and Hirakawa’s method (Sekine and Hirakawa, 1986):
In order to get a better understanding of the electrochemical process on the metal surface, adsorption characteristics are also studied for 2MBI. This process is closely related to the adsorption of the inhibitor molecules (Hackerman and Sudbury, 1950; Hackerman, 1962) and adsorption is known to depend on the chemical structure (Cheng et al., 1999; Ateya et al., 1984; Akiyama and Nobe, 1970). Adsorption isotherms are very important in determining the mechanism of organic electrochemical reactions. The most frequently used adsorption isotherms are Langmuir, Temkin and Frumkin.
In hydrochloric acid solution, the organic compound follows the Langmuir adsorption isotherm. This is as follows:
From the values of surface coverage, the linear regressions between C/θ and C are calculated and the parameters (adsorption coefficients, slopes, and linear regression coefficients) are listed in Table 5. Fig. 8 shows the relationship between C/θ and C at various temperatures. These results show that the linear regression coefficients (r) and the slopes are almost equal 1.000, indicating that the adsorption of inhibitor onto steel surface agrees the Langmuir adsorption isotherm.
Temperature (K)
Linear regression coefficient
K
Slope
(kJ mol−1)
(kJ mol−1)
(J mol−1 K−1)
308
1.0000
45537.8
0.998
−37.78
−36.26
4.93
313
0.9999
27426.9
1.007
−37.07
−36.26
2.59
323
0.9999
14316.4
1.006
−36.51
−36.26
0.77
333
1.0000
10838.5
0.994
−36.87
−36.26
1.83
343
0.9999
10830.2
1.004
−37.97
−36.26
4.99
353
0.9899
4616.8
0.899
−36.58
−36.26
0.91
Curve fitting of the corrosion data for mild steel in 1.0 M HCl in the presence of different concentrations of 2MBI to the Langmuir isotherm at different temperatures.
In addition, the equilibrium constant of the adsorption process (K) decreases with increasing temperature (Table 5). It is well known that K designates the adsorption power of inhibitor onto the steel surface; clearly, 2MBI gives higher values of K at lower temperatures, indicating that it was adsorbed strongly onto the steel surface. Thus, the inhibition efficiency decreased slightly with the increase in temperature as the result of the improvement for the desorption of 2MBI from the steel surface.
The corrosion inhibition of 2MBI for steel may be well explained by using thermodynamic model, so, the heat, the free energy and the entropy of adsorption are calculated to elucidate the phenomenon for the inhibition action of 2MBI.
According to the Van’t Hoff equation (Tang et al., 2003; Zhao and Mu, 1999):
To obtain the adsorption heat, the linear regression between lnK and 1/T is dealt with, the relationship between lnK and 1/T is shown in Fig. 9. Under the experimental conditions, the adsorption heat could be approximately regarded as the standard adsorption heat (
) (Mu et al., 2004; Zhao and Mu, 1999). The obtained value of
is −36.26 kJ mol−1. The standard adsorption free energy (
) is obtained according to the following equation (Khamis, 1990):
Variation of ln(K) with 1/T for mild steel in 1.0 M HCl in the presence of 2MBI.
The negative values of (Table 5) ensure the spontaneity of the adsorption process and stability of the adsorbed layer on the steel surface (Ali et al., 2005). Furthermore, it is found that slightly increases with temperature.
The negative values of also show that the adsorption of inhibitor is an exothermic process (Gomma and Wahdan, 1995). Generally, an exothermic process signifies either physi- or chemisorption while endothermic process is attributable unequivocally to chemisorption (Durnie et al., 1999). In an exothermic process, physisorption is distinguished from chemisorption by considering the absolute value of a physisorption process is lower than 40 kJ mol−1 while the adsorption heat of a chemisorption process approaches 100 kJ mol−1 (Martinez and Stern, 2002). In the present case; the standard adsorption heat −36.26 kJ mol−1 shows that a comprehensive adsorption (physical adsorption) might occur (Tang et al., 2003). The same results were obtained in previous studies (Trachli et al., 2002; Li and Mu, 2005). The adsorption of inhibitor molecules is accompanied by positive values of .
= −36.26 kJ mol−1 found by the Van’t Hoff equation, may be also evaluated by the Gibbs–Helmholtz equation, which is defined as follows:
The variation of
/T with 1/T gives a straight line with a slope that equals
(Fig. 10). It can be seen from Fig. 10 that
/T decreases slightly with 1/T in a linear fashion. The values of
is negative (
= −36.26 kJ mol−1), reflecting the exothermic behaviour of adsorption on the steel surface.Variation of
/T with 1/T for mild steel in 1.0 M HCl in the presence of 2MBI.
The value of the enthalpy of adsorption found by the two methods such as Van’t Hoff and Gibbs–Helmholtz relations are in good agreement.
3.6 Electrochemical impedance spectroscopy (EIS)
Fig. 11 shows a typical set of Nyquist plots for mild steel in 1.0 M HCl in the absence and presence of various concentrations of 2-mercapto benzimidazole. It is clear from these plots that the impedance response of mild steel has significantly changed after the addition of 2MBI in 1.0 M HCl. The impedance parameters derived from these plots using equivalent circuit described elsewhere (Khaled, 2003) are given in Table 6.Nyquist diagrams for steel in 1.0 M HCl in the presence of different of 2MBI concentrations.
2MBI concentration (M)
Rt (Ω cm2)
fmax (Hz)
Cdl (μF cm−2)
χ%
0
106.9
15.82
94.1
10−6
112.5
17.86
79.2
5
10−5
196.5
10.01
81.1
46
5 × 10−5
307.6
6.79
76.2
65
10−4
699.1
3.57
63.7
85
5 × 10−4
1302.1
2.23
54.8
92
10−3
2026.7
1.58
49.7
95
The charge transfer resistance, Rt values are calculated from the difference in impedance at lower and higher frequencies, as suggested by Tsuru and Haruyama (Mansfeld et al., 1981a,b; Tsuru et al., 1978). To obtain the double capacitance (Cdl), the frequency at which the imaginary component of the impedance is maximum (−Zmax) is found and Cdl values are obtained from the equation:
Table 6 indicates that by increasing the concentration of 2MBI, the Cdl values tended to decrease and the inhibition efficiency improved. This decrease in the Cdl, which can result from a decrease in local dielectric constant and/or an increase in the thickness of the electrical double layer, signifying that 2MBI molecules act by adsorption at the metal/solution interface (McCafferty and Hackerman, 1972).
The double layer capacitance Cdl is expressed in the Helmotz model by:
The decrease in Cdl values may be interpreted either by a decrease of local dielectric constant ɛ (McCafferty and Hackerman, 1972) or by the thickness of the adsorbate layer of inhibitor at the metal surface (Bastidas et al., 2000).
The data obtained from EIS are in good harmony with that obtained from potentiodynamic polarisation. The impedance spectra obtained on mild steel in 1.0 M HCl solutions consist of one depressed capacitive loop (one time constant in the Bode-phase representation, not shown here). When Nyquist plot contains a ‘‘depressed semicircle with the center under the real axis’’, such behaviour is characteristic for solid electrodes and often referred to as frequency dispersion that has been attributed to roughness and other inhomogeneities of the solid surface (Khaled, 2003).
Deviation of this kind often referred to as frequency dispersion. Therefore, a constant phase element (CPE) instead of a capacitive element is used in Fig. 2 to get a more accurate fit of experimental data sets using generally more complicated equivalent circuits. The impedance, Z, of CPE has the form (Khaled, 2003):
4 Conclusions
From the above results and discussion, the following conclusions are drawn:
-
2MBI inhibits the corrosion of steel in 1.0 M HCl. The inhibition efficiency increases with the inhibitor concentration, but decrease slightly with the temperature.
-
2MBI acts as a mixed-type inhibitor.
-
The adsorption of 2MBI on the steel surface from 1.0 M HCl obeys a Langmuir adsorption isotherm. The adsorption process is a spontaneous and exothermic process.
-
The kinetic and thermodynamic parameters of corrosion and adsorption processes are determined.
-
The results obtained from weight loss, potentiodynamic polarisation and impedance spectroscopy are in good agreement.
References
- Appl. Surf. Sci.. 2003;206:129.
- J. Electrochem. Soc.. 1970;17:999.
- Corros. Sci.. 2005;47:2659.
- Corros. Sci.. 2009;51:1836.
- Corros. Sci.. 2007;49:2936.
- Corros. Sci.. 1984;24:509.
- Corrosion. 1989;45:780.
- J. Appl. Electrochem.. 2000;30:1173.
- Corros. Sci.. 2000;42:127.
- Corros. Sci.. 1999;41:321.
- Corros. Sci.. 2004;46:1613.
- J. Electrochem. Soc.. 1999;146:1751.
- Corrosion. 1996;52:103.
- Mater. Chem. Phys.. 2004;83:129.
- Werks. Korros.. 1972;23:253.
- Mater. Chem. Phys.. 1995;39:211.
- Corrosion. 1962;18:332.
- J. Electrochem. Soc.. 1950;94:4.
- Mater. Chem. Phys.. 1997;51:203.
- Mater. Chem. Phys.. 1997;48:164.
- Appl. Surf. Sci.. 1996;93:59.
- Electrochim. Acta. 2003;48:2493.
- J. Appl. Electrochem.. 2009;39:429.
- Mater. Chem. Phys.. 2009;117:148.
- Corrosion. 1990;46:476.
- Appl. Surf. Sci.. 2005;252:1254.
- Corrosion (Houston). 2001;57(1):3.
- Corrosion. 1981;37:301.
- Corrosion. 1982;38:570.
- Appl. Surf. Sci.. 2002;199:83.
- J. Electrochem. Soc.. 1972;119:146.
- Corros. Sci.. 1998;40:391.
- Mater. Chem. Phys.. 2004;86:59.
- J. Appl. Electrochem.. 2009;39:1075.
- Corros. Sci.. 2003;45:33.
- Riggs Jr., O.L., 1974. in: Nathan, C.C. (Ed.), Corrosion Inhibitors, second ed. NACE, Houston, TX.
- Corrosion. 1986;42:272.
- Orient. J. Chem.. 2008;25:67.
- Electrochim. Acta. 2002;47:4189.
- Br. Corros. J.. 1984;25:175.
- Corros. Sci.. 2003;45:2251.
- Prog. Org. Coat.. 2002;44:17.
- J. Jpn. Soc. Corros. Eng.. 1978;27:573.
- Corros. Sci.. 1975;31:97.
- Corros. Sci.. 2001;43:2281.
- Mater. Chem. Phys.. 2002;77:655.
- Corros. Sci.. 2003;45:677.
- J. Phys. Chem.. 1991;95:7380.
- Corros. Sci.. 2004;46:3031.
- Corros. Sci.. 1999;41:1937.