Freeze-out model of light nuclei formation in heavy-ion collision transport

This paper proposes a hybrid coarse-graining model that combines dynamical transport and thermal cluster production to predict light nuclei yields, spectra, and elliptic flows in semi-peripheral Au+Au collisions at 1.23 A GeV, effectively bridging nucleon and cluster descriptions at freeze-out while accounting for thermal non-uniformity and collective transport.

The Gist

Imagine a heavy-ion collision (smashing two heavy atomic nuclei together) as a chaotic, high-speed crash of two massive trucks. Inside the wreckage, the matter gets so hot and dense that it turns into a "nuclear hellfire," a soup of particles so energetic that you'd expect any small, fragile structures to be instantly vaporized.

Yet, strangely, tiny structures called light nuclei (like deuterons, which are just a proton and a neutron stuck together) survive this explosion and are found in the debris. Scientists have long wondered: How do these fragile things survive the fire?

This paper proposes a new way to understand and predict how these particles form and survive. Here is the breakdown using simple analogies:

The Problem: Two Different Ways of Looking at the Crash

Currently, scientists use two main tools to study these crashes, but they don't always agree:

1. The "Traffic Cam" (Transport Models): This tracks every single particle (protons and neutrons) as they bounce around like billiard balls. It's great for seeing how they move, but it's terrible at predicting when they decide to stick together to form a cluster. It's like trying to predict a traffic jam by watching every car individually; you miss the big picture of the gridlock.
2. The "Weather Report" (Thermal Models): This treats the matter like a gas in a room. It assumes everything has settled down and reached a comfortable temperature. It's great at predicting how many clusters form based on temperature, but it ignores the fact that the "room" is expanding and swirling with currents.

The Solution: The "Hybrid Freeze-Out" Model

The authors propose a new approach called the Hybrid Coarse-Grained Freeze-Out (HCGF) model. Think of it as a smart switch that changes the camera angle at the perfect moment.

1. The Hot Phase (The Traffic Cam): In the beginning, when the crash is hottest and most violent, the model tracks individual particles (protons and neutrons) as they zoom around.
2. The "Freeze-Out" Moment (The Switch): As the explosion expands, the density drops. The authors set a specific "freeze-out" line (a density threshold). Once the matter drops below this density, the model stops tracking individual bounces.
3. The Thermal Phase (The Weather Report): At this exact moment, the model says, "Okay, the chaos has settled enough." It instantly calculates how many clusters form based on the local temperature and pressure, just like a weather report predicts rain based on humidity.

The Key Insight:
The paper argues that when these clusters form, they release a tiny bit of energy (like a magnet snapping shut). This release actually makes the local temperature slightly higher than if the particles had stayed separate. The model accounts for this "heating up" effect, which previous methods often missed.

What Did They Find?

The team tested this model on a specific type of collision (Gold nuclei smashing into Gold nuclei). Here is what they discovered:

• It Matches Reality: The model successfully predicted how many protons, neutrons, and light clusters were produced, matching real-world data from the HADES experiment.
• Clusters are "Late Bloomers": The model shows that light clusters form later in the explosion than free protons. Because they form later, they are carried by the "wind" of the explosion (collective flow) differently.
• Temperature Differences: The model reveals that the free protons come from a wider range of temperatures (some hot, some cooler), while the clusters mostly come from a specific, slightly cooler "zone" where the conditions were just right for them to stick together.

The Big Picture

Think of the explosion as a giant, expanding balloon.

• Old models tried to guess the final contents of the balloon by either watching every rubber molecule bounce (too messy) or by assuming the balloon was a static room (too simple).
• This new model watches the molecules bounce until the balloon stretches enough, then instantly calculates the final contents based on the balloon's current size and temperature.

By combining the movement of the particles with the rules of thermal equilibrium, this new "Hybrid" model gives a much clearer picture of how the universe builds these fragile nuclear structures out of the ashes of a nuclear fire. It helps scientists better understand the "rules of the road" (the Equation of State) that govern how matter behaves under extreme pressure.

Technical Summary

Technical Summary: Freeze-out Model of Light Nuclei Formation in Heavy-Ion Collision Transport

Problem Statement
Light nuclei (clusters) are abundant products in heavy-ion collisions (HIC) at intermediate beam energies and serve as sensitive probes for inferring the nuclear Equation of State (EoS), particularly regarding pressures in the high-density region. However, predicting cluster production remains more challenging than predicting nucleon yields within standard transport models. While some models treat clusters dynamically as independent particles, most rely on post-processing transport results using methods like coalescence or simulated annealing. These approaches often struggle to simultaneously account for the dynamical transport features of the system and the thermal non-uniformities present in the hot, strongly interacting matter created during collisions. Furthermore, the persistent survival of these fragile clusters in a "nuclear hellfire" environment requires a deeper understanding of their formation and evolution.

Methodology
The authors propose a novel hybrid approach termed the Hybrid Coarse-Grained Freeze-out (HCGF) model. This method bridges dynamical transport simulations and thermal cluster production specifically for mid-rapidity particles. The core philosophy involves a "coarse-graining" step where the complexity of the system is reduced at the freeze-out stage to a level consistent with practical scales, similar to techniques used in polymer physics or hydrodynamics.

The model operates through the following steps:

1. Transport Phase: Participant nucleons propagate via a transport model (specifically UrQMD 4.0 with momentum-dependent mean-field potentials) until the system reaches a defined local freeze-out density, $n_{fr} \approx 0.02 \text{ fm}^{-3}$ (approximately 1/8 of saturation density).
2. Chemical Freeze-out Matching: At $n_{fr}$, the transport description is matched to a thermal model. The authors equate the densities of conserved quantities from the transport simulation (energy density $\varepsilon$, baryon density $n_B$, isospin density $n_{I3}$, and strangeness density $n_S$) to their local equilibrium values calculated within an Ideal Nucleon-and-Cluster Gas (INCG) model.
   • The matching conditions are defined by Eqs. (1a)–(1d), solving for temperature ($T$) and chemical potentials ($\mu_B, \mu_{I3}, \mu_S$) that satisfy the conservation laws.
3. Cluster Production: Once $T$ and chemical potentials are determined, yields and momentum distributions for light clusters (protons, neutrons, deuterons, tritons, $^3$He, $^4$He) are calculated assuming chemical equilibrium in the presence of a collective flow field.
4. Post-Processing: The resulting clusters and free nucleons propagate to detectors without further inelastic reactions. The model assumes that scattering after chemical freeze-out does not significantly alter chemical composition.

The approach is designed to be agnostic regarding the specific transport or thermal model used, provided the particle states are fixed. In this study, it is applied to semi-peripheral Au+Au collisions at $E_{lab} = 1.23$ AGeV.

Key Results
The application of the HCGF model yields several significant findings:

• Thermal Effects of Clustering: The formation of clusters releases binding energy, which raises the freeze-out temperature compared to a nucleon-only scenario. For Au+Au collisions at 1.23 AGeV, the average freeze-out temperature for the cluster-inclusive model is approximately 30 MeV, compared to 24 MeV for a nucleon-only model. This "coalescence heating" is a direct consequence of the energy release during cluster formation.
• Yield Predictions: The model successfully reproduces experimental data from the HADES collaboration for proton, neutron, and light cluster yields. The HCGF predictions align well with the grey band representing experimental data, outperforming the raw UrQMD transport simulation without explicit clustering.
• Spectral Shapes: The model predicts that light cluster spectra peak at lower scaled transverse momentum ($p_t/A$) compared to protons. This is consistent with the expectation that clusters are emitted later in the evolution when transverse flow and temperature are lower. The asymptotic slopes of the spectra reflect the effective temperature and collective flow of the emission region.
• Elliptic Flow: The analysis of elliptic flow ($v_2$) reveals that protons are more sensitive to the full temperature distribution of the system than clusters. Protons originate from a wider range of freeze-out temperatures, whereas clusters predominantly emerge from cooler regions ($T \lesssim 40$ MeV). Differences between proton and cluster flow become pronounced at higher transverse momenta ($p_t/A \gtrsim 0.4$ GeV/c).

Significance and Claims
The paper claims that the HCGF framework provides a thermodynamically consistent method for describing cluster production at mid-rapidity. Its primary significance lies in its ability to:

1. Preserve the dynamical features of transport models (such as collective flow and non-uniformity) while incorporating the internal energy of bound states.
2. Offer a simple, instantaneous description of cluster production that avoids the technical difficulties and ambiguities associated with dynamic cluster-production modeling.
3. Successfully describe yields, spectral shapes, and flow observables for both nucleons and light clusters within a unified ensemble.

The authors emphasize that this approach is particularly effective for mid-rapidity particles where local thermal equilibrium is expected to be reached at late stages of evolution. While the current implementation focuses on ground states, the framework is presented as extensible to include excited states, decays, and particles at large rapidity in future work. The model serves as a bridge between hydrodynamic descriptions and transport simulations, offering a robust tool for constraining the nuclear EoS using cluster observables.