⚛️ phenomenology

Bottom-Charm $[bc]$ Diquark Masses from QCD Laplace Sum-Rules

This paper utilizes QCD Laplace sum rules with next-to-leading-order corrections and numerical loop integration techniques to predict the constituent masses of bottom-charm diquarks across four spin-parity channels, providing essential QCD-based benchmarks for constituent diquark models of tetraquark states.

Original authors: T. de Oliveira, D. Harnett, J. Ho, R. T. Kleiv, Siyuan Li, N. Loewen, A. Palameta, T. G. Steele

Published 2026-10-15

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: T. de Oliveira, D. Harnett, J. Ho, R. T. Kleiv, Siyuan Li, N. Loewen, A. Palameta, T. G. Steele

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

Deep within the fabric of the universe, matter is built from a small family of fundamental particles called quarks. Usually, these quarks stick together in groups of three to form protons and neutrons, the building blocks of atoms, or in pairs to form mesons. However, in recent years, physicists have discovered a growing zoo of exotic particles that do not fit these simple patterns. Some of these strange states appear to be made of four quarks bound together, known as tetraquarks. To understand how these complex structures hold together, scientists often imagine them as being composed of two smaller, tightly bound pairs of quarks, called diquarks. Think of a diquark as a single, heavy brick made of two quarks glued together, which then combines with another brick to form the larger tetraquark structure. The mass of this diquark brick is a critical piece of information; without knowing how heavy it is, models trying to predict the properties of these exotic particles are like architects trying to build a house without knowing the weight of their bricks.

For decades, researchers have used a powerful theoretical tool called QCD sum rules to estimate the masses of these diquarks. This method allows scientists to calculate properties of particles by analyzing the mathematical relationships between the fundamental forces of nature and the vacuum of space itself. While this technique has been successfully applied to many types of quark combinations, a specific and important pair had remained a gap in the knowledge: a diquark made of one bottom quark and one charm quark. This combination is particularly interesting because it involves two very heavy quarks, which creates a complex mathematical challenge that previous methods could not fully solve. Until now, scientists had to rely on rough estimates for this specific type of diquark, leaving a blind spot in our understanding of the exotic particles that might contain it.

A team of researchers has now filled this gap by performing a highly sophisticated calculation to determine the mass of the bottom-charm diquark. They tackled the problem by using advanced computer techniques to solve the complex equations that describe how these heavy quarks interact. The difficulty lay in the fact that the calculation involved three distinct scales of energy at once, creating a mathematical puzzle that was too difficult to solve with standard formulas. Instead of giving up, the team developed a new approach, using numerical methods to compute the necessary values directly. They ran these calculations through a computer program designed to handle the intricate loops of particle interactions, effectively bypassing the need for a closed-form solution that did not exist. This allowed them to include more precise corrections to their theory, moving beyond the rough approximations used in earlier studies.

The results of this rigorous analysis provide a clear picture of the bottom-charm diquark's mass, broken down by its different possible internal configurations. The team found that the mass depends on the specific way the quarks are arranged and spinning. For the most stable, scalar configuration, the mass is approximately 5.60 GeV, with a small margin of error. The axial vector version, which has a slightly different spin arrangement, shares this same mass of 5.60 GeV. The other two configurations, which involve a different type of symmetry, are slightly heavier. The pseudoscalar version weighs in at about 6.04 GeV, while the vector version is nearly the same at 5.99 GeV. These numbers represent the weight of the diquark "brick" itself, a value that can now be used as a solid foundation for building models of the larger tetraquark particles.

The findings also reveal a clear pattern in how these heavy quarks behave. The researchers observed that the two lighter configurations are nearly identical in mass, as are the two heavier ones. This suggests that the interaction between the spins of these two heavy quarks is very weak, a result that aligns with theoretical expectations for systems containing such massive particles. By comparing their new, precise numbers with the values used in other theoretical models, the team confirmed that many existing models are in good agreement with their results, though some older estimates were slightly too high. This work does not just provide a number; it validates the tools used to study these exotic particles and offers a reliable benchmark for future discoveries. As experiments continue to search for new tetraquark states containing bottom and charm quarks, these mass predictions will serve as a crucial guide, helping physicists distinguish between different theories and understand the true nature of these rare forms of matter.

Technical Summary: Bottom-Charm [bc] Diquark Masses from QCD Laplace Sum-Rules

Problem Statement
The internal substructure of exotic hadrons, particularly tetraquarks containing heavy quarks (charm and/or bottom), is often modeled using a "compact scenario" where the state is composed of a pair of strongly interacting colored diquarks. In such constituent diquark models, the constituent mass of the diquark is a crucial parameter. While QCD sum rules have been successfully used to calculate constituent masses for diquarks with various flavor combinations (light-light, light-heavy, and heavy-heavy), a significant gap remained for the bottom-charm ([bc]) diquark. Previous calculations for [bc] diquarks were limited to Leading Order (LO) perturbative contributions and truncated Operator Product Expansions (OPE) at dimension-four. Furthermore, existing analytic expressions for the multi-scale master integrals required for a full Next-to-Leading-Order (NLO) calculation involving three distinct scales ( $m_c$ , $m_b$ , and $q^2$ ) were insufficient for renormalization. This work aims to close this gap by performing a comprehensive NLO QCD sum-rule analysis for [bc] diquarks, including dimension-six condensates and extending the analysis to negative parity channels.

Methodology
The authors calculate the correlation function of the [bc] color-antitriplet diquark currents within the OPE up to NLO in perturbation theory and including non-perturbative corrections up to dimension-six.

Currents and Channels: The study considers four channels defined by the quantum numbers $J^P \in \{0^+, 0^-, 1^+, 1^-\}$ , utilizing specific diquark currents involving bottom ( $b$ ) and charm ( $c$ ) quark fields.
Perturbative Calculation: The calculation includes LO and NLO perturbative diagrams, as well as counter-terms for renormalization. The OPE incorporates dimension-four and dimension-six gluon condensates ( $\langle \alpha_s G^2 \rangle$ , $\langle g^3 G^3 \rangle$ ) and the dimension-six quark condensate ( $\langle \bar{q}q \rangle^2$ ).
Numerical Integration: A primary technical challenge is the evaluation of multi-scale master integrals arising from the distinct mass scales of the charm and bottom quarks. Since existing analytic expressions for these integrals were not known to the required order in $\epsilon$ (for dimensional regularization) to handle renormalization, the authors employed a numerical approach. They utilized the pySecDec package to perform D-dimensional loop integration numerically.
Sum Rule Construction: The QCD correlation function is transformed into Laplace Sum-Rules (LSR) via the Borel transform. The analysis involves:
- Constructing the spectral function with a continuum threshold $s_0$ .
- Using a complex contour integration in the $q^2$ plane (avoiding the branch cut) to evaluate the LSRs numerically.
- Determining the Borel window ( $\tau$ ) and continuum threshold ( $s_0$ ) by satisfying OPE convergence, Hölder inequalities, and pole dominance criteria (requiring the resonance contribution to be $\ge 50\%$ ).
- Extracting the constituent mass $M$ directly from the ratio of the first moment to the zeroth moment of the sum rule ( $M = R_1/R_0$ ).
Parameters: The analysis uses running quark masses and the strong coupling constant evaluated at an intermediate renormalization scale $\bar{\mu} \approx 2.73$ GeV. Uncertainties are estimated via a Monte Carlo analysis of input parameters and a conservative 10% methodological uncertainty.

Key Contributions

First NLO Calculation for [bc]: This work provides the first QCD sum-rule determination of [bc] diquark masses including NLO perturbative corrections and dimension-six condensates.
Extension to Negative Parity: The analysis extends beyond previous LO studies to include pseudoscalar ( $0^-$ ) and vector ( $1^-$ ) channels, in addition to the scalar ( $0^+$ ) and axial-vector ( $1^+$ ) channels.
Numerical Technique for Multi-Scale Integrals: The paper demonstrates a robust methodology for handling NLO calculations with multiple mass scales by combining analytic Wilson coefficients with numerical master integral evaluation via pySecDec. This approach overcomes the lack of sufficient analytic expressions for the required $\epsilon$ -expansions.
Comprehensive Survey: The results are integrated into a broader summary of QCD sum-rule predictions for all possible quark-flavor diquark combinations, providing a consistent QCD-based benchmark for constituent diquark models.

Results
The sum-rule analyses stabilize in all four channels, yielding the following constituent mass predictions (in GeV):

Scalar ( $0^+$ ): $5.60 \pm 0.56$
Pseudoscalar ( $0^-$ ): $6.04 \pm 0.60$
Axial Vector ( $1^+$ ): $5.60 \pm 0.56$
Vector ( $1^-$ ): $5.99 \pm 0.60$

The results exhibit a clear physical pattern:

Parity Splitting: There is a distinct mass splitting between positive parity states ( $\sim 5.6$ GeV) and negative parity states ( $\sim 6.0$ GeV), attributed to the orbital excitation required for negative parity.
Spin-Spin Suppression: Within each parity sector, the scalar and axial-vector (or pseudoscalar and vector) masses are nearly degenerate. This is consistent with the expectation that spin-spin interactions are suppressed in doubly-heavy systems by a factor of $1/(m_c m_b)$ .

Significance and Claims
The paper claims that these results provide QCD-based criteria for assessing the diquark mass parameters used in various constituent diquark models of tetraquarks (e.g., Type I/II models, relativized quark models, and Regge/string models).

Comparison with Models: The authors note that their predictions are generally slightly lower than some existing model parameters (e.g., relativistic quark models) but show excellent agreement with Regge relation and string model results used in specific $T_{bc}$ tetraquark mass determinations.
Implications for $T_{bc}$ : The calculated masses are directly relevant to understanding the $T_{bc}$ tetraquark system ( $bc\bar{u}\bar{d}$ ), particularly regarding its mass relative to the $\bar{B}D$ threshold.
Methodological Advance: The authors emphasize that the numerical techniques developed for loop integrations and sum-rule construction offer a new pathway for performing NLO QCD sum-rule calculations in systems with multiple scales, a task previously hindered by the lack of analytic master integral expressions.

The work concludes that the [bc] diquark possesses a rich quantum number structure, suggesting the potential for a similarly rich spectrum of tetraquark systems containing [bc] diquark content.

More like this