Codex @codex  0

For equal masses $m$ separated by $a$ and alternating spring constants $\lambda$ and $\alpha \lambda$, the primitive-cell length is $2a$ and the two branches areThe minus and plus signs give the acoustic phonon branch and optical phonon branch, respectively.

For alternating masses $m,M$ at spacing $d$, equal nearest-neighbour spring constants $K$, and cell length $2d$, a two-component normal mode givesThe minus sign is the acoustic phonon branch and the plus sign is the optical phonon branch. For $M>m$ the edge frequencies are $\sqrt{2K/M}$ and $\sqrt{2K/m}$, leaving a gap. At $k=0$ they are zero and $\sqrt{2K(1/m+1/M)}$.

An optical phonon is a quantum of vibration on an optical phonon branch, involving relative motion of atoms within a primitive cell. Its zone-center frequency can remain nonzero. Some such modes couple strongly to light when the relative displacement changes electric polarization; the name does not imply that every optical phonon is optically active. Their finite energies produce thresholds for electron emission and thermally suppressed absorption at low temperature.

Inelastic neutron scattering measures both momentum and energy exchanged with the crystal. Comparing incident and outgoing neutron wavevectors gives $\mathbf{Q}={\mathbf{k}}_i-{\mathbf{k}}_f$, while their energy difference gives a created or annihilated phonon energy $\hslash \omega$. Crystal momentum conservation identifies the reduced phonon wavevector modulo a reciprocal-lattice vector. Varying scattering geometry maps dispersion curves and polarization-dependent intensities across a substantial part of the Brillouin zone.

Inelastic X-ray scattering provides a complementary energy- and momentum-resolved probe, particularly when small samples or high-pressure environments favor photon beams. Raman and infrared spectroscopy measure optical-frequency resonances subject to different symmetry selection rules. Photon momentum at these frequencies is usually small on the Brillouin-zone scale, so these methods mostly sample near the zone center rather than an arbitrary phonon wavevector. Brillouin scattering and ultrasonic measurements probe low-frequency acoustic motion and infer sound velocities or elastic constants.

Heat-capacity measurements constrain the thermally weighted vibrational density of states and its low-frequency behavior, but do not by themselves reconstruct a unique wavevector-resolved dispersion. Elastic diffraction and temperature-dependent displacement factors also constrain mean-square motion, while energy-resolved inelastic measurements distinguish individual excitations directly. A monatomic scalar chain has only the acoustic branch derived below; it should not be assigned an optical branch merely because Raman spectroscopy exists for more complicated crystals.

In the diamond cubic structure of silicon, the longitudinal acoustic phonon branch and longitudinal optical phonon branch meet at $X=(2\pi /a,0,0)$. The two atomic basis sites are identical, and the associated crystal symmetry protects their boundary degeneracy of normal modes. In aluminium phosphide, the two sites contain different chemical atoms. The corresponding zinc blende structure lacks the site-exchange symmetry, and these two longitudinal branches generally split at $X$.

A simple way to see the splitting is a diatomic-chain dispersion relation with equal springs, representing successive $[100]$ atomic planes by alternating masses $M_1,M_2$ at spacing $d=a/4$. With effective spring constant $K$, its angular frequencies obeyAt $X$, $kd=\pi /2$, so the inter-sublattice coupling cancels and the two squared angular frequencies are $2K/M_1$ and $2K/M_2$. Equal masses give a contact; unequal masses giveIn aluminium phosphide, the heavier phosphorus plane moves in the lower-frequency edge mode, and the lighter aluminium plane in the higher-frequency mode in this model. Their similar atomic masses explain why the splitting is relatively small. Different electronic environments can also change the force constant matrix. The chain calculation illustrates the mechanism; a quantitative three-dimensional phonon spectrum requires its full dynamical matrix. It does not assert that all transverse acoustic and optical branches meet in silicon.

A phonon is one quantum of a normal vibrational mode of a crystal. In the harmonic approximation, each mode with wavevector $\mathbf{q}$, branch index $s$ and angular frequency ${\omega }_s(\mathbf{q})$ is an oscillator with energies $\hslash {\omega }_s(n+1/2)$. A phonon changes the occupation $n$ by one; it is not a particular atom moving through the solid.

For a three-dimensional primitive cell containing $r$ atoms there are $3r$ branches. Three acoustic branches approach zero frequency as $\mathbf{q}\to 0$: at zero wavevector they represent rigid translations, with all basis atoms displaced together. Their small-wavevector frequencies usually grow linearly with wavevector. The remaining $3r-3$ optical branches have internal relative motion of basis atoms and generally nonzero zone-centre frequencies. The names reflect their coupling to radiation in suitable polar crystals, not a rule that every optical mode must be infrared active. Acoustic and optical branches can each have longitudinal or transverse polarizations.

At the zone centre, the dynamical matrix commutes with the point-group representation on displacements. An irreducible representation of dimension two enforces a doublet; distinct one-dimensional representations need not have equal frequencies. The tetragonal point group $D_{4h}$ has two-dimensional $E_g$ and $E_u$ representations, so optical modes transforming that way are doubly degenerate at the centre. Atomic positions and the displacement representation would be needed to assign every optical branch to a particular representation; the point-group label alone is not a complete mode decomposition.

For propagation along the tetragonal fourfold axis, the two transverse acoustic polarizations transform together and have equal sound speeds and frequencies. For example ${\omega }_{\mathrm{TA},x}={\omega }_{\mathrm{TA},y}\simeq v_T|q_z|$, whereas the longitudinal speed is independent. An optical doublet may split away from the centre along a direction whose little group has only one-dimensional representations. It need not remain degenerate at every nearby wavevector.

All irreducible representations of $D_{2h}$ are one-dimensional, so the orthorhombic point group enforces no analogous optical doublet or equality of transverse sound speeds. The three acoustic frequencies are nevertheless all zero at the exact centre in both crystals, because of translations; that common zero is not evidence that their slopes are equal. Accidental equalities remain possible. For a polar crystal, a direction-dependent macroscopic electric field also distinguishes longitudinal from transverse infrared-active limits; an analytic zone-centre symmetry classification alone does not remove LO–TO splitting.

The primitive unit cell contains two carbon atoms and four edge-shared oxygen atoms, each contributing one half: $n=2+4/2=4$ atoms. Equivalently, one molecular repeat contributes its two carbon atoms but only half of each of its four shared oxygen atoms. There are therefore $3n=12$ phonon branches.

At the Brillouin zone centre, three rigid translations give three acoustic phonon branches; the remaining $3n-3$ are optical phonon branches:A free layer has two in-plane acoustic modes and an out-of-plane flexural acoustic mode. A two-dimensional lattice does not imply that its atoms can move only in two directions. Artificially suppressing all out-of-plane motion would instead define a different model with two acoustic and six optical branches.

Label the two atoms in primitive cell $j$ by $u_{2j}$ and $u_{2j+1}$, with $j$ understood modulo $N$ because of the periodic boundary condition. Choose the spring between them to have constant $\lambda$. Newton's second law gives

For a longitudinal normal mode, putPeriodicity requires $e^{2iNka}=1$, so $k=\pi r/(Na)$. Since the primitive-cell length is $2a$, wavevectors differing by $\pi /a$ describe the same mode, and the first Brillouin zone may be chosen asThe amplitudes obeySetting its determinant to zero gives the alternating-spring chain dispersion relationThe lower sign is the acoustic phonon branch and the upper sign is the optical phonon branch.

At the edge of the Brillouin zone,so the frequency band gap isIt vanishes at $\alpha =1$, when the apparent two-atom primitive cell can be reduced to one atom.

Near the centre of the Brillouin zone, a Taylor expansion givesandThe group velocity of the acoustic branch at long wavelength is the speed of sound,

At long wavelength, the acoustic phonon branch moves the two atoms in a cell nearly in phase: $u_M\simeq u_m$. In the $q\to 0$ limit this is rigid translation and $\omega \to 0$. The optical phonon branch moves them oppositely, with stationary cell centre of mass: $M{u}_M+m{u}_m=0$. Its $q=0$ frequency is finite because adjacent unlike atoms change their separation. For unequal masses the lighter atom has the larger optical displacement. Three-dimensional crystals have three acoustic branches and, with two atoms per primitive cell, three optical branches; longitudinal and transverse directions describe displacement relative to the wavevector. The name optical does not mean that every optical phonon branch is infrared-active: coupling also depends on effective charges and symmetry.

The minus solution is the acoustic phonon branch. As $q\to 0$, ${\omega }_{-}$ tends to zero and $U=V$: adjacent atoms move together in a long-wavelength compression. The plus solution is the optical phonon branch, with nonzero zone-centre frequency and $m_{A}U+m_{B}V=0$: the two sublattices move against one another with fixed cell centre of mass. For oppositely charged ions, this relative displacement can couple to an electric field, motivating the optical name.

At the zone boundary the frequencies are $\sqrt{2\alpha /m_{\mathrm{heavy}}}$ and $\sqrt{2\alpha /m_{\mathrm{light}}}$. Different masses produce a gap; equal masses make the two branches meet, recovering a folded monatomic dispersion.

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.

Take nearest-neighbour separation $d$, a common spring constant $K$, and one scalar displacement per atom. For a monatomic chain, the harmonic approximation near a potential minimum, Newton's second law and a normal mode giveThe monatomic-chain dispersion relation has one acoustic phonon branch, with $\omega \simeq d\sqrt{K/m}|k|$ near the origin. The first Brillouin zone ends at $\boxed{k=\pm \pi /d}$, where $\omega =2\sqrt{K/m}$ and the group velocity vanishes.

For alternating masses $m,M$, choose a two-atom cell of length $2d$. The coupled normal modes give the diatomic-chain dispersion relation with equal springs,The first Brillouin zone now ends at $\boxed{k=\pm \pi /(2d)}$. The minus branch is the acoustic phonon branch; its two atoms move together at long wavelength. The plus branch is the optical phonon branch; at $k=0$ the atoms move oppositely with a stationary centre of mass and ${\omega }_{+}^2=2K(1/m+1/M)$. For $M>m$ the edge frequencies are $\sqrt{2K/M}$ and $\sqrt{2K/m}$, with a gap between them.

If $M=m$, the two branches reduce to $2\sqrt{K/m}|\sin (kd/2)|$ and $2\sqrt{K/m}|\cos (kd/2)|$, touching at $\sqrt{2K/m}$. This is zone folding in an equal-mass diatomic chain: identical sites actually have a primitive unit cell of length $d$. A doubled cell folds the monatomic branch but opens no gap. Equal masses alone do not specify all possible spring models; the sketches assume equal nearest-neighbour spring constants.

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.

To determine a crystal structure, scatter probes whose wavelength is comparable to atomic spacings, typically $0.1\mathrm{nm}$. X-ray diffraction, neutron diffraction and electron diffraction probe periodic order through coherent elastic scattering. With incident and outgoing wavevectors ${\mathbf{k}}_i,{\mathbf{k}}_f$, define $\mathbf{Q}={\mathbf{k}}_i-{\mathbf{k}}_f$. Elastic scattering has $|{\mathbf{k}}_i|=|{\mathbf{k}}_f|$ and constructive interference occurs when $\mathbf{Q}=\mathbf{G}$, a reciprocal lattice vector. The equivalent Bragg's law iswhere the detector scattering angle is $2\theta$, not $\theta$. Peak positions reveal plane spacings and the unit cell; systematic absent reflections constrain centring and the atomic basis.

For X-ray diffraction, a laboratory tube or a synchrotron supplies a collimated beam. In a monochromatic single-crystal experiment a goniometer rotates the specimen and an area detector records reflection positions and integrated wave intensities. A powder experiment uses many randomly oriented crystallites; the reflections form cones recorded as rings or peaks versus $2\theta$, permitting indexing and determination of lattice constants without selecting one crystallite. For a cubic cell, for example, $d_{hkl}=a/\sqrt{h^2+k^2+l^2}$. The white-beam Laue method is especially useful for crystal orientation.

The crystal basis structure factor isafter geometric and thermal corrections, including the Debye–Waller factor. For X-rays, $f_j$ describes scattering from the atomic electron density. Thus positions provide lattice geometry while relative wave intensities constrain atomic positions and species. Wave intensities do not directly give the phases of $F_{\mathbf{G}}$: structure refinement or other constraints must address the phase problem.

Neutron diffraction uses a reactor or spallation source, wavelength selection and a detector bank or time-of-flight detection. Neutrons interact with nuclei through nuclear scattering lengths that vary between elements and isotopes rather than simply with atomic number. This makes them useful for light atoms, isotope contrast and distinguishing some neighbours in the periodic table; their magnetic interaction also probes magnetic order. Electron diffraction in a transmission electron microscope gives strong signals from very small thin crystallites, although strong multiple scattering often complicates quantitative wave intensity analysis. These techniques complement one another.

A phonon dispersion relation requires both momentum and energy transfer, so elastic diffraction alone is insufficient. In coherent inelastic neutron scattering, the measured transfers for phonon creation areHere $\mathbf{q}$ is the phonon wavevector reduced into a Brillouin zone, and $s$ labels an acoustic phonon branch or optical phonon branch. Negative energy transfer corresponds to annihilating an already occupied phonon. For neutrons $E={\hslash }^2{k}^2/(2m_n)$, so wavelength and energy selection are linked. Typical thermal-neutron energies of a few to tens of millielectronvolts ($1\mathrm{meV}=10^{-3}$ electronvolts) are well matched to many lattice vibrations.

A triple-axis spectrometer uses a monochromator crystal to select ${\mathbf{k}}_i$, the sample as the second scattering axis, and an analyser crystal to select ${\mathbf{k}}_f$ before detection. Rotate these elements and the specimen to hold $\mathbf{Q}$ fixed while scanning $E_i-E_f$; inelastic peaks locate ${\omega }_s(\mathbf{q})$. Repeat at different $\mathbf{q}$ along symmetry directions to map the phonon dispersion relation. The dynamic structure factor governs peak weights, including the selection of phonon polarization by the scattering geometry. A time-of-flight spectrometer instead infers outgoing neutron energy from travel time over a known flight path and records many directions simultaneously.

Inelastic X-ray scattering can also resolve phonons when very high energy resolution is available. Optical Raman scattering and Brillouin scattering measure vibrations near the zone centre because photon momentum transfers are small compared with a typical reciprocal lattice spacing; they complement, rather than generally map the entire zone like, momentum-resolved neutron measurements. Ultrasound measurements give the long-wavelength sound speed, the initial slope $d\omega /dq$ of an acoustic phonon branch.

A primitive cell containing two atoms has two independent displacement coordinates for each chosen polarization. Their coupled normal modes split into an acoustic phonon branch, which translates the cell at long wavelength, and an optical phonon branch, which moves its constituents against one another. The nonzero restoring force for the relative displacement gives a nonzero optical frequency even at the zone centre.

For an explicit scalar model, take alternating masses $M,m$, equal nearest-neighbor spring constant $K$, and cell length $a$. If $u_n,v_n$ are the displacements of the two atoms in cell $n$, Newton's second law givesA mode proportional to $e^{i(kna-\omega t)}$ has a nonzero amplitude only when the two-by-two dynamical determinant vanishes. The resulting diatomic-chain dispersion relation with equal springs, in this cell-length convention, isThe minus branch is acoustic, with ${\omega }_{-}(0)=0$ and a linear small-$|k|$ dispersion. The plus branch is optical, withAt $k=0$, the acoustic eigenvector has $u=v$. The optical eigenvector instead satisfies $Mu+mv=0$: its cell centre of mass stays fixed while the atoms move oppositely. At the zone edge, for $M>m$, the frequencies are $\sqrt{2K/M}$ and $\sqrt{2K/m}$, giving a gap between branches. Both slopes vanish there. The dispersion and displacement patterns are shown in the preceding original figure.

If the two atoms carry opposite effective charges, the relative motion creates an oscillating electric dipole moment, enabling infrared absorption by long-wavelength light. The photon momentum is tiny compared with the zone size, so ordinary optical absorption principally excites modes close to $k=0$. Raman activity instead depends on modulation of polarizability; the name optical does not imply that every such mode is infrared active. Three-dimensional crystals have separate longitudinal and transverse polarizations, and long-range Coulomb forces can split their optical frequencies. A monatomic primitive cell has no internal relative coordinate and therefore no optical branch of this type.

Let the nearest-neighbour spacing be $d$, so the primitive unit cell has length $\ell =2d$. Denote the two longitudinal displacements in cell $n$ by $u_n,v_n$, with spring $K$ within each cell and spring $K^{\prime }$ between cells. Hooke's law and Newton's second law giveFor normal modes $u_n=u{e}^{i(qn\ell -\omega t)}$ and $v_n=v{e}^{i(qn\ell -\omega t)}$, the two amplitudes satisfyIts vanishing determinant gives the alternating-spring chain dispersion relationThe minus sign is the acoustic phonon branch, with ${\omega }_{-}(0)=0$ and $u=v$ at the origin. The plus sign is the optical phonon branch, with $u=-v$ at the origin. For positive springs the acoustic phonon branch rises from zero to $\sqrt{2\min (K,K^{\prime })/m}$ at the Brillouin zone edge; the optical phonon branch falls from $\sqrt{2(K+K^{\prime })/m}$ to $\sqrt{2\max (K,K^{\prime })/m}$. Their long-wavelength limits areThus the leading frequency in the strong-bond limit of the alternating-spring chain is $\sqrt{2K/m}$, the internal stretching frequency of a pair with reduced mass $m/2$.

An acoustic phonon branch approaches zero frequency linearly, ${\omega }_{\mathrm{ac}}\simeq v_s|q|$, because a long-wavelength translation has no restoring force. An optical phonon branch has nonzero frequency at $q=0$ and represents relative motion of atoms within a primitive unit cell. Such an optical branch requires more than one independent atom per primitive cell; a truly monatomic simple-cubic model has only the acoustic branches.

The Bose occupation is $N_q=[e^{\hslash {\omega }_q/(k_{B}T)}-1{]}^{-1}$. At low temperature, appreciably populated acoustic modes haveAs temperature falls, the available phonon wavevectors concentrate near the zone centre and optical modes freeze out when $k_{B}T\ll \hslash {\omega }_{\mathrm{op}}(0)$. These are thermal-occupation restrictions, not a change in the allowed geometric Brillouin zone.

Let adjacent equal masses $m$ be separated by $a$, with alternating strong and weak spring constants $K$ and $\kappa$, $K\gg \kappa >0$. The pattern repeats after $2a$, even though the masses are identical. Thus the first Brillouin zone isChoose the strong bond inside each two-atom primitive unit cell. With displacement amplitudes $(u,v)e^{i(2naq-\omega t)}$, the normal mode equations areThe alternating-spring chain dispersion relation is consequentlyThe lower acoustic phonon branch has in-phase motion within the strongly bound pair and ${\omega }_{-}(0)=0$. Its small-$q$ slope is the speed of sound,The upper optical phonon branch has opposite motion within each pair, with ${\omega }_{+}(0)=\sqrt{2(K+\kappa )/m}$. Here “optical” names the gapped branch; equal neutral masses need not couple to light through an electric dipole.

In the strong-bond limit of the alternating-spring chain, the pairs act almost as rigid dimers in the low-frequency branch, while their internal stretch gives a high, nearly flat branch:At the zone edge the two frequencies are $\sqrt{2\kappa /m}$ and $\sqrt{2K/m}$, separated by a large forbidden frequency interval. Neighbouring cells are out of phase there. Both branches have horizontal tangents at the edge.

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.