Pith. sign in

REVIEW 4 major objections 5 minor 9 references

Connecting Pre-Thermal and Hydrodynamizing Attractors With Adiabatic Hydrodynamization

T0 review · 4 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The hydrodynamizing attractor of an expanding gluon plasma is the late-time ground state of the same effective Hamiltonian whose early-time ground-state band is the pre-thermal attractor.

desk verdict A plainly written proceedings paper that gives a compact numerical illustration of a plausible unified attractor picture, but the central 'no mixing from higher bands' claim is not yet shown to be independent of the chosen rescaling and basis. read the letter →

arxiv 2504.13332 v1 pith:F6SFIP46 submitted 2025-04-17 hep-ph

Krishna Rajagopal , Bruno Scheihing-Hitschfeld , Rachel Steinhorst This is my paper · ORCID
classification hep-ph
keywords adiabatichydrodynamizationattractorquark-gluonplasmakinetictheorybottom-upthermalizationeffectiveHamiltonianprescalinglongitudinalexpansion
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that the early far-from-equilibrium attractor of an expanding gluon plasma and the later hydrodynamization attractor are not separate phenomena. Using the adiabatic hydrodynamization (AH) framework, it shows that the pre-thermal attractor is a band of nearly degenerate low-energy modes, and that hydrodynamization happens when all but one of these modes rise in effective energy, leaving a unique ground state. The surviving mode is the hydrodynamic mode, built exclusively from the lowest band with no mixing from higher-energy modes. The claim matters because it turns a sequence of empirically observed attractors into a single intuitive mechanism of staged memory loss.

What carries the argument

The central object is the effective Hamiltonian $H_{\mathrm{eff}}$ defined by $H_{\mathrm{eff}} w = -\partial_y w$, where $y \equiv \ln(\tau/\tau_I)$ and the distribution function is rescaled as $f(p,u,\tau) = A(\tau)\, w(p/D(\tau), u, r(\tau), \tau)$ with $u \equiv p_z/p$. Its instantaneous eigenstates have effective energies whose excited modes decay roughly as $e^{-\int \epsilon_i(\tau')\, d\tau'}$ when $H_{\mathrm{eff}}$ evolves slowly. The argument is carried by the specific choices of the rescalings: $A(\tau)$ from number conservation, $D(\tau)$ from a chosen relaxation equation toward the effective temperature, and $r(\tau)$ fixed by requiring the first basis state to satisfy the evolution equation; together with a basis that interpolates between early-time Gaussian and late-time exponential behavior, these choices produce the band structure and its late-time splitting into a unique ground state.

What would settle it

Repeat the calculation with more than the twelve basis states, or with a different time evolution for $D(\tau)$ such as following the full effective temperature directly, and measure the overlap of the late-time ground state with higher bands. If the unique ground state mixes with excited bands, or if the lowest band fails to split into a single surviving mode under a valid alternative choice, the claimed continuous connection is an artifact of the basis rather than a property of the kinetic theory.

Watch

Extended reading notes

Core claim

For a boost-invariant, longitudinally expanding gluon gas described by a simplified QCD kinetic theory with small-angle elastic scattering, the paper demonstrates that the pre-thermal attractor and the hydrodynamizing attractor are continuously connected through the spectrum of a time-dependent effective Hamiltonian. At early times the low-energy eigenstates form a nearly degenerate band that acts as an attractor surface; at later times all but one of these energies lift off, and the lone remaining ground state evolves adiabatically into the hydrodynamic mode. The hydrodynamic mode therefore arises from within the pre-thermal band, with no contribution from higher bands, so the system's memory loss proceeds in two stages: first the excited longitudinal modes are forgotten, then the remaining band collapses to a single hydrodynamic state.

Load-bearing premise

The load-bearing premise is that the chosen coordinate rescalings and basis, including the hand-picked evolution equation for $D(\tau)$ and the condition fixing $r(\tau)$, faithfully represent the dynamics rather than imposing the observed band structure and absence of mixing.

Editorial extensions

If this is right

  • If the central claim is correct, the pre-thermal and hydrodynamizing attractors are the same object at different epochs, so no new mode needs to be invoked for hydrodynamization.
  • Memory loss in the expanding gluon gas occurs in two well-defined stages, with the later stage consisting solely of the lowest band collapsing to one surviving mode.
  • The AH description provides a mechanistic explanation for why attractor behavior occurs: the system settles into the instantaneous ground state of a slowly evolving effective Hamiltonian.
  • The time-dependent basis and scalar choices used here track the full distribution through both attractors, suggesting that a single reduced description can cover the entire pre-hydrodynamic evolution.
  • Including number-non-conserving processes should hasten hydrodynamization while preserving the same band-to-ground-state structure, making the scenario applicable to realistic QCD kinetic theory.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Extending the paper's logic, the same band-splitting mechanism should be visible in full QCD effective kinetic theory once inelastic scattering is added, and one could test this by projecting exact kinetic-theory solutions onto a time-dependent basis and tracking the lowest-band overlaps.
  • The staged memory loss suggests a selection principle: among all possible coordinate rescalings, the physically realized attractor is the one that makes the effective Hamiltonian adiabatic and gapped; this could be used to predict new attractors in other expanding systems.
  • If the hydrodynamic mode carries no admixture from higher bands, then all surviving initial-state information before hydrodynamization is encoded in the lowest band, which may offer a practical reduced description for initial-state modeling in heavy-ion collisions.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper applies the Adiabatic Hydrodynamization (AH) framework to a simplified kinetic theory of a longitudinally expanding gluon gas, using the small-angle elastic collision kernel and no inelastic processes, to argue that the pre-thermal (BMSS/dilute) attractor and the hydrodynamizing attractor are continuously connected. It defines a time-dependent effective Hamiltonian Heff and shows numerically, for two couplings, that the early-time attractor is a band of near-degenerate low-energy modes and the late-time hydrodynamic mode emerges when all but one of the energies in that band lift off, leaving a unique ground state with no mixing from higher bands.

Significance. If correct, the AH framework offers a concrete physical mechanism for the two-stage memory loss in pre-hydrodynamic QGP and explains why the hydrodynamic mode inherits the pre-thermal attractor sector. The paper is largely a proceedings summary of Ref. [3], but it states a clear, falsifiable thesis. Its strengths are the explicit definition of a simplified collision kernel, the self-consistency check between the scaling exponents extracted from f and from the first basis state, and the comparison with the analytic prediction of Ref. [7]. The claim that the hydrodynamic mode is built from the lowest dilute band is a testable structure prediction for any kinetic theory.

major comments (4)
  1. [Sec. 2, Eqs. (2)-(4), (7)-(9)] The eigenvalue spectrum of Heff is only defined after choosing the time-dependent rescalings A(τ), D(τ), and r(τ), and the resulting spectrum is not invariant under changes of these rescalings. The choice ∂_yD/D = 10(1 - D/D_pE) is presented without derivation or robustness study, and r(τ) is fixed by a Galerkin condition on the first basis state. Since the central claim of Sec. 3 ('no mixing from higher-energy modes') is read directly from the spectrum of this Heff, the paper needs to demonstrate that the band structure and the late-time unique ground state survive other admissible rescalings (e.g., D fixed by a moment of f, or a different relaxation rate) and that the conclusion is not an artifact of the chosen frame.
  2. [Sec. 2, condition defining r(τ)] The condition ∫ p u² (∂_y + Heff) ψ10^(R) = 0 is a Galerkin projection that forces the lowest basis state to be an approximate solution of the projected evolution. The subsequent observation that f is well captured by the first basis state (agreement of the scaling exponents in Figs. 1-2) is therefore a self-consistency check, but it does not establish that the physical solution has negligible overlap with higher bands in a frame-independent sense. The authors should report the overlaps of f onto the excited basis states or otherwise show that the 'no mixing' conclusion is not imposed by the basis choice.
  3. [Sec. 3, Figs. 1 and 2] All numerical results are obtained with a 12-state basis, and no convergence test with respect to the basis size is reported; the text refers to Ref. [3] for numerical details. Because the claim that exactly one mode survives to form the hydrodynamic ground state is a spectral degeneracy/lifting statement, the 12-state truncation could in principle alter the late-time gap structure. The manuscript should include a truncation check (at least for the late-time spectrum) or state that such a check was performed in Ref. [3] and summarize its outcome.
  4. [Sec. 3, gs = 10^-3 case] The paper states that the analytical prediction of Eq. (13) 'clearly describes' the effective energies in the dilute regime, but no quantitative comparison is shown. The claim that the evolution is not adiabatic because the gap is 'vanishingly small' and Heff is time-dependent requires an estimate of the adiabaticity parameter (the ratio of the gap to the rate of change of Heff). Without such a comparison, the interpretation of the dilute regime as an 'attractor surface' rather than a non-adiabatic blob is not fully supported. This is a load-bearing point for the claim of continuous connection, so a quantitative measure of adiabaticity is needed.
minor comments (5)
  1. [Sec. 2, Eq. (8)] The variable χ appears in Heff and in the basis of Eq. (9) but is never defined; it should be stated that χ = p/D.
  2. [Sec. 2, Eq. (6)] The text says the explicit Bose enhancement factor is omitted, yet Ia[f] contains (1+f); the distinction between the Bose factor in the integral and the Boltzmann simplification should be clarified.
  3. [Sec. 3, Eq. (13)] The indices n and m are not mapped onto the basis states used in the numerical calculation; a brief explanation of how Eq. (13) is compared with the left panels of Figs. 1-2 would improve readability.
  4. [Outlook] The phrase 'collsion kernel' should be 'collision kernel'.
  5. [General] The paper is a proceedings contribution and refers to Ref. [3] for the numerical implementation; as a standalone manuscript it would benefit from a short summary of the spectral projection method and the initial condition.

Circularity Check

2 steps flagged · score 4.0 of 10

The attractor-as-ground-state identification is definitional, and the no-mixing conclusion is drawn in a hand-fitted rescaling frame without a test of frame independence.

  1. self definitional [Sec. 1, Eqs. (2)-(4) and the text following them.]
    "This inspires the definition of an effective Hamiltonian operator Heff satisfying Heff w =−∂yw (2) ... This is natural in the sense that for A, B, C satisfying ∂yA A =α, ∂yB B =−β, ∂yC C =−γ, (4), when f is prescaling, we will have∂yw = 0. That is, ws is an instantaneous eigenstate of Heff with “effective energy” 0."

    By Eq. (2), Heff is defined so that Heff w = -∂_y w for the rescaled distribution w. By Eq. (4), a prescaling solution has ∂_y w = 0. Therefore the statement that a prescaling attractor is a zero-energy eigenstate of Heff is not a discovered property of the dynamics but a direct consequence of the definitions: Heff is constructed to map w to -∂_y w, and prescaling is defined as ∂_y w = 0. The subsequent adiabatic language (“the attractor is the ground state”) thus partly restates the prescaling condition in new terminology rather than deriving it from the Boltzmann equation.

  2. fitted input called prediction [Sec. 2 (choice of r(y)) and Sec. 3 (conclusion of no mixing).]
    "Finally, we use our choice of r(y) to attempt to maximize the extent to which f is described by ψ(R)10, in particular by setting R p u2(∂y + Heff)ψ(R)10 = 0 and solving for∂yr. ... All but one of the effective energies which form the low-energy dilute “band” lift off at late times y ≳ 15 to reveal a unique ground state, the hydrodynamizing attractor."

    The function r(τ) is not predicted; it is fit by a Galerkin condition that forces the first basis state ψ10 to be an approximate solution of the projected evolution equation. The paper then presents the agreement of the scaling exponents (β1, γ1, α1 with β, γ, α) as evidence that f is “well-captured” by this basis, and reads off the central no-mixing/band-connectivity claim from the same fitted frame. No test is reported of an alternative admissible r(τ) choice or of a larger basis, so the conclusion that the hydrodynamizing mode emerges exclusively from the lowest band is at least partly built into the time-dependent basis ansatz rather than being shown to be a frame-independent property of the kinetic theory.

full rationale

The paper’s core new content — the numerical energy-level calculation showing that all but one of the low-band effective energies lift off and that the hydrodynamic mode is the remaining ground state — does not literally reduce to a single equation: it is obtained by solving the projected kinetic equation in a 12-state basis, and that calculation has independent content. However, two load-bearing elements are partially circular. First, the identification of an attractor with a zero effective energy eigenstate is a definitional restatement of Eqs. (2)-(4): Heff is defined through Heff w = -∂_y w, and prescaling is defined as ∂_y w = 0, so the “ground state” label is attached by construction. Second, the rescaling function r(τ) is chosen by an explicit Galerkin projection designed to make the first basis state an approximate solution; the later claims that f is well described by ψ10 and that the hydrodynamic mode has no mixing from higher bands are made in that same fitted frame, with no test of whether they survive a different r(τ), a different relaxation choice for D(τ), or a larger truncation. The proceedings also delegates numerical implementation details to the same-authors’ Ref. [3], which is a normal citation to the full paper rather than an appeal to authority, so I do not count the self-citation itself as a circular step. On balance, the central spectral observation is not forced by the definitions alone, but the interpretation of it as the adiabatic connection between attractors is entangled with the hand-chosen rescalings; this warrants a moderate circularity score of 4 rather than 0 or 2.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

Free parameters: three hand-chosen numbers (rate coefficient 10, initial condition parameters, basis size 12) affect the numerical demonstration but are not fitted to external data. Axioms: the geometric/kinetic theory assumptions are standard for the setup; the AH-specific axioms (zero-eigenvalue attractor, chosen rescalings, adiabaticity) are the framework's own postulates. The central claim depends on those postulates, so the contribution is an explanatory reinterpretation plus a numerical demonstration within a simplified model, not a parameter-free derivation. No new physical entities are introduced; the effective Hamiltonian and effective energy bands are mathematical constructs, not new degrees of freedom.

free parameters (3)
  • Rate coefficient in the D(τ) evolution ODE = 10
    Chosen by hand in ∂_yD/D = 10(1 - D/D_pE); it sets how fast the rescaling length D approaches its equilibrium value and therefore affects the adiabatic evolution and the effective energy gap structure.
  • Initial condition parameters (r_I, τ_I Q_s, σ_0) = r_I = √3, τ_I Q_s = 1, σ_0 = 1
    The single initial condition used for the numerical demonstration, Eq. (10); attractor claims are meant to be generic, but no variation of initial conditions is shown.
  • Number of basis states in the spectral projection = 12
    Truncation of the basis; no convergence study with respect to basis size is presented in this paper, so the stability of the band structure and no-mixing conclusion is not demonstrated here.
assumptions (5)
  • domain assumption The system is boost-invariant and homogeneous in the transverse plane, reducing the Boltzmann equation to Eq. (5).
    Section 2. This restricts the analysis to the standard longitudinal-expansion setup; it is a modeling assumption about the heavy-ion collision geometry.
  • domain assumption The collision kernel is the small-angle elastic kernel of Eq. (6), without inelastic scatterings and without explicit Bose enhancement in C[f], though (1+f) is retained in Ia.
    Section 2 and acknowledged in Sec. 4. The absence of number-changing processes makes this a simplified model; the paper does not claim to describe the full QCD kinetic theory.
  • ad hoc to paper The attractor can be characterized as the instantaneous zero-eigenvalue state of Heff defined by Heff w = -∂_y w.
    Eqs. (2)-(4). This is the defining construction of the AH framework; it builds the 'attractor equals ground state' identification into the formalism. Its validity is what the paper is testing for the specific kinetic theory.
  • ad hoc to paper The time-dependent scalings A(τ), D(τ), r(τ) and the 12-state polynomial basis adequately represent the true dynamics.
    Sec. 2. A is set by number conservation, D by a hand-chosen ODE with a rate coefficient 10, and r by a projection condition that forces the first basis state to describe f. The finite basis truncation is not tested in this paper.
  • domain assumption The adiabatic approximation holds: Heff evolves slowly enough relative to the effective energy gaps that excited modes decay as exp(-∫ ϵ_i dτ') and the ground state dominates.
    Sec. 1. This is the physical premise that lets the authors interpret the ground state as the attractor; the numerical results are meant to check it, but the paper does not quantify the adiabaticity condition.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Connecting Pre-Thermal and Hydrodynamizing Attractors With Adiabatic Hydrodynamization." pith.science (2026). https://pith.science/paper/F6SFIP46

@misc{pith2026250413332,
  author       = {Pith},
  title        = {Pith review of: Connecting Pre-Thermal and Hydrodynamizing Attractors With Adiabatic Hydrodynamization},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/F6SFIP46}},
  note         = {Machine review of arXiv:2504.13332}
}
read the original abstract

The far-from equilibrium dynamics of the pre-hydrodynamic quark-gluon plasma (QGP) formed in heavy ion collisions can be characterized by distinct stages, during each of which the system loses some memory of its initial condition, until only the hydrodynamic modes remain. However, even though it has been repeatedly observed, finding intuitive physical explanations of how and why attractor behavior occurs has remained a challenge. The Adiabatic Hydrodynamization (AH) framework provides exactly such an explanation, showing that the attractor solution can be thought of as the ground state of an analog to quantum mechanical adiabatic evolution, provided we identify appropriate coordinate rescalings. Using the example of a simplified QCD kinetic theory in the small-angle scattering limit, we show how AH can explain both the early pre-hydrodynamic attractor and the later hydrodynamizing attractor in a longitudinally expanding gluon gas in a unified framework. By doing this, we provide a unified description of, and intuition for, all the stages of what in QCD would be bottom-up thermalization, starting from a pre-hydrodynamic attractor and ending with hydrodynamization.

Figures

Figures reproduced from arXiv: 2504.13332 by the authors.

Figure 1
Figure 1. On the left, scaling exponents, and on the right, effective energy levels, each for gs = 10−3 . Scaling exponents are as defined in Eqns. (11) and (12). The good agreement between the two scaling exponent definitions suggests our choice of basis and scaling describes the evolution of f well. Both scaling regimes are characterized by well-separated bands of nearly degenerate modes, suggesting a partial loss of memory… view at source ↗
Figure 2
Figure 2. On the left, scaling exponents, and on the right, effective energy levels, each for gs = 1. Scaling exponents are again as defined in Eqns. (11) and (12), and as in the weak coupling case agreement suggests a good time-dependent basis choice. Between times approximately y = 3 and y = 15, both scaling exponents and effective energy levels seem to resemble the dilute scaling regime, although not as cleanly as at weake… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

9 extracted references · 7 linked inside Pith

  1. [3]

    Rajagopal, B

    K. Rajagopal, B. Scheihing-Hitschfeld, R. Steinhorst, Adiabatic Hydrodynamization and the Emergence of Attractors: a Unified Description of Hydrodynamization in Kinetic Theory (2024), 2405.17545

    arXiv 2024

  2. [7]

    Brewer, B

    J. Brewer, B. Scheihing-Hitschfeld, Y. Yin, Scaling and adiabaticity in a rapidly expanding gluon plasma , JHEP 05, 145 (2022), 2203.02427 . 10.1007/JHEP05(2022)145

    arXiv 2022

  3. [1]

    Kurkela, W

    A. Kurkela, W. van der Schee, U.A. Wiedemann, B. Wu, Early- and Late-Time Behavior of Attractors in Heavy-Ion Collisions , Phys. Rev. Lett. 124, 102301 (2020), 1907.08101 . 10.1103/PhysRevLett.124.102301

    arXiv 2020

  4. [2]

    Brewer, L

    J. Brewer, L. Yan, Y. Yin, Adiabatic hydrodynamization in rapidly-expanding quark gluon plasma , Phys. Lett. B 816, 136189 (2021), 1910.00021 . 10.1016/j.physletb.2021.136189

    arXiv 2021

  5. [4]

    Berges, K

    J. Berges, K. Boguslavski, S. Schlichting, R. Venugopalan, Turbulent thermalization process in heavy-ion collisions at ultrarelativistic energies , Phys. Rev. D 89, 074011 (2014), 1303.5650 . 10.1103/PhysRevD.89.074011

    arXiv 2014

  6. [5]

    Baier, A.H

    R. Baier, A.H. Mueller, D. Schiff, D.T. Son, 'Bottom up' thermalization in heavy ion collisions , Phys. Lett. B 502, 51 (2001), hep-ph/0009237 . 10.1016/S0370-2693(01)00191-5

    arXiv 2001

  7. [6]

    Mazeliauskas, J

    A. Mazeliauskas, J. Berges, Prescaling and far-from-equilibrium hydrodynamics in the quark-gluon plasma , Phys. Rev. Lett. 122, 122301 (2019), 1810.10554 . 10.1103/PhysRevLett.122.122301

    arXiv 2019

  8. [8]

    , " * write output.state after.block = add.period write newline

    ENTRY address archive author booktitle chapter collaboration edition editor eid eprint howpublished institution isbn issn journal key month note number numpages organization pages publisher school series title type url doi volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION ...

Show all 9 references
  1. [9]

    write FUNCTION output.bibitem newline

    " write FUNCTION output.bibitem newline " " " " * write cite write " " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" *...

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.