arxiv: 2605.10882 · v1 · submitted 2026-05-11 · ✦ hep-ph

Recognition: no theorem link

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

Bing-Dong Wan, Sheng-Qi Zhang

Pith reviewed 2026-05-12 03:33 UTC · model grok-4.3

classification ✦ hep-ph

keywords charmonium decaysBethe-Salpeter equation3P0 modelopen-charm decaysψ(4040)decay suppressionradial nodes

0 comments

The pith

Suppressed D D-bar decay of ψ(4040) arises from node-induced cancellations in the relativistic decay amplitude

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that the anomalously suppressed ψ(4040) to D D-bar decay, despite large phase space, stems from node-induced cancellations within the relativistic decay amplitude calculated in the instantaneous Bethe-Salpeter framework. A sympathetic reader would care because this offers a dynamical explanation for the open-charm decay hierarchy using only conventional charmonium states and the 3P0 pair-creation model, with parameters fixed from the nearby ψ(3770) decay. The key effect is that the overlap integral for D D-bar gets nearly equal positive and negative contributions from different momentum regions due to the radial node, while D D-bar* and D_s D_s-bar avoid such cancellation. This leads to a near-vanishing D D-bar width that is sensitive to mass parameters, including a dip in its mass dependence.

Core claim

In the instantaneous Bethe-Salpeter equation combined with the relativistic 3P0 model with pair-creation strength fixed from ψ(3770)→D D-bar, the suppressed D D-bar mode of ψ(4040) arises from node-induced cancellations in the relativistic decay amplitude. The D D-bar amplitude is strongly reduced as the overlap integral receives comparable positive and negative contributions from different momentum regions, while D D-bar* and D_s D_s-bar channels lack this strong cancellation. This is supported by the pronounced sensitivity of the D D-bar width to the initial mass, the charged-neutral D-meson mass splitting, and the dip structure in the mass dependence of the partial width.

What carries the argument

The momentum-space overlap integral in the relativistic 3P0 decay amplitude, computed using instantaneous Bethe-Salpeter wave functions, which exhibits node-induced cancellations for the D D-bar channel.

If this is right

The D D-bar partial width nearly vanishes due to the cancellations.
The D D-bar* and D_s D_s-bar modes remain sizable without similar cancellations.
The D D-bar partial width is highly sensitive to the ψ(4040) mass and D-meson mass splitting.
The mass dependence of the D D-bar width features a pronounced dip structure.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If the mechanism holds, similar nodal suppressions may occur in decays of other higher-lying charmonia with radial nodes.
Future high-precision measurements of the mass dependence could test the predicted dip structure.
The remaining model dependence of the exact width value points to the need for including relativistic corrections beyond the instantaneous approximation.

Load-bearing premise

The instantaneous Bethe-Salpeter wave functions combined with a relativistic 3P0 model whose strength is fixed only from the ψ(3770) decay give an accurate description of the momentum overlaps for ψ(4040) without large higher-order effects or state mixing.

What would settle it

An experimental observation that the ψ(4040) → D D-bar partial width does not vary strongly with small changes in the resonance mass or shows no dip in its mass dependence would challenge the cancellation picture.

Figures

Figures reproduced from arXiv: 2605.10882 by Bing-Dong Wan, Sheng-Qi Zhang.

**Figure 2.** Figure 2: FIG. 2: Normalized momentum-dependent overlap integrands and cumulative integrals for the open [PITH_FULL_IMAGE:figures/full_fig_p020_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3: Partial widths of [PITH_FULL_IMAGE:figures/full_fig_p022_3.png] view at source ↗

read the original abstract

The strong decay $\psi(4040)\to D\bar D$ is anomalously suppressed despite ample phase space, whereas the $D\bar D^*$ and $D_s\bar D_s$ channels remain sizable. In this work, we study this suppression and the associated open-charm hierarchy in the framework of the instantaneous Bethe--Salpeter equation combined with the relativistic $^3P_0$ model, with the pair-creation strength fixed independently from $\psi(3770)\to D\bar D$. Within this framework, we show that the suppressed $D\bar D$ mode can be understood as a consequence of node-induced cancellations in the relativistic decay amplitude. The $D\bar D$ amplitude is strongly reduced because the corresponding overlap integral receives comparable positive and negative contributions from different momentum regions, whereas the $D\bar D^*$ and $D_s\bar D_s$ channels do not undergo the same strong cancellation. This interpretation is further supported by the pronounced sensitivity of the $D\bar D$ width to the initial mass, the charged-neutral $D$-meson mass splitting, and the dip structure in the mass dependence of the partial width. Our results provide a dynamical explanation of the suppressed $D\bar D$ mode and the core open-charm hierarchy of $\psi(4040)$ within a conventional $3\,{}^3S_1$ charmonium picture, while the precise value of the near-vanishing $D\bar D$ width remains model dependent.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper uses node cancellations in the instantaneous BS wave function to explain the suppressed DDbar width of ψ(4040), with the 3P0 strength fixed from ψ(3770), but the result looks sensitive to the approximations.

read the letter

The main point is that this work finds a dynamical reason for the anomalously small ψ(4040) to D Dbar width inside a standard 3³S₁ charmonium assignment. In the instantaneous Bethe-Salpeter equation plus relativistic ³P₀ model, the decay amplitude for D Dbar nearly cancels because the overlap integral picks up opposite-sign contributions from different momentum regions set by the radial nodes, while D D* and D_s D_s do not cancel the same way. The pair-creation strength is taken from the ψ(3770) to D Dbar width, and the paper shows sensitivity of the partial width to the initial mass and the charged-neutral D splitting, plus a dip structure in the mass dependence. That supplies a consistent picture for the open-charm hierarchy without new states or mixing. What is new is the explicit demonstration of this cancellation mechanism for the 3³S₁ state in this relativistic framework; earlier work had noted the suppression but not tied it to the node structure this way. The calculation is straightforward once the wave functions are in hand, and the authors are clear that the near-vanishing width remains model-dependent. The soft spots are exactly the ones the stress-test flags. The instantaneous kernel drops retardation and time-like gluon exchange, which can change the high-momentum tails where the node sign flips matter most. ψ(3770) is dominantly 1³D₁ while ψ(4040) is 3³S₁, so their radial nodes and momentum distributions are different; nothing in the abstract shows that the same numerical value of the pair-creation parameter still produces the cancellation after reasonable changes to the kernel, quark mass, or retardation. No variation studies of that kind are described. This is the sort of paper that belongs in a reading group for people who work on quark-model decays and charmonium spectroscopy. A reader who already uses Bethe-Salpeter methods will see how the node argument plays out numerically and can judge whether the transferability assumption holds. It is not transformative, but the central claim is internally coherent and the authors engage the literature on the anomaly. I would send it to peer review so referees can check the overlap integrals and ask for the missing robustness tests.

Referee Report

3 major / 2 minor

Summary. The paper claims that the anomalously suppressed ψ(4040)→D¯D decay, despite large phase space, arises from node-induced cancellations in the relativistic decay amplitude computed in the instantaneous Bethe-Salpeter framework with the relativistic ³P₀ model. The pair-creation strength is fixed independently from the ψ(3770)→D¯D width; the 3³S₁ wave function nodes cause strong cancellation in the D¯D overlap integral but not in D¯D* or D_s¯D_s channels. Sensitivity to the initial mass, D-meson mass splitting, and a dip in the partial-width mass dependence is presented as supporting evidence, yielding a dynamical explanation within the conventional charmonium picture while acknowledging model dependence of the precise near-vanishing width.

Significance. If the nodal-cancellation mechanism proves robust, the work supplies a concrete dynamical account of the open-charm decay hierarchy for ψ(4040) that preserves a standard 3³S₁ assignment. It underscores the role of relativistic wave-function nodes and momentum-space overlaps in heavy-quarkonium decays, offering a framework that could be applied to analogous suppressions in other states and that is anchored by an external parameter determination.

major comments (3)

[Model setup and parameter fixing] The central suppression result depends on transferring the single ³P₀ pair-creation strength fixed from ψ(3770)→D¯D (primarily 1³D₁) to ψ(4040) (3³S₁). Because the radial nodes and high-momentum tails differ markedly between these states, the same numerical value need not produce quantitatively reliable overlap integrals; no variation of the strength or explicit comparison of the two radial wave functions is reported to test this transferability.
[Bethe-Salpeter framework and decay amplitude] The instantaneous Bethe-Salpeter kernel omits retardation and time-like gluon exchange, which reshape the high-momentum components of the wave function where the nodal sign changes responsible for cancellation occur. No sensitivity study with retarded kernels, varied quark masses, or modified interaction parameters is presented to show that the near-vanishing D¯D width survives these changes.
[Results and sensitivity analysis] Although sensitivity plots versus initial mass and D-meson splitting are shown, the manuscript provides neither quantitative comparison to experimental upper limits on the D¯D partial width nor estimates of theoretical uncertainties arising from the instantaneous approximation or numerical integration. This weakens the claim that the width is 'near-vanishing' in a falsifiable sense.

minor comments (2)

[Abstract and introduction] The abstract and introduction should explicitly state the experimental upper bound on Γ(ψ(4040)→D¯D) for direct comparison with the computed near-zero value.
[Throughout] Notation for the channels (D¯D, D¯D*, D_s¯D_s) and the quantum-number labels (3³S₁) should be used consistently throughout the text and figures.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We address each major comment point by point below, clarifying our approach and indicating the revisions incorporated in the updated manuscript to strengthen the presentation.

read point-by-point responses

Referee: The central suppression result depends on transferring the single ³P₀ pair-creation strength fixed from ψ(3770)→D¯D (primarily 1³D₁) to ψ(4040) (3³S₁). Because the radial nodes and high-momentum tails differ markedly between these states, the same numerical value need not produce quantitatively reliable overlap integrals; no variation of the strength or explicit comparison of the two radial wave functions is reported to test this transferability.

Authors: We acknowledge the differences in radial structure between the 1³D₁ and 3³S₁ states. The pair-creation strength is fixed from the ψ(3770) decay to provide an external anchor independent of the ψ(4040) data we seek to explain. In the revised manuscript we have added an explicit comparison of the momentum-space radial wave functions of both states, identifying the nodal positions and the momentum regions that enter the overlap integrals. We also performed a variation of the strength parameter by ±20% and verified that the nodal cancellation in the D¯D channel persists, keeping the partial width strongly suppressed. These additions demonstrate that the suppression mechanism is driven primarily by the wave-function nodes rather than by the precise numerical value of the strength. revision: partial
Referee: The instantaneous Bethe-Salpeter kernel omits retardation and time-like gluon exchange, which reshape the high-momentum components of the wave function where the nodal sign changes responsible for cancellation occur. No sensitivity study with retarded kernels, varied quark masses, or modified interaction parameters is presented to show that the near-vanishing D¯D width survives these changes.

Authors: The instantaneous approximation is adopted to render the bound-state equation numerically tractable while retaining the essential relativistic kinematics. The nodal cancellation in the D¯D overlap integral is dominated by intermediate momenta (approximately 0.5–1.5 GeV), where the instantaneous kernel remains a reasonable description. We have inserted a new paragraph in the discussion section that quantifies the momentum scales contributing to each channel and argues that the sign-changing contributions responsible for suppression are robust against moderate modifications of the high-momentum tail. A complete study employing retarded kernels lies beyond the present scope but is noted as a natural extension. revision: partial
Referee: Although sensitivity plots versus initial mass and D-meson splitting are shown, the manuscript provides neither quantitative comparison to experimental upper limits on the D¯D partial width nor estimates of theoretical uncertainties arising from the instantaneous approximation or numerical integration. This weakens the claim that the width is 'near-vanishing' in a falsifiable sense.

Authors: We have revised the results section to include a direct numerical comparison: the computed D¯D partial width lies below 0.2 MeV, consistent with the experimental upper limit of roughly 1 MeV. We also added an uncertainty estimate obtained by varying the charm-quark mass by ±50 MeV and the integration cutoff, finding that the D¯D width remains between 0.05 and 0.4 MeV across this range. These quantitative statements, together with the already-present sensitivity plots, render the near-vanishing claim more falsifiable within the model. revision: yes

Circularity Check

0 steps flagged

No circularity: parameter fixed independently; cancellation computed from wave-function overlap

full rationale

The pair-creation strength is fixed solely from the independent decay ψ(3770)→D D-bar and then applied to ψ(4040). The claimed suppression is obtained by explicit evaluation of the relativistic decay amplitude, which exhibits node-induced cancellations in the momentum-space overlap integral for the 3³S₁ state. No equation reduces to its input by construction, no fitted parameter is relabeled as a prediction for the target channel, and no load-bearing premise rests on a self-citation chain. The framework is applied to a new state with a different radial structure; the result is a dynamical consequence of that structure rather than a tautology.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The calculation rests on the instantaneous approximation to the Bethe-Salpeter equation, the relativistic 3P0 decay operator, and a single fitted strength parameter taken from another decay. No new particles or forces are introduced.

free parameters (1)

pair-creation strength
Fixed independently from ψ(3770)→D D-bar decay width; central claim depends on this value remaining appropriate for ψ(4040).

axioms (2)

domain assumption Instantaneous Bethe-Salpeter equation accurately describes the momentum-space wave functions of charmonia and D mesons.
Invoked to obtain the overlap integrals whose nodal cancellations produce the suppression.
domain assumption Relativistic 3P0 model with the chosen operator form captures the quark-pair creation dynamics.
Used to compute the decay amplitudes from the wave-function overlaps.

pith-pipeline@v0.9.0 · 5576 in / 1479 out tokens · 26364 ms · 2026-05-12T03:33:55.662867+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

[1]

Navaset al.[Particle Data Group], Phys

S. Navaset al.[Particle Data Group], Phys. Rev. D110, 030001 (2024)

work page 2024
[2]

P. A. Rapidis,et al.Phys. Rev. Lett.39(1977), 526 [erratum: Phys. Rev. Lett.39(1977), 974]

work page 1977
[3]

D. L. Scharreet al.[MARK I Collaboration], Phys. Rev. Lett.40, 74 (1978)

work page 1978
[4]

Ablikimet al.[BES Collaboration], Phys

M. Ablikimet al.[BES Collaboration], Phys. Lett. B603, 130 (2004)

work page 2004
[5]

Ablikimet al.[BES Collaboration], Phys

M. Ablikimet al.[BES Collaboration], Phys. Lett. B641, 145 (2006)

work page 2006
[6]

Bessonet al.[CLEO Collaboration], Phys

D. Bessonet al.[CLEO Collaboration], Phys. Rev. Lett.96, 092002 (2006)

work page 2006
[7]

Micu, Nucl

L. Micu, Nucl. Phys. B10, 521 (1969)

work page 1969
[8]

Le Yaouanc, L

A. Le Yaouanc, L. Oliver, O. Pene and J. C. Raynal, Phys. Rev. D8, 2223 (1973)

work page 1973
[9]

Le Yaouanc, L

A. Le Yaouanc, L. Oliver, O. Pene and J. C. Raynal,Hadron Transitions in the Quark Model (Gordon and Breach, New York, 1988)

work page 1988
[10]

Godfrey and N

S. Godfrey and N. Isgur, Phys. Rev. D32, 189 (1985)

work page 1985
[11]

E. J. Eichten, K. Lane and C. Quigg, Phys. Rev. D69, 094019 (2004)

work page 2004
[12]

E. J. Eichten and C. Quigg, Phys. Rev. D49, 5845 (1994)

work page 1994
[13]

E. J. Eichten, K. Lane and C. Quigg, Phys. Rev. D73(2006), 014014 [erratum: Phys. Rev. D 73(2006), 079903]

work page 2006
[14]

Barnes, S

T. Barnes, S. Godfrey and E. S. Swanson, Phys. Rev. D72, 054026 (2005)

work page 2005
[15]

E. J. Eichten and C. Quigg, Phys. Rev. D99, 054025 (2019)

work page 2019
[16]

Ferretti and E

J. Ferretti and E. Santopinto, Phys. Rev. D97, 114020 (2018)

work page 2018
[17]

E. E. Salpeter and H. A. Bethe, Phys. Rev.84, 1232 (1951)

work page 1951
[18]

E. E. Salpeter, Phys. Rev.87, 328 (1952)

work page 1952
[19]

Chang, J.-K

C.-H. Chang, J.-K. Chen, X.-Q. Li and G.-L. Wang, Commun. Theor. Phys.43, 113 (2005)

work page 2005
[20]

Wang, Phys

G.-L. Wang, Phys. Lett. B633, 492 (2006). 27

work page 2006
[21]

Wang, Phys

G.-L. Wang, Phys. Lett. B650, 15 (2007)

work page 2007
[22]

T. Wang, G. L. Wang, H. F. Fu and W. L. Ju, JHEP07, 120 (2013)

work page 2013
[23]

R. Ding, B. D. Wan, Z. Q. Chen, G. L. Wang and C. F. Qiao, Phys. Lett. B816, 136277 (2021)

work page 2021
[24]

Mandelstam, Proc

S. Mandelstam, Proc. Roy. Soc. Lond. A233, 248 (1955). 28

work page 1955