Butler–Volmer equation

In electrochemistry, the Butler–Volmer equation (named after John Alfred Valentine Butler[1] and Max Volmer), also known as Erdey-Grúz–Volmer equation, is one of the most fundamental relationships in electrochemical kinetics.

It describes how the electrical current through an electrode depends on the voltage difference between the electrode and the bulk electrolyte for a simple, unimolecular redox reaction, considering that both a cathodic and an anodic reaction occur on the same electrode:[2] The Butler–Volmer equation is: or in a more compact form: where: The right hand figure shows plots valid for

The theoretical values of the Tafel equation constants are different for the cathodic and anodic processes.

The more general form of the Butler–Volmer equation, applicable to the mass transfer-influenced conditions, can be written as:[3] where: The above form simplifies to the conventional one (shown at the top of the article) when the concentration of the electroactive species at the surface is equal to that in the bulk.

There are two rates which determine the current-voltage relationship for an electrode.

First is the rate of the chemical reaction at the electrode, which consumes reactants and produces products.

The second is the rate at which reactants are provided, and products removed, from the electrode region by various processes including diffusion, migration, and convection.

The simple Butler–Volmer equation assumes that the concentrations at the electrode are practically equal to the concentrations in the bulk electrolyte, allowing the current to be expressed as a function of potential only.

Despite this limitation, the utility of the Butler–Volmer equation in electrochemistry is wide, and it is often considered to be "central in the phenomenological electrode kinetics".

[4] The extended Butler–Volmer equation does not make this assumption, but rather takes the concentrations at the electrode as given, yielding a relationship in which the current is expressed as a function not only of potential, but of the given concentrations as well.

The mass-transfer rate may be relatively small, but its only effect on the chemical reaction is through the altered (given) concentrations.

A full treatment, which yields the current as a function of potential only, will be expressed by the extended Butler–Volmer equation, but will require explicit inclusion of mass transfer effects in order to express the concentrations as functions of the potential.

The following derivation of the extended Butler–Volmer equation is adapted from that of Bard and Faulkner[3] and Newman and Thomas-Alyea.

[5] For a simple unimolecular, one-step reaction of the form: The forward and backward reaction rates (vf and vb) and, from Faraday's laws of electrolysis, the associated electrical current densities (j), may be written as: where kf and kb are the reaction rate constants, with units of frequency (1/time) and co and cr are the surface concentrations (mol/area) of the oxidized and reduced molecules, respectively (written as co(0,t) and cr(0,t) in the previous section).

The net rate of reaction v and net current density j are then:[Note 2] The figure above plots various Gibbs energy curves as a function of the reaction coordinate ξ.

The reaction coordinate is roughly a measure of distance, with the body of the electrode being on the left, the bulk solution being on the right.

The blue energy curve shows the increase in Gibbs energy for an oxidized molecule as it moves closer to the surface of the electrode when no potential is applied.

Applying a potential E to the electrode will move the energy curve downward[Note 3] (to the red curve) by nFE and the intersection point will move to

[Note 4] Assume that the rate constants are well approximated by an Arrhenius equation, where the Af and Ab are constants such that Af co = Ab cr is the "correctly oriented" O-R collision frequency, and the exponential term (Boltzmann factor) is the fraction of those collisions with sufficient energy to overcome the barrier and react.

Assuming that the energy curves are practically linear in the transition region, they may be represented there by: The charge transfer coefficient for this simple case is equivalent to the symmetry factor, and can be expressed in terms of the slopes of the energy curves: It follows that: For conciseness, define: The rate constants can now be expressed as: where the rate constants at zero potential are: The current density j as a function of applied potential E may now be written:[5]: § 8.3 At a certain voltage Ee, equilibrium will attain and the forward and backward rates (vf and vb) will be equal.

Note that the net current density at equilibrium will be zero.

The equilibrium rate constants are then: Solving the above for kfo and kbo in terms of the equilibrium concentrations coe and cre and the exchange current density jo, the current density j as a function of applied potential E may now be written:[5]: § 8.3 Assuming that equilibrium holds in the bulk solution, with concentrations

Note that E-Ee is also known as η, the activation overpotential.

For the simple reaction, the change in Gibbs energy is:[Note 5] where aoe and are are the activities at equilibrium.

is the standard potential Defining the formal potential:[3]: § 2.1.6 the equilibrium potential is then: Substituting this equilibrium potential into the Butler–Volmer equation yields: which may also be written in terms of the standard rate constant ko as:[3]: § 3.3.2 The standard rate constant is an important descriptor of electrode behavior, independent of concentrations.

It is a measure of the rate at which the system will approach equilibrium.

For nearly ideal electrodes with small ko, large changes in the overpotential are required to generate a significant current.

The upper graph shows the current density as function of the overpotential η . The anodic and cathodic current densities are shown as j a and j c , respectively for α=α a c =0.5 and j 0 =1mAcm −2 (close to values for platinum and palladium). The lower graph shows the logarithmic plot for different values of α (Tafel plot).
A plot of various Gibbs energies as a function of reaction coordinate. The reaction will proceed towards the lower energy - reducing for the blue curve, oxidizing for the red curve. The green curve illustrates equilibrium.