Beam energy dependence of identified particle production in heavy-ion collisions using a parton-hadron string dynamics model

Using the parton-hadron string dynamics (PHSD) transport model, this study predicts transverse momentum spectra and yield characteristics of identified particles in Au+Au collisions across various beam energies and centralities, highlighting the critical roles of baryon stopping and strangeness production while providing theoretical context for current and future heavy-ion beam energy scan programs.

The Gist

Imagine you are a chef trying to understand the secret recipe of the universe's most extreme soup. This soup isn't made of vegetables; it's made of the fundamental building blocks of matter (quarks and gluons) smashed together at incredible speeds.

This paper is a "cooking simulation" where the authors use a sophisticated computer program (called PHSD) to predict what happens when they smash two heavy gold atoms together at different speeds. They are trying to figure out how the "ingredients" (particles like pions, kaons, protons) behave when the heat and pressure change.

Here is the breakdown of their findings in everyday language:

1. The Experiment: Smashing Gold at Different Speeds

Think of the collision like two freight trains crashing into each other.

• The Setup: They smashed gold trains together at four different speeds (energies): 6.7, 8, 11, and 25 units of speed.
• The Goal: They wanted to see how the "debris" (the particles flying out) changed depending on how hard the crash was.
• The Context: This speed range is special. It's like the "Goldilocks zone" for high-density matter. It's fast enough to create a hot, dense soup, but slow enough that the matter is packed so tightly it resembles the core of a neutron star. This is exactly the territory that future experiments in Germany (FAIR) and Russia (NICA) are planning to explore.

2. The Main Characters: The Particles

When the gold trains crash, they spit out different types of particles. The authors tracked four main groups:

• Pions ($\pi$): The light, fast "messengers" of the crash.
• Kaons ($K$): The "strange" particles, carrying a special property called "strangeness."
• Protons ($p$): The heavy, sturdy "bricks" of normal matter.
• Antiprotons ($\bar{p}$): The "evil twins" of protons. If they touch a normal proton, they annihilate each other in a flash of energy.

3. The Big Discoveries (The "Flavor" of the Crash)

A. The "Baryon Traffic Jam" (Baryon Stopping)

Imagine a highway where cars (protons) are driving toward a tunnel.

• At high speeds (25 units): The cars zoom right through the tunnel. They don't stop much.
• At low speeds (6.7 units): The cars hit a massive traffic jam right in the middle of the tunnel. They get stuck.
• The Result: At lower speeds, the crash creates a huge pile-up of protons right in the center. This is called baryon stopping. The computer model showed that as the speed went down, the number of protons in the middle actually increased because they were getting stuck there.

B. The "Vanishing Act" (Antimatter Suppression)

Now, imagine the "evil twins" (antiprotons) trying to enter that traffic jam.

• Because the traffic jam is so full of normal protons, the antiprotons get swarmed. They crash into the protons and annihilate (disappear in a puff of energy).
• The Result: At lower speeds, where the traffic jam is worst, the number of antiprotons drops dramatically. The model predicts that the "low-speed" crashes are very hostile to antimatter because there are too many normal protons around to eat them.

C. The "Strange" Shift (Kaons)

Strangeness (kaons) is produced in two ways:

1. Pair Production: Creating a strange particle and an anti-strange particle together (like buying a pair of shoes). This happens more at high speeds.
2. Associated Production: Creating a strange particle alongside a normal particle (like buying a shoe and a sock). This dominates at low speeds.

• The Result: As the speed dropped, the "pair" production stopped, and the "sock" production took over. This changed the ratio of positive to negative kaons, giving the scientists a clue about how the "soup" was being cooked.

D. The "Speed of the Debris" (Transverse Momentum)

When the trains crash, the debris flies out sideways.

• High Speed Crash: The debris flies out fast and hard (high momentum). It's like a high-pressure hose.
• Low Speed Crash: The debris moves slower and more lazily (low momentum).
• The Twist: Interestingly, the antiprotons at low speeds were flying out faster than the protons. Why? Because the slow, lazy antiprotons got eaten by the proton traffic jam. Only the fast, energetic antiprotons survived the crash. This left a "harder" (faster) average speed for the survivors.

4. Why Does This Matter?

This paper is like a practice run for future experiments.

• Scientists at the RHIC (USA), FAIR (Germany), and NICA (Russia) are about to run these exact experiments.
• The authors used their computer model to say, "If you smash gold at these speeds, here is exactly what you should expect to see."
• By comparing their predictions with the real data, scientists can figure out the "Equation of State" of the universe. Think of this as figuring out the rules of how matter behaves when it is squeezed as hard as it can possibly be squeezed.

Summary

In simple terms, this paper says: "We used a super-computer to simulate smashing gold atoms at different speeds. We found that at lower speeds, the matter gets so crowded that protons get stuck in the middle, and antimatter gets eaten up. This helps us understand how the universe behaves under extreme pressure, which is crucial for the big experiments happening right now."

Technical Summary

1. Problem Statement and Motivation

The primary objective of heavy-ion collision experiments is to map the Quantum Chromodynamics (QCD) phase diagram, specifically exploring the transition from hadronic matter to a deconfined Quark-Gluon Plasma (QGP) at high temperatures and high baryon chemical potentials ($\mu_B$).

• The Gap: While high-energy collisions (e.g., at LHC and top RHIC energies) are well-studied, the region of high baryon density and intermediate beam energies ($\sqrt{s_{NN}} \approx 4 - 7$ GeV, corresponding to $E_{lab} \approx 6.7 - 25$ A GeV) remains less understood. This regime is critical for understanding the onset of deconfinement, the nuclear equation of state (EoS), and the behavior of matter at densities comparable to neutron star interiors.
• The Challenge: Existing theoretical models often struggle to consistently describe bulk observables (yields, spectra, ratios) across this wide energy range, particularly regarding the interplay between baryon stopping, strangeness production mechanisms, and final-state interactions.
• Context: This energy range overlaps with the upcoming experimental programs at FAIR (CBM experiment) and NICA (MPD experiment), making theoretical benchmarks essential for interpreting future data.

2. Methodology

The authors utilized the Parton-Hadron String Dynamics (PHSD) transport model (version 4.1) to simulate Au + Au collisions.

• Model Framework: PHSD is a microscopic covariant dynamical approach based on the Dynamical Quasi-Particle Model (DQPM) within the Kadanoff–Baym framework. It unifies the description of partonic and hadronic phases, incorporating off-shell transport equations and scalar mean fields.
• Simulation Setup:
  • Collisions: Au + Au.
  • Beam Energies ($E_{lab}$): 6.7, 8, 11, and 25 A GeV (corresponding to $\sqrt{s_{NN}} \approx 4 - 7$ GeV).
  • Centrality: Nine centrality classes ranging from 0–5% (central) to 70–80% (peripheral), defined by charged particle multiplicity in $|\eta| < 0.5$.
  • Event Sample: 50 million minimum-bias events generated per energy.
• Observables: The study focused on identified hadrons ($\pi^\pm, K^\pm, p, \bar{p}$) at mid-rapidity ($|y| < 0.5$). Key metrics included:
  • Transverse momentum ($p_T$) spectra.
  • Integrated yields ($dN/dy$).
  • Mean transverse momentum ($\langle p_T \rangle$).
  • Particle ratios (e.g., $\bar{p}/p$, $K^-/K^+$, $K/\pi$).

3. Key Contributions

• Systematic Energy Scan: Provided a comprehensive theoretical prediction for identified particle production across the specific energy window relevant to the FAIR and NICA experiments.
• Mechanism Identification: Explicitly quantified the roles of baryon stopping, strangeness production (associated vs. pair production), and baryon-antibaryon annihilation in shaping the final state observables at high baryon density.
• Model Validation: Demonstrated that the PHSD model, which includes off-shell dynamics and in-medium modifications, can qualitatively reproduce experimental trends from AGS and RHIC, validating its use for future predictions.

4. Key Results

A. Transverse Momentum ($p_T$) Spectra

• Spectra decrease with increasing $p_T$ and from central to peripheral collisions.
• Energy Dependence: Invariant yields for pions, kaons, and antiprotons decrease as beam energy decreases. Conversely, proton yields increase as beam energy decreases, indicating the accumulation of initial-state baryons at mid-rapidity due to enhanced baryon stopping.

B. Particle Yields ($dN/dy$)

• Centrality: Yields generally scale with the number of participating nucleons ($N_{part}$), though $K^-$ and $\bar{p}$ show weaker centrality dependence at lower energies.
• Energy Trends:
  • Pions: Less sensitive to beam energy variations.
  • Kaons: $K^-$ yields are significantly suppressed compared to $K^+$ at lower energies, reflecting the dominance of associated production over pair production.
  • Protons/Antiprotons: Proton yields rise at lower energies (baryon stopping), while $\bar{p}$ yields drop sharply due to annihilation in the dense baryonic medium.

C. Mean Transverse Momentum ($\langle p_T \rangle$)

• $\langle p_T \rangle$ increases with particle mass (mass ordering).
• Centrality: $\langle p_T \rangle$ decreases from central to peripheral collisions, reflecting reduced collective transverse expansion.
• Energy: $\langle p_T \rangle$ decreases with decreasing beam energy for all species, indicating weaker collective dynamics at lower energies.
• Proton vs. Antiproton Splitting: A distinct feature is that $\langle p_T \rangle_{\bar{p}} > \langle p_T \rangle_{p}$ across all energies. This is attributed to baryon-antibaryon annihilation, which selectively depletes low-$p_T$ antiprotons, hardening their spectrum, while protons are replenished by transported initial-state baryons in the low-$p_T$ region.

D. Particle Ratios

• $\pi^-/\pi^+$: Shows a slight enhancement at the lowest energy (6.7 A GeV), attributed to isospin effects and resonance decays.
• $K^-/K^+$: Increases with beam energy. The sharp drop at lower energies signals a transition from pair production dominance (high energy) to associated production dominance (low energy).
• $\bar{p}/p$: Increases with beam energy, reflecting reduced baryon stopping and lower net-baryon density at higher energies.
• $K/\pi$ and $\bar{p}/\pi$: Increase with beam energy, signaling enhanced strangeness and pair production.
• $p/\pi^+$: Decreases with beam energy, opposite to $\bar{p}/\pi^-$, highlighting the shift from baryon-stopping dominance to pair-production dominance.

5. Significance and Conclusion

• Theoretical Benchmark: The PHSD results show good qualitative agreement with existing experimental data from AGS and RHIC, establishing the model as a reliable tool for the high-baryon-density regime.
• Physics Insights: The study confirms that at low beam energies ($E_{lab} < 25$ A GeV), the system is dominated by baryon stopping and annihilation processes. The suppression of $\bar{p}$ and $K^-$ relative to their antiparticles/counterparts provides clear signatures of the dense baryonic environment.
• Future Impact: These predictions serve as critical baselines for the CBM (FAIR) and MPD (NICA) experiments. Understanding these mechanisms is essential for:
  • Constraining the Equation of State (EoS) of dense nuclear matter.
  • Identifying signatures of the QCD phase transition and chiral symmetry restoration.
  • Distinguishing between hadronic and partonic degrees of freedom in the early stages of collisions.

In summary, the paper successfully utilizes the PHSD model to decode the complex interplay of transport and production mechanisms in heavy-ion collisions at intermediate energies, offering vital theoretical guidance for the next generation of heavy-ion physics experiments.