Nonperturbative completion of the inclusive $\bar B\to X_s\gamma$ prediction: resolved-photon contributions and the prospects for a clean determination of the chromomagnetic coefficient $C_8$

The 1995 source paper by Kapustin, Ligeti and Politzer established that the $O_8$–$O_8$ contribution to the inclusive $\bar B\to X_s\gamma$ photon spectrum carries large collinear logarithms of the $b$-quark mass, resummed them via quark and gluon fragmentation functions down to an infrared scale $\mu\sim\Lambda_{\rm QCD}$, and posed the open questions that define this problem: the uncalculable nonperturbative decay-function inputs and residual $\mu$-dependence, the then-incomplete full NLO calculation of the rate and spectrum, and the prospect of a theoretically clean determination of the chromomagnetic Wilson coefficient $C_8$ from the photon spectrum around $E_\gamma\sim1$ GeV. The nonperturbative half of the question remains open: the resolved-photon (subleading-power, $O(\Lambda_{\rm QCD}/m_b)$) contributions, the modern form of the source paper's uncalculable fragmentation functions, are now the largest residual theory uncertainty on the rate and CP asymmetry, and the low-energy spectrum where $C_8$ could be determined has never been measured. The surviving open problem is the nonperturbative completion of the inclusive $\bar B\to X_s\gamma$ prediction — a first-principles or rigorously bounded determination of the resolved-photon contributions — together with the predictability and measurability of the $E_\gamma\sim1$ GeV spectrum for a clean $C_8$ determination.

Background

The inclusive radiative decay $\bar B\to X_s\gamma$, in which a $B$ meson decays into any strange hadronic state $X_s$ plus a photon, arises in the Standard Model only through one-loop (so-called penguin) diagrams and is therefore unusually sensitive to hypothetical new particles contributing at the same loop level. To compute it, one integrates out the heavy $W$ boson and top quark and works with an effective weak Hamiltonian built from local operators $O_i$ multiplied by Wilson coefficients $C_i$ (the short-distance couplings obtained by matching at the electroweak scale and evolving down with the renormalization group). Two dipole operators are central here: the electromagnetic dipole operator $O_7$, which emits the photon directly off the $b$ or $s$ quark and dominates the rate, and the chromomagnetic dipole operator $O_8$, which mediates the flavor-changing gluon emission $b\to s g$; their coefficients are $C_7$ and $C_8$. The source paper (Kapustin, Ligeti and Politzer, 1995) analyzed the contribution to the photon spectrum proportional to $C_8^2$, denoted $\Gamma_{88}$, arising from $b\to s g$ followed by photon emission off the light partons. In perturbation theory this piece contains collinear logarithms $\log(m_b^2/m_s^2)$: a photon emitted collinearly from a nearly on-shell massless quark or gluon produces a divergence cut off only by the light-quark mass. The paper showed that the leading logarithms $[\alpha_s\log(m_b^2/\mu^2)]^n$ can be summed to all orders by organizing the emission into quark and gluon to photon fragmentation functions (probability distributions for a parton to produce a photon carrying a given fraction of its energy), whose leading-logarithmic normalization is calculable by an argument due to Witten for the photon structure function (the distribution functions of a high-energy photon, measured in deep-inelastic electron–photon scattering). The price is an infrared matching scale $\mu$ of order $\Lambda_{\rm QCD}$, the intrinsic strong-interaction scale: the nonperturbative parts of the fragmentation functions are absorbed into $\mu$ and must be supplied from experiment or nonperturbative methods — they are not calculable in perturbation theory, leaving a parametrical uncertainty. The paper closed with the observations that define the problem: the full next-to-leading-order (NLO) calculation was then incomplete, with the technically most challenging parts still uncalculated; $\Gamma_{88}$ dominates the spectrum below about 1.2 GeV, so a measurement of the photon spectrum around $E_\gamma\sim1$ GeV — probably impossible at CLEO (the Cornell electron–positron collider experiment then taking the data), perhaps feasible at the then-planned asymmetric $B$ factories (colliders whose two beams of unequal energy boost the produced $B$ mesons in the laboratory) — would provide a theoretically clean determination of $C_8$; and combining the resummed perturbative pieces with nonperturbative inputs quantitatively remained to be done. The contemporaneous CLEO measurement (Alam et al., 1995) had just made the missing NLO calculation the limiting factor in constraining new physics. The subsequent literature transformed this problem. The NLO program was completed (Chetyrkin, Misiak and Münz 1997; Buras, Czarnecki, Misiak and Urban 2002; Buras and Misiak 2002), the last multi-parton NLO pieces were finished in 2025–26 (Brune, Huber and Moos), and the next-to-next-to-leading-order (NNLO) prediction is now complete with exact charm-mass dependence, giving a Standard Model branching ratio of $(3.54\pm0.14)\times10^{-4}$ for $E_\gamma>1.6$ GeV, in agreement with the experimental average $(3.49\pm0.19)\times10^{-4}$ (Misiak et al. 2007, 2020, 2026). The photon spectrum near its endpoint was reformulated in soft-collinear effective theory (SCET, the effective field theory describing the collinear and soft radiation in such decays) as a perturbatively calculable coefficient convoluted with a nonperturbative shape function, which describes the distribution of the residual light-parton momentum inside the $B$ meson (Kagan and Neubert 1999; Ligeti, Luke, Manohar and Wise 1999). At subleading power, $O(\Lambda_{\rm QCD}/m_b)$, Benzke, Lee, Neubert and Paz (2010) identified resolved-photon contributions — photons coupling to light partons rather than directly to the weak interaction — which first appear at order $1/m_b$, are described by new nonperturbative soft and jet functions, cannot be expressed in the local heavy-quark expansion, and now form the largest single component of the theory uncertainty on the total rate (Gunawardana and Paz 2019). The subleading-power $O_8$–$O_8$ channel, the descendant of the source paper's $\Gamma_{88}$, suffers from endpoint divergences that block a fixed-order treatment; Hurth and Szafron (2023) proved a refactorisation theorem permitting its controlled resummation, a program still under construction (Bartocci, Böer and Hurth 2025; Benzke, Garzelli and Hurth 2026). Experimentally, the spectrum has been measured by the Belle and Belle II experiments (the detectors at the KEK asymmetric $B$ factory in Japan) only down to photon energies of 1.7–1.8 GeV; the $E_\gamma\sim1$ GeV region has never been measured, and no determination of $C_8$ from the spectrum exists.

Problem Statement

The problem is the completion of the quantitative prediction of the inclusive $\bar B\to X_s\gamma$ rate and photon spectrum posed by the 1995 source paper. Its perturbative components — the full NLO (and subsequently NNLO) calculation of the rate and spectrum, and the combination of resummed perturbative contributions with nonperturbative inputs at leading power — have been discharged by later literature; they are therefore excluded from the open core. The surviving open objective, whose two components are both required and serve the single goal of a complete, nonperturbatively controlled prediction of this decay, is:

(1) Determine, from first principles or with rigorously controlled uncertainties, the resolved-photon (subleading-power, $O(\Lambda_{\rm QCD}/m_b)$) contributions to the total $\bar B\to X_s\gamma$ rate and photon spectrum — including the subleading-power continuation of the $O_8$–$O_8$ channel whose leading-power collinear logarithms the source paper resummed via fragmentation functions down to the infrared matching scale $\mu\sim\Lambda_{\rm QCD}$. These terms are currently parameterized by nonperturbative soft functions (describing soft, low-momentum gluon radiation), subleading jet functions (describing collinear partonic radiation) and subleading shape functions (the subleading-power generalizations of the $B$-meson momentum-distribution function) that are not determined by other observables, and they constitute the largest residual theory uncertainty on the rate and on its CP asymmetry (the small difference between the rates of the CP-conjugate $B$ and $\bar B$ decays); the goal is a determination, or a rigorously derived bound, of these contributions that removes the resolved-photon uncertainty as the dominant source of theoretical error on the rate and on its CP asymmetry.

(2) Establish whether the photon spectrum in the low-energy region ($E_\gamma\sim1$ GeV in the $B$-meson rest frame), where the $O_8$–$O_8$ contribution dominates, can be predicted and measured with sufficient control — including the resolved-photon contamination that first enters at this order and was unknown in 1995 — to provide the theoretically clean determination of the chromomagnetic Wilson coefficient $C_8$ that the source paper proposed. This component uses the resolved-photon and $O_8$–$O_8$ determination sought in (1) as its theory input, and its deliverable is a definite verdict rather than an extraction: it comprises the experimental feasibility of the measurement (no spectrum below 1.7–1.8 GeV has ever been measured) and the fidelity with which $C_8$, or the $O_8$–$O_8$ contribution, could be extracted from such data.

References

1. A. Kapustin, Z. Ligeti, H. David Politzer, 'Leading logarithms of the b quark mass in inclusive B → Xs γ decay', Phys. Lett. B 357 (1995) 653–658, DOI 10.1016/0370-2693(95)00962-K, arXiv:hep-ph/9507248, https://arxiv.org/abs/hep-ph/9507248
2. A. Ali, C. Greub, 'Photon energy spectrum in $B\to X_s\gamma$ and comparison with data', Phys. Lett. B 361 (1995) 146–154, DOI 10.1016/0370-2693(95)01118-A, arXiv:hep-ph/9506374, https://arxiv.org/abs/hep-ph/9506374
3. M. S. Alam et al. (CLEO Collaboration), 'First Measurement of the Rate for the Inclusive Radiative Penguin Decay b → s γ', Phys. Rev. Lett. 74 (1995) 2885, DOI 10.1103/PhysRevLett.74.2885, https://doi.org/10.1103/PhysRevLett.74.2885
4. K. G. Chetyrkin, M. Misiak, M. Münz, 'Weak radiative B-meson decay beyond leading logarithms', Phys. Lett. B 400 (1997) 206, DOI 10.1016/S0370-2693(97)00324-9, arXiv:hep-ph/9612313, https://arxiv.org/abs/hep-ph/9612313
5. A. L. Kagan, M. Neubert, 'QCD anatomy of B → X_s gamma decays', Eur. Phys. J. C 7 (1999) 5, DOI 10.1007/s100529800959, arXiv:hep-ph/9805303, https://arxiv.org/abs/hep-ph/9805303
6. Z. Ligeti, M. E. Luke, A. V. Manohar, M. B. Wise, 'The B → X_s gamma photon spectrum', Phys. Rev. D 60 (1999) 034019, DOI 10.1103/PhysRevD.60.034019, arXiv:hep-ph/9903305, https://arxiv.org/abs/hep-ph/9903305
7. A. J. Buras, A. Czarnecki, M. Misiak, J. Urban, 'Completing the NLO QCD calculation of B → X_s γ', Nucl. Phys. B 631 (2002) 219–238, DOI 10.1016/S0550-3213(02)00261-4, https://doi.org/10.1016/S0550-3213(02)00261-4
8. A. J. Buras, M. Misiak, 'B → X_s gamma after completion of the NLO QCD calculations', Acta. Phys. Polon. B 33 (2002) 2597, arXiv:hep-ph/0207131, https://arxiv.org/abs/hep-ph/0207131
9. S. J. Lee, M. Neubert, G. Paz, 'Enhanced Non-local Power Corrections to the B̄ → X_s γ Decay Rate', Phys. Rev. D 75 (2007) 114005, DOI 10.1103/PhysRevD.75.114005, arXiv:hep-ph/0609224, https://arxiv.org/abs/hep-ph/0609224
10. M. Misiak, H. M. Asatrian, K. Bieri, M. Czakon, A. Czarnecki, T. Ewerth, A. Ferroglia, P. Gambino, M. Gorbahn, C. Greub, U. Haisch, A. Hovhannisyan, T. Hurth, A. Mitov, V. Poghosyan, M. Ślusarczyk, M. Steinhauser, 'Estimate of B(B̄ → X_s γ) at O(α_s²)', Phys. Rev. Lett. 98 (2007) 022002, DOI 10.1103/PhysRevLett.98.022002, arXiv:hep-ph/0609232, https://arxiv.org/abs/hep-ph/0609232
11. M. Benzke, S. J. Lee, M. Neubert, G. Paz, 'Factorization at Subleading Power and Irreducible Uncertainties in B̄ → X_s γ Decay', JHEP 08 (2010) 099, DOI 10.1007/JHEP08(2010)099, arXiv:1003.5012, https://arxiv.org/abs/1003.5012
12. M. Kamiński, M. Misiak, M. Poradziński, 'Tree-level contributions to B → X_s γ', Phys. Rev. D 86 (2012) 094004, DOI 10.1103/PhysRevD.86.094004, arXiv:1209.0965, https://arxiv.org/abs/1209.0965
13. A. Gunawardana, G. Paz, 'Reevaluating uncertainties in B̄ → X_s γ decay', JHEP 11 (2019) 141, DOI 10.1007/JHEP11(2019)141, arXiv:1908.02812, https://arxiv.org/abs/1908.02812
14. M. Misiak, A. Rehman, M. Steinhauser, 'Towards B̄ → X_s γ at the NNLO in QCD without interpolation in m_c', JHEP 06 (2020) 175, DOI 10.1007/JHEP06(2020)175, arXiv:2002.01548, https://arxiv.org/abs/2002.01548
15. F. U. Bernlochner, H. Lacker, Z. Ligeti, I. W. Stewart, F. J. Tackmann, K. Tackmann (SIMBA collaboration), 'Precision Global Determination of the B → X_s γ Decay Rate', Phys. Rev. Lett. 127 (2021) 102001, DOI 10.1103/PhysRevLett.127.102001, arXiv:2007.04320, https://arxiv.org/abs/2007.04320
16. T. Hurth, R. Szafron, 'Refactorisation in subleading B̄ → X_s γ', Nucl. Phys. B 991 (2023) 116200, DOI 10.1016/j.nuclphysb.2023.116200, arXiv:2301.01739, https://arxiv.org/abs/2301.01739
17. R. Bartocci, P. Böer, T. Hurth, 'Renormalisation group evolution of the shape function g_17 in B̄ → X_s γ and B̄ → X_s ℓ⁺ℓ⁻ at subleading power', JHEP 04 (2025) 066, DOI 10.1007/JHEP04(2025)066, arXiv:2411.16634, https://arxiv.org/abs/2411.16634
18. M. Benzke, M. V. Garzelli, T. Hurth, 'Update of the nonlocal sub-leading $O_1$–$O_7$ contribution to $\bar B\to X_s\gamma$ at LO', Phys. Rev. D 113 (2026) 076006, DOI 10.1103/g6ds-ldx4, arXiv:2512.08902, https://arxiv.org/abs/2512.08902
19. R. Bartocci, P. Böer, T. Hurth, 'NLO analysis of the subleading-power $Q_1$–$Q_7^\gamma$ interference in $\bar B\to X_s\gamma$ at large photon energies', PoS EPS-HEP2025 (2026) 301, DOI 10.22323/1.485.0301, arXiv:2510.18811, https://arxiv.org/abs/2510.18811
20. K. Brune, T. Huber, L.-T. Moos, 'Multi-parton contributions to B̄ → X_s γ at NLO', JHEP 01 (2026) 142, DOI 10.1007/JHEP01(2026)142, arXiv:2509.22564, https://arxiv.org/abs/2509.22564
21. 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, 'The inclusive B̄ → X_s γ decay rate with higher precision', arXiv:2608.29863, https://arxiv.org/abs/2608.29863
22. M. Czaja, M. Czakon, T. Huber, M. Misiak, M. Niggetiedt, A. Rehman, K. Schönwald, M. Steinhauser, 'NNLO QCD corrections to the weak radiative B-meson decay with exact dependence on m_c', arXiv:2608.29864, https://arxiv.org/abs/2608.29864
23. Belle II Collaboration (F. Abudinén et al.), 'Measurement of the photon-energy spectrum in inclusive B → X_s γ decays identified using hadronic decays of the recoil B meson in 2019–2021 Belle II data', arXiv:2210.10220, https://arxiv.org/abs/2210.10220
24. A. Limosani et al. (Belle Collaboration), 'Measurement of Inclusive Radiative B-Meson Decays with a Photon Energy Threshold of 1.7 GeV', Phys. Rev. Lett. 103 (2009) 241801, DOI 10.1103/PhysRevLett.103.241801, https://doi.org/10.1103/PhysRevLett.103.241801
25. C. Schwanda et al. (Belle Collaboration), 'Measurement of the moments of the photon energy spectrum in B → X_s γ decays and determination of |V_cb| and m_b at Belle', Phys. Rev. D 78 (2008) 032016, DOI 10.1103/PhysRevD.78.032016, arXiv:0803.2158, https://arxiv.org/abs/0803.2158