⚛️ lattice

Calibrated correlation between heavy-quark masses and Hadronic Vacuum Polarization observables at the precision frontier

This paper establishes a self-consistent framework using relativistic QCD Sum Rules to simultaneously determine heavy-quark masses and their Hadronic Vacuum Polarization contributions to the muon anomalous magnetic moment, leveraging their intrinsic anticorrelation to achieve unprecedented phenomenological precision and provide a direct diagnostic for residual theoretical uncertainties.

Original authors: Arnau Beltran, Pere Masjuan, Antonio Rivera

Published 2026-08-12

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Arnau Beltran, Pere Masjuan, Antonio Rivera

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to bake the perfect cake, but you don't have a recipe for the most important ingredient: the sugar. Instead, you have to guess how much sugar is in the cake by tasting the crumbs left over. In the world of particle physics, scientists are trying to predict the behavior of a tiny particle called the muon. They are incredibly good at this, but there is a tiny, stubborn mismatch between their predictions and what they actually measure in the lab. The biggest source of this mismatch is a "cloud" of virtual particles that pop in and out of existence around the muon. This cloud, known as Hadronic Vacuum Polarization (HVP), is like the sugar in our cake—it's hard to measure directly because it's made of messy, complex interactions.

For a long time, scientists have tried to calculate this "sugar" using two different methods. One method looks at real-world data from particle collisions (the "tasting crumbs" approach), and the other uses supercomputers to simulate the laws of physics from scratch (the "theoretical recipe" approach). The problem is that these two methods disagree with each other, and no one is sure why. It's like if your taste test said the cake needed two cups of sugar, but your recipe book said it needed three. To solve this mystery, scientists need to find a way to check their math in a place where the ingredients are cleaner and easier to understand. This is where heavy quarks—specifically the "charm" and "bottom" types—come in. They are like the heavy, dense chocolate chips in our cake: they are easier to track than the fluffy sugar, and they offer a cleaner way to test the rules of the game.

In this paper, a team of physicists proposes a clever new way to solve the "sugar" problem by realizing that the weight of the chocolate chips and the amount of sugar they contribute to the cake are actually two sides of the same coin. Traditionally, scientists would first measure the weight of the chocolate chips (the heavy-quark mass) and then use that number as a fixed input to calculate the sugar contribution. The authors argue that this is like measuring the weight of the chips, writing it down, and then forgetting that the chips and the sugar are baked from the exact same batter. They show that the mass of the heavy quark and its contribution to the muon's magnetic moment are not independent; they are both derived from the same underlying "spectral function," which is just a fancy way of describing the energy distribution of the particles involved.

The authors introduce a new mathematical tool called "Generalized Correlated Sum Rules." Think of this as a new kind of scale that can weigh the chocolate chips and the sugar contribution at the exact same time, rather than weighing them separately. By treating the mass and the magnetic moment contribution as a single, linked pair, they can use the fact that if one goes up, the other tends to go down (a relationship called anticorrelation) to sharpen their measurements. This simultaneous approach allows them to cancel out some of the guesswork and reduce the uncertainty in their results.

Using this method, the team calculated the heavy-quark contributions to the muon's magnetic moment with unprecedented precision. For the charm quark, they found a value of 14.46(13) × 10⁻¹⁰ at the leading order, and for the bottom quark, 0.3009(17) × 10⁻¹⁰. They also determined the masses of these quarks with high accuracy: the charm quark mass is 1267.1(6.8) MeV, and the bottom quark mass is 4182.3(7.2) MeV. When they compared their results with the latest computer simulations (lattice QCD), they found excellent agreement, especially for the charm sector. This suggests that their new "linked scale" method is working correctly and that the tension between different ways of calculating the muon's behavior might be resolved by better understanding these correlations.

The paper also points out that this method does more than just give a number; it acts as a diagnostic tool. Because the theory and the data are forced to agree on the mass and the magnetic moment simultaneously, any remaining disagreement at other levels of calculation reveals exactly where the models might be failing. It's like if your cake scale told you that the weight was perfect, but the texture was slightly off; the scale would then point directly to the flour, not the sugar, as the problem. The authors found that the remaining uncertainties are now dominated by the precision of the experimental data on the particles themselves, rather than the theoretical models used to describe them. This is a significant step forward, suggesting that future improvements will come from better measurements of the particles rather than just tweaking the math.

While this work focuses on the heavy quarks, the authors suggest that the strategy could be applied to the more difficult "light quark" sector, which is the main source of the current disagreement in the muon's magnetic moment. They don't claim to have solved the entire mystery of the muon yet, but they have built a powerful new tool that makes the path forward much clearer. By recognizing that the mass and the magnetic contribution are partners in a dance, rather than strangers passing in the night, they have tightened the grip on the numbers and provided a fresh perspective on one of physics' most persistent puzzles.

Technical Summary: Calibrated Correlation between Heavy-Quark Masses and Hadronic Vacuum Polarization Observables

Problem Statement
The theoretical prediction of the muon anomalous magnetic moment, $a_\mu$ , is currently limited by the uncertainty in the Hadronic Vacuum Polarization (HVP) contribution. A significant tension exists between dispersive determinations (based on $e^+e^- \to \text{hadrons}$ data) and lattice QCD calculations. While the light-quark sector is dominated by complex resonance structures and thresholds that make disentangling uncertainties difficult, the heavy-quark (charm and bottom) sectors offer a theoretically cleaner environment. In these sectors, both the timelike (dispersive) and spacelike (perturbative QCD) descriptions of the vector-current correlator are under good perturbative control. However, conventional approaches treat the determination of the heavy-quark mass ( $\hat{m}_q$ ) and its HVP contribution ( $a_\mu^{HVP, q}$ ) as independent problems: either the mass is determined first and used as an external input for $a_\mu$ , or $a_\mu$ is computed from a spectral model without reference to the mass. This paper argues that these are not independent problems, as both quantities are weighted integrals of the same hadronic spectral function, differing only in their integration kernels.

Methodology
The authors introduce a "generalized, kernel-weighted, correlated sum rule" framework within Relativistic QCD Sum Rules (RQSR) to simultaneously determine the $\overline{\text{MS}}$ heavy-quark mass ( $\hat{m}_q$ ) and its HVP contribution ( $a_\mu^{(i)}$ ).

Generalized Moments: The standard moment sum rules, defined as derivatives of the correlator $\Pi(q^2)$ at $q^2=0$ , are generalized to arbitrary admissible kernels $K(s)$ . The standard moments use $K(s)=1$ (weighting by $s^{-(n+1)}$ ), while the HVP contribution uses the physical QED kernel $\hat{K}^{(i)}(s)$ (weighting by $\hat{K}^{(i)}(s)/s^2$ ).
Simultaneous Extraction: By promoting the HVP kernel to a free choice within the sum rule formalism, the authors solve the self-consistency condition $M^{\text{th}}_n[K] = M^{\text{had}}_n[K]$ for a pair of moments (typically the zeroth moment $M_0$ and a higher moment $M_n$ ). This allows $\hat{m}_q$ and the continuum shape parameter $\lambda_q^3$ to be fixed simultaneously with the HVP observable.
Intrinsic Correlation: The framework explicitly exploits the anticorrelation between $\hat{m}_q$ and $a_\mu^{(i)}$ . Since both are derived from the same spectral function, fixing them together reduces the final uncertainty on $a_\mu^{(i)}$ compared to treating the mass as an independent external input.
Self-Calibration and Diagnostics: The method enforces agreement between the perturbative and hadronic descriptions of $a_\mu^{(i)}$ directly at the extraction point (where the kernel matches the observable). The residual spread between theory and hadronic sides at other moment pairs serves as a direct diagnostic for residual theory/model dependence, including duality violations and continuum modeling systematics.
Continuum Ansatz Extensions: The authors refine the continuum ansatz used in previous works:
- Charm: Explicitly include the $\psi(3770)$ and $\psi(4040)$ resonances as narrow states rather than absorbing them into the continuum.
- Bottom: Extend mass corrections in the continuum ansatz from $O(\hat{m}_q^2/s)$ to $O(\hat{m}_q^4/s^2)$ , including associated $O(\alpha_s^3)$ corrections.

Key Contributions

Unified Framework: A novel sum-rule construction that treats the heavy-quark mass and its HVP contribution as correlated observables derived from a single spectral function, rather than sequential steps.
Uncertainty Reduction: By accounting for the intrinsic anticorrelation, the method achieves a significant reduction in the uncertainty of the HVP contribution compared to non-correlated, data-driven determinations that do not account for quark mass or hadronic correlations.
Diagnostic Tool: The residual spread between theory and hadronic evaluations at non-calibration moment pairs provides a specific, observable-dependent measure of duality violation and modeling systematics, which is inaccessible in standard mass-only determinations.
NLO Extension: The framework is extended to Next-to-Leading Order (NLO) kernels, demonstrating applicability to higher-order corrections, with specific handling required for the divergent logarithmic term in the $\hat{K}^{(4a)}$ kernel.

Results
The authors present simultaneous determinations for the charm and bottom sectors at Leading Order (LO) and Next-to-Leading Order (NLO).

Heavy-Quark Masses:
- Charm: $\hat{m}_c(\hat{m}_c) = 1267.1(6.8)$ MeV.
- Bottom: $\hat{m}_b(\hat{m}_b) = 4182.3(7.2)$ MeV.
- These results exhibit unprecedented phenomenological precision.
HVP Contributions (LO):
- Charm: $a_\mu^{c, \text{LO}} = (14.46 \pm 0.13) \times 10^{-10}$ .
- Bottom: $a_\mu^{b, \text{LO}} = (0.3009 \pm 0.0017) \times 10^{-10}$ .
- These values represent an uncertainty reduction of approximately 45% (charm) and 90% (bottom) compared to previous non-correlated, data-driven estimates.
HVP Contributions (NLO):
- Charm: $a_\mu^{c, \text{NLOa}} = (-0.8006 \pm 0.0093) \times 10^{-10}$ ; $a_\mu^{c, \text{NLOb}} = (0.2268 \pm 0.0019) \times 10^{-10}$ .
- Bottom: $a_\mu^{b, \text{NLOa}} = (-0.2295 \pm 0.0013) \times 10^{-11}$ ; $a_\mu^{b, \text{NLOb}} = (0.04730 \pm 0.00027) \times 10^{-11}$ .
- The NLO charm results show good agreement with the first available lattice QCD determination at this order.
Calibration: The continuum shape parameter $\lambda_q^3$ extracted from the sum rules is calibrated against experimental data for the spectral function in the threshold region. The authors find good agreement, validating the continuum ansatz and the treatment of duality.

Significance and Claims
The paper claims that the heavy-quark sector provides a "theoretically clean environment" to quantitatively test sources of uncertainty in HVP determinations. The primary significance lies in the intrinsic construction which:

Sharpens Precision: Exploits the anticorrelation between mass and HVP to reduce final uncertainties without introducing new external inputs.
Diagnoses Systematics: Turns the residual disagreement between perturbative and hadronic descriptions into a direct, observable-specific diagnostic for duality violations and continuum modeling, a feature unavailable in mass-only analyses.
Proof of Principle: Establishes a general strategy applicable to any observable expressible as a kernel-weighted dispersive integral. The authors posit that this approach could eventually be applied to the broader, more difficult light-quark HVP problem to address the current tension between data-driven and lattice determinations, provided the analytic properties of the relevant kernels allow for the necessary subtractions.

The work does not claim to resolve the light-quark HVP tension immediately but offers a refined methodology and a controlled proof of principle for extracting correlated hadronic quantities from a unified sum-rule framework.

Technical Summary: Calibrated Correlation between Heavy-Quark Masses and Hadronic Vacuum Polarization Observables

More like this