The Gibbs-Helmholtz Equation

What the equation does

Free energy is the quantity that decides whether a reaction runs. But the number you care about — ΔG — is not a constant: it changes as you heat the flask. The Gibbs-Helmholtz equation is the rule that governs how ΔG changes with temperature, and its remarkable feature is that the whole change is dictated by a single quantity you can measure independently: the enthalpy ΔH.

Written for a chemical process, it is:

    ⎡ ∂ ⎛ ΔG ⎞ ⎤          ΔH
    ⎢ ── ⎜ ── ⎟ ⎥    =  − ────
    ⎣ ∂T ⎝  T ⎠ ⎦_P        T²

Read it in words: the rate at which the "free energy per kelvin" (ΔG/T) changes as you warm the system is minus the enthalpy, damped by T². The genius is the grouping. You do not track ΔG directly — you track ΔG/T. That one division is what makes the messy entropy term disappear and leaves a slope you can trust.

The derivation, step by step

Everything starts from the definition of Gibbs free energy and one master relation from the fundamental equation of thermodynamics.

1. Definition. For any system at temperature T, G = H − TS.
2. The temperature slope of G. From the fundamental equation dG = −S dT + V dP, holding pressure constant gives (∂G/∂T)_P = −S. This is the piece that carries the inconvenient, temperature-dependent entropy.
3. Differentiate the quotient G/T. Apply the quotient rule to G/T with respect to T at constant P:

     ∂(G/T)/∂T = (1/T)(∂G/∂T)_P − G/T²
               = (1/T)(−S) − G/T²
               = −S/T − (H − TS)/T²

4. Watch the entropy cancel. Expand the last term: −(H − TS)/T² = −H/T² + S/T. The +S/T exactly annihilates the −S/T sitting in front:

     ∂(G/T)/∂T = −S/T − H/T² + S/T = −H/T²

5. Take the difference over a reaction. Subtract the reactant expression from the product expression (Δ means "products minus reactants") and you land on the working form: [∂(ΔG/T)/∂T]_P = −ΔH/T².

That cancellation in step 4 is the entire story. Entropy is temperature-dependent and hard to hold constant; enthalpy is comparatively stable. By choosing to differentiate G/T instead of G, Gibbs and Helmholtz traded the slippery variable for the sturdy one.

A cleaner form: slope against 1/T

The T² in the denominator looks awkward until you change variables. Let x = 1/T, so dx = −(1/T²) dT. Substituting turns the equation into something a chemist can plot with a ruler:

    ⎡ ∂ ⎛ ΔG ⎞ ⎤
    ⎢ ── ⎜ ── ⎟ ⎥    =  ΔH
    ⎣ ∂(1/T) ⎝ T ⎠ ⎦_P

In other words: plot ΔG/T on the y-axis against 1/T on the x-axis, and the slope of that line is exactly ΔH. No calorimeter, no bomb, no ice bath. If ΔH is roughly constant over your temperature range, the plot is a straight line whose steepness hands you the enthalpy directly.

The integrated form

Chemists rarely want the derivative; they want to move a known free energy from one temperature to another. Integrating between T₁ and T₂, treating ΔH as constant, gives the practical two-point form:

    ΔG₂     ΔG₁          ⎛  1     1  ⎞
    ───  −  ───  =  ΔH  ⎜ ──  − ── ⎟
     T₂      T₁          ⎝  T₂    T₁ ⎠

Given ΔG at 298 K and the reaction enthalpy, you can predict ΔG at 350 K, 500 K, or any temperature, as long as ΔH does not drift too much. This is the workhorse expression you actually paste into a spreadsheet.

Worked example: ammonia synthesis at two temperatures

Take the Haber-Bosch reaction, N₂(g) + 3H₂(g) → 2NH₃(g). At 298 K the standard values are ΔH° = −92.2 kJ/mol and ΔG° = −33.0 kJ/mol. What is ΔG° at 500 K, the low end of an industrial reactor?

  ΔG₂/T₂ = ΔG₁/T₁ + ΔH°(1/T₂ − 1/T₁)

  ΔG₁/T₁ = −33 000 J·mol⁻¹ / 298 K = −110.7 J·mol⁻¹·K⁻¹

  1/T₂ − 1/T₁ = 1/500 − 1/298
              = 0.002000 − 0.003356
              = −0.001356 K⁻¹

  ΔH°(1/T₂ − 1/T₁) = (−92 200)(−0.001356)
                    = +125.0 J·mol⁻¹·K⁻¹

  ΔG₂/T₂ = −110.7 + 125.0 = +14.3 J·mol⁻¹·K⁻¹

  ΔG₂ = (+14.3)(500) = +7 200 J·mol⁻¹ ≈ +7.2 kJ/mol

The free energy has flipped sign — from −33 kJ/mol (spontaneous) at 298 K to about +7 kJ/mol (non-spontaneous) at 500 K. That is precisely why ammonia synthesis is fought against by high temperature: the exothermic, entropy-losing reaction becomes thermodynamically unfavourable as you heat it. Industry runs hot anyway (~450 °C) only because the rate is hopeless when cold, then compensates by cranking the pressure to 150–300 atm to drag equilibrium back toward product. The Gibbs-Helmholtz equation is the quantitative statement of that trade-off.

Gibbs-Helmholtz vs its relatives

	Gibbs-Helmholtz	van 't Hoff	Clausius-Clapeyron
Variable tracked	ΔG/T	ln K	ln P
Slope gives	ΔH (vs 1/T)	−ΔH°/R (vs 1/T)	−ΔHvap/R (vs 1/T)
Core statement	∂(ΔG/T)/∂T = −ΔH/T²	d(ln K)/dT = ΔH°/RT²	d(ln P)/dT = ΔHvap/RT²
Input you need	Free energies	Equilibrium constants	Vapor pressures
Applies to	Any process at constant P	Chemical equilibria	Phase transitions
Relationship	Parent form	Gibbs-Helmholtz + ΔG° = −RT ln K	Special case for a phase change

All three are the same T² structure. Once you see that Gibbs-Helmholtz is the general parent — and that substituting the relevant free-energy definition turns it into van 't Hoff (for K) or Clausius-Clapeyron (for vapor pressure) — the family stops being three separate formulas to memorize and becomes one.

The bridge to van 't Hoff

The most important consequence of Gibbs-Helmholtz is the van 't Hoff equation, and the derivation is a two-line substitution. Start from the standard-state form and insert ΔG° = −RT ln K:

  ΔG°/T = −R ln K

  ∂(ΔG°/T)/∂T = −R · d(ln K)/dT      (left side)
  ∂(ΔG°/T)/∂T = −ΔH°/T²              (Gibbs-Helmholtz, right side)

  ⇒  −R · d(ln K)/dT = −ΔH°/T²

  ⇒  d(ln K)/dT = ΔH°/(R T²)          (van 't Hoff)

This is why a van 't Hoff plot — ln K against 1/T — has a slope of −ΔH°/R. It is Gibbs-Helmholtz dressed in equilibrium-constant clothing. Any time you extract a reaction enthalpy from how an equilibrium constant, a solubility, or a binding affinity changes with temperature, you are using the Gibbs-Helmholtz equation whether you name it or not.

Where it earns its keep

• Electrochemistry — cell potentials. Since ΔG = −nFE, the Gibbs-Helmholtz equation converts the measured temperature coefficient of a cell (∂E/∂T) into the reaction enthalpy and entropy. This is how battery and fuel-cell thermodynamics are separated into heat and work contributions without a calorimeter.
• Binding thermodynamics. In biochemistry, plotting the association constant K of a ligand across temperature gives ΔH of binding through the van 't Hoff form — a fast alternative to isothermal titration calorimetry, and a standard cross-check on it.
• Solubility and dissolution. The temperature dependence of a salt's solubility product K_sp yields the enthalpy of dissolution directly, telling you whether a compound dissolves endothermically (dissolves better hot) or exothermically.
• Industrial equilibrium design. Predicting how a synthesis equilibrium (ammonia, methanol, water-gas shift) shifts across a reactor's temperature profile is a Gibbs-Helmholtz calculation at its core.
• Reference-data extrapolation. Thermochemical tables list ΔG°_f at 298 K; Gibbs-Helmholtz (with Kirchhoff's ΔCp correction) propagates those values to the process temperature you actually run at.

Limitations and where it breaks

• ΔH is not truly constant. Kirchhoff's law says (∂ΔH/∂T)_P = ΔCp. Over 100 K or more, ΔH drifts, and the simple integrated two-point form accumulates error. For high-accuracy work, carry ΔCp(T) through the integral rather than treating ΔH as a fixed number.
• Constant pressure only. The equation is derived at constant P (that is why the −S slope comes from dG = −S dT + V dP with dP = 0). The Helmholtz free energy A = U − TS has its own analogous Gibbs-Helmholtz relation, ∂(A/T)/∂T = −U/T², for constant-volume processes.
• Phase changes break the line. If a reactant or product melts, boils, or undergoes a transition within your temperature window, ΔH jumps discontinuously and the plot of ΔG/T versus 1/T kinks. You must split the analysis at each transition.
• It gives the trend, not the sign. Gibbs-Helmholtz tells you which way ΔG moves as T changes — it does not by itself tell you whether ΔG is negative at any given temperature. You still need an absolute value of ΔG at one reference point to anchor the prediction.

Who Gibbs and Helmholtz were

The two names belong to two of the nineteenth century's deepest thinkers, working an ocean apart.

Josiah Willard Gibbs (1839–1903) spent his entire career at Yale, publishing his monumental memoir On the Equilibrium of Heterogeneous Substances in the obscure Transactions of the Connecticut Academy between 1876 and 1878. In it he built the free-energy formalism — chemical potential, the phase rule, and the function we now call G — essentially inventing chemical thermodynamics as a subject. His work was so mathematically dense that it took Europe years to notice; Maxwell was one of the few who read it immediately.

Hermann von Helmholtz (1821–1894) came to thermodynamics from physiology and physics. In his 1882 paper Die Thermodynamik chemischer Vorgänge ("The Thermodynamics of Chemical Processes") he independently drew the sharp distinction between the total energy of a reaction and the portion available to do work — coining the term freie Energie (free energy) — and wrote down the temperature-dependence relation that now bears both men's names. That the two arrived at the same relation within a few years, from an American mathematical physicist and a German physiologist-physicist, is why the equation is doubly named.

The intuition to remember

Strip away the calculus and the equation says one memorable thing: an endothermic process is helped by heat, an exothermic one is hurt by it, and the enthalpy is the exact conversion factor. ΔG/T drops with temperature when ΔH is positive (endothermic), pushing ΔG negative — heating drives the reaction. ΔG/T climbs when ΔH is negative (exothermic), pushing ΔG positive — heating stalls the reaction. Le Chatelier's qualitative rule about temperature is just the sign of ΔH read off the Gibbs-Helmholtz slope. The equation is the number behind the intuition.