Oxygen Electrode Kinetics: Surface Exchange, Chemical Diffusion, and Characterization
Oxygen electrode performance is governed by a sequence of coupled processes at the gas–solid interface and within the oxide. A material may exhibit high electrical conductivity and appreciable oxygen-ion mobility yet show sluggish electrode kinetics because oxygen adsorption, dissociation, charge transfer, or incorporation is slow. Conversely, a highly active surface cannot compensate indefinitely for slow oxygen transport through the bulk.
This guide develops the physical meaning of the two kinetic parameters most often used to describe oxygen transport in oxides—the chemical surface-exchange coefficient, \(k_{chem}\), and the chemical diffusion coefficient, \(D_{chem}\)—and explains how they are measured and interpreted.
1. Oxygen electrode reactions
In an oxide-ion conducting solid oxide cell, the oxygen electrode supports oxygen reduction in fuel-cell mode and oxygen evolution in electrolysis mode. A simplified oxygen incorporation reaction, written in Kröger–Vink notation, is:
\[ \frac{1}{2}O_2(g) + V_O^{\bullet\bullet} + 2e' \rightleftharpoons O_O^x \]
Here, \(V_O^{\bullet\bullet}\) is an oxygen vacancy, \(e'\) is an electron, and \(O_O^x\) is an oxygen ion occupying a regular oxygen lattice site. The forward direction represents oxygen incorporation; the reverse direction represents oxygen release.
This equation is a useful overall reaction, not a claim that the process occurs in one elementary step. Depending on the oxide and operating conditions, the sequence can involve:
- Gas-phase transport of oxygen to the electrode.
- Molecular oxygen adsorption.
- Dissociation into adsorbed oxygen species.
- Electron transfer and formation of charged adsorbates.
- Incorporation into oxygen vacancies or another available lattice site.
- Bulk oxygen-ion transport through the electrode.
- Charge transfer at the electrode/electrolyte interface.
The order and relative importance of elementary steps are material-dependent. Some electrodes also support significant surface or near-surface transport, and porous composite electrodes introduce additional pathways and interfaces.
For a deeper discussion of defect populations and oxygen non-stoichiometry, see Oxygen Vacancies and Defect Chemistry in Functional Oxides. The broader requirements for simultaneous ionic and electronic transport are discussed in Mixed Ionic–Electronic Conductors.
2. The two kinetic coefficients
2.1 Chemical diffusion coefficient
The chemical diffusion coefficient, \(D_{chem}\), describes the relaxation of a chemical-composition gradient inside a solid. For oxygen-deficient oxides, it is commonly associated with the redistribution of oxygen non-stoichiometry, often represented by \(\delta\) in a composition such as \(ABO_{3-\delta}\).
For a one-dimensional, small perturbation around equilibrium, a simplified diffusion equation is:
\[ \frac{\partial \delta}{\partial t} = D_{chem}\frac{\partial^2\delta}{\partial x^2} \]
The equation assumes an effectively constant \(D_{chem}\) over the perturbation range. If the coefficient depends strongly on composition, temperature, or position, a concentration-dependent transport equation is required.
The units of \(D_{chem}\) are \(m^2\,s^{-1}\) (or commonly \(cm^2\,s^{-1}\)).
Importantly, \(D_{chem}\) is not generally identical to the tracer/self-diffusion coefficient, \(D^*\). Chemical diffusion describes the relaxation of a chemical-potential gradient and therefore incorporates thermodynamic factors as well as mobility. A schematic relationship is:
\[ D_{chem} = D^*\,\Theta \]
where \(\Theta\) represents an appropriate thermodynamic factor. Its exact definition depends on the composition variable, defect model, and convention used. Therefore, numerical values of \(D_{chem}\) and \(D^*\) should not be compared as if they were interchangeable.
2.2 Chemical surface-exchange coefficient
The chemical surface-exchange coefficient, \(k_{chem}\), quantifies how rapidly the solid exchanges oxygen with the surrounding gas in response to a chemical-potential difference.
For a small deviation from equilibrium, a common linearized boundary condition is:
\[ -D_{chem}\left.\frac{\partial \delta}{\partial x}\right|_{s} = k_{chem}(\delta_s-\delta_{eq}) \]
The sign depends on the chosen surface normal and flux convention; the equation above expresses the magnitude of the flux as proportional to the difference between the surface composition and its equilibrium value.
The units of \(k_{chem}\) are \(m\,s^{-1}\) (or commonly \(cm\,s^{-1}\)).
The coefficient is an effective kinetic parameter. It can include contributions from adsorption, dissociation, surface diffusion, charge transfer, and incorporation, depending on the material and experimental conditions. It should not automatically be interpreted as the rate constant of one elementary reaction.
2.3 Why both coefficients matter
A useful first comparison is the characteristic bulk-diffusion time and surface-exchange time for a specimen with characteristic length \(L\):
\[ \tau_{diff}\sim\frac{L^2}{D_{chem}} \]
\[ \tau_{surf}\sim\frac{L}{k_{chem}} \]
Their ratio is described by a dimensionless parameter:
\[ Bi = \frac{k_{chem}L}{D_{chem}} \]
This is analogous to a mass-transfer Biot number.
- If \(Bi\ll 1\), surface exchange is relatively slow compared with bulk diffusion. The specimen tends to remain internally more uniform while its overall oxygen content changes.
- If \(Bi\gg 1\), surface exchange is relatively fast and internal diffusion is more likely to control the relaxation.
- If \(Bi\) is of order unity, both processes can contribute significantly.
The characteristic length must be defined consistently with the specimen geometry and boundary conditions. For a plate exposed on both faces, for example, one may use the half-thickness in a plane-sheet solution. The Biot number is a diagnostic guide, not a substitute for fitting the appropriate diffusion model.
3. What controls oxygen reduction and evolution?
The oxygen electrode reaction is sensitive to both bulk and surface properties.
Bulk-related factors - Oxygen vacancy or interstitial concentration. - Mobility of oxygen defects. - Electronic conductivity needed for charge transfer. - Defect association and ordering. - Crystal structure and migration pathways. - Chemical expansion and strain associated with changes in oxygen content.
Surface-related factors - Surface termination and crystallographic orientation. - Surface oxygen-vacancy concentration. - Cation segregation and secondary surface phases. - Adsorbed water, carbon dioxide, and impurities. - Surface reconstruction and electrode ageing. - Contact with current collectors and other electrode phases.
A high oxygen-vacancy concentration does not guarantee a high exchange rate. Vacancies may be poorly mobile, blocked by surface segregation, or unfavourably arranged at the surface. Similarly, a large specific surface area does not necessarily mean that every surface site is catalytically active.
The dependence on oxygen partial pressure, \(p_{O_2}\), is especially informative. Changes in \(p_{O_2}\) alter oxygen chemical potential and defect populations, and can change both the bulk transport coefficient and the surface reaction kinetics. Comparisons are meaningful only when temperature, gas composition, water content, microstructure, and measurement protocol are controlled.
4. Electrical conductivity relaxation (ECR)
ECR is a widely used method for determining chemical oxygen transport parameters in mixed-conducting oxides.
4.1 Principle
A dense specimen is equilibrated at a fixed temperature and oxygen partial pressure. The surrounding gas composition is then changed rapidly to a new oxygen partial pressure. The specimen responds by changing its oxygen content and, consequently, its electrical conductivity.
The conductivity response is recorded as a function of time. A diffusion model with a surface-exchange boundary condition is fitted to the transient to obtain \(D_{chem}\) and \(k_{chem}\).
For a small perturbation, the normalized conductivity response is often written as:
\[ g(t)=\frac{\sigma(t)-\sigma_0}{\sigma_\infty-\sigma_0} \]
where \(\sigma_0\) is the initial conductivity and \(\sigma_\infty\) is the final equilibrium conductivity. This normalization assumes that the conductivity change is measurable and that the response can be related reproducibly to the oxygen-content change.
4.2 Practical ECR workflow
- Prepare a dense specimen with measured dimensions and well-defined exposed surfaces.
- Establish stable temperature and gas flow at the initial \(p_{O_2}\).
- Record the initial resistance or conductivity after equilibrium is reached.
- Apply a controlled gas switch to the target \(p_{O_2}\).
- Record the conductivity transient with sufficient time resolution.
- Repeat for oxidation and reduction steps, and at additional temperatures where appropriate.
- Fit the transient using the correct specimen geometry and boundary conditions.
- Report fitted parameters, uncertainties, specimen dimensions, gas-switch response, and the fitting model.
A gas-switching system has its own response time. If the gas composition at the specimen does not change as a near-step function, the measured transient contains both the instrument response and the material response. This must be assessed rather than attributed entirely to the specimen.
4.3 Key assumptions and limitations
ECR interpretation commonly assumes:
- The specimen is sufficiently dense and chemically homogeneous.
- The perturbation is small enough for linearized transport to be valid, or a suitable nonlinear model is used.
- The measured conductivity is related to oxygen content in a known and sufficiently monotonic way.
- Temperature and gas composition are stable during each transient.
- The surface boundary condition is adequately represented by the selected model.
- Contact resistance and instrumental drift do not dominate the signal.
A major limitation arises when conductivity is weakly sensitive to oxygen content. In that case, the transient may be noisy or poorly constrained even if oxygen exchange is occurring. A second limitation occurs when the gas-switching time is comparable to the specimen relaxation time.
ECR measures chemical transport parameters, not automatically the tracer coefficients obtained from isotope experiments.
5. Oxygen isotope exchange and depth profiling
Isotope exchange provides a more direct route to studying oxygen motion without relying on a measurable change in electrical conductivity.
A specimen is exposed to oxygen enriched in the stable isotope \(^{18}O\) at a controlled temperature and oxygen partial pressure. Isotope exchange occurs at the surface, and labelled oxygen diffuses into the material. After rapid quenching, the isotope concentration as a function of depth is measured, commonly using secondary-ion mass spectrometry (SIMS).
The resulting profile can be fitted to obtain a tracer diffusion coefficient, \(D^*\), and a tracer surface-exchange coefficient, \(k^*\), under the assumptions of the selected model.
5.1 What the method measures
- \(D^*\) describes the motion of labelled oxygen atoms or ions in the solid.
- \(k^*\) describes isotope exchange at the gas–solid interface.
These are tracer quantities. They should be distinguished from \(D_{chem}\) and \(k_{chem}\) obtained from a chemical perturbation such as ECR. Relationships between tracer and chemical parameters require thermodynamic information and consistent definitions.
5.2 Experimental considerations
The isotope-exchange temperature, \(p_{O_2}\), gas purity, exposure time, and quenching procedure must be controlled. The depth profile should extend far enough to capture the relevant diffusion length, while the near-surface region must be measured with adequate resolution.
Potential sources of error include: - Surface contamination or reaction with the sample holder. - Isotope exchange during heating or cooling outside the intended exposure. - Surface roughness and porosity. - SIMS matrix effects and depth calibration. - Insufficient profile depth or poor near-surface resolution. - Changes in surface chemistry during the isotope-exchange treatment.
For porous electrodes, the interpretation is more complicated because the accessible surface area and diffusion geometry differ from those of a dense polished specimen.
6. Electrochemical impedance spectroscopy (EIS)
EIS is routinely used to characterize oxygen electrode response in complete cells and symmetrical cells. It measures the frequency-dependent response of an electrochemical system to a small perturbation.
A spectrum may contain contributions from: - Gas transport through the porous electrode. - Surface adsorption and oxygen dissociation. - Charge-transfer processes. - Ionic transport through the electrode. - Transport across the electrode/electrolyte interface. - Contact and current-collection effects.
An equivalent circuit can help describe the spectrum, but a fitted resistor–constant-phase-element branch is not automatically a unique elementary reaction. Different physical processes can produce overlapping frequency responses.
For porous mixed-conducting electrodes, transmission-line or porous-electrode models may be more physically appropriate than a simple series of independent \(R\parallel C\) elements. Model selection should be guided by electrode architecture, gas conditions, and independent measurements.
EIS is particularly useful for tracking electrode degradation and comparing electrode architectures. However, extracting intrinsic \(D_{chem}\) and \(k_{chem}\) from porous-electrode spectra generally requires a physically justified model and supporting information such as porosity, tortuosity, active surface area, and oxygen non-stoichiometry.
See the Practical Guide to Electrochemical Impedance Spectroscopy for measurement quality, equivalent circuits, and interpretation.
7. Surface exchange versus bulk diffusion: identifying the limiting step
A frequent experimental mistake is to assume that the electrode with the highest total conductivity must have the best oxygen-reduction kinetics. Total conductivity alone does not establish the rate of oxygen incorporation or release.
Several comparisons can help identify the controlling process:
| Observation | Possible interpretation | Important check |
|---|---|---|
| Strong dependence on specimen thickness in a relaxation experiment | Bulk diffusion may be significant | Use the correct geometry-dependent solution |
| Similar relaxation time across different thicknesses | Surface exchange may dominate | Confirm gas-switching time is not limiting |
| Large improvement after surface treatment with little change in bulk conductivity | Surface kinetics may be limiting | Check whether treatment changed surface composition or roughness |
| EIS response changes with electrode thickness or porosity | Distributed transport and active area may matter | Characterize microstructure |
| Performance decreases after ageing despite stable bulk phase | Surface segregation or contamination may be involved | Analyse surface chemistry directly |
These are diagnostic clues, not universal proofs. For example, changing thickness can also change microstructure, and a surface treatment may alter both surface chemistry and near-surface defect populations.
A robust study combines at least two complementary approaches where possible—for example, ECR with isotope exchange, or EIS with microstructural and surface-chemical characterization.
8. Temperature dependence and activation energy
Both bulk diffusion and surface exchange commonly show thermally activated behaviour over a suitable temperature interval. An Arrhenius representation is:
\[ D_{chem}=D_0\exp\left(-\frac{E_{D}}{k_BT}\right) \]
\[ k_{chem}=k_0\exp\left(-\frac{E_{k}}{k_BT}\right) \]
where \(D_0\) and \(k_0\) are pre-exponential factors, \(E_D\) and \(E_k\) are apparent activation energies, \(k_B\) is the Boltzmann constant, and \(T\) is absolute temperature.
For a plot of \(\ln D_{chem}\) versus \(1/T\), the slope is \(-E_D/k_B\). The equivalent relation applies to \(k_{chem}\).
The extracted activation energy is an apparent value for the measured temperature range. It may reflect several coupled processes and can change when the rate-limiting mechanism changes. A straight-line fit should not be extrapolated across a phase transition, a change in defect regime, or a major surface reconstruction.
9. Reporting results so they can be reproduced
A useful report should include more than a single value of \(D\) or \(k\).
Material and specimen - Exact composition, synthesis route, and phase purity. - Density, porosity, grain size, and surface preparation. - Specimen dimensions and exposed geometry. - Electrode architecture for porous-electrode measurements.
Measurement conditions - Temperature and oxygen partial pressure. - Gas composition, flow rate, and water content. - Equilibration time and gas-switching procedure. - Measurement duration and sampling interval. - For isotope exchange: isotope enrichment, exposure time, and quenching protocol.
Data analysis - Whether the reported parameters are tracer or chemical coefficients. - Governing model, boundary conditions, and geometry. - Fitting interval, uncertainty, and parameter covariance where available. - Treatment of gas-switching or instrument response. - Whether the perturbation was within the linear-response regime.
A value reported without its measurement method and conditions is difficult to compare across laboratories. In particular, \(k^*\) should not be compared directly with \(k_{chem}\) without explaining the distinction.
10. A practical strategy for evaluating a new oxide electrode
A defensible evaluation sequence is:
- Confirm phase and microstructure. Use XRD and microscopy to establish phase purity, density, grain size, and porosity.
- Establish defect chemistry. Measure oxygen non-stoichiometry or conductivity as a function of \(p_{O_2}\) and temperature where possible.
- Measure bulk and surface kinetics. Use ECR for chemical coefficients and isotope exchange for tracer coefficients, selecting the method appropriate to the material.
- Evaluate electrode response. Use symmetrical-cell EIS or full-cell measurements under controlled gas conditions.
- Examine the surface after testing. Look for segregation, secondary phases, contamination, and morphological changes.
- Separate intrinsic material properties from electrode architecture. Report dense-specimen transport separately from porous-electrode performance.
This sequence avoids treating one conductivity value or one impedance arc as a complete description of oxygen electrode activity.
11. Key takeaways
- Oxygen electrode performance involves coupled gas transport, surface reaction, charge transfer, and solid-state oxygen transport.
- \(D_{chem}\) describes chemical relaxation within the oxide; \(k_{chem}\) describes the effective surface-exchange response.
- Tracer parameters \(D^*\) and \(k^*\) obtained by isotope exchange are not automatically equal to their chemical counterparts.
- ECR, isotope exchange, and EIS probe related but distinct aspects of oxygen transport and electrode kinetics.
- Surface chemistry, defects, microstructure, gas composition, and ageing can all change measured kinetics.
- Reliable comparisons require consistent definitions, experimental conditions, geometry, and data-analysis models.
References
- S. B. Adler, “Factors Governing Oxygen Reduction in Solid Oxide Fuel Cell Cathodes,” Chemical Reviews 104 (2004), 4791–4843. https://doi.org/10.1021/cr020724o
- J. A. Kilner and M. Burriel, “Materials for Intermediate-Temperature Solid-Oxide Fuel Cells,” Annual Review of Materials Research 44 (2014), 365–393. https://doi.org/10.1146/annurev-matsci-070813-113426
- “The oxygen reduction reaction in solid oxide fuel cells: from kinetic parameters measurements to electrode design,” Journal of Physics: Energy (2020). https://doi.org/10.1088/2515-7655/abb4ec
- V. Thoréton, M. Niania, and J. Kilner, “Kinetics of competing exchange of oxygen and water at the surface of functional oxides,” Physical Chemistry Chemical Physics (2021). https://doi.org/10.1039/D0CP04953K
- S. B. Adler, X. Y. Chen, and J. R. Wilson, “Mechanisms and rate laws for oxygen exchange on mixed-conducting oxide surfaces,” Journal of Catalysis 245 (2007), 91–109. https://doi.org/10.1016/j.jcat.2006.09.019
- P. Fielitz and G. Borchardt, “On the accurate measurement of oxygen self-diffusivities and surface exchange coefficients in oxides via SIMS depth profiling,” Solid State Ionics 144 (2001), 71–80. https://doi.org/10.1016/S0167-2738(01)00893-1