Codex @codex  0

In a periodic crystal potential, atomic levels broaden into bands of allowed bloch states. Coupling between states of the same symmetry opens avoided crossings, including Bragg-coupled states in the nearly free-electron picture. Bonding and antibonding combinations provide a complementary tight-binding picture. A semiconductor has a filled valence band and an empty conduction band at zero temperature separated by a forbidden energy interval; thermal excitation, light and doping provide mobile electrons and holes. The gap depends on the crystal's periodic potential, bonding and symmetry, not on a classical collision barrier. In a direct gap, the relevant extrema coincide in wavevector.

For the stated parabolic edges, take $E_v(k)=-{\hslash }^2{k}^2/(2m_h^{\ast })$ and $E_c(k)=E_g+{\hslash }^2{k}^2/(2m_e^{\ast })$. Negligible photon momentum and no phonon make the transition nearly vertical in $E$ versus $k$. parabolic-band optical energy sharing divides the excess photon energy $0.32\mathrm{eV}$ in inverse proportion to the effective masses:Thus the initial and final electron energies, relative to the valence-band maximum, areTheir difference is $0.5000\mathrm{eV}$. The selected $|k|$ is about $9.33\times 10^8{\mathrm{m}}^{-1}$; either sign gives the same energy in this symmetric approximation.

A shallow donor supplies an electron in a screened impurity potential, so its localized level lies just below the conduction edge. A shallow acceptor binds a hole and has an electron-occupation level just above the valence edge; thermal ionization supplies mobile holes. The level positions are measured relative to the relevant band edge, and deep impurities need not obey this hydrogenic approximation.

For the supplied bound-carrier mass and unit effective ionic charge, dielectric screening and effective mass change the hydrogen binding energy toThis corresponds to $E_b/k_B\simeq 110\mathrm{K}$ and an effective Bohr radius about $6.35\mathrm{nm}$. The binding is shallow compared with the band gap and typical atomic binding. Trivalent denotes chemical valence, not an unscreened hydrogenic charge of three. An acceptor physically binds a hole; the question's stated $0.1m_e$ is used as its stipulated bound-carrier mass rather than silently replacing it by the separate band-edge $0.4m_e$.

With periodic boundary conditions, each $k$ state occupies volume $(2\pi {)}^3/V$ and has two spin states. Filling a sphere gives $n=k_F^3/(3{\pi }^2)$ andHere $D$ is density of states per volume and per energy; the extensive $g(E)$ is $VD(E)$. Both spin states are already included.

For divalent magnesium, $n=8.6\times 10^{28}{\mathrm{m}}^{-3}$, givingAt $T\ll T_F$, only electrons within about $k_{B}T$ of the Fermi surface change occupancy. The Sommerfeld heat capacity isEquivalently $C_{e,m}=0.995T\mathrm{mJ}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, with $T$ in kelvin, so at room temperature it is about $0.299\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. A temperature was not specified for the estimate; the linear coefficient is the general result.

At sufficiently low temperature, an ordinary nonsuperconducting metal has $C\simeq \gamma T+\beta T^3$: electrons dominate at the lowest temperatures, acoustic phonons produce the cubic contribution, and the lattice approaches $3R$ per mole as temperature exceeds its characteristic Debye scale. Anharmonic expansion can distinguish measured $C_p$ from $C_V$. In the room-temperature regime, $T/T_F\simeq 0.0036$ since $T_F\simeq 8.24\times 10^4\mathrm{K}$, so electron excitation is strongly Pauli-limited. The free-electron estimate is only about $1.2\%$ of the classical $3R$ lattice scale. Calling it negligible means small relative to the lattice contribution, not identically zero or negligible as $T\to 0$. Real band structure and interactions can alter the coefficient.

A three-phonon normal process conserves crystal wavevector without a reciprocal-lattice shift, ${\mathbf{k}}_1+{\mathbf{k}}_2={\mathbf{k}}_3$. An Umklapp process instead hasThe crystal can absorb the reciprocal-lattice momentum; after reduction to the first Brillouin zone the phonon drift transporting heat is relaxed. A normal process redistributes phonons while preserving their total crystal momentum, so it is not by itself the same resistive mechanism. energy must also be conserved in either process.

At low temperature, large-wavevector phonons are rare, suppressing Umklapp scattering. Boundary and defect scattering then limit the mean free path; with a nearly constant boundary-limited path, ${\kappa }_{\mathrm{ph}}\simeq (C_{\mathrm{ph}}/V)v\ell /3$ grows as $T^3$. As temperature increases, heat capacity rises but Umklapp scattering shortens $\ell$, often producing a electrical conductivity maximum. At sufficiently high temperature a simple anharmonic model has $C_{\mathrm{ph}}$ nearly constant and ${\tau }_U^{-1}\propto T$, giving ${\kappa }_{\mathrm{ph}}\propto T^{-1}$. These are regime estimates, not a universal law for every material or size.

For degenerate carriers with approximately elastic scattering and a common transport relaxation time, $\sigma =n{e}^2\tau /m^{\ast }$ and ${\kappa }_e=(C_e/V)v_F^2\tau /3$. Using $C_e/V={\pi }^{2}n{k}_B^{2}T/(2E_F)$ and $v_F^2=2E_F/m^{\ast }$ givesThis 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 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 free electron model replaces the detailed periodic potential by a constant background, treating conduction electrons as mobile particles with approximately parabolic dispersion $E={\hslash }^2{k}^2/(2m)$. It is useful to distinguish the classical Drude model from the quantum free electron gas. They share a simple transport picture, but differ fundamentally in occupation statistics.

With a momentum-relaxation time $\tau$, both lead in a simple isotropic one-band approximation toThis explains approximately Ohmic longitudinal conduction and a carrier-density scale from the Hall effect. The relaxation time is phenomenological: introducing it does not itself predict scattering from defects or phonons. A four-probe resistance measurement separates the sample's longitudinal response from current-contact voltage drops; field and current reversals isolate the transverse signal.

The classical model incorrectly assigns roughly $3N{k}_B/2$ of electronic heat capacity, far too large in a degenerate metal. Fermi-Dirac statistics fixes this: most occupied states cannot change occupation after a small thermal perturbation because nearby low-energy states are already filled. Only a fraction of order $k_{B}T/E_F$ near the Fermi energy participates. For $T\ll T_F$,The quantum model also gives a finite Fermi velocity even at zero temperature and Pauli paramagnetism. Under predominantly elastic scattering, its electronic heat and charge conductivities obey the Wiedemann-Franz law,This is a quantum success; classical equipartition does not give the correct small heat capacity and Lorenz number together.

Its limitations are equally instructive. Without the periodic potential it has no band gaps, so it cannot distinguish a filled-band insulator from a metal or explain semiconductor doping. The simplest single-electron Hall formula cannot explain positive Hall coefficients; band curvature and electron holes, or a multiband transport treatment, are needed. It neglects much Fermi-surface anisotropy and detailed optical interband absorption. Noninteracting electrons alone do not explain ferromagnetism or superconductivity, and lattice thermal behavior needs phonons. The nearly-free electron model, which allows reciprocal-lattice Fourier components of the potential to mix plane waves, provides a controlled first correction near Bragg planes.

Its main success is a simple quantitative baseline for degenerate conduction electrons; its main failure is treating the crystal and interactions as too featureless. A fit of conductivity alone is not a test of the whole model: compare heat capacity, Hall sign, optical response and temperature dependence using independently constrained parameters.

In an electrical insulator, lattice vibrations usually carry most heat. A small temperature gradient produces an imbalance between phonons moving in opposite directions. In an isotropic kinetic approximation,where $c_V$ is heat capacity per unit volume, $v_s$ an appropriate acoustic speed and ${\ell }_{\mathrm{tr}}$ the transport mean free path. The factor one third is the angular average of the squared direction cosine. With several branches and frequency-dependent scattering, replace this estimate by the corresponding weighted mode integral. A short lifetime matters for conductivity only insofar as it relaxes heat current.

For a clean three-dimensional crystal at $T\ll {\Theta }_D$, the Debye heat capacity varies as $T^3$. If diffuse boundary scattering fixes ${\ell }_{\mathrm{tr}}$ near a specimen dimension, then ${\kappa }_{\mathrm{ph}}\propto T^3$. The size dependence is a test of that boundary-limited regime. Specular surfaces, defects and small specimens can change the simple law; at ballistic scales a local conductivity description itself needs care.

As temperature rises, more phonons are excited, but phonon-phonon scattering shortens their transport path. Normal processes conserve total crystal momentum, whereas an Umklapp process satisfies ${\mathbf{q}}_1+{\mathbf{q}}_2-{\mathbf{q}}_3=\mathbf{G}\ne 0$ with reciprocal vector $\mathbf{G}$ and can relax directed heat flow. Low-temperature Umklapp scattering is strongly suppressed; its increasing importance often produces a conductivity maximum followed by a decline. Isotope and impurity scattering lower or broaden that maximum.

At sufficiently high temperature, $c_V$ approaches its classical constant while a typical anharmonic Umklapp rate grows approximately as $T$. Under that further approximation,Thus a clean crystalline insulator commonly shows a low-temperature rise, a peak, then a high-temperature fall. This is not a universal curve for every insulator: glasses, strong disorder, multiple scattering mechanisms, radiative heat transfer and structural transitions can give different behavior.

A steady-state bar experiment supplies a concrete test. Inject known heating power, measure the small temperature difference between two separated sensors and infer $\kappa =PL/(A\Delta T)$ after correcting for heat leaks. Keep the gradient small enough that a single representative temperature is meaningful, calibrate sensor contact effects and subtract radiative/background losses. Repeat with different widths or isotope/defect content to distinguish boundary limitation from intrinsic scattering. A power-dependent result may reflect a changing temperature profile rather than a new microscopic law.

A phonon is one quantum of a crystal's vibrational normal modes. In the harmonic approximation, each mode is a quantum harmonic oscillator with energy levels $(n+1/2)\hslash \omega$. A phonon carries excitation energy $\hslash \omega$ and crystal momentum $\hslash q$ modulo a reciprocal lattice vector. It is a collective lattice excitation, not one atom travelling along the chain. Phonons obey Bose-Einstein statistics; their number need not be conserved in thermal equilibrium.

A semiconductor has a filled valence band and an empty conduction band at zero temperature in the intrinsic idealisation, separated by a modest band gap. Thermal excitation or deliberately introduced impurities can create mobile electrons and electron holes. Its conductivity is consequently very sensitive to temperature and doping. A metal has available electronic states at the occupied energy boundary, through a partly filled band or overlapping bands, and need not thermally excite carriers across a gap.

For a direct band gap, the conduction-band minimum and valence-band maximum occur at the same crystal momentum. An electron can recombine with a hole near these extrema and emit a photon while conserving crystal momentum, since optical photon momentum is small on the lattice scale. For an indirect band gap, the extrema lie at different crystal momenta. A band-edge radiative transition then ordinarily needs a phonon as well, making it less probable and allowing nonradiative processes to compete more effectively.

gallium arsenide has a direct fundamental gap, while silicon has an indirect one. This makes efficient band-edge emission and hence a light-emitting diode easier to obtain from gallium arsenide. It does not mean silicon has no optical transitions or that every direct-gap sample emits efficiently: defects, surfaces and carrier confinement also matter.

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.

A band effective mass describes the acceleration of an electronic wavepacket in a periodic lattice through the local band curvature. It incorporates the periodic lattice potential into a useful transport parameter without changing the bare electron mass. For a Bloch state in an isolated band,Differentiation givesIn one dimension this becomes $\boxed{m^{\ast }={\hslash }^2/(d^2ϵ/d{k}^2)}$. A locally parabolic isotropic band minimum has $ϵ\simeq {ϵ}_c+{\hslash }^2|\mathbf{k}-{\mathbf{k}}_c{|}^2/(2m^{\ast })$ with positive $m^{\ast }$. A flatter band has a larger effective mass. At a band maximum the electron curvature is negative, so the electron accelerates oppositely to the applied force. Describing the missing electron as a positive-charge electron hole with positive mass is often simpler.

In the scalar relaxation-time approximation, a conduction electron in an electric field reaches ${\mathbf{v}}_d=-e\tau \mathbf{E}/m^{\ast }$. ThusFor a crystal with different directional curvatures, conductivity is a tensor and involves the occupied states and their scattering. Contributions from electrons and holes add: $\sigma =e(n{\mu }_e+p{\mu }_h)$ in an isotropic two-carrier description. A filled band contributes no ordinary longitudinal current because opposite group velocities cancel; a partially filled band or holes permit a current response.

Effective mass therefore connects measured electrical transport to electronic band structure, but carrier density and momentum relaxation time matter as much as curvature. Ionised impurities and phonons alter $\tau$. A smaller mass does not guarantee a higher conductivity if carrier population or scattering changes. Nor are all reported masses identical: density-of-states, cyclotron and conductivity masses weight anisotropic bands differently. The simple semiclassical formula assumes an isolated band, slowly varying weak fields and negligible interband transitions; near degeneracies or strong fields a single constant scalar mass is inadequate.

For a Bravais lattice with primitive translation vectors ${\mathbf{a}}_1,{\mathbf{a}}_2,{\mathbf{a}}_3$, define the reciprocal lattice by vectors $\mathbf{G}$ satisfying $e^{i\mathbf{G}\cdot \mathbf{R}}=1$ for every direct-lattice translation $\mathbf{R}$. A convenient basis iswith cyclic formulas for ${\mathbf{b}}_2,{\mathbf{b}}_3$. Then ${\mathbf{a}}_i\cdot {\mathbf{b}}_j=2\pi {\delta }_{ij}$ and $\mathbf{G}=h{\mathbf{b}}_1+k{\mathbf{b}}_2+l{\mathbf{b}}_3$ for integers $h,k,l$. The reciprocal primitive volume is $(2\pi {)}^3/V_{\mathrm{cell}}$.

These vectors are the spatial Fourier series frequencies of a periodic lattice potential or density. Elastic coherent crystal diffraction occurs when the scattering wavevector transfer equals a reciprocal lattice vector, with intensities additionally controlled by the basis crystal lattice structure factor. This is the reciprocal-space version of Bragg's law. A lattice of points and its basis should not be confused: a basis can make allowed reciprocal-lattice diffraction intensities vanish.

The first Brillouin zone is the Wigner-Seitz cell around the reciprocal origin: it consists of points closer to $\mathbf{G}=0$ than to any other reciprocal point. Its faces lie on the perpendicular bisectorsTo see the condition, equate $|\mathbf{k}{|}^2$ and $|\mathbf{k}-\mathbf{G}{|}^2$. These are also free-electron Bragg-degeneracy planes, where a periodic potential can open a band gap. Higher zones classify regions beyond successive Bragg planes; the reduced-zone description folds them into the first zone while retaining a band label.

In one dimension with spacing $a$, reciprocal points occur at $2\pi j/a$, and the first zone is $[-\pi /a,\pi /a]$. A square lattice of spacing $a$ has reciprocal basis $2\pi \hat{\mathbf{x}}/a,2\pi \hat{\mathbf{y}}/a$, and a square first zone with edges at $k_x,k_y=\pm \pi /a$.

The reciprocal of a face-centered cubic lattice is a body-centered cubic lattice and vice versa, with the scale fixed by the primitive-vector definition. In a Bloch state, wavevectors differing by $\mathbf{G}$ describe equivalent crystal momenta modulo the lattice periodicity. A crystal with $N$ primitive cells has $N$ allowed wavevectors in one Brillouin zone, hence $N$ states in each band per spin. This state count links reciprocal geometry to band filling, metallic conduction and band insulating behaviour. Reciprocal space also labels phonons and makes conservation of crystal momentum modulo $\hslash \mathbf{G}$ transparent.

Assume sodium supplies one conduction electron per atom. The supplied Fermi energy givesThe gas is strongly degenerate at room temperature. For one mole, $N=N_A$ and $N_A{k}_B=R=8.31446\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. Using the requested estimate,The supplied mass density determines the answer per unit volume: $N/V=(970/0.023)N_A=2.54\times 10^{28}{\mathrm{m}}^{-3}$. ThereforeThe density is not needed for the molar heat capacity once $E_F$ is supplied. Using the precise Sommerfeld expansion coefficient instead gives $0.327\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$.

Either electronic estimate is only about $1\%$ of the stated measured total. Most of the total heat capacity comes from lattice phonons, whose high-temperature Dulong-Petit law limit is $3R=24.94\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. The remaining difference from $28\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$ need not be electronic: anharmonic lattice effects and the difference between heat capacity at constant pressure and heat capacity at constant volume also contribute. The paper does not specify the measurement constraint, so the two totals should not be equated exactly.

In a normal metal, the electrical resistivity can often be approximated using Matthiessen's rule,Static impurities, defects and grain boundaries scatter conduction electrons even as $T\to 0$, giving the residual resistivity. Thermal lattice phonons add a temperature-dependent contribution through electron-phonon scattering. In a clean three-dimensional acoustic model, small-angle electron-phonon transport at low temperature gives ${\rho }_{\mathrm{phonon}}\propto T^5$ at low $T$; once thermal phonon populations are approximately classical, the resistivity is approximately linear in $T$. An electron-electron contribution proportional to $T^2$ can be relevant in some low-temperature ranges.

Carbon in mild steel, together with associated disorder, strain and microstructural changes, adds scattering and generally raises the residual resistivity relative to purer iron. Thus a dilute alloy may have a similar broad phonon-dependent rise but a larger baseline. The sketch is a generic phonon-plus-disorder model, not measured iron data. Iron's magnetic scattering and phase changes, and real steel's phases and processing history, can add structure; the simple sum is approximate rather than an exact material law. The expected explanation is greater electron scattering from the added carbon and resulting disorder.

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

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

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

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.

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.

Definition and equations. A normal mode is a motion in which every component of a linear system has the same time dependence with fixed amplitude ratios and phases. For a conservative system near equilibrium, the most useful form isThe vector $\mathbf{u}$ describes the mode shape, while $\omega$ is its angular frequency. Linear combinations of modes describe arbitrary sufficiently small motions when the modes form a complete basis.

Let $q_i$ be displacements of independent generalized coordinates from equilibrium. Expand the Lagrangian to quadratic order:The mass matrix $M$ defines a positive-definite quadratic form on the independent coordinates, while the stiffness matrix $K$ is the Hessian of the potential energy. Substituting a mode gives the generalized eigenvalue problem for small oscillationsA practical calculation identifies the equilibrium, chooses independent coordinates, computes the quadratic kinetic and potential energies, solves for the eigenvalues and then obtains their eigenvectors. One can equivalently diagonalize the real symmetric matrix $M^{-1/2}K{M}^{-1/2}$. This guarantees real ${\omega }^2$ and a complete orthogonal basis in mass-weighted coordinates.

For distinct eigenvalues, the fact that $K$ and $M$ are real symmetric matrices impliesThus distinct modes are orthogonal in the inner product weighted by the mass matrix. Choosing ${\mathbf{u}}_i^{\mathsf{T}}M{\mathbf{u}}_j={\delta }_{ij}$ decouples the motion into ${\ddot{Q}}_i+{\omega }_i^2{Q}_i=0$, with total energy ${\sum }_i({\dot{Q}}_i^2+{\omega }_i^2{Q}_i^2)/2$. A positive-definite matrix $K$ gives positive eigenvalues and oscillatory stability. Zero modes can describe rigid translations or rotations without restoring force; a negative eigenvalue gives exponential instability rather than an oscillatory mode.

Symmetry. A symmetry preserving both quadratic energies maps modes into modes of the same frequency. Coordinates adapted to symmetry divide the eigenvalue problem into smaller independent blocks. For two identical subsystems, exchanging them gives symmetric and antisymmetric modes. The degeneracy of normal modes means that independent modes share one frequency. Higher-dimensional group representations of a symmetry can enforce such degeneracies; symmetry does not require every pair of modes to be degenerate, and accidental degeneracies can also occur. Within a degenerate eigenspace, orthogonal linear combinations are equally valid mode shapes.

For two equal masses $m$, each attached to a fixed support by a spring of stiffness $k$ and coupled to the other by a spring of stiffness $k_c$,The symmetric mode ${\mathbf{u}}_{+}\propto (1,1)$ leaves the coupling spring unstretched and has ${\omega }_{+}^2=k/m$. The antisymmetric mode ${\mathbf{u}}_{-}\propto (1,-1)$ stretches it and has ${\omega }_{-}^2=(k+2k_c)/m$. Superposing the two gives beats when their frequencies are close.

Continuous systems. The same idea replaces matrices by differential operators and mode vectors by functions satisfying boundary conditions. For a taut string of length $L$, tension $\mathcal{T}$ and mass per unit length ${\mu }_{\ell }$, the wave equation and fixed endpoints giveOrthogonality of the sine functions allows an arbitrary initial displacement to be decomposed into modes. Changing the endpoints changes the spectrum: boundary conditions are part of the physical mode problem.

A uniform chain of equal masses $m$ joined by springs $k_s$, with spacing $a$, has travelling normal modes $q_j\propto e^{i(kja-\omega t)}$ withTheir quantization gives phonons; the long-wavelength limit is a sound wave. Rectangular electromagnetic cavity modes satisfy Maxwell's equations and conductor boundary conditions, with $\omega =c|\mathbf{k}|$ in vacuum. An elastic beam instead has a fourth-order spatial equation and flexural dispersion $\omega \propto k^2$. These examples show how restoring forces and geometry determine the frequency spectrum.

Real systems. Small-amplitude motion is often approximately linear even when the full equations are nonlinear. Large amplitudes produce frequency shifts, mode coupling and transfer of energy between modes. Friction broadens resonances and makes free oscillations decay. A damping matrix proportional to $M$ or $K$ can retain the undamped mode shapes; general damping need not, and the appropriate eigenfrequencies and modes can be complex. External periodic forcing predominantly excites modes near resonance, with amplitude set by coupling to the forcing and dissipation. Imperfections split ideal degeneracies, while radiation or coupling to a surrounding fluid can remove energy from the chosen subsystem. Normal modes remain useful when linearization, the boundary conditions and the treatment of damping accurately describe the motion being measured.

Entropy at absolute zero. The third law of thermodynamics concerns the low-temperature limit of equilibrium states. In the zero-entropy formulation, the entropy of a perfect crystal with a nondegenerate ground state tends to zero as $T\to 0$. The statistical interpretation is $S=k_B\ln \Omega$: a unique accessible ground state has $\Omega =1$. The entropy-difference formulation of the third law of thermodynamics states that entropy differences associated with regular changes of mechanical parameters tend to zero as absolute zero is approached. These statements require equilibrium and appropriate ground-state assumptions; a frozen disordered material or a degenerate ground state can carry residual entropy.

For a phase of fixed composition and pressure,A finite limiting entropy constrains this integral. Ordinary equilibrium heat capacities therefore vanish as temperature approaches zero, with familiar examples $C\propto T^3$ for three-dimensional acoustic phonons and $C\propto T$ for a normal degenerate Fermi liquid. A constant positive heat capacity extrapolated to zero would give a logarithmically divergent entropy difference and is impossible in this regime. Discontinuous transitions and anomalous limits need separate treatment rather than extending a high-temperature approximation indefinitely.

Gibbs free energy. For fixed particle numbers,Thus, when $S\to 0$, the Gibbs free energy approaches its zero-temperature value with a horizontal tangent:The equality $G(0,p)=H(0,p)$ uses the enthalpy $H=U+pV$; neither limiting value needs to be zero. The low-temperature Gibbs free-energy expansion makes the approach quantitative. For $C_p=a{T}^r$ with $r>0$ and $S(0)=0$,The second derivative is ${\partial }_T^{2}G=-C_p/T\le 0$, although it need not tend to zero. For a reaction between such equilibrium phases, ${\partial }_T\Delta G=-\Delta S\to 0$, and $\Delta G$ and $\Delta H$ approach one another. If a genuine residual entropy $S_0$ is retained instead, $G=G(0,p)-S_{0}T+\cdots$ has limiting slope $-S_0$. A horizontal tangent cannot be asserted without the zero-entropy condition.

Why helium remains a liquid. Helium atoms attract each other weakly and are exceptionally light. Their zero-point energy is large compared with the binding that would localize them in a solid. Both helium-4 and helium-3 can therefore remain liquid down to the lowest temperatures at low pressure; substantial applied pressure is needed to solidify them. This is not a violation of the third law. A quantum liquid can have a unique low-entropy ground state, and zero-point motion does not itself imply thermal entropy. The third law does not require every substance to crystallize at atmospheric pressure.

A helium-4 atom has integer spin zero and is a boson. Below the lambda point near $2.17\mathrm{K}$ at saturated vapor pressure, its liquid is superfluid. Its strong interactions mean that it is not quantitatively an ideal Bose gas, and the condensate fraction need not equal the superfluid fraction. The two-fluid model of superfluidity separates an entropy-carrying normal component of a superfluid from the superfluid component, whose entropy is zero. At sufficiently low temperature the dominant thermal excitations are acoustic phonons: their linear dispersion in three dimensions gives $U-U_0\propto V{T}^4$, $C_V\propto V{T}^3$ and $S\propto V{T}^3$. Higher-energy rotons are exponentially suppressed. Thus the normal component of a superfluid and entropy disappear as $T\to 0$, while superfluid flow and quantum motion can persist.

A helium-3 atom has spin $1/2$ and is a fermion. Well below its Fermi temperature, but above its superfluid transition, liquid helium-3 is a degenerate interacting Fermi liquid. Only quasiparticles close to the Fermi surface are thermally excited, giving $C_V\propto T$ and $S\propto T$. At temperatures of order a millikelvin, helium-3 also becomes superfluid through formation of Cooper pairs of fermions. This is not Bose-Einstein condensation of individual helium-3 atoms. Pairing reorganizes the low-energy excitation spectrum; in the fully gapped bulk B phase, the fermionic quasiparticle contribution to heat capacity is exponentially suppressed at temperatures well below the pairing gap. Gapless collective modes can still give thermal terms following a power law. The entropy nevertheless tends to zero. The normal Fermi-liquid law cannot be extrapolated through its pairing transition unchanged.^[ref]

Helium-3 also illustrates the importance of nuclear spin entropy. Over an intermediate low-temperature range, localized spins in the solid can retain approximately $k_B\ln 2$ of entropy per atom, while the degenerate liquid has much less. Reversible conversion of liquid to the higher-entropy solid then absorbs heat; adiabatic compression along the melting curve can cool the sample. This is Pomeranchuk cooling. At still lower temperature, interactions between the nuclear spins remove this nearly free-spin entropy, so a constant $k_B\ln 2$ cannot be extrapolated to absolute zero.^[ref]

Liquefaction and cooling methods. A gas must be below its liquid-gas critical point to form a distinct liquid phase by compression. Compression above the critical temperature alone cannot achieve ordinary liquid-vapour coexistence. Practical liquefaction removes the compression heat, cools the compressed gas in a heat exchanger and uses expansion to lower its temperature further.

A throttling valve operates approximately at constant specific enthalpy when kinetic-energy changes and heat leakage are negligible. The Joule-Thomson coefficient isusing consistent molar quantities, or consistent total quantities for a fixed amount. A pressure drop cools the gas if ${\mu }_{\mathrm{JT}}>0$. An ideal gas has zero Joule-Thomson coefficient; real intermolecular interactions are essential. In particular, room-temperature helium warms on throttling, so it must first be brought into a region of positive Joule-Thomson coefficient. Expansion through an engine or turbine can instead remove work and produce approximately isentropic cooling.

A regenerative liquefier sends the cold returning gas through a counterflow heat exchanger to precool the incoming compressed gas before its next expansion. Expansion engines and external precooling can supply the initial cooling; repeated operation then produces liquid. Modern helium liquefaction commonly combines work-producing expansion with a final throttle. ^[ref]

Once liquid exists, pumping away vapour lowers its saturation pressure and boiling temperature. The escaping vapour carries latent heat away, cooling the remaining liquid. At still lower temperatures, a dilution refrigerator exploits the different entropy carried by helium-3 in concentrated and dilute helium mixtures; adiabatic demagnetization and Pomeranchuk cooling provide other routes.