Pith. sign in

REVIEW 3 major objections 6 minor 28 references

Elliptic and triangular flows in dAu collisions at 200 GeV in the fusing color string model

T0 review · 3 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read In the fusing color string model, path-length-dependent quenching inside fused strings reproduces the measured elliptic flow of p-Au and d-Au collisions at 200 GeV and the observed ordering v2(d-Au)>v2(p-Au), while triangular flow comes…

desk verdict A modest but honest extension of the fusing string model to small-system flow; the v2 ordering is plausible, but the result leans entirely on an unproven QED quenching law. read the letter →

arxiv 1909.02131 v1 pith:XZXOLPBP submitted 2019-09-04 hep-ph nucl-th

M.A. Braun , C. Pajares This is my paper · ORCID
classification hep-phnucl-th
keywords colorstringmodelfusionpartonquenchingellipticflowtriangularsmallcollisionsystemsp-Aucollisionsd-Au
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 claims that the elliptic flow measured in deuteron-gold and proton-gold collisions at 200 GeV can be explained entirely by the fusing color string model, without invoking hydrodynamic expansion of a quark-gluon plasma. The mechanism is geometric: a produced particle loses energy as it crosses the gluonic fields of the strings it passes through, and because the string distribution in each event is azimuthally anisotropic, this path-length-dependent quenching generates a nonzero v2. With a single quenching parameter κ fixed by mid-central Pb-Pb data, the model reproduces the PHENIX v2(pT) points for central d-Au and p-Au collisions and the ordering v2(d-Au) > v2(p-Au). The triangular flow v3 comes out larger than the data, in the same way as models based on Color Glass Condensate initial conditions with hydrodynamic evolution, which the authors trace to their simplified parton-hadron duality treatment of hadronization.

What carries the argument

The central machinery is the fusing color string model with percolation, in which each string is a droplet of gluonic field of finite transverse size; the load-bearing formula is the QED-inspired quenching law p0(p,l) = p(1+$κp^{{-1/3}}$$T^{{2/3}}$l)^3, which converts path length l inside strings into an enhanced initial momentum that suppresses emission along directions with more string matter. The Monte Carlo implementation, with the Hulthen deuteron wave function, Gaussian nucleon density, string fusion, and multiplicity-based centrality selection, converts this geometry into the flow coefficients vn via event-by-event Fourier decomposition.

What would settle it

A measurement that breaks the predicted pattern would falsify the model: for instance, if PHENIX or STAR data at 200 GeV showed v2(p-Au) ≥ v2(d-Au) for central collisions, or if a first-principles QCD calculation showed that the Nikishov QED quenching formula is not applicable to partons in color strings, the mechanism would be ruled out.

Watch

Extended reading notes

Core claim

In the color string picture with fusion and percolation, azimuthal anisotropies arise because a produced parton must traverse the gluonic fields of the strings it crosses, and its initial transverse momentum is enhanced by the quenching factor p0(p,l) = p(1+$κp^{{-1/3}}$$T^{{2/3}}$l)^3, where l is the path length inside each string and T its tension. Running a Monte Carlo that places interacting nucleons, distributes strings, fuses overlapping ones, and computes the quenching factor for each emission direction, the authors reproduce the experimentally observed v2(pT) for central d-Au and p-Au collisions at 200 GeV and the ordering v2(d-Au)>v2(p-Au), with κ=0.6 fixed by mid-central Pb-Pb data. The triangular flow v3(pT) overshoots the PHENIX data, in the same way as initial-state CGC plus hydrodynamics models. The authors conclude that the number of emitting sources in d-Au being roughly twice that in p-Au does not suppress v2, because the strings communicate through the common gluonic field.

Load-bearing premise

The result rests on applying the QED formula for energy loss of a charged particle in an external electromagnetic field, Eq. (7), to partons crossing fusing color strings in QCD; the paper does not derive this quenching from QCD, so if that transfer fails the predicted v2, v3 and their ordering do not follow.

Editorial extensions

If this is right

  • If the model is right, the observed v2 in p-Au and d-Au collisions does not require a hydrodynamic description of a quark-gluon plasma; a two-stage string-emission-with-quenching scenario suffices.
  • The ordering v2(d-Au)>v2(p-Au) for central collisions is natural in this picture despite d-Au having roughly twice as many sources, because strings in the overlap region communicate through a common gluonic field rather than remaining independent emitters.
  • The predicted decrease of the d-Au versus p-Au flow difference with centrality, and its near-vanishing in peripheral collisions, is a testable consequence of the geometry.
  • Because v3 overshoots the data in both this model and CGC-plus-hydro approaches, the discrepancy likely lies in the hadronization or final-state stage rather than the initial-state dynamics, pointing to where the models should be improved.

Reading between the lines

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

  • A natural testable extension would be to include a more detailed hadronization stage, such as fluctuations in parton-hadron conversion, and check whether v3 drops toward the data while v2 remains stable; if so, the v3 overshoot would be diagnosed as a final-state effect common to all such models.
  • The same geometric-quenching mechanism should produce predictions for other small systems such as He-Au or O-O collisions at RHIC and LHC energies, and the model's scaling with ε p^{2/3} T^{1/3} l could be checked across these systems.
  • The paper's reliance on a QED quenching formula suggests a possible weakness: if the formula is not valid in QCD, the apparent success of the model might be accidental; deriving the equivalent quenching law from perturbative QCD or from string dynamics would either firm up the mechanism or reveal its limits.
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 / 6 minor

Summary. The paper applies the fusing color string model to p-Au and d-Au collisions at 200 GeV, computing elliptic flow v2 and triangular flow v3 as functions of transverse momentum. The model represents each nucleon-nucleon interaction by color strings that fuse when overlapping, and generates azimuthal anisotropy through a path-length-dependent quenching of the final particle momentum, Eq. (7). The quenching coefficient kappa is fixed to 0.6 by fitting the integrated v2 in mid-central Pb-Pb collisions. The authors compare their results with PHENIX central-collision data, report that the observed ordering v_n(d-Au) > v_n(p-Au) is reproduced, that v2 agrees satisfactorily with data, and that v3 overshoots the measured values. They also present predictions for mid-central and peripheral collisions.

Significance. If the results are robust, the paper offers a simple, non-hydrodynamic explanation of the observed small-system flow hierarchy, with only one tunable parameter (kappa). The central comparison is not circular: kappa is fixed to Pb-Pb data, while the p-Au and d-Au v2 and v3 values are predictions. The paper also honestly notes the v3 overshoot, which it shares with CGC-plus-hydrodynamics approaches. Strengths include a concrete falsifiable prediction for the centrality dependence of the d-Au versus p-Au difference and a transparent physical mechanism. The main limitations are the largely undefended transfer of the QED quenching formula to QCD strings and the incomplete specification of the Monte Carlo implementation, both of which affect the confidence one can place in the quantitative comparison.

major comments (3)
  1. [Section 2, Eq. (7)] The entire anisotropy mechanism rests on Eq. (7), p0(p,l) = p (1 + kappa p^{-1/3} T^{2/3} l)^3, which is taken from the QED treatment of a charged particle in an external electromagnetic field [24]. The paper does not derive this law for a parton traversing fused color strings, does not specify how the QED formulas map to QCD, and does not discuss non-Abelian corrections or the regime of validity. Since the functional form and the p^{-1/3} T^{2/3} l scaling determine the pT dependence of v2 and the magnitude of v3, a different energy-loss law (linear, BDMPS-Z, or different power of p) would change the predictions. Fitting kappa to one Pb-Pb value fixes the overall normalization but does not validate the functional form. The authors should either provide a derivation or a clear phenomenological argument for transferring Eq. (7) to color strings, and should quantify how sensitive the reported v2 and v3 are to the assumed functional form.
  2. [Section 3, Monte Carlo procedure] The simulation details are not fully specified: the number of exchanged strings is described only as taken from previous calculations, the distribution of string positions and the fusion algorithm are described qualitatively, and the statistical error is quoted as 'around 5%' without stating the number of events, the binning, or how the error was estimated. The paper should provide enough information for the calculation to be reproduced, or release the code/data. Without this, the reader cannot judge whether the reported agreement with data and the d-Au/p-Au ordering are stable with respect to the Monte Carlo implementation.
  3. [Section 3, centrality selection] The centrality classes are defined by multiplicity windows at fixed impact parameter b with 0.9 mu_max < mu < mu_max for central collisions, but the comparison is made to experimental 0-5% centrality data. The text does not demonstrate that this multiplicity cut actually corresponds to the experimental 0-5% centrality bin, nor how the impact parameter is sampled or averaged in the Monte Carlo. Because v2 and v3 are known to be centrality-dependent, a mismatch in the centrality definition could shift the curves in Fig. 1 relative to data. The authors should clarify the relationship between their multiplicity windows and the experimental centrality classes.
minor comments (6)
  1. [Section 2, Eq. (1)] Equation (1) appears to contain an extra parenthesis: C(phi) = A(1 + (1 + 2 sum ...)) has an unbalanced bracket; it should likely read C(phi) = A (1 + 2 sum_{n>=1} v_n cos(n phi)).
  2. [Abstract and Section 1] The abstract states the paper studies flows 'for p-Au and d-Au collisions', which is clear, but the phrase 'at 5-13 TeV GeV' in Section 2 is an awkward typo; the energy should be 5.02 TeV or similar, and 'GeV' is redundant.
  3. [Section 2] The text says 'anisotropy od string distribution'; this should be 'of the string distribution'.
  4. [Section 2] The scaling property attributed to Eq. (7) is stated without showing the derivation or a reference to where it is derived; the authors should add a short explanation of how the product epsilon p^{2/3} T^{1/3} l arises.
  5. [Section 3] In Figure 1 the caption says 'Experimental data ... are from [13]' but the data points for p-Au and d-Au are not distinguished in the caption; the figure legend should state which points correspond to which system.
  6. [References] A few references contain typographical errors, e.g., [6] '1804.029442' appears to have an extra digit, and [13] 'ncl-ex' should be 'nucl-ex'. These should be corrected.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the p-Au/d-Au v2 and v3 predictions are external to the fitted parameter and benchmarked against PHENIX data.

full rationale

The derivation chain is not circular in the sense defined by the protocol. The only fitted quantity is the quenching coefficient kappa = 0.6, adjusted to the integrated mid-central Pb-Pb v2. The paper then computes v2(pT), v3(pT), and the d-Au/p-Au ordering for 200 GeV p-Au and d-Au collisions from the string-fusion geometry and the path-length-dependent quenching formula (Eq. (7)), comparing these outputs with external PHENIX measurements [13]. None of the compared small-system observables is used to fix kappa or any other parameter, so the central predictions are not equal to inputs by construction. The quenching law Eq. (7) is an explicitly stated modeling assumption borrowed from QED [24] and the authors' earlier string-model application [28]; whether that transfer is physically justified is a correctness question, not a circularity one. The routine self-citations ([20], [25], [26], [28]) support the model framework and previous scaling tests, but they do not define the small-system flow coefficients in terms of the fitted kappa. The v3 overshoot is a genuine failed prediction, further indicating that the outputs are not constructed to match the data. Therefore no circular step of the specified kinds is present.

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

The model introduces no new entities beyond the existing color string, fusion, and percolation framework. The central claim rests on calibrated inputs: the quenching coefficient kappa, the inherited string numbers, and the manually chosen centrality windows. These are all stated in the text, though the string numbers are not quantified.

free parameters (3)
  • kappa = 0.6
    Universal quenching coefficient adjusted to reproduce the experimental integrated v2 in mid-central Pb-Pb collisions, then applied to p-Au and d-Au. This calibration is stated explicitly in Section 2.
  • number of exchanged strings = not specified in this paper
    The number of strings per interacting nucleon pair is taken from previous calculations [20]. This sets the source density and therefore the magnitude of the quenching effect and the flow coefficients.
  • centrality multiplicity windows = 0.9 mu_max, 0.45-0.55 mu_max, below 0.1 mu_max
    Centrality classes are defined by arbitrary multiplicity windows at fixed impact parameter. Changing these thresholds changes which events are averaged and thus the predicted flow values.
assumptions (6)
  • domain assumption The QED-inspired quenching formula Eq. (7) applies to partons traversing fusing color strings.
    This is the central mechanism producing flow. The formula is borrowed from [24] and not re-derived for QCD strings.
  • domain assumption String decay follows the Schwinger mechanism, with fluctuating string tension producing a thermal transverse momentum distribution.
    Used in Eqs. (5) and (6) to set the initial parton transverse momentum spectrum before quenching.
  • domain assumption The color string fusion and percolation picture, including the number of strings and their tensions, is taken from the authors' earlier review [20].
    The entire simulation is built on this framework, which is not re-derived here.
  • domain assumption Parton-hadron duality holds for the final particle distribution.
    Mentioned in Section 3 as the simplified treatment of hadronization; the authors suggest this may be the reason v3 is overpredicted.
  • ad hoc to paper Multiplicity windows at fixed impact parameter define centrality classes.
    The 10 percent, 45-55 percent, and sub-10 percent windows of mu_max are chosen by hand to emulate experimental central, mid-central, and peripheral collisions.
  • domain assumption The deuteron wavefunction is the Hulthen wavefunction.
    Standard input for distributing the two nucleons in the deuteron projectile.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Elliptic and triangular flows in dAu collisions at 200 GeV in the fusing color string model." pith.science (2026). https://pith.science/paper/XZXOLPBP

@misc{pith2026190902131,
  author       = {Pith},
  title        = {Pith review of: Elliptic and triangular flows in dAu collisions at 200 GeV in the fusing color string model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XZXOLPBP}},
  note         = {Machine review of arXiv:1909.02131}
}
abstract

In the color string picture with fusion and percolation the elliptic and triangular flows are studied for p-Au and d-Au collisions at 200 GeV. The ordering $v_n(d-Au)>v_n(p-Au)$ observed experimentally for central collisions is reproduced.The calculated elliptic flow $v_2$ at central collisions agrees satisfactorily with thedata. The triangular flow $v_3$ is found to be greater than the experimental values, similar to the resultsobtained in the approach based on the Color Glass Condensate initial conditions with subsequenthydrodynamical evolution.

Figures

Figures reproduced from arXiv: 1909.02131 by the authors.

Figure 1
Figure 1. The calculated flow coefficients v2 (left panel) and v3 (right panel) as function of transverse momenta pT for d-Au (upper curves) and p-Au central collisions at 200 GeV. Experimental data for d − Au (upper points) and p − Au collisions at 0-5% centrality are from [13]. events in which µ lies within a certain part of the total interval µ < µmax where µmax is the largest multiplicity for all b and runs. In such an ap… view at source ↗
Figure 2
Figure 2. The calculated flow coefficients v2 (left panel) and v3 (right panel) as function of transverse momenta pT for d-Au (upper curves) and p-Au mid-central collisions at 200 GeV. 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0 0.5 1 1.5 2 2.5 3 v2 pT GeV/c peripheral ’paupd.res’ u 1:2 ’paupp.res’ 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0 0.5 1 1.5 2 2.5 3 v3 pT GeV/c peripheral ’paupd.res’ u 1:3 ’paupp.res’ u 1:3 [PITH_FULL_IM… view at source ↗
Figure 3
Figure 3. The calculated flow coefficients v2 (left panel) and v3 (right panel) as function of transverse momenta pT . for d-Au and p-Au peripheral collisions at 200 GeV. For d-Au collisions v2 corresponds to the lower curve and v3 corresponds to the upper curve. 6 [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

28 extracted references · 20 canonical work pages

  1. [24]

    A.I.Nikishov, Nucl. Phys. B 21 (1970) 346

    work page 1970

  2. [1]

    C 80 (2009) 024909 [nucl-ex/ 0905.1070]

    S.Afanasiev et al , PHENIX collab., Phys.Rev. C 80 (2009) 024909 [nucl-ex/ 0905.1070]

    arXiv 2009

  3. [2]

    R.Aamodt et al , ALICE collab., Phys. Rev. Lett. 107 (2011) 032301 [nucl.ex/1105.3865]

    arXiv 2011

  4. [3]

    A.Adare et al , PHENIX collab., Phys. Rev. Lett. 107 (2011) 252301 [nucl- ex/1105.3928]

    arXiv 2011

  5. [4]

    B.Alver et al PHOBOS collab,Phys. REv. Lett. 98 (2007) 242302 [nucl-ex/061003 7]

    work page 2007

  6. [5]

    L.Adamczyk et al (STAR collab.), Phys. Rev. Lett. 115 (2015) 22 2301 [nucl- ex/1505.07812]

    arXiv 2015

  7. [6]

    S Acharya et al (ALICE collab.) JHEP 1807 (2018) 103 [nucl-ex/18 04.029442]

    work page 2018

  8. [7]

    JHEP 09 (2010) 091

    V.Khachatryan et al. JHEP 09 (2010) 091. 7

    work page 2010

Show all 28 references
  1. [8]

    Phys.Lett

    S.Chatrchyan et al. Phys.Lett. B 718 (2013) 795

    2013

  2. [9]

    Phys.Lett

    B.Abelev et al. Phys.Lett. B 719 (2013) 29

    2013

  3. [10]

    G.Aad et al. Phys. Rev. Lett. 110 (2013) 182302

    2013

  4. [11]

    A.Adare et al. Phys. Rev. Lett. 111 (2013) 212301

    2013

  5. [12]

    A.Adare et al. Phys. Rev. Lett. 114 (2015) 192301

    2015

  6. [13]

    PHENIX collab

    C.Aidala et al. PHENIX collab. Nature Physics 15 (2019) 214 [ncl-ex/1805.02973]

    2019

  7. [14]

    Bozek,Phys

    P. Bozek,Phys. Rev. C 85 (2012) 014911 [hep-ph/1112.0915]

    2012 arXiv

  8. [15]

    M.Habich, J.L.Nagle, P.Romatschke Eur. Phys. J C 75 (2015) 15

    2015

  9. [16]

    C.Shen,J-F.Paquet,G.S.Denicol, S,Jeon, C.Gale, Phys. Rev. C95 (2017) 014906

    2017

  10. [17]

    M.Mace, V.Skokov,P.Tribedy, R.Venugopalan, Phys. Lett. 788 (2019) 161 [hep- ph/1807.00825]

    2019 arXiv

  11. [18]

    M.Mace, V.Skokov,P.Tribedy, R.Venugopalan, Phys. Rev. Lett 121 (2018) 052301 [hep/ph/1805.09342]; erratum ibid 123 (2019) 039901

    2018 arXiv

  12. [19]

    B.Schenke, C,Shen, P.Tribedy, arXiv:1908.06212[nucl-th]

    1908 arXiv

  13. [20]

    Rep., 599 (2015) 1-50

    M.A.Braun, J.Dias de Deus, A.S.Hirsch, C.Pajares, R.P.Sharenberg , B.K.Srivastava, Phys. Rep., 599 (2015) 1-50

    2015

  14. [21]

    J.Dias de Deus and C.Pajares, Phys. Lett. B 695 (2011) 211

    2011

  15. [22]

    A.Bialas, Phys. Lett. B 466 (1999) 301

    1999

  16. [23]

    J.Dias de Deus and C.Pajares, Phys. Lett. B 642 (2006) 455

    2006

  17. [25]

    C.Andres,J.Dias de Deus, A.Moscoso,C.Pajares, C.Salgado, Phys. Rev. C 92 (2015) 03496

    2015

  18. [26]

    C.Andres, M.A.Braun, C.Pajares, Eur. Phys. J. A 53 (2017) 41

    2017

  19. [27]

    PHENIX collab

    A.Adare et al. PHENIX collab. Phys. Rev. C 90 (2014) 034902

    2014

  20. [28]

    M.A.Braun, C.Pajares, V.V.Vechernin, Nucl. Phys. A 906 (2013) 14. 8

    2013

Pith tools

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