High-Momentum Tails in Quantum Systems

• High-momentum tails (HMTs) are asymptotic power-law decays (often k⁻⁴) in particle momentum distributions reflecting universal short-range correlations in many-body systems.
• They are derived using renormalization group and operator product expansion methods, isolating two-body physics into a contact parameter that links theory to observables.
• HMTs are observed in diverse systems—from ultracold atomic gases to finite nuclei and high-energy collisions—with experimental probes and theoretical models addressing equilibrium and dynamical phenomena.

High-momentum tails (HMTs) refer to the asymptotic, algebraic decay at large momentum in single-particle momentum distributions and related correlation functions, observed across a range of many-body quantum systems, including ultracold atomic gases, nuclear systems, strongly correlated electron gases, and effective quantum field theories. In such systems, the momentum distribution $n(k)$ at high momentum $k$, far above the typical Fermi or thermal scale, generically exhibits a power-law decay, often of the form $n(k)\sim C/k^4$, where $C$ is a state-dependent "contact" parameter encoding short-range correlations. This universal tail reflects the underlying short-distance operator structure, emergent from scale separation and contact interactions, and is robust across dimensions and statistics.

1. Universal Factorization and Operator Structure

At the core of the HMTs is an operator factorization principle formally derived via renormalization group (RG) or operator product expansion (OPE) techniques. For any many-body state $|\Psi\rangle$ of a system with Hamiltonian $H$ constructed with (or RG-evolved to) a low-energy resolution scale $\Lambda$ and high momenta $q\gg\Lambda$, the single-particle momentum distribution satisfies

$$n(q) = \langle\Psi|a_q^\dagger a_q|\Psi\rangle \simeq F_n(q;\Lambda)\,\langle\Psi|{\hat C}|\Psi\rangle_\Lambda,$$

where $F_n(q;\Lambda)$ is a universal function of $q$ controlled entirely by two-body physics, while $\langle\Psi|{\hat C}|\Psi\rangle_\Lambda$ is a matrix element of a "contact operator" depending only on the low-momentum structure of $|\Psi\rangle$ (Bogner et al., 2012). Analogous factorization applies to the static structure factor and other local observables: $$S_{\uparrow\downarrow}(q) \simeq F_S(q;\Lambda) \langle\Psi|{\hat C}|\Psi\rangle_\Lambda,$$ where $F_S(q;\Lambda)$ is constructed from universal high-momentum coefficients.

The universal functions are calculable, e.g., via

$$F_n(q;\Lambda) = \gamma^2(q;\Lambda), \quad \gamma(q;\Lambda) = -\int_{|\vec{q'}|>\Lambda} d^3q' \langle q|[Q_\Lambda H Q_\Lambda]^{-1}|q'\rangle V(q',0),$$

for two-body interactions $V$. In paradigmatic cases:

• Unitary Fermi Gas: $n(q)\simeq C/q^4$, with $C$ identified as Tan’s contact.
• Electron Gas: $n(q)\sim 1/q^8$ (Kimball tail), reflecting Coulomb interactions.

All many-body dependence is isolated in the "contact" $\langle\Psi|{\hat C}|\Psi\rangle_\Lambda$, paralleling the Wilson coefficient decomposition of local operators in field theory OPE and enabling direct inference of high-momentum observables from low-energy many-body theory (Bogner et al., 2012).

2. HMTs in Cold Atoms and Strongly Correlated Gases

For systems of quantum particles with short-range (contact) interactions, such as ultracold atomic Fermi gases and 1D mixtures, HMTs acquire a direct signature in the measured $n(k)$: $n(k) \simeq \frac{C}{k^4}\quad (k \gg k_F),$ where $C$ is Tan’s contact, encoding the probability of finding two particles at vanishing separation (Decamp et al., 2016). The universality of this $k^{-4}$ decay has been observed in both bosonic and fermionic systems and is linked to the non-analytic cusp condition the many-body wavefunction obeys as two coordinates approach, an exact consequence of the zero-range boundary condition.

In multicomponent 1D fermionic mixtures with $SU(\kappa)$ symmetry, the component-resolved contacts $C_\sigma$ characterize the magnetic symmetry of the many-body state. In the fermionized limit, the value of $C_\sigma$ is controlled by the Young tableau of the permutation symmetry, making $k^{-4}$ tail measurements a probe of quantum magnetism (Decamp et al., 2016). Analytical expressions for $C_\sigma$ reveal $N^{5/2}$ scaling with particle number and explicitly encode temperature and interaction-strength dependence. Under local density approximation and Bethe Ansatz, LDA contacts agree with matrix-product-state results to high precision, supporting their universality.

Out-of-equilibrium dynamics can generate anomalously large HMT amplitudes, as demonstrated for Bose-Einstein condensates with dilute impurities—where $1/k^4$ tails appear with amplitudes exceeding the two-body contact equilibrium value by orders of magnitude due to impurity-bath scattering during time-of-flight expansion (Cayla et al., 2022). This indicates that HMTs can encode both equilibrium and dynamical information.

3. HMTs in Finite Nuclei and Nuclear Matter

In nuclear many-body systems, HMTs in $n(k)$ stem predominantly from the dominant neutron-proton tensor component of the nucleon-nucleon interaction, supplemented by a short-range repulsive core. State-of-the-art ab initio approaches, such as extended Brueckner-Hartree-Fock (EBHF) and lowest order cluster (LOC) expansions, demonstrate that, irrespective of isospin asymmetry, approx. $17\%-18\%$ of nucleons in finite nuclei reside in the HMT ($p>300$ MeV/$c$), with

$$n_A^\tau(k) = n_{\text{SD}}^\tau(k) + \delta n_A^\tau(k), \qquad \tau = n,p$$

where $n_{\text{SD}}^\tau$ is mean-field Slater determinant and $\delta n_A^\tau$ encodes correlated SRC strength (Fan et al., 2022).

Universal features include:

• Equal high-momentum neutron and proton occupation: $(N/Z)_{\text{high}}\simeq1$ even for large $N/Z$ nuclei.
• HMT fraction is robust across nuclear mass and isospin.
• Tensor correlations responsible for HMTs are enhanced at the nuclear surface.

The resulting $k^{-4}$ (or exponential) tail at high $k$ is in agreement with high-precision electron $(e,e'pN)$ knockout measurements and heavy-ion photon production observables. Isospin-dependent correction factors are required to enforce quantitative agreement with data, primarily ensuring high-$k$ proton and neutron fractions are equal.

4. Experimental Probes and Transport Phenomenology

HMTs in nuclei and nuclear matter are probed both in electron scattering and in heavy-ion collision observables:

• Photon Production: Energetic photons above 30 MeV in central heavy-ion reactions ($^{86}$Kr+$^{124}$Sn at 25 MeV/u) are sensitive to the HMT fraction, as high-momentum nucleon pairs ($np$) dominate photon-bremsstrahlung. Comparison of data to IBUU transport model simulations incorporating adjustable HMT fractions yields quantitative constraints, e.g., an HMT fraction of $\sim$15% is favored, with zero-HMT models statistically disfavored (Qin et al., 2023).
• Observables: The shape of photon energy and transverse-momentum spectra (especially yield ratios at $E_\gamma>100$ MeV or $p_T>150$ MeV) is highly sensitive to the HMT shape/exponent (e.g., $1/k^4$, $1/k^6$, $1/k^9$), while angular spectra are not (Guo et al., 2021). Additionally, the shape of the $\pi^-/\pi^+$ ratio and the momentum-resolved $n/p$ ratio in heavy-ion collision products provide independent sensitivity to HMT parameterizations (Yang et al., 2018).
• Strategy: Extraction of HMT fraction proceeds by maximizing likelihood between data and transport modeling as a function of $R_{\mathrm{HMT}}$, filtering the output with realistic detector response.

5. Finite-Size and Trap Effects on HMTs

In confined geometries, e.g., 1D gases in a box potential, the HMT can display additional structure beyond the canonical $1/k^4$ decay. For strongly interacting systems, the presence of rigid boundaries introduces two finite-size corrections: $$n(k) \sim \frac{{\mathcal C}^\mathrm{mix}}{k^4} + \frac{\mathcal{B}_N}{k^4} + \frac{(-1)^{N+1}\mathcal{A}_N}{k^4} \left[\sum_\sigma \frac{N_\sigma}{N} c_\sigma^{(1,N)}\right] \cos(kL + \phi)$$ where $\mathcal{B}_N$ is a non-oscillatory shift due to wall-induced "half-cusps," and the oscillatory component (frequency $\sim kL$) encodes the long-range spin correlations of the many-body state (Aupetit-Diallo et al., 2023). The phase and amplitude of these oscillations are directly connected to the ground-state symmetry and spin coherence. Experimental observation of these structured HMTs permits direct nonlocal probing of quantum order and is robust against shot-to-shot atom number fluctuations in the bosonic ground state.

6. Theoretical Methods for Removing or Controlling HMTs

In effective field theories and low-momentum Hamiltonians, unphysical HMTs associated with model-dependent high-momentum form factors may contaminate numerical results but can be systematically integrated out. The SRG (similarity renormalization group) with a block-diagonal generator applied to a coarse grid in $p$-space decouples high-momentum ($p>\Lambda$) and low-momentum ($p<\Lambda$) components, yielding a Hamiltonian that is strictly scattering-equivalent to the original but manifestly devoid of HMTs: $$H_s = \begin{pmatrix} H^{PP}_s & 0 \ 0 & H^{QQ}_s \end{pmatrix}$$ with only the $P$ block relevant for physical observables at $p<\Lambda$. All phase shifts and observables dependent on the low-momentum block are strictly unchanged, and the computational cost is greatly reduced (Gómez-Rocha et al., 2020).

7. HMTs in High-Energy Collider and Plasma Physics

High-momentum tails also appear as signatures in collider and plasma settings:

• Collider EFT: In $pp\to\ell^+\ell^-$ Drell–Yan production, the high-$p_T$ tail of the dilepton invariant-mass spectrum encodes contact interactions, with dimension-six SMEFT operators generating new contributions in the tail region well above the electroweak scale. Precision fits to ATLAS LHC data place stringent limits ($v^2/\Lambda^2\sim10^{-4}$–$10^{-2}$) on the associated four-fermion operator coefficients (Greljo et al., 2017).
• Kappa Velocities in Plasmas: In non-Maxwellian driven gases and space-plasma contexts, stochastic acceleration and velocity-dependent friction yield stationary distributions with heavy power-law tails, $\mu_{\mathrm{st}}(v) \propto v^{-\alpha}$ with $\alpha$ set by microscopic parameters (Banerjee et al., 2019). These "kappa-distributions" are the velocity-space analog of HMTs.

---

The study of high-momentum tails thus provides a unifying probe of short-range, dynamical, and symmetry-driven correlations in quantum many-body physics, with implications for cold atom systems, nuclear structure, heavy-ion collisions, condensed matter, and high-energy phenomena. The universal $k^{-4}$ or related algebraic tails, fixed by two-body physics but with amplitude controlled by many-body state contacts (or their analogs), have become central diagnostic tools in the analysis and interpretation of strongly correlated systems.