Beyond the second-order valence approximation: a systematically improvable, first-principles light-front Hamiltonian QCD calculation of heavy-quarkonium spectra

In 1997, Brisudová, Perry and Wilson computed the renormalized effective Hamiltonian of light-front QCD to second order in the coupling using the similarity renormalization group with coupling coherence and, within a valence $q\bar{q}$ approximation with a nonrelativistic reduction, obtained a Coulomb-plus-logarithmic confining potential that fits the low-lying charmonium and bottomonium levels and predicts the charmonium hyperfine splitting in reasonable agreement with experiment. They framed the result as an initial crude step toward an accurate first-principles light-front QCD calculation. The open problem is to complete that program: include the mixing of Fock sectors beyond the valence $q\bar{q}$ component (which requires the nonperturbative light-front three-body problem), carry the renormalization-group expansion beyond second order with running couplings and higher-order counterterms, demonstrate restoration of rotational and Lorentz symmetry in observables, and quantify the sensitivity of the fitted parameters to truncations and cutoff choices.

Background

Quantum chromodynamics (QCD), the $SU(3)$ gauge theory of quarks and gluons, has no known analytic solution in the confinement regime, and computing hadron masses and wave functions directly from the theory is a long-standing challenge. One route is Hamiltonian: instead of the usual equal-time quantization, one uses Dirac's front form (light-front or light-cone quantization), with coordinates $x^\pm = x^0 \pm x^3$ and fields quantized at equal $x^+$ . In this form the free-theory vacuum remains the vacuum of the full theory in the presence of cutoffs, so a hadron state can be expanded in Fock sectors — components with a definite number of partons, such as the valence $q\bar{q}$ component of a meson, the $q\bar{q}g$ component with one gluon, and so on. The price is that symmetries such as rotations are no longer kinematic: rotations mix transverse and longitudinal momenta and must emerge dynamically in observables.

Renormalizing a light-front Hamiltonian is subtle because many symmetries are not manifest and locality is violated in the longitudinal direction. Głazek and Wilson (Phys. Rev. D 48, 5863 (1993)) introduced the similarity renormalization group: a unitary similarity transformation that renders the Hamiltonian band-diagonal in the free light-front energy, so that matrix elements couple only states of nearby energy and ultraviolet divergences can be absorbed into cutoff-independent counterterms. Together with coupling coherence — the requirement that the Hamiltonian reproduce its own form under the renormalization group, which fixes the finite parts of the counterterms — this was proposed by Wilson, Walhout, Harindranath, Zhang, Perry and Głazek (Phys. Rev. D 49, 6720 (1994)) as a route to deriving a constituent-parton picture of hadrons from weak-coupling QCD.

Brisudová and Perry (Phys. Rev. D 54, 1831 (1996)) developed the second-order effective Hamiltonian in detail, and Brisudová, Perry and Wilson (Phys. Rev. Lett. 78, 1227 (1997)) applied it to heavy quarkonia — the charmonium ( $c\bar{c}$ ) and bottomonium ( $b\bar{b}$ ) bound states. The effective Hamiltonian to $\mathcal{O}(g^2)$ contains a short-range Coulomb potential and a long-range logarithmic confining potential. For the bound-state calculation it was split as $H_0 + V$ : $H_0$ contains the nonrelativistic reduction of the kinetic energy, the effective one-body terms, the Coulomb potential, and the rotationally symmetric part of the confining potential, and is solved nonperturbatively; the remainder $V$ — including spin-dependent terms and the emission and absorption of low-energy gluons — is treated in bound-state perturbation theory. With the parameters $M$ (quark mass), $\alpha$ (coupling) and $\Lambda$ (cutoff) fit to the 1S, 1P and 2S levels (for charmonium, $M_c = 1.5$ GeV, $\alpha = 0.5$ , $\Lambda = 1.7$ GeV), the calculation predicted the charmonium hyperfine splitting of 0.13 GeV against the experimental $0.118 \pm 0.002$ GeV, and estimated the corrections from rotationally non-invariant terms at a few percent for the ground state and about 10% for the lowest excited states.

The calculation rests on three approximations: the valence $q\bar{q}$ approximation (excluding low-energy gluon emission and absorption from $H_0$ , so that mixing with higher Fock sectors enters only at second order in $V$ and requires solving the nonperturbative three-body light-front problem), the nonrelativistic reduction appropriate for heavy quarks, and the truncation of the similarity renormalization to $\mathcal{O}(g^2)$ . The paper describes its result as 'an initial crude step toward an accurate first principles, light-front QCD calculation'. The open problem is what must be computed to close the gap between this second-order valence approximation and a systematically improvable first-principles calculation. As of 2026 the direct continuation of this program — the renormalization group procedure for effective particles — still works at second order, with an effective-gluon-mass ansatz standing in for higher Fock components, so the problem remains open.

Problem Statement

Starting from the canonical light-front QCD Hamiltonian and its similarity-renormalization-group renormalization (in its original form or its successor, the renormalization group procedure for effective particles), construct a systematically improvable, first-principles light-front Hamiltonian calculation of the heavy-quarkonium (charmonium and bottomonium) spectrum that goes beyond the second-order, valence $q\bar{q}$ , nonrelativistic approximation of Brisudová, Perry and Wilson. The target is a calculation that, within one coherent framework, provides all four of the following, jointly:

(a) Fock-sector mixing beyond valence: a treatment of the higher Fock components (at minimum the $q\bar{q}g$ sector) that computes the corrections from particle-number-changing interactions and states an explicit bound on the residual error from the higher sectors left out, quantified from the contribution of the neglected sectors, including the required nonperturbative three-body light-front bound-state calculation, instead of absorbing those components into an uncontrolled ansatz;

(b) higher orders: the effective Hamiltonian carried beyond second order to at least the order at which the renormalization-group running of the effective coupling first appears (the fourth order of the weak-coupling expansion), with the finite counterterms beyond second order established, so that quantitative predictions carry a truncation error that is estimated rather than assumed;

(c) symmetry restoration: control of rotational and full Lorentz symmetry breaking induced by the regulator and the truncations, demonstrating restoration of these symmetries in observables — rotational multiplets degenerate to within the truncation error estimated in (b) — rather than only estimating the size of their violation; and

(d) parameter stability: a quantified account, over a stated range of the similarity cutoff $\Lambda$ and scale, of how the fitted parameters ( $M$ , $\alpha$ ) and the predicted observables change under higher-order corrections and under variation of that cutoff and scale, with $M$ and $\alpha$ re-fitted at each value of the cutoff.

The unifying objective is a prediction of the heavy-quarkonium spectrum (the charmonium hyperfine splitting being the source calculation's benchmark prediction) whose approximations can be improved systematically and whose input parameters are fixed consistently with QCD rather than tuned to absorb uncomputed effects.

References

M. M. Brisudová, R. J. Perry and K. G. Wilson, 'Quarkonia in Hamiltonian Light-Front QCD', Phys. Rev. Lett. 78, 1227 (1997), DOI 10.1103/PhysRevLett.78.1227, https://doi.org/10.1103/PhysRevLett.78.1227 (preprint arXiv:hep-ph/9607280, https://arxiv.org/abs/hep-ph/9607280)
S. D. Głazek and K. G. Wilson, 'Renormalization of Hamiltonians', Phys. Rev. D 48, 5863 (1993), DOI 10.1103/PhysRevD.48.5863, https://doi.org/10.1103/PhysRevD.48.5863
K. G. Wilson, T. S. Walhout, A. Harindranath, Wei-Min Zhang, R. J. Perry and S. D. Głazek, 'Nonperturbative QCD: A weak-coupling treatment on the light front', Phys. Rev. D 49, 6720 (1994), DOI 10.1103/PhysRevD.49.6720, https://doi.org/10.1103/PhysRevD.49.6720
M. Brisudová and R. Perry, 'Initial studies of bound states in light-front QCD', Phys. Rev. D 54, 1831 (1996), DOI 10.1103/PhysRevD.54.1831, https://doi.org/10.1103/PhysRevD.54.1831
M. M. Brisudová, S. Szpigel and R. J. Perry, 'Effects of massive gluons on quarkonia in light-front QCD', Phys. Lett. B 421, 334 (1998), DOI 10.1016/S0370-2693(97)01556-6, https://doi.org/10.1016/S0370-2693(97)01556-6 (preprint arXiv:hep-ph/9709479, https://arxiv.org/abs/hep-ph/9709479)
R. J. Perry, 'Using Wilson's renormalization group to repair symmetries', Phys. Rep. 348, 33 (2001), DOI 10.1016/S0370-1573(00)00135-6, https://doi.org/10.1016/S0370-1573(00)00135-6
M. Brisudová, 'Small x divergences in a heavy quark-antiquark state', Nucl. Phys. B Proc. Suppl. 108, 173 (2002), DOI 10.1016/S0920-5632(02)01324-5, https://doi.org/10.1016/S0920-5632(02)01324-5
S. D. Głazek and J. Młynik, 'Boost-invariant Hamiltonian approach to heavy quarkonia', Phys. Rev. D 74, 105015 (2006), DOI 10.1103/PhysRevD.74.105015, https://doi.org/10.1103/PhysRevD.74.105015 (preprint arXiv:hep-th/0606235, https://arxiv.org/abs/hep-th/0606235)
S. D. Głazek, M. Gómez-Rocha, J. More and K. Serafin, 'Renormalized quark-antiquark Hamiltonian induced by a gluon mass ansatz in heavy-flavor QCD', Phys. Lett. B 773, 172 (2017), DOI 10.1016/j.physletb.2017.08.018, https://doi.org/10.1016/j.physletb.2017.08.018 (preprint arXiv:1705.07629, https://arxiv.org/abs/1705.07629)
K. Serafin, M. Gómez-Rocha, J. More and S. D. Głazek, 'Dynamics of heavy quarks in the Fock space', Phys. Rev. D 109, 016017 (2024), DOI 10.1103/PhysRevD.109.016017, https://doi.org/10.1103/PhysRevD.109.016017 (preprint arXiv:2310.00365, https://arxiv.org/abs/2310.00365)
S. D. Głazek, 'Renormalization group procedure for effective particles', Eur. Phys. J. Spec. Top. 235, 1513 (2026), DOI 10.1140/epjs/s11734-026-02127-y, https://doi.org/10.1140/epjs/s11734-026-02127-y
S. D. Głazek, 'Gluon mass and small-x dynamics in hadrons', preprint arXiv:2606.29455 (2026), https://arxiv.org/abs/2606.29455
S. J. Brodsky, H.-C. Pauli and S. S. Pinsky, 'Quantum chromodynamics and other field theories on the light cone', Phys. Rep. 301, 299 (1998), DOI 10.1016/S0370-1573(97)00089-6, https://doi.org/10.1016/S0370-1573(97)00089-6 (preprint arXiv:hep-ph/9705477, https://arxiv.org/abs/hep-ph/9705477)
B. L. G. Bakker, A. Bassetto, S. J. Brodsky, W. Broniowski, S. Dalley, T. Frederico, S. D. Głazek, J. R. Hiller, C.-R. Ji, V. Karmanov, D. Kulshreshtha, J.-F. Mathiot, W. Melnitchouk, G. A. Miller, J. Papavassiliou, W. N. Polyzou, N. G. Stefanis, J. P. Vary, A. Ilderton and T. Heinzl, 'Light-front quantum chromodynamics: a framework for the analysis of hadron physics', Nucl. Phys. B Proc. Suppl. 251-252, 165 (2014), DOI 10.1016/j.nuclphysbps.2014.05.004, https://doi.org/10.1016/j.nuclphysbps.2014.05.004 (preprint arXiv:1309.6333, https://arxiv.org/abs/1309.6333)
J. R. Hiller, 'Nonperturbative light-front Hamiltonian methods', Prog. Part. Nucl. Phys. 90, 75 (2016), DOI 10.1016/j.ppnp.2016.06.002, https://doi.org/10.1016/j.ppnp.2016.06.002 (preprint arXiv:1606.08348, https://arxiv.org/abs/1606.08348)
Y. Li, P. Maris, X. Zhao and J. P. Vary, 'Heavy quarkonium in a holographic basis', Phys. Lett. B 758, 118 (2016), DOI 10.1016/j.physletb.2016.04.065, https://doi.org/10.1016/j.physletb.2016.04.065 (preprint arXiv:1509.07212, https://arxiv.org/abs/1509.07212)
S. Xu, Y. Liu, C. Mondal, J. Lan, X. Zhao, Y. Li and J. P. Vary, 'Towards a first principles light-front Hamiltonian for the nucleon', Phys. Lett. B 867, 139599 (2025), DOI 10.1016/j.physletb.2025.139599, https://doi.org/10.1016/j.physletb.2025.139599 (preprint arXiv:2408.11298, https://arxiv.org/abs/2408.11298)