Study of Form Factors and Observables in $B_c^- \rightarrow \bar{D}^{(*)0}\ell^-\barν_{\ell}$ and $B_c^- \rightarrow D^{(*)-}\ell^+\ell^-$ decays

This paper investigates the Standard Model predictions for $B_c^-$ decays into charmed mesons and leptons by employing perturbative QCD form factors constrained by lattice QCD inputs and heavy quark spin symmetry to calculate branching fractions, lepton flavor violating observables, and detailed angular distributions.

The Gist

Imagine the universe as a giant, bustling construction site. At this site, there are massive, heavy machines called $B_c$ mesons. These machines are unique because they are built from two very heavy parts stuck together: a "bottom" quark and a "charm" quark. Unlike other machines in the family that are built from one heavy part and one light part, these two heavy parts make the $B_c$ meson behave differently.

This paper is a detailed blueprint and a set of predictions about how these $B_c$ machines fall apart (decay) into smaller, simpler machines. Specifically, the authors are looking at two types of breakdowns:

1. The "Standard" Breakdown: Where the machine splits into a lighter car ($\bar{D}$ or $D^*$) and a pair of particles (a lepton and a neutrino).
2. The "Rare" Breakdown: A much more unusual event where the machine splits into a lighter car and a pair of charged particles (like an electron and a positron) without a neutrino. This is rare because it's like a car spontaneously turning into two other cars and a pair of twins without any external help—it only happens through complex, hidden loops in the laws of physics.

Here is a simple breakdown of what the authors did and found:

1. The Problem: We Didn't Know the "Shape" of the Machine

To predict how these machines break, you need to know exactly how the parts inside are arranged. In physics, this arrangement is described by something called a wave function (or Light Cone Distribution Amplitude). Think of this as the "blueprint" or the "DNA" of the machine.

In previous studies, scientists just guessed what this blueprint looked like, picking a random shape and hoping it was right. It was like trying to predict how a car crashes without knowing if it's a sedan or a truck.

The Innovation:
The authors of this paper decided to stop guessing. They used a "data-driven" approach. They took existing, high-precision measurements from other experiments (like the HPQCD lattice data) and worked backward. They asked: "What shape of the blueprint would make our math match the real-world data?"

They treated the shape of the blueprint as a mystery variable and used a statistical method (like a super-advanced curve-fitting game) to find the exact numbers that fit the data best. This allowed them to create a much more accurate blueprint for the $B_c$ and $D$ mesons.

2. The Bridge: Connecting the Known to the Unknown

The authors had a lot of data on how a $B$ meson (a different machine) breaks down, but they needed to know about the $B_c$ meson. They used a set of rules called Heavy Quark Spin Symmetry.

Think of this like a translator. If you know how a heavy truck ($B$) behaves, and you know the rules of the road (symmetry), you can predict how a slightly different heavy truck ($B_c$) will behave, even if you haven't seen it crash yet. They used these rules to translate their new, accurate blueprints from the known machines to the unknown ones, filling in the gaps for the entire range of possible outcomes.

3. The Predictions: What Happens When They Break?

Once they had the correct blueprints and the translation rules, they ran the numbers to predict what happens when these machines break. They calculated:

• Branching Fractions: How often does a specific type of breakdown happen? (e.g., "Out of 10,000 $B_c$ machines, how many will turn into a $D^*$ and a tau particle?")
• Lepton Flavor Universality: The Standard Model says that electrons, muons, and taus should behave exactly the same way, except for their weight. The authors calculated the ratio of heavy tau decays to light electron/muon decays to see if nature follows the rules perfectly.
• Angular Observables: This is the most detailed part. When the machine breaks, the pieces fly off in specific directions. The authors predicted the angles at which these pieces would fly. Imagine a pinball machine where the ball bounces off flippers; they predicted exactly where the ball would land. These angles are very sensitive to "New Physics"—if the ball lands somewhere unexpected, it might mean there are new, unknown forces at play.

4. The Results

• Precision: Their predictions are much more precise than previous guesses because they used real data to fix the blueprints.
• The "Clean" Observables: They identified specific angles and ratios that are "clean," meaning they are less affected by the messy internal details of the machine and more likely to show us if the Standard Model is wrong.
• CP Asymmetry: They predicted a tiny difference between how a machine breaks and how its "mirror image" (antimatter) breaks. This difference is very small but non-zero, which is a standard prediction of the current laws of physics.

Summary

In short, this paper is like a team of engineers who stopped guessing how a complex machine works. Instead, they measured the machine's vibrations to reverse-engineer its exact internal design. With this new, accurate design, they simulated thousands of crash scenarios to predict exactly how often the machine breaks, what pieces fly off, and in which direction.

Their goal isn't to build a new car, but to provide a baseline. If future experiments (like those at the LHCb detector) see these machines breaking in a way that doesn't match these precise predictions, it will be a huge signal that there is "New Physics" hiding in the shadows, waiting to be discovered.

Technical Summary

1. Problem Statement

The paper addresses the lack of precise theoretical predictions for the semileptonic and rare decays of the $B_c$ meson, specifically the channels $B_c^- \to \bar{D}^{(*)0}\ell^-\bar{\nu}_\ell$ and $B_c^- \to D^{(*)-}\ell^+\ell^-$.

• Context: While the $B_c$ meson (a bound state of $b$ and $c$ quarks) offers a unique laboratory for testing the Standard Model (SM) and probing New Physics (NP), precise calculations of its decay rates and angular observables are hindered by uncertainties in hadronic form factors over the full physical $q^2$ (momentum transfer) range.
• Challenge: Existing perturbative QCD (pQCD) calculations often rely on fixed, model-dependent parameters for Light Cone Distribution Amplitudes (LCDAs), leading to large theoretical uncertainties. Furthermore, pQCD is reliable only at low $q^2$, while lattice QCD is reliable only near the zero-recoil point ($q^2_{max}$). Bridging these two regimes to predict observables across the entire kinematic range is a significant challenge.

2. Methodology

The authors employ a hybrid, data-driven approach combining perturbative QCD (pQCD), Heavy Quark Spin Symmetry (HQSS), and Lattice QCD inputs.

A. Framework and Wave Functions

• Modified pQCD: The analysis uses a modified pQCD framework tailored for the $B_c$ meson. Unlike standard $B$ meson decays, the $B_c$ involves two heavy quarks. The authors modify the Sudakov factors to account for the finite charm quark mass scale, improving the treatment of $k_T$ resummation.
• LCDA Extraction: Instead of fixing the shape parameters of the meson wave functions (LCDAs) to arbitrary values, the authors treat them as free parameters. They perform a global $\chi^2$ minimization using:
  • Lattice QCD data from HPQCD and MILC for $B \to D^{(*)}$ and $B_c \to D$ form factors at $q^2=0$.
  • Light-Cone Sum Rule (LCSR) inputs.
  • This allows for a unified, data-driven determination of the shape parameters ($\omega_{B_c}, \omega_B, C_{D^{(*)}}, \omega_{D^{(*)}}$).

B. Extending Form Factors to Full $q^2$ Range

Since pQCD is only valid at low $q^2$, the authors use Heavy Quark Effective Theory (HQET) and Heavy Quark Spin Symmetry (HQSS) to connect the $B_c \to D$ and $B_c \to D^*$ form factors.

• Universal Functions: They express all 10 relevant form factors (3 for $B_c \to D$ and 7 for $B_c \to D^*$) in terms of two universal soft functions, $\Sigma_1(w)$ and $\Sigma_2(w)$, where $w$ is the recoil parameter.
• Taylor Expansion: These universal functions are expanded in a Taylor series around the zero-recoil point ($w=1$) using HPQCD lattice data to extract the expansion coefficients.
• BCL Parametrization: To cover the full physical $q^2$ region, the authors fit the form factors using the Boyd-Grinstein-Lebed (BCL) $z$-expansion. They combine:
  1. pQCD predictions at $q^2=0$ (constraining the low-$q^2$ behavior).
  2. HQSS-derived values at high $q^2$ (constraining the high-$q^2$ behavior).
  3. Lattice data points.
  • Systematic uncertainties from truncating the $z$-expansion and higher-twist effects are conservatively estimated (assigned a 30% uncertainty to LO predictions).

C. Observable Prediction

With the form factors determined over the full $q^2$ range, the authors calculate:

• Branching Fractions: For charged current ($B_c \to \bar{D}^{(*)0}\ell\bar{\nu}$) and Flavor Changing Neutral Current (FCNC) rare decays ($B_c \to D^{(*)-}\ell^+\ell^-$ and $B_c \to D^{(*)-}\nu\bar{\nu}$).
• Angular Observables: A detailed analysis of the cascade decay $B_c^- \to D^{*-}(\to D^0\pi^-)\ell^+\ell^-$, including CP-averaged angular coefficients ($S_i$), forward-backward asymmetry ($A_{FB}$), polarization fractions ($F_L, F_T$), and "clean" angular observables ($P_i, P'_i$) that are less sensitive to hadronic uncertainties.

3. Key Contributions

1. Data-Driven LCDA Parameters: The paper provides the first extraction of LCDA shape parameters for $B_c$ and $D^{(*)}$ mesons using a global fit to lattice and LCSR data, rather than relying on fixed model assumptions. This yields precise values for $\omega_{B_c} \approx 1.06$ GeV and $\omega_B \approx 0.50$ GeV.
2. Unified Form Factor Description: By leveraging HQSS relations, the authors successfully link $B_c \to D$ and $B_c \to D^*$ transitions, allowing the extraction of $B_c \to D^*$ form factors (which lack direct lattice data) using $B_c \to D$ lattice inputs.
3. Comprehensive SM Predictions: The work provides the most precise Standard Model predictions to date for a wide array of observables in $B_c$ decays, including branching ratios, lepton flavor universality (LFU) ratios ($R_{D^{(*)}}$), and detailed angular distributions for both light ($e, \mu$) and heavy ($\tau$) leptons.
4. Systematic Uncertainty Handling: The authors rigorously account for theoretical uncertainties, including higher-order loop corrections, higher-twist effects, and truncation errors in the $z$-expansion, providing robust error estimates.

4. Key Results

• Form Factors at $q^2=0$:
  • $F_+^{B_c \to D}(0) = 0.179(32)$
  • $A_0^{B_c \to D^*}(0) = 0.201(70)$
  • These values show good agreement with previous pQCD and Lattice results but with improved error propagation.
• Branching Fractions:
  • Predicted $\mathcal{B}(B_c^- \to \bar{D}^{*0}\ell^-\bar{\nu}_\ell) \approx 1.25 \times 10^{-4}$ (significantly more precise than previous pQCD estimates).
  • Predicted $\mathcal{B}(B_c^- \to D^{*-}\ell^+\ell^-) \approx 1.51 \times 10^{-8}$.
  • The $\tau$-mode branching fractions are suppressed due to phase space, consistent with SM expectations.
• Lepton Flavor Universality (LFU):
  • Calculated ratios $R_{D^*} = \mathcal{B}(B_c \to D^*\tau\nu)/\mathcal{B}(B_c \to D^*\ell\nu) = 0.632(22)$ and $R_{D^*}^{\ell^+\ell^-} = 0.166(10)$. These serve as SM baselines for future NP searches.
• Angular Observables:
  • Identified zero-crossing points for angular coefficients $S_5$ ($q^2_0 \approx 2.21$ GeV$^2$) and $S_{6s}$ ($q^2_0 \approx 3.49$ GeV$^2$).
  • Provided $q^2$-binned predictions for clean observables $P_2$ and $P'_5$, which are sensitive to NP. The authors note that current SM predictions for these in the $B_c$ channel are crucial for interpreting future LHCb data, similar to the anomalies seen in $B \to K^*\mu^+\mu^-$.
• CP Asymmetry: Predicted small but non-zero CP asymmetries ($A_{CP}$) in rare decays, primarily driven by the phase in CKM elements $V_{ub}$ and $V_{td}$.

5. Significance

This work is significant for the upcoming high-luminosity era at the LHC (specifically LHCb), where $B_c$ meson production rates are expected to increase by orders of magnitude.

• Precision Baseline: By reducing theoretical uncertainties through data-driven LCDA extraction and rigorous error propagation, the paper establishes a robust Standard Model baseline.
• New Physics Probe: The precise predictions for angular observables (like $P_2$ and $P'_5$) and LFU ratios allow experimentalists to distinguish between SM effects and potential New Physics contributions (e.g., leptoquarks or $Z'$ bosons) with greater confidence.
• Methodological Advancement: The approach of combining pQCD, HQSS, and lattice inputs to span the full $q^2$ range offers a template for analyzing other heavy-heavy or heavy-light meson decays where direct lattice calculations are sparse.

In summary, the paper transforms the theoretical understanding of $B_c$ decays from model-dependent estimates to precision predictions, directly enabling the next generation of experimental tests of the Standard Model.