{"id":"92b6b4d2-5046-437a-9187-a9119b7d6f37","arxiv_id":"1908.04956","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A Wigner-function-based quark coalescence algorithm in the AMPT model improves the energy dependence of proton directed flow and reproduces baryon-baryon anti-correlations.","lead":"This paper changes how quarks combine into hadrons in the AMPT model, using Wigner functions so that partons close in both position and momentum are favored. The new algorithm preserves partonic flow information better and produces a qualitatively correct proton directed flow sign at RHIC beam energies, plus baryon anti-correlations seen in data.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central v1 validation compares model protons at 20–30% centrality with STAR data from the 10–40% bin; because directed flow is centrality dependent, the claimed sign sequence may be a comparison artifact.","rationale":"The reader's weakest_assumption focuses on the vacuum-fixed Wigner widths and statistical factors gM and gB, and that is a legitimate concern: Eq. (7) makes the meson-versus-baryon choice by comparing absolute Wigner-function values, so these coefficients are not physically innocuous, and no sensitivity study is provided. I regard this as a secondary issue rather than the most load-bearing one, because the paper's headline claim is specifically the qualitative reproduction of the proton directed-flow slope, and that claim rests on Fig. 8. There, the model centrality (20–30%) differs from the STAR bin (10–40%) shown in the same figure, and the parton cross section is simultaneously recalibrated to v2. A sign sequence that is robust only to one of these choices would not validate the hadronization mechanism. The proposed centrality-matched rerun is a single, decisive check that does not require new theory input. Since the reader already assigned CONDITIONAL and flagged the centrality issue in the rationale, my read does not change the verdict; it sharpens the condition: match the experimental centrality and test sensitivity of the v1 result to the Wigner parameters before accepting the central claim.","tokens_in":19360,"tokens_out":8058,"duration_ms":95825,"concrete_test":"Recompute panel (b) of Fig. 8 using the improved AMPT model in the same centrality bin as the STAR data (10–40%), keeping the same parton cross section (3 mb), Lund parameters, and event-plane analysis. Then extract dv1/dy at y=0 for final protons at 7.7, 11.5, 19.6, 27, and 39 GeV, and compare the sign sequence with STAR. If the positive-to-negative transition moves or disappears in the matched bin, the central claim loses its main empirical support; if the sign sequence is unchanged, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's strongest claim is anchored by Fig. 8 in Sec. III.D: the improved AMPT model is said to preserve parton directed flow through hadronization and to give a qualitatively correct collision-energy dependence of the proton v1 slope. But the model curves are computed for midcentral 20–30% Au+Au collisions, while the STAR points shown for comparison are for midcentral 10–40% collisions, as stated in the figure caption. Directed flow is known to be sensitive to centrality: the magnitude and even the sign of dv1/dy at a fixed energy can change with the centrality selection. The comparison is further confounded because the improved model uses a parton scattering cross section (3 mb) that is refitted to reproduce v2, whereas the original model uses 1.5 mb. Thus the headline test has two uncontrolled degrees of freedom: a mismatched centrality bin and a recalibrated cross section. If the positive-to-negative sign transition in Fig. 8 appears only because of the narrower 20–30% selection, then the claim that the Wigner-function coalescence 'preserves better the parton dynamics' is not established. The Wigner-parameter sensitivity identified by the reader is a real secondary issue, but even a complete sensitivity scan of gM/gB and the widths would not resolve the centrality mismatch in the decisive comparison.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper modifies the hadronization step of the string-melting AMPT model by replacing spatial-only nearest-neighbor coalescence with a greedy algorithm that forms mesons, baryons, and antibaryons from parton combinations having the largest Wigner function in phase space. The Wigner widths are fixed using vacuum pion and proton charge radii, the statistical weights are set to gM=1/36 and gB=1/108, and net baryon, electric, and strangeness charges are conserved by consuming all partons. The authors compare the original and improved versions for elliptic flow and vn scaling at RHIC and LHC energies, the directed-flow slope at RHIC beam energy scan energies, particle yield ratios, and di-hadron correlations in p+p collisions. The main claims are that the improved algorithm preserves parton dynamics better through hadronization, reproduces the qualitative collision-energy dependence of the proton v1 slope, gives reasonable yield ratios at top RHIC and LHC energies, and produces the near-side anti-correlation in baryon-baryon and antibaryon-antibaryon correlations seen by ALICE.","tokens_in":19610,"tokens_out":5482,"duration_ms":54477,"significance":"Hadronization is a major uncertainty in interpreting transport model results, and the proposed algorithm is a concrete, reproducible modification of a widely used code. The paper's strengths are that the Wigner-function inputs are tied to measured charge radii, the algorithm removes the coalescence-ordering ambiguity of the original model, and several predictions (the v1 sign sequence, di-hadron correlations, and yield ratios) are not used to fit parameters. If the centrality and parameter-sensitivity concerns below are resolved, the paper would provide a useful benchmark for future AMPT studies. However, the headline v1 validation is weakened by a centrality mismatch and by a simultaneous rescaling of the parton cross section, so the significance of the central claim is currently conditional.","major_comments":[{"comment":"The central validation of the v1 claim compares model protons in midcentral 20-30% Au+Au collisions with STAR data from the 10-40% centrality bin, as stated in the figure caption, while the improved model also uses a parton scattering cross section of 3 mb rather than the 1.5 mb used by the original model. Since the directed-flow slope is known to depend on centrality, the apparent reproduction of the positive-to-negative sign change may be an artifact of the narrower centrality selection or of the recalibrated cross section. The authors should either compute the same centrality bin as the data or demonstrate that the sign sequence is stable across centralities and for a fixed parton cross section.","section":"Section III.D, Fig. 8"},{"comment":"The ranking of meson versus baryon formation in Eq. (7) depends on the relative magnitudes of fM, fB, and fbarB, which are controlled by the Gaussian widths fixed to vacuum pion and proton charge radii and by the statistical weights gM=1/36 and gB=1/108. These choices are asserted without a sensitivity test, and if in-medium widths or relative weights differ, the preferred parton combinations and all downstream hadron observables change. The paper should include a sensitivity scan over these widths and weights, or provide a physical derivation for them, before claiming that the Wigner-function approach is parameter-free.","section":"Section II, Eqs. (1)-(7)"},{"comment":"The algorithm selects the globally largest Wigner function and consumes all partons, rather than sampling from the Wigner function as a formation probability, so the relation between the stated 'formation probabilities' and the actual greedy maximum selection is not established. The final paragraph acknowledges that partons largely separated in phase space are still forced to coalesce; this limitation should be quantified by comparing the greedy selection with a probabilistic sampling version or with a version that leaves non-coalescing partons to fragmentation.","section":"Section II, Eq. (7); Section IV"}],"minor_comments":[{"comment":"There are several typographical errors: the title has 'mult iphase' instead of 'multiphase', Section III.A has 'scarifices' instead of 'sacrifices', and 'presumedly' should be 'presumably'.","section":"Title and Section III.A"},{"comment":"In the paragraph describing the Lund string fragmentation parameters, 'from 7.7 GeV to 39 TeV' should read 'from 7.7 GeV to 39 GeV', since the paper otherwise discusses RHIC BES energies.","section":"Section III"},{"comment":"The text states that the improved model's near-side anti-correlation is qualitatively consistent with ALICE data, but Fig. 10 contains only model curves and no data points; the comparison would be more transparent if the ALICE data were included in the figure or if the specific ALICE reference were cited directly in the caption.","section":"Section III.F, Fig. 10"}],"recommendation":"major_revision","confidential_remarks":"The centrality mismatch in Fig. 8 is the main obstacle to accepting the central v1 claim. I would want either a recomputation at the STAR centrality bin or a demonstration that the sign sequence is centrality-independent. The relation to Ref. [49] should also be clarified: the authors state that similar results are obtained there, but the novelty of the present work relative to Ref. [49] is not made explicit."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid, useful model-improvement paper, and the main new idea is worth taking seriously, but the headline directed-flow comparison has a binning mismatch that should be fixed before I would trust the sign-change claim.\n\nWhat's new: the authors replace the AMPT string-melting spatial coalescence with a Wigner-function-weighted phase-space coalescence. The algorithm is clearly specified: pick the maximal Wigner function among meson, baryon, and antibaryon candidates, conserve net charges, and remove the old ordering ambiguity. They then show the improved version preserves parton momentum-space information better and connect it to three non-fit observables: the proton v1 sign sequence at RHIC BES, near-side baryon-baryon and antibaryon-antibaryon anti-correlations in p+p at 7 TeV, and several antiparticle-to-particle yield ratios at top RHIC/LHC energies. Those are real checks, not calibration, and the anti-correlation result in particular is a nice qualitative hit that the original spatial coalescence misses.\n\nCredit where due: the mechanism is described with enough equations to reproduce. The v2 agreement is honestly presented as calibration, since the parton cross section is refitted for each model version. The paper also explicitly admits the BES yield ratios are not reproduced and that additional physics is needed. That is good scholarly behavior.\n\nSoft spots, in order of severity. First, the central v1 claim in Fig. 8 compares model protons at 20-30% centrality with STAR data from the 10-40% bin. Directed flow is centrality dependent, and the sign and magnitude of dv1/dy can shift with centrality class. With the cross section also changed from 1.5 to 3 mb, the comparison has two uncontrolled degrees of freedom. The qualitative sign sequence may survive, but the paper does not currently show that. This is the main comparison, so it needs a same-centrality figure or a centrality robustness test. Second, the Wigner widths are fixed by vacuum pion and proton charge radii, and the statistical factors gM and gB are set by hand; there is no sensitivity scan. The ranking in Eq. (7) depends on those choices. The reader's concern is real, though secondary. Third, the v2 agreement is calibration, not validation, so it should not be counted as independent support.\n\nMinor point: the greedy maximum-selection algorithm is a heuristic approximation to dynamical coalescence. That is acceptable for a transport model, but it should be acknowledged as such.\n\nVerdict: conditional. This paper deserves a serious referee, and with a same-centrality v1 plot and a Wigner-parameter sensitivity study it could be a publishable contribution. I would bring it to group and would cite it if I were working on coalescence or AMPT. Recommend peer review, with those two revisions requested.","headline":"A useful, clearly specified Wigner-coalescence upgrade for AMPT, but the headline v1 comparison uses mismatched centralities and needs a sensitivity check before the sign-change claim is settled.","tokens_in":20158,"tokens_out":2700,"would_cite":true,"duration_ms":27853,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The improved hadronization algorithm preserves partonic directed flow, reproducing the measured proton flow sign change with collision energy.","keywords":["Wigner function","quark coalescence","hadronization","AMPT model","directed flow","di-hadron correlations","anisotropic flow","beam energy scan"],"falsifier":"Vary the Gaussian widths in Eqs. (3) and (6) by $\\pm 30\\%$ or change the $g_M/g_B$ ratio while rerunning the hadronization on the same freeze-out parton distribution; if the proton $\\mathrm{d}v_1/\\mathrm{d}y$ sign change between 11.5 and 19.6 GeV or the near-side baryon-antibaryon anticorrelation disappears, the central claim is not robust. Also compare the resulting $v_1$ slope magnitude quantitatively with the measured values, since the paper only claims qualitative agreement.","tokens_in":19098,"feed_emoji":"⚛️","tokens_out":12558,"duration_ms":101820,"temperature":0.7,"pith_summary":"The paper rewrites the hadronization step of the string-melting multiphase transport model so that partons are combined into hadrons not by spatial proximity alone but by the highest value of a Wigner function, a Gaussian in both relative coordinate and relative momentum. The widths are fixed by the measured charge radii of pions and protons, and the relative meson/baryon weights are fixed constants, so the algorithm introduces no free parameters. The payoff is that hadronization no longer scrambles the partonic dynamics: the proton directed-flow slope now changes sign with collision energy in the pattern seen in beam-energy-scan data, and baryon-baryon and antibaryon-antibaryon correlations in p+p collisions acquire the near-side anticorrelation measured in high-energy hadron collisions. These features were absent in the original spatial-coalescence algorithm. If the claim holds, hadronization is a controllable element of transport simulations rather than an adjustable black box.","feed_headline":"New hadronization rule reproduces proton flow sign change","feed_subtitle":"Picking partons close in momentum as well as position reproduces the proton flow sign change seen in heavy-ion data.","key_machinery":"The Wigner function of the valence parton combination. For a quark-antiquark pair it reads $f_M(\\rho,k_\\rho) = 8g_M\\exp(-\\rho^2/\\sigma_\\rho^2 - k_\\rho^2\\sigma_\\rho^2)$, with $\\rho$ and $k_\\rho$ the relative coordinate and momentum in the pair's center-of-mass frame; for baryons an analogous product over two Jacobi coordinates. The selection rule is to form the hadron from the combination whose Wigner function is largest, repeated until all partons are used, which simultaneously removes the ordering ambiguity between meson-first and baryon-first coalescence and favors partons that are close in momentum as well as space. The widths $\\sigma_\\rho,\\sigma_\\lambda$ are fixed from the root-mean-square charge radii of pions (0.61 fm) and protons (0.877 fm) through Eqs. (3) and (6).","core_discovery":"The central discovery is that the choice of which partons coalesce controls whether the hadron phase inherits the partonic directed flow. With the improved algorithm, the slope $\\mathrm{d}v_1/\\mathrm{d}y|_{y=0}$ for protons is negative at collision energies above 7.7 GeV and flips sign between 11.5 and 19.6 GeV in agreement with measured beam-energy-scan results, whereas the original nearest-in-space coalescence turns these slopes positive at all energies. The same algorithm produces a dip instead of a peak on the near side of baryon-baryon and antibaryon-antibaryon azimuthal correlations in p+p collisions at 7 TeV, consistent with measured data. This is achieved by ranking all possible two- and three-parton combinations by the value of the corresponding Wigner function, with Gaussian widths fixed by vacuum pion and proton radii.","pith_inferences":["Because the widths are fixed by vacuum charge radii, a natural test is to rerun the algorithm with in-medium or temperature-dependent widths; if the $v_1$ sign change or the p+p anticorrelation is sensitive to those widths, then hadronization itself is encoding in-medium information that current transport models discard.","The same Wigner-function ranking could be applied to leftover partons that are currently forced into hadrons despite being far apart in phase space; fragmenting those instead would provide a direct handle on the transition between coalescence and fragmentation at high transverse momentum.","If the improved hadronization preserves partonic dynamics this well, other observables that encode early-time flow, such as the mass ordering of $v_2$ at low $p_T$ or HBT radii, may shift as well; the paper does not examine these.","The larger cross sections required to fit the elliptic flow after the hadronization change suggest that earlier extractions of the specific shear viscosity from this transport model would need to be revisited with the new algorithm."],"forward_implications":["The hadronization step acts as a dynamical filter: from the same parton freeze-out distribution, spatial coalescence overestimates hadron elliptic flow, so the parton-scattering cross section needed to match data rises from 1.5 mb to 3 mb at 200 GeV and to 2 mb at 2.76 TeV.","The near-side anticorrelation between two baryons or two antibaryons in p+p collisions emerges naturally from momentum-favored coalescence rather than from a separate fragmentation mechanism.","The sign change of the proton directed-flow slope with collision energy survives hadronization only if coalescence preserves the partonic flow; the original spatial coalescence destroys the sign change entirely.","The relative production of strange baryons and antibaryons depends strongly on whether mesons or baryons are formed first; the Wigner ranking removes this ordering ambiguity and brings these ratios closer to measured values at top collision energies.","Both the original and improved algorithms preserve the constituent-quark-number scaling of anisotropic flows, but the improved model yields better $v_n/v_2^{n/2}$ scaling, especially at the higher collision energy."],"supporting_citations":[{"why":"Defines the string-melting multiphase transport model whose hadronization is being improved.","marker":"[3]"},{"why":"Introduces the dynamical-coalescence Wigner-function method adopted here.","marker":"[12]"},{"why":"Reports the measured proton directed-flow slopes at beam-energy-scan energies that the improved model reproduces.","marker":"[38]"},{"why":"Shows hadronization significantly affects proton directed flow, motivating the improved algorithm.","marker":"[43]"},{"why":"Provides the p+p di-hadron correlations in meson, baryon, and antibaryon channels that the improved model matches.","marker":"[48]"},{"why":"Introduced a meson/baryon preference parameter in AMPT coalescence, the direct predecessor of the Wigner-function ranking.","marker":"[53]"},{"why":"Supplies the pion and proton charge radii that fix the Gaussian widths of the Wigner functions.","marker":"[55]"},{"why":"Independently shows an improved hadronization can reproduce the baryon-antibaryon anticorrelation.","marker":"[49]"}],"fun_headline_variants":["Momentum-aware coalescence flips proton flow sign","Wigner function hadronization matches proton flow reversal","Proton flow sign change from momentum-space pairing","New coalescence rule reproduces proton flow sign jump"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The ranking of meson versus baryon formation depends on fixed Gaussian widths taken from vacuum pion and proton radii and on hand-set statistical weights of $1/36$ and $1/108$; if these do not represent the in-medium coalescence probability for freeze-out partons, the preferred parton combinations and all downstream observables change, yet the paper does not test this sensitivity.","fun_headline_variants_meta":{"raw":{"variants":["Momentum-aware coalescence flips proton flow sign","Wigner function hadronization matches proton flow reversal","Proton flow sign change from momentum-space pairing","New coalescence rule reproduces proton flow sign jump"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000247,"raw_usage":{"total_tokens":1493,"prompt_tokens":846,"completion_tokens":647,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":462,"completion_tokens_details":{"reasoning_tokens":584}},"tokens_in":462,"tokens_out":647,"duration_ms":6727,"temperature":1.0,"reasoning_tokens":584,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:27:50.864444+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Vary the Gaussian widths in Eqs. (3) and (6) by $\\pm 30\\%$ or change the $g_M/g_B$ ratio while rerunning the hadronization on the same freeze-out parton distribution; if the proton $\\mathrm{d}v_1/\\mathrm{d}y$ sign change between 11.5 and 19.6 GeV or the near-side baryon-antibaryon anticorrelation disappears, the central claim is not robust. Also compare the resulting $v_1$ slope magnitude quantitatively with the measured values, since the paper only claims qualitative agreement.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the string-melting multiphase transport model whose hadronization is being improved."},{"cited_title":"Greco, C","cited_arxiv_id":null,"evidence_quote":"Introduces the dynamical-coalescence Wigner-function method adopted here."},{"cited_title":"Adamczyk et al","cited_arxiv_id":null,"evidence_quote":"Reports the measured proton directed-flow slopes at beam-energy-scan energies that the improved model reproduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows hadronization significantly affects proton directed flow, motivating the improved algorithm."},{"cited_title":"Adam et al","cited_arxiv_id":null,"evidence_quote":"Provides the p+p di-hadron correlations in meson, baryon, and antibaryon channels that the improved model matches."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduced a meson/baryon preference parameter in AMPT coalescence, the direct predecessor of the Wigner-function ranking."},{"cited_title":"Beringer et al","cited_arxiv_id":null,"evidence_quote":"Supplies the pion and proton charge radii that fix the Gaussian widths of the Wigner functions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Independently shows an improved hadronization can reproduce the baryon-antibaryon anticorrelation."}],"review_version":1}