Thermodynamic Equilibrium and the Mass Action Law in Solids: Reaction Gibbs Energy, Equilibrium Constants, and Oxide Reactions
Thermodynamic Equilibrium and the Mass Action Law in Solids: Reaction Gibbs Energy, Equilibrium Constants, and Oxide Reactions
The mass-action law is used throughout solid-state chemistry to describe oxidation, reduction, hydration, defect formation, and phase transformations. In materials research, it is often introduced as a formula to memorize. A more reliable approach is to derive it from the Gibbs free-energy change of a reaction.
This connection clarifies what an equilibrium constant means, why activities rather than raw concentrations appear in thermodynamic expressions, how to determine the direction of a reaction away from equilibrium, and how temperature changes the equilibrium state.
The focus here is the thermodynamic structure of chemical equilibrium. Detailed Kröger–Vink bookkeeping, charge-neutrality models, and Brouwer-diagram construction are treated separately in Defect Equilibria in Functional Oxides.
1. Defining a chemical reaction
Consider a general reaction involving species (A_i):
[ _i _i A_i=0 ]
The stoichiometric coefficient (_i) is positive for products and negative for reactants. For example,
[ ]
can be written as
[ -A-B+C=0 ]
The reaction extent, (), describes how far the reaction has proceeded. The change in the amount of each species is
[ dn_i=_i,d ]
At constant temperature and pressure, the differential Gibbs free energy is
[ dG=_i_i,dn_i ]
Substituting the stoichiometric relation gives
[ dG=(_i_i_i)d ]
The reaction Gibbs energy is therefore
[ ]
It is the change in Gibbs free energy per unit extent of reaction, under the stated thermodynamic conditions.
2. Reaction direction and equilibrium
The sign of (_rG) determines the thermodynamic direction in which the reaction can proceed spontaneously at constant temperature and pressure:
- (_rG<0): the forward reaction lowers Gibbs free energy.
- (_rG>0): the reverse reaction lowers Gibbs free energy.
- (_rG=0): the reaction is at equilibrium.
These statements describe thermodynamic driving force, not reaction rate. A reaction with a strongly negative (_rG) can still be slow if it has a large activation barrier or is limited by diffusion, nucleation, or interfacial processes.
For a system containing several independent reactions, equilibrium requires the reaction Gibbs energy of each independent reaction to be zero, subject to the conservation constraints.
3. Deriving the reaction quotient
The chemical potential of species (i) can be expressed in terms of its activity:
[ _i=_i^+RTa_i ]
Substitute this into the reaction Gibbs energy:
[ _rG =_i_i_i^ +RT_i_ia_i ]
Define the standard reaction Gibbs energy:
[ _rG=_i_i_i ]
Using the logarithm rule,
[ _i_ia_i =(_i a_i^{_i}) ]
The reaction quotient is
[ Q=_i a_i^{_i} ]
and therefore
[ ]
This is the central equation connecting reaction thermodynamics to composition. It applies away from equilibrium as well as at equilibrium.
The activity of each species must be defined using a consistent standard-state convention. For gases, activities are commonly expressed using fugacity divided by standard pressure; for ideal gases this is approximated by partial pressure divided by standard pressure. For solutions and solid solutions, activity coefficients may be required.
For a detailed discussion of activity and standard states, see Chemical Potential and Activity in Solid Materials.
4. The equilibrium constant
At equilibrium,
[ _rG=0 ]
and the reaction quotient takes its equilibrium value, denoted (K):
[ Q=K ]
Thus,
[ 0=_rG^+RTK ]
or
[ ]
Equivalently,
[ ]
The equilibrium constant is dimensionless when the activities are dimensionless and consistently defined.
The magnitude of (K) indicates the relative thermodynamic preference under the specified standard-state convention:
- (K): products are favoured at equilibrium.
- (K): reactants are favoured at equilibrium.
- (K): neither side is overwhelmingly favoured.
A large equilibrium constant does not mean the reaction is fast. Nor does it mean the products are pure: the equilibrium composition depends on the complete reaction stoichiometry, mass balance, and activities.
5. A quantitative example: reaction direction away from equilibrium
Consider a hypothetical reaction
[ ]
with dimensionless activities (a_A), (a_B), and (a_C). Its reaction quotient is
[ Q= ]
Suppose the equilibrium constant at a selected temperature is (K=100), while the current activities give (Q=10).
The reaction Gibbs energy is
[ _rG=RT() ]
At (T=800 K),
[ _rG =(8.314)(800)(0.1) ]
[ _rG ]
The forward reaction is thermodynamically favourable under these conditions because (Q<K). The system tends to consume reactants and form product until the reaction quotient reaches its equilibrium value, provided no kinetic or transport constraint prevents equilibration.
The example is illustrative: the activities are dimensionless, and the stated (K) is assumed to apply at 800 K.
6. How pure solids and liquids enter the equilibrium expression
A pure solid or pure liquid in its standard state is assigned unit activity:
[ a_{}=1 ]
This means its activity does not appear explicitly in the reaction quotient while that pure phase remains present and its standard-state approximation is valid.
For example, consider the reduction reaction
[ ]
If (MO) and (M) are treated as pure solid phases, then
[ Q= {a_{MO}a_{}} ]
and with (a_M=a_{MO}=1),
[ Q= ]
For ideal gases,
[ Q= {p_{}/p^} = ]
The pressure ratio is dimensionless. If either solid is actually a solid solution or a non-stoichiometric phase with composition-dependent chemical potential, assigning unit activity may be inappropriate.
A phase absent from the system also requires care: the unit-activity convention for a present pure phase should not be used to imply that an absent phase has a finite amount.
7. Equilibrium constants for defect reactions
The same thermodynamic derivation applies to point-defect reactions in oxides. Consider oxygen loss:
[ ]
The reaction quotient is formally
[ Q_{}= ]
The thermodynamic equilibrium constant is the value of this quotient at equilibrium.
For a dilute-defect model, the regular oxygen-site activity may be approximated as constant and absorbed into an effective equilibrium constant. Defect activities may also be approximated by normalized site fractions. These approximations lead to the familiar concentration-based mass-action expressions.
The distinction is important:
- The thermodynamic expression is written in activities.
- The concentration form follows only after a model for the activities is selected.
- Charge neutrality and site balance are additional constraints; they are not supplied by the mass-action equation alone.
For the detailed defect species, charge-neutrality relations, and oxygen-pressure regimes, see Defect Equilibria in Functional Oxides.
8. Hydration equilibrium in proton-conducting oxides
A common hydration reaction in acceptor-doped perovskite oxides is
[ ]
The corresponding activity-based equilibrium constant is
[ K_{}= ]
If the regular oxygen-site activity is treated as approximately constant, it can be incorporated into an effective equilibrium constant. In a dilute-defect approximation, normalized defect concentrations may then be used.
The reaction shows why proton concentration depends on water activity and vacancy availability. However, this equilibrium relation alone does not determine the complete defect population. Charge neutrality, dopant concentration, site occupancy, and competing reactions must also be considered.
The thermodynamic equilibrium composition should also be distinguished from proton conductivity: mobility and transport kinetics are separate material properties.
9. Combining multiple equilibria
A useful feature of thermodynamic equilibrium constants is that they can be combined algebraically.
If reaction 1 has equilibrium constant (K_1), and reaction 2 has equilibrium constant (K_2), then:
- Reversing a reaction changes its equilibrium constant to (1/K).
- Multiplying every stoichiometric coefficient by (n) changes the equilibrium constant to (K^n).
- Adding independent reactions multiplies their equilibrium constants.
For example, if
[ ,K_1 ]
and
[ ,K_2 ]
then the overall reaction
[ ]
has
[ K_{}=K_1K_2 ]
This follows because standard reaction Gibbs energies add:
[ rG{}^ = _rG_1+_rG_2 ]
Combining reactions is useful when constructing defect models involving oxygen exchange, hydration, electronic compensation, and phase formation. It also provides a consistency check: a closed cycle of reactions must produce a consistent net standard Gibbs energy.
10. Temperature dependence of the equilibrium constant
The equilibrium constant depends on temperature because the standard reaction Gibbs energy does:
[ _rG=_rH-T_rS^ ]
Combining this with
[ _rG^=-RTK ]
gives
[ K= - + ]
If the reaction enthalpy is approximately constant over the temperature interval, differentiation yields the van ’t Hoff equation:
[ ]
For an endothermic forward reaction, (_rH^>0), so (K) increases with temperature under this approximation. For an exothermic forward reaction, (K) decreases with temperature.
The assumption of constant reaction enthalpy is not always adequate over a broad temperature range. Heat-capacity differences between products and reactants can make (_rH^) temperature-dependent.
The van ’t Hoff equation is developed in more detail in the next article of this series.
11. Practical workflow for using the mass-action law
When applying a mass-action expression to an oxide or ceramic reaction:
- Write the balanced reaction. Check elemental balance, charge balance where applicable, and the correct stoichiometric coefficients.
- Choose the standard states. Identify the reference state for gases, pure solids, liquids, and solution components.
- Write activities first. Construct (Q) from the stoichiometric coefficients before introducing concentration approximations.
- Determine the reaction Gibbs energy. Use (_rG=_rG^+RTQ) to assess the reaction direction.
- Set the equilibrium condition. At equilibrium, (Q=K) and (_rG=0).
- Apply material balances. Conserve the total amount of each chemical component.
- Apply charge neutrality and site balance. These constraints are essential in ionic solids.
- Check the model assumptions. Assess ideality, dilute-defect limits, fixed site populations, and the presence of competing phases.
- Separate equilibrium from kinetics. Confirm whether the experimental time and temperature permit equilibration.
12. Common mistakes
- Using concentrations as if they were always activities. This is valid only under an appropriate ideal or dilute approximation.
- Forgetting standard-state normalization. Gas activities are dimensionless pressure or fugacity ratios.
- Using the equilibrium constant to predict reaction rate. (K) describes equilibrium, not kinetics.
- **Assuming (_rG^) at equilibrium.** The actual reaction Gibbs energy is zero at equilibrium; the standard reaction Gibbs energy generally is not.
- Omitting charge neutrality in defect chemistry. Mass action alone does not determine every defect concentration.
- Including pure solid activities unnecessarily. A present pure solid in its standard state has unit activity, but solid solutions and non-stoichiometric phases need a suitable model.
- Changing reaction normalization without changing (K). Multiplying a reaction by (n) changes its equilibrium constant to (K^n).
13. Summary
The mass-action law follows directly from the dependence of chemical potential on activity. For a reaction with quotient (Q),
[ _rG=_rG^+RTQ ]
At equilibrium, (Q=K), giving
[ _rG^=-RTK ]
This framework identifies the thermodynamic direction of a reaction, defines the equilibrium constant, and connects standard reaction energies to measurable or modelled activities.
In functional oxides, the same relations describe oxygen exchange, reduction, hydration, and defect formation. Reliable use of the mass-action law requires consistent standard states and activity models, together with material balance, charge neutrality, and a clear distinction between equilibrium thermodynamics and reaction kinetics.
Further reading
- D. R. Gaskell and D. E. Laughlin, Introduction to the Thermodynamics of Materials, 6th ed., CRC Press.
- H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, 2nd ed., Wiley.
- P. Atkins, J. de Paula, and J. Keeler, Atkins’ Physical Chemistry, Oxford University Press.
- A. R. West, Solid State Chemistry and its Applications, 2nd ed., Wiley.