Degenerate Fermi Gas at Zero Temperature

Cool an ideal gas of fermions toward absolute zero and something strange happens: it refuses to sit still. A classical gas would freeze into a single low-energy heap, but fermions cannot — the Pauli principle forbids two of them from sharing a single-particle state. So even at T=0T=0 the particles are forced to stack into successively higher momentum states, filling every level from the bottom up to a sharp ceiling. That ceiling is the Fermi energy εF\varepsilon_F, the filled region in momentum space is the Fermi sphere, and the whole construction is the ground state of the degenerate Fermi gas. This article establishes that ground state and its exact zero-temperature thermodynamics. We take the T→0T\to0 limit of the mean occupation from aq-sm-fermi-dirac-distribution, whose Fermi function nˉ=1/(eβ(ε−μ)+1)\bar n=1/(e^{\beta(\varepsilon-\mu)}+1) collapses to a step Θ(εF−ε)\Theta(\varepsilon_F-\varepsilon): every state below εF\varepsilon_F singly occupied, every state above empty. Then a single counting argument — how many single-particle states fit inside a sphere of radius kFk_F in kk-space — fixes everything. The density of states g(ε)∝εg(\varepsilon)\propto\sqrt\varepsilon that we plug into the moments is owned by aq-sm-density-of-states-free-particles, and the resulting energy scale εF(n)\varepsilon_F(n) is developed in aq-sm-fermi-energy-momentum-temperature; here we take both as given and read off the two owned results of this article: the filled-sphere count N=gsVkF3/6π2N=g_sVk_F^3/6\pi^2 and the exact ground-state energy E0=35NεFE_0=\tfrac35 N\varepsilon_F. The payoff is enormous. The Fermi temperature TF=εF/kBT_F=\varepsilon_F/k_B of metallic electrons reaches 10410^4–10510^5 K, so room-temperature electrons are deeply degenerate — the ground-state picture, not thermal physics, governs them. The gas stores a colossal zero-point kinetic energy and pushes back with a degeneracy pressure that supports white dwarfs against gravity, even when completely cold.

Key takeaways
  • At T=0T=0 the Fermi-Dirac occupation is a sharp step nˉ(ε)=Θ(εF−ε)\bar n(\varepsilon)=\Theta(\varepsilon_F-\varepsilon): all states below εF\varepsilon_F filled, all above empty — the Fermi sea.
  • In kk-space the filled states form the Fermi sphere of radius kFk_F; counting states inside gives N=gsVkF3/6π2N=g_sVk_F^3/6\pi^2, hence kF=(3π2n)1/3k_F=(3\pi^2 n)^{1/3} for spin-12\tfrac12 fermions (gs=2g_s=2).
  • The energy per particle is a pure number times the Fermi energy: E0/N=35εFE_0/N=\tfrac35\varepsilon_F, an entirely quantum zero-point energy that survives at T=0T=0.
  • The cold gas exerts a finite degeneracy pressure P0=25nεF=23E0/VP_0=\tfrac25 n\varepsilon_F=\tfrac23 E_0/V — a consequence of exclusion, not of temperature.
  • Key idea: the T=0T=0 ground state is the filled Fermi sphere; counting states inside it (N=gsVkF3/6π2N=g_sVk_F^3/6\pi^2) and taking moments of the ε\sqrt\varepsilon density of states gives E0=35NεFE_0=\tfrac35 N\varepsilon_F exactly — every ground-state property follows from this one geometric picture.
The core insight

At T=0T=0 the Fermi function is a step Θ(εF−ε)\Theta(\varepsilon_F-\varepsilon), so the particles fill a sphere of radius kFk_F in momentum space. Requiring the sphere to hold exactly NN particles fixes kF=(3π2n)1/3k_F=(3\pi^2 n)^{1/3} (for gs=2g_s=2), and integrating the free-particle energy ℏ2k2/2m\hbar^2k^2/2m over the filled sphere gives the exact ground-state energy E0=35NεFE_0=\tfrac35 N\varepsilon_F. Everything cold about a Fermi gas — its enormous mean energy, its degeneracy pressure, its stiffness — is the arithmetic of a filled ball.

Where this lives today

The degenerate Fermi gas is the working model for the electrons in every metal and the matter in every compact star. In a metal the conduction electrons form a Fermi sea with TF∼104T_F\sim10^4–10510^5 K, so at ordinary temperatures T/TF≪1T/T_F\ll1 and only the thin shell of width ∼kBT\sim k_BT around the Fermi surface can be thermally excited — the reason the electronic heat capacity is a small linear-in-TT term rather than the classical 32NkB\tfrac32 Nk_B, computed in aq-sm-electronic-heat-capacity-metals via the aq-sm-sommerfeld-expansion built on this exact ground state. In astrophysics, electron degeneracy pressure supports white dwarfs and neutron degeneracy supports neutron stars; the softening of the equation of state from n5/3n^{5/3} to the relativistic n4/3n^{4/3} produces the Chandrasekhar mass ≈1.4 M⊙\approx1.4\,M_\odot. And in the laboratory, trapped 40^{40}K and 6^6Li atoms cooled below their Fermi temperature realize the textbook ground state directly: experiments image the sharp Fermi surface in time-of-flight, measure the 35εF\tfrac35\varepsilon_F energy per particle, and observe Pauli blocking of collisions.

The physical picture

Pour fermions into a box one at a time. The first drops into the lowest single-particle state; the next cannot join it (beyond the gsg_s-fold spin degeneracy), so it takes the next state up; and so on. Like water filling a basin, the particles settle into the lowest available levels, and when all NN are in, the "water line" is the Fermi energy εF\varepsilon_F. Because a free particle's energy depends only on the magnitude of its wavevector, ε=ℏ2k2/2m\varepsilon=\hbar^2k^2/2m, the occupied states fill a ball in kk-space — the Fermi sphere — bounded by the Fermi surface at ∣k∣=kF|\mathbf k|=k_F. The step-function occupation is the smoking gun of degeneracy. At T=0T=0 there is no thermal energy to promote a particle across the finite gap to an empty state above the surface, so the distribution is perfectly sharp. The first figure is the hero: the T=0T=0 step, and beside it the filled Fermi sphere it describes.

Figure
The T=0 Fermi step and the filled Fermi sphere it describes
T = 0: a sharp step in energy is a filled sphere in k-space energy ε → ⟨n⟩ 1 0 ε = ε_F filled ⟨n⟩ = 1 empty k_F Fermi sphere: filled to |k| = k_F every k inside radius k_F is occupied

The punchline: a sharp step at εF\varepsilon_F in energy is exactly a solidly filled ball of radius kFk_F in momentum space. There is no gradient, no tail — the ground state is a filled sphere.

The ground-state energy is purely kinetic and purely quantum. Even with no interactions and no temperature, antisymmetry forces most particles into high-momentum states. Averaging the energy over the filled sphere gives each particle 35εF\tfrac35\varepsilon_F on average — a huge number for a "cold" gas. The second figure shows why the average is 35\tfrac35 and not 12\tfrac12: the states pile up as ε\sqrt\varepsilon, so there are far more high-energy states than low-energy ones, dragging the mean upward.

Figure
Why the mean energy is three-fifths of the Fermi energy: the sqrt-epsilon density of states weights the top heavily
g(ε) ∝ √ε: the top of the sea is crowded, so ⟨ε⟩ = (3/5) ε_F energy ε → g(ε) ∝ √ε ε_F filled at T = 0 ⟨ε⟩ = (3/5)ε_F empty above ε_F

Because g(ε)∝εg(\varepsilon)\propto\sqrt\varepsilon rises toward the surface, the mean occupied energy is pulled well above the midpoint — the exact moment ratio gives 35\tfrac35.

Density controls everything. Since εF∝n2/3\varepsilon_F\propto n^{2/3}, compressing the gas raises the Fermi energy — and hence the pressure — steeply. The cold Fermi gas resists compression far more stiffly than a classical gas at the same temperature, and that stiffness is what holds up a white dwarf.

The core equations

The starting point is the zero-temperature occupation. Taking T→0T\to0 in the Fermi function (owned by aq-sm-fermi-dirac-distribution) sends the smooth shelf to a sharp step at the chemical potential, which at T=0T=0 is the Fermi energy.

nˉ(ε)=1e(ε−μ)/kBT+1 → T→0  Θ(εF−ε),εF≡μ(T=0)(1)\bar n(\varepsilon)=\frac{1}{e^{(\varepsilon-\mu)/k_BT}+1}\ \xrightarrow{\ T\to0\ }\ \Theta(\varepsilon_F-\varepsilon),\qquad \varepsilon_F\equiv\mu(T{=}0) \qquad (1)

where: nˉ(ε)\bar n(\varepsilon) — mean occupation of a state of energy ε\varepsilon; Θ\Theta — Heaviside step; εF\varepsilon_F — Fermi energy, the highest occupied energy at T=0T=0. Physical reading: at T=0T=0 occupation is all-or-nothing, so counting occupied states is a pure geometry problem — the volume of a sphere in kk-space. Counting the states inside the Fermi sphere fixes its radius. Each single-particle kk-state occupies a cell of volume (2π/L)3(2\pi/L)^3, so the number of states with ∣k∣<kF|\mathbf k|<k_F, times the spin degeneracy, must equal NN.

N=gs V(2π)3∫∣k∣<kF ⁣d3k=gs V(2π)3⋅43πkF3=gsV kF36π2(2)N=g_s\,\frac{V}{(2\pi)^3}\int_{|\mathbf k|<k_F}\!\mathrm d^3k=g_s\,\frac{V}{(2\pi)^3}\cdot\frac{4}{3}\pi k_F^3=\frac{g_sV\,k_F^3}{6\pi^2} \qquad (2)

where: gs=2s+1g_s=2s+1 — spin degeneracy (=2=2 for electrons); VV — volume; kFk_F — Fermi wavevector. Physical reading: this is the owned counting result of the article — the filled sphere holds exactly NN particles, and inverting gives the Fermi wavevector purely from the density. Solving Eq. (2) for kFk_F gives the Fermi wavevector as a function of density alone — no temperature appears.

kF=(6π2ngs)1/3 → gs=2  (3π2n)1/3,n=NV(3)k_F=\left(\frac{6\pi^2 n}{g_s}\right)^{1/3}\ \xrightarrow{\ g_s=2\ }\ (3\pi^2 n)^{1/3},\qquad n=\frac{N}{V} \qquad (3)

where: nn — number density. Physical reading: doubling the density raises kFk_F by 21/32^{1/3}; the Fermi energy εF=ℏ2kF2/2m\varepsilon_F=\hbar^2k_F^2/2m (its density scaling developed in aq-sm-fermi-energy-momentum-temperature) then rises by 22/32^{2/3}. The ground-state energy is the sum of the kinetic energies of all filled states. Integrating ℏ2k2/2m\hbar^2k^2/2m over the filled sphere — equivalently, taking the first moment of the ε\sqrt\varepsilon density of states owned by aq-sm-density-of-states-free-particles — gives a strikingly clean result.

E0=gsV(2π)3∫∣k∣<kF ⁣ℏ2k22m d3k=∫0εF ⁣ε g(ε) dε=35NεF(4)E_0=g_s\frac{V}{(2\pi)^3}\int_{|\mathbf k|<k_F}\!\frac{\hbar^2k^2}{2m}\,\mathrm d^3k=\int_0^{\varepsilon_F}\!\varepsilon\,g(\varepsilon)\,\mathrm d\varepsilon=\frac{3}{5}N\varepsilon_F \qquad (4)

where: E0E_0 — total (zero-point) kinetic energy of the filled Fermi sphere; g(ε)=3N2ε/εF3/2g(\varepsilon)=\tfrac{3N}{2}\sqrt\varepsilon/\varepsilon_F^{3/2}. Physical reading: this is the second owned result — the energy per particle is 35εF\tfrac35\varepsilon_F, a pure number set by the geometry of a filled ball, independent of gsg_s, VV, or the density of states prefactor. Because εF\varepsilon_F carries all the volume dependence at fixed NN, the cold gas has a finite pressure. Any non-relativistic ideal gas obeys PV=23EPV=\tfrac23 E, so the ground-state energy immediately yields the degeneracy pressure.

P0=23E0V=25n εF=(3π2)2/35ℏ2m n5/3(5)P_0=\frac{2}{3}\frac{E_0}{V}=\frac{2}{5}n\,\varepsilon_F=\frac{(3\pi^2)^{2/3}}{5}\frac{\hbar^2}{m}\,n^{5/3} \qquad (5)

where: P0P_0 — zero-temperature pressure, finite despite T=0T=0. Physical reading: this "degeneracy pressure" is quantum, not thermal — its steep n5/3n^{5/3} growth is the mechanical backbone of white-dwarf stability, developed fully in aq-sm-degeneracy-pressure-fermions.

Derivations

Derivation ACounting the Fermi sphere and the exact ground-state energy

The goal is to derive the two owned results of this article — the filled-sphere particle count N=gsVkF3/6π2N=g_sVk_F^3/6\pi^2 (hence kF=(3π2n)1/3k_F=(3\pi^2 n)^{1/3} for gs=2g_s=2) and the exact ground-state energy E0=35NεFE_0=\tfrac35 N\varepsilon_F — from the single fact that at T=0T=0 the occupation is a step. The strategy is pure counting: convert the kk-space sum to an integral using the uniform mode density V/(2π)3V/(2\pi)^3 (owned by aq-sm-density-of-states-k-space), fill the sphere, and take the ratio of the energy moment to the number moment so that all prefactors cancel. Watch the one subtlety: the answer 35\tfrac35 must be independent of gsg_s, VV, and mm — it is a property of the geometry, not the material.

Step 1. Convert the sum over occupied states to an integral. At T=0T=0 every kk-state with ∣k∣<kF|\mathbf k|<k_F is filled once per spin, so N=gs∑∣k∣<kF1→gsV(2π)3∫∣k∣<kFd3kN=g_s\sum_{|\mathbf k|<k_F}1\to g_s\dfrac{V}{(2\pi)^3}\displaystyle\int_{|\mathbf k|<k_F}\mathrm d^3k. Why: each allowed k\mathbf k occupies a cell of volume (2π/L)3(2\pi/L)^3 in a box of side LL, so the number of states in a kk-space region is its volume times V/(2π)3V/(2\pi)^3.

Step 2. Do the sphere integral. ∫∣k∣<kFd3k=43πkF3\displaystyle\int_{|\mathbf k|<k_F}\mathrm d^3k=\tfrac43\pi k_F^3, so N=gsV(2π)343πkF3=gsVkF36π2N=g_s\dfrac{V}{(2\pi)^3}\tfrac43\pi k_F^3=\dfrac{g_sV k_F^3}{6\pi^2}. Inverting, kF=(6π2n/gs)1/3k_F=(6\pi^2 n/g_s)^{1/3}, which for gs=2g_s=2 is (3π2n)1/3(3\pi^2 n)^{1/3}. Why: the 43πkF3\tfrac43\pi k_F^3 is just the volume of the filled ball; the factor (2π)3=8π3(2\pi)^3=8\pi^3 combines with 43π\tfrac43\pi to give the clean 6π26\pi^2.

Step 3. Write the energy as the same integral weighted by ε=ℏ2k2/2m\varepsilon=\hbar^2k^2/2m. E0=gsV(2π)3∫0kFℏ2k22m 4πk2 dk=gsV(2π)34πℏ22mkF55=gsVℏ2kF520π2mE_0=g_s\dfrac{V}{(2\pi)^3}\displaystyle\int_0^{k_F}\dfrac{\hbar^2k^2}{2m}\,4\pi k^2\,\mathrm dk=g_s\dfrac{V}{(2\pi)^3}\dfrac{4\pi\hbar^2}{2m}\dfrac{k_F^5}{5}=\dfrac{g_sV\hbar^2 k_F^5}{20\pi^2 m}. Why: the shell at radius kk has volume 4πk2 dk4\pi k^2\,\mathrm dk and each state there carries energy ℏ2k2/2m\hbar^2k^2/2m; the integral ∫0kFk4 dk=kF5/5\int_0^{k_F}k^4\,\mathrm dk=k_F^5/5 supplies the 15\tfrac15.

Step 4. Take the ratio to eliminate every prefactor. E0N=gsVℏ2kF5/20π2mgsVkF3/6π2=620ℏ2kF2m=310ℏ2kF2m=35ℏ2kF22m=35εF\dfrac{E_0}{N}=\dfrac{g_sV\hbar^2 k_F^5/20\pi^2 m}{g_sV k_F^3/6\pi^2}=\dfrac{6}{20}\dfrac{\hbar^2 k_F^2}{m}=\dfrac{3}{10}\dfrac{\hbar^2k_F^2}{m}=\dfrac{3}{5}\dfrac{\hbar^2k_F^2}{2m}=\dfrac35\varepsilon_F. Why: gsg_s, VV, mm, and the numerical prefactors all cancel in the ratio, leaving the pure geometric number 35\tfrac35 times εF=ℏ2kF2/2m\varepsilon_F=\hbar^2k_F^2/2m. Equivalently, the same 35\tfrac35 is the moment ratio ∫0εFε3/2dε / (εF∫0εFε1/2dε)=2/52/3=35\int_0^{\varepsilon_F}\varepsilon^{3/2}\mathrm d\varepsilon\,/\,(\varepsilon_F\int_0^{\varepsilon_F}\varepsilon^{1/2}\mathrm d\varepsilon)=\tfrac{2/5}{2/3}=\tfrac35 of the ε\sqrt\varepsilon density of states. Verified: the sphere count N=gsVkF3/6π2N=g_sVk_F^3/6\pi^2, the inversion kF=(3π2n)1/3k_F=(3\pi^2 n)^{1/3} at gs=2g_s=2, the energy integral E0=gsVℏ2kF5/20π2mE_0=g_sV\hbar^2k_F^5/20\pi^2 m, and both the kk-space and DOS moment ratios giving E0/N=35εFE_0/N=\tfrac35\varepsilon_F are checked symbolically in verify_degenerate-fermi-gas-zero-temperature.py.

Go deeper · Derivation B · 1 minDegeneracy pressure at absolute zero
Derivation BDegeneracy pressure at absolute zero

The goal is to show that the cold gas exerts a finite pressure P0=25nεF=23E0/VP_0=\tfrac25 n\varepsilon_F=\tfrac23 E_0/V, a purely quantum "degeneracy pressure" with no thermal component. The route is mechanical: since S=0S=0 at T=0T=0, the pressure is the adiabatic derivative P0=−(∂E0/∂V)NP_0=-(\partial E_0/\partial V)_N, and the only volume dependence sits in εF∝V−2/3\varepsilon_F\propto V^{-2/3}. Watch the interpretation: this pressure does not vanish as T→0T\to0 — it is the mechanical face of the Pauli principle, and its steep density scaling is what makes it astrophysically decisive. The full stellar treatment is forward-teased to aq-sm-degeneracy-pressure-fermions; here we establish only the equation of state.

Step 1. Identify the pressure with the adiabatic energy derivative. At T=0T=0 the entropy is fixed at S=0S=0, so P0=−(∂E0∂V)N,SP_0=-\left(\dfrac{\partial E_0}{\partial V}\right)_{N,S}. Why: the first law dE=T dS−P dV\mathrm dE=T\,\mathrm dS-P\,\mathrm dV with dS=0\mathrm dS=0 gives P=−(∂E/∂V)SP=-(\partial E/\partial V)_S directly.

Step 2. Extract the volume dependence. With n=N/Vn=N/V, εF=ℏ22m(3π2n)2/3∝n2/3∝V−2/3\varepsilon_F=\dfrac{\hbar^2}{2m}(3\pi^2 n)^{2/3}\propto n^{2/3}\propto V^{-2/3}, so E0=35NεF=C V−2/3E_0=\tfrac35 N\varepsilon_F=C\,V^{-2/3} at fixed NN, where CC collects the VV-independent constants. Why: only εF\varepsilon_F depends on VV at fixed NN; the number of particles and the 35\tfrac35 do not.

Step 3. Differentiate. P0=−∂∂V ⁣(CV−2/3)=23 CV−5/3=23E0VP_0=-\dfrac{\partial}{\partial V}\!\left(CV^{-2/3}\right)=\tfrac23\,CV^{-5/3}=\dfrac23\dfrac{E_0}{V}. Substituting E0/V=35nεFE_0/V=\tfrac35 n\varepsilon_F gives P0=25nεFP_0=\tfrac25 n\varepsilon_F, and with εF∝n2/3\varepsilon_F\propto n^{2/3} this is P0=(3π2)2/35ℏ2mn5/3P_0=\dfrac{(3\pi^2)^{2/3}}{5}\dfrac{\hbar^2}{m}n^{5/3}. Why: the exponent −23-\tfrac23 produces the universal PV=23EPV=\tfrac23 E of a non-relativistic gas; the steep n5/3n^{5/3} growth is what resists gravitational compression. Verified: the mechanical relation P=−(∂E/∂V)N=23E/VP=-(\partial E/\partial V)_N=\tfrac23 E/V for an ε∝k2\varepsilon\propto k^2 gas, the identity P0=25nεFP_0=\tfrac25 n\varepsilon_F, and the n5/3n^{5/3} form (with the ultrarelativistic 34εF\tfrac34\varepsilon_F, n4/3n^{4/3} cross-check) are verified in verify_degenerate-fermi-gas-zero-temperature.py.

Limiting cases

The T=0T=0 ground state is itself a limit, and it is the reference point for several regimes.

CaseResultMeaning
Strong degeneracy T≪TFT\ll T_Fnˉ→Θ(εF−ε)\bar n\to\Theta(\varepsilon_F-\varepsilon)Occupation is a sharp step; finite-TT corrections are analytic in (T/TF)2(T/T_F)^2 via the Sommerfeld expansion. Metals at 300 K live here.
Classical T≫TFT\gg T_Fnˉ→e(μ−ε)/kBT≪1\bar n\to e^{(\mu-\varepsilon)/k_BT}\ll1Occupations are tiny; Fermi-Dirac reduces to Maxwell-Boltzmann and the sphere picture is irrelevant (nλT3≪1n\lambda_T^3\ll1).
Ultra-relativistic pFc≫mc2p_Fc\gg mc^2E0=34NεF, P0=14nεFE_0=\tfrac34 N\varepsilon_F,\ P_0=\tfrac14 n\varepsilon_FLinear dispersion ε=ℏkc\varepsilon=\hbar kc softens the equation of state to P∝n4/3P\propto n^{4/3} — the origin of the Chandrasekhar mass.
Full spin polarization gs=1g_s=1kF=(6π2n)1/3k_F=(6\pi^2 n)^{1/3}One spin species holds all NN; kFk_F and εF\varepsilon_F rise by 21/32^{1/3} and 22/3≈1.592^{2/3}\approx1.59 — polarizing costs kinetic energy.

Worked examples

Example 1 — Fill four fermions into oscillator-like levels and read off the ground state.

Step 1. Setup. Take single-particle levels with energies εk=k\varepsilon_k=k (in units of some ε0\varepsilon_0), each spin-degenerate (gs=2g_s=2), and place N=4N=4 spin-12\tfrac12 fermions. Fill from the bottom, at most 22 per level.

Step 2. Fill. Level ε=0\varepsilon=0 takes 22 (spin up, spin down); level ε=1\varepsilon=1 takes the remaining 22. Levels ε≥2\varepsilon\ge2 are empty. The highest occupied energy is εF=1\varepsilon_F=\mathbf{1} (in units of ε0\varepsilon_0).

Step 3. Ground-state energy by direct sum. E0=2(0)+2(1)=2 ε0E_0=2(0)+2(1)=\mathbf{2\,\varepsilon_0}, so E0/N=2/4=0.5 ε0E_0/N=2/4=\mathbf{0.5\,\varepsilon_0}. This discrete toy has E0/N=12εFE_0/N=\tfrac12\varepsilon_F, not 35\tfrac35 — the 35\tfrac35 is the continuum result for a ε\sqrt\varepsilon density of states, which this tiny system does not yet approximate.

Step 4. Interpret. The exclusion principle alone forces E0>0E_0>0: four distinguishable particles could all sit at ε=0\varepsilon=0 for E0=0E_0=0, but fermions cannot. This is the microcosm of the Fermi sea — filling forced upward by antisymmetry, giving a finite zero-point energy even with no interactions and no temperature.

Example 2 — Fermi parameters of copper's conduction electrons.

Step 1. Setup. Copper contributes one 4s4s electron per atom, giving n=8.47×1028 m−3n=8.47\times10^{28}\,\mathrm{m^{-3}}; with gs=2g_s=2, use kF=(3π2n)1/3k_F=(3\pi^2 n)^{1/3}.

Step 2. Fermi wavevector. 3π2n=3π2(8.47×1028)=2.508×1030 m−33\pi^2 n=3\pi^2(8.47\times10^{28})=2.508\times10^{30}\,\mathrm{m^{-3}}, so kF=(2.508×1030)1/3=1.36×1010 m−1k_F=(2.508\times10^{30})^{1/3}=\mathbf{1.36\times10^{10}\,m^{-1}}.

Step 3. Fermi energy and temperature. εF=ℏ2kF22me=(1.0546×10−34)2(1.36×1010)22(9.109×10−31)=1.127×10−18 J=7.03 eV\varepsilon_F=\dfrac{\hbar^2k_F^2}{2m_e}=\dfrac{(1.0546\times10^{-34})^2(1.36\times10^{10})^2}{2(9.109\times10^{-31})}=1.127\times10^{-18}\,\mathrm J=\mathbf{7.03\ eV}, and TF=εF/kB=8.16×104 KT_F=\varepsilon_F/k_B=\mathbf{8.16\times10^{4}\,K}.

Step 4. The degeneracy verdict. Room temperature is 300 K300\,\mathrm K, so T/TF=300/(8.16×104)=3.7×10−3T/T_F=300/(8.16\times10^4)=\mathbf{3.7\times10^{-3}} — the electron gas is deeply degenerate, sitting essentially in its ground state. Each electron carries E0/N=35εF=4.2 eVE_0/N=\tfrac35\varepsilon_F=\mathbf{4.2\ eV} of zero-point kinetic energy, and the Fermi velocity is vF=ℏkF/me=1.57×106 m/sv_F=\hbar k_F/m_e=\mathbf{1.57\times10^{6}\,m/s} — about 0.5%0.5\% of the speed of light, all of it quantum zero-point motion, none of it thermal.

Origins · 1 minsubrahmanyan chandrasekhar and how this was discovered
SC
subrahmanyan chandrasekhar
1910–1995 · The Indian-American astrophysicist who, on the voyage to England in 1930, applied the cold degenerate Fermi gas to white dwarfs. Combining the relativistic equation of state P₀∝n^{4/3} with hydrostatic equilibrium, he found a maximum mass — the Chandrasekhar limit ≈1.4 solar masses — above which electron degeneracy pressure cannot halt collapse. Ridiculed by Eddington at the time, the result won him the 1983 Nobel Prize and underlies the theory of neutron stars and Type Ia supernovae.

History
1926
— Fermi and, independently, Dirac derive the statistics of exclusion-obeying particles
1927
— Sommerfeld applies Fermi-Dirac statistics to metallic electrons, explaining the small electronic heat capacity
1928
— Bloch extends the free-electron gas to periodic lattices, founding band theory
1930
— Chandrasekhar computes the white-dwarf mass limit on his voyage to Cambridge
1932
— Landau independently proposes a maximum mass for degenerate stars
1967
— Pulsars are discovered, confirming neutron stars supported by neutron degeneracy pressure
1999
— DeMarco and Jin produce the first degenerate atomic Fermi gas of potassium-40

Going further

  • Texts: Pathria & Beale, Statistical Mechanics (3rd ed.), §8.1; Kardar, Statistical Physics of Particles, Ch. 7; Reif, Fundamentals of Statistical and Thermal Physics, §9.16–9.17; Landau & Lifshitz, Statistical Physics Part 1, §57; Ashcroft & Mermin, Solid State Physics, Ch. 2.
  • Builds on: the T=0T=0 step from Fermi-Dirac Distribution; the master sum N=∫g(ε)nˉ dεN=\int g(\varepsilon)\bar n\,\mathrm d\varepsilon from Quantum Statistics of Ideal Gases: Occupation Numbers; the ε\sqrt\varepsilon density of states from Density of States for Free Particles; the energy scale from Fermi Energy, Momentum and Fermi Temperature.
  • Leads to: the Sommerfeld Expansion (finite-TT corrections about this step) · Electronic Heat Capacity of Metals · Degeneracy Pressure of Fermions (the stellar equation of state) · Chemical Potential of a Fermi Gas vs Temperature.
  • Provenance: AllQFT `sm` export (Pathria & Beale, §8.1), aq-sm-degenerate-fermi-gas-zero-temperature, ref STAT.IQGAS.DFG0.1218. All counting, energy-moment, pressure, and Fermi-parameter claims python-verified (22 checks, all PASS). No base error found — the migrated formulae were algebra- and sign-consistent; the copper numbers εF=7.03\varepsilon_F=7.03 eV, TF=8.16×104T_F=8.16\times10^4 K, vF=1.57×106v_F=1.57\times10^6 m/s reproduce the Ashcroft-Mermin values.
Why isn't the ground-state energy just zero at T=0?

For distinguishable classical particles it would be — all NN could pile into the lowest state. But fermions cannot: the Pauli principle allows at most gsg_s per single-particle state, so filling NN of them forces particles into successively higher momentum states. The result is a huge zero-point kinetic energy E0=35NεFE_0=\tfrac35 N\varepsilon_F that exists with no interactions and no temperature. It is a direct, quantitative consequence of antisymmetry — exclusion converted into energy.

Is the degeneracy pressure a kind of thermal pressure?

No — it is the opposite. Thermal pressure P=nkBTP=nk_BT vanishes at T=0T=0; the degeneracy pressure P0=25nεFP_0=\tfrac25 n\varepsilon_F is maximal there and is entirely quantum. It comes from the volume dependence of the zero-point energy: squeezing the box raises εF∝V−2/3\varepsilon_F\propto V^{-2/3}, so the filled sphere pushes back. This is why a white dwarf, though "cold" in the thermal sense, resists gravity: its support is exclusion, not heat.

How can occupation be a sharp step when real metals are at 300 K?

Because 300 K300\,\mathrm K is tiny compared to the Fermi temperature TF∼8×104 KT_F\sim8\times10^4\,\mathrm K. The step softens only over an energy width ∼kBT\sim k_BT, which for copper is ∼0.5%\sim0.5\% of εF\varepsilon_F. To the eye the occupation is still a step; only a thin shell near the surface is "active." That is exactly why the ground-state treatment works so well for metals and why finite-TT effects are captured by a systematic expansion in (T/TF)2(T/T_F)^2.

Where does the three-fifths come from, intuitively?

The density of states rises as ε\sqrt\varepsilon, so there are far more states near the top of the sea than near the bottom. The mean occupied energy is therefore pulled well above the midpoint 12εF\tfrac12\varepsilon_F. Precisely, ⟨ε⟩=∫0εFε⋅ε1/2 dε/∫0εFε1/2 dε=εF5/2/(5/2)εF3/2/(3/2)⋅εF−1⋅εF=35εF\langle\varepsilon\rangle=\int_0^{\varepsilon_F}\varepsilon\cdot\varepsilon^{1/2}\,\mathrm d\varepsilon/\int_0^{\varepsilon_F}\varepsilon^{1/2}\,\mathrm d\varepsilon=\tfrac{\varepsilon_F^{5/2}/(5/2)}{\varepsilon_F^{3/2}/(3/2)}\cdot\varepsilon_F^{-1}\cdot\varepsilon_F=\tfrac35\varepsilon_F. The crowded top does the work.

What changes when the fermions become relativistic?

In dense stellar matter the Fermi momentum can be large enough that pFc≫mc2p_Fc\gg mc^2, and the dispersion switches from ε=ℏ2k2/2m\varepsilon=\hbar^2k^2/2m to ε=ℏkc\varepsilon=\hbar kc. The energy moment then gives E0=34NεFE_0=\tfrac34 N\varepsilon_F and the pressure P0=14nεF∝n4/3P_0=\tfrac14 n\varepsilon_F\propto n^{4/3} — a softer equation of state. That softening is decisive: a softer P(n)P(n) cannot always beat gravity, which is why there is a maximum mass, the Chandrasekhar limit ≈1.4 M⊙\approx1.4\,M_\odot.

Does this ground state assume the fermions don't interact?

Yes — it is the ideal-gas ground state, with independent modes and a filled Fermi sphere. Remarkably, Landau's Fermi-liquid theory shows the picture survives interactions: the low-energy excitations of a real metal are long-lived quasiparticles that still fill a Fermi sea up to a (renormalized) Fermi surface, which is why the free-electron model predicts the right εF\varepsilon_F, heat capacity, and Fermi velocity for copper. Strong correlations (Mott insulators, heavy fermions) can break this, but the ideal ground state remains the reference.

Why is the Fermi velocity so large for a cold gas?

Because vF=ℏkF/mv_F=\hbar k_F/m and kFk_F is fixed by density, not temperature. For copper vF≈1.57×106 m/sv_F\approx1.57\times10^6\,\mathrm{m/s} — about 0.5%0.5\% of cc — even at T=0T=0. This is entirely zero-point motion forced by exclusion: the electrons at the Fermi surface are moving fast not because they are hot but because antisymmetry has pushed them to high momentum. It sets the scale for electrical conduction, since only these fast surface electrons carry current.

Before this

Where this leads