Impact of the $a_1(1260) \pi$ cascade contribution on $D^0 \to \pi^+ \pi^- \ell^+ \ell^-$ decays

This paper revises the Standard Model description of the rare decays $D^0 \to \pi^+ \pi^- \ell^+ \ell^-$ by incorporating the previously overlooked $a_1(1260)\pi$ cascade contribution, which significantly enhances the predicted decay rate and achieves unprecedented agreement with LHCb data while maintaining consistency with hadronic parameters from analogous four-body decays.

The Gist

Imagine the universe as a giant, chaotic dance floor where tiny particles called "quarks" are constantly swapping partners and spinning around. Most of the time, we know the rules of this dance very well (this is called the Standard Model). But sometimes, the dancers do something unexpected, and physicists get excited because it might mean there's a new, hidden dancer on the floor (called "New Physics") that we haven't seen yet.

This paper is about a specific dance move involving a particle called the D0 meson. Scientists have been watching this particle decay (break apart) into four smaller pieces: two pions (like tiny balls) and two leptons (like electrons or muons).

Here is the simple breakdown of what the authors did and found:

1. The Missing Piece of the Puzzle

For a long time, scientists tried to predict how this D0 meson breaks apart. They had a good map of the dance floor, but when they compared their map to the actual footage from the LHCb experiment (a giant particle detector), the numbers didn't quite match. It was like trying to predict the path of a bouncing ball, but the ball kept landing in a spot your math said was impossible.

The authors realized they were missing a specific "cascade" move.

• The Old Way (The Direct Jump): They thought the D0 meson just broke apart directly into the final pieces.
• The New Way (The Cascade): They realized the D0 meson actually takes a two-step detour. It first turns into a heavy, unstable "middleman" particle called the a1(1260). This middleman then quickly breaks into a rho particle and a pion. Finally, the rho particle breaks into the two leptons we see.

Think of it like a relay race. The old model thought the runner just sprinted from start to finish. The new model realizes the runner actually passed a baton to a teammate (the a1), who then passed it to another teammate (the rho), who finally crossed the finish line.

2. Why This Matters

When the authors added this "relay race" (cascade) move to their calculations, everything clicked.

• The Fit: Their new prediction matched the experimental data almost perfectly. It was like finally solving a jigsaw puzzle where the last piece was the one you were holding upside down.
• The Size: This "relay race" isn't a tiny, rare side effect. It turns out to be one of the biggest contributors to the whole process. It's as important as the main events everyone was already watching.

3. The "Hidden" Signals

The most exciting part is what this new move does to the "angles" of the dance.

• In the old model, certain angles of the particles' movement were predicted to be perfectly flat or zero. It was like saying, "No matter how you spin, you will always face North."
• With the new cascade move, the authors predict that these angles will now tilt. They will point in specific, non-zero directions.
• Why is this cool? If future experiments see these angles tilting exactly as predicted, it confirms our understanding of the Standard Model is solid. If the angles tilt in a different way than predicted, that would be a smoking gun for "New Physics"—a sign that a new, unknown force is interfering with the dance.

4. Checking the Dance Floor

To make sure they weren't just making up numbers, the authors compared their results to other types of particle decays (where the D0 meson breaks into four pions instead of leptons).

• They found that the "relay race" move (the cascade) is just as popular in those other dances as it is in the one they studied.
• This consistency suggests their model is robust and that they are correctly describing how these particles interact, even when they are doing something complex and messy.

The Bottom Line

The authors didn't discover a new particle. Instead, they realized they were ignoring a very common, complex step in the dance routine. By adding this step back in, they fixed the math, matched the data perfectly, and created a new, more sensitive tool (the angular observables) to catch any future "New Physics" that might try to sneak onto the dance floor.

In short: They found the missing step in the dance, fixed the choreography, and now have a better way to spot if a ghost is dancing with them.

Technical Summary

Technical Summary: Impact of the $a_1(1260)\pi$ Cascade Contribution on $D^0 \to \pi^+\pi^-\ell^+\ell^-$ Decays

Problem Statement
Recent experimental measurements of the rare charm-meson decay $D^0 \to \pi^+\pi^-\ell^+\ell^-$ by LHCb have revealed tensions with Standard Model (SM) predictions, particularly regarding decay rate distributions over the dipion mass and certain CP-symmetric angular observables. Previous theoretical analyses, relying on quasi-two-body (Q2B) topologies mediated by vector ($\rho^0, \omega$) and scalar ($f_0$) resonances, provided a qualitative description but failed to fully account for systematic deviations in the data. Furthermore, specific angular observables that theoretically vanish in the SM (null tests) were measured with non-zero values, suggesting either new physics (NP) or unaccounted-for long-distance SM contributions. The authors identify a missing component in the theoretical framework: the cascade-type topology involving the axial-vector resonance $a_1(1260)$.

Methodology
The authors extend their previous SM calculation by incorporating a cascade topology where the $D^0$ meson decays weakly into an intermediate $a_1(1260)^+$ and a $\pi^-$, followed by the strong decay $a_1^+ \to \rho^0 \pi^+$ and the electromagnetic decay $\rho^0 \to \ell^+\ell^-$.

• Theoretical Framework: The effective Hamiltonian for $c \to u \mu^+\mu^-$ is utilized, focusing on the four-quark operators $Q_1$ and $Q_2$. The calculation is performed in the large-$N_C$ limit, where the cascade amplitude is dominated by the insertion of the $Q_1$ operator (associated with Wilson coefficient $C_1$), which is approximately three times larger than $C_2$ governing the leading Q2B topologies.
• Amplitude Construction: The decay amplitude is modeled using an isobar-like approach. The $a_1 \to \rho \pi$ transition is parameterized with a constant coupling, utilizing a data-driven lineshape for the $a_1$ resonance based on CLEO-c amplitude analyses. The $D \to \pi$ form factors are taken from lattice QCD calculations.
• Kinematics: Unlike Q2B topologies, the cascade topology introduces a distinct kinematic dependence where the squared momentum of the axial resonance ($k^2$) depends on the dipion invariant mass ($p^2$), the dilepton invariant mass ($q^2$), and the angle $\theta_\pi$. This leads to a non-trivial dependence of transversity form factors on $\theta_\pi$.
• Fitting Procedure: The model parameters, including normalization factors and relative strong phases, are fitted to the binned differential branching ratios of $D^0 \to \pi^+\pi^-\mu^+\mu^-$ provided by LHCb. The fit covers the $p^2$ range from threshold to $1.0 \text{ GeV}^2$ and $q^2$ from $0.5$ to $1.2 \text{ GeV}^2$, excluding regions with potential contamination or poorly known heavy resonances.

Key Contributions and Results

1. Improved Data Description: The inclusion of the $a_1\pi$ cascade contribution significantly improves the agreement between the SM prediction and LHCb data. The $\chi^2_{\text{min}}$ per degree of freedom drops from approximately 2 in previous analyses to roughly 1 (specifically, $\chi^2_{\text{min}} \approx 82$ for $\sim 82$ degrees of freedom, corresponding to a $p$-value of $\sim 48\%$). The improvement is most pronounced in the description of the $p^2$ differential branching ratio, particularly in the region above the $\rho^0$ resonance where the cascade contribution exhibits a long tail.
2. Fit Fractions: The fitted contribution of the cascade topology is substantial, accounting for approximately 40% of the total decay rate, comparable to the contributions from the $\sigma$ and $\rho^0$ resonances. The interference between the cascade and other components (specifically $\sigma \to \pi\pi$ and $\rho^0 \to \pi\pi$) is found to be destructive in the low-$p^2$ region.
3. Angular Observables: The presence of the cascade topology modifies the angular structure of the decay.
   • It introduces non-zero values for observables that previously vanished in the SM, such as $\langle I_7 \rangle_-$, $\langle I_3 \rangle_-$, and $\langle I_9 \rangle_-$.
   • Specifically, the interference between the cascade and Q2B topologies allows for non-vanishing relative phases between transversity form factors of the same partial wave, enabling non-zero values for $S_8$ and $S_9$.
   • The authors note that while the SM prediction for $\langle I_7 \rangle_-$ remains small ($<1\%$), the observable is sensitive to NP-induced $C_{10}$ and could be enhanced in specific kinematic regions.
4. Consistency with Non-Leptonic Decays: The normalization factor extracted for the cascade contribution in the rare decay ($ga_{1\rho\pi}f_{a_1}B_{\text{casc}}$) is found to be consistent with values extracted from amplitude analyses of non-leptonic four-body decays ($D^0 \to 4\pi$ and $D^0 \to 2\pi 2K$). This supports the universality of hadronic effects in charm decays and validates the inclusion of this topology.

Significance
The paper claims that the $a_1(1260)\pi$ cascade contribution is one of the largest components of the $D^0 \to \pi^+\pi^-\ell^+\ell^-$ decay rate, a feature previously unaccounted for in rare decay analyses. By incorporating this topology, the authors achieve an unprecedented agreement between SM predictions and experimental data, resolving previous tensions in decay rate distributions.

Crucially, the work demonstrates that the non-zero values observed in certain angular observables can be explained within the SM through the interference of the cascade topology with Q2B amplitudes, rather than necessarily implying New Physics. The authors argue that this cascade component is phenomenologically relevant and essential for setting meaningful bounds on NP-driven Wilson coefficients. They conclude that future experimental analyses, particularly those measuring the $d\Gamma/d\cos\theta_\pi$ distribution and optimized angular observables, are necessary to further test the hadronic dynamics and constrain potential new physics contributions.