{"id":"13ded92c-eb8b-4cf1-b85f-770691402a34","arxiv_id":"2501.15670","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"New lattice QCD measurements of hybrid static potentials, static-light thresholds, and gluelump masses on N_f=3+1 ensembles with pions near 420 MeV, using Laplace trial states.","lead":"Lattice QCD simulations with nearly physical pions are used to measure the energy of a quark-antiquark pair connected by an excited gluon field. The resulting hybrid static potentials and gluelump masses are inputs for models of exotic mesons such as the XYZ states.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The hybrid potentials rely on a truncated Laplacian eigenbasis (N_vec=100-200) with no convergence test; if these modes are insufficient, every extracted hybrid potential and its R->0 match to gluelumps is biased.","rationale":"The reader's weakest_assumption, the completeness of the truncated Laplacian eigenbasis, is also the most load-bearing premise I find. The abstract's central claim is an actual measurement; if the operator basis is too small, the variational estimates of the hybrid potentials are systematically shifted, and the claimed potentials, along with their apparent agreement with the R->0 gluelump limits, could be artifacts of the truncation rather than properties of QCD. The paper provides no convergence evidence, and because the method is nonstandard (Laplace trial states replacing Wilson loops), it cannot fall back on prior validation in the same context. I considered the self-energy mismatch of Eq. (3) as an alternative concern; however, the paper acknowledges it explicitly and it affects only the common additive normalization when overlaying gluelumps, not the potential shape or the gluelump mass splittings. The eigenbasis truncation is unacknowledged and threatens all hybrid potential data. The proposed check, doubling N_vec on A1, would settle it. Since this is an addressable systematic rather than a demonstrated error, the CONDITIONAL verdict stands.","tokens_in":8532,"tokens_out":6305,"duration_ms":58784,"concrete_test":"Recompute the hybrid static potentials in the Sigma_g and Pi_u channels on ensemble A1 (N_vec=100) with N_vec=200 Laplacian eigenvectors, keeping all other parameters fixed, and compare the effective masses at the same R/r0 values; if the shifts exceed the statistical errors, the truncated basis is not converged and the presented potentials cannot be considered final. If resources allow, repeat on A1h (200 to 300).","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section 2 the Laplace trial-state correlators C_Sigma and C_Pi sum over i,j running over the N_vec lowest eigenvectors of the 3D covariant Laplacian (Table 1: 100 on A1/A2, 200 on A1h). The seven Gaussian profiles and subsequent GEVP only reweight these fixed modes; they cannot add support from higher eigenvectors. The central result, the hybrid static potentials on N_f=3+1 ensembles, therefore inherits a variational bias if the truncated product space v_i(x) v_j^dagger(y) is not saturated. No N_vec-dependence study and no cross-check against standard Wilson-loop hybrid potentials (e.g., Ref. [6]) is presented. This is the weakest link because the gluelump masses themselves rest on an independent loop-shape basis with internal E/T2 consistency; the hybrid potentials have no comparable safeguard. The paper's own caveat in Section 4 that the gluelump self-energy 'does not completely remove' the mismatch is a second, acknowledged systematic, but it affects only the overlay, not the potential values.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings paper presents lattice QCD computations of hybrid static quark-antiquark potentials using a novel \"Laplace trial state\" construction, in which gluonic excitations are introduced through covariant derivatives of the lowest Laplacian eigenvectors, with Gaussian profiles and a GEVP used to improve overlaps. The same paper computes gluelump masses on an extended basis of 35 spatial loop shapes for ensembles with N_f=0, 2, and 3+1 dynamical quarks (pion masses down to 406 MeV). The hybrid potentials are compared with static-light S- and P-wave thresholds, and the lowest gluelump masses are overlaid at R->0 as the expected limit of the hybrid potentials. Statistical errors are estimated with the Gamma method via pyerrors.","tokens_in":8739,"tokens_out":7076,"duration_ms":63602,"significance":"If the results are correct, they provide new lattice-QCD inputs for Born-Oppenheimer effective models of hybrid mesons, including excited hybrid channels and string-breaking thresholds for N_f=3+1 QCD with physical-ish pions. The Laplace trial-state method offers practical advantages: off-axis separations can be computed easily, and perambulator-based contractions open the way to studies of string breaking. The internal cross-checks are a genuine strength: the E and T2 gluelump representations agree for the J=2 channel, and the A1/A2 volumes give consistent gluelump masses. No machine-checked proofs or released code accompany the paper, so reproducibility rests on the textual description. The central new systematic—truncation of the Laplacian eigenbasis—is not tested, which limits the strength of the claims at this stage.","major_comments":[{"comment":"The hybrid potentials are constructed from the lowest N_vec Laplacian eigenvectors (N_vec=100 on A1/A2 and 200 on A1h, Table 1), and the correlation functions sum only over these modes and their covariant derivatives. The paper provides no convergence test with respect to N_vec and no comparison of the resulting potentials with standard Wilson-loop hybrid potentials (e.g., Ref. [6]). Since the trial-state space is confined to these modes, an incomplete basis could bias every extracted V(R) and the R->0 comparison with gluelumps. The authors should add a controlled check, for example by repeating one ensemble with a larger N_vec, or by comparing at least the ground-state Sigma_g^+ potential with a Wilson-loop calculation on the same ensembles, and discuss the expected size of the truncation effect.","section":"Section 2, Eqs. for C_Sigma and C_Pi, Table 1"},{"comment":"The overlay of gluelump masses on the hybrid potentials at R->0 relies on the approximate self-energy relation m_self = (N^2/(N^2-1)) V_self, and the text states that subtracting twice the static-light S-wave mass \"does not completely remove\" the self-energy contribution. The figure shows no systematic uncertainty for this residual self-energy difference. Because the claim that gluelump masses are the R->0 limit of the hybrid potentials is central, the authors should either quantify the residual mismatch (for instance by estimating V_self on the same ensembles through a Wilson-loop computation) or explicitly state that the R->0 correspondence is qualitative and carries an unquantified self-energy uncertainty.","section":"Section 4, Figure 2 and Eq. (3)"}],"minor_comments":[{"comment":"The symbol \\bar{\\tau}_{ij} is introduced in the text but is not used in the displayed equations, where only \\tau_{ij} and \\tau_{ji} appear; please clarify the notation for the conjugate perambulator.","section":"Section 2, first paragraph"},{"comment":"The labels for excited channels such as \\Sigma_g' and \\Pi_u' are not defined in the caption or the text; please add a brief explanation, for example that these are the first excited states obtained from the GEVP.","section":"Figure 2 and caption"},{"comment":"The text refers to \"A1 heavy\" while the table uses \"A1h\"; to avoid confusion, a single consistent ensemble name should be used throughout.","section":"Section 3, Table 1 and Section 4"},{"comment":"The statement that \"the individual energies consistently increase with decreasing pion mass\" is ambiguous, since the plotted quantities are mass differences with respect to the T_1^{+-} state; please specify whether this refers to the absolute gluelump masses or to the relative splittings.","section":"Section 4, first paragraph"}],"recommendation":"major_revision","confidential_remarks":"This is a LATTICE2024 proceedings paper, so some brevity is expected, but the absent N_vec convergence test is a load-bearing issue for the hybrid-potential claim. The authors could address it in a revision by adding a convergence study or a direct Wilson-loop cross-check; this is feasible within the scope of a proceedings update. The self-energy caveat is acknowledged but should be quantified to support the R->0 correspondence. I see no circularity or misuse of prior results; the main concerns are technical validation and presentation clarity."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a credible proceedings report that delivers the first hybrid static potentials and gluelump masses on N_f=3+1 ensembles with pions around 410 MeV. The methods are not new—the Laplace trial states come from the authors' earlier work, and the gluelump loop basis from ref [13]—but the combination and the data are new, and the internal consistency checks are genuinely reassuring: A1 and A2 agree across volumes, and E and T2 representations agree for the J=2 gluelump channel.\n\nWhat the paper does well: it is transparent about what is measured and what is not. The gamma-method errors are standard, the ensembles are well described, and the authors flag the self-energy caveat in the gluelump/hybrid overlay instead of hiding it. The observation that string-breaking thresholds are accessible on A2 is a nice, concrete outlook.\n\nThe soft spots, in order of importance. First, there is no convergence test for the Laplacian eigenbasis truncation. The trial states use only the lowest 100–200 eigenvectors per time slice; if those modes do not saturate the gluonic excitations, the hybrid potentials inherit a variational bias. The paper shows no N_vec-dependence and no cross-check against the Wilson-loop hybrid potentials from Capitani et al. This is the load-bearing technical premise, and it is unchecked. Second, all results are at a single lattice spacing; no continuum extrapolation is attempted, which limits the quantitative value. Third, the gluelump masses are only plotted as relative splittings and the self-energy mismatch in Eq. (3) is acknowledged but not quantified. These are addressable limitations, not fatal flaws, and they are normal for a proceedings contribution.\n\nOne thing I want to push back on: the stress-test note worries that the truncated basis is a 'weakest link' with no safeguard. On reading the paper, the internal E/T2 agreement and the volume agreement do provide some safeguard that the method is not wildly off, but they do not protect against a common bias from the eigenbasis truncation. So the concern stands, but its severity is moderate, not critical. The paper is honest about what is preliminary.\n\nWho is this for? Lattice practitioners and hadron phenomenologists using Born-Oppenheimer models for hybrid and exotic mesons. It gives them a new data point at nearly physical pion mass, with the caveat that it is not yet continuum-extrapolated.\n\nRecommendation: if this crossed my desk as a journal submission, I would send it to review. The missing N_vec study and continuum extrapolation would be the referee's main asks, but the paper is a legitimate and useful contribution. I would not desk-reject it.","headline":"New N_f=3+1 hybrid potential and gluelump data with honest caveats, but the truncated Laplace basis needs a convergence check.","tokens_in":9320,"tokens_out":2521,"would_cite":true,"duration_ms":21944,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81T25","81V05"],"pacs":["12.38.Gc"],"model":"deepseek-v4-flash","headline":"This paper reports new lattice QCD measurements of hybrid static potentials and gluelump masses on 3+1-flavor ensembles with 420 MeV pions, including string-breaking thresholds.","keywords":["hybrid static potentials","gluelumps","lattice QCD","Laplace trial states","static-light thresholds","string breaking","Born-Oppenheimer effective theory","N_f=3+1 ensembles"],"falsifier":"Compute the same hybrid potentials and gluelump masses on the A2 ensemble with twice the number of kept Laplace eigenvectors, or with standard Wilson-loop hybrid operators, and look for differences beyond statistical errors; agreement would validate the truncated basis, and disagreement would show it biases the results.","tokens_in":8318,"feed_emoji":"⚛️","tokens_out":9051,"duration_ms":72507,"temperature":0.7,"pith_summary":"The paper works out, on dynamical QCD ensembles with three light quark flavors plus charm and pion masses down to about 420 MeV, the static potentials that bind a heavy quark and antiquark into ordinary and exotic hybrid mesons. Its central result is a spectrum of hybrid static potentials in the $\\Sigma_g$, $\\Sigma_u$, $\\Pi_g$, and $\\Pi_u$ channels, computed with operators built from eigenmodes of the lattice Laplace operator instead of Wilson loops, together with the thresholds at which the confining string breaks into static-light mesons. It also delivers gluelump masses — the energy levels of gluonic excitations around a static color-octet source — which are the $R\\to0$ limit of these hybrid potentials and serve as input to Born-Oppenheimer models of hybrid mesons. On the largest-volume ensemble the authors observe the onset of string breaking for $\\Sigma_g$ and $\\Pi_u$, and they show that the $E$ and $T_2$ gluelump channels, which must coincide in the continuum, agree on each ensemble. A sympathetic reader would take this as evidence that Laplace-trial-state operators can replace Wilson loops for hybrid spectroscopy on realistic ensembles and give physics consistent with the established gluelump picture.","feed_headline":"New lattice QCD data pin down hybrid-meson forces and gluelump masses","feed_subtitle":"Quark-gluon hybrid potentials and their string-breaking thresholds measured with 420 MeV pions.","key_machinery":"The load-bearing object is the Laplace trial state: a correlator in which the spatial Wilson line of a Wilson loop is replaced by products of Laplace eigenvectors $v_i(\\vec x)$ and static perambulators $\\tau_{ij}=v_i^\\dagger U_t v_j$, with gluonic excitations inserted as symmetric covariant derivatives $\\nabla_{\\vec k} v(\\vec x)=\\frac12[U_k(\\vec x)v(\\vec x+\\hat k)-U_k^\\dagger(\\vec x-\\hat k)v(\\vec x-\\hat k)]$ of the eigenvectors. Gaussian profile weights $\\rho_k(\\lambda_i)=e^{-\\lambda_i^2/4\\sigma_k^2}$ assign different weights to eigenmodes, and a generalized eigenvalue problem selects the optimal profiles. This machinery lets one compute correlations at arbitrary on- and off-axis separations without expensive loop re-evaluations. For gluelumps the paper uses a complementary tool: a 35-shape basis of spatial Wilson loops, projected to cubic-group irreps $A_1,A_2,E,T_1,T_2$ with parity and charge conjugation, connected by adjoint static lines, with the adjoint self-energy related to the fundamental one by the Casimir ratio $C_A/C_F=N^2/(N^2-1)$.","core_discovery":"The paper claims that hybrid static potentials extracted from Laplace trial states built from covariant derivatives of lattice Laplace eigenmodes reproduce the expected gluelump multiplets in the $R\\to0$ limit and, at pion masses near 420 MeV, show the string-breaking thresholds where $\\Sigma_g$ and $\\Pi_u$ hybrids meet static-light meson states. On the $N_f=3+1$ ensemble A2 ($m_\\pi\\simeq 409$ MeV, $48^3$ spatial volume) this threshold lies within the lattice, so the potentials are measured across the crossing region. The gluelump spectrum is obtained from a large basis of spatial loop shapes projected onto irreducible representations of the cubic group, with masses given relative to the lowest $T_1^{+-}$ state; the $E^{PC}$ and $T_2^{PC}$ channels are consistent for each $PC$, which is the expected continuum degeneracy for $J=2$. The paper therefore presents these results as new lattice QCD inputs for effective-theory descriptions of hybrid and exotic mesons.","pith_inferences":["Because the Laplace-trial-state machinery computes correlations at arbitrary separations from precomputed perambulators, the same construction should transfer to tetraquark and multiquark static potentials more cheaply than Wilson-loop operators; the paper lists this as a plan, so the transfer is an expectation rather than a demonstrated result.","The trend shown for gluelump masses — larger masses at smaller pion mass — suggests that Born-Oppenheimer hybrid models calibrated on quenched or heavy-pion gluelumps may misplace hybrid levels at physical pion mass; this extrapolation is ours, not the paper's.","A decisive test of the operator replacement would be to compare these Laplace-trial-state potentials with Wilson-loop hybrid potentials on the same A2 ensemble; the paper does not show such a comparison."],"forward_implications":["On ensemble A2 the $\\Sigma_g$ and $\\Pi_u$ hybrid potentials cross the static-light thresholds within the measured range, so hybrid string breaking can be studied directly on that ensemble.","The $E^{PC}$ and $T_2^{PC}$ gluelump channels agree within errors on each ensemble, confirming that the loop-shape basis keeps cutoff effects small enough for the continuum $J=2$ degeneracy to be visible.","Gluelump masses rise systematically as the pion mass is lowered from the quenched value through 2.2 GeV and 788 MeV to 406 MeV, so dynamical fermion effects matter for these masses.","Because Laplace-trial-state correlators give easy access to off-axis separations, the hybrid spectrum can be resolved at many more $R$ values than with standard Wilson-loop operators."],"supporting_citations":[{"why":"Introduces the Laplace trial-state correlator, the eigenmode-based replacement for the spatial Wilson line that the hybrid potential extraction uses.","marker":"[14]"},{"why":"Supplies the extended basis of spatial loop shapes and projection methods used for the gluelump operators.","marker":"[13]"},{"why":"Establishes the classical method for computing gluelump masses from adjoint color sources, which this work refines.","marker":"[9]"},{"why":"Provides the SU(3) gluelump mass ordering and splittings that the present measurements extend to dynamical ensembles.","marker":"[11]"},{"why":"Sets the scale for the N_f=3+1 ensembles, giving the lattice spacings used to convert results to physical units.","marker":"[26]"},{"why":"Gives the prior high-precision Wilson-loop hybrid static potentials in SU(3) that the Laplace-trial-state results are compared with.","marker":"[6]"},{"why":"Computes static-light meson spectroscopy with optimal profiles, providing the S- and P-wave thresholds used for string breaking.","marker":"[42]"},{"why":"Supplies the Born-Oppenheimer effective theory for hybrids, the application for which these gluelump masses are essential input.","marker":"[5]"}],"fun_headline_variants":["Lattice QCD measures hybrid meson forces and gluelump masses","Hybrid static potentials and gluelump masses from lattice QCD","Exotic meson clues: hybrid potentials and gluelumps on the lattice","String-breaking thresholds for hybrid mesons via lattice QCD","Lattice QCD yields hybrid-meson potentials and gluelump spectrum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole extraction rests on the assumption that keeping only the lowest 100 or 200 Laplace-eigenvector modes on each time slice gives a complete enough basis for the hybrid and gluelump states, a point the paper does not test by varying that number.","fun_headline_variants_meta":{"raw":{"variants":["Lattice QCD measures hybrid meson forces and gluelump masses","Hybrid static potentials and gluelump masses from lattice QCD","Exotic meson clues: hybrid potentials and gluelumps on the lattice","String-breaking thresholds for hybrid mesons via lattice QCD","Lattice QCD yields hybrid-meson potentials and gluelump spectrum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00022,"raw_usage":{"total_tokens":1404,"prompt_tokens":858,"completion_tokens":546,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":474,"completion_tokens_details":{"reasoning_tokens":453}},"tokens_in":474,"tokens_out":546,"duration_ms":4843,"temperature":1.0,"reasoning_tokens":453,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:04:28.316850+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same hybrid potentials and gluelump masses on the A2 ensemble with twice the number of kept Laplace eigenvectors, or with standard Wilson-loop hybrid operators, and look for differences beyond statistical errors; agreement would validate the truncated basis, and disagreement would show it biases the results.","supporting_citations":[{"cited_title":"Constructing static quark-anti-quark creation operators from Laplacian eigenmodes","cited_arxiv_id":"2212.08485","evidence_quote":"Introduces the Laplace trial-state correlator, the eigenmode-based replacement for the spatial Wilson line that the hybrid potential extraction uses."},{"cited_title":"Static-light meson spectroscopy with optimal distillation profiles","cited_arxiv_id":"2501.12863","evidence_quote":"Computes static-light meson spectroscopy with optimal profiles, providing the S- and P-wave thresholds used for string breaking."}],"review_version":1}