{"id":"c19aa7d4-5a8d-4fd4-82fc-852bce7cf9ec","arxiv_id":"1909.03901","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"If X(3872) is a charm-meson molecule, e+e- -> X gamma has a narrow triangle-singularity peak at 2.2 MeV above the D*0 Dbar*0 threshold, with a peak cross section near 0.5 pb.","lead":"This paper calculates how often electron-positron collisions should produce the X(3872) particle together with a photon, assuming X(3872) is a loose two-meson molecule. It predicts a narrow bump in the rate about 2.2 MeV above the D*0 Dbar*0 threshold, which the BESIII experiment could look for.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The predicted ~0.5 pb peak is only the triangle-diagram contribution; the short-distance e+e-→Xγ background is assumed negligible without quantitative estimate, though the authors concede a comparable background would turn the peak into a bump or dip.","rationale":"The paper's central claim, conditional on X(3872) being a weakly bound D*0 Dbar*0 molecule, is that e+e-→Xγ has a narrow triangle-singularity peak at W about 2.2 MeV with σ about 0.51 pb, and that observing it would support the molecular picture. The triangle loop amplitude F(W) is derived explicitly (Eqs. (13)-(22)), the singularity position follows from k=(µ/M0)q, and the comparison with the Dubynskiy-Voloshin absorptive contribution is internally consistent; I do not find a technical flaw in the loop calculation itself. The least secure condition is the normalization premise: the predicted peak is the triangle-diagram contribution only, and the short-distance amplitudes (creation of D Dbar, D* Dbar, or D Dbar* followed by binding into X) are assumed negligible. The Summary concedes that if these amplitudes are comparable, the cross section shows only a small peak or a dip on a smooth background, and it calls for quantitative estimates that are not provided. This premise is load-bearing because it decides whether the observable has a peak at all, not merely its height; the 0.2-0.9 pb band from varying A2/A0 is a band for the triangle contribution alone. The alternative concerns are real but secondary: the |A0|,|A2| fit in Section II carries no propagated errors, and the isospin-0 approximation is plausible given the broad ψ(4040) dominance. I agree with the reader's weakest assumption. A decisive check exists: BESIII Ref. [24] measured e+e-→Xγ at sqrt(s)=4.009 GeV, below the D*0 Dbar*0 threshold, so the triangle mechanism is kinematically absent and that point bounds the short-distance background. Since the paper already flags the limitation and states its prediction conditionally, the verdict CONDITIONAL stands; if the 4.009-GeV point or a targeted scan shows a background comparable to the predicted peak, the paper would need to be reframed as a shape-only statement or the observability claim withdrawn.","tokens_in":21006,"tokens_out":28770,"duration_ms":285253,"concrete_test":"Use the BESIII measurement at sqrt(s)=4.009 GeV (Ref. [24]), which lies about 4.7 MeV below the D*0 Dbar*0 threshold (4.0137 GeV); there the triangle mechanism is kinematically shut off, so the measured σ(e+e-→Xγ)×Br(J/ψπ+π-) bounds the short-distance background. Divide this product by the allowed branching fraction Br in [4%,33%] (Ref. [33]) to obtain the implied background σ_bg, then compare with the predicted peak σ of about 0.51 pb, minding the ψ(4040)-tail growth between 4.009 and 4.016 GeV. If the implied σ_bg band overlaps 0.51 pb, the neglect of short-distance amplitudes is falsified and the peak would be obscured; if the band sits far below 0.51 pb, the assumption is supported. A targeted BESIII scan across 4.010-4.020 GeV would settle the question definitively.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that, if X(3872) is a D*0 Dbar*0 molecule, e+e-→Xγ has a narrow triangle-singularity peak at 2.2 MeV above the D*0 Dbar*0 threshold with height about 0.51 pb (0.2-0.9 pb with the A2/A0 phase varied), and that its observation would strongly support the molecular identification. The predicted peak is the contribution of the charm-meson triangle diagrams in Fig. 4 only. As the Summary (Section VII) states, there are also short-distance amplitudes that begin with the creation of D Dbar, D* Dbar, or D Dbar* pairs (pointed out in Ref. [21]); these are essentially constant across the peak region. The paper assumes these are negligible ('We have assumed the short-distance amplitudes are negligible compared to the amplitude from the triangle diagrams. Quantitative estimates of the short-distance amplitudes would be useful.') and gives no estimate. The authors concede that if they are larger than the triangle amplitude, the cross section near the D*0 Dbar*0 threshold shows only a small peak or even a dip on a smooth background. Thus the quantitative prediction in the abstract and the BESIII observability conclusion stand or fall with this unquantified premise; the 0.2-0.9 pb band is a range for the triangle contribution, not for the full observable cross section. A direct check is available: BESIII (Ref. [24]) measured e+e-→Xγ at sqrt(s)=4.009 GeV, below the D*0 Dbar*0 threshold, where the triangle singularity is kinematically absent; this point bounds the short-distance background. Secondary limitations are the unpropagated errors in |A0|,|A2| from the Eq. (6) fit and the isospin-0 assumption, but the background question is the more load-bearing gap: it decides whether the observable has a peak at all.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper calculates the cross section for e+e- -> X(3872) gamma in the near-threshold region, assuming that X(3872) is a weakly bound D*0 Dbar0 molecule. The production proceeds through e+e- -> D*0 Dbar*0 followed by rescattering of the charm-meson pair into X gamma. The loop amplitude is reduced to an analytic scalar integral, and a triangle singularity produces a narrow peak 2.2 MeV above the D*0 Dbar*0 threshold. Using amplitudes fitted to Belle data on e+e- -> D*+D*-, the peak height is predicted to be about 0.51 pb, with a range of 0.2 to 0.9 pb when the A2/A0 phase is varied. The paper also compares the full amplitude with the absorptive contribution computed previously by Dubynskiy and Voloshin, and concludes that the absorptive contribution is not a good approximation near the peak.","tokens_in":21536,"tokens_out":5502,"duration_ms":55999,"significance":"If the result holds, it provides a distinctive, falsifiable line-shape prediction for X(3872) as a charm-meson molecule. The calculation is carried through analytically, with explicit cross-section formulas, and the peak position is derived from kinematics and is robust to the binding-energy scan. The paper is careful to disclose the A2/A0 phase ambiguity, to scan the binding energy, and to compare with the earlier Dubynskiy-Voloshin wavefunction prescription, including an instructive comparison between the full amplitude and its absorptive part. The main weakness, acknowledged in the text, is the unquantified short-distance background, which is load-bearing for the absolute normalization claim.","major_comments":[{"comment":"The predicted peak height is only the triangle-diagram contribution. In the Summary the authors state, \"We have assumed the short-distance amplitudes are negligible compared to the amplitude from the triangle diagrams. Quantitative estimates of the short-distance amplitudes would be useful.\" Because the abstract and the Summary present 0.51 pb and the 0.2-0.9 pb range as the prediction for the cross section, the central quantitative claim depends on an unquantified premise. This needs to be either quantified, for example by using the BESIII measurement at sqrt(s)=4.009 GeV from Ref. [24] where the triangle singularity is absent, or explicitly reframed as the triangle-only contribution.","section":"Section VII (Summary)"},{"comment":"The normalization of the cross section rests on |A0| and |A2| extracted from the Uglov et al. fit to e+e- -> D*+D*-, but the paper's own comparison with the alternative analysis of Du, Meissner, and Wang (Ref. [28]) shows a significantly different ratio of spin-2 to spin-0 cross sections (0.81 vs 2.92 at sqrt(s)=4.040 GeV). The quoted normalization range 0.47-1.80 only scans the A2/A0 phase with |A0|^2+|A2|^2 fixed; it does not include the uncertainties in that sum or in the amplitude ratio, so the uncertainty in the predicted peak height is underestimated.","section":"Section II, Eqs. (6)-(7)"},{"comment":"The identification of the neutral amplitudes A0 and A2 with the charged ones assumes that isospin-1 amplitudes are negligible. The paper argues from psi(4040) dominance, but no quantitative bound on the isospin-1 contamination is provided. Since the neutral amplitude is the difference of isospin-0 and isospin-1 amplitudes, a moderate isospin-1 amplitude could change the normalization beyond the quoted range, making this a load-bearing assumption for the absolute prediction.","section":"Section II, isospin assumption"}],"minor_comments":[{"comment":"There is a typo: \"triangle singularitiy\" should read \"triangle singularity\".","section":"Section VII"},{"comment":"In the sentence \"They did not measure the cross section at energies between 4.009 MeV and 4.178 MeV,\" the units should be GeV, not MeV.","section":"Section IV"},{"comment":"The phrase \"restricted to a singe quadrant\" should read \"restricted to a single quadrant\".","section":"Section VII"},{"comment":"The ratio in Eq. (7) is typeset as \"A 2/A0 = +/-1.9i\"; the spacing in the numerator is a formatting error and should be \"A2/A0\".","section":"Eq. (7)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a detailed companion to the authors' earlier letter (Ref. [22]) and is appropriate for a specialized hadron-physics journal. The main revision needed is to qualify or quantify the normalization claims, since the unquantified short-distance background is acknowledged in the text but the abstract and Summary present the triangle-only peak height as the predicted cross section. A quantitative bound from the existing BESIII data at 4.009 GeV would substantially strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know: this is the long version of a result that already appeared in the authors' PRD letter (Ref. [22]), and it is worth reading for the technical content even if you have seen the letter. The genuinely new pieces are the explicit analytic reduction of the triangle loop, the demonstration that the Dubynskiy-Voloshin absorptive approximation is quantitatively bad near the peak (height off by roughly 40%, position shifted by about 1.3 MeV), and a cleaner moving-frame wavefunction. The prediction itself is what it was: a narrow peak in e+e- -> X gamma at about 2.2 MeV above the D*0 Dbar*0 threshold, with height around 0.5 pb for the preferred amplitudes and 0.2-0.9 pb if the A2/A0 phase is allowed to roam.\n\nThe derivation is careful and the paper is unusually honest about its own limits. The loop integral is reduced to a scalar integral and evaluated analytically; the peak position follows from kinematics and known masses, so circularity is not a problem. The self-cited inputs (XEFT, wavefunction, branching-fraction bounds) are background tools, not the target result.\n\nThe main soft spot is exactly the one the authors flag in the Summary: the short-distance amplitudes that start with D Dbar, D* Dbar, or D Dbar* are assumed negligible, and no estimate is given. If they are comparable to the triangle amplitude, the peak can become a bump or a dip on a smooth background. That is load-bearing, and the stress-test note is right. There is a possible handle they do not use: BESIII measured e+e- -> X gamma at 4.009 GeV, below threshold, where the triangle singularity is kinematically absent, and that point could bound the short-distance background. A referee should ask for that or for an estimate based on the same short-distance production mechanism.\n\nThe other issues are more minor. The normalization from the Uglov et al. fit has no propagated errors; the authors' comparison with Du-Meissner-Wang suggests the total |A0|^2+|A2|^2 is stable at the few-percent level, but the A2/A0 ratio is less well determined and drives the 0.2-0.9 pb band. The isospin-1 suppression is plausible given psi(4040) dominance, but it is an assumption. Neither of these undermines the shape prediction; they only affect the absolute normalization.\n\nWho is this for? Anyone working on X(3872), triangle singularities, or e+e- charm-threshold physics. It deserves a serious referee. I would send it out and ask for the short-distance background estimate (or a shape-only framing) plus propagated errors, but I would not desk-reject it.","headline":"A careful, honest working-out of a triangle-singularity prediction for e+e- -> X gamma; the central caveat is the unquantified short-distance background, and the paper itself says so.","tokens_in":21951,"tokens_out":2592,"would_cite":true,"duration_ms":29609,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["14.80.Va","67.85.Bc","31.15.bt"],"model":"deepseek-v4-flash","headline":"Under the charm-meson-molecule picture of $X(3872)$, $e^+e^-$ annihilation into $X\\gamma$ develops a narrow, normalized peak from a triangle singularity, $2.2$ MeV above the $D^{*0}\\bar D^{*0}$ threshold at $\\sqrt{s}\\simeq 4.016$ GeV.","keywords":["X(3872)","charm-meson molecule","triangle singularity","e+e- annihilation","cross section","exotic hadrons","effective field theory","D* Dbar* threshold"],"falsifier":"Scan $e^+e^- \\to X(3872)\\gamma$ in steps of about 0.5 MeV across $\\sqrt{s}=4.010$ to $4.030$ GeV with enough integrated luminosity to resolve a cross section near 0.5 pb; absence of the predicted narrow peak 2.2 MeV above the $D^{*0}\\bar D^{*0}$ threshold would falsify the central claim, unless the short-distance background amplitudes turn out to dominate.","tokens_in":20811,"feed_emoji":"⚛️","tokens_out":12181,"duration_ms":113634,"temperature":0.7,"pith_summary":"The paper aims to prove a sharp production prediction: if the $X(3872)$ is a weakly bound charm-meson molecule, then $e^+e^-$ annihilation produces $X(3872)\\gamma$ with a narrow peak in the cross section, located $2.2$ MeV above the $D^{*0}\\bar D^{*0}$ threshold at $\\sqrt{s}\\approx 4.016$ GeV. The peak comes from a charm-meson triangle singularity, a kinematic point where all three charm mesons in the rescattering diagram are simultaneously on their mass shells. The predicted peak height is $0.51$ pb for the preferred amplitudes, and between roughly $0.2$ and $0.9$ pb as the ratio of the two P-wave production amplitudes is varied. This matters because the energy region has not yet been measured, so a dedicated scan could confirm or exclude the molecular interpretation of the $X(3872)$. The paper also establishes that the previously computed absorptive contribution alone is not a good approximation to the peak, since it sits about $1.3$ MeV higher and reaches only about $58\\%$ of the full peak height.","feed_headline":"Triangle singularity predicts a sharp 0.5 pb X(3872) gamma peak","feed_subtitle":"The peak sits 2.2 MeV above the D*0 Dbar*0 threshold, where no e+e- experiment has looked yet.","key_machinery":"The load-bearing object is the charm-meson triangle singularity, a kinematic singularity that arises when the three virtual particles forming a triangle diagram can all be on their mass shells at the same energy. Here it is realized in the two rescattering diagrams for $e^+e^- \\to X\\gamma$, and it enters through the loop amplitude $F(W)$, a scalar integral over the undetermined loop energy and momentum. The amplitude is reduced analytically to an integral over one variable and then evaluated in closed form; the singularity sits in the argument of a logarithm whose denominator vanishes at the triangle energy $W_\\Delta$. A secondary piece of machinery is the bound-state wavefunction for the $X$ in a moving frame, expressed in terms of the relative velocity of its constituents, which is used to compare the full amplitude with the earlier absorptive approximation.","core_discovery":"The central claim is that a triangle singularity governs $e^+e^- \\to X(3872)\\gamma$ near the $D^{*0}\\bar D^{*0}$ threshold whenever $X(3872)$ is a weakly bound charm-meson molecule. The virtual photon creates a P-wave $D^{*0}\\bar D^{*0}$ pair at short distances; one of the charm mesons radiates a photon; and the remaining $D^0$ or $\\bar D^0$ combines with the other charm meson to form the $X$. At one specific energy all three charm-meson lines in the triangle go on shell simultaneously, producing a logarithmic singularity in the loop amplitude $F(W)$. Evaluating the scalar loop integral analytically, the cross section $\\sigma[e^+e^- \\to X\\gamma]$ acquires a narrow peak at $W=2.2$ MeV, i.e. $\\sqrt{s}\\approx 4.016$ GeV, whose height is insensitive to the binding energy for $|E_X|$ between $0.10$ and $0.30$ MeV. The paper further shows that the absorptive part of the amplitude, which corresponds to on-shell $D^{*0}\\bar D^{*0}$ intermediate states, is not an adequate replacement for the full amplitude: it peaks about $1.3$ MeV higher and at only about $58\\%$ of the full peak height.","pith_inferences":["A dedicated energy scan in steps of order 0.5 MeV across $\\sqrt{s}=4.010$ to $4.030$ GeV would settle the prediction; the existing data gap makes this an immediately available experimental test.","The triangle-singularity mechanism is not channel-specific: any short-distance production of $D^{*0}\\bar D^{*0}$ can feed it, so analogous narrow peaks may appear in other final states, and the analytic loop amplitude here is a template for such predictions.","Precise line-shape data could also constrain the spin composition of P-wave charm-meson-pair production, because the normalization of the peak depends on the ratio $A_2/A_0$.","If no peak appears, the short-distance amplitudes that bypass the triangle would have to be large; quantifying those amplitudes is then necessary before ruling out the molecule picture, since the paper leaves that input unquantified."],"forward_implications":["The predicted peak sits at $\\sqrt{s}\\approx 4.016$ GeV, 2.2 MeV above the $D^{*0}\\bar D^{*0}$ threshold, in an energy gap that existing $e^+e^-$ measurements skipped; a fine scan there can test the molecular hypothesis.","The peak height is $0.51$ pb for the preferred $A_0,A_2$ amplitudes and ranges from about $0.2$ to $0.9$ pb under all complex values of $A_2/A_0$ consistent with the same $|A_0|^2+|A_2|^2$; its position is insensitive to the $X$ binding energy for $|E_X|$ between 0.10 and 0.30 MeV.","After multiplication by the $X\\to J/\\psi\\pi^+\\pi^-$ branching fraction, loosely bounded between 4% and 33%, the visible peak can be a sizable fraction of the $X\\gamma$ cross sections already measured at higher energies.","The absorptive contribution alone is not a good approximation in the peak region, so any extraction of the peak from data must use the full dispersive loop amplitude.","The same triangle-singularity peak appears with roughly the same shape for a zero-energy resonance or a virtual state, so observing the peak supports the molecular or resonant interpretation but would not by itself distinguish a narrow bound state from those alternatives."],"supporting_citations":[{"why":"Identifies the charm-meson triangle singularity that produces a narrow $X\\gamma$ peak near the $D^{*0}\\bar D^{*0}$ threshold in high-energy production.","marker":"[20]"},{"why":"Supplies the prior absorptive calculation whose peak this paper shows is not a good approximation to the full triangle-singularity line shape.","marker":"[21]"},{"why":"The earlier short report of this calculation, whose results are here derived in full detail.","marker":"[22]"},{"why":"Supplies the fit to measured $e^+e^- \\to D^{*+}D^{*-}$ cross sections used to fix the $A_0$ and $A_2$ amplitudes that set the peak normalization.","marker":"[23]"},{"why":"Gives existing $e^+e^- \\to \\gamma X(3872)$ measurements that skipped the predicted threshold region.","marker":"[24]"},{"why":"Provides the later measurements at coarser energies above the threshold, defining the gap where the peak would appear.","marker":"[25]"},{"why":"Fixes the measured $X$ mass relative to the $D^{*0}\\bar D^0$ threshold, determining the binding momentum $\\gamma_X$.","marker":"[13]"},{"why":"Provides the $D^*$ width and radiative width used for $\\Gamma_{*0}$ and the transition magnetic moment $\\nu$.","marker":"[29]"},{"why":"Gives the branching-fraction bounds on $X\\to J/\\psi\\pi^+\\pi^-$ used to estimate the observable height of the peak.","marker":"[33]"}],"fun_headline_variants":["Triangle singularity pins X(3872) gamma peak at 4.016 GeV","Narrow X(3872) gamma peak from charm-meson triangle loop","e+e- to X gamma: triangle singularity shapes the peak","Sharp X(3872) gamma resonance predicted in e+e- collisions","X(3872) gamma cross section peaks via triangle singularity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that short-distance amplitudes for $e^+e^- \\to X\\gamma$ that do not pass through the $D^{*0}\\bar D^{*0}$ triangle are negligible; if they are large, the predicted sharp peak becomes a small bump or even a dip on a smooth background.","fun_headline_variants_meta":{"raw":{"variants":["Triangle singularity pins X(3872) gamma peak at 4.016 GeV","Narrow X(3872) gamma peak from charm-meson triangle loop","e+e- to X gamma: triangle singularity shapes the peak","Sharp X(3872) gamma resonance predicted in e+e- collisions","X(3872) gamma cross section peaks via triangle singularity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000838,"raw_usage":{"total_tokens":3686,"prompt_tokens":1013,"completion_tokens":2673,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":629,"completion_tokens_details":{"reasoning_tokens":2578}},"tokens_in":629,"tokens_out":2673,"duration_ms":19879,"temperature":1.0,"reasoning_tokens":2578,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:45:21.991926+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Scan $e^+e^- \\to X(3872)\\gamma$ in steps of about 0.5 MeV across $\\sqrt{s}=4.010$ to $4.030$ GeV with enough integrated luminosity to resolve a cross section near 0.5 pb; absence of the predicted narrow peak 2.2 MeV above the $D^{*0}\\bar D^{*0}$ threshold would falsify the central claim, unless the short-distance background amplitudes turn out to dominate.","supporting_citations":[{"cited_title":"Tanabashi et al","cited_arxiv_id":null,"evidence_quote":"Supplies the prior absorptive calculation whose peak this paper shows is not a good approximation to the full triangle-singularity line shape."},{"cited_title":"Production of $X(3872)$ Accompanied by a Soft Pion at Hadron Colliders","cited_arxiv_id":"1903.04355","evidence_quote":"The earlier short report of this calculation, whose results are here derived in full detail."},{"cited_title":"Triangle Singularity in the Production of $X(3872)$ and a Photon in $e^+e^-$ Annihilation","cited_arxiv_id":"1904.12915","evidence_quote":"Provides the $D^*$ width and radiative width used for $\\Gamma_{*0}$ and the transition magnetic moment $\\nu$."}],"review_version":1}