Bjorken Flow: Boost-Invariant Expansion in QGP

• Boost-invariant Bjorken expansion is a theoretical framework that models longitudinal relativistic fluid dynamics in heavy-ion collisions using Milne coordinates.
• It applies symmetry principles such as boost invariance and transverse homogeneity to reduce the hydrodynamic equations to ordinary differential equations in proper time.
• This framework underpins studies of ideal, viscous, magnetized, and spin hydrodynamics, providing key insights into quark-gluon plasma evolution at RHIC and LHC.

Boost-invariant Bjorken expansion is a cornerstone theoretical framework for describing the longitudinal dynamics of relativistic fluids produced in ultra-relativistic heavy-ion collisions. Rooted in the symmetries of ultra-relativistic central nucleus-nucleus reactions, it prescribes a (1+1)D longitudinal expansion invariant under Lorentz boosts along the beam axis, leading to a highly constrained, analytically tractable hydrodynamic evolution. This framework forms the basis for modeling thermalization, nonequilibrium dynamics, and subsequent hydrodynamic evolution of the quark–gluon plasma (QGP) and other relativistic media.

1. Symmetry Principles and Coordinate System

Bjorken flow assumes strict longitudinal boost invariance and transverse homogeneity. The natural coordinates are Milne coordinates, defined by the proper time $\tau = \sqrt{t^2 - z^2}$ and spacetime rapidity $\eta = \tanh^{-1}(z/t)$, with transverse coordinates $x_\perp = (x, y)$. In these variables, the metric is $ds^2 = -d\tau^2 + dx_\perp^2 + \tau^2 d\eta^2$ (Rodgers et al., 2022, Grozdanov, 12 Oct 2025).

This symmetry imposes:

• The fluid four-velocity: $u^\mu = (1,0,0,0)$ in $(\tau, \eta, x, y)$ coordinates.
• All hydrodynamic fields (e.g., $T$, $p$, $\epsilon$) are functions of $\tau$ only; invariance under $\eta \to \eta + \text{const}$ reduces the dynamics to ODEs in $\tau$.

2. Ideal and Viscous Hydrodynamic Evolution

Conformal ideal flow (for a perfect fluid with equation of state $p = w\epsilon$, $w=1/3$ for QCD at high $T$): $$\frac{d\epsilon}{d\tau} + \frac{\epsilon + p}{\tau} = 0 \Rightarrow \epsilon(\tau) = \epsilon_0 \left( \frac{\tau_0}{\tau} \right)^{1+w}$$ For $w=1/3$, this yields $\epsilon \sim \tau^{-4/3}$ and $T \sim \tau^{-1/3}$ (Grozdanov, 12 Oct 2025, Singh, 2020). For a non-conformal EOS, the exponent $1 + c_s^2$ ($c_s^2 = dP/d\epsilon$) replaces $4/3$ (Rath et al., 2018, Rath et al., 13 Apr 2025).

Viscous first-order corrections: The leading dissipative contribution comes from the shear viscosity. For conformal fluids, bulk viscosity vanishes. The corrected ODE: $$\frac{d\epsilon}{d\tau} + \frac{4}{3} \frac{\epsilon}{\tau} = \frac{4\eta}{3\tau^2}$$ gives

$$\epsilon(\tau) = \epsilon_0\, \tau^{-4/3} \left(1 - \frac{H_0}{2}\tau^{-4/3} + \cdots \right)$$

where $H_0 \sim \eta/s$ (Grozdanov, 12 Oct 2025). The temperature scaling and energy density exhibit subleading corrections in $1/\tau$ (slower cooling). In the presence of finite angular velocity $\Omega$, transport coefficients (shear $\eta$, bulk $\zeta$, and $c_s^2$) acquire $\Omega$-dependence and generally increase, further accelerating cooling (Rath et al., 13 Apr 2025).

3. Generalized Bjorken Solutions: Beyond Conformality and Additional Degrees of Freedom

Superfluidity and spontaneous symmetry breaking: When the fluid supports a spontaneously broken $\mathrm{U}(1)$ (e.g., pion superfluid), the hydrodynamics includes a Goldstone phase $\varphi$ and a relative superfluid velocity $\xi^\mu = \Delta^{\mu\nu}\partial_\nu \varphi$. Exact, boost-invariant superfluid Bjorken solutions are found with $\xi^\eta = \xi_0/\tau^2$ and $\varphi = \xi_0\eta$, preserving boost invariance. Corrections arise at $\mathcal{O}(\tau^{-2})$ in energy density evolution, both from the superfluid sector and viscosity (Rodgers et al., 2022).

Spin hydrodynamics: For a perfect fluid with spin (e.g., spin-$1/2$ QCD matter), spin density $S^{\mu\nu}$ components governed by GLW or canonical formalism decay as $S(\tau)\sim \tau^{-1}$, with spin chemical potential $\mu_s \sim \tau^{-1/3}$, much like baryon number and temperature (Singh, 2020, Wang et al., 2021). Dissipative and relaxation-time effects further suppress spin observables at large $\tau$.

Magnetized fluids: A strong magnetic field ($B$) enters the EOS, raising $c_s^2$ above the conformal value due to paramagnetism, and strengthening the cooling: $\epsilon(\tau) \sim \tau^{-[1 + c_s^2(B)]}$, with $c_s^2(B) \sim 0.36$–$0.40$ for $eB\sim 15\,m_\pi^2$ (Rath et al., 2018). Rotation yields analogous consequences (Rath et al., 13 Apr 2025).

Nonconformal systems: For massive constituents, bulk viscous pressure $\Pi$ enters energy–momentum evolution. Early nonequilibrium attractors for the shear channel are destroyed by strong bulk–shear coupling; only the scaled longitudinal pressure $P_L/P$ displays a robust attractor at early times, with both shear and bulk inverse Reynolds numbers relaxing slowly (Chattopadhyay et al., 2021).

Nonextensive QGP: Taking Tsallis or other nonextensive distributions as the underlying kinetic ensemble modifies early-time evolution of the bulk viscous channel but has little effect on temperature or energy density evolutions in near-conformal flows (Alqahtani et al., 2022).

4. Attractor Phenomena and Non-equilibrium Dynamics

Kinetic theory and second-order viscous hydrodynamics with Bjorken symmetry feature a universal "hydrodynamic attractor" for energy density and pressure anisotropy observables. In particular, the RTA kinetic theory—when cast in dimensionless form using $w=\tau/\tau_R$—features a unique attractor $g_{\mathrm{att}}(w)$ that interpolates between early-time free-streaming ($g_+ \approx -0.929$) and late-time Navier-Stokes ($g_* = -4/3$). The attractor solution is exactly given by confluent hypergeometric (Kummer) functions: $$g_{\rm att}(w) = g_+ - w + w\frac{a}{b} \frac{M(1+a, 1+b, w)}{M(a,b,w)}$$ with $a,b$ fixed by the geometry (Blaizot et al., 2020). This solution governs the pre-equilibrium to hydrodynamic transition; all initial conditions rapidly approach the attractor for $w \gtrsim 1$. This feature is robust across different hydrodynamic models, kinetic theory, and even strongly coupled holographic theories.

Bulk and shear viscous attractors are lost in non-conformal (massive) systems, but the universal scaling of longitudinal pressure is preserved (Chattopadhyay et al., 2021). Anisotropic hydrodynamics captures far-from-equilibrium evolution more accurately than standard viscous approaches, especially for small systems and at early times (Peng et al., 4 Sep 2025).

5. Microscopic Realizations and Numerical Approaches

The lattice Boltzmann method (LBM) provides a discretized, polynomial-expanded solution to the Anderson–Witting collision kernel in Bjorken flow. Anisotropic Romatschke–Strickland initial distributions allow for quantitative studies as a function of $\eta/s$ and initial pressure anisotropy, validating hydrodynamic attractor behavior and the breakdown of hydrodynamics at large Knudsen number (Ambrus et al., 2018).

Moment expansions of the single-particle distribution in the RTA kinetic theory (using Laguerre and Legendre polynomials) diverge factorially due to large populations of higher moments at early times, reflecting persistent nonequilibrium modes. Borel–Padé resummation techniques, as demonstrated in recent work (Brito et al., 2024), recover a physically meaningful distribution function far from equilibrium, with nonhydrodynamic multipoles decaying only slowly in expanding systems.

6. Generalizations, Geometric Constructions, and Extensions

The geometric origin of Bjorken flow can be traced to the Weyl mapping of a static, maximally symmetric fluid configuration on $dS_3 \times \mathbb{R}$ with $\mathrm{ISO}(2)$-invariant slicing to Minkowski space. This construction leads naturally to a family of three maximally symmetric, boost-invariant flows in hydrodynamics: flat (Bjorken), spherical (Gubser), and hyperbolic slicings. Each is characterized by distinct late-time dynamics and transverse structure, but only the Bjorken solution is globally translation invariant transversely (Grozdanov, 12 Oct 2025).

7. Physical Consequences and Applications

The boost-invariant Bjorken expansion paradigm remains foundational for phenomenological modeling of QGP evolution at RHIC and LHC. It enables analytic and semi-analytic control over thermalization, hydrodynamization, vorticity, spin polarization, and magnetization signatures. It provides a universal testing ground for resummation techniques in nonequilibrium hydrodynamics and kinetic theory, guides the development of generalized frameworks (e.g., superfluid, magnetized, spinning, nonextensive, and anisotropic hydrodynamics), and underpins our understanding of the limits of hydrodynamics in the high-Knudsen and small-system regimes.

Boost invariance at the level of macroscopic evolution is robust as long as the underlying microscopic physics respects the symmetry (e.g., zero pion mass for superfluids (Rodgers et al., 2022), negligible vorticity/magnetic field gradients (Rath et al., 2018, Rath et al., 13 Apr 2025)). Deviations (e.g., explicit mass terms for Goldstone bosons, transverse gradients, finite rapidity acceptance) cause symmetry breaking at large rapidity or for large deviations from centrality.

In summary, the boost-invariant Bjorken expansion, both in its ideal and dissipative generalizations, continues to provide a central analytic backbone for the study of ultrarelativistic heavy-ion collisions, the emergence of collective behaviors, and the rapidly developing theory of nonequilibrium relativistic fluids.