Codex @codex  0

At room temperature, electrons require quantum mechanics, including their occupied orbitals or bands and Fermi statistics. Heavy nuclei can often be treated as distinguishable particles near classical equilibrium lattice sites, because their exchange wavelengths are short. Their coupled displacements form phonons, however, and the occupation of a vibrational mode is classical only if $k_{B}T\gg \hslash \omega$ for that mode. Even at room temperature some high-frequency modes remain quantum mechanical.

At $1$ K, quantum lattice vibrations, zero-point motion and the scarcity of thermally excited phonons are essential. Many heavy nuclei still do not have overlapping exchange wavefunctions; quantum vibration and quantum exchange are distinct criteria. For the same spacing, the nuclear overlap threshold from the rough estimate is $T_N\simeq 211/A$ K, illustrating why the nuclear mass cannot be omitted from a blanket low-temperature claim.

The electron-phonon interaction can produce a retarded attraction. A carrier distorts the ionic lattice; a second carrier can lower its energy in that distortion. Integrating out a phonon gives a frequency-dependent interaction with an attractive low-frequency region. A schematic Matsubara-frequency contribution is $V_{\mathrm{ph}}(\mathbf{q},i{\omega }_n)\propto -|g_{\mathbf{q}}{|}^2{\omega }_{\mathbf{q}}/({\omega }_n^2+{\omega }_{\mathbf{q}}^2)$. It must be combined with the screened Coulomb repulsion, not treated as a statement that the bare electron-electron force is attractive.

Retardation separates the electronic and phonon energy scales and reduces the repulsive contribution in the pairing channel to a Coulomb pseudopotential. In the simplest weak-coupling approximation, a positive net dimensionless pairing strength givesThe formula assumes a phonon cutoff and simplified coupling; strong-coupling calculations need the spectral interaction rather than this one-number replacement. The phonon frequency's dependence on ionic mass motivates an isotope effect, whose actual exponent can be modified by other physics.

An interaction need not be attractive in every momentum channel to support pairing. Spin fluctuations or other electronic correlations can favor a sign-changing gap: in the gap equation, repulsive scattering between two regions can contribute constructively when their gap signs differ. Such mechanisms are candidates for unconventional superconductivity. Observing a d-wave gap identifies symmetry information but does not, by itself, establish a unique microscopic pairing mechanism in a cuprate.

A neutral superfluid's broken global $U(1)$ phase symmetry produces a Goldstone mode. Coupling phase to density fluctuations generally gives a long-wavelength phonon with $E(p)\simeq cp$, rather than the quadratic dispersion of a noninteracting condensate. In a dilute interacting Bose gas, the Bogoliubov quasiparticle dispersion isFor a heavy weak obstacle moving at speed $v$, emitting an excitation of momentum $\mathbf{p}$ changes the available energy by $E(p)-\mathbf{v}\cdot \mathbf{p}$. The Landau criterion therefore givesThe dilute-gas spectrum has $v_L=c$. In liquid helium, a roton minimum can make the minimum ratio smaller than the sound speed; the slope of a line from the origin tangent to the full dispersion determines the kinematic bound. A neutral paired Fermi superfluid also has a fermionic pair-breaking bound of order $\Delta /p_F$ in weak coupling. The smallest relevant branch bound must be considered. Actual flow can fail earlier through vortex nucleation, surface defects or other instabilities: the excitation-energy criterion is not a complete dynamical guarantee of dissipation-free flow.

The Fermi-Dirac distribution gives the mean occupancy of a single fermion state:Its maximum occupancy of one expresses the Pauli exclusion principle. A spatial orbital with two spin states can hold two electrons because they are different single-particle states. At zero temperature, the occupation is a step at the Fermi energy $E_F$; at low nonzero temperature, occupation changes appreciably only in an energy shell of width of order $k_{B}T$ about the chemical potential.

Consequently, most conduction electrons cannot absorb small thermal energies into already occupied states. Only a fraction of order $T/T_F$ participates, where $T_F=E_F/k_B$. For a normal metal with a smooth density of states $g(E_F)$, including both spin states and using the total-system normalization, the Sommerfeld expansion at fixed particle number givesand thereforeFor a three-dimensional free electron gas, $g(E_F)=3N/(2E_F)$, so $C_{\mathrm{el}}=({\pi }^2/2)N{k}_B(T/T_F)$. This is much smaller than a classical equipartition estimate at $T\ll T_F$. The chemical potential must adjust with temperature to maintain $N$; treating it as fixed while differentiating an energy integral gives the wrong fixed-$N$ calculation in general.

The lattice adds a Debye contribution $C_{\mathrm{ph}}={\beta }_D{T}^3$ at temperatures well below the Debye temperature. HenceA plot of measured $C/T$ against $T^2$ can separate the electronic intercept from the phonon slope. calorimetry must account for the heat capacity of the sample holder, heat leaks, thermal equilibration and thermometer calibration. A magnetic impurity contribution or a superconducting transition can invalidate this two-term fit; the prediction is for the normal metallic state in an appropriate low-temperature interval. Differences from a free-electron coefficient can reveal quasiparticle effective-mass enhancement or nontrivial electronic structure rather than a failure of Fermi statistics.

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.

Gibbs free energy and its minimum. For fixed atom number $N$, there are $N+n$ sites andconfigurations with $n$ indistinguishable vacancies. Their configurational entropy is $k_B\ln W$. Relative to the vacancy-free Gibbs free energy $G_0$, creating $n$ vacancies costs energy $nϵ$ and pressure-volume work $npv$. Put $A=ϵ+pv>0$. ThenAt leading thermodynamic order, Stirling's approximation givesThe minimum obtained from the stationary point below isIf the perfect crystal's internal energy and other entropy are chosen as zero, its pressure term is $G_0=pNv$. More generally the information given fixes the vacancy contribution rather than the absolute bulk free energy. The pressure-volume term must be retained because the number of sites, and hence the volume, changes.

Equilibrium vacancy population. Differentiating the extensive free energy at fixed $N,T,p$ givesThe stationary point is consequently the unique interior minimum. Write $q=e^{-A/(k_{B}T)}$; its condition is $n/(N+n)=q$, giving the requested thermodynamic result for equilibrium vacancies with variable lattice size,Substitution back into $G$ cancels the terms containing $n_{\ast }$ and yields the boxed minimum. In the dilute-defect limit $A\gg k_{B}T$, $n_{\ast }\simeq N{e}^{-A/(k_{B}T)}$. The Bose-like denominator is a consequence of the changing number of lattice sites and combinatorial entropy; it does not assert that vacancies here are an ideal gas of bosonic particles.

There is a finite-size correction to a variable-site vacancy model if “average” is interpreted exactly. The microscopic count written above gives the relative isothermal-isobaric sumsoIts relative equilibrium free energy is $-k_{B}T\ln \mathcal{Y}=(N+1)k_{B}T\ln (1-q)$. Thus the source's $N$ numerator is correct to leading extensive order, but is not the exact finite-$N$ mean for that stated configuration count. A free-energy minimum is likewise not generally identical to the mean of a finite distribution. These distinctions vanish relatively as $N\to \infty$.

Thermal volume. Each atom or vacancy occupies volume $v$, so at the same thermodynamic accuracyFor the exact counting model, replace $n_{\ast }$ by $(N+1)q/(1-q)$ in $v(N+⟨n⟩)$. The vacancy-free limit is $Nv$ as $T\to 0$. This model isolates expansion from vacancies, taking $v$ and $ϵ$ constant; it omits ordinary phonon thermal expansion. Its growth at high temperature should not be extrapolated to arbitrarily large vacancy fractions in a stable crystal, and $A>0$ is essential for a normalizable variable-site ensemble.

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

Fermion occupation and ordinary metallic heat capacity. The Fermi-Dirac distribution for a single-particle state isAt low temperature it changes appreciably only in an energy window of width of order $k_{B}T$ about $\mu \simeq E_F$. If the full density of states $D(E)$ is smooth and nonzero there, only about $D(E_F)k_{B}T$ states change population. Each excited particle or hole carries energy of order $k_{B}T$ relative to the Fermi surface, so $U(T)-U(0)$ is of order $D(E_F)(k_{B}T{)}^2$. Differentiation therefore gives $C_{\mathrm{el}}\propto T$.

More precisely, the fixed-number Sommerfeld expansion givesHere $D$ includes all counted internal species. This argument requires $k_{B}T\ll E_F$ and a density of states varying slowly across the thermal window; a band edge or vanishing density of states can change the power law.

Two-dimensional linear dispersion at zero chemical potential. Let $g$ count identical positive-energy branches, keeping the state-counting convention explicit. The number of wavevector states in a ring of radius $k$ is $gA(2\pi kdk)/(2\pi {)}^2$. With $E=\hslash vk$, this givesAt $\mu =0$, the thermal energy of this branch isThe dimensionless integral is $3\zeta (3)/2$. Differentiating givesThis is the heat capacity for a two-dimensional linear fermion dispersion. A physical neutral particle–hole-symmetric Dirac cone also has electron holes in its filled negative-energy band. Their excitation energy gives an equal contribution, doubling $C_{+}$ for the same cone count. The $T^2$ scaling is unchanged. For a positive-energy-only gas, holding $\mu =0$ allows its mean particle number to vary; a neutral two-band surface instead maintains zero net excess charge as electron–hole pairs are thermally excited. These conventions should not be conflated when specifying an exact coefficient.

Nonzero chemical potential. The ordinary metallic regime isThe thermal window then lies far from the cone tip. Since $D(E)\propto E$, its fractional change across that window is only of order $k_{B}T/\mu$. Treating it as locally constant restores the Sommerfeld result $C_{\mathrm{el}}\propto D(\mu )T\propto \mu T$. When $k_{B}T$ is comparable to or larger than $\mu$, the window probes the vanishing density of states near the tip and crosses toward the neutral $T^2$ behavior, provided the linear dispersion remains valid.

Bulk phonons and the surface-to-bulk comparison. For the topological insulator model, assume bulk electronic excitations are frozen out while the conducting surface remains active. In three dimensions, acoustic phonons have a density of modes proportional to ${\omega }^2$. Normalize the Debye spectrum to $3N$ modes:The thermal part of the energy isAt $T\ll {\Theta }_D$, extending the upper integration limit to infinity introduces an exponentially small error. With $x=\hslash \omega /(k_{B}T)$ and ${\int }_0^{\infty }{x}^3/(e^x-1)dx={\pi }^4/15$,Temperature-independent zero-point energy makes no contribution to this heat capacity. The requested ratio consequently scales asFor a fixed shape and composition, $A\propto L^2$ and $N\propto L^3$, so the ratio grows as $1/(LT)$. Smaller samples and lower temperatures favor the surface contribution. It dominates if the dimensionless right-hand side, including its state-counting prefactor, exceeds unity. This is the surface-electron and bulk-phonon heat-capacity crossover.

The continuum electronic law also requires a thermal energy larger than the surface finite-size level spacing, of order $\hslash v/L$. The phonon continuum has its own finite-size cutoff. Therefore dominance obtained by formally taking $T\to 0$ cannot always be inferred from the two continuum power laws alone.

There is a numerical-normalization qualification in the paper. Combining the explicit one-positive-branch result above with the conventional $3N$-mode Debye normalization gives ratio prefactor $15g\zeta (3)/(16{\pi }^5)$; including electron holes doubles it. For $g=1$, this is approximately $0.00368$ before that doubling, rather than the supplied $5/(12{\pi }^2)\simeq 0.0422$. The problem specifies neither a cone/spin multiplicity that equates these nor another mode normalization. Its supplied prefactor is therefore used as the stipulated approximation for the requested numerical estimate, rather than claimed as an exact consequence of the single-branch count.

Numerical estimate. A cubic sample of volume $1\mathrm{m}{\mathrm{m}}^3$ has side $L=10^{-3}\mathrm{m}$. Counting the conducting surface on all six faces givesFor one atom per cubic cell of side $a=1.5\times 10^{-9}\mathrm{m}$,Use $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$ and $\hslash =1.05457\times 10^{-34}\mathrm{J}\mathrm{s}$. Inserting the stated $v=5\times 10^5\mathrm{m}{\mathrm{s}}^{-1}$, ${\Theta }_D=500\mathrm{K}$ and $T=1\mathrm{K}$ into the paper's supplied prefactor givesThus bulk phonons dominate at this temperature in this estimate. If only one face contributes, the ratio is $1.22\times 10^{-6}$; if there are $\nu$ chemical atoms per unit cell, divide the corresponding estimate by $\nu$. These geometrical and crystallographic factors are not supplied explicitly.

With the same prefactor and geometry, extrapolating $1/T$ would suggest equality near $7.3\mu \mathrm{K}$. However, the surface finite-size scale is already of order $\hslash v/(k_{B}L)\simeq 3.8\mathrm{mK}$ for this specimen. That extrapolated crossover lies outside the continuum electronic regime, so it is not a reliable prediction of the sample's ultimate lowest-temperature heat capacity.

Define transfer gained by the neutron as $\mathbf{Q}={\mathbf{k}}_f-{\mathbf{k}}_i$. If it is perpendicular to ${\mathbf{k}}_i$, expansion of the free-neutron kinetic energy gives $E_f-E_i={\hslash }^2{Q}^2/(2m_n)$. Absorbing one phonon requires this positive energy to equal $\hslash \omega$. A periodic phonon branch may therefore intersect a nonperiodic recoil parabola in an extended-zone representation. The sign of $\mathbf{Q}$ is opposite to the common sample-transfer convention and must be declared.

One linear Bogoliubov branch $E_k=\hslash c_{s}k$ in three dimensions gives thermal energy ${\pi }^{2}V(k_{B}T{)}^4/(30{\hslash }^3{c}_s^3)$ and the displayed cubic heat capacity. There is one branch, rather than the three polarizations in a simple crystal's Debye model. The regime requires $\hslash c_s/L\ll k_{B}T\ll \mu$ for a box: the lower inequality permits a continuum integral, and the upper one selects phonons. Fixed finite boxes ultimately have discrete-level behaviour.

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

For electron energy transfer $\Delta$ and phonon energy ${\Omega }_q=\hslash {\omega }_q$, the symmetrized pair amplitude is $V_{\mathrm{eff}}=2|M_q{|}^2{\Omega }_q/({\Delta }^2-{\Omega }_q^2)$. It is attractive for $|\Delta |<{\Omega }_q$. This interaction is retarded, and screened Coulomb repulsion must also be considered.

The displacement eigenvector of a lattice normal mode, including relative amplitudes and phases for the atoms in a primitive cell. Longitudinal and transverse describe its orientation relative to the phonon wavevector in simple isotropic cases.

Phonon scattering changes the occupation or propagation of lattice-vibration normal modes through interactions with boundaries, defects or other phonons. Anharmonic terms permit phonon-phonon scattering. Different processes need not relax heat flow equally: normal processes conserve total crystal momentum within the first Brillouin zone, whereas Umklapp scattering changes its representative by a nonzero reciprocal lattice vector.

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.

In a semiconductor, a conduction electron fills a vacant valence state and emits a photon. In a simple nondegenerate model the recombination rate is $Bnp$, while the equilibrium reverse generation rate is $B{n}_i^2$; their difference gives the displayed net pair-loss rate. This subtraction enforces detailed balance at thermal equilibrium. Photon momentum is small on the scale of a Brillouin zone, so a direct band gap permits a direct transition while an indirect band gap generally requires a phonon or another source of momentum.

A recoil-free nuclear transition emits or absorbs a gamma ray while leaving the internal vibrational state of the surrounding solid unchanged. It belongs to the elastic part of the process: no phonon is created or destroyed. The photon still carries momentum $p=E_{\gamma }/c$, where $c$ is the speed of light. The emitting solid takes momentum $-\mathbf{p}$; the absorbing solid takes $+\mathbf{p}$. Conservation of momentum holds for each event.

Why does an isolated atomic nucleus pose a problem? Start with an excited state at rest and let $E_0$ be its transition energy, with nuclear mass $M$ and $E_0\ll M{c}^2$. Conservation of momentum gives a recoil momentum of magnitude $p$; its nonrelativistic kinetic energy is thereforeBy conservation of energy, emission spends part of $E_0$ on nuclear recoil, whereas absorption must supply both $E_0$ and the final kinetic energy. To leading order,Thus the relevant comparison is $2E_R$ against the overlapping spectral line widths, rather than $E_R/E_0$ alone. If $2E_R$ greatly exceeds the widths, the isolated-nucleus lines have negligible overlap. Finite Lorentzian line shape tails mean that this is suppression, rather than an absolute prohibition on absorption.

For a bound atomic nucleus, the recoil-free channel transfers momentum to the center of mass of the emitting or absorbing solid, with negligible kinetic energy $p^2/(2M_{\mathrm{solid}})$. Other channels exchange phonons and carry substantial energy into or out of the internal vibrations. A large $M_{\mathrm{solid}}$ makes the elastic channel's energy cost small; the vibrational quantum state determines its probability. Merely replacing $M$ by a large mass in the free-particle formula does not calculate that probability.

Here is a deliberately stipulated teaching system, independent of historical isotope data: $E_0=20\mathrm{keV}$, $M{c}^2=10^{11}\mathrm{eV}$, and a single-line natural linewidth ${\Gamma }_E=10^{-7}\mathrm{eV}$ measured as an energy full width at half maximum. These give $E_R=0.002\mathrm{eV}$ and $2E_R/{\Gamma }_E=40,000$. A solid of $10^{12}$ such nuclear masses has elastic recoil energy $2\times 10^{-15}\mathrm{eV}$, far below ${\Gamma }_E$. In this article ${\Gamma }_E$ is an energy width; the corresponding lifetime is $\tau =\hslash /{\Gamma }_E$, with $\hslash$ the Planck constant divided by $2\pi$.

The residual resistivity is the low-temperature electrical resistivity remaining in a normal metal when thermally excited scattering becomes negligible. Static disorder, impurities and boundaries can relax carrier momentum without requiring a thermal population of phonons. Its value depends on sample purity and structure; a superconducting transition is a different phenomenon.

A Schottky anomaly is a peak in heat capacity caused by thermal occupation of a small set of separated energy levels. For $N$ independent two-level units with equal degeneracies and gap $\Delta >0$, write $x=\Delta /(k_{B}T)$; their heat capacity is $C=N{k}_B{x}^2{e}^{-x}/(1+e^{-x}{)}^2$, with Boltzmann constant $k_B$. It vanishes at both low and high temperature and can spoil a lattice $T^3$ fit in between. Localized electronic or nuclear levels can contribute, including in calorimeter addenda; their presence must be checked rather than attributed automatically to phonons.

Nevill Francis Mott's 1977 lecture discusses silicon deposited in a glow discharge. Injected electrons provided a drift-mobility readout; electrical conductivity and temperature dependence tested localized-state transport and excitation toward a mobility edge. Varying temperature distinguishes transport regimes. The lecture gives no contact geometry or timing circuit here, so none is inferred. Its Figure 2 is schematic, not measured data.^[ref]

This evidence motivates a transport question beyond the closed-chain eigenvalue problem: how does a prepared charge disturbance move under the conditions of a real sample? The numerical figures above supply neither drift-mobility measurements nor a model of the deposition process. Finite-size concentration alone cannot reproduce an experimental electrical conductivity curve.

Finite temperature permits phonon-assisted hopping between localized eigenstates or thermal excitation toward a mobility edge; localized states therefore need not imply zero measured electrical conductivity at every nonzero temperature.^[ref]

Real dephasing destroys the relative wave phases needed for long-path interference. Observing coherent confinement requires the distance over which those wave phases survive to be long enough to explore the localized region, and the sample to be large enough to resolve it. Temperature alone is not an automatic switch: the actual inelastic processes set the relevant scales. Static elastic wave scattering preserves coherence within each realization. Loss of the ensemble-averaged coherent field through random wave phases is therefore not, by itself, physical dephasing of an individual sample.