⚛️ phenomenology

The inclusive $\bar B \to X_s \gamma$ decay rate with higher precision

This paper presents a high-precision Standard Model prediction for the inclusive $\bar B \to X_s \gamma$ decay branching ratio, $(3.54 \pm 0.14)\times 10^{-4}$ , achieved by incorporating complete $\mathcal{O}(\alpha_s)$ corrections and exact $\mathcal{O}(\alpha_s^2)$ charm-mass dependence, which aligns with experimental data and sets a 95% confidence level lower bound of 670 GeV on the charged Higgs boson mass in the Two-Higgs-Doublet Model II.

Original authors: M. Misiak, H. M. Asatrian, K. Brune, M. Czaja, M. Czakon, M. Fael, C. Greub, T. Huber, F. Lange, L. T. Moos, M. Niggetiedt, A. Rehman, K. Schönwald, M. Steinhauser

Published 2026-09-01

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: M. Misiak, H. M. Asatrian, K. Brune, M. Czaja, M. Czakon, M. Fael, C. Greub, T. Huber, F. Lange, L. T. Moos, 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 subatomic world, matter is built from a small family of particles called quarks. Among these, the bottom quark is one of the heaviest, and like all heavy things, it is unstable. It does not last long before it transforms into lighter particles. This transformation is not random; it follows strict rules set by the Standard Model, the best theory we have for how the universe works at its most fundamental level. One of the most interesting ways a bottom quark can change is by emitting a flash of light—a photon—and turning into a strange quark. This process is rare, occurring in only a few out of every hundred thousand decays, but it is a powerful tool for scientists. Because the process is so rare and so sensitive to the underlying laws of physics, even a tiny deviation from the predicted rate could signal the presence of new, undiscovered particles or forces lurking in the background.

A team of researchers has now performed a new, highly refined calculation of how often this specific transformation happens. They focused on the inclusive decay of the B meson, a particle containing a bottom quark, into a strange quark and a photon. By improving the precision of their theoretical prediction, they have tightened the net around the Standard Model, making it much harder for new physics to hide in the shadows. Their work confirms that the current experimental measurements match the theory almost perfectly, but it also sets a new, stricter limit on the mass of a hypothetical particle called the charged Higgs boson, which appears in many theories that extend beyond our current understanding.

The researchers began by tackling the messy details of the calculation. In the quantum world, particles do not just interact in simple pairs; they can briefly create clouds of other particles that pop in and out of existence. To predict the decay rate accurately, the team had to account for these complex interactions, including scenarios where the bottom quark splits into four or five particles before recombining. Previous calculations had to make educated guesses about the mass of the charm quark, a heavy particle that often appears in these virtual clouds. The team replaced these guesses with an exact calculation that accounts for the true mass of the charm quark. This was a massive computational task, requiring them to solve equations that describe four loops of particle interactions simultaneously. The result was a significant shift in the predicted rate, increasing it by about 4.2 percent compared to older estimates, and removing a large source of uncertainty that had plagued previous studies.

With the theoretical framework now more precise, the team combined their new numbers with the most up-to-date experimental data. They found that the predicted rate for the decay is (3.54 ± 0.14) × 10⁻⁴. This number represents the probability of the event occurring, and it aligns almost perfectly with the current experimental average of (3.49 ± 0.19) × 10⁻⁴. The agreement is so close that it suggests the Standard Model is holding up remarkably well under scrutiny. However, this agreement is not a dead end; it is a powerful constraint. Because the theory and experiment match so closely, any new physics that tries to alter this decay rate must be very subtle or involve very heavy particles that are difficult to produce.

The researchers used this tight agreement to test a popular extension of the Standard Model known as the two-Higgs-doublet model. This theory proposes that there is not just one Higgs boson, as we have observed, but a whole family of them, including a charged version. The team calculated how the presence of this charged Higgs boson would change the decay rate. They found that if this particle were too light, it would have pushed the decay rate higher than what is actually observed. By comparing their precise prediction with the experimental data, they determined that the charged Higgs boson must be heavier than 670 GeV at a 95 percent confidence level. This is a significant finding because it rules out a wide range of masses for this hypothetical particle, narrowing the search for new physics.

This work represents a milestone in precision physics. By moving from approximations to exact calculations and incorporating the latest experimental data, the team has sharpened our view of the subatomic world. They have shown that the Standard Model remains robust, but they have also drawn a clearer boundary around where new physics might still be hiding. The fact that the charged Higgs boson, if it exists, must be heavier than 670 GeV means that future experiments will need to reach even higher energies to find it. For now, the universe continues to behave exactly as the Standard Model predicts, but with a level of precision that leaves very little room for error.

Technical Summary: The Inclusive $\bar{B} \to X_s\gamma$ Decay Rate with Higher Precision

Problem and Motivation
The inclusive weak radiative decay of the $B$ meson, specifically the CP- and isospin-averaged branching ratio $\mathcal{B}_{s\gamma}$ for $\bar{B} \to X_s\gamma$ , serves as a critical probe for constraining physics beyond the Standard Model (SM). While the current experimental average is $(3.49 \pm 0.19) \times 10^{-4}$ (with a photon energy cut $E_\gamma > 1.6$ GeV), achieving a comparable theoretical precision requires a rigorous treatment of perturbative QCD corrections and non-perturbative effects. Previous theoretical analyses faced limitations, particularly regarding the treatment of the charm quark mass ( $m_c$ ) dependence in next-to-next-to-leading order (NNLO) corrections and the formal completeness of the leading-order (LO) perturbative series.

Methodology
The authors perform an upgraded analysis within the Standard Model using the low-energy effective theory (LEFT) framework. The calculation proceeds through the following methodological steps:

Perturbative Corrections:
- $\mathcal{O}(\alpha_s)$ Completeness: The analysis incorporates four- and five-body final state contributions arising from four-quark LEFT operators with $(\bar{s}b)(\bar{q}q)$ flavor content (where $q \in \{u, d, s\}$ ). This inclusion renders the $\mathcal{O}(\alpha_s)$ calculation formally complete.
- $\mathcal{O}(\alpha_s^2)$ Exactness: A primary methodological advancement is the exact calculation of the $\mathcal{O}(\alpha_s^2)$ corrections for the dominant interference terms $\hat{G}^{(2)}_{17}$ and $\hat{G}^{(2)}_{27}$ . Unlike previous works that interpolated between the limits $m_c=0$ and $m_c \gg m_b$ , this study utilizes four-loop propagator diagrams with unitarity cuts to determine the exact dependence on the charm quark mass.
- Renormalization: The calculation accounts for the renormalization group evolution of Wilson coefficients from the electroweak scale ( $\mu_0 \sim M_W, m_t$ ) down to the low-energy scale ( $\mu_b \sim m_b/2$ ).
Non-Perturbative Effects:
- Heavy Quark Expansion (HQE): Non-perturbative corrections are parameterized using HQE matrix elements. These parameters are adopted from the most recent kinetic-scheme fits of semileptonic $B$ -meson decays, which include $\mathcal{O}(\alpha_s^3)$ corrections and up-to-date experimental data.
- Resolved Photons and Collinear Logarithms: The analysis estimates non-perturbative effects arising from "resolved photons" and collinear logarithms. Constraints are derived from measured isospin asymmetry and Soft-Collinear Effective Theory (SCET) analyses. The uncertainty associated with these effects is carefully propagated, particularly for the $Q_7$ - $Q_7$ interference and the photon energy cut extrapolation.

Key Contributions

Exact $m_c$ Dependence: The paper provides the first exact calculation of the $\mathcal{O}(\alpha_s^2)$ corrections for $\hat{G}^{(2)}_{17}$ and $\hat{G}^{(2)}_{27}$ with arbitrary $m_c$ and $m_b$ . This removes a significant source of uncertainty previously addressed only by interpolation.
Formal Completeness: By including the previously missing four- and five-body final states, the $\mathcal{O}(\alpha_s)$ calculation is now formally complete.
Updated Parameters: The study utilizes the latest kinetic-scheme fits for the $b$ -quark mass and HQE matrix elements, ensuring consistency with modern semileptonic decay data.

Results
The authors report the following results for the Standard Model prediction of the branching ratio with a photon energy cut $E_\gamma > 1.6$ GeV:
$\mathcal{B}_{s\gamma}^{\text{SM}} = (3.54 \pm 0.14) \times 10^{-4}$
The total uncertainty of approximately $\pm 4\%$ is a combination of higher-order theoretical uncertainties ( $\pm 3\%$ ) and parametric uncertainties ( $\pm 2.7\%$ ).

Impact of Exact $m_c$ Calculation: The exact treatment of the charm mass leads to an enhancement of the branching ratio by approximately $4.2\%$ compared to previous interpolated results. This shift exceeds the previously estimated interpolation uncertainty by a factor of 1.4.
Scale Dependence: The inclusion of $\mathcal{O}(\alpha_s^2)$ corrections significantly reduces the dependence of the result on the renormalization scale $\mu_b$ . Varying $\mu_b$ in the range $[1.5, 10]$ GeV results in a variation of only $\pm 1.3\%$ .
Agreement with Experiment: The theoretical result is in perfect agreement with the experimental average of $(3.49 \pm 0.19) \times 10^{-4}$ .

Significance and Constraints on New Physics
The paper claims that the improved precision allows for robust constraints on Beyond Standard Model (BSM) physics, despite the lack of a direct signal for new physics in the current data.

Two-Higgs-Doublet Model (2HDM-II): In the context of the 2HDM-II, where contributions depend on $\tan\beta$ and the charged Higgs mass $M_{H^\pm}$ , the authors derive an absolute lower bound. For any value of $\tan\beta$ , the charged Higgs mass must satisfy:
$M_{H^\pm} > 670 \text{ GeV} \quad (95\% \text{ C.L.})$
The authors note that this bound is stronger than those available from direct LHC searches for the considered model.
Generic BSM Extensions: The paper provides a phenomenological formula relating deviations in $\mathcal{B}_{s\gamma}$ to modifications of the Wilson coefficients $C_7$ and $C_8$ at the electroweak scale:
$\mathcal{B}_{s\gamma} \times 10^4 \approx (3.54 \pm 0.14) - 8.51 \Delta C_7 - 2.16 \Delta C_8$
This linearized relation facilitates constraints on generic extensions of the SM, assuming new physics effects are not large enough to invalidate the perturbative expansion.

In conclusion, the paper establishes a new benchmark for the theoretical precision of inclusive radiative $B$ decays, removing a major source of theoretical uncertainty and tightening constraints on charged Higgs bosons in the 2HDM-II.

More like this