{"id":"ed3dd0c0-5df9-45eb-8496-67a115982381","arxiv_id":"1908.02113","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":2.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"This review organizes the weak-coupling theory of early-time thermalization in heavy-ion collisions around over- and under-occupied plasma stages and quotes a hydrodynamization time of about 1 fm/c.","lead":"Two heavy-ion theorists review what happens in the first instant after two nuclei collide, when a dense soup of quarks and gluons forms. The review explains the main weak-coupling picture, bottom-up thermalization, and gives a concrete estimate for when the soup starts behaving as a fluid.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Quantitative hydrodynamization time (eq. 81) lacks an uncertainty estimate for the leading-order-to-realistic-coupling extrapolation; the λ-scaling universality in Fig. 5 does not test the absolute relaxation rate.","rationale":"Good-faith reading: this is an accurate and well-cited review of the weak-coupling thermalization program. It explicitly flags the two main caveats—the extrapolation to realistic coupling (Sec. 5.1) and the incomplete understanding of plasma instabilities (Sec. 4.3)—and the numerical results (Kurkela-Zhu, Kurkela-Mazeliauskas) provide real, citable evidence for the bottom-up attractor picture. The review's synthesis is faithful to the literature and its stated limitations. My concern targets the quantitative headline claim, eq. (81), rather than the qualitative bottom-up picture. The paper does not attempt to quantify the systematic error from truncating the weak-coupling expansion at LO while pushing λ to 10-25; this is the least secure link between ab initio kinetic theory and the quoted 1.1 fm number. I agree with the reader's weakest_assumption in spirit, but would isolate the coupling extrapolation as the single most load-bearing element, since plasma instabilities are at least argued (and partly supported by classical-statistical simulations) to be bounded and not to change the qualitative picture. The proposed check—an NLO estimate of the relaxation rate and its effect on eq. (81)—would decide whether the extrapolation lands. Because the review is a synthesis rather than a claim of new results, and because it already discloses the need for extrapolation, this concern does not overturn the ACCEPT verdict; it argues for treating the 1.1 fm number as an order-of-magnitude estimate unless the NLO check is performed.","tokens_in":25771,"tokens_out":17940,"duration_ms":186635,"concrete_test":"Recompute eq. (81) with next-to-leading-order input for η/s and for the relaxation rate, using the NLO kinetic-theory results of Ghiglieri, Moore, and Teaney (JHEP 03:179 (2018)) at the same λ values (10 and 25) and with the same T_Id normalization. Alternatively, rerun the Kurkela-Mazeliauskas kinetic code with an NLO-corrected collision kernel (or with the NLO η/s) and extract τhydro from the matching to viscous hydrodynamics. If τhydro changes by more than 0.3 fm at η/s≈0.16, or by more than 30% at either coupling, the extrapolation is not controlled and eq. (81) should carry a large systematic uncertainty; if the shift is smaller, the quantitative claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is that weak-coupling kinetic theory yields the hydrodynamization time in eq. (81), τhydro≈1.1 fm for central Pb+Pb. The underlying simulations solve leading-order (LO) QCD effective kinetic theory at 't Hooft couplings λ=10-25, i.e. α_s≈0.27-0.66. At these couplings the LO quasiparticle expansion is not under perturbative control: collision kernels are LO, screening is an isotropic Debye mass, and the η/s used to map λ is itself computed in LO kinetic theory. Section 5.1 concedes that 'an extrapolation to sizeable coupling strength is required,' but the review provides no estimate of the resulting systematic error. The apparent coupling-independence of Fig. 5 does not settle this: scaling time by the theory's own τR^eq=4πη/s/T_Id removes the coupling dependence of the overall normalization, so it cannot validate the absolute magnitude of τR^eq. An O(1) error in the LO rate enters eq. (81) as (2)^{3/2}≈2.8, shifting τhydro from ~1 fm to ~3 fm and moving the small-system threshold in eq. (84) accordingly. Thus the quantitative bridge to data rests on an unquantified extrapolation of LO weak-coupling kinetics.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript is an invited review of the weak-coupling description of the earliest stage of ultra-relativistic heavy-ion collisions. After a concise introduction to the CGC initial state and the BMSS bottom-up scenario, the authors review the elements of leading-order QCD kinetic theory (elastic scattering, Fokker-Planck transport coefficients, collinear radiation with LPM suppression), derive the scaling exponents for overoccupied and underoccupied plasmas, and summarize numerical simulations of the approach to hydrodynamics. The review concludes with quantitative estimates for the hydrodynamization time tau_hydro ~ 1.1 fm (Eq. 81), for pre-flow response functions (Eqs. 82-83), and for the small-system onset of hydrodynamics (Eq. 84).","tokens_in":26059,"tokens_out":7296,"duration_ms":73197,"significance":"The review is valuable as a pedagogical and critical synthesis. It correctly identifies the universal scaling results (alpha = -4/7, beta = -1/7; Q_s tau ~ alpha_s^(-13/5)) and presents them with transparent derivations that a reader can follow. The combination of parametric estimates, kinetic-theory simulations, and phenomenological constraints in one place will be useful to the community. The most important caveat is that the quantitative bridge to phenomenology is built on leading-order kinetic theory at couplings where the expansion is not controlled; this caveat is acknowledged in places but not quantified in the central formulas.","major_comments":[{"comment":"The hydrodynamization time in Eq. (81) and the small-system threshold in Eq. (84) are load-bearing quantitative claims, but the manuscript does not provide an uncertainty estimate for the extrapolation from leading-order kinetic theory at lambda = 10-25 to the couplings relevant to Pb+Pb. Rescaling time by tau_R^eq = 4 pi eta/s / T_Id in Fig. 5(a) removes the coupling dependence of the overall normalization and therefore does not test the absolute value of tau_R^eq. Because the same LO rate enters the translation from lambda to eta/s, an O(1) error in the rate changes Eq. (81) by a factor (2)^(3/2) and shifts Eq. (84) by the corresponding amount. The authors should either supply a quantitative systematic-error estimate for Eqs. (81) and (84) or state explicitly that these are parametric, order-of-magnitude estimates.","section":"Section 5.1, Eq. (81), Fig. 5(a)"},{"comment":"The quantitative statements in Section 5.1 and Section 5.3 assume that plasma instabilities do not qualitatively modify the bottom-up evolution, but Section 4.3 states that the role of instabilities in underoccupied systems \"has not been fully clarified\" and that current kinetic theory implementations ignore them. Since Eq. (81) and the pre-flow response functions of Eq. (82) inherit this assumption, the review should either quantify the expected sensitivity of tau_hydro to instabilities, for instance through the modified exponents of Ref. 63, or carry a prominent caveat in Section 5.1. As written, the main quantitative claims are stronger than the stated theoretical basis.","section":"Sections 4.3 and 5.1"}],"minor_comments":[{"comment":"The phrase \"a brief expose\" should be \"a brief exposé\".","section":"Abstract"},{"comment":"The word \"mircoseconds\" in the introductory paragraph is a typo for \"microseconds\".","section":"Section 1"},{"comment":"The text \"produced by direction radiation by the bath\" should read \"produced by direct radiation from the bath\".","section":"Section 4.3"},{"comment":"The sentence containing \"an be transformed\" should read \"can be transformed\".","section":"Section 4.2.2, discussion after Eq. (66)"},{"comment":"The equation for e_soft(tau) in the second phase is again numbered (58), and the equation for e_soft(tau) in the final phase is again numbered (64), duplicating the numbers already used in Section 4.2; these should be renumbered or cited with reference to the original equations.","section":"Section 4.3"}],"recommendation":"major_revision","confidential_remarks":"The review is competent and likely useful, but the strongest quantitative claims rest on simulations and analyses from the authors' own group (Refs. 30, 71, 72). This is not a reason to reject a review, but the editor may wish to ask the authors to state explicitly, where relevant, which numerical results have been independently reproduced."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a competent, honest review, not a research paper. If you need a single accessible account of the weak-coupling bottom-up program with the modern kinetic-theory numbers, this is the go-to reference. The catch is that its headline number — τhydro ≈ 1.1 fm, eq. (81) — inherits an unquantified extrapolation from leading-order kinetic theory at λ=10–25 to realistic coupling, and the paper says so but never tells you how uncertain the number is.\n\nWhat is genuinely good: the review organizes the material better than the original papers. The Fokker-Planck and LPM rate discussion in Sec. 3 is clear and internally consistent. The distinction between overoccupied and underoccupied systems, with the scaling exponents α=-4/7 and β=-1/7 and the bottom-up time Qsτ∼α_s^-13/5, is presented accurately and tied to the right sources. Section 5.3's description of pre-flow response functions is a genuinely useful package for phenomenology. The paper also flags its own limitations: Sec. 5.1 concedes the extrapolation to sizeable coupling, and Sec. 4.3 says plasma instabilities are not fully clarified. That honesty is real, not cosmetic.\n\nWhere I part company with the reader's report only slightly: the stress-test note is right, and the reader's soundness=8 is fair but maybe generous for the quantitative claim. Fig. 5's apparent coupling independence does not test the absolute relaxation rate, because time is scaled by the theory's own τR^eq=4πη/s/T_Id. That scaling removes the normalization issue. An O(1) error in the LO rate changes eq. (81) by (2)^{3/2}≈2.8, moving τhydro from ~1 fm to ~3 fm and shifting the small-system threshold in eq. (84). The review acknowledges the need for extrapolation but gives no error bar. For a review that is not fatal, but it should be said explicitly, and the paper should add a sentence or two quantifying this.\n\nThe self-citation point is a non-issue here. Several anchors come from overlapping-author papers, but this is a review by people who did much of the work; the underlying simulations are published and separately checkable. I would not downgrade for that.\n\nBottom line: for a graduate student or a phenomenologist wanting the state of the weak-coupling early-time program, this paper is worth reading. It is not original research, and the central number needs a caveat. A serious referee should engage with it; the right outcome is acceptance after minor revision, with the uncertainty on eq. (81) made explicit.","headline":"A well-executed review of weak-coupling early-time dynamics whose headline hydrodynamization time is solid as a synthesis but carries an unstated LO-to-realistic-coupling uncertainty.","tokens_in":26559,"tokens_out":2059,"would_cite":true,"duration_ms":21866,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q","12.38.Mh"],"model":"deepseek-v4-flash","headline":"Weak-coupling QCD kinetic theory puts the first fm/c of a heavy-ion collision within quantitative reach.","keywords":["heavy-ion collisions","quark-gluon plasma","QCD kinetic theory","thermalization","bottom-up scenario","hydrodynamization","pre-equilibrium dynamics","color-glass condensate"],"falsifier":"Take the strongest form: start an expanding, anisotropic SU(3) plasma at $\\lambda\\approx10$ from the same CGC-motivated initial condition and evolve it with a method that includes plasma instabilities, such as classical-statistical field theory at early times; if the normalized pressure components $P_L/e$ and $P_T/e$ deviate from the kinetic-theory curves in the review's Fig. 5(a) by more than the expected viscous-hydro spread at $\\tau/\\tau_R^{\\mathrm{eq}} \\approx 1$, then eq. (81) and the claimed hydrodynamization time are not supported.","tokens_in":25587,"feed_emoji":"⚛️","tokens_out":9726,"duration_ms":91076,"temperature":0.7,"pith_summary":"This review argues that the first fm/c of a heavy-ion collision—the stretch of time before a quark-gluon plasma becomes a fluid—can be described from first principles by weakly coupled QCD kinetic theory. The bottom-up thermalization picture, which began as parametric estimates, is now realized in numerical simulations that start from Color-Glass-Condensate initial conditions and evolve the full energy-momentum tensor through overoccupied and underoccupied stages. The quantitative payoff is a hydrodynamization time for central lead-lead collisions, $\\tau_{\\mathrm{hydro}} \\approx 1.1$ fm, with the dependence on shear viscosity and multiplicity given in closed form. If the picture is right, the pre-equilibrium phase is not a free parameter but a calculable bridge from the initial state to the hydrodynamic initial conditions.","feed_headline":"The first fm/c of a heavy-ion collision is calculable","feed_subtitle":"Kinetic theory bridges the initial gluon state to viscous hydrodynamics at τ≈1.1 fm.","key_machinery":"The central object is the bottom-up thermalization scenario, the three-stage weak-coupling mechanism by which a highly anisotropic, overoccupied gluon plasma relaxes: momentum diffusion and longitudinal expansion compete in the first stage, radiation populates a soft bath in the second, and an inverse energy cascade transfers energy from hard to soft modes in the third. The quantitative machinery is the leading-order QCD Boltzmann equation with $2\\leftrightarrow 2$ elastic scattering and $1\\leftrightarrow 2$ LPM-suppressed collinear radiation, solved numerically in the expanding boost-invariant geometry. The account is carried by a chain of parametric and numerical results: the non-thermal fixed point with scaling exponents $\\alpha = -4/7$ and $\\beta = -1/7$ for overoccupied plasmas; the turbulent inverse-cascade spectrum with spectral index $\\kappa = 7/2$ for the underoccupied cascade; and the criterion that hydrodynamics applies once the proper time is of order the equilibrium relaxation time, $\\tau_{\\mathrm{hydro}} \\approx \\tau_R^{\\mathrm{eq}}$.","core_discovery":"The review's central claim is that weak-coupling QCD kinetic theory quantitatively describes the first fm/c of a heavy-ion collision, making viscous hydrodynamics applicable by the time quoted in eq. (81): $\\tau_{\\mathrm{hydro}} \\approx 1.1$ fm $\\left(4\\pi(\\eta/s)/2\\right)^{3/2} \\left(\\langle\\tau s\\rangle/(4.1~\\mathrm{GeV}^2)\\right)^{-1/2} \\left(\\nu_{\\mathrm{eff}}/40\\right)^{1/2}$. The same framework converts the original bottom-up scenario into a three-stage sequence—an overoccupied gluon plasma that first becomes more anisotropic, a phase in which a soft thermal bath builds up through LPM-suppressed radiation, and a final inverse energy cascade that transfers energy from hard modes to the bath. Numerical solutions of the leading-order QCD Boltzmann equation with elastic and inelastic processes reproduce this sequence and show that the energy-momentum tensor approaches the hydrodynamic asymptotic behavior even while the pressure anisotropy remains of order one. The paper packages the result as a pre-flow linear-response kernel that connects classical-field initial conditions, kinetic theory, and hydrodynamic initial conditions in event-by-event simulations.","pith_inferences":["Editorial extension: if eq. (81) survives contact with data, the same kinetic-theory machinery should determine the pre-equilibrium electromagnetic emissivity (photons and dileptons) once the quark sector is treated at full leading order; the review notes that complete calculation is still missing.","Editorial extension: the multiplicity threshold $dN_{\\mathrm{ch}}/dy \\approx 70$ is a sharp, falsifiable bookmark: flow-like correlations in small collision systems should turn off near that value if the weak-coupling picture is the whole story, whereas a much lower threshold would suggest strong-coupling or non-kinetic mechanisms.","Editorial extension: because the review brackets plasma instabilities as an open issue, a direct test would be to compare the kinetic-theory pressure curves with classical-statistical simulations of expanding anisotropic plasmas at $\\lambda \\approx 10$ that include instabilities; the two should agree at the 10–20% level at $\\tau/\\tau_R^{\\mathrm{eq}} \\sim 1$ if the extrapolation is safe."],"forward_implications":["Hydrodynamic initial conditions become calculable: eq. (81) fixes when and how the pre-equilibrium $T^{\\mu\\nu}$ is matched to viscous hydrodynamics for large collision systems.","Hydrodynamization precedes local equilibrium: the fluid equations become reliable while the longitudinal and transverse pressures still differ by an order-one amount, so the success of hydrodynamics does not imply rapid isotropization.","Pre-equilibrium production is quantified: the kinetic-theory evolution generates a factor of roughly two to three in entropy before hydrodynamics, which directly links initial-state models to measured charged-particle multiplicities.","Event-by-event pre-flow can be computed: the linear response in eq. (82) gives the pre-equilibrium buildup of transverse flow from initial geometry fluctuations, replacing ad hoc hydrodynamic initialization.","Small systems have a concrete threshold: using eq. (84), hydrodynamics becomes marginal when the charged-particle multiplicity per unit rapidity is of order $dN_{\\mathrm{ch}}/dy \\approx 70$, giving a target for small-system flow studies."],"supporting_citations":[{"why":"Seminal bottom-up thermalization scenario: provides the three-stage parametric cascade that the review quantifies.","marker":"(6)"},{"why":"Analytic matching of classical color fields across the collision, giving the CGC initial conditions for the stress tensor.","marker":"(9, 10, 11)"},{"why":"Leading-order QCD Boltzmann equation and splitting rates used as the kinetic theory's collision kernel.","marker":"(33)"},{"why":"Parametric analysis of overoccupied and underoccupied thermalization and the direct-radiation estimate for the soft bath.","marker":"(36)"},{"why":"Classical-statistical simulations establishing the non-thermal fixed point and scaling exponents in overoccupied Yang-Mills plasmas.","marker":"(48)"},{"why":"Kinetic-theory simulations giving concrete thermalization times and confirming the cascade picture.","marker":"(49)"},{"why":"Numerical simulation of the full bottom-up scenario from CGC-motivated initial conditions that underlies the hydrodynamization curves.","marker":"(67)"},{"why":"Linear-response pre-flow formalism and the package connecting kinetic theory to hydrodynamics shown in the review's figures.","marker":"(30)"},{"why":"Constrains $\\langle\\tau s\\rangle$ for central lead-lead collisions and demonstrates hydrodynamic description of the pre-equilibrium $T^{\\mu\\nu}$; basis for eq. (81).","marker":"(71)"}],"fun_headline_variants":["Kinetic theory quantifies the first fm/c of heavy-ion collisions","First fm/c of heavy-ion collisions now quantitatively understood","From glasma to hydrodynamics: the first fm/c decoded","Weak-coupling QCD predicts when hydrodynamics kicks in","Kinetic theory bridges the first fm/c to hydrodynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that kinetic-theory results computed at weak coupling remain quantitatively valid when the coupling is extrapolated to realistic values ($\\lambda\\approx10$–$25$, i.e. $\\eta/s\\approx0.62$–$0.16$) and that plasma instabilities do not change the bottom-up dynamics; the paper itself flags both as unproven.","fun_headline_variants_meta":{"raw":{"variants":["Kinetic theory quantifies the first fm/c of heavy-ion collisions","First fm/c of heavy-ion collisions now quantitatively understood","From glasma to hydrodynamics: the first fm/c decoded","Weak-coupling QCD predicts when hydrodynamics kicks in","Kinetic theory bridges the first fm/c to hydrodynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00089,"raw_usage":{"total_tokens":3823,"prompt_tokens":910,"completion_tokens":2913,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":2830}},"tokens_in":526,"tokens_out":2913,"duration_ms":38726,"temperature":1.0,"reasoning_tokens":2830,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:54:06.640091+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the strongest form: start an expanding, anisotropic SU(3) plasma at $\\lambda\\approx10$ from the same CGC-motivated initial condition and evolve it with a method that includes plasma instabilities, such as classical-statistical field theory at early times; if the normalized pressure components $P_L/e$ and $P_T/e$ deviate from the kinetic-theory curves in the review's Fig. 5(a) by more than the expected viscous-hydro spread at $\\tau/\\tau_R^{\\mathrm{eq}} \\approx 1$, then eq. (81) and the claimed hydrodynamization time are not supported.","supporting_citations":[],"review_version":1}