⚛️ phenomenology

Multi-parton contributions to $\bar B \to X_s γ$ at NLO

This paper presents the first complete calculation of the remaining formally NLO multi-parton contributions to the inclusive radiative decay $\bar{B} \to X_s \gamma$ , providing fully analytic results in terms of multiple polylogarithms and demonstrating that their numerical impact on the decay rate is small due to a partial cancellation with leading-order terms.

Original authors: Kevin Brune, Tobias Huber, Lars-Thorben Moos

Published 2026-01-30

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Kevin Brune, Tobias Huber, Lars-Thorben Moos

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 the universe as a giant, incredibly complex machine. For decades, physicists have been trying to build a perfect blueprint of how this machine works, called the Standard Model. One of the most sensitive parts of this machine is a specific type of particle decay: a heavy "bottom" quark turning into a "strange" quark while shooting out a flash of light (a photon). This process is like a rare, high-stakes magic trick that happens inside a particle accelerator.

For a long time, scientists have been measuring how often this trick happens. They are getting very precise measurements. However, to know if the machine is working exactly as the blueprint predicts—or if there's a hidden "new physics" glitch—they need an equally precise theoretical prediction.

The Puzzle of the Missing Pieces

Think of the theoretical prediction as a massive jigsaw puzzle. Over the last few decades, scientists have successfully placed most of the pieces. They have calculated the main events (the "two-body" decay, where the bottom quark just turns into a strange quark and a photon) with incredible detail, even accounting for tiny quantum fluctuations up to the third level of complexity (NNLO).

But, there was a small, stubborn gap in the puzzle.

While the main event is simple, nature sometimes gets messy. Occasionally, the bottom quark doesn't just produce a photon; it also accidentally creates extra particles, like a pair of light quarks (a "multi-parton" state) or even a gluon (the particle that holds quarks together). These are like the "extra crumbs" that fall off the table during the magic trick.

Previously, scientists had calculated the "four-body" version of these messy events (bottom → strange + 2 light quarks + photon) and the "five-body" version (adding a gluon). However, they were missing the Next-to-Leading Order (NLO) corrections for the four-body events.

Think of it this way: You have calculated the cost of a meal (the basic decay). You have also calculated the cost of the meal plus a side dish (the extra particles). But you hadn't yet calculated the "tip" or the "service charge" (the quantum corrections) specifically for the meal with the side dish. Without this tip, the total bill wasn't mathematically complete, even if the missing amount was small.

What This Paper Does

This paper by Kevin Brune, Tobias Huber, and Lars-Thorben Moos is the act of calculating that missing "tip." They computed the final, missing mathematical pieces required to make the theoretical prediction for this decay formally complete at the NLO level.

Here is how they tackled the challenge, using some creative analogies:

1. The "Reading Point" Rule (Handling the $\gamma_5$ Problem)
In the math of particle physics, there is a tricky object called $\gamma_5$ (gamma-five). It's like a special compass that only works in a 4-dimensional world. When scientists try to do calculations in the "fuzzy" mathematical space used for quantum mechanics (which has a tiny bit of extra dimension, $D=4-2\epsilon$ ), this compass starts spinning wildly.
The authors used a specific set of rules (the "KKS scheme") to handle this. Imagine you are reading a book, but the pages are slightly transparent. To make sure you don't get confused, they decided to always start reading from the same specific page (the "reading point") and never flip the book around. This keeps the math consistent, even if it feels a bit rigid.

2. The "Reverse Unitarity" Trick
The paper involves calculating the probability of particles flying out in all directions. This is like trying to count every possible way a handful of marbles can scatter after hitting a wall.
Usually, this is done by integrating over "phase space" (all possible angles and speeds). The authors used a clever trick called "reverse unitarity." Imagine taking a movie of the particles scattering and playing it backwards. By doing this, they could turn the messy problem of "scattering particles" into a cleaner problem of "particles moving in loops" (which is a type of math problem they know how to solve very well). This allowed them to reduce thousands of complex equations down to a manageable list of about 60 "master integrals" (the fundamental building blocks of the answer).

3. The "Collinear Logarithms"
When a photon is shot out, it sometimes flies almost perfectly parallel to a light quark, like two cars driving bumper-to-bumper at high speed. In the math, this creates a "singularity" (a number that tries to become infinite).
To fix this, the authors pretend the light quarks have a tiny bit of mass (like adding a tiny bit of weight to a feather). This stops the infinity. However, this introduces a new term in the equation called a "collinear logarithm." It's like a penalty fee that depends on how light the quark is. The authors calculated exactly how big this fee is.

The Result: A Small but Necessary Correction

After all this heavy lifting, what did they find?

The Size of the Effect: The missing pieces they calculated turned out to be very small. The numerical impact on the total decay rate is less than 1% (specifically, in the "per mille" range, or parts per thousand).
Why so small? There was a partial cancellation. The "LO" (Leading Order) contribution and the "NLO" (Next-to-Leading Order) contribution were fighting each other, canceling out much of the effect. It's like two people pushing a heavy box in opposite directions; the box barely moves.
The Importance: Even though the number is small, the calculation is crucial. In high-precision physics, you cannot have a "rough draft" of a prediction. If you want to claim that an experiment has discovered "New Physics" (something outside the Standard Model), you must be 100% sure that your theoretical prediction is complete. This paper provides the final piece of the puzzle, ensuring that the theoretical prediction is as sharp as the experimental measurements.

Summary

In short, this paper is the final polish on a very high-precision calculation. The authors didn't discover a new particle or a new force. Instead, they did the unglamorous but essential work of filling in the last few missing numbers in the theoretical blueprint. They used advanced mathematical tricks to handle complex quantum rules and confirmed that, while these specific multi-particle corrections are tiny, they are now fully accounted for. This allows physicists to look at experimental data with confidence, knowing that any remaining discrepancy is likely due to new physics, not a missing math term.

Technical Summary: Multi-parton contributions to $\bar{B} \to X_s\gamma$ at NLO

Problem Statement
The inclusive radiative decay $\bar{B} \to X_s\gamma$ serves as a precision test of the Standard Model (SM) and a sensitive probe for New Physics. While the theoretical prediction for this decay has been pushed to Next-to-Next-to-Leading Order (NNLO) in QCD for the dominant two-body partonic process ( $b \to s\gamma$ ), the Next-to-Leading Order (NLO) corrections for multi-parton final states remain incomplete. Specifically, the NLO contributions involving four-body final states ( $b \to s q \bar{q} \gamma$ ) and their associated five-body bremsstrahlung corrections ( $b \to s q \bar{q} g \gamma$ ) were not fully known. Previous work [17] calculated the subset of NLO contributions where four-body virtual corrections renormalize among themselves. However, the specific one-loop four-particle diagrams that must be supplemented by five-particle tree-level cuts (gluon bremsstrahlung) to ensure infrared (IR) finiteness and formal completeness at NLO were missing. This paper addresses this gap to formally complete the purely perturbative NLO calculation for $\bar{B} \to X_s\gamma$ at leading power.

Methodology
The authors perform a complete calculation of the interference terms between current-current and QCD penguin operators contributing to the decay rate. The methodology involves several advanced technical steps:

Diagram Generation and Algebra: Using QGRAF, the authors generated 176 diagrams for the four-body contribution and 400 for the five-body contribution per operator insertion. Dirac and color algebra were performed using FORM.
Treatment of $\gamma_5$ : Given the presence of axial currents, the authors employed the KKS scheme (Kreimer, Körner, Schilcher). This scheme maintains the anti-commutativity $\{\gamma_5, \gamma_\mu\} = 0$ but forbids the cyclicity of traces involving an odd number of $\gamma_5$ matrices. They utilized a "reading point" strategy and averaged over axial starting points to handle trace products consistently in $D = 4 - 2\epsilon$ dimensions.
Integral Reduction: The phase-space integrals were converted into four-loop propagator-type integrals via the method of reverse unitarity. These were reduced to a set of master integrals using integration-by-parts (IBP) identities and Laporta's algorithm (implemented in FIRE). The reduction yielded seven distinct integral topologies.
Master Integral Evaluation: The master integrals were evaluated analytically as functions of the photon energy cut variable $z$ $z$ and the dimensional regulator $\epsilon$ $ϵ$ .
- Simple integrals were computed via explicit phase-space integration.
- Complex integrals were solved using the method of differential equations, transformed into $\epsilon$ -form, and expressed in terms of Harmonic Polylogarithms (HPLs) and Goncharov polylogarithms.
- Boundary conditions were determined using asymptotic expansions (via hypergeometric functions or Mellin-Barnes representations) and by integrating over the full phase space (without cuts) to relate cut and cut-less integrals.
Renormalization and IR Subtraction:
- UV Renormalization: Performed in the on-shell scheme for the $b$ -quark mass and fields. Evanescent operators were included to handle the operator basis in $D$ dimensions.
- IR Subtraction: To handle collinear divergences arising from massless light quarks, the authors employed a subtraction scheme converting dimensional regulators into mass-regulated collinear logarithms ( $\ln(m_b/m_q)$ ). This involved convoluting tree-level and one-loop matrix elements with LO and NLO IR-subtraction kernels ( $S_0$ and $S_1$ ).

Key Contributions

Formal Completion: The paper provides the final missing piece required to formally complete the NLO QCD calculation for $\bar{B} \to X_s\gamma$ at leading power, specifically the $b \to s q \bar{q} \gamma$ one-loop diagrams and their associated $b \to s q \bar{q} g \gamma$ real radiation.
Analytic Results: The authors derive fully analytic results for the relevant matrix elements in terms of multiple polylogarithms, including the dependence on collinear logarithms.
Technical Implementation: The work demonstrates the application of modern integral reduction and differential equation techniques to complex four- and five-particle phase-space integrals with energy cuts.

Results
The numerical impact of these multi-parton corrections is found to be small.

Magnitude: The total correction $\Delta\Gamma$ (sum of four- and five-body contributions) is in the per-mille range relative to the total decay width. Relative to the leading-order two-body width, the correction is approximately $0.1\% - 0.2\%$ depending on the choice of the light quark mass parameter ( $m_q$ ).
Cancellation: There is a partial cancellation between the Leading Order (LO) and NLO contributions. Both the previously calculated NLO four-body terms (from ref. [17]) and the new five-body terms calculated here are negative, reducing the LO prediction.
Dependence: The results show a dependence on the renormalization scale $\mu$ and the light quark mass $m_q$ (via collinear logarithms). The inclusion of these terms improves the renormalization group invariance of the prediction compared to the LO result.
Scheme Dependence: The authors investigated the dependence on the heavy quark mass scheme (pole, 1S, $\overline{\text{MS}}$ ). While the pole and 1S schemes yield similar results, the $\overline{\text{MS}}$ scheme shows larger variations, a behavior attributed to the specific handling of the charm quark mass and the order of the perturbative expansion.

Significance
The paper claims that its primary significance lies in the formal completion of the NLO perturbative calculation for $\bar{B} \to X_s\gamma$ . By calculating the last missing multi-parton contributions, the authors ensure that the theoretical prediction is complete at the NLO level, which is a necessary prerequisite for matching the anticipated experimental precision of the Belle II experiment (projected uncertainty of $\pm 2.6\%$ ). While the numerical impact of these specific corrections is small due to CKM suppression and partial cancellations, their inclusion is essential for a rigorous theoretical framework. The authors note that a thorough investigation of scheme dependencies and higher-order power corrections requires embedding these results into the broader context of existing NNLO corrections, which is beyond the scope of this specific work.

The Puzzle of the Missing Pieces

What This Paper Does

The Result: A Small but Necessary Correction

Summary

More like this