Rate Constant for Ion Transfer in Inhomogeneous Media at the

Rate Constant for Ion Transfer in Inhomogeneous Media at the Interface of Immiscible ... Chem. B , 1997, 101 (50), pp 10801–10806 ... B 101, 50, 108...
0 downloads 0 Views 152KB Size
J. Phys. Chem. B 1997, 101, 10801-10806

10801

Rate Constant for Ion Transfer in Inhomogeneous Media at the Interface of Immiscible Electrolytes Kyo1 sti Kontturi,*,† Jose´ A. Manzanares,‡ Lasse Murtoma1 ki,† and David J. Schiffrin§ Laboratory of Physical Chemistry and Electrochemistry, Helsinki UniVersity of Technology, P.O. Box 6100, FIN-02015 HUT, Finland, Department of Thermodynamics, Faculty of Physics, UniVersity of Valencia, E-46100 Burjassot, Spain, and Chemistry Department, UniVersity of LiVerpool, P.O. Box 147, LiVerpool, U.K. L69 3BX ReceiVed: June 24, 1997; In Final Form: September 16, 1997X

An expression for the impedance due to ion transfer across the interface of immiscible electrolyte solutions (ITIES) is derived, assuming a continuous change of the ionic resolvation energy across an interfacial layer. The model is analogous to the Goldman constant field approximation. The result is practically identical to the one obtained assuming an activated step occurs during the transfer, i.e., the classical Butler-Volmer relationship, so that the standard rate constant k0 in the Butler-Volmer theory is replaced by the quantity D/δ, where D is the average diffusion coefficient of the transferring ion inside the interfacial layer and δ is the thickness of the interfacial layer. At potentials close to the standard transfer potential, the apparent chargetransfer coefficient R takes the value of 1/2. Combination of the charge-transfer impedance and an inner layer capacitance in parallel gives a complete semicircle in an impedance plot, characteristic of a parallel RC circuit.

Introduction The mechanism of the ion transfer across molecular interfaces is still unclear partly because both the structure of these interfaces and the Galvani potential distribution present are still poorly understood. Different interfacial models have been proposed. In the original work of Gavach,1 Buck et al.,2 and Samec et al.,3 a sharp boundary with two ion-free inner layers, one at each side of the interface, was proposed, similar to a modified Verwey-Niessen model. This was a direct transposition of the classical electrochemical formalisms for electrontransfer reactions at metallic electrodes. This interfacial model has been more recently applied to the study of interfacial electron-transfer reactions using the scanning electrochemical microscope.4 In contrast to treating interfaces as molecularly sharp boundaries, a mixed solvent layer was suggested by Girault and Schiffrin.5 More recently, comprehensive molecular dynamics simulations by Benjamin6 indicate that these boundaries are highly disordered, with “fingering” phenomena prevailing. The possible influence of these interfacial structures on ion transfer has been evaluated.6b The main question being discussed at present is whether ion transfer is a process containing an activated step corresponding, for instance, to resolvation or is it just a transport problem. If the former approach is accepted, the potential dependence of the transfer current will follow a classical Butler-Volmer relationship.7 The difficulty with this model is the uncertainty regarding the potential distribution across the interfacial region. This is particularly important if it is considered that most of the potential drop normally occurs in the organic phase. When the Butler-Volmer model is applied to the analysis of experimental data, the standard rate constant, k0, takes rather high values, on the order of 0.1 cm s-1. Recently, Girault et al. measured the impedance of a microhole ITIES which facilitates the mass transfer on both sides of the interface, but could not †

Helsinki University of Technology. University of Valencia. § University of Liverpool. X Abstract published in AdVance ACS Abstracts, December 1, 1997. ‡

S1089-5647(97)02059-2 CCC: $14.00

find any charge-transfer resistance due to ion transfer, and concluded that k0 has to take values on the order of at least 10 cm s-1.8 A proposed alternative model considers the interfacial region as an inhomogeneous phase through which the trace ion is transferred. The solvent properties and salt concentrations of the adjoining phases are considered to vary continuously through this interfacial region. This model has been considered by Kakiuchi,9 Kontturi et al.,10 and Senda,11 who integrated the Nernst-Planck equation across the interfacial region under steady-state conditions. The current-voltage characteristics thus obtained could be expressed as that of a first-order heterogeneous reaction, which enabled the comparison with the ButlerVolmer approach. The fundamental result of these analyses is that the Nernst-Planck model explains the salient features of the kinetic data, such as the nonlinear current-voltage characteristics and a value close to 1/2 for the transfer coefficient R. Furthermore, the formal standard rate constant can be related to the diffusion coefficient of the transferring ion across the interfacial region and to its thickness. Quite recently, Aoki presented a more general theory of ion transfer12 which includes these two approaches as its limiting cases. The comparison of the above two formalisms is not obvious for fast ion-transfer reactions. Rate constants are conveniently measured by perturbation techniques, such as ac impedance spectroscopy.7a Unfortunately, the classical analysis based, for example, on a Randlestype equivalent circuit, includes explicitly a particular potential dependence of the rate constant, in all cases, giVen as the Butler-Volmer equation. The question is whether the rate constant derived from transient techniques is uniquely related to the use of a specific kinetic model or, indeed, no particular kinetic equation for ion transfer needs to be assumed and all the kinetic information obtained from these techniques is simply a consequence of the diffusion and migration characteristics of an ion in an inhomogeneous region as described by the NernstPlanck equation. Previous theoretical ac impedance analysis in systems related to the present one has been carried out by © 1997 American Chemical Society

10802 J. Phys. Chem. B, Vol. 101, No. 50, 1997

Kontturi et al. and considering that the activity coefficient of the tracer ion is independent of concentration, the flux is given by

[

J ) -D

Figure 1. Schematic presentation of the interfacial region.

Buck et al.13a,b More recently, this analysis has been extended by the same authors by the use of transmission lines to the study of ion- and electron-transfer reactions in membranes.13c,d The numerical analysis by Buck13a referred to the membrane case, in which activated ion transfer at the membrane-solution boundary takes place. The purpose of the present work is to analyse the ac response for ion transfer across the interface of immiscible electrolytes and compare the results derived from the solution of the diffusion and migration problem with those for which an activated step process (the Butler-Volmer approach) was considered. In particular, it was regarded important to be able to ascertain if the observation of an ac response containing a classical charge-transfer resistance would be equally well derived from the phenomenological migration and diffusion description of ion transport. Theory The case considered here is that of a cation from the aqueous phase transferring across the interface. The concentration of the cation is cw in the aqueous phase and co in the organic phase. Concentration polarization is neglected so that interfacial concentrations are equal to those in the bulk; the effect of the concentration polarization was studied in ref 10. The reason for doing this is to highlight the nature of interfacial ion transfer. The diffusional components (Warburg impedance) in both adjoining phases can be simply taken into account in much the same way as is usually done in the impedance theory for metal electrodes.14 The Warburg impedance can then be subtracted from the experimental data.15 Furthermore, it is assumed that the entire changes of the Galvani potential and standard chemical potential of the ion take place in the interfacial region, the thickness of which is denoted by δ; see Figure 1. The latter assumption may not be quite correct; that is, the change of Galvani potential probably extends over a wider distance than δ. Refining the model accordingly would, first, bring about unnecessary mathematical complications and, second, have only a minor effect on the final result, as we have shown in ref 10. The analysis is carried out assuming the presence of supporting electrolytes in both phases. Two different cases will be considered: (i) equilibrium impedance and (ii) overpotential impedance. Solution of the General Problem. The ion transfer is described as an electrodiffusion problem, and the flux of the transferring ion, J, follows the phenomenological equation16

Dc dµ˜ J)RT dx

(1)

where D and c are the diffusion coefficient and concentration of the ion being transferred and µ˜ is its electrochemical potential; R and T have their usual significance. For a cation transfer

]

∂c(x) ∂(Φ(x) - Φ0(x)) + c(x)f ∂x ∂x

(2)

where c(x) is the concentration of the tracer ion a distance x from the aqueous phase boundary (taken as the origin of coordinates), f ) zF/RT, and Φ0(x) ) -µ0(x)/zF. µ0(x) is the standard chemical potential of the tracer ion for the solvent composition at point x. The quantity µ0(x) recognizes that ion transfer occurs across an inhomogeneous interfacial region, in which not only potential changes but also the chemical solvation energy gradient must be considered.17 The inner potentials given in eq 2 are absolute quantities, but integration across the interfacial region leads to more convenient relative quantities expressed with reference to the either of the phases. Here, we use the convention ∆wo Φ ) Φ(0) - Φ(δ) and ∆wo Φ0 ) Φ0(0) - Φ0(δ). Hence, the more positive ∆wo Φ - ∆wo Φ0 is made, the more the cation flux from the aqueous to the organic phase is accelerated. The solution of eq 2 under steady-state conditions (ss) is dependent on the assumptions made on the electrical and standard chemical potentials gradient in the interfacial region. For membrane or confined regions, the Goldman constant field assumption is commonly employed to describe the potential distribution.18 This approximation implies that no space charges are present in the interfacial phase. However, even in cases where the potential distribution deviates significantly from linearity, there is little difference in the calculated values of the current-voltage response.10,19 The latter observation is very important, since it shows that the fine details of the interfacial potential distribution have a small effect on the calculated values of the rate constants for ion transfer. Similar arguments apply to the chemical potential contribution to the driving force for interfacial ion transfer,10,19 and therefore, (2) becomes

[

]

dcss(x) + css(x)Gss dx

Jss ) -D

(3)

where Jss and css(x) are the steady-state flux and ion concentration at x. Gss is defined as

f(∆wo Φss - ∆wo Φ0) δ

Gss ) -

(4)

Direct integration of (3) between x ) 0 and x ) δ leads to

Jss ) DGss

cw - coeGssδ eGssδ - 1

(5)

From (3) and (5) the concentration distribution of the cation is in the interfacial region accordingly given by

css(x) )

cw(eGss(δ-x) - 1) - co(eGss(δ-x) - eGssδ) eGssδ - 1

(6)

Note that eq 6 leads to correct limits css(0) ) cw and css(δ) ) co. The system is perturbed by applying an ac potential difference to the interfacial region in addition to the steadystate value: ∆wo Φ ) ∆wo Φss + ∆wo Φac. This results in a perturbation of the steady-state cation distribution within the interfacial region, c(x,t) ) css(x) + cac(x,t), which leads to an ac current. In the present analysis, however, the boundary concentrations remain unchanged since no mass-transfer limitation outside the interfacial region is considered here.10 Hence,

Ion Transfer in Inhomogeneous Media

J. Phys. Chem. B, Vol. 101, No. 50, 1997 10803

cac(0,t) ) cac(δ,t) ) 0. In what follows, for the sake of simplicity, css(x) and cac(x,t) are denoted as css and cac, respectively. Let us rewrite eq 2, in a form that makes explicit reference to the deviation of the flux from its steady-state values, Jss, by the applied ac signal, Jac:

[dxd (c

J ) Jss + Jac ) -D

Hence, inserting (14) into (13), the impedance of the system can be calculated for different boundary conditions. (i) Equilibrium Impedance. At equilibrium, no steady-state current flows through the interface and Jss ) 0. Therefore, the steady-state concentration is equal to the equilibrium concentration, i.e., css ) ceq, which can be calculated from (3):

]

ceq(x) ) cwe-Geqx

(7)

where Geq corresponds to (4) in the absence of current. Taking, as discussed above, cac(0) ) cac(δ) ) 0, the solution of (10) is

+ cac) + (css + cac)(Gss + Gac)

ss

where Gac ) -f∆wo Φac/δ. From (3) and (7), and neglecting the second-order ac term, the ac component of the flux is given by

[

]

dcac -Jac ) D + cssGac + cacGss dx

(8)

The ac component of the concentration in the case of a sinusoidal potential perturbation, ∆wo Φac ) ∆wo Φ ˆ ejωt, can be w jω(t+θ) expressed as cac ) cˆe , where ∆o Φ ˆ is the peak potential perturbation, θ the phase angle with respect to the exciting potential function, and ω its angular frequency; j is the imaginary unit. As it will be seen, it is not necessary to calculate θ to obtain the impedance. The continuity condition of the flux is then given by

∂Jac ∂cac ) ) jωcac ∂x ∂t

(9)

-

cac ) D GssGaccw jω er1x + r2δ(1 - er1δ) - er1δ + r2x(1 - er2δ) er2δ - er1δ

∂cac jω ∂css c + G ) -G ss ac ac ∂x D ∂x ∂x2

(10)

where it has been assumed that ∂Gac /∂x ) ∂Gss /∂x ) 0 according to the Goldman constant field approximation. The total electric current is given by13a

(

)

∂2Φ I ) A FJ -  ∂x ∂t

Geq ( 2

r1,2 ) -

(∫

Iδ ) A F

∂ J dx +  ∆wo Φ ∂t

δ

0

∫0

δ

cac dx )

(12)

 w ∆o Φac ∫0δJac dx + jωδF

(13)

where Iac is the ac component of the current. From (8)

∫0 Jac dx ) -D[Gss∫0 cac dx + Gac∫0 css dx] δ

δ

δ

(14)

(17)

[

]

Gaccwδ2e-Geqδ/2 cosh(r1δ) - 1 cosh(r2δ) - 1 )sinh((r2 - r1)δ/2) (r1δ)2 (r2δ)2

Expanding cosh x ≈ 1 + x2/2 + x4/24, sinh x ≈ x + x3/6 and using the binomial expansion (1 + x)-1 ≈ 1 - x, after lengthy algebra (18) can be approximated to

GacGeqcwδ3e-Geqδ/2 12

∫0δcac dx ≈ -

(19)

Also, the second integral term on the right-hand side of (14) can be expanded to give

∫0

δ

(

(Geqδ)2 1 - e-Geqδ ) cwδe-Geqδ/2 1 + + ... Geqδ 24

ceq dx ) cwδ

)

(20) Therefore (see eq 14),

∫0δceq dx + Geq∫0δcac dx ≈

(

Gacδcwe-Geqδ/2 1 -

so that in the case of a sinusoidal perturbation it must be satisfied that

Iac 1 ) AF δ

Geq 2 jω + 2 D

DGacGeqcw r12(er2δ - 1)2er1δ - r22(er1δ - 1)2er2δ jω r1r2(r1 + r2)(er2δ - er1δ) (18)

Gac

)

x( )

From (16) and (17) and, after some rearrangements,

(11)

where the first term inside the parentheses is the current due to cation transfer (faradaic current) and the second term is the Maxwell’s displacement current (capacitive current). A is the interfacial area and  the mean permittivity of the interfacial region. Thus, even though J is position dependent, the total current I is position independent. Since no information of the time variation of the electric field inside the interfacial region can be obtained from our theoretical model, an alternative approach must be followed. Integration of (11) with respect to the spatial coordinate from x ) 0 to x ) δ gives

- e(r1+r2)x (16)

where r1,2 are the roots of the characteristic function

Applying (9) into (8) and after minor arrangements,

∂2cac

(15)

)

(Geqδ)2 + ... (21) 24

Finally, from (13), (14), and (21) it is found that

Iac )

[

]

A AF2Dcw -Geqδ/2 e + jω ∆wo Φac RTδ δ

(22)

The two contributions inside the brackets are admittance components corresponding to a parallel RC combination. In this case, the resistance corresponds to the charge-transfer

10804 J. Phys. Chem. B, Vol. 101, No. 50, 1997

Kontturi et al.

resistance, Rct,14 and the capacitance to that of the interfacial region of thickness δ. Therefore, at equilibrium

Rct )

RTδ f exp - (∆wo Φeq - ∆wo Φ0) 2 w 2 AF Dc

[

]

(23)

and the interfacial layer capacitance takes the form C ) A/δ;  ) r0, where r is the relative permittivity of the interfacial layer and 0 ) 8.854 × 10-14 F cm-1. (ii) Overpotential Impedance. In this case a steady-state current is flowing across the interface. Considering external potential polarization in terms of an overpotential, η, it is found that Gss ) Geq + Gov where Gov ) -η/δ ) -f(∆wo Φ ∆wo Φeq)/δ. In order to solve (14), css (eq. 6) has to be integrated across the interfacial region between 0 and δ. This gives

∫0

δ

given by

(cw - co)δ (cw - coeGssδ)δ css dx ) Gssδ eGssδ - 1

[

) cwδe-Geqδ/2

Iac w ∆o Φac

]

2 sinh(Geqδ/2) sinh(Govδ/2) + Gssδ sinh(Gssδ/2)

eq

[

(

( (

[

≈ cwδe-Geqδ/2 1 +

Geq(Geq - 2Gov)δ2 + ... 24

]

]

Rct )

(25)

The solution of (10) is here omitted, but the corresponding integral in (14) takes the form

G (cw - co)δ2

∫0δcac dx ) cosh(rac δ) - cosh(r δ) × 2

[

1

cosh(r1δ) - 1 2

(r1δ)

]

cosh(r2δ) - 1 -

(r2δ)2

(26)

Gac(cw - co)δ2 1 1 + ((r1δ)2 + (r2δ)2) + ... ≈12 20

[

[

]

The form of the admittance in this case is similar to that for the equilibrium case, as shown in (22), and therefore, the chargetransfer resistance is given by

) )

(Govδ)2 + ... 24 (Gssδ)2 + ... Gssδ 1 + 24

Govδ 1 +

AF2Dcwe-Geqδ/2 1 1 - Geq(Geq + 2Gov)δ2 + ... + RTδ 24 A (28) δ

)

(Geqδ)2 + ... Geqδ 1 + 24 + Gssδ

)



Expanding the sinh term, this can be approximated to

∫0δcss dx ≈ cwδe-G δ/2

Figure 2. Numerical simulation of the equilibrium impedance. r ) 40, δ/A ) 10-7 cm-1 and Dcw ) 10-15 mol cm-1 s-1. f ) 0.01-100 kHz, 5 points per decade. Numerical values indicate values of ∆wo Φeq - ∆wo Φ0 used.

]

where r1 and r2 correspond to (17) replacing Geq with Gss. Neglecting the last term inside the brackets, (26) simplifies to

Gaccwδ2e-Geqδ/2 2 sinh(Geqδ/2) ≈ 12 GacGeqcwδ3e-Geqδ/2 (27) 12

∫0δcac dx ≈ -

From (13) and (14) it follows that the interfacial admittance is

[

{

]}

f RTδ exp - (∆wo Φeq - ∆wo Φ0) × 2 AF2Dcw 1 1 + Geq(Geq + 2Gov)δ2 + ... (29) 24

[

]

Results and Discussion The approximations implicit in the series expansions made above may seem to restrict the present analysis to potentials close to the standard potential, while in practice the equilibrium potential can be varied in a range of several hundred millivolts by varying the concentration ratio of the potential-determining ion between the two phases. Therefore, the validity of the approximations made have been examined with numerical simulations. Computations were carried out with Matlab for Windows in a Pentium processor personal computer. A simulation of the equilibrium impedance is presented in Figure 2, using the complete forms of the integrals of (18) and (20). As can be seen, the impedance plot is similar to that of a parallel RC circuit. The potential dependence of Rct for this case is shown in Figure 3; Rct was obtained from the intercept of the impedance diagram with the real impedance axis at the lowfrequency limit. A similar plot for Rct in the overpotential case ii is shown in Figure 4, using the corresponding complete integrals given in (24) and (26). Equations 23 and 29 describing the charge-transfer resistance are important results. Firstly, since this is only a resistive component, its combination in parallel with the interfacial capacitance leads to a semicircle in the impedance plots indistinguishable from an activated charge-transfer reaction.14 Secondly, the potential dependence of Rct is identical to that predicted from the Butler-Volmer (B-V) formalism:14

) RB-V ct

RT exp[-Rf(∆wo Φeq - ∆wo Φ0)] 2 0 w AF k c

(30)

with a transfer coefficient R ) 1/2. The apparent charge-transfer

Ion Transfer in Inhomogeneous Media

J. Phys. Chem. B, Vol. 101, No. 50, 1997 10805

Figure 3. Potential dependence of the equilibrium charge-transfer resistance. Straight line corresponds to (23). Parameters are the same as in Figure 2.

Figure 4. Potential dependence of the charge-transfer resistance at the concentration ratios co/cw ) 0.01 (s) and 0.1 (- - -), corresponding to equilibrium potentials ∆wo Φeq - ∆wo Φ0 ) -118 and -59 mV, respectively. Parameters are the same as in Figure 2.

coefficient derived from the present work is given by

RT d ln Rct F d∆wΦ

R)-

o

(31)

eq

and from the potential dependence of the simulation results shown in Figure 3, a value of R ) 1/2 can be calculated at potentials close the equilibrium potential. This is strictly in agreement with the theoretical calculations presented in (23) and (29). It should be stressed, though, that the physical meaning of R from the Butler-Volmer and from the solution of the Nernst-Planck equations is quite different. In the former case, the transfer involves a single activated step, which could be, for instance, a potential dependent ionic resolvation process.17 However, in the Nernst-Planck case, R ) 1/2 occurs as a natural consequence of the form of the diffusion and migration equation combined with the symmetry of the Goldman approximation. The similarity between (23) and (30) holds as long as the equilibrium potential is close to the standard transfer potential. However, since the potential dependence of Rct is different in the two formalisms, the comparison is not straightforward at other potentials, while keeping the apparent standard rate constant defined as k0 ) D/δ. Although the two formalisms lead to identical relationships between the charge-transfer resistance and the potential, the physical models underlying the two approaches are very

different. As mentioned above, the idea behind the ButlerVolmer equation is that in the transition from water to the organic phase a potential dependent rate-determining step is present. In contrast, the treatment of ion transfer as a diffusion and migration problem does not make any assumptions regarding the individual steps occurring during transfer, but rather treats it as a continuous migration and diffusion through an interfacial region.9-12 The average diffusion coefficient for this process contains, in this case, the information on resolvation phenomena, through the incidence of an average friction coefficient that incorporates resolvation energy considerations. To give an idea of the order of magnitude of the quantities discussed, for k0 ) 0.22 cm s-1,15 (for example, for the tetraethylammonium transfer across the water-nitrobenzene interface), D ) 10-7 cm2 s-1 for δ ) 4.5 nm. The value of D estimated in the interfacial region for this value of k0 is lower than that in the adjoining phases, as expected if a resolvation process occurs during the transfer. The difference with the activated step approach is mainly in the possibility that resolvation occurs continuously during transfer. However, the recent results of Girault et al.8 for ion transfer across the water/ 1,2-dichloroethane interface give values of D much closer to those in the bulk phases. This gives additional support to the idea of treating the structure of these liquid/liquid interfaces as a mixed solvent layer.5 Of course, the value of its thickness, δ, is arbitrarily chosen, and future work is required to understand the dynamics of resolvation. A lead in this direction is the advanced molecular dynamics simulations of Benjamin et al.6 These nonequilibrium calculations for ion transfer indicate that surface roughness due to the phenomena of “fingering”6 and the associated transfer of water of hydration into the organic phase are responsible for the increased friction coefficient for ionic movement. This implies that the solvent exchange process occurs at a much faster rate than ionic migration through the thickness δ. Conclusions It has been shown that the ac impedance technique, commonly employed to obtain rate constants for charge-transfer reactions across immiscible interfaces, gives values that are indistinguishable for an activated step, or from a migration and diffusion model. The apparent charge-transfer coefficient takes the value of 1/2 at potentials close to the standard transfer potential. This shows that it is not necessary to interpret these rate constants as resulting from a single rate-determining step, but a resolvation model affecting the average diffusion coefficient (or friction coefficient) of the transferred ion across the interfacial region can account for the measured rate constants. Acknowledgment. This work has been carried out with support of the Training and Mobility of Researchers (TMR) of the European Union, Contract No. ERB-FMRX-CT96-0078. References and Notes (1) (a) Gavach, C.; Seta, P.; D’Epenoux, B. J. Electroanal. Chem. 1977, 83, 225. (b) Gros, M.; Gromb, C.; Gavach, C. J. Electroanal. Chem. 1978, 89, 29. (2) Reid, J.D.; Melroy, O. R.; Buck, R. O. J. Electroanal. Chem. 1983, 147, 71. (3) Samec, Z.; Marecˇek, V.; Homolka, D. J. Electroanal. Chem. 1985, 187, 31. (4) Tsionsky, M.; Bard, A. J.; Mirkin, M. V. J. Phys. Chem. 1996, 100, 17881. (5) (a) Girault, H. H.; Schiffrin, D.J. J. Electroanal. Chem. 1983, 190, 43. (b) Girault, H. H. In Modern Aspects of Electrochemistry; Bockris, J. O’M., Conway, B. E., White, R. E., Eds.; Plenum: New York, 1993; Vol. 25, p 1.

10806 J. Phys. Chem. B, Vol. 101, No. 50, 1997 (6) (a) Benjamin, I. J. Phys. Chem. 1991, 95, 6675. (b) Benjamin, I. J. Phys. Chem. 1992, 97, 1432. (c) Benjamin, I. Chem. ReV. 1996, 96, 1449. (d) Scheiwghofer, K. J.; Benjamin, I. J. Electroanal. Chem. 1995, 391, 1. (7) (a) Wandlowski, T.; Marecˇek, V.; Holub, K.; Samec, Z. J. Phys. Chem. 1989, 93, 8204. (b) Kakiuchi, T.; Noguchi, J.; Kotani, M.; Senda, M. J. Electroanal. Chem. 1990, 296, 517. (c) Shao, Y.; Campbell, J. A.; Girault, H. H. J. Electroanal. Chem. 1991, 300, 415. (d) Girault, H. H.; Schiffrin, D. J. J. Electroanal. Chem. 1985, 195, 213. (8) Beattie, P.D.; Delahay, A.; Girault, H.H. Electrochim. Acta 1995, 40, 2961. (9) Kakiuchi, T. J. Electroanal. Chem. 1992, 322, 55. (10) Kontturi, K.; Manzanares, J. A.; Murtoma¨ki, L. Electrochim. Acta 1995, 40, 2979. (11) Senda, M. Electrochim. Acta 1995, 40, 2993. (12) Aoki, K. Electrochim. Acta 1996, 41, 2321. (13) (a) Brumleve, T. R.; Buck, R. P. J. Electroanal. Chem. 1978, 90, 1. (b) Brumleve, T. R.; Buck, R. P. J. Electroanal. Chem. 1981, 126, 73.

Kontturi et al. (c) Buck, R. P.; Mundt, C. J. Chem. Soc., Faraday Trans. 1996, 92, 3947. (d) Buck, R. P.; Mundt, C. J. Chem. Soc., Faraday Trans. 1996, 92, 4987. (14) Bard, A. J.; Faulkner, L.R. Electrochemical Methods; John Wiley & Sons: New York, 1980; Chapter 9. (15) Lhotsky´, A.; Holub, K.; Neuzˇil, P.; Marecˇek, V. J. Chem. Soc., Faraday Trans. 1996, 96, 3851. (16) Haase, R. Thermodynamics of IrreVersible Processes; Addison-Wesley: Reading, 1969; Chapter 4. (17) Girault, H. H.; Schiffrin, D. J. J. Electroanal. Chem. 1985, 195, 213. (18) (a) Aguilella, V.; Belaya, M.; Levadny, V. Langmuir 1996, 12, 4817. (b) Taskinen, P.; Kontturi, K.; Sipila¨, A. Finn. Chem. Lett. 1980, 97. (c) Goldman, D. J. Gen. Physiol. 1943, 27, 37. (d) Hodgkin, A.; Katz, B. J. Physiol. 1949, 108, 37. (19) Pellicer, J.; Mafe´, S.; Aguilella, V. Ber. Bunsen-Ges. Phys. Chem. 1986, 90, 867.