Recognition: unknown
Hybrid hadrons at rest and on the light front
Pith reviewed 2026-05-07 11:43 UTC · model grok-4.3
The pith
Heavy hybrid hadrons are described by treating the gluon as a massive quasiparticle in a Born-Oppenheimer framework, yielding light-front wave functions and gluon PDFs for ccg and qqqg systems.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We present a unified description of heavy hybrid hadrons based on a constituent-gluon picture embedded in the Born-Oppenheimer (BO) framework. In this approach, the gluonic excitation is treated as a dynamical quasiparticle with a mass generated by instanton-induced interactions. We propose a simple variational derivation of the BO potentials. The main focus of the paper is the derivation of light-front wave functions for hybrid systems, specifically for the ccg and qqqg cases. We employ both variational methods and numerical solutions of the Schrödinger equation in momentum representation. Using the resulting wave functions, we compute the gluon PDFs for these systems.
What carries the argument
The Born-Oppenheimer potentials for heavy quark-gluon systems, derived variationally with a massive constituent gluon generated by instanton interactions, which are then used to construct light-front wave functions by solving the Schrödinger equation in momentum space.
If this is right
- Light-front wave functions for the ccg and qqqg hybrid systems follow from solving the momentum-space Schrödinger equation on the derived Born-Oppenheimer potentials.
- Gluon parton distribution functions are obtained directly by integrating the light-front wave functions over the transverse momenta.
- The same framework describes the hybrids both in their rest frame and in light-front kinematics without additional assumptions.
- Variational methods provide a simple approximation to the full numerical solution for the wave functions.
Where Pith is reading between the lines
- The derived PDFs could be inserted into models of high-energy scattering to predict cross sections involving hybrid states.
- The approach suggests a route to compute additional light-front observables such as form factors once the wave functions are in hand.
- If the instanton-generated mass proves robust, the same quasiparticle treatment might apply to other gluonic excitations beyond the lowest hybrid states.
Load-bearing premise
The gluonic excitation can be treated as a dynamical quasiparticle whose mass is generated by instanton-induced interactions, and that a variational approach to the Born-Oppenheimer potentials is adequate for deriving reliable light-front wave functions.
What would settle it
A significant mismatch between the gluon PDFs computed from these light-front wave functions and PDFs extracted from lattice QCD calculations or from experimental data on hybrid hadron production would falsify the central claim.
Figures
read the original abstract
We present a unified description of heavy hybrid hadrons based on a constituent-gluon picture embedded in the Born-Oppenheimer (BO) framework. In this approach, the gluonic excitation is treated as a dynamical quasiparticle with a mass generated by instanton-induced interactions. We propose a simple variational derivation of the BO potentials. The main focus of the paper is the derivation of light-front wave functions for hybrid systems, specifically for the $ccg$ and $qqqg$ cases. We employ both variational methods and numerical solutions of the Schr\"odinger equation in momentum representation. Using the resulting wave functions, we compute the gluon PDFs for these systems.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims to present a unified description of heavy hybrid hadrons based on a constituent-gluon picture embedded in the Born-Oppenheimer (BO) framework, treating the gluonic excitation as a dynamical quasiparticle with a mass generated by instanton-induced interactions. It proposes a simple variational derivation of the BO potentials. The main focus is the derivation of light-front wave functions for hybrid systems, specifically for the ccg and qqqg cases, using both variational methods and numerical solutions of the Schrödinger equation in momentum representation, followed by computation of the gluon PDFs for these systems.
Significance. If validated, this framework could bridge non-relativistic effective models with light-front dynamics for exotic hybrid states, offering a practical route to gluon PDFs that are otherwise challenging to compute. The combination of variational and numerical techniques in momentum space is a methodological strength, and extending the BO approach to light-front wave functions addresses a relevant gap in hybrid hadron phenomenology.
major comments (2)
- Abstract: the description of variational and numerical methods supplies no error estimates, comparisons to data or lattice results, or validation of the quasiparticle mass assumption, leaving the central claim of computing usable gluon PDFs without demonstrated quantitative support.
- The section on the constituent-gluon picture and variational derivation of BO potentials: the gluon is modeled as a massive quasiparticle from instanton-induced interactions without explicit separation of this free parameter from predictions or sensitivity tests, which risks circularity when the resulting wave functions are used to extract PDFs.
Simulated Author's Rebuttal
We thank the referee for the positive evaluation of the significance of our work and for the detailed comments. We address the major concerns point by point below, indicating the revisions we will implement.
read point-by-point responses
-
Referee: Abstract: the description of variational and numerical methods supplies no error estimates, comparisons to data or lattice results, or validation of the quasiparticle mass assumption, leaving the central claim of computing usable gluon PDFs without demonstrated quantitative support.
Authors: We agree that the abstract should more clearly indicate the scope of the calculation. In the revised manuscript we will modify the abstract to state that the gluon PDFs are model predictions obtained within the constituent-gluon Born-Oppenheimer framework, that error control is provided by the agreement between the variational and numerical solutions of the Schrödinger equation (detailed in Sections 3 and 4), and that direct comparisons with lattice QCD or experimental data lie outside the present study but are a natural extension. The quasiparticle mass is treated as an input parameter whose value is taken from the instanton literature; its role is discussed in the main text. revision: yes
-
Referee: The section on the constituent-gluon picture and variational derivation of BO potentials: the gluon is modeled as a massive quasiparticle from instanton-induced interactions without explicit separation of this free parameter from predictions or sensitivity tests, which risks circularity when the resulting wave functions are used to extract PDFs.
Authors: The gluon mass enters the model as a single free parameter whose origin is the instanton-induced effective interaction. The variational derivation of the BO potentials itself is performed for arbitrary mass and does not presuppose a numerical value. To remove any appearance of circularity we will add an explicit sensitivity analysis (new subsection or appendix) in which the mass is varied over a physically motivated interval and the resulting changes in the light-front wave functions and PDFs are quantified. This will demonstrate which features are robust and which depend on the precise mass choice. revision: yes
Circularity Check
No significant circularity: derivation uses standard variational/numerical methods on assumed inputs without reducing predictions to fits by construction
full rationale
The paper's chain starts from an effective constituent-gluon model with instanton-generated mass (an external assumption), applies a variational ansatz to generate BO potentials, solves the Schrödinger equation variationally or numerically in momentum space to obtain light-front wave functions, and finally computes gluon PDFs from those wave functions. None of these steps is shown to be equivalent to its inputs by definition or by fitting the target PDFs themselves; the PDFs are downstream outputs of the wave functions rather than parameters adjusted to match data. Self-citations or prior instanton work are not load-bearing for the central LF wave-function derivation, which remains self-contained against external benchmarks like lattice or other models. This is the normal non-circular outcome for a model-building paper.
Axiom & Free-Parameter Ledger
free parameters (1)
- gluon quasiparticle mass
axioms (2)
- domain assumption Born-Oppenheimer approximation applies to heavy hybrid hadrons with gluonic excitations
- domain assumption Instanton effects generate an effective mass for the constituent gluon
invented entities (1)
-
constituent gluon quasiparticle
no independent evidence
Reference graph
Works this paper leans on
-
[1]
,[21] ,[111] . (116) The color-singlet hybrid is obtained by coupling a color-octetqqqcore to the octet gluon, (qqq)8 ⊗g 8 →1 C,(117) so the three-quark core must carry mixed permuta- tion symmetry [21]C ∼ . A. LF kinematics and Jacobi variables We use the same LF kinematics and Jacobi mo- menta as in the text. The parton fractions satisfy xi ≥0, x 1 +x 2...
-
[2]
These are the natural fixed-mq spin-flavor basis states for the quark core. H. Antisymmetricqqqcore on the LF The fully antisymmetric quark core is obtained by taking the[111]component of[21] C ⊗[21] SF Ψcore qqq (N, mq) = 1√ 2 Cρ Φλ SF (N, mq)− C λ Φρ SF (N, mq) . (150) This is the unique lowest quark-core structure com- patible with the Pauli principle,...
-
[3]
Depending on the interaction, parity, and additional dynamical assumptions, notalloftheseneedcontributeequally. Similarly, for fixed Jz = +3 2 (164) one may use |1⟩= Φ LF +1 Ψcore qqq (N,+ 1 2)|1,0⟩,(165) |2⟩= Φ LF +1 Ψcore qqq (N,− 1 2)|1,+1⟩,(166) |3⟩= Φ LF −1 Ψcore qqq (N,+ 1 2)|1,+1⟩,(167) if a constituent-gluon basis including theλ g = 0 polarization...
-
[4]
bridging scale
ItisagainusedtocalculatethegluonPDF,shown in Fig. 7. It is important to note that LF wave functions provide a complete description of the correspond- ing states. While in this work we have focused on single-particle momentum distributions (PDFs), these wave functions can be used to compute a wide range of observables, including TMDs, GPDs, form factors, a...
-
[5]
One finds ⟨⃗ p2⟩Π = 2 a2 ⊥ + 1 2a2 ∥ ,⟨⃗ p 2⟩Σ = 1 a2 ⊥ + 3 2a2 ∥ ,(D6) and therefore TΠ =m g + 1 2mg 2 a2 ⊥ + 1 2a2 ∥ ! , TΣ =m g + 1 2mg 1 a2 ⊥ + 3 2a2 ∥ ! .(D7)
Kinetic The kinetic expectation values are most easily evaluated from⟨⃗ p2⟩= R d3x|∇ψ| 2. One finds ⟨⃗ p2⟩Π = 2 a2 ⊥ + 1 2a2 ∥ ,⟨⃗ p 2⟩Σ = 1 a2 ⊥ + 3 2a2 ∥ ,(D6) and therefore TΠ =m g + 1 2mg 2 a2 ⊥ + 1 2a2 ∥ ! , TΣ =m g + 1 2mg 1 a2 ⊥ + 3 2a2 ∥ ! .(D7)
-
[6]
Coulomb The Coulomb matrix elements are I(Π) C (r;a ⊥, a∥) = Z d3x|ψ Π(⃗ x)|2 1 |⃗ x−⃗ r/2|+ 1 |⃗ x+⃗ r/2| , (D8) I(Σ) C (r;a ⊥, a∥) = Z d3x|ψ Σ(⃗ x)|2 1 |⃗ x−⃗ r/2|+ 1 |⃗ x+⃗ r/2| , (D9) 23 which, after the trivialϕintegration, reduce to the two-dimensional forms I(Π) C (r;a ⊥, a∥) = 2√π a4 ⊥a∥ Z ∞ 0 dρ ρ3e−ρ2/a2 ⊥ Z ∞ −∞ dz e−z2/a2 ∥ 1p ρ2 + (z−r/2) 2 +...
-
[7]
Confining The confining matrix elements are defined analogously I(Π) σ (r;a ⊥, a∥) = Z d3x|ψ Π(⃗ x)|2 ⃗ x−⃗ r 2 + ⃗ x+⃗ r 2 ,(D12) I(Σ) σ (r;a ⊥, a∥) = Z d3x|ψ Σ(⃗ x)|2 ⃗ x−⃗ r 2 + ⃗ x+⃗ r 2 .(D13) After the angular integration these become I(Π) σ (r;a ⊥, a∥) = 2√π a4 ⊥a∥ Z ∞ 0 dρ ρ3e−ρ2/a2 ⊥ Z ∞ −∞ dz e−z2/a2 ∥ p ρ2 + (z−r/2) 2 + p ρ2 + (z+r/2) 2 , (D14)...
-
[8]
Instanton For a Gaussian instanton-induced interaction, Vinst(⃗ x) =−GΓe−x2/ρ2 0 ,(D16) withρ 0 the instanton size, the expectation values are elementary. For theΠu ansatz one obtains IΠ(a⊥, a∥) = Z d3x|ψ Π(⃗ x)|2e−x2/ρ2 0 = a−2 ⊥ a−2 ⊥ +ρ −2 0 2 a−2 ∥ a−2 ∥ +ρ −2 0 !1/2 ,(D17) while for theΣ− u ansatz one finds IΣ(a⊥, a∥) = Z d3x|ψ Σ(⃗ x)|2e−x2/ρ2 0 = a−...
-
[9]
III is recovered in the isotropic limit
Connection to isotropic ansatz It is essential to verify that the spherical ansatz of Sec. III is recovered in the isotropic limit. Setting a⊥ =a ∥ = 1 β ,(D22) one finds ψ(m=±1) Π ∝ρ e imϕe−β2(ρ2+z2)/2 ∝r e −β2r2/2Y1,±1(ˆx), (D23) and ψΣ ∝z e −β2(ρ2+z2)/2 ∝r e −β2r2/2Y10(ˆx).(D24) The kinetic terms reduce correctly to ⟨⃗ p2⟩Π = 5 2 β2,⟨⃗ p 2⟩Σ = 5 2 β2,(...
-
[10]
Glueballs, Constituent Gluons and Instantons
E. Shuryak and I. Zahed, (2026), arXiv:2604.04803 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2026
- [11]
-
[12]
M. Berwein, N. Brambilla, A. Mohapatra, and A. Vairo, Phys. Rev. D110, 094040 (2024), arXiv:2408.04719 [hep-ph]
- [13]
-
[14]
E. Braaten, C. Langmack, and D. H. Smith, Phys. Rev. D90, 014044 (2014), arXiv:1402.0438 [hep-ph]
- [15]
- [16]
-
[17]
Schlosser and M
C. Schlosser and M. Wagner, Phys. Rev. D105, 054503 (2022)
2022
-
[18]
E. Shuryak and I. Zahed, Phys. Rev. D107, 034026 (2023), arXiv:2202.00167 [hep-ph]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.