Use of a small-voltage electrical response to measure sodium motion in rf sputtered SiO2 films

Use of a small-voltage electrical response to measure sodium motion in rf sputtered SiO2 films

Journal of Non-Crystalline Solids 33 (1979) 391-403 © North-Holland Publishing Company USE OF A SMALL-VOLTAGE ELECTRICAL RESPONSE TO MEASURE SODIUM M...

513KB Sizes 0 Downloads 40 Views

Journal of Non-Crystalline Solids 33 (1979) 391-403 © North-Holland Publishing Company

USE OF A SMALL-VOLTAGE ELECTRICAL RESPONSE TO MEASURE SODIUM MOTION IN RF SPUTTERED SiO2 FILMS M. MEAUDRE, R. MEAUDRE and A. DEGUIN Laboratoire de Physique Electronique, Universit~ de Lyon 1, 69621 Villeurbanne, France

Received 5 January 1979 Revised manuscript received 12 March 1979

Low voltage transient and ac responses of rf sputtered SiO2 films are studied in the temperature range 473-593 K. Experimental results can be explained in terms of a model with trapped carriers and with only one type of recombining mobile charge carriers. Two distinct. phenomena clearly appear in transient current curves; the rapid one is attributed to the space charge effect, the other to generation recombination. Theoretical and experimental curves are in good agreement. The concentrations and mobilities are deducted of carriers identified as Na÷ ions. Concentrations values are compared (for glass and silica substrates) with ion microprobe analysis results.

1. Introduction There are a considerable number o f publications concerning thermally grown SiO: films [1] but fewer projects have been carried out on rf sputtered SiO2 films, and those done have generally studied the influence o f experimental conditions on the thickness, density, structure and composition o f films [ 2 - 1 7 ] . Within the scope o f electrical measurements on rf sputtered SiO2 films we can quote: the investigation o f dielectric behaviour in the frequency range 103-107 Hz and dc leakage currents by Franklin [18]; the study o f thermally stimulated currents by Hickmott [19] to measure additional Na motion and b y Hino [20] and Hino et al., [21] to determine the relaxation time o f dipoles; the use o f conventional C - V method by Kumar and Gregor [ 14] and Schreiber and Fr6schle [ 16] to obtain the ion concentration or interface state density; the measurements o f breakdown fields by Franklin [18] Davidse and Maissel [2] and Schreiber and Fr6schle [16]. None o f these projects gives information on the motion o f the intrinsic charge carrier existing in the films. Therefore it was o f interest to study polarization processes likely to clarify this point. This investigation o f dielectric behaviour was performed b y transient and alternating current measurements. Such an analysis has never before been made neither in rf sputtered SiO2 films nor in thermally grown SiO: films. 391

392

M. Meaudre et al. / Sodium motion in R F sputtered SiO 2 films

2. Samples The rf sputtering equipment used for the deposition of SiO2 films was a diode sputtering system [22]. Typical data for the experiments were given in a previous paper [23]. The films obtained were analyzed by different methods. Electron diffraction has shown their amorphous character. It is important to state the amorphous character, as the position of absorption bands in the IR spectra can be influenced by crystalline zones present in the fihn. We obtained the infrared spectrum represented in fig. 1. The Si-O stretching band is shifted to 9.43/~m; for analogous displacements of the S i - O stretching band,Pliskin et al. [24] have calculated that these films are only slightly oxygen deficient, having a formula of SiOx with 1.95 < x < 1.99. Stoichiometry of our SiO2 films was confirmed by d, p nuclear reactions. Electrical measurements were made on A u - S i O 2 - A u samples. Au was pre-evaporated on glass (type I sample) or silica (type II sample) substrates. After sputtering of SiO2 an Au film was deposited to establish a capacitor Au-SiO2-Au. The contact areas were 3.14 mm 2. SiO2 film thickness varied between 500 and 3000 A. Before electrical measurements, in order to eliminate moisture, which generally masks true conductivity and frequency spectral characteristics, the films had been heated at 200°C in the vacuum chamber with a pressure lower than 10 -s Torr for about five days. 9

I0

WAVEL£lqCTH (MICRONS)

15

20

25

100

I

80

A o

i

40

20

I

I

1300

1200

I

1100

I

1000

I

I

1

I

900

800

700

600

~W~L~ER

I

500

(CM-I }

Fig. 1. Infrared spectra for a rf sputtered SiO2 film, about 3000 A thick.

400

393

M. Meaudre et al. /Sodium motion in RF sputtered SiO 2 films

3. Experimental 3.1. M e t h o d

Electrical measurements have been made in a vacuum lower than 5 X 10 -6 Torr. The transient current response due to a step voltage of 50 mV was measured with a Keithley model 610 C electrometer and recorded. Capacitances Cp and conductances G o were inferred from the following method: a low frequency sinusoidal voltage u(t) = Re[V exp(/cot)] was applied to the sample and the ac current i(t) passing through it was measured. One has: i(t) = Re[(Gp +/Cpco) V(cos cot +/sin co0]

= GpVcos cot - CpcoV sin cot. The resulting i(t) = f(o) curve is an ellipse, i(t) was measured in the frequency range 5 X 10-4-0.2 Hz with a Keithley model 610 C electrometer and i(t) =rio)

1

DETAIL

I ~ 10-

~1~

.

~'~260°C° •

*

--.

lo-1

,~ \

l

~

"~o

o

280"C

319,5oc lOO 200

300

400

10_2

~ ~ ~ 10-3

~ 2 8 0 0 c -''~'e'-~v o 299 5°C

319, 5°C IO=I,8XI0-9A

i¢7,4x10-11A IO=2'3XI0-IOA

I0=6,9XI0-i0 A I

I0000

I

20000

I

30000

T I E (SEC) Fig. 2. Normalized current versus time isotherms for type I f'dm (P17C1 sample). Applied voltage: 50 mV; SiO2 f'dm thickness: 3000 A; relative dielectric constant er = 3.8; experimental points and theoretical curves calculated with parameters given in table 1.

394

M. Meaudre et al. /Sodium motion in R F sputtered SiO 2 film* e~

o

d

_E

>

o t~

Io

o

o

110

C~

< E

o

,~"a

l

.~

0 O Q

o 0I -I

o /~T68~D (I~ZI~A~ON

? o

M. Meaudre et al. / S o d i u m motion in R F sputtered SiO 2 films

395

10 5

Io

104

22o 200,5 °C "

~

\

\ N

\. N

~

.

103

l

10.4

10.3

1

I

I

10. 2

10-1

1

(Hz)

(a}

103

|





259,5"C



239,5"C



220"C

,

200,5"C

10 2

10

10 -4

I

I

10-3

10.2

I

10 -1

I

I

FRDQtJ~Y (ltz) (b) Fig. 4. (a) Parallel capacitance, (b) parallel conductance versus frequency, for the P2oC 6 sample at various temperatures. Experimental points and theoretical curves calculated with parameters given in table 1.

396

M. Meaudre et al. / Sodium motion in RF sputtered SiO 2 films

curves were plotted using a Hewlett Packard model 7000 A / X - Y recorder. A Tektronix model 547 oscilloscope was used in the frequency range 0 . 2 - 2 0 Hz. For ~otl = 2k~r ,

v(tl) = V ,

i q l ) = GpVand Gp : i ( t l ) / V .

For cot2 = (rr/2) + 2klr,

o(t2) = 0 ,

i(t2) = -CpeoVand Cp = - i q 2 ) / V w .

3.2. Transient current

The time response of the current due to a step voltage was studied in the temperature range 473-593 K. The current-time curves obtained for samples P17C~ (type I sample) and P2oC 6 (type II sample) are shown in figs. 2 and 3. Depolarization currents were also studied. Polarization and depolarization currents followed Ohm's law closely for applied electric fields up to 10 6 V m -I. The principle of superposition was obeyed. 3.3. L o w frequency behaviour

Dielectric behaviour of sample P2oC6 (type II sample) was investigated from 5 × 10-4-20 Hz in the temperature range 473-530 K. Figs. 4(a) and 4(b) show a dielectric dispersion at low frequencies, this dispersion is characterized by large variations in capacitance. Alternating currents were not measured in sample P~7C1 but the dielectric behaviour of analog type I films was previously given [23]. The static capacitances of type I samples was found to be larger than those of type II samples.

4. Results

The large experimental values of static capacitance and the fact that they do not depend practically on electrode separation L (for type I film having L = 600 A, 1200 A, 2500 A, we measured respectively 1.8/aF, 2/~F, 1.6/IF) lead us to exclude a dipolar model or the Friauf or Macdonald explanation [25,26]. Jaffe and Lemay [27] have given an expression for the transient current in the case of non-recombining positive and negative carriers with equal mobilities and the same boundary conditions. They assumed a uniform electric field so that Poisson's equation is not obeyed and the voltage is at least as large as 8KT/e. The theoretical time constant ratios were not observed experimentally. In the same case of non-recombining carders with equal mobilities Macdonald [28] gave an expression for the current which only applies to interfaces, and the corresponding (non-blocking situation) static capacitance is proportional to L. Fatuzzo and Coppo [29] have taken into account

M. Meaudre et al. / Sodium motion in RF sputtered Si02 films

397

the space charge layers spontaneously formed at the crystal surfaces; they assumed completely blocking electrodes and equal carrier mobilities. The transient current is governed by only one voltage dependent time constant; this fact and the complete blocking of carriers were not observed by us. Fonash et al., [30-31] have investigated the case of Schottky barriers but the relation between initial currents and time constants were not satisfied, furthermore an electronic conduction seemed unlikely in our low voltage low frequency situation. Our model [32] taking into account generation recombination at all the states of the gap can give theoretical dielectric dispersion curves similar to experimental ones but leading to a high and unacceptable density of states. Finally when only one type of carriers can recombine at a single level we have shown that the transient response is made up of the sum of two exponential decays [33-35]; when certain experimental curves are fitted to this model it is impossible to satisfy the theoretical relations existing between the different coefficients. For all the previous reasons we studied a slightly more sophisticated model [36]. We considered a sample with concentration M d of immobile neutral centres before dissociation. After dissociation n(x, t) carriers of mobility g and M a trapped carriers appear, nd(x, t) centres are not dissociated. The number of carriers produced per second at position x is proportional to the product of na(x, t) and a generation coefficient k, the number of recombinations per second is ~n(x, t) [Ma - na(x, t)]. When a low step voltage V is applied at time t = 0, the transient current for electrodes of unit area may be stated in the following form:

n oep V J(t) : Z[l+(~Z/2O'] [(1 +

~(Ma +n°)~2po] e x p ( - P ° ( 1 + ~(Ma+n°)~t]2po /]

2/3(Ma+no) fl(Ma+no)(~ p~ /3no exp(-/3not) ÷ Po

~nO]exp(_flnot)+ ~D1 Po /

(1) no, D and ~ represent the equilibrium concentration, the diffusivity and the discharge parameter of the carriers respectively.

= Fnoe, ( 2

-'

is the space charge relaxation time, L is the sample length, LD -- (e.D/noeP) 112 the Debye length and e the static dielectric constant.

¢ = ¢[½, 2, - ~(no +M.) t] is the degenerate hypergeometric function. The complex dielectric constant e' - ] e " is then:

e'

LBrl

6LC

1 +¢o=r~

po(l+co=r22)

- - - +

LCT" 3

LCT"4

l+co2ra:

1+¢o2r]

+ - - + - - + e

(2)

398

M. Meaudre et al.

/ Sodium

motion in R F sputtered SiO 2 films

X X X X X

m

~T

T XX

X X X X X

XX

~

V

X

X X X X X

~

X

V

X X X X X

X X X X X

~T

X

T

X X X X X o

I

~>~_,-..

~T

77

X X X X X

X

M. Meaudre et al. /Sodium motion in RF sputtered Si02 films e"-

LBwr~

6LCcor2

LC~z]

LCwr~

399

F

1 + co2r~ + po(1 + ~o2r~) t- 1 +W27"~ ÷ 1- +w2r~ + co

(3)

with

B:L[1

F =

noela ( +(~L/2D)] 1

noela ~L L [1 + (~L/2D)] 2 D '

r31 =/3no ÷ --0"14Crll B +4C '

~(Ma + no)~ 2p~

]'

~(m a + no)

noela

C=L[I+(~L/2D)]

7"11 =p0(1 + fl(Ma +no)] 2po

4po

T~I =7"; 1 =flno,

/ '

1.74Crl 1 r41 =/3no + - B +4C

When using a least squares computer program our experimental results were successfully fitted with eq. (1). Second order terms such as 13(Ma + no)~no/P2o are not taken into account during the fitting of transient current curves, nevertheless they were given in eq. (1) because they lead to first order terms in the expression of the complex dielectric constant. Theoretical curves and experimental points are shown in figs. 2 and 3; a very good agreement is obtained and this was for transient current durations as long as eight hours. Values of parameters are listed in table I, they have also been used in eqs. (2) and (3) to calculate the theoretical Cp and Gp curves of fig. 4. One can note from the Gp curves that for frequencies larger than 0.1 Hz, another phenomenon appears, a work in progress shows that a hopping process is very likely.

5. Discussion Electronic conduction may be due to the motion of free carriers or alternatively to the motion of quasi-localized carriers. Mobility values o f free carriers [ 3 7 - 3 9 ] are considerably larger than the ones given in table 1. Mobility values o f carriers hopping between localized sites can be very low and only the smallness of the activation energy remains a plausible argument for determining the nature of the carriers. As a general rule, activation energies in excess o f 0.6 eV indicate an ionic transport [40], as activation energies deducted from figs. 5(a) and 5(b) are 1 eV we conclude electronic transport is not responsible for the observed effects. Furthermore, ionic conduction is characterized by low mobilities, indeed, in fused silica Lee [41] reported a mobility of 2.8 × 10 -12 m 2 V -I s -1 for H ÷ at 140°C, Owen and Douglas [42] reported one of 2.8 × 10 -18 m s V -1 s -1 for Na ÷ at 140°C; in thermally grown SiO2 films Hofstein [43] measured a mobility of 1.4 × 10 -13 m 2 V -1 s -1 for H ÷ at 140°C, Sugano et al., [44] reported one of 6 × 10 -16 m s V -1 s -1 for Na ÷ at 200°C and Snow et al., [45] used a mobility value found by Owen and Douglas to explain their results. Though slightly lower, our mobility values may be compared with the previous

400

M. Meaudre et al.

/ Sodium

motion in R F sputtered SiO 2 films

10-16 A !

%

10-17 *

1 eV.

10-18 1,7

1,8

1,9

F~CIPBOCAL TE/~ERA'IURE (10 -3 "K-1)

(a)

!

%

10-16

,.o

10-17 I

I

1,9

2

1

~'~

2,1

I~CIPBOCAL 'I'I~ERATURE (10 -3 *K-1)

(b) Fig. 5. Plots of carrier mobility as a function of reciprocal temperature. (a) Type I film (P17C1 sample), (b) type II film (P2oC6 sample).

ones for Na ÷. Activation energy f o r Na ÷ mobility in SiO2 was found 0 . 9 - 1 . 4 eV [ 4 4 - 4 5 ] . According to Kriegler and Bartnikas [46] values of about 1 eV, which is in good agreement with ours, would really characterize the mobility, whilst higher values would characterize trapping effects at the interfaces. So observed conduction mechanisms seem t o be due to Na ions. Analysis of films by an ion probe microanalyzer confirmed the presence of sodium. Profiles of Na concentration are given in figs. 6(a) and 6(b) for type I and type II films. Na concentrations deducted from electrical measurements are higher than Na concentrations measured by analysis, but both methods give the same variation of content between films o f type I and type II. Consequently polarization effects are undoubtedly due to Na ion motion. In the case of ionic conduction, according to some authors, [ 4 7 - 5 0 ] , the Ein-

401

M. Meaudre et al. / Sodium motion in R F sputtered SiO 2 films

-A 1019

1021

I

I

7

1o2O1~

~ 10 1

r.

iii o o u

zo I

1019h 0

,

I

,

100o

I 2000

Depth (al



,

3000

(A)

.



.

|

5OO

Depth

|

10oo

(~)

(b)

Fig. 6. Na profiles derived from ion microprobe analysis for (a) type I film (PI 7C1 sample), (b) type II film (P20C 6 sample).

stein relation may be replaced by I~/D = e/fKT where f is a correlation factor smaller than unity; this leads to smaller concentration values ( × f ) a n d thus to a better agreement with microprobe analysis results. One has now to justify the set of equations [36] leading to eqs. (1), (21) and (3). In SiO2 films, some S i - O - S i linkages are broken, resulting in oxygens attached to only one silicon, or non-bridging oxygens. According to Charles [51 ] there are positions of stability for Na ÷ ions about non-bridging oxygens separated from each other by energy barriers. We then consider that the M a immobile neutral centres existing per unit volume before dissociation, consist of non-bridging oxygen, that is a negatively charged oxygen, with an Na ÷ ion in the neighbourhood [52]. After dissociation n(x, t) positive carriers of mobility/a (Na ÷ ions) and M a trapped carriers appear, ha(X, t) centres are not dissociated. Charge motion is due to the applied electric field and diffusion process. Na ÷ ions are partially blocked at the electrodes. Eqs. (1), (2) and (3) have been obtained in the case of negative mobile and trapped carriers but it is easily shown that changing the sign of these carriers does not change anything in the final result. The presence of trapped ions has often been observed [53-56]. It must be underlined that our/3no values are of the same order of magnitude as the k's derived from Hickmott [54] in the case of rf sputtered SiO2 films for 1-1.65 eV trap depth.

402

M. Meaudre et al. / Sodium motion in R F sputtered SiO 2 films

6. Conclusion Low voltage electrical behaviour of rf sputtered SiO2 films has been studied by t r a n s i e n t and alternating Current measurements. The interpretation o f experimental results has necessitated the change o f a previously published model by adding trapped carriers which do not take part in conduction. Good fitting o f experimental curves has allowed us to obtain the concentration and mobility o f carriers considered as Na ÷ ions. When glass substrates with a high Na content are used, SiO2 films are polluted by sodium contained in substrates. When silicon dioxide substrates are used, very clean films are obtained since the Na content was o f the same order as that estimated by Ridley [57] for the purest thermally grown SiO2 films. Differences between the Na concentrations obtained by microprobe analysis and deduced from electrical measurements probably arise from our model which is too simple. Different boundary conditions leading to smaller concentrations could be envisaged [58,59], unfortunately the corresponding expressions for the transient current are not yet available. Furthermore it is very likely that the main discrepancy arises from non-uniform ionic distributions such as those appearing in figs. 6(a) and 6(b) and not taken into account by the theory.

Acknowledgements The authors wish to thank Mrs A. Joly and Mr R. Chifflet for computer analysis, and Dr B. Blanchard for ion microprobe analysis.

References [1] A.H. Agajanian, Solid St. Technol. 20 (1977) 36. [2] P.D. Davidse and L.I. Maissel, J. Appl. Phys. 37 (1966) 574. [3] W.A. Pliskin, P.D. Davidse, H.S. Lehman and L.I. Maissel, I.B.M.J. Res. Develop. U.S.A. 11(1967) 461. [4] W.A. Pliskin, D.R. Kerr and J.A. Perri, Phys. Thin Films U.S.A. 4 (1967) 257. [5] W.A. Pliskin, Thin Solid Films 2 (1968) 1. [6] R.E. Jones, H.F. Winters and L.I. Maissel, J. Vac. Sei. Technol. U.S.A. 5 (1968) 84. [7] J.G. Carver and B,J. Franklin, Signals Research and Development Establishment Report no. 69018, Ministry of Technology, Christchurch, Hants (1969). [8] W. Hoffmeister and B.J. Zuegel, Thin Solid Films 3 (1969) 35. [9] J: Gyulai, Acta Phys. Chem. 16 (1970) 119. [10] D.H. Grantham, E.L. Paradis and D.J: Quinn, J. Vac. Sci. Technol. 7 (1970) 343. [11] G.C. Schwartz and R.E. Jones, I.B.M.J. Res. Develop. U.S.A. 14 (1970) 52. [12] T.W. Hiekmott, J. Appl. Phys. 42 (1971) 2543. [13] S. Petersson, G. Linker and O. Meyer, Phys. Stat. Solidi (a)14 (1972) 605.

M. Meaudre et al. /Sodium motion in R F sputtered SiO 2 films

[14] [15] [16] [17] ] 18]

403

S. Kumar and L.V. Gregor, J. Electrochem. Soc. 120 (1973) 1285. R.S. Butlin and M.K. Mac Phun, I.E.E.E. Conf. Publ. (1975). H.U. Schreiber and E. Froschle, J. Electrochem. Soc. 123 (1976) 30. K. Thiemt and R. Kaltofen, Thin Solid Films 38 (1976) 193. B.J. Franklin, Signals Research and Development Establishment, report no. 69020, Ministry of Technology Christchurch, Hants (1969). [19] T.W. Hickmott, Appl. Phys. Letters 22 (1973) 267. [20] T. Hino, Japan J. Appl. Phys. 11 (1972) 1573. [21] T. Hino, K. Suzuki and K. Yamashita, Japan J. Appl. Phys. 12 (1973) 651. [22] L. Holland, T. Putner and G.N. Jackson, J. Sci. lnstrum. 1 (1963) 32. [23] M. Meaudre and R. Meaudre, Solid St. Electron, 16 (1973) 1205. [24] W.A. Pliskin, P.D. Davidse, H.S. Lehman and L.I. Maissel, I.B.M.J. Res. Develop. U.S.A. 11 (1967) 461. [25] R.J. Friauf, J. Chem. Phys. 22 (1954) 1329. [26] J.R. Macdonald, J. Chem. Phys. 58 (1973) 4982. [27] G. Jaffe and C.Z. Lemay, J. Chem. Phys. 21 (1973) 920. [28] J.R. Macdonald, J. Chem. Phys. 54 (1971) 2026. [29] E. Fatuzzo and S. Coppo, J. Appl. Phys. 43 (1972) 1457. [30] S.J. Fonash, J.A. Roger, J. Pivot and A. Cachard, J. Appl. Phys. 45 (1974) 1223. [31 ] S.J. Fonash, J.A. Roger and C.H.S. Dupuy, J. Appl. Phys. 45 (1974) 2907. [32] R. Meaudre and M. Meaudre, Phys. Stat. Solidi (a)37 (1976) 633. [33] R. Meaudre and G. Mesnard, J. Phys. C: Solid St. Phys. 7 (1974) 1271.. [34] R. Meaudre, J. Phys. C: Solid St. Phys. 9 (1976) 599. [35] J.R. Macdonald, D.R. Franceschetti and R. Meaudre, J. Phys. C: Solid St. Phys. 10 (1977) 1459. [36] R. Meaudre and M. Meaudre, Phys. Rev. (1979) in press. [37] R.C. Hughes, Phys. Rev. Letters 30 (1973) 1333. [38] R. Williams, Phys. Rev. 140 2A (1965) A569. [39] O.L. Curtis, Jr, J.R. Srour and K.Y. Chiu, J. Appl. Phys. 45 (1974) 4506. [40] A.K. Jonscher, Thin Solid Films 1 (1967) 213. [41 ] R.W. Lee, J. Phys. Chem. Glasses 5 (1964) 35. [42] A.E. Owen and R.W. Douglas, J. Soc. Glass. Technol. 43 (1959) 159. [43] S.H. Hofstein, I.E.E.E. Trans. Electron Devices ED 13 (1966) 222. [44] T. Sugano, K. Hoh, K. Kudo and N. Hishinuma, Japan J. Appl. Phys. 7 (1968) 715. [45] E.H. Snow, A.S. Grove, B.E. Deal and C.T. Sah, J. Appl. Phys. 36 (1965) 1664. [46] R.J. Kriegler and R. Bartnikas, I.E.E.E. Trans. Electron. Devices 20 (1973) 722. [47] R.H. Doremus, Modern aspects of the vitreous state (Butterworths, Washington, 1962) ch. 2, p. 1. [48] Y. Haven and B. Verkerk, Phys. Chem. Glasses 6 (1965) 38. [49] M.E. Milberg, Fast Ion Transport in Solids (North Holland, Amsterdam and London, 1973) p. 373. [50] A.D. Le Claire, ibid. (1973) 51. [ 51 ] R.J. Charles, J. Appl. Phys. 32 ( 1961) 1115. [52] N.F. Mott, Rev. Phys. Appl. 12 (1977) 619. [53] E. Yon, W.H. Ko and A.B. Kupper, I.E.E.E. Trans. Electron. Devices ED 13 (1966) 276. [54] T.W. Hickmott, Appl. Phys. Letters 22 (1973) 267. [55] T.W. Hickmott, J. Appl. Phys. 46 (1975) 2583. [561 C.M. Osburn, J. Solid St. Chem. 12 (1975) 232. [57] B.K. Ridley, J. Appl. Phys. 46 (1975) 998. [581 S. Lanyi, J. Phys. Chem. Solids 36 (1975) 775. [59] D.R. Franceschetti and J.R. Macdonald, J. Electr. Chem. 82 (1977) 271.