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.

Electronic heat capacity is the contribution from electronic excitations to heat capacity. A simple normal metal with a smooth density of states near its Fermi energy has $C_{\mathrm{m}}^{\mathrm{el}}\simeq \gamma T$ well below its electron Fermi temperature. Only a narrow energy interval near the Fermi energy changes occupation on heating, explaining the small fraction of active Electrons. The Debye T-cubed test separates this contribution from lattice phonons. Interactions, a superconductivity transition or other electronic structure can invalidate a constant $\gamma$ description.

Low-temperature calorimetry measures a small sample's heat capacity by resolving a small temperature change after a known electrical energy input. The experimental challenge is that the support, heater and thermometer also warm up. Their calorimeter addenda must be removed before interpreting the sample's phonons or electronic heat capacity. The Debye T-cubed test turns the corrected response into a test of the low-temperature vibrational spectrum.

Manuel Brando's 2009 experiment used a sapphire platform suspended by four alloyed platinum wires from a copper ring fixed to a cryostat stage. Grease attached the sample; a chip heater and ruthenium-oxide resistance thermometer served platform PLS-2. A controller stabilized the bath; a lock-in amplifier supplied heating pulses and read the calibrated sensor. The steady rise calibrated thermal conductance as $G=P/\Delta T$. A separate $59.1\mathrm{mg}$, $99.99\%$ pure gold sample, compared with reference heat capacities, established the calorimeter addenda. Small pulses limited property variation; imperfect sample contact required a second relaxation time. The platform contribution rose below $1.5\mathrm{K}$, making simple extrapolation unsafe.^[ref]

For an original idealization, assume sample and platform share one temperature $T$, the bath remains at $T_b$, and both the total heat capacity and the weak-link thermal conductance $G$ are approximately constant across one small pulse. Ignore other heat transfer and any phase transition. The first law of thermodynamics becomesHere $C_s$ is sample heat capacity, $C_a$ is calorimeter addenda, $P$ is deposited electrical power from Joule heating, and $\tau$ is a time constant. This is the same balance as calorimeter heat-loss correction, used here with a deliberately weak cryogenic link and a separately subtracted platform. For a constant pulse $P_0$ beginning at bath temperature, then switched off at $t_h$,The cooling exponential decay determines $\tau$, while the heating plateau determines $G=P_0/(T_{\infty }-T_b)$. Together with calorimeter addenda, they give $C_s=G\tau -C_a$. Cooling alone determines a ratio, not an absolute heat capacity. With finite $t_h$, the observed end rise is smaller than the plateau; fit the full response or apply the displayed correction rather than identifying the two.

For the stipulated figure parameters, $G=5.00\times 10^{-8}\mathrm{W}{\mathrm{K}}^{-1}$, $\tau =20.0\mathrm{s}$, and $C_a=2.00\times 10^{-7}\mathrm{J}{\mathrm{K}}^{-1}$ give $C_s=8.00\times 10^{-7}\mathrm{J}{\mathrm{K}}^{-1}$. The asymptotic rise is $0.0400\mathrm{K}$, only $2\%$ of bath temperature. At $100\mathrm{s}$, the rise is $0.0397305\mathrm{K}$. Of $2.00\times 10^{-7}\mathrm{J}$ input, only $3.97305\times 10^{-8}\mathrm{J}$ remains stored: the other $1.60270\times 10^{-7}\mathrm{J}$ has flowed to the bath. Treating this long pulse as an isolated $Q/C$ experiment would overestimate total heat capacity by a factor of $5.034$.

Suppose independent standard uncertainties are $u(G)=0.05\times 10^{-8}\mathrm{W}{\mathrm{K}}^{-1}$, $u(\tau )=0.2\mathrm{s}$ and $u(C_a)=0.10\times 10^{-7}\mathrm{J}{\mathrm{K}}^{-1}$. First-order propagation of uncertainty givesThe uncertainty is one standard deviation under those assumptions, about $2.2\%$ of sample heat capacity. Shared thermometer calibration errors can correlate fitted $G$ and $\tau$; then their covariance contributes too. Smaller samples amplify the relative influence of calorimeter addenda. Proposed checks include varying pulse amplitude, repeating with a reference sample, checking sensor self-heating, and comparing single- and two-time-constant fits. A drifting bath, slow internal equilibration or a transition invalidates this simple reduction even if a cooling curve looks smooth.

The Debye model treats a solid's low-frequency lattice vibrations as acoustic phonons in an elastic continuum. There are three polarizations per wavevector, one longitudinal and two transverse, with approximately linear dispersion relations. A cutoff is chosen to preserve the total $3N$ vibrational modes of $N$ atoms. The phonon number itself is not fixed; mode occupations follow the Bose-Einstein distribution with zero chemical potential.

For one polarization of speed $v_s$, counting wavevector states in a spherical shell gives $g_s(\omega )=V{\omega }^2/(2{\pi }^2{v}_s^3)$. Define an average Debye speed by $3/v_D^3=1/v_L^3+2/v_T^3$. The total density of states can then be writtenThe Debye wavevector is $k_D=(6{\pi }^{2}N/V{)}^{1/3}$, the Debye frequency is ${\omega }_D=v_D{k}_D$, and the Debye temperature is ${\Theta }_D=\hslash {\omega }_D/k_B$.

Omitting the temperature-independent zero-point energy, the vibrational energy and resulting Debye heat capacity areAt low temperature the upper limit can be extended to infinity. The integral is $4{\pi }^4/15$, givingOnly the long-wavelength acoustic modes are appreciably populated, explaining the Debye T-cubed test. At high temperature every mode contributes $k_B$ by equipartition theorem, so $\boxed{C_V\to 3N{k}_B}$, the Dulong-Petit law.

The model succeeds especially at low temperature, where long-wavelength elasticity is reliable. It replaces real anisotropic branches by a simple continuum and cannot reproduce detailed optical phonon branches or every intermediate-temperature feature. Metals also have an electronic heat capacity approximately proportional to $T$; defects, magnetic excitations or a phase transition can add other terms. Experiments often measure heat capacity at constant pressure, which need not equal the harmonic model's $C_V$ exactly.

Only a shell of energy width of order $k_{B}T$ near the Fermi energy can change occupation: much lower states are blocked by the Pauli exclusion principle. The thermal-shell estimate of electronic heat capacity countsTo recover the rough coefficient requested in the paper, assign each of these active particles an excitation energy of order $3k_{B}T/2$. ThenThe differentiation matters: the active shell itself widens as temperature rises. This argument explains the small linear-in-$T$ electronic heat capacity; its shell width and average excitation energy are heuristic, so it does not establish a precise numerical prefactor.

A controlled calculation from the Fermi-Dirac distribution gives a nearby but different result. At fixed particle number, the Sommerfeld expansion yieldsThus $9/2$ is a rough estimate of the exact free-electron coefficient ${\pi }^2/2\simeq 4.935$, not an exact consequence of the stated density of states. Both require $T\ll T_F$.

The quantum free electron model treats conduction Electrons as independent particles moving in an approximately constant potential, with parabolic dispersion relation $E={\hslash }^2{k}^2/(2m)$. Fermi-Dirac statistics successfully explains important features such as a Fermi energy and a small electronic heat capacity, but a constant potential gives no band gaps and cannot by itself explain why different band fillings produce metals, semiconductors and band insulators. It also misses lattice-dependent effective masses, detailed directional responses and the topology of real Fermi surfaces. Scattering and finite electrical resistivity require additional physics.

The nearly-free electron model introduces a weak periodic lattice potential. Bloch's theorem replaces individual free plane waves by Bloch states, and Fourier components of the potential mix plane waves whose wavevectors differ by a reciprocal lattice vector. Near degeneracy at a Brillouin zone boundary, this coupling opens a band gap; for a two-wave crossing its size is $2|V_G|$. Away from the crossing, energies remain close to free-particle values. Band curvature produces band effective masses, and band filling determines whether the Fermi energy lies inside an allowed band or a gap.

The improvement is the inclusion of lattice periodicity, rather than a switch from classical to quantum statistics: the quantum free electron model already uses Fermi-Dirac statistics. The strong delta barriers in the later numerical example are treated exactly by the Kronig-Penney model, rather than by assuming that their strength is perturbatively small.

In a metal, mobile Electrons transport heat as well as electric charge. A temperature gradient creates an imbalance between carriers arriving from warmer and colder regions. The resulting heat current obeys Fourier's law ${\mathbf{j}}_Q=-\kappa \nabla T$ for an isotropic material, where thermal conductivity $\kappa$ has units $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$. For a bar of area $A$ carrying steady heat flow $P$, $P=-\kappa AdT/dx$ when parasitic losses are negligible.

A kinetic estimate for the electronic contribution iswhere $C_e$ is the electronic heat capacity per volume, $v_F$ the Fermi velocity, and ${\ell }_{\kappa }=v_F{\tau }_{\kappa }$ the heat-transport mean free path. Fermi-Dirac statistics and the Pauli exclusion principle restrict thermal rearrangements to states near the Fermi surface, rather than allowing every Electron to contribute a classical heat capacity. A degenerate parabolic band givesUsing the Drude model form $\sigma =n{e}^2{\tau }_{\sigma }/m$ for electrical transport, the relation between thermal and electrical relaxation times in metals becomesFor scattering that relaxes heat and charge currents similarly, ${\tau }_{\kappa }\simeq {\tau }_{\sigma }$ and the Wiedemann-Franz law follows:Thus good electrical conductors are usually good thermal conductors. The correct quantum Lorenz number differs from a naive classical electron-gas estimate.

Scattering from impurities, boundaries and phonons limits the heat-transport path. In an impurity-dominated low-temperature regime with approximately constant relaxation time, ${\kappa }_e$ is approximately proportional to $T$. Where phonon-limited electrical resistivity is proportional to $T$ and the Wiedemann-Franz approximation is applicable, ${\kappa }_e$ can be roughly temperature independent. Inelastic processes can relax heat and charge currents differently, so the Lorenz ratio is not universal at all temperatures. Lattice phonons add a separate positive thermal contribution; the measured total is $\kappa ={\kappa }_e+{\kappa }_{\mathrm{ph}}$, and the electronic Wiedemann-Franz expression should not be applied blindly to that total.

In a metal, a partially filled energy band or overlapping energy bands supplies Electrons with nearby unoccupied states. At thermal equilibrium, oppositely directed group velocities are equally populated and the net electric current vanishes. An applied electric field distorts the occupation distribution, producing a small mean drift velocity. A filled energy band alone has no such low-field intraband response because its velocity contributions cancel over the entire Brillouin zone.

In the Drude model, let $n$ be the mobile Electron number density, let $e>0$ be the elementary charge, and let $\tau$ be the momentum relaxation time. The mean motion obeysThe steady state is ${\mathbf{v}}_d=-e\tau \mathbf{E}/m$. The current density is opposite to the Electron drift velocity:This gives microscopic Ohm's law, with electrical conductivity $\sigma$ in $\mathrm{S}{\mathrm{m}}^{-1}$ and electrical resistivity ${\rho }_{\mathrm{el}}$ in $\Omega \mathrm{m}$. For a uniform wire of length $\ell$ and cross-sectional area $A$, its electrical resistance is $R={\rho }_{\mathrm{el}}\ell /A$. The dissipated power per volume is $\mathbf{J}\cdot \mathbf{E}=\sigma E^2$, the Joule heating supplied by the external source.

The free electron model with Fermi-Dirac statistics improves the microscopic picture. The Pauli exclusion principle fills states up to the Fermi energy at low temperature. For a spherical, parabolic energy band with two spin states,Only states near the Fermi surface change their occupation in a weak-field response; the large equilibrium Fermi velocity is not the small drift velocity. The same $n{e}^2\tau /m$ formula results for this isotropic parabolic model when $\tau$ is evaluated for transport near the Fermi surface. In a crystal, the energy band group velocity, effective mass and possibly direction-dependent scattering replace the simplest bare-particle picture. The mean free path is approximately ${\ell }_{\mathrm{mfp}}=v_F\tau$.

A perfectly periodic static lattice is already included in Bloch states and does not by itself give a finite electrical resistivity. Momentum relaxation comes from defects, impurities, boundaries, and electron-phonon scattering. Independent scattering mechanisms approximately add their rates, giving Matthiessen's rule,The approximately temperature-independent ${\rho }_0$ is the residual resistivity. In ordinary simple metals, the phonon contribution is approximately proportional to temperature at high temperature; a three-dimensional acoustic-phonon transport model gives a $T^5$ contribution well below the relevant Debye temperature. Other mechanisms can dominate in particular materials or regimes, so neither power law is universal.

The classical Drude model incorrectly treats the equilibrium carrier statistics. Fermi-Dirac statistics explains the small electronic heat capacity, while electrons carrying both charge and heat give the approximate Wiedemann-Franz law ${\kappa }_e/(\sigma T)={\pi }^2{k}_B^2/(3e^2)$ for degenerate carriers with appropriate nearly elastic scattering. Both the simple electrical conductivity formula and this thermal relation depend on the assumptions about energy bands and scattering, rather than being exact for all metals.

A solid stores thermal energy in lattice vibrations and, when available, electronic and other excitations. For a nonmagnetic insulating crystal at sufficiently low temperature, the gap suppresses mobile electronic excitations and the leading bulk contribution normally comes from acoustic phonons. The three-dimensional Debye heat capacity isThe vibrational density of states scales as ${\omega }^2$, so only a fraction proportional to $T^3$ of the modes is thermally active at low temperature. The two limiting forms areThe high-temperature value is the Dulong-Petit law, or $3R$ per mole of atoms. Optical phonons add thermally activated contributions as their characteristic energies become accessible. A multiatom cell has more branches, but still three vibrational degrees of freedom per atom.

A normal metal has an additional low-temperature electronic heat capacity. The Pauli exclusion principle restricts excitations to a shell of width about $k_{B}T$ around the Fermi energy. The number of thermally active Electrons is therefore proportional to $T$, rather than the total Electron number contributing a classical heat capacity. With spin-summed electronic density of states $g(E_F)$ for the whole specimen,For a three-dimensional free-electron gas, $g(E_F)=3N_e/(2E_F)$ and $C_{\mathrm{el}}=({\pi }^2/2)N_e{k}_B(T/T_F)$. Only a small fraction of the classical electronic value is accessible when $T\ll T_F$, even though the electronic term dominates the phonon term as $T\to 0$ in a normal gapless metal.

Plotting $C/T$ against $T^2$ gives a straight line with intercept $\gamma$ and slope $\beta$ in the appropriate low-temperature range. This separates the electronic density-of-states information from the acoustic-phonon scale. In a Fermi liquid, the interacting quasiparticle density of states replaces the free value, so an enhanced effective mass enhances $\gamma$.

The qualifications matter experimentally. Superconductivity can gap electronic excitations, magnetic order adds magnetic excitations, defects can generate additional low-temperature terms, and low-dimensional or very small specimens need not obey the bulk cubic law. The Debye expression is for constant-volume heat capacity; constant-pressure heat capacity differs through thermal expansion, although the difference is often small at low temperature. Thus the normal-metal and ordinary insulating limits are useful leading behaviors, not universal formulas for every solid in every phase.

In the three-dimensional Debye model, the lattice has $3N$ normal modes for $N$ atoms. Replacing their dispersion by a linear acoustic dispersion up to a cutoff gives the Debye heat capacityAt low temperature the upper limit can be replaced by infinity. Since the integral then equals $4{\pi }^4/15$,At high temperature each harmonic mode contributes $k_B$ by the equipartition theorem, giving the Dulong-Petit lawFor one mole replace $N{k}_B$ by $R$. The curve starts with zero slope, rises as $T^3$, and levels off toward the classical limit. This is the lattice contribution; a metal's electronic heat capacity is a separate low-temperature term.

The Debye wavevector follows by counting three polarizations inside the cutoff sphere:With the supplied speed of sound treated as the effective Debye speed, ${\omega }_D=v{k}_D$ andHere $k_D\simeq 1.52\times 10^{10}{\mathrm{m}}^{-1}$. Real crystals have different longitudinal and transverse sound speeds; the one-speed estimate uses the stated value, rather than inferring those unprovided branches.

For a degenerate Fermi gas, the Pauli exclusion principle restricts thermally changed occupations to an energy shell of width $k_{B}T$ around the Fermi energy. The number of active Electrons is of order $g(E_F)k_{B}T\sim N{k}_{B}T/E_F$; each carries an excitation energy of order $k_{B}T$. Thus the thermal energy grows as $N(k_{B}T{)}^2/E_F$ and the electronic heat capacity is proportional to $T$.

More precisely, the fixed-particle-number Sommerfeld expansion givesHere $g$ is the total spin-summed density of states for the sample. For a three-dimensional parabolic band, $g(E_F)=3N/(2E_F)$, yieldingThe linear law concerns the electronic contribution; the low-temperature lattice heat capacity usually adds a $T^3$ term.

The Wiedemann-Franz law relates the electronic thermal conductivity of a normal degenerate metal to its electrical conductivity:An isotropic kinetic theory estimate gives $K_e=c_e{v}_F^2{\tau }_Q/3$. With electronic heat capacity density $c_e={\pi }^{2}n{k}_B^{2}T/(2E_F)$, $v_F^2=2E_F/m^{\ast }$ and $\sigma =n{e}^2{\tau }_{\sigma }/m^{\ast }$, the ratio is $L_0{\tau }_Q/{\tau }_{\sigma }$. Approximately elastic scattering with a smooth energy dependence gives equal heat and charge relaxation times. Inelastic electron-phonon scattering can make them different; the law also concerns the electronic contribution, not additional lattice heat conduction.^[ref]