{"id":"b29001e7-a185-4125-8ecc-86b4fce6fc6e","arxiv_id":"2501.19350","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The ETMC collaboration presents preliminary electro-quenched valence quark-connected isospin-breaking corrections to the leading-order hadronic vacuum polarization contribution to the muon anomalous magnetic moment on two lattice volumes at a single lattice spacing.","lead":"This paper reports preliminary lattice QCD calculations of the valence isospin-breaking corrections to the muon's anomalous magnetic moment, a key ingredient in the hadronic vacuum polarization contribution. It is part of the ETMC collaboration's effort to independently compute a quantity that helps test the Standard Model against the Fermilab muon g-2 measurement.","discovery_kind":"replication","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The linear r_m extrapolation in Eq. (18) from r_m = 3–9 to r_m = 1 is the softest point; subset fits do not test curvature below the lowest simulated quark mass, so the Table 3 light values may carry an unquantified chiral bias.","rationale":"This is an explicit feasibility/proceedings report with a well-established RM123 framework, and its claims are scoped to electro-quenched, valence quark-connected corrections. I found no internal inconsistency in the counterterm procedure: counterterms are fixed by meson masses, not by a_mu, so there is no circularity. The main numerical dependence of the headline light correction goes through the chiral extrapolation of Eq. (18). The reader's weakest assumption is essentially correct, but the statement that the paper provides 'no curvature test or alternative ansatz' is too strong: the subset fits (3,5,7) versus (5,7,9) are a partial linearity check. They cannot, however, constrain curvature below r_m = 3, and the quoted errors are statistical only. Therefore the reader's CONDITIONAL verdict remains appropriate; the concern is real but already identified, and it is not severe enough to move the verdict to REJECT.","tokens_in":9453,"tokens_out":4184,"duration_ms":45496,"concrete_test":"Reanalyse the B64 data of Fig. 3 with a quadratic ansatz Δa_mu^HVP(l; t_cut, r_m) = A(t_cut) + c1 r_m + c2 r_m^2 and evaluate at r_m = 1; if c2 is compatible with zero while A changes by less than the quoted statistical error, the linear ansatz is supported. As a stronger check, compute the light correction at r_m = 2 (or, ideally, at r_m = 1) and compare with the extrapolated value; agreement within one statistical sigma would settle the concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central light-quark numbers in Table 3 are obtained by Eq. (18), a linear chiral extrapolation in r_m evaluated at r_m = 1 using data only at r_m = 3, 5, 7, 9 (Fig. 3, bottom-left). None of the simulated points lies below three times the physical light mass. The alternative fits quoted in the text (r_m in [3,7] and [5,9]) check the stability of the slope over the fitted window, but they cannot detect curvature that sets in below r_m = 3. The QED and SIB pieces of ΔC_JJ(light) are separately shown in Fig. 3 and both vary steeply with r_m; there is no independent constraint (e.g. from chiral perturbation theory or a direct r_m ≈ 1 computation) on the functional form. The paper explicitly reports only statistical errors in Sec. 5, and the spread of the subset fits is used as a systematic estimate but is not propagated into the quoted errors in Table 3. If the true r_m dependence has downward curvature between 1 and 3, the intercept Δa_mu^HVP(l; t_cut) would shift by more than the 0.44/0.86 statistical errors, directly changing the light correction and therefore the headline result.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Lattice 2024 proceedings paper applies the RM123 expansion to compute valence, quark-connected, leading isospin-breaking corrections to the light, strange and charm contributions to a_mu^HVP-LO on two ETMC ensembles (B48 and B64) with lattice spacing a ~ 0.08 fm and linear sizes L ~ 3.8 fm and L ~ 5.1 fm. The counterterms are fixed by parity restoration for the critical mass shift and by the pi+, K+, K0 and D_s meson masses for the quark mass shifts, with QED finite-size effects on meson masses corrected via Eq. (10); the HVP correction is obtained from the time-momentum representation. For the light quark, the correction is computed at r_m = 3, 5, 7, 9 times the physical light mass and extrapolated linearly (Eq. (18)) to r_m = 1. The main results are Table 3: Delta a_mu^HVP(l) = 3.41(44) x 10^-10 (B48) and 4.79(86) x 10^-10 (B64); Delta a_mu^HVP(s) = 0.0049(10) and 0.0059(7) x 10^-10; Delta a_mu^HVP(c) = 0.1369(12) and 0.1363(11) x 10^-10. The paper explicitly states that all quoted errors are statistical only and that the results are preliminary.","tokens_in":9740,"tokens_out":9372,"duration_ms":89362,"significance":"The manuscript is a clear status report from an ongoing ETMC calculation. Its strengths are the use of the standard RM123 framework, an internally consistent counterterm setup that does not use a_mu itself as an input (so there is no circularity), two volumes for a first finite-size check, and precise and t_cut-stable strange and charm results. The authors are appropriately cautious in marking the results as preliminary. The main limitation is the light-quark chiral extrapolation: the central light values in Table 3 rest on a linear fit over a factor of three in quark mass with no data below r_m = 3, so the quoted statistical errors do not include the leading systematic risk. If the result holds up, the strange and charm numbers are useful intermediate checks, but the light values are not yet competitive with the final precision targets until the chiral form, continuum limit, and finite-size effects are addressed.","major_comments":[{"comment":"The central light-quark values in Table 3 are obtained by the linear ansatz Delta a_mu^HVP(l; t_cut, r_m) = Delta a_mu^HVP(l; t_cut) + c1 r_m, fitted to data at r_m = 3, 5, 7, 9 and evaluated at r_m = 1. No simulated point lies below three times the physical light mass, and the alternative fits with r_m in [3,7] and [5,9] only probe the slope inside the fitted window; they cannot detect curvature setting in below r_m = 3. Because the QED and SIB components of the integrand both vary steeply with r_m (Fig. 3, top panels), the linear form is a nontrivial assumption. Please add a curvature test (for example, a c2 r_m^2 term or a data point at r_m close to 1) or explicitly state that the Table 3 light values are model-dependent estimates pending such a check.","section":"§4.1, Eq. (18), Fig. 3"},{"comment":"The errors quoted in Table 3 are statistical only, yet the spread among the three chiral fits shown in Fig. 3 (bottom-left: 4.34(76), 4.84(95), 4.75(74) at t_cut = 2.54 on B64) is described as the estimated systematic error of the extrapolation but is not propagated into the final results. This makes it difficult to interpret the comparison in Sec. 5 of the paper's accuracy with that of other collaborations. Please either include the extrapolation systematic in the quoted uncertainties or clearly identify the Table 3 light values as central values of a preliminary analysis with an additional unquantified systematic.","section":"§5 and Table 3"},{"comment":"The two light results, 3.41(44) (B48) and 4.79(86) (B64), differ by 1.38 x 10^-10 against a combined statistical error of about 0.97 x 10^-10, i.e., only 1.4 sigma, not 'about two standard deviations' as stated. Since no finite-volume correction is applied to the light HVP correlator (Eq. (10) is used only for the meson masses in the counterterm determination), the volume dependence of the light value is unresolved. Please report the actual significance and comment on the implications for the final error budget.","section":"§4.1 and Table 3"}],"minor_comments":[{"comment":"The heading 'The LIB correntions' contains a typo; it should be 'corrections', and later in the same section 'respecively' should be 'respectively'.","section":"§4.1"},{"comment":"The heading 'Acknowlogments' is misspelled; it should be 'Acknowledgments'.","section":"Acknowledgments"},{"comment":"The electro-quenched approximation is never explicitly defined in the text; please add a sentence clarifying that it neglects QED effects on sea quarks and the associated determinant expansion.","section":"§2"}],"recommendation":"major_revision","confidential_remarks":"For the editor: this is a proceedings paper, and the requested changes are within the scope of a revision. The chiral extrapolation test is the key issue; if the authors prefer to keep the paper as a short status report, they should at least remove the comparison of their accuracy with final results of other collaborations and state explicitly that the light values are model-dependent. I see no reason to reject."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a Lattice 2024 proceedings from ETMC reporting preliminary valence-connected leading isospin breaking corrections to a_mu^HVP, split into light, strange, and charm pieces. The RM123 method is applied to two volumes at a single lattice spacing (a ~ 0.08 fm). The actual numbers in Table 3 are new for ETMC, but the quantity was already computed by BMW and by Giusti et al., so the value here is as an independent cross-check within an ongoing program, not a first result.\n\nWhat the paper does well: the counterterm setup is clean and non-circular. Equation (9) fixes the bare quark mass shifts from pion, kaon, and D_s masses; the critical mass shift comes from the parity condition; neither uses a_mu itself. The Z_V determination is shown with a clear plateaux. The paper is also honest about its scope: electro-quenched valence-connected only, statistical errors only, systematic assessment deferred, no direct comparison with BMW/Giusti yet. For a proceedings, that's the right level of candor.\n\nThe softest point is the light-quark chiral extrapolation. Eq. (18) is a linear ansatz in r_m, fit to r_m = 3, 5, 7, 9 and evaluated at r_m = 1. The subset fits (r_m in [3,7] and [5,9]) check the slope over the fitted window, but they cannot detect curvature between r_m = 1 and 3. The stress-test note is right about that. The two light results, 3.41(44) and 4.79(86) on B48 and B64, differ by about 1.4 sigma, which could be FSE or curvature, but with statistical errors only it's not diagnostic. The authors flag all of this explicitly, so it's a known limitation, not a hidden one.\n\nOther soft spots: one lattice spacing, no FSE correction for the light HVP, and the chiral intercept as quoted carries only the linear-fit error, with the subset-fit spread reported but not propagated into Table 3. For a journal submission, that would be a real barrier; for a proceedings, it's acceptable as long as the reader treats the numbers as preliminary.\n\nWho gets value: lattice HVP practitioners and anyone tracking the g-2 tension who wants a status update on ETMC's LIB program. It's a well-written, honest progress report. I'd be happy to see it refereed for the proceedings; the referee should check the internal consistency of the counterterm and extrapolation procedure, not expect a final number.","headline":"Honest, well-executed ETMC progress report on valence LIB to HVP; the light chiral extrapolation is the main soft spot, but the paper is transparent about it.","tokens_in":10313,"tokens_out":3150,"would_cite":true,"duration_ms":29527,"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":"This lattice calculation reports the valence-connected leading isospin-breaking correction to the muon HVP as $3.41(44)$ and $4.79(86)\\times10^{-10}$ for the light quarks, with strange and charm contributions about one and three orders of…","keywords":["lattice QCD","hadronic vacuum polarization","muon g-2","isospin breaking","QED corrections","RM123 approach","twisted mass fermions","chiral extrapolation"],"falsifier":"Compute $\\Delta a_\\mu^{\\mathrm{HVP}}(\\ell)$ directly at $r_m = 1$, the physical light-quark mass, with enough statistics to avoid the chiral extrapolation; if the direct value falls outside the linear-fit band reported in Table 3, the light correction is not what the paper quotes.","tokens_in":9223,"feed_emoji":"🧲","tokens_out":10528,"duration_ms":91218,"temperature":0.7,"pith_summary":"This paper reports a lattice QCD+QED calculation of the leading isospin-breaking corrections to the light, strange, and charm quark-connected contributions to the leading-order hadronic vacuum polarization (HVP) part of the muon anomalous magnetic moment, $a_\\mu^{\\mathrm{HVP}}$. The calculation uses the RM123 expansion, treating electromagnetic effects and the up-down quark mass difference as small perturbations around isospin-symmetric QCD to first order in $\\alpha_{\\mathrm{em}}$ and $(\\mu_d-\\mu_u)/\\Lambda_{\\mathrm{QCD}}$. At the current, explicitly preliminary stage, the light-quark correction is the dominant one, $\\Delta a_\\mu^{\\mathrm{HVP}}(\\ell) = 3.41(44)\\times 10^{-10}$ on the smaller volume and $4.79(86)\\times 10^{-10}$ on the larger, while the strange and charm corrections are one and three orders of magnitude smaller, respectively. The results are obtained at a single lattice spacing ($a \\sim 0.08$ fm) with statistical errors only, so the paper's contribution is a feasibility demonstration that this method can reach the sub-percent accuracy needed for the comparison with the Fermilab measurement of the muon $g-2$.","feed_headline":"Light-quark isospin correction to muon g-2: 3.4-4.8 ×10⁻¹⁰","feed_subtitle":"Valence-connected QED and up-down mass effects now match other lattice groups; systematics still to come.","key_machinery":"The machinery is the RM123 expansion, a first-order perturbation theory in $\\alpha_{\\mathrm{em}}$ and $\\mu_u-\\mu_d$ around an isospin-symmetric QCD ensemble: all QED and strong-isospin effects are obtained by differentiating correlation functions with respect to $e^2$, bare quark masses, and critical masses. The HVP integral uses the time-momentum representation with the analytic kernel $K(m_\\mu t)$, and the counterterms (the critical-mass shift $\\Delta \\bar m_{cr}$ and the bare quark-mass shifts $\\Delta \\bar\\mu$, $\\Delta \\mu_{ud}$, $\\Delta \\mu_s$, $\\Delta \\mu_c$) are fixed by parity-restoration conditions and by matching the $\\pi^+$, $K^+$, $K^0$, and $D_s$ masses, with QED finite-size effects removed by a known $1/L$ formula. The light-quark result requires the chiral extrapolation ansatz $\\Delta a_\\mu^{\\mathrm{HVP}}(\\ell; t_{\\mathrm{cut}}, r_m) = \\Delta a_\\mu^{\\mathrm{HVP}}(\\ell; t_{\\mathrm{cut}}) + c_1 r_m$, applied at each time cutoff $t_{\\mathrm{cut}}$.","core_discovery":"The central result is the set of valence, quark-connected leading isospin-breaking corrections computed in the electro-quenched approximation, in which sea quarks carry no electric charge and the lattice spacing is kept fixed: $\\Delta a_\\mu^{\\mathrm{HVP}}(\\ell) = 3.41(44)\\times 10^{-10}$ (B48) and $4.79(86)\\times 10^{-10}$ (B64), $\\Delta a_\\mu^{\\mathrm{HVP}}(s) = 0.0049(10)\\times 10^{-10}$ and $0.0059(7)\\times 10^{-10}$, and $\\Delta a_\\mu^{\\mathrm{HVP}}(c) = 0.1369(12)\\times 10^{-10}$ and $0.1363(11)\\times 10^{-10}$. The light correction is obtained by a linear chiral extrapolation in the quark-mass factor $r_m$ from $r_m = 3,5,7,9$ down to the physical point $r_m = 1$; the strange and charm corrections come from plateaux in $t_{\\mathrm{cut}}$ with no significant signal-to-noise degradation. The two volumes agree within about two standard deviations, and the quoted uncertainties are statistical only. The paper concludes that the accuracy is in line with earlier results by other collaborations and by the same collaboration, while a full account of systematic errors is deferred.","pith_inferences":["A natural test of the extrapolation is to add a curvature term, e.g. $c_2 r_m^2$, to Eq. (18) or to include data at $r_m = 2$; even with existing points, a quadratic fit would show whether the linear ansatz is safe.","The near-vanishing of the strange correction suggests that the future struggle for the LIB part of $a_\\mu^{\\mathrm{HVP}}$ will concentrate on the light connected contribution and on the disconnected and sea-quark terms, not on heavy quarks.","The time dependence of the light integrand shown in the paper could be used to form short-distance and long-distance window quantities, which would localize the B48-B64 discrepancy in Euclidean time and sharpen the comparison with the BMW result."],"forward_implications":["If these numbers hold, the light-quark connected correction shifts $a_\\mu^{\\mathrm{HVP}}$ by roughly $4\\times 10^{-10}$, a few permille effect that must be included in the theory prediction.","The strange correction is negligible at current precision (about $0.005\\times 10^{-10}$), while the charm correction is about $0.137\\times 10^{-10}$ and becomes relevant only at sub-permille total accuracy.","The two-volume agreement within about two standard deviations means residual finite-size effects are not yet controlled at the quoted statistical precision.","Because the calculation is electro-quenched and uses one lattice spacing, the present numbers are not final LIB corrections; sea-quark QED, disconnected, and continuum-extrapolated contributions are all still missing."],"supporting_citations":[{"why":"Proposed the $r_m$ chiral-extrapolation procedure and provides the reference BMW HVP value whose LIB corrections this work computes.","marker":"[1]"},{"why":"RM123 method: supplies the first-order expansion in $\\alpha_{\\mathrm{em}}$ and $(\\mu_d-\\mu_u)/\\Lambda_{\\mathrm{QCD}}$ that defines the calculation.","marker":"[4]"},{"why":"Parity-restoration conditions used to determine the critical-mass counterterms.","marker":"[7]"},{"why":"QED finite-size-effect formula applied to meson masses in the counterterm system.","marker":"[8]"},{"why":"Time-momentum representation that turns the vector-vector correlator into $a_\\mu^{\\mathrm{HVP}}$.","marker":"[9]"},{"why":"Provides the $Z_V$ renormalization procedure adapted to compute the QED correction $\\Delta Z_V/Z_V$.","marker":"[10]"},{"why":"Earlier ETMC computation of the same LIB corrections whose accuracy this work matches; the direct comparison baseline.","marker":"[11]"}],"fun_headline_variants":["Valence isospin breaking in muon g-2: 3.4-4.8e-10","Light-quark valence isospin correction to g-2: 3.4-4.8e-10","Valence quark isospin effect on muon g-2: 3.4-4.8e-10","Preliminary valence isospin correction to muon g-2: 3.4-4.8e-10","Muon g-2 isospin breaking: valence quarks give 3.4-4.8e-10"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the light-quark correction follows a straight line in the quark-mass factor $r_m$ over the range from $r_m = 3$ down to $r_m = 1$; if the true curve bends in that interval, the quoted central value shifts by more than its statistical error.","fun_headline_variants_meta":{"raw":{"variants":["Valence isospin breaking in muon g-2: 3.4-4.8e-10","Light-quark valence isospin correction to g-2: 3.4-4.8e-10","Valence quark isospin effect on muon g-2: 3.4-4.8e-10","Preliminary valence isospin correction to muon g-2: 3.4-4.8e-10","Muon g-2 isospin breaking: valence quarks give 3.4-4.8e-10"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001051,"raw_usage":{"total_tokens":4444,"prompt_tokens":1002,"completion_tokens":3442,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":618,"completion_tokens_details":{"reasoning_tokens":3292}},"tokens_in":618,"tokens_out":3442,"duration_ms":26795,"temperature":1.0,"reasoning_tokens":3292,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T20:26:49.313906+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $\\Delta a_\\mu^{\\mathrm{HVP}}(\\ell)$ directly at $r_m = 1$, the physical light-quark mass, with enough statistics to avoid the chiral extrapolation; if the direct value falls outside the linear-fit band reported in Table 3, the light correction is not what the paper quotes.","supporting_citations":[{"cited_title":"Sea quark QED effects and twisted mass fermions","cited_arxiv_id":"1612.02265","evidence_quote":"Parity-restoration conditions used to determine the critical-mass counterterms."}],"review_version":1}