Examination of the $c\bar{c}+n+^{10}$Be bound-state problem within three cluster models based on QCD charmonium-nucleon interactions

Using the hyperspherical harmonics method and effective potentials derived from HAL QCD lattice QCD results, this study predicts that the $c\bar{c}+n+^{10}$Be three-cluster system forms bound states with central binding energies ranging from 1.91 to 3.55 MeV and root-mean-square radii of approximately 2.5 fm.

The Gist

Imagine the atomic nucleus not just as a solid ball of protons and neutrons, but as a tiny, bustling dance floor. Usually, the dancers are the familiar particles that make up ordinary matter. But what happens if you invite a very heavy, exotic guest to the party?

This paper explores a hypothetical scenario where a "charmonium" particle (a heavy pair of quarks, like a tiny, dense weight) joins a specific dance floor made of a Beryllium-10 nucleus and a single neutron. The researchers are asking: Will this heavy guest stick to the dance floor, or will it bounce right off?

Here is a breakdown of their investigation using simple analogies:

1. The Exotic Guest: The "Heavy Quark"

In the world of subatomic physics, most particles are made of "light" ingredients. But this study focuses on charmonium ($c\bar{c}$), which is like a heavy, dense weight made of "charm" quarks. Think of it as a bowling ball in a room full of ping-pong balls. The paper looks at two types of these heavy guests: the $J/\psi$ and the $\eta_c$.

2. The Dance Floor: The Beryllium-10 Nucleus

The "stage" for this experiment is a specific type of atomic nucleus called Beryllium-10, plus one extra neutron.

• The Setup: The researchers treat this system as a three-part team: The heavy guest (charmonium), the extra neutron, and the Beryllium-10 core.
• The Halo Effect: The Beryllium-10 nucleus is described as having a "halo" nature. Imagine a tight core (the Beryllium) with a loose, fuzzy cloud of a neutron orbiting it, like a fuzzy halo around a planet. The heavy guest is expected to interact with this whole fuzzy system.

3. The Invisible Glue: QCD Forces

How does the heavy guest stick to the dance floor?

• The Problem: Usually, particles stick together by swapping lighter particles (like mesons). But because the heavy guest is made of heavy quarks, this usual "glue" is very weak or blocked by rules of physics (called the OZI rule).
• The Solution: The paper suggests the glue comes from QCD van der Waals forces. You can think of this like a very subtle, invisible magnetic pull generated by the exchange of multiple "gluons" (the particles that hold quarks together). It's a weak force, but if it's strong enough, it could hold the heavy guest in place.

4. The Method: The "Folding" Recipe

To figure out if the guest sticks, the researchers had to calculate the strength of this invisible glue.

• Step 1: They started with the most accurate "recipe" available for how a single heavy guest interacts with a single neutron. This recipe comes from supercomputer simulations (Lattice QCD) run by the HAL QCD Collaboration.
• Step 2: Since the dance floor is a whole nucleus (Beryllium-10), not just one neutron, they used a method called single-folding. Imagine taking the "glue" recipe for one neutron and spreading it out over the entire shape of the Beryllium nucleus, averaging it out to see how the whole nucleus feels to the guest.

5. The Results: A Successful "Hug"

Using a sophisticated mathematical tool called the hyperspherical harmonics method (which is like a high-tech way to map the movements of three dancing partners), they solved the equations to see if a stable "bound state" forms.

The findings are positive:

• It Sticks: The calculations show that the heavy guest does get trapped by the Beryllium-10 and the neutron. It forms a stable, bound state.
• How Strong? The "hug" isn't incredibly tight, but it is real.
  • The strongest "hug" (binding energy) is about 4.28 MeV (or 3.55 MeV if you average out the spin details).
  • The weakest "hug" is about 1.91 MeV.
  • Analogy: In the world of nuclear physics, these are small but significant energies, meaning the system is stable enough to exist for a measurable amount of time.
• Size: The resulting "dance trio" is slightly larger than the original nucleus, with a radius of about 2.5 femtometers (a femtometer is one-quadrillionth of a meter).

6. The Big Picture

The paper concludes that while we haven't seen this specific "charmonium-nucleus" system in a lab yet, the math says it should exist. It's a theoretical prediction that the heavy guest can find a comfortable spot inside this specific nuclear arrangement, held there by the subtle, multi-gluon forces of the strong interaction.

The authors note that spotting this in the real world is hard because creating these heavy particles and getting them to stick to a nucleus requires very specific, high-energy conditions, likely to be found in major particle accelerators like those at Jefferson Lab or FAIR. But for now, the math says the party is possible.

Technical Summary

1. Problem Statement

The paper investigates the theoretical possibility of forming a bound state in a hypothetical charmonium–nucleus system composed of a charmonium meson ($c\bar{c}$, specifically $J/\psi$ or $\eta_c$), a neutron ($n$), and a $^{10}\text{Be}$ nucleus.

• Context: While the existence of charmonium-nucleus bound states has been a subject of theoretical interest since Brodsky's proposal regarding QCD van der Waals forces, experimental detection remains elusive.
• Specific Challenge: Previous studies often relied on phenomenological potentials or lattice QCD results at unphysical quark masses. This study aims to utilize state-of-the-art lattice QCD interactions (from the HAL QCD Collaboration) calculated at nearly physical pion masses to determine if a three-body system ($c\bar{c} + n + ^{10}\text{Be}$) can form a bound state.
• Goal: To calculate the binding energies and root-mean-square (RMS) matter radii of this system for various spin channels ($J/\psi N$ and $\eta_c N$) and assess the viability of such "charmed hypernuclei."

2. Methodology

The study employs a rigorous three-cluster model solved using the Hyperspherical Harmonics (HH) method.

A. Interaction Potentials

1. $c\bar{c}$-Nucleon ($c\bar{c}$-N) Interactions:

   • The fundamental inputs are the S-wave potentials between charmonium ($J/\psi, \eta_c$) and nucleons ($N$).
   • These are derived from HAL QCD Collaboration lattice QCD simulations at $m_\pi \approx 146.4$ MeV (nearly physical).
   • Four specific channels are considered:
     • Spin-3/2 $J/\psi N$ ($^4S_{3/2}$)
     • Spin-1/2 $J/\psi N$ ($^2S_{1/2}$)
     • Spin-1/2 $\eta_c N$ ($^2S_{1/2}$)
     • Spin-averaged $J/\psi N$
   • The potentials are parameterized using three-range Gaussian functions fitted to the lattice data.

2. $^{10}\text{Be}$-$c\bar{c}$ Effective Potential:

   • Since $^{10}\text{Be}$ is a tightly bound nucleus, the interaction between the charmonium and the nucleus is approximated using the Single-Folding Procedure (SFP).
   • The potential $U_{^{10}\text{Be}-c\bar{c}}(r)$ is calculated by folding the $c\bar{c}$-N potential with the matter density distribution of $^{10}\text{Be}$ (modeled as a Woods-Saxon distribution).
   • The resulting effective potentials are fitted to a Woods-Saxon (WS) form for computational efficiency in the three-body calculation.

3. Neutron-$^{10}\text{Be}$ ($n$-$^{10}\text{Be}$) Interaction:

   • Modeled using a standard Woods-Saxon potential including central and spin-orbit terms.
   • Parameters are chosen to reproduce the properties of the $^{11}\text{Be}$ halo nucleus, treating $^{10}\text{Be}$ as a core capable of excitation.
   • Phase-equivalent shallow potentials are used to remove Pauli-forbidden states.

B. Three-Body Formalism

• Framework: The system is treated as a three-body problem ($c\bar{c} + n + \text{core}$) using Jacobi coordinates.
• Equations: The total wave function is expanded using Hyperspherical Harmonics (HH). The problem is reduced to solving a set of coupled differential equations for the hyperradial functions $R_\beta(\rho)$, derived from the Faddeev equations.
• Core Assumption: The $^{10}\text{Be}$ core is treated as an explicit degree of freedom with intrinsic structure, while the $c\bar{c}$ and $n$ are treated as clusters relative to this core.

3. Key Contributions

• First Application of Near-Physical Lattice QCD to $c\bar{c} + n + A$ Systems: This is one of the first studies to apply HAL QCD potentials calculated at nearly physical pion masses to a three-body charmonium-nucleus system.
• Spin-Dependence Analysis: The study explicitly investigates how the binding energy depends on the spin configuration of the underlying $c\bar{c}$-N interaction (comparing spin-3/2, spin-1/2, and spin-averaged channels).
• Validation of Cluster Models: It demonstrates the applicability of the hyperspherical harmonics method combined with single-folding potentials for heavy-flavor nuclear systems.

4. Numerical Results

The calculations predict that the $c\bar{c} + n + ^{10}\text{Be}$ system forms bound states across all considered channels.

Interaction Channel	Two-Body Binding ($^{10}\text{Be}$-$c\bar{c}$) [MeV]	Three-Body Binding ($B_3$) [MeV]	RMS Matter Radius [fm]
Spin-3/2 $J/\psi N$	2.18	3.12 (3.47)*	2.49 (2.49)*
Spin-1/2 $J/\psi N$	3.24	4.28 (3.55)*	2.45 (2.48)*
Spin-1/2 $\eta_c N$	1.11	1.91	2.60
Spin-Averaged $J/\psi N$	2.52	3.55	2.48

*Values in parentheses correspond to calculations using the spin-averaged $^{10}\text{Be}$-$J/\psi$ potential derived from Eq. (2), while the primary values use spin-polarized potentials.

• Binding Energies: The system is bound in all cases. The strongest binding is found for the Spin-1/2 $J/\psi N$ channel ($\approx 4.28$ MeV), while the weakest is for the $\eta_c N$ channel ($\approx 1.91$ MeV).
• Radii: The predicted RMS matter radii are approximately 2.45–2.60 fm, indicating a compact but extended structure compared to the $^{10}\text{Be}$ core alone ($2.30$ fm).
• Potential Depth: The effective $^{10}\text{Be}$-$c\bar{c}$ potentials are attractive, with depths ($U_0$) ranging from roughly 9.6 MeV ($\eta_c$) to 14.3 MeV ($J/\psi$ spin-1/2).

5. Significance and Outlook

• Experimental Feasibility: The predicted binding energies (1.9–4.3 MeV) are within the range potentially detectable by current or near-future facilities. The paper suggests that experiments at Jefferson Lab (JLab), J-PARC, FAIR, or via relativistic heavy-ion collisions (RHIC/LHC) could search for these states.
• QCD Insights: The existence of these bound states would provide direct evidence for the attractive nature of QCD van der Waals forces (multi-gluon exchange) in nuclear environments and offer insights into the in-medium modification of charmonium properties.
• Future Directions: The authors note that while the current Faddeev-based cluster model is precise, future work should incorporate:
  • Calculations at the exact physical point by HAL QCD.
  • Three-body forces (beyond binary interactions).
  • Comparison with ab initio methods (e.g., No-Core Shell Model) and experimental data as it becomes available.

In conclusion, the paper provides a robust theoretical prediction that charmed hypernuclei of the type $c\bar{c} + n + ^{10}\text{Be}$ are likely to exist, with binding energies sufficient to warrant dedicated experimental searches.