Skip to content

The Fermi momentum pF=ℏkFp_{\mathrm F}=\hbar k_{\mathrm F} is the momentum at the boundary of the filled zero-temperature states of a uniform ideal Fermi gas. The Fermi energy is the one-particle energy at that boundary:

EF≡ϵ(kF).E_{\mathrm F} \equiv \epsilon(k_{\mathrm F}).

For a nonrelativistic particle of mass mm,

EF=ℏ2kF22m=pF22m.E_{\mathrm F} = \frac{\hbar^2k_{\mathrm F}^2}{2m} = \frac{p_{\mathrm F}^2}{2m}.

These quantities are not independent thermodynamic parameters. Once the spatial dimension, number density, internal degeneracy, and dispersion are specified, state counting fixes kFk_{\mathrm F} and hence every other Fermi scale.

The most common source of wrong factors is a silent change between:

  • total density, summed over all internal components;
  • density per spin or hyperfine component;
  • spinless or fully polarized conventions;
  • continuum momentum and crystal momentum;
  • kinetic energy and total relativistic energy.

The formulas below keep those conventions explicit. The ideal Fermi gas develops the thermodynamic model; degenerate Fermi gases explains the low-temperature physics.

Helpful background. Pauli exclusion explains the filling rule; plane-wave normalization explains the reciprocal-space state density.

Unless stated otherwise, assume:

  • a uniform continuum system in dd spatial dimensions;
  • a large periodic box of dd-dimensional measure Vd=Ld\mathcal V_d=L^d;
  • NN particles and total density n=N/Vdn=N/\mathcal V_d;
  • gg degenerate internal modes at every wave vector;
  • an isotropic, monotonically increasing dispersion ϵ(k)\epsilon(k);
  • zero temperature when defining the occupied Fermi region;
  • one fermion at most in each complete one-particle mode.

Internal modes may be spin, hyperfine, or flavor states sharing the same dispersion and Fermi boundary. Equivalent valleys can also supply a multiplicity when each identical pocket is counted once using momentum relative to its own minimum; they are generally centered at different crystal momenta. Unequal populations or inequivalent dispersions require separate counts.

This page writes the Fermi energy as EFE_{\mathrm F}. Many texts and neighboring pages write ϵF\epsilon_{\mathrm F} for the same quantity. The notation does not imply a different physical definition:

EF=ϵF=ϵ(kF).E_{\mathrm F} = \epsilon_{\mathrm F} = \epsilon(k_{\mathrm F}).

For the nonrelativistic formulas, the energy zero is the bottom of the band or free-particle spectrum. The physical degeneracy scale is the energy above that bottom. If an offset ϵ0\epsilon_0 is retained in ϵ(k)\epsilon(k), use [ϵ(kF)−ϵ0]/kB[\epsilon(k_{\mathrm F})-\epsilon_0]/k_{\mathrm B} for the Fermi temperature and the corresponding energy difference for a microscopic Fermi time. At finite temperature, this scale remains defined from zero-temperature filling. It is generally not equal to the temperature-dependent chemical potential μ(T)\mu(T).

For a quadratic dispersion,

ϵ(k)=ℏ2k22m,\epsilon(k) = \frac{\hbar^2k^2}{2m},

the standard formulas are:

Dimension Total density Fermi wave number
1D n=gkF/πn=gk_{\mathrm F}/\pi kF=πn/gk_{\mathrm F}=\pi n/g
2D n=gkF2/(4π)n=gk_{\mathrm F}^2/(4\pi) kF=4πn/gk_{\mathrm F}=\sqrt{4\pi n/g}
3D n=gkF3/(6π2)n=gk_{\mathrm F}^3/(6\pi^2) kF=(6π2n/g)1/3k_{\mathrm F}=(6\pi^2n/g)^{1/3}

The corresponding Fermi energies are:

Dimension Fermi energy
1D EF=ℏ2π2n2/(2mg2)E_{\mathrm F}=\hbar^2\pi^2n^2/(2mg^2)
2D EF=2πℏ2n/(mg)E_{\mathrm F}=2\pi\hbar^2n/(mg)
3D EF=ℏ22m(6π2n/g)2/3E_{\mathrm F}=\dfrac{\hbar^2}{2m}(6\pi^2n/g)^{2/3}

For a balanced spin-1/21/2 gas, g=2g=2 and nn is the total density:

dkF1πn/222πn3(3π2n)1/3\begin{array}{c|c} d & k_{\mathrm F} \\ \hline 1 & \pi n/2 \\ 2 & \sqrt{2\pi n} \\ 3 & (3\pi^2n)^{1/3} \end{array}

For a spinless or fully polarized gas, set g=1g=1.

The 1D Fermi boundary consists of the two points ±kF\pm k_{\mathrm F}; in 2D it is a circle; in 3D it is a sphere. These formulas require monotone isotropic filling from the minimum. The corresponding densities have units of inverse length, inverse area, and inverse volume, respectively. The geometric derivation below explains all three coefficients at once.

Periodic boundary conditions quantize each wave-vector component in steps of 2π/L2\pi/L. One allowed k\mathbf k point per internal mode therefore occupies dd-dimensional wave-vector volume

(2πL)d.\left( \frac{2\pi}{L} \right)^d.

Equivalently, the number of wave-vector states per internal mode in an element ddkd^dk is

Vd(2π)d ddk.\frac{\mathcal V_d}{(2\pi)^d} \,d^dk.

Let Ωd\Omega_d denote the volume of the unit ball in dd Euclidean dimensions:

Ωd=πd/2Γ(d/2+1).\Omega_d = \frac{\pi^{d/2}} {\Gamma(d/2+1)}.

The ball of radius kk has volume Ωdkd\Omega_d k^d. Including gg internal modes, the continuum approximation to the integrated state count is

Nd(k)=gVd(2π)dΩdkd.\mathcal N_d(k) = g \frac{\mathcal V_d}{(2\pi)^d} \Omega_d k^d.

At T=0T=0, the ground state fills the NN lowest one-particle modes. For a monotone isotropic dispersion, this occupied region is the ball ∣k∣≤kF|\mathbf k|\leq k_{\mathrm F}, so

N=gVd(2π)dΩdkFd.N = g \frac{\mathcal V_d}{(2\pi)^d} \Omega_d k_{\mathrm F}^d.

Dividing by Vd\mathcal V_d gives the master relation

n=gΩd(2π)dkFd.n = g \frac{\Omega_d}{(2\pi)^d} k_{\mathrm F}^d.

Solving for the Fermi wave number,

kF=[(2π)dgΩdn]1/d.k_{\mathrm F} = \left[ \frac{(2\pi)^d}{g\Omega_d} n \right]^{1/d}.

The numerical coefficients in 1D, 2D, and 3D are therefore geometric coefficients, not additional dynamical assumptions.

Filled Fermi regions shown as an interval in one dimension, a disk in two dimensions, and a ball in three dimensions

For an isotropic monotone dispersion, the zero-temperature occupied region is the dd-ball Bd(kF)\mathcal B_d(k_{\mathrm F}). Its geometric measure is 2kF2k_{\mathrm F} in 1D, πkF2\pi k_{\mathrm F}^2 in 2D, and 4πkF3/34\pi k_{\mathrm F}^3/3 in 3D. Multiplication by gVd/(2π)dg\mathcal V_d/(2\pi)^d converts that measure into a state count.

The degeneracy factor counts independent one-particle modes at the same k\mathbf k. It is part of the state count, not a correction applied after the count.

If spin is conserved and all 2s+12s+1 spin projections are degenerate and equally populated, then

gs=2s+1.g_s = 2s+1.

Examples include:

  • spinless fermions or a fully polarized sample: g=1g=1;
  • balanced spin-1/21/2 fermions: g=2g=2;
  • an idealized balanced spin-3/23/2 multiplet: g=4g=4.

This formula does not say that every spin multiplet is physically populated. Zeeman splitting, spin–orbit coupling, preparation constraints, interactions, or selection rules may remove the degeneracy.

If spin, valley, and another internal label are all independent and exactly degenerate, their multiplicities multiply:

g=gsgvgother.g = g_s g_v g_{\mathrm{other}}.

One must not multiply by a label that is already counted as a separate band or included explicitly in the sum over states.

Suppose gg components are balanced, each with density nan_a. Then

n=∑a=1gna=gna.n = \sum_{a=1}^{g}n_a = g n_a.

The per-component state count contains no additional degeneracy factor:

na=Ωd(2π)dkFd.n_a = \frac{\Omega_d}{(2\pi)^d} k_{\mathrm F}^d.

Substituting na=n/gn_a=n/g recovers the total-density formula.

If components have unequal densities, there is generally no single common Fermi wave number. Each component has

kF,a=[(2π)dΩdna]1/d.k_{{\mathrm F},a} = \left[ \frac{(2\pi)^d}{\Omega_d} n_a \right]^{1/d}.

For two spin components in three dimensions,

kF,↑=(6π2n↑)1/3,kF,↓=(6π2n↓)1/3.k_{{\mathrm F},\uparrow} = (6\pi^2n_\uparrow)^{1/3}, \qquad k_{{\mathrm F},\downarrow} = (6\pi^2n_\downarrow)^{1/3}.

Writing one g=2g=2 formula for an imbalanced gas hides the two distinct Fermi surfaces.

At fixed total density,

kF∝g−1/d,k_{\mathrm F} \propto g^{-1/d},

and, for a quadratic dispersion,

EF∝g−2/d.E_{\mathrm F} \propto g^{-2/d}.

Reducing the number of available internal components forces particles to occupy a larger region of momentum space. In three dimensions, fully polarizing a previously balanced spin-1/21/2 ideal gas changes the scales by

kF(g=1)kF(g=2)=21/3,\frac{k_{\mathrm F}^{(g=1)}} {k_{\mathrm F}^{(g=2)}} = 2^{1/3},

and

EF(g=1)EF(g=2)=22/3.\frac{E_{\mathrm F}^{(g=1)}} {E_{\mathrm F}^{(g=2)}} = 2^{2/3}.

Once kFk_{\mathrm F} is known, several useful scales follow.

Wave number and momentum differ by ℏ\hbar:

pF=ℏkF.p_{\mathrm F} = \hbar k_{\mathrm F}.

Thus kFk_{\mathrm F} has units of inverse length, while pFp_{\mathrm F} has units of momentum. The phrase Fermi momentum is often used informally for either quantity; a careful formula should display which one is meant.

For a general isotropic dispersion, the group velocity at the Fermi boundary is

vF=1ℏdϵdk∣k=kF.v_{\mathrm F} = \frac{1}{\hbar} \left. \frac{d\epsilon}{dk} \right|_{k=k_{\mathrm F}}.

For a quadratic dispersion,

vF=ℏkFm=2EFm.v_{\mathrm F} = \frac{\hbar k_{\mathrm F}}{m} = \sqrt{ \frac{2E_{\mathrm F}}{m} }.

The Fermi temperature is the energy scale expressed in kelvin:

TF≡EFkB.T_{\mathrm F} \equiv \frac{E_{\mathrm F}}{k_{\mathrm B}}.

It is not the temperature of the zero-temperature gas. It defines the reduced temperature

θ=TTF,\theta = \frac{T}{T_{\mathrm F}},

which distinguishes the degenerate regime θ≪1\theta\ll1 from the classical regime θ≫1\theta\gg1 under the appropriate density conditions.

The de Broglie wavelength associated with pFp_{\mathrm F} is

λF=hpF=2πkF.\lambda_{\mathrm F} = \frac{h}{p_{\mathrm F}} = \frac{2\pi}{k_{\mathrm F}}.

Many-body estimates also use the shorter length kF−1k_{\mathrm F}^{-1}. These differ by 2π2\pi and should not be interchanged silently.

A natural microscopic time is

tF=ℏEF.t_{\mathrm F} = \frac{\hbar}{E_{\mathrm F}}.

For a quadratic dispersion,

tF=2mℏkF2.t_{\mathrm F} = \frac{2m}{\hbar k_{\mathrm F}^2}.

This is a scale, not automatically a collision time or equilibration time.

If the nonrelativistic Fermi energy is known, then

kF=2mEFℏ.k_{\mathrm F} = \frac{\sqrt{2mE_{\mathrm F}}}{\hbar}.

The corresponding density in dd dimensions is

n=gΩd(2π)d(2mEFℏ2)d/2.n = g \frac{\Omega_d}{(2\pi)^d} \left( \frac{2mE_{\mathrm F}}{\hbar^2} \right)^{d/2}.

Explicitly,

n1D=gπ2mEFℏ,n2D=gmEF2πℏ2,n3D=g6π2(2mEFℏ2)3/2.\begin{aligned} n_{1\mathrm D} &= \frac{g}{\pi} \frac{\sqrt{2mE_{\mathrm F}}}{\hbar}, \\ n_{2\mathrm D} &= \frac{gmE_{\mathrm F}} {2\pi\hbar^2}, \\ n_{3\mathrm D} &= \frac{g}{6\pi^2} \left( \frac{2mE_{\mathrm F}}{\hbar^2} \right)^{3/2}. \end{aligned}

The densities have different dimensions: inverse length in 1D, inverse area in 2D, and inverse volume in 3D.

The integrated number of states below energy EE for a quadratic dispersion is

Nd(E)Vd=gΩd(2π)d(2mEℏ2)d/2.\frac{\mathcal N_d(E)}{\mathcal V_d} = g \frac{\Omega_d}{(2\pi)^d} \left( \frac{2mE}{\hbar^2} \right)^{d/2}.

Differentiation gives the density of states per dd-dimensional measure:

ρd(E)=1VddNddE.\rho_d(E) = \frac{1}{\mathcal V_d} \frac{d\mathcal N_d}{dE}.

At the Fermi energy,

ρd(EF)=d2nEF.\rho_d(E_{\mathrm F}) = \frac{d}{2} \frac{n}{E_{\mathrm F}}.

This identity is a fast consistency check. It uses a density of states that includes the same internal factor gg as the density nn.

The energy dependence is

ρd(E)∝Ed/2−1.\rho_d(E) \propto E^{d/2-1}.

Thus the free-particle density of states diverges as E−1/2E^{-1/2} in 1D, is constant in 2D, and grows as E1/2E^{1/2} in 3D. The general density-of-states normalization is explained in Density of States: First Encounter.

For a quadratic dispersion, integrating all occupied energies gives

U0N=dd+2EF.\frac{U_0}{N} = \frac{d}{d+2} E_{\mathrm F}.

Scale invariance of the nonrelativistic kinetic energy gives

P0=2dU0Vd=2d+2nEF.P_0 = \frac{2}{d} \frac{U_0}{\mathcal V_d} = \frac{2}{d+2} nE_{\mathrm F}.

In one dimension, P0P_0 is a force; in two dimensions, it is force per length; in three dimensions, it is the usual force per area. These identities are useful checks, but the full equation-of-state derivation belongs to the Ideal Fermi Gas page.

Worked example: a three-dimensional electron density

Section titled “Worked example: a three-dimensional electron density”

Consider nonrelativistic electrons with total density

n=8.5×1028 m−3n = 8.5\times10^{28}\ \mathrm{m}^{-3}

and spin degeneracy g=2g=2. Then

kF=(3π2n)1/3≈1.36×1010 m−1.k_{\mathrm F} = (3\pi^2n)^{1/3} \approx 1.36\times10^{10}\ \mathrm{m}^{-1}.

The momentum scale is

pF=ℏkF≈1.43×10−24 kg m s−1.p_{\mathrm F} = \hbar k_{\mathrm F} \approx 1.43\times10^{-24}\ \mathrm{kg\,m\,s^{-1}}.

Using the electron mass,

EF=ℏ2kF22me≈7.0 eV.E_{\mathrm F} = \frac{\hbar^2k_{\mathrm F}^2}{2m_e} \approx 7.0\ \mathrm{eV}.

Consequently,

vF≈1.6×106 m s−1,v_{\mathrm F} \approx 1.6\times10^6\ \mathrm{m\,s^{-1}},

and

TF≈8.1×104 K.T_{\mathrm F} \approx 8.1\times10^4\ \mathrm K.

An ordinary laboratory temperature can therefore be much smaller than TFT_{\mathrm F} even when it is hundreds of kelvin. This is why conduction electrons are often strongly degenerate. A real metal may require band-dependent effective masses and a nonspherical Fermi surface, so the numerical calculation is an ideal free-electron estimate rather than a complete material model.

Worked example: a two-dimensional electron gas

Section titled “Worked example: a two-dimensional electron gas”

Take a two-dimensional spin-degenerate gas with sheet density

n2D=1.0×1015 m−2,n_{2\mathrm D} = 1.0\times10^{15}\ \mathrm{m}^{-2},

g=2g=2, and effective mass m∗=0.067mem^*=0.067m_e. The Fermi wave number is

kF=2πn2D≈7.93×107 m−1.k_{\mathrm F} = \sqrt{2\pi n_{2\mathrm D}} \approx 7.93\times10^7\ \mathrm{m}^{-1}.

The quadratic-band Fermi energy is

EF=πℏ2n2Dm∗≈3.6 meV.E_{\mathrm F} = \frac{\pi\hbar^2n_{2\mathrm D}}{m^*} \approx 3.6\ \mathrm{meV}.

Equivalently,

TF≈41 K.T_{\mathrm F} \approx 41\ \mathrm K.

The use of m∗m^* rather than the bare electron mass is part of the model specification. Valley degeneracy, spin splitting, or nonparabolicity would modify the count or the dispersion.

The state-counting relation for kFk_{\mathrm F} depends on the occupied kk-space volume, not on the detailed energy dispersion. For a continuum dispersion, or an isotropic pocket counted within its domain of validity, the same density relation holds:

n=gΩd(2π)dkFd.n = g \frac{\Omega_d}{(2\pi)^d} k_{\mathrm F}^d.

Only the map from kFk_{\mathrm F} to energy and velocity changes:

EF=ϵ(kF),E_{\mathrm F} = \epsilon(k_{\mathrm F}),

and

vF=1ℏϵ′(kF).v_{\mathrm F} = \frac{1}{\hbar} \epsilon'(k_{\mathrm F}).

For a power-law dispersion

ϵ(k)=Aks,\epsilon(k) = A k^s,

one obtains

EF∝ns/d,E_{\mathrm F} \propto n^{s/d},

and

ρd(EF)=dsnEF.\rho_d(E_{\mathrm F}) = \frac{d}{s} \frac{n}{E_{\mathrm F}}.

The familiar factor d/2d/2 is the special case s=2s=2.

Relativistic state counting gives the same kFk_{\mathrm F} for a specified density and degeneracy. What changes is the energy relation. With

pF=ℏkF,p_{\mathrm F} = \hbar k_{\mathrm F},

the total one-particle energy at the Fermi boundary is

EF=pF2c2+m2c4.\mathcal E_{\mathrm F} = \sqrt{ p_{\mathrm F}^2c^2 + m^2c^4 }.

If Fermi energy means kinetic energy measured from the rest energy, then

EFkin=pF2c2+m2c4−mc2.E_{\mathrm F}^{\mathrm{kin}} = \sqrt{ p_{\mathrm F}^2c^2 + m^2c^4 } - mc^2.

At zero temperature, a relativistic chemical potential that includes rest energy is

μ0=EF,\mu_0 = \mathcal E_{\mathrm F},

whereas a convention that subtracts rest energy gives

μ0kin=EFkin.\mu_0^{\mathrm{kin}} = E_{\mathrm F}^{\mathrm{kin}}.

Both conventions are legitimate; mixing them in one calculation is not. The nonrelativistic approximation requires

pF≪mc.p_{\mathrm F} \ll mc.

The scalar kFk_{\mathrm F} is exact for a spherical Fermi boundary. In a crystal band, the occupied region is determined by

En(k)≤μ0,E_n(\mathbf k) \leq \mu_0,

inside the Brillouin zone. The boundary can be warped, disconnected, open across a zone boundary, or split among bands. Then one should speak of:

  • the Fermi surface En(k)=μ0E_n(\mathbf k)=\mu_0;
  • direction-dependent Fermi wave vectors;
  • extremal radii or orbit areas;
  • an effective spherical kFk_{\mathrm F} only when explicitly defined.

For an anisotropic quadratic band,

ϵ(k)=ℏ22(kx2mx+ky2my+kz2mz),\epsilon(\mathbf k) = \frac{\hbar^2}{2} \left( \frac{k_x^2}{m_x} + \frac{k_y^2}{m_y} + \frac{k_z^2}{m_z} \right),

constant-energy surfaces are ellipsoids with semiaxes kF,i=2miEF/ℏk_{{\mathrm F},i}=\sqrt{2m_iE_{\mathrm F}}/\hbar. Their volume gives

n=g6π2(2EFℏ2)3/2mxmymz.n=\frac{g}{6\pi^2}\left(\frac{2E_{\mathrm F}}{\hbar^2}\right)^{3/2} \sqrt{m_xm_ym_z}.

The density-of-states mass mDOS=(mxmymz)1/3m_{\mathrm{DOS}}=(m_xm_ym_z)^{1/3} reproduces this total count in the spherical formula. It does not reproduce the directional Fermi radii or velocities; those still need the individual masses.

Fermi Surface develops the general many-body concept and low-energy excitations near nonspherical boundaries. The linked band-structure discussions explain Bloch-band topology and material-specific surfaces.

A finite trap does not have translational invariance, so momentum states are not labeled by a uniform continuum density. Exact state counting uses trap eigenlevels.

In a slowly varying potential, a local-density approximation may define

n(r)=gΩd(2π)dkF(r)d,n(\mathbf r) = g \frac{\Omega_d}{(2\pi)^d} k_{\mathrm F}(\mathbf r)^d,

with a local Fermi energy

EF(r)=ℏ2kF(r)22m.E_{\mathrm F}(\mathbf r) = \frac{\hbar^2k_{\mathrm F}(\mathbf r)^2}{2m}.

These are local scales. They need not equal the global Fermi energy defined by the highest occupied trap level. Quantum gases in traps derives the trap conventions and local-density approximation.

In a finite box, allowed momenta are discrete. The continuum formula replaces a lattice-point count by a geometric volume, so it becomes accurate when many modes are occupied and boundary corrections are small.

For a finite noninteracting system, several quantities may be called the Fermi energy:

  • the highest occupied one-particle energy;
  • the lowest unoccupied one-particle energy;
  • the midpoint between them;
  • an addition or removal chemical potential;
  • the continuum estimate from the mean density.

For the gapless continuum gas, these differences vanish in the ordinary thermodynamic limit, but can be a shell spacing in a small system. A band insulator with a persistent gap does not have the same coincidence. A partially filled degenerate shell also makes the zero-temperature boundary less sharp than the continuum picture suggests.

Translate between total and component density

Section titled “Translate between total and component density”

A balanced two-component gas in two dimensions has total sheet density nn. Compute kFk_{\mathrm F} first from the total-density formula with g=2g=2, then from the density nσ=n/2n_\sigma=n/2 of one component. Verify agreement.

Solution

Using total density and g=2g=2,

kF=4πn2=2πn.k_{\mathrm F} = \sqrt{ \frac{4\pi n}{2} } = \sqrt{2\pi n}.

For one component, no extra degeneracy factor is included:

nσ=kF24π.n_\sigma = \frac{k_{\mathrm F}^2}{4\pi}.

Since nσ=n/2n_\sigma=n/2,

kF=4πnσ=2πn.k_{\mathrm F} = \sqrt{4\pi n_\sigma} = \sqrt{2\pi n}.

Both expressions count each spin component once and therefore agree.

Starting from Nd(E)/Vd∝Ed/2\mathcal N_d(E)/\mathcal V_d\propto E^{d/2} for a quadratic dispersion, show that

ρd(EF)=d2nEF.\rho_d(E_{\mathrm F}) = \frac{d}{2} \frac{n}{E_{\mathrm F}}.
Solution

Write

Nd(E)Vd=CEd/2,\frac{\mathcal N_d(E)}{\mathcal V_d} = C E^{d/2},

where CC contains the mass, degeneracy, and geometric factors. Differentiation gives

ρd(E)=d2CEd/2−1.\rho_d(E) = \frac{d}{2} C E^{d/2-1}.

At E=EFE=E_{\mathrm F},

n=CEFd/2.n = C E_{\mathrm F}^{d/2}.

Eliminating CC yields

ρd(EF)=d2nEF.\rho_d(E_{\mathrm F}) = \frac{d}{2} \frac{n}{E_{\mathrm F}}.

Use ρd(E)∝Ed/2−1\rho_d(E)\propto E^{d/2-1} to prove

U0N=dd+2EF.\frac{U_0}{N} = \frac{d}{d+2} E_{\mathrm F}.

Check the results in 1D, 2D, and 3D.

Solution

The particle number and energy are proportional to

N∝∫0EFEd/2−1 dE=2dEFd/2,N \propto \int_0^{E_{\mathrm F}} E^{d/2-1}\,dE = \frac{2}{d} E_{\mathrm F}^{d/2},

and

U0∝∫0EFEd/2 dE=2d+2EFd/2+1.U_0 \propto \int_0^{E_{\mathrm F}} E^{d/2}\,dE = \frac{2}{d+2} E_{\mathrm F}^{d/2+1}.

Taking the ratio gives

U0N=dd+2EF.\frac{U_0}{N} = \frac{d}{d+2} E_{\mathrm F}.

Thus

U0N={EF/3,d=1,EF/2,d=2,3EF/5,d=3.\frac{U_0}{N} = \begin{cases} E_{\mathrm F}/3, & d=1,\\ E_{\mathrm F}/2, & d=2,\\ 3E_{\mathrm F}/5, & d=3. \end{cases}

Show that the relativistic kinetic Fermi energy approaches pF2/(2m)p_{\mathrm F}^2/(2m) when pF≪mcp_{\mathrm F}\ll mc, and approaches pFcp_{\mathrm F}c when pF≫mcp_{\mathrm F}\gg mc.

Solution

The kinetic Fermi energy is

EFkin=mc2[1+(pFmc)2−1].E_{\mathrm F}^{\mathrm{kin}} = mc^2 \left[ \sqrt{ 1+ \left( \frac{p_{\mathrm F}}{mc} \right)^2 } -1 \right].

For x=pF/(mc)≪1x=p_{\mathrm F}/(mc)\ll1,

1+x2=1+x22+O(x4),\sqrt{1+x^2} = 1+ \frac{x^2}{2} + O(x^4),

so

EFkin=pF22m+O ⁣(pF4m3c2).E_{\mathrm F}^{\mathrm{kin}} = \frac{p_{\mathrm F}^2}{2m} + O\!\left( \frac{p_{\mathrm F}^4}{m^3c^2} \right).

For x≫1x\gg1,

1+x2=x+O(x−1),\sqrt{1+x^2} = x+ O(x^{-1}),

and therefore

EFkin=pFc−mc2+O ⁣(m2c3pF).E_{\mathrm F}^{\mathrm{kin}} = p_{\mathrm F}c - mc^2 + O\!\left( \frac{m^2c^3}{p_{\mathrm F}} \right).

The leading density-dependent term is pFcp_{\mathrm F}c. The subtractive rest-energy constant remains because the kinetic-energy convention measures from mc2mc^2.

At fixed total density and degeneracy, compare the quadratic band with masses (m,4m,9m)(m,4m,9m) to the isotropic band with mass mm. Find the energy ratio, the three radius ratios, and the velocity ratios.

Solution

The density-of-states mass is mDOS=361/3mm_{\mathrm{DOS}}=36^{1/3}m, so EFell/EFsph=36−1/3E_{\mathrm F}^{\mathrm{ell}}/E_{\mathrm F}^{\mathrm{sph}}=36^{-1/3}. Since each radius is proportional to miEF\sqrt{m_iE_{\mathrm F}}, the ratios are 36−1/6(1,2,3)36^{-1/6}(1,2,3). Their product is one, as fixed occupied volume requires. The velocities are ℏki/mi\hbar k_i/m_i, so their ratios are 36−1/6(1,1/2,1/3)36^{-1/6}(1,1/2,1/3). One scalar mass preserves the state count but not these directional properties.

  • N. W. Ashcroft and N. D. Mermin, Solid State Physics, Holt, Rinehart and Winston (1976) — free-electron and band-theory Fermi scales.
  • A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, Dover (2003) — Fermi momentum, density of states, and many-body conventions.
  • S. Giorgini, L. P. Pitaevskii, and S. Stringari, “Theory of ultracold atomic Fermi gases,” Reviews of Modern Physics 80, 1215–1274 (2008) — balanced and imbalanced atomic-gas conventions.
  • G. F. Giuliani and G. Vignale, Quantum Theory of the Electron Liquid, Cambridge University Press (2005) — electron-gas conventions in two and three dimensions.
  • X.-W. Guan, M. T. Batchelor, and C. Lee, “Fermi gases in one dimension: From Bethe ansatz to experiments,” Reviews of Modern Physics 85, 1633–1691 (2013) — one-dimensional Fermi scales and component conventions.
  • K. Huang, Statistical Mechanics, 2nd ed., Wiley (1987) — ideal Fermi-gas thermodynamics and dimensional state counting.
  • M. Kardar, Statistical Physics of Particles, Cambridge University Press (2007) — phase-space counting and quantum statistics.
  • R. K. Pathria and P. D. Beale, Statistical Mechanics, 3rd ed., Elsevier (2011) — continuum state counting and ideal quantum gases.