{"id":"d44457dc-59a7-49e8-b0b0-ee6732101a6b","arxiv_id":"2505.03411","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Forward-backward multiplicity correlations in UrQMD-hydro Au+Au collisions at 10 AGeV are long-range (correlation length above one rapidity unit) and depend on the Cooper-Frye freeze-out prescription.","lead":"Simulations of 10 AGeV gold-gold collisions show that particle counts in forward and backward rapidity windows are correlated, with the correlation weakening as the gap between windows grows. The study also finds that turning off the Cooper-Frye freeze-out step removes the correlation signal, pointing to freeze-out dynamics as the source.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Flat σ² for the no-CF runs does not establish lack of F-B correlation; the paper never shows bcorr(ygap) for that sample, so the central Cooper-Frye conclusion is unsupported.","rationale":"The reader's Gaussian-ansatz concern is legitimate but secondary: it affects the quantitative λ>1 LRC claim, whereas the paper's distinctive new statement—that Cooper-Frye freeze-out is necessary to generate F-B correlations—rests entirely on the σ² plot for the no-CF sample. As printed, Eq. (5) is dimensionally inconsistent, so the plotted σ² must come from an unstated corrected formula. More importantly, σ² measures the scaled variance of nF−nB, not the correlation strength bcorr; flatness in ygap does not imply vanishing correlation. Since bcorr(ygap) is never shown for the no-CF run, the decisive comparison is missing. A simple re-analysis of the stored events would settle the issue. I keep the reader's conditional verdict, but for a different reason: the paper needs the bcorr(ygap) curve for the no-CF sample before the central claim can be accepted.","tokens_in":4240,"tokens_out":10004,"duration_ms":103055,"concrete_test":"Re-analyze the stored 60k UrQMD-hydro-without-CF events: compute bcorr(ygap) using Eq. (2) with the same window width and gap definition as the other samples, and compute σ² with the corrected variance-of-difference formula. If bcorr is not flat in ygap or is not consistent with zero, the 'lack of correlation' statement is false and the Cooper-Frye conclusion fails. If bcorr is flat and the paper reports the exact formula and window width used, the conclusion is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing weakness is not the Gaussian fit but the interpretation of σ², which carries the paper's distinctive claim about Cooper-Frye freeze-out. In §4 and Fig. 5, the no-CF run is declared to show 'lack of correlation' solely because σ² is flat versus ygap. Even after correcting the apparent typo in Eq. (5) (as printed it contains V² terms and cannot yield the plotted values of order 1–10; the standard variance-of-difference form is σ²=(V_FF+V_BB−2V_FB)/⟨nF+nB⟩), this observable is the scaled variance of the multiplicity difference nF−nB, not the F-B correlation strength bcorr. A flat σ² is compatible with a nonzero, gap-independent bcorr, and conversely bcorr can fall with ygap while σ² stays flat. The authors never show bcorr(ygap) for the no-CF sample, so the conclusion that Cooper-Frye freeze-out is necessary for generating F-B correlations is unsupported. The λ>1 LRC claim also rests on the self-cited Gaussian ansatz, but the CF claim is the more consequential and less supported step.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies forward-backward (F-B) multiplicity correlations in Au+Au collisions at 10 AGeV using three variants of UrQMD-generated events: default UrQMD, UrQMD with an intermediate hydrodynamic stage (UrQMD-hydro), and UrQMD-hydro with the Cooper-Frye (C-F) freeze-out step disabled. The authors define the correlation strength bcorr from the slope of ⟨nB⟩ versus ⟨nF⟩ and from the covariance/variance ratio, fit bcorr as a function of the rapidity gap ygap to a Gaussian ansatz exp(−ygap^2/δ^2), and convert the fitted width δ into a correlation length λ = 2δ/√π. They report λ > 1 for both UrQMD-default and UrQMD-hydro, interpret σ² > 1 versus ygap as evidence for long-range correlations, and observe that σ² is flat for the no-CF sample, which they interpret as a lack of correlation and hence as evidence that the Cooper-Frye mechanism is necessary for generating F-B correlations.","tokens_in":4509,"tokens_out":3676,"duration_ms":38141,"significance":"If established, the central claims would be of interest to the heavy-ion community planning SIS100 measurements, because they provide a concrete model prediction for how F-B correlations and their correlation length behave at lower beam energies and how the hydrodynamic freeze-out stage affects these observables. The paper's strength is that it compares three model variants with fixed collision parameters and presents a quantitative observable that could be checked against future SIS100 data; the analysis is also straightforward to reproduce once the model settings are specified. However, the key quantitative conclusions rest on an unvalidated Gaussian parameterization, a dimensionally inconsistent formula for σ², and an interpretive step that equates a flat σ² with the absence of F-B correlation. These issues currently prevent the results from being regarded as established. The paper does not provide machine-checked proofs or released analysis code; its contribution is an exploratory model comparison.","major_comments":[{"comment":"Equation (5) as printed is dimensionally inconsistent: if V_FF, V_BB, and V_FB are variances/covariances of multiplicities (dimension n^2), the numerator contains n^4 terms divided by ⟨n_F + n_B⟩ (dimension n), which cannot yield the order-unity values plotted in Figs. 4 and 5. The standard variance-of-difference form is σ² = (V_FF + V_BB − 2 V_FB)/⟨n_F + n_B⟩, and this is presumably what was computed. Because the σ² > 1 interpretation and the flat-σ² conclusion for the no-CF sample depend on this quantity, the formula must be corrected and the figures re-checked against the corrected expression.","section":"Sec. 2, Eq. (5)"},{"comment":"A flat σ² as a function of ygap does not imply the absence of F-B correlation. σ² is a scaled variance of the forward-backward multiplicity difference, not the correlation strength bcorr; a ygap-independent σ² is compatible with a nonzero and even ygap-independent bcorr, and conversely a falling bcorr can coexist with a flat σ². The paper never shows bcorr(ygap) for the UrQMD-hydro without C-F sample, so the conclusion that the Cooper-Frye approach is necessary for generating F-B correlations is unsupported by the presented evidence.","section":"Sec. 4, Fig. 5"},{"comment":"The Gaussian parameterization bcorr ∝ exp(−ygap^2/δ^2) and the conversion λ = 2δ/√π are attributed only to the authors' own proceedings reference [3] without derivation or external validation. The reported correlation lengths λ ≈ 1.35–1.50 are therefore fit parameters of an assumed functional form; if the true decay were exponential or otherwise non-Gaussian, the extracted λ values and the statement that λ > 1 indicates long-range correlation would lose their meaning. The authors should either derive the Gaussian form from a model, compare it with alternative fits, or benchmark the extracted λ against a known observable with an established source of correlation.","section":"Sec. 2, Eqs. (3)-(4) and Table 1"},{"comment":"The manuscript never states the value of the rapidity window width δy used in the analysis, although the definition of ygap and the placement of the forward and backward windows depend on it; without this parameter the results in Table 1 and Figs. 3-5 cannot be reproduced or compared with other analyses. Additionally, the figures show no statistical uncertainties or systematic variations, so the claimed differences between UrQMD-default and UrQMD-hydro (and the flatness of the no-CF σ²) cannot be quantitatively assessed.","section":"Sec. 2 and Figs. 3-5"}],"minor_comments":[{"comment":"The sentence beginning \"Let consider the first moment\" should read \"Let us consider\" and the notation ⟨n_B⟩ n_F in the same section is confusing and should be cleaned up.","section":"Sec. 1, Introduction"},{"comment":"The regression form ⟨n_B⟩ = a + bcorr ⟨n_F⟩ should be stated with the random variable n_F (or its conditional average) made explicit, since Eq. (2) then defines bcorr through the covariance/variance ratio; a brief derivation would clarify why the two definitions coincide.","section":"Sec. 2, Eq. (1)"},{"comment":"Typo: \"UrQMD-defaut\" should be \"UrQMD-default\"; also the axis label \"corrb\" should be \"b_corr\".","section":"Fig. 3 caption"},{"comment":"Several typographical errors occur, including \"eqaution\", \"diffrent\", and the inconsistent use of \"C-F\" versus \"CF\" for the Cooper-Frye approach; these should be corrected.","section":"Throughout"},{"comment":"The color-bar scale in Fig. 2 (event counts up to 24000) makes the scatter plot difficult to read; a density plot with a clear color scale or contours would better support the claim of a forward-backward correlation.","section":"Sec. 4, Fig. 2"},{"comment":"The statement that λ > 1 means \"particles are still correlated at the rapidity gap more than one unit\" is imprecise, since λ is a parameter of a Gaussian fit and not an operator-defined threshold; the wording should be revised to avoid implying a direct measurement of a physical length.","section":"Sec. 5, Summary"}],"recommendation":"major_revision","confidential_remarks":"The paper reads as a extended proceedings contribution. The main results hinge on self-cited formulas and on a claim about Cooper-Frye that is not supported by the shown observable. I would ask for a substantive revision rather than a quick acceptance. If the authors can supply the corrected σ² formula, show bcorr for the no-CF sample, and justify or replace the Gaussian λ extraction, the work could become a useful model comparison for the SIS100 program."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a small model study with one clean qualitative result and one unsupported headline claim. Read it for the UrQMD-hydro comparison; don't believe the Cooper-Frye conclusion as written.\n\nWhat's new: the group compares bcorr and sigma2 for UrQMD-default vs UrQMD-hydro at 10 AGeV Au+Au, and adds a third run with Cooper-Frye disabled. The bcorr(ygap) trend is clean: monotonically decreasing, higher for hydro. That's a legitimate extension of ref. [6]'s UrQMD-only study, and the SIS100 community can use it as a baseline prediction.\n\nCredit where due: they generated 150k/100k/60k events, show the raw 2D correlation, and give fit parameters with errors. That's reproducible enough for a model study.\n\nSoft spots, in order of seriousness:\n\n1. The Cooper-Frye claim. The paper says the no-CF run shows 'lack of correlation' because sigma2 is flat vs ygap. But sigma2, even after fixing the typo in Eq. (5), is the scaled variance of the multiplicity difference, not bcorr. A flat sigma2 doesn't imply zero F-B correlation; bcorr could be nonzero and constant. They never show bcorr(ygap) for the no-CF sample. So the central conclusion – that Cooper-Frye freeze-out is necessary for F-B correlations – is unsupported. That's not a minor caveat; it's the paper's distinctive claim.\n\n2. The Gaussian ansatz. Eqs. (3)-(4) come from their own proceedings ref. [3], with no derivation or external motivation. The reported lambda ~1.35-1.50 are fit parameters of that assumed shape. If the decay is not Gaussian, lambda loses meaning. The 'lambda>1 implies LRC' step is therefore conditional.\n\n3. Missing settings and a typo. The rapidity window width is never stated. Eq. (5) as printed has squared variances/covariances, which is dimensionally wrong and can't reproduce the plotted values. I assume it's a typo for the standard form, but it should be fixed. Figures 3-5 have no statistical uncertainties on the points, which matters when interpreting flatness.\n\nThe stress-test note lands. The bcorr result itself is fine; it's the interpretation of sigma2 that doesn't hold up.\n\nWho this is for: people working on FB correlations at SIS100 energies, and anyone using UrQMD-hydro. The paper is short, clear, and could be a useful proceedings-level contribution after major revision.\n\nRecommendation: send to peer review, but the referees should push for a bcorr(ygap) plot for the no-CF run, a statement of the rapidity window, and a derivation or external reference for the Gaussian form. If the authors can't show bcorr for the no-CF sample, they need to soften the Cooper-Frye claim.","headline":"A clean qualitative model study whose central statistical claim about Cooper-Frye freeze-out does not follow from the plotted observable.","tokens_in":5049,"tokens_out":2518,"would_cite":false,"duration_ms":24035,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q"],"model":"deepseek-v4-flash","headline":"In 10 AGeV Au+Au collisions, a transport-plus-hydro hybrid simulation produces forward-backward multiplicity correlations that survive rapidity gaps above one unit, and the Cooper-Frye freeze-out step is needed to generate them.","keywords":["forward-backward correlation","long-range correlation","UrQMD-hydro","dynamical fluctuations","rapidity gap","Cooper-Frye freeze-out","correlation length","Au+Au collisions at 10 AGeV"],"falsifier":"Refit the same $b_{\\mathrm{corr}}$ versus $y_{\\mathrm{gap}}$ points with an exponential and with a power law using the same number of parameters; if either fits comparably or better, the extracted $\\lambda$ is a shape artifact. A second, model-independent check is to measure the rapidity gap at which $b_{\\mathrm{corr}}$ crosses a fixed small threshold without assuming any functional form; if that crossing sits below one unit of rapidity, the $\\lambda > 1$ long-range conclusion collapses.","tokens_in":4016,"feed_emoji":"⚛️","tokens_out":9633,"duration_ms":81840,"temperature":0.7,"pith_summary":"This paper tries to establish that hybrid transport-plus-hydro simulations of 10 AGeV Au+Au collisions contain long-range forward-backward multiplicity correlations, not just short-range noise, and that the hydrodynamic stage strengthens them relative to the pure transport calculation. If true, collisions at fixed-target energies around 10 AGeV should show measurable correlations between particles emitted in forward and backward rapidity windows separated by more than one unit of rapidity, with the correlations tracing back to the collective evolution of the fireball. The paper further claims that removing the Cooper-Frye freeze-out step makes the fluctuation measure $\\sigma^2$ flat in the rapidity gap, so the freeze-out procedure is what converts hydrodynamic flow into correlated final-state particles. A reader should care because this identifies a specific dynamical mechanism — not just statistical noise — as the source of event-by-event fluctuations at energies planned for new fixed-target experiments.","feed_headline":"Rapidity correlations survive a one-unit gap at 10 AGeV","feed_subtitle":"Hybrid simulation finds long-range correlations only when Cooper-Frye freeze-out is on.","key_machinery":"The load-bearing quantities are the forward-backward correlation strength, defined as the covariance of forward and backward multiplicities divided by the forward multiplicity variance, and the assumed Gaussian decay $b_{\\mathrm{corr}} \\propto \\exp(-y_{\\mathrm{gap}}^2/\\delta^2)$ with the conversion $\\lambda = 2\\delta/\\sqrt{\\pi}$. A second observable, $\\sigma^2$, is built from the forward and backward variances and their covariance and is compared against the Poisson baseline of 1 to classify correlations as short- or long-range. The argument is carried by comparing three simulated data sets — default transport, hybrid with Cooper-Frye freeze-out, and hybrid with freeze-out disabled — where the third run isolates the freeze-out mechanism.","core_discovery":"The central claim is that for minimum-bias Au+Au at 10 AGeV, the hybrid UrQMD-hydro model yields a forward-backward correlation strength $b_{\\mathrm{corr}}$ that decreases monotonically with the rapidity gap, and the fitted Gaussian parameter $\\delta$ corresponds to correlation lengths $\\lambda \\approx 1.35$ for the hybrid run and $\\lambda \\approx 1.50$ for the default transport run. Because both extracted lengths exceed unity in rapidity units, the paper reads them as the signature of long-range correlation: particles separated by more than one unit of rapidity are still correlated. The companion result is that the hybrid run with the Cooper-Frye particle-production step disabled shows $\\sigma^2$ flat against $y_{\\mathrm{gap}}$, which the paper interprets as an absence of correlation, implying that the Cooper-Frye freeze-out stage is the active ingredient producing the observed long-range correlations in the hybrid calculation.","pith_inferences":["Editorial inference: The Gaussian-derived $\\lambda$ could be checked against a model-independent length scale, such as the rapidity gap at which $b_{\\mathrm{corr}}$ falls to half its peak; a large disagreement would show that the functional-form assumption, not the data, drives the $\\lambda > 1$ conclusion.","Editorial inference: Because the paper uses only minimum-bias events, a natural next test is centrality-binned runs; if the correlation strengthens with centrality as it does at higher beam energies, the 10 AGeV result would connect cleanly to the established energy pattern.","Editorial inference: The flat $\\sigma^2$ in the run without Cooper-Frye could reflect the absence of collective velocity smearing rather than a true loss of correlation; computing $b_{\\mathrm{corr}}$ directly for that run, rather than only $\\sigma^2$, would separate the two explanations.","Editorial inference: A direct extension would be to scan beam energies from 10 to 40 AGeV in the same hybrid setup and test whether the extracted correlation length grows with collision energy; this paper reports a single energy."],"forward_implications":["At fixed-target energies near 10 AGeV, minimum-bias Au+Au events should show a measurable forward-backward correlation that survives rapidity gaps above one unit, so detector acceptances should be designed to cover forward and backward windows with such gaps.","The hydrodynamic stage raises $b_{\\mathrm{corr}}$ relative to the default transport calculation, making the correlation strength a model discriminator sensitive to the equation of state and collective expansion.","Switching off the Cooper-Frye freeze-out flattens $\\sigma^2$ versus $y_{\\mathrm{gap}}$, so the freeze-out treatment is not merely a detail of particle spectra; it determines whether long-range multiplicity correlations appear in the hybrid calculation at all.","The $\\sigma^2 > 1$ values seen in the two runs with active freeze-out indicate deviations from Poisson emission across all measured gaps, so independent single-particle production cannot describe these events by itself."],"supporting_citations":[{"why":"Supplies the Gaussian parameterization of $b_{\\mathrm{corr}}$ versus $y_{\\mathrm{gap}}$ and the relation $\\lambda = 2\\delta/\\sqrt{\\pi}$ used to extract correlation lengths.","marker":"[3]"},{"why":"Provides the prior UrQMD forward-backward correlation and fluctuation calculation at the same fixed-target energies that this work extends with a hydrodynamic stage.","marker":"[6]"},{"why":"Gives the experimental measurement and the covariance-to-variance definition of $b_{\\mathrm{corr}}$ that the analysis adopts.","marker":"[2]"},{"why":"Establishes long-range forward-backward multiplicity correlations at high energy, the benchmark that motivates looking for the same effect at fixed-target energies.","marker":"[5]"},{"why":"Defines the microscopic transport model used to generate the default data set.","marker":"[7]"},{"why":"Describes the integrated transport-plus-hydrodynamic approach used to generate the hybrid UrQMD-hydro data.","marker":"[10]"},{"why":"Supplies the chiral equation-of-state hydrodynamics with Cooper-Frye freeze-out used in the hybrid runs.","marker":"[12]"},{"why":"Documents the UrQMD version details needed to reproduce the default simulations.","marker":"[9]"}],"fun_headline_variants":["Long-range rapidity correlations vanish without Cooper-Frye at 10 GeV","Forward-backward correlation length ~1.35 rapidity units in hybrid run","Hybrid model needs Cooper-Frye to see correlations across rapidity gap","Gap-dependent forward-backward fluctuations at 10 AGeV Au+Au","Cooper-Frye freeze-out is key for long-range rapidity correlations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire long-range-correlation conclusion rests on assuming that the correlation strength falls off as a bell-shaped Gaussian in the rapidity gap, an assumption quoted from the authors' own earlier work without derivation; if the true falloff is exponential or some other shape, the reported correlation lengths and the 'greater than one means long-range' reading would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Long-range rapidity correlations vanish without Cooper-Frye at 10 GeV","Forward-backward correlation length ~1.35 rapidity units in hybrid run","Hybrid model needs Cooper-Frye to see correlations across rapidity gap","Gap-dependent forward-backward fluctuations at 10 AGeV Au+Au","Cooper-Frye freeze-out is key for long-range rapidity correlations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000748,"raw_usage":{"total_tokens":3273,"prompt_tokens":828,"completion_tokens":2445,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":444,"completion_tokens_details":{"reasoning_tokens":2344}},"tokens_in":444,"tokens_out":2445,"duration_ms":17677,"temperature":1.0,"reasoning_tokens":2344,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:52:01.014438+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Refit the same $b_{\\mathrm{corr}}$ versus $y_{\\mathrm{gap}}$ points with an exponential and with a power law using the same number of parameters; if either fits comparably or better, the extracted $\\lambda$ is a shape artifact. A second, model-independent check is to measure the rapidity gap at which $b_{\\mathrm{corr}}$ crosses a fixed small threshold without assuming any functional form; if that crossing sits below one unit of rapidity, the $\\lambda > 1$ long-range conclusion collapses.","supporting_citations":[{"cited_title":"Gope et al., Investigation of forward-backward correlation in hybrid UrQMD-hydro at 10 AGeV Au+Au collisions, Proceedings of the DAE Symp","cited_arxiv_id":null,"evidence_quote":"Supplies the Gaussian parameterization of $b_{\\mathrm{corr}}$ versus $y_{\\mathrm{gap}}$ and the relation $\\lambda = 2\\delta/\\sqrt{\\pi}$ used to extract correlation lengths."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the prior UrQMD forward-backward correlation and fluctuation calculation at the same fixed-target energies that this work extends with a hydrodynamic stage."},{"cited_title":"Adam et al","cited_arxiv_id":null,"evidence_quote":"Gives the experimental measurement and the covariance-to-variance definition of $b_{\\mathrm{corr}}$ that the analysis adopts."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes long-range forward-backward multiplicity correlations at high energy, the benchmark that motivates looking for the same effect at fixed-target energies."},{"cited_title":"Bass et al., Microscopic models for ultrarelativistic heavy ion colli- sions, Prog.Part.Nucl.Phys","cited_arxiv_id":null,"evidence_quote":"Defines the microscopic transport model used to generate the default data set."},{"cited_title":"Steinheimer et al., Hydrodynamics with a chiral hadronic equation of state including quark degrees of freedom, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the chiral equation-of-state hydrodynamics with Cooper-Frye freeze-out used in the hybrid runs."}],"review_version":1}