 Debye T-cubed test

The Debye T-cubed test asks whether a bulk three-dimensional crystal's lattice molar heat capacity approaches a cubic law well below its Debye temperature $T_D$. It tests the low-frequency acoustic phonon branches in the Debye model, rather than asserting that every contribution to a real sample's heat capacity is cubic. Peter Debye's continuum description makes a striking connection: a thermal measurement can reveal the low-frequency vibrational density of states without resolving individual vibrations.^[ref]

For a monatomic bulk model, the Debye model has vibrational density of states $g(\omega )=9N{\omega }^2/{\omega }_D^3$ up to the Debye frequency ${\omega }_D$, containing $3N$ normal modes. With $x=\hslash \omega /(k_{B}T)$, the thermal vibrational internal energy becomesThe reduced Planck constant is $\hslash$, the Boltzmann constant is $k_B$, and $N$ counts chemical atoms. For $T\ll T_D$, replace the upper limit by infinity. The integral is ${\pi }^4/15$, as in the Stefan–Boltzmann law derivation. Differentiating at fixed volume gives, per mole of chemical atoms,where $R$ is the molar gas constant. The cubic exponent comes from three-dimensional acoustic phonon branches with linear low-wavevector dispersion relations. The coefficient encodes their propagation speeds. Formula-unit normalization in a compound requires its number of chemical atoms per formula unit; do not silently substitute moles of formula units for moles of chemical atoms.

For a simple metal in its nonsuperconducting state, well below both Debye temperature and electron Fermi temperature, include electronic heat capacity:This approximation also requires negligible difference between constant-pressure and constant-volume heat capacities at the measurement precision, and no significant magnetic, nuclear or Schottky anomaly contribution in the fit interval. Thus linear regression of $C_{p,\mathrm{m}}/T$ against $T^2$ separates an electronic intercept from a lattice slope. For an insulating crystal with no relevant low-energy magnetic excitations, the electronic term can instead be negligible; a superconductivity transition changes the electronic response.^[ref]

After subtracting calorimeter addenda, Brando's second gold sample, $72.6\mathrm{mg}$, gave a linear $C_{p,\mathrm{m}}/T$ versus $T^2$ plot in the nonsuperconducting state. The reported electronic intercept was $7.0\times 10^{-4}\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$; deviations from reference heat capacities were below $2\%$.^[ref] This supports the combined linear-plus-cubic description over the plotted interval. Because the platform subtraction used reference gold, it is also an instrument validation with a shared material standard, rather than a wholly independent determination of the Debye model from an uncalibrated instrument.

For a worked original inference, stipulate values of $y=C_{p,\mathrm{m}}/T$ equal to $(0.001100,0.002300,0.004300)\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$ at $T=(1,2,3)\mathrm{K}$. Assign independent $u(y)=0.000020\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$ and negligible temperature error. Equal-weight linear regression on $x=T^2$ givesThe points sit exactly on the stipulated line; these standard uncertainties come from the assigned input precision, not an estimate from zero residual scatter. With $R=8.314462618\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, inversion givesAt $2\mathrm{K}$, the electronic and lattice molar heat capacities are respectively $0.0014$ and $0.0032\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. Their equality occurs at $\sqrt{\gamma /\beta }=1.323\mathrm{K}$: even in this clean model, the total heat capacity approaches a linear, rather than cubic, law as $T$ tends to zero in the normal state.

A straight line over a narrow interval is a useful diagnostic, not a universal proof. Inspect residuals and repeat the fit over several upper and lower temperature cutoffs; curvature can signal higher lattice terms, an unremoved Schottky anomaly, or a faulty calorimeter addenda subtraction. Correlated thermometer calibration and sample-amount errors affect both axes and can dominate the illustrative linear regression uncertainty. Low-dimensional samples, additional low-energy excitations and proximity to a phase transition need other models. The test is valuable precisely because its failures can reveal physics beyond the bulk acoustic approximation.