{"id":"66a34a25-dc9b-4c42-85d7-132142397622","arxiv_id":"2510.24507","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"In an effective chiral model retaining only the zero field mode, finite volume smooths the transition, regulates critical fluctuations, and shifts the baryon-kurtosis peak to lower collision energies.","lead":"This paper studies how the small size of a heavy-ion fireball changes predicted signals for a possible critical point in the QCD phase diagram. It finds that finite-size effects can shift the phase boundary and reshape fluctuation observables, so a kurtosis peak along the freeze-out line does not by itself prove a critical point.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The zero-mode-only reduction places all size dependence into the βV factor; omitted non-zero modes and gradient terms could control the L≈10–40 fm thresholds and the R42 shift, so the model's quantitative claims are not yet secured.","rationale":"The reader identified the same load-bearing assumption: the entire finite-size effect enters only through βV multiplying a size-independent Ueff, with non-zero modes, momentum discretization, vacuum fluctuations, and gradient terms omitted. My read agrees with that assessment. The concern is genuinely load-bearing because the paper's own quantitative statements—the L values at which the phase boundary changes, the L≈40 fm sensitivity of R42, and the shift of the dip-peak structure—are all controlled by this simplification. That said, the paper repeatedly frames itself as a qualitative illustration, explicitly states the omission in Sec. II A, and the internal derivations, including the double-Gaussian coexistence analysis and the Binder-cumulant construction, are internally consistent. The first-order coexistence scaling χ∝L^d is a robust consequence of two-state fluctuations and would survive in a more complete treatment. For these reasons I would not move the verdict to REJECT or UNVERDICTED. The appropriate action remains CONDITIONAL: the manuscript should either correct the over-promising abstract or include a quantitative test of the zero-mode-only truncation before its thresholds and phenomenological shifts are cited as finite-volume predictions. No code or data are provided, but the methods are sufficiently detailed that the proposed momentum-mode truncation check can be implemented directly.","tokens_in":18005,"tokens_out":9627,"duration_ms":96188,"concrete_test":"Recompute ⟨φ⟩, the chiral susceptibility, and R42 in a cubic box of side L using a controlled truncation that keeps the zero mode plus the N lowest non-zero Fourier modes, with the gradient term (2π n/L)²|φ_n|² and the fermionic determinant evaluated with the same discrete momenta as in Ref. [37]. Run this for L = 2, 4, 6, 10, 20, 40, and 100 fm and compare the phase-boundary shift and the R42 dip-peak location with the zero-mode-only results. If including just the first few momentum modes changes these quantities by more than about 20% at L ≤ 40 fm, the central simplification is not robust and the reported thresholds are model artifacts.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central construction is Eq. (4), Z = ∫ dφ exp(−βV Ueff(φ)), obtained by keeping only a single spacetime-independent mode. Section II A explicitly says: 'we will completely omit the momentum space constraints for simplicity, which eliminates the size dependence of Ueff.' Thus every finite-size effect in the paper—rounding of the transition, the reported phase-boundary shift near L≈10–20 fm, the R42 sensitivity below L≈40 fm, and the dip-peak shift in Fig. 12—comes solely from the multiplicative factor βV in front of a size-independent effective potential. The abstract, by contrast, advertises 'momentum-space discretization, and gradient effects modeled via a prescribed finite-volume profile,' which never appears in the body. If the discarded non-zero modes contribute at these physically relevant sizes—for instance through a gradient cost of order (2π/L)² or through modified momentum sums in the fermionic matter term of Eq. (10)—then the thresholds and the shift are not robust predictions of a finite-volume model but properties of the chosen zero-mode ansatz. The paper's assertion that including momentum constraints does not modify the scaling behavior is neither shown nor accompanied by data, and it does not address the R42 shift or the phase-boundary comparison with lattice results in Fig. 7.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a finite-volume extension of a quark-meson-like effective model in which only the spacetime-independent zero mode of the scalar field is retained in the partition function, giving Z = ∫ dφ exp(−βV U_eff(φ)). Finite-size effects then enter through the βV prefactor multiplying a size-independent effective potential. The authors study the chiral condensate, chiral susceptibilities, Binder cumulant, and net-baryon number cumulant ratio R42 as functions of system size L. They report finite-size scaling collapse near the critical endpoint (using an effective exponent \\tildeν = 2/3), an L^d scaling of the susceptibility at the first-order transition due to phase coexistence, a size-independent Binder cumulant crossing, and a finite-size shift of the R42 dip-peak along a parametrized freeze-out line. The phase boundary is also compared with lattice results obtained by extrapolation from imaginary chemical potentials.","tokens_in":18530,"tokens_out":23189,"duration_ms":190297,"significance":"If accepted as a qualitative finite-volume model, the construction is useful: it gives a transparent analytic mechanism for finite-size rounding, a double-Gaussian derivation of coexistence scaling, and an illustration of why susceptibility ratios such as R42 can be volume-dependent near a critical point. The paper honestly states that vacuum fluctuations and momentum-space constraints are neglected. Its main value is as a controlled toy model for interpreting finite-size effects in heavy-ion phenomenology. However, the quantitative thresholds (L ≈ 10–20 fm for the phase boundary, L ≈ 40 fm for R42) are properties of the zero-mode ansatz and are not secured against the omitted non-zero modes and gradient terms; the abstract, as supplied for review, overstates the scope.","major_comments":[{"comment":"The abstract supplied for this manuscript promises momentum-space discretization and gradient effects via a prescribed finite-volume profile, but the body states the opposite: 'we will completely omit the momentum space constraints for simplicity, which eliminates the size dependence of Ueff.' Consequently every finite-size effect in the paper—the rounding of the transition, the L ≲ 10–20 fm phase-boundary shift in Fig. 7, and the R42 dip-peak shift in Fig. 12—comes solely from the multiplicative factor βV in Eq. (4) in front of a size-independent Ueff. The claim in Sec. III A that including momentum constraints 'does not modify the scaling behavior' is not supported by any displayed numerical check and does not address the sub-20 fm thresholds or the R42 shift. Please either implement and quantify the promised momentum/gradient treatment, or revise the abstract and present the quantitative thresholds as properties of the zero-mode model only.","section":"Sec. II A, Eq. (4)"},{"comment":"The statement 'generally χ_k ∝ L^{kd}' is not correct for the susceptibility definitions in Eq. (A2). From χ_k = (βV)^{k−1} κ_k and the two-peaked coexistence distribution with κ_k = O(1), one obtains χ_k ∝ L^{d(k−1)}; for example, χ_2 ∝ L^d (as the same sentence states) and χ_4 ∝ L^{3d}. The later statement in Sec. IV that R42 scales as L^6 at the first-order transition is consistent with the corrected law, so the displayed power L^{kd} appears to be a factor-L^d typo, but it should be fixed because it is an explicit quantitative claim used in the discussion of the first-order peak.","section":"Sec. III A"},{"comment":"The text states that the finite-size phase boundaries 'match surprisingly well' with the lattice results of Refs. [19,20]. This is a quantitative-sounding claim made with a model that keeps only the zero mode, neglects vacuum fluctuations, and omits momentum-space constraints; the authors themselves note in Sec. III A that the direction and magnitude of the vacuum contribution depend on the type of constraint, boundary conditions, and treatment. Please either provide a quantitative comparison (e.g., a band of lattice results with a definition of agreement) or explicitly downgrade this statement to 'qualitative similarity' throughout the text and the conclusion.","section":"Sec. III C, Fig. 7"}],"minor_comments":[{"comment":"Once the momentum-space constraints are omitted, U_eff no longer depends on L; the L argument in U_eff(φ,T,μ,L) is misleading and should be removed.","section":"Sec. II A, Eq. (8)"},{"comment":"The factor N_c is never defined in the text; please state that N_c = 3 is used.","section":"Eq. (10)"},{"comment":"Please specify the units of α and β in the quartic freeze-out parametrization; from the subsequent text ('fixing β = 0.06') one infers they are dimensionful, but the dimensions should be stated explicitly.","section":"Eq. (25)"},{"comment":"The statement that \\tildeν is introduced 'merely to shorten our notation' understates its role; \\tildeν = (γ + 2β)/d is the effective finite-size exponent required when hyperscaling is violated in the mean-field approximation, and this should be said more precisely.","section":"Sec. II B"},{"comment":"The crossing of the Binder cumulant curves 'recovers the L→∞ CEP' because τ is defined with respect to that same CEP; the sentence should be phrased as a consistency check rather than as an independent determination of the CEP.","section":"Sec. III E, Fig. 9 inset"},{"comment":"The title in the manuscript header differs from the title in the submission metadata; please ensure they match. There are also numerous typographical errors (e.g., 'resuts', 'comming', 'te minima', 'sligthly') and a careful proofreading pass is needed.","section":"Title and typos"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the zero-mode construction is a reasonable qualitative contribution. The main issue is the gap between the advertised scope and the actual calculation, plus the unsupported assertion that momentum constraints do not affect the results. The L^{kd} scaling statement is a clear error but easily corrected. I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nThe paper is worth reading for one idea: near a critical endpoint in a finite system, a non-monotonic kurtosis along the freeze-out line is not a clean CEP signature, and finite-size rounding sets in around L≈10–40 fm. The authors get there by promoting the usual mean-field minimization to an integral over one homogeneous mode, Z=∫dφ exp(−βV U_eff(φ)). That zero-mode trick is not new (they credit Ref. [43]), but applying it to a quark-meson-type model and mapping the Binder cumulant and R42 freeze-out behavior is a concrete, useful exercise.\n\nWhat they do well: the finite-size scaling collapse for the chiral condensate and susceptibility is a genuine check, and the coexistence-driven L^d susceptibility at the first-order transition is derived analytically via a double-Gaussian approximation. They also state limitations honestly: the vacuum fluctuations are dropped, momentum constraints are omitted, and the model is called schematic. The Binder cumulant crossing recovers the infinite-volume CEP as a consistency check, not a fit.\n\nThe soft spot is real and has two parts. First, the abstract promises momentum-space discretization and gradient effects that never appear in the body; Section II A says the opposite—momentum constraints are 'completely omitted.' That mismatch should be fixed before publication. Second, because U_eff is size-independent, every finite-size effect in the paper, including the R42 dip-peak shift and the L≈10–40 fm thresholds, comes solely from the βV factor. If non-zero modes or gradient terms contribute at those sizes, the thresholds are not robust. The authors assert that adding momentum constraints does not change the scaling behavior, but they don't show the check, and it doesn't address the R42 shift. This is a limitation, not a fatal flaw: the paper is explicitly qualitative, and the zero-mode truncation is a legitimate way to regulate the critical divergence. But the quantitative thresholds should be advertised as model-dependent.\n\nVerdict: the central caution—dip-peak structure does not imply CEP proximity—holds up as a qualitative statement. The Binder cumulant analysis is solid. The paper deserves a serious referee, with a request to align the abstract with the body and to state the zero-mode-only caveat in the conclusions. I'd bring it to a reading group but wouldn't cite the quantitative thresholds without qualification.","headline":"A useful qualitative caution about non-monotonic kurtosis and finite-size rounding, but the abstract overpromises and the quantitative thresholds rest on a single-mode truncation.","tokens_in":18835,"tokens_out":2693,"would_cite":true,"duration_ms":23134,"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":"The paper claims that in an effective quark-meson model, keeping the integration over a single constant order-parameter mode rather than imposing momentum-space constraints carries the finite-size physics, rounding the transition…","keywords":["finite-size effects","chiral order parameter","critical endpoint","net-baryon number fluctuations","Binder cumulant","kurtosis ratio","QCD phase diagram","quark-meson model"],"falsifier":"Recompute the central observables with the non-zero modes restored (finite-box fermion-mode sums or a low-momentum cutoff in the matter potential, plus the scalar gradient term) at $L \\approx 10$ fm and $T$ near the would-be critical endpoint; the paper's own premise is that these are negligible, so if the chiral susceptibility or the $R_{42}$ dip-peak position shifts by more than the quoted few-percent level, the central claim fails. Alternatively, a lattice QCD calculation of the chiral susceptibility at finite volume near the critical point should show the predicted rounding with $\\chi_2 \\propto L^{3/2}$; a persistent divergence or a different exponent would falsify the framework.","tokens_in":17810,"feed_emoji":"⚛️","tokens_out":9864,"duration_ms":85823,"temperature":0.7,"pith_summary":"The paper proposes that finite-size effects in effective models of the QCD phase diagram can be captured not by modifying the momentum spectrum but by keeping the functional integral over a single spacetime-independent order-parameter mode: $Z = \\int d\\varphi \\, e^{-\\beta V U_{\\mathrm{eff}}(\\varphi,T,\\mu)}$. Because the factor $\\beta V$ multiplies a size-independent effective potential, the integral acts as a statistical average over mean-field configurations, rounding the transition at finite volume and producing the correct finite-size scaling at the critical point. The authors apply this to a quark-meson-like model and show that at a first-order transition the coexistence of two minima makes susceptibilities grow as $L^d$, so the maximum of the susceptibility cannot locate the critical endpoint; the Binder cumulant remains size-independent and can. They then compute the net-baryon kurtosis ratio $R_{42} = \\chi_4^B/\\chi_2^B$ along the freeze-out line and find that for $L \\lesssim 10$ fm the dip-peak structure shifts to lower collision energies, while for $L \\gtrsim 40$ fm it is indistinguishable from the infinite-volume result. A sympathetic reader would care because heavy-ion fireballs are only a few to tens of fermis across, and this gives a qualitative, controlled way to estimate what those finite sizes do to fluctuation observables used in critical-endpoint searches.","feed_headline":"One zero-mode integral reproduces finite-size QCD scaling","feed_subtitle":"Volume-weighted averaging over the order parameter rounds the transition and shifts fluctuation peaks at fireball sizes.","key_machinery":"The machinery is the zero-mode, volume-weighted partition function $Z = \\int d\\varphi \\, e^{-\\beta V U_{\\mathrm{eff}}(\\varphi,T,\\mu)}$, built from a size-independent effective potential for a quark-meson-like model, with all finite-size dependence entering only through the factor $\\beta V = L^3/T$. The integral turns the mean-field free-energy landscape into a probability distribution over the constant field, so a double-well potential becomes a two-peaked probability density whose relative weights are controlled by $V$; this is what produces the melting of the transition at small $L$ and the coexistence-driven $L^d$ enhancement of susceptibilities at a first-order boundary. The double-Gaussian approximation for that density supplies the analytic $\\chi(\\eta)$ formula used to subtract the coexistence peak and verify the scaling. The Binder cumulant $\\kappa_B = \\chi_4/(\\beta L^d \\chi_2^2)$, being volume independent by construction, provides the size-independent diagnostic that locates the infinite-volume critical endpoint.","core_discovery":"The central claim is that the volume dependence of the grand-canonical partition function itself, rather than the discretization of momenta, is the essential finite-size effect in a mean-field chiral model. Keeping only the zero mode $\\varphi$ of the scalar field, the authors write $Z = \\int_{-\\infty}^{\\infty} d\\varphi \\, e^{-\\beta V U_{\\mathrm{eff}}(\\varphi,T,\\mu)}$ with $U_{\\mathrm{eff}}$ independent of $L$. In the infinite-volume limit the weight becomes a sum of Dirac deltas at the minima, recovering ordinary mean-field theory; at finite $L$ the peaks melt. This yields a rounded crossover instead of a sharp first-order transition, a decrease of the pseudocritical temperature at small $L$ (about 6.5 percent at $L = 3$ fm, in line with lattice extrapolations), and controlled scaling: near the critical point the singular parts collapse with an effective exponent $\\tilde{\\nu} = 2/3$, i.e. $\\chi_2 \\propto L^{3/2}$, while at the first-order transition the two-peaked probability distribution gives $\\chi_k \\propto L^{kd}$. Along the freeze-out line, $R_{42} = \\chi_4^B/\\chi_2^B$ is unchanged for $L \\gtrsim 40$ fm and becomes only slightly smaller, but shifts to lower $\\sqrt{s}$, for $L \\approx 10$ fm. The paper concludes that the finite-volume dip-peak structure does not by itself indicate proximity of a critical endpoint; it can simply reflect the separation between the freeze-out line and the phase boundary.","pith_inferences":["Because the finite-size effect is purely through $\\beta V$, the model implies a universality the paper does not state: any two systems with the same $L^3/T$ but different $L$ and $T$ would give identical $R_{42}$, which could be checked by comparing fluctuation data across collision systems of different sizes.","A direct extension would be to compute Lee-Yang zeros of the partition function in the complex temperature or chemical-potential plane; the zero-mode weight suggests they move with $L$ in a calculable way, so finite-size data could be extrapolated to locate the infinite-volume critical endpoint more sharply than by peak positions.","The coexistence argument implies that in a genuine first-order transition at finite volume, the measured baryon-number susceptibility should grow roughly as $L^3$ while the condensate stays flat; comparing peripheral and central collisions across a first-order boundary would distinguish coexistence broadening from critical rounding.","The Binder-cumulant crossing technique, if applied to net-baryon cumulants rather than chiral ones, would give a finite-size robust estimator of the critical-endpoint location in experiment; the paper only demonstrates it for the chiral susceptibility."],"forward_implications":["At any finite $L$, the first-order transition is smoothed into a crossover; no true phase transition and no true divergence of $\\chi_2$ occurs in a finite fireball.","Near the critical endpoint, singular quantities obey finite-size scaling with an effective exponent $\\tilde{\\nu} = 2/3$; for $L \\gtrsim 20$ fm the data collapse onto universal curves, and for $L \\lesssim 10$--$20$ fm genuine finite-size corrections take over.","Along a first-order transition, the two-peaked structure of the probability distribution makes susceptibilities scale as $\\chi_k \\propto L^{kd}$; hence the maximum of the chiral or baryon susceptibility cannot be used to locate the critical endpoint at finite volume.","The Binder cumulant $\\kappa_B = \\chi_4/(\\beta L^d \\chi_2^2)$ is volume independent at both second- and first-order transitions and crosses at the infinite-volume critical endpoint, providing a size-robust locator.","Along the freeze-out line, $R_{42} = \\chi_4^B/\\chi_2^B$ is unchanged for $L \\gtrsim 40$ fm and only slightly reduced for $L \\approx 10$ fm, but its dip-peak feature shifts to lower $\\sqrt{s}$; because susceptibilities grow monotonically along the phase boundary, such a non-monotonic structure does not necessarily signal proximity to the critical endpoint."],"supporting_citations":[{"why":"Defines the mean-field model with momentum-space constraints that the present work replaces, and supplies the fermionic-matter potential and the observation that constraints alone do not give finite-size scaling.","marker":"[37]"},{"why":"Provides the finite-size scaling ansatz and the double-Gaussian treatment of coexistence used to derive the $L^d$ susceptibility growth at first-order transitions.","marker":"[7]"},{"why":"Gives the double-Gaussian probability density and the coexistence contribution to the susceptibility used in the subtraction and sector projection.","marker":"[8]"},{"why":"Introduces the effective exponent $\\tilde{\\nu} = (\\gamma+2\\beta)/d = 2/3$ used to collapse the critical-point scaling data in the mean-field approximation.","marker":"[44]"},{"why":"Lattice results for the size dependence of the phase boundary at $\\mu = 0$ that the model's downward shift is compared with.","marker":"[19]"},{"why":"Follow-up lattice extrapolation from imaginary chemical potentials used as the comparison for the finite-size phase boundary.","marker":"[20]"},{"why":"Freeze-out parametrization $T(\\sqrt{s})$, $\\mu(\\sqrt{s})$ that converts the model results into $\\sqrt{s}$-dependent $R_{42}$ predictions.","marker":"[53]"}],"fun_headline_variants":["Zero-mode weight melts QCD transition peaks","Finite size from a single integral, not momentum grids","Volume rounding: one zero-mode integral drives scaling","Zero-mode integral captures finite-size QCD scaling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"All finite-size effects are assumed to enter only through the factor $\\beta V$ multiplying a size-independent potential: the non-zero Fourier modes of the field, the momentum discretization of the fermion fluctuations, the gradient term, and the fermionic vacuum contribution are all dropped, so the numbers below $L \\sim 10$--$20$ fm are not under quantitative control.","fun_headline_variants_meta":{"raw":{"variants":["Zero-mode weight melts QCD transition peaks","Finite size from a single integral, not momentum grids","Volume rounding: one zero-mode integral drives scaling","Zero-mode integral captures finite-size QCD scaling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000202,"raw_usage":{"total_tokens":1394,"prompt_tokens":971,"completion_tokens":423,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":587,"completion_tokens_details":{"reasoning_tokens":363}},"tokens_in":587,"tokens_out":423,"duration_ms":4086,"temperature":1.0,"reasoning_tokens":363,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:41:26.468856+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the central observables with the non-zero modes restored (finite-box fermion-mode sums or a low-momentum cutoff in the matter potential, plus the scalar gradient term) at $L \\approx 10$ fm and $T$ near the would-be critical endpoint; the paper's own premise is that these are negligible, so if the chiral susceptibility or the $R_{42}$ dip-peak position shifts by more than the quoted few-percent level, the central claim fails. Alternatively, a lattice QCD calculation of the chiral susceptibility at finite volume near the critical point should show the predicted rounding with $\\chi_2 \\propto L^{3/2}$; a persistent divergence or a different exponent would falsify the framework.","supporting_citations":[{"cited_title":"Binder, M","cited_arxiv_id":null,"evidence_quote":"Provides the finite-size scaling ansatz and the double-Gaussian treatment of coexistence used to derive the $L^d$ susceptibility growth at first-order transitions."},{"cited_title":"Binder, Z","cited_arxiv_id":null,"evidence_quote":"Introduces the effective exponent $\\tilde{\\nu} = (\\gamma+2\\beta)/d = 2/3$ used to collapse the critical-point scaling data in the mean-field approximation."},{"cited_title":"Finite volume effects near the chiral crossover","cited_arxiv_id":"2401.01169","evidence_quote":"Lattice results for the size dependence of the phase boundary at $\\mu = 0$ that the model's downward shift is compared with."}],"review_version":2}