{"id":"1e345e36-90bb-4d27-b2c6-7d660d4008d5","arxiv_id":"1908.05710","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":2.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A coupled-channel extension of the unquenched quark model with one per-multiplet subtraction parameter predicts that χc(2P) states, especially the X(3872), have large molecular components while χb(3P) states are almost pure bottomonia.","lead":"This paper summarizes a quark-model approach that adds meson-meson loops, or threshold corrections, to the wave functions of heavy quarkonium states, and applies it to charmonium and bottomonium multiplets. The approach is presented as a fix for the lack of convergence in older unquenched quark-model calculations and as a way to tell ordinary quarkonia from molecules and other exotic states.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Renormalization factor in Eq. (7) appears to have the wrong sign, making the Table II continuum probability inconsistent with the Table I mass shifts.","rationale":"The reader's weakest-assumption concerned the experimental identification of χc0(3915) as the 2P state. That is a legitimate issue for the mass-pattern claim, but it is an external assignment ambiguity. The concern raised here is internal to the formalism: Eq. (7), as printed, rescales the pair-creation coupling in the wrong direction relative to the subtraction in Eq. (5). Because Table II's P_cont = 0.853 is the central quantitative evidence for the abstract's claim that χc(2P) states have non-negligible molecular components, this inconsistency is more load-bearing than the assignment question for the paper's headline result. The concrete test is a direct recomputation with a clearly specified sign correction, so it can settle the matter without new experimental input. If the corrected P_cont remains above 0.5, the qualitative conclusion stands but the quantitative value and the text need revision; hence CONDITIONAL rather than REJECT. The paper's other limitations (proceedings-style summary, derivation deferred to Ref. [47], use of a per-multiplet subtraction parameter) are already reflected in the reader's verdict and do not change this assessment.","tokens_in":9584,"tokens_out":15789,"duration_ms":153883,"concrete_test":"Recompute R_i from Eq. (7) as printed and with the corrected sign, using the raw self-energies implied by Table I and Δ = +16 MeV. For each multiplet member, evaluate Eq. (6) with γ̃_i = γ√R_i using the same vertices and form factors as Ref. [47], and compare the resulting P_cont with Table II. Then substitute the renormalized coupling back into Eq. (4) and check whether the resulting self-energy equals the Table I value: the printed sign gives −97 MeV for χc1 instead of −65 MeV, while the corrected sign gives −65 MeV. If recomputation confirms the mismatch, Table II's 0.853 value is not the continuum probability consistent with the reported mass shifts, and the numbers must be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Eq. (7) defines R_i = (Σ(M_i) − Δ)/Σ(M_i), and Eq. (8) sets the renormalized pair-creation strength γ̃_i = γ√R_i, which is then used in Eq. (6) to compute P_cont. Under the text's definition of Δ as the smallest |Σ| among multiplet members, Table I implies raw self-energies Σ(hc) = −32, Σ(χc0) = −16, Σ(χc1) = −81, Σ(χc2) = −46 MeV, with Δ = +16 MeV. For χc1, Eq. (7) gives R = (−81 − 16)/(−81) = 1.198, so γ̃²/γ² = 1.198. The self-energy corresponding to this stronger coupling is R·Σ = −97 MeV, not the −65 MeV shift reported in Table I. The intended renormalization should reduce the coupling: R = (Σ + Δ)/Σ = 0.802 reproduces −65 MeV. As written, Table II's P_cont = 0.853 is therefore computed with a coupling inconsistent with the mass-shift calculation. With the corrected sign, the continuum probability would be roughly 0.853 × 0.802 / 1.198 ≈ 0.57. The qualitative claim of a dominant molecular component survives, but the headline number in Table II is suspect and needs recomputation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a UQM-based coupled-channel approach to heavy quarkonium-like states, in which bare quark-model masses E_A are corrected by meson-meson loop self-energies Σ(M_A) plus a per-multiplet subtraction Δ. It applies the method to the χc(2P) and χb(3P) multiplets, reporting mass shifts that move χc1(2P) toward X(3872) and describing χc(2P) states as having significant molecular components (P_cont=0.853 for X(3872)) while χb(3P) states remain almost pure. It also sketches the calculation of open-flavor and hidden-flavor decay widths, quoting a J/ψω to J/ψρ ratio compatible with data. The paper is a concise, reference-heavy summary of a published companion calculation rather than a full self-contained derivation.","tokens_in":9889,"tokens_out":11627,"duration_ms":109938,"significance":"If the model and the numerical results are correct, the paper would provide a simple renormalization prescription for UQM loop calculations and a concrete diagnostic for quarkonium-like versus molecular states. The central idea is valuable, and the comparison with experiment at the 30–50 MeV level is a genuine strength. However, the quantitative claim about the X(3872) continuum probability and the methodological claim about convergence are currently undermined by an algebraic inconsistency in Eq. (7) and by unresolved normalization and assignment issues. The paper provides no machine-checked proofs or code; its numerical input is transparently traced to cited models and to Ref. [47].","major_comments":[{"comment":"The sign convention in Eq. (7) is inconsistent with Eq. (5) and Table I. For the χc(2P) multiplet, Table I implies raw self-energies Σ ≈ (−32, −16, −81, −46) MeV for hc, χc0, χc1, χc2, with Δ = +16 MeV if Σ + Δ is the tabulated quantity. Then for χc1, R = (Σ−Δ)/Σ = (−81−16)/(−81) ≈ 1.198, so Eq. (8) increases the pair-creation strength and would produce a self-energy of about −97 MeV, not the −65 MeV reported. If instead one takes Δ = −16 MeV, Eq. (7) does give R ≈ 0.802, but then Σ + Δ for χc0 would be −32 MeV, contradicting the zero entry in Table I. Thus the equations and the table are mutually inconsistent. Recomputing P_cont with the rescaling that reproduces the −65 MeV shift gives P_cont ≈ 0.82 rather than 0.853, so the qualitative molecular-dominance conclusion survives, but the headline value in Table II must be recalculated.","section":"III B, Eqs. (7)–(8) and Tables I–II"},{"comment":"The probability definition is ambiguous: Eq. (6) writes P_cont as an unnormalized sum over 1/(M_A−E_B−E_C)^2, while the wave function in Eq. (1) contains an explicit normalization factor N. The paper does not state whether Table II includes N or whether the renormalized γ̃ of Eq. (8) is intended to absorb it. Since P_val = 1 − P_cont is used, the normalization convention is consequential and should be stated explicitly.","section":"III B, Eq. (6) and Table II"},{"comment":"The interpretation of the χc(2P) multiplet rests on a debated assignment. The table labels the 0++ candidate as 'χc0(3915) or χc0(2P)', but the argument that threshold effects break the χ-multiplet pattern, and the choice of Δ itself, depend on identifying that state as the 2P partner. If χc0(3915) is not χc0(2P), the smallest self-energy, Δ, and all R_i would change, and the pattern-breaking claim would not be established. The paper should either justify the assignment or present the conclusions as conditional on it.","section":"III A, Table I"},{"comment":"The renormalization prescription is asserted but not derived. The text states that a single Δ removes the convergence problem of UQM, but no argument or numerical demonstration shows that the physical results are stable as the tower of intermediate states is extended. Since the novelty of the paper is precisely this prescription, a convergence test (e.g., M_A and P_cont as functions of the number of included channels) is needed to support the central methodological claim.","section":"III B, Eqs. (7)–(8)"}],"minor_comments":[{"comment":"Figure 1 is referred to in the text, but only the caption is present in the manuscript; the plot itself is missing.","section":"III A"},{"comment":"There are a few grammatical slips, e.g., 'It is worth to observe' should be 'It is worth observing'.","section":"III A"},{"comment":"Equation (12) is called a ratio of amplitudes, but the written object is a ratio of decay widths; the text should say widths.","section":"IV A, Eq. (12)"},{"comment":"The definitions of γ0^eff and the vertex form factors are only given by reference to Ref. [59]; for a self-contained journal submission, these definitions should be reproduced or at least summarized.","section":"III B"}],"recommendation":"major_revision","confidential_remarks":"A practical check for the revision: compare Eq. (7) with the published version (Ferretti and Santopinto, PLB 789, 550 (2019)). If the sign in Eq. (7) is a transcription error, the fix is straightforward but Table II must be recomputed; if the published version has the same sign, the published result itself needs scrutiny. I also note that this manuscript is essentially a proceedings summary of that PLB paper; the incremental content for a refereed journal should be clarified with the editor."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Jacopo,\n\nYou can read this in half an hour. It's a proceedings write-up of the UQM-based coupled-channel model from Ferretti and Santopinto (PLB 789, 550). The genuinely useful thing here is the renormalization idea: subtract the smallest self-energy in each multiplet to get a convergent expansion, then apply the same subtraction to the continuum probability. That is a practical fix to a known problem, and the paper explains it clearly enough to reproduce.\n\nThe comparison between chi_c(2P) and chi_b(3P) is also interesting – the former getting large molecular components, the latter staying essentially pure bottomonium. The table of mass shifts is consistent with experiment at the 30-50 MeV precision claimed. And the ratio of J/psi-omega to J/psi-rho decay widths for the X(3872) is nicely close to the measured 0.8±0.3.\n\nNow the soft spots. The obvious one is Eq. (7). As written, R_i = (Sigma_i – Delta)/Sigma_i. With Sigma_i negative and Delta positive, this gives R_i > 1, which means the renormalized pair-creation strength gamma-tilde is larger than the bare gamma. That would make the self-energy more negative, not less. The whole point of Delta is to reduce the mass shift, so the effective self-energy is Sigma_i + Delta. The ratio that reproduces that shift is (Sigma_i + Delta)/Sigma_i, which is < 1. For the numbers in Table I, R for chi_c1 would be about 0.80, not 1.20. So the continuum probability in Table II, 0.853, is computed with the wrong renormalization. Re-scaling gives P_cont ≈ 0.57. The qualitative conclusion – X(3872) has a sizable molecular component – survives, but the headline number is off by a factor of ~1.5.\n\nThe other soft spots are milder. This is explicitly a summary of Ref. [47]; there are no new derivations or calculations. The assignment of chi_c0(3915) to the 2P multiplet is assumed without discussing the alternative. There are no uncertainties on the predicted masses, and the decay width section is a quote of previous work, with the actual coupled-channel decay amplitudes left to a future paper.\n\nWho should read this? Anyone who wants a quick, guided tour of the unquenched quark model and this particular coupled-channel trick. It's a fine proceedings piece. But as a standalone research paper it doesn't carry its weight – the sign error alone means the quantitative claim in Table II can't be trusted. If this lands on a referee's desk, I'd ask for the derivation and the corrected renormalization factor, and for an explicit statement that this is a summary, not new research. I would not cite it in place of the original.\n\nCheers.","headline":"Useful proceedings-style summary of a published coupled-channel model, but the renormalization factor in Eq. (7) looks sign-inconsistent and the paper adds no new results.","tokens_in":10422,"tokens_out":6176,"would_cite":false,"duration_ms":55131,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Threshold meson–meson loops make X(3872) a mostly molecular state","keywords":["heavy quarkonium spectroscopy","threshold effects","unquenched quark model","coupled-channel model","X(3872)","charmonium","bottomonium","hadronic molecules"],"falsifier":"Measure the hc(2P) state: the model explicitly assumes its physical mass equals its bare mass of 3956 MeV, so finding hc(2P) shifted by tens of MeV would directly falsify the multiplet-subtraction assumption; alternatively, a lattice QCD calculation of the χc1(2P) self-energy that gives a shift much smaller than the claimed −65 MeV would rule out the threshold-dominance picture.","tokens_in":2144,"feed_emoji":"⚛️","tokens_out":6693,"duration_ms":101381,"temperature":0.7,"pith_summary":"This paper tries to establish that threshold effects—virtual meson–meson loops dressing a bare quark–antiquark state—matter quantitatively for heavy quarkonium, and that a renormalized coupled-channel version of the Unquenched Quark Model can handle them without the usual convergence trouble. Applied to the χc(2P) multiplet, the model identifies X(3872) as a charmonium-like state with a dominant molecular-type component (85.3% continuum probability) and produces a self-energy shift of −65 MeV for its χc1 member, moving the bare mass of 3953 MeV toward the experimental 3871.69 MeV. In the bottomonium χb(3P) multiplet, the same corrections are only a few MeV, so those states stay almost pure quarkonia. The paper also uses the formalism for hidden-flavor decays, predicting a J/ψω to J/ψρ ratio for X(3872) of 0.6, compatible with the measured 0.8 ± 0.3.","feed_headline":"Meson-meson loops turn X(3872) into a mostly molecular state","feed_subtitle":"Threshold shifts of up to −65 MeV reshape the χc(2P) multiplet but leave bottomonium nearly pure.","key_machinery":"The central object is the self-energy $\\Sigma(M_A)=\\sum_{BC\\ell J}\\int_0^\\infty k^2dk\\,|\\langle BCk\\ell J|T^\\dagger|A\\rangle|^2/(M_A-E_B-E_C)$, computed with a $^3P_0$ pair-creation operator $T^\\dagger$ that couples the bare quarkonium state $|A\\rangle$ to meson–meson continuum states $|BC\\rangle$. The paper's renormalization prescription replaces $M_A=E_A+\\Sigma(M_A)$ by $M_A=E_A+\\Sigma(M_A)+\\Delta$, where $\\Delta$ is the smallest self-energy correction among the multiplet members, which removes the convergence problem of the standard UQM. The molecular fraction is controlled by the continuum probability $P^{\\rm cont}_A=\\sum_{BC\\ell J}\\int q^2dq\\,|\\langle BCq\\ell J|T^\\dagger|A\\rangle|^2/(M_A-E_B-E_C)^2$, evaluated with a rescaled pair-creation strength $\\tilde\\gamma^{\\rm eff}_{0,i}=\\gamma^{\\rm eff}_0\\sqrt{R_i}$, where $R_i=(\\Sigma(M_i)-\\Delta)/\\Sigma(M_i)$.","core_discovery":"The central claim is that a UQM-based coupled-channel model, with one subtraction constant per multiplet, can reproduce the masses of the χc(2P) and χb(3P) quarkonia including threshold corrections. The key quantitative results are that the χc1(2P), identified with X(3872), receives a −65 MeV self-energy shift and acquires a 0.853 continuum (D D̄, D D̄*, ...) probability, while the other χc(2P) members shift by −30 MeV or 0, breaking the usual χ-multiplet mass pattern Mχ0 < Mχ1 ≈ Mh < Mχ2. In contrast, the χb(3P) members shift by only −2 to −7 MeV, leaving the bottomonium multiplet pattern intact and the states almost pure bb̄. The paper interprets this contrast as evidence that threshold effects selectively turn some charmonium states into molecular-like objects, while bottomonium states remain essentially unquenched.","pith_inferences":["If the 85.3% molecular fraction is correct, X(3872) should sit very close to the D0D̄*0 threshold and have a large D0D̄*0 scattering length; a high-precision measurement of its line shape near threshold could test this picture directly.","Applying the same multiplet-by-multiplet subtraction scheme to other quarkonia, such as χc(3P) or ψ(4S)-like states, would make a concrete prediction about which states should show strong threshold distortions and which should stay almost pure quarkonia.","The single free parameter Δ per multiplet is a testable modelling choice: replacing it with a fitted subtraction constant determined from more states would show whether the scheme is a general renormalization procedure or an ad hoc calibration.","The assumed equality of the unobserved hc(2P) physical mass with its bare mass is a strong constraint; a future measurement of hc(2P) that shows a sizable shift would require modifying the multiplet pattern and the interpretation of the threshold effects."],"forward_implications":["X(3872) is predicted to be predominantly a D0D̄*0-like molecular state, with only a 14.7% valence charmonium core, so its internal structure is far from a plain c c̄ state.","Threshold corrections of up to −65 MeV are large enough to break the usual spin-mass ordering of the χc(2P) multiplet, meaning simple quark-model mass patterns can be unreliable when open-flavor thresholds are nearby.","The χb(3P) multiplet remains essentially a pure bottomonium multiplet, with threshold shifts of only a few MeV, so bottomonium spectroscopy is much less affected by meson–meson loops.","The same formalism predicts the hidden-flavor decay ratio Γ(X → J/ψω)/Γ(X → J/ψρ) = 0.6, which agrees with the measured 0.8 ± 0.3, suggesting the molecular component can simultaneously explain the mass and the decay pattern of X(3872).","The renormalized coupled-channel prescription, subtracting the smallest self-energy of a multiplet, provides a route to computing other observables, such as open-flavor strong decay amplitudes, without the convergence issues that plagued earlier UQM calculations."],"supporting_citations":[{"why":"Supplies the UQM-based coupled-channel model, the subtraction prescription, and all the numerical results for the χc(2P) and χb(3P) mass shifts and continuum probabilities.","marker":"[47]"},{"why":"Provides the standard UQM formalism, including the pair-creation operator and the self-energy calculation that the coupled-channel approach builds on.","marker":"[42, 43]"},{"why":"Supplies the bare quark-model masses E_A for the charmonium and bottomonium states used as input to the coupled-channel calculation.","marker":"[64]"},{"why":"Provides the experimental masses and the measured decay ratio against which the model predictions are compared.","marker":"[4]"},{"why":"Defines the effective pair-creation strength and the continuum-probability formula used to compute the 0.853 molecular fraction of X(3872).","marker":"[59, 60]"},{"why":"Provides the 3P0 pair-creation model used for strong and hidden-flavor decay amplitudes in the UQM framework.","marker":"[63]"},{"why":"Gives the experimental X(3872) → J/ψω to J/ψρ ratio used to validate the decay-rate calculation.","marker":"[73]"}],"fun_headline_variants":["X(3872) goes molecular from 65 MeV threshold shift","Threshold loops make X(3872) 85% molecular, bottomonia pure","Coupled channels: X(3872) mostly molecule, bottomonium not","Unquenching splits quarkonia: X(3872) molecular, χb pure","How threshold effects turn X(3872) into a molecule"],"cache_read_input_tokens":12416,"weakest_assumption_plain":"The calculation treats the observed χc0(3915), χc1(3872), and χc2(3930) as the three spin partners of the χc(2P) multiplet and takes the unobserved hc(2P) mass to be its bare value; if that experimental classification is wrong, the claimed threshold-induced breaking of the multiplet mass pattern is not established.","fun_headline_variants_meta":{"raw":{"variants":["X(3872) goes molecular from 65 MeV threshold shift","Threshold loops make X(3872) 85% molecular, bottomonia pure","Coupled channels: X(3872) mostly molecule, bottomonium not","Unquenching splits quarkonia: X(3872) molecular, χb pure","How threshold effects turn X(3872) into a molecule"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00013,"raw_usage":{"total_tokens":1135,"prompt_tokens":967,"completion_tokens":168,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":69}},"tokens_in":583,"tokens_out":168,"duration_ms":2556,"temperature":1.0,"reasoning_tokens":69,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:06:45.158465+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the hc(2P) state: the model explicitly assumes its physical mass equals its bare mass of 3956 MeV, so finding hc(2P) shifted by tens of MeV would directly falsify the multiplet-subtraction assumption; alternatively, a lattice QCD calculation of the χc1(2P) self-energy that gives a shift much smaller than the claimed −65 MeV would rule out the threshold-dominance picture.","supporting_citations":[{"cited_title":"Ferretti and E","cited_arxiv_id":null,"evidence_quote":"Supplies the UQM-based coupled-channel model, the subtraction prescription, and all the numerical results for the χc(2P) and χb(3P) mass shifts and continuum probabilities."},{"cited_title":"del Amo Sanchez et al","cited_arxiv_id":null,"evidence_quote":"Gives the experimental X(3872) → J/ψω to J/ψρ ratio used to validate the decay-rate calculation."}],"review_version":1}