{"id":"4604d311-7f94-4d41-9e41-8fcafdeecced","arxiv_id":"2501.16920","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The authors derive and apply RG equations for the lattice spacing in 2+1 clover fermion QCD, finding only a small quark mass dependence and smoother lattice spacing ratios.","lead":"This paper derives renormalisation group equations for 2+1 clover fermions and uses them, together with pion mass and gradient flow data at five lattice spacings, to study how the lattice spacing scales with coupling and quark mass. The value is procedural: a principled way to reduce one systematic error in lattice QCD continuum extrapolations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The mass-axis proxy linking X_pi^lat^2/X_t0^lat^2 to a fixed renormalised quark mass across beta is the load-bearing assumption; its direct evidence is shown at a single beta and O(a^2) artifacts in the ratio are untested.","rationale":"I read the central claim as: the lattice spacing ratios obtained by matching the proxy y follow the RG equations, with a mass correction of only a few percent. The paper's derivation in Sec. 2 is clean, and the near-flatness of the ratios in y (Fig. 2, right) is visible even before the beta-function fit. The beta-function model (Eq. (27)) and the one-parameter D fit are secondary: the observed flatness in y is essentially a comparison of slopes in the fits of Fig. 2 (left). The weakest link is the calibration of the mass axis. Eq. (21) expresses y as a function of mbar_rgi only for a single beta; using it across beta presupposes that lattice artifacts cancel in the ratio to the claimed accuracy. The stationarity of X_s^lat^2 under flavour-symmetry breaking is demonstrated at beta = 5.50; it stabilises the singlet quantities when the average mass is fixed, but it does not test whether the ratio's relation to mbar_rgi is beta-independent. Because the reconstructed a^2 ratios and the D term are all anchored to this proxy, a few-percent violation would directly bias the final scaling plot. The paper is appropriately cautious, labels results as preliminary, and proposes more data and improved fits; my concern does not change the verdict, but it does justify the CONDITIONAL status and motivates the specific cross-proxy check described above.","tokens_in":8664,"tokens_out":20334,"duration_ms":175131,"concrete_test":"Test the proxy by repeating the Section 4-5 analysis with an independent flavour-singlet mass proxy, X_N^lat^2/X_t0^lat^2, using the X_N data already shown in Fig. 1. Apply the same fits x = y(A + B y) and Eq. (24) to extract a^2(beta, mbar_rgi)/a^2(beta_ref, mbar_rgi) for the nucleon proxy. If the resulting ratios agree with the pion-proxy results within the quoted few-percent error at beta = 5.40 and ~1-2% at finer beta, the mass-axis calibration is validated; if they disagree, the proxy assumption fails and the inferred scaling ratios are biased.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central comparison across beta relies on interpreting lines of constant y = X_pi^lat^2/X_t0^lat^2 as lines of constant renormalised quark mass mbar_rgi (Eq. (21)). This requires y to be a beta-independent function of mbar_rgi at the few-percent accuracy claimed. The paper's support is the assertion that X_s^lat^2 is constant along a mbar = const lines (Sec. 3); however, the direct evidence in Fig. 1 (right) is at a single beta = 5.50. That stationarity under flavour-symmetry breaking at fixed beta is necessary but not sufficient: the ratio X_pi^lat^2/X_t0^lat^2 can still carry O(a^2) discretisation errors that differ by beta, so equal y at different beta need not correspond to equal mbar_rgi. If the proxy is biased at the few-percent level, the interpolated lattice spacing ratios extracted from x = y(A + B y), the coefficients s^2(beta)/s^2(beta_ref) and D_pi/t0(beta,beta_ref), and the final scaling plot in Fig. 4 all inherit the bias. This is a calibration concern about the mass axis, distinct from the beta-function model in Eq. (27), and it is the most load-bearing uncertainty in the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Lattice 2024 proceedings contribution derives renormalisation-group equations for the lattice spacing of 2+1 clover fermions, including the leading O(am_q) mass-dependence terms, and applies them to five lattice spacings (beta = 5.40, 5.50, 5.65, 5.80, 5.95) using pion mass and gradient-flow observables. The central construction is the RG solution s(g0) for the beta function and r(g0) for the mass dependence of the lattice spacing, Eqs. (12)-(16). On the data side, the authors use y = X_pi^lat^2 / X_t0^lat^2 as a proxy for the RGI quark mass, fit x = X_pi^lat^2 = y(A+By) at each beta, extract s^2(beta)/s^2(beta_ref) and the mass-slope coefficient D_{pi/t0}, fit the beta function with a [2/2] Pade, and reconstruct a^2(beta,mbar*)/a^2(beta_ref,mbar*) at two values of mbar*. They find only a small mass dependence, mainly at beta = 5.40, and a smoother lattice-spacing curve than their previous separate determinations.","tokens_in":8882,"tokens_out":4622,"duration_ms":41206,"significance":"The paper proposes a conceptually attractive way to use RG equations to constrain lattice-spacing ratios across beta, which could lead to smoother continuum extrapolations. The derivation in Section 2 is clean and standard: the equations for u(g0), v(g0), s(g0) and r(g0) are internally consistent, and the reduction to the usual beta-function form is explicit. The authors are also honest that the results are preliminary and that the mass proxy is a calibration assumption. If the cross-beta calibration of y with mbar_rgi is confirmed at the few-percent level, the claimed scaling result is a useful methodological step. The paper does not overclaim an independent prediction, and the data are standard for the collaboration, though no reproducibility package is provided for this proceedings article.","major_comments":[{"comment":"The identification of y = X_pi^lat^2 / X_t0^lat^2 with a beta-independent function of the renormalised quark mass is the load-bearing assumption of the analysis, and the evidence offered for it is not sufficient for the few-percent accuracy claimed. The stationarity of X_s^lat under flavour-symmetry breaking (Fig. 1, left) is demonstrated at a single beta = 5.50 and only shows that X_s^lat is constant along lines of constant a mbar at that beta; it does not control O(a^2) discretisation errors in the ratio X_pi^lat^2 / X_t0^lat^2 that can differ from one beta to another. Consequently, equal y at different beta values need not correspond to equal mbar_rgi, and the interpolated values of x, the ratios s^2(beta)/s^2(beta_ref), the coefficient D_{pi/t0}, and the final scaling plot in Fig. 4 would all inherit the resulting bias. I would like to see a quantitative estimate of this effect, for example by comparing y against a directly determined mbar_rgi at two or more beta values, or by repeating the analysis after dropping one beta and checking the stability of the extracted ratios.","section":"Section 3, Eq. (21), Fig. 2"},{"comment":"The final a^2 ratios are reconstructed from the same fits that define the parameters entering Eq. (24): A and B per beta from the x = y(A+By) fits, b2^eff from the beta-function fit, and the proportionality constant in Eq. (28). The agreement between the RG-guided curve and the magenta points from QCDSF15 is therefore a consistency check of the fit model rather than an independent determination, and the manuscript should say so explicitly. To make the claim of a smoother lattice-spacing set quantitative, the authors should report the uncertainty on the reconstructed ratios that propagates from all fit parameters, and/or perform a leave-one-beta-out exercise.","section":"Section 5, Eq. (24), Fig. 4"},{"comment":"The [2/2] Pade form for the beta function is introduced as the simplest choice, but the sensitivity of s^2(beta)/s^2(beta_ref) and D_{pi/t0} to this model choice is not quantified. Since b2^eff is an effective parameter and the two-loop result is visibly different in Fig. 3, the analysis should include an estimate of the model uncertainty, for instance by using alternative Pade orders or by allowing b2 to vary within a reasonable range. Without such an estimate, the fit to B0(g0) cannot be distinguished from a simple empirical interpolation.","section":"Section 5, Eq. (27)"}],"minor_comments":[{"comment":"In the last expression of Eq. (28), the integration variable is written as 1/(2 b0 g0^2) without stating the substitution; please define it explicitly so that the incomplete-Gamma form is transparent.","section":"Section 5, Eq. (28)"},{"comment":"The quantity X_pi^2(mbar_rgi)/X_t0^2(mbar_rgi) is used in Eq. (24) and in Fig. 4 before being defined; a sentence introducing this ratio as the continuum counterpart of the proxy y would improve readability.","section":"Section 4, Eq. (22)"},{"comment":"The statement that X_s^lat^2 is constant along any a mbar = const line is made without an error budget; please quote the size of the observed scatter at beta = 5.50 and, if available, at other beta values.","section":"Section 3, paragraph after Fig. 1"},{"comment":"Reference [8] is given as an arXiv preprint; if a published version exists, it should be cited instead or in addition.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is an explicitly preliminary proceedings contribution, and the referee report should not be read as demanding a full journal-scale study. The main unresolved issue, the cross-beta calibration of the mass proxy, is genuine and load-bearing for the central scaling claim; a targeted test would materially strengthen the paper. The relationship to reference [1] (in preparation) is not assessable from this text, but the present contribution is self-contained enough for a proceedings paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful core is the leading-order RG solution for 2+1 clover fermions and the first application to QCDSF data at five beta values. The derivation is standard and clean, the data analysis is direct, and the authors are honest about what is preliminary.\n\nWhat is actually new: the RG equations and the b_g/r relation are known and properly cited; the novelty is the explicit leading-order solution specialized to 2+1 clover fermions and the first construction of RG-guided lattice spacing ratios from five beta values. That is a legitimate extension, not a breakthrough, but a useful one for scale setting.\n\nWhat the paper does well: the right panel of Fig. 2 shows directly that the mass dependence of the lattice spacing ratio is small, with the largest deviation at the coarsest beta. The RG-guided ratios in Fig. 4 look smoother than the separate determinations, which is exactly what the method is for. The beta-function Padé and the D coefficient are clearly presented as fits, with the two-loop curve shown for comparison. The paper does not overclaim; it labels the results preliminary and lists concrete next steps.\n\nThe soft spot is the mass-axis proxy. The central comparison across beta treats y = X_pi^lat^2/X_t0^lat^2 as a beta-independent proxy for the renormalised quark mass. The support in Sec. 3 is the claim that X_s^lat^2 is constant along lines of constant a mbar, but the evidence shown is at a single beta (5.50). Stationarity under flavour-symmetry breaking at fixed beta is necessary but not sufficient: the ratio can still carry O(a^2) discretisation errors that differ by beta, so equal y at different beta need not correspond to equal mbar_rgi. If that proxy is biased at the few-percent level, the interpolated lattice spacing ratios, the D_pi/t0 term, and the final scaling plot all inherit the bias. This is a calibration concern about the mass axis, distinct from the beta-function model in Eq. (27). The beta-function model is a real modeling choice, but less worrying because it is explicitly a fit and the alternative two-loop result is shown.\n\nThere are no critical red flags. The final ratios are reconstructions of fitted parameters rather than independent predictions, but the paper does not claim otherwise. It ships no code or data, which is normal for a proceedings talk.\n\nWho this is for: lattice QCD practitioners working on scale setting and continuum extrapolations. It deserves a serious referee as a proceedings contribution, and I would accept it for peer review. Before the method is used to set the scale, the mass-axis proxy should be tested at multiple beta values, ideally with a direct determination of mbar from the hopping parameter critical line. I would cite it for the RG-guided ratios, but not yet rely on the absolute scale.","headline":"A clean, honest RG treatment of lattice spacing scaling for 2+1 clover fermions, with a load-bearing assumption about the mass-axis proxy that needs testing at more beta values.","tokens_in":9534,"tokens_out":2080,"would_cite":true,"duration_ms":17835,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Gc"],"model":"deepseek-v4-flash","headline":"The paper derives RG equations for 2+1 clover fermions and shows the lattice spacing scales with little dependence on the renormalized quark mass.","keywords":["lattice QCD","renormalisation group","clover fermions","lattice spacing scaling","SU(3) flavour symmetry","gradient flow","pion mass","beta function"],"falsifier":"On existing ensembles, compute the ratio $a^2(\\beta,\\bar m_{\\rm rgi})/a^2(\\beta_{\\rm ref},\\bar m_{\\rm rgi})$ directly from separate lattice-spacing determinations at two values of $\\bar m_{\\rm rgi}$ that differ by about 20% in the proxy $X_\\pi^2/X_{t_0}^2$; the RG claim predicts the ratio changes by roughly 2–4%, so a change of 10% or more, or visible curvature in the x-versus-y plots, would contradict the scaling conclusion.","tokens_in":8376,"feed_emoji":"⚛️","tokens_out":9525,"duration_ms":74873,"temperature":0.7,"pith_summary":"This paper asks whether the lattice spacing of 2+1 flavor clover-fermion QCD follows the renormalization group as the bare coupling and quark mass are varied. The authors derive the RG equation for the lattice spacing along lines of constant physics, solve it in terms of the beta function and mass anomalous dimensions, and test it against pion mass and gradient-flow data at five lattice spacings. They find scaling: the lattice-spacing ratio between different beta values depends on the renormalized quark mass only at the few-percent level, so an error in choosing the physical quark mass mostly shifts the overall scale. If the result holds, lattice-spacing ratios can be fixed more precisely from accurate pion and t0 data, giving smoother continuum extrapolations.","feed_headline":"RG scaling holds for lattice spacings at five beta values","feed_subtitle":"Pion mass and gradient-flow data show quark-mass effects of a few percent, smoothing continuum extrapolations.","key_machinery":"The central object is the renormalization-group operator acting at fixed physics, $a\\,\\partial/\\partial a|_{\\rm physics}$, expressed in the bare coupling $g_0$ and lattice quark masses $a m_q$. At leading order in $a m_q$ it is controlled by the $\\beta$-function $B_0(g_0)=-b_0 g_0^3-\\cdots$, the mass anomalous dimension combination $G_0 = 1-\\gamma_m^{\\rm NS}$, the singlet–non-singlet difference $H_0=\\gamma_m^{\\rm NS}-\\gamma_m^{\\rm S}$, and the improvement coefficient $B_1 = b^{\\rm lat}_{10}g_0^3+\\cdots$ used to define the mass correction $r(g_0)$. Solving the RG equations gives $a = s(g_0)\\{1+a\\bar m\\, r(g_0)\\}$ with $s(g_0)$ fixed by the $\\beta$-function; the data analysis then rides on the identity that at common $\\bar m_{\\rm rgi}$ the measured ratio $x(\\beta)/x(\\beta_{\\rm ref})$ equals the lattice-spacing ratio $a^2(\\beta)/a^2(\\beta_{\\rm ref})$. The practical machinery is a two-parameter fit $x=y(A+By)$ in the mass-proxy $y=X_\\pi^{{\\rm lat}\\,2}/X_{t_0}^{{\\rm lat}\\,2}$, a $[2/2]$ Padé fit to the $\\beta$-function, and a one-parameter integral fit for the mass-correction coefficient $D_{\\pi/t_0}$.","core_discovery":"The central claim is that the lattice spacing in 2+1 flavor clover QCD obeys the renormalization group along lines of constant physics. Concretely, the spacing satisfies $a = s(g_0)\\{1 + a\\bar m\\, r(g_0)\\}$, where $s(g_0)$ is the standard $\\beta$-function solution and the mass correction $r(g_0)$ is small, of order $-b^{\\rm lat}_{10}g_0^2/b_0$ at weak coupling. Using the ratio $X_\\pi^{{\\rm lat}\\,2}/X_{t_0}^{{\\rm lat}\\,2}$ as a proxy for the renormalized quark mass $\\bar m_{\\rm rgi}$, the authors show that the ratio of lattice spacings at fixed $\\bar m_{\\rm rgi}$, $a^2(\\beta,\\bar m_{\\rm rgi})/a^2(\\beta_{\\rm ref},\\bar m_{\\rm rgi})$, is nearly flat: a 10% change in the mass proxy moves the coarsest ratio by only about 1–2%. Reconstructing the spacing ratio from $s^2(\\beta)/s^2(\\beta_{\\rm ref})$ plus a fitted linear term $D_{\\pi/t_0}$ in the mass proxy gives a smoother set of lattice spacings across $\\beta = 5.40, 5.50, 5.65, 5.80, 5.95$ than earlier separate determinations, leading the authors to conclude that scaling holds with little dependence on $\\bar m_{\\rm rgi}^*$.","pith_inferences":["If the scaling claim holds at finer couplings, the same RG-guided interpolation could set lattice spacings at intermediate beta values from accurate pion and t0 data alone, reducing the number of dedicated scale-setting runs.","The same machinery, with appropriate anomalous dimensions, could be applied to hadron matrix elements rather than just the lattice spacing, giving RG-guided joint continuum extrapolations.","A direct test of the constancy assumption would be to check whether $X_s^{{\\rm lat}\\,2}$ shows curvature along $a\\bar m = \\text{const}$ lines once more SU(3)-symmetric data at finer lattices are included; the paper's own plots show little evidence for it now.","Repeating the analysis with an independent third value of the physical mass proxy (beyond the two averages used here) would test whether the reported insensitivity to $\\bar m_{\\rm rgi}^*$ persists at the percent level."],"forward_implications":["Lattice spacing ratios across many beta values can be determined more precisely from accurate pion mass and gradient-flow data than from separate fits, since the RG equation fixes the beta-dependence.","An error in choosing the physical quark mass $\\bar m_{\\rm rgi}^*$ shifts the whole lattice-spacing set by a nearly constant overall factor, so absolute scale errors decouple from relative spacing errors.","The RG-guided ratios give a smoother lattice spacing curve across $\\beta = 5.40$–$5.95$, which is expected to reduce the noise in $a^2$ continuum extrapolations of hadron masses and matrix elements.","The numerically small $r(g_0)$ and $b_g(g_0)$ coefficients confirm the standard expectation that the two-loop $\\beta$-function describes the coupling dependence of the lattice spacing in this range, with corrections of a few percent at the coarsest lattice."],"supporting_citations":[{"why":"introduces the redefinition of the coupling used to derive the relation between b_g and r.","marker":"[2]"},{"why":"supports the expectation that b_g and r are numerically small.","marker":"[3]"},{"why":"updates that expectation with modern data.","marker":"[4]"},{"why":"introduces the SU(3)-flavour-singlet quantities whose constancy defines lines of constant physics.","marker":"[5]"},{"why":"provides the lattice data for the flavour singlet combinations plotted in Fig. 1.","marker":"[6]"},{"why":"provides the non-perturbatively O(a)-improved clover action used for the ensembles.","marker":"[7]"},{"why":"supplies the physical value of the mass proxy and the previous separate lattice-spacing determinations used for comparison.","marker":"[8]"}],"fun_headline_variants":["RG equations tame 2+1 clover lattice spacings","Clover fermion RG: scaling across five beta values","Five spacings show RG scaling for clover QCD","Mass-corrected RG smooths clover lattice spacing","Clover RG fixes scale across five beta values"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result rests on taking the flavour-singlet combination $X_\\pi^{{\\rm lat}\\,2}/X_{t_0}^{{\\rm lat}\\,2}$ as a clean stand-in for the renormalised quark mass: if that quantity bends along lines of constant $a\\bar m$ at the few-percent level, the derived spacing ratios and the scaling plot would be biased.","fun_headline_variants_meta":{"raw":{"variants":["RG equations tame 2+1 clover lattice spacings","Clover fermion RG: scaling across five beta values","Five spacings show RG scaling for clover QCD","Mass-corrected RG smooths clover lattice spacing","Clover RG fixes scale across five beta values"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000252,"raw_usage":{"total_tokens":1544,"prompt_tokens":914,"completion_tokens":630,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":530,"completion_tokens_details":{"reasoning_tokens":550}},"tokens_in":530,"tokens_out":630,"duration_ms":5249,"temperature":1.0,"reasoning_tokens":550,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T05:36:23.302301+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"On existing ensembles, compute the ratio $a^2(\\beta,\\bar m_{\\rm rgi})/a^2(\\beta_{\\rm ref},\\bar m_{\\rm rgi})$ directly from separate lattice-spacing determinations at two values of $\\bar m_{\\rm rgi}$ that differ by about 20% in the proxy $X_\\pi^2/X_{t_0}^2$; the RG claim predicts the ratio changes by roughly 2–4%, so a change of 10% or more, or visible curvature in the x-versus-y plots, would contradict the scaling conclusion.","supporting_citations":[],"review_version":1}