Digital simulation of chronopotentiometric and steady-state voltammetric curves at microelectrodes in the presence of a low concentration of supporting electrolyte

A simulation scheme for the calculation of theoretical chronopotentiograms at microelectrodes in solutions containing low amounts of supporting electrolyte is presented. The scheme allows computation of the changes in the concentration profiles of the substrates, products and the supporting electrolyte ions with time. The electrode potentials that are established after reaching the steady-state, together with the appropriate current intensities, can be used for constructing the steady-state voltammograms. The simulation of the mixed diffusional and migrational transport is based on the Crank-Nicolson method with an exponentially expanding time and space grids. The scheme does not impose any limitations on diffusion coefficients and it can be applied both to simple electrode reactions (one reactant-one product) and more complicated reactions under the assumption that the double-layer thickness is small in comparison to the diffusion layer. Five simple types of electrode reactions and an example of a more complicated scheme were considered. The results obtained demonstrate that the dependence of the steady-state limiting current on the support ratio (csupp.el./csubst) depends not only on the charge of the reactant and the product, but also on the diffusion coefficient ratio of the substrate and product. If the difference between diffusion coefficients is large, the predictions based on simpler theories available in literature can become invalid.