Hidden-charm $uds\,c\bar c$ pentaquarks as flavor eigenstates in a constituent quark model

Using a diffusion Monte Carlo algorithm within a constituent quark model, this study demonstrates that explaining the observed $P_{cs}(4338)$ and $P_{cs}(4459)$ pentaquarks requires the total wavefunction to be an eigenvector of the SU(3) flavor operator, a condition that yields two mass-compatible structures and predicts additional states below the $J/\psi\Lambda$ threshold, whereas imposing only isospin symmetry results in a single structure.

The Gist

Imagine the universe is built out of tiny, invisible LEGO bricks called quarks. Usually, these bricks snap together in two standard ways:

• Mesons: Two bricks stuck together (a quark and an anti-quark).
• Baryons: Three bricks stuck together (like the protons and neutrons in your body).

But for decades, physicists have wondered: What if we tried to snap five bricks together? These hypothetical five-brick structures are called pentaquarks.

Recently, a giant particle collider (LHCb) found evidence of two strange, heavy pentaquarks containing a "hidden charm" (a charm quark and its anti-particle). They named them Pcs(4338) and Pcs(4459). The big mystery was: How are these five bricks actually arranged inside? Are they a tight ball? A loose molecule? A specific pattern?

This paper is like a digital simulation lab where three physicists tried to figure out the answer. Here is how they did it, explained simply:

1. The Simulation: A "Diffusion" Game

The authors used a super-smart computer algorithm called Diffusion Monte Carlo (DMC).

• The Analogy: Imagine you are trying to find the deepest valley in a foggy mountain range. You can't see the whole map, so you drop thousands of "walkers" (virtual particles) that bounce around. Over time, they naturally "diffuse" and settle into the lowest point (the ground state energy).
• The Goal: They wanted to see where these five quarks would settle to find the most stable, lowest-energy arrangement.

2. The Rules of the Game

To make the simulation work, they had to follow the strict rules of quantum mechanics:

• The Bricks: They used up, down, strange, charm, and anti-charm quarks.
• The Glue: They used a specific "glue" formula (the AL1 potential) that mimics how quarks attract and repel each other.
• The Identity Crisis: This is the most important part. In quantum mechanics, identical particles (like the up, down, and strange quarks) are like identical twins. You can't tell them apart. If you swap them, the whole system has to look the same (or exactly opposite).

3. The Big Discovery: The "Flavor" Key

The researchers tried two different ways to set up the rules for these twins:

Attempt A: The "Loose" Rule (Only Isospin)
They treated the quarks as if only their electric charge mattered (Isospin).

• The Result: The simulation only found one stable pentaquark.
• The Problem: This one pentaquark was too light. It would fall apart before it could be seen in the experiments. It didn't match the two real-world discoveries (Pcs(4338) and Pcs(4459)).
• The Lesson: This approach was too simple. It was like trying to sort a deck of cards by only looking at the red suits, ignoring the numbers.

Attempt B: The "Strict" Rule (Full Flavor Symmetry)
They realized they needed to treat the quarks as a complete family (SU(3) flavor symmetry). They forced the simulation to respect the deep symmetry between the up, down, and strange quarks.

• The Result: Suddenly, the simulation found two distinct, stable structures!
  1. Structure 1: A heavy, compact ball of quarks. Its mass matched the Pcs(4459) discovery.
  2. Structure 2: A slightly lighter, but still compact ball. Its mass matched the Pcs(4338) discovery.
• The "Aha!" Moment: By respecting the full "flavor" rules, the math naturally split into two different solutions that perfectly matched the two real-world particles.

4. What Do They Look Like?

The simulation also gave them a "snapshot" of the internal structure:

• Not Loose Molecules: They expected these pentaquarks to be like a loose molecule (a baryon floating next to a meson, like two magnets barely touching).
• The Reality: Instead, the simulation showed they are tight, compact balls. All five quarks are huddled close together, like a tightly packed group hug, rather than two separate groups floating near each other.
  • One looks a bit like a distorted Lambda baryon hugging a charm-anticharm pair.
  • The other looks like a tight cluster of three quarks hugging a charm quark, which is hugging an anti-charm quark.

5. The Hidden Twins

The simulation also predicted two other pentaquarks that are slightly lighter than the ones we found.

• Why haven't we seen them? They are too light to break apart into the "J/psi + Lambda" pieces that the detectors are currently looking for. They would likely decay into a different, harder-to-see combination (eta-c + Lambda).
• The Metaphor: It's like finding two bright stars in the sky, but the simulation says there are two dimmer stars right next to them that are too faint for our current telescopes to see.

Summary

The paper tells us that to understand these exotic five-quark particles, we can't just look at them simply. We have to respect the complex "family rules" (flavor symmetry) of the quarks. When we do that, the math reveals two distinct, compact structures that perfectly explain the two mysterious particles LHCb found. It turns out nature didn't just make one type of pentaquark; it made two different "flavors" of the same family, and we needed the right mathematical lens to see them both.

Technical Summary

1. Problem Statement

The paper addresses the theoretical description of hidden-charm pentaquarks with the quark content $uds\bar{c}c$. Specifically, it aims to explain the existence and properties of two recently observed experimental candidates:

• $P_{cs}(4338)$: Observed in $B^- \to J/\psi \Lambda \bar{p}$, with $J^P = 1/2^-$ and $I=0$.
• $P_{cs}(4459)$: Observed in $\Xi_b^- \to J/\psi \Lambda K^-$, with unknown spin-parity but favored $I=0$.

Previous theoretical works often treated the strange quark ($s$) as distinguishable from the up ($u$) and down ($d$) quarks at the wavefunction level, effectively restricting the flavor symmetry to isospin ($I$) alone. The authors hypothesize that this restriction may be insufficient to reproduce the experimental spectrum, which shows two distinct structures. They propose that the total wavefunction must be an eigenvector of the full $SU(3)$ flavor operator, not just the isospin operator, to correctly model these states.

2. Methodology

The authors employ a non-relativistic constituent quark model solved via the Diffusion Monte Carlo (DMC) algorithm.

• Hamiltonian: The system is governed by the AL1 potential, which includes:
  • A Coulomb-like term (modeling one-gluon exchange).
  • A linear confining potential.
  • A hyperfine spin-spin interaction.
  • Parameters are fixed based on fits to known baryon properties and are not re-tuned for the pentaquarks.
• DMC Technique:
  • The Schrödinger equation is solved in imaginary time to project out the ground state.
  • The method uses "importance sampling" with a trial wavefunction ($\Psi_T$) that includes spatial, spin, color, and flavor degrees of freedom.
  • The trial function is constructed as a linear combination of spin-color-flavor basis states multiplied by a Jastrow-type spatial ansatz (symmetric under particle exchange for $L=0$).
• Symmetry Constraints:
  • Isospin: Fixed at $I=0$ to match experimental observations.
  • Flavor: The authors test two scenarios:
    1. Restricted Flavor: Only $I=0$ is imposed (treating $s$ as distinct).
    2. Full $SU(3)$ Flavor: The wavefunction is required to be an eigenvector of the flavor Casimir operator ($F^2$).
  • Antisymmetry: The wavefunction must be antisymmetric under the exchange of identical fermions ($u, d, s$). The authors construct two distinct symmetry-constrained families of wavefunctions:
    • Family I: Treats $u, d, s$ on equal footing (full antisymmetry among light quarks).
    • Family II: Treats two light quarks symmetrically while correlating them with the $c$ quark ($qqc$) and the third with $\bar{c}$ ($q'\bar{c}$).

3. Key Contributions

• Flavor Symmetry Requirement: The paper demonstrates that reproducing the two experimental pentaquark signals ($P_{cs}(4338)$ and $P_{cs}(4459)$) requires the wavefunction to be an eigenstate of the $SU(3)$ flavor operator. If only isospin ($I=0$) is imposed, the model yields only a single state, failing to explain the two observed structures.
• Four Distinct Configurations: By combining the two symmetry constructions (Family I and II) with the two flavor sectors ($F^2=0$ singlet and $F^2=3$ octet), the authors identify four distinct pentaquark configurations ($Ia, Ib, IIa, IIb$).
• Internal Structure Analysis: The study provides a detailed characterization of the spatial structure (interquark distances) of these states, distinguishing between compact multiquark configurations and loosely bound molecular states.

4. Results

Mass Spectrum

Using the DMC algorithm, the authors calculated the masses of the four configurations:

State	Flavor Symmetry ($F^2$)	Mass (MeV)	Relation to Thresholds
**$Ia$** |$0$ (Singlet)	$4473 \pm 5$	Above $J/\psi \Lambda$ ($4257$ MeV)
**$Ib$** |$3$ (Octet)	$4237 \pm 3$	Below $J/\psi \Lambda$, Above $\eta_c \Lambda$ ($4161$ MeV)
**$IIa$** |$3$ (Octet)	$4350 \pm 6$	Above $J/\psi \Lambda$
**$IIb$** |$0$ (Singlet)	$4245 \pm 7$	Below $J/\psi \Lambda$, Above $\eta_c \Lambda$

• Experimental Matching:
  • **$Ia$** ($4473$ MeV) is compatible with $P_{cs}(4459)$.
  • **$IIa$** ($4350$ MeV) is compatible with $P_{cs}(4338)$.
  • Both states lie above the $J/\psi \Lambda$ threshold, allowing decay into the observed channel.
• Unobserved States:
  • States $Ib$ and $IIb$ lie below the $J/\psi \Lambda$ threshold but above the $\eta_c \Lambda$ threshold. The authors predict these would decay predominantly into $\eta_c \Lambda$, a channel not yet experimentally explored.
• Isospin-Only Scenario:
  • If the flavor symmetry is broken (only $I=0$ imposed), the model yields a single state ($III$) with mass $4200 \pm 4$ MeV. This state is below the $J/\psi \Lambda$ threshold and cannot explain the experimental data, proving the necessity of the full $SU(3)$ flavor treatment.

Internal Structure

Analysis of mean interquark distances ($\langle r_{ij} \rangle$) reveals:

• $Ia$ and $IIa$ (The Candidates): Both are compact structures.
  • $Ia$ resembles a distorted $\Lambda + c\bar{c}$ configuration where the $c\bar{c}$ separation is significantly larger than in charmonium.
  • $IIa$ resembles a compact $qqc + q'\bar{c}$ arrangement (baryon-meson like but tightly bound).
  • Neither shows the large spatial separation characteristic of loosely bound hadronic molecules.
• $Ib$ and $IIb$: These states exhibit larger separations, resembling a baryon plus a $\bar{c}c$ meson with constituents fairly separated.

5. Significance

• Theoretical Necessity of Flavor Symmetry: The work establishes that treating the strange quark as merely distinct from $u/d$ (isospin only) is insufficient for describing hidden-charm pentaquarks. The $SU(3)$ flavor symmetry must be explicitly enforced in the wavefunction to generate the multiplicity of states observed experimentally.
• Nature of Pentaquarks: The results suggest that the observed $P_{cs}$ states are likely compact multiquark states rather than simple hadronic molecules, as indicated by the small interquark distances and the specific mass ordering derived from the constituent quark model.
• Predictive Power: The model predicts the existence of two additional states ($Ib, IIb$) that should be searched for in the $\eta_c \Lambda$ decay channel, providing a clear target for future experimental verification.
• Methodological Robustness: The successful application of DMC with full spin-color-flavor antisymmetrization demonstrates the capability of this technique to handle complex multiquark systems with high precision (estimated theoretical uncertainty of 10-20%).