Observation of $D_s^+ \to a_0(980)^+f_0(500)$ in the Amplitude Analysis of $D_s^+ \to π^+ π^0 π^0 η$

Using 7.33 fb$^{-1}$ of data collected by the BESIII detector, this study reports the first observation of the decay $D_s^+ \to \pi^+\pi^0\pi^0\eta$ and measures an unexpectedly large branching fraction for the intermediate process $D_s^+ \to a_0(980)^+f_0(500)$ with over $10\sigma$ significance, providing new constraints on the internal structure of light scalar mesons.

The Gist

Imagine the subatomic world as a bustling, chaotic construction site. In this site, particles are constantly being built, broken apart, and rebuilt. For a long time, physicists have been trying to understand the blueprints of some very specific, mysterious "bricks" called scalar mesons (specifically the $a_0(980)$ and $f_500$).

The big question has been: Are these bricks made of just two smaller pieces (a quark and an antiquark), or are they actually complex, four-piece structures (tetraquarks)?

This paper is like a team of master detectives (the BESIII Collaboration) finally solving a cold case that reveals the true nature of these mysterious bricks.

The Crime Scene: The $D_s$ Particle

The detectives focused on a specific event: the decay of a particle called the $D_s^+$ (a charmed meson). Think of the $D_s^+$ as a heavy, unstable delivery truck that crashes and breaks apart into smaller pieces.

Usually, when this truck crashes, it spits out a chaotic mess of particles: a positive pion ($\pi^+$), two neutral pions ($\pi^0$), and an eta particle ($\eta$).

• The Goal: The team wanted to see how this crash happened. Did the pieces fly off randomly, or did they form temporary "clumps" (intermediate particles) before breaking apart completely?

The Investigation: Amplitude Analysis

To solve this, the team didn't just count the debris; they performed a high-tech "amplitude analysis."

• The Analogy: Imagine listening to a symphony orchestra. You can hear the final song, but to understand the music, you need to separate the sound of the violins, the trumpets, and the drums.
• The Method: The team used a massive dataset (equivalent to 7.33 "inverse femtobarns" of data—think of this as a giant pile of crash reports) to separate the different "instruments" (decay paths) playing in the background.

The Big Discovery: Two New "Clumps"

The analysis revealed two major ways the $D_s^+$ truck breaks apart:

1. The Dominant Path (The "Standard" Crash):
   Most of the time (about 60%), the truck breaks into a specific clump called $a_1(1260)$ and an $\eta$ particle. This was expected, like finding a familiar car model in a junkyard. It confirmed that the laws of physics (specifically "isospin symmetry") are working correctly, just like a perfectly balanced scale.

2. The Shocking Discovery (The "Exotic" Crash):
   The real surprise was finding a second major path where the truck breaks into $a_0(980)$ and $f_0(500)$.

   • The Surprise: This is a "Double Scalar" decay ($D_s \to S + S$). In the standard model of physics, if these particles were simple two-piece bricks, this crash should be extremely rare—almost impossible.
   • The Result: The team found this happens 1% of the time. That is huge for this type of reaction. It's like finding a rare, exotic flower growing in a desert where nothing should grow.

What Does This Mean? (The Tetraquark Theory)

Why is this 1% so important?

• The Two-Piece Theory: If the $a_0(980)$ and $f_0(500)$ were simple two-piece bricks, the "engine" driving this crash would be very weak, and we shouldn't see it often.
• The Four-Piece Theory: If these particles are actually tetraquarks (four-piece structures), the "engine" is much stronger. The fact that they saw this crash so frequently suggests these particles are likely tetraquarks (or at least have a strong tetraquark component).

The Metaphor:
Imagine you are trying to open a locked box.

• If the lock is a simple keyhole (two-quark), you need a very specific, tiny key.
• If the lock is a complex combination dial (tetraquark), you can open it with a much larger, more common tool.
  The fact that the physicists found the "large tool" working so often suggests the lock is actually a complex combination dial.

The "Ghost" in the Machine

The team also noticed a strange bump in the data around 0.8 GeV (a specific energy level) that didn't match any known particle.

• The Analogy: It's like hearing a ghostly whisper in the recording that doesn't match any of the known instruments.
• Significance: This suggests there might be even more mysterious structures or interactions happening that we haven't discovered yet. It's a hint that our map of the subatomic world is still incomplete.

The Bottom Line

This paper is a landmark discovery because:

1. First Observation: It is the first time scientists have ever seen the $D_s^+$ decay into two scalar mesons ($a_0$ and $f_0$).
2. Structural Clue: The high frequency of this decay provides strong evidence that light scalar mesons are likely tetraquarks (four-quark states), solving a decades-old mystery about their internal structure.
3. New Physics: It opens the door to understanding how the strong force (the glue holding matter together) behaves in complex, non-perturbative ways.

In short, the BESIII team didn't just find a new particle; they found a new way of building matter, proving that the subatomic world is more complex and fascinating than we previously imagined.

Technical Summary

1. Problem and Motivation

The paper addresses two primary challenges in hadronic physics:

• Nature of Light Scalar Mesons: The internal structures of light scalar mesons below 1 GeV/$c^2$ (specifically $a_0(980)$ and $f_0(500)$) remain poorly understood. Theoretical models debate whether they are conventional $q\bar{q}$ states, compact tetraquarks, or hadronic molecules.
• Decay Dynamics of $D_s$ Mesons: While many $D_{(s)} \to SP$ (Scalar-Pseudoscalar) and $SV$ (Scalar-Vector) decays have been studied, no $D_s \to SS$ (Scalar-Scalar) decays had been observed.
  • If scalar mesons are conventional $q\bar{q}$ states, the decay $D^+_s \to a_0(980)^+f_0(500)$ is expected to proceed via pure weak annihilation (WA), which is highly suppressed.
  • If scalar mesons are tetraquarks, the decay could receive contributions from external $W$-emission diagrams, potentially enhancing the branching fraction (BF) significantly.
• Isospin Symmetry Validation: Previous studies observed $D^+_s \to a_1(1260)^+\eta$ with $a_1(1260)^+ \to \rho(770)^0\pi^+$. This paper aims to observe the charge-conjugate channel $a_1(1260)^+ \to \rho(770)^+\pi^0$ in the $D^+_s \to \pi^+\pi^0\pi^0\eta$ decay to validate isospin symmetry.

2. Methodology

The analysis utilizes data collected by the BESIII detector at the BEPCII collider.

• Dataset: An integrated luminosity of 7.33 fb$^{-1}$ collected in $e^+e^-$ collisions at center-of-mass energies between 4.128 and 4.226 GeV.
• Event Selection (Double-Tag Method):
  • Single Tag (ST): $D^-_s$ mesons are reconstructed in nine hadronic decay modes (e.g., $K^-_S K^-$, $K^+K^-\pi^-$, etc.).
  • Double Tag (DT): The signal $D^+_s \to \pi^+\pi^0\pi^0\eta$ is reconstructed in the recoil against the tagged $D^-_s$.
  • Kinematic Fits: A 9-constraint kinematic fit is applied to improve momentum resolution and reduce combinatorial background.
  • Background Suppression: Events with invariant masses corresponding to $\eta'$ ($M_{\pi^0\pi^0\eta} \in [0.946, 0.970]$ GeV/$c^2$) are separated to isolate the non-$\eta'$ signal. Fake $\pi^0$ combinations are vetoed.
• Amplitude Analysis:
  • An unbinned maximum likelihood fit is performed using the isobar model in the covariant tensor formalism.
  • Signal Amplitudes: Modeled as a coherent sum of intermediate resonances. Key components include:
    • $a_1(1260)^+$ (modeled with Relativistic Breit-Wigner).
    • $\rho(770)$ (modeled with Gounaris-Sakurai line shape).
    • $a_0(980)$ (modeled with a coupled Flatté formula).
    • $f_0(500)$ (modeled using a formula derived from $\pi\pi$ elastic scattering).
    • $S$-wave $\pi^0\pi^0$ (modeled using K-matrix formalism).
  • Background: Modeled using inclusive Monte Carlo (MC) samples weighted by an XGBoost algorithm to match data distributions.
• Branching Fraction Measurement:
  • Calculated using the formula: $\mathcal{B} = \frac{Y_{sig}}{\sum Y_{tag} \cdot \epsilon_{tag,sig}/\epsilon_{tag}}$, where $Y$ represents yields and $\epsilon$ represents efficiencies derived from MC simulations.

3. Key Contributions

1. First Observation of $D^+_s \to \pi^+\pi^0\pi^0\eta$: This is the first reported observation of this specific decay mode.
2. First Observation of $D^+_s \to SS$: The paper reports the first observation of a $D_s \to$ Scalar-Scalar decay, specifically $D^+_s \to a_0(980)^+f_0(500)$.
3. Amplitude Analysis of $D^+_s \to \pi^+\pi^0\pi^0\eta$: A comprehensive amplitude analysis is performed to disentangle the complex multi-body final state into its intermediate resonant components.
4. Isospin Symmetry Test: The study provides a direct comparison between $a_1(1260)^+ \to \rho(770)^0\pi^+$ and $a_1(1260)^+ \to \rho(770)^+\pi^0$ decays.

4. Results

• Branching Fractions (BF):

  • Total Non-$\eta'$ BF: $\mathcal{B}(D^+_s \to \pi^+\pi^0\pi^0\eta|_{\text{non-}\eta'}) = (2.97 \pm 0.23_{\text{stat}} \pm 0.14_{\text{syst}})\%$.
  • $\eta'$ Contribution: $\mathcal{B}(D^+_s \to \pi^+\eta', \eta' \to \pi^0\pi^0\eta) = (0.82 \pm 0.08_{\text{stat}} \pm 0.02_{\text{syst}})\%$.
  • Derived $\mathcal{B}(D^+_s \to \pi^+\eta')$: $(3.66 \pm 0.36_{\text{stat}} \pm 0.09_{\text{syst}} \pm 0.08_{\text{PDG}})\%$.
  • Impact: These results reduce the inferred branching fraction of unobserved $D^+_s \to \eta X$ decays to $(0.1 \pm 3.1)\%$.

• Intermediate States (from Amplitude Analysis):

  • $D^+_s \to a_1(1260)^+\eta$ ($a_1 \to \rho^+\pi^0$):
    • Fraction: $59.7 \pm 5.2 \pm 2.9\%$.
    • BF: $(1.77 \pm 0.21_{\text{stat}} \pm 0.12_{\text{syst}})\%$.
    • Significance: $>10\sigma$.
    • Significance: This validates isospin symmetry when compared to the previously measured $a_1 \to \rho^0\pi^+$ channel ($1.73 \pm 0.14\%$).
  • $D^+_s \to a_0(980)^+f_0(500)$:
    • Fraction: $33.1 \pm 4.9 \pm 6.8\%$.
    • BF: $(0.98 \pm 0.16_{\text{stat}} \pm 0.22_{\text{syst}})\%$.
    • Significance: $>10\sigma$.
  • $D^+_s \to \pi^+(\pi^0\pi^0)_{S\text{-wave}}\eta$:
    • BF: $(0.28 \pm 0.06_{\text{stat}} \pm 0.13_{\text{syst}})\%$ (Significance $4.5\sigma$).
  • $D^+_s \to a_0(980)^0\rho(770)^+$: Not observed with current statistics (Significance $2.4\sigma$).

5. Significance and Implications

• Evidence for Tetraquark Structure or Final State Interactions: The measured branching fraction for $D^+_s \to a_0(980)^+f_0(500)$ is unexpectedly large (order of $10^{-2}$).
  • Under the conventional $q\bar{q}$ meson hypothesis, this decay should be suppressed by the small $a_0(980)^+$ decay constant due to the pure weak annihilation mechanism.
  • The large observed BF suggests either a substantial tetraquark component in the $a_0(980)$ and $f_0(500)$ (allowing for external $W$-emission contributions) or significant enhancements from final state interactions (e.g., triangle rescattering).
• Constraints on Scalar Meson Models: These results provide crucial new constraints for theoretical models attempting to explain the nature of light scalar mesons, favoring interpretations beyond simple quark-antiquark states.
• Validation of Isospin Symmetry: The consistency between the $\rho^0\pi^+$ and $\rho^+\pi^0$ decay modes of the $a_1(1260)^+$ reinforces the validity of isospin symmetry in charmed meson decays.
• Completeness of $D_s$ Decays: The measurement significantly reduces the "missing" branching fraction for $D_s \to \eta X$ decays, improving the global understanding of charmed meson decay dynamics.

In conclusion, this paper presents a landmark observation of a $D_s \to SS$ decay with a branching fraction that challenges the conventional quark model, offering compelling evidence for exotic hadron structures or complex strong interaction dynamics.