Semileptonic decay form factors of $\Xi_b^0 \rightarrow \Xi_c^+\ell\bar{\nu}_{\ell}$ in HQET

This paper investigates the semileptonic decay $\Xi_b^0 \rightarrow \Xi_c^+\ell\bar{\nu}_{\ell}$ using a phenomenological quark model combined with Heavy Quark Effective Theory to calculate form factors, decay rates, and the lepton flavor universality ratio, finding results consistent with existing theoretical predictions.

The Gist

Imagine the universe as a giant, bustling construction site. In this site, there are tiny, fundamental building blocks called quarks. Usually, they stick together in groups of three to form larger structures called baryons (which include protons and neutrons).

This paper is like a detailed architectural blueprint for a very specific, rare renovation project happening on that construction site. Here is the story of what the scientists did, explained simply:

1. The Cast of Characters

• The Heavy Movers: The main characters are two heavy baryons. One is the $\Xi^0_b$ (Xi-zero-bottom), which contains a very heavy "bottom" quark. The other is the $\Xi^+_c$ (Xi-plus-charm), which contains a slightly lighter "charm" quark.
• The Transformation: The scientists are studying a "semileptonic decay." Think of this as the heavy-bottom baryon deciding to downsize. It sheds its heavy bottom quark, swaps it for a lighter charm quark, and in the process, spits out a pair of particles: a lepton (like an electron or a tau particle) and a neutrino (a ghost-like particle that barely interacts with anything).
• The Goal: They want to calculate exactly how this transformation happens. How fast does it go? How much energy is released? And does it matter if the lepton is an electron or a heavy tau?

2. The Toolkit: The "Hyper-Central" Model

To predict how this dance happens, the authors used a specific mathematical tool called the Hypercentral Constituent Quark Model (HCQM).

• The Analogy: Imagine the three quarks inside the baryon are three dancers holding hands in a circle, spinning around a common center.
• The Problem: Calculating the exact movement of three dancers who are constantly pulling and pushing each other is incredibly hard.
• The Solution: The scientists used a special coordinate system (Jacobi coordinates) that treats the whole trio as a single, unified "super-dancer" moving in a 6-dimensional space. They gave this super-dancer a "potential energy" map (a mix of a spring-like force and a gravitational-like pull) to simulate how the quarks hold each other together.

3. The "Form Factors": The Shape-Shifting Rules

When the heavy baryon changes into the lighter one, it doesn't just snap instantly. It goes through a transition. The scientists needed to describe the "shape" of this transition.

• The Analogy: Think of the baryon as a piece of clay. When it changes from a heavy shape to a light shape, it stretches and squishes. The Form Factors are the mathematical rules that describe exactly how the clay stretches.
• The Heavy Quark Trick (HQET): Because the bottom quark is so heavy, it acts like a heavy anchor. The scientists used a theory called Heavy Quark Effective Theory (HQET).
  • Simple version: If you have a giant anchor (the heavy quark) and a tiny boat (the light quarks), the boat's movement is mostly dictated by the anchor. This simplifies the math, allowing them to predict the "stretching rules" (form factors) with high precision.
• The Result: They found that two specific rules (called $f_1$ and $g_1$) are the "bosses" of the process. They do almost all the work, while the other four rules are just minor corrections, like tiny ripples on a wave.

4. The Race: Electrons vs. Taus

The decay can happen in two ways:

1. The Sprint: The baryon turns into a light electron and a neutrino.
2. The Marathon: The baryon turns into a heavy tau particle and a neutrino.

The scientists calculated the "speed" (decay rate) for both.

• The Finding: The "Sprint" (electron) happens much more often than the "Marathon" (tau).
• The Ratio: They calculated a ratio called $R(\Xi_c)$, which is roughly 0.3. This means for every 10 times the heavy baryon decays into an electron, it only decays into a tau about 3 times.

5. Why Does This Matter?

• Testing the Rules: The Standard Model is the "rulebook" of physics. By calculating these numbers so precisely, the scientists are checking if the universe is following the rulebook. If future experiments find different numbers, it could mean there are new particles or forces we don't know about yet.
• The "R-Puzzle": There is a mystery in physics called the "R-puzzle," where some heavy particle decays seem to happen more often than the rulebook predicts. This paper provides a new, very careful calculation to see if the $\Xi_b$ baryon fits the pattern or if it's part of the mystery.
• Future Proofing: Since we haven't measured this specific decay in a lab yet, this paper acts as a target. When experimentalists at big machines (like the LHC) finally catch this decay in action, they will compare their numbers to this paper. If they match, our understanding of the universe is solid. If they don't, we might be on the verge of a huge discovery!

Summary

In short, these scientists built a sophisticated mathematical simulation of a heavy particle changing into a lighter one. They used a "heavy anchor" theory to simplify the complex dance of quarks, calculated the exact "stretching rules" of the transformation, and predicted how often it happens with different types of particles. Their results agree with other theories, giving us confidence that our current understanding of the subatomic world is on the right track.

Technical Summary

1. Problem Statement

The paper addresses the theoretical description of the semileptonic decay of the bottom-strange baryon $\Xi^0_b$ into the charmed-strange baryon $\Xi^+_c$ ($\Xi^0_b \to \Xi^+_c \ell \bar{\nu}_\ell$, where $\ell = e, \tau$).

• Context: Heavy-to-heavy semileptonic decays are crucial for testing the Standard Model (SM) and extracting Cabibbo-Kobayashi-Maskawa (CKM) matrix elements (specifically $|V_{cb}|$).
• Gap: While the $\Lambda^0_b$ baryon has been extensively studied, other bottom baryons like $\Xi^0_b$ and $\Xi^-_b$ are less explored theoretically. Recent experimental observations of related hadronic decays suggest that semileptonic transitions may soon be observable.
• Objective: To provide a consistent theoretical analysis of the $\Xi^0_b \to \Xi^+_c$ transition, calculating the transition form factors, decay rates, branching ratios, and the Lepton Flavour Universality (LFU) ratio $R(\Xi_c)$, while incorporating subleading corrections in the Heavy Quark Effective Theory (HQET) framework.

2. Methodology

The authors employ a hybrid theoretical approach combining the Hypercentral Constituent Quark Model (HCQM) with Heavy Quark Effective Theory (HQET).

A. Wave Function and Mass Calculation (HCQM)

• Framework: The internal dynamics of the baryons are modeled using Jacobi coordinates and a six-dimensional hyperradial Schrödinger equation.
• Potential: The interaction potential consists of a hyperCoulomb term plus a linear confining term ($V(x) = \tau/x + \beta x + V_0$) and a perturbative spin-dependent hyperfine interaction.
• Parameters: Model parameters (quark masses $m_q$, potential strengths $\beta, V_0$, and hyperfine parameters $A, x_0$) are tuned to reproduce the experimental ground-state masses of $\Xi_b$ and $\Xi_c$ baryons.
• Output: This yields the ground-state wavefunctions $\psi_{\nu\gamma}(x)$, which are used to compute the baryon masses and the Isgur-Wise function.

B. Form Factor Calculation (HQET)

• Isgur-Wise Function: The leading-order form factors are reduced to a single universal function, $\xi(\omega)$, where $\omega$ is the velocity transfer. The authors expand $\xi(\omega)$ in a Taylor series around the zero-recoil point ($\omega=1$) up to the quadratic order:
  $$\xi(\omega) = 1 - \rho^2(\omega - 1) + c(\omega - 1)^2$$
  The slope ($\rho^2$) and curvature ($c$) parameters are calculated directly from the HCQM wavefunctions.
• $1/m_Q$ Corrections: To improve accuracy beyond the infinite mass limit, subleading corrections of order $1/m_Q$ are included. These corrections introduce additional functions ($B_1, B_2$) and depend on the parameter $\bar{\Lambda} = m_{\Xi_b} - m_b$.
• Form Factor Mapping: The six invariant form factors ($F_{1,2,3}$ and $G_{1,2,3}$) are derived from the HQET expansion. These are then converted into the helicity basis form factors ($f_{1,2,3}$ and $g_{1,2,3}$) required for numerical evaluation of decay rates.

C. Decay Rate and Observables

• Helicity Formalism: The differential decay rate $d\Gamma/dq^2$ is calculated using the helicity amplitudes derived from the form factors.
• Integration: The total decay width ($\Gamma$) is obtained by integrating over the kinematically allowed $q^2$ range.
• Branching Ratio & LFU: The branching ratio is calculated using the known lifetime of $\Xi^0_b$. The LFU ratio $R(\Xi_c) = \mathcal{B}(\Xi^0_b \to \Xi^+_c \tau \bar{\nu}_\tau) / \mathcal{B}(\Xi^0_b \to \Xi^+_c e \bar{\nu}_e)$ is computed to test lepton universality.

3. Key Contributions

1. Consistent Framework: The paper successfully integrates HCQM wavefunctions with HQET expansions to calculate form factors including $1/m_Q$ corrections, providing a more rigorous treatment than models relying solely on phenomenological fits.
2. First-Principles Parameters: The slope ($\rho^2 = 1.83$) and curvature ($c = 0.84$) of the Isgur-Wise function are derived directly from the calculated wavefunctions, rather than being fitted to data.
3. Comprehensive Observable Set: The study provides a complete set of predictions for:
   • Six transition form factors ($f_i, g_i$) as a function of $q^2$.
   • Total decay widths for both electronic and tauonic channels.
   • Branching ratios.
   • The LFU ratio $R(\Xi_c)$.
4. Uncertainty Analysis: The authors perform a systematic error analysis by varying quark mass parameters ($\pm 0.05$ GeV), quantifying the theoretical uncertainties in the decay rates and branching ratios.

4. Key Results

• Masses: The model predicts $M_{\Xi_b} = 5.795$ GeV and $M_{\Xi_c} = 2.468$ GeV, in excellent agreement with Particle Data Group (PDG) values.
• Form Factors:
  • The dominant form factors are $f_1$ and $g_1$, which are nearly identical numerically and increase with $q^2$.
  • Subleading form factors ($f_2, f_3, g_2, g_3$) are suppressed by $1/m_Q$ corrections but are non-zero.
  • Results at maximum and minimum recoil points agree well with previous theoretical studies (NRQM, Lattice QCD, RQM).
• Decay Widths:
  • $\Gamma(\Xi^0_b \to \Xi^+_c e \bar{\nu}_e) = 3.81^{+1.36}_{-1.17} \times 10^{10} \, \text{s}^{-1}$.
  • $\Gamma(\Xi^0_b \to \Xi^+_c \tau \bar{\nu}_\tau) = 1.24^{+0.17}_{-0.17} \times 10^{10} \, \text{s}^{-1}$.
  • These values lie at the lower end of the range predicted by other theoretical models but are consistent with relativistic quark models.
• Branching Ratios:
  • $\mathcal{B}(\Xi^0_b \to \Xi^+_c e \bar{\nu}_e) = 5.60^{+1.99}_{-1.72}\%$.
  • $\mathcal{B}(\Xi^0_b \to \Xi^+_c \tau \bar{\nu}_\tau) = 1.83^{+0.25}_{-0.25}\%$.
• LFU Ratio: The calculated ratio is $R(\Xi_c) \approx 0.325$, consistent with existing theoretical predictions (e.g., Ref [20]).

5. Significance

• Experimental Benchmark: With the LHCb and CMS collaborations observing excited bottom-strange states, precise theoretical predictions for $\Xi^0_b$ decays are timely. This work provides a benchmark for future experimental measurements of $\Xi^0_b \to \Xi^+_c \ell \bar{\nu}_\ell$.
• Standard Model Tests: The calculation of $R(\Xi_c)$ contributes to the broader effort to understand Lepton Flavour Universality. While no experimental data exists yet for this specific ratio, the theoretical prediction serves as a baseline for detecting potential New Physics deviations (similar to the $R(D^{(*)})$ anomaly).
• Methodological Validation: The agreement of the calculated slope and curvature parameters with other studies validates the HCQM+HQET approach for singly heavy baryons, suggesting it is a reliable tool for studying non-perturbative QCD dynamics in heavy baryon systems.
• Future Directions: The authors note that reducing the spread in theoretical predictions will require both precise experimental data and refined calculations (e.g., higher-order HQET corrections and Lattice QCD).