{"id":"84f3d6c6-7f53-4586-89d1-5de345d1556f","arxiv_id":"2502.08072","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Coupled-channel effects can induce a 35 degree 4S-3D mixing angle in charmonium, reproducing the ψ(4220) properties and predicting a ψ(4380) partner with distinctive decays.","lead":"This paper uses a coupled-channel model to show that interactions with open-charm meson pairs can generate a large mixing between the 4S and 3D charmonium configurations, explaining the mass and width of the ψ(4220) and predicting a partner state around 4380 MeV. It matters because it offers a dynamical mechanism, rather than a simple parametrization, for why these charmonium states mix, and it suggests decay channels where the predicted partner could be found.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 35° mixing angle rests on an unvalidated off-diagonal mass shift: the once-subtracted dispersion relation and truncated loop set in Eqs. (16)–(17) are not tested for sensitivity, so the central coupled-channel mixing claim is not yet robust.","rationale":"The paper is a competent application of a standard coupled-channel framework. The internal logic is coherent: bare GI masses, QPC vertices, once-subtracted self-energies, and a 2×2 diagonalization. The resulting ψ(4220) mass and width are in the right ballpark, and the partner at 4387 MeV is a falsifiable prediction. However, the central new result—that the coupled channels produce θ=35°—is controlled entirely by the off-diagonal ΔM_SD. That quantity is not shown channel-by-channel, and the subtraction scheme plus truncated channel list are exactly the places where model dependence enters. The required off-diagonal is ~70 MeV, so even a 20–30% change in the dominant loop contributions, or a different subtraction point, could move θ substantially. The paper presents no such checks. This is not an internal inconsistency, but it is a missing validation step for the headline conclusion. The reader's conditional verdict is appropriate; no stronger action is warranted.","tokens_in":19363,"tokens_out":13988,"duration_ms":117086,"concrete_test":"Recompute the effective 2×2 mass matrix and the resulting mixing angle θ with the same inputs but set the subtraction point in Eqs. (16)–(17) to M_ψ(2S)=3.686 GeV and then to M_ψ(4260)=4.230 GeV; if |θ| varies by more than 10° across this range, the large-mixing conclusion is an artifact of the subtraction scheme rather than a robust prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical result, θ=35°, is produced by the off-diagonal mass shift ΔM_SD in Eqs. (16)–(17), evaluated with the J/ψ mass as subtraction point and a truncated set of hadronic loops. The subtraction in Eq. (16) is an algebraic identity at fixed amplitude, so the physical result depends on the subtraction point whenever the loop sum is truncated; high-mass channels omitted from Table II contribute an unknown, M_J/ψ-dependent constant. For θ=35° with diagonal masses of 4279 and 4339 MeV, the required off-diagonal element is about 70 MeV, larger than the 0.5° tensor contribution and comparable to the individual diagonal shifts. The paper reports neither the channel-by-channel decomposition of ΔM_SD nor any variation of the subtraction point, channel set, or QPC strength γ=0.44. In particular, the D*D1(2430)0 and D*D1(2420) loops, whose thresholds lie near or above the bare masses, are among the largest entries in Table II and are the most sensitive to the D1 mixing angle and to the broad D1(2430) width; their off-diagonal contributions could dominate or cancel. Without a sensitivity study, the large mixing angle is not established as a robust dynamical prediction, so the identification of ψ(4220) as the 4S–3D mixed state rests on an unquantified model choice.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":19702,"tokens_out":6633,"duration_ms":58810,"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":[{"comment":"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":"Section III.B, Eqs. (13)–(17)"},{"comment":"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":"Section III.B and Table III"},{"comment":"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.","section":"Section III.A and Table II"},{"comment":"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.","section":"Eqs. (16)–(17) and Table II"}],"minor_comments":[{"comment":"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.","section":"Section III.A, Fig. 2"},{"comment":"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.","section":"Abstract and Section IV"},{"comment":"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.","section":"Eq. (19)"},{"comment":"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.","section":"Notation in Eqs. (11)–(13)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope and builds on a recognizable line of work by the same group. The main concern is not the overall scheme but the lack of any sensitivity analysis for the off-diagonal mass shift that produces the headline mixing angle; this is fixable within the manuscript's scope and should be the primary request for revision. I would also encourage the editor to ask the authors to state clearly which of their results are robust under reasonable variations of the truncation and subtraction choices."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nRead 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.\n\nWhat 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.\n\nNow 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.\n\nNone 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.\n\nThis 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.","headline":"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.","tokens_in":20271,"tokens_out":4984,"would_cite":true,"duration_ms":39067,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["charmonium","ψ(4220)","ψ(4380)","4S-3D mixing","coupled-channel model","hadronic loops","OZI-allowed strong decays","XYZ states"],"falsifier":"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.","tokens_in":19145,"feed_emoji":"⚛️","tokens_out":9824,"duration_ms":79591,"temperature":0.7,"pith_summary":"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.","feed_headline":"35° mix places ψ(4220) at 4235.9 MeV, predicts ψ(4380)","feed_subtitle":"Coupled-channel loops, not the tensor force, drive the large mixing angle and match the observed 32 MeV width.","key_machinery":"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).","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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°."],"supporting_citations":[{"why":"Introduces the 4S-3D mixing scheme for ψ(4220) and the phenomenological 30°–36° mixing angle that this paper aims to explain dynamically.","marker":"[50]"},{"why":"Supplies the once-subtracted dispersion relation that tames the infinite hadronic-loop sums and underlies the off-diagonal mass shift ΔM_SD.","marker":"[82]"},{"why":"Supplies the potential model defining the bare ψ(4S) and ψ(3D) masses and the tensor-force baseline of 0.5° that the coupled-channel result must exceed.","marker":"[57]"},{"why":"Fixes the potential-model parameters and basis size used for all wave functions, masses, and decay amplitudes.","marker":"[63]"},{"why":"Provides the quark-pair-creation model whose transition amplitudes give the meson-loop couplings and OZI-allowed decay widths.","marker":"[85, 86]"},{"why":"Supplies the experimental ψ(4220) mass and width, 4222.1±2.3 MeV and 49±7 MeV, against which the predictions are compared.","marker":"[87]"},{"why":"Establishes the coupled-channel-induced S-D mixing framework for charmonia that the present calculation applies to the 4S-3D pair.","marker":"[79]"},{"why":"Offers the screened-potential prediction of the ψ(4S) mass that agrees with the diagonal coupled-channel result and anchors the comparison.","marker":"[42]"}],"fun_headline_variants":["35° mix from loops: ψ(4220) width 32 MeV, ψ(4380) predicted","Hadronic loops, not tensor force, set 35° charmonium mixing","Coupled channels yield 35° mixing for ψ(4220) and ψ(4380)","Loops lower charmonia masses, 35° mix explains ψ(4220) and ψ(4380)","Predicting ψ(4380): 35° mixing via coupled-channel loops"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["35° mix from loops: ψ(4220) width 32 MeV, ψ(4380) predicted","Hadronic loops, not tensor force, set 35° charmonium mixing","Coupled channels yield 35° mixing for ψ(4220) and ψ(4380)","Loops lower charmonia masses, 35° mix explains ψ(4220) and ψ(4380)","Predicting ψ(4380): 35° mixing via coupled-channel loops"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000839,"raw_usage":{"total_tokens":3672,"prompt_tokens":972,"completion_tokens":2700,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":588,"completion_tokens_details":{"reasoning_tokens":2579}},"tokens_in":588,"tokens_out":2700,"duration_ms":17333,"temperature":1.0,"reasoning_tokens":2579,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T10:55:54.115267+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}