Description:
The attached file is an updated version of a simple, easy-to-use cyclic voltammetry (CV) simulator written in Microsoft Excel (for details, see Brown, J. H.; J. Chem. Educ., 2016, 93 (7), 1326–1329. DOI: 10.1021/acs.jchemed.6b00052). The spreadsheet supports up to four species undergoing up to four redox couples each and assumes the first process is a charge transfer.
Experimental CVs with up to 10,000 data points in text format can be easily copied/pasted into the spreadsheet for comparison to the simulation. The spreadsheet can be operated in protected mode or unprotected and modified as needed. The spreadsheet is provided with no guarantee and it may not be suitable for all applications.
System requirements: Microsoft Windows 8.1 Office software or later. New version (v_12) uploaded to website on 12/10/18.
Important notes:
You must list your potential values in the proper order for your simulations to work. List your formal redox potential (E^{o}') values in decreasing order if reductions occur during the forward sweep. Reverse this order if oxidations occur during the forward sweep. You can then coalesce two adjacent peaks by exchanging their E^{o}' values. Set the E^{o}' values of unused peaks beyond the switching potential (E_{2}) in the proper order noted above to exclude them from the model. Set the starting potential (E_{1}) value greater than the switching potential (E_{2}) if reductions occur during the forward sweep. Reverse this order if oxidations occur during the forward sweep.
Updates:
Added first-order intervening rate constants (k_{I1}) between electron transfers on 12/20/16. The smaller you make a k_{I1} value, the slower the corresponding species is populated. Symbols k_{R1} and k_{R2} are now used to indicate first- and second-order apparent forward homogeneous rate constants that remove a corresponding species from the cycle.
Separated the equations involving the number of electrons transferred on 12/22/16. The number of electrons in the rate determining step (n_{a}) now controls the peak width. The total number of electrons transferred (n) now controls the total current. For details, see Bard, A. J.; Faulkner, L. R. in Electrochemical Methods: Fundamentals and Applications, Wiley: New York, 1980.
Incorporated the Saul’yev RL variant to simulate Fick’s Second Law of Diffusion on 1/8/17. This allows the simulation of thin reaction layers. For details, see Britz, D. in Digital Simulation in Electrochemistry 3^{ed} ed., Springer: Berlin, 2010. Reducing the variable x_{scale }reduces the length of the diffusion grids and places the concentration points closer to the working electrode surface. However, using diffusion grids that are too short or long will distort your simulations, so adjust this parameter only as needed.
Moved the intermediate calculations for the Saul’yev RL variant to their own rows to avoid confusion with the redox variables on 12/9/18. Also removed the page for the simulation of adsorbed species for simplicity.
Variables:
Several variables from the original spreadsheet were renamed to incorporate new features. The following table describes the variables that are used in the newest version of the spreadsheet.
Variable | Definition | Value and/or unit(s) |
Gain | Amplification factor for experimental background | Unitless multiplier |
Offset | Constant added to the experimental current | A |
E^{o}ˊ | Formal potential for a redox couple | V |
n | Total number of electrons transferred per species (controls total peak area) | Unitless |
n_{a} | Number of electrons transferred in the rate determining step (controls peak width) | Unitless |
D | Common diffusion coefficient | ≈ 1 x 10^{-5} cm^{2}/s |
k^{o} | Standard heterogeneous rate constant for an electron transfer | cm/s |
a | Transfer coefficient for a redox couple | ≈ 0.5 (unitless) |
C_{bulk} | Bulk concentration of the analyte | mol/cm^{3} |
E_{1} | Starting potential of the CV experiment | V |
E_{2} | Switching potential of the CV experiment | V |
v | Scan rate used during the CV experiment | V/s |
A | Surface area of the working electrode | cm^{2} |
T | Temperature of the sample during experiment | K |
x_{scale} | Compress or expand x-axes of diffusion grids | Percentage (%) |
k_{I1} | Rate constant for first-order intervening homogenous reaction between electron transfers | s^{-1} |
k_{R1} | Rate constant for a first-order homogeneous reaction that removes a species from electron transfer cycle | s^{-1} |
k_{R2} | Rate constant for second-order homogeneous reaction that removes a species from electron transfer cycle | cm^{3}/mol•s |
Bug report:
Bug in the time increment calculations was corrected on 11/16/15. Bug in the concentration calculations of oxidized form of the 4^{th} redox couple was corrected on 12/8/15. Bug in the Import data page graphics was corrected on 12/31/15. Bug in the data selection process was corrected on 4/25/16. Bug in the intervening homogeneous reactions was corrected on 1/12/17. Corrected a minor bug in the stacked plot on 12/11/18.