Chemical Thermodynamics

Chemical thermodynamics is the branch of physical chemistry and macroscopic physics that applies the four laws of thermodynamics and statistical mechanics to analyze energy transformations, phase transitions, directionality (spontaneity), and equilibrium compositions in chemical reactions and multicomponent mixtures. Formulated in its modern mathematical architecture by Josiah Willard Gibbs between 1873 and 1878 (On the Equilibrium of Heterogeneous Substances), and expanded by Hermann von Helmholtz, Jacobus Henricus van 't Hoff, Walther Nernst, and Gilbert N. Lewis, chemical thermodynamics connects calorimetric heat changes (enthalpy, H) and configurational disorder (entropy, S) to the Gibbs free energy (G) and chemical potential (μ_i).

Within chemistry, chemical engineering, geochemistry, and cellular bioenergetics (such as ATP-coupled metabolism), chemical thermodynamics answers two foundational questions about any proposed transformation of chemical bonds: (1) whether a reaction can proceed spontaneously in the forward direction under specified temperature T and pressure P (Δ_rG < 0), and (2) what the final equilibrium ratio of products to reactants (K_eq) will be once net macroscopic reaction ceases (Δ_rG = 0).

Internal Energy (U), Enthalpy (H), and Entropy (S)

A thermodynamic state function is a macroscopic property (U, H, S, A, G, P, V, T) whose change Δ*X* = X_final − X_initial depends only on the initial and final equilibrium states of the system, independent of the path taken, whereas heat (q) and work (w) are path-dependent inexact differentials (δq, δw). The foundation of chemical thermodynamics rests on three core laws (plus the Zeroth Law defining temperature T):

1. First Law (Conservation of Energy): The change in internal energy (U)—the total microscopic thermal kinetic energy and potential energy of all molecules in a closed system—equals the heat absorbed by the system plus the work done on the system: dU = δq + δw = δq − P_ext dV + δw_non-PV. Because laboratory and biological reactions typically occur at constant atmospheric pressure (P) rather than constant volume (V), defining the enthalpy (H = U + PV) ensures that the enthalpy change at constant pressure with zero non-expansion work equals the heat transferred: Δ*H* = q_P (exothermic for Δ*H* < 0; endothermic for Δ*H* > 0, governed by Hess's Law of constant heat summation).
2. Second Law (Clausius Inequality & Boltzmann Entropy): Every spontaneous process increases the total entropy of the universe (Δ*S*_univ = Δ*S*_sys + Δ*S*_surr ≥ 0), where the differential of entropy (S) for a reversible transfer of heat is dS = δq_rev / T, connected microscopically to the number of accessible quantum microstates W via Ludwig Boltzmann's statistical formula S = k_B ln W (and canonical partition function Q in molecular dynamics).
3. Third Law (Nernst–Planck Theorem, 1906–1912): The entropy of any pure, defect-free crystalline substance approaches zero as absolute temperature approaches absolute zero (lim_{*T*→0 K} S = 0), allowing calorimetric determination of standard third-law molar entropies (*S*°₂₉₈ > 0) for all elements and compounds.

Legendre Transforms: Helmholtz (A) and Gibbs (G) Free Energies

To evaluate spontaneity (dS_univ ≥ 0) using system variables alone without tracking the external surroundings, Helmholtz and Gibbs performed Legendre transforms on U(S, V, {*n*_i}) to replace natural variables (S, V) with experimentally controlled intensive variables (T, P):

• Helmholtz Free Energy (A or F = U − TS, natural variables T, V, {*n*_i}): At constant T and V, spontaneity requires dA_T,V ≤ 0, and −Δ*A* equals the maximum total reversible work extractable from the system.
• Gibbs Free Energy (G = H − TS = U + PV − TS, natural variables T, P, {*n*_i}): At constant temperature T and pressure P, subtracting *T*Δ*S* from Δ*H* yields the master Gibbs–Helmholtz equation:

Because −Δ*H*/T equals the entropy change of the constant-temperature surroundings (Δ*S*_surr), the identity Δ*G* = −*T* Δ*S*_univ holds identically at constant T and P. Furthermore, −Δ*G* equals the maximum reversible non-expansion work (w_max,non-PV)—such as electrical work w_elec = −*nFE*_cell in an electrochemical galvanic cell or mechanical/biosynthetic work driven by ATP in living cells.

Gibbs Chemical Potential (μ_i) and Thermodynamic Activity (a_i)

For an open or reactive multicomponent system containing n_i moles of chemical species i, J. Willard Gibbs (1876) defined the chemical potential (μ_i) as the partial molar Gibbs free energy at constant T, P, and moles of all other species n_j≠i:

where μ_i°(T) is the standard chemical potential of pure species i at standard pressure *P*° = 1 bar, R = 8.3144626 J mol⁻¹ K⁻¹ is the molar gas constant, and a_i is the dimensionless thermodynamic activity (a_i = f_i/*P*° ≈ P_i/1 bar for an ideal gas; a_i = γ_i [i] / 1 M for a solute with activity coefficient γ_i; and a_i = 1 for a pure solid or pure liquid solvent). By Euler's theorem for homogeneous functions of degree 1, the total Gibbs energy of a mixture is simply G = Σ_i n_i μ_i, and matter flows spontaneously from regions of higher chemical potential to regions of lower chemical potential until μ_i^(α) = μ_i^(β) across all coexisting phases.

Extent of Reaction (ξ), Reaction Quotient (Q), and Δ_r*G*° = −*RT* ln K_eq

For a general balanced chemical reaction written with signed stoichiometric coefficients ν_i (positive for products, negative for reactants; e.g., ν_A A + ν_B B ⇌ ν_C C + ν_D D) and molar extent of reaction ξ (dn_i = ν_i *d*ξ), the slope of the Gibbs energy curve with respect to reaction progress at instantaneous composition is the reaction Gibbs energy (Δ_rG):

Defining the standard reaction Gibbs energy Δ_r*G*° = Σ_i ν_i μ_i° = Σ ν_i Δ_f*G*°_i and the instantaneous reaction quotient Q = Π_i a_i^νi = (a_C^νC a_D^νD) / (a_A^|νA| a_B^|νB|) yields the van 't Hoff reaction isotherm:

Notice the essential distinction between Δ_r*G*° (a fixed constant at temperature T describing a hypothetical state where every reactant and product is at unit activity a_i = 1, Q = 1) and Δ_rG (the actual instantaneous thermodynamic driving force at current concentrations Q):

• If Q < K_eq, then Δ_rG < 0: the forward reaction (reactants → products) is exergonic (spontaneous).
• If Q > K_eq, then Δ_rG > 0: the forward reaction is endergonic, and the reverse reaction (products → reactants) proceeds spontaneously.
• At chemical equilibrium, the Gibbs energy reaches its minimum with respect to ξ (Δ_rG = 0) and the reaction quotient equals the thermodynamic equilibrium constant (Q = K_eq), yielding the foundational identity of chemical thermodynamics:

Because K_eq depends exponentially on −Δ_r*G*°/RT, at 298.15 K (RT = 2.479 kJ mol⁻¹, 2.303*RT* = 5.708 kJ mol⁻¹) a modest standard free-energy change of just −5.71 kJ mol⁻¹ shifts K_eq by a factor of 10¹, while −34.2 kJ mol⁻¹ drives the equilibrium constant to K_eq = 10⁶.

The van 't Hoff Equation and Le Chatelier's Principle

Substituting Δ_r*G*° = Δ_r*H*° − *T*Δ_r*S*° into ln K_eq = −Δ_r*G*°/(RT) and differentiating with respect to temperature T via the Gibbs–Helmholtz relation [∂(Δ*G*°/T)/∂*T*]_P = −Δ*H*°/T² yields the van 't Hoff equation:

The van 't Hoff equation provides the quantitative basis for Le Chatelier's principle regarding temperature: for an endothermic reaction (Δ_r*H*° > 0), d(ln K_eq)/dT > 0, so raising temperature increases K_eq and shifts equilibrium toward products; for an exothermic reaction (Δ_r*H*° < 0) (such as ammonia synthesis in the Haber–Bosch process, N₂ + 3 H₂ ⇌ 2 NH₃, Δ_r*H*° = −92.2 kJ mol⁻¹), raising temperature decreases K_eq and shifts equilibrium back toward reactants.

Additivity of Δ*G* and Thermodynamic Coupling in Metallurgy and Bioenergetics

Because G is a state function, an endergonic reaction (Δ*G*₁ > 0) that cannot proceed on its own can be driven forward in the same phase by coupling it through a shared chemical intermediate to a strongly exergonic reaction (Δ*G*₂ << 0) such that the net sum satisfies Δ*G*_net = Δ*G*₁ + Δ*G*₂ < 0 (and K_net = K₁ × K₂ > 1). In extractive metallurgy (Ellingham diagrams), endergonic reduction of metal oxides (such as Fe₂O₃ → 2 Fe + 3/2 O₂, Δ*G*° > 0) is coupled to the strongly exergonic oxidation of carbon to CO and CO₂. In living systems, endergonic biosynthetic reactions—such as polypeptide translation, DNA replication, and active ion pumping—are coupled via phosphorylated enzyme intermediates to the exergonic hydrolysis of ATP (Δ*G*°′ = −30.5 kJ mol⁻¹; physiological Δ*G*_p ≈ −55 kJ mol⁻¹).

What is the difference between Δ*G*° (standard Gibbs energy change) and Δ*G* (actual Gibbs energy change)?

Δ_r*G*° is a fixed reference constant at a given temperature T that specifies the free-energy change when every reactant and product is held at its standard state of unit activity (a_i = 1, i.e., 1 bar partial pressure or 1 M effective concentration, where Q = 1); it determines the equilibrium constant via Δ_r*G*° = −*RT* ln K_eq. Δ_rG = Δ_r*G*° + RT ln Q is the instantaneous driving force at the actual concentrations currently present in the reaction vessel (Q). Even if a reaction has a positive Δ_r*G*° > 0 (K_eq < 1), it will proceed spontaneously in the forward direction (Δ_rG < 0) whenever the starting product-to-reactant ratio Q is smaller than K_eq.

How can an endothermic reaction (Δ*H* > 0) be spontaneous at room temperature?

Spontaneity at constant temperature and pressure is governed not by enthalpy Δ*H* alone, but by the Gibbs free energy Δ*G* = Δ*H* − *T*Δ*S* (which reflects the total entropy change of the system plus surroundings, Δ*S*_univ = −Δ*G*/T). Whenever a reaction produces a large positive increase in system entropy (Δ*S* > 0)—for example, by dissolving an ionic salt like ammonium nitrate (NH₄NO₃) into mobile hydrated ions in water, or generating gas molecules from solids (such as Ba(OH)₂·8H₂O + 2 NH₄SCN)—the entropic term −*T*Δ*S* can be more negative than the positive +Δ*H* term once T > Δ*H*/Δ*S*, yielding Δ*G* < 0 even as the reaction absorbs heat and cools the beaker.

Does a negative Δ*G* mean that a chemical reaction will happen rapidly?

No. Thermodynamics (Δ*G*) determines whether a reaction is energetically downhill (spontaneous) and where its equilibrium lies (K_eq), whereas chemical kinetics (the activation free energy barrier Δ*G*^‡) determines how fast the reaction occurs. For example, the conversion of diamond into graphite at 298 K and 1 bar is thermodynamically spontaneous (Δ*G*° = −2.90 kJ mol⁻¹), and the oxidation of glucose or paper in air (O₂) is strongly exergonic (Δ*G*° = −2,879 kJ mol⁻¹), yet both persist indefinitely at room temperature in a kinetically metastable state because breaking their strong covalent chemical bonds requires overcoming a massive activation barrier Δ*G*^‡ according to the Eyring/Arrhenius exponential rate law k ∝ exp(−Δ*G*^‡/RT).

1. Gibbs, J. W. (1876–1878). "On the equilibrium of heterogeneous substances." Transactions of the Connecticut Academy of Arts and Sciences, 3, 108–248, 343–524. https://doi.org/10.5479/sil.421748.39088007099781
2. van 't Hoff, J. H. (1884). Études de dynamique chimique. Amsterdam: Frederik Muller & Co. https://doi.org/10.1002/bbpc.18850901902
3. Lewis, G. N., & Randall, M. (1923). Thermodynamics and the Free Energy of Chemical Substances. McGraw-Hill. https://archive.org/details/thermodynamicsfr00lewiuoft
4. Alberty, R. A. (2003). Thermodynamics of Biochemical Reactions. John Wiley & Sons. https://doi.org/10.1002/0471332607