Proton and kaon production in Au+Au collisions at $\sqrt{s_{\rm NN}}=3$ GeV

Using an extended isospin- and momentum-dependent Boltzmann-Uehling-Uhlenbeck transport model, this study demonstrates that incorporating momentum dependence in the nuclear mean field is essential for accurately reproducing STAR experimental data on proton, kaon, and $\Lambda$ production in Au+Au collisions at $\sqrt{s_{\rm NN}}=3$ GeV, whereas momentum-independent models only partially describe the results.

The Gist

Imagine you are trying to understand how a giant, super-dense crowd behaves when two massive groups of people crash into each other at high speed. In the world of physics, this "crowd" is made of protons and neutrons (nucleons), and the crash is a heavy-ion collision. Scientists smash gold atoms together to recreate conditions similar to the inside of a neutron star or the universe just after the Big Bang.

This paper is like a detective story where the authors are trying to figure out the "rules of the road" that govern how these particles push and pull on each other during the crash. Specifically, they are looking at collisions happening at a specific energy level (3 GeV) and asking: What is the best way to describe the force field that holds these particles together?

Here is the breakdown of their investigation using simple analogies:

1. The Two Competing Theories

The scientists tested two different ideas about how particles interact:

• The "Static" Rulebook (Momentum-Independent): Imagine a crowd where everyone pushes back with the same force, no matter how fast they are running. Whether you are jogging or sprinting, the resistance you feel is exactly the same. This is the "Momentum-Independent" (MID) model.
• The "Dynamic" Rulebook (Momentum-Dependent): Imagine a crowd where the resistance changes based on your speed. If you run fast, the crowd pushes back harder; if you walk slowly, they push back less. This is the "Momentum-Dependent" (MDI) model.

They also tested two different "stiffness" settings for the crowd:

• Soft: The crowd is like a sponge; it squishes down easily.
• Stiff: The crowd is like a steel wall; it resists squishing.

2. The Experiment: The Gold Crash

The researchers used a super-computer simulation (a digital crash test) to smash gold atoms together. They looked at three main things the particles did after the crash:

• The Speed: How fast were the protons moving sideways?
• The Direction: Did the particles flow straight ahead or bounce off to the sides in an oval shape?
• The New Guests: Did the crash create new particles like Kaons (a type of meson) and Lambda particles (a type of hyperon)?

They compared their computer simulations against real data collected by the STAR experiment at the Relativistic Heavy Ion Collider (RHIC).

3. The Findings: Who Won the Race?

The results were very clear:

• The "Static" Rulebook Failed: Whether they used the "Soft" or "Stiff" version of the Static Rulebook, the computer simulations could only partially match the real-world data. It was like trying to describe a complex dance using only a simple, rigid robot. It got some steps right, but missed the nuance.
• The "Dynamic" Rulebook Won: The model that accounted for speed (the Momentum-Dependent model) with a "Soft" crowd setting (specifically, an incompressibility value of 230 MeV) matched the real data almost perfectly.

The Analogy: Think of the particles as cars on a highway.

• The Static model assumes the air resistance is the same whether the car is going 10 mph or 100 mph. When the scientists tried to predict the crash, the cars didn't scatter the way real cars do.
• The Dynamic model knows that air resistance increases with speed. When they used this rule, the simulated cars scattered and flowed exactly like the real data showed.

4. Why Does This Matter?

The paper concludes that to understand the "Equation of State" (the rulebook for how dense matter behaves), you must account for the fact that the force between particles changes depending on how fast they are moving.

Even though the "Static" models could explain some parts of the crash (like the sideways flow of protons), they failed to explain the oval-shaped flow (elliptic flow) and the production of new particles (Kaons and Lambdas) as accurately as the "Dynamic" model.

The Bottom Line

The authors found that the universe, at least in these high-energy crashes, behaves like a dynamic crowd where the push-and-pull forces change based on speed. Ignoring this speed-dependence leads to an incomplete picture. By using the model that respects this speed-dependence, they were able to accurately describe the behavior of protons, Kaons, and Lambda particles in gold-gold collisions, helping us better understand the fundamental properties of dense nuclear matter.

Technical Summary

1. Problem Statement

The study addresses the significant uncertainty surrounding the Equation of State (EoS) of nuclear matter at high densities (3–4 times the saturation density, $\rho_0$). While heavy-ion collisions (HICs) provide a laboratory to study these conditions, previous analyses of data from the STAR Collaboration at $\sqrt{s_{NN}} = 3$ GeV have yielded conflicting conclusions:

• Some studies suggest that a momentum-independent (MID) nuclear mean field with either a soft ($K_0 = 230$ MeV) or stiff ($K_0 = 380$ MeV) EoS can reproduce experimental data.
• Other studies indicate that only a momentum-dependent (MDI) mean field with a soft EoS provides a good fit.

The core problem is determining whether the momentum dependence of the nuclear mean field is a critical factor in describing collective flows and particle production at this specific energy, and resolving the ambiguity between soft and stiff EoS scenarios when momentum dependence is neglected.

2. Methodology

The authors utilized an extended Isospin- and Momentum-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model to simulate Au+Au collisions at $\sqrt{s_{NN}} = 3$ GeV.

• Model Upgrades: The model was upgraded to include higher-energy reaction channels and new particles/resonances essential for this energy regime, specifically:
  • Resonances: $N^*(1535)$, $\Delta(1232)$, $N^*(1440)$.
  • Mesons and Hyperons: $\eta$, $K$, $\Lambda$, $\Sigma$.
  • Reaction channels: $NN \to NN + \eta$, $NN \to NN + K$, and associated production channels.
• Cross Sections: Cross sections for resonance production were determined using parameterization formulas (VerWest and Arndt) and empirical data (e.g., for $N^*(1535)$ via $\eta$ production).
• Mean Field Scenarios: Three distinct scenarios were compared against STAR experimental data:
  1. MDI (Momentum-Dependent): Incompressibility $K_0 = 230$ MeV (Soft EoS).
  2. MID (Momentum-Independent): Soft EoS with $K_0 = 230$ MeV.
  3. MID (Momentum-Independent): Stiff EoS with $K_0 = 380$ MeV.
• Observables: The study focused on:
  • 0–10% Centrality: Proton transverse momentum ($p_t$) spectra and rapidity-dependent mean transverse momentum ($\langle p_t \rangle$).
  • 10–40% Centrality: Directed flow ($v_1$) and elliptic flow ($v_2$) for protons, $K^+$ mesons, and associated $\Lambda$ hyperons.
• Kaon Potential: The simulation incorporated the empirical kaon-nucleon scattering length scenario, including medium modifications to the kaon mass and a repulsive kaon potential.

3. Key Contributions

• Systematic Comparison: The paper provides a rigorous comparison of momentum-dependent vs. momentum-independent mean fields specifically at $\sqrt{s_{NN}} = 3$ GeV, a region where the transition from hadronic to partonic degrees of freedom is being investigated.
• Resolution of EoS Ambiguity: It clarifies that while different models can fit certain observables (like directed flow), they fail to simultaneously describe the full set of flow observables (particularly elliptic flow) without the correct momentum dependence.
• Extension to Strangeness: The study extends the analysis of flow dynamics to strange particles ($K^+$ and $\Lambda$), demonstrating that the conclusions drawn from protons hold true for strangeness production.

4. Results

The simulation results revealed distinct differences between the three scenarios:

• Central Density Evolution: The MDI scenario ($K_0=230$ MeV) produced a central compression density intermediate between the two MID scenarios. This suggests that momentum dependence and bulk incompressibility ($K_0$) have competing but compensating effects on nuclear compression.
• Proton Transverse Momentum ($\langle p_t \rangle$):
  • The MDI ($K_0=230$) and MID ($K_0=380$) scenarios reasonably fit the experimental $\langle p_t \rangle$ data.
  • The MID ($K_0=230$) scenario significantly underestimated the proton transverse momentum.
• Directed Flow ($v_1$):
  • Both MDI ($K_0=230$) and MID ($K_0=380$) successfully reproduced the proton directed flow data.
  • MID ($K_0=230$) underestimated the amplitude of the directed flow.
• Elliptic Flow ($v_2$) - The Decisive Observable:
  • MID ($K_0=380$, Stiff): Overestimated the magnitude of proton elliptic flow.
  • MID ($K_0=230$, Soft): Underestimated the magnitude of proton elliptic flow.
  • MDI ($K_0=230$): Successfully reproduced the experimental elliptic flow data. The momentum dependence enhanced the pressure gradient sufficiently to match the data without requiring a stiff EoS.
• Kaon and $\Lambda$ Flows:
  • The $K^+$ and $\Lambda$ flows exhibited similar trends to protons.
  • Only the MDI ($K_0=230$) scenario provided a fair simultaneous fit to both the directed and elliptic flows of $K^+$ and $\Lambda$. The MID scenarios failed to describe the elliptic flow data accurately.

5. Significance

• Necessity of Momentum Dependence: The study conclusively demonstrates that the momentum dependence of the nuclear mean field is essential for accurately describing nuclear matter properties at $\sqrt{s_{NN}} = 3$ GeV. Models ignoring this dependence (MID) cannot simultaneously reproduce directed and elliptic flow data, regardless of whether a soft or stiff EoS is chosen.
• EoS Constraints: The findings suggest that when momentum dependence is properly accounted for, the experimental data favors a soft nuclear EoS ($K_0 \approx 230$ MeV). Previous conclusions favoring a stiff EoS were likely artifacts of neglecting momentum dependence.
• Validation of Transport Models: The successful reproduction of STAR data using the extended IBUU model validates the inclusion of specific resonance channels ($N^*(1535)$) and strangeness production mechanisms for this energy regime.
• Implications for Neutron Stars: By constraining the EoS at 3–4 times saturation density, these results contribute to a better understanding of the internal structure of neutron stars and the behavior of dense matter in the early universe.

In summary, the paper establishes that the momentum dependence of the nuclear mean field is not a minor correction but a fundamental component required to interpret heavy-ion collision data at $\sqrt{s_{NN}} = 3$ GeV, pointing toward a soft equation of state for nuclear matter.