⚛️ lattice

Renormalon subtracted nonrelativistic QCD for heavy hadron systems

This paper presents a renormalon-subtracted potential nonrelativistic QCD framework that combines variational and Green's function Monte Carlo methods with NNLO potentials to predict the masses of fully-heavy baryons and the binding properties of fully-heavy tetraquarks, achieving improved perturbative stability while identifying discrepancies with lattice QCD that are consistent with neglected $1/m_Q$ corrections.

Original authors: Benoît Assi, Andreas S. Kronfeld, Simon Vaiva, Michael L. Wagman

Published 2026-07-13

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Benoît Assi, Andreas S. Kronfeld, Simon Vaiva, Michael L. Wagman

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to build a house out of invisible, super-heavy Lego bricks. These bricks are "heavy quarks," the tiny building blocks that make up particles like protons and neutrons. Physicists have a set of blueprints called QCD (Quantum Chromodynamics) to figure out how these bricks stick together. But there's a catch: when the physicists try to calculate how heavy the final house will be, their math starts to go haywire.

The numbers in their equations start growing so fast—like a snowball rolling down a hill that suddenly turns into a mountain—that the answer becomes useless. It's as if every time they tried to add a new layer of detail to the blueprint, the estimated weight of the house jumped by 100 units, then 200, then 1,000, making it impossible to know if the house weighs 10 tons or 10,000 tons. This mathematical glitch is called a renormalon, and it's been a major headache for decades.

The Magic Eraser: Minimal Renormalon Subtraction (MRS)

In this paper, the authors introduce a clever trick called Minimal Renormalon Subtraction (MRS). Think of the renormalon glitch as a sticky, greasy spot on the Lego bricks that makes them impossible to stack neatly. The MRS method acts like a specialized solvent that dissolves that specific grease. It doesn't just ignore the problem; it isolates the "greasy" part of the math, sums it up into a single, clean correction, and removes it from the main calculation.

Once they use this "magic eraser," the math suddenly behaves. Instead of the numbers jumping wildly, they settle down. The authors found that by using MRS, the uncertainty in their weight predictions shrank dramatically. For some heavy particles, the range of possible weights went from a wild swing of 700 MeV (a huge difference in particle physics) down to a tiny 10 MeV. It's the difference between guessing a person's weight by saying "somewhere between a bicycle and a car" versus saying "somewhere between 150 and 151 pounds."

Building the Heavy Houses

With their new, stable math, the team built models of three types of heavy "houses":

Mesons: Pairs of heavy quarks (like a charm quark and an anti-charm quark).
Baryons: Trios of heavy quarks (like three charm quarks or three bottom quarks).
Tetraquarks: Four-quark structures (like two heavy quarks and two heavy anti-quarks).

They tuned their model using known particles (like the $J/\psi$ and $\Upsilon$ ) to get the "weight" of the individual quarks just right. Then, they predicted the masses of brand-new, super-heavy combinations that haven't been measured yet, such as the $\Omega_{ccc}$ (three charm quarks) and the $\Omega_{bbb}$ (three bottom quarks).

The "Missing Piece" Mystery

Here is where it gets interesting. When the authors compared their MRS-predicted masses to the "gold standard" measurements from Lattice QCD (a different, super-powerful computer simulation method), they found a small but consistent gap.

Their predictions were 125–175 MeV lighter than the Lattice QCD results.
This gap wasn't random noise; it shrank as the quarks got heavier. For the lightest heavy quarks (charm), the gap was about 3.8%, but for the heaviest (bottom), it dropped to 1.0%.

The authors are confident this isn't because their math is broken. Instead, they suggest they are missing a specific "furniture" in their house. Their current model treats the quarks as if they are perfectly still (static). However, in reality, these quarks wiggle and spin. The authors argue that the missing 125–175 MeV is likely due to these missing "wiggles" (relativistic corrections) and spin effects that their static model doesn't include yet. It's like predicting the weight of a car by only weighing the metal frame and forgetting to add the engine and the seats. The math is solid, but the model is slightly incomplete.

The Top Quark: The Ghost in the Machine

The team also looked at particles containing top quarks. Top quarks are so heavy and unstable that they decay (fall apart) before they can ever form a real, lasting particle. You can't actually build a "top house" in a lab.

However, the authors used their MRS method to predict what these "ghost houses" would weigh if they could exist. They calculated the masses of top-containing baryons (like $\Omega_{ttt}$ ) and tetraquarks. Even though these particles vanish instantly, their calculations suggest that the "ghost" effects of top quarks could leave a fingerprint in particle colliders like the LHC, perhaps showing up as a slight bump in the data when top quarks are created near a specific energy threshold.

The Four-Quark Puzzle

Finally, the team tackled the tricky question of tetraquarks (four-quark particles). Previous studies suggested that if you have four heavy quarks of the same weight (like four bottom quarks), they won't stick together to form a bound state; they'd just drift apart.

The authors confirmed this: equal-mass tetraquarks like $bb\bar{b}\bar{b}$ do not bind. However, they found that if you mix the weights—making two quarks very heavy and two slightly lighter—the math changes. They discovered a "critical ratio" for binding. If the heavy quarks are more than about 12 times heavier than the lighter ones (in their fixed-order math) or 11 times heavier (in their MRS math), the four quarks snap together into a stable molecule.

Using this, they predicted that combinations like $bb\bar{t}\bar{t}$ (two bottom quarks and two top antiquarks) and $cc\bar{t}\bar{t}$ would indeed form bound states, with binding energies ranging from 476 MeV to 995 MeV.

The Bottom Line

This paper doesn't claim to have solved the entire mystery of heavy particles. Instead, it suggests that by using the MRS method, physicists can finally get a stable, reliable view of how these heavy quarks interact. The math is now stable enough to see clearly that there is a small, consistent gap between their static models and reality, likely caused by the quarks' motion and spin.

The authors propose that this framework is a powerful new tool. It acts like a high-precision telescope that has finally been focused, allowing them to predict the properties of exotic, heavy particles with a level of confidence that was previously impossible, paving the way for future experiments to hunt down these heavy "ghosts" and "molecules" in the data.

Technical Summary: Renormalon Subtracted Nonrelativistic QCD for Heavy Hadron Systems

Problem Statement
Precision spectroscopy of heavy hadron systems (quarkonia, triply-heavy baryons, and fully-heavy tetraquarks) using Potential Nonrelativistic QCD (pNRQCD) has historically faced a significant obstacle: the renormalon ambiguity inherent in the static potential. In fixed-order perturbative calculations, the coefficients of the static potential exhibit factorial growth, leading to renormalization-scale ( $\mu$ ) variations that induce shifts of order 100 MeV in binding energies. These shifts are comparable to the target precision required to match lattice QCD benchmarks. Furthermore, successive perturbative orders (LO, NLO, NNLO) often fail to converge, with uncertainty bands failing to overlap, particularly for systems involving multiple heavy quarks where delicate cancellations among color channels amplify sensitivity to the static potential.

Methodology
The authors implement a Minimal Renormalon Subtraction (MRS) scheme within a pNRQCD framework combined with Variational Monte Carlo (VMC) and Green's Function Monte Carlo (GFMC) methods.

MRS Formalism: The MRS procedure reorganizes the perturbative series by isolating the leading renormalon ambiguity into an explicit nonperturbative power correction ( $\Lambda/Q$ ). This is achieved by:
1. Subtracting the asymptotic renormalon contributions from fixed-order potential coefficients.
2. Resumming the remaining tail of the series using a Borel integral.
3. Absorbing the resulting constant into the quark masses (for singlet potentials) or neglecting it in other cases.
  The MRS static potential is evaluated at each spatial point $r$ using a dynamical scale $\mu'(r)$ that interpolates smoothly between the natural scale $1/r$ and a fixed cutoff $\mu_{cut}$ to avoid derivative discontinuities and ensure stable stochastic sampling.
Computational Approach:
- Hamiltonian: The system is described by a Schrödinger equation with a Hamiltonian containing nonrelativistic kinetic energy and static potentials (two-body and, at NNLO, three-body). Spin-dependent operators are neglected, consistent with the static approximation ( $O(m_Q^0)$ ).
- Wavefunctions: Ground-state energies are extracted using trial wavefunctions. For quarkonia and baryons, Coulombic ansätze are used. For tetraquarks, Hylleraas-type molecular trial wavefunctions are employed to capture correlations at different length scales.
- Extraction: GFMC projects the exact ground state from trial wavefunctions via imaginary-time evolution. Statistical uncertainties are estimated via bootstrap resampling, while systematic errors from finite projection time and discretization are assessed by varying parameters.
Parameter Tuning: Charm ( $m_c$ ) and bottom ( $m_b$ ) quark masses are tuned at each perturbative order to reproduce the experimentally measured spin-averaged $1S$ quarkonium masses ( $J/\psi, \eta_c$ and $\Upsilon, \eta_b$ ). These tuned masses serve as fixed inputs for all subsequent predictions, rendering the mixed-flavor and multi-quark predictions parameter-free.

Key Contributions and Results

Stabilization of Perturbative Series:
The application of MRS dramatically stabilizes the perturbative expansion.
- Scale Dependence: Scale variation bands for binding energies are reduced by factors of 3–10. For example, in $\Omega_{bbb}$ , the NNLO scale variation drops from $\sim 700$ MeV (fixed-order) to $\sim 10$ MeV (MRS).
- Convergence: Order-by-order shifts shrink significantly. In the MRS scheme, shifts between NLO and NNLO for baryon masses are reduced to 10–30 MeV (compared to 50–80 MeV in fixed-order), and the uncertainty bands for LO, NLO, and NNLO overlap, signaling a well-behaved perturbative treatment.
Quark Mass Determination:
MRS-tuned quark masses exhibit sub-percent variation across perturbative orders (e.g., $m_b$ varies from 4.81 GeV to 4.80 GeV), whereas fixed-order masses show variations of 4–11%. This confirms that MRS successfully isolates and resums the leading renormalon.
Triply-Heavy Baryon Spectroscopy:
Using the tuned masses, the authors predict masses for $\Omega_{ccc}$ , $\Omega_{ccb}$ , $\Omega_{cbb}$ , and $\Omega_{bbb}$ .
- Comparison with Lattice QCD: At NNLO with MRS, predictions undershoot lattice QCD results by 125–175 MeV (approx. 1–2% of the baryon mass).
- Interpretation: The fractional discrepancy decreases as $\sim 1/m_Q$ (from 3.8% for $\Omega_{ccc}$ to 1.0% for $\Omega_{bbb}$ ). The authors attribute this residual difference to neglected $O(1/m_Q)$ velocity-dependent corrections and $O(1/m_Q^2)$ spin-dependent terms (spin-orbit, spin-spin, Darwin) present in lattice-NRQCD but absent in their static potential. The reduced scale variation in MRS makes this discrepancy unambiguous, distinguishing it from perturbative instability.
- Top-Containing Baryons: Predictions are provided for $\Omega_{tcc}$ , $\Omega_{tbb}$ , and $\Omega_{ttt}$ , serving as benchmarks for the $m_Q \to \infty$ limit and near-threshold dynamics at the LHC.
Fully-Heavy Tetraquarks:
The framework is extended to $QQ\bar{Q}'\bar{Q}'$ systems to map the binding landscape as a function of the heavy-to-light mass ratio $m_Q/m_{Q'}$ .
- Critical Ratio: The authors determine the critical mass ratio for binding. In the MRS scheme, the threshold for binding shifts to $m_Q/m_{Q'} \approx 0.085$ (NLO), lower than the fixed-order threshold of $\approx 0.12$ .
- Bound States: Bound states are predicted for $bb\bar{t}\bar{t}$ , $bc\bar{t}\bar{t}$ , and $cc\bar{t}\bar{t}$ . No binding is observed for equal-mass or near-equal-mass combinations (e.g., $bb\bar{b}\bar{b}$ , $cc\bar{b}\bar{b}$ ) where the mass ratio exceeds the critical threshold.
- Binding Energies: For bound top-containing tetraquarks at fixed-order NLO, binding energies range from 400–600 MeV below the two-meson threshold.

Significance
The paper claims that this work represents the first application of Minimal Renormalion Subtraction (MRS) to bound-state problems in pNRQCD. The primary significance lies in the demonstration that MRS can systematically remove the leading renormalon ambiguity from the static potential, thereby:

Restoring stable power counting and improving perturbative convergence for multi-heavy-quark systems.
Reducing renormalization-scale dependence by roughly two orders of magnitude in some cases, allowing for a clear interpretation of discrepancies with lattice QCD as missing physical corrections (relativistic/spin effects) rather than theoretical artifacts.
Providing a parameter-free, controlled expansion for predicting the spectra of triply-heavy baryons and fully-heavy tetraquarks, including those containing top quarks, which serve as theoretical benchmarks for heavy-quark dynamics in the infinite mass limit.

The authors note that while the static potential treatment is successful, full quantitative control requires future inclusion of $O(1/m_Q)$ relativistic corrections, MRS for non-singlet colored channels, and resummation of three-body potentials.