Coupled-channel study of $4S$-$3D$ mixing dynamics in $\psi(4220)$ and $\psi(4380)$

Si-Qiang Luo; Xiang Liu; Zi-Long Man; Zi-Yue Bai

REVIEW 4 major objections 4 minor 1 cited by

Coupled-channel study of $4S$-$3D$ mixing dynamics in $\psi(4220)$ and $\psi(4380)$

T0 review · 4 major / 4 minor · reviewed 2026-08-08 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A 35° 4S-3D mixing angle generated by coupled-channel effects identifies ψ(4220) as a mixed charmonium state with mass 4235.9 MeV and width 32.0 MeV, and predicts a ψ(4380) partner at 4387.1 MeV.

desk verdict Coupled-channel S–D mixing for ψ(4220) is a plausible mechanism, but the headline 35° angle is not yet demonstrated because the off-diagonal shift comes with no sensitivity analysis. read the letter →

arxiv 2502.08072 v2 pith:2WNP2MXF submitted 2025-02-12 hep-ph hep-ex

Zi-Long Man , Si-Qiang Luo , Zi-Yue Bai , Xiang Liu This is my paper · ORCID

classification hep-phhep-ex

keywords charmoniumψ(4220)ψ(4380)4S-3Dmixingcoupled-channelmodelhadronicloopsOZI-allowedstrongdecaysXYZstates

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper tries to establish that ψ(4220) is the lower member of a 4S-3D mixed charmonium pair, with the large mixing angle generated dynamically by coupled-channel effects rather than by the conventional tensor force. If true, this resolves the residual mass discrepancy left by a pure 4S state, which still sits about 60 MeV too high after loop corrections: the mixed state lands at 4235.9 MeV with a width of 32.0 MeV, close to the measured ψ(4220). The same scheme predicts a partner ψ(4380) at 4387.1 MeV with a 44.1 MeV width and distinctive decay channels, so the whole picture is experimentally testable.

What carries the argument

The load-bearing object is the off-diagonal self-energy ΔM_SD(M), a loop integral that connects the ψ(4S) and ψ(3D) bare states through intermediate meson-antimeson pairs. Since the full sum over loops is infinite in principle, the paper computes it with a once-subtracted dispersion relation anchored at the J/ψ mass, and evaluates the loop couplings with a quark-pair-creation model. The argument turns on ΔM_SD being comparable in magnitude to the diagonal mass shifts, which is what elevates the potential-model mixing angle from 0.5° to the large 35° angle needed to identify ψ(4220).

What would settle it

Recompute the mixing angle using a different method of handling the infinite loop sum or a different set of intermediate meson channels; if the angle drops well below 30°, the mechanism fails. Alternatively, search for a narrow vector state near 4387 MeV with a width near 44 MeV in e+e− data; its non-appearance would falsify the partner prediction.

Watch

Extended reading notes

Core claim

The paper's central claim is that the charmonia near 4.2–4.4 GeV are not pure ψ(4S) and ψ(3D) states but mixtures produced by hadronic loops. Starting from bare masses 4433.0 MeV and 4491.3 MeV, the coupled-channel self-energies lower the diagonal masses by roughly 154 MeV each, and the off-diagonal shift becomes comparable in size, so solving the 2×2 mass equation gives a mixing angle θ = 35°. The lower mixed state then sits at 4235.9 MeV with a two-body OZI-allowed strong decay width of 32.0 MeV, matching the measured ψ(4220); the higher mixed state is predicted at 4387.1 MeV with width 44.1 MeV and identified with the narrow ψ(4380) structure.

Load-bearing premise

The large mixing angle rests on a particular way of taming the infinite sum over intermediate meson pairs and on which pairs are included; if that way is arbitrary, the angle itself may be an artifact.

Editorial extensions

If this is right

The ψ(4220) should be regarded as the lower 4S-3D mixed charmonium state, which removes the roughly 60 MeV gap between the corrected ψ(4S) mass and the measured ψ(4220) mass.
A pure ψ(4S) would have a total width of 52.7 MeV, whereas the mixed ψ(4220) is predicted at 32.0 MeV, so the narrower measured width is direct evidence for mixing under this scheme.
The dominant decay mode of ψ(4220) is predicted to be D*D* at 91.4%, so a measurement of e+e− → D*D* near 4.24 GeV can either confirm or exclude the assignment.
A partner state ψ(4380) should exist at 4387.1 MeV with width 44.1 MeV and dominant decays into DD2*(2460), D*D*, DD1(2430), and DD1(2420), making it searchable in existing e+e− data.
The decay pattern of ψ(4380) discriminates the mixed assignment from a pure ψ(3D): the DD2*(2460) branching fraction jumps from 0.9% to 37.3% when the mixing angle is 35°.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Our inference: the same once-subtracted loop machinery should generate sizable mixing angles for other high-radial charmonium pairs, such as 5S-4D, so the mechanism can be tested by whether ψ(4500)-type partners appear with the predicted masses and widths.
Our inference: because the ψ(4220) width of 32.0 MeV relies on the D*D* channel opening just above threshold, precise measurements of the D*D* line shape around 4.24 GeV would provide a sharper test of the 91.4% branching fraction than total-width comparisons alone.
Our inference: if a future calculation with a different subtraction prescription or an expanded loop set yields a much smaller mixing angle, the identification of ψ(4220) as the lower mixed state would lose its dynamical support despite the mass and width agreement.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Dear Colleague,

Read the coupled-channel 4S–3D mixing paper for ψ(4220)/ψ(4380). The headline: it does a credible job of showing that hadronic loops can generate a large S–D mixing angle, and the resulting mass (4235.9 MeV) and width (32.0 MeV) land near the ψ(4220) data. But the central angle θ = 35° is produced by an off-diagonal mass shift ΔM_SD that is neither decomposed channel-by-channel nor tested for sensitivity to the model choices that define it. That is a real gap between claim and evidence.

What is genuinely useful: the calculation uses standard machinery—GI potential for bare masses, QPC for couplings, once-subtracted dispersion relations for self-energies—and the parameters are mostly external. The mixing angle is not fitted; it emerges from the equations, which is methodologically honest. The prediction of a partner state at 4387 MeV with width 44 MeV and a distinctive DD*_2(2460) decay is concrete and testable at BESIII.

Now the soft spot, and it is exactly where the stress-test note points. A 35° angle with diagonal masses of 4279 and 4339 MeV requires an off-diagonal element of roughly 70–80 MeV, comparable to the large diagonal shifts. That element comes from Eq. (17) with the subtraction point fixed at M_J/ψ and a truncated loop set. The paper reports no decomposition of ΔM_SD by channel, so we cannot see whether it is dominated by the D*D1 loops whose thresholds sit near the bare masses, and no variation of the subtraction point, channel set, or γ. The abstract's claim that DD1 controls the lower state while D*D1 affects the higher state is not supported by Table II: D*D1(2430)0 is the largest single D1 contribution to both diagonal shifts. Both channels clearly matter, and the attribution in the abstract oversimplifies.

None of this makes the paper unserious. The framework is coherent, the literature is handled honestly (including a note added acknowledging the parallel bottomonium work), and the derived mass and width of the lower state are in the right ballpark. The issue is that the central result—the large mixing angle—rests on model choices that could easily dominate the outcome, and the paper calls it robust without a robustness check.

This paper is for hadron spectroscopists working on charmonium above 4 GeV. It deserves a serious referee, but the referee should insist on a sensitivity study: channel-by-channel ΔM_SD, variation of the subtraction point, and variation of γ. Without that, the paper presents a plausible model, not yet a demonstration.

Referee Report

4 major / 4 minor

Summary. This paper proposes that the ψ(4220) and ψ(4380) are the lower and upper members of a 4S–3D charmonium mixing pair, with the mixing induced by coupled-channel (hadronic loop) effects rather than by the tensor force. The authors compute bare 4S and 3D masses with the Godfrey–Isgur potential, then evaluate diagonal and off-diagonal mass shifts from a truncated set of open-charm loops using a once-subtracted dispersion relation with the J/ψ mass as subtraction point and QPC vertices. The resulting lower mixed state has mass 4235.9 MeV, width 32.0 MeV, and mixing angle θ=35°, while the upper state has mass 4387.1 MeV and width 44.1 MeV. The paper argues that DD1 loops dominate the lower state and D*D1 loops the upper state, and it gives partial-width predictions for experimental searches.

Significance. If the central prediction is robust, the paper would provide a dynamical explanation for the unusually large 4S–3D mixing angle that previous phenomenological studies had to insert by hand, and it would sharpen the case for a narrow partner state near 4387 MeV. The framework is coherent and the paper is explicit about its parameter choices (γ=0.44 fixed from P-wave charmonia, GI parameters from Ref. [63], φ from the heavy-quark limit). The main weakness is that the decisive off-diagonal mass shift, which generates θ=35°, is computed with a truncated loop set and a fixed subtraction point, and the paper offers no sensitivity analysis. Because the mixing angle is the load-bearing quantity for the ψ(4220) identification, the result is not yet robust as a dynamical prediction.

major comments (4)

[Section III.B, Eqs. (13)–(17)] The central quantity ΔM_SD(M) is not characterized enough to support the claimed θ=35°. With the diagonal masses 4279.1 and 4338.5 MeV, the relation tan(2θ)=2V/(M_D−M_S) requires an off-diagonal element V≈80 MeV, yet the paper reports no channel-by-channel decomposition of ΔM_SD and no variation of the subtraction point, the channel set, or γ=0.44. Eq. (17) vanishes identically at the subtraction point M_J/ψ, so any omitted high-mass channel contributes a subtraction-point-dependent constant at M≈4.2 GeV; the statement that the once-subtracted dispersion relation 'effectively limits the number of loops' is not a substitute for a convergence test. Without such a test, the large mixing angle could be an artifact of the truncation and subtraction scheme rather than a robust prediction.
[Section III.B and Table III] The claimed agreement with experiment is not quantified. The predicted mass 4235.9 MeV differs from the PDG value 4222.1±2.3 MeV by 13.8 MeV, which is about six times the quoted experimental uncertainty, and the computed two-body OZI width 32.0 MeV differs from the measured total width 49±7 MeV by about 17 MeV. The comparison is presented without any estimate of theoretical uncertainty from the GI parameters, the QPC strength, the truncation of the loop sum, or neglected decay modes. Before identifying ψ(4220) with ψ'_{4S-3D}, an error budget is needed.
[Section III.A and Table II] Several of the largest loop contributions, especially D*D1(2430)0 with ΔM_i=−22.8 MeV for ψ(4S) and −21.4 MeV for ψ(3D), are computed by treating D1(2430) as a stable two-body state even though its physical width is about 314 MeV. The threshold position for such a broad state is ill-defined, and convolution over the D1 spectral function could significantly change both the diagonal shifts and the off-diagonal ΔM_SD. The same channels also depend on the D1 mixing angle φ, which is fixed to −54.7° from the heavy-quark limit without any variation. Given that these channels are among the dominant entries, the sensitivity of the mixing angle to this approximation should be tested.
[Eqs. (16)–(17) and Table II] The paper's criterion for including loops is that thresholds lie below the bare masses, but this is a model cut that is not justified by the once-subtracted dispersion relation. The subtracted expression still receives contributions from all channels at M≈4.2 GeV, and there is no argument that channels with thresholds just above 4.5 GeV are negligible. A practical test would be to vary the cutoff by adding the next few expected channels (for example, D*D*0(2550) or D_sD_s1(2536)) and to vary the subtraction point within a reasonable range; the stability of θ under these changes should be reported.

minor comments (4)

[Section III.A, Fig. 2] The text says ψ(4S) has three discontinuities corresponding to DD1(2430)0, D*D1(2430)0, and D*D0(2300), but Table II also assigns nonzero mass shifts to DD1(2420) and D*D1(2420) for ψ(4S); please clarify which channels produce cusps and why D*D1(2420) is not counted.
[Abstract and Section IV] The statement that the results 'align with experimental observations' is stronger than the numbers warrant, given the 13.8 MeV mass gap and 17 MeV width gap; a softer wording with an explicit caveat about model uncertainty would be more accurate.
[Eq. (19)] The quantity Γtotal is defined as a sum of two-body OZI-allowed widths, but it is later compared with the measured total width of ψ(4220); the paper should state explicitly that contributions from three-body channels, radiative decays, and hidden-charm final states are assumed negligible.
[Notation in Eqs. (11)–(13)] The symbols M0_S and M0_D are used for bare masses while M denotes the physical mass in Eq. (24); this notational clash is confusing and should be clarified, for example by using M_bare and M_phys.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the coupled-channel θ=35° mixing angle is a computed output, not a fit to the ψ(4220) mass or width; self-citations are contextual rather than load-bearing.

full rationale

The central numerical claims—θ=35°, m=4235.9 MeV, and Γ=32.0 MeV—are outputs of diagonalizing the 2×2 effective Hamiltonian in Eq. (13), with diagonal mass shifts from Eq. (16) and the off-diagonal loop integral ΔM_SD from Eq. (17). The inputs are fixed before the ψ(4220) data enter: the GI-model bare masses (4433.0 and 4491.3 MeV), the QPC strength γ=0.44 chosen to reproduce the ψ(3915) and ψ(3930) widths, and the once-subtracted loop scheme with subtraction point M_J/ψ. The PDG values for ψ(4220) appear only as post-hoc comparisons in Tables II and III, not as fitted constraints. No equation or parameter is defined in terms of the target mass, width, or mixing angle, so no self-definitional or fitted-input-called-prediction circularity is present. The paper does rely on Ref. [50], a same-group paper, for the prior 4S–3D mixing scheme and on Ref. [63] for GI-model parameters, but these are not load-bearing in a circular sense: Ref. [50] supplies a phenomenological fit that this work independently tests with a distinct dynamical mechanism, and Ref. [63] supplies parameters fitted to low-lying charmonia rather than to ψ(4220). The real weakness—sensitivity of ΔM_SD to the truncated hadronic-loop set and to the subtraction point—is a robustness or correctness concern, not circularity, because the paper does not tune these choices to reproduce the quoted output. Accordingly, no circular reduction can be exhibited; the score reflects only the presence of minor self-citations that are not load-bearing.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central result depends on the GI bare masses, the QPC coupling strength, the chosen channel set, and the subtraction scheme. None of these are fitted to the ψ(4220) or ψ(4380) parameters, so the calculation is not circular, but the many modeling choices mean the numerical output should be interpreted with caution.

free parameters (4)

QPC pair creation strength γ = 0.44
Fitted to reproduce the decay widths of ψ(3915) and ψ(3930); controls all charmonium to meson-meson coupling vertices and thus the magnitude of ΔM and ΔM_SD.
GI model parameters (Table I) = see Table I (from Ref [63])
Fitted to low-lying charmonium masses in prior work; these determine the bare masses and wave functions used as inputs.
Subtraction point M_J/ψ = 3.0969 GeV (J/ψ mass)
Chosen by hand as the subtraction point in the once-subtracted dispersion relation; the dependence of results on this choice is not studied.
D1 mixing angle φ = -54.7 degrees
Taken from the heavy quark limit to define D1(2420) and D1(2430); not varied or fitted.

assumptions (4)

domain assumption The once-subtracted dispersion relation with subtraction point at M_J/ψ renders the mass-shift integrals finite and physical.
Invoked from Ref [82] without derivation; the choice of subtraction point is a modeling assumption that could affect the large off-diagonal mixing.
domain assumption The quark pair creation model describes the coupling between charmonium and meson-meson channels.
The vertex strength and spin/flavor/color factors are assumed to follow the QPC Hamiltonian of Refs [85,86].
ad hoc to paper Only hadronic loops with thresholds below the bare masses are included, and interactions between the final mesons B and C are neglected.
Sec. III states this choice without justification; omitting B-C interactions is a significant simplification that could alter the mass shifts.
domain assumption HQSS is used to treat D/D* and D1(2420)/D1(2430) as mixed states with fixed mixing angle.
Footnote and Eq. (22) rely on the heavy quark limit; corrections are not estimated.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Coupled-channel study of $4S$-$3D$ mixing dynamics in $\psi(4220)$ and $\psi(4380)$." pith.science (2026). https://pith.science/paper/2WNP2MXF

@misc{pith2026250208072,
  author       = {Pith},
  title        = {Pith review of: Coupled-channel study of $4S$-$3D$ mixing dynamics in $\psi(4220)$ and $\psi(4380)$},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2WNP2MXF}},
  note         = {Machine review of arXiv:2502.08072}
}

abstract

Among charmoniumlike $XYZ$ states, the $\psi(4220)$ and $\psi(4380)$ states have emerged as key candidates for exploring the charmonium spectrum. In this work, we propose a $4S$-$3D$ charmonium mixing scheme for the $\psi(4220)$ and $\psi(4380)$, induced by coupled-channel effects. By constructing a coupled-channel model, we identify the dynamical mechanism responsible for the large mixing angle observed in previous studies, which cannot be explained by conventional potential models alone. Our analysis reveals that the $DD_1$ channel significantly influences the lower state ($\psi(4220)$), while the $D^*D_1$ channel primarily affects the higher state ($\psi(4380)$). Furthermore, we investigate the two-body Okubo-Zweig-Iizuka (OZI)-allowed strong decay behaviors of these states, providing insights into their total widths. This study not only supports the $4S$-$3D$ mixing scheme but also offers a deeper understanding of the role of coupled channels in shaping the charmonium spectrum above 4 GeV. Our results align with experimental observations and provide a framework for interpreting future data on charmonium states.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

Discovery potential of charmonium $2P$ states through the $e^+e^- \to \gamma D\bar{D}$ processes
hep-ph 2024-12 conditional novelty 5.0 of 10

Using the measured e+e- -> gamma X(3872) cross section as a reference, the authors predict measurable e+e- -> gamma chi_c0(2P) and e+e- -> gamma chi_c2(2P) production and propose the gamma D Dbar final state as a sear...

Pith tools

Reviewed August 8, 2026 · model on record in the stance chip above.