{"id":"efb6ae73-c59e-4921-a7d6-dbe0c574dd7b","arxiv_id":"2608.03104","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Using one-boson-exchange Bethe-Salpeter equations, the authors find possible S-wave bound states in isoscalar D*Dbar*/B*Bbar* and in the 0(1+) and 1(2+) doubly heavy systems, with the isovector hidden-heavy channels unbound within the parameter scan.","lead":"This paper uses the Bethe-Salpeter equation with meson exchange to ask whether pairs of heavy vector mesons can bind together. It finds candidate bound states in several channels, but the predictions depend heavily on a free cutoff parameter that the model does not fix.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"D*D*/B*B* kernel is t-channel-only for identical bosons; if the u-channel exchange is omitted, the doubly-heavy bound-state solutions are not Bose-symmetric.","rationale":"The reader's weakest_assumption correctly flags the uncalibrated cutoff alpha as limiting the physical interpretation. I agree that the alpha dependence is a real limitation. However, the more load-bearing issue for the paper's internal correctness is the treatment of identical particles in the D*D* and B*B* sectors. The BS equation for two identical bosons must be solved with a Bose-symmetric kernel; the paper's Eqs. (14)-(16) contain only the t-channel one-boson-exchange diagram and no u-channel exchange. Because the constituent mesons are identical, the exchange diagram contributes at the same order and is not negligible. Omitting it means the equation does not respect the interchange symmetry of the wavefunction, so the reported solutions in Sec. III.B may not be physical. This affects the core claim for the doubly heavy systems. The hidden-heavy sector is not affected by this issue. The concrete test of adding the u-channel term (or comparing with a published BS calculation that includes it, e.g., Ref. [45]) would settle the matter. The paper is otherwise transparent about the model assumptions and the alpha dependence. Therefore I keep the reader's CONDITIONAL verdict, but the condition should include a demonstration that the symmetrized kernel produces the same doubly-heavy bound states.","tokens_in":13637,"tokens_out":17524,"duration_ms":171208,"concrete_test":"Recompute the D*D* (and B*B*) spectrum using the symmetrized kernel K_sym(p,q) = K(p,q) + K(p,-q), where the second term is the u-channel diagram with the appropriate permutation of the vector and isospin indices, keeping all masses, couplings, and the form factor (18) unchanged. If the I(JP) = 0(1+) and 1(2+) bound states disappear or need a substantially different alpha, the doubly-heavy claim is not robust. A simpler check is to reproduce the D*D* results of Ref. [45], which is a BS calculation that includes crossed diagrams.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim for the doubly heavy systems rests on the BS equation for two identical vector mesons, Eqs. (14)-(16). For identical bosons, the irreducible two-body kernel must include both the direct (t-channel) and crossed (u-channel) one-boson-exchange diagrams; the crossed diagram has a different momentum transfer (p+q versus p-q) and a different contraction of vector indices. The paper presents only the t-channel kernels and does not state that the kernel was symmetrized or that a factor of 1/2 was used. As written, the equation is not invariant under p -> -p with the vector indices exchanged, which is required for a Bose-symmetric bound-state amplitude (the flavor wavefunctions in Eq. (3) establish the symmetry of the state, but the kernel must be symmetrized too). Without the u-channel, the solutions in Figs. 3 and 4 for I(JP) = 0(1+), 1(0+), 1(2+) are not guaranteed to be physical two-boson bound states. This is an internal completeness issue, independent of the uncalibrated cutoff alpha; it affects the entire doubly-heavy sector of the paper's conclusions.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies S-wave molecular bound states in the hidden-heavy systems D*Dbar* and B*Bbar* and in the doubly heavy systems D*D* and Bbar*Bbar* within the Bethe-Salpeter formalism. The authors construct one-boson-exchange kernels (sigma, pi, eta, rho, omega) in the ladder and instantaneous approximations, introduce a monopole form factor with cutoff Lambda = m + alpha Lambda_QCD, reduce the four-dimensional BS equations to one-dimensional eigenvalue problems, and solve them numerically. They report bound-state solutions for the isoscalar hidden-heavy channels with J^PC = 0^{++}, 1^{+-}, and 2^{++}, no isovector solutions in the scanned parameter range, and doubly heavy solutions in I(J^P) = 0(1^+), 1(0^+), and 1(2^+), with the 1(0^+) channel requiring a comparatively large cutoff. The bottom analogues are found to bind at smaller cutoff values because of their larger reduced masses.","tokens_in":13833,"tokens_out":10389,"duration_ms":105399,"significance":"If the results survive the technical issues below, the paper would provide a systematic one-boson-exchange BS survey of the D*/B* vector-meson pair systems and a useful comparison with earlier meson-exchange, EFT, complex-scaling, and lattice studies. The authors deserve credit for being explicit that alpha is not determined from first principles and that large-alpha solutions should be interpreted cautiously; the comparison with the existing literature is broad and helpful, and the enumeration of Bose-allowed quantum numbers for the identical-particle systems is correct. The main barriers are the unsymmetrized kernel for identical bosons, the lack of quantitative cutoff thresholds, and the absence of the reduced equations needed to reproduce the numerical solutions.","major_comments":[{"comment":"The doubly heavy sector is not Bose-symmetric as written. For two identical D* or B* mesons, the irreducible two-body kernel in Eq. (14) must be symmetric under particle exchange. Equation (16) contains only t-channel one-boson-exchange diagrams; the crossed u-channel diagrams, obtained by exchanging the final-state momenta and interchanging the vector indices, are neither included nor stated to be negligible. The flavor wave functions in Eq. (3) determine the symmetry of the state, but the kernel itself must also be symmetrized, so that the equation is invariant under p -> -p together with an exchange of the vector indices. As it stands, the solutions shown in Figs. 3 and 4 cannot be identified with physical two-boson bound states, and the conclusion is independent of the regulator alpha. The hidden-heavy sector is unaffected because D* and Dbar* are distinguishable, but the doubly heavy analysis in Secs. III B and IV must be redone with a Bose-symmetrized kernel or with an explicit argument that the u-channel contribution vanishes.","section":"Sec. II B, Eqs. (14)-(16)"},{"comment":"The predictive content of the paper is weakened by the fact that alpha is a free scan parameter (0.5 <= alpha <= 10) and that threshold values are not reported. Each claimed bound state appears only above some alpha, and the distinction between 'moderate' and 'substantially large' cutoffs is never quantified. Please add a table, or explicit in-text values, of alpha_min and Lambda_min for every channel shown in Figs. 1-4, and state the criterion used to judge a cutoff as reasonable. Without these numbers, the reader cannot quantitatively compare the present results with the Lambda around 0.5 GeV used in Refs. [15,54,62], and the claim that the bottom systems bind more favorably is not quantitatively testable.","section":"Sec. II B, Eq. (18), and Sec. III"},{"comment":"The reduction from the four-dimensional BS equation to the one-dimensional integral equations is not shown. The text states that the p_l integration is performed by contour integration and the azimuthal integration analytically, but the resulting projected kernels and the discrete matrix eigenvalue problem are not given. Without these expressions, or a clear pointer to a companion derivation, the numerical solutions in Figs. 1-4 cannot be checked. Please provide an appendix or supplemental material with the reduced equations, the quadrature rule, and a convergence test.","section":"Sec. III, paragraph on numerical reduction"}],"minor_comments":[{"comment":"The text contains grammatical errors: 'is express in terms' and 'can be express as' should be 'is expressed in terms' and 'can be expressed as'.","section":"Sec. II A, around Eq. (8)"},{"comment":"The allowed quantum numbers are listed as 'I(JP) = 0(1+), 1(0+) 1(2+)'; a comma is missing between 1(0+) and 1(2+).","section":"Sec. III B"},{"comment":"The phrase 'the 0(1+) and 1(2+) solutions' is ambiguous; please write 'the I(J^P) = 0(1^+) and I(J^P) = 1(2^+) channels'.","section":"Sec. IV"},{"comment":"The text says the Lorentz tensor structures are 'exactly the same' as in Eq. (11); while the same tensor basis is used, the charge-conjugation properties differ between the hidden-heavy J = 1 channel and the doubly heavy J = 1 channel. Please clarify that the identity is at the level of the tensor basis only.","section":"Eq. (17)"},{"comment":"The phrase 'the bottom systems are bounded more favorably' should be reworded, for example as 'the bottom systems bind more readily' or 'form bound states at smaller cutoff values'.","section":"Abstract and Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The u-channel omission for the identical-boson systems is the decisive technical issue; if the authors can demonstrate that the u-channel contributions cancel in the adopted approximation, or if they redo the calculation with a symmetrized kernel, the doubly heavy sector could be made sound. The free-cutoff issue is not fatal by itself because the dependence is displayed, but the paper should state the alpha_min values and define what counts as a reasonable cutoff. I would not recommend acceptance before these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"One thing to know: this paper is a standard BS calculation with OBE kernels, and the qualitative results are not new. The hidden-heavy sector—isoscalar 0++, 1+-, 2++ bound states, no isovector—matches earlier OBE and BS studies. The doubly-heavy sector also reproduces previously reported 0(1+) and 1(2+) candidates. What the paper does well: the derivation is clean, the numerical reduction to one-dimensional eigenvalue equations is standard, and the authors are transparent about the free cutoff alpha. They show binding energy vs alpha curves, note that the 1(0+) state needs a large cutoff, and explicitly say alpha cannot be fixed from first principles. That honesty is worth something.\n\nThe soft spots are real. The entire bound-state inventory is a scan over alpha, with no external calibration and no error bars. Since every channel binds only above some alpha threshold, the predictions are conditional on short-distance physics that the OBE model does not constrain. That is not fatal—this is how most OBE molecule papers operate—but it does mean the paper is not a sharp prediction. The authors half-say this themselves.\n\nThe bigger technical concern is in the D*D*/B*B* section. For identical bosons, the two-body kernel must be symmetrized; the direct t-channel diagram alone is not enough. The paper presents only t-channel kernels (Eq. 16) and never mentions a u-channel exchange or a symmetrization procedure. If that is really all they solved, the doubly-heavy bound-state solutions in Figs. 3 and 4 are not obviously physical two-boson states. This is an internal completeness problem, independent of the cutoff, and it deserves a direct answer from the authors.\n\nThe citation pattern is okay—they cite the relevant OBE, BS, lattice, and EFT papers—but the novelty claim is modest. Refs [42], [46], [65], [17], and [68] already report the same qualitative results. I would not cite this paper for a new result; I might cite it as one more data point in a model comparison.\n\nRecommendation: send it to peer review, but with a request to fix/symmetrize the identical-boson kernel and to discuss cutoff sensitivity more quantitatively. For a reading group, only if someone is working on BS molecules.","headline":"A competent but incremental Bethe-Salpeter scan of four heavy-vector-meson systems; the hidden-heavy results are consistent with earlier OBE work, while the doubly-heavy section may be missing the u-channel exchange required for identical bosons.","tokens_in":14423,"tokens_out":3230,"would_cite":false,"duration_ms":30242,"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":"Heavy vector-meson pairs can bind as S-wave hadronic molecules, a Bethe-Salpeter one-boson-exchange calculation concludes.","keywords":["hadronic molecules","Bethe-Salpeter equation","one-boson-exchange interaction","heavy vector mesons","S-wave bound states","exotic hadrons","heavy-quark spin symmetry","molecular states near threshold"],"falsifier":"A lattice QCD calculation of $S$-wave $D^*\\bar D^*$ scattering with $I=0$ and $J^{PC}=0^{++}$, $1^{+-}$, $2^{++}$ that finds no near-threshold pole would directly contradict the central prediction; likewise, a coupled-channel meson-exchange calculation that removes these bound states once $D\\bar D$ and $D\\bar D^*$ channels are included would show the one-boson-exchange-only result is an artifact of the truncated kernel.","tokens_in":13436,"feed_emoji":"⚛️","tokens_out":14361,"duration_ms":110796,"temperature":0.7,"pith_summary":"The paper asks whether two heavy vector mesons—$D^*\\bar D^*$, $B^*\\bar B^*$ (one heavy and one anti-heavy partner) and $D^*D^*$, $\\bar B^*\\bar B^*$ (doubly heavy)—can form zero-angular-momentum ($S$-wave) molecular bound states. Treating their interaction through exchanges of $\\sigma$, $\\pi$, $\\eta$, $\\rho$, and $\\omega$ mesons and solving the Bethe-Salpeter equation in ladder and instantaneous approximations, it finds bound-state solutions in the total-isospin-zero (isoscalar) hidden-heavy channels $J^{PC}=0^{++}$, $1^{+-}$, and $2^{++}$, and no isovector solutions in the parameter range considered. In the doubly heavy systems, solutions appear for $I(J^P)=0(1^+)$, $1(0^+)$, and $1(2^+)$, with the $1(0^+)$ state requiring an unusually large cutoff and therefore judged less certain. Bottom systems bind more easily than charmed ones because their larger reduced masses lower the kinetic cost of binding. If these states exist, they would be near-threshold molecular partners of states like $X(3872)$ and the $Z_b$ resonances, and would broaden the experimental search for exotic hadrons.","feed_headline":"Heavy vector meson pairs bind as molecules, model says","feed_subtitle":"A Bethe-Salpeter calculation predicts near-threshold molecular bound states, with bottom analogues binding most readily.","key_machinery":"The load-bearing machinery is the Bethe-Salpeter equation for two vector mesons, with a covariant $S$-wave wave function for each allowed $J^{PC}$ or $J^P$, and a kernel built from $t$-channel one-boson exchange of $\\sigma$, $\\pi$, $\\eta$, $\\rho$, and $\\omega$ mesons. The equation is taken in the ladder and instantaneous approximations, reduced by contour and azimuthal integrations to a one-dimensional eigenvalue problem, and solved for the binding energy; a bound state is declared when the eigenvalue reaches unity. Short-distance physics enters through a monopole form factor $F(k^2)=(\\Lambda^2-m^2)/(\\Lambda^2-k^2)$ at each vertex, with $\\Lambda=m+\\alpha\\Lambda_{\\rm QCD}$ and $\\Lambda_{\\rm QCD}=220$ MeV, so the single free parameter $\\alpha$ controls how much short-range attraction survives.","core_discovery":"The paper's central claim is that one-boson-exchange dynamics, treated covariantly through the Bethe-Salpeter equation, is strong enough in certain channels to bind pairs of charmed or bottom vector mesons. Specifically, the isoscalar $D^*\\bar D^*$ and $B^*\\bar B^*$ systems admit $S$-wave bound solutions for $J^{PC}=0^{++}$, $1^{+-}$, and $2^{++}$ once the monopole cutoff is large enough, while the corresponding isovector systems never bind within $\\alpha\\in[0.5,10]$. For the identical-boson $D^*D^*$ and $\\bar B^*\\bar B^*$ systems, Bose symmetry leaves the allowed channels $I(J^P)=0(1^+)$, $1(0^+)$, and $1(2^+)$; the first two bind at moderate cutoffs, whereas $1(0^+)$ requires a substantially larger cutoff and is flagged as less reliable. Bottom analogues bind with smaller cutoffs for the same binding energy and are therefore more favorable candidates. Every reported bound state is a solution of a one-dimensional integral equation with eigenvalue unity for binding energies in the range 1 to 50 MeV.","pith_inferences":["Beyond the paper, the isoscalar $1^{+-}$ $D^*\\bar D^*$ state would complete a heavy-quark spin-symmetry multiplet around $X(3872)$; discovering it would support the molecular picture, while not discovering it would mainly constrain the allowed cutoff range.","A sharper test of the model would be to compute near-threshold scattering lengths and line shapes, which are more sensitive to the interaction than the existence of a pole.","A natural extension not pursued here is the inclusion of coupled channels such as $D\\bar D$ and $D\\bar D^*$ alongside $D^*\\bar D^*$; the fragile $1(0^+)$ doubly heavy state may disappear or be stabilized once those channels are added.","Because bottom states bind at smaller $\\alpha$, future searches in the bottom sector are the most cost-effective way to confirm or rule out the predicted molecular spectrum."],"forward_implications":["Isoscalar hidden-heavy channels with $J^{PC}=0^{++}$, $1^{+-}$, and $2^{++}$ are predicted to host $S$-wave molecular states below the $D^*\\bar D^*$ and $B^*\\bar B^*$ thresholds.","No isovector $D^*\\bar D^*$ or $B^*\\bar B^*$ bound states are expected within $\\alpha\\le 10$, so isovector candidates, if confirmed, would require dynamics beyond this one-boson-exchange kernel.","In the doubly heavy sector, the $I(J^P)=0(1^+)$ channel is the most favorable molecular configuration, followed by $1(2^+)$, while the $1(0^+)$ state is the least trustworthy because it appears only at large cutoff.","Bottom-sector analogues require smaller cutoffs than charmed counterparts for the same binding energy, making $B^*\\bar B^*$ and $\\bar B^*\\bar B^*$ channels the most favorable experimental targets.","Binding energies rise monotonically with $\\alpha$, so the model's quantitative masses are regulator dependent; searches should target near-threshold states rather than a single predicted mass."],"supporting_citations":[{"why":"Supplies the effective Lagrangians and coupling constants used to build the one-boson-exchange vertices for charmed vector mesons.","marker":"[61]"},{"why":"Establishes the Bethe-Salpeter method for $D\\bar D^*$/$B\\bar B^*$ systems that this paper extends to vector-meson pairs.","marker":"[59]"},{"why":"Applies the same BS formalism to heavy baryonium and dibaryon systems, providing the methodological template.","marker":"[60]"},{"why":"Gives the heavy-quark spin-symmetry prediction of $D^*\\bar D^*$ and $B^*\\bar B^*$ partners of $X(3872)$, the baseline this calculation is compared with.","marker":"[54]"},{"why":"A meson-exchange model that obtained isoscalar $D^*\\bar D^*$ bound states, providing the channel-by-channel comparison.","marker":"[62]"},{"why":"A recent one-boson-exchange analysis of short-range dynamics in the same systems, used to benchmark the present results.","marker":"[46]"},{"why":"A quasipotential Bethe-Salpeter study of $D^{(*)}D^{(*)}$ and $B^{(*)}B^{(*)}$ interactions, the closest alternative calculation of the same channels.","marker":"[42]"},{"why":"A complex-scaling calculation of bound and resonant $D^{(*)}D^{(*)}$ and $D^{(*)}\\bar D^{(*)}$ states, used to compare isovector and doubly heavy findings.","marker":"[17]"},{"why":"A lattice QCD study of isovector $DD$, $DD^*$, and $D^*D^*$ scattering that found no pole, used as a check on the isovector predictions.","marker":"[67]"},{"why":"An effective field theory coupled-channel study identifying $I(J^P)=0(1^+)$ as the most favorable doubly heavy configuration, supporting the same conclusion here.","marker":"[58]"}],"fun_headline_variants":["Bottom vector meson pairs bind more easily than charmed","Bethe-Salpeter predicts molecular bound states of heavy meson pairs","Charmed and bottom meson pairs form S-wave bound states","Heavy meson molecules bind most easily for bottom pairs","Hidden-heavy vector meson pairs bind in isoscalar channels"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results stand on the assumption that tuning the single cutoff parameter $\\alpha$ over a wide range can stand in for the true short-distance strong interaction; if that interaction differs from the one-boson-exchange form, the predicted bound states could disappear.","fun_headline_variants_meta":{"raw":{"variants":["Bottom vector meson pairs bind more easily than charmed","Bethe-Salpeter predicts molecular bound states of heavy meson pairs","Charmed and bottom meson pairs form S-wave bound states","Heavy meson molecules bind most easily for bottom pairs","Hidden-heavy vector meson pairs bind in isoscalar channels"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001387,"raw_usage":{"total_tokens":5645,"prompt_tokens":1003,"completion_tokens":4642,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":619,"completion_tokens_details":{"reasoning_tokens":4555}},"tokens_in":619,"tokens_out":4642,"duration_ms":29864,"temperature":1.0,"reasoning_tokens":4555,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:51:39.827797+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice QCD calculation of $S$-wave $D^*\\bar D^*$ scattering with $I=0$ and $J^{PC}=0^{++}$, $1^{+-}$, $2^{++}$ that finds no near-threshold pole would directly contradict the central prediction; likewise, a coupled-channel meson-exchange calculation that removes these bound states once $D\\bar D$ and $D\\bar D^*$ channels are included would show the one-boson-exchange-only result is an artifact of the truncated kernel.","supporting_citations":[{"cited_title":"Possible bound states of Heavy Baryonium and Heavy Dibaryon systems","cited_arxiv_id":"2409.03315","evidence_quote":"Supplies the effective Lagrangians and coupling constants used to build the one-boson-exchange vertices for charmed vector mesons."},{"cited_title":"Systematic Study of Coupled-Channel Dynamics in Doubly Heavy Hadronic Molecules","cited_arxiv_id":"2606.01177","evidence_quote":"Establishes the Bethe-Salpeter method for $D\\bar D^*$/$B\\bar B^*$ systems that this paper extends to vector-meson pairs."},{"cited_title":"Radiative decays of the neutral $Z_c(3900)$ and $Z_c(4020)$","cited_arxiv_id":"2210.06783","evidence_quote":"Gives the heavy-quark spin-symmetry prediction of $D^*\\bar D^*$ and $B^*\\bar B^*$ partners of $X(3872)$, the baseline this calculation is compared with."},{"cited_title":"Bound States of the Heavy Flavor Vector Mesons and Y(4008) and $Z^{+}_1(4050)$","cited_arxiv_id":"0905.1188","evidence_quote":"A meson-exchange model that obtained isoscalar $D^*\\bar D^*$ bound states, providing the channel-by-channel comparison."},{"cited_title":"X(3872) as a molecular $D\\bar{D}^*$ state in the Bethe-Salpeter equation approach","cited_arxiv_id":"1710.07424","evidence_quote":"A complex-scaling calculation of bound and resonant $D^{(*)}D^{(*)}$ and $D^{(*)}\\bar D^{(*)}$ states, used to compare isovector and doubly heavy findings."},{"cited_title":"Zhang, M.-Z","cited_arxiv_id":null,"evidence_quote":"An effective field theory coupled-channel study identifying $I(J^P)=0(1^+)$ as the most favorable doubly heavy configuration, supporting the same conclusion here."}],"review_version":2}