Reaction equilibria

You already know how to write KcK_c or KpK_p for a reaction and predict which way Le Chatelier's principle pushes it — that is CBSE/JEE-level equilibrium. What you have not yet seen is why KK takes the form it does, why it depends on temperature alone and never on pressure or concentration, or why heating an exothermic reaction reliably lowers KK. CY11001 answers all three from a single starting point: the chemical potential expression μi=μi∘+RTln⁡ai\mu_i=\mu_i^\circ+RT\ln a_i carried over from the previous article. Writing the equilibrium condition in terms of species activities aia_i produces the general equilibrium constant K=∏iaiνiK=\prod_i a_i^{\nu_i}, and differentiating ln⁡K\ln K with respect to TT produces the van't Hoff equation — a genuine derivative that replaces the qualitative "system opposes the change" heuristic with an equation you can actually compute.

The core insight

The thermodynamic equilibrium constant K=∏iaiνiK=\prod_i a_i^{\nu_i} follows from ΔrG∘=−RTln⁡K\Delta_rG^\circ=-RT\ln K and depends on temperature alone, and differentiating this relation gives the van't Hoff equation dln⁡K/dT=ΔrH∘/(RT2)d\ln K/dT=\Delta_rH^\circ/(RT^2), which turns Le Chatelier's temperature rule into a computable derivative.

The core idea

From chemical potential to KK. The equilibrium condition for a reaction is that the weighted sum of chemical potentials of products equals that of reactants, each weighted by its stoichiometric coefficient νi\nu_i (negative for reactants, positive for products). Substituting the general form μi=μi∘+RTln⁡ai\mu_i=\mu_i^\circ+RT\ln a_i — where aia_i is the species activity, a dimensionless "effective concentration" referenced to a standard state (ai=pi/p∘a_i=p_i/p^\circ for an ideal gas, ai≈ci/c∘a_i\approx c_i/c^\circ for a dilute solute) — and rearranging gives ΔrG∘=−RTln⁡K\Delta_rG^\circ=-RT\ln K with

K=∏iaiνiK=\prod_i a_i^{\nu_i}

where: aia_i is the activity of species ii and νi\nu_i its stoichiometric coefficient. Because every μi∘\mu_i^\circ on the right of ΔrG∘=−RTln⁡K\Delta_rG^\circ=-RT\ln K is a standard chemical potential — fixed once TT is fixed, independent of how much of each species is actually present — KK itself is a function of TT alone. This is the fact CBSE states but CY11001 proves: changing pressure, volume or concentration shifts the equilibrium composition to keep KK satisfied, but never changes KK itself.

The van't Hoff equation. Differentiate ln⁡K=−ΔrG∘/(RT)\ln K=-\Delta_rG^\circ/(RT) with respect to TT and use the Gibbs-Helmholtz relation (∂(ΔG∘/T)/∂T)P=−ΔH∘/T2\left(\partial(\Delta G^\circ/T)/\partial T\right)_P=-\Delta H^\circ/T^2:

dln⁡KdT=ΔrH∘RT2\frac{d\ln K}{dT}=\frac{\Delta_rH^\circ}{RT^2}

where: ΔrH∘\Delta_rH^\circ is the standard reaction enthalpy (assumed roughly constant over the temperature range considered) and R=8.314 J mol−1K−1R=8.314\ \text{J mol}^{-1}\text{K}^{-1}. Integrating between two temperatures gives the working form used in calculations:

ln⁡K2K1=−ΔrH∘R(1T2−1T1)\ln\frac{K_2}{K_1}=-\frac{\Delta_rH^\circ}{R}\left(\frac{1}{T_2}-\frac{1}{T_1}\right)

This is structurally identical to the integrated Clausius-Clapeyron equation from phase equilibrium, because both descend from the same d(ΔG∘/T)/dT=−ΔH∘/T2d(\Delta G^\circ/T)/dT=-\Delta H^\circ/T^2 mathematics.

Le Chatelier's rule, now a sign. For an exothermic reaction (ΔrH∘<0\Delta_rH^\circ<0), the van't Hoff equation gives dln⁡K/dT<0d\ln K/dT<0, so KK falls as TT rises — heating disfavours products. For an endothermic reaction (ΔrH∘>0\Delta_rH^\circ>0), KK rises with TT. The qualitative CBSE rule is exactly this sign, now derived rather than asserted.

Figure
A graph of the natural log of K plotted against one over T showing a straight line whose downward slope equals minus the standard reaction enthalpy divided by R for an exothermic reaction
van't Hoff plot: ln K versus 1/T 1/T (higher toward the right = lower T) ln K slope = -(delta H standard)/R K at 500 K K at 600 K, smaller Negative slope means an exothermic reaction: K drops as T rises
Derivation ADeriving the van't Hoff equation from the equilibrium constant

The goal is to show, using only the chain rule and the Gibbs-Helmholtz relation already established for ΔG∘/T\Delta G^\circ/T, that the temperature dependence of the equilibrium constant KK is governed by a single derivative dln⁡K/dT=ΔrH∘/(RT2)d\ln K/dT=\Delta_rH^\circ/(RT^2), and that this derivative — not a qualitative appeal to "the system resists change" — is the correct thermodynamic content behind Le Chatelier's temperature rule.

Step 1. Start from the relation between KK and standard Gibbs energy, ΔrG∘=−RTln⁡K\Delta_rG^\circ=-RT\ln K, and rearrange it to isolate ln⁡K\ln K: ln⁡K=−ΔrG∘RT=−1R(ΔrG∘T)\ln K=-\dfrac{\Delta_rG^\circ}{RT}=-\dfrac{1}{R}\left(\dfrac{\Delta_rG^\circ}{T}\right). Written this way, ln⁡K\ln K is a constant (−1/R-1/R) times the very quantity ΔrG∘/T\Delta_rG^\circ/T whose temperature derivative the Gibbs-Helmholtz relation already gives.

Step 2. Differentiate both sides with respect to TT at constant pressure. dln⁡KdT=−1RddT(ΔrG∘T)\dfrac{d\ln K}{dT}=-\dfrac{1}{R}\dfrac{d}{dT}\left(\dfrac{\Delta_rG^\circ}{T}\right). This step is purely mechanical differentiation of Step 1; no new physics has entered yet.

Step 3. Substitute the Gibbs-Helmholtz relation, (∂(ΔG∘/T)/∂T)P=−ΔH∘/T2\left(\partial(\Delta G^\circ/T)/\partial T\right)_P=-\Delta H^\circ/T^2, applied to the reaction quantities ΔrG∘\Delta_rG^\circ and ΔrH∘\Delta_rH^\circ: dln⁡KdT=−1R(−ΔrH∘T2)=ΔrH∘RT2\dfrac{d\ln K}{dT}=-\dfrac{1}{R}\left(-\dfrac{\Delta_rH^\circ}{T^2}\right)=\dfrac{\Delta_rH^\circ}{RT^2}. This is the van't Hoff equation.

Step 4. Integrate assuming ΔrH∘\Delta_rH^\circ is roughly constant over T1T_1 to T2T_2. ∫ln⁡K1ln⁡K2dln⁡K=ΔrH∘R∫T1T2dTT2\displaystyle\int_{\ln K_1}^{\ln K_2}d\ln K=\frac{\Delta_rH^\circ}{R}\int_{T_1}^{T_2}\frac{dT}{T^2} gives ln⁡K2K1=−ΔrH∘R(1T2−1T1)\ln\dfrac{K_2}{K_1}=-\dfrac{\Delta_rH^\circ}{R}\left(\dfrac{1}{T_2}-\dfrac{1}{T_1}\right), the working integrated form. Reading the sign directly: if ΔrH∘<0\Delta_rH^\circ<0 (exothermic) and T2>T1T_2>T_1, the bracket (1/T2−1/T1)(1/T_2-1/T_1) is negative, so the two negatives cancel and ln⁡(K2/K1)\ln(K_2/K_1) is negative — KK falls on heating, exactly the exothermic Le Chatelier rule, now obtained as an algebraic consequence rather than assumed.

Worked example

For N2(g)+3H2(g)⇌2NH3(g)\text{N}_2(g)+3\text{H}_2(g)\rightleftharpoons 2\text{NH}_3(g), Kp=1.45×10−5K_p=1.45\times10^{-5} at 500 K and ΔrH∘=−92.4 kJ mol−1\Delta_rH^\circ=-92.4\ \text{kJ mol}^{-1} (exothermic, roughly constant over the range). Estimate KpK_p at 600 K, and use the answer to find the degree of dissociation regime it implies is favoured.

Step 1. Apply the integrated van't Hoff equation with T1=500T_1=500 K, K1=1.45×10−5K_1=1.45\times10^{-5}, T2=600T_2=600 K. ln⁡K2K1=−ΔrH∘R(1T2−1T1)=−(−92400)8.314(1600−1500)\ln\dfrac{K_2}{K_1}=-\dfrac{\Delta_rH^\circ}{R}\left(\dfrac{1}{T_2}-\dfrac{1}{T_1}\right)=-\dfrac{(-92400)}{8.314}\left(\dfrac{1}{600}-\dfrac{1}{500}\right).

Step 2. Evaluate the bracket. 1600=1.6667×10−3\frac{1}{600}=1.6667\times10^{-3}, 1500=2.0000×10−3\frac{1}{500}=2.0000\times10^{-3}, so the difference is −3.333×10−4 K−1-3.333\times10^{-4}\ \text{K}^{-1}.

Step 3. Combine. ln⁡K2K1=−−924008.314×(−3.333×10−4)=11114.3×(−3.333×10−4)=−3.705\ln\dfrac{K_2}{K_1}=-\dfrac{-92400}{8.314}\times(-3.333\times10^{-4})=11114.3\times(-3.333\times10^{-4})=-3.705, so K2K1=e−3.705=0.02461\dfrac{K_2}{K_1}=e^{-3.705}=0.02461 and K2=(1.45×10−5)(0.02461)=3.57×10−7K_2=(1.45\times10^{-5})(0.02461)=3.57\times10^{-7}.

Step 4. Interpret via the dissociation link. Since KpK_p falls by roughly two orders of magnitude on heating, the equilibrium extent of the reverse reaction (ammonia decomposing back to N2\text{N}_2 and H2\text{H}_2) is strongly favoured at 600 K compared with 500 K — consistent with the general degree-of-dissociation relation Kp=α21−α2⋅Pp∘K_p=\dfrac{\alpha^2}{1-\alpha^2}\cdot\dfrac{P}{p^\circ} derived for a dissociating gas-phase equilibrium: a smaller KpK_p at fixed PP forces a smaller equilibrium α\alpha for the forward process, i.e. less product. Answer: Kp(600 K)≈3.57×10−7K_p(600\ \text{K})\approx3.57\times10^{-7}, roughly 1/401/40 of its value at 500 K, obtained purely from ΔrH∘\Delta_rH^\circ — the quantitative version of "heating an exothermic reaction disfavours products."

Competitive stretch

Suppose instead you were given KpK_p at both 500 K and 600 K but not ΔrH∘\Delta_rH^\circ — could you recover it? Yes: the same integrated van't Hoff equation is linear in the unknown ΔrH∘\Delta_rH^\circ once K1K_1, K2K_2, T1T_1 and T2T_2 are known, so it can be solved by rearrangement rather than integration. Using the two values already found, K1=1.45×10−5K_1=1.45\times10^{-5} at T1=500T_1=500 K and K2=3.57×10−7K_2=3.57\times10^{-7} at T2=600T_2=600 K, rearrange ln⁡(K2/K1)=−ΔrH∘R(1T2−1T1)\ln(K_2/K_1)=-\dfrac{\Delta_rH^\circ}{R}\left(\dfrac{1}{T_2}-\dfrac{1}{T_1}\right) for ΔrH∘\Delta_rH^\circ: ΔrH∘=−Rln⁡(K2/K1)/(1T2−1T1)=−(8.314)(−3.705)/(−3.333×10−4)=−92,400 J mol−1\Delta_rH^\circ=-R\ln(K_2/K_1)\Big/\left(\dfrac{1}{T_2}-\dfrac{1}{T_1}\right)=-(8.314)(-3.705)/(-3.333\times10^{-4})=-92{,}400\ \text{J mol}^{-1}, recovering the same −92.4 kJ mol−1-92.4\ \text{kJ mol}^{-1} used to construct the problem — the check that the van't Hoff equation is a genuinely invertible relation between ΔrH∘\Delta_rH^\circ and the pair of equilibrium constants, exactly the way two points on a straight line fix its slope.

Common traps

  • Believing KK itself changes when reactant is added or pressure is changed — KK is a function of TT only; what changes is QQ relative to KK, or (for pressure and volume changes) the equilibrium composition/extent that still satisfies the fixed KpK_p, never KpK_p itself.
  • Applying the integrated van't Hoff equation across a temperature range where ΔrH∘\Delta_rH^\circ is not roughly constant, such as one spanning a phase change of a reactant — the constant-ΔH\Delta H assumption can fail over wide ranges, exactly as in Clausius-Clapeyron.
  • Reversing the sign rule under exam pressure: an increase in TT increases KK for endothermic reactions (ΔrH∘>0\Delta_rH^\circ>0) and decreases KK for exothermic reactions — students commonly remember it backwards.
  • Forgetting that KK built from activities is dimensionless by construction (each activity already divides by p∘p^\circ or c∘c^\circ); plugging a KpK_p quoted with explicit pressure units directly into ΔrG∘=−RTln⁡K\Delta_rG^\circ=-RT\ln K produces a dimensionally meaningless logarithm.

Ask the tutor: If you want Le Chatelier's pressure and concentration rules derived thermodynamically the same way the temperature rule was derived here, ask the AI tutor to work through the QQ-versus-KK argument for those cases too.

If K only depends on temperature, why does adding more reactant shift the equilibrium at all?

Adding reactant does not change KK, but it does change the reaction quotient QQ, which is computed the same way as KK but from the current (non-equilibrium) activities. Immediately after adding reactant, Q<KQ<K, so the reaction proceeds forward until the new composition again satisfies Q=KQ=K. The system moves to a new equilibrium composition; the value of KK itself never moves, because KK is fixed by the standard chemical potentials at that temperature only.

Why does the van't Hoff equation look so similar to Clausius-Clapeyron?

Both are integrated forms of the same underlying mathematical statement, d(ΔG∘/T)/dT=−ΔH∘/T2d(\Delta G^\circ/T)/dT=-\Delta H^\circ/T^2, applied to two different equilibria. For phase equilibrium, ΔH∘\Delta H^\circ is the enthalpy of the phase change and the "equilibrium constant" is effectively the vapour pressure; for reaction equilibrium, ΔH∘\Delta H^\circ is the reaction enthalpy and the constant is KK. The shared Gibbs-Helmholtz origin is why both integrate to the same 1/T2−1/T11/T_2-1/T_1 functional form.

How do I know when the constant-delta-H assumption in the integrated van't Hoff equation is safe to use?

It is safe when the temperature range is narrow enough, and free of any phase change of a reactant or product, that ΔrH∘\Delta_rH^\circ does not vary appreciably — typically a few tens of kelvin for problems at this level. Over wider ranges, ΔrH∘\Delta_rH^\circ itself depends on TT through the heat capacities of reactants and products, and the simple integrated form used in the worked example above would need a correction term. Unless a problem explicitly signals a wide range or a phase change, treat ΔrH∘\Delta_rH^\circ as constant, exactly as with Clausius-Clapeyron.

Before this

Where this leads