{"id":"0fdca96b-285a-4bf4-8856-6e9bef076f85","arxiv_id":"2412.06400","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"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 search channel.","lead":"This paper predicts that the charmonium states chi_c0(2P) and chi_c2(2P) can be produced in electron-positron collisions through a radiative decay chain that leaves a photon and a D Dbar pair. It uses measured gamma X(3872) data as a calibrator and recommends that BESIII and Belle II search for these states in the gamma D Dbar invariant mass spectrum.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Predicted measurability of chi_c2(2P) hinges on ad hoc g_psiD1D = gmax/4; Table II shows the ratio Gamma(gamma chi_c2)/Gamma(gamma chi_c1) drops by an order of magnitude at gmax or gmax/2, so the quoted 0.47-0.64 range is not robust.","rationale":"The reader's weakest_assumption and my independent pass converge on g_psiD1D. I checked the normalization chain: Eq. (1) is a standard Breit-Wigner cross section, Eq. (29) uses BR(X(3872)->pi+pi-J/psi) = 3.5% and assumes psi(4230) dominance, and none of these is the dominant fragility because an overall normalization error shifts all three chi_cJ cross sections together and does not destroy the angular separation test. By contrast, Gamma(gamma chi_c2)/Gamma(gamma chi_c1) is the only place where an unconstrained parameter multiplies the central prediction by an order of magnitude, and it directly controls the measurability of chi_c2(2P). Fig. 2(a) shows only mild alpha dependence, but Table II exposes the much larger g_psiD1D dependence. Because the abstract asserts that e+e- -> gamma D Dbar is an ideal process without the gmax/4 caveat, the quantitative claim outstrips its input constraint. This is fully consistent with the reader's conditional verdict, so I recommend no change to that verdict.","tokens_in":19046,"tokens_out":4798,"duration_ms":50248,"concrete_test":"Recompute the hadronic-loop ratios of Sec. II with g_psiD1D fixed by the coupled-channel psi(4220)/psi(4380) solution of Ref. [34] (or by a direct determination of psi(4230) -> D1(2420)Dbar), for alpha = 3, 4, and 5, and compare with Table II. If the resulting Gamma(gamma chi_c2):Gamma(gamma chi_c1) is at or below the gmax/2 endpoint of about 0.17, then the Fig. 2(b) sigma[gamma chi_c2(2P)] prediction is at least a factor of 3-10 too large, and the 'ideal process' conclusion should be downgraded to an upper bound pending a precise g_psiD1D.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is that sigma[e+e- -> gamma chi_c2(2P)] is large enough to make e+e- -> gamma D Dbar an ideal discovery channel. That claim is controlled by the computed ratio Gamma(gamma chi_c2)/Gamma(gamma chi_c1), which is dominated by the D1(2420)Dbar loop and therefore by the coupling g_psiD1D. In Sec. III the authors derive gmax = 10.19 GeV^-1/2 from the measured R[psi(4230)] using a D1(2420)-mediated model of e+e- -> pi+ D*− D0, then arbitrarily set g_psiD1D = gmax/4 to absorb contributions from other intermediate channels. Table II shows the consequence: with gmax the ratio Gamma(gamma chi_c2):Gamma(gamma chi_c1) is 0.029-0.042, with gmax/2 it is 0.118-0.168, and the quoted 0.47-0.64 is reached only for the unshown gmax/4 case. Since sigma[gamma chi_c2(2P)] is obtained by multiplying the BESIII-extracted sigma[gamma chi_c1(2P)] by this ratio, the predicted chi_c2 signal and hence the 'ideal process' claim vary by more than an order of magnitude under a plausible change of a single unmeasured coupling. No uncertainty band for g_psiD1D is propagated into the cross-section figures. The authors note that a precise determination awaits additional measurements, but the abstract and conclusion state the discovery potential without that caveat.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes e+e− → γ D Dbar at √s = 4.23 GeV as a discovery channel for the charmonium 2P states χc0(2P) and χc2(2P). The production is assumed to proceed through ψ(4230) → γχcJ(2P), with the radiative transition computed in a hadronic-loop model that includes D(∗)D(∗) and D1(2420)Dbar loops plus gauge-invariant contact terms. Using the measured BESIII cross section for e+e− → γX(3872) as a reference and identifying X(3872) with χc1(2P), the authors derive cross sections for e+e− → γχc0(2P) and e+e− → γχc2(2P), then predict the D Dbar invariant mass spectrum of e+e− → γD Dbar and the photon angular distribution for the two states. The central quantitative outputs are the partial-width ratios in Eq. (28), the scaled cross sections in Fig. 2(b), the invariant-mass spectra in Fig. 3, and the angular parameter β2 in Table III.","tokens_in":19481,"tokens_out":4277,"duration_ms":44163,"significance":"If the predictions are robust, the paper gives a concrete, testable production mechanism for two poorly established charmonium states and provides phenomenological guidance for BESIII and Belle II. The hadronic-loop amplitudes are presented in enough detail to be reproduced, and the contact terms are fixed by the Ward-Takahashi identity rather than fitted, which is a genuine strength. The predicted β0 = 1 for ψ → γχc0(2P) is essentially model-independent, and the differential spectra and angular distributions are falsifiable. However, the headline conclusion that e+e− → γD Dbar is an ideal discovery channel is not yet supported with quantified uncertainties, because the predicted χc2(2P) signal is controlled by an unconstrained coupling choice.","major_comments":[{"comment":"The stress-test concern is borne out. The central claim that σ[γχc2(2P)] is measurable rests on the ratio Γγχc2/Γγχc1, and the quoted range 0.47–0.64 in Eq. (28) is obtained only for gψD1D = gmax/4. Table II shows that at gψD1D = gmax/2 the ratio is 0.118–0.168 and at gψD1D = gmax it is 0.029–0.042, i.e. a drop by more than an order of magnitude at the upper end of the plausible interval. Since σ[γχc2(2P)] is obtained in Fig. 2(b) by multiplying the BESIII-extracted σ[γχc1(2P)] by exactly this ratio (Eq. (2)), the predicted χc2 signal—and therefore the \"ideal process\" conclusion in the abstract—is not robust to the choice of gψD1D. The authors acknowledge the dependence in the text, but no uncertainty band is propagated into Fig. 2(b), Fig. 3, or the conclusions. This is load-bearing and should be fixed by either constraining gψD1D with a data-driven error estimate or by presenting the cross-section predictions as functions of gψD1D and softening the abstract and conclusions accordingly.","section":"§III, Eq. (28) and Table II"},{"comment":"The derivation of gmax assumes that the measured quantity Rψ(4230) = 2.70 eV is saturated by a single D1(2420)Dbar loop, and the subsequent reduction gψD1D = gmax/4 is introduced as an order-of-magnitude guess to absorb other intermediate channels. The factor 1/4 is not obtained from any fit or independent constraint, and it is the single parameter that most strongly controls the final cross sections. The manuscript should either derive a range for this factor from the coupled-channel results it cites, or explicitly present the central predictions as conditional on gψD1D = gmax/4 rather than as unconditional discovery potentials.","section":"§III, Eq. (24)"},{"comment":"The absolute normalization of σ[γχc0(2P)] and σ[γχc2(2P)] is obtained by dividing the measured e+e− → γX(3872) → γπ+π−J/ψ cross section by BR(X(3872) → π+π−J/ψ) = 3.5% and identifying X(3872) with χc1(2P). This introduces an external systematic uncertainty that is not propagated into the figures, and the comparison with data is made \"ignoring background contributions\" without reporting a fit quality or an uncertainty on the extracted σ[γχc1(2P)]. The absolute scale of the predictions therefore carries an additional unquantified systematic error beyond the α and gψD1D variations. The authors should quote this uncertainty or explicitly state in the conclusions that the absolute cross sections are uncertain at this level.","section":"§III, Eq. (29) and Fig. 2(b)"}],"minor_comments":[{"comment":"The notation \"dN/Nd cosθ\" should be written as dN/(N d cosθ) for clarity.","section":"§III, Eq. (30)"},{"comment":"The caption says \"The data points are obtained by dividing the original BESIII experiment data [23] by the branching fraction of X(3872)→π+π−J/ψ,\" but the curve labeled σ[γX(3872)] appears to be the reference cross section before this division; the relation between the displayed curves and the data should be stated explicitly.","section":"§III, Fig. 2(b)"},{"comment":"The abstract and concluding paragraph state that e+e− → γD Dbar is an ideal process without repeating the caveat, stated in §III, that a precise determination of gψD1D requires additional experimental measurements; a one-sentence qualification in the conclusions would make the paper internally consistent.","section":"§IV"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The paper gives concrete, new predictions for e+e- -> gamma chi_c0/2(2P) at sqrt(s)=4.23 GeV and a gamma D Dbar mass spectrum, using psi(4230) as the intermediate state. The genuinely useful piece is the angular distribution: beta0=1 for chi_c0 and beta2~0.11 for chi_c2, nearly independent of the cutoff alpha. That is a real handle for separating two states whose masses differ by only 0.4 MeV.\n\nThe soft spot, and it is a big one, is the normalization. The cross sections are obtained by rescaling the BESIII e+e- -> gamma X(3872) data with computed width ratios. The ratio Gamma(gamma chi_c2)/Gamma(gamma chi_c1) is quoted as 0.47-0.64 in Eq. (28), but Table II shows that with g_psiD1D = gmax (the value derived from R_psi(4230)) that ratio is 0.029-0.042, and with gmax/2 it is 0.118-0.168. The 0.47-0.64 range appears to come from the unshown gmax/4 case, which the authors adopt ad hoc to absorb other channels. So the predicted chi_c2 signal, and hence the claim that e+e- -> gamma D Dbar is 'ideal' for identifying it, changes by more than an order of magnitude under a plausible variation of a single coupling. The authors do flag that g_psiD1D is not well determined, but the abstract and conclusion do not carry that caveat.\n\nOther, lesser issues: the gamma D Dbar spectrum is plotted without any continuum background, and at 4.23 GeV the open-charm background is not negligible. And because the input sigma[gamma chi_c1(2P)] is itself extracted from the measured e+e- -> gamma X(3872) cross section, this is not an independent prediction of the production rate; it is a model-dependent ratio times an experimental number. The small mass gap means the D Dbar spectrum alone will not separate the two states — the angular distribution is the better discriminator, and that part survives.\n\nCredit where due: the hadronic loop calculation is careful, with gauge-invariant contact terms and the D1(2420)bar-D channel included. Table II is honest — it shows what happens when you vary g_psiD1D. The beta2 result being alpha-insensitive is a nice, robust output.\n\nWho is this for: hadron-spectroscopy phenomenologists and BESIII/Belle II experimentalists. It deserves a serious referee: the process proposal and angular observable are worth publishing, but the cross-section claims need a proper uncertainty treatment for g_psiD1D and a background estimate. Send it to peer review, and ask for those two things.","headline":"A testable proposal with one clean observable (the chi_c0 vs chi_c2 angular distribution), but the preprint's central cross-section claim swings by an order of magnitude with the ad hoc coupling g_psiD1D.","tokens_in":20062,"tokens_out":2909,"would_cite":true,"duration_ms":28851,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper predicts that $e^+e^-\\to \\gamma D\\bar D$ at $\\sqrt{s}=4.23$ GeV is a discovery channel for both $\\chi_{c0}(2P)$ and $\\chi_{c2}(2P)$, with a measurable $\\chi_{c2}$ signal and a photon angular distribution that separates the two…","keywords":["charmonium 2P states","e+e− annihilation","γ D\\bar D production","hadronic loop mechanism","χc0(2P)","χc2(2P)","ψ(4230)","radiative charmonium transitions"],"falsifier":"Measure the $e^+e^-\\to\\gamma D\\bar D$ cross section near $\\sqrt{s}=4.23$ GeV and look for a $D\\bar D$ invariant-mass enhancement around 3.922 GeV whose photon polar-angle distribution in the $e^+e^-$ center-of-mass frame matches $1+\\cos^2\\theta$ (for $\\chi_{c0}$) or the nearly flat shape with $\\beta_2\\approx 0.11$ (for $\\chi_{c2}$); absence of such an enhancement would falsify the prediction. A separate decisive check is a precise measurement of $\\psi(4230)\\to\\pi^+ D^{*-}D^0$, which fixes $g_{\\psi D_1 D}$ and determines whether the claimed cross sections survive.","tokens_in":18789,"feed_emoji":"🔭","tokens_out":12574,"duration_ms":102377,"temperature":0.7,"pith_summary":"This paper argues that the radiative process $e^+e^-\\to \\gamma D\\bar D$ at $\\sqrt{s}=4.23$ GeV should be a discovery channel for the two missing charmonium $2P$ states, $\\chi_{c0}(2P)$ and $\\chi_{c2}(2P)$. Its strategy is to scale from the measured $e^+e^-\\to \\gamma X(3872)$ cross section, treating $X(3872)$ as the $\\chi_{c1}(2P)$ partner, and to fix the relative production rates with a hadronic-loop calculation of $\\psi(4230)\\to \\gamma \\chi_{cJ}(2P)$. The predicted $e^+e^-\\to \\gamma \\chi_{c2}(2P)$ cross section is comparable in size to the measured $\\gamma X(3872)$ reference and therefore within reach of BESIII and Belle II. Because both $2P$ states decay predominantly to $D\\bar D$, their production would show up as a bump in the $D\\bar D$ invariant mass spectrum, with a photon angular distribution that separates the two states even though their masses differ by less than 1 MeV.","feed_headline":"One radiative process should expose both missing 2P charmonia","feed_subtitle":"Hadronic-loop ratios scale from measured γX(3872) data to a visible χc2(2P) signal, separable by photon angle.","key_machinery":"The central object is the hadronic-loop amplitude for $\\psi(4230)\\to \\gamma \\chi_{cJ}(2P)$: triangle diagrams with virtual $D$, $D^*$, and $D_1(2420)$ meson loops plus contact diagrams whose strength is fixed by the Ward-Takahashi identity. A dipole form factor $F(q^2)=\\left((m_E^2-\\Lambda^2)/(q^2-\\Lambda^2)\\right)^2$ with $\\Lambda=m_E+\\alpha\\Lambda_{\\rm QCD}$ regularises the loops, and the couplings of the $\\chi_{cJ}(2P)$ multiplet to charmed mesons share a universal constant $g_p$ that cancels in the width ratios. The result is a prediction that depends only on the cutoff parameter $\\alpha\\in[3,5]$ and on the $\\psi(4230)\\to D_1\\bar D$ coupling $g_{\\psi D_1 D}$.","core_discovery":"Treating $X(3872)$ as $\\chi_{c1}(2P)$, the BESIII data on $e^+e^-\\to \\gamma X(3872)$ fix the cross section for $e^+e^-\\to \\psi(4230)\\to \\gamma \\chi_{c1}(2P)$, and the paper's hadronic-loop calculation yields the quoted width ratios $\\Gamma_{\\gamma\\chi_{c0}}:\\Gamma_{\\gamma\\chi_{c1}}:\\Gamma_{\\gamma\\chi_{c2}} = (0.07\\text{--}0.13):1:(0.47\\text{--}0.64)$ for cutoff $\\alpha=3\\text{--}5$. From these ratios the paper predicts that $\\sigma[e^+e^-\\to \\gamma \\chi_{c2}(2P)]$ is slightly below but of the same order as $\\sigma[\\gamma \\chi_{c1}(2P)]$, while $\\sigma[\\gamma \\chi_{c0}(2P)]$ is smaller but still favourably detectable; both are large enough for the next round of BESIII and Belle II data. With $\\chi_{c0}(2P)\\to D\\bar D$ treated as 100\\% and $\\chi_{c2}(2P)\\to D\\bar D$ as about 60\\%, the process $e^+e^-\\to \\gamma D\\bar D$ should show two overlapping 3.9 GeV structures in the $D\\bar D$ mass distribution, and the angular distribution of the photon ($\\beta_0=1$ for $\\chi_{c0}$ versus $\\beta_2\\approx 0.11$ for $\\chi_{c2}$) offers a kinematic handle to tell them apart.","pith_inferences":["If the prediction holds, the same scaling method could be applied to other hidden-charm states whose $\\gamma X$ production through $\\psi(4230)$ is kinematically open and whose branching ratios are known.","The near-flat photon distribution predicted for $\\chi_{c2}$ means a modest number of events may distinguish the two states by fitting the photon polar angle alone, without resolving their sub-MeV mass gap.","A precise measurement of $\\psi(4230)\\to \\pi^+ D^{*-}D^0$ would pin down the one coupling that currently dominates the uncertainty and would turn the prediction from a range into a sharp number."],"forward_implications":["The $e^+e^-\\to \\gamma D\\bar D$ channel would let $\\chi_{c0}(2P)$ and $\\chi_{c2}(2P)$ be observed simultaneously in one production process at BESIII or Belle II.","The $D\\bar D$ invariant mass spectrum should contain a two-peak structure near 3.9 GeV whose relative size is set by the computed width ratios.","The photon angular distribution distinguishes $\\chi_{c0}(2P)$ and $\\chi_{c2}(2P)$ even though their masses are almost equal.","A $\\chi_{c2}(2P)$ signal at the predicted level would corroborate the treatment of $\\psi(4230)$ as an unquenched charmonium state and support the hadronic-loop description of its radiative decays."],"supporting_citations":[{"why":"supplies the measured $e^+e^-\\to\\gamma X(3872)$ cross-section data that are scaled to $\\sigma[\\gamma\\chi_{c1}(2P)]$","marker":"[23]"},{"why":"first observation of $e^+e^-\\to\\gamma X(3872)$ and of the resonance-like structure near 4.2 GeV that motivates $\\psi(4230)$ as the intermediate state","marker":"[22]"},{"why":"provides the $\\psi(4230)$ couplings $g_{\\psi DD}$, $g_{\\psi D^*D}$, $g_{\\psi D^*D^*}$ and the dilepton width used as input","marker":"[8]"},{"why":"LHCb amplitude analysis that established $\\chi_{c0}(2P)$ and $\\chi_{c2}(2P)$ in $B^+\\to D^+D^-K^+$ and gives their masses","marker":"[16]"},{"why":"model-independent confirmation of the resonant structure near 3.9 GeV in $B^+\\to D^+D^-K^+$ that fixes the two $2P$ states","marker":"[17]"},{"why":"establishes the $D\\bar D$ decay pattern and the two-substructure picture of the 3.9 GeV enhancement used for the branching ratios","marker":"[14]"},{"why":"provides the hadronic-loop technique and the gauge-invariance-fixed contact terms used in the width calculation","marker":"[40]"},{"why":"supplies the effective Lagrangians and Feynman rules for $\\psi\\to\\gamma\\chi_{cJ}(2P)$ used in the loop amplitudes","marker":"[42]"},{"why":"gives the measured $\\Gamma_{e^+e^-}^{\\psi(4230)}\\cdot BR(\\psi(4230)\\to\\pi^+D^{*-}D^0)$ used to estimate the $g_{\\psi D_1 D}$ coupling","marker":"[46]"}],"fun_headline_variants":["Two missing charmonia could show up in one photon channel","Radiative process may reveal both 2P charmonia at BESIII","Photon-tagged D mesons could unmask χc0 and χc2","Angular cut separates χc0 and χc2 in γ D Dbar events"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The size of the predicted cross sections rests on the assumed value $g_{\\psi D_1 D}=g_{\\rm max}/4\\approx 2.55$ GeV$^{-1/2}$ for the $\\psi(4230)$ coupling to the $D_1(2420)\\bar D$ channel, chosen because $g_{\\rm max}$ overestimates the coupling; if the true coupling were closer to $g_{\\rm max}$, the ratio $\\Gamma_{\\gamma\\chi_{c2}}:\\Gamma_{\\gamma\\chi_{c1}}$ would shrink from about 0.5 to below 0.05 and the measurable signal would largely disappear.","fun_headline_variants_meta":{"raw":{"variants":["Two missing charmonia could show up in one photon channel","Radiative process may reveal both 2P charmonia at BESIII","Photon-tagged D mesons could unmask χc0 and χc2","Angular cut separates χc0 and χc2 in γ D Dbar events"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000245,"raw_usage":{"total_tokens":1657,"prompt_tokens":1186,"completion_tokens":471,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":802,"completion_tokens_details":{"reasoning_tokens":385}},"tokens_in":802,"tokens_out":471,"duration_ms":5123,"temperature":1.0,"reasoning_tokens":385,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T19:42:39.251839+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $e^+e^-\\to\\gamma D\\bar D$ cross section near $\\sqrt{s}=4.23$ GeV and look for a $D\\bar D$ invariant-mass enhancement around 3.922 GeV whose photon polar-angle distribution in the $e^+e^-$ center-of-mass frame matches $1+\\cos^2\\theta$ (for $\\chi_{c0}$) or the nearly flat shape with $\\beta_2\\approx 0.11$ (for $\\chi_{c2}$); absence of such an enhancement would falsify the prediction. A separate decisive check is a precise measurement of $\\psi(4230)\\to\\pi^+ D^{*-}D^0$, which fixes $g_{\\psi D_1 D}$ and determines whether the claimed cross sections survive.","supporting_citations":[{"cited_title":"Is the X(3915) the chi_{c0}(2P)?","cited_arxiv_id":"1410.6534","evidence_quote":"establishes the $D\\bar D$ decay pattern and the two-substructure picture of the 3.9 GeV enhancement used for the branching ratios"},{"cited_title":"Understanding $B^-\\rightarrow X(3823)K^-$ via rescattering mechanism and predicting $B^-\\to \\eta_{c2} (^1D_2)/\\psi_3(^3D_3)K^-$","cited_arxiv_id":"1605.04776","evidence_quote":"gives the measured $\\Gamma_{e^+e^-}^{\\psi(4230)}\\cdot BR(\\psi(4230)\\to\\pi^+D^{*-}D^0)$ used to estimate the $g_{\\psi D_1 D}$ coupling"}],"review_version":1}