Study of the $Ω_{ccc}Ω_{ccc}$ and $Ω_{bbb}Ω_{bbb}$ dibaryons in QCD Sum Rules

Using QCD sum rules with the iterative dispersion relation method to handle complex loop diagrams, this study predicts that scalar $\Omega_{ccc}\Omega_{ccc}$ and $\Omega_{bbb}\Omega_{bbb}$ dibaryons lie lower in energy than their tensor counterparts, with the former appearing as a resonance just above the threshold and the latter potentially forming bound states.

The Gist

The Big Picture: Building "Super-Heavy" Lego Castles

Imagine the universe is built out of tiny, fundamental Lego bricks called quarks. Usually, these bricks snap together in small groups:

• 3 bricks make a baryon (like a proton or neutron, which are the building blocks of our atoms).
• 4 bricks make a tetraquark (a rare, exotic shape recently spotted by scientists).

But what happens if you try to snap 6 bricks together? That's a dibaryon. It's like trying to build a massive, double-sized castle out of Lego.

For decades, scientists have only confirmed one stable "double castle" in nature: the deuteron (two protons/neutrons stuck together). But recently, experiments found a whole family of "super-heavy" 4-brick structures made entirely of charm quarks. This discovery made physicists wonder: Could there be even bigger, super-heavy 6-brick castles made entirely of charm or bottom quarks?

This paper is a theoretical investigation to see if these "Super-Heavy Castles" (specifically called $\Omega_{ccc}\Omega_{ccc}$ and $\Omega_{bbb}\Omega_{bbb}$) can actually exist without falling apart.

---

The Characters: The "Heavyweights"

To understand the study, you need to know the players:

1. Charm Quarks ($c$): These are heavy, like a bowling ball.
2. Bottom Quarks ($b$): These are super heavy, like a grand piano.
3. The $\Omega$ Baryon: Think of this as a specific type of Lego structure made of three heavy quarks.
   • $\Omega_{ccc}$ is a structure made of three charm quarks.
   • $\Omega_{bbb}$ is a structure made of three bottom quarks.

The paper asks: If you take two of these heavy structures and try to stick them together, do they form a stable bond, or do they bounce off each other?

The Method: The "Quantum Calculator"

The authors didn't build these in a lab (it's too expensive and hard to make). Instead, they used a mathematical tool called QCD Sum Rules.

Think of this like trying to guess the weight of a sealed, heavy box without opening it. You can't see inside, but you can:

1. Shake the box (mathematical equations).
2. Listen to how it rattles (spectral density).
3. Use the laws of physics (Quantum Chromodynamics) to calculate what must be inside based on how it moves.

The authors used a very specific type of "mathematical Lego" (called an interpolating current) to represent these 6-quark structures. They then ran a massive calculation to see if the math predicts a stable, heavy object.

The Technical Hurdle: The "Five-Loop Banana"

The math involved here is notoriously difficult. The authors had to calculate something called a "five-loop banana diagram."

• The Analogy: Imagine trying to calculate the path of a ball bouncing through a maze. Now imagine the ball is bouncing through a maze that has five layers of walls, and every time it hits a wall, it splits into more paths.
• The Problem: Standard math tools get stuck in these mazes. They often hit a "divergence," which is like a calculator trying to divide by zero and crashing.
• The Solution: The authors used a clever new trick called the Iterative Dispersion Relation (IDR) method. Instead of trying to solve the whole maze at once, they solved it layer by layer, like peeling an onion. They also fixed a specific "glitch" (the small-circle divergence) that usually breaks these calculations, ensuring their results are accurate.

The Results: The "Heavy" Verdict

After running their super-complex calculations, here is what they found:

1. The Shape Matters (Scalar vs. Tensor)

The "castles" can be built in two different shapes:

• Scalar (The Flat Stack): Like stacking bricks flat.
• Tensor (The Tangled Knot): Like twisting the bricks together.
• Finding: In both the charm and bottom worlds, the Flat Stack (Scalar) is always lighter and more stable than the Tangled Knot (Tensor). Nature prefers the simpler shape.

2. The Charm Castle ($\Omega_{ccc}\Omega_{ccc}$)

• The Result: The mass of this double-charm castle is calculated to be slightly heavier than two separate charm castles floating apart.
• The Analogy: Imagine two magnets. If you push them together, they repel slightly. They are almost stuck, but not quite.
• Conclusion: It likely does not form a bound state. It's a "resonance" that might exist for a split second before flying apart.

3. The Bottom Castle ($\Omega_{bbb}\Omega_{bbb}$)

• The Result: The mass of this double-bottom castle is significantly lighter than two separate bottom castles.
• The Analogy: Imagine two magnets that snap together with a loud click. They are stuck tight.
• Conclusion: This does form a bound state. It is a stable particle that could theoretically exist in the universe.

Why Does This Matter?

1. Testing the Rules of Physics: Since these particles are so heavy, we can ignore some of the messy "relativistic" effects (like time slowing down) that complicate lighter particles. This gives us a "clean" test of how the Strong Force (the glue holding atoms together) works.
2. No Light Quarks: These particles are made only of heavy quarks. There are no light quarks (like up or down) to mess things up by swapping around. It's a pristine environment to study how heavy things interact.
3. Future Discovery: If the Large Hadron Collider (LHC) or future colliders start smashing particles together, they might finally spot these "Bottom Castles." This paper tells experimentalists exactly what mass to look for (around 26.6 GeV for the bottom one).

Summary in One Sentence

Using advanced math to solve a "five-layered maze," this paper predicts that while a double-charm particle is too heavy to stick together, a double-bottom particle is light enough to form a stable, new kind of "super-atom" waiting to be discovered.

Technical Summary

1. Problem Statement

The paper addresses the theoretical investigation of fully-heavy dibaryons, specifically the $\Omega_{ccc}\Omega_{ccc}$ (six charm quarks) and $\Omega_{bbb}\Omega_{bbb}$ (six bottom quarks) systems.

• Context: Following the experimental observation of fully-charm tetraquarks ($X(6600)$, $X(6900)$, etc.) by LHCb, ATLAS, and CMS, there is growing interest in whether bound states of two heavy baryons (dibaryons) exist.
• Theoretical Challenge: Previous studies using Lattice QCD, Quark Models, and One-Boson Exchange (OBE) models have yielded conflicting results regarding the existence and binding energy of these states. Some predict bound states, while others suggest they are unbound or only loosely bound.
• Specific Technical Hurdle: Calculating the correlation functions for these six-quark systems involves massive five-loop Feynman diagrams (banana diagrams). Standard techniques like Integration-by-Parts (IBP) reduction become computationally prohibitive for such high-loop orders. Furthermore, the inclusion of non-perturbative gluon condensate terms introduces a "small-circle divergence" problem in the dispersion relations, which has historically been mishandled in similar calculations, leading to inaccurate results.

2. Methodology

The authors employ QCD Sum Rules to analyze the scalar ($J^P = 0^+$) and tensor ($J^P = 2^+$) channels of the $\Omega_{ccc}\Omega_{ccc}$ and $\Omega_{bbb}\Omega_{bbb}$ systems.

• Interpolating Current: They construct a symmetric tensor interpolating current in a hadronic molecular configuration:
  $$J_{\mu\nu}(x) = \epsilon_{i_1j_1k_1}\epsilon_{i_2j_2k_2}[Q^T_{i_1}(x)C\gamma_\mu Q_{j_1}(x)]Q^T_{k_1}(x) \cdot C\gamma_5 \cdot Q_{k_2}(x)[Q^T_{i_2}(x)C\gamma_\nu Q_{j_2}(x)]$$
  This current couples to both scalar and tensor dibaryon states.
• Operator Product Expansion (OPE): The correlation function is calculated up to dimension-four gluon condensates ($\langle g_s^2 G^2 \rangle$). Dimension-six triple-gluon condensates are neglected as they are negligible in fully-heavy systems.
• Computational Techniques:
  • Symmetry Reduction: To handle the 720 Wick contraction terms, the authors implemented a simplification algorithm based on relabeling propagators, color factors, and Dirac matrices, reducing the final distinct terms to 36.
  • Iterative Dispersion Relation (IDR): To compute the massive five-loop banana diagrams, the authors utilize the IDR method. This approach constructs multi-loop integrals recursively by treating the dispersion factor $(s-q^2)^{-1}$ as a propagator for the next loop, avoiding the intractable IBP reduction.
  • Handling Small-Circle Divergence: A key innovation is the rigorous treatment of the small-circle divergence arising from high-power propagators in the gluon condensate terms. Instead of using standard dispersion relations (which fail here) or complex Generalized Dispersion Relations (GDR), they adopt a strategy (referenced from Ref. [47]) that avoids GDR calculations entirely while ensuring numerical stability and correctness.
• Dimensionless Formulation: The spectral densities are reduced to dimensionless forms by factoring out the heavy quark mass ($m_Q$). This allows the same OPE series to be used for both charm and bottom systems, with differences arising only from the explicit value of $m_Q$.

3. Key Contributions

1. Advanced Computational Framework: The successful application of the IDR method to massive five-loop diagrams in the context of QCD sum rules, significantly improving computational efficiency over traditional loop-by-loop Feynman parameterization.
2. Resolution of Divergence Issues: A robust solution to the small-circle divergence problem in non-perturbative terms, correcting previous inaccuracies found in similar literature.
3. Systematic Comparison: A unified analysis of both charm ($\Omega_{ccc}\Omega_{ccc}$) and bottom ($\Omega_{bbb}\Omega_{bbb}$) sectors within the same formalism, allowing for a direct comparison of binding trends across heavy quark flavors.
4. Spin Dependence: Simultaneous investigation of scalar ($0^+$) and tensor ($2^+$) channels to determine spin-dependent mass splitting.

4. Results

The numerical analysis yields the following mass predictions (using $\overline{MS}$ quark masses):

System	$J^P$	Mass (GeV)	Status vs. Threshold
$\Omega_{ccc}\Omega_{ccc}$	$0^+$ (Scalar)	$9.77 \pm 0.04$	Slightly above $2m_{\Omega_{ccc}}$ (Unbound/Resonance)
$\Omega_{ccc}\Omega_{ccc}$	$2^+$ (Tensor)	$10.45 \pm 0.04$	Above threshold
$\Omega_{bbb}\Omega_{bbb}$	$0^+$ (Scalar)	$26.60 \pm 0.05$	Below $2m_{\Omega_{bbb}}$ (Bound State)
$\Omega_{bbb}\Omega_{bbb}$	$2^+$ (Tensor)	$27.05 \pm 0.06$	Below threshold

• Spin Ordering: For both charm and bottom systems, the scalar dibaryon lies lower in mass than its tensor counterpart.
• Binding Nature:
  • The $\Omega_{ccc}\Omega_{ccc}$ scalar state is predicted to be slightly above the two-baryon threshold, suggesting it does not form a stable bound state but may exist as a resonance.
  • The $\Omega_{bbb}\Omega_{bbb}$ system is predicted to form a bound state, as its mass is significantly lower than the $2\Omega_{bbb}$ threshold.
• Comparison with Literature: The scalar $\Omega_{ccc}\Omega_{ccc}$ mass aligns well with some quark model predictions but differs from others. The $\Omega_{bbb}\Omega_{bbb}$ mass is notably lower than many previous theoretical predictions, reinforcing the likelihood of a bound state in the bottom sector.

5. Significance

• Validation of Heavy Quark Dynamics: The results support the hypothesis that fully-heavy systems, free from light quark exchanges, are governed primarily by heavy quarkonia interactions and Coulomb-like forces, making them ideal testing grounds for QCD.
• Experimental Guidance: The prediction of a bound $\Omega_{bbb}\Omega_{bbb}$ state provides a concrete target for future high-energy experiments (e.g., LHCb, FCC) to search for fully-heavy dibaryons.
• Methodological Benchmark: The paper establishes a reliable computational pipeline for high-loop QCD sum rule calculations involving heavy quarks, resolving long-standing technical issues regarding multi-loop integrals and dispersion relation divergences. This methodology can be extended to other exotic hadronic states.