Three-dimensional sizes and shapes of pion emission in heavy-ion collisions

This paper presents a three-dimensional analysis of the two-pion source in Au+Au collisions at 200 GeV using Monte-Carlo simulations and compares these results with recent centrality-dependent measurements from the PHENIX Collaboration to refine the description of pion emission shapes using Levy-stable distributions.

The Gist

Imagine you are trying to figure out the shape of a cloud of smoke that was just released from a firework, but you can't see the smoke itself. Instead, you only have a high-speed camera that takes pictures of the smoke particles as they fly away from each other. By looking at how close or far apart these particles are when they hit your camera, you can try to reconstruct what the original cloud looked like.

This is essentially what physicists do in heavy-ion collisions. They smash gold atoms together at nearly the speed of light to create a tiny, super-hot "soup" of particles (called a quark-gluon plasma). As this soup cools down, it freezes into a cloud of pions (a type of subatomic particle). The scientists want to know: What is the 3D shape and size of this pion cloud?

Here is a breakdown of the paper's story, using simple analogies:

1. The Detective Work: "Femtoscopy"

The scientists use a technique called femtoscopy. Think of it like trying to guess the size of a room by listening to how echoes bounce off the walls.

• The Echo: When two identical pions are created close together, they behave like twins. Because of a quantum rule (Bose-Einstein statistics), they tend to stick together or avoid each other in specific ways depending on how far apart they are.
• The Clue: By measuring the "relative momentum" (how fast they are moving away from each other), the scientists can work backward to figure out the distance between them when they were born.

2. The Shape of the Cloud: The "Levy" Distribution

For a long time, scientists thought these clouds were simple, round, and smooth, like a perfectly fluffy marshmallow (a Gaussian distribution).

• The Twist: New, more precise data showed that these clouds aren't perfect marshmallows. They have "tails." Imagine a marshmallow that has some long, stringy bits stretching out far from the center.
• The Solution: The paper uses a mathematical shape called a Lévy-stable distribution. Think of this as a "super-marshmallow" that can stretch out into long, thin tendrils. This shape fits the data much better than a simple round ball.

3. The Simulation vs. Reality: The "Video Game" Test

The authors used a super-computer simulation called EPOS3. You can think of EPOS3 as a incredibly detailed video game that simulates the laws of physics for these collisions.

• The Goal: They wanted to see if their video game could reproduce the real-life "photos" taken by the PHENIX experiment (a real detector at the Relativistic Heavy Ion Collider).
• The Setup: They simulated 300,000 collisions of gold atoms and measured the "pion clouds" in the game.

4. The Results: A Tale of Two Collisions

When they compared the "Video Game" (EPOS3) to the "Real Photos" (PHENIX data), they found a split result:

• The Peripheral Collisions (The Glancing Blows):
  Imagine two cars barely brushing past each other. The resulting "pion cloud" is small and messy.

  • Result: The video game matched the real photos perfectly. The shape, size, and "stretchiness" (the Lévy parameters) were spot on.

• The Central Collisions (The Head-On Crash):
  Imagine two cars smashing head-on. This creates a massive, hot, dense explosion.

  • Result: The video game started to drift apart from reality.
    • The Shape: The game predicted the cloud was too "round" or didn't have the right "stretchiness" (the Lévy index $\alpha$ was off).
    • The Size: The game got the general size right, but there were small errors, especially in the direction of the beam.
    • The "Halo" Effect: The game struggled to explain the "correlation strength" (how many pions are "twins" vs. how many are just random noise) in these big crashes.

5. Why the Discrepancy? (The Missing Ingredients)

Why did the video game fail for the big crashes? The authors suggest a few missing ingredients in their recipe:

1. The "Static Electricity" Effect (Coulomb Scattering): Pions are electrically charged. In a massive, dense crash, they push and pull on each other like magnets. The simulation might not be accounting for this "static electricity" push perfectly.
2. The "Heavy" Particles: The simulation might be missing some specific heavy particles that decay into pions later, creating those long "tails" in the cloud.
3. The "Medium" Effect: Inside the hot soup, particles might act heavier or different than they do in a vacuum. The simulation assumes they act normally, but maybe they don't.

The Big Takeaway

The paper concludes that while our current "video games" (EPOS3) are amazing at predicting what happens in glancing collisions, they need a software update to handle the massive, head-on crashes perfectly.

However, there is good news: The simulation did correctly predict the "scaled" behavior (how the shape changes as you look at different energy levels). This means the basic physics engine is working; we just need to add a few more "patches" (like better handling of electric forces or particle decays) to make the simulation perfect for the most extreme collisions.

In short: We are getting closer to understanding the 3D shape of the universe's smallest explosions, but we still have a little bit of "glitch" to fix in our simulations for the biggest crashes.

Technical Summary

1. Problem Statement

In high-energy heavy-ion physics, femtoscopy (specifically two-pion Bose-Einstein correlation measurements) is used to reconstruct the space-time geometry of the particle-emitting source created in collisions. While experimental data from the PHENIX collaboration and others have confirmed that the pion pair source is well-described by Lévy-stable distributions (rather than simple Gaussians) due to long tails caused by resonance decays and hadronic rescattering, a direct quantitative comparison between state-of-the-art phenomenological models and experimental data remains necessary.

The specific problem addressed is the lack of a comprehensive, three-dimensional comparison between the EPOS3 Monte Carlo simulation model and recent PHENIX centrality-dependent measurements. Previous studies (e.g., Ref. [11]) established the viability of Lévy distributions in EPOS3 on an event-by-event basis, but a detailed analysis of centrality dependence and a direct comparison with experimental source parameters (Lévy index $\alpha$, scale $R$, and correlation strength $\lambda$) were needed to validate the model's ability to describe the underlying physics of pion emission.

2. Methodology

The authors performed a rigorous analysis using the following steps:

• Simulation Setup:

  • Model: EPOS3 (version 359), a hybrid model combining hydrodynamic evolution with a hadronic rescattering phase.
  • System: Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV.
  • Sample: 300,000 minimum bias events.
  • Kinematics: Pions selected with $0.15 < p_T < 1.0$ GeV/c and $|\eta| < 1$. Pairs were selected based on relative momentum $|q_{LCMS}| < \sqrt{0.15 \cdot m_T}$ to match experimental fit limits.
  • Centrality: Divided into six 10% wide classes (0–60%).
  • Transverse Mass ($m_T$): Analyzed in 11 bins ranging from 0.175 to 0.725 GeV/c.

• Source Reconstruction:

  • The three-dimensional spatial separation vector $\rho$ was calculated in the Longitudinal Co-Moving System (LCMS) using "out-side-long" coordinates.
  • The pair source distribution $D(\rho)$ was constructed from freeze-out coordinates.
  • Core-Halo Picture: The analysis distinguishes between the "core" (primordial pions and short-lived resonances like $\rho, \Delta$) and the "halo" (long-lived resonances like $\eta, K_S^0, \Lambda$). The correlation strength $\lambda$ is defined as the integral of the core-core contribution $D^{(c,c)}$.

• Fitting Procedure:

  • One-dimensional projections of $D(\rho)$ along the out, side, and long axes were fitted simultaneously using 1D Lévy-stable distributions.
  • Fit Function: $D^{(c,c)}(\rho_\nu) = \lambda \cdot L_{1D}(\rho_\nu, \alpha, R_\nu)$, where $\alpha$ is the Lévy exponent and $R$ is the scale parameter.
  • Optimization: Log-likelihood minimization was used (preferred over $\chi^2$ due to low statistics in tail bins).
  • Event Averaging: To ensure stability, the number of events summed ($N_{events}$) was optimized per centrality and $k_T$ bin (ranging from 50 to 3000) to reach convergence of fit parameters.

• Systematic Uncertainties:

  • Variations were tested for $N_{events}$, the pair kinematic limit ($Q_{max}$), the fit range ($\rho_{fit}^{max}$), and the normalization limit ($\rho_{\lambda}^{max}$).
  • The normalization limit $\rho_{\lambda}^{max}$ was varied from 1 mm to 50 cm to simulate different experimental tracking resolutions (DCA cuts).

3. Key Contributions

• 3D Lévy Analysis in EPOS3: Provided the first detailed 3D analysis of the pion source in EPOS3 simulations, confirming that the source shape is well-described by an elliptically contoured Lévy-stable distribution.
• Direct Data-Model Comparison: Conducted a systematic, centrality-dependent comparison of EPOS3 results against the latest PHENIX experimental data (Ref. [7]).
• Resolution of Scaling Ambiguities: Demonstrated that while absolute correlation strengths ($\lambda$) depend heavily on experimental spatial resolution (tracking cuts), the scaled correlation strength ($\lambda/\lambda_{max}$) is robust and well-described by the model across all centralities without requiring in-medium modifications.
• Identification of Discrepancies: Pinpointed specific regimes (central collisions) where the model fails to reproduce the Lévy index ($\alpha$) and source sizes, suggesting missing physics in the current EPOS3 implementation.

4. Results

• Lévy Index ($\alpha$):

  • EPOS3 captures the $m_T$ trend but fails to reproduce the centrality dependence observed in PHENIX data.
  • Discrepancy: There is an increasing deviation in the absolute magnitude of $\alpha$ for central collisions (0–20%), where the simulation predicts values that differ significantly from data.
  • Possible Causes: The authors suggest missing Coulomb scattering effects (which affect quantum-statistical correlations) or in-medium mass/width modifications of hadrons.

• Lévy Scale Parameters ($R_{out}, R_{side}, R_{long}$):

  • The simulations reproduce the qualitative trends and absolute magnitudes of the radii.
  • Discrepancy: The most central classes show deviations, particularly in the longitudinal direction. The average scale $R_{avg}$ is well described for centralities >20%, but discrepancies persist in the most central collisions.

• Correlation Strength ($\lambda$):

  • The absolute value of $\lambda$ is highly sensitive to the normalization range $\rho_{\lambda}^{max}$ (simulating detector resolution).
  • Different $\rho_{\lambda}^{max}$ values are required to match data for different centralities (e.g., 1 mm for 0–10%, 50 cm for peripheral), suggesting that experimental pair resolutions or hadronic dynamics (neglected resonances) may vary or be incomplete in the model.

• Scaled Correlation Strength ($\lambda/\lambda_{max}$):

  • Major Success: The EPOS3 model reproduces the $\lambda/\lambda_{max}$ data across all centralities and $m_T$ bins with high statistical compatibility (Confidence Level > 0.1%).
  • This indicates that the core-halo fraction and its kinematic variation are correctly modeled, removing the need for ad-hoc in-medium modifications to explain this specific observable.

5. Significance

• Validation of Transport Models: The study confirms that hydrodynamics-based transport models coupled with hadronic rescattering (like EPOS3) can capture key femtoscopic observables, specifically the Lévy-stable nature of the source and the scaled correlation strength.
• Physics Insights: The results challenge the interpretation in previous experimental analyses (Ref. [7]) that suggested in-medium modifications were necessary to explain $\lambda/\lambda_{max}$. The authors show that standard hadronic dynamics in EPOS3 suffice for the scaled strength.
• Future Directions: The persistent discrepancies in central collisions for the Lévy index ($\alpha$) and absolute sizes highlight specific areas for model improvement. The authors propose that incorporating Coulomb scattering and in-medium modifications of hadronic properties are the next critical steps to bridge the gap between simulation and precision experimental data.

In conclusion, the paper establishes EPOS3 as a robust tool for describing pion emission geometries in peripheral and mid-central collisions, while identifying central collisions as a regime where additional physical effects must be incorporated to achieve precision agreement with data.