Analysis of the semileptonic decays of $Ξ_{cc}$ and $Ω_{cc}$ baryons in QCD sum rules

This paper employs three-point QCD sum rules to systematically analyze the spin-1/2 to spin-3/2 weak transition form factors of doubly charmed baryons, subsequently predicting the semileptonic decay rates of $\Xi_{cc}$ and $\Omega_{cc}$ to singly charmed baryons to aid in understanding heavy baryon dynamics and searching for new physics.

The Gist

Imagine the universe as a giant, bustling construction site. Most of the buildings we know (like protons and neutrons) are made of three bricks called quarks. Usually, these bricks come in pairs or triplets of light, fast-moving types (like up and down quarks).

But sometimes, nature builds something special: a "heavy-duty" building made of two heavy bricks (charm quarks) and one light one. These are called doubly charmed baryons (specifically $\Xi_{cc}$ and $\Omega_{cc}$). For a long time, these were just theoretical ghosts. It wasn't until 2017 that the LHCb experiment actually found one, like spotting a rare, mythical creature in the wild.

This paper is a detailed theoretical study of how these rare creatures decay (fall apart) and transform into other particles. Here is the breakdown of what the authors did, using simple analogies:

1. The Mission: Predicting the "Fall"

When a heavy particle decays, it doesn't just vanish; it changes into a lighter particle, a lepton (like an electron or muon), and a neutrino. This is called a semileptonic decay.

To predict exactly how fast this happens and how likely it is, physicists need to calculate something called Form Factors.

• The Analogy: Imagine you are trying to predict how a specific type of clay pot will shatter when dropped. You need to know the pot's internal structure, its density, and how the clay is bonded. The "Form Factors" are the mathematical blueprint of that internal structure. Without them, you can't predict the shatter pattern (the decay rate).

2. The Tool: QCD Sum Rules (The "X-Ray Machine")

The authors used a method called QCD Sum Rules.

• The Analogy: Quantum Chromodynamics (QCD) is the rulebook for how quarks stick together. But at low energies (inside a baryon), the rules get messy and "fuzzy" (non-perturbative). You can't just do simple math; it's like trying to calculate the exact path of a single drop of ink spreading in a glass of water.
• The Method: The authors used a "three-point" correlation function. Think of this as a three-way conversation between:
  1. The initial heavy baryon (the parent).
  2. The final lighter baryon (the child).
  3. The weak force current (the messenger causing the change).
• They calculated this conversation from two sides:
  • The Phenomenological Side: What we see in the real world (the particles).
  • The QCD Side: What's happening deep down with the quarks and gluons.
• By matching these two sides, they could extract the "Form Factors" (the blueprint).

3. The Challenge: The "Noise" Problem

The paper mentions a specific difficulty: Spin and Parity.

• The Analogy: Imagine you are trying to listen to a specific singer (the particle you want to study) in a crowded stadium. But there are other singers (excited states) and people shouting (noise) that sound similar.
• The authors had to filter out the "noise" from particles with the wrong spin (how they rotate) or parity (how they look in a mirror). They did this by setting up a system of 16 linear equations.
• Think of it like a complex puzzle where they had to solve for 16 different variables simultaneously to isolate the exact signal of the particle they cared about, ignoring the "ghosts" that looked similar but weren't the target.

4. The Calculation: From "Space" to "Time"

The math works best in "space-like" regions (a mathematical zone where the energy transfer is imaginary). But real decays happen in "time-like" regions (real energy).

• The Analogy: It's like measuring the shape of a shadow to figure out the shape of the object casting it. You measure the shadow (space-like), but you need to know the object (time-like).
• The authors calculated the numbers in the "shadow" zone and then used a fitting function (a mathematical curve) to stretch and extrapolate those numbers into the real world. It's like taking a photo of a shadow and using software to reconstruct the 3D object.

5. The Results: What Did They Find?

They analyzed four specific decay paths:

1. $\Xi_{cc}^{++} \to \Sigma_c^{*+}$
2. $\Xi_{cc}^{++} \to \Xi_c^{\prime *+}$
3. $\Omega_{cc}^{+} \to \Xi_c^{\prime *0}$
4. $\Omega_{cc}^{+} \to \Omega_c^{*0}$

Key Findings:

• Branching Fractions: They calculated the "branching ratio," which is the percentage of the time a particle chooses a specific decay path. They found that some paths (like $\Omega_{cc} \to \Omega_c^*$) are much more likely than others.
• Comparison: Their results were slightly different from other models (like the "Quark Model"), but when they applied their results to the final decay rates, the numbers were surprisingly close to other predictions. This suggests their "blueprint" is reliable.
• Symmetry Breaking: They found that the universe doesn't treat all quarks perfectly equally (SU(3) symmetry breaking). The mass of the strange quark makes a difference, causing the decay rates to deviate from perfect mathematical symmetry.

Why Does This Matter?

• Testing the Standard Model: By predicting these decays accurately, scientists can compare their predictions with future experiments. If the real world disagrees with the prediction, it might mean New Physics exists—something beyond our current understanding of the universe.
• Mapping the Heavy World: Since doubly heavy baryons are rare and hard to find, having accurate theoretical maps helps experimentalists know exactly where to look and what to expect when they finally catch one.

In Summary:
This paper is like a master architect drawing up the blueprints for how rare, heavy-duty particle buildings collapse. They used a sophisticated mathematical "X-ray" to see inside the particles, filtered out the noise of similar-looking particles, and predicted exactly how they will fall apart. This helps experimentalists at places like CERN know what to look for, potentially leading to discoveries that could rewrite the laws of physics.

Technical Summary

1. Problem Statement

The paper addresses the theoretical understanding of weak decay processes involving doubly charmed baryons ($\Xi_{cc}$ and $\Omega_{cc}$). Specifically, it focuses on the semileptonic decays where a spin-$1/2^+$ doubly charmed baryon transitions to a spin-$3/2^+$ singly charmed baryon ($1/2^+ \to 3/2^+$).

• Context: While single heavy baryons have been extensively studied, doubly heavy baryons are a newer frontier (with $\Xi_{cc}^{++}$ discovered by LHCb in 2017).
• Gap: Existing literature often focuses on $1/2^+ \to 1/2^+$ transitions or relies heavily on quark models which depend on specific spatial wavefunctions. There is a lack of systematic, non-perturbative calculations for $1/2^+ \to 3/2^+$ transitions involving doubly heavy baryons.
• Goal: To calculate the transition form factors for $\Xi_{cc}^{++} \to \Sigma_c^{*+}$, $\Xi_{cc}^{++} \to \Xi_c^{\prime *+}$, $\Omega_{cc}^{+} \to \Xi_c^{\prime *0}$, and $\Omega_{cc}^{+} \to \Omega_c^{*0}$, and subsequently predict their semileptonic decay widths and branching ratios.

2. Methodology

The authors employ the Three-Point QCD Sum Rules (QCDSR) framework, a non-perturbative method that bridges QCD and hadronic physics via the Operator Product Expansion (OPE).

A. Correlation Function and Interpolating Currents

• A three-point correlation function is constructed involving the initial baryon current ($J_{B_1}$), the final baryon current ($J_{B_2}^\mu$), and the weak transition current ($J_\nu^{V-A}$).
• Currents: Specific interpolating currents are defined for the doubly charmed baryons ($\Xi_{cc}, \Omega_{cc}$) and the singly charmed spin-3/2 baryons ($\Sigma_c^*, \Xi_c^{\prime *}, \Omega_c^*$).
• Phenomenological Side: The correlation function is saturated with a complete set of hadronic states. Crucially, the authors consider all possible couplings of the interpolating currents to hadronic states, including:
  • Ground states ($1/2^+, 3/2^+$).
  • Excited/negative parity states ($1/2^-, 3/2^-$).
  • Noise Elimination: To isolate the desired $1/2^+ \to 3/2^+$ transition, the authors systematically eliminate Dirac structures proportional to $\gamma_\mu$ and $p'_\mu$ which correspond to contaminating $1/2^\pm$ and $3/2^-$ states. This results in a system of 16 linear equations derived from 16 distinct Dirac structures.

B. QCD Side and OPE

• The correlation function is calculated at the quark-gluon level using the OPE.
• Contributions Included: The calculation goes beyond leading order, incorporating:
  • Perturbative term (dimension 0).
  • Quark condensate $\langle \bar{q}q \rangle$ (dimension 3).
  • Gluon condensate $\langle g_s^2 GG \rangle$ (dimension 4).
  • Mixed quark-gluon condensate $\langle \bar{q}g_s \sigma G q \rangle$ (dimension 5).
  • Four-quark condensate $g_s^2 \langle \bar{q}q \rangle^2$ (dimension 6).
• Spectral Density: The spectral densities are derived using Cutkosky's rules, handling the integration of delta functions for on-shell quarks.

C. Extraction of Form Factors

• Borel Transformation: Double Borel transformations are applied to both sides to suppress higher excited states and continuum contributions.
• Stability: The "Borel platform" is determined by ensuring pole dominance (>40%) and OPE convergence (high-dimension condensates < 2%).
• Extrapolation: Form factors are initially calculated in the space-like region ($Q^2 > 0$). To access the physical time-like region ($Q^2 < 0$) required for decay rates, the results are fitted using a $z$-series expansion (a model-independent parameterization) and extrapolated.

3. Key Contributions

1. Systematic Treatment of Spin-3/2 Transitions: Unlike previous works that often simplified the spin structure, this paper rigorously handles the $1/2^+ \to 3/2^+$ transition by solving a coupled system of 16 linear equations to disentangle the 8 form factors ($F_{1-4}, G_{1-4}$) from the complex Dirac structure of the correlation function.
2. High-Dimensional Condensates: The inclusion of dimension-6 condensates ($g_s^2 \langle \bar{q}q \rangle^2$) and dimension-5 mixed condensates improves the precision and reliability of the sum rules compared to calculations truncating at lower dimensions.
3. Comprehensive Decay Analysis: The study provides the first QCDSR predictions for the semileptonic decays of $\Omega_{cc}$ and $\Xi_{cc}$ into spin-3/2 states, covering both electron and muon channels.
4. SU(3) Symmetry Breaking Analysis: The authors explicitly analyze the deviation from flavor SU(3) symmetry relations in the decay widths, attributing the differences to the non-negligible mass of the strange quark compared to up/down quarks.

4. Key Results

• Form Factors: The momentum-dependent form factors $F_i(Q^2)$ and $G_i(Q^2)$ were successfully extracted.
  • The authors note that their values at $Q^2=0$ differ from Light-Front Quark Model (LFQM) predictions, but the discrepancy diminishes after extrapolation to the physical region ($q^2_{max}$).
  • The fitted form factors show that the absolute values generally increase when extrapolated from space-like to time-like regions (contrary to some LFQM results where they decrease).
• Decay Widths and Branching Ratios:
  • Calculated partial decay widths ($\Gamma$) and branching fractions ($\mathcal{B}$) for $l = e, \mu$.
  • Dominant Channels: The decays $\Xi_{cc}^{++} \to \Xi_c^{\prime *+} l^+ \nu_l$ and $\Omega_{cc}^{+} \to \Omega_c^{*0} l^+ \nu_l$ are predicted to have significantly larger branching fractions (~3.7% and ~4.7% respectively) compared to the $\Sigma_c^*$ and $\Xi_c^{\prime *}$ modes involving $\Omega_{cc}$ or $\Xi_{cc}$ to $\Sigma_c^*$.
  • Comparison: The predicted widths are slightly lower than those from covariant confined quark models (Refs [14, 23]) but are of the same order of magnitude ($10^{-14}$ GeV).
  • Longitudinal vs. Transverse: The ratio of longitudinal to transverse decay widths ($\Gamma_L/\Gamma_T$) is found to be approximately 1.5–1.65.

5. Significance

• Experimental Guidance: The predictions provide crucial input for experimentalists (e.g., at LHCb, Belle II, and future facilities) searching for doubly heavy baryons. The identified dominant decay channels offer the most promising avenues for reconstruction.
• New Physics Potential: Precise knowledge of these form factors and decay rates is essential for extracting CKM matrix elements ($V_{cs}, V_{cd}$) and searching for deviations that could signal physics beyond the Standard Model.
• Theoretical Validation: The work demonstrates the robustness of QCDSR in handling complex spin transitions and high-dimensional condensates, offering a complementary perspective to quark models and lattice QCD (which is currently challenging for these specific heavy-light-heavy systems).
• Future Directions: The authors highlight that their framework can be extended to heavy-to-heavy transitions and tensor form factors (relevant for radiative and dilepton decays), paving the way for a more complete understanding of heavy baryon dynamics.