Pith. sign in

REVIEW 3 major objections 5 minor 12 references

The (3+1)D structure of the dilute Glasma

T0 review · 3 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash

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

Pith's one-line read The dilute Glasma's large-rapidity energy spike is a coordinate artifact

desk verdict Useful summary of the dilute Glasma, but the two headline claims are either a restatement of the IR regulator or an under-derived coordinate argument. read the letter →

arxiv 2508.20725 v1 pith:YFVJT5M2 submitted 2025-08-28 hep-ph nucl-th

Andreas Ipp , Markus Leuthner , David I. Müller , Sören Schlichting , Kayran Schmidt , Pragya Singh This is my paper
classification hep-phnucl-th
keywords colorglasscondensateGlasmaheavy-ioncollisionsclassicalYang-MillsMcLerran-Venugopalanmodelenergy-momentumtensorMilnecoordinatessaturationscale
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 proceedings paper claims that in the dilute, leading-order approximation of the color glass condensate, the full (3+1)-dimensional Glasma field after a heavy-ion collision can be computed event by event, and that two features of the result overturn naive readings of standard Milne-coordinate plots. The large enhancement of the T^{ττ} component at high spacetime rapidity is not deposited energy; it is a coordinate effect caused by large longitudinal pressures in the co-moving Bjorken frame. In the transverse plane, the size of the fluctuating energy-density domains is set by the infrared regulator 1/m rather than by a saturation scale, because the dilute approximation never develops Q_s. If these claims are right, boost-invariant Milne-coordinate treatments misplace forward-rapidity energy and cannot capture the dynamics of initial conditions with extended longitudinal structure.

What carries the argument

The dilute (3+1)D Glasma ansatz: in covariant gauge the total field is A^μ_A + A^μ_B + a^μ, with the non-linear remainder a^μ carrying the collision-induced Glasma field; at leading order in ρ_A ρ_B, f^{μν} reduces to an abelian wave-equation solution expressed as integrals of commutators of the single-nucleus field strengths over the backward light cone. This expression is evaluated numerically with Monte Carlo integration on a (τ, η_s, x) grid. The nuclear model is a generalized McLerran-Venugopalan model: Gaussian color charges with a finite (1+1)-dimensional envelope T(x^±, x) and a longitudinal correlation function U_ξ; the IR regulator m in the momentum-space solution for A^∓ sets the

What would settle it

Run the same three-dimensional nuclear model through the full classical Yang-Mills equations without the ρ_A ρ_B truncation, at RHIC/LHC values of g^2 μ, and compare the rapidity dependence of T^{ττ} and ε_LRF with the dilute result; if the large-rapidity enhancement vanishes or changes character, the dilute coordinate-artifact interpretation does not extend to the realistic case. A second, cheaper check: vary the IR regulator m and measure the transverse correlation width of ε_LRF — if the width does not follow 1/m, the dilute domain-size claim is wrong.

Watch

Extended reading notes

Core claim

Working to leading order m=n=1 in an expansion in the color charge densities ρ_A^m ρ_B^n of the two nuclei, the authors write the Glasma field strength as an effectively abelian tensor f^{μν}=∂^μ a^ν−∂^ν a^μ whose sources are the known single-nucleus field strengths, integrated over the backward light cone of the observer. Using a three-dimensional generalization of the McLerran-Venugopalan model with finite longitudinal extent, a tunable longitudinal correlation scale ξ, and an infrared regulator m, they compute T^{μν} and the local rest frame energy density ε_LRF event by event. They find that T^{ττ} grows with rapidity while ε_LRF and T^{xx}+T^{yy} do not, and they attribute the growth to

Load-bearing premise

The whole calculation assumes the leading-order dilute limit — where the gluon field is effectively abelian and no saturation scale forms — still captures the longitudinal dynamics that matter for real heavy-ion collisions; if nonlinear color-charge terms are important at RHIC/LHC energies, both the Milne-frame artifact diagnosis and the 1/m domain-size result would not carry over.

Editorial extensions

If this is right

  • Milne-coordinate (Bjorken co-moving frame) plots of T^{ττ} should not be interpreted as physical energy deposition at forward rapidity for collisions with extended longitudinal structure; a more suitable frame or observable is needed.
  • In the dilute approximation, transverse domain sizes are controlled by the ad hoc IR regulator 1/m, not by a saturation scale Q_s∝g^2 μ, so dilute-Glasma initial conditions cannot claim to predict correlation lengths from first principles.
  • The gauge-invariant local rest frame energy density ε_LRF and the transverse pressure sum T^{xx}+T^{yy} are free of the large-rapidity artifact and can serve as robust diagnostics of longitudinal structure.
  • Event by event, ε_LRF is deposited only in the transverse overlap region of the nuclear profiles, so three-dimensional nuclear structure leaves an imprint in the initial-state energy-density pattern.
  • Varying the longitudinal correlation scale ξ changes the rapidity width of the profiles, providing a direct handle on how nucleon longitudinal structure enters the Glasma.

Reading between the lines

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

  • If the coordinate-artifact diagnosis survives in the full nonlinear Glasma, then existing boost-invariant initial condition models may systematically misassign energy to large rapidities, and switching to Cartesian-frame or lab-frame observables could alter the starting profiles used in hydrodynamics.
  • A direct test: run full classical Yang-Mills simulations (no ρ_A ρ_B truncation) with the same three-dimensional nuclear model; if the T^{ττ} enhancement disappears once a saturation scale develops, the dilute result is specific to the abelian limit, and if it persists, the Milne-frame caveat applies more broadly.
  • The 1/m domain-size law is a falsifiable prediction of the dilute limit: varying m should move the transverse correlation width of ε_LRF linearly, whereas a saturation-scale picture would make the width insensitive to m at fixed g^2 μ.
  • The same backward-light-cone integral machinery could be applied to other observables, such as the Glasma's chromo-electric and chromo-magnetic fields separately, to see whether the coordinate artifact shows up in field-strength components or only in the energy-momentum tensor.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. This proceedings paper presents a (3+1)-dimensional dilute Glasma computation within a generalized McLerran-Venugopalan model with tunable longitudinal and transverse correlation scales. The field strength tensor is computed to leading order in the color-charge expansion (m=n=1), yielding an effectively abelian Glasma field, and the energy-momentum tensor and local rest-frame energy density are evaluated event-by-event on a Milne grid. The main reported findings are: (i) the rapidity profile of T^{ττ} grows substantially at large spacetime rapidity η_s while ε_LRF and T^{xx}+T^{yy} do not; this enhancement is attributed to a coordinate effect in the Milne frame, leading the authors to argue that Milne coordinates are impractical for initial conditions with extended longitudinal structure; and (ii) the transverse size of energy-density domains is set by the infrared regulator m, since no saturation scale is generated in the dilute limit. The paper summarizes and interprets results previously derived in [5,6].

Significance. If the central interpretive claim is correct, the paper has practical implications for the choice of coordinate frame in (3+1)-dimensional heavy-ion initial conditions and early-time hydrodynamics: standard Milne-coordinate treatments may misrepresent the dynamics when the initial state has finite longitudinal structure. The numerical event-by-event framework and the use of a tunable nuclear model are useful, and the underlying leading-order field-strength computation is already published and internally consistent. The paper is honest about the dilute approximation and the absence of a dynamically generated saturation scale. However, the distinct new claim that the T^{ττ} enhancement is a coordinate artifact and that Milne coordinates are unsuitable is asserted rather than demonstrated; this is the load-bearing point that needs strengthening.

major comments (3)
  1. [Sec. 4 / Fig. 2 / Conclusion] The central claim that the large-η_s enhancement of T^{ττ} is a coordinate effect and that Milne coordinates are 'too impractical' is not derived. The text states: 'We explain this enhancement as a coordinate effect arising from large longitudinal pressures in the co-moving Bjorken frame. These pressures lead to growing T^{ττ} because T^{μν} is traceless in any frame.' Tracelessness alone does not imply this behavior; one must show how the local rest-frame quantities combine, e.g. T^{ττ} = (ε+p_L) cosh^2(η_s - y_flow) - p_L (with appropriate transverse averaging), and verify from the computed field configuration that y_flow deviates from η_s at large rapidity. Without such a derivation, the conclusion rests on comparing a coordinate component with a scalar (ε_LRF), which is not by itself evidence of practical inadequacy. This is the main novel point beyond [5,6] and needs either a deriva
  2. [Sec. 4 / Fig. 2] The rapidity profiles are averaged over only 10 central collision events, and no statistical uncertainty (standard error, event-by-event spread, or centrality bin) is shown. Since the T^{ττ}/ε_LRF separation is the primary effect, the reader cannot judge whether the enhancement is robust or driven by rare configurations. Please include error bars or at least a statement of the event-by-event variance, and specify the exact centrality/impact-parameter averaging procedure.
  3. [Sec. 4 / Eq. (4) / Conclusion] The statement that 'the size of the transverse domains is governed by 1/m' is a structural consequence of introducing m as the IR regulator in Eq. (4) and working in the dilute abelian limit where no other transverse scale exists; the paper itself notes that g^2 μ is only an overall prefactor. This is not a dynamical prediction. The concluding phrase characterizing 1/m as 'closely related to a saturation scale' is potentially misleading, since a saturation scale is dynamically generated while m is an input parameter. Please reframe this as a model property and avoid implying a dynamical relation that the dilute approximation explicitly excludes.
minor comments (5)
  1. [Eq. (1)] The square-root symbols over T(x^±, x) and T(y^±, y) appear to be missing or mis-typeset in the current text; the expression should read √T(x^±, x) √T(y^±, y), not 'pT'.
  2. [Eq. (3)] The formula for U_ξ is hard to parse as printed: the factors 1/(√(2π) ξ^2) and the exponents lose their divisions. Please typeset the equation with unambiguous parentheses and division bars.
  3. [Sec. 3] The symbol γ is defined as √s_NN/(2m_u) using 'm_u' for the nucleon mass, which conflicts notationally with the MV parameter μ used in Eq. (1). Use a distinct symbol, e.g. m_N.
  4. [Sec. 4 / Fig. 2 caption] The caption describes the T^{ττ} enhancement as 'unphysical'; if the point is that it is a coordinate artifact, the enhancement is a real component of the tensor in Milne coordinates, not an unphysical value. Rephrasing to 'coordinate-dependent' would be more precise.
  5. [General] The paper states that the f^{μν} integrals were presented in [6] but does not reproduce them. For a standalone proceedings this is acceptable, but a brief reminder of the integral form or a pointer to the equation number in [6] would improve readability.

Circularity Check

1 steps flagged · score 6.0 of 10

The 1/m transverse-domain result is a restatement of the IR regulator introduced in Eq. (4); the Milne-frame discussion is an underived interpretation, not a circular reduction.

  1. self definitional [Sec. 3, Eq. (4); Sec. 4, Fig. 1 caption and paragraph after Fig. 1]
    "The IR regulator m = 0.2 GeV sets the size of the fluctuating domains to∼ 1/m. ... In contrast to Glauber, MC-KLN and IP-Glasma results [10], the size of the transverse domains is governed by 1/m, as it was the case for the single-nucleus fields. This is because in the dilute approximation, there does not develop a saturation scale Qs∝ g^2 μ as it does for non-perturbative calculations. Instead, the IR regulator introduced in the nuclear model remains as the only transverse scale in the dilute Glasma and g^2 μ becomes an overall prefactor."

    The single-nucleus field is defined in Eq. (4) as A^∓,a = ∫ d^2k ρ̃^a(x±,k)/(k^2+m^2) e^{-k^2/(2Λ^2_UV)} e^{-ik·x}; the denominator introduces m as the transverse correlation scale by construction. The paper then reports as a result that transverse Glasma domains are sized by 1/m, acknowledging "as it was the case for the single-nucleus fields" and that m "remains as the only transverse scale." At leading order m=n=1 the field is abelian and no new scale such as Q_s can appear, so the outcome is fixed by the input propagator. The claim is therefore a restatement of the model input, not a prediction independent of it.

full rationale

Only one load-bearing step exhibits a genuine by-construction reduction: the transverse domain size 1/m is put in via the IR regulator in Eq. (4) and then returned as a finding. The other headline claim—that the large-|η_s| growth of T^{ττ} is a coordinate effect that makes Milne coordinates impractical—is not circular in the formal sense: T^{ττ} is computed from the field-strength tensor, and the interpretation requires additional physical reasoning (e.g., transforming to the local rest frame), which the paper asserts rather than derives. That is a justification weakness, not a circular reduction. The dilute-Glasma expressions are cited from prior work by the same group [5–7], but those are published external references and the present paper does not need a self-citation to force its conclusion; the 1/m result would follow from the model equations even without the citation. Self-citation is present but not the vehicle of the circularity. Given that one of the two central results reduces by construction, the circularity score is 6 (partial circularity).

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

The central results rest on a dilute (leading-order) expansion with a Gaussian MV-model input, two tunable regulators (m, xi) that set the very scales the paper reports, and a standard Landau-frame definition. No new entities are introduced.

free parameters (4)
  • IR regulator m = m = 0.2 GeV (Fig. 1); m = 0.2 mu and 2.0 mu (Fig. 2)
    Introduced in Eq. (4) as the infrared regulator that sets transverse domain size; the paper's central statement that domain size is ~1/m is a direct consequence of this choice, making it a load-bearing free parameter.
  • Longitudinal correlation scale xi = xi = 0.1 R_l and 2.0 R_l (Fig. 2)
    Introduced in Eq. (3) to tune longitudinal correlations; the rapidity profiles shown depend on this parameter, so the observed longitudinal structure is controlled by the input, not discovered.
  • MV parameter mu = not specified; sets g^2 mu^2 in Eq. (1)
    Overall strength of color charge fluctuations; the paper notes g^2 mu becomes an overall prefactor in the dilute Glasma, so it does not set shape but is a free input scale.
  • UV regulator Lambda_UV = Lambda_UV = 10 GeV
    Numerical regulator chosen in Eq. (4) to match lattice spacing; not a physical prediction, but an arbitrary choice.
assumptions (5)
  • domain assumption The nuclear color currents are recoilless and frozen in light-cone time: J_mu^{A/B} = delta^{mu pm} rho_{A/B}, with rho depending on x+- and transverse coordinate.
    Statement in Sec. 2; standard CGC approximation, but it fixes the kinematics and is needed for the wave-equation setup.
  • domain assumption Color charges are sampled from a Gaussian functional with correlation function Eq. (1) determined by g^2 mu^2 and the envelope T.
    Sec. 3; the MV model generalization is assumed, not derived; all event-by-event results inherit this distribution.
  • domain assumption The dilute expansion in powers of rho_A^m rho_B^n is truncated at leading order m=n=1, so the Glasma field strength is effectively abelian: f^mu nu = partial^mu a^nu - partial^nu a^mu.
    Sec. 2; this is the load-bearing approximation that removes nonlinear saturation effects and makes the field strength expressions tractable.
  • domain assumption The local rest frame energy density is defined by the Landau condition T^mu_nu u^nu = epsilon_LRF u^mu.
    Eq. (5); standard definition, but needed to interpret the rapidity profiles and identify the Milne-frame artifact.
  • domain assumption The longitudinal extent of the nucleus is described by the stretched Woods-Saxon envelope T(x+- , x) with gamma = sqrt(s_NN)/(2 m_u), Eq. (2).
    A model choice for how the nuclear thickness extends in x+-; it sets the physical rapidity range and the scale R_l used for xi.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The (3+1)D structure of the dilute Glasma." pith.science (2026). https://pith.science/paper/YFVJT5M2

@misc{pith2026250820725,
  author       = {Pith},
  title        = {Pith review of: The (3+1)D structure of the dilute Glasma},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YFVJT5M2}},
  note         = {Machine review of arXiv:2508.20725}
}
read the original abstract

We study the (3+1)D structure of the Glasma in the dilute approximation, which allows us to describe the longitudinal dynamics that arise from the three-dimensional nuclear structure. We employ a nuclear model with tunable longitudinal and transverse fluctuation scales that generalizes the McLerran-Venugopalan model. We discuss the longitudinal profiles of the energy-momentum tensor and the transverse structure of the local rest frame energy density.

Figures

Figures reproduced from arXiv: 2508.20725 by the authors.

Figure 1
Figure 1. Local rest frame energy density ϵLRF in the transverse x = (x, y) plane at mid-rapidity ηs = 0 and proper time τ = 0.4 fm/c for a single event with impact parameter equal to R and √ sNN = 200 GeV corresponding to Au+Au collisions at RHIC. The IR regulator m = 0.2 GeV sets the size of the fluctuating domains to ∼ 1/m. (upper signs) and right-moving B (lower signs). The color charge distributions ρA/B(x ± , x) of the … view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

12 extracted references · 2 canonical work pages

  1. [1]

    Acharya et al

    S. Acharya et al. (ALICE), The ALICE experiment: a journey through QCD , Eur. Phys. J. C 84, 813 (2024), 2211.04384 . 10.1140/epjc/s10052-024-12935-y

    arXiv 2024

  2. [2]

    Gelis, E

    F. Gelis, E. Iancu, J. Jalilian-Marian, R. Venugopalan, The Color Glass Condensate , Ann. Rev. Nucl. Part. Sci. 60, 463 (2010), 1002.0333 . 10.1146/annurev.nucl.010909.083629

    arXiv 2010

  3. [3]

    Gelis, Color Glass Condensate and Glasma , Int

    F. Gelis, Color Glass Condensate and Glasma , Int. J. Mod. Phys. A 28, 1330001 (2013), 1211.3327 . 10.1142/S0217751X13300019

    arXiv 2013

  4. [4]

    Lappi, L

    T. Lappi, L. McLerran, Some features of the glasma , Nucl. Phys. A 772, 200 (2006), hep-ph/0602189 . 10.1016/j.nuclphysa.2006.04.001

    arXiv 2006

  5. [5]

    Ipp, D.I

    A. Ipp, D.I. M \"u ller, S. Schlichting, P. Singh, Spacetime structure of (3+1)D color fields in high energy nuclear collisions , Phys. Rev. D 104, 114040 (2021), 2109.05028 . 10.1103/PhysRevD.104.114040

    arXiv 2021

  6. [6]

    A. Ipp, M. Leuthner, D.I. M \"u ller, S. Schlichting, K. Schmidt, P. Singh, Energy-momentum tensor of the dilute (3+1)D glasma , Phys. Rev. D 109, 094040 (2024), 2401.10320 . 10.1103/PhysRevD.109.094040

    arXiv 2024

  7. [7]

    Aspects of the dilute Glasma

    M. Leuthner, Aspects of the dilute Glasma , Ph.D. thesis, Vienna, Tech. U. (2025), 2501.16216 . 10.34726/hss.2025.99382

    work page Pith review arXiv 2025

  8. [8]

    McLerran, R

    L.D. McLerran, R. Venugopalan, Gluon distribution functions for very large nuclei at small transverse momentum , Phys. Rev. D 49, 3352 (1994), hep-ph/9311205 . 10.1103/PhysRevD.49.3352

    arXiv 1994

Show all 12 references
  1. [9]

    McLerran, R

    L.D. McLerran, R. Venugopalan, Computing quark and gluon distribution functions for very large nuclei , Phys. Rev. D 49, 2233 (1994), hep-ph/9309289 . 10.1103/PhysRevD.49.2233

    1994 arXiv

  2. [10]

    Schenke, P

    B. Schenke, P. Tribedy, R. Venugopalan, Fluctuating Glasma initial conditions and flow in heavy ion collisions , Phys. Rev. Lett. 108, 252301 (2012), 1202.6646 . 10.1103/PhysRevLett.108.252301

    2012 arXiv

  3. [11]

    , " * 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 IN...

  4. [12]

    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 5, 2026 · model on record in the stance chip above.