CHEM465/865, 2004-3, Lecture 20, 27 th Sep., 2004

Hydrodynamic Electrodes and Microelectrodes

So far we have been considering processes at planar electrodes. We have focused on the interplay of diffusion and kinetics (i.e. charge transfer as described for instance by the different formulations of the Butler-Volmer equation). In most cases, diffusion is the most significant transport limitation. Diffusion limitations arise inevitably, since any reaction consumes reactant molecules. This consumption depletes reactant (the so-called electroactive species) in the vicinity of the electrode, which leads to a non-uniform distribution (see the previous notes).

______Note: In principle, we would have to consider the accumulation of product species in the vicinity of the electrode as well. This would not change the basic phenomenology, i.e. the interplay between kinetics and transport would remain the same. But it would make the mathematical formalism considerably more complicated. In order to simplify things, we, thus, focus entirely on the reactant distribution, as the species being consumed. ______

In this part, we are considering a semiinfinite system: The planar electrode is assumed to have a huge surface area and the solution is considered to be an infinite reservoir of reactant. This simple system has only one characteristic length scale: the thickness of the diffusion layer (or mean free path) δδδ. Sometimes the diffusion layer is referred to as the “Nernst layer” . Now: let’s consider again the interplay of kinetics and diffusion limitations. Kinetic limitations are represented by the rate constant k 0 (or equivalently by the 0=== 0bα b 1 −−− α exchange current density j nFkcred c ox ). Diffusion limitations are represented by the diffusion constant D and by the diffusion layer thickness δδδ . D We can define a diffusion rate in the following way: k === . The corresponding diff δδδ

nFDc b diffusion-limited current is j= nFk c b = ox . diff diff ox δδδ

0 The two rates, k and kdiff , determine the interplay between kinetics and 0 >>>>>> mass transport. Reactions for which k k diff are called reversible reactions. 0 <<<<<< Reactions for which k k diff are called irreversible. In the chapter on cyclic we will consider this distinction in more detail. The rates k 0 can vary over wide ranges (from 10 1 cm/s for facile reactions down to 10 -14 cm/s for rather slow reactions). They are determined by the electronic structure of the metal (the Fermi level), the structure of the solution and the LUMO (lowest unoccupied molecular orbital) for reduction or the HOMO (highest occupied molecular orbital) for oxidation of the species in solution undergoing reduction or oxidation. In general, k 0 depends on the type of metal, its surface structure, the type of electrolyte and the redox species. The values of diffusion coefficients in aqueous solution are usually in the range of 10 -5 cm 2/s, with only a small range of variation. This parameter cannot be controlled in an experiment. In this part we will learn ways to control the interplay between kinetics and mass transport. Which options exist to influence this interplay? As you know already, the rate of electron transfer in an electrochemical system is controlled by the electrode potential E. This is the most important variable in an electrochemical experiment. We have studied its effects already in detail. What are the other options of experimental control? Subsequently, experimental techniques will be discussed, which allow to control the rates of mass transport, i.e. allow to control kdiff . There are two principal ways to achieve that: Hydrodynamic devices – forced convection. They help to confine concentration variations to a thin region near the electrode surface. Control the electrode geometry, i.e. use (ultra-)microelectrodes .

Overall, these measures raise the rates of mass transport. Fast rates of mass transport make it possible to study the kinetics of fast electron transfer reactions. Hydrodynamic Devices These devices use convection to enhance and control the rate of mass transport to the electrode surface. Detectable currents are increased and the sensitivity of voltammetric measurements is enhanced.

Two approaches are possible: The electrode is held in a fixed position and solution is flowed over the electrode surface by an applied force, usually an applied pressure gradient (e.g. wall-jet electrode) The electrode is designed to move which acts to mix the solution via convection.

In order to be able to perform a quantitative analysis of the electrode processes, the introduced convection must be predictable. The flow of the solution must be laminar rather than turbulent in order to lead to well-defined, reproducible results. The figure below compares turbulent and laminar flow.

Dropping Mercury Electrode (DME) Historically, this is the first used hydrodynamic technique. The electrochemical cell with potentiostat, (mercury), counter electrode and reference electrode is shown below.

mercury drop electrode surface

A large reservoir of mercury is connected to a capillary. In the capillary, mercury flows under the influence of gravitation. The drop at the opening of the capillary grows in time until it reaches a critical size. At some point, the mercury drop from detaches from the tip and falls down. The surface of the working electrode, which is the surface of the drop, is thus refreshed in regular cycles. A big advantage of this electrode is, that the continuous refreshing minimizes problems of electrode poisoning. Clearly, the measured current at will be a function of the surface area of the drop. It increases continuously with drop size. When the drop falls of, the current drops rapidly. The following picture shows the cyclic current at a DME as a function of time (in this plot: for a fixed potential).

Evidently, the surface area is an important property of this electrode. We, thus, have to consider current and NOT current densities at this electrode. Consider the limiting current due to diffusion as determined by COTTRELL- equation, D I((()( t))) === nFA c b L Hg ox πππ t Assume that the surface of the drop is an ideal sphere (this is of course a === πππ 2 simplification!), i.e. AHg 4 r 0 , where r0 is the radius of the drop. We assume a

-1 constant mass flux mHg [mg s ] of liquid Hg in the capillary. This mass flux determines the time-dependence of the drop radius,

1/ 3 3m t  r === Hg  0 πρπρπρ  4 Hg  ρρρ where Hg is the density of Hg. An effective diffusion coefficient has to be used in 7 the Cottrell equation, D=== D . Using all these definitions in the Cottrell- eff 3 equation, the so-called ILKOVIC-equation is obtained, ((( ))) === 1/2b 2/31/6 ItL708 nDcmt ox Hg

2 where IL is in , D in cm /s, mHg in mg/s, and t in s. The dropping mercury electrode can be used for voltammetric measurements, i.e. measure the current as a function of applied electrode potential E. For historic reasons, this method is called . It was first introduced by Heyrovsky in the 1920’s. The following picture shows a linear scan polarograph (obtained by ) for two reactions:

Curve A: Cd 2+ + 2 e - + Hg  Cd(Hg) (reduction of Cd ions) + - Curve B: 2H + 2 e  H2(g) (hydrogen evolution)

Diffusion current

IL

Note: The spikes visible in this plot represent the cycles of the drop lifetime. The drop lifetime is a constant (determined by the height of the mercury column in the reservoir and by the mass flow rate mHg ). At large cathodic overpotentials, the current reaches a plateau. The height of this plateau is determined by the Ilkovic equation, as specified above. In principle, this electrode is operated at steady state.

Rotating Disc Electrode (RDE) The rotating disc electrode is the most widely used is the first used hydrodynamic electrode. A disc electrode is embedded into the bottom face of an insulating rod (e.g. Teflon). The rod rotates at a constant angular velocity ωωω . The rotation drags solution to the electrode surface, resulting in a vortex, as shown below. Due to this drag and the steady laminar flow to the electrode surface, solution is continuously replaced. This electrode can be operated at steady state.

At angular velocities ωωω <<< 60 s -1 the flow profile will be laminar. Moreover, the electrode disc is considered small compared to the surface area of the insulating rod. This provides uniform conditions at the electrode surface. The figure below shows potential sweep voltammograms (for a cathodic reaction, i.e. negative current, negative electrode potential) measured at an RDE at various frequencies.

Obviously, the current in the mass transfer limited region (large | E|) is independent of time, but it is controlled by the rotation speed. In order to understand this behaviour, we have to return to the concept of the diffusion layer.

What controls the thickness of the diffusion layer?

The transport equations are given by the following modified form of Fick’s equations: vvv = − ∇ molar flux: J vcox D c ox ∂∂∂c vvv t-variation: ox =−∇vcb + Dc ∆ ∂∂∂t ox ox convection + diffusion

The results from linear sweep voltammetry indicate, that δδδ is now controlled by convection and not by diffusion! Stationary operation is possible.

∂∂∂c Consider the steady state limiting current density, i.e. consider the case ox === 0 ∂∂∂t The theory for this electrode was developed by Levich: Levich theory ! Details of this theoretical solution will be skipped here. They can be found for instance in the book: – Principles, Methods and Applications , C.M.A. Brett, A.M.O. Brett, Oxford University Press, Oxford, 1993, section 5.9). The thickness of the diffusion layer is given by:

1/3 1/6−−− 1/2 δ=== 1.61 D ν ω where ν is the kinematic viscosity [cm 2 s -1]. The mass transport limited current density is, thus, given as a function of the rotation speed as

cb j=== nFD ox L δδδ

=== b 2/3ν−−− 1/6 ω 1/2 0.602 nFcox D

ωωω 1/ 2 A plot of jL vs. will , thus result in a straight line.

The slope of this straight line is determined by the bulk concentration of reactant, by the diffusion coefficient and by the kinematic viscosity, i.e. === b 2/3ννν −−− 1/6 slope0.602 nFcox D .

Some notes on using the rotating disc electrode: The angular velocity ωωω should be small enough so that the flow profile will be laminar. The so-called Reynolds number Re, a characteristic of the flow profile, has to be smaller than a critical value. The Reynolds number is defined by vl Re === , ννν where v is a characteristic velocity of the fluid relative to the solid surface of the electrode, l is a characteristic length of the electrode, and ν is the kinematic viscosity (a measure of the inner friction within the fluid). For an RDE, v is the linear velocity at the outer edge of the disc electrode, given by v =ωωω ⋅ r and l=== r is the radius of the disc. The critical Reynolds- 5 Number is Re crit =10 . The criterion for laminar flow is, thus ωωωr⋅⋅⋅ r ωωω A Re= =disc << 10 5 ν πν Kinematic viscosities for dilute aqueous solutions are typically in the range

− ννν === 10 −−2 cm2 s -1 and we obtain the criterion ωωω << ⋅ 3 2 -1 Adisc 3 10 cm s . Consider for example the following disc areas and corresponding upper limits on angular velocity: A === 1 cm2 ⇒⇒⇒ ωωω << 3 ⋅ 103 s -1 , disc === 2⇒⇒⇒ ωωω << ⋅ 5 -1 Adisc 1 mm 3 10 s Smaller electrode surface area – operation at larger angular velicities possible. On the other hand, the angular velocity controls the thickness of the diffusion layer (see expression for δδδ in Levich theory), which controls the limiting current density due to diffusion. The larger ωωω is, the smaller is δ ωωω thinner is and the is jL . Therefore, should not be too small. The typical frequency range for operating an RDE with 0.3 cm radius is 10 s-1 << ωωω << 310 ⋅ 3 s -1 The electrode and rod need to be perfect cylinders to avoid wobbling around the axis. This is not a problem for mm sized electrodes. It could be a problem for smaller electrodes.

Capabilities of RDEs: ηηη Operation in the diffusion limited regime (at large ): From a plot of jL

vs. ωωω 1/ 2 it is possible to determine the diffusion coefficient of reactant species in the solution or determine the concentration of reactant. Operation in the kinetic or mixed regime (small to medium ηηη ): The total

current density j is determined by the following relation, involving an

activation controlled current density jac and the diffusion-limited current

density jdiff

11= + 1 = 1 + 1 ωωω1/ 2 jjac j diff j ac B This can be thought of as a series equivalent circuit of mass transport and activation barrier that electroactive species have to overcome in order to react. On the right hand side, the frequency dependence of the diffusion

1 − limited current has been inserted. It is evident, that a plot of over ωωω −−1/ 2 j (at a fixed potential) will result in a straight line. The intercept with the ωωω → ∞ ordinate, i.e. the limit will give jac , as shown in the following plot.

Summary: Three types of effects that you have to be able to distinguish Potential E Time scales, rates of processes Geometry, length scales! The latter two are interrelated!

Voltammetry: electrochemistry techniques based on current measurement as a function of applied electrode potential – what do you need for that? Electrochemical cell!