{"id":"7b2d3d58-2808-4f2e-8c91-a55d92456e4f","arxiv_id":"2501.03509","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"New lattice QCD results for second- and fourth-order quark number susceptibilities and charge fluctuation ratios using Möbius domain wall fermions show a possible discrepancy with staggered fermions in the electric charge susceptibility below 160 MeV.","lead":"Lattice QCD simulations with Möbius domain wall fermions produce new results for quark number susceptibilities and conserved charge fluctuations at a physical pion mass. The paper compares these first-principles results with hadron resonance gas models and perturbation theory, and finds a discrepancy in the electric charge susceptibility with staggered fermion results at low temperatures.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Single-spacing result: MDWF chi_Q^2 closer to HRG than staggered is not yet robust, since only Ntau=12 is used and the discrepancy could be an MDWF discretization artifact.","rationale":"The paper is a clearly labeled preliminary proceedings report, and the observables and lattice setup are described following standard methods. The central claim, however, is that the MDWF value of chi_Q^2 below Tpc is closer to HRG than staggered results, implying that staggered fermion results may be affected by discretization effects. The most load-bearing assumption behind this claim is that the Ntau=12 MDWF result is already close to the continuum. The paper provides no continuum extrapolation and only one lattice spacing, and it explicitly identifies the need for additional spacings as future work. Agreement with HRG models does not replace a continuum limit, because HRG is a model and not a first-principles target; moreover, the staggered comparison itself is at Ntau=16, so the observed difference could reflect MDWF coarseness rather than staggered error. This concern is the same one identified by the reader, and it justifies a CONDITIONAL verdict rather than a rejection: the paper is honest and methodologically sound, but the headline physics interpretation is not yet supported. No code or data artifacts are provided, which is typical for proceedings and limits reproducibility but is not a load-bearing concern. The residual-mass correction described in Section 3 is reasonable but secondary to the lattice-spacing issue. The concrete test of adding Ntau=16 (and ideally Ntau=8) would settle whether the discrepancy persists as an O(a^2) effect or reflects a real staggered-systematic difference.","tokens_in":8727,"tokens_out":2970,"duration_ms":28183,"concrete_test":"Calculate chi_Q^2 on a new ensemble with Ntau=16 (and ideally also Ntau=8) along the same LCP with ml/ms=1/27.4 and physical strange quark mass, matching the physical volume and pion mass. Compare the values at T=149.7 and 154.6 MeV to the Ntau=12 MDWF points and the staggered Ntau=16 data in Figure 2. If the finer MDWF point moves toward the staggered result by more than the combined uncertainty, the discrepancy is an O(a^2) discretization effect; if it remains at the HRG value, the claim of a staggered discretization effect is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central observation is the discrepancy between MDWF and staggered chi_Q^2 below about 160 MeV, with MDWF closer to HRG (Section 4.1, Figure 2), and the suggestion that MDWF may expose staggered discretization effects. This interpretation requires the MDWF point to be near the continuum limit. However, all MDWF results are at a single lattice spacing, Ntau=12 (Section 3), while the staggered comparison data are at Ntau=16, i.e. a finer spacing. At Ntau=12, O(a^2) errors can be substantial, and MDWF also has residual chiral symmetry breaking controlled only by Ls=12. Without a second MDWF spacing (Ntau=16 or Ntau=8), the observed gap could equally be an MDWF lattice artifact; agreement with HRG is a model comparison, not a continuum check. The authors themselves defer this check to future work, stating in Section 4.1 and Section 5 that additional lattice spacings are needed, so the load-bearing assumption is explicitly untested.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This LATTICE2024 proceedings paper from the JLQCD collaboration presents second- and fourth-order quark number susceptibilities and the derived second-order conserved charge (B, Q, S) fluctuations for (2+1)-flavor QCD with Möbius domain wall fermions (MDWF) at a pion mass of 135 MeV. The simulations use a 36^3×12 lattice with L_s=12, M_5=1, three-level stout smearing, and m_l/m_s=1/27.4 along a line of constant physics. The results are compared with the PDGHRG and QMHRG2020 hadron resonance gas models below T_pc, with O(g^2) perturbation theory at high temperature, and with staggered fermion results at Nτ=16 (HISQ and stout). The paper's central new observation is that χ_Q^2 from MDWF lies above the staggered values and closer to the HRG predictions below about 160 MeV (Fig. 2), along with preliminary leading-order kurtosis ratios R^Q_42 and R^S_42 (Eq. 20). The paper explicitly labels the results preliminary and defers additional lattice spacings and continuum extrapolation to future work.","tokens_in":8943,"tokens_out":10286,"duration_ms":99490,"significance":"If the χ_Q^2 discrepancy between MDWF and staggered fermions survives a continuum check, the paper would provide a genuinely valuable independent, chirally symmetric determination of these fluctuation observables and a concrete hint of discretization effects in staggered calculations. The manuscript has real strengths: the cumulant formalism in Eqs. (1)–(18) is standard, correctly stated, and correctly reduced to the K-term expectation values; the statistical errors are shown; the statistics are substantial at about 20,000 trajectories per temperature; and the authors are honest that these are preliminary results. The main limitation is that the load-bearing comparison rests on a single MDWF lattice spacing, so the strength of the central claim is currently capped; the requested changes are targeted at making that limitation explicit and quantitative rather than at any error in the derivation.","major_comments":[{"comment":"The central observation that χ_Q^2 from MDWF is closer to the HRG predictions than the staggered results below T≈160 MeV rests on a single MDWF lattice spacing (Nτ=12, a≈0.11 fm at T≈150 MeV) compared with staggered data at the finer Nτ=16. Without a second MDWF spacing, the gap is as easily explained by an O(a^2) artifact in the MDWF data as by a discretization effect in the staggered data, and agreement with an HRG model is not a substitute for a continuum check. The paper's own statement that additional spacings will be studied in the future should be moved to the point of the claim and strengthened: the claim should be explicitly qualified as a single-spacing observation, or the comparison should include the continuum-extrapolated staggered results from Ref. [20], which would test whether the Nτ=12 MDWF point is actually closer to the continuum than the Nτ=16 staggered points.","section":"§4.1, Fig. 2"},{"comment":"The statement that \"HRG model calculations of R^Q_42 overshoot the lattice data\" is not quantitatively supported as written: the quoted values are R^Q_42 = 1.05±0.46 at T=149.7 MeV and 1.00±0.53 at T=154.6 MeV, while the HRG predictions are not given in the text. The paper should quote the PDGHRG and QMHRG2020 values used and state the separation in units of the combined uncertainty; with errors of this size, the difference is at most 1–2σ and the word \"overshoot\" overstates the case. The analogous claim that R^S_42 is \"consistent with both HRG models\" should also be documented with the model numbers rather than left implicit in the figure.","section":"§4.3, Eq. (20), Fig. 4"},{"comment":"Since χ_B^2 and χ_S^2 in Fig. 2 show no obvious MDWF–staggered discrepancy, the χ_Q^2 difference must originate in a specific combination of the diagonal and off-diagonal light-quark terms, for example χ_u^2 and χ_ud^11. The paper does not identify which contribution drives the effect. A table or panel showing the individual MDWF components alongside the staggered ones would let the reader judge whether the discrepancy is a light-quark or disconnected-sector effect, which is directly relevant to the paper's main claim and would make the preliminary result substantially more informative.","section":"§4.1, Eqs. (17)–(18)"}],"minor_comments":[{"comment":"The operator is called the \"MDMF Dirac operator\" but should read \"MDWF\", and the display equation for the link substitution contains a stray period after the arrow.","section":"§2, text after Eq. (3)"},{"comment":"The phrase \"weakly interacting gas of hadrons and gluons\" is imprecise: the ideal-gas limit discussed here is a gas of quarks and gluons, not hadrons.","section":"§4.1, first paragraph"},{"comment":"The residual mass m_res is used to correct the bare quark masses and to define the line of constant physics, but its value in lattice units is not reported; since it is taken from a different mass-ratio ensemble, the value and its assumed β-dependence should be stated so the reader can assess the LCP definition.","section":"§3"},{"comment":"The paper reports statistical errors but does not discuss any systematic error budget; at a minimum, the size of the scale-setting uncertainty, the m_res uncertainty, and finite-volume effects should be estimated or explicitly argued to be negligible relative to the statistical errors.","section":"§3 and §4"},{"comment":"The Fig. 3 caption has a stray \"and.\" after \"free quark gas\", and the Fig. 2 caption should say \"at finite lattice spacing\" rather than \"at finite lattice\".","section":"Fig. 3 caption and Fig. 2 caption"},{"comment":"The noise discussion in §4.2 would be more useful if the number of stochastic sources and the dilution scheme used for the D_n^f estimates were reported, and Eq. (19) should define the shorthand X=Q,S before it is used in the text.","section":"§4.2 and Eq. (19)"}],"recommendation":"major_revision","confidential_remarks":"This is a LATTICE2024 proceedings contribution, so the preliminary status of the results is appropriate for the venue; nevertheless, the single-spacing issue is exactly the point that should gate the interpretation, and the authors' own hedging does not fully neutralize the rhetorical framing of the χ_Q^2 result in §4.1 and §5. The requested changes are modest, local, and within the manuscript's scope: re-qualify the central claim at the point where it is made, compare with the continuum-extrapolated staggered results where available, and quantify the HRG comparisons for the kurtosis ratios. I see no issues with the citation pattern; the self-citations are to the collaboration's own scale-setting and residual-mass calibrations, which is standard practice."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nShort version: this is a solid, clearly labeled preliminary proceedings paper from JLQCD. The new thing is the first set of second- and fourth-order quark number susceptibilities with Möbius domain wall fermions at a physical pion mass, including leading kurtosis ratios. The formalism is standard and correct, the statistical errors are shown, and the comparisons to HRG and O(g^2) perturbation theory are sensible. The authors are honest about noise: the fourth-order light-quark quantities have large errors, and R_Q42 at T≈150-155 MeV carries ±0.5 uncertainties, so the kurtosis results are not yet discriminative.\n\nThe interesting physics claim is in Figure 2: below about 160 MeV the MDWF chi_Q^2 lies between the staggered results and the HRG curves, closer to HRG. The authors float the idea that MDWF might expose staggered discretization effects. That is a legitimate hypothesis, but at this stage it is an observation, not a demonstrated discrepancy. All MDWF points are on a single lattice spacing, Ntau=12, while the staggered comparison data are at Ntau=16. At Ntau=12, O(a^2) effects plus residual chiral symmetry breaking with Ls=12 could plausibly shift chi_Q^2 by the amount in question. The paper says in Section 5 that additional lattice spacings are needed; I agree, and I would not cite the discrepancy as established until Ntau=16 or 8 results exist.\n\nMinor points: the residual mass is taken from the m_l/m_s=1/10 series and assumed mass-independent; plausible, but a test or a reference would have been nice. No code or data are provided, which is typical for proceedings and limits reproducibility but is not a fatal flaw here.\n\nWho should read this: people working on lattice QCD thermodynamics, especially on charge fluctuations, will want to know these numbers exist and that JLQCD is heading toward a continuum MDWF equation of state. It is not a paper that settles anything by itself.\n\nRecommendation: as a proceedings contribution, it deserves a serious referee and is acceptable in that context. As a journal paper it would need at least one more lattice spacing and a continuum extrapolation before the chi_Q^2 claim can be evaluated. I would accept it for review rather than desk-reject.","headline":"Solid preliminary MDWF susceptibility results; the chi_Q^2 vs staggered gap is a single-spacing observation, not a demonstrated discrepancy.","tokens_in":9513,"tokens_out":3353,"would_cite":true,"duration_ms":33343,"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":"Möbius domain wall fermion results put low-temperature electric charge fluctuations closer to hadron resonance gas predictions than staggered fermion calculations do.","keywords":["quark number susceptibility","conserved charge fluctuation","Möbius domain wall fermions","lattice QCD","hadron resonance gas","chemical potential","kurtosis"],"falsifier":"Compute $\\chi^Q_2$ with the same Möbius domain wall action on a finer lattice, for example $N_\\tau = 16$, at the same physical pion mass and temperatures between roughly 140 and 160 MeV. If the values move toward the staggered results as the spacing decreases, the low-temperature discrepancy is a discretization effect; if they stay near the hadron resonance gas curves, the staggered calculations are the ones being displaced.","tokens_in":8531,"feed_emoji":"⚛️","tokens_out":9609,"duration_ms":81810,"temperature":0.7,"pith_summary":"This paper works out second- and fourth-order quark number susceptibilities and the corresponding fluctuations of baryon number, electric charge, and strangeness for (2+1)-flavor QCD, using Möbius domain wall fermions with a 135 MeV pion mass. The central result is that below the QCD crossover temperature, the electric charge susceptibility $\\chi^Q_2$ computed with this chirally symmetric fermion formulation is larger than published staggered-fermion results and closer to hadron resonance gas model predictions. If the calculation is right, the discrepancy is a signal that staggered discretizations may distort this observable at finite lattice spacing, and the domain-wall data become a benchmark for a continuum-extrapolated determination. At high temperature the same susceptibilities approach the $\\mathcal{O}(g^2)$ perturbative band, and the leading-order kurtosis ratios $\\chi^Q_4/\\chi^Q_2$ and $\\chi^S_4/\\chi^S_2$ are reported for comparison with heavy-ion freeze-out analyses.","feed_headline":"Chiral QCD lattice moves charge fluctuation toward hadron gas","feed_subtitle":"At a 135 MeV pion mass, the chirally symmetric calculation lands closer to hadron resonance gas than staggered data below 160 MeV.","key_machinery":"The organizing object is the set of generalized quark number susceptibilities $\\chi^{uds}_{ijk}$, the Taylor coefficients of the QCD pressure in the quark chemical potentials $\\hat\\mu_u,\\hat\\mu_d,\\hat\\mu_s$. These are computed as expectation values of derivatives $D^f_n = \\partial^n \\ln \\det M_f / \\partial \\hat\\mu_f^n$, where the Möbius domain wall fermion determinant $M_f$ carries the chemical potential through modified temporal gauge links $U_{\\pm 4} \\to e^{\\pm \\hat\\mu_f} U_{\\pm 4}$. Möbius domain wall fermions are a lattice fermion formulation that keeps chiral symmetry to good accuracy at finite lattice spacing. The same coefficients are linearly transformed to conserved charge susceptibilities $\\chi_2^B,\\chi_2^Q,\\chi_2^S$ and the fourth-order combinations. A useful structural fact exploited here is that, with degenerate $u$ and $d$ quarks, $(D^u_1 D^d_1) = (D^u_1)^2$ in the expectation values, which cancels part of the noise and makes $\\chi^Q_2$ the cleanest of the second-order observables.","core_discovery":"The paper claims that a (2+1)-flavor lattice QCD calculation with Möbius domain wall fermions, at a single lattice spacing set by $N_\\tau = 12$ and a physical light quark mass ($m_l/m_s = 1/27.4$, pion mass 135 MeV), produces second-order conserved charge fluctuations that behave differently from staggered-fermion results in the hadronic phase. In particular, for temperatures below about 160 MeV the electric charge susceptibility $\\chi^Q_2$ lies above the staggered values and in better agreement with both the 3/4-star resonance gas model and the QMHRG2020 hadron resonance gas model, while the baryon, strangeness, and mixed susceptibilities show no such tension. Fourth-order susceptibilities are noisier but consistent with the free-quark gas and $\\mathcal{O}(g^2)$ perturbation theory above $T_{pc}$; the resulting leading-order kurtosis ratios are $R^Q_{42} = 1.05 \\pm 0.46$ at $T = 149.7$ MeV and $R^S_{42} = 1.38 \\pm 0.09$ at the same temperature, quantities relevant for locating the QCD critical point in heavy-ion data.","pith_inferences":["Beyond the paper: if the $\\chi^Q_2$ enhancement survives a continuum extrapolation, freeze-out parameters extracted from staggered $\\chi^Q_2$ in heavy-ion analyses would need to be re-examined, because the electric charge susceptibility is the observable most directly tied to the net-charge cumulants measured in experiment.","Beyond the paper: the fact that only $\\chi^Q_2$ shows the discrepancy points toward taste-breaking or rooting artifacts that couple to electric charge rather than to baryon number; testing the ratio $\\chi^Q_2/\\chi^B_2$ on the same ensembles would sharpen this.","Beyond the paper: a natural next step the authors do not report is to compute the same observables with two lattice spacings and perform a continuum extrapolation using the $N_\\tau = 12$ data together with a planned $N_\\tau = 16$ ensemble, which would directly test the weakest assumption.","Beyond the paper: the noise pattern suggests that $\\chi^Q_2$ is the observable where domain-wall fermions have the best chance of beating staggered results; future high-statistics comparisons should focus there rather than on $\\chi^B_2$, where the stochastic error is larger."],"forward_implications":["A continuum-extrapolated MDWF calculation would decide whether the staggered-fermion $\\chi^Q_2$ below $T_{pc}$ is shifted by discretization effects.","The $\\chi^Q_2$ agreement with hadron resonance gas models strengthens the hadronic degrees-of-freedom interpretation of the transition region, without requiring extra states beyond the resonance list for the non-strange sector.","The leading-order kurtosis ratios $R^Q_{42}$ and $R^S_{42}$ give baseline lattice predictions that can be confronted with heavy-ion freeze-out analyses once experimental errors shrink.","The high-temperature approach to the $\\mathcal{O}(g^2)$ band confirms that degrees of freedom in the quark-gluon plasma are weakly interacting at $T \\gtrsim 180$ MeV, consistent with previous staggered results for diagonal susceptibilities."],"supporting_citations":[{"why":"Supplies the Möbius domain wall fermion algorithm used to define the lattice Dirac operator and the action for all ensembles.","marker":"[3]"},{"why":"Fixes the line of constant physics and the mass renormalization that set the pion mass to 135 MeV.","marker":"[5]"},{"why":"Motivates the leading-order kurtosis ratios as observables sensitive to the QCD critical point.","marker":"[7]"},{"why":"Documents the earlier work on noise in these susceptibilities and the u/d degeneracy cancellation that makes $\\chi^Q_2$ clean.","marker":"[9]"},{"why":"Provides the $\\mathcal{O}(g^2)$ perturbative QCD pressure and susceptibility result used as the high-temperature reference band.","marker":"[12]"},{"why":"Supplies the Taylor-expansion framework and running-coupling scale choice for the perturbative comparison.","marker":"[17]"},{"why":"Gives previous staggered-fermion results for conserved charge fluctuations used as the first comparison set.","marker":"[18]"},{"why":"Gives additional staggered results on the QCD transition, used to check consistency of diagonal susceptibilities.","marker":"[19]"},{"why":"Provides the QMHRG2020 hadron resonance gas model and the HISQ staggered data shown in the comparison plots.","marker":"[20]"},{"why":"Provides the stout-staggered data used as the second staggered comparison in the conserved-charge plots.","marker":"[21]"}],"fun_headline_variants":["Möbius domain wall QCD brings charge fluctuation to hadron gas","Charge fluctuation from Möbius fermions favors hadron gas","Möbius fermion lattice QCD: charge fluctuation matches HRG","At 135 MeV pion, Möbius fermions push charge susceptibility to HRG"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole comparison with hadron resonance gas models presumes that the single coarse lattice spacing used here ($N_\\tau = 12$) does not introduce a significant discretization error in $\\chi^Q_2$; without a second lattice spacing or a continuum extrapolation, the low-temperature difference from staggered results could just as easily be a lattice artifact.","fun_headline_variants_meta":{"raw":{"variants":["Möbius domain wall QCD brings charge fluctuation to hadron gas","Charge fluctuation from Möbius fermions favors hadron gas","Möbius fermion lattice QCD: charge fluctuation matches HRG","At 135 MeV pion, Möbius fermions push charge susceptibility to HRG"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000948,"raw_usage":{"total_tokens":4027,"prompt_tokens":909,"completion_tokens":3118,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":525,"completion_tokens_details":{"reasoning_tokens":3034}},"tokens_in":525,"tokens_out":3118,"duration_ms":21009,"temperature":1.0,"reasoning_tokens":3034,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:52:14.152062+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $\\chi^Q_2$ with the same Möbius domain wall action on a finer lattice, for example $N_\\tau = 16$, at the same physical pion mass and temperatures between roughly 140 and 160 MeV. If the values move toward the staggered results as the spacing decreases, the low-temperature discrepancy is a discretization effect; if they stay near the hadron resonance gas curves, the staggered calculations are the ones being displaced.","supporting_citations":[{"cited_title":"8 (2+1)-Flavor QCD with MDW Fermions Jishnu Goswami","cited_arxiv_id":null,"evidence_quote":"Documents the earlier work on noise in these susceptibilities and the u/d degeneracy cancellation that makes $\\chi^Q_2$ clean."},{"cited_title":"The pressure of QCD at finite temperatures and chemical potentials","cited_arxiv_id":"hep-ph/0305183","evidence_quote":"Provides the $\\mathcal{O}(g^2)$ perturbative QCD pressure and susceptibility result used as the high-temperature reference band."}],"review_version":1}