Butler–Volmer with Concentration Terms

Introduction

The simple Butler–Volmer equation treats concentrations at the reaction plane as fixed. This is a good starting point for small currents or vigorous transport, but a reacting electrode consumes one species and produces another. Their concentrations at the surface may then differ markedly from values sampled in the bulk solution. A current–potential curve can bend or flatten because of that difference even when the underlying elementary electron-transfer kinetics retain the same rate constants. This page adds surface composition to the two partial reaction rates and shows why concentration and potential must be solved together.

Core explanation

For the one-electron couple O + e⁻ ⇌ R, define positive current as oxidation and E as the electrode's reduction potential. Let c O,s and c R,s be concentrations at the reaction plane, while c O,b and c R,b are bulk values. A simple mass-action interfacial law is j = F[k ox(E)c R,s − k red(E)c O,s]. The first term is positive anodic oxidation; the second is negative cathodic reduction. The rate constants are heterogeneous, with dimensions of length per time if concentrations are in amount per volume and j is current per area. This form is a deliberately idealized one-electron model; activities, adsorption or coupled chemical reactions can require a richer law.

Take E°′ as a formal potential for the specified medium. A common parameterization is k red(E) = k⁰ exp[−αF(E−E°′)/RT] and k ox(E) = k⁰ exp[(1−α)F(E−E°′)/RT], with 0 < α < 1. The formal potential incorporates a chosen concentration convention and solution environment. IUPAC notes that formal electrode potential depends on electrolyte composition and uses normalized concentrations, so it must not be confused automatically with a thermodynamic standard potential based on activities. A carefully written equation treats logarithm arguments and concentration ratios as dimensionless.

At open-circuit equilibrium with no concentration gradient, surface and bulk values agree. The Nernst relation for this idealized one-electron pair is E eq = E°′ + (RT/F)ln(c O,b/c R,b). Inserting this E eq into the rate law makes the two partial currents equal. Their common magnitude is j₀ = Fk⁰(c O,b)^(1−α)(c R,b)^α when a consistent concentration normalization is understood. This formula makes an important point: the observed exchange current density can change with composition even if the intrinsic k⁰ of a chosen surface is unchanged.

Now define η = E−E eq using that bulk equilibrium state as the reference and normalize each surface concentration by its bulk value. Algebra gives the concentration-aware form

j = j₀[(c R,s/c R,b)exp((1−α)Fη/RT) − (c O,s/c O,b)exp(−αFη/RT)].

If both ratios equal one, this reduces to the equation from the preceding page. If reduction consumes O so that c O,s/c O,b falls, its cathodic partial current is smaller than the fixed-concentration prediction. R may accumulate at the surface, increasing the opposing oxidation term. The bulk-referenced η here is convenient for a transport calculation. One can instead use a local surface-equilibrium potential, but then the local Nernst shift and prefactors must be transformed consistently. Mixing a bulk η with a surface-based j₀ without the corresponding concentration factors double-counts or omits effects.

The surface concentrations are not arbitrary fitting decorations. Diffusion, migration and convection determine how reactants arrive and products depart. At steady state, the faradaic consumption rate at the interface must match the supply flux through the adjoining electrolyte. Under a simple planar diffusion-layer approximation, a reactant O supplied from the bulk has flux approximately D O(c O,b−c O,s)/δ, where δ is an effective layer thickness. For a one-electron reduction, the magnitude of its transport contribution is F times that flux. As c O,s approaches zero, further potential drive cannot make negative concentration; the current tends toward a transport-controlled limit. IUPAC's description of limiting diffusion current explicitly links the limiting condition to near-zero surface concentration in its polarographic setting. The exact limiting-current formula depends on geometry and flow.

This coupling is why the apparent Tafel slope may curve at high current. A plot of measured log j against η contains both exponential electron-transfer response and changes in surface composition. Solution resistance can further shift the actual interfacial potential from the potentiostat's displayed value. Rotating-disk speed, stirring, concentration, temperature and electrode size are useful independent controls. A primary educational treatment of voltammetric kinetics expresses forward and reverse rates using surface rather than bulk concentrations for precisely this reason.

Step-by-step reasoning

1. Write O + e⁻ ⇌ R and declare the positive current direction. 2. Distinguish c O,s and c R,s from the measured or prepared bulk concentrations. 3. Choose the formal or standard potential appropriate to the concentration or activity convention. 4. Find E eq from the bulk ratio and define η = E−E eq for that chosen state. 5. Evaluate the two potential-dependent partial currents multiplied by their surface concentrations. 6. Couple the resulting j to a transport flux or independently measured surface composition; check that no predicted surface concentration becomes negative.

Visual explanation

Plot concentration versus distance away from a cathodic electrode. Show c O,s at the surface lower than c O,b far away; show product R rising from c R,b to c R,s near the surface. Beside it, draw a current–potential curve: the kinetic-only cathodic branch grows exponentially, while the transport-coupled branch bends toward a horizontal limiting-current region. Label the diffusion layer and note that its thickness changes with convection. The surface value, not the bulk label on a bottle, enters the interfacial reaction law.

Real-world analogy

A restaurant may list abundant ingredients in its warehouse, but a cook can use only what reaches the counter. Faster cooking lowers stock at the counter unless delivery speeds up. Bulk concentration resembles warehouse stock, surface concentration resembles the counter, and interfacial kinetics resembles the cook. Unlike a restaurant, molecules move by diffusion, migration and flow; their concentration gradients and reaction rates can be expressed quantitatively.

Real-world example

An analyst reduces a dissolved metal ion at an electrode. At low drive, a small change in potential causes a large kinetic change. At stronger cathodic drive, the ion concentration immediately at the electrode falls while the beaker's bulk concentration is still almost unchanged. Increasing rotation speed raises the measured plateau current. That dependence supports a transport interpretation; it does not mean the electron-transfer barrier became lower when the electrode rotated. To estimate intrinsic k⁰, the analyst must fit a model that accounts for local concentrations and resistive losses.

Why?

Why can current stop growing even though the applied potential keeps moving in the reduction direction? The potential may continue making charge transfer faster for any O molecule that reaches the surface, but supply of O can become the bottleneck. Once c O,s is near zero, the interface cannot consume more O than diffusion or convection delivers. Further drive may promote a different reaction instead, such as solvent reduction, so a plateau belongs to a specified product and potential window rather than to all possible current.

Common misconception

“The beaker concentration is the concentration in the Butler–Volmer rate law.” The rate law applies at the reaction plane, and bulk and surface values separate at finite current. Another error is interpreting a flattened current curve as proof that electron transfer has stopped responding to potential; transport can hide the kinetic increase. A third is treating a limiting current as universal for a species. It changes with bulk concentration, electrode geometry, diffusion coefficient and flow conditions.

Worked example

For a one-electron couple at 298 K, take α = 0.50, j₀ = 1.00 mA cm⁻² and η = −0.059 V relative to the bulk equilibrium potential. Suppose independent transport analysis gives c O,s/c O,b = 0.40 and c R,s/c R,b = 1.20. Since F η /(RT) ≈ 2.30, the anodic exponential is exp(−1.15) ≈ 0.317 and the cathodic exponential is exp(+1.15) ≈ 3.16. Then j ≈ 1.00[(1.20)(0.317) − (0.40)(3.16)] = −0.884 mA cm⁻². Incorrectly using both concentration ratios as one would give about −2.84 mA cm⁻². The difference is not a change in α; it follows from reactant depletion and product accumulation at the surface. A self-consistent model must also check that the chosen surface ratios satisfy the corresponding transport fluxes.

Quick check

1. Which concentration enters the reduction term of the interfacial rate law for O + e⁻ ⇌ R? Answer: The concentration of O immediately at the electrode surface, c O,s, enters the reduction term. 2. What physical change often causes a cathodic limiting current for dissolved O? Answer: Its surface concentration approaches zero because transport cannot deliver O as rapidly as the electrode could consume it.

Exam focus

Write the two partial terms with correct species and signs before simplifying. Keep bulk and surface subscripts visible throughout a derivation. State whether E°′ and E eq use concentration or activity conventions. In a numerical calculation, check dimensionless concentration ratios and the sign of η. Explain how increasing stirring changes a transport-limited current without necessarily changing the intrinsic charge-transfer constant. A high-current Tafel fit that ignores depletion should be treated skeptically.

Advanced insight

Real electrochemical interfaces can have an electrical migration term, activity gradients, specific adsorption and a compact-layer potential drop. In concentrated electrolytes, species activities and transport coefficients can vary with composition. Porous electrodes add spatial gradients inside their pores, so the “surface” values differ from particle to particle. The concentration-aware Butler–Volmer law becomes a local boundary condition coupled to continuum mass-transport equations, rather than a stand-alone formula fitted to terminal current. This is the bridge from a simple voltammetry model to full battery and fuel-cell simulations.

Summary

Interfacial electron-transfer rates depend on concentrations at the reaction plane. A concentration-aware Butler–Volmer equation multiplies the oxidation and reduction exponentials by the appropriate surface-to-bulk ratios when overpotential is referenced to the bulk equilibrium state. Depletion or buildup changes the measured current and can create a limiting-current region. Transport and kinetics must be combined consistently before assigning a curve shape to an intrinsic electron-transfer parameter.

Practice questions

1. During strong reduction of O, is c O,s usually greater or less than c O,b? Answer: It is usually less, because O is consumed at the electrode faster than the bulk can fully replenish it. 2. If c R,s/c R,b increases while all other factors are fixed, which partial term grows in the stated sign convention? Answer: The positive anodic oxidation term grows because more R is available at the interface. 3. Why can a measured exchange current density change when bulk composition changes even if k⁰ stays fixed? Answer: j₀ contains concentration factors for O and R, so changing the bulk ratio or levels changes the equilibrium partial-current magnitudes. 4. What experiment could help distinguish transport limitation from a purely kinetic current plateau? Answer: Change controlled convection, such as rotating-disk speed; a transport limit should change with the supply rate.