Elliptic flow of deuterons from simulations with hybrid model

This study utilizes a hybrid hydrodynamic-transport model to simulate Pb+Pb collisions at 2.76 TeV, demonstrating that the coalescence mechanism for deuteron production provides a better description of experimental elliptic flow data than direct thermal production.

The Gist

Imagine a giant, high-speed collision between two heavy atomic nuclei (like lead atoms) smashing together at nearly the speed of light. This creates a tiny, super-hot fireball of energy that behaves like a perfect fluid. As this fireball cools down, it freezes into a soup of individual particles called hadrons (mostly protons and neutrons).

Sometimes, these protons and neutrons stick together to form deuterons (a simple nucleus made of one proton and one neutron). The big mystery in physics is: How do they stick together?

There are two main theories about how this happens:

1. The "Statistical Wedding" (Thermal Production): Imagine the fireball is a crowded dance floor. As the music stops (the system freezes), everyone pairs up randomly based on how crowded the room is. If the room is big enough, deuterons just "appear" out of the crowd.
2. The "High-Five" (Coalescence): Imagine the dancers are running off the floor. If two specific dancers (a proton and a neutron) happen to be running very close to each other and moving in the same direction, they high-five and stick together after they leave the main party.

The Experiment: The "Shape" of the Flow

The scientists in this paper wanted to figure out which theory is correct. They didn't just count how many deuterons were made; they looked at Elliptic Flow.

The Analogy:
Imagine the fireball isn't a perfect circle; it's shaped like a rugby ball (oval) because the two lead nuclei didn't hit head-on. As this oval fireball expands, particles flying out in the "long" direction of the oval have an easier path than those flying out the "short" way. This creates a flow pattern that is stronger in one direction than the other.

The question is: Does the way deuterons form change this flow pattern?

• If they form via Coalescence, they are like two people who just high-fived. Their flow depends on how close they were to each other at the moment they left.
• If they form via Thermal Production, they are like a couple that was already married before leaving the party. Their flow depends on the size of the whole room they were in.

The Simulation: A Digital Time Machine

The authors built a super-complex computer simulation (a "hybrid model") to act as a time machine.

1. The Start: They used a model called TRENTo to simulate the initial smash-up.
2. The Fluid: They used vHLLE to simulate the hot, fluid-like expansion of the fireball.
3. The Aftermath: They used SMASH to simulate the chaotic traffic of particles as they cool down and scatter off each other.

They ran this simulation twice:

• Run A: Deuterons were made by the "Statistical Wedding" method (appearing at the start of the cooling phase).
• Run B: Deuterons were made by the "High-Five" method (sticking together later based on proximity).

The Results: Who Won?

When they compared their digital results to real data from the ALICE experiment at CERN (the Large Hadron Collider):

• The "Statistical Wedding" (Thermal) failed.
  The simulation predicted that deuterons formed this way would have a stronger directional flow than what was actually observed. It was like predicting that a couple married in a huge ballroom would run faster than they actually did. Also, this method slightly underestimated the total number of deuterons in less central collisions.

• The "High-Five" (Coalescence) won.
  The simulation where protons and neutrons stuck together based on their proximity matched the real-world data almost perfectly. The flow patterns and the total numbers lined up with what the ALICE experiment saw.

The Twist: Why is this different from before?

The authors noted that in a previous study (by two of the same authors), the results were the opposite. In that older study, the "Statistical Wedding" seemed to work better.

Why the change?
The difference is in the timing.

• In the old study, the "Statistical Wedding" happened instantly at the very end. The deuterons were born and left immediately.
• In this new, more realistic study, the "Statistical Wedding" deuterons were born at the start of the cooling phase, but then they had to survive a chaotic traffic jam (scattering with other particles) for a long time before escaping. This extra traffic jam messed up their flow pattern, making them look different from reality.

The Bottom Line

Using a more realistic, modern simulation, the authors concluded that deuterons are likely formed by the "High-Five" method (Coalescence). They don't just magically appear from the heat; they are formed when protons and neutrons happen to be running close together as the fireball cools down.

This study shows that looking at the shape of how particles flow (elliptic flow) is a powerful tool to distinguish between different theories of how matter is built, even though it requires very sophisticated computer models to get the answer right.

Technical Summary

Here is a detailed technical summary of the paper "Elliptic flow of deuterons from simulations with hybrid model" by Tomáš Poledníček, Radka Vozábová, and Boris Tomášik.

1. Problem Statement

The production mechanism of light nuclear clusters (specifically deuterons) in relativistic heavy-ion collisions remains an open question in high-energy nuclear physics. Two primary theoretical frameworks compete:

• Statistical (Thermal) Model: Assumes deuterons are formed in chemical equilibrium at the hadronization hypersurface, with abundances determined by the grand-canonical ensemble and a freeze-out temperature (~150 MeV).
• Coalescence Model: Assumes deuterons form dynamically after nucleons decouple from the system, binding together based on their proximity in phase space (momentum and position).

While the statistical model successfully describes cluster multiplicities, it faces conceptual challenges regarding the binding energy of deuterons and the screening effects in a dense hadronic medium. Previous studies suggested that differential elliptic flow ($v_2$) could act as a discriminator between these mechanisms, as coalescence depends on the size of the homogeneity region relative to the cluster wave function, whereas thermal production scales with volume. However, a previous study by the authors (Ref. [30]) using a blast-wave model found that coalescence yielded higher $v_2$ than thermal production. This paper aims to re-evaluate this hypothesis using a more realistic hybrid model to determine if $v_2$ can definitively distinguish the production mechanisms.

2. Methodology

The authors employed a state-of-the-art hybrid simulation chain to model Pb+Pb collisions at $\sqrt{s_{NN}} = 2.76$ TeV. The simulation framework consists of:

• Initial Conditions: Generated using TRENTo 3D, a 3D model for initial energy and entropy density distributions.
• Deconfined Phase: Evolved using vHLLE (a viscous hydrodynamics solver) with shear viscosity $\eta/s = 0.12$ and bulk viscosity $\zeta/s$ peaking at 0.4 near the critical temperature.
• Particlization: The fluid evolves until a critical energy density ($\epsilon_{crit} = 0.32$ GeV/fm$^3$), where a Hadron Sampler converts fluid elements into hadrons.
• Hadronic Afterburner: The hadronic phase is simulated using SMASH v2.2 (Simulating Many Accelerated Strongly-interacting Hadrons). This stage handles rescattering, resonance decays, and annihilation.

Production Mechanisms Tested:

1. Coalescence: Nucleons (protons and neutrons) are propagated through the SMASH phase until their last scattering or resonance decay. A coalescence algorithm then pairs nucleons based on relative momentum ($\Delta p < 0.310$ GeV/c) and distance ($\Delta r < 3.5$ fm) in the pair rest frame to form deuterons.
2. Direct Thermal Production: Deuterons are generated statistically alongside other hadrons at the particlization hypersurface. They are then treated as point-like particles entering the SMASH transport phase, where they undergo rescattering and potential destruction/creation.

Analysis:

• Simulations covered centrality bins from 0% to 70%.
• Approximately $1.5 \times 10^6$events were generated per centrality bin (3000 hydro events$\times$ 500 SMASH runs).
• Differential elliptic flow ($v_2$) was calculated using the scalar product method.
• Results were compared against experimental data from the ALICE collaboration.

3. Key Contributions

• Implementation of Thermal Deuterons in Hybrid Models: The authors implemented a rigorous test of the "thermal" production scenario within a full hybrid framework (Hydro + Transport), where deuterons are generated at hadronization and subsequently evolve through the hadronic phase. This contrasts with simpler blast-wave models where particles are frozen out immediately.
• Re-evaluation of $v_2$ as a Discriminator: The study challenges the universality of the previous finding that coalescence always yields higher $v_2$. It demonstrates that the choice of model (blast-wave vs. hybrid transport) fundamentally alters the relationship between production mechanism and elliptic flow.
• Quantitative Comparison: The paper provides a direct, quantitative comparison of both mechanisms against ALICE data for both transverse momentum ($p_T$) spectra and differential $v_2$.

4. Results

• Baseline Validation: The hybrid model successfully reproduced the $p_T$ spectra and elliptic flow ($v_2$) of protons, pions, and kaons across various centrality bins, validating the underlying hydrodynamic and transport parameters.
• Transverse Momentum Spectra: Both the coalescence and direct thermal production mechanisms agreed reasonably well with experimental deuteron $p_T$ spectra within error bars. However, the thermal model slightly underpredicted the yield in the most non-central bin (40–60%).
• Elliptic Flow ($v_2$) - The Critical Finding:
  • Thermal Production: In the hybrid model, deuterons produced thermally and evolved through the hadronic phase exhibited higher elliptic flow ($v_2$) than those formed via coalescence. This is attributed to the fact that thermally produced deuterons exist early in the hadronic phase, allowing them to accumulate significant $v_2$ through rescattering with the abundant pions (large cross-sections) over a long duration (up to 100 fm/c).
  • Coalescence: Deuterons formed via coalescence (at the end of the evolution) showed lower $v_2$.
  • Contrast with Previous Work: This result is the opposite of the authors' previous study (Ref. [30]) using a blast-wave model, where coalescence yielded higher $v_2$. The difference arises because the blast-wave model froze out clusters immediately, whereas the hybrid model allows thermal clusters to interact and build up flow in the hadronic phase.
• Comparison to Data:
  • Coalescence: The $v_2$ predictions from the coalescence model matched the experimental ALICE data well, particularly in central collisions.
  • Thermal Production: The direct thermal production model overshot the measured $v_2$ values significantly.

5. Significance and Conclusion

• Preference for Coalescence: Despite the reversal in the theoretical trend of $v_2$ compared to the previous blast-wave study, the coalescence model remains consistent with experimental data, while the direct thermal production model fails to reproduce the observed elliptic flow.
• Limitations of $v_2$ as a Simple Probe: The study concludes that differential elliptic flow cannot be used as a simple, universal discriminator between production mechanisms without accounting for the specific dynamics of the hadronic phase. The "thermal" $v_2$ is highly sensitive to the duration and cross-sections of the hadronic rescattering phase.
• Future Directions: The authors suggest that the discrepancy in thermal models highlights the need for a deeper understanding of deuteron transport in the hot hadronic medium. They propose that future studies should focus on larger clusters (e.g., tritons, $^4$He), though this is computationally demanding due to the suppression of yields and the complexity of modeling unbound cluster dynamics.

In summary, the paper demonstrates that while the specific behavior of elliptic flow depends heavily on the simulation framework, coalescence remains the preferred mechanism for deuteron production in Pb+Pb collisions at LHC energies when compared against high-precision experimental data.