Hadronic vacuum polarization contribution to $a_\mu$ from functional methods with strong and electromagnetic isospin breaking

Adnan Bashir; Angel S. Miramontes; Christian S. Fischer; Pablo Roig

Continuum DSE/BSE QCD with pion back-reaction, a dynamical ρ, and self-consistent isospin breaking gives a_μ^{HVP,LO} = (710.0 ± 14.5)×10⁻¹⁰, matching lattice results.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-15 02:19 UTC pith:U3OMKC2Q

load-bearing objection Continuum DSE/BSE HVP with self-consistent ISB lands near lattice numbers, but the truncation controlling the intermediate-Q^{2} region remains the open question and we only have the abstract. the 3 major comments →

arxiv 2607.12918 v1 pith:U3OMKC2Q submitted 2026-07-14 hep-ph hep-latnucl-th

Hadronic vacuum polarization contribution to a_μ from functional methods with strong and electromagnetic isospin breaking

Angel S. Miramontes , Adnan Bashir , Christian S. Fischer , Pablo Roig This is my paper

classification hep-ph hep-latnucl-th

keywords hadronic vacuum polarizationmuon anomalous magnetic momentDyson-Schwinger equationsBethe-Salpeter equationsisospin breakingquark-photon vertexcontinuum QCDrho resonance

verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper sets out to compute the leading-order hadronic vacuum-polarization contribution to the muon's anomalous magnetic moment from continuum quantum chromodynamics, not lattice simulations. Working in the Dyson–Schwinger and Bethe–Salpeter framework, the authors include pion back-reaction, a fully dressed quark–photon vertex that dynamically generates the ρ resonance in the timelike region, and both strong and electromagnetic isospin breaking treated self-consistently at the quark level. They report a_μ^{HVP,LO}(u+d+s+c) with isospin breaking equal to 709.7×10⁻¹⁰ and a final result including the bottom quark of (710.0 ± 14.5)×10⁻¹⁰, with an isospin-breaking shift of only 4.5×10⁻¹⁰ (0.6%). A sympathetic reader cares because hadronic vacuum polarization dominates the theory uncertainty on the muon magnetic moment, a precision probe of the Standard Model; an independent continuum determination that agrees with lattice QCD and quantifies isospin breaking strengthens that comparison and shows that the shift, though modest, is not negligible at the target precision.

Core claim

With pion back-reaction, a fully dressed quark–photon vertex containing a dynamically generated ρ resonance, and self-consistent strong plus electromagnetic isospin breaking, continuum Dyson–Schwinger/Bethe–Salpeter QCD yields a_μ^{HVP,LO}(u+d+s+c)|_ISB = 709.7×10⁻¹⁰ and a final a_μ^{HVP,LO}(u+d+s+c+b)|_ISB = (710.0 ± 14.5)×10⁻¹⁰, with Δa_μ^{HVP,LO} = 4.5×10⁻¹⁰ (0.6%) from isospin breaking, in good agreement with recent lattice-QCD determinations.

What carries the argument

The load-bearing mechanism is the continuum Dyson–Schwinger and Bethe–Salpeter system built around a fully dressed quark–photon vertex that produces a dynamical ρ resonance in the timelike region, together with pion back-reaction and self-consistent strong and electromagnetic isospin breaking at the quark level; this generates the hadronic vacuum polarization without lattice discretization.

Load-bearing premise

The chosen truncation of the Dyson–Schwinger and Bethe–Salpeter equations—the ansatz for the quark–gluon interaction and the construction of the dressed quark–photon vertex—captures the nonperturbative dynamics of hadronic vacuum polarization at the few-percent precision claimed.

What would settle it

A controlled lattice-QCD determination of a_μ^{HVP,LO}(u+d+s+c+b) with isospin breaking that lies outside the (710.0 ± 14.5)×10⁻¹⁰ window, or a continuum recalculation with a substantially different quark–gluon interaction that moves the central value by more than the quoted systematic uncertainty.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

Continuum functional methods can reach lattice-level agreement on a_μ^{HVP,LO} once resonance structure and isospin breaking are included.
Isospin-breaking corrections of order 0.6% should be retained rather than dropped in high-precision HVP determinations.
The bottom-quark piece is small; the u+d+s+c result already carries the bulk of the leading-order HVP.
Independent continuum and lattice HVP determinations now agree at the few-percent level, tightening the theory side of the muon g−2 comparison.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If the truncation holds, the same DSE/BSE setup with controlled isospin breaking could be extended to higher-order HVP or light-by-light scattering.
A 0.6% ISB shift alone is unlikely to close any residual theory–experiment tension in a_μ; the bulk of that tension would have to sit elsewhere.
Varying the quark–gluon interaction ansatz while freezing the vertex and ISB treatment would map the systematic band more tightly than the quoted ±14.5.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

Continuum DSE/BSE HVP with self-consistent ISB lands near lattice numbers, but the truncation controlling the intermediate-Q^{2} region remains the open question and we only have the abstract.

read the letter

The one thing you need to know is that they report a_μ^{HVP,LO}(u+d+s+c+b)|_ISB = (710.0 ± 14.5)×10^{-10} from a continuum DSE/BSE calculation that includes pion back-reaction, a fully dressed quark-photon vertex with a dynamical ρ, and self-consistent strong plus electromagnetic isospin breaking at the quark level. The ISB shift itself is 4.5×10^{-10} (0.6 %). That number sits close to recent lattice results and is the piece that matters for the g-2 theory side.

What is actually new is the simultaneous packaging of those three ingredients inside one continuum framework. Earlier DSE/BSE HVP calculations existed and lattice groups already have ISB, but the combination with a dynamically generated ρ in the vertex and quark-level ISB treated self-consistently is a genuine step for functional methods. They produce a concrete, comparable number rather than a qualitative statement, and they correctly flag that the ISB piece, while modest, is not negligible. That is real progress and deserves credit.

The soft spot is exactly the one the stress-test flags and it is load-bearing: the result lives or dies by the quark-gluon interaction ansatz and the construction of the dressed vertex. Those choices dominate the intermediate-Q^{2} region that drives the HVP integral. With only the abstract we have no truncation tests, no free-parameter list, and only an “indicative” systematic. The ±14.5 band cannot be audited yet. Circularity looks moderate—they appear to constrain with meson observables rather than a_μ itself—but that cannot be verified here either. No code or data artifacts are visible.

This paper is for the g-2 theory community and for people who push continuum functional methods toward precision observables. A lattice reader will want the systematics section; a DSE reader will want the vertex details. The work is coherent on its own terms and engages the right literature, so it is serious. I would send it to peer review. A competent referee can demand the missing error budget and truncation studies; the central claim is important enough to warrant that time rather than a desk reject.

Referee Report

3 major / 1 minor

Summary. The manuscript reports a continuum-QCD evaluation of the leading-order hadronic vacuum-polarization contribution to the muon anomalous magnetic moment in the Dyson–Schwinger/Bethe–Salpeter framework. The calculation includes pion back-reaction, a fully dressed quark–photon vertex that dynamically generates a ρ-resonance structure in the timelike region, and self-consistent strong plus electromagnetic isospin breaking at the quark level. The central results quoted are a_μ^{HVP,LO}(u+d+s+c)|_ISB = 709.7 × 10^{-10}, an isospin-breaking shift Δa_μ^{HVP,LO} = 4.5 × 10^{-10} (0.6 %), and a final value including the bottom quark of (710.0 ± 14.5) × 10^{-10}, stated to be in good agreement with recent lattice-QCD determinations.

Significance. If the numerical result and its error budget survive full scrutiny of the underlying truncation and systematics, the work would supply an independent continuum functional determination of a_μ^{HVP,LO} that incorporates a dynamical ρ in the quark–photon vertex and a self-consistent treatment of both strong and electromagnetic isospin breaking. The explicit 0.6 % ISB shift is of direct phenomenological interest for the muon g-2 program and would complement lattice evaluations. The abstract alone, however, does not yet establish that the claimed few-percent precision is under control.

major comments (3)

The central numerical claim (a_μ^{HVP,LO}(u+d+s+c)|_ISB = 709.7 × 10^{-10} and the final (710.0 ± 14.5) × 10^{-10}) rests on a specific DSE/BSE truncation: the ansatz for the quark–gluon interaction and the construction of the fully dressed quark–photon vertex that produces the dynamical ρ. These choices control the intermediate-Q^{2} region that dominates the HVP integral. The abstract asserts that the ingredients are present and that the result agrees with lattice QCD, but supplies no independent verification (e.g., truncation-variation tests or comparison to known intermediate observables) that the truncation error lies inside the quoted ±14.5 band. This premise is load-bearing for any few-percent claim.
The final error ±14.5 is described only as “an indicative estimate of the systematic uncertainties.” No breakdown of the error budget (interaction-parameter variation, vertex truncation, continuum extrapolation, missing higher-order effects, etc.) is given in the abstract. Without an auditable derivation of this band, the quantitative agreement with lattice determinations and the significance of the 4.5 × 10^{-10} ISB shift cannot be assessed.
Only the abstract is available for review. Intermediate results (quark propagators, vertex form factors, the HVP integrand itself, and any sensitivity studies) that would allow a referee to test the load-bearing modeling choices are therefore inaccessible. A full assessment of soundness is not possible on the present material.

minor comments (1)

The abstract is clearly written and the numerical claims are stated unambiguously; no presentation issues can be identified from the abstract alone.

Circularity Check

0 steps flagged

No circularity identifiable from the abstract alone; HVP result is a continuum DSE/BSE calculation benchmarked externally against lattice QCD.

full rationale

Only the abstract is available; no equations, parameter tables, or derivation chain can be inspected. The abstract reports a first-principles-style continuum calculation of a_μ^{HVP,LO} that includes pion back-reaction, a dressed quark–photon vertex with dynamical ρ structure, and self-consistent strong+EM isospin breaking, then quotes a numerical value (709.7×10^{-10}, final 710.0±14.5×10^{-10}) that is compared to independent lattice-QCD determinations. Nothing in the abstract text shows that the HVP integral or the isospin-breaking shift is obtained by fitting to a_μ itself, by renaming a known empirical pattern, or by a load-bearing self-citation that forces the result by construction. Standard practice in this framework is to constrain the quark–gluon interaction and vertex ansätze by other hadronic observables (masses, decay constants, etc.) rather than by the muon anomaly; that modeling choice is a truncation/systematic issue, not circularity under the stated criteria. With no quotable reduction of output to input, the honest finding is score 0 and empty steps.

Axiom & Free-Parameter Ledger

1 free parameters · 4 axioms · 0 invented entities

Abstract-only review: free parameters of the DSE interaction and vertex ansätze are not listed. The calculation rests on standard continuum QCD functional equations plus domain assumptions about truncations that close the infinite tower of Green’s functions. No new particles or forces are invented; the ρ is dynamically generated, not postulated as an elementary field. Counts below are the minimal set visible from the abstract.

free parameters (1)

DSE/BSE interaction and vertex truncation parameters (unspecified)
Continuum functional HVP calculations almost always introduce a finite set of scales/strengths in the effective quark–gluon interaction and in the quark–photon vertex construction. Values are not given in the abstract; they are free parameters relative to the HVP claim until fixed by independent data.

axioms (4)

domain assumption Truncated Dyson–Schwinger and Bethe–Salpeter equations with pion back-reaction adequately represent continuum QCD for LO HVP at the few-percent level.
Invoked as the computational framework of the whole paper; not proved from first principles in the abstract.
domain assumption A fully dressed quark–photon vertex that dynamically generates a ρ-resonance structure in the timelike region is sufficient for the HVP integral.
Stated as an incorporated ingredient; the completeness of that vertex construction is assumed.
domain assumption Strong and electromagnetic isospin breaking can be treated self-consistently at the quark level within the same truncation.
Central methodological claim of the abstract; validity depends on how mass and charge differences are inserted into the kernels.
standard math Standard continuum QCD and QED (quark masses, charges, electromagnetic coupling) as external inputs.
Background field theory assumed throughout.

pith-pipeline@v1.1.0-grok45 · 6159 in / 2910 out tokens · 28766 ms · 2026-07-15T02:19:49.558975+00:00 · methodology

0 comments

read the original abstract

We present a continuum-QCD determination of the leading-order hadronic vacuum-polarization contribution to the anomalous magnetic moment of the muon within the Dyson-Schwinger and Bethe-Salpeter equation framework. The calculation incorporates pion back-reaction, a fully dressed quark-photon vertex with a dynamically generated $\rho$-resonance structure in the timelike region, and both strong and electromagnetic isospin breaking treated self-consistently at the quark level. We obtain $a_\mu^{\mathrm{HVP,LO}}(u+d+s+c)|_{\mathrm{ISB}} = 709.7 \times 10^{-10}$, in good agreement with recent lattice-QCD determinations, and find an isospin-breaking shift of $\Delta a_\mu^{\mathrm{HVP,LO}} = 4.5 \times 10^{-10}$ ($0.6\%$), demonstrating that isospin-breaking effects, while quantitatively modest, are not negligible. Including the bottom-quark contribution and an indicative estimate of the systematic uncertainties, we obtain our final result, $a_\mu^{\mathrm{HVP,LO}}(u+d+s+c+b)|_{\mathrm{ISB}} = (710.0 \pm 14.5) \times 10^{-10}$.