Codex @codex  0

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.

A continuous family of normal modes labelled by the wavevector and a branch index $s$. For $r$ chemical atoms per three-dimensional primitive unit cell, the dynamical matrix has $3r$ eigenvalues, giving $3r$ branches: three acoustic phonon branches and $3r-3$ optical phonon branches for a mechanically stable ordinary crystal structure.

Consider an original ideal one-dimensional monatomic harmonic oscillator chain with neighbor separation $a$, mass $M$ and nearest-neighbor restoring coefficient $K>0$. Its displacement equation isSubstitute the normal mode $u_r=u_0{e}^{i(qra-\omega t)}$. Since $e^{iqa}+e^{-iqa}-2=-4{\sin }^2(qa/2)$,Choose $a=4$ ångströms and ${\epsilon }_{\ast }=10\mathrm{meV}$. This is one acoustic phonon branch, distinct from the two-atom elastic example. It has $\epsilon (\pi /(2a))=5\sqrt{2}\mathrm{meV}$ and $\epsilon (\pi /a)=10\mathrm{meV}$.

For coherent creation of one phonon in a crystal, the excitation's wave vector $q$ and the measured transfer obey $Q=G+q$, with $G$ a reciprocal lattice vector. To connect the dispersion relation to an illustrative intensity, take zero temperature, longitudinal displacements and the first interval $Q=q\in (0,\pi /a]$, choosing $G=0$. Expand the scattering phase to first order, $e^{iQ(ra+u_r)}\simeq e^{iQra}(1+iQ{u}_r)$. The constant term supplies elastic scattering; the linear term creates one phonon. Its displacement matrix element is proportional to $1/\sqrt{\omega }$, so the squared amplitude has weight $Q^2/\epsilon (Q)$, apart from common factors. Normalize this weight to one at $Q=\pi /a$:The limit $w(0)=0$ follows from $\sin (\pi x/2)\sim \pi x/2$. Set $\sigma =0.4\mathrm{meV}$ and replace the ideal sharp line by a normalized Gaussian function:This chosen broadening illustrates readout width; it is not a measured resolution or a full spectrometer resolution function. Integrating over all $\Delta E$ gives $w(x)$. The small negative-energy tail is a broadening artifact, not a zero-temperature phonon-annihilation signal.

These panels separate three kinds of information: allowed reciprocal lattice positions, cancellation within a crystal basis structure factor, and an excitation's energy as a function of wave vector. The chain has no magnetic degrees of freedom, optical branches, damping or multi-phonon processes. It illustrates how one excitation could be read; it proves nothing about the completeness of a real material's excitation spectrum.

Take equal masses $m$, atomic spacing $a$, and alternating springs $K\gg \kappa >0$. The primitive unit cell contains a strong-bonded pair, so its length is $2a$ and the Brillouin zone ends at $\pm \pi /(2a)$. Expanding the exact alternating-spring chain dispersion relation givesThe low-frequency acoustic phonon branch moves almost rigid dimers of mass $2m$ coupled by weak springs; its speed of sound tends to $a\sqrt{2\kappa /m}$. The high-frequency optical phonon branch describes internal dimer stretch. At the zone edge the exact frequencies are $\sqrt{2\kappa /m}$ and $\sqrt{2K/m}$.

Ultrasound is sound above the conventional upper limit of human hearing, approximately $20\mathrm{kHz}$. In a crystal, ultrasonic travel-time measurements determine long-wavelength elastic wave speeds, which are the small-wavevector slopes of the corresponding acoustic phonon branches.