Skip to content

Low-Dimensional Quantum Gases

A low-dimensional quantum gas is a many-particle gas whose thermally and dynamically accessible motion is restricted to one or two spatial directions. The phrase includes genuinely one- or two-dimensional models and experimentally realized quasi-one-dimensional or quasi-two-dimensional gases in which strong confinement freezes the remaining directions.

Dimensionality is not a cosmetic change of coordinates. For particles with quadratic dispersion,

ϵk=ℏ2k22m,\epsilon_{\mathbf k} = \frac{\hbar^2k^2}{2m},

the density of one-particle states per active volume behaves as

Dd(ϵ)Vd∝ϵd/2−1.\frac{D_d(\epsilon)}{V_d} \propto \epsilon^{d/2-1}.

It diverges at the band edge in one dimension, is constant in two dimensions, and vanishes as ϵ\sqrt{\epsilon} in three dimensions. That single change reorganizes ideal-gas thermodynamics, the infrared capacity of a Bose gas, the geometry of a Fermi sea, and the importance of long-wavelength fluctuations.

Low dimension does not merely reduce the number of momentum components. It increases the relative weight of low-energy modes, where statistics, interactions, and collective fluctuations are hardest to ignore.

This page owns uniform one- and two-dimensional state counting, ideal Bose and Fermi consequences, the distinction between strict and effective dimensionality, and the fluctuation logic that points toward low-dimensional many-body physics. Bose–Einstein Condensation owns the full condensation criterion. Sommerfeld Expansion owns systematic low-temperature Fermi integrals. Lieb–Liniger Model Preview owns the exact finite-coupling contact gas, while Tonks–Girardeau Gas Preview owns the impenetrable Bose–Fermi map. This page previews the Berezinskii–Kosterlitz–Thouless transition; Superfluidity in Condensed Matter applies that logic to neutral films and their material evidence, while the Luttinger-liquid page owns its dedicated theory.

Ultracold Atoms owns the platform-level freeze-out, trap-calibration, and interaction-tuning checks used to establish that a laboratory gas realizes this regime.

Three different statements are often compressed into the same phrase.

  1. A strictly dd-dimensional model has configuration space Rd\mathbb R^d, a dd-dimensional measure, and interactions defined in that dimension.
  2. A kinematically quasi-dd-dimensional gas lives in three-dimensional space but occupies only the transverse ground mode over the energies being probed.
  3. A geometrically anisotropic gas is merely much longer or wider in some directions. It need not be dynamically low dimensional if excited transverse modes are populated.

For example, a long cigar-shaped cloud is not automatically one dimensional. It becomes quasi-one-dimensional only when transverse excitations are energetically inaccessible and interactions do not strongly admix them.

If a transverse direction is harmonically confined with frequency ω⊥\omega_\perp, adjacent transverse levels are separated by

Δ⊥=ℏω⊥.\Delta_\perp = \hbar\omega_\perp.

A practical low-energy requirement is

Eactive≪ℏω⊥,E_{\mathrm{active}} \ll \hbar\omega_\perp,

where EactiveE_{\mathrm{active}} must include every relevant scale: kBTk_{\mathrm B}T, a Fermi energy, interaction energy per particle, collective-mode frequencies, drive bandwidths, and any chemical-potential offset measured from the transverse ground threshold.

The word effective matters. A sample can be one dimensional for equilibrium thermodynamics and three dimensional during a violent quench that excites transverse motion.

Consider a large periodic box with active dd-dimensional volume VdV_d and internal multiplicity gg. One momentum state occupies volume (2π)d/Vd(2\pi)^d/V_d in k\mathbf k-space. The number of states inside a sphere of radius kk is therefore

Nd(k)=gVd(2π)dSd−1dkd,\mathcal N_d(k) = g\frac{V_d}{(2\pi)^d} \frac{S_{d-1}}{d}k^d,

where

Sd−1=2πd/2Γ(d/2)S_{d-1} = \frac{2\pi^{d/2}}{\Gamma(d/2)}

is the area of the unit (d−1)(d-1)-sphere.

For quadratic dispersion,

k=2mϵℏ.k = \frac{\sqrt{2m\epsilon}}{\hbar}.

The integrated density of states becomes

Nd(ϵ)Vd=gΓ(d/2+1)(m2πℏ2)d/2ϵd/2.\frac{\mathcal N_d(\epsilon)}{V_d} = \frac{g}{\Gamma(d/2+1)} \left( \frac{m}{2\pi\hbar^2} \right)^{d/2} \epsilon^{d/2}.

Differentiating with respect to energy gives the density of states

Dd(ϵ)Vd=gΓ(d/2)(m2πℏ2)d/2ϵd/2−1.\frac{D_d(\epsilon)}{V_d} = \frac{g}{\Gamma(d/2)} \left( \frac{m}{2\pi\hbar^2} \right)^{d/2} \epsilon^{d/2-1}.

Here V1=LV_1=L, V2=AV_2=A, and V3=VV_3=V. The units change with dimension: D1/LD_1/L is a density per length per energy, while D2/AD_2/A is a density per area per energy.

The explicit results are

D1(ϵ)L=gπℏm2ϵ,\frac{D_1(\epsilon)}{L} = \frac{g}{\pi\hbar} \sqrt{\frac{m}{2\epsilon}},

and

D2(ϵ)A=gm2πℏ2.\frac{D_2(\epsilon)}{A} = \frac{gm}{2\pi\hbar^2}.

For comparison,

D3(ϵ)V=g4π2(2mℏ2)3/2ϵ.\frac{D_3(\epsilon)}{V} = \frac{g}{4\pi^2} \left( \frac{2m}{\hbar^2} \right)^{3/2} \sqrt{\epsilon}.

The one-dimensional divergence at ϵ=0\epsilon=0 is integrable by itself: the integrated number of states still vanishes as ϵ\sqrt{\epsilon}. The more severe Bose infrared problem appears only after this density of states is multiplied by the divergent low-energy occupation factor.

Confinement geometries, low-dimensional densities of states, and Bose infrared scaling

Strong confinement freezes transverse motion and leaves a tube-like quasi-1D gas or a plane-like quasi-2D gas. For quadratic dispersion, D1(ϵ)∝ϵ−1/2D_1(\epsilon)\propto\epsilon^{-1/2} and D2(ϵ)D_2(\epsilon) is constant. Multiplication by the low-energy Bose factor nB(ϵ)∼(kBT)/ϵn_{\mathrm B}(\epsilon)\sim(k_{\mathrm B}T)/\epsilon makes the excited-state integral infrared divergent for d=1,2d=1,2, whereas the d=3d=3 singularity remains integrable.

Thermal Wavelength and Phase-Space Density

Section titled “Thermal Wavelength and Phase-Space Density”

The thermal de Broglie wavelength is dimension independent as a length,

λT=2πℏ2mkBT.\lambda_T = \sqrt{ \frac{2\pi\hbar^2}{mk_{\mathrm B}T} }.

What changes is the dimensionless degeneracy parameter. If nd=N/Vdn_d=N/V_d, then

Dd=ndλTdg.\mathcal D_d = \frac{n_d\lambda_T^d}{g}.

Thus the relevant combinations are n1λT/gn_1\lambda_T/g in one dimension and n2λT2/gn_2\lambda_T^2/g in two dimensions. Quantum statistics becomes important when Dd\mathcal D_d is not small. This criterion diagnoses degeneracy, not condensation: a uniform two-dimensional Bose gas can be deeply degenerate even though its ideal thermodynamic limit has no finite-temperature Bose–Einstein transition.

For fugacity

z=eβμ,β=1kBT,z = e^{\beta\mu}, \qquad \beta = \frac{1}{k_{\mathrm B}T},

the uniform ideal-gas number equations are

nd(B)=gλTdLi⁡d/2(z),0<z≤1,n_d^{(\mathrm B)} = \frac{g}{\lambda_T^d} \operatorname{Li}_{d/2}(z), \qquad 0<z\le 1,

and

nd(F)=−gλTdLi⁡d/2(−z),z>0.n_d^{(\mathrm F)} = -\frac{g}{\lambda_T^d} \operatorname{Li}_{d/2}(-z), \qquad z>0.

These equations follow by integrating the Bose–Einstein or Fermi–Dirac occupation against Dd(ϵ)D_d(\epsilon). Their order, d/2d/2, is a direct imprint of spatial dimension.

At a putative ideal-gas condensation threshold, the chemical potential approaches the one-particle ground energy. Taking that energy to be zero,

μ→0−.\mu \to 0^-.

For ϵ≪kBT\epsilon\ll k_{\mathrm B}T,

1eβϵ−1≃kBTϵ.\frac{1}{e^{\beta\epsilon}-1} \simeq \frac{k_{\mathrm B}T}{\epsilon}.

The low-energy contribution to the excited density therefore scales as

nex∼∫0dϵ ϵd/2−2.n_{\mathrm{ex}} \sim \int_0 d\epsilon\, \epsilon^{d/2-2}.

The lower limit converges only if

d>2.d > 2.

For a uniform ideal gas with quadratic dispersion:

  • in d=1d=1, the integrand behaves as ϵ−3/2\epsilon^{-3/2};
  • in d=2d=2, it behaves as ϵ−1\epsilon^{-1} and diverges logarithmically;
  • in d=3d=3, it behaves as ϵ−1/2\epsilon^{-1/2} and is integrable.

The failure in one and two dimensions is not a shortage of low-energy states. It is the opposite: the combined low-energy state density and Bose enhancement allow the excited modes to absorb arbitrarily many particles as the system grows. There is no finite excited-state capacity that forces a macroscopic population into one mode.

The Bose–Einstein Condensation page gives the general density-of-states criterion and owns the full transition derivation. Here the result serves as the first signal that low-dimensional thermodynamic limits are infrared sensitive.

Two dimensions provide an especially transparent example because

Li⁡1(z)=−ln⁡(1−z).\operatorname{Li}_1(z) = -\ln(1-z).

The number equation is

n2λT2g=−ln⁡(1−z).\frac{n_2\lambda_T^2}{g} = -\ln(1-z).

At any finite density and nonzero temperature,

z=1−exp⁡(−n2λT2g)<1.z = 1 - \exp\left( -\frac{n_2\lambda_T^2}{g} \right) < 1.

Consequently,

μ=kBTln⁡[1−exp⁡(−n2λT2g)]<0.\mu = k_{\mathrm B}T \ln\left[ 1 - \exp\left( -\frac{n_2\lambda_T^2}{g} \right) \right] < 0.

As the gas becomes highly degenerate, μ\mu approaches zero exponentially closely but does not reach it at finite TT in the uniform thermodynamic limit. A finite experimental resolution can therefore make a strongly degenerate ideal gas look threshold-like even when no exact phase transition exists.

In one dimension,

n1=gλTLi⁡1/2(z).n_1 = \frac{g}{\lambda_T} \operatorname{Li}_{1/2}(z).

Write

z=e−α,α=−βμ>0.z = e^{-\alpha}, \qquad \alpha = -\beta\mu > 0.

As α→0+\alpha\to0^+,

Li⁡1/2(e−α)≃πα+ζ(1/2)+O(α).\operatorname{Li}_{1/2}(e^{-\alpha}) \simeq \sqrt{\frac{\pi}{\alpha}} + \zeta(1/2) + O(\sqrt{\alpha}).

Keeping the leading term gives

α≃πg2n12λT2,\alpha \simeq \frac{\pi g^2}{n_1^2\lambda_T^2},

and hence

μ≃−g2m(kBT)22ℏ2n12.\mu \simeq -\frac{g^2m(k_{\mathrm B}T)^2}{2\hbar^2n_1^2}.

The chemical potential again remains negative at every nonzero temperature. Its algebraic approach to zero is different from the exponentially small two-dimensional result, reflecting the stronger one-dimensional infrared divergence.

This asymptotic formula is for a uniform, ideal continuum gas at fixed line density. It is not an equation of state for an interacting one-dimensional Bose fluid.

Finite Size Is Not a Thermodynamic Transition

Section titled “Finite Size Is Not a Thermodynamic Transition”

A finite box has a lowest nonzero excitation energy. For a periodic length LL,

ϵL∼ℏ22m(2πL)2.\epsilon_L \sim \frac{\hbar^2}{2m} \left( \frac{2\pi}{L} \right)^2.

This energy acts as an infrared cutoff. The excited population is then finite at μ=0\mu=0, and a rapid increase of ground-mode occupation can occur as temperature or particle number is varied.

The cutoff disappears as L→∞L\to\infty. In two dimensions the maximum excited density grows logarithmically with system size,

nex,max(2)∼gλT2ln⁡(kBTϵL),n_{\mathrm{ex,max}}^{(2)} \sim \frac{g}{\lambda_T^2} \ln\left( \frac{k_{\mathrm B}T}{\epsilon_L} \right),

up to geometry-dependent constants. In one dimension it grows more strongly. Therefore the apparent threshold drifts with size rather than approaching a nonzero bulk critical temperature.

Three observations must be kept separate:

  • a single mode can be highly occupied in a finite sample;
  • a finite sample can have a long coherence length compared with its size;
  • a thermodynamic phase transition requires a specified infinite-system limit and nonanalytic bulk behavior.

The Thermodynamic Limit page develops that distinction systematically.

Fermions do not have the same excited-state-capacity problem. At zero temperature, all one-particle states are filled up to kFk_{\mathrm F}, with at most one fermion per internal state and momentum mode.

In one dimension, the occupied region is the interval

−kF≤k≤kF.-k_{\mathrm F} \le k \le k_{\mathrm F}.

Counting both Fermi points gives

n1=gkFπ,n_1 = \frac{gk_{\mathrm F}}{\pi},

so

kF(1)=πn1g,k_{\mathrm F}^{(1)} = \frac{\pi n_1}{g},

and

ϵF(1)=ℏ2π2n122mg2.\epsilon_{\mathrm F}^{(1)} = \frac{\hbar^2\pi^2n_1^2}{2mg^2}.

In two dimensions, the occupied region is a disk. Its area is πkF2\pi k_{\mathrm F}^2, so

n2=gkF24π,n_2 = \frac{gk_{\mathrm F}^2}{4\pi},

and

kF(2)=4πn2g.k_{\mathrm F}^{(2)} = \sqrt{\frac{4\pi n_2}{g}}.

The corresponding Fermi energy is

ϵF(2)=2πℏ2n2mg.\epsilon_{\mathrm F}^{(2)} = \frac{2\pi\hbar^2n_2}{mg}.

The one-dimensional Fermi surface consists of the two points ±kF\pm k_{\mathrm F}. The two-dimensional Fermi surface is a circle for an isotropic quadratic dispersion. These are boundaries in momentum space, not material surfaces in real space.

The Fermi Momentum and Fermi Energy page owns the general dd-dimensional counting framework, while Fermi Surface develops the geometric and low-energy meaning of that boundary.

For a quadratic dispersion in dd dimensions,

E0N=dd+2ϵF,\frac{E_0}{N} = \frac{d}{d+2} \epsilon_{\mathrm F},

and scale invariance gives

P0Vd=2dE0.P_0V_d = \frac{2}{d}E_0.

Thus a one-dimensional ideal Fermi gas obeys

E0N=ϵF3,P0=23n1ϵF,\frac{E_0}{N} = \frac{\epsilon_{\mathrm F}}{3}, \qquad P_0 = \frac{2}{3}n_1\epsilon_{\mathrm F},

where one-dimensional pressure has units of force. In two dimensions,

E0N=ϵF2,P0=12n2ϵF.\frac{E_0}{N} = \frac{\epsilon_{\mathrm F}}{2}, \qquad P_0 = \frac{1}{2}n_2\epsilon_{\mathrm F}.

At low temperature, the ideal-gas heat capacity is

CVNkB=π2d6TTF+O(T3TF3),\frac{C_V}{Nk_{\mathrm B}} = \frac{\pi^2d}{6} \frac{T}{T_{\mathrm F}} + O\left( \frac{T^3}{T_{\mathrm F}^3} \right),

provided the Fermi energy lies well away from a singular band edge and the density of states is smooth over the thermal window. In one dimension this condition matters because D1(ϵ)D_1(\epsilon) diverges at zero energy.

Because the two-dimensional quadratic density of states is constant, the number equation can be inverted exactly:

n2=gλT2ln⁡(1+z).n_2 = \frac{g}{\lambda_T^2} \ln(1+z).

Using

n2λT2g=ϵFkBT,\frac{n_2\lambda_T^2}{g} = \frac{\epsilon_{\mathrm F}}{k_{\mathrm B}T},

one obtains

μ(T)=kBTln⁡(eϵF/kBT−1).\mu(T) = k_{\mathrm B}T \ln\left( e^{\epsilon_{\mathrm F}/k_{\mathrm B}T} - 1 \right).

Equivalently,

μ(T)=ϵF+kBTln⁡(1−e−ϵF/kBT).\mu(T) = \epsilon_{\mathrm F} + k_{\mathrm B}T \ln\left( 1-e^{-\epsilon_{\mathrm F}/k_{\mathrm B}T} \right).

The low-temperature correction is exponentially small. A naive power-series Sommerfeld shift vanishes because D2′(ϵ)=0D_2'(\epsilon)=0. The Sommerfeld Expansion page explains why this does not mean that every thermal correction vanishes.

In the grand canonical ensemble,

Var⁡(N)=kBT(∂N∂μ)T,Vd.\operatorname{Var}(N) = k_{\mathrm B}T \left( \frac{\partial N}{\partial\mu} \right)_{T,V_d}.

For a uniform ideal Bose gas above any finite-size crossover,

Var⁡(N)N=Li⁡d/2−1(z)Li⁡d/2(z).\frac{\operatorname{Var}(N)}{N} = \frac{ \operatorname{Li}_{d/2-1}(z) }{ \operatorname{Li}_{d/2}(z) }.

As z→1−z\to1^-, the numerator weights the infrared modes even more strongly than the density does. In two dimensions it contains

Li⁡0(z)=z1−z,\operatorname{Li}_0(z) = \frac{z}{1-z},

and in one dimension it contains Li⁡−1/2(z)\operatorname{Li}_{-1/2}(z). Both diverge strongly. Low-dimensional ideal Bose gases are therefore highly compressible and fluctuation sensitive in the degenerate regime.

This formula describes total-number fluctuations in a grand canonical ideal gas. Fluctuations in a finite imaging cell, a canonical sample, or an interacting fluid require the appropriate correlation function and ensemble. It would be a mistake to infer a literal laboratory variance without matching those conditions.

Fermi occupation fluctuations instead contain f(1−f)f(1-f). Pauli blocking suppresses fluctuations deep inside the filled sea, leaving a thermal shell near the Fermi boundary. Low dimension still changes the available phase space, but it does not create the Bose infrared pileup.

State counting explains the ideal gas. Interacting low-dimensional fluids add a second infrared mechanism: collective phase fluctuations.

Suppose amplitude fluctuations are weak enough that a bosonic field can be written locally as

ψ(r)≃nseiθ(r).\psi(\mathbf r) \simeq \sqrt{n_s} e^{i\theta(\mathbf r)}.

At long wavelengths, a phase-only free energy has the form

Fθ=ℏ2ns2m∫ddr ∣∇θ∣2,F_\theta = \frac{\hbar^2n_s}{2m} \int d^dr\, |\boldsymbol\nabla\theta|^2,

where nsn_s is the superfluid number density in the active dimension. Thermal equipartition then gives the schematic variance

⟨[θ(r)−θ(0)]2⟩=2mkBTℏ2ns∫ddk(2π)d1−cos⁡(k⋅r)k2.\left\langle [\theta(\mathbf r)-\theta(\mathbf 0)]^2 \right\rangle = \frac{2mk_{\mathrm B}T}{\hbar^2n_s} \int\frac{d^dk}{(2\pi)^d} \frac{1-\cos(\mathbf k\cdot\mathbf r)}{k^2}.

The integral has different infrared behavior in each dimension:

DimensionLong-distance phase variancePhase-only coherence
d=1d=1proportional to ∣x∣\lvert x\rvert at T>0T>0exponential decay
d=2d=2proportional to ln⁡r\ln r in the vortex-bound regimealgebraic decay
d=3d=3approaches a finite infrared contributiontrue long-range phase order can survive

In one dimension, the Gaussian phase theory gives

⟨[θ(x)−θ(0)]2⟩≃mkBTℏ2ns∣x∣,\left\langle [\theta(x)-\theta(0)]^2 \right\rangle \simeq \frac{mk_{\mathrm B}T}{\hbar^2n_s}|x|,

and therefore

g(1)(x)∝exp⁡(−∣x∣ℓϕ),ℓϕ=2ℏ2nsmkBT,g^{(1)}(x) \propto \exp\left( -\frac{|x|}{\ell_\phi} \right), \qquad \ell_\phi = \frac{2\hbar^2n_s}{mk_{\mathrm B}T},

within this simple phase-only approximation.

In two dimensions,

⟨[θ(r)−θ(0)]2⟩≃mkBTπℏ2nsln⁡(ra),\left\langle [\theta(\mathbf r)-\theta(\mathbf 0)]^2 \right\rangle \simeq \frac{mk_{\mathrm B}T}{\pi\hbar^2n_s} \ln\left( \frac{r}{a} \right),

so

g(1)(r)∝(ar)η,η=mkBT2πℏ2ns.g^{(1)}(r) \propto \left( \frac{a}{r} \right)^\eta, \qquad \eta = \frac{mk_{\mathrm B}T}{2\pi\hbar^2n_s}.

Here aa is a short-distance cutoff. Vortices, density fluctuations, scale dependence of nsn_s, and finite geometry must be added for a full theory. The Correlation Functions Overview page owns the general interpretation of exponential, algebraic, and long-range correlation behavior.

Hohenberg used a Bogoliubov inequality to exclude ordinary finite-temperature long-range order for broad classes of homogeneous one- and two-dimensional Bose and Fermi systems with sufficiently short-range interactions. For a Bose fluid, the result formalizes the destructive role of long-wavelength fluctuations.

The hypotheses matter. The result does not say that:

  • every finite low-dimensional cloud lacks a large coherent mode;
  • two-dimensional superfluidity is impossible;
  • discrete symmetries cannot order in two dimensions;
  • long-range interactions obey the same conclusion without further analysis;
  • a harmonic trap has the same density of states as a uniform box;
  • zero-temperature quantum phases are covered by a finite-temperature theorem.

It is therefore safer to state the result with its physical setting than to invoke “the Mermin–Wagner theorem” as a universal slogan. Hohenberg’s continuum-fluid argument, the Mermin–Wagner spin result, and later generalizations are related but not interchangeable statements.

An interacting two-dimensional Bose fluid can be superfluid without possessing a nonzero infinite-distance limit of g(1)(r)g^{(1)}(r). In the low-temperature phase, vortex–antivortex pairs are bound and correlations decay algebraically. At the Berezinskii–Kosterlitz–Thouless transition, free vortices proliferate and correlations become short ranged.

For a homogeneous bosonic fluid, the universal jump can be written as

ℏ2ns(TBKT−)m=2kBTBKTπ.\frac{\hbar^2n_s(T_{\mathrm{BKT}}^-)}{m} = \frac{2k_{\mathrm B}T_{\mathrm{BKT}}}{\pi}.

Equivalently,

ns(TBKT−)λTBKT2=4.n_s(T_{\mathrm{BKT}}^-) \lambda_{T_{\mathrm{BKT}}}^2 = 4.

This is a statement about the renormalized long-wavelength superfluid density just below the transition, not the total density and not a generic finite trapped cloud. In an inhomogeneous trap, local critical behavior is rounded into a crossover and interpreted with finite-size and local-density analysis.

The key distinction is

Bose–Einstein condensation≠superfluidity≠BKT quasi-long-range order.\text{Bose–Einstein condensation} \ne \text{superfluidity} \ne \text{BKT quasi-long-range order}.

These notions can coexist in some systems, but none is a definition of the others.

One Dimension and the Luttinger-Liquid Preview

Section titled “One Dimension and the Luttinger-Liquid Preview”

In one dimension, particles cannot pass one another without meeting. This kinematic fact magnifies the consequences of short-range interactions and makes collective descriptions unusually powerful.

A canonical continuum model is the repulsive contact Bose gas,

H=∫dx [ℏ22m(∂xψ†)(∂xψ)+g1D2ψ†ψ†ψψ].H = \int dx\, \left[ \frac{\hbar^2}{2m} (\partial_x\psi^\dagger) (\partial_x\psi) + \frac{g_{1\mathrm D}}{2} \psi^\dagger\psi^\dagger\psi\psi \right].

Its dimensionless interaction parameter is

γ=mg1Dℏ2n1.\gamma = \frac{mg_{1\mathrm D}}{\hbar^2n_1}.

At fixed coupling, lowering the line density increases γ\gamma. The dilute limit can therefore be strongly correlated rather than more ideal. The weakly interacting regime and the strongly repulsive Tonks–Girardeau regime are separated by a continuous crossover in this model; Lieb–Liniger Model Preview owns the exact finite-coupling spectrum, thermodynamics, observables, and benchmark.

For broad classes of gapless one-dimensional interacting systems, the low-energy degrees of freedom are collective density and phase waves rather than long-lived particle-like excitations at a Fermi surface. A Luttinger-liquid description organizes correlation exponents through a sound velocity and a dimensionless interaction parameter. It also predicts phenomena such as spin–charge separation in suitable multicomponent fermion systems.

This is an orientation, not a derivation. Luttinger Liquid Preview owns bosonization, the effective Hamiltonian, correlation exponents, perturbations, and breakdown conditions.

Freezing transverse motion does not justify carrying a three-dimensional contact coupling into a lower-dimensional Hamiltonian unchanged.

Let the three-dimensional field approximately factorize as

Ψ(r)≃ϕ⊥(r⊥)ψ(x)\Psi(\mathbf r) \simeq \phi_\perp(\mathbf r_\perp) \psi(x)

for a quasi-one-dimensional gas. A simple projection of a weak three-dimensional contact interaction gives

g1D(proj)=g3D∫d2r⊥ ∣ϕ⊥(r⊥)∣4.g_{1\mathrm D}^{(\mathrm{proj})} = g_{3\mathrm D} \int d^2r_\perp\, |\phi_\perp(\mathbf r_\perp)|^4.

For a harmonic transverse ground state with oscillator length

a⊥=ℏmω⊥,a_\perp = \sqrt{\frac{\hbar}{m\omega_\perp}},

this projection yields

g1D(proj)=g3D2πa⊥2.g_{1\mathrm D}^{(\mathrm{proj})} = \frac{g_{3\mathrm D}}{2\pi a_\perp^2}.

Similarly, a quasi-two-dimensional projection gives

g2D(proj)=g3D∫dz ∣ϕ0(z)∣4=g3D2πazg_{2\mathrm D}^{(\mathrm{proj})} = g_{3\mathrm D} \int dz\, |\phi_0(z)|^4 = \frac{g_{3\mathrm D}}{\sqrt{2\pi}a_z}

for a harmonic ground mode with length aza_z.

These projected formulas are only leading weak-coupling estimates. Exact low-energy scattering includes virtual transverse excitations. In quasi-one-dimensional geometry this produces a confinement-induced resonance; in two dimensions the effective scattering amplitude depends logarithmically on collision energy or density. The coupling must therefore be matched to the actual confinement and scattering convention before quantitative use.

Uniform and Trapped Dimensionality Are Different Questions

Section titled “Uniform and Trapped Dimensionality Are Different Questions”

Changing the external potential changes the density of states independently of changing spatial dimension.

For a dd-dimensional harmonic trap in the semiclassical regime,

Dtrap,d(ϵ)∝ϵd−1.D_{\mathrm{trap},d}(\epsilon) \propto \epsilon^{d-1}.

The ideal Bose excited-state integral then has low-energy behavior

∫0dϵ ϵd−2.\int_0 d\epsilon\, \epsilon^{d-2}.

A two-dimensional harmonic trap has a finite semiclassical excited-state capacity and can support an ideal-gas condensation transition in its appropriate trap thermodynamic limit. A one-dimensional harmonic trap is marginal and retains a logarithmic divergence in the standard limit.

This does not contradict the uniform result. The box and the trap are different sequences of spectra. A finite trapped gas also rounds any sharp limiting singularity.

The Quantum Gases in Traps page owns harmonic spectra, local-density profiles, imaging observables, and the trap thermodynamic limit.

Suppose a three-dimensional harmonic trap has frequencies

ωx≪ωy,ωz.\omega_x \ll \omega_y, \omega_z.

If

kBT,ϵF,μint,ℏωprobe≪ℏωy,ℏωz,k_{\mathrm B}T, \epsilon_{\mathrm F}, \mu_{\mathrm{int}}, \hbar\omega_{\mathrm{probe}} \ll \hbar\omega_y, \hbar\omega_z,

then motion is effectively one dimensional along xx. For a quasi-two-dimensional gas, only one direction needs to be frozen.

The inequalities are not sharp phase boundaries. As an active scale approaches a transverse gap:

  • excited transverse subbands acquire thermal population;
  • the density of states develops step-like or threshold structure;
  • scattering can couple open and closed transverse channels;
  • thermodynamic observables cross over continuously between dimensions;
  • local dimensionality can vary across an inhomogeneous cloud.

A reliable dimensionality claim should state the confinement frequencies and compare them with all relevant energy scales. A cloud aspect ratio alone is insufficient evidence.

Low-dimensional atomic gases are commonly produced with anisotropic magnetic or optical confinement, atom-chip waveguides, or optical lattices that form arrays of tubes or planes.

Useful observables include:

  • transverse populations: spectroscopy or thermometry can test whether excited frozen modes are occupied;
  • in-situ density profiles: an equation of state reveals the dimension-dependent density of states and compressibility;
  • momentum distributions: narrow central features diagnose long coherence lengths but do not alone prove true condensation;
  • interference correlations: algebraic or exponential coherence can be distinguished over an accessible range;
  • density correlations: bunching, antibunching, and local pair suppression probe statistics and interaction strength;
  • collective modes: dimensional crossover and scale anomalies modify breathing and sound modes;
  • vortex statistics: bound pairs and free vortices are central signatures of two-dimensional BKT physics.

Landmark cold-atom experiments realized strongly correlated one-dimensional bosons in optical tube arrays and observed BKT-type behavior in trapped two-dimensional Bose gases. Such experiments are finite, inhomogeneous, and measured over finite dynamic range. Their comparison with homogeneous thermodynamic statements therefore requires the physical audit in Finite-Size Effects together with trap-specific and imaging analysis.

Calling every elongated trap one dimensional

Section titled “Calling every elongated trap one dimensional”

Geometry is not enough. Compare temperature, interaction, Fermi, and drive energies with the transverse gap.

A line density, area density, and volume density have different units. So do their pressures and couplings.

Concluding that no BEC means no degeneracy

Section titled “Concluding that no BEC means no degeneracy”

The phase-space density ndλTd/gn_d\lambda_T^d/g can be large in one or two dimensions even without a uniform ideal-gas condensation transition.

Concluding that no BEC means no superfluidity

Section titled “Concluding that no BEC means no superfluidity”

An interacting two-dimensional gas can have BKT superfluidity and algebraic order without ordinary long-range one-body order.

Treating a finite-size crossover as a bulk transition

Section titled “Treating a finite-size crossover as a bulk transition”

Large ground-mode occupation or system-wide coherence in a finite cloud does not by itself establish a nonanalytic thermodynamic limit.

Applying the Hohenberg or Mermin–Wagner name without hypotheses

Section titled “Applying the Hohenberg or Mermin–Wagner name without hypotheses”

Dimension, temperature, interaction range, symmetry, homogeneity, and the observable called “order” all matter.

Assuming weak density means weak interactions in one dimension

Section titled “Assuming weak density means weak interactions in one dimension”

For the Lieb–Liniger parameter γ=mg1D/(ℏ2n1)\gamma=mg_{1\mathrm D}/(\hbar^2n_1), lowering density strengthens the dimensionless interaction.

Using a projected coupling as an exact coupling

Section titled “Using a projected coupling as an exact coupling”

Virtual transverse excitations renormalize low-dimensional scattering and can generate confinement-induced resonances.

Treating Fermi points as ordinary quasiparticle surfaces

Section titled “Treating Fermi points as ordinary quasiparticle surfaces”

The ideal one-dimensional gas has Fermi points, but interactions can replace the Landau-quasiparticle picture with Luttinger-liquid collective modes.

For quadratic continuum particles,

Dd(ϵ)Vd=gΓ(d/2)(m2πℏ2)d/2ϵd/2−1.\frac{D_d(\epsilon)}{V_d} = \frac{g}{\Gamma(d/2)} \left( \frac{m}{2\pi\hbar^2} \right)^{d/2} \epsilon^{d/2-1}.

The main consequences are:

  • D1(ϵ)∝ϵ−1/2D_1(\epsilon)\propto\epsilon^{-1/2}, while D2(ϵ)D_2(\epsilon) is constant;
  • a uniform ideal Bose gas with quadratic dispersion has no finite-temperature condensation transition for d≤2d\le2;
  • finite systems and traps can show strong occupation crossovers, and a two-dimensional harmonic trap changes the ideal condensation criterion;
  • one-dimensional Fermi seas terminate at two Fermi points, while two-dimensional isotropic seas have circular boundaries;
  • thermal and phase fluctuations are enhanced at long wavelength;
  • interacting two-dimensional Bose fluids can undergo a BKT transition without ordinary long-range order;
  • generic gapless one-dimensional interacting fluids lead toward Luttinger-liquid behavior;
  • effective dimension and effective coupling must both be derived from the confinement scales.

Low-Dimensional Quantum Matter places these gas results in the broader materials ledger of subband occupancy, environmental dimension, electronic platforms, disorder, and dimensionality evidence.

Starting from periodic momentum quantization, derive Dd(ϵ)/VdD_d(\epsilon)/V_d for quadratic dispersion. Verify the explicit one- and two-dimensional formulas, including all factors of 22, π\pi, and ℏ\hbar.

Solution

The number of states inside a dd-dimensional momentum sphere is

Nd(k)=gVd(2π)dSd−1dkd.\mathcal N_d(k) = g\frac{V_d}{(2\pi)^d} \frac{S_{d-1}}{d}k^d.

Using Sd−1=2πd/2/Γ(d/2)S_{d-1}=2\pi^{d/2}/\Gamma(d/2) and k=2mϵ/ℏk=\sqrt{2m\epsilon}/\hbar gives

Nd(ϵ)Vd=gΓ(d/2+1)(m2πℏ2)d/2ϵd/2.\frac{\mathcal N_d(\epsilon)}{V_d} = \frac{g}{\Gamma(d/2+1)} \left( \frac{m}{2\pi\hbar^2} \right)^{d/2} \epsilon^{d/2}.

Differentiation yields

Dd(ϵ)Vd=gΓ(d/2)(m2πℏ2)d/2ϵd/2−1.\frac{D_d(\epsilon)}{V_d} = \frac{g}{\Gamma(d/2)} \left( \frac{m}{2\pi\hbar^2} \right)^{d/2} \epsilon^{d/2-1}.

For d=1d=1, Γ(1/2)=π\Gamma(1/2)=\sqrt\pi, so

D1L=gπℏm2ϵ.\frac{D_1}{L} = \frac{g}{\pi\hbar} \sqrt{\frac{m}{2\epsilon}}.

For d=2d=2, Γ(1)=1\Gamma(1)=1, giving

D2A=gm2πℏ2.\frac{D_2}{A} = \frac{gm}{2\pi\hbar^2}.

Exercise 2: Infrared convergence for a general dispersion

Section titled “Exercise 2: Infrared convergence for a general dispersion”

Let ϵk=Aks\epsilon_{\mathbf k}=Ak^s at small kk, with A>0A>0. Determine the condition on dd and ss for a uniform ideal Bose gas to have finite excited-state capacity at nonzero temperature.

Solution

The number of states below kk scales as kdk^d. Since k∝ϵ1/sk\propto\epsilon^{1/s},

D(ϵ)∝ϵd/s−1.D(\epsilon) \propto \epsilon^{d/s-1}.

At μ=0\mu=0, the Bose factor is proportional to 1/ϵ1/\epsilon at low energy. The excited-state integrand is therefore

D(ϵ)nB(ϵ)∝ϵd/s−2.D(\epsilon)n_{\mathrm B}(\epsilon) \propto \epsilon^{d/s-2}.

The integral at zero converges when its exponent is greater than −1-1:

ds−2>−1.\frac{d}{s}-2 > -1.

Hence

d>s.d > s.

For s=2s=2, this reduces to d>2d>2.

Exercise 3: Invert the two-dimensional Bose number equation

Section titled “Exercise 3: Invert the two-dimensional Bose number equation”

For a uniform ideal Bose gas in two dimensions, derive z(n2,T)z(n_2,T) and μ(n2,T)\mu(n_2,T). Find the leading highly degenerate asymptotic form of μ\mu.

Solution

Using Li⁡1(z)=−ln⁡(1−z)\operatorname{Li}_1(z)=-\ln(1-z),

n2λT2g=−ln⁡(1−z).\frac{n_2\lambda_T^2}{g} = -\ln(1-z).

Exponentiating gives

z=1−e−n2λT2/g.z = 1 - e^{-n_2\lambda_T^2/g}.

Therefore

μ=kBTln⁡(1−e−n2λT2/g).\mu = k_{\mathrm B}T \ln\left( 1-e^{-n_2\lambda_T^2/g} \right).

When n2λT2/g≫1n_2\lambda_T^2/g\gg1, use ln⁡(1−x)≃−x\ln(1-x)\simeq-x:

μ≃−kBTexp⁡(−n2λT2g).\mu \simeq -k_{\mathrm B}T \exp\left( -\frac{n_2\lambda_T^2}{g} \right).

It is exponentially close to zero but remains negative at finite temperature.

Exercise 4: One- and two-dimensional Fermi scales

Section titled “Exercise 4: One- and two-dimensional Fermi scales”

Derive kFk_{\mathrm F}, ϵF\epsilon_{\mathrm F}, and E0/NE_0/N for uniform ideal Fermi gases in one and two dimensions with internal multiplicity gg.

Solution

In one dimension, the occupied interval has length 2kF2k_{\mathrm F}:

n1=g2kF2π=gkFπ.n_1 = g\frac{2k_{\mathrm F}}{2\pi} = \frac{gk_{\mathrm F}}{\pi}.

Thus

kF=πn1g,ϵF=ℏ2π2n122mg2.k_{\mathrm F} = \frac{\pi n_1}{g}, \qquad \epsilon_{\mathrm F} = \frac{\hbar^2\pi^2n_1^2}{2mg^2}.

Integrating the quadratic energy symmetrically from −kF-k_{\mathrm F} to kFk_{\mathrm F} gives

E0N=ϵF3.\frac{E_0}{N} = \frac{\epsilon_{\mathrm F}}{3}.

In two dimensions,

n2=gπkF2(2π)2=gkF24π.n_2 = g\frac{\pi k_{\mathrm F}^2}{(2\pi)^2} = \frac{gk_{\mathrm F}^2}{4\pi}.

Hence

kF=4πn2g,ϵF=2πℏ2n2mg.k_{\mathrm F} = \sqrt{\frac{4\pi n_2}{g}}, \qquad \epsilon_{\mathrm F} = \frac{2\pi\hbar^2n_2}{mg}.

The constant density of states gives

E0N=ϵF2.\frac{E_0}{N} = \frac{\epsilon_{\mathrm F}}{2}.

Exercise 5: Phase fluctuations in one dimension

Section titled “Exercise 5: Phase fluctuations in one dimension”

Evaluate

I(x)=∫−∞∞dk2π1−cos⁡(kx)k2.I(x) = \int_{-\infty}^{\infty} \frac{dk}{2\pi} \frac{1-\cos(kx)}{k^2}.

Use it to recover the phase-only coherence length ℓϕ\ell_\phi.

Solution

Differentiate with respect to ∣x∣|x|, or use the standard Fourier integral:

I(x)=∣x∣2.I(x) = \frac{|x|}{2}.

Substitution into

⟨[θ(x)−θ(0)]2⟩=2mkBTℏ2nsI(x)\left\langle [\theta(x)-\theta(0)]^2 \right\rangle = \frac{2mk_{\mathrm B}T}{\hbar^2n_s}I(x)

gives

⟨[θ(x)−θ(0)]2⟩=mkBTℏ2ns∣x∣.\left\langle [\theta(x)-\theta(0)]^2 \right\rangle = \frac{mk_{\mathrm B}T}{\hbar^2n_s}|x|.

For Gaussian phase fluctuations,

⟨ei[θ(x)−θ(0)]⟩=exp⁡[−12⟨[θ(x)−θ(0)]2⟩].\left\langle e^{i[\theta(x)-\theta(0)]} \right\rangle = \exp\left[ -\frac{1}{2} \left\langle [\theta(x)-\theta(0)]^2 \right\rangle \right].

Therefore

g(1)(x)∝e−∣x∣/ℓϕ,ℓϕ=2ℏ2nsmkBT.g^{(1)}(x) \propto e^{-|x|/\ell_\phi}, \qquad \ell_\phi = \frac{2\hbar^2n_s}{mk_{\mathrm B}T}.

Exercise 6: Uniform plane versus harmonic plane

Section titled “Exercise 6: Uniform plane versus harmonic plane”

Explain why a uniform two-dimensional ideal Bose gas lacks a finite-temperature condensation transition while a semiclassical two-dimensional harmonic trap can have finite excited-state capacity.

Solution

For the uniform quadratic gas,

Dbox,2(ϵ)∝ϵ0.D_{\mathrm{box},2}(\epsilon) \propto \epsilon^0.

At μ=0\mu=0, multiplication by the Bose factor gives an integrand proportional to ϵ−1\epsilon^{-1}, which diverges logarithmically.

For a two-dimensional harmonic trap, semiclassical phase-space counting gives

Dtrap,2(ϵ)∝ϵ.D_{\mathrm{trap},2}(\epsilon) \propto \epsilon.

The low-energy Bose integrand is then proportional to a constant, so its lower limit converges. The trap changes the spectrum and therefore the thermodynamic limit; it is not merely a boundary around the uniform gas.

A gas has transverse frequency ω⊥\omega_\perp and active scales kBTk_{\mathrm B}T, ϵF\epsilon_{\mathrm F}, μint\mu_{\mathrm{int}}, and ℏωdrive\hbar\omega_{\mathrm{drive}}. Formulate a conservative dimensionless test for quasi-one-dimensional behavior. What conclusion follows if one ratio is of order unity?

Solution

Define

ηdim=max⁡(kBT,ϵF,∣μint∣,ℏωdrive)ℏω⊥.\eta_{\mathrm{dim}} = \frac{ \max( k_{\mathrm B}T, \epsilon_{\mathrm F}, |\mu_{\mathrm{int}}|, \hbar\omega_{\mathrm{drive}} ) }{ \hbar\omega_\perp }.

A conservative frozen-mode regime requires

ηdim≪1.\eta_{\mathrm{dim}} \ll 1.

If one ratio is of order unity, the corresponding process can populate or virtually mix transverse excitations. The system lies in a dimensional crossover for that observable or protocol, even if its equilibrium image remains highly elongated.

Classify each statement as correct, incorrect, or incomplete.

  1. “A uniform two-dimensional ideal Bose gas cannot be quantum degenerate.”
  2. “No finite-temperature BEC in two dimensions implies no two-dimensional superfluidity.”
  3. “A finite one-dimensional gas may have coherence across the whole sample.”
  4. “Projecting onto a transverse ground state always gives the exact effective coupling.”
Solution
  1. Incorrect. The phase-space density n2λT2/gn_2\lambda_T^2/g can be arbitrarily large. The absent feature is a uniform ideal-gas thermodynamic condensation transition, not degeneracy.
  2. Incorrect. An interacting two-dimensional Bose fluid can have BKT superfluidity with algebraic correlations and no ordinary infinite-distance condensate order.
  3. Correct with finite-size qualification. If the coherence length exceeds the system size, the sample can appear globally coherent even though the corresponding infinite system has no finite-temperature long-range order.
  4. Incorrect. Simple projection is a weak-coupling estimate. Virtual transverse excitation and low-dimensional scattering renormalize the coupling and can produce a confinement-induced resonance.
  1. K. Huang, Statistical Mechanics, 2nd ed. (Wiley, 1987) — ideal quantum gases, density-of-states methods, and thermodynamic limits.
  2. R. K. Pathria and P. D. Beale, Statistical Mechanics, 4th ed. (Academic Press, 2021) — Bose and Fermi functions in general dimension.
  3. L. Pitaevskii and S. Stringari, Bose–Einstein Condensation and Superfluidity (Oxford University Press, 2016) — trapped gases, low-dimensional fluctuations, and superfluidity.
  4. C. J. Pethick and H. Smith, Bose–Einstein Condensation in Dilute Gases, 2nd ed. (Cambridge University Press, 2008) — effective dimensionality and dilute-gas scales.
  5. P. C. Hohenberg, “Existence of Long-Range Order in One and Two Dimensions”, Physical Review 158, 383–386 (1967) — Bogoliubov-inequality restrictions on finite-temperature long-range order.
  6. V. L. Berezinskii, “Destruction of Long-Range Order in One-Dimensional and Two-Dimensional Systems Having a Continuous Symmetry Group,” Soviet Physics JETP 32, 493–500 (1971); 34, 610–616 (1972) — topological-defect mechanism and low-dimensional correlations.
  7. J. M. Kosterlitz and D. J. Thouless, “Ordering, metastability and phase transitions in two-dimensional systems”, Journal of Physics C 6, 1181–1203 (1973) — vortex-unbinding transition.
  8. D. R. Nelson and J. M. Kosterlitz, “Universal Jump in the Superfluid Density of Two-Dimensional Superfluids”, Physical Review Letters 39, 1201–1205 (1977) — universal superfluid-density jump.
  9. E. H. Lieb and W. Liniger, “Exact Analysis of an Interacting Bose Gas. I. The General Solution and the Ground State”, Physical Review 130, 1605–1616 (1963) — exact repulsive one-dimensional contact Bose gas.
  10. M. Girardeau, “Relationship between Systems of Impenetrable Bosons and Fermions in One Dimension”, Journal of Mathematical Physics 1, 516–523 (1960) — Bose–Fermi mapping in the hard-core limit.
  11. M. Olshanii, “Atomic Scattering in the Presence of an External Confinement and a Gas of Impenetrable Bosons”, Physical Review Letters 81, 938–941 (1998) — confinement-induced renormalization and resonance in quasi-one-dimensional gases.
  12. M. A. Cazalilla, R. Citro, T. Giamarchi, E. Orignac, and M. Rigol, “One dimensional bosons: From condensed matter systems to ultracold gases”, Reviews of Modern Physics 83, 1405–1466 (2011) — exact models, Luttinger liquids, and experiments.
  13. V. Bagnato and D. Kleppner, “Bose–Einstein condensation in low-dimensional traps”, Physical Review A 44, 7439–7441 (1991) — how trapping modifies the dimensional condensation criterion.
  14. B. Paredes et al., “Tonks–Girardeau gas of ultracold atoms in an optical lattice”, Nature 429, 277–281 (2004) — experimental realization of strongly correlated bosons in tube arrays.
  15. T. Kinoshita, T. Wenger, and D. S. Weiss, “Observation of a One-Dimensional Tonks–Girardeau Gas”, Science 305, 1125–1128 (2004) — complementary one-dimensional hard-core-boson experiment.
  16. Z. Hadzibabic, P. Krüger, M. Cheneau, B. Battelier, and J. Dalibard, “Berezinskii–Kosterlitz–Thouless crossover in a trapped atomic gas”, Nature 441, 1118–1121 (2006) — coherence and vortex signatures in a trapped two-dimensional Bose gas.
  17. P. Krüger, Z. Hadzibabic, and J. Dalibard, “Critical Point of an Interacting Two-Dimensional Atomic Bose Gas”, Physical Review Letters 99, 040402 (2007) — critical atom number and the failure of an ideal-gas BEC description for an interacting 2D gas.