Debye T-cubed test 2026-10-18
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 . 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 up to the Debye frequency , containing normal modes. With , the thermal vibrational internal energy becomes
The reduced Planck constant is , the Boltzmann constant is , and counts chemical atoms. For , replace the upper limit by infinity. The integral is , as in the Stefan–Boltzmann law derivation. Differentiating at fixed volume gives, per mole of chemical atoms,
where 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 against 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, , gave a linear versus plot in the nonsuperconducting state. The reported electronic intercept was ; deviations from reference heat capacities were below .[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.
Figure 1.
The Debye T-cubed test in a normal metal: lattice slope and electronic intercept
. Original constructed model with and . It represents a stipulated monatomic normal metal, not digitized gold measurements. Solid curves apply only over the chosen – fit interval, with negligible other terms and . Dashed portions illustrate extrapolation; error bars are assigned one-standard deviation uncertainties, not scatter measured in the experiment. Details.
For a worked original inference, stipulate values of equal to at . Assign independent and negligible temperature error. Equal-weight linear regression on gives
The points sit exactly on the stipulated line; these standard uncertainties come from the assigned input precision, not an estimate from zero residual scatter. With , inversion gives
At , the electronic and lattice molar heat capacities are respectively and . Their equality occurs at : even in this clean model, the total heat capacity approaches a linear, rather than cubic, law as 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 Created 2026-10-18 Updated 2026-10-20
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 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 description.
For quasiparticles with density of states including both spins, the leading electronic heat capacity is with . A narrow heavy band makes this coefficient large. The density of states and heat capacity must use consistent extensive, molar or volumetric normalization.
A Kondo lattice develops coherent low-energy quasiparticle bands below a material-dependent coherence scale. This lattice scale need not equal the single-impurity Kondo screening scale. Coherence can produce a heavy Fermi liquid with narrow bands, large electronic heat capacity and low-temperature quadratic resistivity.
Lattice qualification of an electron-only heat-capacity ratio Created 2026-10-20 Updated 2026-11-05
An electron-only comparison omits the lattice contribution to a metal's measured heat capacity. In a simple low-temperature normal metal, the same sample normalization gives
where the linear term is electronic heat capacity and the cubic term is the leading Debye model lattice contribution. Both coefficients must use compatible units and amounts of material.
If a reference gas has heat capacity , its ratio to the total metal value is
The electron-only ratio is therefore reduced when lattice phonons contribute appreciably. The crossover estimate is , where the two terms are equal. Superconductivity, magnetic excitations or other low-energy structure can add further terms. A ratio that is valid for electronic degrees of freedom alone must consequently be distinguished from a ratio of whole-sample calorimetric measurements.
Low-temperature heat capacity measurement Created 2026-10-18 Updated 2026-10-20
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 . A separate , 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 , making simple extrapolation unsafe.[ref]
For an original idealization, assume sample and platform share one temperature , the bath remains at , and both the total heat capacity and the weak-link thermal conductance are approximately constant across one small pulse. Ignore other heat transfer and any phase transition. The first law of thermodynamics becomes
Here is sample heat capacity, is calorimeter addenda, is deposited electrical power from Joule heating, and 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 beginning at bath temperature, then switched off at ,
The cooling exponential decay determines , while the heating plateau determines . Together with calorimeter addenda, they give . Cooling alone determines a ratio, not an absolute heat capacity. With finite , 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, , , and give . The asymptotic rise is , only of bath temperature. At , the rise is . Of input, only remains stored: the other has flowed to the bath. Treating this long pulse as an isolated experiment would overestimate total heat capacity by a factor of .
Suppose independent standard uncertainties are , and . First-order propagation of uncertainty gives
The uncertainty is one standard deviation under those assumptions, about of sample heat capacity. Shared thermometer calibration errors can correlate fitted and ; 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.
Pauli exclusion principle prevents most filled electronic states from changing occupation at . Only a shell of width of order around the Fermi energy participates. Since the density of states near that energy is of order per volume and energy, the fraction of active electrons is of order . Their typical change of energy is of order , so the thermal energy increment per volume scales as
Differentiating with respect to temperature gives the electronic heat capacity
The simple shell argument fixes the linear law, not its coefficient. A Sommerfeld expansion using the spin-summed free-electron density of states gives , hence , consistent with the supplied value. The result excludes strong band-structure or interaction corrections and requires a degenerate electron population.
The density of states counts one-electron states per energy interval, including both electron spin orientations. It may be a total or a per-volume quantity; must have the matching normalization. The Fermi-Dirac distribution
is the occupation probability of such a state. In the fixed-volume free electron model, the integration variable and the density of states do not depend on . The occupation depends explicitly on and implicitly through the chemical potential needed to preserve electron number. Thermal expansion or a changing electronic band structure would require additional effects outside this model.
At , only states within order of the Fermi energy change their occupation. There are order such states, each carrying an excitation energy of order . Thus and the electronic heat capacity is proportional to . Using the supplied coefficient,
The shell argument fixes the scaling, not the numerical coefficient; the controlled Sommerfeld expansion fixes that coefficient at fixed electron number. Filled states far below cannot accept arbitrary thermal excitations because of the Pauli exclusion principle.
At fixed volume and particle number, reversible heating obeys . Integrating the electronic heat capacity from zero temperature gives
With no extensive residual ground-state degeneracy, the third law of thermodynamics sets the electronic . Therefore
The entropy normalization follows that of : total, per volume or per mole. Lattice contributions are additional and are not part of this electronic result.
For degenerate carriers with approximately elastic scattering and a common transport relaxation time, and . Using and gives
This is the Wiedemann–Franz law. The same mobile electrons carry charge and heat, but Pauli statistics set the electronic heat capacity and hence the Lorenz number. The relation applies to electronic thermal conductivity, not automatically to total electrical conductivity including phonons. Inelastic electron-phonon scattering can relax energy and charge currents differently; nondegeneracy and strong energy dependence of scattering also limit the simple form.
A meaningful test combines a four-terminal electrical measurement with a steady thermal-gradient measurement on the same sample and temperature range. Contact heat resistance, parasitic heat leaks, phonon conduction and sample geometry must be controlled before interpreting a measured Lorenz ratio as an electronic departure from the law.
The Pauli exclusion principle confines thermally excited conduction electrons to an energy interval of order around the Fermi surface. Their excess energy is therefore proportional to , giving an electronic heat capacity . More precisely, . In a three-dimensional Debye model, acoustic phonons have density of states proportional to ; their thermal energy scales as , giving
per mole of atoms. This assumes a normal metal without a lower-temperature superconducting or magnetic anomaly.
For , the lattice approaches the Dulong-Petit value, so , usually dominated by . Each classical vibrational mode contributes through kinetic and potential equipartition. This is a high-temperature solid-state regime, not a claim that the solid persists to arbitrarily high temperature; anharmonicity and can matter.
The body-centred cubic cell has two atoms, so . Counting three acoustic branches gives . By the Debye-temperature definition,
This is an effective Debye sound speed; separate longitudinal and transverse speeds are not determined. For isotropic branches, .
The free electron model replaces the microscopic ionic potential and interacting electronic motion by independent mobile electrons in a uniform neutralizing background, with approximately parabolic dispersion . Its quantum version fills states according to Fermi-Dirac statistics and the Pauli exclusion principle. Scattering is represented by a momentum relaxation time , rather than derived from the model's ideal free trajectories. This model explains a Fermi energy much larger than room-temperature thermal energy, an electronic heat capacity proportional to temperature at low temperature, and a Drude electrical conductivity. With appropriate scattering assumptions it also yields the Wiedemann-Franz law. It cannot by itself explain band gaps, detailed band anisotropy, strong correlations, or the microscopic temperature dependence of .
The inertial term is mass times the rate of change of mean drift velocity. The electric force is , taking the printed as the signed electron charge, and describes loss of directed momentum to scattering. electrons retain large microscopic velocities; this damping is of their mean directed motion. In a steady DC field,
Both charge factors are negative for electrons, so the conductivity is positive. If instead denotes a positive charge magnitude, the force and current equations each require a minus sign and give the same result.
For AC response choose physical fields as real parts of phasors proportional to . Substitution gives
Thus
The real part decreases from to zero. The imaginary part starts at zero, reaches at , and falls as . The real part falls as at high frequency. With an convention, only the sign of the imaginary part changes. Positive gives positive cycle-averaged dissipation .
Figure 1.
Drude conductivity and the dispersion and density of states of a monatomic harmonic chain
. Details.
The supplied DC conductivity determines the transport time:
Inertia matters when , giving , or . Quoting the former number as hertz would introduce an erroneous factor of .
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 vibrational modes of atoms. The phonon number itself is not fixed; mode occupations follow the Bose-Einstein distribution with zero chemical potential.
For one polarization of speed , counting wavevector states in a spherical shell gives . Define an average Debye speed by . The total density of states can then be written
The Debye wavevector is , the Debye frequency is , and the Debye temperature is .
Omitting the temperature-independent zero-point energy, the vibrational energy and resulting Debye heat capacity are
At low temperature the upper limit can be extended to infinity. The integral is , giving
Only the long-wavelength acoustic modes are appreciably populated, explaining the Debye T-cubed test. At high temperature every mode contributes by equipartition theorem, so , 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 ; 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 exactly.
Only a shell of energy width of order 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 counts
To recover the rough coefficient requested in the paper, assign each of these active particles an excitation energy of order . Then
The differentiation matters: the active shell itself widens as temperature rises. This argument explains the small linear-in- 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 yields
Thus is a rough estimate of the exact free-electron coefficient , not an exact consequence of the stated density of states. Both require .
The quantum free electron model treats conduction electrons as independent particles moving in an approximately constant potential, with parabolic dispersion relation . 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 . 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 for an isotropic material, where thermal conductivity has units . For a bar of area carrying steady heat flow , when parasitic losses are negligible.
A kinetic estimate for the electronic contribution is
where is the electronic heat capacity per volume, the Fermi velocity, and 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 gives
Using the Drude model form for electrical transport, the relation between thermal and electrical relaxation times in metals becomes
For scattering that relaxes heat and charge currents similarly, 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, is approximately proportional to . Where phonon-limited electrical resistivity is proportional to and the Wiedemann-Franz approximation is applicable, 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 , 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 be the mobile electron number density, let be the elementary charge, and let be the momentum relaxation time. The mean motion obeys
The steady state is . The current density is opposite to the electron drift velocity:
This gives microscopic Ohm's law, with electrical conductivity in and electrical resistivity in . For a uniform wire of length and cross-sectional area , its electrical resistance is . The dissipated power per volume is , 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 formula results for this isotropic parabolic model when 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 .
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 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 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 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 straight line in against means . Reading its intercept and slope gives approximately and ; the graphical precision is only a few percent. Hence the molar heat capacity is approximately
The linear electronic heat capacity comes from conduction electrons within energy of order of the Fermi energy; the Pauli exclusion principle prevents deeply occupied states from contributing freely. The cubic lattice term is the low-temperature Debye heat capacity of three-dimensional acoustic phonons, whose low-frequency density of states is proportional to frequency squared. These are molar coefficients, so the temperature powers carry the extra kelvin factors indicated.
Angle-resolved photoemission spectroscopy maps occupied energy bands: ultraviolet or x-ray photons eject electrons, and an energy- and angle-resolving analyzer measures and emission direction. The photoelectric effect gives the binding energy ; the parallel momentum follows from , up to surface reciprocal vectors. Clean surfaces and controlled photon energy are needed, and an occupied-state measurement alone does not map an empty conduction band.
Scanning tunneling spectroscopy instead brings a conducting tip close to a surface and modulates the bias while recording the differential tunnel current. At low temperature with a smooth tip density of states and tunnelling matrix element, probes the local density of states near energy . The signal is thermally broadened and depends on tunnelling selection, so it is not an unqualified map of the whole bulk band structure.
Bulk probes complement these surface measurements. In the de Haas-van Alphen effect, a low-temperature magnetization or torque measurement versus magnetic field gives quantum-oscillation frequencies , measuring extremal Fermi surface sections; temperature damping gives a cyclotron effective mass. The linear electronic heat capacity gives the density of low-energy excitations at the Fermi level. Optical absorption measures transition thresholds and matrix-element-weighted joint densities, rather than a unique full dispersion by itself.
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 is
The vibrational density of states scales as , so only a fraction proportional to of the modes is thermally active at low temperature. The two limiting forms are
The high-temperature value is the Dulong-Petit law, or 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 around the Fermi energy. The number of thermally active electrons is therefore proportional to , rather than the total electron number contributing a classical heat capacity. With spin-summed electronic density of states for the whole specimen,
For a three-dimensional free-electron gas, and . Only a small fraction of the classical electronic value is accessible when , even though the electronic term dominates the phonon term as in a normal gapless metal.
Figure 1.
Debye lattice heat capacity and the low-temperature electronic intercept in a normal metal
. Details.
Plotting against gives a straight line with intercept and slope 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 .
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.