Pith. sign in

REVIEW 4 major objections 6 minor 37 references

Existence of dynamical fluctuation in AMPT generated data for Au+Au collisions at 10 AGeV

T0 review · 4 major / 6 minor · reviewed 2026-08-10 · deepseek-v4-flash

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

Pith's one-line read Intermittent, self-similar density fluctuations appear in both modes of AMPT-generated 10 AGeV Au+Au collisions, with string melting producing the stronger high-order signal.

desk verdict Useful low-energy AMPT intermittency point, but the missing χ-transformation equation and lack of controls keep it from being fully convincing as a dynamical signal. read the letter →

arxiv 2501.05175 v1 pith:REPO2Z3K submitted 2025-01-09 hep-ph

Somen Gope , Supriya Das , Saikat Biswas This is my paper · ORCID
classification hep-ph PACS 25.75.-q
keywords ScaledfactorialmomentIntermittencyAMPTmodelStringmeltingmechanismHeavy-ioncollisionsDynamicalfluctuationsCriticalpoint
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 tries to establish that particle production in AMPT simulations of 10 AGeV Au+Au collisions contains genuine dynamical, self-similar fluctuations, or intermittency, rather than only statistical noise. Using scaled factorial moments in a two-dimensional momentum-space coordinate system, it finds power-law growth of the moments with the number of bins in both the default and string-melting modes of AMPT. The string-melting mode shows stronger intermittency at higher moment orders, which the authors connect to a more fully partonic final state. The result matters because intermittent fluctuations are a proposed observable for locating a critical endpoint in the QCD phase diagram at FAIR/SIS100 energies, so confirming that the observable works in a transport model at 10 AGeV is a step toward that search.

What carries the argument

The Scaled Factorial Moment (SFM) of order $q$, $\langle F_q\rangle$, is computed from bin counts $K_m$ in a partition of phase space; the factorial structure removes Poisson statistical fluctuations and isolates dynamical density fluctuations. The power-law dependence $\langle F_q\rangle \propto M^{\alpha_q}$ on the number of bins $M$ is the intermittency signature, with $\alpha_q$ the intermittency index and $d_q = \alpha_q/(q-1)$ the anomalous fractal dimension. Before forming moments, the analysis maps each particle's $(p_x,p_y)$ into cumulant coordinates $\chi(p_x),\chi(p_y)$ in $[0,1]$ to flatten the single-particle density, and it averages using the horizontal averaging form (Eq. 5).

What would settle it

Recomputing $\langle F_q\rangle$ in raw $(p_x,p_y)$ bins without the $\chi$ transformation, or on synthetic Poisson-random events with the same single-particle density, should make the power-law slopes disappear if the reported intermittency is truly dynamical; if the slopes persist, the paper's central claim would be falsified.

Watch

Extended reading notes

Core claim

In 10 AGeV Au+Au collisions simulated with AMPT, the horizontally averaged scaled factorial moments in two-dimensional $\chi(p_x)$--$\chi(p_y)$ space follow $\langle F_q\rangle \propto M^{\alpha_q}$ for moment orders $q=2$ through $6$, in both the default and string-melting modes. The paper interprets this power law as a sign of intermittency, meaning self-similar, fractal-like density fluctuations in the emitted charged particles. The intermittency indices $\alpha_q$ are larger in the string-melting mode, especially at $q=5$ and $q=6$ ($\alpha_6 \approx 0.081$ versus $0.047$), which the paper reads as stronger dynamical fluctuations arising from the more partonic dynamics of string melting. The anomalous fractal dimension $d_q = \alpha_q/(q-1)$ increases with $q$, which the authors link to multifractal, branching-like particle production, and they note that $\langle F_q\rangle$ values at 200 GeV are similar to those at 10 AGeV.

Load-bearing premise

The analysis assumes that the $\chi(p_x)$--$\chi(p_y)$ mapping, whose explicit formula is deferred to an earlier paper, fully flattens the single-particle density so that the measured factorial-moment slopes come from dynamical fluctuations rather than from the shape of the spectrum.

Editorial extensions

If this is right

  • AMPT default and string-melting modes both reproduce the expected self-similar fluctuation signature at 10 AGeV, so the SFM method can be applied to FAIR-energy experimental data with a known model baseline.
  • The larger high-order intermittency indices in string-melting mode make $\alpha_q$ at $q \geq 5$ a sensitive probe of partonic degrees of freedom in heavy-ion collisions.
  • The increasing $d_q$ with $q$ indicates non-uniform, multifractal fluctuations rather than a single fractal dimension, consistent with branching particle production.
  • The similar $\langle F_q\rangle$ behavior at 10 AGeV and 200 GeV suggests that within the AMPT model the intermittency pattern is roughly energy-independent over this range.

Reading between the lines

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

  • A natural testable extension is to repeat the analysis at intermediate beam energies (roughly 5 to 30 AGeV) and map $\alpha_q$ versus energy; if the paper's picture is right, the separation between string-melting and default modes should grow as the partonic phase becomes more dominant.
  • The explicit formula for the $\chi(p_x)$--$\chi(p_y)$ transformation is deferred to an earlier paper, so the result is not reproducible from this paper alone; supplying that formula and testing how the slopes respond to its parameters would settle whether the flattening is complete.
  • If the intermittency indices remain robust under alternative flattening maps, the high-order $\alpha_q$ difference between modes could serve as a model discriminator for whether a QGP-like partonic stage is formed at SIS100 energies.
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

4 major / 6 minor

Summary. The paper analyzes AMPT model-generated Au+Au collisions at 10 AGeV, in both default and string-melting modes, using the scaled factorial moment (SFM) technique in two-dimensional χ(px)−χ(py) space. The authors report a monotonic rise of ⟨F_q⟩ with M² for q = 2–6, fit intermittency indices α_q, and compute the anomalous fractal dimension d_q. They claim that both AMPT modes exhibit power-law behavior characteristic of intermittency, with stronger, more complex fluctuations in the string-melting mode at higher orders, and they interpret this as evidence relevant to the QCD critical-point search at FAIR energies. The paper also compares π+ transverse-mass spectra with E895 data and shows an analogous 200 GeV AMPT(default) result.

Significance. If the result is correct, it would provide model-based evidence that both hadronic and partonic scenarios encoded in AMPT produce self-similar density fluctuations at a FAIR-relevant beam energy, with a quantitative difference between the two modes at q = 5, 6. This would be a useful input for planning intermittency measurements in CBM. The authors follow the standard SFM recipe and present their figures clearly. The significance is, however, presently limited by the absence of the explicit χ-transformation equation, the lack of event/centrality/version/control details, and the unquantified error treatment; these omissions prevent the reader from verifying that the reported slopes measure dynamical fluctuations rather than single-particle spectral shape or binning artifacts.

major comments (4)
  1. [Section 4, 'These transformations are defined as follows'] The displayed equation for χ(px) is missing: the text reads 'Replace px by py in the following equation for χ(py)' and then no equation follows, with all parameters deferred to reference [6]. Since Eqs. (3)–(5) cancel Poisson noise only after the single-particle density is flattened, the whole extraction of α_q depends on this undocumented coordinate map. Provide the explicit χ transformation, its inverse/Jacobian if relevant, and a quantitative flatness test (e.g., χ²/ndf or KS statistic) for the distribution shown in Fig. 2.
  2. [Section 4, Table 1 and Figs. 3–4] The manuscript does not state the number of events, centrality selection, impact-parameter range, AMPT version, or model parameters used to generate the data, nor does it show a mixed-event or Poisson-pseudo-event baseline. Without such a baseline, the monotonic rise of ⟨F_q⟩ with ln M² could reflect residual density nonuniformity from the χ map, binning effects, or statistical artifacts rather than genuine dynamical intermittency. Add these reproducibility details and a control analysis; this is load-bearing for the claim that the slopes in Table 1 measure dynamical fluctuations.
  3. [Section 4, Table 1 and Figs. 3–4] The fitting procedure for the intermittency indices α_q is not described. The table reports slopes with statistical errors but no fit range, no χ²/ndf, and no statement of whether the fit is a weighted least-squares over all M² points or over a restricted region. Because the central claim '⟨F_q⟩ ∝ M^{α_q}' rests entirely on these slopes, specify the fit procedure and show the fitted lines on the ⟨F_q⟩ versus ln M² plots, so the reader can judge the quality of the power-law description.
  4. [Section 4, paragraph 'The errors in this analysis...'] The statement that the errors are 'only independent statistical errors not considered for different bin sizes' is asserted without quantitative support, and the text that excluding bin-size-dependent error 'does not change the outcomes appreciably' cites earlier work [35–37] but does not demonstrate this for the present data. Since SFM analyses are known to be sensitive to bin-size correlations, either include the bin-size-dependent errors or provide a numerical comparison showing that the extracted α_q are stable under their inclusion.
minor comments (6)
  1. [Section 2.1, Eq. (1)] The parameters a and b in the Lund fragmentation function are called 'free parameters' but their numerical values used in the AMPT runs are not given; please specify them (or state that the AMPT defaults were used).
  2. [Section 4, Fig. 1 and surrounding text] The text and Fig. 1 refer to '8 AGeV energy' while the rest of the paper analyzes 10 AGeV; clarify whether the E895 comparison is at 8 AGeV and, if so, state the centrality and rapidity acceptance used for the mT spectra.
  3. [Equations (3)–(5)] The normalization notation is inconsistent: Eq. (3) uses K, Eq. (4) uses ⟨K⟩, and Eq. (5) uses ⟨n⟩. Define each quantity and ensure the definitions are consistent across the equations.
  4. [Section 4, Fig. 7] The 200 GeV AMPT(default) comparison is presented without any event details, multiplicity ranges, or a corresponding table of α_q values, so the claimed similarity to the 10 AGeV results cannot be verified.
  5. [Section 4, Figs. 2–4 captions and text] There are minor typos and notation issues: 'cahrged' in the Fig. 6 caption, 'Px−Py' in the Summary, and 'χ(px−py)' in Section 4 should be 'χ(px)−χ(py)'.
  6. [References, [6]] Reference [6] is incomplete; please provide the full bibliographic details (volume, article number, page range) for 'S. Gope and B. Bhattacharjee, Eur. Phys. J. A (2021)'.

Circularity Check

1 steps flagged · score 2.0 of 10

No true circularity: the SFM slopes are measured observables, but the flattening map that turns them into dynamical-fluctuation claims is left unspecified and imported from a same-author citation.

  1. other [Section 4, paragraph following Fig. 1 (χ(px)–χ(py) transformation)]
    "To eliminate the initial shape dependence of the two-dimensional density distribution, the px and py values of each primary charged particle are converted into new variables, χ(px) and χ(py). These transformations are defined as follows: Replace px by py in the following equation for χ(py). ... Other parameters are defined in one of our earlier published papers [6]."

    The load-bearing premise is that the χ(px)–χ(py) coordinate map removes the single-particle spectral shape so that the fitted SFM slopes α_q in Table 1 measure dynamical fluctuations. That premise cannot be checked in this paper: the defining equation for χ(py) is omitted, the parameters are deferred to the authors' earlier paper [6], and the 'flatness' of Fig. 2 is offered visually with no quantitative test. This is an omitted-support/self-citation step in the derivation chain, not a fitted-input-called-prediction; the α_q values themselves are honestly reported fitted slopes. It does not make the central claim definitionally equivalent to its inputs, but it raises the burden on ref. [6] as the only articulated justification for interpreting the slopes as intermittency.

full rationale

The paper is an empirical SFM analysis of AMPT events, not a first-principles derivation. Eqs. (2)–(5) are the standard Bialas–Peschanski factorial-moment definitions, and the intermittency indices in Table 1 are the straight-line slopes of the <F_q> vs ln M^2 plots shown in Figs. 3–4. Reporting these fitted slopes as 'intermittency indices' is a definitional summary of the measurement, not a circular prediction. The central concern is the χ(px)–χ(py) flattening map: its formula is absent, its parameters are in the same-author reference [6], and its success is asserted from the visual flatness of Fig. 2 without a χ²/KS or bin-size test. This is a missing-support/self-citation weakness rather than a demonstrated equivalence between input and output, so it does not force a high circularity score. I therefore score 2: some self-citation in a load-bearing methodological step, but the central claim still has independent empirical content in the AMPT-generated data and the fitted moments.

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

The paper introduces no new entities. Its claims rest on the standard Bialas-Peschanski SFM framework, on the chi-transformation (defined only by reference to the authors' earlier paper [6]), and on the assumption that bin-size correlated errors can be ignored. The reported intermittency indices are fitted slopes, not independent inputs.

free parameters (1)
  • Intermittency index alpha_q (slopes of ln<F_q> vs ln M^2) = Default: 3.8e-3, 1.1e-2, 2.0e-2, 3.3e-2, 4.7e-2; SM: 3.1e-3, 9.6e-3, 2.1e-2, 4.4e-2, 8.1e-2
    Slopes from linear fits of ln<F_q> vs ln M^2 in Figs. 3 and 4; these are the paper's reported result, and the central claim that intermittency exists rests on them.
assumptions (4)
  • domain assumption Bialas-Peschanski scaling F_q proportional to M^alpha_q constitutes a signature of intermittency and self-similar particle production.
    Section 3, theoretical formulation; this is the interpretive bridge from a fitted power law to a physical claim about density fluctuations and possible critical behavior.
  • domain assumption The chi(px)-chi(py) transformation removes single-particle spectral shape dependence so that measured SFM fluctuations are dynamical.
    Section 4, paragraph on the transformation; the equation itself is not given, but the flat distribution in Fig. 2 is used as evidence.
  • domain assumption Ignoring bin-size correlated errors does not materially change the intermittency indices.
    Section 4, after Figs. 3-4; the paper relies on refs [35-37] rather than on an own estimate of the bin-correlated covariance.
  • domain assumption AMPT default and string melting modes bracket the partonic and hadronic physics relevant for low-energy heavy-ion collisions.
    Sections 2 and 5; the comparison of the two modes is meaningful only if string melting genuinely creates a partonic (QGP-like) phase and default does not.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Existence of dynamical fluctuation in AMPT generated data for Au+Au collisions at 10 AGeV." pith.science (2026). https://pith.science/paper/REPO2Z3K

@misc{pith2026250105175,
  author       = {Pith},
  title        = {Pith review of: Existence of dynamical fluctuation in AMPT generated data for Au+Au collisions at 10 AGeV},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/REPO2Z3K}},
  note         = {Machine review of arXiv:2501.05175}
}
abstract

In this study, the intermittency behavior of emitted particles produced in heavy ion collisions has been studied using both modes (default $\&$ string melting) of A Multi Phase Transport (AMPT) model-generated data. We adopted one of the most conventional and successful techniques, the Scaled Factorial Moment (SFM) method, using Monte Carlo (MC) data for 10 AGeV Au+Au collisions in search of intermittency in the model-generated data. Our interest is to search for intermittency behavior of particles that leads to multiplicity fluctuations and that would reveal a phase transition from hadronic matter to QGP. In this article, the intermittency values for both modes of AMPT data are presented. The results obtain some insight into the dynamics of heavy ion collisions and the formation of QGP.

Figures

Figures reproduced from arXiv: 2501.05175 by the authors.

Figure 1
Figure 1. mT spectra of π + at 8 AGeV energy using (a) AMPT generated data, and (b) ratio of mT spectra of experimental (E895 experiment[34]) data to AMPT generated data Plotted the mT -spectra of π + shown in [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. χ(px)vsχ(py) distribution for charged particles produced in 10 AGeV Au+Au collisions using AMPT event generator. px and py values of each primary charged particle are converted into new variables, χ(px) and χ(py). These transformations are defined as follows: Replace px by py in the following equation for χ(py). where px,min and py,min are the minimum values of px and py, respectively. px,max and py,max are the maxi… view at source ↗
Figure 3
Figure 3. Variations of hFqi against lnM2 for charged particles produced in Au+Au collisions at 10 AGeV energy using AMPT generator with default mode. The sizes of the bins follow the sequence of 1, 4, 9, 16, 25, 36, and so on up to 100, and the uniform width of each bin is dχ(px) × dχ(py). To compute SFM for a bin, we carefully count the number of particles in each bin and use these counts for that particular bin. The SFMs o… view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Variations of hFqi against lnM2 for charged particles produced in Au+Au collisions at 10 AGeV energy using AMPT generator with string melting mode. The hFqi for different values of q are then calculated for χ(px − py) space using both the sets of AMPT (default) and AMP…
Figure 5
Figure 5. Figure 5: αq vs. q plot for charged particle produced uing both the modes of AMPT generator. 2 3 4 5 6 q 0.005 0.01 0.015 q0.02 d AMPT (default) AMPT (SM) [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]
Figure 6
Figure 6. Figure 6: The variation of dq with q plot for all cahrged particles produced in Au+Au collisions at 10 AGeV energy using both modes of AMPT generator. 11 [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]
Figure 7
Figure 7. Figure 7: The hFq i vs. lnM2 plots for all charged particles produced in Au+Au collisions at 200 GeV energy using AMPT(default) generator [PITH_FULL_IMAGE:figures/full_fig_p012_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

37 extracted references · 36 canonical work pages

  1. [6]

    Gope and B

    S. Gope and B. Bhattacharjee, Eur. Phys. J. A , (2021)

    work page 2021

  2. [1]

    Busza et al., arXiv:1802.04801v2 [hep-ph] 23 Feb 2018

    W. Busza et al., arXiv:1802.04801v2 [hep-ph] 23 Feb 2018. 13

    arXiv 2018

  3. [2]

    QCD Phase Diagram: Phase Transition, Critical Point and Fluctuations

    B. Mohanty, arXiv: 0907.4476v3 [nucl-ex] 5 Oct 2009

    work page Pith review arXiv 2009

  4. [3]

    Odyniec, J

    G. Odyniec, J. Phys.: Conf. Ser. 455 012037 (2013)

    work page 2013

  5. [4]

    V. V. Braguta, Symmetry 15(7), 1466 (2023)

    work page 2023

  6. [5]

    Strodthoff et al., Phys

    N. Strodthoff et al., Phys. Rev. D 85, 074007 (2012)

    work page 2012

  7. [7]

    J. B. Kogut et al., Phys. Lett. B, 464, 1–2, 1999

    work page 1999

  8. [8]

    Anticic et al., (NA49 Collaboration), Eur

    T. Anticic et al., (NA49 Collaboration), Eur. Phys. J. C 75, 587 (20 15)

    work page

Show all 37 references
  1. [9]

    Bhattacharjee et al., Fractals, 26 1850015 (2018)

    S. Bhattacharjee et al., Fractals, 26 1850015 (2018)

    2018

  2. [10]

    Gope and B

    S. Gope and B. Bhattacharjee, Acta. Phys. Pol. B 53, 9-A3 (2 022)

  3. [11]

    Zhu et al., Physical Review C, 72(5), 051901 (2005)

    L. Zhu et al., Physical Review C, 72(5), 051901 (2005)

    2005

  4. [12]

    Zhang et al., Physical Review C, 61(1), 014901 (2000)

    B. Zhang et al., Physical Review C, 61(1), 014901 (2000)

    2000

  5. [13]

    Z. W. Lin et al., Physical Review C, 69(3), 034904 (2004)

    2004

  6. [14]

    Z. W. Lin et al., Physical Review C, 65(4), 044902 (2002)

    2002

  7. [15]

    Seddiki (CBM Collaboration), J

    S. Seddiki (CBM Collaboration), J. Phys. Conference Series 50 3 012027 (2014)

    2014

  8. [16]

    Satszel et al

    P. Satszel et al. (CBM Collaboration), Acta Physica Polonica B 41 341 (2010)

    2010

  9. [17]

    Bleicher et al., Journal of Physics G, 25(7), 1859-1874 (199 9)

    M. Bleicher et al., Journal of Physics G, 25(7), 1859-1874 (199 9)

  10. [18]

    Z. W. Lin et al.,Physical Review C, 72(6), 064901 (2005)

    2005

  11. [19]

    J. Wu et. al., Phys. Lett. B 801, 135186 (2020)

    2020

  12. [20]

    W., Nuclear Physics A, 774, 469-472 (2007)

    Xu, N., Lin, Z. W., Nuclear Physics A, 774, 469-472 (2007)

    2007

  13. [21]

    Zhang et al.,Physical Review C, 65(5), 054909 (2002)

    B. Zhang et al.,Physical Review C, 65(5), 054909 (2002)

    2002

  14. [22]

    Bialas and R

    A. Bialas and R. Peschanski, Nuclear Physics B, 273(2), 703-72 0 (1986). 14

    1986

  15. [23]

    Bialas and R

    A. Bialas and R. Peschanski, Physical Review D, 35(10), 3327-3 334 (1987)

    1987

  16. [24]

    Bialas et al., Nucl

    A. Bialas et al., Nucl. Phys. B 308 857-867 (1988)

    1988

  17. [25]

    R. C. Hwa and Q. Zhang, Phys. Rev. D 62 014003 (2000)

    2000

  18. [26]

    Davidson and L

    N. Davidson and L. Stumpf, Journal of Physics G, 22(4), 661- 674 (1996)

    1996

  19. [27]

    Frye and R

    M. Frye and R. Smith, European Physical Journal C, 35(3), 40 3-415 (2004)

    2004

  20. [28]

    V. Koch, V. and R. C. Hwa, Physical Review D, 44(5), 1422-143 5 (1991)

    1991

  21. [29]

    Bialas, R

    A. Bialas, R. Peschanski, Physical Review C, 60(4), 045206 (19 99)

  22. [30]

    Ferreira et al., Nuclear Physics A, 641(1), 120-138 (1998)

    A. Ferreira et al., Nuclear Physics A, 641(1), 120-138 (1998)

    1998

  23. [31]

    Aichelin et al., Physics Reports, 370(5), 245-288 (2002)

    J. Aichelin et al., Physics Reports, 370(5), 245-288 (2002)

    2002

  24. [32]

    Stephanov et al., Reviews of Modern Physics, 81(2), 387-40 7 (2008)

    M. Stephanov et al., Reviews of Modern Physics, 81(2), 387-40 7 (2008)

    2008

  25. [33]

    T. S. Biro et al., Journal of Physics G: Nuclear and Particle Physic s, 30(5), 927-934 (2004)

    2004

  26. [34]

    J. L. Klay et al., E895 Collaboration, Phys. Rev. C 68 054905 (200 3)

  27. [35]

    Ajinenko et al.), Phys

    NA22 Collab., (I.V. Ajinenko et al.), Phys. Lett. B 222 306 (1989)

    1989

  28. [36]

    C. B. Chiu et al, Mod. Phys. Lett. A 5, 2651 (1990)

    1990

  29. [37]

    N. M. Agarbabyan et al., Phys. Lett. B 382, 305 (1996). 15

    1996

Pith tools

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