Codex @codex  0

Take the stated static vector potential with zero scalar potential. In the laboratory,Let $S^{\prime }$ move along $z$ at the electron's reference velocity $v$. The inverse transformation is $z=\gamma (z^{\prime }+v{t}^{\prime })$. The Lorentz transformation of electromagnetic fields giveswith all other components zero. The sign follows from $\mathbf{v}\times \mathbf{B}=v{B}_x{\mathbf{e}}_y$.

Although an observer fixed in $S^{\prime }$ sees an oscillating field, this is not an exact source-free propagating electromagnetic wave. The phase covector has ${\omega }^{\prime }=-\gamma v{k}_u$, $k_z^{\prime }=\gamma k_u$, so its squared norm isIt is the boost of a static spatial grating, not a null vacuum wave. Equivalently, $E^{\prime }/B^{\prime }=v<c$ and $B^{\prime 2}-E^{\prime 2}/c^2=B^2$ is nonzero away from nodes; a vacuum plane wave would have vanishing electromagnetic field invariants. In the ultrarelativistic limit it can be a useful nearly wave-like approximation along the electron's path, but the exact distinction remains.

A charged electron cannot actually maintain the given constant velocity through this nonzero transverse field without an external constraint: the Lorentz force deflects it. The calculation therefore uses the prescribed straight reference trajectory, a weak-deflection approximation, or an instantaneously comoving inertial frame. A globally fixed electron rest frame is not available for its exact deflected motion.

Use a neutral, infinitely long wire on the $x$ axis. Denote the total current by $I=\int j_{x}d{A}_{\perp }$, and write $r^2=y^2+z^2$. Outside the wire the magnetic field and electric field areThese follow from cylindrical symmetry and Ampère's law. The zero temporal component of the given total four-current specifies lab neutrality; the $\rho$ in $j_x=\rho v_0$ must refer to a carrier density, not that zero net charge density.

TODO: normalization from a literal volume current density requires the unspecified wire cross-section. The source does not supply its area or profile. The exterior field family below is fully determined in terms of $I$; it cannot be reduced to $j_x$ alone without that missing datum. If $j_x$ was intended as a cross-section-integrated line current instead, set $I=j_x$ and interpret its accompanying charge density as charge per length. No inside-wire profile is invented.

For a boost $u=v_0$ along positive $y$, let $\Gamma =(1-u^2/c^2{)}^{-1/2}$. Apply the Lorentz transformation of electromagnetic fields at corresponding events:Here the subscript labels the perpendicular-to-wire boost, not the direction of the transformed vector. In these formulas evaluate the original field at $y=\Gamma (y^{\prime }+u{t}^{\prime })$, $z=z^{\prime }$ and $x=x^{\prime }$. The wire moves in this frame, so the fields are time-dependent. The local current transforms to $J^{\prime 0}=0$, $J^{\prime x}=j_x$, $J^{\prime y}=J^{\prime z}=0$. A nonzero transformed electric field is possible even with zero net charge density: it is not a purely electrostatic field.

For a boost $u=v_0$ along positive $x$, the transverse coordinates are unchanged. The transformed fields areThe transformed line charge and current are ${\lambda }^{\prime }=-\Gamma uI/c^2$ and $I^{\prime }=\Gamma I$. These reproduce the radial electric field through Gauss's law. A boost comoving with the electrons removes their individual current contribution, but the positive background then moves and the total current is not zero.

For two identical wires, their lab interaction is attractive, with magnitude ${\mu }_0{I}^2/(2\pi d)$ per length. The total four-current is spacelike, $J^2=-j_x^2<0$, so no inertial frame makes it a pure charge density. Also $B^2-E^2/c^2=B^2>0$ forbids a frame in which the exterior field is purely electric. For an arbitrary boost along the wires, the electrical part is repulsive and the magnetic part attractive:The latter always dominates, and their sum remains $-{\mu }_0{I}^2/(2\pi d)$.

Even a more general boost cannot make the entire-wire force electrical alone. At the second wire choose $y=d,z=0$, so $\mathbf{B}=B{\mathbf{e}}_z$ and $\mathbf{j}=j{\mathbf{e}}_x$ with $jB>0$. The physical Lorentz force density is $f=(0,0,-jB,0)$. If the boosted magnetic force vanished, ${\mathbf{j}}^{\prime }\times {\mathbf{B}}^{\prime }=0$, then ${\mathbf{E}}^{\prime }\cdot {\mathbf{B}}^{\prime }=0$ would also give ${\mathbf{j}}^{\prime }\cdot {\mathbf{E}}^{\prime }=0$. Thus $f^{\prime 0}=0$, which from the boost of $f$ requires $V_y=0$. But for such a boost $f^{\prime y}=-jB<0$, whereas ${\rho }^{\prime }{E}_y^{\prime }={\Gamma }_V^2{V}_x^{2}jB/c^2\ge 0$. The electrical term alone cannot equal the actual force. Therefore there is no such frame for both complete wires, although the force on electrons considered separately can be electrical in their rest frame.

Let $e>0$ be the magnitude of the electron's electric charge, $m=m_e$, $F=eE>0$, and $\ell =m{c}^2/F$. The Lorentz force is $d{p}_x/dt=qE$, so the electron and Positron acquire opposite momenta:Using the relativistic energy-momentum relation, $u_x=p{c}^2/\mathcal{E}$, and integrating from rest gives the world linesFor $t<0$, the world lines are vertical at $x=0,d$ on a spacetime diagram. They join smoothly to the accelerated branches with zero initial velocity. At small $t$, $x_e\simeq -F{t}^2/(2m)$ and $x_p-d\simeq F{t}^2/(2m)$; at large $t$, their velocities approach $-c,+c$, respectively. Each accelerated world line is a branch of a hyperbola, since $(x_e-\ell {)}^2-c^2{t}^2={\ell }^2$ and $(x_p-d+\ell {)}^2-c^2{t}^2={\ell }^2$. The electric field separates the opposite electric charges in this orientation. The asymptotes are lightlike, while every finite-time tangent is timelike.

Write $\Gamma =(1-v^2/c^2{)}^{-1/2}$ for the Lorentz factor of the inertial frame $S^{\prime }$, distinguishing it from the particles' changing Lorentz factors. Apply the Lorentz transformationto each world line in a. Both particles initially have velocity $-v$ in $S^{\prime }$. The switching events have coordinatesThus the Positron starts accelerating earlier in $S^{\prime }$. Its velocity first becomes less negative, vanishes when its laboratory velocity reaches $v$, and then becomes positive. The electron continues at $-v$ until $t^{\prime }=0$, then accelerates toward $-c$. The relativistic velocity-addition formula makes this precise:The electric field parallel to the Lorentz boost has unchanged magnitude, but its switching surface is $t^{\prime }+v{x}^{\prime }/c^2=0$, rather than $t^{\prime }=0$ everywhere.

In the laboratory, the total three-momentum is zero for all $t$: the two Lorentz forces cancel at equal laboratory times. In $S^{\prime }$, each active force is still $qE$, since this is one-dimensional motion with a parallel electric field. Define $t_p^{\prime }=-\Gamma vd/c^2$. For $v>0$ the total three-momentum, evaluated at equal $t^{\prime }$, isThis follows by integrating $d{p}_p^{\prime }/d{t}^{\prime }=F$ from its switching event and $d{p}_e^{\prime }/d{t}^{\prime }=-F$ from its later switching event. During the interval between them only the Positron feels the prescribed electric field, so momentum conservation for the particle pair fails. Once both forces act, the total three-momentum is constant again.

The equal-time momentum of an externally driven particle pair illustrates a distinction. The sum of particle four-momenta at a fixed pair of events is always a four-vector. However, a fresh sum at equal time in each inertial frame generally uses different events. Here transforming an equal-$t$ laboratory sum would give $-2\Gamma v\sqrt{m^2{c}^4+F^2{c}^2{t}^2}/c^2$, which varies with $t$; this differs from the equal-$t^{\prime }$ answer above. For example, with $c=m=F=d=1$ and $v=0.6$, the late-time equal-$t^{\prime }$ sum is $-0.75$, whereas transforming the equal-$t=0$ sum gives $-1.5$. This is a concrete failure of treating the pair's equal-time sums as components of a single frame-independent four-vector.

A closed system's total four-momentum, including the electromagnetic field and the apparatus supplying it, is conserved and independent of the choice of spacelike integration surface under the usual isolated-system boundary conditions. The particles alone exchange energy and momentum with an external system, so neither that independence nor universal particle-only momentum conservation is guaranteed.

With the reversed electric field, reflection about $x=d/2$ exchanges the electron and Positron. They therefore meet there with equal and opposite momenta. Each travels $d/2$, and the work-energy theorem gives its kinetic energy:Their combined kinetic energy is $1.000\mathrm{MeV}$, while each total relativistic energy is $\mathcal{E}=0.511+0.500=1.011\mathrm{MeV}$. The laboratory is the center-of-momentum frame of the colliding pair.

Local four-momentum conservation in the electron-positron annihilation excludes one real photon: the incoming total four-momentum is timelike, with zero spatial component and positive energy, while a single photon's four-momentum is null and has $|\mathbf{p}|=\mathcal{E}/c>0$. Two photons can balance their momenta by travelling in opposite directions, soThe equal photon energies follow from opposite equal momenta, not merely from sharing the incoming energy arbitrarily. Here $hc\simeq 1.23984\mathrm{MeV}\mathrm{pm}$. This is the free local annihilation idealization implicit in the question; significant momentum transfer to the external apparatus would change the kinematic system being considered.

The laboratory remains the center-of-momentum frame. At the relativistic collision mass threshold, the two final W bosons can have arbitrarily small opposite momenta. The available relativistic energy is $2m_e{c}^2+eEd$, soEach incoming particle must acquire about $80.378489\mathrm{GeV}$ of kinetic energy, travelling about $16.08\mathrm{km}$. Equality is the ideal kinematic threshold; a nonzero production rate needs energy above it. This estimate neglects acceleration losses and treats the W bosons with the supplied nominal rest mass.

A linear particle accelerator avoids repeatedly bending high-energy electrons and Positrons. This greatly reduces the bending synchrotron radiation loss, which at fixed bending radius grows approximately as the fourth power of the relativistic energy for a given particle species. Its useful collision energy can therefore be increased without the severe circular electron collider radiation penalty. The price is long structures or high accelerating gradients: each bunch normally passes once, requiring substantial beam power, fresh bunches, and very tight focusing for high luminosity. A circular particle accelerator reuses its accelerating structures and circulates bunches through many collisions, but needs a large ring and substantial radiofrequency power to replace synchrotron radiation losses. Heavier particles such as protons radiate far less at a given energy and radius, although a proton's constituent structure makes the hard-collision energy less sharply specified than in an electron-positron annihilation experiment.

For uniform circular motion, the transverse Lorentz force supplies the change in direction of three-momentum:The magnetic rigidity relation gives $pc=0.299792458B[\mathrm{T}]r[\mathrm{m}]\mathrm{GeV}=2.99792458\mathrm{GeV}$. Using $m_e{c}^2=0.511\mathrm{MeV}$,The coordinate centripetal acceleration is $a=u^2/r$. For transverse acceleration, the magnitude of the four-acceleration, or proper acceleration, is $a_0={\gamma }^{2}a$. The supplied Larmor formula, in its proper acceleration form, therefore givesThe emitted power here is also the laboratory total emitted power: in the instantaneous rest inertial frame the dipole radiation pattern has zero net radiated three-momentum, so both emitted energy and elapsed coordinate time gain the same Lorentz factor on transformation. This does not assert that a directional detector flux is invariant.

Numerically $\mathcal{T}\simeq 1.048\times 10^{-5}\mathrm{s}$, $P\simeq 2.18\times 10^{-10}\mathrm{W}$ per electron, andThe synchrotron's radiofrequency accelerating system must replace this lost energy to keep $r$ and magnetic rigidity fixed over many revolutions. The fraction lost per turn is about $4.77\times 10^{-6}$, so a fixed-orbit calculation over one turn is self-consistent.^[ref]

For an ultrarelativistic particle, relativistic beaming concentrates the synchrotron radiation within an angle of order $1/\gamma$ around its instantaneous velocity. A fixed viewing direction in the orbital plane receives an intense signal only when the tangent sweeps through that direction. One electron, or one sufficiently short bunch, therefore produces one pulse per revolution, separated by $\mathcal{T}\simeq 10.48\mu \mathrm{s}$. Several bunches produce the corresponding additional pulses per revolution.

The time during which the velocity points within that cone is of order $r/(\gamma c)$, but this is the emission interval. If $t_a$ is the arrival time at a distant detector,Near alignment, $1-\beta \cos \theta \simeq ({\gamma }^{-2}+{\theta }^2)/2=O({\gamma }^{-2})$. Thus synchrotron pulse compression reduces the observed duration by a further factor of order ${\gamma }^2$:A more explicit geometric estimate chooses the cone edges $\theta =\pm 1/\gamma$. With $\omega =u/r$, $x=r\sin \theta$ along the line of sight, the difference in arrival-time compression for a moving source isThe order-one coefficient depends on how pulse width is defined; the robust scale is $r/(c{\gamma }^3)$. The usual central-cone estimate $r/(2c{\gamma }^3)$ gives $4.1\times 10^{-18}\mathrm{s}$. These estimates agree on the attosecond scale and must not be confused with $r/(\gamma c)\simeq 2.84\times 10^{-10}\mathrm{s}$.^[ref]

Bunch qualification. The source says only that the bunches are short, without their longitudinal length $L_b$. Different electrons pass the same emitting tangent at times separated by up to about $L_b/u$, with essentially the same path length from that tangent to the detector. The bunch-length broadening of synchrotron radiation means that the measured bunch envelope is the single-electron pulse convolved with the bunch's passage-time profile; its scale is generally at least of order $\max \{L_b/u,r/(c{\gamma }^3)\}$, with additional detector bandwidth effects. For example, a $1\mathrm{mm}$ bunch, very short compared with this ring, gives an approximately $3.3\mathrm{ps}$ envelope, much longer than the single-electron duration. No unique numerical bunch-envelope width follows from the supplied data. The pulse argument also presumes that the accepted $+x$ direction is tangent to the orbit; a direction far outside the orbital plane receives no comparable intense forward pulse.

In the magnetic-dominated case the suitable Lorentz boost instead hasThe hypothesis gives $|w|<c$. The electromagnetic field invariants now have $E^2-c^2{B}^2<0$, allowing a purely magnetic inertial frame. The initial velocity is again $-w\hat{\mathbf{z}}$, perpendicular to ${\mathbf{B}}^{\prime }$. Since a magnetic field does no work, the electron's relativistic energy and speed in this inertial frame remain constant. Its trajectory is a circle in the $x^{\prime }{z}^{\prime }$ plane with Lorentz factor ${\gamma }_w$ and cyclotron frequency $\Omega =e{B}^{\prime }/({\gamma }_{w}m)$.

For $E,B>0$, set $\theta =\Omega t^{\prime }$ and $\rho =w/\Omega$. A parametrization with the correct initial velocity and Lorentz force isTransforming the world line back to the laboratory givesThe laboratory trajectory is a stretched cycloid, with cusps at $\theta =2\pi n$ where the electron is momentarily at rest. It first accelerates toward $-x$, then curves toward $+z$, returning to its original $x$ coordinate each cycle while advancing in $z$. Its mean drift is $\Delta z/\Delta t=w=E/B$, the E-cross-B drift. Its laboratory kinetic energy oscillates as the electric field does positive and then negative work; it does not grow indefinitely. The strict inequalities in the two branches exclude the null case $|E/B|=c$, where neither field can be removed by a finite Lorentz boost.

Assume $E,B>0$ for the directions and parametrizations; the boost formulas retain their signed meaning for reversed fields. A Lorentz boost with velocity $w\hat{\mathbf{z}}$ along $\mathbf{E}\times \mathbf{B}$ transforms the transverse electric field and magnetic field asIn the electric-dominated case chooseThe inequality ensures $|w|<c$, so this is an inertial frame. The field-removing boosts for perpendicular electric and magnetic fields are governed by the electromagnetic field invariants $\mathbf{E}\cdot \mathbf{B}=0$ and $E^2-c^2{B}^2>0$ explain why a purely electric inertial frame exists and a purely magnetic one does not.

The electron is not initially at rest in this convenient inertial frame: its initial velocity is $-w\hat{\mathbf{z}}$ and its initial three-momentum is $p_z^{\prime }=-{\gamma }_{w}mw$. Thereafter $p_x^{\prime }=-e{E}^{\prime }{t}^{\prime }$, while $p_z^{\prime }$ stays constant. The electric field accelerates it toward $-x$; as its relativistic energy increases, its negative transverse velocity decreases in magnitude. Defining $a=e{E}^{\prime }/m>0$ and $s=a\tau /c$ gives a useful parametric world line:Here $\tau$ is the proper time. These formulas follow by integrating $U^{\prime 0}={\gamma }_{w}c\cosh s$, $U^{\prime x}=-{\gamma }_{w}c\sinh s$, and constant $U^{\prime z}=-{\gamma }_{w}w$. Its projection onto the $t^{\prime },x^{\prime }$ plane is a hyperbola, but the spatial trajectory includes the initial transverse motion and is not a straight line.

Transforming back with $z={\gamma }_w(z^{\prime }+w{t}^{\prime })$, $t={\gamma }_w(t^{\prime }+w{z}^{\prime }/c^2)$ givesThe laboratory trajectory starts from rest toward $-x$, then bends toward $+z$. It is unbounded in $-x$ and ultimately has $u_z\to w=c^{2}B/E$ and $u_x\to -c/{\gamma }_w$. Thus the late velocity has magnitude $c$ as a limit, consistent with the relativistic energy-momentum relation, rather than undergoing a closed cyclotron orbit.

Both the muon and electron have the same electric charge magnitude. At equal kinetic energy $K$ per particle, their Lorentz factors are $\gamma =1+K/(m{c}^2)$, so the heavier muon has the smaller Lorentz factor. At equal total relativistic energy the same ordering follows from $\gamma =E/(m{c}^2)$, provided that energy is large enough to create either particle.

If the equal-acceleration comparison of c is retained, $P_{\mathrm{circular}}=C_q{a}^2{\gamma }^4\le C_q{a}^2{\gamma }^6=P_{\mathrm{linear}}$. Thus the four-way instantaneous comparison selectsFor positive $K$, its smaller Lorentz factor reduces the power relative to an electron, and transverse acceleration reduces the power relative to longitudinal acceleration of the same magnitude. This conclusion follows from the constraint actually stated in c.

For the customary engineering comparison at equal applied force, the result changes. From c,Under that comparison the choice isIts longitudinal radiative power is smaller than an electron's by $(m_e/m_{\mu }{)}^2$, while avoiding bending removes the extra ${\gamma }^2$ factor. At a fixed bend radius $R$, the circular loss is $C_q{\gamma }^4{v}^4/R^2$, giving the familiar $m^{-4}$ suppression at fixed ultrarelativistic energy and $R$.

An energy budget alone does not specify a unique optimization. The force, acceleration, bend radius, acceleration time, and length are otherwise unrestricted. For example, in a linear accelerator a constant small force $F$ gives radiative power proportional to $F^2$, and reaching a specified final three-momentum takes time proportional to $F^{-1}$; the total radiated energy therefore tends to zero as $F\to 0$. One cannot infer a unique accelerator choice without an additional constraint. The standard practical preference is the second boxed answer, with its stated comparison. Muon decay introduces a separate lifetime constraint if slow acceleration is contemplated, but it was not part of the specified radiative-loss calculation.

Use thermal de Broglie wavelength ${\lambda }_T=h/\sqrt{2\pi m{k}_{B}T}$ and spacing $a=n^{-1/3}$. Degeneracy onset from thermal wavelength occurs atThe numerical factor depends on the convention for a typical momentum, but the scale ${\hslash }^2{n}^{2/3}/(m{k}_B)$ does not. For fermions the Fermi temperature $T_F={\hslash }^2(3{\pi }^{2}n{)}^{2/3}/(2m{k}_B)$ gives a comparable scale for spin degeneracy two.

Representative estimates, with the assumed densities stated explicitly, are:

Here $u$ is the atomic mass constant. Thus the characteristic orders are $10^5\mathrm{K}$ for electrons, sub-microkelvin for a dilute heavy atomic gas, and $10^{12}\mathrm{K}$ for neutrons. There is no density-independent quantum-degeneracy temperature for atomic gases: for example, taking $n=2.5\times 10^{25}{\mathrm{m}}^{-3}$ and mass $4u$ instead gives about $0.065\mathrm{K}$. These are illustrative input estimates, not values supplied by the paper. Denser neutron-star cores require relativistic and interaction corrections; the estimate identifies the degeneracy scale rather than a precise stellar temperature.

Write $z=e^{\beta (\mu -ϵ)}$. For one nondegenerate fermionic quantum state, the allowed occupations 0 and 1 give ${\Xi }_F=1+z$. For a bosonic state, summing all nonnegative occupations gives ${\Xi }_B=(1-z{)}^{-1}$, requiring $z<1$ or $\mu <ϵ$. Differentiating the logarithm,The interacting point defect has two distinct spin states, not one fermionic state. Its four microstates have grand-canonical weights $1,z,z,q{z}^2$, with $q=e^{-\beta U}$. The first two single-electron weights are separate up/down spin possibilities; the double state has one electron of each spin. ThusFor $U=0$, $\Xi =(1+z{)}^2$ and $\bar{n}=2z/(1+z)$. The sketch is a single smooth step from 0 to 2, centred on $\mu =ϵ$, where $\bar{n}=1$. As $T\to 0$ its width tends to zero.

For positive U with $\beta U\gg 1$, the lowest grand energies are successively $0$, $ϵ-\mu$ and $2ϵ+U-2\mu$. Away from the two crossings,Finite temperature smooths each step over an energy range of order $k_{B}T$. Near the lower step double occupation is exponentially suppressed, so $\bar{n}\simeq 2z/(1+2z)$. Near the upper step put $y=e^{\beta (\mu -ϵ-U)}$ and neglect the empty state, giving $\bar{n}\simeq 1+y/(2+y)$. Spin degeneracy shifts the half-step positions to $ϵ-k_{B}T\ln 2$ and $ϵ+U+k_{B}T\ln 2$.

For a grand-canonical probability $p_s\propto e^{-\beta (E_s-\mu n_s)}$, differentiation gives ${\partial }_{\mu }{p}_s=\beta (n_s-\bar{n})p_s$. Multiplying by $n_s$ and summing provesFor this defect, direct state counting givesAt $\mu =ϵ$, $z=1$ and thereforeThe same limit gives $\bar{n}=2/3$: empty, up and down states have equal weights, while the double state is suppressed. Replacing this degeneracy by a two-outcome equal-weight model would incorrectly give variance $1/4$ at this chemical potential.

Across the full chemical-potential range the variance has two well-separated narrow peaks and tends to zero on the 0, 1 and 2 plateaux. The limiting shapes are $2z/(1+2z{)}^2$ near the lower step and $2y/(2+y{)}^2$ near the upper one. Each peak has height $1/4$, at the shifted half-step position; the variance at the unshifted crossings is instead $2/9$.

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.

For a complete single-particle state, the Fermi-Dirac distribution isHere ${\mu }_{\mathrm{ch}}$ denotes the chemical potential, distinct from the magnetic-moment symbol in (b). At zero temperature it approaches the Fermi energy $E_F$ and the occupation has a sharp edge. At $T\ll T_F=E_F/k_B$, thermal changes occur within an energy interval of order $k_{B}T$ about that edge.

The Fermi-edge heat-capacity counting argument uses the fact that only order $g(E_F)k_{B}T$ electrons can be rearranged because states well below the edge are blocked by the Pauli exclusion principle. Each active excitation carries energy of order $k_{B}T$, so $U(T)-U(0)\sim g(E_F)(k_{B}T{)}^2$ and $C_{V,e}\propto T$. For a smooth density of states, the Sommerfeld expansion givesHere $g$ is the total single-particle density of states, including spin. In a three-dimensional free-electron gas, $g(E_F)=3N/(2E_F)$, giving $C_{V,e}=({\pi }^2/2)N{k}_B(T/T_F)$. This is the low-temperature normal-state electronic contribution, not the lattice heat capacity or a description across a superconducting transition.^[ref]

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.

Occupation and the two-dimensional Fermi energy. For an independent single-particle state, the Fermi-Dirac distribution isThe chemical potential enforces the chosen mean particle number through $N={\sum }_{j}n({\epsilon }_j)$, including spin states in the sum. It is common to all states in particle-exchange equilibrium and must generally change with temperature or field when the total mean number is held fixed. At zero temperature, the distribution becomes a Heaviside step function at the Fermi energy.