{"id":"b2633518-428b-48cd-a5df-5a8c14abb33f","arxiv_id":"1909.01962","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A first-principles lattice calculation with Nf=2+1+1 QCD plus quenched QED gives the leading isospin-breaking correction to the muon HVP as 7.1 (2.9) x 10^-10, the most precise such estimate to date.","lead":"This lattice QCD paper computes the electromagnetic and strong isospin-breaking corrections to the hadronic vacuum polarization contribution to the muon's anomalous magnetic moment, finding delta a_mu^HVP(udsc) = 7.1 (2.9) x 10^-10. The value is currently the most precise lattice estimate of this correction, which matters because the muon g-2 experiment may soon expose new physics beyond the Standard Model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central 'most accurate' claim depends on an estimated, not computed, qQED+disconnected systematic; a factor-3 error would overturn it.","rationale":"The paper is methodologically sound and the central value is supported by a detailed error budget, three lattice spacings, multiple volumes, and physical extrapolations. The RM123 approach is standard, and the local current avoids disconnected contributions by construction, but physical HVP requires those contributions. The single most load-bearing premise is that the qQED approximation plus neglect of disconnected diagrams introduces an uncertainty no larger than the QED connected contribution (1.0-1.2e-10). The paper provides no derivation for this estimate; it is a guess motivated by one RBC/UKQCD diagram. If the true effect is only a few times larger, the quoted total error of 2.9e-10 would still beat the 5.1e-10 dispersive estimate, but a threshold around 3.5e-10 exists where the 'most accurate' claim stops being true. The strange-sector 2-sigma deviation from RBC/UKQCD reinforces the possibility of unaccounted systematics, though it has little impact on the total since the strange contribution is tiny. The verdict should remain conditional: the result is credible and likely correct, but the headline precision and 'most accurate' status depend on a systematic that has not been computed. A dedicated calculation of the neglected QED sea and disconnected diagrams would settle this concern directly. We therefore see no reason to change the reader's conditional verdict.","tokens_in":14041,"tokens_out":5973,"duration_ms":62415,"concrete_test":"Compute the leading O(alpha_em) sea-quark QED corrections and the full set of quark-disconnected diagrams for delta a_mu^HVP(ud) on the ETMC B55.32 ensemble (a approx 0.082 fm, M_pi approx 375 MeV), using the same RM123 expansion and QED_L prescription. Compare the sum of these neglected contributions with the assigned 1.2e-10 systematic in Eq. (3.15); if the sum is larger than about 3.5e-10, the total uncertainty would exceed the BMW error and the 'most accurate' claim would not survive.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that Eq. (3.15) gives the most accurate determination of IB corrections—hinges on the 1.2e-10 uncertainty assigned to the qQED approximation and neglect of quark-disconnected diagrams. This is an order-of-magnitude estimate, not a calculation: the paper cites RBC/UKQCD's observation that one disconnected diagram is comparable to the QED connected term, and then assumes the combined qQED sea effect and all disconnected diagrams are similar in size. The qQED sea effect (treating dynamical sea quarks as electrically neutral) is not evaluated at all. If the true sum of these neglected contributions were about 3.5e-10, the total uncertainty in quadrature with the 2.6e-10 would exceed the BMW dispersive error of 5.1e-10, and the 'most accurate' claim would fail. The paper itself reports a 2-sigma deviation in the strange contribution from RBC/UKQCD (Sec. 3), which hints that unaccounted systematics may exist. This is not an internal inconsistency, but an unquantified systematic that directly controls the precision headline.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings contribution reports a lattice QCD calculation of the leading electromagnetic and strong isospin-breaking (IB) corrections to the quark-connected hadronic-vacuum-polarization (HVP) contribution to the muon magnetic anomaly. The calculation uses the RM123 expansion in the quenched-QED (qQED) approximation, ETMC N_f=2+1+1 configurations at three lattice spacings, several volumes, and pion masses between about 210 and 450 MeV, with QED_L zero-mode removal and nonperturbative renormalization constants. After chiral, continuum, and infinite-volume extrapolations the authors obtain delta a_mu^HVP(ud)=7.1(2.5) x 10^-10, delta a_mu^HVP(s)=-0.0053(33) x 10^-10, and delta a_mu^HVP(c)=0.0182(36) x 10^-10. Adding an estimated 1.2 x 10^-10 uncertainty for the qQED approximation and neglected quark-disconnected diagrams, they quote delta a_mu^HVP(udsc)=7.1(2.9) x 10^-10 and claim this is the most accurate determination of the IB corrections to a_mu^HVP to date.","tokens_in":14300,"tokens_out":8138,"duration_ms":82736,"significance":"If the quoted uncertainty is reliable, this is a valuable first-principles cross-check of both the dispersive estimates and other lattice determinations, improving on the BMW value 7.8(5.1) x 10^-10 and the RBC/UKQCD value 9.5(10.2) x 10^-10. The paper is methodologically transparent in several respects: it uses the time-momentum representation, an O(a)-improved twisted-mass setup, dedicated finite-volume ensembles, and nonperturbative renormalization constants from Ref. [27]. It also gives a detailed error budget separating statistical, input, chiral, finite-volume, and discretization effects. The main caveat is that the 'most accurate' claim depends on an estimated, rather than computed, uncertainty for the qQED approximation and disconnected diagrams; this is a load-bearing point that needs to be addressed before the headline claim can be accepted.","major_comments":[{"comment":"The central quantitative claim rests on an estimate rather than a calculation. The text states that the uncertainty from the qQED approximation and the neglect of quark-disconnected diagrams is 'approximately equal to our QED contribution (3.8)', i.e. about 1.1 x 10^-10, quoted as 1.2 x 10^-10 in Eq. (3.15). This estimate is not derived: the qQED sea-quark effect is not evaluated, and the cited RBC/UKQCD observation concerns a single disconnected diagram, not the full set of disconnected contractions nor the qQED sea effect. The 'most accurate' claim fails if the neglected effect is larger than estimated; for example, if the combined qQED/disconnected systematic were about 4.4 x 10^-10 (a factor of roughly four larger), the total uncertainty in Eq. (3.15) would become comparable to the BMW error of 5.1 x 10^-10, and the claimed superiority over BMW would disappear. I do not see an internal inconsistency, but the precision headline needs either a computed qQED/disconnected contribution, a conservative upper bound, or a revised wording that restricts the claim to the connected qQED contribution.","section":"Section 3, Eq. (3.15)"},{"comment":"The paper reports a roughly 2-sigma deviation from RBC/UKQCD for the strange-quark contribution, delta a_mu^HVP(s). Although the strange term is numerically small and does not by itself change the total central value, this tension between two lattice determinations is exactly the kind of signal that could indicate a missing systematic in one of the calculations, and it is directly relevant to the credibility of the estimated qQED/disconnected systematics. The revised manuscript should discuss possible sources of this discrepancy, such as renormalization constants, qQED sea effects, finite-volume effects, or the different treatment of disconnected diagrams.","section":"Section 3, Eq. (3.13)"},{"comment":"The comparison with BMW and RBC/UKQCD is not fully apples-to-apples. Eq. (3.15) is a connected-only result within the qQED approximation plus an estimated uncertainty, whereas the quoted BMW estimate is a full dispersive QCD+QED estimate and the RBC/UKQCD result includes at least one disconnected diagram. If the qQED/disconnected uncertainty remains an estimate, the phrase 'most accurate determination of the IB contribution to a_mu^HVP' should be either qualified to 'most accurate connected qQED determination' or backed by a calculation that places the neglected effects on the same footing as the other determinations.","section":"Section 3, Eqs. (3.15)-(3.16)"}],"minor_comments":[{"comment":"The text between Section 2 and Section 3 contains a large duplicated passage from Ref. [14], beginning with 'where mud 1/4...' and repeated several times. This appears to be a copy-paste error and should be removed from the submitted version.","section":"Section 2 (near Figure 1)"},{"comment":"The fit coefficients delta A_l1l, delta A_l2, delta D_l, and delta F_l in Eq. (3.4) are introduced only in prose; they should be defined explicitly in the equation or immediately after it so the fit ansatz is self-contained.","section":"Section 3, Eq. (3.4)"},{"comment":"The symbol 'qQED+disc' is used for an uncertainty, but the text says only that it is 'approximately equal to our QED contribution (3.8)'. Please state explicitly whether this is meant as a 1-sigma systematic uncertainty and justify the rounding from 1.1 x 10^-10 to 1.2 x 10^-10.","section":"Section 3, Eq. (3.15)"},{"comment":"The figure caption states that the errors are statistical only, but the final result in Eq. (3.6) includes several systematic errors. It would be useful to show the fully combined uncertainty on the physical-point marker, or to explain in the caption why the systematic errors are not displayed there.","section":"Figure 3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a proceedings contribution whose technical core is largely already published in Ref. [12]; the genuinely new element is the 'most accurate' claim and the associated qQED/disconnected uncertainty estimate. For the journal version I would ask for either a derivation or a more cautious claim. If the authors can provide a conservative bound or recompute the central quantity with the neglected effects included, the result could be acceptable. I have no concerns about citation practice or scope for a hep-lat venue."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a well-executed proceedings summary of a lattice result that already appeared in the authors' PRD paper. The genuinely new piece is the improved strange and charm corrections using nonperturbative renormalization from Ref. [27]. The paper is honest about its main weakness, but that weakness is real: the 'most accurate' claim hangs on an estimated, not computed, qQED and disconnected-diagram uncertainty.\n\nWhat's good: the lattice setup is standard and carefully controlled—three lattice spacings, several volumes, pion masses down to ~210 MeV, and a detailed error budget separating statistical, input, chiral, finite-volume, and discretization effects. The ratio method for the light quark contribution is sensible and reduces scale-setting systematics. The strange/charm update is a genuine improvement, cutting uncertainties by a factor of ~3. The paper also flags the ~2 sigma discrepancy with RBC/UKQCD in the strange channel rather than hiding it. That is exactly the right way to present a result that might indicate residual systematics.\n\nWhere the soft spots are: the qQED approximation and neglect of disconnected diagrams are assigned an uncertainty 'approximately equal to our QED contribution' (1.0e-10). That is an order-of-magnitude estimate, not a calculation. It might be right, and it is not crazy—RBC/UKQCD saw one disconnected diagram comparable to the QED connected term—but if the true combined effect is a few times larger, the 'most accurate determination' claim falls apart. The paper's own 2-sigma strange discrepancy is a small hint that something is not fully understood. Also, relative to Ref. [12] the core result is not new; this proceedings is a summary plus an incremental update. That is normal for a proceedings, but it means the novelty is modest.\n\nOverall: the central number is plausible and the methodology is solid, but the headline precision depends on an untested systematic assumption. The paper is transparent about this, which is to its credit.\n\nWho it's for: lattice QCD practitioners and people working on the SM prediction of g-2. It deserves a serious referee if submitted as a journal article, though the authors' own PRD is the primary reference. I would not cite the proceedings in preference to the PRD, but I would read it to understand the updated strange/charm numbers. Recommend: send to peer review; the qQED estimate should be stressed in the report.","headline":"A solid, transparent proceedings summary of an already-published lattice result; the 'most accurate' claim rests on an estimated qQED+disconnected uncertainty, but the paper is honest about it.","tokens_in":14870,"tokens_out":2252,"would_cite":false,"duration_ms":23051,"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 paper reports a lattice QCD calculation of the leading electromagnetic and strong isospin-breaking corrections to the hadronic-vacuum-polarization contribution to the muon anomalous magnetic moment, finding $\\delta a_\\mu^{\\rm…","keywords":["muon anomalous magnetic moment","hadronic vacuum polarization","isospin breaking","lattice QCD","RM123 method","quenched QED approximation","muon g-2","strong isospin breaking"],"falsifier":"Computing the quark-disconnected QED diagram for the light-quark contribution on the same gauge ensembles and comparing its size with the connected QED piece $1.1\\,(1.0)\\times10^{-10}$ would test the central systematic assumption; a disconnected contribution substantially larger than $1\\times10^{-10}$ would move $\\delta a_\\mu^{\\rm HVP}(udsc)$ outside the quoted total uncertainty of $2.9\\times10^{-10}$.","tokens_in":13803,"feed_emoji":"🧲","tokens_out":17231,"duration_ms":143131,"temperature":0.7,"pith_summary":"This paper reports a lattice QCD calculation of the leading electromagnetic and strong isospin-breaking corrections to the quark-connected hadronic-vacuum-polarization (HVP) contribution to the muon’s anomalous magnetic moment. Working in the RM123 expansion in powers of $\\alpha_{\\rm em}$ and $(m_d-m_u)/\\Lambda_{\\rm QCD}$, in the quenched-QED approximation, and using gauge ensembles at three lattice spacings, several volumes, and pion masses between about 220 and 490 MeV, the authors obtain $\\delta a_\\mu^{\\rm HVP}(udsc)=7.1\\,(2.9)\\times10^{-10}$ after continuum, infinite-volume, and physical-mass extrapolations. The light-quark correction dominates, $\\delta a_\\mu^{\\rm HVP}(ud)=7.1\\,(2.5)\\times10^{-10}$, about 85% of it from strong isospin breaking rather than QED; the strange and charm corrections are tiny. The paper states this is currently the most accurate determination of the isospin-breaking corrections to $a_\\mu^{\\rm HVP}$, and when added to leading-order HVP estimates it gives $a_\\mu^{\\rm HVP}=682\\,(19)\\times10^{-10}$, consistent with dispersive determinations. If the result holds, it tightens the first-principles HVP prediction at the precision target of the new muon $g-2$ experiments.","feed_headline":"Lattice QCD computes muon g-2 isospin shift: 7.1(2.9) x 10^-10","feed_subtitle":"Direct lattice QCD yields the leading electromagnetic and up-down mass corrections to hadronic vacuum polarization.","key_machinery":"The argument is carried by the RM123 method: an expansion of the QCD path integral in powers of the quark mass difference and of $\\alpha_{\\rm em}$ around an isosymmetric QCD baseline, with a mass and coupling matching prescription in the $\\overline{\\rm MS}$ scheme at 2 GeV. The vector-current correlator $V(t)$ is expanded as $V^{(0)}(t)+\\delta V(t)$, and $\\delta V(t)$ is built from self-energy, exchange, tadpole, pseudoscalar, and scalar insertion diagrams; the photon zero mode is removed with the QEDL prescription. The quenched-QED approximation treats sea quarks as electrically neutral, and quark-disconnected diagrams are neglected. The HVP integral uses the time-momentum representation, and for the light quark the ratio $\\delta a_\\mu^{\\rm HVP}(ud)/a_\\mu^{{\\rm HVP},(0)}(ud)$ is extrapolated with a fit that includes chiral logarithms or quadratic terms, an $a^2$ discretization term, and exponential or power-law finite-volume corrections; strange and charm use linear fits with $1/L^3$ terms. Non-perturbative QCD+QED renormalization constants for the vector current and quark masses are the input that improves the strange and charm errors by roughly a factor of three relative to the earlier calculation.","core_discovery":"At leading order in $\\alpha_{\\rm em}$ and $(m_d-m_u)/\\Lambda_{\\rm QCD}$, the paper’s central claim is that the total quark-connected isospin-breaking correction to the HVP contribution is $\\delta a_\\mu^{\\rm HVP}(udsc)=7.1\\,(2.9)\\times10^{-10}$, the most accurate determination available. The individual flavor contributions after extrapolation are $\\delta a_\\mu^{\\rm HVP}(ud)=7.1\\,(2.5)\\times10^{-10}$, $\\delta a_\\mu^{\\rm HVP}(s)=-0.0053\\,(33)\\times10^{-10}$, and $\\delta a_\\mu^{\\rm HVP}(c)=0.0182\\,(36)\\times10^{-10}$. The light-quark term splits into a QED part $1.1\\,(1.0)\\times10^{-10}$ and a strong-isospin part $6.0\\,(2.3)\\times10^{-10}$, so the correction is dominated by the up-down mass difference. Adding an estimated uncertainty for the quenched-QED and disconnected-diagram approximation, the authors quote $\\delta a_\\mu^{\\rm HVP}(udsc)=7.1\\,(2.6)\\,(1.2)\\times10^{-10}$ and combine it with leading-order connected and disconnected HVP results to obtain $a_\\mu^{\\rm HVP}=682\\,(19)\\times10^{-10}$, which agrees with dispersive determinations. The result agrees within errors with earlier dispersive and lattice estimates and is more precise than them.","pith_inferences":["The assigned qQED-plus-disconnected uncertainty is about one-third of the total error; if the disconnected QED light-quark diagram is as large as one lattice group’s estimate, that uncertainty is a floor, and a future unquenched-QED simulation could shift the central value by roughly $1\\times10^{-10}$ or more.","Because the paper extrapolates the ratio $\\delta a_\\mu^{\\rm HVP}(ud)/a_\\mu^{{\\rm HVP},(0)}(ud)$, its absolute IB correction is tied to the leading-order HVP normalization; another group with a more precise $a_\\mu^{{\\rm HVP},(0)}$ could rescale the ratio to obtain an updated absolute value immediately.","The same RM123 machinery could be applied to the hadronic light-by-light contribution to $g-2$, whose leading isospin-breaking corrections are so far unknown; a disconnected diagram issue there may be more severe than in HVP.","The published strange-quark value differs from a separate lattice result by about two standard deviations; a combined reanalysis of the two raw datasets would clarify whether this is a statistical fluctuation or a systematic difference."],"forward_implications":["At the precision of the new muon $g-2$ experiments, the isospin-breaking correction of about $7\\times10^{-10}$ is not negligible and must be included in a first-principles HVP prediction.","Because the light-quark correction is dominated by the strong up-down mass difference term $6.0\\,(2.3)\\times10^{-10}$, improving the input value of $m_d-m_u$ directly shrinks the total error.","The strange and charm corrections, both below $10^{-10}$, can be neglected at current precision without materially changing the result.","Adding the IB correction to leading-order HVP estimates gives $a_\\mu^{\\rm HVP}=682\\,(19)\\times10^{-10}$, compatible with dispersive determinations, so the IB corrections do not by themselves alter the muon $g-2$ discrepancy."],"supporting_citations":[{"why":"the companion full calculation that supplies the lattice ensembles, the extrapolation ansatz, and the detailed error budget.","marker":"[12]"},{"why":"introduces the RM123 expansion of the path integral in powers of the quark mass difference and $\\alpha_{\\rm em}$.","marker":"[13]"},{"why":"derives the RM123 insertion diagrams, including tadpole and pseudoscalar counterterms, used for the IB correlator.","marker":"[14]"},{"why":"defines the matching prescription that sets the separation between isosymmetric QCD and isospin-breaking corrections.","marker":"[16]"},{"why":"provides the input value of $m_d-m_u$ used for the strong-isospin contribution.","marker":"[17]"},{"why":"gives the earlier strange and charm IB results and the leading-order strange and charm HVP values used here.","marker":"[23]"},{"why":"supplies the non-perturbative QCD+QED renormalization constants that improve the strange and charm errors by about a factor of three.","marker":"[27]"},{"why":"provides the leading-order light HVP value and lattice finite-volume estimates used in the light-quark combination.","marker":"[28]"},{"why":"gives the dispersive-based estimate that the new result is compared against and that supports the disconnected contribution estimate.","marker":"[36]"},{"why":"provides the lattice result including one disconnected QED diagram, used as comparison and as the basis for the qQED-plus-disconnected uncertainty.","marker":"[38]"}],"fun_headline_variants":["Lattice QCD gives most precise isospin shift for muon g-2: 7.1(2.9)e-10","Isospin correction to muon g-2: lattice QCD yields 7.1(2.9)e-10","Muon g-2 isospin shift: 7.1(2.9)e-10 from lattice QCD","Lattice QCD sets isospin correction for muon g-2 at 7.1(2.9)e-10"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the uncertainty from the quenched-QED approximation together with the neglect of quark-disconnected diagrams is about $1.0\\times10^{-10}$, the size of the connected QED term; that estimate is not derived, so if either effect is larger the central value shifts beyond the quoted error.","fun_headline_variants_meta":{"raw":{"variants":["Lattice QCD gives most precise isospin shift for muon g-2: 7.1(2.9)e-10","Isospin correction to muon g-2: lattice QCD yields 7.1(2.9)e-10","Muon g-2 isospin shift: 7.1(2.9)e-10 from lattice QCD","Lattice QCD sets isospin correction for muon g-2 at 7.1(2.9)e-10"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000776,"raw_usage":{"total_tokens":3565,"prompt_tokens":1208,"completion_tokens":2357,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":824,"completion_tokens_details":{"reasoning_tokens":2227}},"tokens_in":824,"tokens_out":2357,"duration_ms":15589,"temperature":1.0,"reasoning_tokens":2227,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:03:35.843661+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Computing the quark-disconnected QED diagram for the light-quark contribution on the same gauge ensembles and comparing its size with the connected QED piece $1.1\\,(1.0)\\times10^{-10}$ would test the central systematic assumption; a disconnected contribution substantially larger than $1\\times10^{-10}$ would move $\\delta a_\\mu^{\\rm HVP}(udsc)$ outside the quoted total uncertainty of $2.9\\times10^{-10}$.","supporting_citations":[],"review_version":1}