Kinetic analysis based on

working surface interpolation

In order to deduce unknown values of kinetic parameters (such as rate constants and equilibrium constants) from experimental data, a 'fitting' (optimisation) procedure must be used. For a given trial set of parameters, the experimental response may be predicted by the model. In order to minimise the deviation between the experimental and theoretically-predicted response, the trial parameters must be adjusted and the process repeated iteratively.

The repeated computation of physical problems in more than one dimension can be very resource-intensive in terms of CPU time and memory^{1,2}. However, the partial differential equations describing many physical problems can be reduced to a dimensionless form so that the experimental observable may be represented as a sole function of a small number of dimensionless parameters which may be represented graphically as a working curve^{3} or surface^{4,5}. Once this surface has been generated using a powerful workstation or supercomputer, the simulated response may be reconstructed from the stored surface using an appropriate interpolation algorithm.

In this chapter, the merits of bilinear, cubic convolution, bicubic spline and Artificial Neural Network interpolation methods are assessed using the working surfaces generated in Chapter 6. Working surface interpolation is shown to be a viable cost-effective alternative to optimisation by repeated simulation, especially for problems with mass transport in more than one dimension.

By using the Internet, the surfaces need not even be distributed and stored on experimentalists' personal computers - the data can be stored on a central World Wide Web server which can compute the dimensionless parameters, perform the interpolation and return the results. In the last section of this chapter this is realised. Initially, the dimensionless parameters which may be used to describe an electrochemical reaction mechanism are summarised.

ElectrodeGeometry | Characteristictime scale (t _{c} = d^{2}/D) |

Spherical | |

Microdisc | |

Rotating disc^{*} | |

Wall-Jet^{*} | |

Micro-Jet^{*} | |

Channel^{*} |

Table 10.1 shows the characteristic time scale parameters for a number of common electrode geometries which reach a steady-state. The quadratic dependence of the time scale on the electrode radius is the reason why microelectrodes can be used to study significantly faster processes than electrodes of traditional (mm) dimensions. Note the analogy between the hydrodynamic electrodes (under first-order convection approximations) - the Peclet number, Pe, is the ratio of mass transport by convection to that by diffusion (defined symbolically for each geometry in section 10.3).

In Appendices 1-4, the dimensionless parameters on which the current depends are derived for a number of electrode geometries. These, summarised in Tables 10.2 - 10.5, define the working curve, surface or hypersurface that completely characterises a particular system.

Hemispherical/spherical electrode | - |

Microdisc Electrode | - |

Rotating disc electrode | Sc |

Rotating disc electrode (high Schmidt n^{o}) | - |

Wall-jet electrode 1 ^{st} order convection approximation | Pe |

Wall-jet electrode (no radial diffusion) 1 ^{st} order convection approximation | - |

Channel (microband) electrode | p_{1}, p_{2} |

Channel (microband) electrode Lévêque approximation | Pe |

Channel electrode (no axial diffusion) Lévêque approximation | - |

Table 10.3: Dimensionless parameters due to chemical species

Different diffusion coefficients (N species) [N-1 dimensionless parameters] | D_{B}/D_{A} ... D_{N}/D_{A} |

Reagent concentrations (N reagents) [N-1 dimensionless parameters] | [B]_{bulk}/[A]_{bulk} .. [N]_{bulk}/[A]_{bulk} |

Table 10.4: Dimensionless parameters due to kinetics

Reversible electron transfer | Q |

Quasi-reversible electron transfer | Q, K_{o} |

Irreversible homogeneous chemical reaction | K |

Irreversible heterogeneous chemical reaction | K_{het} |

Reversible homogeneous chemical reaction | K, K_{eq} |

Table 10.5: Dimensionless parameters due to time dependence

Chronoamperometry | t |

Linear-sweep/Cyclic voltammetry | s |

Dimensionless time and rate constants may be formulated using the time scale and average diffusion layer thickness:

K_{hom} = k.t_{c} = kd^{2}/D | (10.1) |

K_{het} = k(t_{c}/D)^{1/2} = kd/D | (10.2) |

K_{o} = k(t_{o}_{c}/D)^{1/2} = kd/D_{o} | (10.3) |

t = t/t_{c} = Dt/d^{2} | (10.4) |

(10.5) |

The dimensionless parameters in the tables above may be summed to predict on which dimensionless parameters the measured parameter (usually I_{lim}, N_{eff} or DE_{1/2}) depends. For example, consider the chronoamperometric study of a chemically irreversible, electrochemically reversible EC' reaction:

(10.6) |

at a wall-jet electrode (at low Peclet numbers where radial diffusion cannot be ignored), assuming the oxidised and reduced forms of the catalyst have equal diffusion coefficients which are different from that of the substrate. If the potential was not high enough to give the limiting current at steady-state, the dimensionless parameters would be:

t, Pe, K, D_{Z}/D_{A}, [Z]_{bulk}/[A]_{bulk}, Q | (10.7) |

Clearly such a high-dimensional hypersurface would require a large number of simulated values to define it. Suppose 20 points were enough to define each co-ordinate - 64 million points would be required! Fortunately, many experiments are conducted under conditions where there are fewer degrees of freedom. A more likely wall-jet experiment would be measuring steady-state limiting current at flow rates where radial diffusion effects were negligible. The current would then be a function of just 3 parameters:

K, D_{Z}/D_{A}, [Z]_{bulk}/[A]_{bulk} | (10.8) |

requiring a more acceptable 8000 points to be simulated and stored.

From Table 10.1, it is evident that the dimensionless parameters for hydrodynamic electrodes, under a 1^{st} order convection approximation, are all related to the (shear rate) Peclet number, Pe. Throughout this thesis, a shear-rate Peclet number, denoted P_{s}, has been used which is proportional, but not necessarily equal to, Pe. The constants of proportionality arise from normalisation by electrode diameter vs. radius, factors of p, empirical hydrodynamic constants etc. The definitions of P_{s} (and K) employed were chosen to maintain consistency with other definitions in the literature^{3,6}. Here these constants are summarised, allowing the characteristic timescale and all the associated parameters (10.1)-(10.5), to be calculated from the P_{s} values reported in this work. To begin, definitions of P_{s} are summarised in Table 10.6 together with P - the 'raw' Peclet number in terms of the experimental variables, with any constants of proportionality removed. This provides a 'baseline' from which to relate Pe and P_{s}. Note (comparing Table 10.1 with Table 10.6) that the time scale for the uniformly accessible RDE and MJE electrodes is independent of the electrode radius.

ElectrodeGeometry | Peclet number used | 'Raw' Peclet number |

Rotating disc^{a} | ||

Wall-Jet^{b} | ||

Micro-Jet^{c} | ||

Channel |

Next some of the equations from the previous section are re-written to allow for the constants of proportionality:

P_{s} = k_{p} P | Pe = k P_{s} | (10.9) |

d = k_{d} x_{e}P^{-1/3} | d = k_{dp} x_{e}P_{s-1/3} | (10.10) |

K = k_{k} (kx_{e2}/D) P^{-2/3} | K = k_{kp} (kx_{e}^{2}/D) P_{s-2/3} | (10.11) |

Nu = k_{n} P^{1/3} | Nu = k_{np} P_{s}^{1/3} | (10.12) |

Substituting P_{s} into the above expressions gives the relationships:

k_{np }= k_{n}k_{p-1/3} and k_{dp }= k_{d}k_{p}^{1/3} and k_{k} = k_{kp}k_{p}^{-2/3} | (10.13) |

And:

Nu = x_{e}/dNu = Pe ^{1/3} | k_{n}=1/k_{d}k _{ }= k_{np}^{3} | (10.14) |

Hence a 'conversion table' may be created for each electrode geometry, shown in Table 10.7.

Table 10.7: Constants for each electrode geometry

By equating the (average) diffusion layer thickness of various electrode geometries, a set of equivalence relationships can be generated^{3,6}. The 'equivalence' factors between the geometries are given in Table 10.8 and Table 10.9.

Geometry relationship | Peclet Number Difference | Value |

Channel/RDE | log_{10}Ps_{CHE} - log_{10}Ps_{RDE} | 0.854 |

Wall-jet/RDE | log_{10}Ps_{WJE} - log_{10}Ps_{RDE} | 0.060 |

Micro-jet/RDE | log_{10}Ps_{MJE} - log_{10}Ps_{RDE} | 0.932 |

Table 10.9: 'Equivalence' between dimensionless rate constants at various electrode geometries

Geometry relationship | Difference in K | Value |

Disc/Sphere | log_{10}K_{Disc} - log_{10}K_{sphere} | -3 log_{10}(4/p) |

RDE/Sphere | log_{10}K_{RDE} - log_{10}K_{sphere} | -0.384 |

ChE/Sphere | log_{10}K_{ChE} - log_{10}K_{sphere} | 0.185 |

WJE/Sphere | log_{10}K_{WJE} - log_{10}K_{sphere} | -0.310 |

MJE/Sphere | log_{10}K_{MJE} - log_{10}K_{sphere} | 0.237 |

Channel/RDE | log_{10}K_{CHE} - log_{10}K_{RDE} | 0.569 |

Wall-jet/RDE | log_{10}K_{WJE} - log_{10}K_{RDE} | 0.074 |

Micro-jet/RDE | log_{10}K_{MJE} - log_{10}K_{RDE} | 0.621 |

For non-uniformly accessible electrodes, the 'equivalence' is only approximate since it is based on an average diffusion layer thickness (therefore the deviation gives a measure of the degree of non-uniformity). For uniformly accessible hydrodynamic electrodes this equivalence is valid within the 1^{st}-order convection approximation. Since the effects of radial diffusion are negligible for the MJE above log_{10}P_{s} of 4.5 the RDE analogy, originally proposed by Albery and Bruckenstein^{10}, may be employed. This allows the analysis of steady-state voltammograms at the MJE using the working curves computed for the RDE in Chapter 7.

In order to generate comprehensive working curves/surfaces, simulations must be conducted over a wide range of conditions (both in terms of mass-transport and kinetics). If the values interpolated from the surface are to be accurate, a reasonably high resolution mesh of points is required to define the surface. Therefore the simulation methods used to generate such a surface should ideally be able to run without manual optimisation of the parameters at each point. Below is a summary of the key features of the automated simulations used to generate working curves/surfaces throughout this thesis.

- For each geometry a conformal mapping was used, ensuring good convergence with relatively few nodes at lower rate constants.
- At high rate constants the thinning of the reaction layer necessitated an increased number of nodes in the co-ordinate normal to the electrode surface.

For the channel microband electrode, the simulation geometry needs to be modulated as a function of Peclet number. At high Peclet numbers, the diffusion layer is constrained close to the electrode surface and there is little propagation of axial diffusion upstream or downstream of the electrode. The opposite is true at low Peclet numbers: a large space upstream and downstream is necessary to accommodate axial diffusion. The 2-D channel simulation is slow to converge due to the boundary singularity on the upstream edge of the electrode, so it is not possible simply to adopt the maximum space values and push up the number of nodes to remedy this problem.

To allow this simulation to run automatically, subroutines were written to 'search' for the edge of the diffusion layer. A coarse mesh problem was run first across the entire channel height. The resulting concentration profile was then searched, locating the point at which the normalised concentration departed from unity. A percentage error bound was added, and this was fixed as the fraction of the channel to be simulated. A similar search was conducted upstream of the electrode, starting from a maximum value of 200 electrode lengths. Since it was found in Chapter 6 that the upstream and downstream contributions from axial diffusion were approximately equal, it is possible to allow the same amount of propagation space downstream as upstream, thus only a single parameter needs to be optimised. This was also found to be optimal empirically.

For the wall-jet, two similar parameters require optimisation - the space above the electrode and the space downstream. The former is addressed by the conformal space, since the diffusion layer thickness scales with the h co-ordinate - hence a value of h_{max }= 1 appears to be suitable down to Ps=1. In the case of the space downstream, l_{d}, an alternative strategy was adopted. Simulations were manually optimised across a broad range of Peclet numbers (1 - 10^{10}). A simple quadratic fitted as a function of log(P_{s}) over the range 0-4 above which l_{d}=0.5 was sufficient, namely:

(10.15) |

so that intermediate values could be interpolated.

For a channel macroelectrode, the effects of axial diffusion can be ignored and the current may be expressed as a sole function of dimensionless kinetic parameters, as shown in Table 10.2. Therefore the working surfaces for common reaction mechanisms (such as those simulated in sections 6.2.1 and 6.2.2) collapse to working curves which may be simulated efficiently using a frontal strategy such as the BI method.

Compton et al.^{11-14} have published working curves for a number of common mechanisms using the BI method with a Cartesian mesh. Inevitably working curves published graphically are inconvenient for quantitative analysis, since a value must be interpolated for each measured current (the 'traditional' remedy to this problem is to use approximate functions to the curve or surface, based on regression^{15} or asymptotic formulae^{16}). These were therefore re-simulated using the BI method and stored in an electronic format at a high enough resolution (steps of 0.05 in log K) to allow accurate linear interpolation. In order to improve on the accuracy of the published data, and generate working curves over a wider range of rate constants, the following measures were taken:

- An exponentially expanding y-grid was used. This makes a significant difference to the convergence rate in this co-ordinate, especially at high rate constants where the reaction layer is very much thinner than the diffusion layer.
- Implicit kinetic terms were used where possible (e.g. for an ECE reaction). This improves the stability and accuracy at high rate constants (see section 3.7).
- For mechanisms where explicit kinetic terms could not be avoided, the simulation code contained checks for negative concentrations which result from breakdown of the explicit approximation. If negative concentrations were detected, the number of nodes in the x co-ordinate was doubled and the simulation was restarted. This allowed working curves to be generated up to rate constants where the amount of CPU time associated with a very large N
_{GX}value became prohibitive. - Where heterogeneous kinetics were simulated, back-to-back simulation grids
were used to facilitate coupling without resorting to an iterative scheme
(see section 3.6.1). This gave considerable
improvement in accuracy over the published working curves
^{17}(in which errors accumulated with increasing rate constant), and the simulation time was reduced from hours to seconds.

Simulations were conducted under conditions where the Lévêque approximation is valid (small electrode length of 0.1mm, high flow rate of 0.1cm^{3}s^{-1}) using a minimum of 1000x1000 nodes and rectangular integration.

The rate constants in both the EC and EC_{2} mechanisms may be characterised by measuring the shift in half-wave potential as a function of mass transport, then comparing this with a working curve of dimensionless shift vs. dimensionless rate constant. At high rate constants this behaviour limits to that predicted by simple reaction-layer theory - the shift in half-wave potential becomes a linear function of the log of the dimensionless rate constant. Reaction layer expressions have been derived for spherical^{18,19} and rotating disc electrodes (see Table 7.1) and are summarised Table 10.10.

The 'equivalence' factors, given in the previous section, allow reaction layer limits to be applied to other electrode geometries, though as shown in Chapter 9, the 'equivalence' relationship will not necessarily work for the full working curve, especially at lower rate constants. The transformations are applied to the EC and EC_{2} working curves in Figure 10.1 which were computed in the previous chapters. Note that the small difference between the WJE and RDE is about the same for both EC and EC_{2} processes - contradicting earlier work by Compton et al.^{6}, where a marked difference was observed for the EC_{2} process. One might expect second-order kinetics to pronounce the non-uniform concentration distribution at the electrode surface. This effect does not appear to be 'felt' in the shift in the half-wave potential. The effect on the limiting current (measured as N_{eff}) is discussed in the conclusions chapter.

Table 10.10: Reaction layer expressions for EC and EC_{2} reactions

Mechanism | Shift in dimensionless HWP |

EC | |

EC_{2} |

Due to the thin reaction layer at high rate constants, the (stiff) systems are computationally expensive to simulate accurately. Once the reaction layer limit has been reached, a linear extrapolation may be used to predict the behaviour at high rate constants. However one must be confident that one is not extrapolating using simulation data before the fully linear region, or that the spatial convergence is not slipping towards high rate constants (due to the thinning reaction layer). The safest way to perform the extrapolation is to use one of the reaction layer expressions from Table 10.10.

The working surfaces for E, ECE and EC_{2}E reactions at a (microband) channel electrode were used as test cases. These are representative of many of the log-scale working curves and surfaces found in voltammetry - the response curve is either sigmoidal or an increasing gradient up to a linear limit.

Interpolation was conducted on the decadic logarithm of the current so that the percentage error from interpolation would be approximately uniform across the surface. For ECE and EC_{2}E mechanisms, the interpolation was conducted on the N_{eff} value rather than on its logarithm, since N_{eff} spans a narrow range of magnitude (1-2).

Four methods were investigated. The first, bilinear interpolation (using the INTERPOLATE routine in IDL 4^{(r)20}) was chosen because it is the simplest and most easily programmed two-dimensional interpolation method. Interpolation using bicubic convolution was also assessed (using the /CUBIC option of the IDL INTERPOLATE routine).

Bicubic spline interpolation was conducted using the NAG library routines E02DAF to fit the spline surface and E02DBF to evaluate it at the sought co-ordinates. These routines were used with the parameters suggested in the NAG FORTRAN library manual^{21} to produce 'interpolating splines'. These give more accurate reconstruction of values near the ends of the range than a natural bicubic spline interpolation such as the routine E01ACF^{21} (where the second derivatives are zero at the ends of the ranges).

Figure 10.2: An Artificial Neural Network with an input layer of 2 nodes, 1 or more hidden layers of internal nodes, and an output layer with a single node.

The last method investigated was an Artificial Neural Network (ANN) based on sigmoidal transfer functions^{22}. A neural network with 2 inputs, 1 output and a number of internal nodes (shown in Figure 10.2) serves as a non-linear regression method in which the variable coefficients (connection weightings) are optimised for the data set via the simple back-propagating training procedure. For a fully-connected network (which was used in all cases) with 2 inputs and 1 output, the number of connections are given by:

3l_{1} | for one hidden layer |

2l_{1} + l_{1}.l_{2} + l_{2} | for two hidden layers |

2l_{1} + l_{2}(l_{1}+l_{3}) + l_{3} | for three hidden layers |

where l_{1}, l_{2} and l_{3} are the numbers of nodes in the first, second and third hidden layers, respectively.

The commercial package ARD Propagator(r) V1.0^{23} was used for ANN interpolation. In all cases a validation data set (the points which were to be interpolated) was used when training the neural network in order to prevent over-training. In order to be used with the neural network, the data was normalised to values suitable for the network transfer functions. Sigmoidal transfer functions were used, taking values between zero and 1, so the data was normalised to values between 0.05 and 0.95 before use with the ANN.

In order to evaluate the performance of various interpolation methods, a pair of surface data sets was simulated for each mechanism, the second (validation) surface being offset by half an interval in both the x and y directions from the first (training) data set, as shown in Table 10.11. The observable at the co-ordinates corresponding to the validation data set could then be interpolated using the training data set and the error measured. The parameters used in ARD Propagator^{(r)} for Neural Network fitting are shown in Table 10.12. Suitable learning rate and momentum factors were deduced empirically.

Surface | x | x_{min} | x_{max} | Dx | y | y_{min} | y_{max} | Dy |

ET training | log_{10}p_{1} | -2.4 | 2 | 0.2 | log_{10}p_{2} | 2 | 6 | 0.2 |

ET validation | log_{10}p_{1} | -2.3 | 1.9 | 0.2 | log_{10}p_{2} | 2.1 | 5.9 | 0.2 |

ECE training | log_{10}P_{s} | -4 | 4.4 | 0.4 | log_{10}K_{ECE} | -3 | 2 | 0.2 |

ECE validation | log_{10}P_{s} | -3.8 | 4.2 | 0.4 | log_{10}K_{ECE} | -2.9 | 1.9 | 0.2 |

EC_{2}E training | log_{10}P_{s} | -4 | 4.4 | 0.4 | log_{10}K_{EC2E} | -3 | 2 | 0.2 |

EC_{2}E validation | log_{10}P_{s} | -3.8 | 4.2 | 0.4 | log_{10}K_{EC2E} | -2.9 | 1.9 | 0.2 |

Table 10.12: Parameters for ARD Propagator(r) neural network

Parameter | Value |

Transfer functions | Sigmoidal |

Random seed | 0 |

Connectivity | Full |

Learning Rule | Generalised Delta |

Learning Rate | 0.4 |

Momentum Factor | 0.5 |

Training patterns order | Random |

Training Cycles | 100,000 |

Initial Weights | -1 to 1 |

Input Noise | none |

Figure 10.3: Number of connections in an Artificial Neural network with an equal number of nodes per layer (where the number of nodes was not exactly divisible by the number of layers, nodes were dropped to maximise the number of connections).

The interpolation errors for various network architectures are shown in Table 10.13. The first column corresponds to the number of nodes in each layer - two in the input layer and one in the output layer in all cases, but a variation in the number of hidden layers and nodes in the hidden layers. As one would expect, fitting improves and training time increases with the overall number of hidden-layer nodes as the number of variable parameters (weights) increases. Also perhaps as would be expected, the mean error decreases as the number of connections (and hence 'fitting' parameters) increases. The optimal architecture appears to be with two hidden layers as this gives rise to the maximum number of connections for a given number of nodes, as shown in Figure 10.3, though three hidden layers give a lower maximum error with a large number of nodes.

E | ECE | EC_{2}E | |||||

Architecture | Connections | Mean | Max | Mean | Max | Mean | Max |

2 10 1 | 30 | 0.554% | 3.485% | 0.119% | 0.405% | 0.036% | 0.187% |

2 5 5 1 | 40 | 0.371% | 2.464% | 0.067% | 0.297% | 0.023% | 0.162% |

2 3 4 3 1 | 33 | 0.433% | 4.171% | 0.082% | 0.404% | 0.034% | 0.192% |

2 20 1 | 60 | 0.549% | 3.457% | 0.099% | 0.464% | 0.034% | 0.197% |

2 10 10 1 | 130 | 0.371% | 2.200% | 0.049% | 0.273% | 0.022% | 0.161% |

2 6 7 7 1 | 111 | 0.370% | 2.414% | 0.048% | 0.311% | 0.019% | 0.172% |

2 40 1 | 120 | 0.664% | 3.811% | 0.105% | 0.419% | 0.044% | 0.158% |

2 20 20 1 | 460 | 0.316% | 2.18% | 0.0382% | 0.265% | 0.021% | 0.150% |

2 13 14 13 1 | 403 | 0.327% | 2.03% | 0.0429% | 0.236% | 0.021% | 0.142% |

Table 10.14 shows the percentage error in current for interpolation on each of the three surfaces using the three methods. The success of simple bilinear interpolation on the electron transfer surface can be rationalised by the planarity of the surface in the region where the Levich equation holds.

E | ECE | EC_{2}E | ||||

Interpolation method | Mean | Max | Mean | Max | Mean | Max |

Linear | 0.134% | 0.927% | 0.060% | 0.228% | 0.027% | 0.134% |

Cubic convolution | 0.525% | 5.925% | 0.066% | 0.465% | 0.030% | 0.279% |

Bicubic Spline | 0.043% | 0.894% | 0.012% | 0.225% | 0.007% | 0.134% |

Neural Network^{*} | 0.316% | 2.18% | 0.0382% | 0.265% | 0.021% | 0.142% |

Bicubic spline interpolation is clearly the superior method, though linear and ANN interpolation both give highly respectable results. For the ECE and EC_{2}E surface the 'maximum' error occurs at very low P_{s} (<-3.8) where there is a minor turning point in the data due to slightly imperfect convergence (maximum error of 3-4% in this region compared with <1% elsewhere). This is due to the extremely high number of nodes required to converge the Cartesian grid in the x direction under these conditions where there is a very large amount of axial diffusion relative to the size of the electrode. With 'perfect' (infinitely converged) data the results from bicubic spline and linear interpolation would be considerably better, however rather than rejecting data from this region, the data with a realistic imperfection was chosen evaluate the interpolation methods.

Working surfaces, together with an appropriate interpolation algorithm, allow experimental responses to be predicted accurately and instantaneously by the non-specialist. Of the interpolation methods investigated, bicubic spline interpolation was shown to be superior to both bilinear interpolation and Artificial Neural Networks. Cubic convolution interpolation performed relatively poorly in all cases. Bilinear interpolation was found to give highly satisfactory results for minimal programming effort and CPU time and is therefore the choice method for data sets of this resolution. For lower resolution surfaces bicubic spline interpolation would probably be the optimal method. For problems where the experimental observable must be represented by a higher-dimension (hyper) surface, the choice is more restricted - only linear interpolation and ANNs can readily be extended to higher dimensions.

The working curves and surfaces generated in this thesis form a mechanism 'library' for a number of common geometries. These were used as the basis of a data analysis program^{24,25}, available via the World Wide Web (WWW) at http://physchem.ox.ac.uk:8000/wwwda. The user simply selects a geometry and mechanism, enters the cell parameters along with the limiting current/half-wave potential as a function of radius/rotation speed/flow rate. Data may be entered directly via the WWW, or uploaded as a file by FTP. For each experimental data set, the theoretical response is predicted for the chosen mechanism and plotted as a function of mass transport appropriate to the electrode geometry.

In order to reduce the number of dimensionless parameters to a manageable number (Table 10.15):

- All species were assumed to have equal diffusion coefficients.
- For the EC' mechanism at the rotating disc electrode, a 1
^{st}order convection approximation was used to reduce the number of dimensionless parameters to 2. - For homogeneous processes at the channel electrode, the Lévêque approximation was used to reduce the number of dimensionless parameters to 2.

Mechanisms: | E | EC,EC_{2},ECE,EC_{2}EDISP1,DISP2 | EC' |

Spherical electrode | const. | K | K, r |

Microdisc electrode | const. | K | K, r |

Rotating disc electrode (1 ^{st} order convection) | const. | K | K, r |

Rotating disc electrode (9 ^{th} order convection) | Sc | Sc, K | Sc, K, r |

Wall-jet electrode (No radial diffusion) | const. | K | K, r |

Wall-jet electrode | P_{s} | P_{s}, K^{1} | Pr_{s}, K, |

Channel electrode (Lévêque approximation) | const. | K^{2} | K, r |

Channel (microband) electrode (Lévêque approximation) | P_{s} | P_{s}, K^{3} | Pr_{s}, K, |

Channel (microband) electrode (parabolic flow) | p_{1}, p_{2} | p_{1}, p_{2}, K | pr_{1}, p_{2}, K, |

A measure of the 'goodness of fit' is computed between the experimental and theoretical sets of data. The least squares error, given by^{26}:

(10.16) |

(where x denotes the observable) would seem a natural choice to measure the goodness of fit, but it is biased by outlying experimental data points. Therefore the more robust mean absolute deviation (MAD)^{26} is used:

(10.17) |

This is scaled when treating limiting currents or N_{eff} values, so that the relative error is minimised:

(10.18) |

in order not to bias the error by the magnitude of the current. For half-wave potentials, where the theoretical E_{1/2} value may be very close to zero, the absolute error is minimised.

Figure 10.4: Error surface fitting a rate constant and diffusion coefficient for an ECE mechanism. The experimental N

This error may be minimised to fit (i.e. optimise) unknown parameters such as rate constants or diffusion coefficients. A Golden Section Search^{26} is used to find a local minimum in one dimension. Downhill Simplex minimisation^{26} is used for multidimensional searches. One-, two- and multi-dimensional sampling routines are also available so that the error curve or surface may be visualised, for example see Figure 10.4, allowing the global minimum to be discerned from local minima. This service has been used as part of a number of mechanistic investigations to analyse data gathered at an array of channel electrodes, in order to distinguish between candidate mechanisms and fit diffusion coefficient and rate constant values to experimental data^{27-29}.

- 1 D. Britz, J. Electroanal. Chem., 406, (1996), 15.
- 2 M. Rudolph, in I. Rubinstein (Ed.) Physical Electrochemistry, Marcel Dekker, New York, (1995), p81.
- 3 P.R. Unwin, R.G. Compton, J. Electroanal. Chem., 245, (1988), 287.
- 4 J.A. Alden, R.G. Compton, R.A.W. Dryfe, J. Electroanal. Chem., 397, (1995), 11.
- 5 M.J. Bidwell, J.A. Alden, R.G. Compton, J. Electroanal. Chem, 414, (1996), 247.
- 6 R.G. Compton, A.C. Fisher, G.P. Tyley, J. Appl. Electrochem., 21, (1991), 295.
- 7 Yu.V. Pleskov, V.Yu. Filinovski, The rotating disc electrode, Plenum, New York, (1976), p17.
- 8 R.G. Compton, C.R. Greaves, A.M. Waller, J. Appl. Electrochem., 20, (1990), 575.
- 9 J.V.Macpherson, S.Marcar, P.R.Unwin, Anal.Chem., 66, (1994), 2175.
- 10 W.J. Albery, S. Bruckenstein, J. Electroanal. Chem., 144, (1983), 105.
- 11 R.G. Compton, M.B.G. Pilkington, G.M. Stearn, J.Chem.Soc. Faraday Trans. 1, 84, (1988), 2155.
- 12 R.G. Compton, M.B.G. Pilkington, J. Chem. Soc. Faraday Trans. I, 85, (1989), 2255.
- 13 R.G. Compton, J.C. Eklund, L. Nei, A.M. Bond, R. Colton, Y.A. Mah, J. Electroanal. Chem., 385, (1995), 249.
- 14 J.A. Cooper, R.G. Compton, Electroanalysis, 10, (1998), 141.
- 15 J.A. Alden, R.G. Compton, J. Electroanal. Chem., 404, (1996), 27.
- 16 W.M. Leslie, J.A. Alden, R.G. Compton, T. Silk, J. Phys. Chem., 100, (1996), 14130.
- 17 R.G. Compton, J.C. Eklund, L. Nei, A.M. Bond, R. Colton, Y.A. Mah, J. Electroanal. Chem., 385, (1995), 249.
- 18 Q. Zhuang, H. Chen, J. Electroanal. Chem., 346, (1993), 29.
- 19 Q. Zhuang, D. Sun, J. Electroanal. Chem., 440, (1997), 103.
- 20 Interactive Data Language, Floating point systems; Berks, UK. http://www.floating.co.uk/idl/index.html.
- 21 NAG FORTRAN Library, E01 Introduction, section 3.3.1.
- 22 H.M. Cartwright, Applications of Artificial Intellegence in Chemistry, Oxford University Press, (1993).
- 23 ARD Corporation; Columbia, USA. See http://www.cs.umbc.edu/~zwa/Gator/Description.html.
- 24 J.A Alden, R.G. Compton, Electroanalysis, 10, (1998), 207.
- 25 Announcement: J. Electroanal. Chem., 443, (1998), 227.
- 26 W.H. Press, S.A. Teukolsky, W.T. Vetterling and B.P. Flannery, Numerical Recipes in C: The Art of Scientific Computing, 2nd Edition, Cambridge University Press, (1992).
- 27 J.A .Alden, M.A. Feldman, E. Hill, F. Prieto, M. Oyama, B.A. Coles and R.G. Compton, Anal. Chem., 70, (1998), 1707.
- 28 F. Prieto, M. Oyama, B.A. Coles, J.A. Alden, R.G. Compton, S.Okazaki, J. Electroanal. Chem., (1998), in press.
- 29 J.A Cooper, R.G. Compton, Electroanalysis, submitted.