⚛️ phenomenology

NNLO QCD corrections to the weak radiative $B$ -meson decay with exact dependence on $m_c$

This paper presents the first complete Next-to-Next-to-Leading-Order QCD calculation for the inclusive weak radiative $B$ -meson decay with exact dependence on the physical charm quark mass, yielding a Standard Model branching ratio of $(3.54 \pm 0.14)\times 10^{-4}$ that is in excellent agreement with experimental measurements.

Original authors: M. Czaja, M. Czakon, T. Huber, M. Misiak, M. Niggetiedt, A. Rehman, K. Schönwald, M. Steinhauser

Published 2026-09-01

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: M. Czaja, M. Czakon, T. Huber, M. Misiak, M. Niggetiedt, A. Rehman, K. Schönwald, M. Steinhauser

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

In the vast landscape of modern physics, the Standard Model serves as our most reliable map, charting how the fundamental particles of the universe interact. Yet, physicists suspect this map is incomplete, likely just a low-energy sketch of a deeper, more complex reality. To find the hidden terrain beyond this map, researchers do not always need to build larger machines to smash particles together; sometimes, they need only to look more closely at the subtle, rare behaviors of heavy particles that already exist. One such particle is the B-meson, a short-lived object containing a heavy bottom quark. Occasionally, this particle decays by emitting a photon, a particle of light, in a process that is forbidden by the simplest rules of physics but allowed by the intricate machinery of the Standard Model. Because this decay is so rare and sensitive to the underlying laws of nature, even a tiny deviation between the predicted rate and the observed rate could signal the presence of new, undiscovered particles or forces.

For decades, scientists have been refining their predictions for how often this specific decay should occur. The challenge lies in the fact that the calculation involves the interplay of many forces and particles, including the charm quark, which is significantly lighter than the bottom quark but heavy enough to complicate the math. Previous attempts to predict the decay rate had to rely on a mathematical shortcut: they calculated the result for two extreme scenarios—one where the charm quark was massless and another where it was infinitely heavy—and then guessed the answer for the real, physical mass of the charm quark by interpolating between those two extremes. This guesswork introduced a significant uncertainty, leaving a gap in our understanding that could have hidden new physics or obscured the true nature of the known laws.

A team of theoretical physicists has now closed that gap by performing a calculation that removes the need for any guessing. Instead of relying on an interpolation between extremes, they computed the decay rate using the exact, physical mass of the charm quark. This was a monumental task, requiring the evaluation of hundreds of thousands of complex diagrams that represent the quantum interactions occurring during the decay. These diagrams, which involve four loops of virtual particles, had to be solved with a precision that accounts for the specific mass of the charm quark relative to the bottom quark. The researchers used advanced mathematical techniques to reduce these massive calculations into a manageable set of core integrals, which they then solved with extreme numerical precision.

The result of this rigorous effort is a new, highly precise prediction for the branching ratio of the B-meson decay. The team found that the decay occurs with a frequency of approximately 3.54 in every 10,000 instances, with a very small margin of error. This new theoretical value aligns remarkably well with the current experimental average measured by particle detectors around the world, which stands at 3.49 with a slightly larger uncertainty. The agreement between the new, exact calculation and the experimental data is a strong confirmation that the Standard Model holds up under this specific, high-precision test.

Crucially, the new calculation reveals that the previous method of interpolation had slightly underestimated the decay rate. The difference, while small in absolute terms, was significant enough to shift the theoretical prediction by about four percent. This shift brings the theory even closer to the experimental measurements, reducing the room for potential discrepancies that might have hinted at new physics. The researchers also noted that the uncertainties in their prediction are now dominated by the precision of the input parameters and the experimental measurements, rather than by the theoretical approximations used in the past.

This work does not discover new particles, but it solidifies the foundation upon which the search for them rests. By removing the interpolation uncertainty, the team has provided a sharper tool for future investigations. If new physics exists in this sector, it must now be even more subtle than previously thought, as the Standard Model's prediction has been refined to a level that matches the experimental data with striking accuracy. The study stands as a testament to the power of precise theoretical calculation, showing that even without building a new collider, we can push the boundaries of our knowledge by solving the complex equations that govern the subatomic world.

Technical Summary: NNLO QCD Corrections to the Weak Radiative B-meson Decay with Exact Dependence on $m_c$

Problem Statement
The inclusive weak radiative decay of the $B$ meson ( $B \to X_s \gamma$ ) serves as a critical constraint on Beyond-the-Standard-Model (BSM) physics. Precise Standard Model (SM) predictions for its branching ratio are essential to match current experimental precision. While the Next-to-Next-to-Leading Order (NNLO) QCD corrections to the Wilson coefficients were established previously, the calculation of the corresponding interference matrices $\hat{G}_{ij}$ (specifically $\hat{G}_{17}$ and $\hat{G}_{27}$ ) remained incomplete. Previous evaluations relied on interpolating between the limits of zero charm quark mass ( $m_c = 0$ ) and infinite charm mass ( $m_c \gg m_b$ ). This interpolation introduced a theoretical uncertainty estimated at $\pm 3\%$ , which was a significant component of the total error budget. The primary challenge addressed in this work is the determination of these NNLO corrections for the physical value of the charm quark mass, $m_c$ , without resorting to interpolation.

Methodology
The authors calculate the NNLO QCD contributions to the matrices $\hat{G}_{17}$ and $\hat{G}_{27}$ , which arise from the interference of the electromagnetic dipole operator $Q_7$ with the four-quark current-current operators $Q_1$ and $Q_2$ . The calculation involves the following technical steps:

Diagram Generation and Setup: The computation requires evaluating hundreds of four-loop propagator diagrams with unitarity cuts and two mass scales ( $m_b$ and $m_c$ ). The authors utilize FeynArts and QGRAF for diagram generation. The calculation is performed in the Feynman-'t Hooft gauge ( $\xi=1$ ) using Naive Dimensional Regularization (NDR) with fully anticommuting $\gamma_5$ .
Integration and Reduction: The unrenormalized amplitudes are expressed as linear combinations of approximately $5 \times 10^5$ scalar integrals. Using the "reverse unitarity" method, cut propagators are treated to allow the application of Integration-By-Parts (IBP) identities. The code Kira 2.0 is employed to reduce these integrals to a set of roughly 500 Master Integrals (MIs). This reduction phase was computationally intensive, requiring weeks of CPU time and terabytes of RAM for certain families.
Evaluation of Master Integrals: The MIs are evaluated using two complementary approaches:
- Differential Equations (DEs): A system of first-order linear DEs with respect to the mass ratio $z = m_c^2/m_b^2$ is derived. Initial conditions are determined via asymptotic expansions in the limit $z \gg 1$ (where the charm loop contracts to a point, reducing the problem to three-loop integrals). The DEs are then solved numerically along complex contours to reach the physical value of $z \approx 0.05$ , avoiding branch points at $z=1/4$ .
- Expand and Match: As a cross-check, the "expand and match" technique is used to construct semi-analytic results around specific values of $z$ .
- AMFlow: The numerical tool AMFlow is utilized to compute boundary values and cross-check results, particularly for two- and four-particle cuts.
Renormalization and Combination: After UV renormalization, the contributions from two-, three-, and four-particle cuts are summed. The cancellation of infrared and ultraviolet poles is verified to an accuracy better than $10^{-55}$ . The authors also address specific counterterm contributions involving charm loops that were previously skipped in the literature.

Key Contributions

Exact $m_c$ Dependence: The paper provides the first determination of the NNLO corrections to $\hat{G}_{17}$ and $\hat{G}_{27}$ for the physical charm mass, eliminating the need for the $m_c=0$ to $m_c \gg m_b$ interpolation.
Analytical and Numerical Results: The authors present numerical fits for the functions $F_1(z, 1)$ and $F_2(z, 1)$ (which parameterize the NNLO contributions) across the physical region of $z$ . They also provide power-log expansions for small $z$ and large- $z$ asymptotic expansions.
Resolution of Interpolation Uncertainty: By calculating the exact dependence, the paper removes the $\pm 3\%$ uncertainty associated with the previous interpolation method.
Updated SM Prediction: The authors update the global SM prediction for the branching ratio $\mathcal{B}(B \to X_s \gamma)$ , incorporating the new NNLO results, updated input parameters (including recent HQE fits), and a revised treatment of non-perturbative effects.

Results
The final SM prediction for the CP- and isospin-averaged branching ratio with a photon energy cut $E_\gamma > 1.6$ GeV is:
$\mathcal{B}^{SM}_{s\gamma} = (3.54 \pm 0.14) \times 10^{-4}$
The total uncertainty of approximately $\pm 4\%$ is derived by combining parametric uncertainties ( $\pm 2.7\%$ ) and higher-order theoretical uncertainties ( $\pm 3\%$ ).

Key observations regarding the results include:

The new central value is approximately $1.7\%$ lower than the previous NLO result but significantly higher than intermediate NNLO estimates that relied on interpolation.
The removal of the interpolation uncertainty improves the precision of the theoretical prediction.
The QCD perturbation series exhibits good behavior, particularly when the renormalization scales $\mu_b$ and $\mu_c$ are set to $\frac{1}{2}m_b$ .
The result is in good agreement with the current Particle Data Group (PDG) experimental average of $(3.49 \pm 0.19) \times 10^{-4}$ .

Significance and Claims
The authors claim that this work represents a significant improvement in the theoretical precision of the $B \to X_s \gamma$ branching ratio. By replacing the interpolation with an exact calculation, the theoretical uncertainty is reduced, tightening the constraints on BSM physics. Specifically, in the context of the Two-Higgs-Doublet Model II, the updated prediction yields a 95% confidence level lower bound on the charged scalar mass of $670$ GeV.

The paper emphasizes that while the numerically dominant NNLO corrections are now included, the calculation is not yet "complete" in all aspects. Remaining approximations include:

The assumption of no restriction on photon energy for the calculated NNLO terms (though the physical cut is applied in the final result).
The treatment of three- and four-body final state interferences involving $Q_8$ only in the BLM approximation.
The neglect of NNLO corrections from the QCD penguin operators $Q_3$ – $Q_6$ .

The authors note that the smallness of the overall NNLO effect at the central scale is due to cancellations between BLM and non-BLM terms, highlighting the necessity of the exact calculation to avoid misleading approximations. The work establishes a new baseline for phenomenological analyses of flavor-changing neutral currents.

More like this