{"id":"79f07ed2-f802-4425-9cbd-4a40b258d371","arxiv_id":"2509.10940","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Ground-state masses of fully heavy tetraquarks are predicted with the hyperspherical hyperradial approximation, yielding e.g. 5.86 GeV for the lightest all-charm state.","lead":"This paper calculates the masses of heavy tetraquarks, particles made of four quarks or antiquarks, using a quantum-mechanics approximation called the hyperspherical method. It predicts specific masses for three tetraquark systems, which can be compared with particle-collider data to help identify these exotic states.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The confinement potential in Eq. (6) includes a constant B=-0.8 GeV that is absent from the solved hyperradial equation (17); until this is explained, Table I does not follow from the stated Hamiltonian.","rationale":"The reader's K=0 concern is reasonable and is supported by the paper's own admission that the real tetraquark wave function may have a small angular dependence. However, the K=0 approximation is at least partially tested by the numerical/variational agreement and by the quoted 0.03 GeV difference. The most load-bearing issue I find is not the K=0 truncation but the inconsistent treatment of the constant B in the confinement potential. Eq. (6) explicitly includes B=-0.8 GeV in ∆V^conf, yet Eq. (17) and Eq. (19) omit it. Since the paper's error estimate is only 0.1 GeV while B is a 0.8 GeV shift, this is a first-order consistency problem. The concrete test I propose would settle the issue: if the masses shift by the full B value, the central numerical claim is not reproducible from the stated Hamiltonian; if the authors clarify that B is already absorbed into the quark masses or the zero of energy, the concern disappears. I keep the verdict conditional rather than rejecting outright because the fix may be purely presentational, but the manuscript as written is not self-consistent. This is why my agreement with the reader's weakest assumption is only partial: the reader identified a legitimate approximation issue, but the B-term inconsistency is a more basic and more easily testable threat to the central claim.","tokens_in":12402,"tokens_out":24687,"duration_ms":309412,"concrete_test":"Redo the variational/numerical solution of Eq. (19) with E replaced by E+B (i.e., add the constant B from Eq. (6) to the averaged potential in Eq. (18)), keeping all other parameters fixed, and recompute the three (cc\\bar c\\bar c) masses in Table I. If the 0++ mass moves by about 0.8 GeV, the table is not a solution of the stated Hamiltonian; if it does not, the paper must state explicitly how B is removed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II defines ∆V^conf as Σ A_ij|r_i-r_j| + B and fixes B = -0.8 GeV. Section III constructs the averaged potential V(R) = <V^C + V^conf> (Eq. 18), but the resulting hyperradial equation (17) and its scaled form (19) contain only +a/R - bR - 6/(µR^2), with no constant term. A constant in the potential shifts every eigenvalue by B. If B is treated literally, the (cc\\bar c\\bar c) 0++ mass in Table I would shift downward by ~0.8 GeV, to about 5.06 GeV, far outside the claimed ±0.1 GeV error. This is an internal inconsistency of the model as written, independent of the K=0 approximation and of the size of the Breit corrections. The paper does not state that B is absorbed into the quark masses or otherwise cancelled.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes masses of fully heavy tetraquarks (cc\\bar c\\bar c), (bb\\bar b\\bar b), and (cc\\bar b\\bar b) in the ground states 0++, 1+-, 2++ using the hyperspherical approach. The four-body Schrödinger equation is reduced to a one-dimensional hyperradial equation by assuming K=0 and averaging over the nine-dimensional angular variables. The hyperradial equation is solved numerically and variationally, with good mutual agreement (E0=0.380 vs 0.382 GeV for (cc\\bar c\\bar c)). Hyperfine splitting is computed from one-gluon-exchange and confinement spin-spin terms, and relativistic and recoil corrections are added. The final masses are listed in Table I, e.g., 5.86, 6.02, and 6.35 GeV for (cc\\bar c\\bar c).","tokens_in":12694,"tokens_out":5578,"duration_ms":70721,"significance":"If the calculation is correct, it provides a simple analytical framework for fully heavy tetraquark masses and wave functions, with predictions that can be compared with LHCb, CMS, and ATLAS searches. A strength of the paper is that all model parameters are taken from earlier meson calculations rather than fitted to tetraquark data, so the tetraquark masses are genuine predictions. The numerical and variational solutions agree well, and the analytical wave function (26) allows transparent computation of hyperfine and relativistic corrections. The main weakness is that the K=0 hyperradial approximation is uncontrolled, and there is an apparent internal inconsistency in the treatment of the constant term B in the confinement potential.","major_comments":[{"comment":"The confinement potential in Eq. (6) contains a constant term B=-0.8 GeV, but the hyperradial Schrödinger equation (17) contains no such constant. Averaging a constant over angles gives the same constant, so if B is part of the Hamiltonian it shifts every eigenvalue by -0.8 GeV. With B included literally, the (cc\\bar c\\bar c) 0++ mass in Table I would be about 5.06 GeV, far outside the claimed ±0.1 GeV error. The manuscript does not state that B is absorbed into the quark masses or cancelled by another term. This is a load-bearing internal inconsistency that must be resolved.","section":"Section II, Eq. (6) and Section III, Eq. (17)"},{"comment":"The K=0 hyperradial approximation replaces the potential by its angular average and assumes the wave function depends only on R. No convergence check with K>0 components is provided. The paper itself notes a 0.03 GeV difference from the variational result, which it attributes to a possible small dependence of the wave function on angles. This approximation is central not only for the binding energy but also for the delta-function matrix elements used in Section IV for hyperfine splitting. Without a quantitative estimate of the angular dependence, the stated theoretical error of 0.1 GeV appears optimistic.","section":"Section III, Eqs. (14)-(18)"},{"comment":"The relativistic, recoil, and contact corrections are presented as results of analytical calculations, but no derivation or intermediate steps are given. The text states that these corrections are 'numerically large' and enter with negative sign, so they materially affect Table I. The reader cannot verify the formulas, and Eq. (42) contains an unusual coefficient (20337) that may indicate a typo. The authors should provide at least a derivation outline and numerical cross-checks of these matrix elements.","section":"Section V, Eqs. (42), (44), (46), (48)"},{"comment":"The effective constants a and b are given explicitly only for tetraquarks with identical quark masses, and the hyperfine coefficient κ is also written for the equal-mass case. However, Table I includes (cc\\bar b\\bar b), where quark masses differ. The paper does not state how the angular averages, the reduced masses, and the delta-function matrix elements are generalized to unequal masses. Without this, the (cc\\bar b\\bar b) entries in Table I are not supported by the equations presented.","section":"Section III, Eq. (20) and Section IV, Eq. (34)"}],"minor_comments":[{"comment":"Typo: 'Schroedinger' should be 'Schrödinger'.","section":"Abstract"},{"comment":"The notation '<δ(r12)>=<δ(ρ)>= ... =δ=0.0846257426 GeV^3' is confusing because δ denotes both the delta function and the numerical value. Use a separate symbol for the matrix element.","section":"Section III, Eq. (33)"},{"comment":"The spin wave functions are labeled χ^{11}_{00}, χ^{11}_{11}, χ^{11}_{22}; the second subscript should denote the total spin projection, not the total spin. The notation is inconsistent and should be clarified.","section":"Section IV, Eqs. (30)-(32)"},{"comment":"The coefficient 20337 inside the parenthesis appears suspicious; please verify the arithmetic. It would also help to state the units of p explicitly in Eq. (42).","section":"Section V, Eq. (42)"},{"comment":"The statement that a and b do not depend on μ is made only for the equal-mass case. For the unequal-mass tetraquark, the derivation should be spelled out.","section":"Section III, Eq. (20)"}],"recommendation":"major_revision","confidential_remarks":"The constant-B issue in Eqs. (6) and (17) is the most serious concern: if it is a literal inconsistency, Table I is not supported; if B is absorbed elsewhere, that must be stated explicitly. The paper also relies heavily on self-citations [23], [42], [43]; while this is not improper, the present manuscript should be self-contained in its derivations. The topic fits the journal's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: Martynenko et al. apply the hyperspherical method to fully heavy tetraquarks and get a set of ground-state masses for (cc\\bar c\\bar c), (bb\\bar b\\bar b), and (cc\\bar b\\bar b) in the 0++, 1+-, 2++ channels. The numbers are plausible and fall within the spread of existing quark-model predictions. What's new here is the combination of the hyperradial approximation with a variational wave function and the QCD-Breit corrections; the numerical and variational solutions agree nicely (0.380 vs 0.382 GeV). The parameters are taken from earlier meson fits rather than fit to tetraquark data, so the circularity burden is low. The paper is honest about the K=0 approximation and notes the 0.03 GeV difference with their previous variational result.\n\nThe soft spot that matters is the constant B = -0.8 GeV in the confinement potential. Equation (6) puts it in the Hamiltonian, but the hyperradial equation (17) and its scaled form (19) have no constant term. Averaging a constant over angles just gives the constant, so it should appear in the potential. The paper never says it is absorbed into the quark masses or redefined away. This is not a minor issue: a constant shift of -0.8 GeV changes every eigenvalue by that amount, so the (cc\\bar c\\bar c) 0++ mass in Table I would drop to about 5.06 GeV if B were included literally. That is far outside the claimed 0.1 GeV error. This is an internal inconsistency in the model as written, and the paper needs to clarify it.\n\nThe K=0 approximation is a second, milder concern; the paper acknowledges it, so it's not a hidden assumption. Section V's corrections are given without derivations, which makes the calculation hard to audit. The uncertainty estimate is a single global 0.1 GeV, not a propagation of parameter errors.\n\nWho should read this: people interested in quark-model predictions for fully heavy tetraquarks, especially in the context of LHCb/CMS/ATLAS searches. The method gives a useful independent check, but the B inconsistency must be fixed before the absolute masses are usable. I'd send it to peer review; a referee should insist on a clear explanation of the missing constant and a sensitivity analysis. If that gets resolved, the paper is a solid contribution.","headline":"A useful hyperspherical cross-check on fully heavy tetraquark masses, but the missing B = -0.8 GeV constant in the solved equation makes the absolute scale ambiguous.","tokens_in":13173,"tokens_out":10643,"would_cite":false,"duration_ms":116949,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.39.Ki","14.40.Pq"],"model":"deepseek-v4-flash","headline":"This paper claims that a one-dimensional hyperspherical equation, obtained by averaging over angles, predicts ground-state masses of fully heavy tetraquarks with about 0.1 GeV accuracy.","keywords":["tetraquarks","hyperspherical method","hyperradial approximation","heavy quarks","quark model","Breit Hamiltonian","hyperfine splitting","mass spectrum"],"falsifier":"Measure the mass of the 0++ fully charmed tetraquark; if it differs from 5.86 GeV by more than the stated ~0.1 GeV error, the angular-averaging step fails. Alternatively, a hyperspherical calculation that retains $K>0$ harmonics and shifts the binding energy by more than ~0.03 GeV would falsify the $K=0$ truncation.","tokens_in":12319,"feed_emoji":"⚛️","tokens_out":4109,"duration_ms":45952,"temperature":0.7,"texified_at":"2026-08-05T20:30:29.695244+00:00","pith_summary":"The paper tries to establish that the hyperradial approximation—reducing the four-quark problem to a single radial Schrödinger equation after angular averaging—can reproduce the masses of heavy tetraquarks such as (cc\\bar c\\bar c), (bb\\bar b\\bar b), and (cc\\bar b\\bar c). It predicts ground-state masses of 5.86, 6.02, and 6.35 GeV for the all-charm tetraquark in the 0++, 1+-, and 2++ channels, with similar predictions for the all-beauty and mixed states. The authors argue that the resulting analytic wave function is simple enough to compute hyperfine splitting and relativistic corrections, and that the approach gives results comparable to other quark-model calculations while offering a transparent physical picture.","texify_model":"deepseek-v4-flash","texify_usage":{"total_tokens":7153,"prompt_tokens":800,"completion_tokens":6353,"prompt_tokens_details":{"cached_tokens":0},"prompt_cache_hit_tokens":0,"prompt_cache_miss_tokens":800,"completion_tokens_details":{"reasoning_tokens":5584}},"feed_headline":"One radial equation predicts tetraquark masses to 0.1 GeV","feed_subtitle":"Angular averaging reduces four-quark binding to one dimension, yielding ground-state masses for all-charm, all-beauty, and mixed tetraquarks","key_machinery":"The central object is the hyperradius $R$, defined by $R^2 = \\rho^2 + \\lambda^2 + \\sigma^2$ in Jacobi coordinates, and the hyperspherical angular momentum $K$. In the hyperradial approximation one sets $K=0$, so the wave function depends only on $R$ and the potential is replaced by its angular average, giving $\\langle 1/\\rho \\rangle = \\frac{35}{16R}$ and $\\langle \\rho \\rangle = \\frac{35R}{64}$. This leads to a one-dimensional equation with effective Coulomb and linear constants $a = \\frac{175}{12\\sqrt{2}} \\alpha_s \\sqrt{m}$ and $b = \\frac{175}{32\\sqrt{2}} A \\sqrt{m}$. The solution is obtained with a variational wave function of the form $\\chi(x) = \\left[ \\frac{2}{9} q p^{9/q} / \\Gamma(9/q) \\right]^{1/2} x^4 e^{-p x^q}$, whose Airy-function asymptotics match the confining potential. This wave function is then used to compute spin-spin hyper","core_discovery":"Within the quark model, the paper reduces the four-body Schrödinger equation for heavy tetraquarks to a one-dimensional equation for a wave function that depends only on the hyperradius $R$, after averaging the potential over the angles of a nine-dimensional hypersphere. The equation contains an effective Coulomb attraction $a/R$ and a linear confining term $bR$, with $a$ and $b$ expressed in terms of quark masses and the strong coupling constant. Solving this equation numerically and with a trial function of the form $\\chi(x) \\propto x^4 \\exp(-p x^q)$, the authors find a nonrelativistic binding energy of 0.380 GeV for (cc\\bar c\\bar c), matching the numerical solution. Adding relativistic kinetic corrections, rec","pith_inferences":["Inference: If the hyperradial approximation is accurate, the same reduction could be applied to excited tetraquark states by retaining K>0 hyperspherical harmonics; the centrifugal term -6/(μR^2) already present would then shift and split the spectrum in a predictable way.","Inference: The predicted 0++ all-charm mass near 5.86 GeV lies below the J/ψ pair threshold, suggesting a narrow state, whereas the 2++ at 6.35 GeV lies above it; this threshold crossing is a testable line-shape prediction not explicitly stated in the paper.","Inference: The value of the tetraquark wave function at zero separation (Ψ_T(0)=0.09 GeV^{9/2}) directly controls production rates; a future measurement of exotic tetraquark production in Higgs decays would provide a quantitative test of both the wave function and the hyperradial approach."],"forward_implications":["If the hyperradial predictions are correct, the 0++ all-charm tetraquark should appear near 5.86 GeV, with a hyperfine splitting of about 0.16 GeV between the 0++ and 1+- states.","The simple analytic wave function (26) can be used to estimate production and decay matrix elements, such as tetraquark production in rare Higgs decays.","The ordering of states (0++ below 1+- below 2++) is a direct consequence of the spin-spin Hamiltonian, providing a clear experimental signature.","The 0.03 GeV difference between the hyperradial and variational binding energies suggests that including small angular dependence in the wave function would shift the masses by at most a few tens of MeV.","The method, with its fixed quark-model parameters, yields masses that fall within the range of other quark-model predictions, supporting the search for these states at colliders."],"fun_headline_variants":["Hyperspherical trick predicts tetraquark masses","Four-quark puzzle solved with one radial equation","Tetraquark spectrum from a single wave function","All-heavy tetraquarks bound by one radial equation","Hyperspherical averaging yields precise tetraquark masses"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The prediction rests on assuming the four-quark wave function is independent of the angular variables in nine-dimensional space, so the real interaction is replaced by its angular average; the paper itself notes a 0.03 GeV binding-energy difference with its variational calculation, which may reflect small angular dependence.","fun_headline_variants_meta":{"raw":{"variants":["Hyperspherical trick predicts tetraquark masses","Four-quark puzzle solved with one radial equation","Tetraquark spectrum from a single wave function","All-heavy tetraquarks bound by one radial equation","Hyperspherical averaging yields precise tetraquark masses"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000475,"raw_usage":{"total_tokens":2132,"prompt_tokens":621,"completion_tokens":1511,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":365,"completion_tokens_details":{"reasoning_tokens":1435}},"tokens_in":365,"tokens_out":1511,"duration_ms":13597,"temperature":1.0,"reasoning_tokens":1435,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T17:21:55.788609+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the mass of the 0++ fully charmed tetraquark; if it differs from 5.86 GeV by more than the stated ~0.1 GeV error, the angular-averaging step fails. Alternatively, a hyperspherical calculation that retains $K>0$ harmonics and shifts the binding energy by more than ~0.03 GeV would falsify the $K=0$ truncation.","supporting_citations":[],"review_version":1}