REVIEW 4 major objections 6 minor 1 cited by
A new open-source clusterization library reproduces established heavy-ion cluster finders and adds a stable-cluster tracking step that recovers nuclei lost by transport models.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-03 19:10 UTC pith:M2UABM45
load-bearing objection A genuinely useful open-source clusterization library with real algorithmic improvements (DSU MST, two-pass SA, sMST tracking), but the sMST validation is only circumstantial and the binding-energy criterion is never defined. the 4 major comments →
Modifications in clusterization procedures for heavy-ion collisions: minimum spanning tree, simulated annealing, coalescence
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that all major clusterization approaches used in heavy-ion transport simulations—minimum spanning tree, simulated annealing, coalescence, and a hybrid of the last two—can be implemented in one open library and validated on exactly the same events. The DSU-based MST gives identical results to the full-graph MST at lower cost; the modified two-pass simulated annealing chain reaches SACA-like agreement with data while reducing computational time; and the new sMST tracking algorithm repairs the artificial dissolution of previously stable clusters during QMD evolution, substantially improving mid-rapidity deuteron and helium-3 yields. For helium-3, the sMST result matches the
What carries the argument
The load-bearing object is the binding-energy-per-nucleon filter, E_bind(C) < E_cut with E_cut = -4 MeV, computed from the four-momenta of constituents in the cluster rest frame, which is reused inside the SA chain, the sMST tracking step, and the mixed coalescence procedure. Alongside it, the sMST tracking algorithm maintains a persistent list of stable clusters, removing a cluster only if a constituent disappears or collides after registration, and applying a preferential-growth / controlled-dissociation update rule to decide whether a new candidate replaces overlapping stored clusters.
Load-bearing premise
The whole filtering and tracking chain assumes that the cluster rest-frame binding energy per nucleon, with the chosen threshold E_cut = -4 MeV, genuinely separates physically bound nuclei from accidental groupings; the paper does not write the formula or validate it against known nuclear binding energies.
What would settle it
Compute the binding-energy-per-nucleon spectrum for identified deuterons, tritons, helium-3, and helium-4 using the same four-momentum formula and compare the peak positions to the empirical binding energies. If the formula labels obviously bound nuclei as unbound (or vice versa) by more than a few MeV, the E_cut filter and every conclusion built on it lose their physical basis.
If this is right
- Cluster finding for heavy-ion simulations can be done without building the full graph, giving identical MST results at lower computational cost.
- A two-stage simulated annealing chain that treats unstable MST clusters one at a time reaches SACA-like agreement with data while reducing runtime.
- Enabling the sMST tracking step restores mid-rapidity deuteron and helium-3 yields that otherwise drop, with the helium-3 yield matching measured data better than spin-isospin coalescence on the same baryons.
- At early times, cluster-finding parameters should be chosen according to the system evolution rather than fixed for the whole reaction.
- Applying the binding-energy cut to coalescence candidates suppresses yields too strongly, but a mixed coalescence-plus-binding-energy procedure offers a way to tune parameters for MST-based approaches.
Where Pith is reading between the lines
- The paper's entire filtering chain rests on an E_bind formula that is never written out; validating that formula against empirical nuclear binding energies would be a direct check of whether the 'stable' clusters are truly bound.
- The same-baryon-set comparison between MST and coalescence for helium-3 suggests a clean experimental discriminator between formation mechanisms, provided both procedures are calibrated to identical acceptance.
- The sMST tracking rule could be adapted to other transport models that do not propagate correlations, and its sensitivity to the 'collision after registration' criterion could be studied by varying that rule.
- The momentum-proximity checks used in the mixed procedure hint that the default MST chain could be extended with relative-momentum cuts, which the paper suggests but does not build into the default algorithm.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a new open-source clusterization library (CCL) for heavy-ion transport simulations, implementing MST-based cluster finding, a two-pass simulated annealing chain, a stable-cluster tracking algorithm ('sMST'), a 'box-like' coalescence procedure, and a 'mixed' coalescence+binding-energy method. The library is applied to PHQMD events and compared with established implementations (PHQMD MST, SACA, UrQMD coalescence) and with ALADiN and NA49 data. The central claims are: (i) CCL reproduces established MST/SACA results on the same events; (ii) sMST tracking 'drastically improves' the stability and yield of mid-rapidity clusters; and (iii) for 3He, sMST matches NA49 data better than coalescence applied to the same baryon set. The paper also reports a first comparison of MST-based versus coalescence 3He yields on identical events.
Significance. If the claims hold, CCL would be a valuable open resource: the field suffers from closed-source, model-specific clusterization tools, and a modular, well-tested library with a direct comparison on identical events is genuinely useful. The authors also provide several concrete algorithmic contributions (DSU-based MST, a two-pass SA chain, an online tracking method that does not require freeze-out time). The paper is honest about some limitations, such as the need for time-dependent SA parameters at low energies. However, the central validation currently rests on two fragile pillars: the binding-energy formula is never defined, and the sMST tracking is validated only by comparison to data on a single event sample without any ground-truth test. The 'better than coalescence' conclusion is likewise based on a single coalescence parameter set. These are correctness-risk issues that can be addressed with additional definitions and controlled tests, but they currently prevent acceptance.
major comments (4)
- [Sec. 3.1.1] The binding energy per nucleon E_bind(C) is the central physical quantity of the paper: it drives the SA acceptance criterion, the E_cut=-4 MeV stability cut in Sec. 3.2, the sMST stability filter in Sec. 4, and the 'mixed' procedure in Sec. 5. Yet it is never written down. The reader cannot reproduce any result, nor can the physical content of 'bound' clusters be assessed. Please provide the explicit formula, including how the cluster rest frame is defined and how the four-momenta are combined (e.g., sum of four-momenta, invariant mass, subtraction of constituent masses). This is a load-bearing omission.
- [Sec. 4, Fig. 4] The claim that sMST 'drastically improves mid-rapidity cluster stability' and 'recovers lost clusters' is not validated against any ground truth. The only evidence is the comparison of d and 3He rapidity spectra with tracking on/off on one PHQMD event sample. The update rule has an inherent asymmetry: a candidate p replaces overlapping stored clusters iff A(p)>A(s_max), otherwise it is discarded. In a dense mid-rapidity environment this rule can merge distinct objects and inflate cluster yields; the only physical filter is the negative-binding-energy test, which is undefined (see above). No false-positive rate, no synthetic-event test, and no comparison with an independent tracking method (e.g., aMST or collision-history-based reconstruction) are provided. Please add at least one controlled test with known cluster content and quantify over-merging and under-recovery.
- [Sec. 5.2, Fig. 7] The conclusion that sMST describes NA49 3He data 'better than coalescence' is based on a single coalescence parameter set (M1, as stated for Fig. 6). Table 3 lists other sets, and M4 is explicitly described as chosen to better describe the data; indeed, in the mixed procedure of Fig. 9, the M4 set with reduced Δr=2.5 fm reproduces the 3He shape much better than M1. Thus the claim is not robust to reasonable parameter choices. Please show the coalescence results for all parameter sets in the same comparison, or clearly state that the conclusion refers only to the default M1 set. Also include statistical uncertainties on the model curves, especially given the limited 50k-event statistics noted in the text.
- [Sec. 6, first bullet] The summary states that the MST-based clusterization overestimates ALADiN data even at t=200 fm/c, while the SA procedures with time-dependent parameters (S1 at t=50 fm/c, S2 at t=100 fm/c) achieve 'perfect' agreement. This means the ALADiN comparison is not a predictive, parameter-free test: the control parameters are adjusted per time step to match the data. The paper is partially transparent about this in Sec. 3.3, but the abstract's wording ('tested with complex observables') does not convey that the agreement involves parameter tuning. Please state this limitation explicitly in the abstract or in the conclusions, and specify which data points were used for tuning and which are genuinely predictions.
minor comments (6)
- [Abstract] There are several formatting/typo issues: '⟨Mimf⟩' and '⟨Zmax⟩' are rendered with missing spaces; 'same set of events' is fine but the abstract does not mention that the sMST tracking and the coalescence comparison use 50k events only.
- [Sec. 3.2] The threshold E_cut=-4 MeV is introduced without justification or a sensitivity study. Please provide a reference or a brief rationale, and consider adding a scan over E_cut to show that the conclusions are robust.
- [Sec. 4] The notation E^cut_bind (used in the text) is inconsistent with E_cut used in Sec. 3.2. Please unify.
- [Fig. 2] The bottom panel label reads '⟨∑ E_bind⟩' but the text refers to '⟨E_bind⟩'; the figure should match the text.
- [Sec. 3.1.3] The steps mode formula in Sec. 3.1.2 uses 'ℓ' for temperature level index, but later in Sec. 3.1.3 the stagnation counter and step definitions would benefit from a concrete pseudo-code or a reference to the code repository.
- [Sec. 5.2] The statement 'the final result must be the sum of kinetic deuterons and those found by (s)MST or formed via coalescence' is important, but Fig. 6 does not show the sum of coalescence and kinetic deuterons (only coalescence alone, and sMST+kinetic). A plot of the coalescence+kinetic sum would make the comparison more direct.
Circularity Check
No significant circularity: the paper's central comparisons are external benchmarks against PHQMD MST, SACA, UrQMD coalescence, and ALADIN/NA49 data.
full rationale
The validity of the main claims rests on comparisons with existing closed-source implementations (PHQMD MST, SACA, UrQMD built-in coalescence) on the same events and with experimental data (ALADIN, NA49). The CCL MST results are checked directly against the PHQMD full-graph MST, the CCL SA routine against SACA, and the CCL coalescence against the UrQMD built-in implementation; these are independent external benchmarks, not restatements of CCL parameters. The explicit statement that the M4 set was 'chosen to better description of available data' (Sec. 5.2) is transparent fitting, not a disguised prediction, and does not support the paper's central sMST-vs-coalescence comparison, which is anchored to NA49 data. Self-citations to refs. [9], [10], and [12] are contextual (code provenance and background) and are not load-bearing for the derivation. The unstated binding-energy formula and user-chosen E_cut threshold in Sec. 3.2/4 are a physical-validation gap, but not circularity, since the sMST stability filter's output is ultimately tested against external NA49 yields. No equation was found that is equal to its input by construction, and no fitted parameter is renamed as a prediction.
Axiom & Free-Parameter Ledger
free parameters (5)
- Clusterization radii r_clust =
2.5, 3.5, 4.0, 9.5 fm
- Momentum thresholds delta_p =
0.285-0.55 GeV/c
- SA control parameters (Tmax, Tmin, steps, Pnew, smin, sden) =
Tmax=1, Tmin=1e-15, steps 500/4000, Pnew 0.1-0.9, smin=50, sden=2
- Binding energy cut E_cut =
-4 MeV
- Coalescence parameter set M4 =
delta_r/delta_p combinations
axioms (4)
- domain assumption The binding energy per nucleon computed from four-momenta in the cluster rest frame correctly identifies physically bound nuclei.
- domain assumption QMD-based transport (PHQMD) propagates the phase-space correlations needed to form final-state clusters.
- domain assumption Clusters found by coordinate-space proximity (MST) at late times correspond to final-state fragments.
- domain assumption The 'lost cluster' phenomenon in QMD is an artifact that can be corrected by tracking stable clusters.
read the original abstract
A new open-source cluster finding library ''Common Clusterization Library'' (CCL) is proposed to describe the clusters production when applied to the transport codes. The new library was applied to the Parton-Hadron-Quantum-Molecular Dynamics approach to be comparable with the established and well known MST and SACA implementations on the exactly same set of the events. The presented procedures were tested with rapidity density distributions and complex observables like $\langle M_{imf} \rangle$ and $\langle Z_{\text{max}} \rangle$ vs $Z_{bound}2$ at the energies of ALADiN and NA49 experiments. An improvement to the simulated annealing chain is proposed to reduce the computational time. The new stable clusters tracking algorithm for the ''lost'' cluster recovering is presented. The results of the ''box-like'' coalescence procedure and a mixed method which incorporates the coalescence algorithm steps and the negative binding energy criteria are presented. The difference between the Helium-3 found by the MST-based cluster finding procedure and the coalescence algorithm applied on the exactly same set of events and baryons is shown for the first time.
Figures
Forward citations
Cited by 1 Pith paper
-
Freeze-out model of light nuclei formation in heavy-ion collision transport
A hybrid freeze-out model combining transport dynamics and thermal cluster production is proposed for light nuclei in intermediate-energy heavy-ion collisions.
Reference graph
Works this paper leans on
-
[1]
J. Aichelin, ’Quantum’ molecular dynamics: A Dynam- ical microscopic n body approach to investigate frag- ment formation and the nuclear equation of state in heavy ion collisions, Phys. Rept. 202 (1991) 233–360. doi:10.1016/0370-1573(91)90094-3
-
[2]
Bor ˇuvka, Über ein Minimalproblem., Práce moravské pˇrirodovˇedecké spoleˇcnosti 3 (1926), 37–58
O. Bor ˇuvka, Über ein Minimalproblem., Práce moravské pˇrirodovˇedecké spoleˇcnosti 3 (1926), 37–58. (1926)
1926
-
[3]
C. Dorso, J. Randrup, Early recognition of clusters in molecular dynamics, Phys. Lett. B 301 (1993) 328–333. doi:10.1016/0370-2693(93)91158-J
-
[4]
Pincus, Letter to the editor—a monte carlo method for the approximate solution of cer- tain types of constrained optimization prob- lems, Operations Research 18 (6) (1970) 1225–
M. Pincus, Letter to the editor—a monte carlo method for the approximate solution of cer- tain types of constrained optimization prob- lems, Operations Research 18 (6) (1970) 1225–
1970
-
[5]
N. Metropolis, A. W. Rosenbluth, M. N. Rosen- bluth, A. H. Teller, E. Teller, Equation of state calculations by fast computing machines, The Jour- nal of Chemical Physics 21 (6) (1953) 1087–1092. doi:10.1063/1.1699114
-
[6]
R. K. Puri, J. Aichelin, Simulated annealing clusterization algorithm for studying the multifragmentation, J. Com- put. Phys. 162 (2000) 245–266. arXiv:nucl-th/9811018, doi:10.1006/jcph.2000.6534
Pith/arXiv arXiv 2000
-
[7]
A. Le Fèvre, J. Aichelin, C. Hartnack, Y . Leifels, FRIGA: A new approach to identify isotopes and hypernuclei inn-body transport models, Phys. Rev. C 100 (3) (2019) 034904. arXiv:1906.06162, doi:10.1103/PhysRevC.100.034904
Pith/arXiv arXiv 2019
-
[8]
Kireyeu, https://github.com/vkireyeu/ccl (2025)
V . Kireyeu, https://github.com/vkireyeu/ccl (2025)
2025
-
[9]
Kireyeu, Cluster dynamics studied with the phase- space Minimum Spanning Tree approach, Phys
V . Kireyeu, Cluster dynamics studied with the phase- space Minimum Spanning Tree approach, Phys. Rev. C 103 (5) (2021) 054905. arXiv:2103.10542, doi:10.1103/PhysRevC.103.054905
Pith/arXiv arXiv 2021
-
[10]
V . Kireyeu, G. Coci, S. Glaessel, J. Aichelin, C. Blume, E. Bratkovskaya, Cluster formation near midrapidity: How the production mechanisms can be identified ex- perimentally, Phys. Rev. C 109 (4) (2024) 044906. arXiv:2304.12019, doi:10.1103/PhysRevC.109.044906
Pith/arXiv arXiv 2024
-
[11]
B. A. Galler, M. J. Fisher, An improved equiva- lence algorithm, Commun. ACM 7 (5) (1964) 301–303. doi:10.1145/364099.364331
arXiv 1964
-
[12]
G. Coci, S. Gläßel, V . Kireyeu, J. Aichelin, C. Blume, E. Bratkovskaya, V . Kolesnikov, V . V oronyuk, Dy- namical mechanisms for deuteron production at mid- rapidity in relativistic heavy-ion collisions from ener- gies available at the GSI Schwerionensynchrotron to those at the BNL Relativistic Heavy Ion Collider, Phys. Rev. C 108 (1) (2023) 014902. arX...
Pith/arXiv arXiv 2023
-
[13]
T. Anticic, et al., Production of deuterium, tritium, and He3 in central Pb+Pb collisions at 20A,30A,40A,80A , and 158A GeV at the CERN Super Proton Synchrotron, Phys. Rev. C 94 (4) (2016) 044906. arXiv:1606.04234, doi:10.1103/PhysRevC.94.044906
Pith/arXiv arXiv 2016
-
[14]
S. T. Butler, C. A. Pearson, Deuterons from High-Energy Proton Bombardment of Matter, Phys. Rev. Lett. 7 (1961) 69–71. doi:10.1103/PhysRevLett.7.69
-
[15]
H. H. Gutbrod, A. Sandoval, P. J. Johansen, A. M. Poskanzer, J. Gosset, W. G. Meyer, G. D. Westfall, R. Stock, Final State Interactions in the Production of Hydrogen and Helium Isotopes by Relativistic Heavy Ions on Uranium, Phys. Rev. Lett. 37 (1976) 667–670. doi:10.1103/PhysRevLett.37.667
-
[16]
S. A. Bass, et al., Microscopic models for ultrarelativis- tic heavy ion collisions, Prog. Part. Nucl. Phys. 41 (1998) 255–369. arXiv:nucl-th/9803035, doi:10.1016/S0146- 6410(98)00058-1
Pith/arXiv arXiv 1998
-
[17]
M. Bleicher, et al., Relativistic hadron hadron collisions in the ultrarelativistic quantum molecular dynamics model, J. Phys. G 25 (1999) 1859–1896. arXiv:hep-ph/9909407, doi:10.1088/0954-3899/25/9/308
Pith/arXiv arXiv 1999
-
[18]
S. Sombun, K. Tomuang, A. Limphirat, P. Hill- mann, C. Herold, J. Steinheimer, Y . Yan, M. Bleicher, Deuteron production from phase-space coalescence in the UrQMD approach, Phys. Rev. C 99 (1) (2019) 014901. arXiv:1805.11509, doi:10.1103/PhysRevC.99.014901
Pith/arXiv arXiv 2019
-
[19]
T. Reichert, J. Steinheimer, V . V ovchenko, B. Döni- gus, M. Bleicher, Energy dependence of light hypernu- clei production in heavy-ion collisions from a coales- cence and statistical-thermal model perspective, Phys. Rev. C 107 (1) (2023) 014912. arXiv:2210.11876, doi:10.1103/PhysRevC.107.014912
Pith/arXiv arXiv 2023
-
[20]
T. Reichert, M. Omana Kuttan, A. Kittiratpattana, N. Buyukcizmeci, A. Botvina, J. Steinheimer, M. Ble- icher, Directed and elliptic flow of light nuclei and hy- pernuclei in Au+Au collisions at √sNN =3 GeV: Coalescence vs. Statistical Fragmentation (4 2025). arXiv:2504.17389. 11
Pith/arXiv arXiv 2025
-
[1228]
arXiv:https://doi.org/10.1287/opre.18.6.1225, doi:10.1287/opre.18.6.1225
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.