Heat Capacities & Temperature Dependence

SECTION 1

Historical Context & Motivation

The concept of heat capacity arose from early attempts to quantify the relationship between heat supplied to a body and the resulting temperature change. Before thermodynamics was formalized as a discipline, natural philosophers struggled with a deceptively simple question: why do different substances warm at different rates when identical amounts of heat are applied? The answer to this question would ultimately require centuries of experimental refinement, leading from qualitative observations about fire and warmth to the precise mathematical framework that underpins modern thermochemistry and chemical engineering. Understanding heat capacity proved essential not only for explaining calorimetric measurements but also for predicting how enthalpy changes with temperature — a capability central to reactor design, atmospheric modeling, and materials science.

1760

Black's Specific Heat

Joseph Black distinguishes between heat and temperature, introducing the concept of specific heat and demonstrating that equal masses of different substances require different amounts of heat to produce the same temperature rise.

1842

Mayer & Mechanical Equivalent

Julius Robert Mayer recognizes the relationship between C_P and C_V for gases, linking the difference to the work of expansion and providing early evidence for energy conservation.

1907

Einstein's Quantum Model

Albert Einstein applies Planck's quantization to the vibrational modes of a solid, explaining why heat capacities of solids decrease at low temperatures — a phenomenon classical equipartition could not account for.

1912

Debye's Improved Model

Peter Debye refines Einstein's approach by treating lattice vibrations as a spectrum of phonon frequencies, yielding the famous T³ dependence of heat capacity at very low temperatures and providing excellent agreement with experiment.

1932

Kirchhoff's Law in Standard Tables

Empirical polynomial expressions for C_P(T) become standardized in thermodynamic data compilations, enabling engineers and chemists to compute enthalpy changes at arbitrary temperatures using Kirchhoff's equation.

The central question that these historical developments converge upon is both practical and profound: given that chemical reactions rarely occur at the standard reference temperature of 298.15 K, how can we reliably predict the enthalpy of reaction at any temperature? The answer lies in the systematic use of heat capacity as a bridge between thermodynamic data measured at one temperature and the conditions of actual interest.

SECTION 2

Core Principles & Definitions

Heat capacity quantifies how much energy a system must absorb to raise its temperature by a given amount. At the molecular level, this energy distributes among translational, rotational, and vibrational degrees of freedom, and the relative contributions of these modes change with temperature. To develop a rigorous treatment, we must distinguish between heat capacity measured at constant volume and at constant pressure, because these two conditions yield fundamentally different thermodynamic quantities.

Heat Capacity at Constant Volume (C_V)

Defined as (∂U/∂T)_V, measuring how internal energy changes with temperature when no expansion work is performed. For an ideal gas, C_V depends only on the molecular degrees of freedom.

Heat Capacity at Constant Pressure (C_P)

Defined as (∂H/∂T)_P, measuring how enthalpy changes with temperature when pressure is held fixed. Since most laboratory and industrial processes occur at constant pressure, C_P is the quantity most frequently tabulated.

C_P − C_V Relationship

For an ideal gas, C_P − C_V = nR, because the extra energy at constant pressure goes into PV work of expansion. For real gases and condensed phases, a more general expression involving the thermal expansion coefficient and compressibility applies.

Molar vs. Specific Heat Capacity

Molar heat capacity (units: J mol⁻¹ K⁻¹) is the capacity per mole of substance, while specific heat capacity (J g⁻¹ K⁻¹) is per unit mass. In thermochemistry, molar quantities are standard; in engineering, mass-based quantities are often more practical.

Temperature Dependence

Heat capacities are not constants — they vary with temperature as higher-energy molecular modes become accessible. This dependence is commonly modeled by polynomial expressions of the form C_P(T) = a + bT + cT² or with inverse-square terms.

✦ KEY TAKEAWAY

Think of heat capacity as a thermal "stiffness" of a substance. A material with a high heat capacity is like a thick sponge that absorbs a large quantity of water (energy) before it feels saturated (changes temperature appreciably). When you need to know the enthalpy at a different temperature, the heat capacity tells you how much enthalpy accumulates per degree — just as knowing the sponge's absorption rate tells you how much water it holds after being under the tap for a given time.

SECTION 3

Visual Explanation — C_P(T) Curves for Common Gases

Molar heat capacity at constant pressure as a function of temperature for argon (monatomic, cyan), nitrogen (diatomic, violet), and carbon dioxide (triatomic, pink). Monatomic gases exhibit a temperature-independent C_P of 5/2 R. Polyatomic molecules show increasing C_P with temperature as vibrational modes become thermally activated.

The diagram above illustrates a central result of statistical thermodynamics: the number and nature of molecular degrees of freedom directly control the heat capacity of a gas. A monatomic ideal gas such as argon possesses only three translational degrees of freedom, each contributing ½R to C_V via the equipartition theorem, yielding C_V = 3/2 R and C_P = 5/2 R, both independent of temperature. Diatomic molecules such as N₂ add two rotational modes at room temperature (contributing R to C_V), and their vibrational mode begins to contribute appreciably above roughly 1000 K. For polyatomic molecules like CO₂, with its multiple vibrational modes, the rise in C_P with temperature is pronounced, reflecting progressive activation of bending and stretching vibrations. This temperature dependence is precisely what necessitates the use of polynomial expressions for C_P(T) in thermochemical calculations.

SECTION 4

Mathematical Framework

The mathematical treatment of heat capacity and its connection to enthalpy begins with the definition of enthalpy and proceeds through the fundamental relation between heat capacity at constant pressure and the temperature dependence of H. We then derive the key working equations — Kirchhoff's equation — that allow computation of reaction enthalpies at temperatures other than the standard reference.

DEFINITION OF C_P

C_P = (∂H / ∂T)_P

C_P is the heat capacity at constant pressure (J mol⁻¹ K⁻¹); H is enthalpy; T is absolute temperature (K). This definition implies that at constant pressure, dH = C_P dT.

ENTHALPY AT TEMPERATURE T₂

H(T₂) = H(T₁) + ∫[T₁ → T₂] C_P dT

By integrating the differential relation dH = C_P dT from temperature T₁ to T₂, one obtains the enthalpy at any new temperature given a known reference value H(T₁). If C_P is constant, the integral reduces to C_P(T₂ − T₁).

EMPIRICAL C_P POLYNOMIAL

C_P(T) = a + bT + cT⁻²

Standard thermodynamic tables report the empirical coefficients a (J mol⁻¹ K⁻¹), b (J mol⁻¹ K⁻²), and c (J mol⁻¹ K) for each substance. Some sources use a variant form with a cT² term or additional terms. The Shomate equation (used by NIST) employs C_P = A + BT + CT² + DT³ + E/T² with T in units of T/1000.

KIRCHHOFF'S EQUATION

ΔH°(T₂) = ΔH°(T₁) + ∫[T₁ → T₂] ΔC_P dT

Here ΔC_P = Σ ν_i C_P,i is the difference in heat capacities between products and reactants, weighted by stoichiometric coefficients ν_i (positive for products, negative for reactants). This is the master equation for translating standard enthalpies of reaction from the tabulated temperature to any other temperature.

The derivation of Kirchhoff's equation follows directly from the linearity of integration and the fact that enthalpy is a state function. Since ΔH° = Σ ν_i H_i°, differentiating with respect to T at constant P gives d(ΔH°)/dT = Σ ν_i C_P,i = ΔC_P. Integration from T₁ to T₂ then yields the result above. When the polynomial form of C_P(T) is substituted, the integration is straightforward and produces terms involving ΔaT, ½Δb T², and −Δc T⁻¹ evaluated at the limits.

SECTION 5

Detailed Breakdown — Polynomial Coefficients & Kirchhoff Integration

To use Kirchhoff's equation in practice, one must evaluate the integral of ΔC_P(T) analytically. When C_P is expressed in the standard polynomial form for each species, the difference ΔC_P inherits the same polynomial structure with coefficients Δa, Δb, and Δc. The following table lists representative polynomial parameters for several common substances, enabling manual computation of enthalpy changes across temperature ranges.

Representative C_P polynomial coefficients (valid approximately 298–2000 K). Data from Atkins' Physical Chemistry.

Substance	a / J mol⁻¹ K⁻¹	b × 10³ / J mol⁻¹ K⁻²	c × 10⁻⁵ / J mol⁻¹ K
H₂(g)	27.28	3.26	0.50
O₂(g)	29.96	4.18	−1.67
N₂(g)	28.58	3.77	−0.50
CO₂(g)	44.22	8.79	−8.62
H₂O(g)	30.54	10.29	0.00

A Hess's law cycle illustrating Kirchhoff's equation. The unknown ΔH°(T₂) (bottom horizontal arrow) equals ΔH°(T₁) minus the heating integral for reactants plus the heating integral for products. The net effect is adding ∫ΔC_P dT to the known ΔH°(T₁).

The pathway diagram above captures the essential logic of Kirchhoff's equation as a thermodynamic cycle. Since enthalpy is a state function, the change ΔH°(T₂) can be computed by any convenient path connecting the same initial and final states. The cycle proceeds as follows: cool reactants from T₂ to T₁ (subtracting ∫C_P(reactants) dT), carry out the reaction at T₁ using the known ΔH°(T₁), then heat the products from T₁ to T₂ (adding ∫C_P(products) dT). Combining these three legs recovers ΔH°(T₂) = ΔH°(T₁) + ∫ΔC_P dT, where ΔC_P = C_P(products) − C_P(reactants). This result is remarkably powerful: it means that a single tabulated ΔH° at 298 K, combined with polynomial C_P data, suffices to predict the enthalpy of reaction at any temperature within the valid range of the polynomials.

SECTION 6

Worked Example — Applying Kirchhoff's Equation

Consider the combustion of carbon monoxide: CO(g) + ½O₂(g) → CO₂(g). The standard enthalpy of reaction at 298 K is ΔH°(298) = −283.0 kJ mol⁻¹. We wish to calculate ΔH° at 500 K using the following simplified C_P data (in J mol⁻¹ K⁻¹): CO: a = 28.16, b = 1.675 × 10⁻³; O₂: a = 29.96, b = 4.18 × 10⁻³; CO₂: a = 44.22, b = 8.79 × 10⁻³. For simplicity, we neglect the c/T² terms.

ΔH° of CO Combustion at 500 K

Step 1 — Compute Δa and ΔbUsing ΔC_P = Σ ν_i C_P,i, compute the Δ coefficients: Δa = a(CO₂) − a(CO) − ½ a(O₂) = 44.22 − 28.16 − ½(29.96) = 44.22 − 28.16 − 14.98 = 1.08 J mol⁻¹ K⁻¹. Similarly, Δb = 8.79 × 10⁻³ − 1.675 × 10⁻³ − ½(4.18 × 10⁻³) = 8.79 − 1.675 − 2.09 = 5.025 × 10⁻³ J mol⁻¹ K⁻².

Δa = 1.08 J mol⁻¹ K⁻¹; Δb = 5.025 × 10⁻³ J mol⁻¹ K⁻²

Step 2 — Set Up the IntegralThe integral of ΔC_P dT from T₁ = 298 K to T₂ = 500 K is: ∫(Δa + Δb T) dT = Δa(T₂ − T₁) + ½ Δb(T₂² − T₁²). Substituting: ΔT = 500 − 298 = 202 K, and T₂² − T₁² = 250000 − 88804 = 161196 K².

ΔT = 202 K; T₂² − T₁² = 161 196 K²

Step 3 — Evaluate the IntegralFirst term: Δa × ΔT = 1.08 × 202 = 218.2 J mol⁻¹. Second term: ½ × 5.025 × 10⁻³ × 161196 = 404.8 J mol⁻¹. Total integral = 218.2 + 404.8 = 623.0 J mol⁻¹ = 0.623 kJ mol⁻¹.

∫ΔC_P dT = +0.623 kJ mol⁻¹

Step 4 — Apply Kirchhoff's EquationΔH°(500) = ΔH°(298) + ∫ΔC_P dT = −283.0 + 0.623 = −282.4 kJ mol⁻¹. Notice that ΔC_P is positive (products have higher heat capacity than reactants for this reaction), so the exothermic enthalpy becomes slightly less negative at higher temperature — the reaction releases slightly less heat at 500 K than at 298 K.

ΔH°(500 K) = −282.4 kJ mol⁻¹

💡 Physical Interpretation

When ΔC_P > 0, products absorb heat more effectively than reactants. As you go to higher T, the products' enthalpy rises faster than the reactants', making the exothermic ΔH° less negative (or an endothermic ΔH° more positive). Conversely, when ΔC_P < 0, the reaction becomes more exothermic with increasing temperature.

SECTION 7

Strengths, Limitations & Approximation Methods

Kirchhoff's equation is exact in principle, but its practical application involves several approximation choices. The accuracy of the final result depends on the quality of the C_P(T) data, the validity of the polynomial fit over the temperature range of interest, and whether phase transitions occur within that range. The table below contrasts three common levels of approximation used in thermochemical calculations.

Levels of approximation in applying Kirchhoff's equation

Approximation Level	Method	When Appropriate
Zeroth-order	Assume ΔC_P ≈ 0, so ΔH° is temperature-independent.	Very small ΔT (< 50 K), or when ΔC_P is genuinely negligible (e.g., reactions with similar molecular complexity on both sides).
Constant ΔC_P	Use C_P values at 298 K and treat ΔC_P as a constant: ΔH°(T₂) = ΔH°(298) + ΔC_P × ΔT.	Moderate temperature ranges (≈ 200 K) and back-of-the-envelope estimates. Quick checks before committing to polynomial integration.
Full polynomial	Integrate C_P(T) = a + bT + cT⁻² analytically or use Shomate parameters with 5-term polynomials.	Large temperature ranges (> 200 K), high-accuracy work, or when species have strongly temperature-dependent C_P (e.g., polyatomic gases).

✦ KEY TAKEAWAY

Choosing an approximation level is analogous to choosing a map resolution for navigation. A zeroth-order approximation is like a globe — it shows you the general direction but lacks detail. The constant ΔC_P approach is a road atlas — adequate for regional travel but potentially misleading on winding mountain passes. The full polynomial integration is a GPS with real-time data — computationally heavier but reliably accurate across the entire route. In research and engineering, the cost of using the wrong approximation can be a reactor designed for the wrong temperature profile.

⚠️ Phase Transitions

If a phase transition occurs between T₁ and T₂, the integral must be split at the transition temperature T_tr, and the enthalpy of transition (ΔH_tr) must be added as a separate term. Different C_P polynomial coefficients apply to each phase.

SECTION 8

Connection to Advanced Theory — Statistical Mechanics & Beyond

The empirical polynomial expressions for C_P(T) used throughout classical thermochemistry have a deeper theoretical foundation in statistical mechanics. From the partition function of a molecular system, one can derive the heat capacity directly, providing not just numerical values but physical insight into why C_P increases with temperature and why it eventually levels off as all modes become fully excited. This connection bridges classical thermodynamics and molecular-level theory, and is explored in depth in statistical thermodynamics courses.

Classical vs. statistical mechanical perspectives on heat capacity

Feature	Classical Thermodynamics	Statistical Mechanics
C_P source	Empirical polynomial fit to experimental calorimetric data	Derived from the molecular partition function: C_V = T(∂²A/∂T²)_V or from ⟨E²⟩ − ⟨E⟩²
Temperature dependence	Captured via fitted coefficients (a, b, c) — phenomenological	Emerges naturally from Boltzmann population of quantized energy levels; explains why modes freeze out at low T
Predictive power	Interpolative within fitted range; extrapolation risky	Fully predictive if molecular spectroscopic constants (vibrational frequencies, moments of inertia) are known
Solids at low T	Dulong–Petit limit (3R per atom) assumed; deviations noted empirically	Einstein and Debye models quantitatively explain C → 0 as T → 0, including the T³ Debye law

Looking forward, the concept of temperature-dependent heat capacity connects to several advanced topics. In chemical equilibrium, the van 't Hoff equation d(ln K)/dT = ΔH°/RT² couples with Kirchhoff's equation to predict how equilibrium constants change with temperature when ΔH° itself is temperature-dependent. In computational chemistry, ab initio calculation of vibrational frequencies enables prediction of C_P(T) for species that have never been experimentally characterized. In materials science, understanding heat capacity is essential for modeling thermal conductivity, phase stability, and the thermodynamic properties of alloys and ceramics at extreme temperatures.

SECTION 9

Practice Problems

PROBLEM 1 — CONCEPTUAL

For an ideal monatomic gas, C_P = 5/2 R and is independent of temperature. Explain, using the equipartition theorem, why polyatomic gases have both a higher C_P and a temperature-dependent C_P, while monatomic gases do not.

PROBLEM 2 — BASIC CALCULATION

The molar C_P of liquid water is approximately 75.3 J mol⁻¹ K⁻¹ and can be treated as constant over small temperature ranges. Calculate the molar enthalpy change when 1 mol of liquid water is heated from 25 °C to 80 °C at constant pressure.

PROBLEM 3 — INTERMEDIATE

For gaseous SO₂, the heat capacity polynomial is C_P(T) = 25.72 + 57.95 × 10⁻³ T − 38.05 × 10⁻⁶ T² (J mol⁻¹ K⁻¹, T in K). Calculate the enthalpy change H(600) − H(298) for one mole of SO₂.

PROBLEM 4 — APPLIED

The water-gas shift reaction is CO(g) + H₂O(g) → CO₂(g) + H₂(g). Given ΔH°(298) = −41.2 kJ mol⁻¹ and the following 298 K molar C_P values (J mol⁻¹ K⁻¹): CO = 29.1, H₂O(g) = 33.6, CO₂ = 37.1, H₂ = 28.8, estimate ΔH° at 700 K using the constant-ΔC_P approximation. Discuss whether this approximation is likely to be accurate over a 400 K range.

PROBLEM 5 — CRITICAL THINKING

Consider a hypothetical reaction where ΔC_P changes sign over the temperature range of interest — positive at low T and negative at high T. Describe qualitatively how ΔH°(T) would behave as a function of temperature. Could there be a temperature at which ΔH° reaches an extremum? If so, what thermodynamic condition defines that temperature? Relate your answer to Le Chatelier's principle.

SUMMARY

Summary — Heat Capacities & Temperature Dependence of Enthalpy

Heat capacity measures the energy required to raise a system's temperature by one degree, with C_P = (∂H/∂T)_P being the quantity most relevant to constant-pressure processes in chemistry. For ideal gases, C_P − C_V = nR, reflecting the work of expansion. The temperature dependence of C_P arises from the progressive activation of molecular rotational and vibrational degrees of freedom, as described by the equipartition theorem and its quantum corrections. Practically, C_P(T) is expressed as an empirical polynomial of the form a + bT + cT⁻².

The central working equation is Kirchhoff's equation: ΔH°(T₂) = ΔH°(T₁) + ∫ΔC_P dT, which exploits the state-function nature of enthalpy to translate reaction enthalpies from a reference temperature to any temperature of interest. Three levels of approximation are available: the zeroth-order (ΔC_P ≈ 0), constant ΔC_P, and full polynomial integration, each appropriate for progressively larger temperature ranges and accuracy requirements. Mastery of these tools is foundational for chemical engineering design, atmospheric chemistry, and the statistical mechanical framework that provides their molecular interpretation.