⚛️ phenomenology

Precision Study of Semileptonic and Non-Leptonic $B_c$ Decays to $\eta_c$ and P Wave Charmonia

This paper presents a precision analysis of semileptonic and non-leptonic $B_c$ meson decays into P-wave charmonium states, utilizing NRQCD-based form factors constrained by experimental radiative decay data to provide updated predictions for branching fractions and lepton flavor universality ratios.

Original authors: Aritra Biswas, Soumitra Nandi, Shantanu Sahoo

Published 2026-07-08

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Aritra Biswas, Soumitra Nandi, Shantanu Sahoo

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 massive, bustling construction site. At the very bottom of the hierarchy, we have tiny workers called quarks. Usually, these workers pair up in twos or threes to build larger structures called mesons.

This paper focuses on a very specific, rare, and heavy-duty construction project: the $B_c$ meson. Think of the $B_c$ as a unique, heavy-duty crane built from two different types of heavy workers: a "bottom" quark and a "charm" quark. Because it's so heavy and made of two different heavy materials, it's hard to build (produce) and even harder to study.

The scientists in this paper are trying to understand exactly how this heavy crane falls apart (decays) into other structures, specifically into a family of "charmonium" buildings. These buildings come in different shapes and sizes. Some are simple, round spheres (called S-wave), while others are more complex, spinning structures with a "P-wave" shape (like a dumbbell or a figure-eight).

Here is a breakdown of what the paper does, using simple analogies:

1. The Problem: Missing Blueprints

To predict exactly how the $B_c$ crane falls apart, the scientists need "blueprints" called Form Factors. These are mathematical maps that tell you the probability of the crane turning into a specific shape.

The Issue: Until now, no one had a reliable, direct blueprint for the complex "P-wave" shapes. Previous attempts were like guessing the shape of a building based on a blurry photo or a rough sketch (model-dependent). There were no clear measurements from experiments or supercomputers (Lattice QCD) to guide them.

2. The Solution: Using "Radiative Decay" as a Mirror

Instead of guessing, the authors decided to use a clever trick. They looked at how these complex P-wave buildings emit light (radiative decay) when they change shape.

The Analogy: Imagine you want to know the exact shape of a spinning top, but you can't see it directly. However, you can see the light it reflects. By measuring the light (the photons) it emits, you can work backward to figure out the top's shape.
The Method: The team took real experimental data on how these particles emit light (like $\chi_{c0} \to \gamma\gamma$ or $\chi_{c1} \to J/\psi\gamma$ ). They used a theoretical framework called NRQCD (Non-Relativistic Quantum Chromodynamics) to translate that light data into the "shape" of the wave functions (the internal structure of the particles).
The Result: They successfully extracted the "derivatives of the wave functions." In our analogy, this is like measuring the exact curvature of the spinning top's surface. This gave them the missing numbers needed to build their blueprints.

3. Drawing the Map (The Form Factors)

Once they had the numbers from the light data, they used them to draw the full map (the Form Factors) for the $B_c$ decays.

The Technique: They used a mathematical tool called the z-expansion (BCL parametrization). Think of this as drawing a smooth, curved road on a map. Instead of guessing the whole road, they pinned down a few key points (using their new data) and used a mathematical rule to connect them smoothly.
The Innovation: They didn't just guess the road; they constrained the shape using real data, making their map much more reliable than previous "sketches."

4. The Predictions: What Happens Next?

With their new, data-driven blueprints, the authors made several predictions:

Semileptonic Decays (The "Leak"): They predicted how often the $B_c$ crane leaks a particle (a lepton) while turning into a P-wave building. They calculated the "Branching Fractions" (the odds of this happening).
- Key Finding: They calculated a ratio called R, which compares how often the crane leaks a heavy particle (tau) versus a light one (muon). This is a test for "Lepton Flavor Universality." If the ratio is different from the Standard Model prediction, it might mean new physics is at play. Their predictions for these ratios are:
  - $R(\chi_{c0}) \approx 0.185$
  - $R(\chi_{c1}) \approx 0.147$
  - $R(h_c) \approx 0.068$
Non-Leptonic Decays (The "Crash"): They also predicted what happens when the $B_c$ crashes into a light particle (like a pion) and turns into a P-wave building. They provided estimates for these rare events, which are currently very hard to measure but important for understanding the background noise in other experiments.
Production in Collisions: Finally, they looked at how these P-wave buildings are created in other scenarios, like when a Z-boson (a heavy particle from the early universe) decays, or when electrons and positrons smash together. They updated the expected rates for these events based on their new wave function data.

Summary

In short, this paper is like a team of architects who finally got a clear, high-resolution photo of a rare, complex building. Instead of guessing what the building looks like, they used that photo to draw a precise blueprint. With this new blueprint, they can now accurately predict how a heavy, rare crane ( $B_c$ ) will fall apart into these complex shapes. This helps other scientists know exactly what to look for in their experiments and ensures they aren't mistaking a "complex building" for a "simple sphere" when searching for new physics.

What the paper does NOT claim:

It does not claim to have discovered new physics (like a new force).
It does not claim these results have immediate medical or technological applications.
It does not claim to have measured these decays directly; it provides the theoretical predictions based on existing data that future experiments (like those at the LHC) can test.

Technical Summary: Precision Study of Semileptonic and Non-Leptonic $B_c$ Decays to $\eta_c$ and P-Wave Charmonia

Problem Statement
The $B_c$ meson, composed of two heavy quarks ( $b\bar{c}$ ), offers a unique laboratory for testing heavy-quark dynamics and searching for physics beyond the Standard Model (SM), particularly through Lepton Flavour Universality (LFU) violation. While semileptonic decays to S-wave charmonia (like $J/\psi$ ) are well-studied, decays to P-wave charmonia ( $\chi_{c0}, \chi_{c1}, h_c$ ) remain poorly understood. These P-wave modes are critical not only as potential backgrounds to S-wave analyses but also as independent probes of non-perturbative QCD effects and spin symmetry. A major obstacle in predicting these decays is the lack of reliable lattice QCD or direct experimental determinations of the relevant transition form factors over the full kinematic $q^2$ range. Existing predictions rely heavily on model-dependent approaches (e.g., perturbative QCD, quark models, QCD sum rules), leading to significant theoretical uncertainties.

Methodology
The authors employ a data-driven approach within the Non-Relativistic QCD (NRQCD) factorization framework to constrain the transition form factors for $B_c \to P$ -wave charmonia decays. The methodology proceeds in several key steps:

Extraction of Wave Function Parameters: Instead of relying on theoretical models for the non-perturbative inputs, the authors extract the derivatives of the radial wave functions at the origin ( $\psi'_{M}(0)$ ) for $\chi_{c0}, \chi_{c1},$ and $h_c$ , as well as the wave functions for $J/\psi$ and $\eta_c$ . This is achieved by performing a global $\chi^2$ fit to experimental data on radiative decays (e.g., $\chi_{c1} \to V\gamma$ , $\chi_{c0} \to \gamma\gamma$ , $h_c \to \eta_c\gamma$ ) and leptonic decays, combined with lattice QCD inputs for decay constants. Nuisance parameters, including relativistic corrections and higher-order operators, are treated with Gaussian priors.
Form Factor Calculation at Maximum Recoil: Using the extracted wave function parameters and the NRQCD analytic expressions (including Leading Order and Relativistic Corrections up to $O(v^2)$ ) from Ref. [45], the authors calculate the transition form factors at $q^2=0$ (maximum recoil). For $B_c \to \chi_{c1}$ and $h_c$ , where Next-to-Leading Order (NLO) corrections are unavailable, a conservative 20% uncertainty is assigned to account for missing higher-order perturbative effects.
Shape Determination and Extrapolation: To describe the form factors over the entire physical $q^2$ $q^{2}$ region, the authors utilize two complementary methods:
- Pole Expansion: Parameters are constrained using spin symmetry arguments, assuming that form factors for transitions to states with the same total spin $J$ share the same shape parameters ( $\alpha_i, \beta$ ). These parameters are first determined by fitting synthetic data generated for $B_c \to J/\psi$ and $B_c \to \eta_c$ (based on previous lattice and HQSS studies).
- $z$ -Expansion (BCL Parametrization): The pole-expansion results are used to generate pseudo-data points, which are then fitted using the model-independent Boyd-Caprini-Lellouch (BCL) parametrization. This approach ensures analyticity and unitarity constraints are satisfied, providing a robust description of the $q^2$ dependence.
Phenomenological Predictions: The constrained form factors are used to compute semileptonic branching fractions, LFU ratios ( $R(H)$ ), exclusive non-leptonic decay rates (using QCD factorization with non-factorizable corrections), and production cross-sections for P-wave charmonia in $Z$ boson decays and $e^+e^-$ annihilation.

Key Contributions and Results

Data-Driven Form Factors: The paper provides the first data-constrained estimates for $B_c \to \chi_{c0,1}, h_c$ form factors. By extracting wave function derivatives directly from experimental radiative decay rates, the authors minimize reliance on purely theoretical models.
LFU Ratios: The authors present precise Standard Model predictions for the LFU ratios:
- $R(\chi_{c0}) = 0.185(3)$
- $R(\chi_{c1}) = 0.147(26)$
- $R(h_c) = 0.068(2)$
  These ratios are theoretically clean, as hadronic uncertainties largely cancel between the numerator and denominator. The precision is limited primarily by the form factor uncertainties, which are estimated to be around 5% for $R(\chi_{c0})$ and $R(h_c)$ .
Semileptonic Branching Fractions: Updated predictions for semileptonic branching fractions are provided, e.g., $\mathcal{B}(B_c \to \eta_c \tau^- \bar{\nu}_\tau) = (1.90 \pm 0.49) \times 10^{-3}$ and $\mathcal{B}(B_c \to \eta_c \mu^- \bar{\nu}_\mu) = (6.51 \pm 1.97) \times 10^{-3}$ . The uncertainties remain relatively large (25–30%) due to the current lack of precise lattice inputs for P-wave form factors.
Non-Leptonic Decays: The study extends to exclusive non-leptonic two-body decays ( $B_c \to \chi_{cJ}, h_c + \pi, K, \rho, K^*$ ). Predictions are provided for branching fractions, including those with Next-to-Leading Order (NLO) factorizable corrections and Leading Order (LO) non-factorizable corrections. For instance, $\mathcal{B}(B_c \to \chi_{c0}\pi)$ is predicted to be $5.3(60) \times 10^{-3}$ (with corrections).
Production Modes: The paper updates predictions for P-wave charmonium production in $Z$ boson decays (e.g., $Z \to \chi_{c0}\gamma$ ) and double-charmonium production in $e^+e^-$ annihilation (e.g., $e^+e^- \to J/\psi + \chi_{c0}$ ), comparing these with previous literature estimates.

Significance and Claims
The authors claim that this work represents one of the earliest efforts to constrain $B_c \to P$ -wave charmonium form factors using experimental data rather than relying solely on theoretical models. By combining NRQCD formalism with a global fit to radiative decay data, they provide a "controlled and systematic basis" for evaluating decay widths.

The paper emphasizes that while the current predictions have sizable uncertainties (driven by the poorly constrained hadronic form factors), the derived LFU ratios are theoretically clean probes. The authors state that their predictions are ready to be tested by future experiments at the LHC and Belle II. They explicitly note that the current analysis is a stepping stone, anticipating that future lattice QCD and Light-Cone Sum Rule (LCSR) calculations will allow for a model-independent treatment of these form factors, thereby enabling more precise SM estimates and more sensitive searches for New Physics. The work does not claim to resolve existing tensions in $R(J/\psi)$ but highlights the necessity of understanding P-wave backgrounds and providing complementary precision tests in the heavy-heavy sector.

1. The Problem: Missing Blueprints

2. The Solution: Using "Radiative Decay" as a Mirror

3. Drawing the Map (The Form Factors)

4. The Predictions: What Happens Next?

Summary

More like this