Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

5.3. Ammonia Formation

Course-wide Conventions & Notation

Overview and Learning Objectives

Ammonia synthesis is one of the most important industrial chemical processes: the Haber–Bosch reaction produces the nitrogen feedstock for most synthetic fertilizers, supporting roughly half the world’s food production. It is also a textbook example of the tension between thermodynamics and kinetics. At room temperature the reaction is thermodynamically favorable (ΔGr∘<0\Delta G_r^\circ < 0) but kinetically inert. Raising the temperature makes the kinetics manageable but can push ΔGr∘\Delta G_r^\circ positive. Increasing pressure can restore favorability for reactions that decrease the number of gas molecules. Understanding and quantifying these tradeoffs is the goal of this section.

We use the free-energy machinery from Section 5.1 and the absolute-entropy framework from Section 5.2 to track how ΔGr∘\Delta G_r^\circ changes with temperature and pressure for the reaction

N2(g)+3 H2(g)→2 NH3(g).\mathrm{N_2(g) + 3\,H_2(g) \rightarrow 2\,NH_3(g)}.

The analysis relies on two tools: heat-capacity-based temperature corrections to ΔH\Delta H and ΔS\Delta S, and the ideal-gas pressure dependence of GG.


Learning objectives:

Core Ideas and Derivations

Thermodynamics vs. kinetics

Standard thermodynamic data at 298.15 K

The standard reaction quantities for ammonia synthesis at T=298.15 KT = 298.15\ \mathrm{K} are:

ΔHr∘(298.15 K)≈−92 kJ mol−1\Delta H_r^\circ(298.15\,\mathrm{K}) \approx -92\ \mathrm{kJ\,mol^{-1}}
ΔSr∘(298.15 K)≈−198 J K−1 mol−1\Delta S_r^\circ(298.15\,\mathrm{K}) \approx -198\ \mathrm{J\,K^{-1}\,mol^{-1}}
ΔGr∘(298.15 K)≈−33 kJ mol−1\Delta G_r^\circ(298.15\,\mathrm{K}) \approx -33\ \mathrm{kJ\,mol^{-1}}

These are consistent via ΔGr∘=ΔHr∘−T ΔSr∘\Delta G_r^\circ = \Delta H_r^\circ - T\,\Delta S_r^\circ. Since ΔGr∘<0\Delta G_r^\circ < 0, the reaction is thermodynamically spontaneous at 298 K and 1 bar.

Why it is too slow at room temperature

Even with ΔGr∘<0\Delta G_r^\circ < 0, the reaction can be extremely slow if the activation barrier is large. For ammonia formation, the N≡N\mathrm{N \equiv N} triple bond (bond dissociation energy ≈945 kJ mol−1\approx 945\ \mathrm{kJ\,mol^{-1}}) makes the uncatalyzed barrier prohibitively high at 298 K.

A catalyst (e.g., an Fe-based catalyst developed by Haber and Bosch, ca. 1909; Nobel Prize 1918) provides an alternative reaction pathway with a lower activation barrier. The catalyst accelerates both the forward and reverse reactions equally, so it does not change ΔGr∘\Delta G_r^\circ or the equilibrium constant — it only changes how fast equilibrium is reached. Raising the temperature further increases the rate (Arrhenius behavior), so industrial reactors operate at 400–500 °C400\text{–}500\,\text{°C} where the barrier is surmountable.

But raising TT introduces a thermodynamic penalty, which we now quantify.


Qualitative effect of temperature

Since ΔSr∘<0\Delta S_r^\circ < 0 for this reaction, the −T ΔSr∘-T\,\Delta S_r^\circ term in

ΔGr∘=ΔHr∘−T ΔSr∘\Delta G_r^\circ = \Delta H_r^\circ - T\,\Delta S_r^\circ

is positive and grows with TT. At sufficiently high temperature, this positive term overwhelms the negative ΔHr∘\Delta H_r^\circ, and ΔGr∘\Delta G_r^\circ changes sign — the reaction becomes non-spontaneous at 1 bar.

The qualitative rule is: for reactions with ΔSr∘<0\Delta S_r^\circ < 0, higher temperatures work against thermodynamic favorability.


Temperature corrections from CPC_P

To quantify the temperature effect, we need ΔHr∘(T)\Delta H_r^\circ(T) and ΔSr∘(T)\Delta S_r^\circ(T) at the operating temperature. The strategy is to correct each species’ enthalpy and entropy from the reference temperature Ti=298.15 KT_i = 298.15\ \mathrm{K} to the target temperature TfT_f using heat-capacity data.

Enthalpy correction

At constant pressure (Section 3.3):

H(Tf)−H(Ti)=∫TiTfCP(T) dTH(T_f) - H(T_i) = \int_{T_i}^{T_f} C_P(T)\,dT

Entropy correction

From the enthalpy differential dH=T dS+V dPdH = T\,dS + V\,dP at constant PP (dP=0dP = 0):

S(Tf)−S(Ti)=∫TiTfCP(T)T dTS(T_f) - S(T_i) = \int_{T_i}^{T_f} \frac{C_P(T)}{T}\,dT

This is the same integral that appears in the absolute-entropy formula (Section 5.2, Eq. (11)), applied over a finite temperature range rather than from 0.


Ideal-gas heat capacities from degrees of freedom

For a simple estimate we treat N2_2, H2_2, and NH3_3 as ideal gases and approximate CPC_P as temperature-independent, using the number of active (translational + rotational) degrees of freedom ff:

CP=(f2+1)RC_P = \left(\frac{f}{2} + 1\right)R

For linear molecules (N2_2, H2_2): f=3 (trans)+2 (rot)=5f = 3\ (\text{trans}) + 2\ (\text{rot}) = 5, giving

CP=72R=29.1 J mol−1 K−1.C_P = \tfrac{7}{2}R = 29.1\ \mathrm{J\,mol^{-1}\,K^{-1}}.

For nonlinear NH3_3: f=3 (trans)+3 (rot)=6f = 3\ (\text{trans}) + 3\ (\text{rot}) = 6, giving

CP=4R=33.3 J mol−1 K−1.C_P = 4R = 33.3\ \mathrm{J\,mol^{-1}\,K^{-1}}.

With constant CPC_P, the integrals simplify to:

H(Tf)−H(Ti)=CP(Tf−Ti)H(T_f) - H(T_i) = C_P(T_f - T_i)
S(Tf)−S(Ti)=CPln⁡ ⁣(TfTi)S(T_f) - S(T_i) = C_P\ln\!\left(\frac{T_f}{T_i}\right)

Temperature effect: 298.15 K to 773.15 K (500°C)

Set Ti=298.15 KT_i = 298.15\ \mathrm{K} and Tf=773.15 KT_f = 773.15\ \mathrm{K}.

Per-species corrections

N2_2 and H2_2 (CP=72RC_P = \tfrac{7}{2}R):

ΔHspecies=72R (773.15−298.15)≈13.8 kJ mol−1\Delta H_{\text{species}} = \tfrac{7}{2}R\,(773.15 - 298.15) \approx 13.8\ \mathrm{kJ\,mol^{-1}}
ΔSspecies=72R ln⁡ ⁣(773.15298.15)≈27.7 J K−1 mol−1\Delta S_{\text{species}} = \tfrac{7}{2}R\,\ln\!\left(\frac{773.15}{298.15}\right) \approx 27.7\ \mathrm{J\,K^{-1}\,mol^{-1}}

NH3_3 (CP=4RC_P = 4R):

ΔHspecies=4R (773.15−298.15)≈15.8 kJ mol−1\Delta H_{\text{species}} = 4R\,(773.15 - 298.15) \approx 15.8\ \mathrm{kJ\,mol^{-1}}
ΔSspecies=4R ln⁡ ⁣(773.15298.15)≈31.7 J K−1 mol−1\Delta S_{\text{species}} = 4R\,\ln\!\left(\frac{773.15}{298.15}\right) \approx 31.7\ \mathrm{J\,K^{-1}\,mol^{-1}}

Reaction-level corrections

Apply stoichiometric weights (2×NH3−1×N2−3×H22 \times \text{NH}_3 - 1 \times \text{N}_2 - 3 \times \text{H}_2):

ΔHr∘(Tf)≈ΔHr∘(Ti)+[2(15.8)−(13.8)−3(13.8)]=ΔHr∘(Ti)−23.6 kJ mol−1\Delta H_r^\circ(T_f) \approx \Delta H_r^\circ(T_i) + \big[2(15.8) - (13.8) - 3(13.8)\big] = \Delta H_r^\circ(T_i) - 23.6\ \mathrm{kJ\,mol^{-1}}
ΔSr∘(Tf)≈ΔSr∘(Ti)+[2(31.7)−(27.7)−3(27.7)]=ΔSr∘(Ti)−47.4 J K−1 mol−1\Delta S_r^\circ(T_f) \approx \Delta S_r^\circ(T_i) + \big[2(31.7) - (27.7) - 3(27.7)\big] = \Delta S_r^\circ(T_i) - 47.4\ \mathrm{J\,K^{-1}\,mol^{-1}}

The reaction-level ΔCP\Delta C_P is negative (products have fewer degrees of freedom than reactants), so both ΔHr∘\Delta H_r^\circ and ΔSr∘\Delta S_r^\circ become more negative with increasing temperature.

ΔGr∘\Delta G_r^\circ at 773.15 K

Combining the corrected values:

ΔHr∘(773)=−92.0+(−23.6)=−115.6 kJ mol−1\Delta H_r^\circ(773) = -92.0 + (-23.6) = -115.6\ \mathrm{kJ\,mol^{-1}}
ΔSr∘(773)=−198+(−47.4)=−245.4 J K−1 mol−1\Delta S_r^\circ(773) = -198 + (-47.4) = -245.4\ \mathrm{J\,K^{-1}\,mol^{-1}}
ΔGr∘(773)=−115.6−(773.15)(−0.2454)≈+74.1 kJ mol−1\Delta G_r^\circ(773) = -115.6 - (773.15)(-0.2454) \approx +74.1\ \mathrm{kJ\,mol^{-1}}

At 500 °C500\,\text{°C} and 1 bar, ammonia formation is non-spontaneous (ΔGr∘>0\Delta G_r^\circ > 0), even though it would be much faster kinetically than at room temperature. This is the central dilemma of the Haber–Bosch process.


Pressure dependence of GG for an ideal gas

Temperature is not the only knob available. For gas-phase reactions, pressure can strongly affect ΔG\Delta G when the reaction changes the number of gas molecules.

From the Gibbs differential dG=−S dT+V dPdG = -S\,dT + V\,dP (Section 5.1) at constant temperature:

dG∣T=V dPdG\big|_T = V\,dP

For one mole of an ideal gas, V=RT/PV = RT/P. Integrating from P∘P^\circ to PP:

G(T,P)−G∘(T)=RTln⁡ ⁣(PP∘)G(T, P) - G^\circ(T) = RT\ln\!\left(\frac{P}{P^\circ}\right)

where G∘(T)G^\circ(T) is the standard-state molar Gibbs energy at pressure P∘=1 barP^\circ = 1\ \mathrm{bar}.

Applying to ammonia formation

If each species is compressed from P∘P^\circ to PP at fixed TT, the pressure correction to the reaction Gibbs energy is

ΔGr(T,P)≈ΔGr∘(T)+Δν RTln⁡ ⁣(PP∘),\Delta G_r(T, P) \approx \Delta G_r^\circ(T) + \Delta\nu\,RT\ln\!\left(\frac{P}{P^\circ}\right),

where Δν\Delta\nu is the change in the number of moles of gas:

Δν=2−(1+3)=−2.\Delta\nu = 2 - (1 + 3) = -2.

Setting ΔGr=0\Delta G_r = 0 to find the pressure at which the reaction is just spontaneous:

P=P∘exp⁡ ⁣(−ΔGr∘(T)Δν RT)P = P^\circ \exp\!\left(-\frac{\Delta G_r^\circ(T)}{\Delta\nu\,RT}\right)

Inserting ΔGr∘(773)≈+74.1 kJ mol−1\Delta G_r^\circ(773) \approx +74.1\ \mathrm{kJ\,mol^{-1}}, Δν=−2\Delta\nu = -2, and T=773.15 KT = 773.15\ \mathrm{K}:

ln⁡ ⁣(PP∘)=−74.1×103(−2)(8.314)(773.15)=5.77⇒P≈320 bar\ln\!\left(\frac{P}{P^\circ}\right) = -\frac{74.1 \times 10^3}{(-2)(8.314)(773.15)} = 5.77 \quad\Rightarrow\quad P \approx 320\ \mathrm{bar}

Because Δν<0\Delta\nu < 0, increasing pressure favors the side with fewer gas molecules (Le Châtelier’s principle). At 500 °C500\,\text{°C}, pressures on the order of several hundred bar can restore thermodynamic favorability.


Summary: the engineering knobs

The Haber–Bosch process operates at 400–500 °C400\text{–}500\,\text{°C} and 150–350 bar150\text{–}350\ \mathrm{bar} over an iron-based catalyst. Each operating parameter addresses a different aspect of the thermodynamics/kinetics tradeoff:

Industrial conditions represent a compromise: high enough temperature for acceptable rates, high enough pressure for acceptable equilibrium yield, with a catalyst to make the whole process feasible.


Worked Examples

Example 1: Crossover temperature at 1 bar

Problem. At what temperature does ΔGr∘\Delta G_r^\circ for ammonia formation change sign (from negative to positive) at 1 bar? Estimate this (a) neglecting CPC_P corrections, and (b) including the constant-CPC_P corrections from this section.

Solution.

(a) Without CPC_P corrections. If we treat ΔHr∘\Delta H_r^\circ and ΔSr∘\Delta S_r^\circ as temperature-independent, the crossover occurs when ΔGr∘=ΔHr∘−T ΔSr∘=0\Delta G_r^\circ = \Delta H_r^\circ - T\,\Delta S_r^\circ = 0, giving

Tcross=ΔHr∘ΔSr∘=−92,000 J mol−1−198 J K−1 mol−1=465 K≈192 °C.T_{\mathrm{cross}} = \frac{\Delta H_r^\circ}{\Delta S_r^\circ} = \frac{-92{,}000\ \mathrm{J\,mol^{-1}}}{-198\ \mathrm{J\,K^{-1}\,mol^{-1}}} = 465\ \mathrm{K} \approx 192\,\text{°C}.

(b) With CPC_P corrections. Including the temperature-dependent corrections from this section, the reaction-level ΔCP=2CP,NH3−CP,N2−3CP,H2=2(4R)−(7R/2)−3(7R/2)=−6R≈−49.9 J K−1 mol−1\Delta C_P = 2C_{P,\mathrm{NH_3}} - C_{P,\mathrm{N_2}} - 3C_{P,\mathrm{H_2}} = 2(4R) - (7R/2) - 3(7R/2) = -6R \approx -49.9\ \mathrm{J\,K^{-1}\,mol^{-1}}. We need to solve

ΔHr∘(Ti)+ΔCP(T−Ti)−T ⁣[ΔSr∘(Ti)+ΔCPln⁡ ⁣(TTi)]=0,\Delta H_r^\circ(T_i) + \Delta C_P(T - T_i) - T\!\left[\Delta S_r^\circ(T_i) + \Delta C_P\ln\!\left(\frac{T}{T_i}\right)\right] = 0,

which cannot be solved in closed form but can be solved numerically. The result is

Tcross≈456 K≈183 °C.T_{\mathrm{cross}} \approx 456\ \mathrm{K} \approx 183\,\text{°C}.

Result. The reaction becomes non-spontaneous (at 1 bar) around 456–465 K456\text{–}465\ \mathrm{K}, depending on whether CPC_P corrections are included. Either way, the crossover is well below the 400–500 °C400\text{–}500\,\text{°C} operating range of Haber–Bosch, confirming that pressure is essential to restore thermodynamic favorability at operating temperatures.


Example 2: Pressure needed at 400°C

Problem. Repeat the pressure-for-spontaneity estimate at T=673.15 KT = 673.15\ \mathrm{K} (400 °C400\,\text{°C}) instead of 500 °C500\,\text{°C}.

Solution.

First, compute ΔGr∘(673 K)\Delta G_r^\circ(673\ \mathrm{K}) using the same constant-CPC_P framework. The temperature shift is ΔT=673.15−298.15=375.0 K\Delta T = 673.15 - 298.15 = 375.0\ \mathrm{K}:

ΔHr∘(673)=−92.0+(−6R)(375.0)=−92.0−18.7=−110.7 kJ mol−1\Delta H_r^\circ(673) = -92.0 + (-6R)(375.0) = -92.0 - 18.7 = -110.7\ \mathrm{kJ\,mol^{-1}}
ΔSr∘(673)=−198+(−6R)ln⁡ ⁣(673.15298.15)=−198−40.7=−238.7 J K−1 mol−1\Delta S_r^\circ(673) = -198 + (-6R)\ln\!\left(\frac{673.15}{298.15}\right) = -198 - 40.7 = -238.7\ \mathrm{J\,K^{-1}\,mol^{-1}}
ΔGr∘(673)=−110.7−(673.15)(−0.2387)=−110.7+160.7=+50.0 kJ mol−1\Delta G_r^\circ(673) = -110.7 - (673.15)(-0.2387) = -110.7 + 160.7 = +50.0\ \mathrm{kJ\,mol^{-1}}

Now use Eq. (26):

ln⁡ ⁣(P1 bar)=−50.0×103(−2)(8.314)(673.15)=4.47⇒P≈87 bar.\ln\!\left(\frac{P}{1\ \mathrm{bar}}\right) = -\frac{50.0 \times 10^3}{(-2)(8.314)(673.15)} = 4.47 \quad\Rightarrow\quad P \approx 87\ \mathrm{bar}.

Result. At 400 °C400\,\text{°C}, only about 87 bar is needed for spontaneity, compared with ∼\sim320 bar at 500 °C500\,\text{°C}. This illustrates the tradeoff: a lower operating temperature demands less pressure for thermodynamic favorability but gives a slower reaction rate. The industrial compromise (400–500 °C400\text{–}500\,\text{°C}, 150–350 bar150\text{–}350\ \mathrm{bar}) sits squarely in the range our estimates predict.

Concept Checks

  1. Why does ΔSr∘<0\Delta S_r^\circ < 0 make higher temperatures less favorable for ammonia formation at equilibrium?

  2. How does the sign of Δν\Delta\nu determine whether increasing pressure helps or hurts product formation?

  3. Which approximations are embedded in using G(P)−G∘=RTln⁡(P/P∘)G(P) - G^\circ = RT\ln(P/P^\circ) for each species? How would the analysis change for a non-ideal gas?

  4. A catalyst accelerates both the forward and reverse reactions. Why does this not shift the equilibrium?

  5. An engineer proposes running at 300°C (where ΔGr∘<0\Delta G_r^\circ < 0 at 1 bar) to avoid the need for high pressure. What practical problem does this create?

Key Takeaways