Skip to main content

Entropy and Gibbs Free Energy

Enthalpy alone is insufficient to predict whether a reaction will occur spontaneously. The dissolution of ammonium nitrate in water is endothermic yet proceeds spontaneously. The decomposition of calcium carbonate requires continuous heating. To resolve these cases, we need a second state function: entropy, and its combination with enthalpy in the Gibbs free energy.

Entropy (SS)

Statistical Definition

Entropy is a measure of the number of microstates (WW) accessible to a system -- that is, the number of ways the energy of the system can be distributed among its particles:

S=kBlnWS = k_B \ln W

where kB=1.381×1023J/Kk_B = 1.381 \times 10^{-23}\,\mathrm{J/K} is the Boltzmann constant. A system with more accessible microstates has higher entropy. This is the Boltzmann equation, engraved on his tombstone.

Thermodynamic Definition

The entropy change is defined as:

ΔS=LBδqrevRB◆◆LBTRB\Delta S = \int \frac◆LB◆\delta q_\mathrm{rev}◆RB◆◆LB◆T◆RB◆

For an isothermal process:

ΔS=LBqrevRB◆◆LBTRB\Delta S = \frac◆LB◆q_\mathrm{rev}◆RB◆◆LB◆T◆RB◆

This shows that entropy is heat divided by temperature. Transferring a given amount of heat at a lower temperature produces a larger entropy change than at a higher temperature.

Factors Affecting Entropy

FactorEffect on SSRationale
More gas molecules producedIncreaseMore ways to distribute energy among more particles
Change from solid to liquid to gasIncreaseGases have the most microstates; solids the fewest
Higher temperatureIncreaseMore energy available to distribute; more accessible microstates
Dissolution of a solid ionic compoundIncrease (usually)Ions dispersed in solution have more freedom than in the lattice
Fewer moles of gasDecreaseFewer particles, fewer microstates
Increased pressure of a gasDecreaseReduced volume constrains the available positions

Standard Entropy Values (SS^\circ)

Standard entropies are absolute values (not relative to a reference, unlike enthalpy). At 0K0\,\mathrm{K}, a perfect crystal has S=0S = 0 (third law of thermodynamics). Standard entropies are always positive.

Typical values (Jmol1K1\mathrm{J\,mol^{-1}\,K^{-1}}):

SubstanceSS^\circ
C(diamond)\mathrm{C}(\mathrm{diamond})2.4
C(graphite)\mathrm{C}(\mathrm{graphite})5.7
Fe(s)\mathrm{Fe}(s)27.3
NaCl(s)\mathrm{NaCl}(s)72.1
H2O(l)\mathrm{H}_2\mathrm{O}(l)69.9
H2O(g)\mathrm{H}_2\mathrm{O}(g)188.8
CO2(g)\mathrm{CO}_2(g)213.7
N2(g)\mathrm{N}_2(g)191.6

Gases have much higher standard entropies than liquids and solids. Diamond has a lower entropy than graphite due to its more rigid, ordered structure.

Standard Entropy Change

ΔS=S(products)S(reactants)\Delta S^\circ = \sum S^\circ(\mathrm{products}) - \sum S^\circ(\mathrm{reactants})

Worked Example. Calculate ΔS\Delta S^\circ for the thermal decomposition of calcium carbonate:

CaCO3(s)CaO(s)+CO2(g)\mathrm{CaCO}_3(s) \to \mathrm{CaO}(s) + \mathrm{CO}_2(g) ΔS=S(CaO)+S(CO2)S(CaCO3)=38.1+213.792.9=+158.9Jmol1K1\Delta S^\circ = S^\circ(\mathrm{CaO}) + S^\circ(\mathrm{CO}_2) - S^\circ(\mathrm{CaCO}_3) = 38.1 + 213.7 - 92.9 = +158.9\,\mathrm{J\,mol^{-1}\,K^{-1}}

The entropy increases because a gas is produced from a solid, creating many more microstates.

The Second Law of Thermodynamics

The total entropy of the universe increases in any spontaneous process:

ΔSuniverse=ΔSsystem+ΔSsurroundings>0\Delta S_\mathrm{universe} = \Delta S_\mathrm{system} + \Delta S_\mathrm{surroundings} \gt 0

At equilibrium: ΔSuniverse=0\Delta S_\mathrm{universe} = 0.

Entropy Change of the Surroundings

The surroundings gain or lose heat as a result of the reaction. The entropy change of the surroundings is:

ΔSsurroundings=LBΔHsystemRB◆◆LBTRB\Delta S_\mathrm{surroundings} = -\frac◆LB◆\Delta H_\mathrm{system}◆RB◆◆LB◆T◆RB◆

The negative sign arises because when the system releases heat (exothermic, ΔH<0\Delta H \lt 0), the surroundings gain that heat and their entropy increases. When the system absorbs heat (endothermic, ΔH>0\Delta H \gt 0), the surroundings lose heat and their entropy decreases.

Combining with the second law:

ΔSuniverse=ΔSsystemLBΔHsystemRB◆◆LBTRB>0\Delta S_\mathrm{universe} = \Delta S_\mathrm{system} - \frac◆LB◆\Delta H_\mathrm{system}◆RB◆◆LB◆T◆RB◆ \gt 0

Multiplying through by TT (which is always positive):

TΔSsystemΔHsystem>0T\Delta S_\mathrm{system} - \Delta H_\mathrm{system} \gt 0

Rearranging:

ΔHsystemTΔSsystem<0\Delta H_\mathrm{system} - T\Delta S_\mathrm{system} \lt 0

Gibbs Free Energy (GG)

Derivation

The inequality above defines the Gibbs free energy:

ΔG=ΔHTΔS\Delta G = \Delta H - T\Delta S

This is the most important equation in chemical thermodynamics. It combines the enthalpy and entropy contributions into a single quantity that determines spontaneity.

Spontaneity Criterion

ΔG\Delta GProcess
ΔG<0\Delta G \lt 0Spontaneous (thermodynamically favourable)
ΔG=0\Delta G = 0At equilibrium
ΔG>0\Delta G \gt 0Non-spontaneous (thermodynamically unfavourable)

A spontaneous process is one that proceeds without external intervention once initiated. "Spontaneous" does not mean "fast" -- kinetics determines the rate; thermodynamics determines the direction.

Standard Gibbs Free Energy Change

ΔG=ΔHTΔS\Delta G^\circ = \Delta H^\circ - T\Delta S^\circ

This uses standard enthalpy and entropy data at 298K298\,\mathrm{K}.

Relationship to the Equilibrium Constant

ΔG=RTlnK\Delta G^\circ = -RT\ln K

where R=8.314Jmol1K1R = 8.314\,\mathrm{J\,mol^{-1}\,K^{-1}}, TT is in Kelvin, and KK is the equilibrium constant (dimensionless, using activities).

This equation is one of the most powerful in chemistry because it connects thermodynamic data to measurable equilibrium constants.

ΔG\Delta G^\circKKEquilibrium Position
ΔG<0\Delta G^\circ \lt 0K>1K \gt 1Products favoured
ΔG=0\Delta G^\circ = 0K=1K = 1Neither favoured
ΔG>0\Delta G^\circ \gt 0K<1K \lt 1Reactants favoured

Worked Example. For the Haber process at 298K298\,\mathrm{K}:

N2(g)+3H2(g)2NH3(g)\mathrm{N}_2(g) + 3\mathrm{H}_2(g) \rightleftharpoons 2\mathrm{NH}_3(g)

ΔH=92kJ/mol\Delta H^\circ = -92\,\mathrm{kJ/mol}, ΔS=199Jmol1K1\Delta S^\circ = -199\,\mathrm{J\,mol^{-1}\,K^{-1}}.

ΔG=92000298×(199)=92000+59302=32698J/mol=32.7kJ/mol\Delta G^\circ = -92000 - 298 \times (-199) = -92000 + 59302 = -32698\,\mathrm{J/mol} = -32.7\,\mathrm{kJ/mol}

Since ΔG<0\Delta G^\circ \lt 0, the reaction is spontaneous at 298K298\,\mathrm{K} (but kinetically extremely slow without a catalyst).

K=exp(LBΔGRB◆◆LBRTRB)=exp(LB32698RB◆◆LB8.314×298RB)=exp(13.19)=5.3×105K = \exp\left(\frac◆LB◆-\Delta G^\circ◆RB◆◆LB◆RT◆RB◆\right) = \exp\left(\frac◆LB◆32698◆RB◆◆LB◆8.314 \times 298◆RB◆\right) = \exp(13.19) = 5.3 \times 10^5

The equilibrium constant is very large, confirming that products are strongly favoured at 298K298\,\mathrm{K}.

Temperature Dependence of Feasibility

The ΔG=ΔHTΔS\Delta G = \Delta H - T\Delta S equation reveals that the spontaneity of a reaction can change with temperature. The four cases:

ΔH\Delta HΔS\Delta SLow TTHigh TTExample
-++SpontaneousSpontaneousCombustion of hydrogen
--SpontaneousNon-spontaneousFreezing of water
++++Non-spontaneousSpontaneousThermal decomposition of CaCO3\mathrm{CaCO}_3
++-Non-spontaneousNon-spontaneousDecomposition of N2O\mathrm{N}_2\mathrm{O}

The Temperature of Equilibrium

The temperature at which a reaction changes from spontaneous to non-spontaneous (or vice versa) is the temperature at which ΔG=0\Delta G = 0:

T=LBΔHRB◆◆LBΔSRBT = \frac◆LB◆\Delta H^\circ◆RB◆◆LB◆\Delta S^\circ◆RB◆

Worked Example. At what temperature does the thermal decomposition of calcium carbonate become spontaneous?

CaCO3(s)CaO(s)+CO2(g)\mathrm{CaCO}_3(s) \to \mathrm{CaO}(s) + \mathrm{CO}_2(g)

ΔH=+178kJ/mol\Delta H^\circ = +178\,\mathrm{kJ/mol}, ΔS=+160Jmol1K1\Delta S^\circ = +160\,\mathrm{J\,mol^{-1}\,K^{-1}}.

T=LB178×103RB◆◆LB160RB=1113K=840CT = \frac◆LB◆178 \times 10^3◆RB◆◆LB◆160◆RB◆ = 1113\,\mathrm{K} = 840^\circ\mathrm{C}

Above 1113K1113\,\mathrm{K}, ΔG<0\Delta G^\circ \lt 0 and decomposition is spontaneous. In a lime kiln, temperatures of approximately 900900--1000C1000^\circ\mathrm{C} are used to ensure thermodynamic feasibility while maintaining a practical rate.

Graphical Interpretation

A plot of ΔG\Delta G^\circ vs TT is a straight line with slope =ΔS= -\Delta S^\circ and y-intercept =ΔH= \Delta H^\circ (at T=0T = 0).

  • When ΔS>0\Delta S^\circ \gt 0: the line slopes downward. The reaction becomes more spontaneous as temperature increases.
  • When ΔS<0\Delta S^\circ \lt 0: the line slopes upward. The reaction becomes less spontaneous as temperature increases.
  • The x-intercept (ΔG=0\Delta G^\circ = 0) gives the equilibrium temperature.

Important caveat: This linear extrapolation assumes that ΔH\Delta H^\circ and ΔS\Delta S^\circ are independent of temperature (Kirchhoff's approximation). This is a reasonable approximation over small temperature ranges but fails over large ranges where heat capacities change significantly.

Non-Standard Conditions: ΔG\Delta G vs ΔG\Delta G^\circ

The standard free energy change (ΔG\Delta G^\circ) applies when all reactants and products are in their standard states (1 mol/dm3^3 for solutions, 1 bar for gases). Under non-standard conditions:

ΔG=ΔG+RTlnQ\Delta G = \Delta G^\circ + RT\ln Q

where QQ is the reaction quotient (the same expression as KK but with current, non-equilibrium concentrations or partial pressures).

At equilibrium, Q=KQ = K and ΔG=0\Delta G = 0, recovering ΔG=RTlnK\Delta G^\circ = -RT\ln K.

Industrial Applications

The Haber Process

N2(g)+3H2(g)2NH3(g)ΔH=92kJ/mol\mathrm{N}_2(g) + 3\mathrm{H}_2(g) \rightleftharpoons 2\mathrm{NH}_3(g) \quad \Delta H^\circ = -92\,\mathrm{kJ/mol}

ΔH<0\Delta H^\circ \lt 0 (exothermic), ΔS<0\Delta S^\circ \lt 0 (4 moles of gas to 2 moles). By Le Chatelier's principle and the Gibbs equation:

  • Low temperature favours the forward reaction (exothermic). But low temperature gives a slow rate.
  • High pressure favours the forward reaction (fewer gas moles on the product side).
  • Compromise: 450C450^\circ\mathrm{C}, 200atm200\,\mathrm{atm}, with an iron catalyst.

The Contact Process

2SO2(g)+O2(g)2SO3(g)ΔH=198kJ/mol2\mathrm{SO}_2(g) + \mathrm{O}_2(g) \rightleftharpoons 2\mathrm{SO}_3(g) \quad \Delta H^\circ = -198\,\mathrm{kJ/mol}

ΔH<0\Delta H^\circ \lt 0, ΔS<0\Delta S^\circ \lt 0 (3 moles to 2 moles). Low temperature favours the product but slows the rate. Industrial conditions: 450C450^\circ\mathrm{C}, 11--2atm2\,\mathrm{atm}, V2O5\mathrm{V}_2\mathrm{O}_5 catalyst.

Extraction of Iron

Fe2O3(s)+3CO(g)2Fe(l)+3CO2(g)ΔH=23kJ/mol\mathrm{Fe}_2\mathrm{O}_3(s) + 3\mathrm{CO}(g) \rightleftharpoons 2\mathrm{Fe}(l) + 3\mathrm{CO}_2(g) \quad \Delta H^\circ = -23\,\mathrm{kJ/mol}

ΔH<0\Delta H^\circ \lt 0 but small. The entropy change is favourable (3 moles of gas reactants to 3 moles of gas products, but the solid is consumed). At the blast furnace temperature (1500C\approx 1500^\circ\mathrm{C}), the reaction is thermodynamically feasible.

Common Pitfalls

  1. Unit mismatch in the Gibbs equation. ΔH\Delta H is typically in kJ/mol\mathrm{kJ/mol}, while ΔS\Delta S is in Jmol1K1\mathrm{J\,mol^{-1}\,K^{-1}}. Always convert to consistent units before combining: either convert ΔH\Delta H to J/mol\mathrm{J/mol} or ΔS\Delta S to kJmol1K1\mathrm{kJ\,mol^{-1}\,K^{-1}}.

  2. Using ΔG<0\Delta G \lt 0 to predict rate. Thermodynamics says nothing about kinetics. A reaction with ΔG0\Delta G \ll 0 may be immeasurably slow (e.g. diamond conversion to graphite at room temperature: ΔG2.9kJ/mol\Delta G^\circ \approx -2.9\,\mathrm{kJ/mol}, but the half-life is effectively infinite).

  3. Forgetting that ΔS\Delta S^\circ values are absolute. Unlike ΔHf\Delta H_f^\circ (which is relative to elements in standard states), SS^\circ values are absolute entropies. SS^\circ of an element in its standard state is not zero (except at 0K0\,\mathrm{K}).

  4. Confusing ΔG\Delta G with ΔG\Delta G^\circ. ΔG\Delta G^\circ is the free energy change under standard conditions. The actual free energy change ΔG\Delta G depends on the specific concentrations/pressures and is given by ΔG=ΔG+RTlnQ\Delta G = \Delta G^\circ + RT\ln Q.

  5. Assuming the linear ΔG\Delta G vs TT relationship holds indefinitely. The equation ΔG=ΔHTΔS\Delta G = \Delta H - T\Delta S assumes ΔH\Delta H and ΔS\Delta S are temperature-independent. Over large temperature ranges, this approximation fails.

Practice Problems

Problem 1

For the reaction C(s)+H2O(g)CO(g)+H2(g)\mathrm{C}(s) + \mathrm{H}_2\mathrm{O}(g) \to \mathrm{CO}(g) + \mathrm{H}_2(g):

ΔH=+131kJ/mol\Delta H^\circ = +131\,\mathrm{kJ/mol}, ΔS=+134Jmol1K1\Delta S^\circ = +134\,\mathrm{J\,mol^{-1}\,K^{-1}}.

(a) Calculate ΔG\Delta G^\circ at 298K298\,\mathrm{K}. Is the reaction spontaneous? (b) At what temperature does the reaction become spontaneous? (c) Calculate KK at 1000K1000\,\mathrm{K}.

Solution:

(a)

ΔG=131000298×134=13100039932=+91068J/mol=+91.1kJ/mol\Delta G^\circ = 131000 - 298 \times 134 = 131000 - 39932 = +91068\,\mathrm{J/mol} = +91.1\,\mathrm{kJ/mol}

ΔG>0\Delta G^\circ \gt 0: not spontaneous at 298K298\,\mathrm{K}.

(b)

T=LBΔHRB◆◆LBΔSRB=131000134=978K=705CT = \frac◆LB◆\Delta H^\circ◆RB◆◆LB◆\Delta S^\circ◆RB◆ = \frac{131000}{134} = 978\,\mathrm{K} = 705^\circ\mathrm{C}

Above 978K978\,\mathrm{K}, the reaction becomes spontaneous.

(c) At 1000K1000\,\mathrm{K}:

ΔG=1310001000×134=131000134000=3000J/mol\Delta G^\circ = 131000 - 1000 \times 134 = 131000 - 134000 = -3000\,\mathrm{J/mol}K=exp(LBΔGRB◆◆LBRTRB)=exp(LB3000RB◆◆LB8.314×1000RB)=exp(0.361)=1.43K = \exp\left(\frac◆LB◆-\Delta G^\circ◆RB◆◆LB◆RT◆RB◆\right) = \exp\left(\frac◆LB◆3000◆RB◆◆LB◆8.314 \times 1000◆RB◆\right) = \exp(0.361) = 1.43

K=1.43K = 1.43, so products are slightly favoured at 1000K1000\,\mathrm{K}.

Problem 2

The melting of ice: H2O(s)H2O(l)\mathrm{H}_2\mathrm{O}(s) \to \mathrm{H}_2\mathrm{O}(l) has ΔH=+6.01kJ/mol\Delta H^\circ = +6.01\,\mathrm{kJ/mol} and ΔS=+22.0Jmol1K1\Delta S^\circ = +22.0\,\mathrm{J\,mol^{-1}\,K^{-1}}.

(a) Calculate the normal melting point of ice. (b) Explain why ice melts spontaneously at room temperature despite the process being endothermic.

Solution:

(a) At the melting point, ΔG=0\Delta G = 0:

T=LBΔHRB◆◆LBΔSRB=601022.0=273K=0CT = \frac◆LB◆\Delta H^\circ◆RB◆◆LB◆\Delta S^\circ◆RB◆ = \frac{6010}{22.0} = 273\,\mathrm{K} = 0^\circ\mathrm{C}

(b) Although ΔH>0\Delta H^\circ \gt 0 (endothermic), ΔS>0\Delta S^\circ \gt 0 (entropy increases). At temperatures above 273K273\,\mathrm{K}, the TΔST\Delta S term exceeds ΔH\Delta H, making ΔG<0\Delta G \lt 0. The entropy gain from the increased disorder of the liquid phase more than compensates for the enthalpy cost of breaking the hydrogen-bonded lattice.

Problem 3

For the reaction 2NO2(g)N2O4(g)2\mathrm{NO}_2(g) \rightleftharpoons \mathrm{N}_2\mathrm{O}_4(g) at 298K298\,\mathrm{K}:

ΔH=57.2kJ/mol\Delta H^\circ = -57.2\,\mathrm{kJ/mol}, ΔS=175.8Jmol1K1\Delta S^\circ = -175.8\,\mathrm{J\,mol^{-1}\,K^{-1}}.

(a) Calculate ΔG\Delta G^\circ and KpK_p at 298K298\,\mathrm{K}. (b) At what temperature does Kp=1K_p = 1? (c) Explain qualitatively whether increasing temperature increases or decreases the yield of N2O4\mathrm{N}_2\mathrm{O}_4.

Solution:

(a)

ΔG=57200298×(175.8)=57200+52388=4812J/mol=4.81kJ/mol\Delta G^\circ = -57200 - 298 \times (-175.8) = -57200 + 52388 = -4812\,\mathrm{J/mol} = -4.81\,\mathrm{kJ/mol}Kp=exp(LB4812RB◆◆LB8.314×298RB)=exp(1.943)=6.98K_p = \exp\left(\frac◆LB◆4812◆RB◆◆LB◆8.314 \times 298◆RB◆\right) = \exp(1.943) = 6.98

(b) Kp=1K_p = 1 when ΔG=0\Delta G^\circ = 0:

T=LBΔHRB◆◆LBΔSRB=57200175.8=325KT = \frac◆LB◆\Delta H^\circ◆RB◆◆LB◆\Delta S^\circ◆RB◆ = \frac{-57200}{-175.8} = 325\,\mathrm{K}

(c) The forward reaction is exothermic (ΔH<0\Delta H^\circ \lt 0) and decreases entropy (ΔS<0\Delta S^\circ \lt 0). Increasing temperature makes ΔG\Delta G^\circ less negative (eventually positive), so KpK_p decreases. The yield of N2O4\mathrm{N}_2\mathrm{O}_4 decreases with increasing temperature. This is consistent with Le Chatelier's principle.

Entropy in Chemical Processes

Entropy of Phase Changes

At a phase transition, ΔG=0\Delta G = 0, so ΔS=ΔH/T\Delta S = \Delta H / T.

For vaporisation: ΔSvap=ΔHvap/Tboil\Delta S_\mathrm{vap} = \Delta H_\mathrm{vap} / T_\mathrm{boil}

Trouton's rule: ΔSvap85Jmol1K1\Delta S_\mathrm{vap} \approx 85\,\mathrm{J\,mol^{-1}\,K^{-1}} for most non-polar liquids. Deviations indicate hydrogen bonding (e.g. water: ΔSvap=109Jmol1K1\Delta S_\mathrm{vap} = 109\,\mathrm{J\,mol^{-1}\,K^{-1}}, due to extra ordering in the liquid from H-bonds).

Entropy of Mixing

When two ideal gases (or two ideal solutions) mix, the entropy always increases:

ΔSmix=nR(x1lnx1+x2lnx2)\Delta S_\mathrm{mix} = -nR(x_1 \ln x_1 + x_2 \ln x_2)

where x1x_1 and x2x_2 are the mole fractions. For equal amounts (x1=x2=0.5x_1 = x_2 = 0.5):

ΔSmix=nR(0.5ln0.5+0.5ln0.5)=nR(0.693)=5.76nJ/K\Delta S_\mathrm{mix} = -nR(0.5 \ln 0.5 + 0.5 \ln 0.5) = -nR(-0.693) = 5.76n\,\mathrm{J/K}

This is the thermodynamic basis for diffusion: gases spontaneously mix because the mixed state has higher entropy.

The Third Law of Thermodynamics

The entropy of a perfect crystal at absolute zero is zero:

S(0K)=0S(0\,\mathrm{K}) = 0

This provides the reference point for absolute entropies (SS^\circ values tabulated in data books). Unlike enthalpy, entropy has an absolute scale.

Worked Example: Calculating ΔS\Delta S^\circ from Absolute Entropies

Calculate ΔS\Delta S^\circ for the combustion of methane:

CH4(g)+2O2(g)CO2(g)+2H2O(l)\mathrm{CH}_4(g) + 2\mathrm{O}_2(g) \to \mathrm{CO}_2(g) + 2\mathrm{H}_2\mathrm{O}(l)

SS^\circ values: CH4(g)=186.3\mathrm{CH}_4(g) = 186.3, O2(g)=205.1\mathrm{O}_2(g) = 205.1, CO2(g)=213.7\mathrm{CO}_2(g) = 213.7, H2O(l)=69.9Jmol1K1\mathrm{H}_2\mathrm{O}(l) = 69.9\,\mathrm{J\,mol^{-1}\,K^{-1}}.

ΔS=[213.7+2(69.9)][186.3+2(205.1)]=353.5596.5=243.0Jmol1K1\Delta S^\circ = [213.7 + 2(69.9)] - [186.3 + 2(205.1)] = 353.5 - 596.5 = -243.0\,\mathrm{J\,mol^{-1}\,K^{-1}}

The entropy decreases because 3 moles of gas produce 1 mole of gas + 2 moles of liquid. The decrease in the number of gaseous molecules dominates.

Gibbs Free Energy in Biological Systems

ATP hydrolysis is the energy currency of cells:

ATP+H2OADP+PiΔG=30.5kJ/mol\mathrm{ATP} + \mathrm{H}_2\mathrm{O} \to \mathrm{ADP} + \mathrm{P}_i \quad \Delta G^\circ = -30.5\,\mathrm{kJ/mol}

The large negative ΔG\Delta G^\circ makes this reaction thermodynamically favourable, and it is coupled to endergonic (unfavourable) reactions in the cell. For example, the synthesis of glucose-6-phosphate from glucose and phosphate (ΔG=+13.8kJ/mol\Delta G^\circ = +13.8\,\mathrm{kJ/mol}) is driven by coupling with ATP hydrolysis:

Glucose+Pi+ATPGlucose-6-phosphate+ADP\mathrm{Glucose} + \mathrm{P}_i + \mathrm{ATP} \to \mathrm{Glucose\text{-}6\text{-}phosphate} + \mathrm{ADP} ΔGoverall=13.8+(30.5)=16.7kJ/mol\Delta G^\circ_\mathrm{overall} = 13.8 + (-30.5) = -16.7\,\mathrm{kJ/mol}

The coupled reaction is spontaneous because the overall ΔG\Delta G^\circ is negative.

Thermodynamic Cycles

Born-Haber Cycles (Recap)

Born-Haber cycles apply Hess's Law to ionic compound formation. They are covered in detail in Born-Haber Cycles.

Enthalpy-Entropy Compensation

Some reactions show enthalpy-entropy compensation: a more exothermic ΔH\Delta H is offset by a more negative ΔS\Delta S, so ΔG\Delta G changes less than expected. This is common in:

  • Solvent reorganisation around dissolved species.
  • Protein folding (hydrophobic effect).
  • Ligand binding.

Worked Examples: Comprehensive Gibbs Free Energy Problems

Problem: Predicting Spontaneity at Different Temperatures

For the reaction C(graphite)+CO2(g)2CO(g)\mathrm{C}(\text{graphite}) + \mathrm{CO}_2(g) \rightleftharpoons 2\mathrm{CO}(g):

ΔH=+173kJ/mol,ΔS=+176Jmol1K1\Delta H^\circ = +173\,\mathrm{kJ/mol}, \quad \Delta S^\circ = +176\,\mathrm{J\,mol^{-1}\,K^{-1}}

(a) Calculate ΔG\Delta G^\circ at 298K298\,\mathrm{K} and 1000K1000\,\mathrm{K}.

(b) At what temperature does the reaction become spontaneous?

(c) Calculate KK at 1000K1000\,\mathrm{K}.

(a) At 298K298\,\mathrm{K}:

ΔG=173×103298×176=17300052448=+120552J/mol=+120.6kJ/mol\Delta G^\circ = 173 \times 10^3 - 298 \times 176 = 173000 - 52448 = +120\,552\,\mathrm{J/mol} = +120.6\,\mathrm{kJ/mol}

Not spontaneous at room temperature.

At 1000K1000\,\mathrm{K}:

ΔG=173×1031000×176=173000176000=3000J/mol=3.0kJ/mol\Delta G^\circ = 173 \times 10^3 - 1000 \times 176 = 173000 - 176000 = -3000\,\mathrm{J/mol} = -3.0\,\mathrm{kJ/mol}

Spontaneous at 1000K1000\,\mathrm{K}.

(b) ΔG=0\Delta G^\circ = 0 when:

T=LBΔHRB◆◆LBΔSRB=LB173×103RB◆◆LB176RB=983KT = \frac◆LB◆\Delta H^\circ◆RB◆◆LB◆\Delta S^\circ◆RB◆ = \frac◆LB◆173 \times 10^3◆RB◆◆LB◆176◆RB◆ = 983\,\mathrm{K}

(c) At 1000K1000\,\mathrm{K}:

K=exp(LBΔGRB◆◆LBRTRB)=exp(LB3000RB◆◆LB8.314×1000RB)=exp(0.361)=1.43K = \exp\left(\frac◆LB◆-\Delta G^\circ◆RB◆◆LB◆RT◆RB◆\right) = \exp\left(\frac◆LB◆3000◆RB◆◆LB◆8.314 \times 1000◆RB◆\right) = \exp(0.361) = 1.43

K>1K > 1, confirming that products are favoured.

Problem: Using Gibbs Energy to Predict Decomposition

Will Ag2CO3(s)\mathrm{Ag}_2\mathrm{CO}_3(s) decompose at 500K500\,\mathrm{K}?

Ag2CO3(s)Ag2O(s)+CO2(g)\mathrm{Ag}_2\mathrm{CO}_3(s) \to \mathrm{Ag}_2\mathrm{O}(s) + \mathrm{CO}_2(g)

ΔH=+82kJ/mol\Delta H^\circ = +82\,\mathrm{kJ/mol}, ΔS=+170Jmol1K1\Delta S^\circ = +170\,\mathrm{J\,mol^{-1}\,K^{-1}}.

ΔG=82000500×170=8200085000=3000J/mol\Delta G^\circ = 82000 - 500 \times 170 = 82000 - 85000 = -3000\,\mathrm{J/mol}

ΔG<0\Delta G^\circ < 0, so the decomposition is spontaneous at 500K500\,\mathrm{K}. (The threshold temperature is T=82000/170=482KT = 82000/170 = 482\,\mathrm{K}.)

Advanced Entropy and Gibbs Energy

Entropy Changes of Mixing

When two ideal gases mix, the entropy always increases because there are more ways to arrange the molecules in the larger volume.

Worked Example: 1.0mol1.0\,\mathrm{mol} of He\mathrm{He} and 1.0mol1.0\,\mathrm{mol} of Ne\mathrm{Ne}, both initially in separate 10.0dm310.0\,\mathrm{dm}^3 containers at 298K298\,\mathrm{K}, are allowed to mix in a combined volume of 20.0dm320.0\,\mathrm{dm}^3. Calculate ΔSmix\Delta S_\text{mix}.

ΔSmix=nR(xHelnxHe+xNelnxNe)\Delta S_\text{mix} = -nR\left(x_\mathrm{He}\ln x_\mathrm{He} + x_\mathrm{Ne}\ln x_\mathrm{Ne}\right)

where xHe=xNe=0.5x_\mathrm{He} = x_\mathrm{Ne} = 0.5:

ΔSmix=(1.0+1.0)×8.314×(0.5ln0.5+0.5ln0.5)\Delta S_\text{mix} = -(1.0 + 1.0) \times 8.314 \times (0.5 \ln 0.5 + 0.5 \ln 0.5)

=2.0×8.314×ln0.5=16.63×(0.693)=+11.5JK1= -2.0 \times 8.314 \times \ln 0.5 = -16.63 \times (-0.693) = +11.5\,\mathrm{J\,K^{-1}}

Phase Transitions and Entropy

At a phase transition, the system is at equilibrium so ΔG=0\Delta G = 0, giving:

ΔStransition=LBΔHtransitionRB◆◆LBTtransitionRB\Delta S_\text{transition} = \frac◆LB◆\Delta H_\text{transition}◆RB◆◆LB◆T_\text{transition}◆RB◆

TransitionΔH\Delta HΔS\Delta SΔG\Delta G
Melting (fusion)Positive (endothermic)Positive (disorder increases)0 (at TmT_\text{m})
Boiling (vaporisation)Positive (endothermic)Positive (large increase in disorder)0 (at TbT_\text{b})
FreezingNegative (exothermic)Negative (order increases)0 (at TmT_\text{m})
SublimationPositive (endothermic)Positive0 (at TsubT_\text{sub})

Worked Example: Calculate the entropy of vaporisation of water at 373K373\,\mathrm{K} given ΔHvap=+40.7kJ/mol\Delta H_\text{vap} = +40.7\,\mathrm{kJ/mol}.

ΔSvap=40700373=+109Jmol1K1\Delta S_\text{vap} = \frac{40700}{373} = +109\,\mathrm{J\,mol^{-1}\,K^{-1}}

This is close to Trouton's rule (ΔSvap88Jmol1K1\Delta S_\text{vap} \approx 88\,\mathrm{J\,mol^{-1}\,K^{-1}} for non-hydrogen-bonding liquids). Water is higher because of extensive hydrogen bonding in the liquid phase.

Born-Haber Cycles and Gibbs Energy

Gibbs energy of formation can be calculated from Born-Haber cycles by using ΔG=ΔHTΔS\Delta G = \Delta H - T\Delta S for each step.

Worked Example: Calculate ΔGf\Delta G_f^\circ for NaCl(s)\mathrm{NaCl}(s) at 298K298\,\mathrm{K}.

Using the Born-Haber cycle values:

  • ΔHf(NaCl)=411kJ/mol\Delta H_f^\circ(\mathrm{NaCl}) = -411\,\mathrm{kJ/mol}
  • ΔSf=72.1Jmol1K1\Delta S_f^\circ = -72.1\,\mathrm{J\,mol^{-1}\,K^{-1}} (system becomes more ordered: solid from gas atoms)

ΔGf=411000298×(72.1)=411000+21486=389514J/mol=389.5kJ/mol\Delta G_f^\circ = -411000 - 298 \times (-72.1) = -411000 + 21486 = -389514\,\mathrm{J/mol} = -389.5\,\mathrm{kJ/mol}

Coupled Reactions in Biochemistry

A thermodynamically unfavourable reaction (ΔG>0\Delta G > 0) can be driven by coupling it to a thermodynamically favourable one (ΔG<0\Delta G < 0), provided the overall ΔG<0\Delta G < 0.

Example: Hydrolysis of ATP: ATP+H2OADP+PiΔG=30.5kJ/mol\mathrm{ATP} + \mathrm{H}_2\mathrm{O} \to \mathrm{ADP} + \mathrm{P}_i \quad \Delta G^\circ = -30.5\,\mathrm{kJ/mol}

This strongly exergonic reaction drives many endergonic processes in cells. If a reaction requires +20kJ/mol+20\,\mathrm{kJ/mol}, coupling with ATP hydrolysis gives:

ΔGoverall=+20.0+(30.5)=10.5kJ/mol\Delta G_\text{overall} = +20.0 + (-30.5) = -10.5\,\mathrm{kJ/mol}

Gibbs Energy and Equilibrium: Quantitative Treatment

The relationship between ΔG\Delta G, the reaction quotient QQ, and the equilibrium constant KK:

ΔG=ΔG+RTlnQ\Delta G = \Delta G^\circ + RT\ln Q

At equilibrium, ΔG=0\Delta G = 0 and Q=KQ = K, giving:

ΔG=RTlnK\Delta G^\circ = -RT\ln K

Worked Example: For the reaction N2O4(g)2NO2(g)\mathrm{N}_2\mathrm{O}_4(g) \rightleftharpoons 2\mathrm{NO}_2(g) at 298K298\,\mathrm{K}:

ΔG=+4.72kJ/mol\Delta G^\circ = +4.72\,\mathrm{kJ/mol}. Calculate KK.

K=exp(LBΔGRB◆◆LBRTRB)=exp(LB4720RB◆◆LB8.314×298RB)=exp(1.905)=0.149K = \exp\left(\frac◆LB◆-\Delta G^\circ◆RB◆◆LB◆RT◆RB◆\right) = \exp\left(\frac◆LB◆-4720◆RB◆◆LB◆8.314 \times 298◆RB◆\right) = \exp(-1.905) = 0.149

If the initial pressure of N2O4\mathrm{N}_2\mathrm{O}_4 is 1.00atm1.00\,\mathrm{atm} and no NO2\mathrm{NO}_2 is present:

Q=LB(pNO2)2RB◆◆LBpN2O4RB=01=0Q = \frac◆LB◆(p_{\mathrm{NO}_2})^2◆RB◆◆LB◆p_{\mathrm{N}_2\mathrm{O}_4}◆RB◆ = \frac{0}{1} = 0

ΔG=4720+8.314×298×ln0=4720+()\Delta G = 4720 + 8.314 \times 298 \times \ln 0 = 4720 + (-\infty) \to -\infty

ΔG\Delta G is very negative, so the forward reaction is strongly favoured initially (the reaction proceeds to the right until equilibrium is reached).

Common Pitfalls

  1. Sign errors in ΔG=ΔHTΔS\Delta G = \Delta H - T\Delta S: Remember the minus sign. A positive ΔH\Delta H and positive ΔS\Delta S means the reaction is spontaneous at high TT (the TΔS-T\Delta S term dominates). Students often incorrectly write ΔG=ΔH+TΔS\Delta G = \Delta H + T\Delta S.

  2. Units of ΔS\Delta S: Always use Jmol1K1\mathrm{J\,mol^{-1}\,K^{-1}} for entropy and kJ/mol\mathrm{kJ/mol} for enthalpy. You must convert one of them before combining. Forgetting to convert ΔS\Delta S from J\mathrm{J} to kJ\mathrm{kJ} (divide by 1000) is the single most common arithmetic error.

  3. Standard vs non-standard conditions: ΔG\Delta G^\circ applies only when all components are in their standard states (1 mol/dm3^3 for solutions, 100kPa100\,\mathrm{kPa} for gases, pure solids/liquids). Under non-standard conditions, use ΔG=ΔG+RTlnQ\Delta G = \Delta G^\circ + RT\ln Q.

  4. Assuming ΔG<0\Delta G^\circ < 0 means the reaction happens quickly: Thermodynamic feasibility does not imply kinetic feasibility. Diamond converting to graphite has ΔG<0\Delta G^\circ < 0 but the rate is essentially zero at room temperature.

  5. Entropy of a pure element: The absolute entropy of a pure element in its standard state at 298K298\,\mathrm{K} is not zero (only SS^\circ at 0K0\,\mathrm{K} is zero, by the third law).

Practical Applications: Gibbs Energy in Industry

Industrial ProcessReactionΔH\Delta H^\circ (kJ/mol)ΔS\Delta S^\circ (J/mol/K)TthresholdT_\text{threshold} (K)
Haber processN2+3H22NH3\mathrm{N}_2 + 3\mathrm{H}_2 \to 2\mathrm{NH}_392-92199-199Not applicable (ΔH<0\Delta H < 0, ΔS<0\Delta S < 0: spontaneous at low TT)
Contact process2SO2+O22SO32\mathrm{SO}_2 + \mathrm{O}_2 \to 2\mathrm{SO}_3198-198190-190Not applicable
Thermal decomposition of CaCO3\mathrm{CaCO}_3CaCO3CaO+CO2\mathrm{CaCO}_3 \to \mathrm{CaO} + \mathrm{CO}_2+178+178+161+16111061106
Roasting of ZnS\mathrm{ZnS}2ZnS+3O22ZnO+2SO22\mathrm{ZnS} + 3\mathrm{O}_2 \to 2\mathrm{ZnO} + 2\mathrm{SO}_2880-880+100+100Spontaneous at all TT

Exam-Style Questions with Full Mark Schemes

Q1 (5 marks)

For the reaction C(s)+H2O(g)CO(g)+H2(g)\mathrm{C}(s) + \mathrm{H}_2\mathrm{O}(g) \to \mathrm{CO}(g) + \mathrm{H}_2(g):

ΔH=+131kJ/mol\Delta H^\circ = +131\,\mathrm{kJ/mol}, ΔS=+134Jmol1K1\Delta S^\circ = +134\,\mathrm{J\,mol^{-1}\,K^{-1}}.

(a) Calculate ΔG\Delta G^\circ at 298K298\,\mathrm{K} and state whether the reaction is feasible. (3 marks)

(b) Calculate the minimum temperature at which the reaction becomes feasible. (2 marks)

Mark Scheme:

(a) ΔG=131000298×134=13100039932=+91068J/mol=+91.1kJ/mol\Delta G^\circ = 131000 - 298 \times 134 = 131000 - 39932 = +91068\,\mathrm{J/mol} = +91.1\,\mathrm{kJ/mol} (2 marks for calculation).

ΔG>0\Delta G^\circ > 0, so the reaction is not feasible at 298K298\,\mathrm{K} (1 mark).

(b) ΔG=0\Delta G^\circ = 0 when T=ΔH/ΔS=131000/134=978KT = \Delta H^\circ / \Delta S^\circ = 131000 / 134 = 978\,\mathrm{K} (2 marks).

The reaction becomes feasible above 978K978\,\mathrm{K}.

Q2 (6 marks)

Explain why the entropy change for the reaction N2O4(g)2NO2(g)\mathrm{N}_2\mathrm{O}_4(g) \rightleftharpoons 2\mathrm{NO}_2(g) is positive, and calculate ΔS\Delta S^\circ given the following standard entropies:

S(N2O4)=304Jmol1K1S^\circ(\mathrm{N}_2\mathrm{O}_4) = 304\,\mathrm{J\,mol^{-1}\,K^{-1}}, S(NO2)=240Jmol1K1S^\circ(\mathrm{NO}_2) = 240\,\mathrm{J\,mol^{-1}\,K^{-1}}.

Mark Scheme:

The entropy change is positive because one mole of gas produces two moles of gas (1 mark). There are more ways to arrange the molecules and more microstates when there are more gas particles (1 mark). The products have greater positional disorder than the reactants (1 mark).

ΔS=2×240304=480304=+176Jmol1K1\Delta S^\circ = 2 \times 240 - 304 = 480 - 304 = +176\,\mathrm{J\,mol^{-1}\,K^{-1}} (3 marks for calculation with correct units).

Q3 (5 marks)

A student claims that because the combustion of methane is highly exothermic (ΔHc=890kJ/mol\Delta H_c^\circ = -890\,\mathrm{kJ/mol}), it must be thermodynamically feasible at all temperatures. Evaluate this claim.

Mark Scheme:

CH4(g)+2O2(g)CO2(g)+2H2O(l)\mathrm{CH}_4(g) + 2\mathrm{O}_2(g) \to \mathrm{CO}_2(g) + 2\mathrm{H}_2\mathrm{O}(l)

ΔS=[213+2(70)][186+2(205)]=353596=243Jmol1K1\Delta S^\circ = [213 + 2(70)] - [186 + 2(205)] = 353 - 596 = -243\,\mathrm{J\,mol^{-1}\,K^{-1}} (2 marks).

Since ΔH<0\Delta H < 0 and ΔS<0\Delta S < 0, the reaction is feasible only when ΔG=ΔHTΔS<0\Delta G = \Delta H - T\Delta S < 0, i.e. when T<ΔH/ΔS=890000/243=3663KT < |\Delta H / \Delta S| = 890000/243 = 3663\,\mathrm{K} (2 marks).

The claim is correct in practice (combustion is feasible at all reasonable temperatures), but incorrect in principle -- at sufficiently high temperatures (above 3663K3663\,\mathrm{K}), the reaction would not be thermodynamically feasible (1 mark).

Q4 (4 marks)

The melting point of sodium is 371K371\,\mathrm{K} and ΔHfus=+2.60kJ/mol\Delta H_\text{fus} = +2.60\,\mathrm{kJ/mol}. Calculate the entropy change of fusion and explain its sign.

Mark Scheme:

ΔSfus=LBΔHfusRB◆◆LBTmRB=2600371=+7.01Jmol1K1\Delta S_\text{fus} = \frac◆LB◆\Delta H_\text{fus}◆RB◆◆LB◆T_\text{m}◆RB◆ = \frac{2600}{371} = +7.01\,\mathrm{J\,mol^{-1}\,K^{-1}} (2 marks).

The entropy change is positive because the solid sodium becomes a liquid, which has greater disorder and more ways to arrange the particles (1 mark). The ions in the liquid are no longer fixed in a lattice and have greater freedom of movement (1 mark).