Pith. Sign in

REVIEW

The Boltzmann equation for gluons at early times after a heavy ion collision

High Energy Physics - Phenomenology · Submitted 15 Sep 1999 · 7 references

Pith's verdict: A clean, honest framework note that reduces the gluon Boltzmann equation to Fokker-Planck form, but delivers parametric scaling rather than evaluated numbers because the MV normalization c and the closure of η^{-1} are left undetermined.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on November 13, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-ph/9909388 (version 1)
Pith IDYWNXSNO5, Pith's permanent record of this paper and its review
CitationA.H. Mueller (Columbia University). "The Boltzmann equation for gluons at early times after a heavy ion collision." arXiv:hep-ph/9909388 (1999). https://pith.science/paper/YWNXSNO5
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Mueller rewrites the first instants of a heavy-ion collision as a Boltzmann equation for gluons, reduces its collision term to a Fokker-Planck drift-and-diffusion form, and reads off how transverse momentum shrinks while longitudinal动量...

The reading

Mueller asks what happens to the gluon cloud in the first instants after a head-on heavy-ion collision, before anything like equilibrium has set in. His answer is to stop tracking individual scatterings and instead write a Boltzmann equation for the gluon occupation number $f(t,p_​z,p_​z)$ whose collision term comes from gluon-gluon scattering in the one-gluon-exchange approximation. In the small-angle logarithmic approximation that collision term collapses into a Fokker-Planck operator with a diffusion piece and a drift piece, Eq. (15), and the drift strength is what actually pushes momentum out of the transverse plane and into the longitudinal direction. Without solving the equation, he extracts its first moments: the typical transverse momentum falls as $1-(cN_c/4\pi)\alpha\xi^2$ in the logarithm $\xi=\ln(t/t_0)$, while $\langle p_z^2\rangle$ rises linearly in $\xi$, so the system is driven toward the isotropic condition $\langle p_z^2\rangle\sim\langle p_\perp\rangle^2$. The payoff is a definite time, $\xi_1=\sqrt{2\pi/(c\alpha N_c)}$, at which transverse and longitudinal momenta become comparable and the early-time treatment stops being usable — which he reads as the onset of equilibration.

What carries the argument

The load-bearing object is the Fokker-Planck reduction of the Boltzmann collision term, Eq. (15): $(\partial/\partial\xi - \partial/\partial p_z\,& p_z)\tilde f = \lambda\eta\xi\,¦^2_{\vec p}\tilde f + 2\lambda\eta^{-1}\xi\,¦_{\vec p}\cdot(\vec v\tilde f)$, with $\lambda = 2\pi\alpha^2 N_c^2/(3(N_c^2-1))$ and $\xi = \ln(t/t_0)/3$ coming from the small-angle integral cut off at $\theta_m^3=t_0/t$. The first term is diffusion in momentum; the second is the drift that converts transverse into longitudinal momentum and is the mechanism the whole paper turns on. The second piece of machinery is the moment hierarchy obtained by multiplying Eq. (15) by $p_\perp^2$ and $p_z^2$ and integrating: Eqs. (18) and (19) are exact consequences of the Fokker-Planck form, and Eqs. (25), (30), (31) and (33) follow once the moments of $\tilde f$ are replaced by their initial-state values, which is where the McLerran-Venugopalan model enters through $\eta$ and $\eta^{-1}$.

Extended reading notes

Core claim

The paper's central claim is that the early evolution of the gluon system produced in a heavy-ion collision is governed by a Boltzmann equation whose collision term, taken from elastic gluon-gluon scattering in one-gluon exchange and integrated over small scattering angles in a logarithmic approximation, reduces to a Fokker-Planck equation for the occupation number, Eq. (15). Two parameters enter: a time-independent $\eta$ fixed by the initial conditions, and a drift coefficient $\eta^{-1}$ that in principle should come self-consistently from the solution. Integrating Eq. (15) over momentum shows the gluon number density falls as $1/t$. Taking moments gives $\langle p_\perp\rangle/\bar Q \simeq 1-(c\alpha N_c/4\pi)\xi^2$ (refined to $\langle p_\perp\rangle^2/\bar Q^2 = 1-(c\alpha N_c/2\pi)\xi^2$ when the time variation of $\eta^{-1}$ is modeled) and $\langle (p_z/\bar Q)^2\rangle \simeq (c\alpha N_c/4\pi)\xi$ with $\xi=\ln(t/t_0)$, so the two become comparable at $\xi_1=\sqrt{2\pi/(c\alpha N_c)}$. Beyond $\xi_1$ the author cannot follow the system analytically, but presumes transverse and longitudinal momenta stay equal at $z=0$ and the system continues toward equilibrium.

Load-bearing premise

The strength of the drift term — how fast scattering moves gluons from transverse into longitudinal momentum — is not determined from the solution itself but is imported from the initial-state model and then extrapolated in time by an ansatz, so the numerical coefficients and the predicted breakdown time rest on that closure rather than on the equation alone.

Editorial extensions

If this is right

  • The gluon number density at $z=0$ falls as $1/t$ while the system is in this regime, a direct consequence of integrating Eq. (15) over momentum.
  • Transverse momentum is depleted quadratically in $\xi=\ln(t/t_0)$ and longitudinal momentum squared grows linearly in $\xi$, giving a definite, calculable shape for the approach to isotropy at early times.
  • The early-time description has a built-in expiry: at $\xi_1=\sqrt{2\pi/(c\alpha N_c)}$ the transverse and longitudinal momenta become comparable, and the approximations behind Eqs. (30), (31) and (33) fail.
  • Because $\xi\approx\ln(1/\alpha)$ already at $t\sim 1/Q_s$ in the McLerran-Venugopalan model, the logarithmic approximation that justifies the Fokker-Planck form is parametrically valid from the moment the gluons are freed.
  • The framework gives a concrete target for a numerical solution of Eq. (15): the coefficient of $\xi^2$ in $\langle p_\perp\rangle$ and the crossing time $\xi_1$ are the quantities a full solution should reproduce or correct.

Reading between the lines

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

  • If the drift coefficient is genuinely self-consistent rather than imported, the crossing time $\xi_1$ becomes a prediction of the theory rather than a parameter, which would turn the onset of isotropy into a calculable number for comparison with transport models of the early collision stage.
  • The same moment method applied to higher moments of Eq. (15) would give the full shape of the distribution at $\xi_1$, not just its width, and would test whether the assumption that $\langle p_z^2\rangle$ and $\langle p_\perp\rangle^2$ stay equal past the crossing is consistent with the equation itself.
  • Because the equation is boost invariant and written at $z=0$, the results translate into a rapidity-space statement: the drift term acts as a diffusion in rapidity, which suggests comparing the predicted spread with measured hadron spectra at low transverse momentum.
  • The author's treatment of $t_0$ in the infrared cutoff as constant is the natural next thing to relax; a time-dependent cutoff would feed back into $\xi$ and shift the coefficients in Eqs. (30) and (33) without changing their functional form.
Original abstract

A Boltzmann equation is given for the early stages of evolution of the gluon system produced in a head-on collision. The collision term is taken from gluon- gluon approximation. < p_\perp > and < p_z^2 > are evaluated as a function of time using initial conditions taken from the McLerran-Venugopalan model.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 2,311 input tokens, $0.0001 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.30 · human participants 0.01 · controlled comparison 0.03 · randomized assignment 0.01 · states quantitative result 0.13 · compares with prior results 0.06 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

A Boltzmann equation is given for the early stages of evolution of the gluon system produced in a head-on collision. The collision term is taken from gluon- gluon approximation. < p_\perp > and < p_z^2 > are evaluated as a function of time using initial conditions taken from the McLerran-Venugopalan model.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 A Boltzmann equation is given for the early stages of evolution of the gluon system produced in a head-on collision. contribution 0.97 0.82 0.17 asserted 1.00
s1 The collision term is taken from gluon- gluon approximation. < p_\perp > and < p_z^2 > are evaluated as a function of time using initial conditions taken from the McLerran-Venugopalan model. method 0.97 0.81 0.10 asserted 1.00

What kind of paper the classifier reports this is: theoretical result 0.96 · evidence from simulation or computation 0.60.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

References
How the literature uses it

33 arXiv papers cite this in 57 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 37 sentences without a confident label

Builds on it 1 paper
Uses its method or tool 7 papers
  • As recently detailed in [ref] , the $2↔2$ kernel is simplified using the diffusion approximation [ref] : [equation] where $ϵ_a=1$ for bosons and $ϵ_a=-1$ for fermions and the source terms [ref] take the form [equation] The space-time-dependent quantities in the above equations, which govern the time evolution of the system, are defined as:

    Quark production in the bottom-up thermalization · The BEDA for longitudinally boost-invariant systems · 0.92

  • Using the same setup as in [ref] , we describe the plasma in terms of quasiparticle distribution functions $f_p≡f(τ,p)$, whose evolution in proper time $τ$ is given by the Boltzmann equation [equation] Here $$ encodes collisions resulting in elastic scatterings, $$ describes effective splitting terms relevant for gluon radiation, and $= - p_z/τ [f_p]p_z$ accounts for the longitudinal expansion of the plasma in the beam direction [ref] .

    Jet momentum broadening during initial stages in heavy-ion collisions · Kinetic theory · 0.97

  • We shall demonstrate the correspondence between the AH and QNM frameworks by making use of the Fokker-Planck-type Boltzmann equation [ref] [equation]

    Adiabatic hydrodynamization and quasinormal modes of nonthermal attractors · Introduction · 0.89

Compares with it 1 paper
Cites it as support 1 paper
Cites it as background 8 papers
  • For an initially over-occupied system of gluons one expects the so-called “bottom-up” equilibration scenario [ref] that indicates a rapid growth of the momentum broadening of hard partons followed by a gradual decay [ref] .

    Sensitivity of jet quenching to the initial state in heavy-ion collisions · Introduction · 0.80

  • The pre-thermal quark-gluon plasma generated in heavy-ion collisions follows the bottom-up thermalization scenario [ref] , which consists of several stages: The first stage is dominated by a large number of hard gluons, and the anisotropy increases due to the longitudinal expansion along the beam axis.

    Jet momentum broadening during initial stages in heavy-ion collisions · Bottom-up thermalization · 0.99

  • Significant progress has been made in the last decade in developing a microscopic picture of the rapid hydrodynamization of the quark-gluon plasma (QGP) formed in ultrarelativistic heavy ion collisions, both in strongly-coupled holographic descriptions [ref] , and in weakly-coupled descriptions employing QCD classical field simulations [ref] , QCD effective kinetic theory (QCD EKT) [ref] , and related kinetic theory descriptions of QCD [ref] .

    Attractors Without Scaling: Adiabatic Hydrodynamization With and Without Inelastic Scattering · Introduction · 1.00

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.