Precise theoretical prediction on branching fractions and polarizations of $D \to V V$ decays

This paper presents a precise theoretical analysis of $D \to VV$ decays using the factorization-assisted topological-amplitude (FAT) approach, where 10 nonfactorizable parameters are globally fitted to 36 experimental data points to reveal that large strong phases in the $E$ amplitude cause destructive interference with the $C$ component, leading to a polarization hierarchy ($f_\parallel > f_L$) that contradicts naive factorization predictions.

The Gist

Imagine the subatomic world as a bustling, chaotic dance floor. In this paper, physicists are trying to predict exactly how certain dancers (particles called D mesons) will spin and pair up when they break apart into two new dancers (particles called vector mesons).

Here is the story of their discovery, explained without the heavy math.

1. The Mystery: The "Spin" Puzzle

In the world of particle physics, when a heavy particle decays (breaks apart), the two new particles don't just fly off randomly. They have a specific "spin" or orientation, much like how a figure skater can spin forward, sideways, or on their side.

For a long time, physicists had a simple rule of thumb (called "naive factorization") that predicted: "Most of the time, the dancers will spin forward (longitudinal)."

But experiments kept showing something weird. Sometimes, the dancers were spinning sideways (transverse) just as much as, or even more than, forward. It was like predicting a soccer ball would always roll straight, but watching it curve wildly to the left instead. This was the "Polarization Anomaly."

2. The Old Map vs. The New GPS

To understand these dances, physicists use a "Topological Map." They draw diagrams showing how the particles interact:

• T (Tree): The main, easy path.
• C (Color-Suppressed): A slightly harder path.
• E (Exchange): A path where particles swap partners.
• A (Annihilation): A path where they destroy each other.

The Problem: The old way of using this map assumed all the dancers were identical twins (a symmetry called SU(3)). But in reality, they are cousins with different sizes and weights. The old map ignored these differences, leading to wrong predictions.

The Solution (The FAT Approach): The authors used a new method called FAT (Factorization-Assisted Topological-Amplitude).

• The Analogy: Imagine you are trying to predict the outcome of a cooking competition. The old method assumed all chefs used the exact same ingredients and tools. The FAT method says, "Okay, let's acknowledge that Chef A has a fancy oven and Chef B has a cast-iron skillet. We will factor out the specific tools (the 'form factors') and focus only on the secret sauce (the non-factorizable parts) that makes the dish taste unique."
• By separating the "tools" from the "secret sauce," they found that the secret sauce is actually universal. They could describe the messy, complex interactions of all these different decays using just a handful of universal numbers.

3. The Discovery: Why the Dancers Spin Sideways

The team took 36 different experimental data points (measurements from real experiments) and used their new "GPS" to fit the 10 most important "secret sauce" numbers.

They found two main reasons why the dancers spin sideways more than expected:

1. The "Clash" of Amplitudes: In some dances, the "Exchange" move (E) and the "Tree" move (T) happen at the same time. Because of a specific "phase shift" (think of it as a timing delay or a step out of sync), these two moves cancel each other out when trying to spin forward, but they boost each other when spinning sideways.

   • Result: The forward spin gets crushed, and the sideways spin wins. This explains why some decays have more sideways spin than the old rules predicted.

2. The "D-Wave" Surprise: Physicists also look at the "shape" of the dance (S-wave, P-wave, D-wave). The old rule said the "S-wave" (a simple, round dance) should always be the biggest.

   • The Twist: The team found that for dances involving the "Exchange" (E) move, the "D-wave" (a more complex, figure-eight dance) actually becomes the dominant move.
   • Why? The "Exchange" move has a huge timing difference between its forward and sideways components. This causes the simple round dance to cancel out, leaving the complex figure-eight dance as the winner. This explains recent experiments that showed the "D-wave" was winning, contradicting old theories.

4. The Prediction: What's Next?

The authors didn't just explain the past; they predicted the future.

• They calculated the exact probabilities (branching fractions) and spin orientations for 28 different decay modes.
• For the 28 modes they could check, their predictions matched the existing experimental data perfectly.
• The Treasure Map: They identified several "unobserved" dances (decays that haven't been measured yet). They predict these will happen frequently enough to be seen by major detectors like BESIII, STCF, Belle II, and LHCb.

Summary

Think of this paper as a master choreographer who finally figured out the secret rhythm of a chaotic dance floor.

• The Problem: The dancers were spinning sideways when they were supposed to spin forward.
• The Fix: They realized the dancers weren't identical; they had different "tools" and "timing."
• The Result: By accounting for these differences, they created a new set of rules that perfectly explains the weird spins and predicts exactly which new dances we should look for next.

This work bridges the gap between messy experimental data and clean theoretical predictions, helping us understand the fundamental forces that govern how matter transforms.

Technical Summary

1. Problem Statement

The paper addresses the theoretical challenges in describing $D \to VV$ decays (where $D$ represents charmed mesons $D^0, D^+, D_s^+$ and $V$ represents vector mesons $\rho, K^*, \omega, \phi$).

• The "Polarization Anomaly": Experimental data reveals significant deviations from naive factorization predictions. Specifically, transverse polarization fractions ($f_\parallel, f_\perp$) are often larger than expected, and in some modes, the D-wave contribution exceeds the S-wave contribution, contradicting the naive hierarchy $|S|^2 > |P|^2 > |D|^2$.
• Theoretical Limitations: Unlike $B$-meson decays, the charm quark mass ($m_c$) is not large enough to justify a heavy-quark expansion ($1/m_c$), rendering standard QCD factorization, perturbative QCD, and SCET less reliable.
• Data Scarcity: While $D \to PP$ and $D \to PV$ decays have been systematically studied, $D \to VV$ data is sparse and lacks precision, making it difficult to constrain phenomenological models using traditional topological diagram approaches which rely heavily on exact $SU(3)$ flavor symmetry.

2. Methodology: Factorization-Assisted Topological-Amplitude (FAT) Approach

The authors employ the FAT approach, a framework designed to handle $SU(3)$ symmetry breaking effects while minimizing free parameters.

• Topological Decomposition: Decay amplitudes are decomposed into four topological diagrams:
  • $T$ (Color-favored emission): Treated as factorizable.
  • $C$ (Color-suppressed emission): Dominated by non-factorizable contributions.
  • $E$ (W-exchange): Dominated by non-factorizable contributions.
  • $A$ (W-annihilation): Neglected due to negligible contribution and lack of data to constrain it.
• Factorization of Symmetry Breaking:
  • The authors explicitly factor out form factors and decay constants from the topological amplitudes.
  • The remaining non-factorizable contributions are assumed to be universal across all decay modes.
  • This allows the incorporation of $SU(3)$ breaking effects (via distinct masses and decay constants for different mesons) without introducing a new free parameter for every specific decay mode.
• Parameterization:
  • $T$ amplitude: Calculated using naive factorization with a single scale parameter $\mu$ (constrained to be $< m_c/2$).
  • $C$ and $E$ amplitudes: Parameterized by universal magnitudes ($\chi$) and strong phases ($\phi$) for each polarization state ($0, \parallel, \perp$).
  • Reduction: By neglecting the $A$ diagram and the transverse component of the $E$ diagram ($\chi_E^\perp \approx 0$), the model reduces the free parameters to 11 (9 magnitudes/phases for $C$ and $E$, plus the scale $\mu$).

3. Key Contributions

• Global Fit: The authors performed a global fit to 36 experimental data points (branching fractions and partial wave components) from BESIII, LHCb, CLEO-c, FOCUS, and Mark III.
• Precision Extraction: They successfully extracted 10 non-factorizable parameters associated with the $C$ and $E$ diagrams with high precision (yielding $\chi^2/\text{d.o.f.} = 8.43$).
• Resolution of Anomalies: The framework successfully explains:
  • The inversion of polarization fractions ($f_\parallel > f_L$) in specific modes.
  • The dominance of D-waves over S-waves in certain decays.
• Prediction Generation: The paper provides precise theoretical predictions for 28 decay modes, including 28 branching fractions and polarization fractions, many of which are unmeasured or poorly constrained.

4. Key Results and Physical Insights

• Amplitude Hierarchies:
  • The fit reveals $|C^0| > |T^0| > |E^0|$ and $|T^{\parallel(\perp)}| \sim |C^{\parallel(\perp)}| > |E^{\parallel(\perp)}|$.
  • Crucially, the non-factorizable $C$ contribution is comparable in magnitude to the factorizable $T$ contribution, challenging the assumption that $C$ is negligible.
• Polarization Mechanism:
  • $T$ and $C$ dominated modes: Follow the hierarchy $|A^0| > |A^\parallel| > |A^\perp|$, leading to longitudinal dominance ($f_L > f_\parallel$).
  • $E$-diagram interference: A large strong phase in the longitudinal $E$ amplitude ($\phi_E^0 \approx 2.23$) causes destructive interference with the $C$ longitudinal component. This suppresses $f_L$ and enhances $f_\parallel$, explaining modes where $f_\parallel > f_L$ (e.g., $D^0 \to \omega \rho^0$).
• Partial Wave Hierarchy:
  • For $T$ and $C$ diagrams, $|S| > |P| > |D|$ holds.
  • For $E$-diagram dominated modes, the large phase difference $|\phi_E^0 - \phi_E^\parallel|$ reverses the hierarchy to $|S| < |D|$. This explains experimental observations where D-wave branching fractions exceed S-wave fractions (e.g., $D^0 \to \rho^0 \bar{K}^{*0}$), contradicting naive factorization.
• Specific Predictions:
  • The model predicts branching fractions for unobserved modes in the range of $10^{-3} \sim 10^{-2}$.
  • It predicts specific polarization patterns for modes like $D^0 \to \rho^0 \rho^0$ (agreeing with FOCUS data: $f_L \approx 77\%$) and $D^0 \to \omega \rho^0$ (where $f_\parallel > f_L$).

5. Significance

• Phenomenological Framework: This work establishes the FAT approach as a robust, systematic tool for analyzing charmed meson decays into vector mesons, overcoming the limitations of naive factorization and strict $SU(3)$ symmetry.
• Experimental Guidance: The paper provides a comprehensive set of predictions for upcoming experiments (BESIII, STCF, Belle II, LHCb). It specifically highlights modes that are currently unmeasured or where current data contradicts simple models, guiding future searches for D-wave dominance and transverse polarization anomalies.
• Understanding QCD Dynamics: By isolating non-factorizable contributions through a global fit, the study offers deeper insight into the strong interaction dynamics (final state interactions and $SU(3)$ breaking) governing the charm sector, which is crucial for understanding CP violation and mixing in the Standard Model.