Improved branching-fraction measurements of $B^0_{(s)} \to K_S^0 h^+ h^{'-}$ decays and first observation of $B^0_(s) \to K_S^0 K^+ K^-$

Using 9 fb$^{-1}$ of LHCb $pp$ collision data, this study reports the first observation of the $B^0_s \to K^0_S K^+ K^-$ decay and provides improved measurements of branching-fraction ratios for various charmless three-body $B^0_{(s)} \to K^0_S h^+ h^{\prime -}$ decays.

The Gist

Imagine the Large Hadron Collider (LHC) at CERN as a giant, high-speed particle racetrack. Scientists smash protons together at nearly the speed of light, creating a chaotic explosion of new particles. Among the debris, they are looking for very specific, rare "ghosts" called B mesons.

This paper is like a detective report from the LHCb team, a group of scientists who specialize in hunting these ghosts. They have analyzed a massive amount of data (equivalent to 9 years of heavy traffic on the racetrack) to study how these B mesons fall apart.

Here is the story of their discovery, explained simply:

1. The Mystery: How Do the Ghosts Break Up?

When a B meson dies, it doesn't just vanish; it breaks apart into three smaller particles. The scientists are interested in a specific family of breakups where the B meson turns into a neutral Kaon ($K^0_S$) and two other charged particles (which can be pions or kaons).

Think of the B meson as a fragile glass vase. When it shatters, it usually breaks into predictable pieces. However, the laws of physics (specifically the "Standard Model") predict exactly how often this should happen. If the vase breaks in a way that is slightly different from the prediction, it might mean there is a "ghost" (a new, undiscovered particle or force) helping to break it.

2. The New Discovery: Finding the "Missing" Breakup

For a long time, the scientists had a hunch that one specific type of breakup—where the B meson turns into a neutral Kaon and two positive/negative Kaons ($B^0_s \to K^0_S K^+ K^-$)—was happening, but they couldn't prove it. It was like hearing a faint whisper in a noisy stadium; they knew it was there, but they couldn't isolate the voice.

The Big News: With this new, larger dataset, the whisper became a shout. They have officially observed this decay for the first time! It's like finally finding that missing piece of a puzzle that everyone thought was lost.

3. The Method: Counting the Pieces

To measure how often these breakups happen, the scientists used a clever trick. Instead of trying to count every single B meson that ever existed (which is impossible), they used a reference point.

• The Reference: They picked a very common breakup ($B^0 \to K^0_S \pi^+ \pi^-$) as their "standard ruler."
• The Ratio: They counted how many times the rare breakups happened compared to the common one.
  • Analogy: Imagine you are counting cars in a parking lot. You know there are exactly 1,000 red sedans (the common breakup). You count the blue convertibles (the rare breakup) and find there are 578 of them. You don't need to know the total number of cars in the universe; you just know the ratio is 0.578. This makes the measurement much more precise.

4. The Tools: The Digital Net

The LHCb detector is like a giant, high-speed camera that takes millions of pictures of these collisions. But most of the pictures are just background noise (other particles).

• The Filter (Trigger): The computer system acts like a bouncer at a club, instantly deciding which events are interesting enough to save.
• The Detective Work (Selection): The scientists then use a "digital net" (a machine learning algorithm called a BDT) to sift through the saved events. They look for specific patterns, like the flight path of the particles, to separate the real signal from the background noise.
• The Mass Fit: Finally, they plot the "mass" of the particles. If the B meson is there, it shows up as a sharp peak on a graph, rising above a flat line of background noise. In this paper, the peak for the new discovery was so clear it stood out with a significance of 10 standard deviations (in science, 5 is usually enough to claim a discovery; 10 is a slam dunk).

5. The Results: What Did They Learn?

The team measured the "branching fractions" (the probability of these breakups) for several different combinations. Here are the key takeaways:

• First Observation: They confirmed the existence of $B^0_s \to K^0_S K^+ K^-$.
• Precision: They measured the ratios of these breakups with incredible precision. For example, the $B^0_s$ meson is about 1.8 times more likely to break into a Kaon and a Pion than the $B^0$ meson is.
• No New Physics (Yet): The results match the predictions of the Standard Model very well. This is actually good news! It means our current understanding of the universe is solid. However, it also means we haven't found the "new physics" (like new particles) yet, so the hunt continues.

6. Why Does This Matter?

You might ask, "Why do we care about how a B meson breaks apart?"

• Testing the Rules: Every time we measure these decays precisely, we are stress-testing the laws of physics. If the numbers were slightly off, it would be a crack in the foundation of our understanding, pointing us toward new theories.
• Understanding the Universe: These decays involve "loop transitions," which are like quantum shortcuts where particles briefly pop in and out of existence. Studying them helps us understand the fundamental forces that hold the universe together.
• The Future: Now that they have observed this new decay, they can use it to study CP violation (a subtle difference between matter and antimatter). This is crucial for understanding why the universe is made of matter and not just empty space.

Summary

In short, the LHCb team used a massive amount of data to catch a rare particle breakup that had been hiding in the shadows. They proved it exists, measured exactly how often it happens, and confirmed that the universe is behaving exactly as our best theories predict. It's a victory for precision, a step forward in our understanding of the cosmos, and a reminder that even in a world of tiny particles, there are still new stories to be told.

Technical Summary

Here is a detailed technical summary of the paper CERN-EP-2026-029 / LHCb-PAPER-2024-029, titled "Improved branching-fraction measurements of $B^0_{(s)} \to K^0_S h^+ h'^-$ decays and first observation of $B^0_s \to K^0_S K^+ K^-$".

1. Problem and Motivation

The paper addresses the study of charmless three-body decays of neutral $B$ mesons ($B^0$ and $B^0_s$) into final states containing a short-lived neutral kaon ($K^0_S$) and two charged hadrons ($h^+ h'^-$), where $h, h' \in \{\pi, K\}$.

• Theoretical Context: These decays are dominated by loop transitions ($b \to q\bar{q}s$) in the Standard Model (SM). Precise measurements of their branching fractions and CP-violating asymmetries are crucial for testing the SM, constraining New Physics (NP) models, and refining theoretical approaches like QCD factorization.
• Previous Limitations: Previous LHCb studies (using 2011–2012 data, 3.0 fb$^{-1}$) observed $B^0_s \to K^0_S \pi^+ \pi^-$ and $B^0_s \to K^0_S K^\pm \pi^\mp$ but found no conclusive evidence for the rare decay $B^0_s \to K^0_S K^+ K^-$ (significance of only 2.5$\sigma$).
• Goal: To improve the precision of branching fraction ratios and achieve the first definitive observation of the $B^0_s \to K^0_S K^+ K^-$ decay using a significantly larger dataset.

2. Methodology

Data Sample

The analysis utilizes proton-proton collision data collected by the LHCb experiment between 2011 and 2018, corresponding to an integrated luminosity of 9 fb$^{-1}$ at center-of-mass energies of 7, 8, and 13 TeV. This dataset is approximately three times larger than that used in the previous LHCb study.

Event Selection and Reconstruction

• Trigger: A two-stage trigger system (hardware and software) selects events with displaced secondary vertices consistent with $b$-hadron decays.
• $K^0_S$ Reconstruction: Two categories are used based on the decay length of the $K^0_S$:
  • Long (LL): Decays within the Vertex Locator (VELO).
  • Downstream (DD): Decays outside the VELO, reconstructed using track segments in the downstream tracker.
• Particle Identification (PID): Ring-imaging Cherenkov (RICH) detectors distinguish between pions and kaons.
• Multivariate Analysis:
  • Topological BDT: Suppresses combinatorial background using decay topology variables (flight distance, impact parameters, etc.).
  • PID BDT: Suppresses "cross-feed" backgrounds where a signal from one channel is misidentified as another (e.g., a pion misidentified as a kaon).
• Vetoing: Decays involving intermediate charm or charmonium resonances (e.g., $D^0$, $J/\psi$) are explicitly reconstructed and vetoed to prevent contamination.

Fit Model and Yield Extraction

• Mass Fits: Unbinned maximum-likelihood fits are performed on the invariant mass spectra ($m_{K^0_S h^+ h'^-}$) in the range 5100–5750 MeV/$c^2$.
• Signal Modeling: Signal peaks are modeled using a sum of two Crystal Ball functions (sharing peak position and width but with opposite tails).
• Background Modeling:
  • Cross-feed: Modeled with the same shape as the signal but constrained by misidentification probabilities derived from simulation.
  • Partially Reconstructed: Modeled using a convolution of a generalized ARGUS function and a Gaussian, based on specific decay modes (e.g., $B^0 \to K^*(892)^0 \rho(770)^0$ with a missing $\pi^0$).
  • Combinatorial: Modeled with a first-order polynomial.
• Efficiency Correction: Due to the non-uniform distribution of events in the Dalitz plot (phase space), efficiencies are calculated as weighted averages. The analysis uses a square Dalitz plot binned in a $10 \times 10$grid. Weights are derived from data using the **sPlot** technique to isolate signal events. Corrections for data-simulation differences (tracking, trigger, PID) are applied using calibration samples (e.g.,$D^{*+} \to D^0 \pi^+$).

Branching Fraction Calculation

Branching fractions are measured relative to the normalization mode $B^0 \to K^0_S \pi^+ \pi^-$. The ratio is calculated as:
$$\frac{\mathcal{B}(B^0_{(s)} \to K^0_S h^+ h'^-)}{\mathcal{B}(B^0 \to K^0_S \pi^+ \pi^-)} = \frac{\varepsilon_{B^0 \to K^0_S \pi^+ \pi^-}}{\varepsilon_{B^0_{(s)} \to K^0_S h^+ h'^-}} \frac{N_{B^0_{(s)} \to K^0_S h^+ h'^-}}{N_{B^0 \to K^0_S \pi^+ \pi^-}} \frac{f_d}{f_{d(s)}}$$
where $N$ is the fitted yield, $\varepsilon$ is the efficiency, and $f_d/f_s$ are the hadronization fractions of $b$ quarks into $B^0$ and $B^0_s$ mesons.

3. Key Contributions

1. First Observation of $B^0_s \to K^0_S K^+ K^-$: The analysis achieves a statistical significance of 8.6 standard deviations (after systematic uncertainties), marking the first definitive observation of this decay mode.
2. Improved Precision: The branching fraction ratios are measured with significantly reduced statistical uncertainties compared to previous LHCb results (using 3 fb$^{-1}$) and are consistent with world averages.
3. Comprehensive Systematic Treatment: The paper provides a detailed breakdown of systematic uncertainties, including fit model variations, efficiency binning schemes, and kinematic reweighting, demonstrating robustness against potential biases.
4. Dalitz Plot Efficiency Correction: The rigorous application of sPlot weights and square Dalitz plot binning to correct for phase-space dependent efficiency variations sets a high standard for multibody decay analyses.

4. Results

The paper reports the following ratios of branching fractions (statistical $\pm$ systematic $\pm f_s/f_d$):

Decay Mode	Ratio Value
$\mathcal{B}(B^0 \to K^0_S K^+ K^-) / \mathcal{B}(B^0 \to K^0_S \pi^+ \pi^-)$	$0.578 \pm 0.007 \pm 0.017$
$\mathcal{B}(B^0 \to K^0_S K^\pm \pi^\mp) / \mathcal{B}(B^0 \to K^0_S \pi^+ \pi^-)$	$0.1363 \pm 0.0035 \pm 0.0051$
$\mathcal{B}(B^0_s \to K^0_S \pi^+ \pi^-) / \mathcal{B}(B^0 \to K^0_S \pi^+ \pi^-)$	$0.269 \pm 0.011 \pm 0.015 \pm 0.008$
$\mathcal{B}(B^0_s \to K^0_S K^+ K^-) / \mathcal{B}(B^0 \to K^0_S \pi^+ \pi^-)$	$0.0303 \pm 0.0041 \pm 0.0025 \pm 0.0009$
$\mathcal{B}(B^0_s \to K^0_S K^\pm \pi^\mp) / \mathcal{B}(B^0 \to K^0_S \pi^+ \pi^-)$	$1.818 \pm 0.021 \pm 0.031 \pm 0.056$

Using the external value for $\mathcal{B}(B^0 \to K^0 \pi^+ \pi^-)$, the absolute time-integrated branching fractions are derived, e.g., $\mathcal{B}(B^0_s \to K^0_S K^+ K^-) = [1.50 \pm 0.20 (\text{stat}) \pm 0.12 (\text{syst}) \pm 0.08 (\text{ext})] \times 10^{-6}$.

5. Significance

• New Physics Sensitivity: The precise measurement of these branching fractions provides critical inputs for testing theoretical predictions. The $B^0_s \to K^0_S K^+ K^-$ mode, being a loop-dominated process, is particularly sensitive to new particles or interactions that could alter the decay amplitude.
• CP Violation and $\gamma$ Angle: While this paper focuses on branching fractions, the established signal yields and clean selection criteria pave the way for future time-dependent amplitude analyses. These will allow for the measurement of mixing-induced CP-violating phases and the determination of the CKM angle $\gamma$.
• Flavor Symmetry Tests: The results enable rigorous tests of isospin, U-spin, and SU(3) flavor symmetries by comparing the rates of $B^0$ and $B^0_s$ decays into similar final states.
• Validation of SM: The measured values are in good agreement with Standard Model predictions and previous experimental determinations, reinforcing the current understanding of $b$-quark decay dynamics while tightening the constraints on potential deviations.

In conclusion, this paper represents a major advancement in the study of charmless $B$ meson decays, successfully observing a previously elusive decay mode and providing the most precise measurements of these branching fraction ratios to date.