Skip to content

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.

We compare the one- and two-dimensional ideal gases, then explain how collective phase fluctuations change interacting fluids. Ultracold atoms explains how transverse spectroscopy and trap calibration establish the effective dimension of a laboratory gas.

Helpful background. Ideal Bose and Fermi gases give the occupation laws; trapped gases explains transverse levels and continuum limits.

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 energies ℏω\hbar\omega, drive bandwidths expressed as energies, 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, with equally populated identical components when a common chemical potential is used. One momentum state occupies volume (2π)d/Vd(2\pi)^d/V_d in k\mathbf k-space. The leading continuum 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.

Integrated Bose low-energy weights versus the lower energy cutoff for one, two, and three dimensions

The dimensionless low-energy integral Id(δ)=∫δ1ud/2−2 duI_d(\delta)=\int_\delta^1 u^{d/2-2}\,du isolates infrared convergence; overall density and temperature factors are omitted. As the cutoff decreases toward the left, the one-dimensional curve diverges as δ−1/2\delta^{-1/2} and the two-dimensional curve grows logarithmically. The three-dimensional curve approaches the finite value two. Both axes are logarithmic.

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.

For d=1,2d=1,2 at finite density and T>0T>0, the Bose relation has z<1z<1 strictly. For d>2d>2 at z=1z=1, it gives only the saturated excited density; a condensate must be added separately. NN in the grand-canonical formulas is a mean, and the Fermi polylogarithm uses its real Fermi-integral branch beyond z=1z=1. 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 derives the transition when it exists. 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(\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.

The ideal zero-temperature one-dimensional sea is an interval, bounded by ±kF\pm k_{\mathrm F}. The two-dimensional isotropic sea is a disk, bounded by a circle. For equally populated components, the counts are:

Dimension Number density Fermi energy
One n1=gkF/πn_1=gk_{\mathrm F}/\pi ϵF=ℏ2π2n12/(2mg2)\epsilon_{\mathrm F}=\hbar^2\pi^2n_1^2/(2mg^2)
Two n2=gkF2/(4π)n_2=gk_{\mathrm F}^2/(4\pi) ϵF=2πℏ2n2/(mg)\epsilon_{\mathrm F}=2\pi\hbar^2n_2/(mg)

Fermi scales derives the general count; Fermi surfaces explains how nearby particle–hole excitations carry small energy. These boundaries lie in wave-vector space. They are not surfaces bounding the gas in real space.

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.

Exponentially small chemical-potential shift in two dimensions

Section titled “Exponentially small chemical-potential shift in two dimensions”

The constant density of states gives n2=(g/λT2)ln⁡(1+z)n_2=(g/\lambda_T^2)\ln(1+z). At fixed density its exact inversion is

μ(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).

Thus the low-temperature chemical-potential shift is exponentially small. The Sommerfeld expansion explains why a vanishing algebraic shift does not eliminate the nonzero heat capacity.

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)≃nqeiθ(r).\psi(\mathbf r) \simeq \sqrt{n_{\mathrm q}} e^{i\theta(\mathbf r)}.

Here nqn_{\mathrm q} is a local amplitude or quasicondensate density, not automatically a condensate density or a superfluid density. The energetic cost of a slowly varying phase is a separate quantity. Write the phase-only free energy as

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

where ℏ2nstiff/m\hbar^2n_{\mathrm{stiff}}/m is the phase stiffness at the length scale being described. A microscopic cutoff is implicit. For the classical, thermally populated long-wavelength modes, equipartition gives

⟨[θ(r)−θ(0)]2⟩=2mkBTℏ2nstiff∫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_{\mathrm{stiff}}} \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:

Dimension Long-distance phase variance Phase-only coherence
d=1d=1 proportional to ∣x∣\lvert x\rvert at T>0T>0 exponential decay
d=2d=2 proportional to ln⁡r\ln r in the vortex-bound regime algebraic decay
d=3d=3 approaches a finite infrared contribution true long-range phase order can survive

In one dimension, the Gaussian phase theory gives

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

and therefore

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

within this local-stiffness approximation. In one dimension at positive temperature, this parameter must not automatically be identified with a nonzero infinite-system equilibrium helicity modulus: phase slips and winding sectors matter at long scales. Quantum fluctuations at zero temperature require the Luttinger description rather than this thermal equipartition integral.

In two dimensions,

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

so

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

Here aa is a short-distance cutoff. In a two-dimensional vortex-bound fluid, using the renormalized long-distance stiffness identifies nstiffn_{\mathrm{stiff}} with the superfluid density nsn_s in the BKT formula below. Vortices, density fluctuations, scale dependence of nstiffn_{\mathrm{stiff}}, and finite geometry must be added for a full theory. The Correlation Functions Overview page explains exponential, algebraic, and long-range correlation behavior.

Hohenberg used a Bogoliubov inequality to exclude the relevant continuous-symmetry long-range order at positive temperature in broad classes of homogeneous one- and two-dimensional systems under its interaction and regularity assumptions. For fermion fluids the superconducting application concerns pair order, not macroscopic occupation of a single fermion mode. 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. A finite trapped cloud rounds the transition and requires finite-size and local-density analysis. An appropriate trap thermodynamic limit can still support sharp critical behavior.

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 develops the exact finite-coupling model and thermodynamics.

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.

Luttinger liquids develops the effective Hamiltonian and correlation exponents, including when relevant perturbations invalidate the gapless description.

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. For a short-range three-dimensional interaction represented at leading order by g3D=4πℏ2as/mg_{3\mathrm D}=4\pi\hbar^2a_s/m, with scattering length asa_s small compared with the confinement length, simple transverse projection 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 derives harmonic spectra, local-density profiles and imaging formulas.

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 comparison of lengths and energy scales in Finite-Size Effects together with trap-specific and imaging analysis.

Infrared convergence for a general dispersion

Section titled “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.

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ℏ2nstiffI(x)\left\langle [\theta(x)-\theta(0)]^2 \right\rangle = \frac{2mk_{\mathrm B}T}{\hbar^2n_{\mathrm{stiff}}}I(x)

gives

⟨[θ(x)−θ(0)]2⟩=mkBTℏ2nstiff∣x∣.\left\langle [\theta(x)-\theta(0)]^2 \right\rangle = \frac{mk_{\mathrm B}T}{\hbar^2n_{\mathrm{stiff}}}|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ℏ2nstiffmkBT.g^{(1)}(x) \propto e^{-|x|/\ell_\phi}, \qquad \ell_\phi = \frac{2\hbar^2n_{\mathrm{stiff}}}{mk_{\mathrm B}T}.

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 simple energy-separation test no longer guarantees a frozen mode. Actual excitation also depends on coupling matrix elements and selection rules. Even well below the gap, virtual excitations can renormalize the effective interaction without populating a transverse excited state.

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.