Sequential Clusterization of Light Nuclei and Hypernuclei in Heavy-Ion Collisions within a Wigner Function Coalescence Framework

This paper investigates the formation of light nuclei and hypernuclei in Au+Au collisions at $\sqrt{s_{NN}}=3~\mathrm{GeV}$ using a parameter-free coalescence framework based on realistic $N$-body wave functions, revealing species-dependent formation times and improving the description of $A=4$ yields through additional cluster-nucleon channels while providing predictions for heavier hypernuclei.

The Gist

Imagine a high-energy particle collision as a massive, chaotic dance floor where thousands of tiny particles (protons and neutrons) are spinning, bumping, and flying apart at incredible speeds. The scientists in this paper wanted to understand how, amidst this chaos, these tiny particles sometimes stick together to form "dance couples" or even small "groups" (like light nuclei and hypernuclei).

Here is a simple breakdown of what they did and what they found, using everyday analogies:

The Setup: A High-Speed Dance Floor

The researchers simulated a collision between two heavy gold atoms (Au+Au) at a specific energy level. Think of this as two crowds of people rushing into a room and colliding. For a split second, it's a hot, dense mess. Then, the crowd expands and cools down.

Usually, scientists assume that these particles only stick together to form groups at the very end of the dance, when the music stops and everyone freezes in place. This is called "kinetic freeze-out."

The New Tool: A Better Blueprint

In the past, scientists used a rough, generic blueprint to guess how these groups form. It was like assuming every dance group has a perfect, tight circle shape. But the paper argues that some groups are actually loose and floppy (like a stretched-out rubber band), and the old blueprint didn't fit them well.

Instead, the authors used a realistic, custom-made blueprint for each group. They solved complex math equations to get the exact shape and size of these particle groups. This allowed them to see the groups exactly as they are, without guessing.

The Big Discovery: Timing is Everything

The most exciting finding is about when these groups form. The researchers tested different "stop times" for the dance floor to see when the groups were most likely to stick together.

• The Small Groups (Deuterons, Tritons, Helium-3): These are like small pairs or trios. The paper found they form late in the process, when the crowd has already spread out and thinned. They need the space to find each other and settle down.
• The Big Groups (Helium-4 and Hypernuclei): These are larger, tighter groups. Surprisingly, the paper found they form much earlier, while the crowd is still very dense and crowded.

The Analogy: Imagine trying to form a huddle.

• If you are a small group of 2 or 3 people, you can wait until the crowd thins out and find your friends easily.
• If you are a large group of 4 people who need to hold hands tightly, you have to grab each other immediately while the crowd is still packed tight. If you wait until the crowd disperses, it's too hard to get all four of you close enough at the same time.

The "Side Door" Effect

The paper also discovered that for the larger groups (like Helium-4), there isn't just one way to form. Sometimes, a smaller group (like a trio) grabs one extra person to become a larger group. The authors found that including these "side door" formation paths was crucial. Without them, their models couldn't explain how many of these large groups were actually being created in the experiments.

The Results: Matching the Real World

When they compared their new, time-sensitive model with real data from the STAR experiment (which actually observes these collisions), the results lined up perfectly.

• The model correctly predicted how many of each type of particle group was made.
• It confirmed that different groups form at different times.
• It showed that the "tighter" the group (the more strongly bound), the earlier it forms.

Looking Ahead: Predicting the Future

Finally, the paper used their new understanding to make a prediction. They calculated how many even heavier, stranger groups (containing two "strange" particles) might be formed in future experiments. They predicted that while these groups are rare, they should be detectable if scientists look at the right moment in the collision.

Summary

In short, this paper says: "Don't assume all particle groups form at the same time."

• Small, loose groups form late, when things calm down.
• Large, tight groups form early, while things are still chaotic and crowded.
• To understand the universe's building blocks, we need to look at the timing of the collision, not just the final result.

Technical Summary

Technical Summary: Sequential Clusterization of Light Nuclei and Hypernuclei in Heavy-Ion Collisions within a Wigner Function Coalescence Framework

Problem Statement
The production of light nuclei and hypernuclei in heavy-ion collisions remains a complex challenge, particularly at low to intermediate energies ($\sqrt{s_{NN}} \sim 2\text{--}10$ GeV). While the coalescence model is the standard approach for describing cluster formation at kinetic freeze-out, current implementations face significant limitations. Specifically, the model often relies on Gaussian approximations for the Wigner phase-space densities of clusters, which fail to capture the extended spatial structures of loosely bound systems like deuterons and light hypernuclei. Furthermore, standard coalescence models typically assume a universal formation time at the final kinetic freeze-out. This assumption may be inadequate for compact, strongly bound nuclei (e.g., $^4$He) which possess higher binding energies and smaller radii, potentially forming earlier in the evolution of the hadronic phase. Additionally, theoretical approaches frequently underestimate the production yields of $A=4$ clusters, suggesting that important formation channels, such as sequential cluster-nucleon coalescence, are missing.

Methodology
The authors investigate these issues using a microscopic transport framework combining the Parton-Hadron-Quantum-Molecular Dynamics (PHQMD) model with a Wigner-function-based coalescence procedure.

• Transport Model: The study utilizes the PHQMD model with a mean-field mode and a hard equation of state ($\kappa = 380$ MeV) to simulate Au+Au collisions at $\sqrt{s_{NN}} = 3$ GeV. The default clusterization algorithms in PHQMD are disabled to isolate the coalescence mechanism.
• Realistic Wave Functions: Instead of Gaussian ansatzes, the Wigner phase-space densities are constructed from realistic few-body wave functions. These are obtained by numerically solving the Schrödinger equation using the hyperspherical harmonics formalism. This provides a parameter-free description of the intrinsic structure of light nuclei and hypernuclei.
• Time-Dependent Coalescence: The coalescence time ($t_{coal}$) is treated as a variable parameter, implemented by applying the coalescence calculation at different evolution cutoff times (ranging from 20 to 145 fm/c) during the hadronic phase. This allows for the exploration of whether different clusters form at distinct stages of the collision evolution.
• Formation Channels: The study incorporates sequential formation channels for $A=4$ systems. For $^4$He, channels include direct four-body coalescence ($p+p+n+n$) as well as sequential channels ($^3$He + $n$ and $t + p$). Similarly, for $^4_\Lambda$H, the direct channel is supplemented by the $^3_\Lambda$H + $n$ channel.
• Weighting: To correct for discrepancies between the transport model and experimental yields of protons and $\Lambda$ hyperons, an empirical weighting procedure is applied to the constituent particles before coalescence.

Key Contributions

1. Parameter-Free Wigner Densities: The work replaces phenomenological Gaussian approximations with realistic, angle-independent Wigner functions derived from hyperspherical harmonic solutions, offering a more rigorous description of cluster structures.
2. Sequential Formation Mechanisms: The paper explicitly introduces and quantifies the impact of cluster-nucleon formation channels for $A=4$ systems, addressing the long-standing underestimation of $^4$He and $^4_\Lambda$H yields in theoretical models.
3. Species-Dependent Formation Times: By varying the coalescence cutoff time, the authors demonstrate that the optimal formation time is not universal but depends on the specific cluster species and its binding properties.

Results

• Coalescence Times: Comparisons with STAR experimental data for rapidity distributions ($dN/dy$) reveal a clear species-dependent pattern. Light nuclei ($d, t, ^3$He) and the hypernucleus $^3_\Lambda$H are best described by coalescence times in the range of $\sim 80\text{--}125$ fm/c, corresponding to later stages of the hadronic evolution. In contrast, $A=4$ systems ($^4$He and $^4_\Lambda$H) favor significantly earlier coalescence times ($\sim 40\text{--}80$ fm/c).
• Impact of Sequential Channels: The inclusion of sequential channels ($^3$He+$n$, $t+p$ for $^4$He; $^3_\Lambda$H+$n$ for $^4_\Lambda$H) substantially improves the description of experimental yields. For $^4$He, the direct four-body channel alone fails to reproduce the data, whereas the multi-channel approach yields a $\chi^2/NDF$ of 1.5 compared to 18.2 for the direct channel.
• Kinematic Distributions: Using the extracted "matched" coalescence times, the model successfully reproduces the transverse momentum ($p_T$) spectra and mean transverse momentum ($\langle p_T \rangle$) for light nuclei and hypernuclei, particularly in the lower $p_T/A$ region ($\lesssim 1$ GeV/c).
• Binding Energy Correlation: When plotted against binding energy per nucleon, the extracted coalescence times suggest a trend where more strongly bound nuclei form earlier. This implies that heavier clusters retain information from the denser, earlier phases of the reaction, whereas loosely bound systems form later in a dilute environment.
• Predictions: The framework provides predictions for the rapidity distributions of heavier hypernuclei, specifically $^5_\Lambda$He and the double-strangeness system $^5_{\Lambda\Lambda}$He, assuming an early coalescence time ($t_{coal} = 40$ fm/c).

Significance
The paper argues that the formation of light nuclei and hypernuclei is not a single-step process occurring solely at kinetic freeze-out but is sensitive to the specific binding properties and spatial extent of the clusters. The finding that $A=4$ systems require earlier formation times suggests they are sensitive probes of the high-baryon-density region of the collision, distinct from the dilute late-stage environment where $A=2$ and $A=3$ clusters form. By incorporating realistic wave functions and sequential formation channels, the study resolves discrepancies in $A=4$ yields and offers a more robust, parameter-free framework for interpreting cluster production data. The results highlight the necessity of moving beyond universal freeze-out assumptions to understand the dynamical evolution of multi-baryon correlations in heavy-ion collisions.