--- title: Pre-equilibrium Energy Loss url: https://www.emergentmind.com/topics/pre-equilibrium-energy-loss type: topic --- # Pre-equilibrium Energy Loss Pre-equilibrium energy loss denotes energy degradation, dissipation, or redistribution that occurs before a many-body system reaches the equilibrium or hydrodynamic regime. In relativistic heavy-ion physics, the term covers both early-time quenching of hard probes and pre-hydrodynamic transport of the bulk energy-momentum tensor; in low-energy nuclear reactions it refers to particle emission and collective radiation before compound-nucleus equilibration; in pre-supernova matter it refers to neutrino and antineutrino energy loss from thermally excited nuclei before collapse [2005.03330]. Across these settings, the defining feature is temporal ordering: the relevant energy-transfer mechanism acts while the medium is still far from equilibrium, anisotropic, or only partially thermalized. ## 1. Conceptual scope and regimes In the heavy-ion literature, the pre-equilibrium interval is the stage before hydrodynamics is assumed applicable. One heavy-flavor study takes this as \(0<\tau<\tau_0\) with \(\tau_0=0.6~\mathrm{fm}\) in the default setup, and explicitly asks when heavy-quark–medium interaction begins and what temperature profile should be used before \(\tau_0\) [2005.03330]. Weak-coupling kinetic descriptions instead formulate the same interval as evolution from a saturated initial state at \(\tau_0\sim Q_s^{-1}\) to a later \(\tau_{\rm hydro}\sim 1~\mathrm{fm}/c\) when viscous hydrodynamics becomes applicable [1610.09912]. In Bayesian jet-quenching phenomenology, hydrodynamics starts at \(\tau_{\rm hyd}=0.4~\mathrm{fm}\), and pre-equilibrium energy loss is represented by allowing quenching to turn on at an earlier onset time \(\tau_{\min}\) [2509.19430]. | Domain | Carrier of energy loss or redistribution | Representative observables | |---|---|---| | Relativistic heavy-ion hard probes | Heavy quarks, jets, partons | \(R_{AA}\), \(v_2\), jet \(R_{AA}\), jet \(v_2\) | | Relativistic heavy-ion bulk medium | \(T^{\mu\nu}\), entropy, preflow | Green functions, hydrodynamic initial conditions | | Low-energy heavy-ion reactions | Dynamical dipole radiation, nucleon emission | \(\gamma\)-ray spectrum, emitted \(N/Z\), angular distributions | | Pre-supernova stellar matter | Weak processes on hot nuclei | Neutrino spectra, antineutrino spectra, energy-loss rates | The phrase does not denote a single microscopic mechanism. In some works it means collisional or radiative parton energy loss before equilibration [1110.4846]; in others it means bulk redistribution of energy and momentum before hydrodynamic matching rather than loss from the system [1805.00961]. A frequent source of confusion is the use of free streaming as a proxy for pre-equilibrium evolution. One isobar-collision study is explicit that free streaming is an *ad hoc extreme scenario*, that the pre-equilibrium stage is poorly understood, and that free streaming with velocity \(c\) should be interpreted as an upper bound on smearing of the initial geometry [2305.03703]. ## 2. Early-time quenching of hard probes in heavy-ion collisions For heavy flavor, pre-equilibrium energy loss has been analyzed in a modified Langevin framework coupled to hydrodynamics. The heavy-quark momentum obeys \[ \frac{d\vec p}{dt}=-\eta_D(p)\,\vec p+\vec\xi+\vec f_g, \] with drag, white-noise diffusion, \[ \langle \xi^i(t)\xi^j(t')\rangle=\kappa\,\delta^{ij}\delta(t-t'), \] and recoil from medium-induced gluon radiation. The fluctuation-dissipation relation gives \[ \eta_D(p)=\frac{\kappa}{2TE}, \qquad D_s=\frac{T}{M\eta_D(0)}=\frac{2T^2}{\kappa}, \] so the interaction strength is characterized by \(D_s(2\pi T)\). Radiation is incorporated probabilistically using a higher-twist gluon-emission spectrum with dead-cone suppression and a low-energy cutoff \(\omega_0=\pi T\) below which radiation is switched off [2005.03330]. The pre-equilibrium stage enters in two ways: by changing the time interval over which the Langevin evolution acts, and by changing the temperature history that controls drag, diffusion, and radiation rates. Four pre-equilibrium temperature profiles were tested for \(0<\tau<\tau_0\): free streaming, \(T(\tau)=0\); linear rise, \(T(\tau)=T(\tau_0)(\tau/\tau_0)\); constant temperature, \(T(\tau)=T(\tau_0)\); and Bjorken-like evolution, \(T(\tau)=T(\tau_0)(\tau_0/\tau)^{1/3}\). An important implementation point is that the same modified Langevin equation is used before \(\tau_0\) whenever a temperature profile is supplied [2005.03330]. The phenomenological consequence is that suppression and flow do not respond identically. After retuning \(D_s\) so that different scenarios yield similar \(R_{AA}\), delaying the interaction onset from \(0.6\) to \(1.2~\mathrm{fm}\) requires about \(35\%\) smaller \(D_s\), i.e. stronger effective coupling, to reproduce the same \(R_{AA}\). Under the same retuning, low-\(p_T\) \(D\)-meson \(v_2\) increases by \(8\%\) at RHIC and \(24\%\) at the LHC. For pre-equilibrium temperature profiles, the free-streaming assumption yields about \(39\%\) larger low-\(p_T\) \(v_2\) at RHIC and \(19\%\) larger at the LHC than constant-temperature and Bjorken-like evolution. The paper’s broader result is that different modeling of initial conditions, pre-equilibrium evolution, and in-medium interaction can each induce about \(10\%-40\%\) uncertainties in low-\(p_T\) heavy-flavor suppression and flow [2005.03330]. A conceptually distinct early-time mechanism attributes azimuthal anisotropy itself to pre-equilibrium energy loss in the strong color field of overlapping strings or string clusters. In that picture, a parton traverses an anisotropic color field before the medium equilibrates, radiates, and loses energy by an amount that depends on the path length \(l(\phi)\). The resulting elliptic-flow scaling law is \[ \frac{v_2(p_T)}{\epsilon\, Q_s^A\, L}=f(\tau), \qquad \tau=\left(\frac{p_T}{Q_s^A}\right)^2, \] with \(f(\tau)\propto \tau^{1/3}\) as the central theoretical expectation for \(p_T\lesssim Q_s^A\). The fitted scaling curve is \(f(\tau)\sim a\tau^b\) with \(a=0.1264\pm0.0076\) and \(b=0.404\pm0.025\), and the model interprets the near-universality across centralities, species, and RHIC/LHC energies as evidence that early-time quenching can generate much of the observed \(v_2\) scaling [1609.03927]. More recent jet phenomenology introduces pre-equilibrium loss through the hydrodynamic attractor. In that framework, the quenched jet spectrum is a convolution of an nPDF-modified vacuum cross section with an energy-loss probability distribution, quenching is applied only for \(\tau\ge \tau_{\min}\), and the same energy-loss machinery is used in the pre-equilibrium stage with an attractor-extrapolated temperature and preflow field. A Bayesian analysis using jet \(R_{AA}\) and jet \(v_2\) finds that jet \(R_{AA}\) alone leaves \(\tau_{\min}\) poorly constrained because earlier onset can be compensated by smaller coupling, but including jet \(v_2\) constrains \(\tau_{\min}\approx 0.24~\mathrm{fm}\), well before \(\tau_{\rm hyd}=0.4~\mathrm{fm}\). The same study reports good simultaneous description of jet \(R_{AA}\), jet \(v_2\), RHIC jet suppression, and charged-hadron \(R_{AA}\), with the main tension in low-\(p_T\) hadron \(v_2\), which is underpredicted [2509.19430]. ## 3. Initial-value formulations and unstable plasma dynamics A separate line of work formulates pre-equilibrium parton energy loss as an initial-value problem in a weakly coupled, anisotropic, unstable quark-gluon plasma. The energetic parton is treated as a classical colored particle obeying the Wong equations, and the instantaneous energy loss is \[ \frac{dE(t)}{dt}=gQ^a\,{\bf E}_a(t,{\bf r}(t))\cdot{\bf v}. \] Because the medium is not stationary, the analysis uses a one-sided Fourier transform, \[ f(\omega,{\bf k})=\int_0^\infty dt\int d^3r\, e^{i(\omega t-{\bf k}\cdot{\bf r})}f(t,{\bf r}), \] so that the energy transfer retains explicit dependence on initial chromodynamic fields and on unstable collective modes [1110.4846]. In equilibrium or any stable plasma, all collective poles lie in the lower half of the complex \(\omega\)-plane, their contributions decay, and the familiar steady collisional energy-loss formula is recovered. In unstable plasmas, however, poles of the retarded propagator enter the upper half-plane, \[ \omega=\omega_r+i\gamma,\qquad \gamma>0, \] producing factors \(e^{\gamma t}\) and hence explicitly time-dependent energy transfer. This is the origin of exponentially growing contributions to pre-equilibrium loss in the weak-coupling unstable-plasma literature [1110.4846]. Two benchmark anisotropic systems have been studied in detail. For a two-stream plasma with \[ f({\bf p})=(2\pi)^3 n\Big[\delta^{(3)}({\bf p}-{\bf q})+\delta^{(3)}({\bf p}+{\bf q})\Big], \] the instability is longitudinal and chromoelectric. The computed energy loss per unit length oscillates strongly in time, has an amplitude that grows with time, and depends strongly on the angle between the parton velocity and the stream direction [1201.1486]. For an extremely prolate plasma, the momentum distribution is infinitely elongated along one direction, unstable Weibel-like modes appear, and the energy loss is strongly time dependent, strongly directional, and can become much larger than the equilibrium loss; in the numerical setup discussed in the paper, the corresponding equilibrium benchmark is about \(0.18\) for \(k_{\max}=5\mu\), \(g=1\), \(C_R=N_c=3\) [1301.4563]. The ultraviolet structure is also nontrivial. For the extremely prolate case, the soft contribution is logarithmically ultraviolet divergent, and a dedicated analysis argues that a good approximation is obtained by cutting the divergence off with the parton energy, \(k_{\max}\sim E\) [1312.5913]. Another extension shows that the sign of the energy transfer can depend on initial conditions. With uncorrelated initial chromodynamic fields, the parton typically loses energy and the magnitude is comparable to equilibrium; with maximally correlated initial fields induced by the parton, the parton can either gain or lose energy, the transfer grows exponentially in time, and the magnitude can far exceed equilibrium loss. The effect is maximal for \({\bf v}\parallel{\bf n}\) in extremely prolate plasmas and for \({\bf v}\perp{\bf n}\) in extremely oblate plasmas [1506.09082]. Strong-coupling holography yields a different nonequilibrium pattern. In collisions of gravitational shock waves dual to colliding sheets of energy, the drag force on a heavy quark moving through far-from-equilibrium matter is not dramatically enhanced relative to equilibrium estimates with a local effective temperature or pressure, but it turns on only after a delay. For zero rapidity the force rises after the collision with a delay roughly of order \(1/(\pi T_{\rm hydro})\) at low velocity and later becomes semi-quantitatively describable by equilibrium drag based on the local hydrodynamic temperature. At larger rapidity, gradients in the fluid velocity produce qualitatively new effects: the force required to hold the quark can point in the same direction as the quark velocity, and a substantial component perpendicular to the velocity can appear even in the local fluid rest frame [1306.0564]. This establishes an important contrast: unstable weak-coupling plasmas can generate exponential enhancement, whereas a strongly coupled far-from-equilibrium collision does not automatically imply extra drag, but does imply delay and gradient corrections. ## 4. Bulk pre-equilibrium transport, preflow, and medium background In weak-coupling kinetic theory, pre-equilibrium “energy loss” often refers not to hard-probe quenching but to transport and redistribution of the bulk energy-momentum tensor. The central object is the Boltzmann equation with leading-order QCD processes, \[ \partial_\tau f + \frac{\mathbf p}{|p|}\cdot \nabla_x f - \frac{p^z}{\tau}\,\partial_{p^z} f = -\mathcal C_{2\leftrightarrow 2}[f]-\mathcal C_{1\leftrightarrow 2}[f], \] where the collision kernel includes elastic \(2\leftrightarrow 2\) scattering and inelastic \(1\leftrightarrow 2\) splittings with screening and Landau-Pomeranchuk-Migdal suppression. The purpose is to propagate a far-from-equilibrium, gluon-dominated, boost-invariant system from \(\tau_0\sim Q_s^{-1}\) to a time when viscous hydrodynamics becomes valid [1605.04287]. The practical output is a Green-function map from initial perturbations to hydrodynamic initial data. In one formulation, the evolved tensor at \(\tau_{\rm hydro}\) is decomposed into a local homogeneous background plus linear response to initial energy and momentum perturbations, \[ T^{\mu\nu}(\tau_{\rm hydro},x) = T^{\mu\nu}_x(\tau_{\rm hydro}) + \frac{T^{\tau \tau}_x(\tau_{\rm hydro})}{T^{\tau\tau}_x(\tau_0)} \int d^2x'\, G^{\mu\nu}_{\alpha\beta}(x,x',\tau_{\rm hydro},\tau_0)\, \delta T^{\alpha\beta}_x(\tau_0,x') . \] This event-by-event, causal construction propagates local energy density and momentum perturbations only within the causal neighborhood \(|x'-x|