{"id":"60787413-14e7-482c-bfe7-0aa2f5b13750","arxiv_id":"2507.21744","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A finite-volume PNJL model calculation of net-charge moments at RHIC BES energies shows non-monotonic behavior for R=2 fm, which the authors connect to the QCD critical point, though the numerical setup is incomplete.","lead":"This paper computes net-charge fluctuation moments in a finite-volume PNJL quark model at RHIC Beam Energy Scan energies and compares them with STAR data. It reports non-monotonic signals for small system sizes that it reads as evidence for the QCD critical point, but key parts of the calculation are left unspecified.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The finite-volume replacement of discrete momentum modes by a continuous integral with cutoff λ=π/R is uncontrolled; the 2 fm non-monotonic signal could be an artifact of that cutoff, so the critical-point interpretation is not established.","rationale":"The paper attempts to connect finite-volume PNJL calculations to STAR net-charge data and to infer critical-point sensitivity. For that connection to hold, the finite-volume model must be a reliable implementation of a small system, and the energy axis must correspond to the actual freeze-out chemical potentials. The weakest point is the former. Eq. (18) is obtained from Eq. (16) by replacing the momentum sum with an integral from λ=π/R to ∞; the paper explicitly says 'the infinite sum has been regarded as an integration over continuous variation of momentum even though with a lower momentum cutoff and the surface and curvature effects have been disregarded.' That statement is an admission that the finite-volume physics is not modeled: a cubic box of side R has discrete modes (p_i=2π n_i/R for bosons, p_i=(2n_i+1)π/R for fermions), and the mode density near p≈π/R is not captured by a sharp lower cutoff. Since R=2 fm corresponds to λ≈98 MeV, the cutoff is not a small perturbation; it changes the vacuum and medium integrals substantially. The claimed non-monotonicity in the 2 fm moment products could therefore be a direct consequence of the chosen cutoff, not of critical fluctuations. This concern is load-bearing because the entire novelty of the paper is the finite-volume effect. The reader's verdict already identified this as the weakest assumption, and I agree. I would not change the REJECT verdict: the central claim is not falsified, but it is not established by the presented calculation. The paper has genuine supporting structure (an established PNJL parameter set, comparison with STAR/HRG/UrQMD/lattice data, self-consistent gap equations), but those do not rescue the uncontrolled finite-volume step. The missing µ_Q mapping is an additional reproducibility failure, but the finite-volume cutoff is the conceptual core. A targeted recomputation with exact discrete mode sums would settle whether the 2 fm non-monotonic signal survives.","tokens_in":13209,"tokens_out":4675,"duration_ms":54098,"concrete_test":"Recompute the CEP and C1-C4 for R=2 and 4 fm using exact discrete momentum sums with antiperiodic boundary conditions for quarks (p_i=(2n_i+1)π/R) instead of continuous integrals with lower cutoff λ=π/R, keeping all other parameters identical. If the 2 fm non-monotonic features in M/σ², Sσ, and κσ² disappear or move substantially, they are artifacts of the cutoff approximation. As a secondary check, obtain the explicit µ_Q(√sNN) mapping used to produce each energy point and rerun the calculation with that mapping; if no mapping exists or it changes the curves, the energy dependence is not reproducible.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the finite-volume three-flavor PNJL model produces non-monotonic net-charge moment products for R=2 fm and that this signals critical fluctuations. The load-bearing premise is the finite-volume prescription in Sec. II and Eq. (18): the discrete momentum sum is replaced by continuous integration with a lower cutoff λ=π/R, and surface/curvature effects are explicitly disregarded. This is not a controlled approximation. For R=2 fm, λ≈98 MeV, a scale comparable to the quark masses and to the model's vacuum cutoff; a single lower cutoff does not reproduce the density of discrete modes or the difference between periodic and antiperiodic boundary conditions. The shift of the CEP in earlier finite-volume PNJL studies depends on this same ad hoc cutoff, so the energy dependence of C1-C4 reported here may be an artifact of λ rather than finite-size critical physics. The inconsistency of the paper's own comparison statements compounds the problem: Sec. III first says the 2 fm results 'fail to capture' STAR/model behavior at lower energies, then claims 'reasonable agreement' at 7-20 GeV, while Sec. IV concludes 'better agreement with STAR data at lower energies.' A secondary but severe omission is that the quark chemical potential used for net-charge is never specified ('fixed quark chemical potential', no µ_Q(√sNN) relation), making the energy axis non-reproducible. Together these undermine the central interpretation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes the first four cumulants C1-C4 of net-charge fluctuations and the volume-independent moment products M/sigma^2, S*sigma, and kappa*sigma^2 in a three-flavor PNJL model with finite volume implemented through a lower momentum cutoff lambda = pi/R. Results are presented for R = 2 fm and R = 4 fm as functions of sqrt(s_NN) from 2.4 to 200 GeV and are compared with STAR net-charge data, lattice QCD, HRG, UrQMD, and HIJING. The central claim is that the 2 fm system shows pronounced non-monotonic behavior and better agreement with STAR at lower energies, which the authors interpret as sensitivity to critical fluctuations in smaller systems.","tokens_in":13520,"tokens_out":4231,"duration_ms":52252,"significance":"If established, the claim would provide a finite-volume, model-based probe of the QCD critical point and a possible tool for extracting freeze-out parameters from net-charge moments. The paper is not circular: the PNJL couplings and Polyakov potential are taken from earlier fits to vacuum meson properties and to lattice QCD thermodynamics, and the freeze-out parametrization in Eq. (19) is fitted to chemical freeze-out data, so the comparison with STAR net-charge moments is not contaminated by using those data as input. However, the significance is currently limited by the ad hoc finite-volume prescription, the absence of an explicit chemical-potential mapping, and the lack of a transparent derivation of C1-C4 from the thermodynamic potential. These omissions make the central non-monotonic signal difficult to reproduce or to distinguish from a cutoff artifact.","major_comments":[{"comment":"The finite-volume prescription is uncontrolled and is load-bearing for the main claim. Replacing the discrete momentum sum over finite-volume modes by a continuum integral with lower cutoff lambda = pi/R disregards surface and curvature effects and the difference between periodic and antiperiodic boundary conditions. For R = 2 fm, lambda is about 98 MeV, which is comparable to the dynamically generated quark masses and not negligible relative to the three-momentum cutoff Lambda = 631 MeV. The authors explicitly state that surface and curvature effects are disregarded, but they do not provide any validation that this cutoff prescription reproduces the true finite-volume mode density. The 2 fm non-monotonic signal in Figures 2 and 3 could therefore be an artifact of the cutoff rather than a finite-size critical effect. The authors should compare the cutoff-integral results with a discrete momentum sum over the appropriate boundary conditions, or at least demonstrate that the non-monotonic features persist under variations of the cutoff prescription.","section":"Sec. II, Eq. (18)"},{"comment":"The quark chemical potential used for net-charge is never specified. The text states that results are obtained at a fixed quark chemical potential, but no relation mu_Q(sqrt(s_NN)) or values of mu_u, mu_d, mu_s are given. The freeze-out parametrization in Eq. (19) fixes T(mu_B) and mu_B(sqrt(s_NN)), but it does not determine the individual quark chemical potentials or the net-charge chemical potential. Without this mapping, the energy dependence of C1-C4 in Figures 1-3 is not reproducible, and the comparison with STAR data cannot be interpreted. This is a central omission because the entire observable depends on the chemical-potential assignment.","section":"Sec. III"},{"comment":"There is no explicit expression for C1-C4 in terms of derivatives of the thermodynamic potential Omega' in Eq. (18). Equations (2)-(7) define cumulants and their relation to susceptibilities, but the manuscript does not state how the net-charge susceptibilities chi_Q^(n) are obtained from Omega'—for example, as derivatives with respect to mu_Q at fixed T and mu_B. Without these expressions, the central results in Figures 1-3 cannot be checked or reproduced. The authors should provide the explicit derivative relations and, if possible, the numerical procedure used to evaluate them.","section":"Secs. II and III"},{"comment":"The comparison statements are internally inconsistent and affect the central interpretation. After Figure 2, the text says that at lower energies 'the 2 fm data fails to capture the experimental and theoretical behavior observed in STAR and model studies.' Later in Section III, for Figure 3, it says the 2 fm results show 'reasonable agreement with both the experimental and lattice QCD data in the energy range of 7-20 GeV.' Section IV then concludes that the 2 fm system shows 'better agreement with STAR data at lower energies.' These statements cannot all be true as written, and the paper should state clearly whether the 2 fm model agrees or disagrees with STAR at low energies.","section":"Secs. III and IV"}],"minor_comments":[{"comment":"The caption describes 'infinite volume systems with lateral size R = 2fm and R = 4fm,' which contradicts the finite-volume implementation described in Sec. II and Eq. (18). This appears to be a typo, but it should be corrected because the distinction between finite and infinite volume is central to the paper.","section":"Fig. 1 caption"},{"comment":"The cross-reference structure is broken: the text refers to 'Section II summarizes the mathematical definitions,' 'Section III briefly outlines the formalism,' and 'the conclusion is given in the Section ??,' but the actual sections are organized differently. The section numbering and cross-references should be corrected.","section":"Introduction"},{"comment":"There are numerous typographical errors that should be fixed, including 'handron' for 'hadron' in the Introduction, 'anamoly' for 'anomaly' in Sec. II, and inconsistent notation such as 'rce' for the trace in the definition of the Polyakov loop. A careful proofread is needed.","section":"Throughout"},{"comment":"The conclusion states that 'both system size and interaction strength play crucial roles,' but the paper does not present eight-quark interaction results, and the interaction strength is not varied in the figures. The conclusion should be limited to what is actually shown.","section":"Sec. IV"},{"comment":"The logarithmic terms in Eq. (18) have unbalanced parentheses, which makes the expression difficult to parse. The authors should rewrite the equation with unambiguous bracket structure.","section":"Eq. (18)"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a topical question, but its central claim currently rests on an uncontrolled finite-volume approximation and on an unspecified chemical-potential mapping. The comparison with STAR data is not circular, which is a point in the paper's favor. If the authors can supply the missing mu_Q(sqrt(s_NN)) relation, explicit formulas for the cumulants, and a robustness check of the finite-volume cutoff prescription, the manuscript could become publishable. Without those additions, the non-monotonic 2 fm signal cannot be distinguished from a cutoff artifact."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe new thing here is a finite-volume three-flavor PNJL computation of net-charge cumulants C1–C4 and the moment products M/σ², Sσ, κσ² at BES energies plus 2.4 and 3 GeV. That specific calculation is not in the cited papers, so it is a genuine numerical extension. The authors take the model parameters and freeze-out parameterization from earlier work rather than fitting to the STAR moments they compare with, so the comparison is not circular. The paper's ambition is modest: it offers a model-to-data tool, not a CEP location.\n\nWhat I find soft:\n\n1. The finite-volume prescription is a single lower momentum cutoff λ=π/R replacing the discrete sum, with surface and curvature effects dropped. That is a crude approximation, and for R=2 fm, λ≈98 MeV, which is the same scale as the quark masses and not far below the UV cutoff. Non-monotonic features that show up only at 2 fm could easily be artifacts of where the cutoff sits rather than finite-size critical physics. The paper acknowledges the approximation but does not quantify its uncertainty.\n\n2. The paper never specifies the quark chemical potentials for net-charge. It says \"fixed quark chemical potential\" but gives no μ_Q(√sNN) mapping. Since the energy dependence is the entire observable, this makes the calculation non-reproducible. You could not repeat the computation from the text.\n\n3. There is an internal inconsistency in the agreement claims. Section III first says the 2 fm results \"fail to capture\" the lower-energy behavior, then says there is \"reasonable agreement\" at 7–20 GeV, and Section IV concludes \"better agreement with STAR data at lower energies.\" Those statements conflict.\n\n4. The cumulants are presented without the explicit derivatives of Ω' that generate them. A reader can infer they come from standard thermodynamic relations, but for a paper whose whole point is these cumulants, the formulas should be there.\n\nThe stress-test note on the uncontrolled cutoff is fair; it is the load-bearing weakness. The circularity concern, by contrast, is not a real problem here — the model is not fit to the STAR moments.\n\nRecommendation: This deserves a serious referee because it is a legitimate extension in a hot-topic area and the flaws are fixable in revision, but the central interpretation should be flagged as not established. If the authors can provide the missing chemical potential definitions, show the cumulant formulas, and add a robustness check on the cutoff (e.g., varying λ or comparing with a discrete-mode sum), the paper could be publishable after major revision.","headline":"A legitimate finite-volume PNJL extension for net-charge cumulants, but the uncontrolled momentum-cutoff prescription and missing chemical-potential definitions leave the critical-point interpretation unestablished.","tokens_in":14123,"tokens_out":2721,"would_cite":false,"duration_ms":29007,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Aw","12.38.Mh","12.39.-x"],"model":"deepseek-v4-flash","headline":"Small fireballs should flash critical-point signals in net-charge moments","keywords":["QCD critical point","net-charge fluctuations","finite volume PNJL model","higher-order cumulants","skewness","kurtosis","RHIC Beam Energy Scan","freeze-out parameters"],"falsifier":"A direct test is to compute the same net-charge cumulants in the finite-volume PNJL model with the full discrete momentum sum (periodic/antiperiodic boundary conditions) instead of the lower-cutoff continuum approximation; if the discrete-sum results do not reproduce the non-monotonic 2 fm behavior in Sσ and κσ2, the claimed critical-point signal is an artifact of the cutoff approximation.","tokens_in":12944,"feed_emoji":"⚛️","tokens_out":1512,"duration_ms":17502,"temperature":0.7,"pith_summary":"This paper claims that the finite-volume three-flavor PNJL model produces non-monotonic, volume-independent moment products of net-charge as a function of collision energy, and that this behavior is a signature of the QCD critical point. The authors compute cumulants C1 through C4 and the products M/σ2, Sσ, and κσ2 for systems of radius 2 fm and 4 fm across RHIC BES energies from 7.7 to 200 GeV, plus 2.4 and 3 GeV, and compare with STAR net-charge data, HRG, UrQMD, HIJING, and lattice QCD. They find that the 2 fm system shows pronounced non-monotonic fluctuations and better agreement with STAR at lower energies, suggesting that smaller fireballs are more sensitive to critical fluctuations near the critical end point (CEP).","feed_headline":"Small fireballs should flash critical-point signals in net-charge moments","feed_subtitle":"A finite-volume PNJL model predicts non-monotonic net-charge moment products that match STAR data at low energy.","key_machinery":"The central object is the finite-volume Polyakov loop extended Nambu-Jona-Lasinio (PNJL) model with six-quark interactions, where finite volume is implemented by imposing a lower momentum cutoff λ = π/R on both vacuum and medium integrals. The cumulants C1–C4 of net-charge are computed from the thermodynamic potential, and the volume-independent ratios M/σ2 = C1/C2, Sσ = C3/C2, and κσ2 = C4/C2 are built to cancel volume dependence. The non-monotonic energy dependence of these ratios is the signal claimed to be connected to the critical point, via the relation of cumulants to susceptibilities and powers of the correlation length.","core_discovery":"The central claim is that in the finite-volume PNJL model, the moment products of net-charge—especially Sσ and κσ2—develop non-monotonic, volume-independent structures as a function of beam energy, and that the 2 fm system exhibits larger fluctuations and better qualitative agreement with STAR net-charge data at lower energies than the 4 fm system. The authors interpret this as evidence that smaller systems are more sensitive to critical fluctuations near the QCD critical point, and that comparing these moment products with STAR results can serve as a tool for extracting freeze-out parameters.","pith_inferences":["A testable extension is to compute the same moment products for net-baryon and net-strangeness with the identical finite-volume cutoff; if the 2 fm non-monotonicity appears in all three conserved charges at the same energy, the critical-point interpretation is strengthened, whereas charge-only signals would suggest a different mechanism.","The claim implies that existing STAR low-energy data already contain the signature, meaning a re-analysis of net-charge moments with explicit centrality-dependent fireball radii—rather than a single 2 fm or 4 fm value—should reproduce the non-monotonic dip-and-rise structure if the model is right.","A stronger test would be to compute the same observables at 3D Ising universality-class critical exponents within the PNJL model; the current non-monotonicity is qualitative, and matching the predicted power-law growth of Sσ and κσ2 near the CEP would distinguish genuine critical behavior from model artifacts."],"forward_implications":["If the 2 fm non-monotonic moment products are correct, smaller fireballs from lower-energy heavy-ion collisions should carry enhanced critical-fluctuation signals detectable in STAR-like net-charge measurements.","Comparing PNJL moment products with measured STAR values can provide a parameter-free route to extract freeze-out temperature and baryon chemical potential at each beam energy.","The volume independence of the moment products means the same computed ratios can be compared across different centralities and system sizes, removing the main ambiguity of unknown interaction volume.","The shift of the CEP toward higher µB and lower T in smaller volumes, implied by the finite-volume model, predicts where in the T-µB plane to scan for non-monotonic signals.","The discrepancy between 2 fm and 4 fm results at low energies indicates that system-size dependence itself is a diagnostic: matching experimental data may require modeling the actual fireball size rather than assuming a single large volume."],"supporting_citations":[{"why":"Establishes that the CEP shifts toward higher µB and lower T with decreasing volume in the PNJL model, the premise for why a 2 fm system should show stronger critical sensitivity.","marker":"[35]"},{"why":"Provides the re-parametrized three-flavor PNJL parametrization (T0, a0, a1, a2, b3, b4, κ) used for the six-quark interaction in the present calculation.","marker":"[40]"},{"why":"Supplies the freeze-out parameterization T(µB) and µB(√sNN) needed to convert beam energy into temperature and chemical potential for the model.","marker":"[13]"},{"why":"The STAR net-charge experimental data against which both Sσ and κσ2 results are compared to claim agreement at low energies.","marker":"[45]"},{"why":"The UrQMD model predictions used to contrast the PNJL results and to support the claim that the 2 fm behavior differs from non-critical transport models.","marker":"[48]"},{"why":"Defines the Landau-Ginzburg Polyakov loop potential with the Z(3) symmetric form used in the PNJL thermodynamic potential.","marker":"[50]"},{"why":"Supplies lattice QCD data for Sσ3/M ratio used to benchmark the PNJL results in the 7–20 GeV range.","marker":"[56]"}],"fun_headline_variants":["Finite volume sharpens net-charge signals for QCD critical point","Small systems amplify net-charge moments toward critical point","Net-charge moment products expose QCD critical point fingerprints","Non-monotonic net-charge moments match STAR data in finite volume"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that replacing the true finite-volume momentum spectrum with a continuous integral cut off at λ = π/R captures the physics of a 2 fm or 4 fm fireball, while surface and curvature effects are negligible.","fun_headline_variants_meta":{"raw":{"variants":["Finite volume sharpens net-charge signals for QCD critical point","Small systems amplify net-charge moments toward critical point","Net-charge moment products expose QCD critical point fingerprints","Non-monotonic net-charge moments match STAR data in finite volume"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000405,"raw_usage":{"total_tokens":2157,"prompt_tokens":1045,"completion_tokens":1112,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":661,"completion_tokens_details":{"reasoning_tokens":1040}},"tokens_in":661,"tokens_out":1112,"duration_ms":8800,"temperature":1.0,"reasoning_tokens":1040,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T12:26:03.593839+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test is to compute the same net-charge cumulants in the finite-volume PNJL model with the full discrete momentum sum (periodic/antiperiodic boundary conditions) instead of the lower-cutoff continuum approximation; if the discrete-sum results do not reproduce the non-monotonic 2 fm behavior in Sσ and κσ2, the claimed critical-point signal is an artifact of the cutoff approximation.","supporting_citations":[{"cited_title":"Bhattacharyya, R","cited_arxiv_id":null,"evidence_quote":"Establishes that the CEP shifts toward higher µB and lower T with decreasing volume in the PNJL model, the premise for why a 2 fm system should show stronger critical sensitivity."},{"cited_title":"Bhattacharyya, S","cited_arxiv_id":null,"evidence_quote":"Provides the re-parametrized three-flavor PNJL parametrization (T0, a0, a1, a2, b3, b4, κ) used for the six-quark interaction in the present calculation."},{"cited_title":"Cleymans, H","cited_arxiv_id":null,"evidence_quote":"Supplies the freeze-out parameterization T(µB) and µB(√sNN) needed to convert beam energy into temperature and chemical potential for the model."},{"cited_title":"Adamczyk et al","cited_arxiv_id":null,"evidence_quote":"The STAR net-charge experimental data against which both Sσ and κσ2 results are compared to claim agreement at low energies."},{"cited_title":"Cumulants of Net-Proton, Net-Kaon and Net-Charge Multiplicity Distributions in Au+Au Collisions at RHIC BES Energies from UrQMD Model","cited_arxiv_id":"1606.03900","evidence_quote":"The UrQMD model predictions used to contrast the PNJL results and to support the claim that the 2 fm behavior differs from non-critical transport models."},{"cited_title":"Ratti, M","cited_arxiv_id":null,"evidence_quote":"Defines the Landau-Ginzburg Polyakov loop potential with the Z(3) symmetric form used in the PNJL thermodynamic potential."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies lattice QCD data for Sσ3/M ratio used to benchmark the PNJL results in the 7–20 GeV range."}],"review_version":1}