# Diffusivity of Point Defects in the Passive Film on Stainless Steel

Authors: A. Fattah-alhosseini; M. H. Alemi; S. Banaei
Journal: International Journal of Electrochemistry
 (2011)

---

## Abstract

The semiconductor properties of passive films formed on AISI 316 stainless steel in sulfuric acid solution were studied by employing Mott-Schottky analysis in conjunction with the point defect model. The donor density of the passive films, which can be estimated by the Mott-Schottky plots, changes depending on the film formation potentials. Based on the Mott-Schottky analysis, an exponential relationship between donor density and the film formation potentials of the passive films was developed. The results showed that the donor densities evaluated from Mott-Schottky plots are in the range 2-3 × 1021 cm−3 and decreased with the film formation potential. By assuming that the donors are oxygen ion vacancies and/or cation interstitials, the diffusion coefficient of the donors, (DO), is calculated to be approximately 3.12 × 10−16 cm2/s.


---

## Body

## 1. Introduction

Compared with many theories qualitatively describing the passive state, the point defect model (PDM) provides a microscopic description of the growth and breakdown of a passive film and an analytical expression for the flux and the concentration of vacancies within the passive film, which affords an opportunity for quantitative analysis [1–3]. The PDM is based on the migration of point defects (oxygen and metal vacancies) under the influence of the electrostatic field in the passive film. Because the passive film is grown into the metal by the generation of oxygen vacancies at the metal/passive film interface and by their annihilation at the passive film/solution interface, the transport properties and spatial distribution of oxygen vacancies or metal cations within the passive film are of great interest [4]. A key parameter in describing the transport of point defects and hence the kinetics of passive film growth is the diffusivity of the point defects (DO). By using Mott-Schottky analysis in conjunction with the PDM for passive films, Sikora et al. [5] determined the oxygen vacancy diffusion coefficient in the WO3−x passive film formed on tungsten in phosphoric acid (pH 0.96) solution at ambient temperature to be in the range of 10−14-10−15 cm2 s-1. Although it is impossible to distinguish between oxygen vacancy and cation interstitial conducting passive films on the basis of PDM itself, the assumption is still valid in this case due to the fact that WO3 is known as priori to be oxygen vacancy. Also, by this method the diffusion coefficient of the donors in the passive film formed on AISI 316L stainless steel is determined [6].It is well recognized that passive films exhibit semiconductor properties, and numerous attempts have been made to interpret these properties in terms of classical theories for lowly and moderately doped semiconductors. However, direct measurement of passive film composition (e.g., by determining the O : M ratio through the film) and indirect electrochemical studies (e.g., by Mott-Schottky analysis) have shown that the vacancy concentrations may be very large (of the order of 1019–1021cm−3)) [7].However, there is still lack of study on the diffusivities of point defects in the passive film formed on stainless steel, which is the most extensively used engineering material in industries. In this work, the Mott-Schottky analysis of AISI 316 stainless steel (316) in sulfuric acid solution was performed, and the defects concentration was measured as a function of the film formation potential (E). The relationship between the donor density (ND) and the film formation potential is discussed in order to understand the property of the passivation of stainless steel.


## 2. Materials and Methods

Specimens were cut from 1-cm-diameter rods of 316; the nominal composition is given in the Table1. The samples were placed in stainless steel sacks and were annealed in inert environment (Ar gas) to eliminate the cold work effect due to cutting process. The annealing was performed at 1050°C for 90 min followed by water quenching. A mechanically polished 316 sample was used as a reference point (polishing under water with SiC paper up to 1200 grit size). Prior to each experiment, samples were degreased in an ultrasonic acetone/ethanol mixed bath, rinsed with distilled water, and finally dried.Table 1
Nominal chemical composition of 316.
ElementsCrNiMoMnSiCPSFe316/wt%16.512.52.221.420.420.090.040.005BalAll electrochemical measurements were performed in a conventional three-electrode cell. The counter electrode was a Pt plate, and all potentials were measured against a saturated calomel electrode (SCE) connected to the cell via a Luggin probe. All experiments were performed at ambient temperature (23–25°C) in 0.05 M H2SO4, pH 1.5.The polarization curve was obtained using a PARSTAT 2273 advanced electrochemical system with 1 mV/s scan rate. The capacitance of the interface was also recorded with a PARSTAT 2273 advanced electrochemical system. Prior to capacitance measurement, working electrodes were initially reduced potentiostatically at −1.2 V for 5 min and then polarized at a desired formation potential for 1 h to form a steady oxide film. Capacitance measurements were carried out on the anodic films at a frequency of 1 kHz using a 10 mV AC signal and a step rate of 25 mV, in the cathodic direction. An important point to be considered in the Mott-Schottky analysis is that the measured capacitance value is influenced by frequency. In order to obtain an optimized frequency for capacitance measurements, a curve of capacitance as a function of frequency at stable potential was obtained. 1 kHz was chosen because a stable capacitance value was obtained as the frequency varied from 1 to 0.001 kHz. The measurement of the capacitance-potential profile of the passive film formed at different potentials was carried out on different samples.


## 3. Results and Discussion

### 3.1. Polarization Behavior

Figure1 shows the potentiodynamic curve for 316 electrode in 0.05 M H2SO4. From the polarization curve, the passive range was determined to be from about −0.2 to 0.6 VSCE. It was observed that before the electrode was transferred to a passive state, an active current peak occurred at around −0.15 VSCE, which could be attributed to the oxidation of Fe2+ to Fe3+ions in the passive film [6].Figure 1
Potentiodynamic polarization curve of 316 in 0.05 M H2SO4 solution with 1 mV/s scan rate.Films were grown at each potential for 1 h to insure that the system was in steady state. During the formation of the oxide films on working electrode, the current density was monitored. It was observed that the current density diminishes with time until a constant value is reached and a steady-state established. Figure2 shows the values of the steady state passive current density (iSS) versus the film formation potential. The steady-state current density is approximately 6.8 μA cm−2.Figure 2
Steady-state current density obtained during the potentiostatic growth of the passive films at different film formation potentials for 1 h.


### 3.2. Mott-Schottky Measurements

The capacitance of the passive film (CT) is considered to be the combination of two series capacitances: the space charge capacitance (CSC) at the film/solution interface and the classical Helmholtz capacitance (CH) [8]:(1)1CT=1CSC+1CH.
Assuming that the Helmholtz capacitance is much larger than the space charge capacitance, the measured CT is equal to CSC. Therefore, the data points in Csc-2 versus E plots can describe the semiconducting behavior of the passive film. Given a correct value of CSC capacitance, the electronic properties are usually deduced from the Mott-Schottky relation [8, 9] as follows:(2)1CT2=1CSC2=2εε0qNq(E-Efb-kBTq),
where q is the elementary charge (−1.602 × 10−19 C for an electron and +1.602 × 10−19 C for electron hole), Nq is the density of charge carriers (NA for acceptors and ND for donors), (ε is the dielectric constant of the passive film, which is assumed as 15.6 [10]), ε0 is the vacuum permittivity (8.854 × 10−14 F/cm), kB is the Boltzmann constant (1.38 × 10−23 J/K), T is the absolute temperature, E is the applied potential, and Efb  is the flat-band potential. At room temperature, kBT/q is neglected because it is only about 25 mV. For a p-type semiconductor, Csc-2 versus E should be linear with a negative slope that is inversely proportional to the acceptor density. On the other hand, an n-type semiconductor yields a positive slope which is inversely proportional to the donor density.Figure3 shows the Mott-Schottky plots obtained at different frequencies for a passive film formed on surface of 316 in 0.05 M H2SO4. As was expected, the Mott-Schottky plots are frequency dependent. However, one can note that below 1 kHz the experimental data are quite stable. Yet, this allowed us to set a threshold frequency in the experimental protocol before conducting any experiments on the surface. According to (2), the slope of the linear portion of Csc-2 versus E plot, gives the charge carrier density (Nq) from the following equation:(3)Nq=2q⋅ε⋅ε0⋅s,
where s is the slope of the Mott-Schottky plot in the linear-region of interest. A higher frequency was selected to minimize the contribution of Helmholtz capacitance, which results from charge redistributions within the solution at the surface. Helmholtz capacitance neither affects the Csc-2 versus E plot, nor the calculation of charge carrier density within the passive film, as only the linear region of the graph is considered to calculate the slope [10].Figure 3
Mott-Schottky plots at different frequencies of the passive films formed on surface in 0.05 M H2SO4. The electrode is passivated for 1 h at 0.6 VSCE.The Mott-Schottky plots for the passive films formed at different film formation potentials on 316 are shown in Figure4. The film formation potential varies between −0.1 and 0.6 VSCE. Four typical examples of Mott-Schottky plots measured at different film formation potentials are shown in Figure 4. From this figure, a linear part can be easily distinguished for potentials greater than 0.0 VSCE. The decrease of capacitance with the applied potential in Figure 4 is attributed to an increase in the thickness of the electron-depleted layer and a diminishing number of charge carriers (donor for an n-type semiconductor), that is the region close to the film/solution is depleted in electrons.Figure 4
Mott-Schottky curves for passive films formed on 316 at different film formation potentials for 1 h in 0.05 M H2SO4.The slopes of the Mott-Schottky curve of the passive film formed at various formation potentials were different. At a formation potential lower than 0.6 VSCE, the value of the Mott-Schottky curve slopes increased with the increase of the formation potential. Therefore, in the following analysis, the potential range from −0.1 to 0.5 VSCE for 316 was chosen to explain the relationship between donor density and formation potential.The observed values ofND obtained for 316 are in the range of 2-3 × 1021 cm−3, which is comparable, and in the similar range reported for the passive film of stainless steels. Di Paola reported a donor density of 3.8 × 1021 cm−3 for 304 passivated in 0.5 kmol/m3 H2SO4 at 0.2 VSCE, which agrees with the present data [11]. The density of charge carriers has been observed to be a strong function of the passive film formation potential. The donor density of passive film of 316L was observed to decrease with an increase in the formation potential.In general, the changes inND do not correspond to impurity concentration but rather to the nonstoichiometry defects in the space charge region or due to the disordered character of the passive film. The donors in semiconducting passive layers are defects themselves, including oxygen ion vacancies and/or cation interstitials. These defects act as dopants, that is, oxygen ion vacancies and/or cation interstitials imparting n-type character. Both oxygen ion vacancies and cation interstitials are electron donors. The presence of such dopant prevents migration of cations from alloys and penetration of harmful anions, such as SO42-, from the solution, thereby improving the corrosion resistance. The decrease in ND means impoverishment of these species, which neutralizes the dopant. Higher ND values are strong indication of a non-stoichiometric or highly disordered passive film [12].


## 3.1. Polarization Behavior

Figure1 shows the potentiodynamic curve for 316 electrode in 0.05 M H2SO4. From the polarization curve, the passive range was determined to be from about −0.2 to 0.6 VSCE. It was observed that before the electrode was transferred to a passive state, an active current peak occurred at around −0.15 VSCE, which could be attributed to the oxidation of Fe2+ to Fe3+ions in the passive film [6].Figure 1
Potentiodynamic polarization curve of 316 in 0.05 M H2SO4 solution with 1 mV/s scan rate.Films were grown at each potential for 1 h to insure that the system was in steady state. During the formation of the oxide films on working electrode, the current density was monitored. It was observed that the current density diminishes with time until a constant value is reached and a steady-state established. Figure2 shows the values of the steady state passive current density (iSS) versus the film formation potential. The steady-state current density is approximately 6.8 μA cm−2.Figure 2
Steady-state current density obtained during the potentiostatic growth of the passive films at different film formation potentials for 1 h.


## 3.2. Mott-Schottky Measurements

The capacitance of the passive film (CT) is considered to be the combination of two series capacitances: the space charge capacitance (CSC) at the film/solution interface and the classical Helmholtz capacitance (CH) [8]:(1)1CT=1CSC+1CH.
Assuming that the Helmholtz capacitance is much larger than the space charge capacitance, the measured CT is equal to CSC. Therefore, the data points in Csc-2 versus E plots can describe the semiconducting behavior of the passive film. Given a correct value of CSC capacitance, the electronic properties are usually deduced from the Mott-Schottky relation [8, 9] as follows:(2)1CT2=1CSC2=2εε0qNq(E-Efb-kBTq),
where q is the elementary charge (−1.602 × 10−19 C for an electron and +1.602 × 10−19 C for electron hole), Nq is the density of charge carriers (NA for acceptors and ND for donors), (ε is the dielectric constant of the passive film, which is assumed as 15.6 [10]), ε0 is the vacuum permittivity (8.854 × 10−14 F/cm), kB is the Boltzmann constant (1.38 × 10−23 J/K), T is the absolute temperature, E is the applied potential, and Efb  is the flat-band potential. At room temperature, kBT/q is neglected because it is only about 25 mV. For a p-type semiconductor, Csc-2 versus E should be linear with a negative slope that is inversely proportional to the acceptor density. On the other hand, an n-type semiconductor yields a positive slope which is inversely proportional to the donor density.Figure3 shows the Mott-Schottky plots obtained at different frequencies for a passive film formed on surface of 316 in 0.05 M H2SO4. As was expected, the Mott-Schottky plots are frequency dependent. However, one can note that below 1 kHz the experimental data are quite stable. Yet, this allowed us to set a threshold frequency in the experimental protocol before conducting any experiments on the surface. According to (2), the slope of the linear portion of Csc-2 versus E plot, gives the charge carrier density (Nq) from the following equation:(3)Nq=2q⋅ε⋅ε0⋅s,
where s is the slope of the Mott-Schottky plot in the linear-region of interest. A higher frequency was selected to minimize the contribution of Helmholtz capacitance, which results from charge redistributions within the solution at the surface. Helmholtz capacitance neither affects the Csc-2 versus E plot, nor the calculation of charge carrier density within the passive film, as only the linear region of the graph is considered to calculate the slope [10].Figure 3
Mott-Schottky plots at different frequencies of the passive films formed on surface in 0.05 M H2SO4. The electrode is passivated for 1 h at 0.6 VSCE.The Mott-Schottky plots for the passive films formed at different film formation potentials on 316 are shown in Figure4. The film formation potential varies between −0.1 and 0.6 VSCE. Four typical examples of Mott-Schottky plots measured at different film formation potentials are shown in Figure 4. From this figure, a linear part can be easily distinguished for potentials greater than 0.0 VSCE. The decrease of capacitance with the applied potential in Figure 4 is attributed to an increase in the thickness of the electron-depleted layer and a diminishing number of charge carriers (donor for an n-type semiconductor), that is the region close to the film/solution is depleted in electrons.Figure 4
Mott-Schottky curves for passive films formed on 316 at different film formation potentials for 1 h in 0.05 M H2SO4.The slopes of the Mott-Schottky curve of the passive film formed at various formation potentials were different. At a formation potential lower than 0.6 VSCE, the value of the Mott-Schottky curve slopes increased with the increase of the formation potential. Therefore, in the following analysis, the potential range from −0.1 to 0.5 VSCE for 316 was chosen to explain the relationship between donor density and formation potential.The observed values ofND obtained for 316 are in the range of 2-3 × 1021 cm−3, which is comparable, and in the similar range reported for the passive film of stainless steels. Di Paola reported a donor density of 3.8 × 1021 cm−3 for 304 passivated in 0.5 kmol/m3 H2SO4 at 0.2 VSCE, which agrees with the present data [11]. The density of charge carriers has been observed to be a strong function of the passive film formation potential. The donor density of passive film of 316L was observed to decrease with an increase in the formation potential.In general, the changes inND do not correspond to impurity concentration but rather to the nonstoichiometry defects in the space charge region or due to the disordered character of the passive film. The donors in semiconducting passive layers are defects themselves, including oxygen ion vacancies and/or cation interstitials. These defects act as dopants, that is, oxygen ion vacancies and/or cation interstitials imparting n-type character. Both oxygen ion vacancies and cation interstitials are electron donors. The presence of such dopant prevents migration of cations from alloys and penetration of harmful anions, such as SO42-, from the solution, thereby improving the corrosion resistance. The decrease in ND means impoverishment of these species, which neutralizes the dopant. Higher ND values are strong indication of a non-stoichiometric or highly disordered passive film [12].


## 4. Relationship betweenND and E

A number of models have been proposed to explain the kinetic, thermodynamic, structural, and electronic properties of passive film formation. Among all the models, PDM is widely accepted. PDM provides microscopic description of the growth and breakdown of the passive film under steady-state and transient conditions. This model is based on the migration of the point defects under the influence of the electrostatic field in the passive film. A passive film is envisaged to grow into the alloy by generation of the oxygen ion vacancies at the alloy/passive film interface and by their annihilation at the passive film/solution interface. Figure5 shows the dynamic and steady-state properties of a passive film expressed by five electrochemical reactions basing on the PDM [13]. (1)
Metal vacancies are annihilated and produced at the alloy/passive film interface:(4)m+Vm3-→r1Mm+3e,Mm→r2M3++Vm3-.(2)
Oxygen ion vacancies are produced and annihilated at the passive film/solution interface:(5)m→r3Mm+VO2++2e,(6)VO2++H2O→r4OO+2H+.(3)
The film dissolves chemically at the passive film/solution interface:(7)M2O3+6H+→r53H2O+2M3+.m represents metal atom, Mm metal cations in film, OO oxygen anions in film, Vm3-  metal vacancy, VO2+  oxygen ion vacancy, M3+ metal ion in solution, and M2O3 chemical expression of the passive film.Figure 5
Schematic representation of five electrochemical processes which occur within the passive film according to the PDM.Because the passive films shown-type semiconductive behavior, the donor density (oxygen ion vacancy) is much larger than that of the acceptor (metal vacancy). Reactions (5) and (6) are dominant in the model. The flux of oxygen ion vacancy during the passive film can be expressed by Nernst-Planck equation [13, 14]:(8)JO=-DO∂CO∂x+2KDOCO,
where DO is the diffusion coefficient of oxygen ion vacancy, K=keF/RT, and ke is the electric field strength of the passive film. The oxygen ion vacancy density in the passive film was obtained by solving the following equation:(9)CO(x)=(CO(L)-JO2KDO)exp⁡(-2KL)exp⁡(2Kx)+JO2KDO,
where L is the thickness of the passive film. The lower concentration of oxygen ion vacancy close to the passive film/solution interface results in the occurrence of an insulating layer. The almost constant concentration of oxygen ion vacancy close to the alloy/passive film interface should result in a linear Mott-Schottky plot. Accordingly, the donor density calculated from Mott-Schottky curves was the oxygen ion vacancy concentration at the alloy/passive film interface [13, 14].(10)ND=CO(0)=(CO(L)-JO2KDO)exp⁡(-2KL)+JO2kDO,ND=CO(L)exp⁡(-2KL)+JO2KDO(1-exp⁡(-2KL)).
As reaction (6) is dominant, JO equals the rate of production of oxygen ion vacancy at the passive film/solution interface. Accordingly,(11)JO=r4CO(L),r4=r40exp⁡(2α4γϕf/s),ϕf/s=ϕf/s0+αE+βpH,
where r4 is the rate constant for reaction (6). The rate constant depends on the potential drop across the passive film/solution interface, ϕf/s. r40 is the rate constant for ϕf/s0=0, α4 is the transfer coefficient, γ=F/RT, ϕf/s0 is the potential drop at the passive film/solution interface when E=0, pH=0, and α,β are constant parameters.(12)CO(L)=JOr4=JOr40exp⁡(2α4γϕf/s)=JOr40exp⁡(2α4γ(ϕf/s0+αE+βpH)),ND=JOr40exp⁡(2α4γ(ϕf/s0+αE+βpH))exp⁡(-2KL)+JO2KDO(1-exp⁡(-2KL)),ND=JOr40exp⁡(-2α4γ(ϕf/s0+αE+βpH))exp⁡(-2KL)+JO2KDO(1-exp⁡(-2KL)),KL=(keγ)L=(keFRT)L=(keFRT)L=FERT,ND=JOr40exp⁡(-2α4γ(ϕf/s0+αE+βpH))exp⁡(-2FERT)+JO2KDO(1-exp⁡(-2FERT)).
Since for thin passive films (<10 nm thick), E may be as large as 106-107 V/cm [14], the value of FE is much larger than RT. So,(13)1-exp⁡(-2FERT)⟶1,ω1=JOr40exp⁡(-2α4γ(ϕf/s0+βpH)),b=(-2FRT-2α4γα)<0,ω2=JO2KDO.
Therefore,(14)ND=ω1exp⁡(bE)+ω2.
So, the donor density of the passive film should decrease exponentially with the increase of the formation potential [15]. The donor density has been calculated from capacitance measurements for films formed at different formation potentials. The best exponential fit to the experimental data is shown in Figure 6. The data plotted in Figure 6 show that the donor density decreases with increasing the film formation potential, which is in agreement with the prediction of PDM.Figure 6
Curve of donor density as function of film formation potential of the passive films formed on 316 in 0.05 M H2SO4.Sikora and Macdonald have shown that the relationship betweenND and the film formation potential can be developed theoretically on the basis of the PDM, assuming that the defects are the electron donor. This theoretical relationship yields a good fit to the experimental results and allows the diffusivity of the defects in the passive film to be calculated. The relevant equation describing the dependence of the donor density on the film formation potential is shown as (14), where ω1, ω2, and b are unknown constants that are to be determined from the experimental data [5]. Based on the nonlinear fitting of the experimental data by software, the exponential relationship between the donor density (ND) and the film formation potential (E) is derived as(15)ND=0.392×1021exp⁡(-4.37E)+2.332×1021.
According to the PDM, the flux of oxygen ion vacancy through the passive film is essential to the film growth process, which supports the existence of oxygen ion vacancy in the film regardless of its concentration. In this concept, the dominant point defects in the passive film are considered to be oxygen ion vacancies and/or cation interstitials acting as electron donors. However, as it is not possible to separate the contribution of oxygen ion vacancies and cation interstitials on the measured diffusivity value based on the PDM, the diffusivity is considered that for the combination of these two point defects [16]. The diffusion coefficient can be calculated from (15) by [6]:(16)DO=JO2Kω2=JORT2Fω2ke,(17)JO=-iSS2e.ke was determined to be approximately 1.02 × 106 V/cm by Szklarska-Smialowska and Kozlowski [16] for the passive film grown on iron. Therefore, it is assumed that the value of ke for the passive film on stainless steel to be approximately 106 V/cm. Substitution of ke (≈106 V/cm),   iSS (≈6.8 × 10−6) A cm−2), e (−1.602 × 10−19 C), ω2 (2.332 × 1021 cm−3), R (8.314 J/mol K), T (298 K), and F (96500 C mol/e) into (16) and (17) yields DO=1.17×10-16 cm2/s. Taking into account the error caused by the assumption of ke in (16), the diffusion coefficient of defects in the passive film formed on 316 in 0.05 M H2SO4 solution is estimated to be approximately in the range of  10−16 cm2/s.


## 5. Summary

Potentiodynamic polarization studies demonstrate that 316 displays a wide passive range in 0.05 M H2SO4 at ambient temperature. Potentiostatic polarization tests revealed that the steady-state current density (iSS) through the passive film formed on 316 for 1 h was independent of formation potential, which is well consistent with the postulation of the PDM. Based on the Mott-Schottky analysis, it was shown that the calculated donor density decreases exponentially with increasing formation potential. The experimental data were interpreted in terms of the PDM for the passivity of 0.05 M H2SO4, assuming that the donors are defects including oxygen vacancies and cation interstitials. The diffusivity of donors in the passive film on 316 was calculated to be approximately in the range of 10−16 cm2/s.



---

*Source: 968512-2011-01-05.xml*