Hydration Thermodynamics in Proton-Conducting Oxides: Water Incorporation, Equilibrium, and Temperature Dependence
Hydration Thermodynamics in Proton-Conducting Oxides: Water Incorporation, Equilibrium, and Temperature Dependence
Proton-conducting oxides take up water from the surrounding atmosphere by incorporating hydrogen into the oxygen sublattice. This process, commonly called hydration, creates protonic defects that can support ionic conduction at intermediate temperatures.
Hydration is not simply a matter of exposing a ceramic to humid gas. The equilibrium proton concentration is governed by the Gibbs energy of the hydration reaction, the concentration of oxygen vacancies, temperature, and water-vapour activity. The resulting proton population can therefore change substantially with temperature and humidity, even when the nominal dopant concentration remains fixed.
This article develops the thermodynamic description of hydration in acceptor-doped oxides, derives the mass-action relationship, explains the roles of hydration enthalpy and entropy, and connects equilibrium water uptake to proton concentration and experimental measurements.
1. The hydration reaction in an oxide lattice
In Kröger–Vink notation, the conventional hydration reaction is
[ ]
An oxygen vacancy is filled by an oxygen atom originating from water, while hydrogen is incorporated as a proton bound to a lattice oxygen. Two hydroxyl-type protonic defects are produced for each water molecule incorporated.
The effective charges are balanced:
- (V_O^{}): charge (+2).
- (O_O^x): charge (0).
- (2OH_O^): total charge (+2).
The reaction is therefore charge balanced in the Kröger–Vink formalism.
For an acceptor-doped perovskite such as a generic (), a lower-valence dopant on the B site is commonly represented as an acceptor defect (M_B’). Under conditions where oxygen vacancies compensate the acceptors, the dry material contains a population of oxygen vacancies. Hydration converts some of those vacancies into protonic defects.
The exact defect chemistry depends on dopant valence, site occupancy, secondary phases, and whether electronic defects also contribute to charge compensation.
2. Mass-action expression for hydration
For the hydration reaction,
[ ]
the thermodynamic equilibrium constant is
[ K_{} = ]
where each (a_i) is a dimensionless activity.
For a dilute defect solution, if the regular oxygen-site activity is approximately constant and defect activities can be approximated by normalized site fractions, the expression becomes
[ K_{}’ = ]
Rearranging,
[ ]
and therefore
[ ]
This relation shows that proton concentration increases with water activity and with the availability of oxygen vacancies, provided the hydration equilibrium constant and the assumptions of the dilute model remain applicable.
For an ideal water-vapour mixture,
[ a_{} ]
where (p^) is the standard pressure. At non-ideal conditions, fugacity should be used instead of partial pressure.
3. Hydration Gibbs energy and the equilibrium constant
The standard Gibbs energy change of hydration is related to the equilibrium constant by
[ ]
The reaction Gibbs energy under non-standard conditions is
[ rG{} = rG{}^ + RTQ_{} ]
where the reaction quotient is
[ Q_{} = ]
At equilibrium,
[ rG{}=0 ]
and therefore
[ Q_{}=K_{} ]
A negative standard hydration Gibbs energy means that hydration is favourable under standard-state conditions. It does not, by itself, mean that all oxygen vacancies will be hydrated: the actual equilibrium depends on temperature, water activity, defect activities, and the available vacancy population.
4. Hydration enthalpy and hydration entropy
The standard Gibbs energy is
[ rG{}^ = rH{}^ - TrS{}^ ]
The hydration enthalpy describes the energetic contribution associated with incorporating water-derived species into the oxide lattice. The hydration entropy includes changes in configurational, vibrational, and gas-phase contributions.
Hydration of many acceptor-doped perovskite oxides is exothermic:
[ rH{}^<0 ]
The entropy change is also often negative because gaseous water is consumed and hydrogen becomes incorporated into the solid lattice:
[ rS{}^<0 ]
These signs are common tendencies, not universal constants. Their values depend on composition, structure, defect association, and the thermodynamic reference states used.
When both enthalpy and entropy are negative, hydration is generally favoured at lower temperature and becomes less favourable as temperature increases, all else being equal. This is one reason proton concentrations in many oxide proton conductors decrease on heating in a fixed humid atmosphere.
5. The van ’t Hoff relationship for hydration
Starting from
[ K_{} = - ]
and using
[ rG{}^ = rH{}^ - TrS{}^ ]
gives
[ K_{} = - + ]
If the reaction enthalpy and entropy are approximately temperature independent over the range considered,
[ ]
For an exothermic hydration reaction, the slope of (K_{}) versus (1/T) is positive.
The integrated form between two temperatures is
[ () = - ( - ) ]
A plot of (K_{}) against (1/T) can therefore be used to estimate the hydration enthalpy, provided the measured or calculated equilibrium constants correspond to the same reaction and standard-state convention.
6. How water-vapour pressure controls proton concentration
For the simplified dilute model,
[ [OH_O^]^2 = K_{}’ [V_O^{}]a_{} ]
If the oxygen-vacancy concentration is approximately constant, then
[ ]
and, for an ideal gas mixture,
[ [OH_O^] p_{}^{1/2} ]
This square-root dependence is a useful limiting result. It assumes that the vacancy population is effectively fixed and that protonic defects remain dilute.
It should not be applied indiscriminately. When hydration substantially depletes the vacancies, the vacancy concentration is no longer constant and the pressure dependence changes.
7. Acceptor-doped material: charge neutrality and hydration extent
Consider a simplified acceptor-doped oxide in which the dominant charged defects are acceptors, oxygen vacancies, and protonic defects. Neglecting electronic carriers and other defects, charge neutrality is
[ 2[V_O^{}] + [OH_O^] = [A’] ]
where ([A’]) is the fixed acceptor concentration.
This relation expresses two limiting compensation mechanisms:
- In a dry state, oxygen vacancies compensate most of the acceptor charge.
- In a highly hydrated state, protonic defects replace a substantial fraction of the vacancy compensation.
Let the total acceptor concentration be fixed and define (v=[V_O^{}]), (h=[OH_O^]), and (a=[A’]). Then
[ 2v+h=a ]
so
[ v= ]
Substituting into the simplified hydration equilibrium,
[ h^2 = K_{}’a_{} ]
or
[ ]
The physically meaningful positive solution is
[ ]
This result describes the transition from a vacancy-dominated dry state to a proton-rich hydrated state within the simplified defect model.
At low water activity, (ha), and the vacancy concentration is approximately (a/2). The square-root pressure dependence is recovered:
[ h ( )^{1/2} ]
At sufficiently high water activity, the proton concentration approaches the acceptor concentration:
[ ha ]
and the oxygen-vacancy concentration approaches zero in this idealized model. Real materials generally retain additional defects and may not reach this mathematical limit.
8. A practical numerical illustration
Suppose an acceptor-doped oxide is described by the simplified relation
[ h^2=K_{}’v a_{} ]
At a selected temperature, take illustrative normalized values:
[ K_{}’=0.20, v=0.020, a_{}=0.10 ]
Then
[ h^2=(0.20)(0.020)(0.10) =4.0^{-4} ]
and
[ ]
The calculated proton site fraction is 0.020 under the stated simplified model.
If water activity is increased by a factor of four while vacancy concentration and temperature are held fixed,
[ h_{} = h_{} ]
so the predicted proton concentration doubles. This scaling applies only while the vacancy concentration remains effectively unchanged.
The values above are illustrative, not reported data for a particular composition.
9. Water activity, humidity, and gas composition
Water-vapour pressure is often controlled by mixing a humidified gas with a dry carrier gas. For an ideal gas mixture,
[ p_{}=y_{}P_{} ]
where (y_{}) is the water-vapour mole fraction and (P_{}) is the total pressure.
Relative humidity is defined with respect to the saturation vapour pressure at the relevant temperature:
[ = {p_{}^{}(T)} ]
Thus, the same relative humidity at two different temperatures does not correspond to the same water-vapour partial pressure. When comparing hydration experiments, report the actual water-vapour pressure or dew point, total pressure, and sample temperature rather than relative humidity alone.
In a gas line, the water-vapour pressure at the sample may differ from the value expected from the humidifier if tubing temperatures, leaks, condensation, or flow conditions are not controlled.
10. Hydration and oxygen nonstoichiometry
For a perovskite written as
[ ]
the oxygen deficiency () is related to the oxygen-vacancy population. Hydration consumes oxygen vacancies and introduces protonic defects.
In the simplified acceptor-compensated model,
[ 2[V_O^{}]+[OH_O^]=[A’] ]
Therefore, increasing proton concentration is accompanied by a decrease in the oxygen-vacancy concentration.
The oxygen content measured under humid conditions may consequently differ from that measured under dry conditions. A mass change during hydration can contain contributions from water incorporation and, depending on the atmosphere and material, changes in oxygen content.
For quantitative interpretation, the composition model must distinguish oxygen vacancies from protonic defects. A single mass change should not automatically be assigned entirely to oxygen uptake or entirely to hydrogen uptake without a balanced reaction and a defined reference state.
See Measuring Oxygen Nonstoichiometry in Functional Oxides for the experimental basis of oxygen-content determination.
11. Hydration, proton mobility, and conductivity
The concentration of protonic defects is not the same as proton conductivity. A simplified expression for protonic conductivity is
[ _{H} = ]
where (n_{H}) is the number density of mobile protons, (q) is the elementary charge, (D_{H}) is the proton diffusion coefficient, and (k_{B}) is Boltzmann’s constant.
Hydration can increase (n_{H}), but increasing temperature can simultaneously increase proton mobility while decreasing the equilibrium proton concentration in an exothermically hydrated oxide.
As a result, proton conductivity can exhibit a temperature dependence that reflects competition between carrier concentration and mobility. It may increase with temperature over one interval even while the equilibrium proton concentration decreases.
The Nernst–Einstein expression is a useful idealized connection between diffusion and conductivity. Correlations between proton motions, defect association, and microstructural effects can cause deviations from this simple relation.
12. Hydration kinetics are not hydration equilibrium
Thermodynamic equilibrium specifies the final defect populations under defined conditions. It does not specify how quickly the material reaches those populations.
The observed hydration rate can be limited by:
- Dissociative adsorption of water at the surface.
- Surface exchange and incorporation of hydrogen.
- Proton diffusion through the bulk.
- Oxygen-vacancy transport and redistribution.
- Grain-boundary resistance.
- Gas transport and humidity equilibration in the measurement system.
A sample may therefore show a slow mass or conductivity response after a humidity step even when the equilibrium state is thermodynamically well defined.
To distinguish equilibrium from kinetics, vary the dwell time, check whether the response reaches a stable plateau, and compare humidification and dehydration paths. Hysteresis can arise from slow equilibration, defect association, phase changes, or irreversible chemical changes.
13. Experimental approaches to hydration thermodynamics
Several complementary methods can be used to investigate hydration:
Thermogravimetric analysis (TGA). Measures mass changes during controlled changes in water-vapour pressure and temperature. With a defined composition and reaction stoichiometry, the mass response can be converted into water uptake.
Electrical conductivity. Tracks changes in transport properties as humidity and temperature are varied. Conductivity alone does not uniquely determine proton concentration because mobility and other carrier contributions also change.
Isotope exchange. Replacing () with () can help distinguish proton-related processes and investigate isotope effects.
Spectroscopic methods. Infrared and related techniques can provide evidence for hydroxyl species and changes in hydrogen bonding or local environments.
Electrochemical measurements. Under suitable cell configurations, transport numbers and concentration-cell responses can help separate ionic and electronic contributions.
A robust thermodynamic analysis combines more than one technique and uses a consistent defect model to interpret the measurements.
14. Common pitfalls
- Using relative humidity without its temperature. Relative humidity does not uniquely specify water-vapour pressure unless the temperature is known.
- Assuming vacancy concentration is always constant. Hydration consumes vacancies; the fixed-vacancy approximation is a limiting regime.
- Confusing proton concentration with proton conductivity. Mobility and microstructure also control conductivity.
- Ignoring electronic defects. At some oxygen chemical potentials, electrons or holes contribute to charge compensation.
- Treating hydration enthalpy as universal. It depends on composition, defect interactions, and the thermodynamic model.
- Extracting equilibrium constants from transient data. Kinetic lag can distort apparent equilibrium values.
- Ignoring gas-line condensation or leaks. The atmosphere at the sample may differ from the nominal humidifier setting.
- Comparing inconsistent standard states. Equilibrium constants and thermodynamic quantities must use compatible activity conventions.
15. Practical workflow for analysing hydration
- Define the oxide composition, acceptor concentration, and relevant defect species.
- Write the balanced hydration reaction in Kröger–Vink notation.
- Specify the activity convention and standard states.
- Derive the mass-action expression.
- Write charge neutrality, including electronic defects if they are relevant.
- Convert measured humidity, dew point, or gas composition into water-vapour activity at the sample.
- Determine equilibrium proton uptake at each temperature and water activity.
- Separate equilibrium carrier concentration from mobility when interpreting conductivity.
- Use temperature-dependent equilibrium data to estimate hydration enthalpy and entropy.
- Check reversibility, equilibration time, and possible phase or chemical changes.
16. Summary
Hydration in proton-conducting oxides is described by
[ ]
and, in a simplified dilute-defect model,
[ [OH_O^]^2 = K_{}’ [V_O^{}]a_{} ]
The equilibrium constant is connected to the standard hydration Gibbs energy by
[ rG{}^ = -RTK_{} ]
Hydration enthalpy and entropy determine how the equilibrium changes with temperature. Water-vapour activity controls the reaction quotient, while charge neutrality determines how protonic defects replace oxygen vacancies.
The central distinction is between thermodynamic proton uptake and proton transport: hydration sets the equilibrium population of protonic defects, whereas conductivity also depends on their mobility, correlations, and the microstructure of the ceramic.
Further reading
- K. D. Kreuer, “Proton-conducting oxides,” Annual Review of Materials Research 33, 333–359 (2003).
- H. Iwahara, “Technological challenges in the application of proton conducting ceramics,” Solid State Ionics 77, 289–298 (1995).
- J. Maier, Physical Chemistry of Ionic Materials: Ions and Electrons in Solids, Wiley.
- A. R. West, Solid State Chemistry and its Applications, 2nd ed., Wiley.