STEMviz
Chemistry · General, Physical and Quantum Chemistry

Thermodynamics of Chemical Reactions and Equations of State

The first law, reversible work for ideal and van der Waals gases, enthalpy and heat capacities, internal pressure and the Joule-Thomson effect, the second law and Gibbs energy, and Kirchhoff's law for reaction enthalpies.

Thermodynamics of Chemical Reactions and Equations of State

Chapter 1 built the equations of state that describe how much volume a gas occupies at given pressure and temperature. This chapter asks a different question: how much energy moves between a chemical system and its surroundings as that state changes, and how is that energy accounted for exactly? Classical thermodynamics answers this with two state functions — internal energy UU and enthalpy HH — and two conservation/direction laws that govern them. The equations of state from Chapter 1 are not left behind: they reappear here as the engine that lets the abstract thermodynamic machinery produce real numbers, from the work done by an expanding van der Waals gas to the cooling of a real gas forced through a throttle.

You should have to hand the van der Waals equation, the critical constants, and the Boyle temperature from Chapter 1. This chapter develops the first law, enthalpy and heat capacities, the internal-pressure concept that measures real-gas non-ideality energetically, the Joule–Thomson effect, a compact but rigorous introduction to the second law and Gibbs energy, and Kirchhoff's law for the temperature dependence of reaction enthalpies. The next chapter uses the Gibbs energy introduced here to define the chemical potential and build the theory of solutions.

1. The First Law of Thermodynamics

The internal energy UU of a system is the total energy contained in it — translational, rotational, vibrational and potential energy of its constituent particles — and it is a state function: its value depends only on the current state of the system, not on the path by which that state was reached. The first law of thermodynamics asserts that UU changes only through heat and work crossing the system boundary:

dU=dq+dwdU = dq + dw

using the IUPAC convention in which qq and ww are both positive when energy flows into the system. Unlike UU, neither qq nor ww is a state function individually — both depend on the path taken between two states — but their sum, dUdU, does not. This is the content of the first law: energy is conserved, and internal energy bookkeeping does not care how the energy arrived.

Key distinction. UU is a state function (exact differential dUdU); qq and ww are path functions (inexact differentials, often written δq,δw\delta q, \delta w). Only sums and differences of state functions between two fixed states are path-independent.

2. Reversible and Irreversible Work of Expansion

For a system doing only expansion work against an external pressure pextp_{\text{ext}}, the work done on the system as its volume changes by dVdV is

dw=−pext dVdw = -p_{\text{ext}}\,dV

the minus sign reflecting that the system loses energy as it expands (dV>0dV>0) against an opposing external pressure. Two limiting cases matter:

  • Free expansion (pext=0p_{\text{ext}}=0, e.g. expansion into a vacuum): w=0w=0 regardless of ΔV\Delta V.
  • Reversible expansion: pextp_{\text{ext}} is adjusted at every instant to equal the system's own pressure pp (quasi-statically, through a continuous sequence of near-equilibrium states). This extracts the maximum possible work from an expansion, since pext=pp_{\text{ext}}=p is the largest external pressure the system can push against without the process becoming impossible (pext>pp_{\text{ext}}>p would compress it instead).

For a reversible isothermal expansion of a perfect gas from V1V_1 to V2V_2 at temperature TT, substitute p=nRT/Vp=nRT/V and integrate:

w=−∫V1V2nRTV dV=−nRTln⁡V2V1w = -\int_{V_1}^{V_2} \frac{nRT}{V}\,dV = -nRT\ln\frac{V_2}{V_1}

For a van der Waals gas, substitute instead p=nRTV−nb−an2V2p = \dfrac{nRT}{V-nb} - \dfrac{an^2}{V^2} (Chapter 1) and integrate term by term:

w=−nRTln⁡V2−nbV1−nb+an2(1V1−1V2)w = -nRT\ln\frac{V_2-nb}{V_1-nb} + an^2\left(\frac{1}{V_1}-\frac{1}{V_2}\right)

Worked example. One mole of CO2_2 (a=3.610 dm6 atm mol−2a=3.610~\text{dm}^6\,\text{atm}\,\text{mol}^{-2}, b=0.0429 dm3 mol−1b=0.0429~\text{dm}^3\,\text{mol}^{-1}) expands reversibly and isothermally at 300 K300~\text{K} from V1=0.500 dm3V_1=0.500~\text{dm}^3 to V2=5.00 dm3V_2=5.00~\text{dm}^3. The perfect-gas prediction is

wideal=−(1.00)(0.08206)(300)ln⁡5.000.500=−24.618×2.3026=−56.69 dm3 atm=−5.74 kJw_{\text{ideal}} = -(1.00)(0.08206)(300)\ln\frac{5.00}{0.500} = -24.618 \times 2.3026 = -56.69~\text{dm}^3\,\text{atm} = -5.74~\text{kJ}

The van der Waals result uses V1−nb=0.4571V_1-nb=0.4571, V2−nb=4.9571V_2-nb=4.9571:

wvdW=−24.618ln⁡4.95710.4571+3.610(10.500−15.00)=−58.68+6.50=−52.18 dm3 atm=−5.29 kJw_{\text{vdW}} = -24.618\ln\frac{4.9571}{0.4571} + 3.610\left(\frac{1}{0.500}-\frac{1}{5.00}\right) = -58.68 + 6.50 = -52.18~\text{dm}^3\,\text{atm} = -5.29~\text{kJ}

The van der Waals gas does less work expanding (5.29 kJ5.29~\text{kJ} against 5.74 kJ5.74~\text{kJ}) than the perfect-gas model predicts — a direct energetic consequence of the lower pressure that intermolecular attraction produces at every volume along the path, exactly the effect Chapter 1 quantified.

3. Enthalpy and Heat Capacities

Enthalpy is defined as H≡U+pVH \equiv U + pV, a state function (being a combination of state functions) chosen so that, at constant pressure, dH=dqpdH = dq_p: enthalpy change equals heat absorbed at constant pressure, the condition under which most chemistry is actually done (open vessels, atmospheric pressure).

The heat capacity at constant volume, CV=(∂U∂T)VC_V = \left(\dfrac{\partial U}{\partial T}\right)_V, and the heat capacity at constant pressure, Cp=(∂H∂T)pC_p = \left(\dfrac{\partial H}{\partial T}\right)_p, measure how much heat is needed per unit temperature rise under each constraint. For a perfect gas, H=U+nRTH = U + nRT, so

(∂H∂T)p=(∂U∂T)p+nR\left(\frac{\partial H}{\partial T}\right)_p = \left(\frac{\partial U}{\partial T}\right)_p + nR

Joule's law (justified in the next section) states that for a perfect gas UU depends on TT alone, so (∂U/∂T)p=(∂U/∂T)V=CV(\partial U/\partial T)_p = (\partial U/\partial T)_V = C_V, giving Mayer's relation:

Cp−CV=nR(perfect gas only)C_p - C_V = nR \qquad \text{(perfect gas only)}

4. Internal Pressure and the Thermodynamic Equation of State

The internal pressure, πT≡(∂U∂V)T\pi_T \equiv \left(\dfrac{\partial U}{\partial V}\right)_T, measures how internal energy changes with volume at constant temperature — physically, the energy cost of pulling molecules apart against their mutual attraction (or pushing them together against repulsion), separate from any change in their kinetic energy. Joule attempted to measure this by letting a gas expand freely (adiabatically, into a vacuum, so q=0q=0 and w=0w=0, hence ΔU=0\Delta U=0) and looking for a temperature change; within the precision available to him, none was found for a near-ideal gas, establishing Joule's law: πT=0\pi_T = 0 for a perfect gas, i.e. U=U(T)U = U(T) only.

A general thermodynamic identity — the thermodynamic equation of state, derivable by combining the first and second laws through a Maxwell relation (the result is quoted here without proof) — expresses the internal pressure in terms of the equation of state alone:

πT=T(∂p∂T)V−p\pi_T = T\left(\frac{\partial p}{\partial T}\right)_V - p

For a perfect gas, p=nRT/Vp=nRT/V gives (∂p∂T)V=nR/V\left(\dfrac{\partial p}{\partial T}\right)_V = nR/V, so T(∂p/∂T)V=nRT/V=pT(\partial p/\partial T)_V = nRT/V = p, and πT=p−p=0\pi_T = p-p=0 — recovering Joule's law exactly.

For a van der Waals gas, p=nRTV−nb−an2V2p = \dfrac{nRT}{V-nb}-\dfrac{an^2}{V^2} gives (∂p∂T)V=nRV−nb\left(\dfrac{\partial p}{\partial T}\right)_V = \dfrac{nR}{V-nb}, so

T(∂p∂T)V=nRTV−nb=p+an2V2T\left(\frac{\partial p}{\partial T}\right)_V = \frac{nRT}{V-nb} = p + \frac{an^2}{V^2}

(the last step uses the van der Waals equation itself to replace nRT/(V−nb)nRT/(V-nb)). Therefore

πT=an2V2(per mole: πT=a/Vm2)\boxed{\pi_T = \frac{an^2}{V^2}} \qquad \text{(per mole: } \pi_T = a/V_m^2\text{)}

Internal energy does depend on volume for a van der Waals gas, growing as the gas expands — exactly as expected, since expanding against attractive forces (without any compensating change in kinetic energy, TT fixed) stores energy as increased separation of attracting molecules.

5. The General Relation Between CpC_p and CVC_V

For any substance, not just a perfect gas, a parallel thermodynamic argument (again using the first and second laws and a Maxwell relation, quoted here without proof) gives the exact relation

Cp−CV=[p+(∂U∂V)T](∂V∂T)p=(πT+p)(∂V∂T)pC_p - C_V = \left[p + \left(\frac{\partial U}{\partial V}\right)_T\right]\left(\frac{\partial V}{\partial T}\right)_p = \left(\pi_T + p\right)\left(\frac{\partial V}{\partial T}\right)_p

For a perfect gas, πT=0\pi_T=0 and p=nRT/Vp = nRT/V combine with (∂V/∂T)p=nR/p(\partial V/\partial T)_p = nR/p to reproduce Cp−CV=p⋅nR/p=nRC_p-C_V = p\cdot nR/p = nR, consistent with Section 3. For a real gas, πT≠0\pi_T \ne 0 and Cp−CVC_p - C_V departs from nRnR — the correction is largest exactly where the gas is least ideal, linking this chapter's thermodynamics directly back to Chapter 1's equation-of-state corrections.

6. The Joule–Thomson Effect

A more sensitive experiment than Joule's free expansion is the Joule–Thomson (throttling) experiment: a gas is forced continuously through a porous plug or valve from high pressure to low pressure, with the whole apparatus thermally insulated. This process occurs at constant enthalpy (ΔH=0\Delta H = 0 for the throttled gas — a standard result of the steady-flow energy balance), and the temperature change observed defines the Joule–Thomson coefficient:

μJT≡(∂T∂p)H\mu_{\text{JT}} \equiv \left(\frac{\partial T}{\partial p}\right)_H

A perfect gas has μJT=0\mu_{\text{JT}}=0 identically. A real gas does not, because H=H(T,p)H=H(T,p) depends on pp through intermolecular forces; using (∂H/∂p)T=Vm−T(∂Vm/∂T)p(\partial H/\partial p)_T = V_m - T(\partial V_m/\partial T)_p (the enthalpy analogue of the internal-pressure identity of Section 4) together with the low-pressure expansion of the van der Waals equation, Vm≈RT/p+B(T)V_m \approx RT/p + B(T) with B(T)=b−a/RTB(T) = b - a/RT (Chapter 1), gives after differentiation and simplification:

μJT≈1Cp,m(2aRT−b)\mu_{\text{JT}} \approx \frac{1}{C_{p,m}}\left(\frac{2a}{RT} - b\right)

This changes sign at the Joule–Thomson inversion temperature, Tinv=2aRbT_{\text{inv}} = \dfrac{2a}{Rb} — exactly twice the Boyle temperature TB=a/RbT_B = a/Rb found in Chapter 1. Below TinvT_{\text{inv}}, μJT>0\mu_{\text{JT}}>0: throttling cools the gas (attraction dominates — expanding against attraction absorbs energy, which must come from the gas's own kinetic energy, lowering TT). Above TinvT_{\text{inv}}, μJT<0\mu_{\text{JT}}<0: throttling heats the gas (the excluded-volume effect dominates).

Worked example. For N2_2 (a=1.352 dm6 atm mol−2a=1.352~\text{dm}^6\,\text{atm}\,\text{mol}^{-2}, b=0.0387 dm3 mol−1b=0.0387~\text{dm}^3\,\text{mol}^{-1}, Cp,m≈29.1 J mol−1K−1=0.2872 dm3 atm mol−1K−1C_{p,m}\approx29.1~\text{J}\,\text{mol}^{-1}\text{K}^{-1} = 0.2872~\text{dm}^3\,\text{atm}\,\text{mol}^{-1}\text{K}^{-1}) at T=300 KT=300~\text{K}:

2aRT=2(1.352)(0.08206)(300)=2.70424.618=0.1098 dm3 mol−1\frac{2a}{RT} = \frac{2(1.352)}{(0.08206)(300)} = \frac{2.704}{24.618} = 0.1098~\text{dm}^3\,\text{mol}^{-1} μJT=0.1098−0.03870.2872=0.07110.2872=0.248 Katm\mu_{\text{JT}} = \frac{0.1098-0.0387}{0.2872} = \frac{0.0711}{0.2872} = 0.248~\dfrac{\text{K}}{\text{atm}}

a positive coefficient, meaning N2_2 cools on throttling at room temperature — consistent with its measured value and with T=300 KT=300~\text{K} being far below its inversion temperature, Tinv=2a/Rb=2(425.7 K)=851.4 KT_{\text{inv}} = 2a/Rb = 2(425.7~\text{K}) = 851.4~\text{K} (twice the Boyle temperature computed in Chapter 1).

7. The Second Law, Entropy and Gibbs Energy

The first law forbids energy from being created or destroyed but says nothing about which processes actually occur spontaneously — a gas never spontaneously compresses itself into one corner of a room, even though doing so would violate no energy conservation law. The second law of thermodynamics supplies the missing direction, through a new state function, entropy SS, defined via dS=dqrevTdS = \dfrac{dq_{\text{rev}}}{T} for a reversible path between two states (entropy, like UU and HH, is path-independent even though qq is not, because the reversible-path definition is shown to be unique to the two end states). The Clausius inequality, in its original system-level form, states that for any real (possibly irreversible) process exchanging heat dqdq with its surroundings at temperature TT,

dSsystem≥dqTdS_{\text{system}} \ge \frac{dq}{T}

with equality only in the reversible limit. Adding the corresponding entropy change of the surroundings, dSsurroundings=−dq/TdS_{\text{surroundings}} = -dq/T for the actual heat transferred, gives the equivalent universe-level statement used throughout this chapter,

dSuniverse=dSsystem+dSsurroundings≥0dS_{\text{universe}} = dS_{\text{system}} + dS_{\text{surroundings}} \ge 0

again with equality only for a reversible process — spontaneous processes always increase the entropy of the universe. For a reversible isothermal expansion of a perfect gas, dqrev=−dw=nRT dV/Vdq_{\text{rev}} = -dw = nRT\,dV/V, so ΔS=nRln⁡(V2/V1)\Delta S = nR\ln(V_2/V_1), positive for an expansion, consistent with the intuitive picture of entropy as a measure of the number of microstates accessible to a system (more volume, more accessible positions, higher entropy).

The second-law criterion dSuniverse≥0dS_{\text{universe}}\ge0 is a statement about the universe (system plus surroundings), which is often inconvenient for predicting the behaviour of a chemical system alone at constant temperature and pressure. The Gibbs energy, G≡H−TSG \equiv H - TS, is constructed precisely to repackage the second law as a statement about the system alone: at constant TT and pp, it can be shown that dSuniverse≥0dS_{\text{universe}} \ge 0 is equivalent to

dG≤0(constant T,p)dG \le 0 \qquad \text{(constant } T, p\text{)}

A process at constant temperature and pressure is spontaneous if and only if it decreases the Gibbs energy of the system; equilibrium is reached when GG is at a minimum (dG=0dG=0). This single criterion — requiring no bookkeeping of the surroundings at all — is the practical workhorse of chemical thermodynamics, and the next chapter builds the chemical potential, μi=(∂G/∂ni)T,p,nj≠i\mu_i = (\partial G/\partial n_i)_{T,p,n_{j\ne i}}, directly on top of it to describe mixtures and solutions.

8. Standard Reaction Enthalpies and Kirchhoff's Law

The standard reaction enthalpy, ΔrH⊖\Delta_rH^{\ominus}, is the enthalpy change when reactants in their standard states (pure, at 1 bar1~\text{bar}, in their stated physical state) are converted completely to products in their standard states, per mole of reaction as written. Because HH is a state function, ΔrH⊖\Delta_rH^{\ominus} can be built from tabulated standard enthalpies of formation, ΔrH⊖=∑νiΔfH⊖(products)−∑νiΔfH⊖(reactants)\Delta_rH^{\ominus} = \sum \nu_i \Delta_fH^{\ominus}(\text{products}) - \sum \nu_i \Delta_fH^{\ominus}(\text{reactants}) (Hess's law), independent of the reaction path.

Reaction enthalpies are tabulated at a reference temperature (conventionally 298 K298~\text{K}) but chemical processes often occur elsewhere. Since HH is a state function, differentiating ΔrH=∑νiHi(products)−∑νiHi(reactants)\Delta_rH = \sum\nu_iH_i(\text{products}) - \sum\nu_iH_i(\text{reactants}) with respect to TT at constant pp, and using (∂Hi/∂T)p=Cp,i(\partial H_i/\partial T)_p = C_{p,i} for each species, gives Kirchhoff's law:

(∂ΔrH∂T)p=ΔrCp≡∑νiCp,i(products)−∑νiCp,i(reactants)\left(\frac{\partial \Delta_rH}{\partial T}\right)_p = \Delta_rC_p \equiv \sum\nu_iC_{p,i}(\text{products}) - \sum\nu_iC_{p,i}(\text{reactants})

If ΔrCp\Delta_rC_p is approximately constant over the temperature range of interest, this integrates to ΔrH(T2)=ΔrH(T1)+ΔrCp(T2−T1)\Delta_rH(T_2) = \Delta_rH(T_1) + \Delta_rC_p(T_2-T_1).

Worked example. For the Haber process, N2(g)+3H2(g)→2NH3(g)\text{N}_2(g) + 3\text{H}_2(g) \rightarrow 2\text{NH}_3(g), ΔrH⊖(298 K)=−92.2 kJ mol−1\Delta_rH^{\ominus}(298~\text{K}) = -92.2~\text{kJ}\,\text{mol}^{-1}. Using Cp,m(N2)=29.1C_{p,m}(\text{N}_2)=29.1, Cp,m(H2)=28.8C_{p,m}(\text{H}_2)=28.8, Cp,m(NH3)=35.1 J mol−1K−1C_{p,m}(\text{NH}_3)=35.1~\text{J}\,\text{mol}^{-1}\text{K}^{-1}:

ΔrCp=2(35.1)−[29.1+3(28.8)]=70.2−115.5=−45.3 J K−1mol−1\Delta_rC_p = 2(35.1) - \left[29.1 + 3(28.8)\right] = 70.2 - 115.5 = -45.3~\text{J}\,\text{K}^{-1}\text{mol}^{-1} ΔrH(500 K)≈−92,200+(−45.3)(500−298)=−92,200−9,151=−101,351 J mol−1≈−101.4 kJ mol−1\Delta_rH(500~\text{K}) \approx -92{,}200 + (-45.3)(500-298) = -92{,}200 - 9{,}151 = -101{,}351~\text{J}\,\text{mol}^{-1} \approx -101.4~\text{kJ}\,\text{mol}^{-1}

The reaction becomes noticeably more exothermic at higher temperature in this (approximate, constant-ΔCp\Delta C_p) treatment, because the products have a lower combined heat capacity than the reactants.

9. Common Exam Pitfalls

  • Dropping the minus sign in dw=−pextdVdw=-p_{\text{ext}}dV, or confusing the IUPAC sign convention (w>0w>0 when work is done on the system) with the older "engineering" convention. State which convention is in use and stay consistent throughout a calculation.
  • Using pp (the system's own pressure) instead of pextp_{\text{ext}} in an irreversible expansion. They coincide only in the reversible limit; an irreversible expansion against a constant external pressure uses w=−pextΔVw=-p_{\text{ext}}\Delta V directly, no integration needed.
  • Assuming Cp−CV=nRC_p - C_V = nR for a real gas. Mayer's relation is a perfect-gas result; the general relation of Section 5 involves the internal pressure πT\pi_T and reduces to nRnR only when πT=0\pi_T=0.
  • Forgetting that μJT\mu_{\text{JT}} can be positive or negative. Below the inversion temperature throttling cools a gas (used industrially to liquefy gases); above it, throttling heats the gas — the sign is not fixed, and depends on TT relative to Tinv=2a/RbT_{\text{inv}}=2a/Rb.
  • Using dSsystem≥0dS_{\text{system}}\ge0 alone as the spontaneity criterion. The Clausius inequality requires dSuniverse=dSsystem+dSsurroundings≥0dS_{\text{universe}}=dS_{\text{system}}+dS_{\text{surroundings}}\ge0; the Gibbs energy criterion dG≤0dG\le0 is what packages this into a system-only statement, but only at constant TT and pp.
  • Applying Kirchhoff's law with reactant/product CpC_p values swapped or with the wrong stoichiometric coefficients. ΔrCp\Delta_rC_p follows the same "products minus reactants, each weighted by its stoichiometric coefficient" pattern as ΔrH\Delta_rH itself.

Summary

  • The first law, dU=dq+dwdU=dq+dw, makes UU a state function; dw=−pextdVdw=-p_{\text{ext}}dV for expansion work, maximised in magnitude for a reversible process where pext=pp_{\text{ext}}=p at every instant.
  • Reversible isothermal work: w=−nRTln⁡(V2/V1)w=-nRT\ln(V_2/V_1) for a perfect gas; w=−nRTln⁡ ⁣[V2−nbV1−nb]+an2 ⁣(1V1−1V2)w=-nRT\ln\!\left[\frac{V_2-nb}{V_1-nb}\right]+an^2\!\left(\frac1{V_1}-\frac1{V_2}\right) for a van der Waals gas.
  • Enthalpy H=U+pVH=U+pV gives dH=dqpdH=dq_p; Mayer's relation Cp−CV=nRC_p-C_V=nR holds for a perfect gas only.
  • Internal pressure πT=(∂U/∂V)T=T(∂p/∂T)V−p\pi_T=(\partial U/\partial V)_T = T(\partial p/\partial T)_V - p is zero for a perfect gas (Joule's law) and equals a/Vm2a/V_m^2 for a van der Waals gas.
  • The general relation Cp−CV=(πT+p)(∂V/∂T)pC_p-C_V=(\pi_T+p)(\partial V/\partial T)_p reduces to Mayer's relation only when πT=0\pi_T=0.
  • The Joule–Thomson coefficient μJT=(∂T/∂p)H≈(2a/RT−b)/Cp,m\mu_{\text{JT}}=(\partial T/\partial p)_H \approx (2a/RT-b)/C_{p,m} changes sign at the inversion temperature Tinv=2a/Rb=2TBT_{\text{inv}}=2a/Rb=2T_B.
  • The second law, dSuniverse≥0dS_{\text{universe}}\ge0, and the Gibbs energy, G=H−TSG=H-TS with spontaneity criterion dG≤0dG\le0 at constant T,pT,p, set up the chemical potential of the next chapter.
  • Kirchhoff's law, (∂ΔrH/∂T)p=ΔrCp(\partial\Delta_rH/\partial T)_p=\Delta_rC_p, gives the temperature dependence of reaction enthalpies from heat-capacity data.