Universal Phase-Space Scaling

• Universal phase-space scaling is the emergence of simple power-law relationships between phase-space quantities and control parameters across diverse systems such as Fermi gases, dark matter halos, and chaotic dynamics.
• It provides a unified framework that collapses complex thermodynamic and dynamical behavior onto dimensionless curves, enabling predictive analyses in quantum criticality, astrophysical simulations, and extended chaos.
• The methodology integrates dimensionless parametrization, numerical simulations, and mappings to universality classes like KPZ, offering insights into self-similarity, critical behavior, and relaxation processes.

Universal phase-space scaling describes the emergence of simple, robust power-law relationships between phase-space quantities—such as the distribution function, energy distribution, or variance of dynamical observables—and system parameters or geometric variables. This universality appears in diverse contexts: the phase diagrams and thermodynamics of imbalanced Fermi superfluids near a quantum multicritical point, the phase-space structure and energy distribution of dark matter halos, and the scaling of Lyapunov-exponent fluctuations in extended chaotic systems. Across these domains, universal phase-space scaling reflects underlying scale invariance, critical behavior, and self-similar relaxation processes.

1. Universal Scaling in Fermi Gases: Multicriticality and Dimensional Analysis

In resonantly interacting two-component Fermi gases, universal phase-space scaling arises near the quantum multicritical point at zero chemical potential $\mu$, effective magnetic field $h$, and inverse scattering length $1/a$ (Frank et al., 2018). In the unitary limit ($1/a=0$), scale invariance restricts the relevant dimensionless parameters to $x = T/\mu$ (relative temperature) and $y = h/\mu$ (relative Zeeman field). All bulk thermodynamics—including pressure $P$, density $n$, spin-imbalance density $s$, and entropy $S$—are encoded in a small set of universal scaling functions (e.g., $h$0) of these ratios. For instance,

$h$1

with other observables obtained by differentiation.

The universality is especially manifest at the multicritical point, which acts as a nontrivial fixed point with three relevant directions ($h$2, $h$3, $h$4). Thermodynamic behavior throughout the phase diagram—transition lines, critical behavior, and even the boundaries of inhomogeneous superfluidity (FFLO phase)—can be collapsed onto universal curves when plotted against the properly rescaled variables. Away from unitarity, an additional scaling variable $h$5 (either $h$6 or $h$7) enters, but the principle of dimensionless parametrization and function collapse remains. This scaling framework provides both predictive power and unification across the BCS, BEC, and unitary regimes.

2. Phase-Space Scaling in Dark Matter Halos: Distribution Functions and Energy Laws

Numerical simulations of collisionless dark matter halos consistently reveal universal phase-space scaling in their isotropic distribution functions $h$8 and differential energy distributions $h$9 (Gross et al., 2024). For halos well fit by a Navarro-Frenk-White (NFW) profile, Gross et al. find that

$1/a$0

across more than a decade in radius, where $1/a$1 is the radius where the gravitational potential equals the energy $1/a$2. The coefficients are fixed by the gravitational constant $1/a$3 and the NFW scale parameters $1/a$4, e.g.,

$1/a$5

This scaling mirrors the exact results for the singular isothermal sphere (SIS), where $1/a$6 and $1/a$7. The empirical exponents very nearly match those predicted by self-similar secondary-infall models and reflect the fact that, at $1/a$8, the NFW halo's density profile approaches the SIS’s $1/a$9 scaling. This universality holds across a wide halo mass range, with deviations only in the innermost (cusp) and outermost (finite-radius) regions.

3. Universal Scaling in Spatially Extended Chaos: Lyapunov Exponent Fluctuations

Universal scaling extends to the fluctuation statistics of dynamical invariants in chaotic extended systems (Pazó et al., 2013). Specifically, Paz, López, and Politi demonstrate that the diffusion coefficient $1/a=0$0 of finite-time Lyapunov exponents (FTLEs) scales non-trivially with system size $1/a=0$1:

$1/a=0$2

The wandering exponent $1/a=0$3 is not model-dependent, but instead set by roughening exponents $1/a=0$4 of the "Lyapunov-surface" through $1/a=0$5. For the maximal FTLE (first Lyapunov vector) in one-dimensional dissipative chaos, universality arises from a mapping to the Kardar-Parisi-Zhang (KPZ) universality class, giving $1/a=0$6, $1/a=0$7, and thus $1/a=0$8. Extensive numerical work confirms $1/a=0$9 holds for diverse systems—including chains of Hénon maps, Lorenz-96, and delay-differential models.

For the bulk of the spectrum (higher FTLEs), the observed exponents differ ($x = T/\mu$0–$x = T/\mu$1 in 1D), with no known stochastic PDE capturing their scaling—a signature of a potentially distinct and still unidentified universality class. This scaling uncovers deeper statistical-geometric properties of phase-space trajectories in high-dimensional, stochastic, or chaotic dynamics.

4. Methodological Frameworks and Derivations

In unequal Fermi gases, the Luttinger-Ward (LW) formalism, truncated at the particle-particle ladder, provides a self-consistent statistical framework for deriving the universal scaling functions in the thermodynamic potentials. Starting from the grand-potential functional $x = T/\mu$2, stationarity with respect to the Green’s functions $x = T/\mu$3 leads to self-consistent Dyson equations. The pressure, density, and other thermodynamic quantities then directly yield the scaling functions $x = T/\mu$4, $x = T/\mu$5, $x = T/\mu$6 via established differentiation (e.g., $x = T/\mu$7).

For self-gravitating halos, direct measurement of $x = T/\mu$8 and $x = T/\mu$9 is performed on simulated particle data binned by $y = h/\mu$0. The empirical power-law forms are validated against known analytical results (SIS, Eddington inversion) and interpreted through the lens of secondary-infall and self-similar relaxation. Full proportionality coefficients are presented in terms of $y = h/\mu$1, $y = h/\mu$2, and $y = h/\mu$3.

For space-time chaos, analysis of Lyapunov exponent statistics combines both the dynamic-scaling theory for surface roughness and extensive simulation. The mapping to KPZ universality leverages the field interpretation $y = h/\mu$4, connecting the fluctuation statistics of dynamical invariants to universality classes of stochastic growth models.

5. Physical and Theoretical Significance

Universal phase-space scaling encodes the system's response to changes in control parameters, independent of microscopic details, so long as the dynamical or statistical regime is properly captured by the scaling variables. In Fermi gases, these scalings unlock unified understanding of the rich phase diagrams and critical phenomena, and connect BCS, unitary, and BEC physics across different detunings and polarizations. In gravitational systems, they condense complex structure-formation histories into succinct, predictive relations that underpin both simulation analysis and semi-analytical modeling.

Similar scaling relations in chaotic systems illuminate how fluctuations, transport, and sensitivity to perturbations are governed by system size and effective field theories (e.g., KPZ), rather than specific dynamical rules. This suggests a deep entanglement between statistical mechanics, dynamical systems, and the geometrization of phase space, encompassing both equilibrium and non-equilibrium contexts.

6. Domains of Applicability, Limitations, and Open Questions

Universal phase-space scaling, while robust, is constrained by symmetry, statistical regime, and the form of interactions:

• Fermi gases: Scaling holds near the unitary point and multicriticality; deviations occur away from the zero-range (scale invariant) limit or when higher-order corrections become relevant.
• Dark matter halos: The empirical scalings are valid for $y = h/\mu$5 in $y = h/\mu$6 for $y = h/\mu$7, $y = h/\mu$8 for $y = h/\mu$9, with predictable deviations near the cusp and at large radii due to changes in the density profile and finite counting statistics.
• Spatiotemporal chaos: Universality for Lyapunov fluctuations is established for the first exponent (KPZ class) and observed for the bulk, but the governing equations for the latter remain unidentified. The persistence of positive roughness exponents ($P$0) and the uniqueness of the bulk universality class are open questions (Pazó et al., 2013).

A plausible implication is that different physical systems—across quantum criticality, self-gravitating relaxation, and chaotic dynamics—may exhibit universal phase-space scaling due to underlying self-similarity, critical phenomena, or universality-class structure. The search for minimal field-theoretic or stochastic models unifying these behaviors (especially for bulk Lyapunov exponents) remains an active area of research.

7. Implications and Analytical Utility

The existence of universal phase-space scaling enables compact parametrizations of complex systems:

• Imbalanced Fermi gases: Rapid, accurate prediction of thermodynamics and phase boundaries across interaction regimes via scaling functions and dimensionless variables (Frank et al., 2018).
• Dark matter halos: Single-line formulas for $P$1 and $P$2, given NFW parameters, bypass iterative Eddington inversions and facilitate both simulation diagnostics and analytical modeling (Gross et al., 2024).
• Chaos theory: Scaling exponents extracted from Lyapunov fluctuations furnish a bridge between dynamical systems theory and statistical field theory, potentially guiding new models for spatiotemporal chaos (Pazó et al., 2013).

More broadly, universal phase-space scaling exemplifies how symmetry, dimensional analysis, and critical-point theory can underpin predictive frameworks throughout many-body physics, astrophysics, and nonlinear dynamics.