⚛️ phenomenology

Addressing the $S$ -wave scalar $f_0(1500)$ -resonance in quasi-four-body FCNC rare $B_s \to f_0(1500) (\to \pi^+ \pi^- ) \ell^+ \ell^- / \nu \bar \nu$ decays

This paper investigates the $f_0(150)$ resonance within the quasi-four-body FCNC rare decays $B_s \to f_0(1500)(\to \pi^+ \pi^-) \ell^+ \ell^- / \nu \bar \nu$ by constructing its twist-2 light-cone distribution amplitudes, calculating the associated transition form factors via QCD light-cone sum rules, and predicting the resulting branching fractions and differential decay widths to provide theoretical guidance for future experimental studies.

Original authors: Xue Zheng, Hai-Bing Fu, Dan-Dan Hu, Jing-Kai Yang, Jian-Qi Chen, Wan-Bing Luo

Published 2026-09-17

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Xue Zheng, Hai-Bing Fu, Dan-Dan Hu, Jing-Kai Yang, Jian-Qi Chen, Wan-Bing Luo

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the subatomic world, matter is built from a small set of fundamental particles called quarks, which are held together by a force known as the strong interaction. This force is so powerful that quarks never appear alone; they are always bound in groups, forming particles like protons and neutrons, or more exotic combinations called mesons. Among these mesons are a family of particles known as scalar resonances. These are short-lived states that exist for only a fleeting moment before decaying into lighter particles. For decades, physicists have struggled to understand the precise internal makeup of these scalar resonances. Are they simple pairs of a quark and an antiquark, or are they more complex structures, perhaps made of four quarks or even temporary molecular bonds? One particular particle, the f0(1500), has been a source of intense debate. It sits in a crowded region of the particle spectrum where several other resonances overlap, making it difficult to isolate its true nature. Its mass is around 1522 million electron volts, and it decays relatively quickly, but its behavior does not fit neatly into the standard models physicists use to describe the universe.

To solve this puzzle, a team of researchers from Guizhou Minzu University in China has turned their attention to a rare and complex type of decay involving the B-meson, a heavy particle containing a bottom quark. They focused on a specific process where a B-meson transforms into an f0(1500) resonance, which then immediately breaks apart into two pions, while simultaneously producing a pair of leptons, such as electrons or muons. This is a "quasi-four-body" decay because the final state consists of four distinct particles, even though the process happens in two steps: the B-meson first creates the intermediate f0(1500), and then that resonance splits. This specific decay channel is governed by a rare phenomenon called a flavor-changing neutral current, a process that is heavily suppressed in the standard model of particle physics. Because it happens so infrequently, any deviation from the predicted rate could signal the presence of new, undiscovered physics. However, to make a precise prediction, the researchers first had to understand the f0(1500) itself, specifically how its internal components share momentum.

The researchers began by constructing a detailed mathematical model of the f0(1500) resonance, treating it as a standard pair of a quark and an antiquark. They developed a new way to describe the distribution of momentum inside this particle, a concept known as a light-cone distribution amplitude. Think of this as a map showing how the speed and direction of the internal quarks are shared between them. Using a theoretical framework called the light-cone harmonic oscillator model, they calculated specific numerical values, known as moments, that characterize this distribution. These calculations were performed at specific energy scales, providing a set of parameters that describe the particle's internal structure with high precision. With this new map in hand, they were able to calculate the transition form factors, which are essentially numbers that describe how easily the heavy B-meson can transform into the lighter f0(1500) resonance. They found that at a specific point of maximum momentum transfer, the vector form factor was approximately 0.390, the scalar form factor was roughly -0.460, and the tensor form factor was about 0.568. These numbers are crucial because they determine the probability of the decay occurring.

Once these internal properties were established, the team simulated the entire decay process to predict how often it would happen. They calculated the branching fraction, which is the percentage of B-mesons that would undergo this specific transformation. For the decay into an electron-positron pair, the predicted rate was about 5.05 times 10 to the power of negative 7. For a muon-antimuon pair, the rate was nearly identical at 5.02 times 10 to the power of negative 7, a similarity that arises from the fundamental symmetry of nature known as lepton universality. The decay into a tau-lepton pair was much rarer, occurring only about 0.096 times 10 to the power of negative 7, because the heavy mass of the tau particle restricts the available space for the decay to happen. They also predicted the rate for the decay into a pair of neutrinos, which was significantly higher at 28.25 times 10 to the power of negative 7.

A critical part of their work involved comparing this complex four-particle outcome with a simpler, older way of thinking about the problem. Traditionally, physicists often use an approximation called the "narrow-width approximation," which assumes the intermediate f0(1500) particle is stable and has a fixed mass, ignoring the fact that it actually has a finite lifetime and a spread of masses. The researchers found that this old method overestimates the likelihood of the decay. When they accounted for the true, finite width of the f0(1500) and the way it interacts with the final particles, the predicted rate for the four-particle decay dropped significantly compared to the simplified three-particle calculation. This difference highlights the importance of treating the resonance as a real, dynamic object rather than a static point. The team also visualized the distribution of the decay energy, showing a clear peak where the f0(1500) resonance forms, followed by a rapid drop-off as the energy moves away from this specific value.

The study concludes that while the simplified models provide a rough estimate, a full understanding of these rare decays requires accounting for the complex, finite nature of the intermediate resonance. The researchers' predictions for the decay rates are now available for experimentalists to test. If future measurements at particle colliders match these numbers, it will confirm our current understanding of the f0(1500) as a quark-antiquark state and validate the standard model's description of these rare processes. If the measurements differ, it could point to new physics or a different internal structure for the f0(1500). For now, the work stands as a precise theoretical reference, offering a clear target for the next generation of experiments to aim for in their search to understand the hidden architecture of matter.

Technical Summary: Addressing the S-wave scalar $f_0(1500)$ -resonance in quasi-four-body FCNC rare $B_s \to f_0(1500)(\to \pi^+\pi^-)\ell^+\ell^-/\nu\bar{\nu}$ decays

Problem Statement
The internal structure of light scalar resonances, particularly the $f_0(1500)$ , remains a central open problem in non-perturbative Quantum Chromodynamics (QCD). While the $f_0(1500)$ is a representative scalar state with $J^{PC}=0^{++}$ , its composition (conventional $q\bar{q}$ , hybrid, molecule, or tetraquark) is heavily debated. A specific challenge arises in phenomenological studies of $B_s$ -meson quasi-four-body rare decays involving this resonance. Existing analyses often rely on the narrow-width approximation (NWA) or constant-width Breit-Wigner propagators, which fail to accurately describe broad resonances like $f_0(1500)$ (width $\sim 108$ MeV) and $f_0(1370)$ (width $200-500$ MeV). These approximations neglect finite-width effects and coupled-channel interference (specifically between $\pi\pi$ and $K\bar{K}$ channels), leading to potential large theoretical errors and a loss of phase information. Furthermore, systematic theoretical investigations of the quasi-four-body rare decay $B_s \to f_0(1500)(\to \pi^+\pi^-)\ell^+\ell^-$ are currently scarce.

Methodology
The authors adopt the conventional quark-antiquark ( $q\bar{q}$ ) picture for the $f_0(1500)$ and employ a multi-step theoretical framework:

Effective Hamiltonian: The study utilizes the Standard Model (SM) effective Hamiltonian for flavor-changing neutral current (FCNC) transitions $b \to s\ell^+\ell^-$ and $b \to s\nu\bar{\nu}$ . The decay amplitude is factorized into a weak transition amplitude, a resonance propagator, and a strong decay vertex.
Resonance Treatment: Instead of the NWA, the authors employ a relativistic Flatté line-shape function to describe the $f_0(1500)$ propagator. This approach incorporates energy-dependent partial decay widths for the $\pi\pi$ and $K\bar{K}$ channels, accounting for finite-width effects and coupled-channel dynamics.
Light-Cone Distribution Amplitudes (LCDA): To handle the non-perturbative dynamics of the $f_0(1500)$ , the authors construct a twist-2 LCDA scheme based on the Light-Cone Harmonic Oscillator (LCHO) model. This model incorporates the Brodsky-Huang-Lepage (BHL) prescription and Wigner-Melosh rotations. The longitudinal distribution function $\phi_{2;f_0(1500)}(x)$ is modeled using a Gegenbauer expansion, and the model parameters are constrained by the first Gegenbauer moment and the root-mean-square transverse momentum.
Form Factor Calculation: The $B_s \to f_0(1500)$ transition form factors (TFFs)—specifically $f_+$ , $f_-$ , and $f_T$ —are calculated using QCD Light-Cone Sum Rules (LCSR). The calculation includes contributions from both twist-2 and twist-3 LCDAs.
Extrapolation and Integration: The TFFs, initially reliable in the low-to-intermediate $q^2$ region ( $0 \le q^2 \le 8$ GeV $^2$ ), are extrapolated to the full physical kinematic region using a simplified $z(q^2)$ -series expansion. The double-differential decay width $d^2\Gamma/dsdq^2$ is then integrated over the physical phase space to obtain branching fractions.

Key Contributions

Novel LCDA Construction: The paper presents a specific twist-2 LCDA scheme for the $f_0(1500)$ derived from the LCHO model, providing explicit values for moments $\langle \xi^n \rangle$ and Gegenbauer moments $a_n$ at scales $\mu_0 = 1$ GeV and $\mu_k = 3$ GeV.
Finite-Width Analysis: The work moves beyond the narrow-width approximation by implementing the Flatté formula, thereby systematically addressing the finite-width effects and coupled-channel interference inherent to the $f_0(1500)$ resonance in quasi-four-body decays.
First Systematic Study: This is the first systematic theoretical investigation of the quasi-four-body rare decays $B_s \to f_0(1500)(\to \pi^+\pi^-)\ell^+\ell^-$ and $B_s \to f_0(1500)(\to \pi^+\pi^-)\nu\bar{\nu}$ .
Comprehensive TFF Calculation: The authors provide the first LCSR calculation of the $B_s \to f_0(1500)$ TFFs that includes both twist-2 and twist-3 contributions, offering a more complete description of the non-perturbative inputs.

Results

Transition Form Factors (at $q^2=0$ ):
- $f_+^{B_s f_0(1500)}(0) = 0.390^{+0.047}_{-0.046}$
- $f_-^{B_s f_0(1500)}(0) = -0.460^{+0.051}_{-0.055}$
- $f_T^{B_s f_0(1500)}(0) = 0.568^{+0.069}_{-0.065}$
  These results show good agreement with previous LCSR predictions (e.g., LCSR'08, LCSR'10, LCSR'23) within uncertainties, though they differ from some pQCD results.
Branching Fractions:
- The branching fractions for the quasi-four-body decays are found to be significantly smaller than those predicted under the narrow-width approximation. For example, $\mathcal{B}(B_s \to f_0(1500)(\to \pi^+\pi^-)\mu^+\mu^-) \approx 5.02 \times 10^{-7}$ , compared to $\approx 6.25 \times 10^{-7}$ for the three-body NWA case.
- The electron and muon channels are nearly equal due to lepton universality.
- The $\tau^+\tau^-$ channel is heavily suppressed by phase space ( $\sim 10^{-9}$ ).
- The neutrino channel $B_s \to f_0(1500)(\to \pi^+\pi^-)\nu\bar{\nu}$ has the largest branching fraction among the modes studied, reported as $28.25^{+6.75}_{-5.84} \times 10^{-7}$ (approximately $2.8 \times 10^{-6}$ ).
Distributions: The double-differential decay width $d^2\Gamma/dsdq^2$ exhibits a clear peak in the $f_0(1500)$ resonance region, with the distribution shape modulated by the finite width of the resonance.

Significance and Claims
The paper claims that its predictions provide a necessary theoretical reference for future experimental measurements of these rare decays, which are sensitive probes for testing the Standard Model and searching for new physics. The authors emphasize that the reduction in branching fractions observed in the quasi-four-body analysis compared to the NWA is a direct consequence of the finite width of the $f_0(1500)$ , final-state strong interactions, and phase-space reduction.

The authors maintain a modest stance, noting that:

There are currently no experimental measurements or prior systematic theoretical calculations for this specific quasi-four-body decay channel to compare against.
The internal structure of the $f_0(1500)$ remains controversial; the results are based on the $q\bar{q}$ assumption.
The primary value of the work lies in establishing a consistent theoretical framework (incorporating Flatté parametrization and LCHO-based LCDAs) to guide future phenomenological research and experimental searches.

More like this