Thermal Expansion in Functional Oxides: Dilatometry, Anisotropy, and Cell Compatibility
Thermal expansion is a basic but consequential property of functional ceramics. In a single-phase ceramic, it determines how dimensions change with temperature. In a multilayer electrochemical cell, differences in thermal expansion between the electrolyte, electrodes, interconnects, and seals can generate residual stresses during fabrication, start-up, operation, and cooling. These stresses may cause cracking, warping, loss of contact, or delamination.
A thermal expansion coefficient (TEC) is therefore more than a tabulated material constant. Its value depends on temperature, composition, crystal symmetry, defect chemistry, atmosphere, and measurement direction. This guide explains how to measure and interpret thermal expansion and how to use it in materials selection for solid oxide and protonic ceramic cells.
1. What does thermal expansion measure?
For a specimen with initial length (L_0), the linear thermal strain is
\[ \varepsilon_{\mathrm{th}}= \frac{L(T)-L_0}{L_0} \]
The average linear thermal expansion coefficient over a temperature interval is
\[ \overline{\alpha}_L= \frac{L(T_2)-L(T_1)} {L(T_1)(T_2-T_1)} \]
For small strains, this is often approximated as
\[ \overline{\alpha}_L\approx \frac{\Delta L/L}{\Delta T} \]
The instantaneous or differential coefficient is
\[ \alpha_L(T)= \frac{1}{L(T)} \left(\frac{dL}{dT}\right)_p \]
where the subscript indicates measurement at the specified pressure. The distinction matters when the coefficient changes appreciably with temperature. A single average value can conceal curvature, phase transitions, or changes in defect chemistry.
Thermal expansion coefficients are commonly reported in () or (10{-6} \mathrm{K{-1}}). Temperature differences in kelvin and degrees Celsius have the same numerical magnitude.
Linear and volumetric expansion
For an isotropic material with small strains,
\[ \beta\approx3\alpha_L \]
where () is the volumetric thermal expansion coefficient. This approximation is not appropriate for strongly anisotropic crystals unless the directional coefficients are accounted for.
For orthorhombic or lower-symmetry materials, the volume coefficient is approximately the sum of the three principal linear coefficients:
\[ \beta=\alpha_a+\alpha_b+\alpha_c \]
This relation is useful when interpreting crystallographic lattice-parameter changes measured by high-temperature X-ray diffraction.
2. Why oxide ceramics expand
Thermal expansion in a crystalline solid arises from the anharmonicity of interatomic interactions. As temperature increases, the average interatomic separation changes because the interatomic potential is not perfectly symmetric around its minimum.
In functional oxides, the observed expansion can also be influenced by:
- changes in lattice vibrations and phonon populations;
- structural phase transitions or octahedral tilting;
- cation substitution and changes in bond stiffness;
- oxygen-vacancy concentration and cation valence;
- hydration and dehydration;
- reduction and oxidation under different oxygen partial pressures;
- microcracking, porosity, and irreversible dimensional changes.
A measured dilatometric curve is consequently a response of the complete specimen under the measurement conditions. It is not always a pure measure of the intrinsic lattice expansion.
3. Thermal expansion anisotropy
In cubic crystals, symmetry requires equal principal linear expansion coefficients in the absence of external fields or texture effects. In tetragonal, orthorhombic, hexagonal, and lower-symmetry crystals, expansion can differ along crystallographic directions.
For a tetragonal crystal,
\[ \alpha_a=\alpha_b\ne\alpha_c \]
For a general orthorhombic crystal,
\[ \alpha_a\ne\alpha_b\ne\alpha_c \]
The volume expansion is related to the sum of the principal coefficients, but a polycrystalline ceramic may show an effective macroscopic response that depends on grain orientation distribution, texture, porosity, and mechanical constraint between grains.
Single crystal versus polycrystalline ceramic
A single-crystal measurement can resolve directional expansion. A randomly oriented, dense polycrystalline specimen often gives an averaged macroscopic response, provided it is sufficiently homogeneous and does not develop microcracks. Strong texture can make the polycrystalline response direction-dependent.
When comparing literature values, check whether the reported value is:
- a single-crystal coefficient along a specified axis;
- a polycrystalline linear coefficient;
- a volumetric coefficient;
- an average over a stated temperature interval;
- measured during heating or cooling;
- obtained in air, inert gas, vacuum, or a controlled oxygen partial pressure.
Values without these details are difficult to compare reliably.
4. Measuring thermal expansion by dilatometry
Push-rod dilatometry is a common technique for measuring the dimensional change of a ceramic specimen as a function of temperature. The specimen is placed in a furnace and its displacement is transmitted to a sensor through a push rod. The instrument records displacement relative to a reference system.
A typical workflow is:
- Prepare a specimen with parallel faces and a well-defined measurement length.
- Measure the initial length and record the specimen geometry.
- Select a holder and push-rod material compatible with the temperature and atmosphere.
- Run a baseline or reference measurement using the instrument’s calibration procedure.
- Apply the selected heating rate and atmosphere.
- Record displacement and temperature during heating and, where relevant, cooling.
- Correct for the thermal expansion of the instrument components.
- Calculate strain and the average or differential expansion coefficient.
- Repeat measurements to assess reproducibility and irreversible changes.
Instrument correction
The raw displacement is not necessarily the specimen displacement. It can include expansion of the push rod, sample holder, and other parts of the measurement train. The instrument correction must be performed according to the dilatometer calibration protocol.
If the corrected specimen displacement is (L_{s}), the linear strain is
\[ \varepsilon_L(T)= \frac{\Delta L_{\mathrm s}(T)}{L_0} \]
A numerical derivative gives an estimate of the instantaneous coefficient:
\[ \alpha_L(T)= \frac{d\varepsilon_L}{dT} \]
Numerical differentiation amplifies noise. Apply a justified smoothing or fitting procedure, report the method, and avoid interpreting small derivative features that are not reproducible in the original displacement data.
Heating-rate effects
A higher heating rate can introduce thermal gradients between the specimen and the thermocouple or within the specimen itself. It can also shift the apparent temperature of a kinetic process. For accurate comparisons, use the same heating rate, sample geometry, atmosphere, and calibration procedure.
For materials that undergo redox or hydration changes, the heating rate and gas-flow conditions can influence the measured curve because chemical equilibration may take time.
5. Thermal expansion and atmosphere-dependent chemistry
Many functional oxides are non-stoichiometric. Their oxygen content and transition-metal oxidation states can change with oxygen partial pressure and temperature. Such changes may alter lattice dimensions in addition to ordinary thermal expansion.
The measured strain can be considered schematically as
\[ \varepsilon_{\mathrm{meas}}= \varepsilon_{\mathrm{thermal}}+ \varepsilon_{\mathrm{chemical}}+ \varepsilon_{\mathrm{irreversible}} \]
This decomposition is conceptual: the contributions can be coupled and are not always independently additive over large strain ranges.
For example, reduction may create oxygen vacancies and change the average cation valence. The associated chemical expansion can either reinforce or oppose the ordinary thermal expansion. If the oxygen content changes during a dilatometry run, the slope of length versus temperature is an apparent coefficient under that atmosphere and protocol, not necessarily the intrinsic TEC of a fixed-composition phase.
To distinguish the contributions, consider:
- measuring under controlled oxygen partial pressure;
- comparing oxidizing and reducing atmospheres;
- performing heating and cooling cycles;
- checking mass change by thermogravimetry;
- measuring lattice parameters by high-temperature XRD;
- allowing sufficient time for chemical equilibration at each temperature.
6. Thermal expansion mismatch in multilayer cells
A solid oxide cell contains multiple bonded or mechanically constrained layers. If two layers have different free thermal strains, their interface must accommodate the mismatch.
For two materials with constant coefficients over a temperature interval,
\[ \Delta\varepsilon_{\mathrm{mis}}= (\alpha_1-\alpha_2)\Delta T \]
For temperature-dependent coefficients, use the integrated difference:
\[ \Delta\varepsilon_{\mathrm{mis}}= \int_{T_1}^{T_2} [\alpha_1(T)-\alpha_2(T)]\,dT \]
This expression estimates the difference in free thermal strain. It is not itself the stress. The actual stress depends on elastic moduli, Poisson’s ratios, layer thicknesses, geometry, porosity, creep, cracking, and mechanical boundary conditions.
A simple mismatch example
Suppose an electrolyte has an average TEC of (10.5{-6} \mathrm{K{-1}}), while an electrode has an average TEC of (14.0{-6} \mathrm{K{-1}}), over a temperature change of 700 K.
The free-strain mismatch is
\[ \Delta\varepsilon_{\mathrm{mis}}= (14.0-10.5)\times10^{-6}\times700 \]
\[ \boxed{\Delta\varepsilon_{\mathrm{mis}}=2.45\times10^{-3}} \]
This corresponds to 0.245% differential free strain. Whether it causes cracking or delamination cannot be decided from the TEC difference alone; the layer architecture, elastic response, interfaces, and stress relaxation must also be considered.
7. From mismatch strain to stress: limits of a simple estimate
For a fully constrained, isotropic elastic material under a simplified uniaxial condition, a stress scale can be estimated as
\[ \sigma\sim E\,\Delta\varepsilon \]
where (E) is Young’s modulus. This is an order-of-magnitude estimate, not a general multilayer solution. A biaxially constrained thin film may require a plane-stress or plane-strain formulation involving Poisson’s ratio. A porous electrode or a high-temperature ceramic that creeps can relax stress substantially.
For an idealized equibiaxial thin film under plane stress, the elastic stress scale is often written
\[ \sigma\approx \frac{E}{1-\nu}\,\Delta\varepsilon \]
where () is Poisson’s ratio. The sign depends on which layer expands more and on whether the temperature is increasing or decreasing. Use finite-element analysis or an appropriate laminate model when designing a real cell stack.
Why thermal cycling is especially demanding
A cell may survive one heating cycle but degrade after repeated cycling. Repeated mismatch stress can initiate or propagate cracks, weaken interfaces, and cause local loss of contact. Start-up and shut-down rates also matter because transient temperature gradients can produce stresses beyond those predicted from uniform-temperature expansion alone.
8. Strategies for improving thermal compatibility
Thermal compatibility can be improved through materials selection and architecture, not only by searching for identical TEC values.
- Composition tuning: substitution can adjust lattice dimensions, phase stability, and thermal expansion, but may also alter conductivity and chemical stability.
- Composite electrodes: mixing an electrode with an electrolyte phase can adjust the effective expansion and improve interfacial contact.
- Graded interfaces: gradual changes in composition or porosity can reduce abrupt property discontinuities.
- Layer-thickness optimization: the stiffness and thickness of each layer influence stress distribution.
- Controlled porosity: porous electrodes can have a lower effective modulus and may accommodate some strain, although excessive porosity harms electrical and mechanical connectivity.
- Compatible processing: co-sintering shrinkage, phase reactions, and residual stresses must be considered alongside TEC.
- Thermal management: controlled start-up and shut-down protocols reduce transient gradients and thermal shock.
A close TEC match is helpful but not sufficient. Chemical compatibility, interdiffusion, mechanical integrity, and electrochemical performance must be evaluated together.
9. Chemical expansion is not the same as thermal expansion
Thermal expansion is the dimensional response to temperature change under defined chemical conditions. Chemical expansion is the dimensional response to a change in composition or defect concentration, such as oxygen non-stoichiometry or hydration.
In a working electrochemical cell, both may occur simultaneously. For example, an electrode can experience a temperature gradient and a spatially varying oxygen chemical potential. Its total strain may then contain thermal, chemical, and mechanical contributions.
Do not report an atmosphere-dependent length change as a pure TEC unless the chemical state is known to remain constant or the chemical contribution has been separated.
10. Complementary characterization
Dilatometry provides a macroscopic length change, but complementary measurements help identify its origin.
- High-temperature XRD: tracks lattice parameters and phase transitions; it measures crystallographic expansion rather than the full macroscopic response of a porous body.
- Thermogravimetric analysis: measures mass changes associated with oxygen exchange, hydration, or other chemical processes.
- DSC/DTA: detects thermal events and phase transitions, but does not directly replace a length-change measurement.
- Elastic-modulus measurements: provide mechanical parameters needed for stress analysis.
- SEM and optical microscopy: reveal cracks, warping, and microstructural changes after thermal cycling.
- In situ diffraction or microscopy: can follow structural evolution under operating temperature and atmosphere.
The most useful combination depends on whether the question concerns intrinsic lattice expansion, chemical expansion, component compatibility, or long-term mechanical degradation.
11. Practical reporting checklist
For a reproducible thermal-expansion study, report:
- Exact composition, phase purity, and specimen processing history.
- Specimen dimensions, density, and porosity where relevant.
- Measurement technique and instrument calibration procedure.
- Heating and cooling rates and the temperature interval.
- Atmosphere, gas flow, oxygen partial pressure, and humidity.
- Whether the result is an average or differential coefficient.
- Whether the value is linear, volumetric, or direction-specific.
- Number of cycles and evidence of reversibility.
- Corrections for holder and push-rod expansion.
- Uncertainty, repeatability, and any evidence of mass or phase changes.
12. Common interpretation errors
- Comparing coefficients measured over different temperature ranges. TEC can vary with temperature.
- Mixing linear and volumetric coefficients. Their numerical values are not interchangeable.
- Ignoring crystal anisotropy. A single directional value may not represent a textured ceramic.
- Ignoring atmosphere. Redox and hydration can change lattice dimensions.
- Treating raw push-rod displacement as specimen expansion. Instrument correction is essential.
- Differentiating noisy data without validation. Numerical derivatives can exaggerate artefacts.
- Assuming TEC mismatch directly gives stress. Elastic constraints and geometry determine the stress state.
- Assuming a close TEC match guarantees compatibility. Chemical reactions, sintering shrinkage, and interdiffusion can still cause failure.
Further reading
- Y. Shen et al., “Bridging the thermal expansion gap in solid oxide fuel cells: towards robust and efficient energy conversion,” Journal of Materials Chemistry A 14 (2026), 5542–5563. https://doi.org/10.1039/D5TA08407E.
- “Thermal properties of solid oxide cells and their governing physicochemical factors: a review,” Journal of Power Sources 667 (2026), 239230. https://doi.org/10.1016/j.jpowsour.2025.239230.
- N. Shah et al., “Mitigating thermal expansion effects in solid oxide fuel cell cathodes: A critical review,” Journal of Power Sources 599 (2024), 234211. https://doi.org/10.1016/j.jpowsour.2024.234211.
- “Solid oxide electrolysis cells – current material development and industrial application,” Journal of Materials Chemistry A (2023). https://doi.org/10.1039/D3TA02161K.
- “Addressing planar solid oxide cell degradation mechanisms: A critical review of selected components,” Electrochemical Science Advances (2022). https://doi.org/10.1002/elsa.202100024.