simulateCV {eChem}

Simulate a Cyclic Voltammetry Experiment

Description

Simulates a cyclic voltammetry experiment as either an E, EC, or CE mechanism, where E is a redox reaction and where C is a chemical reaction that either precedes or follows the redox reaction.

Usage

simulateCV(e.start = 0, e.switch = -0.5, e.form = -0.25,
  mechanism = c("E", "EC", "CE"), ko = 1, kcf = 0, kcr = 0, n = 1,
  alpha = 0.5, d = 1e-05, area = 0.01, temp = 298.15, scan.rate = 1,
  conc.bulk = 0.001, t.units = 2000, x.units = 180, sd.noise = 0)

Arguments

e.start

Initial potential (in volts).

e.switch

Switching potential (in volts).

e.form

Formal potential for the redox reaction (in volts).

mechanism

Mechanism for the electrochemical system; one of E for redox reaction only, EC for redox reaction with a following chemical reaction, or CE for redox reaction with a preceding chemical reaction. Default is E.

ko

Standard heterogeneous electron transfer rate constant for the redox reaction (in cm/s).

kcf

Homogeneous first-order rate constant for the forward chemical reaction (in s^-1).

kcr

Homogeneous first-order rate constant for the reverse chemical reaction (in s^-1).

n

Number of electrons in the redox reaction.

alpha

Transfer coefficient.

d

Diffusion coefficient for Ox and Red (in cm^2 s^-1).

area

Surface area of the electrode (in cm^2).

temp

Temperature (in K).

scan.rate

Rate at which the potential is changed (in V/s).

conc.bulk

Initial bulk concentration of Ox or Red for an E or an EC mechanism, or the combined initial concentrations of Ox and Z, or of Red and Z for a CE mechanism (in mol/L).

t.units

The number of increments in time for the diffusion grids.

x.units

The number of increments in distance for the diffusion grids.

sd.noise

The standard deviation for noise as a percent of maximum current (in \muA).

Value

Returns a list with the following components

expt

type of experiment; defaults to CV for a cyclic voltammetry simulation

mechanism

type of mechanism used for the simulation

file_type

value that indicates whether the output includes all data (full) or a subset of data (reduced); defaults to full for cvSim

current

vector giving the current as a function of time

potential

vector giving the potential as a function of time

time

vector giving the times used for the diffusion grids

distance

vector giving the distances from electrode surface used for the diffusion grids

oxdata

diffusion grid, as a matrix, giving the concentration of Ox

reddata

diffusion grid, as a matrix, giving the concentrations of Red

chemdata

diffusion grid, as a matrix, giving the concentrations of Z

formalE

formal potential for the redox reaction

initialE

initial potential

switchE

switching potential

electrons

number of electrons, n, in the redox reaction

ko

standard heterogeneous electron transfer rate constant

kcf

homogeneous first-order rate constant for forward chemical reaction

kcr

homogeneous first-order rate constant for reverse chemical reaction

alpha

transfer coefficient

diffcoef

diffusion coefficient for Ox and Red

area

surface area for electrode

temperature

temperature

scanrate

scan rate

conc.bulk

initial concentration of Ox or Red for an E or EC mechanism, or the combined initial concentrations of Ox and Z, or of Red and Z for a CE mechanism

tunits

the number of increments in time for the diffusion grids

xunits

the number of increments in distance for the diffusion grids

sdnoise

standard deviation, as percent of maximum current, used to add noise to simulated data

direction

-1 for an initial reduction reaction of Ox to Red; +1 for an initial oxidation reaction of Red to Ox

k_f

vector of forward electron transfer rate constant as a function of potential

k_b

vector of reverse electron transfer rate constant as a function of potential

jox

vector giving the flux of Ox to the electrode surface as a function of potential

jred

vector giving the flux of Red to the electrode surface as a function of potential

Examples

ex_cv = simulateCV(e.start = 0.25, e.switch = -0.25, e.form = 0,
  x.units = 100, t.units = 1000)
str(ex_cv)