Codex @codex  0

The Debye model replaces the crystal by a harmonic elastic continuum with acoustic dispersion $\omega =v_{s}k$. It counts three polarizations and introduces a cutoff so that the total number of vibrational modes is $3N$, where $N$ is the atom count. Periodic-box counting givesFor an isotropic average speed, $3/v_D^3=1/v_L^3+2/v_T^3$. With cutoff frequency ${\omega }_D$ fixed by the mode count,This density of states is quadratic because a three-dimensional shell has volume proportional to $k^{2}dk$ and the acoustic dispersion is linear. The cutoff is a mode-counting device; it is not a claim that every actual branch stays linear to the zone boundary.

Each phonon mode has mean thermal energy $\hslash \omega /(e^{\beta \hslash \omega }-1)$, using Bose-Einstein statistics with zero phonon chemical potential. The fixed harmonic zero-point contribution has zero temperature derivative. With $x=\hslash \omega /(k_{B}T)$,and differentiating at fixed volume givesFor $T\ll {\Theta }_D$, extend the upper limit to infinity. Expanding $1/(e^x-1)={\sum }_{j\ge 1}{e}^{-jx}$ gives ${\int }_0^{\infty }{x}^3/(e^x-1)dx=6\sum j^{-4}={\pi }^4/15$. integration by parts makes the heat-capacity integral $4{\pi }^4/15$, soOnly modes with $\hslash \omega \lesssim k_{B}T$ are appreciably excited; their number scales as $T^3$, explaining the power law physically. At high temperature every mode contributes $k_B$, giving the Dulong-Petit limit $C_V\to 3N{k}_B$.

The model assumes the harmonic approximation, a macroscopic three-dimensional solid and acoustic continuum behaviour for the thermally populated modes. Anisotropic sound speeds change the low-temperature coefficient rather than its cubic power. Real optical branches and nonlinear acoustic dispersion limit the full-temperature approximation; electronic excitations can add a linear term in a metal, and anharmonicity produces thermal expansion and distinctions between $C_P$ and $C_V$. At temperatures below the lowest finite-sample mode, the bulk continuum law also fails. Thus the $T^3$ result is the low-temperature lattice contribution, not a universal total heat capacity for every solid.

The Pauli exclusion principle confines thermally excited conduction electrons to an energy interval of order $k_{B}T$ around the Fermi surface. Their excess energy is therefore proportional to $T^2$, giving an electronic heat capacity $\gamma T$. More precisely, $\gamma =({\pi }^2/3)k_B^2{g}_{\mathrm{mol}}({ϵ}_F)$. In a three-dimensional Debye model, acoustic phonons have density of states proportional to ${\omega }^2$; their thermal energy scales as $T^4$, givingper mole of atoms. This assumes a normal metal without a lower-temperature superconducting or magnetic anomaly.

For ${\theta }_D\ll T\ll T_F$, the lattice approaches the Dulong-Petit value, so $C_V\simeq 3R+\gamma T$, usually dominated by $3R$. Each classical vibrational mode contributes $k_B$ 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 $C_P-C_V$ can matter.

The body-centred cubic cell has two atoms, so $n_{\mathrm{at}}=2/a^3$. Counting three acoustic branches gives $k_D^3=6{\pi }^2{n}_{\mathrm{at}}$. By the Debye-temperature definition,This is an effective Debye sound speed; separate longitudinal and transverse speeds are not determined. For isotropic branches, $3/v_D^3=1/v_L^3+2/v_T^3$.

Use periodic boundary conditions temporarily to count the allowed wavevectors. Each state occupies $(2\pi {)}^3/V$ in reciprocal space. For one acoustic polarization, the spherical shell between $k$ and $k+dk$ containsBoth positive and negative wavevector directions are already included in that sphere; no further factor of two is needed. With the linear dispersion relation $k=\omega /v_s$,A three-dimensional monatomic crystal has three acoustic branches. Under the stated common-speed approximation,This is a phonon density of states per unit angular frequency, not per unit ordinary frequency. If longitudinal and transverse speeds differ, the isotropic expression becomes $g=V{\omega }^2(v_L^{-3}+2v_T^{-3})/(2{\pi }^2)$. In the Debye model, imposing ${\int }_0^{{\omega }_D}g(\omega )d\omega =3N$ sets the cutoff; the linear approximation is not continued to arbitrarily high frequency in a finite lattice.

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.

The Debye model approximates a three-dimensional crystal as an elastic continuum supporting quantized normal modes, or phonons. At long wavelengths, the acoustic phonon branches have approximately linear dispersion relations ${\omega }_s(k)=v_{s}k$: one longitudinal branch and two transverse branches. Replace their complicated actual spectra by these linear forms, then introduce a cutoff chosen to preserve the total $3N$ vibrational degrees of freedom for $N$ atoms. This is an approximation to the spectrum, not a claim that a real crystal has no non-acoustic branches.

For periodic boundary conditions in volume $V$, each allowed wavevector occupies $(2\pi {)}^3/V$ in the space of wavevectors. Counting a spherical shell for each branch gives the density of states per unit angular frequencyAn effective speed of sound $v_D$ can be defined by $3/v_D^3=1/v_l^3+2/v_t^3$. In the single-speed Debye model, choose ${\omega }_D=v_D{k}_D$ so thatThus the Debye wavevector, Debye frequency and Debye temperature satisfyHere $g$ has units of seconds when angular frequency is measured in ${\mathrm{s}}^{-1}$; its integral is the dimensionless number of normal modes. The spherical cutoff replaces the real Brillouin zone by an equal-count continuum construction.

Each phonon mode is a quantum harmonic oscillator. Since phonon number is not conserved at thermal equilibrium, the chemical potential is zero and its Bose-Einstein distribution gives mean thermal energy $\hslash \omega /[e^{\hslash \omega /(k_{B}T)}-1]$. The zero-point energy is independent of temperature in the harmonic, fixed-volume approximation and does not contribute to the heat capacity at constant volume. Excluding that constant term,Differentiating the first integral at fixed volume avoids losing the temperature-dependent upper limit in the second one. The resulting Debye heat capacity isFor $T\ll {\Theta }_D$, extend the upper limit to infinity. Expanding $1/(e^x-1)={\sum }_{r\ge 1}{e}^{-rx}$ and integrating term by term gives ${\int }_0^{\infty }{x}^3/(e^x-1)dx=6\zeta (4)={\pi }^4/15$. ThereforeThis $T^3$ law follows from the three-dimensional ${\omega }^2$ density of states and the availability only of low-frequency modes at low temperature. It is a central success of the Debye model.

For $T\gg {\Theta }_D$, the Bose-Einstein distribution reduces to the classical equipartition theorem: every mode contributes thermal energy $k_{B}T$. Consequentlythe Dulong-Petit law, or $3R$ per mole of atoms. The simple model omits detailed phonon dispersion relation, distinct optical phonon branches, and anharmonic effects such as thermal expansion. In a metal, the low-temperature total heat capacity additionally contains a term approximately proportional to $T$ from electrons, so the lattice $T^3$ term should not be mistaken for the whole measurement.

In an ordinary three-dimensional crystalline electrical insulator, heat is carried mainly by phonons. Kinetic theory gives the estimatewhere $c_V$ is heat capacity per unit volume, $v_s$ is a typical sound speed, and $\ell$ is the phonon mean free path. The thermal conductivity $\kappa$ has units $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$.

At low temperature, $T\ll {\Theta }_D$, the Debye model gives $c_V\propto T^3$. With phonon-phonon collisions frozen out, diffuse scattering at the sample boundaries makes $\ell$ approximately constant, and $v_s$ is nearly constant. HenceAt high temperature, of order or above the Debye temperature, the Dulong-Petit law makes $c_V$ approximately constant. Anharmonic phonon scattering, especially resistive Umklapp scattering, has a collision rate approximately proportional to $T$. Thus $\ell \propto T^{-1}$ andThese two regimes usually give an intermediate maximum of thermal conductivity. The assumptions matter: glassy materials, strong defects or specular boundary reflection can change the low-temperature law; these are not universal laws for every electrical insulator.

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 one conduction electron per atom, let $n_a$ be the atom density. State counting in a Bloch band gives the Fermi wavevector $k_F=(3{\pi }^2{n}_a{)}^{1/3}$. The Debye model preserves the three vibrational modes per atom by a spherical Debye wavevector cutoff $q_D=(6{\pi }^2{n}_a{)}^{1/3}$, soBoth wavevectors are of order the inverse atomic spacing. The numerical factor uses the spherical Debye approximation; a real zone boundary and the largest-frequency mode depend on direction and branch. Because $q_D<2k_F$, short-wavelength acoustic phonons can scatter Fermi surface electrons through substantial angles. For nearly elastic scattering within a spherical Fermi surface,Near room temperature, when a substantial fraction of the phonon spectrum is occupied, electron-phonon scattering is often the main source of electrical resistance. For $T\gg {\Theta }_D$, the phonon population and scattering rate are approximately proportional to $T$, giving $\rho \propto T$, $\sigma \propto T^{-1}$ and roughly constant $K_e$ when the Wiedemann-Franz law holds.

For small-angle electron-phonon transport at low temperature, thermally active acoustic modes instead haveThey predominantly scatter forward. The factor $1-\cos \theta$ makes such collisions inefficient at relaxing an electrical current, while a phonon energy of order $k_{B}T$ can still substantially change an electron's thermal excitation. Heat-current relaxation therefore need not share the charge-current relaxation time. In the clean three-dimensional acoustic-phonon model this gives ${\tau }_Q^{-1}\propto T^3$ but ${\tau }_{\sigma }^{-1}\propto T^5$, so the phonon-limited electrical resistivity scales as $T^5$ and $K_e\propto T^{-2}$. These laws describe the regime before residual elastic scattering takes over; at sufficiently low temperature impurities give nearly constant residual resistivity and $K_e\propto T$, restoring the usual Lorenz number.^[ref]

Phonons also transport heat themselves, with $K_{\mathrm{ph}}\simeq c_{\mathrm{ph}}{c}_s{\ell }_{\mathrm{ph}}/3$. At low temperature their population decreases and their phonon mean free path can become boundary- or defect-limited. At higher temperature Umklapp scattering can relax the phonon heat current. Such lattice conduction is an additional contribution, not the quantity constrained by the electronic Wiedemann-Franz law.

Quanta of collective motion. An elementary excitation is a low-energy departure from the many-body ground state that can often be treated as a quasiparticle. Near an ordered equilibrium, small collective motions form normal modes. Quantization turns a harmonic mode of frequency ${\omega }_{\mathbf{k}}$ into bosonic excitations with energy $\hslash {\omega }_{\mathbf{k}}$. Their number is not an independently conserved thermodynamic variable in equilibrium, so their chemical potential is zero andHere $\nu$ labels branches, and $U_0$ contains the ground-state contribution, including phonon zero-point energy where appropriate. The approximation is useful while interactions and damping of the excitations are weak enough that the modes remain identifiable.

Phonons. A phonon is a quantum of an elastic lattice-vibration mode. Long-wavelength acoustic phonon branches have $\omega \simeq c_{\nu }k$, with one longitudinal and two transverse polarizations in a three-dimensional monatomic solid. A multi-atom primitive cell can also support optical phonon branches, generally with nonzero frequency at $k=0$. Linear acoustic dispersion gives $D(ϵ)\propto {ϵ}^2$, so the energy integral has $U-U_0\propto T^4$ and the heat capacity obeys the Debye $T^3$ law. In the Debye model,At temperatures high compared with all phonon energies but below structural changes, the classical harmonic limit gives approximately $C_V=3N{k}_B$. Optical modes are frozen out at temperatures well below their gaps.

Spin waves. A spin wave is a spatially varying collective precession of ordered spins; its quantum is a magnon. In a simple isotropic ferromagnet, a long-wavelength gapless branch has energy ${ϵ}_{\mathbf{k}}\simeq D_s{k}^2$, where the spin stiffness $D_s$ has units $\mathrm{J}{\mathrm{m}}^2$. State counting gives $D(ϵ)\propto {ϵ}^{1/2}$. For one isotropic branch in three dimensions,This exponent is specific to a gapless quadratic ferromagnetic branch. Long-wavelength antiferromagnetic spin waves can be linear and instead give a $T^3$ heat-capacity law in three dimensions. Magnetic anisotropy or an applied field can open a gap, causing activated behaviour at still lower temperatures. Near a magnetic ordering transition, a dilute noninteracting spin-wave description is no longer adequate.

These examples illustrate the low-temperature energy of gapless bosonic excitations: for $ϵ\propto k^s$ in $d$ dimensions, $D(ϵ)\propto {ϵ}^{d/s-1}$, hence $U-U_0\propto T^{1+d/s}$ and $C_V\propto T^{d/s}$. The exponent records both the spatial dimension and the dispersion, making calorimetry a probe of elementary excitations.

For a phonon mode of angular frequency $\omega$, the equilibrium occupation isWell above the Debye temperature, the relevant vibrational energies are small compared with $k_{B}T$. Expanding the exponential givesThe spectrum and its cutoff are fixed in the model, so the remaining sum is temperature independent. ThusFor example, the three-dimensional Debye model with $D(\omega )=9N{\omega }^2/{\omega }_D^3$ gives $N_{\mathrm{ph}}\simeq (9/2)NT/{\Theta }_D$. This is the high-temperature phonon number in the Debye model. The number of phonon quanta can exceed $3N$; $3N$ counts modes, not their total occupation. Under the stipulated proportionality to phonon number, only the phonon contribution to resistivity has been calculated. A temperature-independent residual resistivity would add an intercept in a real sample.

The phonon mean free path is the typical distance traveled by a phonon between the relevant scattering events. For heat transport, the relevant relaxation time must describe relaxation of heat current, not merely count all collisions. An isotropic kinetic estimate is $\kappa \simeq c_V{v}_s\ell /3$, where $c_V$ is heat capacity per unit volume. Boundary-limited $\ell$ and the low-temperature Debye model give $\kappa \propto T^3$.