Back to Blog•Core Philosophy•September 15, 2026

The Curve We Cannot Read:
On the Structural Indistinguishability of Interfacial Kinetics

Why Butler–Volmer, its resistance-corrected cousins, and Marcus–Hush–Chidsey/CIET kinetics all agree near equilibrium, diverge only where the data cannot yet be trusted, and what breaking that degeneracy would require.

"At the end of my last essay I named the shape of the problem without doing the work: before we can decide whether Butler–Volmer or CIET is correct at high overpotential, we might first need our own Bates–Guggenheim moment, a convention that makes the measurement trustworthy rather than a better guess at the mechanism. I never said what makes the measurement untrustworthy in the first place. That is what this essay is for."

Plot current against overpotential for almost any electrode reaction and you get a shape that looks reassuringly familiar: flat near the origin, rising into a curve, eventually bending over into something like a plateau. Every kinetics textbook shows a version of this curve. The reassuring part is the problem.

Let me put the claim precisely. Over the overpotential range a real cell gives us, several genuinely different hypotheses compete to describe the interface: symmetric Butler–Volmer, asymmetric Butler–Volmer, resistance-corrected ("modified") Butler–Volmer, Marcus–Hush–Chidsey, and CIET. Their agreement runs deeper than a shared ability to fit the same data. To leading order, they are the same equation. The differences meant to distinguish them surface only in a regime where two unrelated artifacts produce the identical signature. Better instruments would not settle it; the obstruction is in the algebra.


The Family We're Actually Choosing Among

It is easy to talk about "BV versus CIET" as if there were two candidates. What we are choosing among is a family, all of it built from the same move: assume a shape for the free-energy landscape of the elementary charge-transfer step, then read off a current from it.

  • Symmetric Butler–Volmer: a fixed 50/50 split of the driving force between the anodic and cathodic activation barriers.
  • Asymmetric Butler–Volmer: the same split, but at a fitted ratio α:(1−α)\alpha : (1-\alpha) instead of 1/2:1/21/2 : 1/2.
  • Modified (resistance-corrected) Butler–Volmer: either of the above, with an assumed series resistance subtracted from the applied bias before fitting.
  • Marcus–Hush–Chidsey (MHC) / Gerischer: a parabolic free-energy landscape integrated over a continuum of electronic states, appropriate for a metal or, with a weighting density of states, a semiconducting oxide like NMC or LCO.
  • CIET: the MHC formalism extended so that the transferring electron and the inserting ion move together, still inheriting MHC's underlying activation picture.

Every member of this family ends in the same place: a current–overpotential relation i(η)i(\eta). The question I want to answer is: how much of that function do we actually need to see before the family members stop agreeing with each other?


Where the Whole Family Agrees

Start with the general (asymmetric) Butler–Volmer equation, since it is the one member of the family we can expand by hand without special functions:

i=i0[exp⁡ ⁣(αFηRT)−exp⁡ ⁣(−(1−α)FηRT)]i = i_0\left[\exp\!\left(\dfrac{\alpha F \eta}{RT}\right) - \exp\!\left(-\dfrac{(1-\alpha) F \eta}{RT}\right)\right]

Let x≡Fη/RTx \equiv F\eta / RT, so that i/i0=eαx−e−(1−α)xi/i_0 = e^{\alpha x} - e^{-(1-\alpha)x}. Expand each exponential to third order:

eαx≈1+αx+(αx)22+(αx)36,e−(1−α)x≈1−(1−α)x+[(1−α)x]22−[(1−α)x]36e^{\alpha x} \approx 1 + \alpha x + \dfrac{(\alpha x)^2}{2} + \dfrac{(\alpha x)^3}{6}, \qquad e^{-(1-\alpha)x} \approx 1 - (1-\alpha)x + \dfrac{[(1-\alpha)x]^2}{2} - \dfrac{[(1-\alpha)x]^3}{6}

Subtracting term by term, the constant "1"s cancel, the linear terms add (αx+(1−α)x=x\alpha x + (1-\alpha)x = x), and the quadratic terms combine as a difference of squares, (αx)2−[(1−α)x]2=(2α−1)x2(\alpha x)^2 - [(1-\alpha)x]^2 = (2\alpha - 1)x^2. What survives is:

ii0=x+2α−12x2+O(x3)\dfrac{i}{i_0} = x + \dfrac{2\alpha - 1}{2}x^2 + O(x^3)

Look at where α\alpha first enters: not in the leading term, only in the second-order one. To leading order,

i≈i0x=Fi0RT η≡ηRct,Rct≡RTFi0i \approx i_0 x = \dfrac{F i_0}{RT}\,\eta \equiv \dfrac{\eta}{R_{ct}}, \qquad R_{ct} \equiv \dfrac{RT}{F i_0}

Key InsightThe charge-transfer resistance RctR_{ct} that a small-signal EIS measurement returns carries no information about α\alpha at all. The expansion never gave it a place to appear: the transfer coefficient — and with it the entire question of how asymmetric the barrier is — enters only at second order, buried under whatever double-layer charging, thermodynamic drift, and ordinary noise already crowd that regime.

The same conclusion holds, model for model, if the exponentials are replaced with the Marcus–Hush–Chidsey rate integral that CIET is built on. Expanded near equilibrium, the MHC current reduces to the identical linear form, i≈i0Fη/RTi \approx i_0 F\eta/RT, with an effective transfer coefficient of exactly 1/21/2 (Chidsey, 1991; Zeng et al., 2014). A fixed asymmetric split and a full Marcus parabola integrated over an electrode's density of states are, in this regime, the same number wearing two different derivations. Neither leaves a trace in the one quantity a linear i-ηi\text{-}\eta or small-amplitude EIS measurement can deliver.


Where the Family Is Supposed to Split

The disagreement is real, but it only shows up far from equilibrium. Butler–Volmer has no ceiling: each exponential term grows without bound as η\eta grows, so the current keeps climbing on a Tafel plot forever. Marcus theory behaves differently by construction. The activation energy for a single barrier is a parabola in the driving force,

ΔG‡(η)=(λ−Fη)24λ\Delta G^{\ddagger}(\eta) = \dfrac{(\lambda - F\eta)^2}{4\lambda}

where λ\lambda is the reorganization energy. As η\eta grows from zero, this barrier first shrinks, reproducing something that looks Tafel-like, until Fη=λF\eta = \lambda, where the reaction becomes barrierless. Push η\eta past that point for an isolated molecular reaction and the barrier grows again: the famous Marcus inverted region, where the rate turns over and falls.

An electrode does not behave this way, which is Chidsey's (1991) central point. A metal or a semiconducting oxide offers a whole continuum of electronic states, not one. Past the barrierless point, there is always some subset of states with the right energy to react without paying the inverted-region penalty, so the net current saturates instead of turning over. CIET inherits this saturating, S-shaped curve because it is built on the same MHC formalism. The one place the family is supposed to disagree with Butler–Volmer is exactly this large-η\eta plateau, a feature set by λ\lambda, and invisible to any EIS spectrum taken at low bias.

Unbounded family

Butler–Volmer, symmetric or asymmetric

A fixed exponential in η\eta, with no built-in mechanism to stop climbing. Any bend in a real curve has to be explained by something added on top, usually transport or resistance.

Saturating family

Marcus–Hush–Chidsey / Gerischer / CIET

Bends over on its own, intrinsically, once FηF\eta approaches λ\lambda. The bend is a real prediction of the mechanism, not a correction applied afterward.


Two Impostors Wearing the Same Bend

A saturating curve is exactly what an intrinsic Marcus-type mechanism predicts. It is unfortunately also what two purely extrinsic effects predict, at similar current densities, in almost any real cell.

The first is mass transport. As η\eta grows, the reactant at the surface depletes, and the current approaches a diffusion-limited value that, for a film of thickness LL, scales roughly as

ilim∼nFDc∗Li_{\text{lim}} \sim \dfrac{nFDc^{*}}{L}

This plateau has nothing to do with the reorganization energy. But taken alone, an i-ηi\text{-}\eta curve that flattens at high η\eta cannot tell you whether it flattened because of λ\lambda or because of LL.

The second is uncompensated resistance. Any series resistance between the reference point and the interface shifts the true local overpotential away from the applied one, ηtrue=ηapplied−iRs\eta_{\text{true}} = \eta_{\text{applied}} - iR_s. At high current, this steals back exactly the kind of curvature that would otherwise look like unbounded Tafel growth. This is precisely why "modified Butler–Volmer" schemes exist: fit RsR_s alongside i0i_0 and α\alpha, and subtract its effect before reporting kinetics.

I used to think of resistance correction as a fix, a way of cleaning the data before the real kinetic question gets asked. I don't think that anymore. Adding RsR_s as a free parameter does not remove the degeneracy between an intrinsic saturation and an artifact; it relocates it. A curve that bends over at high η\eta can now be explained by any combination of a real λ\lambda, a real diffusion length LL, and an assumed RsR_s that trades off against both. "Modified BV" sits inside this family alongside CIET — a third curve that bends like the others, for a reason of its own.


Counting What We Can Actually Solve For

Put the pieces together and the problem is one of counting. A real, finite-range i-ηi\text{-}\eta curve, measured with real noise, realistically yields perhaps three independent features: an initial slope (RctR_{ct}), the overpotential at which curvature becomes visible, and a limiting or asymptotic current. Any model in the family with three or more free parameters can be tuned to reproduce those three features from more than one underlying mechanism, and every model past plain symmetric BV clears that bar the moment you add α\alpha, or λ\lambda, or RsR_s, or LL. The map from "which physical process is happening" to "what the curve looks like" is not one-to-one. Fitting harder does not invert a map that isn't injective.

I don't think this is a coincidence, and I don't think it is unique to kinetics. It is the same shape of problem I keep running into when people claim a DRT deconvolution has cleanly separated charge-transfer resistance from diffusion admittance in an impedance spectrum: an underdetermined inverse problem, dressed up differently depending on which domain (time, frequency, or overpotential) you happen to be squinting at. That comparison deserves its own essay. For now, the point is narrower: no amount of clever regression rescues an argument of the form "the fit is good, therefore the mechanism is right," when several mechanisms are guaranteed to produce equally good fits by construction.


What Would Actually Break the Degeneracy

If the curve itself can't decide, the only way out is architectural: build an experiment where the competing explanations are forced to make different predictions about something other than the shape of one curve.

The diffusion-versus-intrinsic-saturation confound has a clean handle, because the two depend on geometry differently. An intrinsic, λ\lambda-controlled plateau does not care how thick the electrode is. A diffusion-limited plateau scales predictably with the transport length. Vary that length scale independently, rather than inferring it after the fact from the same curve you're trying to interpret, and the two hypotheses stop agreeing.

The resistance-versus-intrinsic-saturation confound needs the same kind of move: a way of reading the true interfacial η\eta directly, close enough to the reaction plane that RsR_s stops being a free parameter to fit around and becomes something the architecture pins down on its own. And underneath both sits the confound my last two essays were circling: none of this length-scale or resistance bookkeeping survives a reference potential that drifts as the reaction proceeds, which is exactly what an intercalation electrode's Ueq(x)U_{eq}(x) does and a catalytic electrode's does not. That has to be settled first; the other two become answerable only once the baseline holds still.

"A kinetic model doesn't earn credibility by fitting the accessible curve. It earns credibility by surviving an architecture built to make the other candidates fail differently. The real question was never which curve fits best; it's which architecture would even let the curves disagree."

Working out what such an architecture looks like for a lithium-intercalation electrode, and being honest this time about the assumptions it would quietly carry, is where I want to go next.

References & Key Literature

  1. Marcus, R. A. "On the Theory of Oxidation–Reduction Reactions Involving Electron Transfer. I." J. Chem. Phys. 24, 966–978 (1956).
  2. Chidsey, C. E. D. "Free Energy and Temperature Dependence of Electron Transfer at the Metal–Electrolyte Interface." Science 251, 919–922 (1991).
  3. Zeng, Y., Smith, R. B., Bai, P. & Bazant, M. Z. "Simple Formula for Marcus–Hush–Chidsey Kinetics." J. Electroanal. Chem. 735, 77–83 (2014).
  4. Zhang, Y. et al. "Lithium-ion intercalation by coupled ion–electron transfer." Science 390, eadq2541 (2025).
  5. Bai, P. & Bazant, M. Z. "Charge transfer kinetics at the solid–solid interface in porous electrodes." Nat. Commun. 5, 3585 (2014).
  6. Bard, A. J. & Faulkner, L. R. Electrochemical Methods: Fundamentals and Applications, 2nd ed. Wiley (2001).