Chemical Potential and Activity in Solid Materials: Oxygen Potential, Gas Equilibria, and Defect Chemistry
Chemical Potential and Activity in Solid Materials: Oxygen Potential, Gas Equilibria, and Defect Chemistry
Chemical potential is the thermodynamic quantity that determines how a component responds to a change in composition and whether it tends to enter, leave, or redistribute within a material. It connects phase equilibrium, gas–solid reactions, defect populations, oxygen nonstoichiometry, hydration, and electrochemical processes within one framework.
For functional oxides, chemical potential is particularly useful because the surrounding atmosphere imposes chemical conditions on the solid. Changing temperature, oxygen partial pressure, or water-vapour pressure changes the chemical potentials of the gas species and can therefore shift the equilibrium concentrations of oxygen vacancies, electronic carriers, and protonic defects.
This article develops the concepts needed to use chemical potential quantitatively, with particular attention to activity, standard states, gas fugacity, and oxygen chemical potential.
1. Definition of chemical potential
For a multicomponent system, the Gibbs free energy is a function of temperature, pressure, and the amounts of each component:
[ G=G(T,P,n_1,n_2,,n_i) ]
The chemical potential of component (i) is its partial molar Gibbs free energy:
[ i= (){T,P,n_{ji}} ]
It measures the change in Gibbs free energy when an infinitesimal amount of component (i) is added while temperature, pressure, and the amounts of all other components are held constant.
For a system with variable composition,
[ dG=-S,dT+V,dP+_i_i,dn_i ]
At constant temperature and pressure,
[ dG=_i_i,dn_i ]
This expression is central to chemical equilibrium. A spontaneous compositional change proceeds in the direction that lowers the total Gibbs free energy. At equilibrium, the free-energy change associated with every allowed independent reaction or transfer is zero.
Chemical potential is an intensive quantity. It is not the same as the total Gibbs energy or the chemical energy stored in an entire sample.
2. Chemical potential in a mixture
For an ideal mixture, the chemical potential of component (i) is
[ _i=_i^+RTx_i ]
where (x_i) is its mole fraction and (_i^) is the chemical potential in the selected standard state at the same temperature and pressure.
For a non-ideal mixture, composition alone does not adequately describe the thermodynamic behaviour. The expression is written using activity:
[ ]
Activity (a_i) is dimensionless. It represents the effective thermodynamic concentration of a species relative to a defined standard state.
For a solution described using mole fractions,
[ a_i=_i x_i ]
where (_i) is the activity coefficient on the chosen mole-fraction standard-state convention. Thus,
[ _i=_i^+RT(_i x_i) ]
The ideal-mixture limit is (_i=1), in which case (a_i=x_i).
The activity coefficient accounts for non-ideal interactions. Depending on the system, these may arise from differences in bonding, local strain, electrostatic interactions, short-range ordering, or changes in the local coordination environment.
Important: activity coefficients depend on the chosen convention and standard state. A numerical activity is meaningful only when its reference convention is specified.
3. Activity and standard states
A standard state is a defined reference state used to express chemical potential. It is not necessarily the most stable state under the actual experimental conditions.
Common conventions include:
| System | Common activity convention |
|---|---|
| Ideal dilute solute | Concentration relative to a reference concentration |
| Liquid or solid solution | (a_i=_i x_i), under a stated mole-fraction convention |
| Pure solid or pure liquid | (a_i=1) in its standard state |
| Ideal gas | (a_i=p_i/p^) |
| Real gas | (a_i=f_i/p^) |
Here (p^) is the standard pressure, conventionally 1 bar, and (f_i) is fugacity.
For a pure crystalline solid, assigning unit activity is often a useful approximation when the solid is treated as a separate pure phase with fixed composition. It is not generally valid for a component dissolved in a solid solution, a non-stoichiometric oxide, or a phase whose composition changes appreciably.
For example, the activity of oxygen in a nearly stoichiometric pure oxide phase may be incorporated into a reaction equilibrium constant. By contrast, the oxygen-vacancy population in a variable-composition oxide must be represented through an appropriate defect activity or a clearly stated concentration approximation.
4. Chemical potential of gases: pressure and fugacity
For an ideal gas, the chemical potential is
[ _i(T,p_i)=_i(T)+RT(\frac{p_i}{p}) ]
The pressure ratio is dimensionless. A logarithm of a dimensional pressure should not be used.
Real gases deviate from ideal-gas behaviour, particularly at high pressure. Fugacity (f_i) replaces partial pressure:
[ _i(T,p_i)=_i(T)+RT(\frac{f_i}{p}) ]
Fugacity is related to pressure through the fugacity coefficient (_i):
[ f_i=_i p_i ]
At sufficiently low pressure, (_i) approaches unity and fugacity approaches partial pressure. For many oxide experiments conducted near atmospheric pressure, the ideal-gas approximation is often adequate, but the assumption should be checked when high-pressure conditions are used.
5. Oxygen chemical potential and oxygen partial pressure
Oxygen chemical potential is one of the most important thermodynamic variables in oxide materials research. For molecular oxygen in an ideal gas atmosphere,
[ {}(T,p{}) = _{}^(T) + RT() ]
Defect reactions are frequently written per oxygen atom rather than per oxygen molecule. The corresponding chemical potential per mole of oxygen atoms is one-half the molecular value:
[ {O} = {} ]
Therefore,
[ {O} = {}^(T) + RT() ]
This distinction matters when comparing expressions for oxygen incorporation or oxygen loss. The factor of one-half comes from the stoichiometry of molecular oxygen.
5.1 How a change in oxygen pressure shifts oxygen potential
At fixed temperature, the change in oxygen chemical potential per mole of oxygen atoms between two atmospheres is
[ _{O} = RT () ]
A decrease in (p_{}) lowers the oxygen chemical potential and makes reducing conditions more favourable. An increase in (p_{}) raises oxygen chemical potential and favours oxidation.
For a tenfold change in oxygen partial pressure,
[ _{O} = RT ]
At (T=1000 K), using (R=8.314 ),
[ _{O} ]
This is the change per mole of oxygen atoms, not per mole of () molecules. The corresponding molecular oxygen chemical-potential change is twice this value.
6. Connecting oxygen chemical potential to defect equilibria
Consider oxygen loss from an oxide lattice:
[ ]
The thermodynamic equilibrium constant is written in terms of activities:
[ K_{} = ]
If regular oxygen sites are treated as the reference population and the oxide is sufficiently close to stoichiometric that (a_{}) is approximately constant, that term can be absorbed into an effective equilibrium constant:
[ K_{}’ a_{} a_{}^2 ()^{1/2} ]
In a dilute-defect approximation, activities may be replaced by appropriately normalized defect concentrations. Under those assumptions,
[ K_{}’ ^2 p_{}^{1/2} ]
This is the thermodynamic basis for oxygen-pressure-dependent defect concentrations and Brouwer-diagram slopes. The resulting concentration–pressure relation still depends on charge neutrality, dopant content, site constraints, and which defects dominate.
For example, in an idealized intrinsic regime where electrons compensate oxygen vacancies,
[ [e’]=2[] ]
Combining this relation with the dilute mass-action expression gives
[ []^3 p_{}^{-1/2} ]
and hence
[ [] p_{}^{-1/6} ]
This slope is not universal. It applies only to the stated defect model and assumptions. Acceptor-doped regimes, hole-dominated regimes, defect association, and non-dilute site occupancy can produce different dependencies.
For a broader treatment of charge neutrality and defect-regime construction, see Defect Equilibria in Functional Oxides.
7. Chemical potential of water vapour and hydration
Water vapour participates directly in the hydration of oxygen vacancies in proton-conducting oxides:
[ ]
The corresponding equilibrium constant is
[ K_{} = ]
For an ideal water-vapour atmosphere,
[ a_{}= ]
and for a real gas, fugacity should be used instead of partial pressure.
Increasing water-vapour activity generally shifts the hydration equilibrium toward protonic defects, provided the temperature, composition, and other relevant chemical potentials are held appropriately. However, the total proton concentration is also constrained by the available oxygen vacancies, dopant concentration, site occupancy, and competing defect reactions.
The distinction between equilibrium proton concentration and proton mobility is essential. Thermodynamics determines the equilibrium state; it does not by itself determine the migration rate or conductivity.
See Proton-Conducting Oxides for the broader defect-chemistry and transport framework.
8. Chemical potential and phase equilibrium
At equilibrium between two phases () and (), the chemical potential of each transferable component must be equal:
[ _i=_i ]
If the chemical potentials differ, transfer of component (i) can lower the total Gibbs free energy. The direction of transfer is from higher chemical potential toward lower chemical potential, subject to the system constraints.
For a reaction,
[ _i _i A_i=0 ]
the reaction Gibbs energy is
[ _r G=_i_i_i ]
Using the activity expression for chemical potential,
[ _r G = _r G^+RTQ ]
where
[ Q=_i a_i^{_i} ]
At equilibrium,
[ _r G=0,Q=K ]
and therefore
[ _r G^=-RTK ]
This relation links standard Gibbs energy, equilibrium constants, and the activities of the reacting species. It is the basis for quantitative analysis of oxidation, reduction, hydration, and phase transformations.
9. Practical interpretation for functional oxide experiments
When interpreting an oxide experiment, it is useful to identify the chemical potentials imposed by the atmosphere before assigning a defect mechanism.
| Experimental variable | Thermodynamic role | Typical consequence |
|---|---|---|
| Temperature, (T) | Changes standard chemical potentials and equilibrium constants | Changes equilibrium defect populations and phase stability |
| Oxygen partial pressure, (p_{}) | Sets oxygen chemical potential | Controls oxidation state and oxygen nonstoichiometry |
| Water-vapour pressure, (p_{}) | Sets water chemical potential | Influences hydration and proton incorporation |
| Hydrogen partial pressure, (p_{}) | Sets hydrogen chemical potential | Influences reduction and gas–solid equilibria |
| Total pressure | Affects gas fugacities and pressure contributions | Important especially at elevated pressures |
A defensible thermodynamic interpretation should state:
- The temperature and gas composition at the sample, not merely the nominal inlet settings.
- Whether ideal-gas partial pressures or fugacities are used.
- The standard-state convention for every activity.
- Whether solid defects are treated as dilute species or through a non-ideal solution model.
- The charge-neutrality condition and the assumed dominant defects.
- Whether the experiment reached equilibrium or may be limited by gas exchange and solid-state kinetics.
In humid atmospheres, the actual oxygen partial pressure and water-vapour pressure should be distinguished. Gas mixtures, leaks, gas humidification, and temperature gradients can make the atmosphere at the sample differ from the nominal supply composition.
10. Common mistakes
- Using pressure inside a logarithm without normalization. Use (p_i/p^), or fugacity divided by the standard pressure.
- Confusing oxygen molecules with oxygen atoms. Expressions written per mole of oxygen atoms contain the factor (1/2) relative to molecular ().
- Treating activity as concentration in every regime. This is an approximation that requires a stated model.
- Assuming every pure-looking oxide has unit activity. Non-stoichiometric solids and solid solutions can have composition-dependent chemical potentials.
- Mixing standard states. Equilibrium constants and standard chemical potentials must use mutually consistent conventions.
- Equating thermodynamic favourability with fast kinetics. A favourable equilibrium does not guarantee that the material reaches it on the experimental timescale.
11. Summary
Chemical potential is the partial molar Gibbs free energy and provides the thermodynamic driving force for component transfer and chemical reactions. Activity expresses the effective concentration of a species relative to a defined standard state, allowing non-ideal behaviour to be incorporated consistently.
For functional oxides, oxygen chemical potential is directly connected to oxygen partial pressure, while water chemical potential is important for hydration and proton incorporation. These quantities enter mass-action expressions and determine equilibrium defect populations together with charge neutrality and site constraints.
The practical lesson is to formulate defect and phase-equilibrium models in terms of chemical potentials and activities first, and only then introduce concentration approximations. This keeps the thermodynamic assumptions explicit and makes the interpretation of oxygen nonstoichiometry, proton uptake, and redox experiments more reliable.