{"id":"4a86723d-1e21-4a1b-8a71-ad4acba7a67e","arxiv_id":"2509.03195","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"PYTHIA and Thermal FIST predict the same proton-deuteron balance scaling but differ sharply in its transverse-momentum dependence, so deuteron-triggered balance functions can in principle tell coalescence from statistical hadronization.","lead":"This paper proposes using deuteron- and proton-triggered balance functions in pp collisions at the LHC to distinguish two competing production mechanisms for light nuclei. The two models make different predictions for how the balance changes with trigger momentum, giving an experimentally testable discriminator that goes beyond just counting nuclei.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thermal FIST's flat pT dependence is built into its uncorrelated sampler, so the claimed mechanism discrimination is not established.","rationale":"The paper is an honest, well-scoped proof-of-concept and the authors explicitly flag that Thermal FIST lacks momentum-space correlations. However, the headline claim that the pT dependence of balance functions discriminates between coalescence and statistical thermal hadronization rests on treating Thermal FIST's flat curve as a thermal-model prediction. Section II shows the flatness is guaranteed by the independent-sampling construction, and Section III confirms that no mechanism couples trigger pT to the balancing-partner distribution. The reader's weakest_assumption identified exactly this point, and I agree it is the load-bearing weakness. The proposed check—adding local momentum-space conservation to the thermal sampler—would determine whether the flatness is robust or an artifact. Pending that check, the verdict should be conditional rather than unconditional acceptance, because the central claim as worded overstates the discriminating power of the observable.","tokens_in":10775,"tokens_out":8069,"duration_ms":82848,"concrete_test":"Re-run the Thermal FIST analysis with local quantum-number conservation imposed in momentum-space subvolumes, using the existing correlation-volume parameter Vc as a momentum-space correlation length, and check whether B(Δy) remains flat in trigger pT. If a pT dependence appears (for example, because high-pT triggers originate from flow cells whose balancing partners are kinematically nearby), the central contrast is an artifact of the uncorrelated sampler rather than a robust thermal-model prediction. An independent, direct check is to measure pT-differential proton- and deuteron-triggered balance functions in ALICE Run 3 pp data; a flat trend would not by itself confirm thermal production unless the correlated thermal model is also tested.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central discriminator (Fig. 5) is the contrast between PYTHIA's narrowing of B(Δy) with trigger pT and Thermal FIST's flat dependence. The flatness is not an emergent prediction of statistical hadronization: Sec. II states that Thermal FIST 'lacks the momentum-space correlations observed in pp collisions' because it 'samples from equilibrium distributions', and Sec. III attributes the flat result to particles being 'produced independently ... with quantum number conservation only imposed globally'. With independent single-particle sampling, the associated yield factorizes with respect to trigger pT, so B(Δy) is pT-independent by construction. The observable therefore tests whether data contain the kind of momentum correlations implemented in PYTHIA's string model, not whether nuclei form by coalescence or thermal hadronization. A thermal/hydrodynamic model with local conservation and flow-induced position-momentum correlations could plausibly produce pT-dependent balance functions, so the paper's conclusion that this observable can discriminate the two production mechanisms goes beyond what the calculation supports.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a new differential observable to distinguish between two production mechanisms of light nuclei in pp collisions: proton- and deuteron-triggered balance functions B(Δy), defined via Eq. (2) as the difference between opposite- and same-baryon-number associated yields. Using PYTHIA 8.3 with deuteron coalescence enabled as a representative of the coalescence picture, and Thermal FIST with a tuned blast-wave boost as a representative of statistical thermal hadronization, the authors compute balance functions for associated antiprotons, Λ baryons, and pions. They find that B_d ≈ 2 B_p in both models, that the deuteron-pion balance vanishes identically in both models, that the Thermal FIST balance function width is controlled by the correlation volume V_c, and, as the central result of Fig. 5, that PYTHIA exhibits a clear narrowing of the balance function with increasing trigger p_T while Thermal FIST shows no p_T dependence. The paper concludes that this observable is a promising discriminator between coalescence and statistical hadronization and motivates a measurement with ALICE Run 3 data.","tokens_in":10970,"tokens_out":6352,"duration_ms":64124,"significance":"If the predicted contrast is robust, the observable would be a valuable new probe of hadronization and nucleus production. The paper has clear strengths: it uses publicly available generators, tests parameter variations (V_c and T), identifies an exact symmetry result for deuteron-pion balance, and proposes an experimentally feasible measurement with existing ALICE capabilities. The principal weakness is that the Thermal FIST p_T independence is essentially built into the model's independent-particle sampling, so the stronger conclusion about discriminating the two production mechanisms needs qualification. The paper is a useful proof-of-concept, but its central interpretive claim goes beyond what the model comparison can establish.","major_comments":[{"comment":"The flat p_T dependence of the Thermal FIST balance functions is a consequence of the model construction, not an emergent prediction of statistical hadronization. The paper itself states in Sec. II that Thermal FIST lacks the momentum-space correlations observed in pp collisions because it samples from equilibrium distributions, and Sec. III attributes the flatness to particles being produced independently with quantum number conservation only imposed globally. With independent single-particle sampling, the associated yield factorizes with respect to the trigger momentum, so B(Δy) is p_T-independent by construction. The contrast in Fig. 5 therefore primarily tests whether the data contain the momentum correlations implemented in PYTHIA's string model, rather than whether nuclei form by coalescence or by thermal hadronization. A thermal or hydrodynamic model that includes local conservation and flow-induced position-momentum correlations could plausibly produce p_T-dependent balance functions. I recommend that the conclusion be reframed to refer to the specific model implementations, and that the claim of discriminating the two production mechanisms be softened unless an additional thermal model with momentum correlations is studied.","section":"Sec. III, Fig. 5 and Sec. II"},{"comment":"No statistical uncertainties or event counts are reported for any of the generator predictions. The central claim is a qualitative comparison of the shapes and their p_T dependence; without statistical errors or at least the number of generated events, it is not possible to assess whether the PYTHIA narrowing and the Thermal FIST flatness are significant relative to Monte Carlo fluctuations. Adding error bars, or a quantitative measure such as the width or RMS of B(Δy) with uncertainties, would also make the predicted experimental discrimination more concrete and would strengthen the proposal for an ALICE measurement.","section":"All figures, especially Fig. 5"}],"minor_comments":[{"comment":"The sentence 'We have tested that the observed small difference agrees with that neutrons are slightly less likely to be balanced by antiprotons than protons' is grammatically unclear; it should read 'agrees with the statement that neutrons are slightly less likely...'.","section":"Sec. III, Fig. 4"},{"comment":"The vanishing deuteron-pion balance is argued in one sentence. Since it is presented as an exact identity, a short formal isospin argument (the deuteron is an isoscalar, so the associated yields of π^+ and π^- are equal) would be more convincing and would clarify the role of baryon-number conservation.","section":"Sec. III, deuteron-pion balance"},{"comment":"The trigger p_T intervals for deuterons are chosen as twice those for protons, but the comparison in Fig. 5 uses these intervals directly. Because coalescence implies p_T(d) ≈ 2 p_T(p), the apparent narrowing pattern for deuterons could be better interpreted if the deuteron intervals were divided by two, or if the choice were explicitly justified in the text.","section":"Sec. III, Fig. 5"},{"comment":"There are minor stylistic inconsistencies, including the hyphenation of Thermal-FIST versus Thermal FIST, and some figure legends omit the variable name in trigger p_T intervals (for example, '1.0 < ... < 8.0 GeV/c'). These should be cleaned up.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely of interest to the journal's audience and the model comparison is clearly presented. The main issue is the overinterpretation of the Thermal FIST result: the flat p_T dependence is built into the uncorrelated sampler, so the observable does not yet establish a discrimination between coalescence and statistical hadronization in general. If the authors qualify the central claim and add statistical information to the figures, the paper would be a solid contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe main thing to know: this is a new, clean observable — deuteron-triggered balance functions — with a genuinely testable prediction. If ALICE Run 3 data show the pT narrowing seen in PYTHIA, that argues for string-like local correlations in hadronization; if flat, it argues against. The paper does that well, and the B_d=2B_p scaling and vanishing d-pion balance are nice cross-checks, both models agree on them.\n\nBut the stress-test note is right, and I think the paper's conclusion overstates it. Sec. II is honest: Thermal FIST samples from equilibrium distributions and only then gets a Blast-Wave boost, so it has no momentum-space correlations by construction. Sec. III says the flat pT dependence 'reflects that particles are produced independently.' That means Fig. 5 is a contrast between a string model with local correlations and an uncorrelated single-particle sampler. It is not a test of statistical hadronization in general. A hydrodynamic thermal model with local conservation and flow-induced correlations could very plausibly produce a pT-dependent balance function. So the paper's claim that this observable discriminates between coalescence and thermal production goes beyond what the calculation supports. The authors flag the cause in the text but don't carry the caveat into the conclusion; they call it a promising discriminator between the two production mechanisms.\n\nThat said, the central observation is still useful. A flat pT dependence in data would be evidence against the kind of color-coherence correlations PYTHIA implements; a narrowing would point toward them. That is a weaker statement but a solid one. The paper would be improved by reframing the conclusion accordingly.\n\nMinor issues: no statistical uncertainties or event counts anywhere, so I can't judge how robust PYTHIA's narrowing is within its own Monte Carlo. The comparison is visual only; a width or ratio metric would sharpen it. The V_c and T scans are helpful but only constrain this particular implementation. The citation pattern is fine — relevant ALICE measurements and the standard phenomenological references are there.\n\nBottom line: a well-scoped phenomenological note, worth sending to a referee. It should be accepted after a revision that softens the interpretation and adds some measure of statistical robustness. I'd bring this to a reading group to discuss the model-comparison logic, but I wouldn't cite it as a definitive discriminator.","headline":"A testable new observable for light-nucleus production that is worth a referee, but the claimed discrimination between thermal and coalescence is weaker than the models actually support because Thermal FIST's flat pT dependence is built in.","tokens_in":11464,"tokens_out":3399,"would_cite":false,"duration_ms":32453,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Proton- and deuteron-triggered balance functions in pp collisions can distinguish deuteron formation by coalescence from formation by statistical thermal hadronization.","keywords":["balance functions","deuteron production","coalescence","statistical hadronization","proton-proton collisions","light nuclei","hadronization mechanism","transverse momentum dependence"],"falsifier":"Measure the width of the proton- and deuteron-triggered balance functions in pp collisions at 13 TeV using LHC Run 3 data across several trigger $p_T$ intervals: if the width narrows with increasing $p_T$, the coalescence/string picture is supported, and if it stays flat, the statistical-thermal picture is supported.","tokens_in":10585,"feed_emoji":"⚛️","tokens_out":8870,"duration_ms":77709,"temperature":0.7,"pith_summary":"This paper proposes a measurement that could decide between the two main explanations for how light nuclei arise in high-energy collisions. It compares proton- and deuteron-triggered balance functions in pp collisions at 13 TeV from a string-fragmentation plus coalescence model and from a statistical thermal model. Both models produce the same simple scaling relation—the deuteron balance function is twice the proton one—and both predict a vanishing deuteron–pion balance. The models disagree sharply when the trigger transverse momentum is varied: the coalescence model shows a narrowing balance function with increasing $p_T$, while the thermal model shows no $p_T$ dependence at all. If confirmed by LHC Run 3 data, this $p_T$ dependence would discriminate the two production scenarios.","feed_headline":"Balance functions' pT narrowing reveals how deuterons form","feed_subtitle":"Narrowing with trigger pT favors coalescence; a flat width favors a thermal fireball.","key_machinery":"The key object is the balance function $B(\\Delta y) = Y_{\\mathrm{opposite}}(\\Delta y) - Y_{\\mathrm{same}}(\\Delta y)$, the difference between the associated yield of particles with opposite baryon number and those with the same baryon number, per trigger particle, as a function of rapidity separation. It isolates the quantum-number-balancing component of hadron correlations. The paper's comparison relies on the relation $B_{\\mathrm{deuteron}}(\\Delta y) \\approx 2 B_{\\mathrm{proton}}(\\Delta y)$, which holds in both models, and on the different response of the width of $B(\\Delta y)$ to trigger $p_T$: the coalescence model produces narrowing while the thermal model stays flat.","core_discovery":"On the paper's own terms, the central discovery is that the transverse-momentum dependence of the balance function separates the two paradigms. In the coalescence-inspired picture, the balancing antibaryon is produced on the same color string as the trigger, so increasing the trigger $p_T$ selects shorter strings and the balance function narrows in rapidity. In the statistical thermal picture, particles are emitted from a common thermal source with baryon number conserved only globally, so the balance function does not depend on trigger $p_T$. This difference is the basis for a new observable that goes beyond nuclei yields, which both models describe equally well.","pith_inferences":["The paper does not test this, but if the $p_T$ narrowing is confirmed, it would suggest that a genuinely thermalized momentum source would erase the narrowing, pushing hadronization pictures toward local color-string balancing rather than global equilibrium.","The deuteron–pion zero balance could serve as a built-in systematic control for detector acceptance and event-mixing corrections; a measured nonzero value would indicate experimental artifacts rather than a new production mechanism.","A testable extension beyond the paper would use heavier clusters such as tritons or helions as triggers, where coalescence predicts balance functions scaling with nucleon number and a correspondingly shifted $p_T$ narrowing, giving a ladder of predictions.","If the discriminating power holds, the observable could also constrain antinucleus production in cosmic-ray sources, since the balance-function shape encodes how antinuclei inherit correlations from their antibaryon parents."],"forward_implications":["A measurement of the $p_T$ dependence of proton- and deuteron-triggered balance functions in existing LHC Run 3 pp data can discriminate coalescence from statistical thermal production of light nuclei.","The scaling relation $B_d \\approx 2 B_p$ holds in both scenarios, so it cannot serve as the discriminator; the $p_T$ dependence is the discriminative handle.","Both models predict a strictly zero deuteron–pion balance, a clean consequence of baryon-number and charge conservation with isospin symmetry that can be checked directly.","The multiplicity dependence seen in the coalescence model is tied to the color-reconnection mechanism; no such dependence appears in the thermal model.","This $p_T$-based discrimination applies to balance functions generally, not only to nuclei, because it probes whether balancing quantum numbers are produced in the same color-coherent process."],"supporting_citations":[{"why":"Supplies the string-fragmentation event generator in which deuteron production via coalescence can be enabled.","marker":"[9]"},{"why":"Supplies the statistical-thermal hadronization package that samples hadrons from an equilibrium distribution with global conservation laws.","marker":"[35]"},{"why":"Provides the empirical cross-section model used to decide which proton-neutron pairs combine into deuterons in the coalescence scenario.","marker":"[34]"},{"why":"Measured Xi-triggered balance functions in pp collisions that demonstrate the observable's sensitivity to hadronization.","marker":"[30]"},{"why":"Earlier comparison of balance functions between the two paradigms; the analysis method here follows it.","marker":"[32]"},{"why":"Data used to tune the blast-wave boost so the thermal model matches measured transverse-momentum spectra in pp collisions.","marker":"[36]"},{"why":"Provides the strangeness-enhancement parameters used for the default correlation volume and temperature in the thermal model.","marker":"[38]"}],"fun_headline_variants":["pT narrowing of balance functions reveals deuteron formation mechanism","Balance function width vs trigger pT distinguishes coalescence and thermal models","Deuteron-triggered balance functions expose production mechanism","Trigger pT dependence of balance functions tells deuteron origin","Narrowing balance functions at high pT favor coalescence over thermal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The flat $p_T$ dependence of the thermal-model balance function is built into the model by construction, because it samples independent particles from an equilibrium source and then applies a tuned boost; if the physical thermal source has intrinsic momentum correlations between a trigger and its balancing partners, the predicted contrast with coalescence would shrink.","fun_headline_variants_meta":{"raw":{"variants":["pT narrowing of balance functions reveals deuteron formation mechanism","Balance function width vs trigger pT distinguishes coalescence and thermal models","Deuteron-triggered balance functions expose production mechanism","Trigger pT dependence of balance functions tells deuteron origin","Narrowing balance functions at high pT favor coalescence over thermal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000624,"raw_usage":{"total_tokens":2848,"prompt_tokens":866,"completion_tokens":1982,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":482,"completion_tokens_details":{"reasoning_tokens":1896}},"tokens_in":482,"tokens_out":1982,"duration_ms":12702,"temperature":1.0,"reasoning_tokens":1896,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:33:12.686139+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the width of the proton- and deuteron-triggered balance functions in pp collisions at 13 TeV using LHC Run 3 data across several trigger $p_T$ intervals: if the width narrows with increasing $p_T$, the coalescence/string picture is supported, and if it stays flat, the statistical-thermal picture is supported.","supporting_citations":[{"cited_title":"An Alternative Formation Model for Antideuterons from Dark Matter","cited_arxiv_id":"1504.07242","evidence_quote":"Provides the empirical cross-section model used to decide which proton-neutron pairs combine into deuterons in the coalescence scenario."},{"cited_title":"Simplifying Strangeness Fluctuations through Balance Functions in Proton-Proton Collisions","cited_arxiv_id":"2506.18375","evidence_quote":"Earlier comparison of balance functions between the two paradigms; the analysis method here follows it."}],"review_version":2}