Nodal mechanism for the suppressed $D\bar D$ decay of $\psi(4040)$ in the Bethe--Salpeter framework

This paper explains the anomalously suppressed $D\bar D$ decay of $\psi(4040)$ within a conventional charmonium framework by demonstrating that node-induced cancellations in the relativistic decay amplitude significantly reduce the $D\bar D$ width while leaving other open-charm channels largely unaffected.

The Gist

Imagine a tiny, heavy particle called $\psi(4040)$. Think of it as a very excited, vibrating drum made of heavy quarks (particles that make up protons and neutrons). When this drum vibrates, it wants to break apart into smaller pieces, specifically pairs of particles called $D$ mesons.

Physicists have a puzzle: When this drum breaks, it almost always splits into a specific combination of pieces ($D\bar{D}^*$ or $D_s\bar{D}_s$). However, there is one combination ($D\bar{D}$) that it should be able to split into easily, but it almost never does. It's like a door that is wide open, but the person inside refuses to walk through it.

This paper explains why that door stays closed.

The "Bumpy Road" Analogy

To understand the solution, imagine the particle's internal structure isn't a smooth ball, but a bumpy road or a wavy ocean.

1. The Wave: Because the $\psi(4040)$ is an excited state (like a drum hit hard), its internal "wave" has a special shape. It goes up, comes down, crosses the zero line, goes down, and comes back up. This crossing point is called a node.
2. The Journey: When the particle decays (breaks apart), it has to "travel" through different speeds (momenta) to create the new pieces.
3. The Cancellation:
   • For the forbidden path ($D\bar{D}$), the journey takes the particle through a section of the wave where the "up" bumps and "down" bumps are perfectly balanced.
   • Imagine walking on a path where every step forward is canceled out by a step backward. You end up going nowhere. In physics terms, the positive and negative parts of the calculation cancel each other out, making the result zero.
   • For the allowed paths ($D\bar{D}^*$ and $D_s\bar{D}_s$), the journey takes them through different parts of the wave where the bumps don't cancel out. They keep moving forward, so these decays happen frequently.

The "Filter" Metaphor

The authors describe this process as a momentum filter.

• The particle's internal structure acts like a sieve.
• The "forbidden" pieces ($D\bar{D}$) are the exact size that fits perfectly into the holes of the sieve, getting filtered out (cancelled).
• The "allowed" pieces are slightly different sizes; they don't fit the holes and pass right through.

The "Tuning Fork" Sensitivity

The paper also points out something very interesting about how sensitive this cancellation is.

• Because the "forbidden" path is so close to being perfectly cancelled out, it is extremely sensitive to tiny changes.
• The authors tested this by slightly changing the "weight" (mass) of the starting particle in their computer model.
• The Result: A tiny shift in weight caused the "forbidden" path to suddenly open up or close down dramatically. It's like a tightrope walker balancing on a wire; a tiny breeze (mass change) makes them fall or stand up.
• In contrast, the "allowed" paths were stable and didn't care much about the tiny weight change. This proves that the suppression isn't just a random accident; it's a specific feature of the wave's shape.

The "Isospin" Twist

The paper also looked at the difference between "charged" and "neutral" versions of the $D$ mesons.

• Normally, the difference between a charged and neutral particle is tiny, like the difference between a red apple and a slightly redder apple.
• However, because the "forbidden" path is already balanced on the edge of zero, this tiny difference gets amplified. It's like a microphone that is turned up so high that even a whisper sounds like a shout. The tiny difference in mass causes a noticeable difference in how often the decay happens, but only because the main signal was already so weak.

The Conclusion

The authors used a sophisticated mathematical framework (the Bethe-Salpeter equation combined with a model called $^3P_0$) to prove that:

1. The $\psi(4040)$ is a standard, conventional particle (a "3S" state).
2. It doesn't need to be a mysterious new type of matter to explain why it doesn't decay into $D\bar{D}$.
3. The reason is purely mathematical geometry: The internal wave of the particle has a "node" (a zero point) that lines up perfectly to cancel out the $D\bar{D}$ decay, while leaving the other decays untouched.

In short, the particle isn't "choosing" not to decay; the laws of physics and the shape of its internal wave make it mathematically impossible for that specific decay to happen, while allowing others to proceed freely.

Technical Summary

Technical Summary: Nodal Mechanism for the Suppressed $\bar{D}D$ Decay of $\psi(4040)$

Problem Statement
The vector charmonium state $\psi(4040)$ exhibits a long-standing anomaly in its strong decay patterns: despite having ample phase space, the decay channel $\psi(4040) \to D\bar{D}$ is anomalously suppressed. Conversely, the $D\bar{D}^*$ and $D_s\bar{D}_s$ channels remain sizable, creating a hierarchy ($\Gamma(D\bar{D}^*) \gg \Gamma(D_s\bar{D}_s) \gg \Gamma(D\bar{D})$) that deviates from naive phase-space expectations. While previous theoretical studies (e.g., nonrelativistic $^3P_0$ models and potential models) have qualitatively attributed this suppression to nodal effects in the radial wave function of a $3\,^3S_1$ assignment, the precise amplitude-level origin remained unclear. Specifically, it was not explicitly demonstrated how the node distributes across relativistic Salpeter components, which quantity changes sign to cause cancellation, and why this cancellation affects $D\bar{D}$ but not $D\bar{D}^*$ or $D_s\bar{D}_s$.

Methodology
The authors investigate this hierarchy within the framework of the instantaneous Bethe–Salpeter (BS) equation combined with the relativistic $^3P_0$ model.

• Bound State Formalism: The initial $\psi(4040)$ state and final heavy-light mesons ($D, D^*, D_s$) are described by relativistic Salpeter wave functions. The BS equation is solved using a screened Cornell potential kernel. For the vector $1^{--}$ state, the wave function comprises eight scalar components, reduced to four independent ones ($f_3, f_4, f_5, f_6$) by constraint equations.
• Decay Mechanism: The OZI-allowed strong decays are modeled by the creation of a light quark-antiquark pair from the vacuum with quantum numbers $J^{PC}=0^{++}$. The transition amplitude is calculated using the relativistic $^3P_0$ interaction Hamiltonian.
• Parameter Fixing: Crucially, the pair-creation strength parameter $\gamma$ is not fitted to $\psi(4040)$ data. Instead, it is determined independently by reproducing the experimental decay width of $\psi(3770) \to D\bar{D}$. This fixed parameter is then applied to $\psi(4040)$ without further adjustment.
• Mass Handling: The calculation utilizes a "diagnostic setup" where the BS eigenstate for the $3\,^3S_1$ charmonium is found at $4051$ MeV, while the physical phase space is evaluated at the experimental pole mass of $4040$ MeV. The authors explicitly scan the initial mass $M_A$ in the range $4.00\text{--}4.07$ GeV to isolate the dynamical nodal mechanism from kinematic mass shifts.

Key Contributions and Results

1. Amplitude-Level Cancellation: The study demonstrates that the suppression of the $D\bar{D}$ mode is not due to the vanishing of a single radial component, but rather a sign-changing cancellation in the complete channel-dependent overlap integrand. In the relativistic BS framework, the node is distributed among multiple Salpeter components. For the $D\bar{D}$ channel, the overlap integral receives comparable positive and negative contributions from different momentum regions, leading to a near-total destructive interference.
2. Channel Selectivity: The same initial-state wave function does not produce a similar cancellation for the $D\bar{D}^*$ or $D_s\bar{D}_s$ channels.
   • For $D_s\bar{D}_s$, the overlap integral probes a different momentum region of the initial wave function, avoiding the nodal zero.
   • For $D\bar{D}^*$, the decay width depends on the incoherent sum $|M_{12}|^2 + |M_{21}|^2$ of two independent amplitude structures. Even if these structures have opposite signs, they do not cancel in the physical rate, allowing the channel to remain sizable.
3. Numerical Hierarchy: Using the fixed $\gamma$ derived from $\psi(3770)$, the calculated partial widths reproduce the experimental hierarchy:
   • $\Gamma(D\bar{D}) \approx 0.028$ MeV (strongly suppressed).
   • $\Gamma(D\bar{D}^* + \text{c.c.}) \approx 26.8$ MeV (dominant).
   • $\Gamma(D_s\bar{D}_s) \approx 7.9$ MeV (sizable).
4. Diagnostic Signatures: The paper identifies two correlated signatures supporting the nodal mechanism:
   • Dip Structure: A mass scan reveals a pronounced dip in the $\Gamma(D\bar{D})$ curve near the physical mass, confirming the physical region lies close to a zero of the overlap amplitude. In contrast, $\Gamma(D_s\bar{D}_s)$ and $\Gamma(D\bar{D}^*)$ vary smoothly.
   • Isospin Sensitivity: The small mass splitting between charged ($D^\pm$) and neutral ($D^0$) mesons is strongly amplified in the $D\bar{D}$ width because the amplitude is already near zero. The authors show that in the isospin limit, the width becomes even smaller, indicating that the physical difference arises from the amplification of small kinematic shifts near a nodal zero, rather than strong isospin violation.

Significance and Limitations
The paper claims to provide a dynamical explanation for the suppressed $D\bar{D}$ mode and the open-charm hierarchy of $\psi(4040)$ within a conventional $3\,^3S_1$ charmonium picture. The significance lies in moving beyond qualitative nodal arguments to a quantitative, amplitude-level demonstration of how relativistic overlap integrals filter specific decay channels.

However, the authors maintain a modest stance regarding quantitative precision:

• The absolute value of the near-vanishing $D\bar{D}$ width is model-dependent and numerically fragile. A shift of only $11$ MeV (between the BS eigenmass and the experimental pole) changes the $D\bar{D}$ width by an order of magnitude.
• The study is restricted to the three primary channels ($D\bar{D}, D\bar{D}^*, D_s\bar{D}_s$) and does not include the experimentally important $D^*\bar{D}^*$ channel or coupled-channel dynamics.
• The results confirm that the suppression mechanism is robust, but the precise numerical value of the $D\bar{D}$ width should be interpreted with caution due to its sensitivity to the mass mismatch.

In conclusion, the work illustrates that the open-flavor decays of radially excited quarkonia are controlled not only by phase space and spin counting but also by the sign structure of relativistic overlap integrals, offering a clear example of "nodal filtering" in charmonium physics.