A Data-Guided Coalescence Model for Light Nuclei and Hypernuclei Production in Relativistic Heavy-Ion Collisions at $\sqrt{s_{\rm{NN}}} = 3$--200 GeV

This paper presents a data-guided coalescence model that successfully predicts light nuclei and hypernuclei production across a wide energy range in relativistic heavy-ion collisions, revealing that hypertriton yields are highly sensitive to wave function assumptions, particularly at low energies and in low-multiplicity environments.

The Gist

The Big Picture: Solving the "Hyperon Puzzle"

Imagine the inside of a neutron star as a crowded, high-pressure dance floor. Usually, this floor is packed with neutrons. But physicists think that if the pressure gets high enough, some neutrons should turn into "hyperons" (a strange type of particle).

Here's the problem: If too many hyperons show up, they act like a soft pillow, making the whole dance floor squishy. Physics says this should cause the star to collapse under its own weight. But astronomers see massive neutron stars that are not collapsing. They are stiff and strong. This contradiction is called the Hyperon Puzzle.

To solve it, we need to understand how these hyperons interact with normal particles. Do they push each other away (stiffening the star) or pull together (softening it)?

The Experiment: A Cosmic "Clay Ball" Maker

Since we can't build a neutron star in a lab, scientists smash heavy atoms (like Gold) together at nearly the speed of light. This creates a tiny, super-hot fireball that mimics the conditions inside a star for a split second.

When the fireball cools down, the particles inside try to stick together to form new "clay balls" (nuclei). Some of these balls contain hyperons. These are called hypernuclei. The specific one this paper studies is the Hypertriton (a tiny ball made of a proton, a neutron, and a hyperon).

The Method: The "Coalescence" Model

The paper uses a method called Coalescence. Think of it like this:

Imagine a giant, chaotic dance floor (the collision) where people (protons, neutrons, and hyperons) are running around. To form a "clay ball" (a nucleus), three people need to run close enough to each other and hold hands.

The paper asks: How likely are they to hold hands?

To answer this, the authors developed a "Data-Guided" approach. Instead of guessing how the dance floor looks, they looked at the actual footprints left behind (experimental data) to figure out the size of the dance floor and how crowded it was.

1. Measuring the Room: They first looked at simple pairs (Deuterons: a proton and neutron holding hands) to measure the size of the "dance floor" (the source size).
2. Predicting the Complex: Once they knew the size of the room, they tried to predict how often three people would hold hands to form a Hypertriton.

The Twist: The Shape of the "Hug"

The most important discovery in this paper is about the shape of the hug.

In physics, particles don't just sit at a fixed distance; they have a "wave function," which is like a probability cloud showing where they are likely to be.

• The Old Guess: Scientists used to guess that the Hypertriton was a very loose, fluffy cloud (like a giant, fluffy dandelion).
• The New Test: The authors tested different shapes: some fluffy, some compact, and some based on specific theories (the "Congleton" models).

The Result:
The shape of the hug matters a lot, especially when the "dance floor" is small.

• At High Energies (Big Dance Floor): The room is huge. Whether the hug is tight or loose, the particles can find each other easily. All the models look similar.
• At Low Energies (Small Dance Floor): The room is tiny. If the hug is too loose (fluffy), the particles can't reach each other in the small space, and they don't stick together. If the hug is tight (compact), they stick easily.

The Analogy: Imagine trying to fit three people into a tiny elevator.

• If they are holding hands loosely with long arms (fluffy wave function), they won't fit.
• If they hug tightly (compact wave function), they fit perfectly.

The paper found that the Hypertriton seems to have a tighter hug than previously thought. The "fluffy" models predicted too few hypernuclei, while the "tighter" models matched the real data much better.

Why This Matters

1. Solving the Puzzle: By figuring out the exact shape of the Hypertriton's hug, scientists can better calculate how hyperons behave in the dense core of a neutron star. This helps explain why those stars don't collapse.
2. The Best Place to Look: The paper suggests that to see these differences clearly, we shouldn't look at the biggest, hottest collisions. We should look at smaller, lower-energy collisions (or collisions that aren't perfectly centered). In these "smaller rooms," the shape of the hug is the deciding factor.
3. Future Experiments: This gives a roadmap for future experiments (like those at the STAR detector or new facilities in China and Germany). They know exactly what to measure to finally pin down the rules of hyperon interactions.

Summary in One Sentence

This paper uses a "data-driven" method to show that the tiny, strange atoms formed in particle collisions are sensitive to the "shape" of their internal structure, proving that studying these particles in smaller, lower-energy collisions is the key to solving the mystery of how neutron stars stay strong.

Technical Summary

1. Problem Statement

The paper addresses the hyperon puzzle in nuclear astrophysics. Theoretical models predict that hyperons (baryons containing strange quarks) should appear in the dense cores of neutron stars, softening the Equation of State (EoS) and reducing the maximum predicted neutron star mass. This contradicts observations of massive neutron stars ($\sim 2 M_\odot$).

• The Gap: Resolving this puzzle requires precise knowledge of hyperon–nucleon (YN) and three-body hyperon–nucleon–nucleon (YNN) interactions.
• The Challenge: Experimental data on these interactions is scarce, particularly at the high densities relevant to neutron stars. While hypernuclear binding energies provide low-density constraints, they do not fully constrain the wave functions or short-range structures necessary for extrapolation to high densities.
• The Opportunity: Relativistic heavy-ion collisions produce light hypernuclei (like the hypertriton, $^3_\Lambda\text{H}$). The production yield of these nuclei is sensitive to their internal wave functions. However, traditional theoretical approaches (thermal models or transport models with coalescence afterburners) often suffer from large uncertainties regarding the source size and phase-space distributions of nucleons, especially at low collision energies.

2. Methodology: A Data-Guided Coalescence Framework

The authors propose a data-driven coalescence approach to minimize model dependencies on the EoS and transport dynamics.

• Core Concept: Instead of relying on transport models to generate nucleon/hyperon spectra, the model uses experimental data (proton and deuteron yields) to extract the source size ($R_{inv}$) and then uses this source size to predict the yields of heavier nuclei ($A=3$) and hypernuclei.
• Source Size Extraction:
  • The coalescence parameter $B_2$ for deuterons is calculated from experimental proton and deuteron spectra using the relation:
    $$B_2 \propto \int |\Phi_d(r)|^2 S_2(r) d^3r$$
    where $\Phi_d$ is the deuteron wave function and $S_2$ is the source function.
  • By assuming a Gaussian source, $R_{inv}$ is extracted from the measured $B_2$.
  • Validation: The method was validated against ALICE data at $\sqrt{s_{NN}} = 5.02$ TeV, showing consistency with femtoscopic correlation measurements (within a $\sim 8\%$ systematic uncertainty).
• Wave Functions:
  • Deuteron: Argonne v18 (default), Hulthen, and Gaussian forms.
  • Triton ($^3$H) and $^3$He: Gaussian wave functions.
  • Hypertriton ($^3_\Lambda\text{H}$): Four different parameterizations were tested:
    1. Congleton (a): Based on 3-body EFT ($Q_\Lambda = 2.5$ fm$^{-1}$).
    2. Congleton (b): 2-body closure approximation ($Q_\Lambda = 1.17$ fm$^{-1}$).
    3. Congleton (c): Alternative parameters near the $\langle r_{d\Lambda} \rangle = 10.9$ fm contour.
    4. Gaussian: A simple Gaussian ansatz.
  • All tested wave functions share a similar Root-Mean-Square (RMS) separation $\langle r_{d\Lambda} \rangle \approx 10.8$ fm, allowing the study of shape sensitivity independent of size.
• Inputs: The model utilizes measured proton and $\Lambda$ spectra (corrected for feed-down) and extrapolates low-$p_T$ regions using blast-wave parameterizations.

3. Key Contributions

1. Novel Framework: Development of a coalescence model that bypasses the uncertainties of transport models by directly constraining the source size via experimental deuteron yields.
2. Systematic Wave Function Study: A comprehensive analysis of how different $^3_\Lambda\text{H}$ wave function shapes affect production yields, keeping the RMS separation constant.
3. Energy Range Coverage: Application of the model across a wide energy range ($\sqrt{s_{NN}} = 3–200$ GeV), covering the Beam Energy Scan (BES) region where the EoS is most uncertain.
4. Multiplicity Scaling: Demonstration that the extracted source radius scales universally with charged-particle multiplicity ($R_{inv} \propto \langle dN_{ch}/d\eta \rangle^{1/3}$), independent of collision energy or system (Au+Au vs. Pb+Pb).

4. Key Results

A. Light Nuclei (Triton and $^3$He)

• Spectra and Yields: The coalescence model using realistic wave functions (Argonne v18 or Hulthen) successfully reproduces experimental triton $p_T$ spectra and yield ratios ($N_t/N_p$) across all centralities and energies.
• Comparison with Models:
  • Thermal Model: Overestimates the $N_t/N_p$ ratio by a factor of $\sim 2$.
  • Blast-Wave Model: Fails to describe the mean transverse momentum ($\langle p_T \rangle$) of tritons, suggesting tritons do not share the same kinetic freeze-out surface as hadrons.
• Wave Function Sensitivity: The choice of deuteron wave function (Argonne vs. Gaussian) significantly impacts predictions, particularly in peripheral collisions. Realistic wave functions are essential for accuracy.

B. Hypertriton ($^3_\Lambda\text{H}$)

• Wave Function Sensitivity: The predicted $^3_\Lambda\text{H}$ yield is highly sensitive to the specific shape of the wave function, even when the RMS separation is fixed.
  • Congleton (a) (most compact in the core) predicts the highest yields.
  • Gaussian (broadest) predicts the lowest yields.
• Energy and Centrality Dependence:
  • Sensitivity to the wave function is strongest at low collision energies ($\sqrt{s_{NN}} = 3$ GeV) and in peripheral collisions (small source sizes).
  • In small systems, the overlap between the source and the wave function varies drastically depending on the wave function's behavior at short distances ($r < 4.5$ fm).
• Comparison with Data:
  • At $\sqrt{s_{NN}} = 3$ GeV, the Congleton (b) and (c) parameterizations agree well with STAR data.
  • The Gaussian wave function tends to underestimate the data, suggesting it is too broad for the actual hypertriton structure.
• Strangeness Population Factor ($S_3$): The ratio $S_3 = (N_{^3_\Lambda\text{H}}/N_\Lambda) / (N_{^3\text{He}}/N_p)$ shows clear dependence on the wave function at low multiplicities. High-energy small systems (Ru+Ru, Zr+Zr at 200 GeV) also offer a promising avenue to probe this structure.

5. Significance and Outlook

• Resolving the Hyperon Puzzle: By constraining the hypertriton wave function through production yields, this work provides indirect constraints on the YN and YNN interactions. This is crucial for understanding the repulsive three-body forces needed to support massive neutron stars.
• Experimental Strategy: The paper establishes that low-energy collisions and small collision systems (peripheral or small nuclei) are the most sensitive probes for hypernuclear structure due to the smaller source size enhancing wave-function overlap effects.
• Future Directions:
  • The framework can be extended to heavier hypernuclei ($^4_\Lambda\text{H}$, $^4_\Lambda\text{He}$).
  • Future precise measurements from the STAR BES-II program and future facilities (FAIR, HIAF) will be critical to distinguish between the tested wave function parameterizations and refine the hypernuclear structure models.

In summary, the paper demonstrates that a data-guided coalescence model is a powerful tool for extracting nuclear structure information from heavy-ion collision data, specifically highlighting the hypertriton as a sensitive probe for the underlying hyperon-nucleon interactions.