Past exam of the thermal and statistical physics course of the University of Cambridge / 2018 / Section B

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.

Classical translational statistics. Write the number density as $n=N/V$ and the thermal de Broglie wavelength as $\lambda =\sqrt{2\pi {\hslash }^2/(m{k}_{B}T)}$. Quantum exchange is negligible when wave packets rarely overlap:The equivalent state-counting statement is that the mean occupation of each thermally relevant one-particle state is much smaller than one. The equations in this question count one internal species; with $g$ equally populated spin states the relevant density per species is $n/g$.

Let $z=e^{\beta \mu }$ be the dimensionless quantum gas fugacity, with the lowest one-particle energy chosen as zero. Both quantum distributions can be writtenIf $z\ll 1$, then $n_{\eta }(E)\simeq z{e}^{-\beta E}\ll 1$ for all $E\ge 0$, recovering Maxwell–Boltzmann statistics. Conversely, low occupation even at the spectrum's bottom requires $z/(1-z)\ll 1$ for bosons or $z/(1+z)\ll 1$ for fermions; either condition impliesA single high-energy state can be sparsely occupied even in a degenerate gas, so the condition must hold across the relevant spectrum, not just for one arbitrarily chosen energy. Also, classical translational statistics does not imply that all internal excitation energies or other quantum effects can be ignored.

Grand partition function, density and pressure. Noninteracting particles have independent occupations of each one-particle state. Put $x_{\mathbf{k}}=z{e}^{-\beta {\epsilon }_{\mathbf{k}}}$. Summing the allowed occupations givesThereforeThe grand potential is ${\Omega }_G=-k_{B}T\ln \Xi =-pV$ for this homogeneous gas. Replacing the sum by $V\int d^{3}k/(2\pi {)}^3$ and using ${\epsilon }_{\mathbf{k}}={\hslash }^2{k}^2/(2m)$ givesThe mean number is $N=z{\partial }_z\ln \Xi$, givingFor bosons, the grand sum requires $z<1$ when the ground energy is zero. A macroscopic condensate at the limiting $z\to 1$ must be separated from the continuum integral in the condensed phase. The following small-$z$ calculation is in the normal, nondegenerate regime, so that complication does not enter.

Density expansion and inversion. Expand the occupation as a geometric series:The Gaussian wavevector integral for its $\ell$th term isIt follows thatFor the inverse series, put $y=n{\lambda }^3$ and try $z=y+c{y}^2+O(y^3)$. Substitution gives $0=(c+\eta /2^{3/2})y^2+O(y^3)$, henceThe leading equality $n{\lambda }^3\simeq z$ directly connects the density and occupancy criteria for classical behavior.

Pressure and its first quantum correction. Expanding the logarithm in the pressure expression gives one extra factor $1/\ell$ compared with the density series:Insert the inverse series. To second order in $y=n{\lambda }^3$,ThusThe first quantum correction at fixed density and temperature isIt is negative for bosons and positive for fermions, and vanishes in the classical density limit.

Virial coefficient. Define the number-density virial expansion byThe coefficient $B_j$ has units ${\mathrm{volume}}^{j-1}$ in this convention. Comparing terms givesThis is the exchange contribution to the second virial coefficient. It arises without a pair-interaction potential; the underlying gas is noninteracting. A molar-density convention would multiply this $B_2$ by the Avogadro constant. For $g$ equally populated internal species, exchange occurs within each species and the coefficient in terms of total density is smaller by a factor $g$.

Why the two signs differ. Bosonic exchange symmetry enhances the probability of nearby identical bosons, whereas fermionic exchange symmetry suppresses nearby identical fermions in the same internal state. In the dilute regime, the normalized equal-time pair detection correlation isThe exchange part follows by squaring the Gaussian one-particle coherence factor $e^{-\pi r^2/{\lambda }^2}$. At zero separation this gives $g^{(2)}(0)=2$ for bosons and $0$ for one-component fermions. Integrating the exchange pair weight provides a direct coefficient check:Bunching therefore acts like an effective statistical attraction in the dilute equation of state, lowering the pressure needed at fixed density. Exclusion acts like an effective statistical repulsion, raising it. Neither description asserts an actual interparticle force. Different internal fermion species do not exclude each other from a common spatial location in the same way, explaining the importance of the species convention.

Order parameter and the free-energy landscape. In Landau theory, take a scalar magnetization component at fixed volume and composition. It is an order parameter: its nonzero spontaneous value distinguishes an ordered phase with ferromagnetism from the symmetric high-temperature phase. One first constructs a constrained order-parameter free energy at specified magnetization by summing over microscopic states compatible with it. An external magnetic field is an independent control, so it can also appear in that constrained function. For a total magnetic moment the conjugate term is $-BM$; if $M$ is a magnetization density, the total-energy term is $-VBM$, or $-BM$ in free-energy density. Minimizing over $M$ then yields the equilibrium thermodynamic potential as a function of the imposed field.

Thus treating both $M$ and $B$ as arguments of a Landau function does not treat them as independent natural variables of the already-minimized equilibrium free energy. The first is the fluctuating or constrained order parameter, and the second its conjugate field.

At zero field, reversal symmetry requires $F(M)=F(-M)$, excluding odd powers. Near a continuous transition, assume the homogeneous constrained free energy is analytic for small $M$, that the quadratic coefficient crosses zero linearly, and that the quartic coefficient is positive. ThenThe temperature $T_c$ is where the zero-field quadratic coefficient vanishes and spontaneous ordering begins in this approximation. Higher powers are negligible sufficiently close to it; spatial fluctuations are neglected in this mean-field treatment.

Stationarity gives $2M[b(T-T_c)+2d{M}^2]=0$. Above $T_c$, the only minimum is at zero. Below it, zero becomes a maximum and two minima appear:The two signs describe spontaneous symmetry breaking: they are the average magnetizations in the symmetry-broken pure phases. A strictly finite system at exactly zero field samples the two wells symmetrically and has ensemble average $⟨M⟩=0$. A spontaneous pure-phase value is defined by selecting a phase, for example taking the thermodynamic limit before removing a small field. This qualification distinguishes a minimum's location from the average over a symmetric two-well distribution.

Minimum free energy and the heat-capacity jump. At a nonzero minimum, $M^2=-b(T-T_c)/(2d)$. Substituting this value givesHencewhile $F_{\min }=F_0$ above the transition. Take the temperature derivatives only after minimizing, so the equilibrium order parameter is allowed to change with temperature. With $C_0=-T{F}_0^{\prime \prime }(T)$,ThereforeThis is the heat-capacity jump in quartic Landau theory. The entropy is continuous at $T_c$, so there is no latent heat. If $F$ denotes free-energy density rather than total free energy, these expressions denote heat capacity per unit volume with the same normalization.

Stabilizing a tricritical expansion. At a tricritical point, the quartic coefficient is tuned to zero at the same point where the quadratic coefficient changes sign. Reversal symmetry still excludes odd powers. The next available stabilizing term is thereforeA negative or vanishing leading sextic coefficient would not stabilize this truncation for large $|M|$; a higher positive term would then be needed and the present tricritical form would be insufficient. The ordinary quartic theory has a continuous transition for $d>0$, while a stabilized theory with $d<0$ can have a first-order transition. Their meeting at zero field is what makes the simultaneous tuning tricritical.

The tricritical minimum. At zero field and $T<T_c$, stationarity givesThe nonzero solutions are minima: their curvature is $2b(T-T_c)+30d^{\prime }{M}^4=8b(T_c-T)>0$. Put $t=T_c-T>0$. Since $d^{\prime }{M}^6=(bt/3)M^2$ at stationarity,ThusThe order parameter now vanishes as $(T_c-T{)}^{1/4}$ rather than the quartic theory's square root.

Heat-capacity singularity. Define $K=2b^{3/2}/\sqrt{27d^{\prime }}>0$, so the ordered-side singular free energy is $F_{\mathrm{sing}}=-K(T_c-T{)}^{3/2}$. Its derivatives areConsequentlyThis is the tricritical heat-capacity amplitude in a sextic Landau expansion. It is positive on the ordered side and diverges as $T\to T_c^{-}$, rather than remaining a finite jump. This value of the heat-capacity critical exponent is the mean-field result along the path keeping the quartic term zero; a detuning to $d>0$ ultimately restores ordinary quartic behavior. Fluctuation corrections to the mean-field theory are not included in the stipulated polynomial minimization.