4.1 Specific heat capacity

Specific heat capacity describes how much heat energy a material needs to absorb for its temperature to rise by a given amount. In solid state physics, understanding specific heat is essential because it reveals how energy is stored in a crystal lattice and how atoms vibrate at different temperatures. Classical predictions work well at high temperatures, but quantum models from Einstein and Debye are needed to explain what happens as temperature drops.

Definition of specific heat capacity

Specific heat capacity quantifies the heat required to raise the temperature of a unit mass of a substance by one degree. The mathematical definition is:

$C = \frac{Q}{m \Delta T}$

where $Q$ is the heat added, $m$ is the mass, and $\Delta T$ is the temperature change.

In solid state physics, you'll more often work with molar specific heat (heat per mole per degree) rather than per-unit-mass, since we're comparing elements and compounds at the atomic level. The distinction matters when you're switching between formulas.

Factors affecting specific heat capacity

Atomic mass and specific heat capacity

Specific heat capacity generally decreases with increasing atomic mass. Heavier atoms vibrate at lower frequencies, so each vibrational mode requires less energy to excite. This is why lead (207 u) has a much lower specific heat per gram than aluminum (27 u).

The Dulong-Petit law captures this trend at high temperatures: the molar specific heat of most elemental solids approaches $3R \approx 24.9 \text{ J mol}^{-1}\text{K}^{-1}$, where $R$ is the universal gas constant. On a per-gram basis, light elements end up with higher specific heats because the same 3R of energy is spread over less mass.

Crystal structure and specific heat capacity

• Crystal structure determines the available vibrational modes (phonon dispersion relations), which directly affect how energy is distributed among atoms.
• Highly symmetric structures like cubic crystals tend to have more equivalent vibrational modes, while less symmetric structures (hexagonal, for example) can show anisotropy, meaning the specific heat may depend on the direction of heat flow in a single crystal.
• These structural effects are most noticeable at intermediate temperatures where not all phonon modes are fully excited.

Classical theory of specific heat

Dulong-Petit law

The classical approach treats each atom in a solid as a harmonic oscillator vibrating in three independent directions (x, y, z). Each direction contributes $k_BT$ of energy on average ($\frac{1}{2}k_BT$ kinetic + $\frac{1}{2}k_BT$ potential), giving a total energy per atom of $3k_BT$. For one mole of atoms, the total thermal energy is $U = 3N_Ak_BT = 3RT$, so the molar heat capacity is:

$C_V = \frac{dU}{dT} = 3R \approx 24.9 \text{ J mol}^{-1}\text{K}^{-1}$

This works surprisingly well for most solids at room temperature and above their Debye temperature.

Limitations of classical theory

• The classical result predicts a constant $C_V$ at all temperatures. Experiments show that specific heat drops toward zero as $T \to 0$, which the classical picture cannot explain.
• The equipartition theorem assigns $k_BT$ of energy to every mode regardless of temperature, but at low temperatures some vibrational modes simply aren't accessible because the available thermal energy is too small to excite them.
• Quantum mechanics is needed to fix this: energy comes in discrete packets (phonons), and modes with energies much larger than $k_BT$ remain "frozen out."

Einstein model of specific heat

Assumptions in Einstein model

Einstein's key simplification was to treat all $N$ atoms in a solid as independent quantum harmonic oscillators, each vibrating at the same frequency $\omega_E$. Each oscillator has quantized energy levels:

$E_n = \left(n + \frac{1}{2}\right)\hbar\omega_E$

where $n = 0, 1, 2, \ldots$ and $\hbar$ is the reduced Planck constant. At low temperatures, most oscillators sit in their ground state because $k_BT$ is too small to promote them to higher levels.

Einstein temperature and specific heat

The Einstein temperature sets the energy scale for these oscillators:

$\Theta_E = \frac{\hbar\omega_E}{k_B}$

When $T \gg \Theta_E$, the oscillators behave classically and you recover the Dulong-Petit value. When $T \ll \Theta_E$, the modes freeze out and $C$ drops exponentially.

The full expression for the molar heat capacity is:

$C_E = 3Nk_B\left(\frac{\Theta_E}{T}\right)^2 \frac{e^{\Theta_E/T}}{\left(e^{\Theta_E/T} - 1\right)^2}$

Successes of Einstein model

• First model to correctly predict that $C \to 0$ as $T \to 0$.
• Introduced the idea that quantized vibrational energy explains the temperature dependence of specific heat.
• Gives a reasonable fit to experimental data across a wide temperature range.

Limitations of Einstein model

• All atoms vibrate at one frequency, which isn't realistic. Real solids have a broad spectrum of vibrational frequencies.
• At low temperatures, the model predicts an exponential decay ($C \propto e^{-\Theta_E/T}$), but experiments show a slower, power-law decrease ($C \propto T^3$). The single-frequency assumption misses the low-frequency acoustic phonons that dominate at low $T$.

Debye model of specific heat

Debye temperature and specific heat

Debye improved on Einstein's approach by treating the solid as an elastic continuum with a spectrum of vibrational frequencies, from zero up to a maximum cutoff frequency $\omega_D$. The Debye temperature is:

$\Theta_D = \frac{\hbar\omega_D}{k_B}$

This temperature marks the boundary between quantum and classical behavior. For $T \gg \Theta_D$, all phonon modes are fully excited and you get the classical $3R$ result. For $T \ll \Theta_D$, only the lowest-frequency modes contribute.

The molar heat capacity in the Debye model is:

$C_D = 9Nk_B\left(\frac{T}{\Theta_D}\right)^3 \int_0^{\Theta_D/T} \frac{x^4 e^x}{(e^x - 1)^2}\, dx$

where $x = \frac{\hbar\omega}{k_BT}$.

Debye $T^3$ law at low temperatures

When $T \ll \Theta_D$, the upper limit of the integral goes to infinity and can be evaluated analytically. The result is:

$C_D \propto T^3$

This Debye $T^3$ law matches experimental low-temperature data extremely well and is one of the major successes of the model. Physically, it reflects the fact that only long-wavelength acoustic phonons (with low frequencies and a density of states that goes as $\omega^2$) are excited at very low temperatures.

Successes of Debye model

• Correctly predicts the $T^3$ behavior at low temperatures, unlike the Einstein model's exponential decay.
• Accounts for the full spectrum of acoustic phonon frequencies rather than assuming a single frequency.
• Provides a good interpolation between the low-temperature quantum regime and the high-temperature classical limit.

Limitations of Debye model

• Assumes a linear (acoustic) dispersion relation $\omega = v|\mathbf{k}|$ for all modes, which breaks down near the Brillouin zone boundary.
• Ignores optical phonon branches, which can contribute significantly in crystals with more than one atom per unit cell.
• The single cutoff frequency $\omega_D$ is an approximation; real phonon densities of states have more complex shapes with Van Hove singularities.

Electronic contribution to specific heat

Free electron model and specific heat

In metals, conduction electrons also carry thermal energy. You might expect each electron to contribute $\frac{3}{2}k_BT$ (classical ideal gas), but this would give a much larger specific heat than observed. The resolution comes from quantum statistics.

Electrons obey Fermi-Dirac statistics, and at typical temperatures only a small fraction of electrons (those within $\sim k_BT$ of the Fermi energy) can absorb thermal energy. The vast majority are deep in the Fermi sea and are blocked by the Pauli exclusion principle. This gives an electronic specific heat:

$C_{el} = \frac{\pi^2}{2} k_B^2 N(E_F) T$

where $N(E_F)$ is the density of states at the Fermi energy.

Fermi energy and electronic specific heat

The Fermi energy $E_F$ is the energy of the highest occupied electron state at absolute zero. Typical values for metals are a few eV (e.g., $E_F \approx 7.0$ eV for copper), corresponding to Fermi temperatures $T_F = E_F/k_B$ of tens of thousands of kelvin.

Because room temperature is far below $T_F$, only a tiny fraction ($\sim T/T_F$) of electrons participate in thermal processes. Metals with a higher $N(E_F)$ have a larger electronic specific heat; transition metals with their narrow d-bands are a good example.

Temperature dependence of electronic specific heat

The electronic contribution is linear in $T$, while the lattice (phonon) contribution goes as $T^3$ at low temperatures. The total low-temperature specific heat of a metal is therefore:

$C = \gamma T + \alpha T^3$

where $\gamma$ (the Sommerfeld coefficient) captures the electronic part and $\alpha$ captures the Debye lattice part.

At very low temperatures (typically below ~10 K), the linear electronic term dominates because $T^3$ shrinks faster than $T$. Plotting $C/T$ versus $T^2$ gives a straight line: the y-intercept is $\gamma$ and the slope is $\alpha$. This is a standard experimental technique for separating the two contributions.

Experimental techniques for measuring specific heat

Adiabatic calorimetry

1. Place the sample in a thermally insulated (adiabatic) container so no heat escapes to the surroundings.
2. Supply a known amount of electrical energy $Q$ to the sample using a heater.
3. Measure the resulting temperature rise $\Delta T$.
4. Calculate specific heat from $C = Q / (m \Delta T)$.

This is the most straightforward method and works well over a wide temperature range, though maintaining truly adiabatic conditions becomes challenging at very low temperatures.

Differential scanning calorimetry (DSC)

1. Heat the sample and an inert reference material at the same constant rate.
2. Measure the difference in heat flow needed to keep both at the same temperature.
3. Extract the specific heat from the heat flow difference and the known heating rate.

DSC is widely used for studying phase transitions because it directly reveals anomalies in heat flow at transition temperatures.

Relaxation calorimetry

1. Apply a small heat pulse to the sample, raising its temperature slightly above the bath temperature.
2. Monitor how the sample temperature decays back to equilibrium over time.
3. Extract the specific heat from the thermal relaxation time constant and the known thermal link conductance between the sample and the bath.

This technique is especially useful for very small samples and at very low temperatures (millikelvin range), where adiabatic methods become impractical.

Applications of specific heat in solids

Thermal energy storage

Materials with high specific heat capacities can absorb and release large amounts of thermal energy with relatively small temperature swings. Phase change materials (PCMs) take this further by exploiting latent heat during melting/solidification. Common PCMs include paraffin wax, salt hydrates, and fatty acids, used in applications from building insulation to spacecraft thermal regulation.

Thermal management in electronics

Electronic components generate significant heat during operation, and materials with high specific heat (combined with high thermal conductivity) are used as heat sinks to absorb transient thermal loads. Copper, aluminum, and diamond are common choices. The specific heat determines how much temperature rise occurs before steady-state heat dissipation takes over.

Specific heat and phase transitions

Specific heat measurements are a powerful probe of phase transitions in solids:

• A first-order transition (like melting) produces a sharp spike or divergence in $C$ because latent heat is absorbed at a single temperature.
• A second-order (continuous) transition shows a discontinuity or lambda-shaped peak in $C$ without latent heat. The classic example is the superconducting transition, where $C$ jumps at $T_c$.
• Magnetic ordering transitions (ferromagnetic, antiferromagnetic) and structural transitions also produce characteristic specific heat anomalies that help identify the transition type and critical temperature.