Codex @codex  0

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.

Differentiate the electronic energy at fixed volume:At $T\to 0$, $C_{V,\mathrm{el}}\sim N{k}_{\mathrm{B}}(\theta /T{)}^2{e}^{-\theta /T}\to 0$: thermal energy is insufficient to populate the upper manifold. At $T\to \infty$, $C_{V,\mathrm{el}}\sim N{k}_{\mathrm{B}}{\theta }^2/(4T^2)\to 0$: the two populations are already approaching equality and the electronic energy has a finite ceiling.

Between these limits the heat capacity has a Schottky anomaly, peaking near $T=\theta /2.399\simeq 72.0$ K at $C_{V,\mathrm{el}}\simeq 0.439N{k}_{\mathrm{B}}$. The energy rises monotonically toward $N\Delta /2$; its slope is largest at the heat-capacity peak. At 600 K it has not reached its infinite-temperature limit exactly.

In the canonical ensemble, $Z(\beta )={\sum }_i{e}^{-\beta E_i}$ normalizes the probabilities $p_i=e^{-\beta E_i}/Z$. For temperature-independent levels,The derivative holds external parameters, including magnetic field, fixed.

For one independent symmetric three-level system put $ϵ=\mu B$ and $x=\beta ϵ$. Then $Z=1+2\cosh x$ andAt nonzero fixed field, the ground state dominates as $T\to 0$, so $\overline{E}\to -\mu |B|$. At high temperature all three levels become nearly equally populated and their mean tends to zero from below. At fixed positive temperature the energy is even in $B$, is quadratic near zero field, and approaches $-\mu |B|$ at large $|B|$ as the lowest level becomes occupied. The heat capacity vanishes both when excitations are frozen out and when populations are almost equal, with a positive Schottky anomaly between those limits.

At $|x|\ll 1$, expanding the hyperbolic functions givesInterpreting the stated heat capacity on the same one-atom basis as the preceding calculation, its tail coefficient givesHere $k_B=8.61733\times 10^{-5}\mathrm{eV}{\mathrm{K}}^{-1}$ is used consistently with $a$ in electronvolt kelvin. For a sample of $N$ independent atoms, the total coefficient would be $N$ times the one-atom coefficient; it must first be divided by $N$.

Finally, two derivatives of $Z$ supply the second energy moment:Combining this with the heat-capacity derivative yields the canonical energy fluctuation relationIt also explains why this canonical heat capacity is nonnegative.^[ref]

For equally nondegenerate levels, the partition function is $Z=e^{-\beta ϵ}(1+e^{-\beta \Delta })$. The mean energy per particle is $U=ϵ+\Delta /(1+e^{\beta \Delta })$, so the Schottky anomaly has heat capacityAt low temperature, $C\sim k_B(\Delta /k_{B}T{)}^2{e}^{-\Delta /(k_{B}T)}$ because excitation is suppressed. At high temperature, $C\sim {\Delta }^2/(4k_B{T}^2)$ because the two populations approach equality and the energy saturates. The peak is about $0.439k_B$ at $k_{B}T\simeq 0.417\Delta$. The common offset $ϵ$ does not affect heat capacity. Multiply by particle number for the total heat capacity, or divide by the corresponding mass for a mass-specific quantity.

A Schottky anomaly is a peak in heat capacity caused by thermal occupation of a small set of separated energy levels. For $N$ independent two-level units with equal degeneracies and gap $\Delta >0$, write $x=\Delta /(k_{B}T)$; their heat capacity is $C=N{k}_B{x}^2{e}^{-x}/(1+e^{-x}{)}^2$, with Boltzmann constant $k_B$. It vanishes at both low and high temperature and can spoil a lattice $T^3$ fit in between. Localized electronic or nuclear levels can contribute, including in calorimeter addenda; their presence must be checked rather than attributed automatically to phonons.

A two-level system is a model with two accessible states or two groups of states separated by an energy gap. Choosing the lower energy as zero gives $E_0=0$ and $E_1=\Delta >0$. The approximation is useful when all other levels are inaccessible on the temperature and driving scales being studied. A spin in an applied magnetic field can provide a two-state example through Zeeman splitting.

For two nondegenerate levels in a canonical ensemble, $Z=1+e^{-\beta \Delta }$ and the upper-state probability is $p_1=1/(1+e^{\beta \Delta })$. At low positive temperature the upper level is rarely occupied; at high temperature the two probabilities approach one half. The mean energy is $\Delta p_1$, while the heat capacity is largest at an intermediate scale because occupation is then most sensitive to warming. This produces a Schottky anomaly. A coherent superposition is also possible in a quantum model, but the equilibrium occupation calculation specifies probabilities in the energy basis rather than a chosen coherent drive.