Search

Alternating current (AC) voltammetric techniques are experimentally powerful as they enable Faradaic current to be isolated from non-Faradaic contributions. Finding the best global fit between experimental voltammetric data and simulations based on reaction models requires searching a substantial parameter space at high resolution. In this paper, we estimate parameters from purely sinusoidal voltammetry (PSV) experiments, investigating the redox reactions of a surface-confined ferrocene derivative. The advantage of PSV is that a complete experiment can be simulated relatively rapidly, compared to other AC voltammetric techniques. In one example involving thermodynamic dispersion, a PSV parameter inference effort requiring 7,500,000 simulations was completed in 7 h, whereas the same process for our previously used technique, ramped Fourier transform AC voltammetry (ramped FTACV), would have taken 4 days. Using both synthetic and experimental data with a surface confined diazonium substituted ferrocene derivative, it is shown that the PSV technique can be used to recover the key chemical and physical parameters. By applying techniques from Bayesian inference and Markov chain Monte Carlo methods, the confidence, distribution, and degree of correlation of the recovered parameters was visualized and quantified.

The redox chemistry of the electron entry/exit site in Escherichia coli hydrogenase-1 is shown to play a vital role in tuning biocatalysis. Inspired by nature, we generate a HyaA-R193L variant to disrupt a proposed Arg-His cation-π interaction in the secondary coordination sphere of the outermost, "distal", iron-sulfur cluster. This rewires the enzyme, enhancing the relative rate of H 2 production and the thermodynamic efficiency of H 2 oxidation catalysis. On the basis of Fourier transformed alternating current voltammetry measurements, we relate these changes in catalysis to a shift in the distal [Fe 4 S 4 ] 2+/1+ redox potential, a previously experimentally inaccessible parameter. Thus, metalloenzyme chemistry is shown to be tuned by the second coordination sphere of an electron transfer site distant from the catalytic center.

Rapid disulfide bond formation and cleavage is an essential mechanism of life. Using large amplitude Fourier transformed alternating current voltammetry (FTacV) we have measured previously uncharacterized disulfide bond redox chemistry in Escherichia coli HypD. This protein is representative of a class of assembly proteins that play an essential role in the biosynthesis of the active site of [NiFe]-hydrogenases, a family of H 2 -activating enzymes. Compared to conventional electrochemical methods, the advantages of the FTacV technique are the high resolution of the faradaic signal in the higher order harmonics and the fact that a single electrochemical experiment contains all the data needed to estimate the (very fast) electron transfer rates (both rate constants ≥ 4000 s -1 ) and quantify the energetics of the cysteine disulfide redox-reaction (reversible potentials for both processes approximately -0.21 ± 0.01 V vs SHE at pH 6). Previously, deriving such data depended on an inefficient manual trial-and-error approach to simulation. As a highly advantageous alternative, we describe herein an automated multiparameter data optimization analysis strategy where the simulated and experimental faradaic current data are compared for both the real and imaginary components in each of the 4th to 12th harmonics after quantifying the charging current data using the time-domain response.

Estimation of thermodynamic and kinetic parameters in electrochemical studies is usually undertaken via comparison of the experimental results with theory based on a model that mimics the experiment. The present study examines the credibility of transient d.c. and a.c. voltammetric theory-experiment comparisons for recovery of the parameters needed to model the ubiquitous mechanism when an electron transfer (E) reaction is followed by a chemical (C) step in the EC process ([Formula: see text]). The data analysis has been undertaken using optimization methods facilitated in some cases by grid computing. These techniques have been applied to the simulated (5% noise added) and experimental (reduction of trans-stilbene) voltammograms to assess the capabilities of parameter recovery of E(0) (reversible potential for the E step), k(0) (heterogeneous electron transfer rate constant at E(0)), α (charge transfer coefficient for the E step), and k(f) and k(b) (forward and backward rate constants for the C step) under different kinetic regimes. The advantages provided by the use of a.c. instead of d.c. voltammetry and data optimization methods over heuristic approaches to "experiment"-theory comparisons are discussed, as are the limitations in the efficient recovery of a unique set of parameters for the EC mechanism. In the particular experimental case examined herein, results for the protonation of the electrochemically generated stilbene dianion demonstrate that, notwithstanding significant advances in experiment and theory of voltammetric analysis, reliable recovery of the parameters for the EC mechanism with a fast chemical process remains a stiff problem.

Commonly, significant discrepancies are reported in theoretical and experimental comparisons of dc voltammograms derived from a monolayer or close to monolayer coverage of redox-active surface-confined molecules. For example, broader-than-predicted voltammetric wave shapes are attributed to the thermodynamic or kinetic dispersion derived from distributions in reversible potentials (E(0)) and electrode kinetics (k(0)), respectively. The recent availability of experimentally estimated distributions of E(0) and k(0) values derived from the analysis of data for small numbers of surface-confined modified azurin metalloprotein molecules now allows more realistic modeling to be undertaken, assuming the same distributions apply under conditions of high surface coverage relevant to voltammetric experiments. In this work, modeling based on conventional and stochastic kinetic theory is considered, and the computationally far more efficient conventional model is shown to be equivalent to the stochastic one when large numbers of molecules are present. Perhaps unexpectedly, when experimentally determined distributions of E(0) and k(0) are input into the model, thermodynamic dispersion is found to be unimportant and only kinetic dispersion contributes significantly to the broadening of dc voltammograms. Simulations of ac voltammetric experiments lead to the conclusion that the ac method, particularly when the analysis of kinetically very sensitive higher-order harmonics is undertaken, are far more sensitive to kinetic dispersion than the dc method. ac methods are therefore concluded to provide a potentially superior strategy for addressing the inverse problem of determining the k(0) distribution that could give rise to the apparent anomalies in surface-confined voltammetry.

Many electrode processes that approach the "reversible" (infinitely fast) limit under voltammetric conditions have been inappropriately analyzed by comparison of experimental data and theory derived from the "quasi-reversible" model. Simulations based on "reversible" and "quasi-reversible" models have been fitted to an extensive series of a.c. voltammetric experiments undertaken at macrodisk glassy carbon (GC) electrodes for oxidation of ferrocene (Fc(0/+)) in CH3CN (0.10 M (n-Bu)4NPF6) and reduction of [Ru(NH3)6](3+) and [Fe(CN)6](3-) in 1 M KCl aqueous electrolyte. The confidence with which parameters such as standard formal potential (E(0)), heterogeneous electron transfer rate constant at E(0) (k(0)), charge transfer coefficient (α), uncompensated resistance (Ru), and double layer capacitance (CDL) can be reported using the "quasi-reversible" model has been assessed using bootstrapping and parameter sweep (contour plot) techniques. Underparameterization, such as that which occurs when modeling CDL with a potential independent value, results in a less than optimal level of experiment-theory agreement. Overparameterization may improve the agreement but easily results in generation of physically meaningful but incorrect values of the recovered parameters, as is the case with the very fast Fc(0/+) and [Ru(NH3)6](3+/2+) processes. In summary, for fast electrode kinetics approaching the "reversible" limit, it is recommended that the "reversible" model be used for theory-experiment comparisons with only E(0), Ru, and CDL being quantified and a lower limit of k(0) being reported; e.g., k(0) ≥ 9 cm s(-1) for the Fc(0/+) process.

Fully automated and computer assisted heuristic data analysis approaches have been applied to a series of AC voltammetric experiments undertaken on the [Fe(CN)6](3-/4-) process at a glassy carbon electrode in 3 M KCl aqueous electrolyte. The recovered parameters in all forms of data analysis encompass E(0) (reversible potential), k(0) (heterogeneous charge transfer rate constant at E(0)), α (charge transfer coefficient), Ru (uncompensated resistance), and Cdl (double layer capacitance). The automated method of analysis employed time domain optimization and Bayesian statistics. This and all other methods assumed the Butler-Volmer model applies for electron transfer kinetics, planar diffusion for mass transport, Ohm's Law for Ru, and a potential-independent Cdl model. Heuristic approaches utilize combinations of Fourier Transform filtering, sensitivity analysis, and simplex-based forms of optimization applied to resolved AC harmonics and rely on experimenter experience to assist in experiment-theory comparisons. Remarkable consistency of parameter evaluation was achieved, although the fully automated time domain method provided consistently higher α values than those based on frequency domain data analysis. The origin of this difference is that the implemented fully automated method requires a perfect model for the double layer capacitance. In contrast, the importance of imperfections in the double layer model is minimized when analysis is performed in the frequency domain. Substantial variation in k(0) values was found by analysis of the 10 data sets for this highly surface-sensitive pathologically variable [Fe(CN)6](3-/4-) process, but remarkably, all fit the quasi-reversible model satisfactorily.

A detailed analysis of the cooperative two-electron transfer of surface-confined cytochrome c peroxidase (CcP) in contact with pH 6.0 phosphate buffer solution has been undertaken. This investigation is prompted by the prospect of achieving a richer understanding of this biologically important system via the employment of kinetically sensitive, but background devoid, higher harmonic components available in the large-amplitude Fourier transform ac voltammetric method. Data obtained from the conventional dc cyclic voltammetric method are also provided for comparison. Theoretical considerations based on both ac and dc approaches are presented for cases where reversible or quasi-reversible cooperative two-electron transfer involves variation in the separation of their reversible potentials, including potential inversion (as described previously for solution phase studies), and reversibility of the electrode processes. Comparison is also made with respect to the case of a simultaneous two-electron transfer process that is unlikely to occur in the physiological situation. Theoretical analysis confirms that the ac higher harmonic components provide greater sensitivity to the various mechanistic nuances that can arise in two-electron surface-confined processes. Experimentally, the ac perturbation with amplitude and frequency of 200 mV and 3.88 Hz, respectively, was employed to detect the electron transfer when CcP is confined to the surface of a graphite electrode. Simulations based on cooperative two-electron transfer with the employment of reversible potentials of 0.745 ± 0.010 V, heterogeneous electron transfer rate constants of between 3 and 10 s(-1) and charge transfer coefficients of 0.5 for both processes fitted experimental data for the fifth to eighth ac harmonics. Imperfections in theory-experiment comparison are consistent with kinetic and thermodynamic dispersion and other nonidealities not included in the theory used to model the voltammetry of surface-confined CcP.

In this article, we review some of our previous work that considers the general problem of numerical simulation of the currents at microelectrodes using an adaptive finite element approach. Microelectrodes typically consist of an electrode embedded (or recessed) in an insulating material. For all such electrodes, numerical simulation is made difficult by the presence of a boundary singularity at the electrode edge (where the electrode meets the insulator), manifested by the large increase in the current density at this point, often referred to as the edge effect. Our approach to overcoming this problem has involved the derivation of an a posteriori bound on the error in the numerical approximation for the current that can be used to drive an adaptive mesh-generation algorithm, allowing calculation of the quantity of interest (the current) to within a prescribed tolerance. We illustrate the generic applicability of the approach by considering a broad range of steady-state applications of the technique.

Large-amplitude sinusoidal ac voltammetric techniques, when analyzed in the frequency domain using the Fourier transform-inverse Fourier transform sequence, produce the expected dc and fundamental harmonic ac responses in addition to very substantial second, third, and higher ac harmonics that arise from the presence of significant nonlinearity. A full numerical simulation of the process, Red right arrow over left arrow Ox + e(-), incorporates terms for the uncompensated resistance (R(u)), capacitance of the double layer (C(dl)), and slow electron transfer kinetics (in particular, the reversible potential (E degrees ), rate constant (k(0)), and charge transfer coefficient (alpha) from the Butler-Volmer model). Identification of intuitively obvious patterns of behavior (with characteristically different sensitivity regimes) in dc, fundamental, and higher harmonic terms enables simple protocols to be developed to estimate R(u), C(dl), E degrees , k(0), and alpha. Thus, if large-amplitude sinusoidal cyclic voltammograms are obtained for two concentrations of the reduced species, data obtained from analysis of the recovered signals provide initial estimates of parameters as follows: (a) the dc cyclic component provides an estimate of E degrees (because the R(u) and k(0) effects are minimized); (b) the fundamental harmonic provides an estimate of C(dl) (because it has a high capacitance-to-faradaic current ratio); and (c) the second harmonic provides an estimate of R(u), k(0), and alpha (because the C(dl) effect is minimized). Methods of refining the initial estimates are then implemented. As a check on the fidelity of the parameters (estimated on the basis of an essentially heuristic approach that solely utilizes the dc, fundamental, and second harmonic voltammograms), comparison of the predicted simulated and experimental third (or higher) harmonic voltammograms can be made to verify that agreement between theory and experiment has been achieved at a predetermined level. The use of the heuristic pattern recognition approach to evaluate the oxidation of ferrocene at a platinum electrode (a reversible process) in the very high resistance solvent dichloromethane (0.1 M Bu(4)NPF(6)) and the reduction of [Fe(CN(6))](3)(-) at a glassy carbon electrode (a quasi-reversible process) in much lower resistance but higher capacitance conditions found in aqueous (0.5 M KCl) media is described and verifies the inherent advantages of employing large-amplitude sinusoidal techniques in quantitative studies of electrode processes.