The Gibbs–Helmholtz Equation

A chemist measures ΔG∘\Delta G^\circ for a reaction at 298 K, but the reactor runs at 350 K. Rerunning the whole free-energy calculation seems necessary. It is not. The Gibbs–Helmholtz equation lets you predict how ΔG\Delta G shifts with temperature using only ΔH\Delta H, a much easier quantity to measure or estimate calorimetrically. It converts a two-variable problem, tracking both energy and disorder, into a one-variable problem governed by enthalpy alone. This article derives the equation from the fundamental properties of GG. It shows why the combination G/TG/T, not GG itself, obeys the clean derivative law. It closes with two worked numerical examples and a preview of how this equation becomes the van't Hoff equation for equilibrium constants.

Key takeaways
  • The Gibbs–Helmholtz equation states (∂(G/T)/∂T)p=−H/T2\left(\partial(G/T)/\partial T\right)_p = -H/T^2.
  • It follows directly from G=H−TSG = H - TS and the relation (∂G/∂T)p=−S(\partial G/\partial T)_p = -S.
  • The trick is dividing by TT first: this eliminates SS from the derivative entirely.
  • Applied to a reaction, it gives (∂(ΔG/T)/∂T)p=−ΔH/T2\left(\partial(\Delta G/T)/\partial T\right)_p = -\Delta H/T^2.
  • Integrating over a temperature range, assuming ΔH\Delta H is constant, predicts ΔG\Delta G at any TT from data at one reference temperature.
  • Key idea: Because SS cancels algebraically, the temperature dependence of ΔG/T\Delta G/T is controlled purely by ΔH\Delta H, making enthalpy data alone sufficient to extrapolate free energies across temperature.
The core insight

Dividing the Gibbs energy by temperature before differentiating causes the entropy term to vanish algebraically, leaving a derivative that depends only on enthalpy — this is the entire content and utility of the Gibbs–Helmholtz equation.

The physical picture

Start from the familiar temperature dependence of Gibbs energy at constant pressure, established from the fundamental equation dG=−S dT+V dpdG = -S\,dT + V\,dp:

(∂G∂T)p=−S\left(\frac{\partial G}{\partial T}\right)_p = -S

This is useful, but it leaves SS sitting in the answer. Entropy is often the harder quantity to pin down experimentally compared with enthalpy, which comes straight from calorimetry. The Gibbs–Helmholtz equation is a rewriting of exactly this same fact, but packaged so that only HH appears, never SS. The key move is to look not at GG itself, but at the ratio G/TG/T. This ratio appears constantly in equilibrium and electrochemical work because ΔG∘/T\Delta G^\circ/T is directly proportional to −Rln⁡K-R\ln K (foreshadowing the van't Hoff equation). When you differentiate G/TG/T with respect to TT at constant pp, the quotient rule produces two terms: one containing (∂G/∂T)p=−S(\partial G/\partial T)_p = -S, and another containing G/T2G/T^2. Because G=H−TSG = H - TS, these two terms combine so that every occurrence of SS cancels exactly, leaving a clean result in HH and TT alone. This cancellation is not a coincidence — it is guaranteed by the definition of GG, and it is the entire reason chemists prefer G/TG/T over GG when studying temperature effects. Physically, this tells you something important: a reaction's free energy changes with temperature for two competing reasons — the entropy term −S dT-S\,dT pulls in one direction, while the changing weight of the TSTS term inside GG pulls in another. The Gibbs–Helmholtz form shows that after these two effects are combined honestly, the net temperature sensitivity of G/TG/T is set only by how much heat the process absorbs or releases, i.e., by HH.

Figure
Slope of G/T versus 1/T reveals ΔH
1/T (K⁻¹) ΔG/T (J/mol·K) slope = ΔH 298 K point 350 K point d(ΔG/T)/d(1/T) = ΔH
Figure
Enthalpy and entropy contributions to G versus T
T (K) G (kJ/mol) G = H − TS, slope −S H (roughly constant) H ≈ −92 kJ/mol (NH₃ synthesis, example)

The core relations

(∂(G/T)∂T)p=−HT2\left(\frac{\partial (G/T)}{\partial T}\right)_p = -\frac{H}{T^2}

where: GG is Gibbs energy, HH is enthalpy, TT is absolute temperature, derivative taken at constant pressure.

(∂(ΔG/T)∂T)p=−ΔHT2\left(\frac{\partial (\Delta G/T)}{\partial T}\right)_p = -\frac{\Delta H}{T^2}

where: applies the same relation to reaction quantities, ΔG\Delta G and ΔH\Delta H for a process at fixed composition/pressure.

ΔG(T2)T2−ΔG(T1)T1=ΔH(1T2−1T1)\frac{\Delta G(T_2)}{T_2} - \frac{\Delta G(T_1)}{T_1} = \Delta H \left(\frac{1}{T_2} - \frac{1}{T_1}\right)

where: integrated form assuming ΔH\Delta H is independent of temperature over T1T_1 to T2T_2.

Derivation AFrom G = H − TS to the Gibbs–Helmholtz equation

Goal. Show that (∂(G/T)/∂T)p=−H/T2\left(\partial(G/T)/\partial T\right)_p = -H/T^2, starting only from the definition of GG and its known first derivative.

Step 1. State the known derivative. From the fundamental equation for Gibbs energy, dG=−S dT+V dpdG = -S\,dT + V\,dp, holding pressure constant gives:

(∂G∂T)p=−S\left(\frac{\partial G}{\partial T}\right)_p = -S

This was established when studying the temperature dependence of GG itself; it is the starting point here, not a new assumption.

Step 2. Differentiate the ratio G/T using the quotient rule. Treat G/TG/T as a function of TT at fixed pp:

(∂(G/T)∂T)p=1T(∂G∂T)p−GT2\left(\frac{\partial (G/T)}{\partial T}\right)_p = \frac{1}{T}\left(\frac{\partial G}{\partial T}\right)_p - \frac{G}{T^2}

This is just calculus — the derivative of a quotient f/gf/g where g=Tg = T.

Step 3. Substitute the known derivative from Step 1. Replace (∂G/∂T)p(\partial G/\partial T)_p with −S-S:

(∂(G/T)∂T)p=−ST−GT2\left(\frac{\partial (G/T)}{\partial T}\right)_p = \frac{-S}{T} - \frac{G}{T^2}

Step 4. Substitute the definition G = H − TS. Write G=H−TSG = H - TS, so G/T2=H/T2−S/TG/T^2 = H/T^2 - S/T:

(∂(G/T)∂T)p=−ST−(HT2−ST)\left(\frac{\partial (G/T)}{\partial T}\right)_p = -\frac{S}{T} - \left(\frac{H}{T^2} - \frac{S}{T}\right)

Step 5. Cancel the entropy terms. The −S/T-S/T from Step 3 and the +S/T+S/T generated in Step 4 cancel exactly:

(∂(G/T)∂T)p=−ST+ST−HT2=−HT2\left(\frac{\partial (G/T)}{\partial T}\right)_p = -\frac{S}{T} + \frac{S}{T} - \frac{H}{T^2} = -\frac{H}{T^2}

This cancellation is exact and general — it holds for any system obeying G=H−TSG = H - TS, i.e., always.

Step 6. Apply to a process (reaction). Since ΔG=Gproducts−Greactants\Delta G = G_{products} - G_{reactants} and each side individually obeys the relation above, subtracting gives the same law for the difference:

(∂(ΔG/T)∂T)p=−ΔHT2\boxed{\left(\frac{\partial (\Delta G/T)}{\partial T}\right)_p = -\frac{\Delta H}{T^2}}

Step 7. Integrate for practical use. If ΔH\Delta H is treated as constant over the temperature interval of interest (a good approximation when ΔCp≈0\Delta C_p \approx 0 or the range is modest), integrate both sides from T1T_1 to T2T_2:

∫T1T2d(ΔGT)=−ΔH∫T1T2dTT2\int_{T_1}^{T_2} d\left(\frac{\Delta G}{T}\right) = -\Delta H \int_{T_1}^{T_2} \frac{dT}{T^2}
ΔG(T2)T2−ΔG(T1)T1=ΔH(1T2−1T1)\frac{\Delta G(T_2)}{T_2} - \frac{\Delta G(T_1)}{T_1} = \Delta H\left(\frac{1}{T_2} - \frac{1}{T_1}\right)

This boxed integrated form is the workhorse equation used to extrapolate free energy data across temperature.

Worked examples

Example 1 — Foundational: extrapolating ΔG for ammonia synthesis. For N2(g)+3H2(g)→2NH3(g)\mathrm{N_2(g) + 3H_2(g) \rightarrow 2NH_3(g)} at 298 K, standard tables give:

ΔH∘=−92.4 kJ mol−1,ΔG298∘=−32.9 kJ mol−1\Delta H^\circ = -92.4\ \text{kJ mol}^{-1}, \qquad \Delta G^\circ_{298} = -32.9\ \text{kJ mol}^{-1}

Find ΔG∘\Delta G^\circ at 500 K, assuming ΔH∘\Delta H^\circ is constant. Step 1. Compute ΔG298∘/T1\Delta G^\circ_{298}/T_1:

ΔG298∘298 K=−32,900 J mol−1298 K=−110.4 J mol−1K−1\frac{\Delta G^\circ_{298}}{298\ \text{K}} = \frac{-32{,}900\ \text{J mol}^{-1}}{298\ \text{K}} = -110.4\ \text{J mol}^{-1}\text{K}^{-1}

Step 2. Apply the integrated Gibbs–Helmholtz relation:

ΔG500∘500=ΔG298∘298+ΔH∘(1500−1298)\frac{\Delta G^\circ_{500}}{500} = \frac{\Delta G^\circ_{298}}{298} + \Delta H^\circ\left(\frac{1}{500} - \frac{1}{298}\right)

Step 3. Compute the bracket:

1500−1298=0.002000−0.003356=−0.001356 K−1\frac{1}{500} - \frac{1}{298} = 0.002000 - 0.003356 = -0.001356\ \text{K}^{-1}

Step 4. Multiply by ΔH∘\Delta H^\circ:

(−92,400 J mol−1)(−0.001356 K−1)=125.3 J mol−1K−1(-92{,}400\ \text{J mol}^{-1})(-0.001356\ \text{K}^{-1}) = 125.3\ \text{J mol}^{-1}\text{K}^{-1}

Step 5. Sum and solve for ΔG500∘\Delta G^\circ_{500}:

ΔG500∘500=−110.4+125.3=14.9 J mol−1K−1\frac{\Delta G^\circ_{500}}{500} = -110.4 + 125.3 = 14.9\ \text{J mol}^{-1}\text{K}^{-1}
ΔG500∘=500×14.9=7450 J mol−1≈+7.5 kJ mol−1\Delta G^\circ_{500} = 500 \times 14.9 = 7450\ \text{J mol}^{-1} \approx +7.5\ \text{kJ mol}^{-1}

Sanity check. ΔG∘\Delta G^\circ flips sign from negative to positive as temperature rises — consistent with Le Chatelier: ammonia synthesis is exothermic with negative ΔS∘\Delta S^\circ (4 mol gas → 2 mol gas), so raising TT should disfavor products. The sign change near 450–500 K matches known industrial behavior (high pressure is needed at operating temperatures near 700 K to compensate).

Example 2 — Applied: comparing predictions using ΔH vs. naive constant-ΔG assumption. A biochemist studies a protein-folding equilibrium with ΔH∘=−180 kJ mol−1\Delta H^\circ = -180\ \text{kJ mol}^{-1} and ΔG310∘=−8.0 kJ mol−1\Delta G^\circ_{310} = -8.0\ \text{kJ mol}^{-1} at physiological temperature (310 K, 37 °C). Estimate ΔG∘\Delta G^\circ at 315 K (mild fever) and compare with the (wrong) assumption that ΔG∘\Delta G^\circ is temperature-independent. Step 1. Compute ΔG310∘/310\Delta G^\circ_{310}/310:

−8000 J mol−1310 K=−25.81 J mol−1K−1\frac{-8000\ \text{J mol}^{-1}}{310\ \text{K}} = -25.81\ \text{J mol}^{-1}\text{K}^{-1}

Step 2. Compute the bracket term:

1315−1310=0.0031746−0.0032258=−5.122×10−5 K−1\frac{1}{315} - \frac{1}{310} = 0.0031746 - 0.0032258 = -5.122\times10^{-5}\ \text{K}^{-1}

Step 3. Multiply by ΔH∘\Delta H^\circ:

(−180,000 J mol−1)(−5.122×10−5 K−1)=9.220 J mol−1K−1(-180{,}000\ \text{J mol}^{-1})(-5.122\times10^{-5}\ \text{K}^{-1}) = 9.220\ \text{J mol}^{-1}\text{K}^{-1}

Step 4. Add:

ΔG315∘315=−25.81+9.22=−16.59 J mol−1K−1\frac{\Delta G^\circ_{315}}{315} = -25.81 + 9.22 = -16.59\ \text{J mol}^{-1}\text{K}^{-1}
ΔG315∘=315×(−16.59)=−5226 J mol−1≈−5.2 kJ mol−1\Delta G^\circ_{315} = 315 \times (-16.59) = -5226\ \text{J mol}^{-1} \approx -5.2\ \text{kJ mol}^{-1}

Comparison. A naive assumption of constant ΔG∘=−8.0 kJ mol−1\Delta G^\circ = -8.0\ \text{kJ mol}^{-1} would badly overstate stability. The Gibbs–Helmholtz correction shows folding becomes noticeably less favorable, a swing of about 2.8 kJ/mol over just 5 K — because the large negative ΔH∘\Delta H^\circ makes ΔG/T\Delta G/T very sensitive to 1/T1/T. This size of shift is exactly why fever-range temperature changes can perturb protein stability biologically.

Limiting cases & common traps

SituationWhat happensWhy
ΔH=0\Delta H = 0ΔG/T\Delta G/T is constant with TTOnly entropy drives the temperature dependence, and ΔG=−TΔS\Delta G = -T\Delta S scales linearly
ΔCp≠0\Delta C_p \neq 0 over a wide rangeIntegrated formula becomes inaccurateAssumption of constant ΔH\Delta H breaks down; need Kirchhoff's law correction
Small T2−T1T_2 - T_1Linear approximation in 1/T1/T is excellentHigher-order terms in Taylor expansion are negligible
Confusing GG and G/TG/T derivativesWrong sign or magnitude(∂G/∂T)p=−S(\partial G/\partial T)_p = -S, NOT −H/T2-H/T^2; the H/T2H/T^2 form applies only to G/TG/T

The chemist behind it

The equation bears two names because it synthesizes work from both scientists on the relationship between heat and available (free) energy in chemical and physical processes during the 1870s–1880s. Josiah Willard Gibbs, working largely in isolation at Yale, developed the rigorous thermodynamic potential GG (originally his "available energy") and the mathematical apparatus of chemical thermodynamics almost single-handedly, publishing in the relatively obscure Transactions of the Connecticut Academy. Hermann von Helmholtz, already a towering figure in German physics and physiology, independently developed the closely related concept of free energy for physical systems and pushed the equation toward practical use in galvanic cells and chemical affinity. Their combined names reflect how a single deep insight, that G/TG/T rather than GG has this simple structure, emerged from two very different scientific traditions converging on the same thermodynamic truth.

Origins · 1 minJosiah Willard Gibbs and how this was discovered
JW
Josiah Willard Gibbs
1839-1903 · American mathematical physicist who founded chemical thermodynamics and defined the Gibbs free energy

Ask the tutor: How does the Gibbs–Helmholtz equation transform into the van't Hoff equation for predicting equilibrium constants at different temperatures?

Why divide by T instead of just using the plain derivative of G

The plain derivative, (∂G/∂T)p=−S(\partial G/\partial T)_p = -S, is correct but requires knowing entropy, which is harder to measure directly than enthalpy. Dividing by TT before differentiating is a deliberate algebraic trick: it makes the entropy terms cancel exactly, leaving a result that depends only on HH and TT. Since ΔH\Delta H is readily available from calorimetry or bond energies, this reformulation turns an entropy-dependent problem into an enthalpy-only problem, which is far more convenient experimentally.

Does the Gibbs–Helmholtz equation assume constant ΔH, or is that only for the integrated form

The differential form, (∂(ΔG/T)/∂T)p=−ΔH/T2\left(\partial(\Delta G/T)/\partial T\right)_p = -\Delta H/T^2, is exact at every temperature with no approximation. The approximation enters only when you integrate this over a finite range and want a simple closed-form answer. There, treating ΔH\Delta H as constant (ignoring ΔCp\Delta C_p) is the standard simplification, valid when the temperature range is modest or heat capacities of reactants and products are similar.

How does this equation connect to the van't Hoff equation

Since ΔG∘=−RTln⁡K\Delta G^\circ = -RT\ln K, dividing by TT gives ΔG∘/T=−Rln⁡K\Delta G^\circ/T = -R\ln K. Substituting this into the Gibbs–Helmholtz relation directly yields (∂ln⁡K/∂T)p=ΔH∘/(RT2)\left(\partial \ln K/\partial T\right)_p = \Delta H^\circ/(RT^2), which is precisely the van't Hoff equation. In this sense, the van't Hoff equation is simply the Gibbs–Helmholtz equation rewritten in terms of the equilibrium constant instead of the free energy itself.

Can this equation be applied to non-reaction processes, like phase transitions

Yes. The derivation never assumed a chemical reaction specifically; it only used G=H−TSG = H - TS and (∂G/∂T)p=−S(\partial G/\partial T)_p = -S, both general thermodynamic facts. For a phase transition such as melting or vaporization, the same equation applies with ΔHfus\Delta H_{fus} or ΔHvap\Delta H_{vap} in place of a reaction enthalpy, letting you predict how the free energy difference between phases shifts with temperature near the transition point.

What goes wrong if I ignore the ΔCp correction over a very wide temperature range

If ΔCp≠0\Delta C_p \neq 0, then ΔH\Delta H itself changes with temperature according to Kirchhoff's law, d(ΔH)/dT=ΔCpd(\Delta H)/dT = \Delta C_p. Using the constant-ΔH\Delta H integrated formula over a wide range then introduces systematic error, typically underestimating or overestimating ΔG\Delta G at the far end of the interval. For precise work spanning hundreds of kelvin, chemists use a temperature-dependent ΔH(T)\Delta H(T) inside the integral rather than treating it as constant.

Before this