{"id":"2656c447-51bc-4147-8690-f8d2d3b9b0ac","arxiv_id":"2508.19337","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"New finite-volume lattice data for the subtracted chiral condensate agree with 3-d O(2) scaling for H<=1/160, while larger H violates the universal curve.","lead":"Lattice QCD simulations of the chiral order parameter show that at small light quark masses, its finite-volume behavior matches the expected O(2) universal scaling, with deviations growing at larger masses. The data give a preliminary chiral transition temperature near 145 MeV on coarse lattices.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Universal O(2)-scaling claim rests on an unjustified 1/L truncation in Eq. (5): terms dropped from the z_L expansion in Eq. (4) could bias the infinite-volume extrapolation at the quoted precision.","rationale":"The reader returned CONDITIONAL because the volume extrapolation and non-universal fits are not fully documented. My stress-test is a sharpened version of that concern. The specific issue is Eq. (5) vs Eq. (4): Eq. (4) is explicitly a Taylor series in z_L, and z_L ∝ 1/L. Therefore the extrapolation ansatz should not arbitrarily start at 1/L^3 unless the lower-order coefficients are known to vanish; the paper does not provide that information. This matters because at the smallest L values the omitted 1/L term is numerically larger than the terms fitted. If that term is present, the extrapolated infinite-volume values used to draw the O(2) line are biased, and the conclusion for H≤1/160 weakens. I did not find an internal inconsistency: the construction of M in Eq. (1) could indeed cancel the naive leading FV terms, and the agreement visible for H=1/160 and 1/240 is suggestive. That is why the verdict remains conditional rather than reject. The paper is a proceedings-style preliminary result; the absence of released data and error bars makes the requested test essential. If the check shows the 1/L coefficients vanish (or that adding them does not move a0 outside errors), the central claim would be substantially supported.","tokens_in":4584,"tokens_out":6341,"duration_ms":68655,"concrete_test":"Re-analyze the H=1/160 (and, if available, H=1/240) data at T=145.1 MeV with the extended ansatz f=a0+a1/L+a2/L^2+a3/L^3+a4/L^4; compare a0 with the Eq. (5) result. Independently, compute the first Taylor coefficient(s) a_{0m} of the O(2) finite-volume scaling function fGχ(z_L) from Ref. [5] with the same conventions; if a01 or a02 is nonzero, Eq. (5) is biased and the verdict should be conditional until fits include those terms. Also report bootstrap errors on a0 for both ansätze.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that M/H^{1/delta} follows the O(2) finite-volume scaling function for H≤1/160—depends on the infinite-volume extrapolation in Fig. 2. That extrapolation uses Eq. (5), a0+a3/L^3+a4/L^4, but the scaling form Eq. (4) is an expansion in z_L = l0 h^{-νc}/L, i.e., generically a power series in 1/L. Eq. (5) simply omits the 1/L and 1/L^2 terms. The paper never states or justifies that the first finite-volume coefficients a_{0m} in Eq. (4) vanish for the subtracted order parameter M. If they do not, the fit is systematically biased: for the smallest volumes used here (Nσ/Nτ=3, i.e. 1/L≈0.33), the omitted terms are orders of magnitude larger than the retained 1/L^3 and 1/L^4 contributions (≈0.037 and 0.012). This would shift the z_L=0 points in Fig. 2 and could alter the apparent agreement with O(2) scaling at H=1/160 and 1/240. No errors, stability checks, or leave-one-volume-out tests are shown, so the extrapolation is not controlled in the sense claimed. Msub is likewise assumed negligible without a quantitative bound.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a finite-volume scaling analysis of the subtracted chiral order parameter M = M_ℓ − H χ_ℓ on N_τ = 8 lattices in (2+1)-flavor QCD with HISQ action, for light-to-strange mass ratios H down to 1/240. The authors compare M/H^{1/δ} with the 3-d O(2) finite-volume scaling function f_{Gχ}(z,z_L) of Ref. [5] and claim that for H ≤ 1/160, the data follow this universal scaling function, with finite-volume effects subdued and sub-leading contributions small. They quantify deviations for larger H and emphasize the preliminary nature of the results, calling for additional data and future N_τ = 12 simulations.","tokens_in":5102,"tokens_out":3263,"duration_ms":35257,"significance":"If the central claim holds, the work provides a potentially important step toward establishing O(2) scaling in (2+1)-flavor QCD at finite lattice spacing, and the proposed method of joint infinite-volume and chiral extrapolations could sharpen determinations of T_c and the chiral critical region. The paper introduces a clean subtracted order parameter, explicitly compares to an independent scaling function, and shows that the large-H data do not fall on the universal curves—so the comparison is not empty. However, the reported results are preliminary: no statistical errors are shown, the universal scale parameters are preliminary, and the infinite-volume extrapolation relies on an ansatz whose validity is not demonstrated. These issues currently prevent the stated conclusion from being fully supported.","major_comments":[{"comment":"The infinite-volume ansatz f = a0 + a3/L^3 + a4/L^4 drops all 1/L and 1/L^2 terms that are generically present in the z_L expansion of Eq. (4). The paper does not justify that the finite-volume coefficients for those orders vanish for the subtracted order parameter M. For the smallest lattices used (N_σ/N_τ = 3), 1/L ≈ 0.33, so an omitted coefficient of order unity would be an order of magnitude larger than the retained 1/L^3 and 1/L^4 terms. Such a bias would shift the z_L = 0 points in Fig. 2 and could artificially create the apparent agreement with O(2) scaling at H = 1/160 and 1/240. Please justify the truncation, show a stability check with 1/L and 1/L^2 terms included, or demonstrate that the corresponding coefficients are consistent with zero.","section":"§3, Eq. (5)"},{"comment":"No statistical error bars are shown on any data point, on the infinite-volume extrapolated values, or on the fitted scaling curves. The statement in the Summary that the analysis uses 'controlled infinite-volume extrapolation' cannot be assessed without uncertainties. Please provide at least typical statistical errors on M, and for the fits in Eq. (5) include χ²/dof, confidence intervals on a0/a3/a4, and ideally a leave-one-volume-out test. The '0.5 MeV error band' on T_c in Fig. 2 also appears without any derivation; its source and meaning should be made explicit.","section":"§3, Figs. 2 and 3"},{"comment":"The sub-leading contribution M_sub is assumed negligible, but no quantitative bound is given. Since the central claim is restricted to H ≤ 1/160, the paper should demonstrate that corrections-to-scaling and regular terms are smaller than the scatter of the data. This could be tested by checking whether the ratio M/H^{1/δ} at fixed z and z_L is independent of H, or by including an M_sub term in the global fit and showing it to be small.","section":"§2, Eq. (2)"},{"comment":"The non-universal scale parameters T_c = 145.1 MeV, z_0 = 1.52, and z_{L,0} = 0.38 are described as 'preliminary', but it is not stated whether they were determined from the data shown here or taken from an external fit. If they are adjusted to maximize the collapse of the small-H data, the comparison in Fig. 3 loses some of its predictive power. Please state clearly how these parameters were obtained, what their uncertainties are, and whether the O(2) agreement persists for fixed, a priori chosen values.","section":"§3, Fig. 3 and text"}],"minor_comments":[{"comment":"There are several typographical and notation issues: 'vale' should be 'value'; 'the to 3-d O(2)' should be 'the 3-d O(2)'; 'e ffects' and '1 /240' have spacing artifacts. The summation ranges in Eq. (4) are garbled ('muX', 'ml') and the indices m, n are never defined; please clarify the Taylor expansion and its limits.","section":"§3, text and figures"},{"comment":"The statement that data for T = 142.8 and 147.4 MeV have been shifted by 'single, constant values' needs a brief explanation of how those constants were determined; otherwise the apparent coincidence across temperatures is not independently checkable.","section":"§3, left of Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"This manuscript reads like a proceedings contribution: it explicitly says results are preliminary and more data are needed. If the journal's standard for publication requires a self-contained, statistically supported analysis, the lack of error bars and the arbitrary 1/L truncation in Eq. (5) are serious. The central idea is interesting, but the paper currently does not provide the 'controlled extrapolation' it claims. The issues are fixable with additional analysis or with a substantially more cautious wording, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read: this is a proceedings-style preprint with genuinely new data at H=1/240 on N_tau=8 lattices, and the authors make a systematic attempt to extract the infinite-volume limit of the subtracted order parameter. The headline claim, that for H<=1/160 the data fall on the 3-d O(2) finite-volume scaling function, is plausible from the figures. Credit where due: the new H=1/240 data are new, the subtracted order parameter is a sensible observable, and the figures clearly show that H>=1/80 deviates from universal scaling, which makes the comparison non-empty. The authors also flag their own caution: they say more 1/240 data at larger volumes are needed.\n\nThe soft spots are real. The infinite-volume extrapolation uses Eq. (5), a polynomial in 1/L^3 and 1/L^4. But the scaling form Eq. (4) is a generic power series in z_L, i.e. in 1/L, and nothing in the paper establishes that the 1/L and 1/L^2 coefficients vanish for the subtracted order parameter. For the smallest volumes shown (N_sigma/N_tau=3) the omitted 1/L term is an order of magnitude larger than the retained 1/L^3 term, so the extrapolation could be systematically biased. That directly affects the claimed agreement for H=1/160 and 1/240. The paper also shows no statistical errors, no fit quality, and no leave-one-volume-out checks. M_sub is assumed negligible without a quantitative bound. The scale parameters (T_c, z0, z_L0) are preliminary and fitted to the same data, which adds another layer of uncertainty, though the fact that H>1/80 data visibly miss the curves keeps the comparison from being circular.\n\nThe paper is honest about its status: it is a proceedings contribution, clearly labeled preliminary, and the authors call for more data. But the phrase 'controlled infinite-volume extrapolation' is stronger than what the analysis supports. I would like to see the authors justify or amend Eq. (5), give errors, and ideally release the data. If they can show that the 1/L and 1/L^2 terms are suppressed or account for them, the scaling conclusion would be much firmer.\n\nWho is this for? Lattice QCD specialists working on the chiral transition and scaling analyses. It's a useful incremental contribution, not a paradigm shift. I'd send it to a referee: it deserves careful review, and a good referee would push for the missing justifications.","headline":"New H=1/240 lattice data for the chiral order parameter look plausibly O(2)-scaling, but the infinite-volume fit drops the leading 1/L terms from the scaling form, so the conclusion is not yet controlled.","tokens_in":5438,"tokens_out":2924,"would_cite":true,"duration_ms":29037,"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 reports that after divergence subtraction and infinite-volume extrapolation, the chiral order parameter in (2+1)-flavor QCD follows the 3-d O(2) finite-volume scaling curve for light-to-strange quark mass ratios at or below 1/160.","keywords":["lattice QCD","chiral phase transition","finite-volume scaling","O(2) universality class","chiral order parameter","subtracted condensate","critical exponents","quark mass ratio"],"falsifier":"A direct test would use the same subtracted order parameter and infinite-volume extrapolation on Nτ=12 (and Nτ=16) lattices at H=1/160 and H=1/240 over the same z_L range; if the extrapolated M/H^{1/δ} no longer collapses onto f_Gχ(z,z_L), the Nτ=8 result is a finite-spacing artifact. Even without new simulations, refitting the existing data with an additional 1/L^5 term in Eq. (5) and checking whether the fitted O(2) parameters shift beyond the quoted 2% would expose sensitivity to the volume ansatz.","tokens_in":4538,"feed_emoji":"⚛️","tokens_out":12706,"duration_ms":122344,"temperature":0.7,"pith_summary":"The paper asks whether the chiral phase transition in (2+1)-flavor QCD shows universal finite-volume scaling when the light quark mass is small enough. Using an improved chiral order parameter that removes additive ultraviolet and linear quark-mass contributions, the authors compare lattice data on Nτ=8 volumes with the finite-volume scaling function of the 3-d O(2) universality class. They find that for light-to-strange mass ratios H ≤ 1/160, the scaled order parameter M/H^{1/δ} follows the universal curve after an infinite-volume extrapolation, while H ≥ 1/80 deviates visibly. This matters because a controlled scaling analysis is the route to extracting critical exponents, identifying the correct universality class, and pinning down the transition temperature and the QCD critical point. The result is preliminary in that additional large-volume H=1/240 data and Nτ=12 simulations are needed to confirm the trend.","feed_headline":"QCD order parameter obeys universal scaling at small quark mass","feed_subtitle":"At mass ratio 1/160 the scaled condensate falls on one universal curve, pinning down the chiral transition temperature.","key_machinery":"The machinery is the subtracted chiral order parameter M = M_l − Hχ_l (Eq. 1), defined so that additive ultraviolet and linear-H regular contributions cancel. The universal part is M = h^{1/δ} f_Gχ(z,z_L) + M_sub, with scaling variables z = t h^{-1/βδ} and z_L = l h^{-ν_c} (Eqs. 2–3). The comparison uses the 3-d O(2) finite-volume scaling function f_Gχ(z,z_L), expanded as a Taylor series in z and z_L (Eq. 4), and an infinite-volume extrapolation ansatz f = a0 + a3/L^3 + a4/L^4 (Eq. 5) that does not assume universality. The key step is that, for fixed T and H, the 1/L^3 and 1/L^4 volume terms are removed, and the remaining z_L dependence is compared directly with the universal O(2) prediction","core_discovery":"The central claim is that a systematic finite-volume analysis of the improved chiral order parameter M = M_l − Hχ_l on Nτ=8 lattices exposes universal 3-d O(2) scaling when the light quark mass is sufficiently small. Extrapolating fixed-temperature, fixed-H data to infinite volume with a 1/L^3 − 1/L^4 ansatz, the ratio M/H^{1/δ} collapses onto the O(2) finite-volume scaling function f_Gχ(z,z_L) for H = 1/160 and H = 1/240 over a temperature window centered on Tc ≃ 145 MeV, with deviations below the 2% level. For H ≥ 1/80, the same ratio departs visibly from the universal curve, and the volume required to reach a given accuracy grows with H. The authors interpret this as evidence that the sca","pith_inferences":["A direct corollary the authors do not spell out is that the same subtracted order parameter could be used to test O(4) or U(2)×U(2) scaling functions directly; whichever function also collapses the Nτ=12 data would identify the continuum universality class without relying on critical-exponent fits.","The apparent temperature independence of finite-volume effects between 142.8 and 147.4 MeV suggests a single z_L scaling curve may describe the whole transition region; one could predict M(T,H,L) at unmeasured temperatures and verify with the existing H=1/240 data.","Because the volume ansatz Eq. (5) is deliberately non-universal, its fitted coefficients a3 and a4 should vanish in the chiral and infinite-volume limits; tracking how they behave as H → 1/240 would provide a consistency check that the observed collapse is not an artifact of the extrapolation."],"forward_implications":["For H ≤ 1/160, controlled infinite-volume extrapolations of M can be performed with modest volumes, making Tc estimates on Nτ=8 lattices more precise and placing a tighter bound on the QCD critical point.","The size of the scaling regime becomes quantitative: reaching about 2% accuracy in M/H^{1/δ} requires an aspect ratio Nσ/Nτ ≈ 7 at H=1/160, versus ≈5 at H=1/27, guiding where future simulations must be run.","The analysis provides a template for extracting the critical exponents β, δ, and ν from the order parameter, and thereby for discriminating among universality classes such as O(2), O(4), and U(2)×U(2) in the continuum limit.","It validates the use of 3-d O(2) finite-volume scaling functions for staggered fermions at finite lattice spacing, a prerequisite for extrapolating to Nτ=12 and for investigating the fate of U(1)_A symmetry."],"supporting_citations":[{"why":"Provides the earlier (2+1)-flavor QCD results and non-universal scale parameters for pseudo-critical-temperature analyses that this work extends to smaller quark masses.","marker":"[1]"},{"why":"Defines the subtracted order parameter M = M_l − Hχ_l that cancels additive ultraviolet and linear-H contributions, the observable on which the whole analysis is based.","marker":"[3, 4]"},{"why":"Gives the 3-d O(2) finite-volume scaling function and its Taylor expansion in z and z_L, the universal curve used for comparison, as well as the preliminary non-universal scale parameters.","marker":"[5]"},{"why":"Supports the absence of conclusive first-order transition evidence in (2+1)-flavor QCD, justifying the use of universal scaling functions.","marker":"[8]"},{"why":"Describes the lattice simulation code used to generate the new H=1/240 and extended Nτ=8 data sets analyzed here.","marker":"[9, 10]"}],"fun_headline_variants":["Universal scaling in QCD at quark mass ratio 1/160","Chiral order parameter collapses to universal curve","Finite-volume analysis confirms O(2) scaling in QCD","Small quark mass yields universal finite-size scaling","QCD chiral transition shows universal scaling at light quark"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The whole conclusion rests on the assumption that the data for H ≤ 1/160 already lie inside the 3-d O(2) scaling window, with the fitted 1/L^3–1/L^4 volume terms and the truncated Taylor expansion of the scaling function absorbing every non-universal correction; if a neglected correction mimics the O(2) curve, the claim fails.","fun_headline_variants_meta":{"raw":{"variants":["Universal scaling in QCD at quark mass ratio 1/160","Chiral order parameter collapses to universal curve","Finite-volume analysis confirms O(2) scaling in QCD","Small quark mass yields universal finite-size scaling","QCD chiral transition shows universal scaling at light quark"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000177,"raw_usage":{"total_tokens":1102,"prompt_tokens":687,"completion_tokens":415,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":431,"completion_tokens_details":{"reasoning_tokens":337}},"tokens_in":431,"tokens_out":415,"duration_ms":4267,"temperature":1.0,"reasoning_tokens":337,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T15:48:07.992859+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would use the same subtracted order parameter and infinite-volume extrapolation on Nτ=12 (and Nτ=16) lattices at H=1/160 and H=1/240 over the same z_L range; if the extrapolated M/H^{1/δ} no longer collapses onto f_Gχ(z,z_L), the Nτ=8 result is a finite-spacing artifact. Even without new simulations, refitting the existing data with an additional 1/L^5 term in Eq. (5) and checking whether the fitted O(2) parameters shift beyond the quoted 2% would expose sensitivity to the volume ansatz.","supporting_citations":[{"cited_title":"Scaling functions of the three-dimensional $Z(2)$, $O(2)$ and $O(4)$ models and their finite size dependence in an external field","cited_arxiv_id":"2304.01710","evidence_quote":"Gives the 3-d O(2) finite-volume scaling function and its Taylor expansion in z and z_L, the universal curve used for comparison, as well as the preliminary non-universal scale parameters."}],"review_version":1}