$\Xi_c \to \Xi$ form factors from lattice QCD with domain-wall quarks: A new piece in the puzzle of $\Xi_c^0$ decay rates

Using lattice QCD with domain-wall and anisotropic clover quark actions, this study calculates the form factors for $\Xi_c \to \Xi$ semileptonic decays and predicts decay rates that are higher than both previous lattice results and experimental measurements, yet consistent with approximate $SU(3)$ flavor symmetry expectations.

The Gist

Imagine the universe is a giant, complex machine made of tiny building blocks called quarks. These quarks stick together to form larger particles, like protons and neutrons, but also heavier, more exotic ones called baryons.

This paper is about a specific, heavy baryon called the $\Xi_c$ (pronounced "Xi-c"). Think of the $\Xi_c$ as a heavy, unstable "cosmic firework." It doesn't last long; it quickly explodes (decays) into lighter particles. One of the ways it explodes is by shooting out a "lepton" (like an electron or a muon) and a ghostly particle called a neutrino. This is called a semileptonic decay.

Scientists want to predict exactly how often this happens and how fast the firework explodes. But there's a problem: the numbers they get from their best theories don't match the numbers they see in real experiments. It's like a weather forecast predicting a sunny day, but it's pouring rain outside.

Here is what the authors of this paper did to solve the mystery, explained simply:

1. The Problem: A Mismatched Puzzle

For a long time, scientists have been trying to figure out the "decay rate" of the $\Xi_c$.

• The Theory: Based on the rules of the Standard Model (the rulebook of particle physics), the $\Xi_c$ should decay at a certain speed.
• The Experiment: When real-life detectors (like those at the Belle and ALICE labs) watch these particles, they see them decaying much slower than the theory predicts.
• The Confusion: A previous computer simulation (a "lattice QCD" calculation) tried to fix this but still got numbers that were too low. The authors of this paper decided to build a better simulation to see if they could finally solve the puzzle.

2. The Method: Building a Digital Universe

To study these particles, you can't just put them in a test tube; they are too small and too fast. Instead, scientists use Lattice QCD.

• The Analogy: Imagine trying to study how water flows. You can't see individual water molecules easily, so you build a giant 3D grid (a lattice) and simulate the water moving through the grid squares.
• The Upgrade: The authors built a much more detailed grid than before.
  • They used four different versions of this grid, some with "coarse" squares and some with very "fine" squares, to make sure their results weren't just an artifact of the grid size.
  • They used a special type of math called Domain-Wall fermions for the lighter quarks (up, down, strange). Think of this as using a high-definition camera for the background scenery.
  • For the heavy "charm" quark, they used a specialized anisotropic clover action. This is like using a heavy-duty, slow-motion camera specifically for the main actor in the scene, ensuring the heavy particle is tracked perfectly without blurring.

3. The Calculation: Watching the "Firework" Explode

The team ran massive supercomputer simulations to calculate something called form factors.

• The Metaphor: Imagine the $\Xi_c$ is a spinning top. As it decays, it changes shape and speed. The "form factors" are like a detailed map of exactly how the top spins, wobbles, and changes shape during that split second before it disappears.
• They calculated these maps for many different scenarios (different energies, different grid sizes) and then used a clever statistical trick called Bayesian Model Averaging.
  • The Analogy: Imagine asking 100 experts to guess the height of a building. Instead of just taking the average, you ask them to guess based on different assumptions (e.g., "What if the ground is soft?" "What if the wind is strong?"). Then, you weigh their answers based on how likely those assumptions are. This gives a much more reliable final answer than just a simple average.

4. The Results: The Theory is Stronger Than We Thought

After crunching the numbers, the authors found something surprising:

• The Prediction: Their new, high-precision simulation predicts that the $\Xi_c$ should decay much faster than the previous computer simulations suggested. In fact, their prediction is about 3.5 times higher than what the current experimental data shows.
• The Comparison:
  • Old Theory: "It decays slowly."
  • New Theory (This Paper): "It should decay very fast!"
  • Current Experiment: "We see it decaying slowly."

5. The Big Mystery Remains (But We Have a Clue)

So, did they solve the puzzle? Not quite, but they narrowed it down.

• Their new calculation is very precise and consistent with a fundamental symmetry of nature called SU(3) flavor symmetry (which is like saying "if you swap ingredients in a recipe, the taste should be similar").
• The fact that their new, more accurate theory still disagrees with the experiment suggests the problem might not be with the theory. The problem might be with the experiment.

The "Smoking Gun" Clue:
The experiment measures the $\Xi_c$ decay by comparing it to another decay (a "normalization mode"). The authors suggest that the measurement of this other decay might be wrong (perhaps underestimated). If the "normalization" measurement is actually higher than we think, then the $\Xi_c$ decay rate would look higher too, and the mystery would be solved!

Summary

This paper is like a team of master architects building a perfect digital model of a collapsing building.

1. They built a better model than anyone before.
2. Their model says the building should collapse very quickly.
3. Real-life observers say the building is collapsing slowly.
4. The authors conclude: "Our model is solid. The observers might be measuring the 'slow collapse' of a different building incorrectly. Let's re-measure that other building!"

This work doesn't just give a number; it gives physicists a new direction to look. It suggests that the "missing piece" of the puzzle isn't in our understanding of the laws of physics, but likely in how we are measuring the data in the real world.

Technical Summary

1. Problem Statement

The paper addresses a significant discrepancy between experimental measurements and theoretical predictions regarding the semileptonic decay rates of the charmed baryon $\Xi_c$ into the $\Xi$ baryon ($\Xi_c \to \Xi \ell^+ \nu_\ell$).

• The Discrepancy: Experimental measurements of the branching fraction $B(\Xi^0_c \to \Xi^- e^+ \nu_e)$ yield values around $1.05\%$ (PDG average), which is significantly lower than predictions based on approximate $SU(3)$ flavor symmetry (expected $\approx 4.0-5.0\%$) and previous theoretical models.
• Previous Theory: A 2021 lattice QCD calculation by Zhang et al. predicted a branching fraction of $2.38\%$, which was higher than experiment but still lower than the $SU(3)$ symmetry expectation. However, that calculation lacked a chiral extrapolation (using pion masses $\sim 290-300$ MeV) and used a different fermion action (clover for all quarks).
• Goal: To resolve this tension by performing a high-precision lattice QCD calculation of the $\Xi_c \to \Xi$ form factors using physical quark masses, multiple lattice spacings, and improved actions to determine Standard Model predictions for the decay rates.

2. Methodology

The authors performed a state-of-the-art lattice QCD simulation with the following technical specifications:

• Ensembles: Four gauge-field ensembles generated by the RBC and UKQCD collaborations, featuring $2+1$ flavors of domain-wall fermions (DWF) for up, down, and strange quarks.
  • Lattice spacings ($a$): Ranging from $0.111$ fm to $0.073$ fm.
  • Pion masses ($m_\pi$): Ranging from $420$ MeV down to $230$ MeV (approaching the physical point).
• Charm Quark Action: An anisotropic clover action was used for the valence charm quark. The parameters (mass, anisotropy, and clover coefficients) were tuned non-perturbatively to reproduce experimental $D_s$ meson properties (rest mass, kinetic mass, and hyperfine splitting).
• Currents and Renormalization:
  • The $c \to s$ vector and axial-vector currents were renormalized using a mostly non-perturbative scheme.
  • Matching factors ($Z_V$) were calculated non-perturbatively.
  • $O(a)$ improvement terms were included, with perturbative coefficients estimated via extrapolation for the finest ensemble.
• Correlation Functions:
  • Calculated two-point and three-point correlation functions for $\Xi_c$ and $\Xi$ baryons.
  • Used 15 different source-sink separations to control excited-state contamination.
  • Employed Bayesian Model Averaging (BMA) using the "perfect model" Akaike Information Criterion (AIC) to extrapolate results to infinite source-sink separation, rather than relying on simple fit-range shifts.
• Extrapolation to Physical Limit:
  • Performed a combined chiral and continuum extrapolation using a BCL $z$-expansion.
  • The expansion included terms for pion mass dependence ($m_\pi$) and lattice spacing dependence ($a$), constrained by heavy-quark symmetry and $SU(3)$ flavor relations.
  • Systematic uncertainties (renormalization, finite volume, isospin breaking, discretization) were incorporated directly into the fit covariance matrix.

3. Key Contributions

• First High-Precision Lattice Calculation with DWF: This is the first lattice calculation of $\Xi_c \to \Xi$ form factors using domain-wall fermions for light/strange quarks, ensuring excellent chiral symmetry properties.
• Non-Perturbative Charm Tuning: The use of non-perturbative tuning for the heavy-quark action parameters reduces systematic errors associated with heavy-quark discretization.
• Advanced Statistical Analysis: The application of Bayesian Model Averaging for excited-state removal provides a more robust and statistically rigorous determination of the ground-state matrix elements compared to previous methods.
• Comprehensive Error Budget: The study includes a detailed treatment of systematic uncertainties, including a conservative $1\%$ uncertainty for missing higher-order corrections in renormalization factors.

4. Key Results

The authors obtained precise Standard Model predictions for the differential and integrated decay rates.

• Form Factors: The calculated vector ($f_+, f_\perp, f_0$) and axial-vector ($g_+, g_\perp, g_0$) form factors are higher than those from the previous lattice calculation (Zhang et al.) and are consistent with expectations from approximate $SU(3)$ flavor symmetry.
• Decay Rates (Normalized by $|V_{cs}|^2$):
  • $\Gamma(\Xi^0_c \to \Xi^- e^+ \nu_e) / |V_{cs}|^2 = 0.2515(73) \text{ ps}^{-1}$
  • $\Gamma(\Xi^0_c \to \Xi^- \mu^+ \nu_\mu) / |V_{cs}|^2 = 0.2437(71) \text{ ps}^{-1}$
• Branching Fractions:
  • $B(\Xi^0_c \to \Xi^- e^+ \nu_e) = \mathbf{3.58(12)\%}$
  • $B(\Xi^0_c \to \Xi^- \mu^+ \nu_\mu) = 3.47(12)\%$
  • $B(\Xi^+_c \to \Xi^0 e^+ \nu_e) = 10.94(34)\%$
• Comparison with Experiment:
  • The predicted branching fraction for $\Xi^0_c \to \Xi^- e^+ \nu_e$ ($3.58\%$) is approximately 3.4 times larger than the current PDG experimental average ($1.05\%$).
  • The discrepancy has a statistical significance of $10.8\sigma$.

5. Significance and Implications

• Confirmation of SU(3) Symmetry: The results strongly support the approximate $SU(3)$ flavor symmetry prediction ($\approx 4\%$), suggesting that the large breaking previously hypothesized to explain the low experimental values is not present in the form factors.
• The "Puzzle" Deepens: The lattice calculation rules out theoretical deficiencies (such as lack of chiral extrapolation or excited-state contamination) as the primary cause of the discrepancy. The tension now lies squarely between the robust Standard Model prediction and the experimental data.
• Potential Experimental Issues: The authors suggest that the experimental measurement of the normalization mode, $B(\Xi^0_c \to \Xi^- \pi^+)$, may be underestimated. Recent global analyses of charm-baryon nonleptonic decays using $SU(3)$ symmetry support the possibility that the normalization mode branching fraction is larger than currently measured.
• Future Directions: The paper calls for new, precise measurements of the $\Xi^0_c \to \Xi^- \pi^+$ normalization mode to resolve the puzzle. If the normalization mode is confirmed to be larger, the extracted semileptonic branching fraction would increase, potentially reconciling experiment with the lattice QCD prediction.

In summary, this work provides a definitive Standard Model calculation that challenges current experimental values, shifting the focus of the "puzzle" from theoretical limitations to potential experimental systematic errors in the normalization channels.