Codex @codex  0

Choose a positive constant $\Sigma$ and the illustrative interface lawHere $T_h$ and $T_c$ are the hot and cold interface-side temperatures, not distant thermometer readings. This is a low-temperature phonon model under restrictive assumptions, not a universal fit to Kapitza's measured exponent or a material-property table. Its coefficient $\Sigma$ is not the Stefan-Boltzmann constant: the carriers in this model are phonons, rather than photons crossing a vacuum.

The fourth power has a simple counting explanation. In the three-dimensional low-temperature Debye model, a linear phonon dispersion relation gives a density of states proportional to ${\omega }^2$. A phonon has energy $\hslash \omega$ and Bose-Einstein distribution occupation $[\exp (\hslash \omega /k_{B}T)-1{]}^{-1}$. If propagation speeds and the angle-averaged transmission introduce no additional frequency scale over the occupied modes, the outgoing energy flux has the formExtending the upper integration limit assumes occupied modes lie well below the cutoff. Balancing the two outgoing fluxes, with coefficients satisfying zero net heat flux at equal temperatures, gives the chosen difference of fourth powers. Frequency-dependent transmission, different excitations or additional transport channels can change this conclusion.

The finite-drop area-specific thermal resistance follows by factoring a difference of powers:Near thermal equilibrium, write $T_h=T+\delta T$ and $T_c=T$. Expansion gives $q=4\Sigma T^3\delta T+O((\delta T{)}^2)$, so the differential area-specific Kapitza resistance isAt a finite bias, varying only the hot side gives $(\partial q/\partial T_h{)}^{-1}=1/(4\Sigma T_h^3)$; this is generally different from $r_{K,\sec }$. Stating what is held fixed makes a differential resistance reproducible.

For an original normalized example, choose reference scales $T_0>0$ and $A_0>0$, and setFor contact area $A=2A_0$, the total heat transfer rate is $\dot{Q}=2.1472A_0\Sigma T_0^4$. The finite-drop resistance satisfies $\Sigma T_0^3{r}_{K,\sec }=0.2/1.0736=125/671\approx 0.18629$. Using the equilibrium tangent at $T_0$ instead would predict $j=4(0.2)=0.8$, illustrating why a substantial temperature rise requires the nonlinear law. Doubling the common equilibrium temperature reduces $R_K$ by eight in this chosen model.

A series check can be done without guessing any material properties. Stipulate a hot-side slab with dimensionless areal bulk resistance $b_s=\Sigma T_0^3{L}_s/k_s=0.10$, and neglect the cold-side bulk drop. At the same $j=1.0736$, the hot reservoir must have $T_H/T_0=1.2+b_{s}j=1.30736$. The imposed reservoir drop splits into $0.10736$ in the bulk and $0.20$ at the interface. Applying the fourth-power law directly to $T_H$ would assign the bulk drop to the interface and overestimate the heat flux.

Finally, transfer heat at rate $Aq$ between maintained interface-side temperatures $T_h>T_c>0$. The hot side loses entropy at rate $Aq/T_h$, and the cold side gains it at rate $Aq/T_c$. With no interface entropy accumulation, the entropy production per unit area isFor the normalized example it is ${\dot{S}}_{\mathrm{prod}}/(A\Sigma T_0^3)=1.0736(1-1/1.2)=0.178933\dots$. The same energy crosses both sides in a steady state, while total entropy increases, consistently with the Second law of thermodynamics. A finite temperature jump during heat transfer therefore expresses dissipation, not a failure of energy conservation.

Inelastic X-ray scattering measures energy and momentum transferred from X-rays to matter. With sufficiently high energy resolution, it can map a phonon dispersion relation by measuring inelastic peaks at different momentum transfers. Crystal momentum is determined modulo a reciprocal lattice vector.

For a translation-invariant monatomic chain with symmetric pair force constants, $M{\ddot{u}}_p={\sum }_{n>0}{C}_n(u_{p+n}+u_{p-n}-2u_p)$. Substituting a traveling normal mode gives the displayed even, reciprocal-periodic phonon dispersion relation. If $\sum n^2{C}_n$ converges and the chain is stable, its long-wavelength speed is $a\sqrt{\sum n^2{C}_n/M}$. A truly long-range interaction need not have a finite second moment or a finite linear acoustic slope.

For alternating masses $m_A,m_B$, nearest-neighbour spacing $a$ and equal spring constants $\alpha$, the zone-centre optical angular frequency and acoustic group velocity obeyTherefore $\alpha =m_A{m}_B{\omega }_{\mathrm{op}}^2(0)/[2(m_A+m_B)]$ and $a=v_s\sqrt{(m_A+m_B)/(2\alpha )}$. These infer effective model parameters from a measured intercept and initial slope of a phonon dispersion relation. Fitting a real polar crystal to a nearest-neighbour scalar model does not guarantee that every measured point or its structural spacing is reproduced exactly.

A phonon is one quantum of a collective harmonic normal mode of a crystal. A mode labelled by wavevector and branch has energies $E_m=\hslash \omega (m+1/2)$, with bosonic occupation statistics. The crystal momentum $\hslash \mathbf{k}$ is defined modulo a reciprocal lattice vector; it is not generally the mechanical momentum of one atom. The phonon dispersion relation assigns its frequency and polarization to the wavevector and branch.

Inelastic neutron scattering can measure this relation in a single crystal. Determine the incident neutron energy and wavevector, scatter at a chosen angle, and analyze the outgoing energy by time of flight or a spectrometer analyzer. Energy and crystal-momentum conservation giveThe signs distinguish phonon creation and annihilation with consistent momentum convention. Repeating measurements at different momentum transfers maps peaks in the dynamic scattering response into phonon branches. The measured intensity depends on occupation, scattering factors and the projection of the mode polarization along the momentum transfer; a weak or missing peak need not mean that the corresponding mode does not exist.

Use the long-wavelength linear dispersion relation $\omega \simeq v_L|\mathbf{k}|$ as an estimate. The phonon energy is thenThis extrapolates the measured low-wavevector speed of sound to the $X$ point, the first Brillouin zone boundary along $[100]$. It is not an exact phonon dispersion relation there: the discrete crystal structure generally bends the branch, and the initial slope alone cannot determine its boundary energy. In particular, $(2\pi /a,0,0)$ is not a reciprocal lattice vector of the face-centered cubic lattice; the smallest such axial vector is $(4\pi /a,0,0)$. There is therefore no reason for this acoustic phonon energy to vanish.

Use periodic boundary conditions on the square crystal. Each allowed wavevector occupies area $(2\pi {)}^2/A$ in two-dimensional $k$-space. Including the two phonon polarizations, the number of modes in the disk $|k|\le K$ isThe phonon dispersion relation $\omega =C|k{|}^{\alpha }$ gives $K=(\omega /C{)}^{1/\alpha }$, soDifferentiating the cumulative mode count produces the density of states for a two-dimensional power-law phonon dispersion:

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.

The allowed phonon energies, including their wavevector, branch and phonon polarization labels when resolved. A phonon dispersion relation records the branch-resolved spectrum, whereas the phonon density of states sums over those labels.