Amplitude Analysis and Branching Fraction Measurement of $D^+ \to \pi^+\pi^0\pi^0$

Using 20.3 fb$^{-1}$ of $e^+e^-$ collision data collected by the BESIII detector, this paper presents the first amplitude analysis of the $D^+ \to \pi^+\pi^0\pi^0$ decay, identifying the $D^+ \to \rho(770)^+\pi^0$ component as dominant and reporting precise measurements of the total branching fraction, intermediate fit fractions, and CP asymmetries.

The Gist

Imagine the subatomic world as a high-stakes dance floor where particles collide, spin, and sometimes break apart into smaller pieces. This paper is a detailed report from the BESIII experiment, a giant "camera" (detector) located in China, which watched millions of these tiny dance moves to understand a specific breakup: a particle called a $D^+$ meson splitting into three pions (a type of particle similar to a proton's lighter cousin).

Here is the story of what they found, explained without the heavy math.

1. The Setup: A Massive Photo Album

The scientists didn't just take a snapshot; they compiled a massive photo album. They collected data from 20.3 billion electron-positron collisions (imagine smashing two tiny magnets together at nearly the speed of light). This huge amount of data allowed them to see rare events that would be invisible in a smaller sample.

Their goal was to study the decay $D^+ \to \pi^+ \pi^0 \pi^0$.

• The $D^+$ meson: The dancer starting the routine.
• The $\pi^+$ and two $\pi^0$s: The three pieces the dancer breaks into.

2. The Mystery: How Did the Breakup Happen?

When a particle breaks into three pieces, it rarely happens all at once. Usually, it's a two-step process. Think of it like a parent ($D^+$) breaking up a toy into three parts.

• Scenario A: The parent breaks the toy into a big chunk and a small chunk, then the big chunk breaks again.
• Scenario B: The parent breaks it into two medium chunks, then one of those breaks again.

In physics, these "chunks" are called intermediate resonances. The scientists wanted to know: Which path did the $D^+$ meson take?

3. The Main Discovery: The "Star" of the Show

Using a technique called Amplitude Analysis (which is like using a super-computer to reverse-engineer the dance steps from the final positions of the dancers), they found that one specific path was the clear winner.

• The Winner: The $D^+$ meson almost always turned into a $\rho(770)^+$ particle and a $\pi^0$ first. Then, the $\rho(770)^+$ quickly fell apart into the remaining $\pi^+$ and $\pi^0$.
• The Analogy: Imagine a magician pulling a rabbit out of a hat, but the rabbit is actually a hat with a smaller rabbit inside. The "big hat" ($\rho$) is the most common way the trick happens.
• The Result: This specific path accounts for about 63.5% of all the breakups. The scientists measured how often this happens (the "Branching Fraction") and found it to be roughly 3 out of every 1,000 $D^+$ mesons.

4. The Supporting Cast

While the $\rho(770)^+$ was the star, there were other, less common ways the breakup could happen:

• A heavier version of the $\rho$ particle ($\rho(1450)$).
• A different particle called $f_2(1270)$.
• A "S-wave" state (a fuzzy, non-resonant cloud of particles).
• The "Interference" Effect: Sometimes, these different paths happen at the same time and mess with each other, like two sound waves canceling each other out or making a louder noise. The scientists measured these "interference fractions" to understand how the different paths mix.

5. The "Mirror" Test: Looking for Differences (CP Violation)

One of the biggest questions in physics is: Does the universe treat matter and antimatter exactly the same?

• The $D^+$ is matter. Its twin, the $D^-$, is antimatter.
• If the laws of physics are perfectly symmetrical, the $D^+$ and $D^-$ should break apart in exactly the same way, at the same rate.
• If they break apart differently, it's called CP Violation (a hint that the universe has a slight preference for matter over antimatter).

The Result: The scientists compared the "dance moves" of the $D^+$ and the $D^-$. They found no significant difference. The rates were identical within the margin of error.

• The Analogy: It's like watching a left-handed dancer and a right-handed dancer perform the exact same routine. They move slightly differently in their hands, but the overall speed and style are the same. No "new physics" (like a hidden force) was found here.

6. Why Does This Matter?

• Testing the Rules: Theoretical physicists have built models (like the "Pole Model" or "Factorization") to predict how often these breakups happen. The BESIII results are like a final exam for these models.
• The Score: The dominant path ($\rho(770)^+$) matches some predictions but disagrees slightly with others. This helps scientists refine their theories about the "strong force" (the glue holding particles together), which is notoriously difficult to calculate.
• Precision: By measuring the exact frequency of these events (about 4.84 out of every 1,000 total decays), they provide a solid benchmark for future experiments.

Summary

The BESIII collaboration took a massive dataset of particle collisions and performed a detailed "forensic analysis" of how a $D^+$ meson breaks into three pions. They discovered that the breakup is dominated by a specific intermediate step involving a $\rho(770)^+$ particle. They also confirmed that matter and antimatter behave identically in this process, finding no evidence of the mysterious "CP violation" that might explain why our universe is made of matter. This work provides precise numbers that help physicists tune their theories of the subatomic world.

Technical Summary

Technical Summary: Amplitude Analysis and Branching Fraction Measurement of $D^+ \to \pi^+\pi^0\pi^0$

Problem and Motivation
Non-leptonic decays of charmed mesons offer a critical platform for investigating non-perturbative Quantum Chromodynamics (QCD) effects and testing weak interaction mechanisms. The three-body decay $D^+ \to \pi^+\pi^0\pi^0$ is theoretically expected to be dominated by the quasi-two-body channel $D^+ \to \rho(770)^+\pi^0$. While theoretical frameworks such as the pole model, topological diagrammatic approaches (TDA), and factorization-assisted topological-amplitude (FAT) methods provide predictions for the branching fraction (BF) of this dominant mode, precise experimental constraints are necessary to refine these models. Furthermore, while direct CP violation has been observed in neutral $D^0$ decays, no clear CP asymmetry has been established in $D^+$ decays. Searching for CP violation in the singly Cabibbo-suppressed decay $D^+ \to \pi^+\pi^0\pi^0$ is significant for probing potential new physics, particularly as CP asymmetries in multibody decays can arise from long-distance interference effects among intermediate resonances, potentially revealing local asymmetries even if the global asymmetry is small.

Methodology
The analysis utilizes $e^+e^-$ collision data collected by the BESIII detector at a center-of-mass energy of $\sqrt{s} = 3.773$ GeV, corresponding to an integrated luminosity of 20.3 fb$^{-1}$. The study employs a "tag" technique where $D^+D^-$ pairs are produced via $\psi(3770)$ decay. One $D^-$ meson is reconstructed in a specific hadronic decay mode (single tag, ST), while the signal $D^+$ is reconstructed in the $D^+ \to \pi^+\pi^0\pi^0$ mode (double tag, DT).

• Event Selection: Charged tracks and photons are selected using the BESIII detector components (MDC, TOF, EMC). Candidates are vetted using beam-constrained mass ($M_{BC}$) and energy difference ($\Delta E$) variables. Specific vetoes are applied to suppress backgrounds, including a $K_S^0$ veto to reject $D^+ \to K_S^0\pi^+$ and kinematic constraints to reject mispartition backgrounds from $D^0\bar{D}^0$ pairs.
• Signal Extraction: An unbinned two-dimensional maximum likelihood fit is performed on the $M_{BC}^{tag}$ vs. $M_{BC}^{sig}$ distribution to extract the signal yield, separating it from combinatorial backgrounds and misreconstructed events.
• Amplitude Analysis: The signal sample undergoes an unbinned amplitude analysis using the isobar model. The total amplitude is a coherent sum of intermediate processes, including $D^+ \to \rho(770)^+\pi^0$, $D^+ \to \rho(1450)^+\pi^0$, $D^+ \to f_2(1270)\pi^+$, and the $\pi^0\pi^0$ S-wave (described via K-matrix parametrization). The analysis accounts for Bose symmetry between the two identical $\pi^0$ mesons. Detection efficiencies and background shapes are modeled using Monte Carlo (MC) simulations, corrected for data-MC differences in tracking, particle identification (PID), and $\pi^0$ reconstruction.
• CP Asymmetry: A simultaneous fit is performed for $D^+$ and $D^-$ samples, allowing complex coefficients for isobar amplitudes and production parameters for the K-matrix S-wave to vary independently between charges.
• Branching Fraction: The BF is calculated using the ratio of double-tag to single-tag yields, corrected for detection efficiencies and sub-decay branching fractions (specifically $\pi^0 \to \gamma\gamma$).

Key Contributions and Results
This work presents the first amplitude analysis of the $D^+ \to \pi^+\pi^0\pi^0$ decay.

1. Branching Fraction Measurement:
   The total branching fraction is measured to be:
   $$\mathcal{B}(D^+ \to \pi^+\pi^0\pi^0) = (4.84 \pm 0.05_{\text{stat}} \pm 0.05_{\text{syst}}) \times 10^{-3}$$
   The CP asymmetry in the total branching fraction is measured as:
   $$A_{CP} = (-1.4 \pm 1.0_{\text{stat}} \pm 0.6_{\text{syst}})\%$$
   No evidence for CP violation is observed in the integrated rate.

2. Amplitude Analysis and Fit Fractions:
   The analysis confirms that the decay is dominated by the $D^+ \to \rho(770)^+\pi^0$ component. The measured fit fractions (FF) and intermediate branching fractions are:

   • $D^+ \to \rho(770)^+\pi^0$: FF = $63.5 \pm 2.0 \pm 1.2\%$; $\mathcal{B} = (3.08 \pm 0.10_{\text{stat}} \pm 0.07_{\text{syst}}) \times 10^{-3}$.
   • $D^+ \to \rho(1450)^+\pi^0$: FF = $5.2 \pm 0.8 \pm 0.7\%$.
   • $D^+ \to f_2(1270)\pi^+$: FF = $4.5 \pm 0.3 \pm 0.2\%$.
   • $D^+ \to (\pi^0\pi^0)_{S\text{-wave}}\pi^+$: FF = $11.6 \pm 0.9 \pm 0.5\%$.
     The sum of fit fractions is less than 100% due to destructive interference, while the net interference is constructive, with significant contributions from $f_2(1270)\pi^+ \times \rho(770)^+\pi^0$ and $\rho(1450)^+\pi^+ \times (\pi^0\pi^0)_{S\text{-wave}}\pi^+$.

3. CP Asymmetry in Intermediate States:
   CP asymmetries were measured for specific intermediate amplitudes. The results are consistent with zero within uncertainties:

   • $A_{CP}(\rho(770)^+\pi^0) = +3.7 \pm 2.9 \pm 1.5\%$
   • $A_{CP}(\rho(1450)^+\pi^0) = -17.5 \pm 9.8 \pm 7.0\%$
   • $A_{CP}(f_2(1270)\pi^+) = -2.2 \pm 6.6 \pm 3.6\%$
   • $A_{CP}((\pi^0\pi^0)_{S\text{-wave}}\pi^+) = +10.2 \pm 6.6 \pm 3.8\%$

Significance
The paper claims that these measurements provide stringent constraints on theoretical models of charmed meson decays. Specifically, the measured branching fraction for the dominant $D^+ \to \rho(770)^+\pi^0$ channel is compared with predictions from the pole model, TDA, and FAT methods. The authors note that while the pole model prediction is consistent with their result, the FAT prediction is approximately 2.3 standard deviations away from the measurement. Additionally, the precise determination of fit fractions and phases for the intermediate states, including the complex S-wave structure, offers new data for understanding the topological amplitude composition in non-leptonic charm decays. The absence of observed CP violation in this channel, despite the sensitivity to local effects in the Dalitz plot, sets limits on potential new physics contributions in singly Cabibbo-suppressed $D^+$ decays.