Sensitivity study of $\bar{K}_1(1270)$ decay dynamics using four $D\to \bar{K}_1(1270)(\to \bar K\pi\pi)e^+\nu$ decay channels

This paper presents a sensitivity study using the BESIII experiment to propose a model-independent method for determining the branching fraction of $\bar{K}_1(1270)$ decays into three-body final states by simultaneously analyzing four specific semileptonic $D$-meson decay channels.

The Gist

Imagine the subatomic world as a giant, chaotic dance floor. In this dance, particles called mesons are the dancers, constantly spinning, colliding, and breaking apart into smaller groups.

This paper is about a specific, somewhat mysterious dancer named $\bar{K}_1(1270)$. It's a "strange" particle (it contains a strange quark) that is very unstable. It doesn't stay together for long; it almost immediately splits into three smaller particles (a kaon and two pions).

The problem is, we don't know exactly how it splits. Does it break apart in one specific way 70% of the time and another way 30%? Or is it a 50/50 split? Currently, our best guesses are fuzzy, like trying to guess the ingredients of a soup by tasting a spoonful from a very large, murky pot.

Here is a simple breakdown of what the scientists in this paper are doing:

1. The Mystery of the "Mixing Angle"

Think of the $\bar{K}_1(1270)$ particle not as a single, solid object, but as a smoothie made of two different "flavors" of particles mixed together.

• In physics, these flavors are called quantum states.
• The "recipe" for this smoothie is determined by something called a mixing angle.
• If we get the recipe wrong, our calculations for how the universe works (specifically a force called the "Strong Force") will be off. To fix the recipe, we need to know exactly how often this particle breaks into specific combinations of smaller particles.

2. The Old Way vs. The New Way

The Old Way (The "Guessing Game"):
Previously, scientists tried to figure out the recipe by looking at just one type of breakup (like $K^- \pi^+ \pi^-$). They had to make a lot of assumptions about the "middle steps" of the decay. It was like trying to figure out a magician's trick by only watching the final reveal, assuming you knew how the cards were shuffled. This led to big uncertainties (about 20% error).

The New Way (The "Simultaneous Audit"):
The authors of this paper propose a smarter, model-independent method. Instead of guessing the middle steps, they look at four different breakup channels at the exact same time:

1. Kaon + 2 charged pions
2. Kaon + 1 charged pion + 1 neutral pion
3. Neutral Kaon + 2 charged pions
4. Neutral Kaon + 1 charged pion + 1 neutral pion

The Analogy:
Imagine you are a detective trying to figure out how much money was stolen from a bank vault.

• Old Method: You look at one security camera angle, but the camera is blurry, so you have to guess how much was taken.
• New Method: You have four different cameras covering the same vault from different angles. You don't need to know how the thief moved (the "model"); you just count the footprints in all four cameras simultaneously. Because the cameras see different angles, the math balances out perfectly, and you can calculate the exact amount stolen without guessing the thief's technique.

3. The "Missing Mass" Trick

To do this, they use data from the BESIII experiment in China. They look at D-mesons (another type of particle) that decay into our mystery particle ($\bar{K}_1$) plus an electron and a neutrino.

The neutrino is a ghost; it passes through everything and can't be seen. However, the scientists use a clever trick called Missing Mass.

• They know exactly how much energy and momentum the D-meson started with.
• They measure everything that is seen (the electron, the pions, the kaon).
• Whatever energy is "missing" must belong to the invisible neutrino.
• By analyzing the "missing" energy, they can reconstruct the properties of the $\bar{K}_1$ particle that decayed.

4. The Simulation (The "Virtual Reality" Test)

Before running the experiment on real data, the scientists ran 2,000 computer simulations (called pseudo-experiments).

• They created a "fake universe" with known rules.
• They applied their new "four-camera" method to this fake data.
• Result: The method worked perfectly. It recovered the correct answers with very little error.

5. Why Does This Matter?

• Precision: This new method reduces the uncertainty from ~20% down to about 5%. That is a massive improvement.
• Model Independence: They don't need to assume how the particle behaves inside; the math does the work for them.
• Future Proofing: While this study uses current data, the method is ready for the Super Tau-Charm Factory (a future, even bigger particle collider). With more data, the precision could get even better, helping us understand the fundamental building blocks of the universe.

The Bottom Line

This paper is a proposal for a better accounting method for particle physics. Instead of guessing how a particle breaks apart based on one clue, the scientists suggest looking at four different clues simultaneously. By cross-referencing them, they can determine the particle's true nature with much higher confidence, helping us solve the mystery of how the "Strong Force" holds our universe together.

Technical Summary

1. Problem Statement

The strange axial-vector mesons $\bar{K}_1(1270)$ and $\bar{K}_1(1400)$ are crucial for understanding non-perturbative Quantum Chromodynamics (QCD) and $SU(3)$ symmetry breaking. Their properties, particularly the mixing angle $\theta_{\bar{K}_1}$, are difficult to determine precisely due to large uncertainties in their branching fractions (BFs) into three-body final states ($\bar{K}\pi\pi$).

• Current Limitations: Existing measurements of $\bar{K}_1(1270)$ decay BFs rely heavily on a 1981 scattering experiment and recent BESIII measurements that assume specific intermediate resonant contributions (e.g., $\bar{K}^*\pi$, $\bar{K}\rho$).
• The Bottleneck: The uncertainties in these assumed BFs propagate to $\sim 20\%$ uncertainty in the total $\bar{K}_1(1270) \to \bar{K}\pi\pi$ branching fraction. This hinders precise measurements of semileptonic $D$-meson decays and limits the ability to determine the photon polarization in $b \to s\gamma$ transitions.
• The Gap: Previous BESIII analyses of $D \to \bar{K}\pi\pi e^+\nu_e$ decays were performed separately for different charge states or relied on amplitude analyses that require detailed modeling of intermediate resonances, introducing model dependence.

2. Methodology

The authors propose a model-independent approach to determine the branching fractions of $\bar{K}_1(1270)$ decays by performing a simultaneous analysis of four specific semileptonic $D$-meson decay channels using data from the BESIII experiment (specifically the $\psi(3770)$ dataset with an integrated luminosity of $20.3 \text{ fb}^{-1}$).

Key Technical Steps:

• Four Decay Channels: The study simultaneously fits signal yields from:
  1. $D^0 \to K^-\pi^+\pi^- e^+\nu_e$
  2. $D^+ \to K^-\pi^+\pi^0 e^+\nu_e$
  3. $D^0 \to K_S^0\pi^-\pi^0 e^+\nu_e$
  4. $D^+ \to K_S^0\pi^+\pi^- e^+\nu_e$
• Variable Definition: Instead of relying on complex amplitude models, the authors define a transition variable $\beta$ based on the ratio of branching fractions:
  $$\beta^{-1} \equiv 1 - \frac{\mathcal{B}(\bar{K}_1^-(1270) \to K^-\pi^+\pi^-)}{\mathcal{B}(\bar{K}_1^0(1270) \to K^-\pi^+\pi^0)}$$
• Model Independence: By expressing the total three-body branching fraction ($\mathcal{B}_{3body}$) and the ratio of specific intermediate modes ($\alpha$) as explicit functions of $\beta$ and external inputs (PDG values for $\bar{K}\omega$ and $\bar{K}f_0$), the method eliminates the need to know the detailed dynamics of intermediate resonances.
• Statistical Analysis:
  • Pseudo-experiments: 2000 pseudo-datasets were generated based on the $M_{miss}^2$ (missing mass squared) distributions expected from the $20.3 \text{ fb}^{-1}$ dataset.
  • Fitting: An unbinned maximum likelihood fit was performed simultaneously across all four channels.
  • Systematic Control: The double-tag method and the use of yield ratios ($\beta$) cancel out common systematic uncertainties (luminosity, tagging efficiency, tracking, PID). Residual systematics (reconstruction efficiency for $\pi^0/K_S^0$ and external BF inputs) were constrained using Gaussian priors.

3. Key Contributions

• Model-Independent Framework: The paper introduces a formalism where the branching fractions of $\bar{K}_1(1270)$ are derived solely from the experimentally determined ratio $\beta$, removing dependence on specific resonance models (like the relative phases or magnitudes of $\bar{K}^*\pi$ vs. $\bar{K}\rho$).
• Simultaneous Multi-Channel Analysis: It demonstrates the feasibility of extracting physics parameters by combining four distinct charge-conjugate decay modes, maximizing statistical power and constraining systematic errors through cross-channel consistency.
• Sensitivity Projection: The study provides a rigorous sensitivity projection for the upcoming BESIII dataset, showing that current data is sufficient to significantly improve precision over previous separate analyses.

4. Results

Based on the sensitivity study using the $20.3 \text{ fb}^{-1}$ dataset:

• Branching Fraction Precision: The method achieves a total systematic uncertainty of approximately 5% for the determination of $\mathcal{B}(\bar{K}_1(1270) \to \bar{K}\pi\pi)$. This is a significant improvement over the $\sim 20\%$ uncertainty in previous inputs.
• Ratio $\alpha$: The ratio $\alpha \equiv \mathcal{B}(\bar{K}_1(1270) \to \bar{K}^*\pi) / \mathcal{B}(\bar{K}_1(1270) \to \bar{K}\rho)$ can be determined with a statistical uncertainty of $\sim 15\%$ (dominated by current statistics), with systematic uncertainties reduced to $\sim 1\%$.
• Comparison with Amplitude Analysis: While the statistical uncertainty on $\alpha$ is larger than that of a full amplitude analysis (which utilizes angular distributions), the systematic uncertainty is significantly lower because it avoids model bias.
• Validation: The "pull" distribution of the fitted parameter $\alpha$ across 2000 pseudo-experiments is consistent with a normal distribution centered at zero, confirming the fit model is unbiased.
• Specific Outputs (Projected):
  • $\mathcal{B}(\bar{K}_1(1270) \to \bar{K}\rho)$: $\sim 72.5 \pm 9.0\%$ (stat) $\pm 0.7\%$ (syst).
  • $\mathcal{B}(\bar{K}_1(1270) \to \bar{K}^*\pi)$: $\sim 16.5 \pm 9.0\%$ (stat) $\pm 0.7\%$ (syst).
  • The uncertainties on the input $\mathcal{B}(D \to \bar{K}_1(1270)e^+\nu_e)$ are considerably reduced compared to previous measurements.

5. Significance

• Fundamental Physics: This approach provides a cleaner path to determining the mixing angle $\theta_{\bar{K}_1}$ and the masses/widths of axial-vector mesons, which are currently poorly constrained.
• Standard Model Tests: Precise knowledge of these branching fractions is essential for reducing theoretical ambiguities in determining photon polarization in $b \to s\gamma$ transitions, a key probe for New Physics.
• Future Prospects: The method is scalable. With the anticipated larger datasets from the Super Tau-Charm Factory, statistical uncertainties are expected to drop by an order of magnitude, making this model-independent technique the gold standard for precision spectroscopy of strange axial-vector mesons.
• Validation of PDG: It offers a robust, data-driven validation of the Particle Data Group's averages, which currently rely on older, less precise scattering data.