Highly accurate, efficient, and automatic computation of reversible cyclic voltammograms, using double exponential formulas for numerical integration

Highly accurate reference values of theoretical cyclic voltammograms are indispensable for those who develop and examine simulation methods or other modelling approaches for electroanalytical chemistry. Reference values having relative errors close to +/- 10(-16) (which is the level of the representation errors of standard double precision variables) are desirable, but they cannot be obtained efficiently using conventional simulations. A class of reversible cyclic voltammograms can be expressed in terms of improper convolution integrals that have to be evaluated numerically. By employing the efficient automatic numerical integration package INTDE of Ooura, based on double exponential formulas, and written in C, it is shown that such voltammograms can be calculated with relative errors close to +/- 10(-16), in a time of ca. 1 ms per single current function value, on a contemporary laptop computer. Furthermore, this method appears resistant to failures that may accompany highly accurate numerical integration by alternative methods and software. This is demonstrated for three kinetic examples under conditions of pure diffusion transport: reversible charge transfer at a planar electrode; reversible charge transfer at a spherical electrode; and ErevCirr' mechanism at a planar electrode.