Non-Relativistic QCD Hamiltonian

• Non-relativistic QCD Hamiltonian is a framework that integrates out high-energy QCD modes to focus on low-energy dynamics of heavy quarks, like bottomonium and charmonium.
• It employs a systematic operator expansion including terms up to O(p^6) and one-loop radiative corrections to improve dispersion matching and mass renormalisation.
• Key lattice techniques such as Fat³ smearing and stability parameter tuning reduce discretisation errors, enabling precise quarkonium spectroscopy and decay prediction.

Non-relativistic QCD (NRQCD) Hamiltonians provide a systematic framework for describing low-energy dynamics of heavy quarks, notably bottomonium and charmonium states, within both continuum and lattice approaches. Central aspects include operator expansions in powers of momentum and velocity, accurate matching of lattice and continuum dispersion relations, and operator improvement through radiative corrections and lattice stabilisation techniques. These structures underpin quantitative predictions in heavy quark physics at the percent-level and are foundational for calculations of quarkonium spectra, transition rates, and non-perturbative QCD phenomena.

1. Hamiltonian Structure and Operator Expansion

The NRQCD Hamiltonian is constructed by integrating out hard QCD modes above the heavy quark mass scale, leading to a non-relativistic theory with quark fields that evolve under an effective Hamiltonian: $$H = H_0 + \delta H_{\text{kin}} + \delta H_{\text{rad}}$$ where $H_0$ contains the leading kinetic term, $\delta H_{\text{kin}}$ incorporates higher-order corrections in the operator expansion up to $\mathcal{O}(p^6)$, and $\delta H_{\text{rad}}$ includes radiative (loop-level) shifts (Davies et al., 2018).

Leading and Subleading Operators

• Leading kinetic: $aH_0 = -\Delta^{(2)} / (2 a m_b)$, where $\Delta^{(2)}$ is the lattice Laplacian with link-pair cancellations to suppress tadpoles.
• $\mathcal{O}(v^4)$ corrections: Operators such as $(\Delta^{(2)})^2$, $i[\nabla \cdot \tilde{E} - \tilde{E} \cdot \nabla]$, $H_0$0, $H_0$1, and $H_0$2 are included with matching coefficients $H_0$3.
• $H_0$4 corrections: Operators like $H_0$5, $H_0$6, and $H_0$7 provide further discretisation error reduction.

2. Radiative Improvement and Dispersion-Relation Matching

One-loop radiative corrections renormalise the heavy quark mass and introduce zero-point energy shifts:

• Mass renormalisation $H_0$8 affects the dispersion relation.
• Zero-point energy shift $H_0$9 enters additively.
• Matching to continuum NRQCD is realized via the on-shell energy expansion: $\delta H_{\text{kin}}$0 and at one-loop,

$\delta H_{\text{kin}}$1

yielding kinetic matching coefficients $\delta H_{\text{kin}}$2 and $\delta H_{\text{kin}}$3, so all $\delta H_{\text{kin}}$4 mismatches between lattice and continuum are eliminated (Davies et al., 2018).

3. Lattice Stabilisation and Operator Improvement Techniques

Stability Parameter $\delta H_{\text{kin}}$5 and Fat$\delta H_{\text{kin}}$6 Smearing

• The evolution operator $\delta H_{\text{kin}}$7 requires $\delta H_{\text{kin}}$8 for stability, with $\delta H_{\text{kin}}$9 the typical choice.
• Fat$\mathcal{O}(p^6)$0 smearing replaces every thin gauge link by a locally averaged link, projected back to $\mathcal{O}(p^6)$1, nearly eliminating unphysical tadpole contributions. This permits setting $\mathcal{O}(p^6)$2 and reduces the magnitude of radiative corrections, making expansion truncations reliable (Davies et al., 2018).

Tadpole Improvement and Matching Scale

Radiative coefficients are either tadpole-improved via mean-field rescaling ($\mathcal{O}(p^6)$3 in Landau gauge) or, equivalently, Fat$\mathcal{O}(p^6)$4 smearing can be used (where all mean-field counterterms vanish). The $\mathcal{O}(p^6)$5 coupling is evaluated at a process-dependent BLM scale $\mathcal{O}(p^6)$6, with $\mathcal{O}(p^6)$7–$\mathcal{O}(p^6)$8 for $\mathcal{O}(p^6)$9, $\delta H_{\text{rad}}$0 and $\delta H_{\text{rad}}$1–$\delta H_{\text{rad}}$2 for $\delta H_{\text{rad}}$3, $\delta H_{\text{rad}}$4 (Davies et al., 2018). This ensures perturbative control and correct matching of kinetic terms.

4. Fully Improved Lattice NRQCD Hamiltonian

Integrating all components, the improved Hamiltonian used in simulations is: $\delta H_{\text{rad}}$5 where each coefficient $\delta H_{\text{rad}}$6 is expanded as $\delta H_{\text{rad}}$7, with all operators discretised on the lattice. Annotations indicate mapping to continuum terms and stabilization mechanisms (Davies et al., 2018).

5. Implications for Spectroscopy and Precision Calculations

By systematically including radiative corrections and higher-order momentum operators, discrete errors are reduced to $\delta H_{\text{rad}}$8, rotational symmetry is restored (as confirmed through momentum tests), and the precision of quarkonium mass predictions is improved to the 1–2% regime (Davies et al., 2018). Mass renormalisation and zero-point energy shifts allow conversion between lattice results and continuum pole masses, crucial for physical observables.

Lattice ensembles incorporating $\delta H_{\text{rad}}$9, $aH_0 = -\Delta^{(2)} / (2 a m_b)$0, $aH_0 = -\Delta^{(2)} / (2 a m_b)$1, and $aH_0 = -\Delta^{(2)} / (2 a m_b)$2 vacuum polarisation are used to non-perturbatively determine states such as $aH_0 = -\Delta^{(2)} / (2 a m_b)$3, $aH_0 = -\Delta^{(2)} / (2 a m_b)$4, and to tune the $aH_0 = -\Delta^{(2)} / (2 a m_b)$5-quark mass. Improvements in action further stabilise time-evolution algorithms and correct for high-momentum artefacts.

6. Relation to Continuum and Effective Field Theory Approaches

The NRQCD Hamiltonian is tightly connected with pNRQCD and continuum effective descriptions. In pNRQCD, the Hamiltonian incorporates matched potentials (static, $aH_0 = -\Delta^{(2)} / (2 a m_b)$6, $aH_0 = -\Delta^{(2)} / (2 a m_b)$7, spin-dependent) with coefficients determined up to two-loop order (Mishima et al., 2024, Segovia et al., 2017), using either two-step or single-step matching procedures, IBP reduction to master integrals, and differential equations in the velocity expansion parameter. These ensure that all heavy-quarkonium binding energies and transitions can be consistently computed up to N$aH_0 = -\Delta^{(2)} / (2 a m_b)$8LO in the weak-coupling regime. The approach also provides a basis for further extensions, such as BCS quasiparticle formation in non-perturbative QCD Hamiltonians (Yepez-Martinez et al., 2021), supporting dynamical generation of constituent quark masses and correlations in low-energy hadron spectra.

7. Summary Table: Key Features of Improved Lattice NRQCD Hamiltonian

Feature	Description	Impact
Operator Expansion	Includes terms up to $aH_0 = -\Delta^{(2)} / (2 a m_b)$9	$\Delta^{(2)}$0 error control
Radiative Corrections	One-loop matching for mass and kinetic terms	1–2% level mass accuracy
Fat$\Delta^{(2)}$1 Smearing	Removes unphysical tadpoles, sets $\Delta^{(2)}$2	Perturbative convergence, stable fits
Stability Parameter $\Delta^{(2)}$3	Ensures $\Delta^{(2)}$4 for stable evolution	No transfer-matrix sign flips
Matching Scale $\Delta^{(2)}$5	BLM scale-setting for $\Delta^{(2)}$6 in $\Delta^{(2)}$7-scheme	Reliable coefficient expansion

The non-relativistic QCD Hamiltonian, in both lattice and continuum forms, is now a rigorously improved tool for heavy quark physics. Its comprehensive operator content, precise radiative matching, and stabilization techniques enable systematic reduction of discretisation and perturbative errors, supporting high-precision spectroscopy, decay, and transition predictions for quarkonium systems.