Centrality Dependence of the Balance Functions for Identified Particles in Pb--Pb Collisions Using Pythia + Angantyr

This study utilizes the PYTHIA 8.3 + Angantyr model to investigate centrality-dependent balance functions for identified pions, kaons, and protons in Pb--Pb collisions at 2.76 TeV, revealing that while the model qualitatively captures peripheral collision trends and resonance-induced correlation dips, it fails to quantitatively reproduce central collision data, thereby highlighting the need for a dedicated heavy-ion tuning of the Angantyr framework.

The Gist

The Big Picture: The "Couples" of the Particle World

Imagine a massive, chaotic dance party inside a tiny, invisible ballroom. This ballroom is created when two giant lead atoms (like heavy bowling balls) smash into each other at nearly the speed of light. This is what happens in Pb-Pb collisions at the Large Hadron Collider.

When these atoms crash, they create a super-hot soup of energy that instantly turns into thousands of new particles. Among these particles are pions (the lightest, most common dancers), kaons (slightly heavier, with a "strange" flavor), and protons (the heavyweights).

The scientists in this paper are trying to answer a specific question: How long does it take for these particles to be born?

To figure this out, they use a tool called a Balance Function. Think of this like a "couples detector." In the universe, particles are often born in pairs: a positive charge and a negative charge (like a positive pion and a negative pion). These pairs are "balanced" because they were created together.

The Core Concept: The "Spread" of the Dance

The Balance Function measures how far apart these "couples" are when they finally stop dancing and fly out of the ballroom.

• The "Early Bird" Theory: If a couple is born at the very beginning of the party (early in the collision), they have a lot of time to get separated. The crowd pushes them apart, and they drift far away from each other before the party ends. This results in a wide balance function (a "broad" spread).
• The "Late Night" Theory: If a couple is born just as the party is ending (late in the collision), they haven't had time to drift apart. They are still standing right next to each other when the lights go out. This results in a narrow balance function (a "tight" cluster).

The Analogy:
Imagine two friends, Alice and Bob, who are born at a festival.

• If they are born at 9:00 AM, they have 12 hours to wander off in different directions. By 9:00 PM, they might be on opposite sides of the city. (Wide spread).
• If they are born at 8:50 PM, just before the festival closes, they are still holding hands when the gates shut. (Narrow spread).

What the Scientists Found

The researchers used a computer program called Pythia8 + Angantyr to simulate this party. They compared their computer simulation to real data from the ALICE experiment at CERN. Here is what they discovered about the different types of particles:

1. Pions (The Lightweights)

• Observation: In real life, pions get closer together (the balance function gets narrower) as the collision gets more "central" (more violent).
• The Analogy: In a huge, violent crash (Central collision), the "party" lasts longer, and the radial flow (the outward push of the crowd) is stronger. This pushes the late-born pion couples closer together before they scatter.
• The Problem: The computer model (Pythia) was good at simulating small, weak crashes (Peripheral collisions), but it failed to get the details right for the big, violent crashes (Central collisions). It couldn't reproduce the "narrowing" effect seen in real life.
• The Dip: The real data shows a tiny "dip" (a hole) right in the center where the couples are born. This is caused by two things:
  1. Resonances: Some pions are born from the decay of unstable "parent" particles (like a $\rho$ meson exploding into two pions).
  2. Bose-Einstein Correlations: A quantum effect where identical particles (like two pions) avoid being in the exact same spot.
     The model could only see this "dip" in small crashes, not big ones, because the computer model has a limit on how big the "room" can be.

2. Kaons (The "Strange" Dancers)

• Observation: Kaons contain "strange" quarks, which are created very early in the collision.
• The Analogy: Because they are born at the very start of the party, they have plenty of time to drift apart. Their "spread" is wide, and it doesn't change much whether the crash is big or small. They are like early arrivals who have already wandered off by the time the main event starts.
• The Result: The computer model worked quite well for kaons! It successfully predicted that their spread stays the same regardless of the crash size.

3. Protons (The Heavyweights)

• Observation: Protons (and their anti-protons) also showed a spread that didn't change much with the size of the crash.
• The Analogy: Like the kaons, protons seem to be born early in the process. They don't get squeezed closer together in big crashes the way pions do.
• The Result: The model worked best for protons when they turned off a specific feature called "Color Reconnection" (a complex rule about how particle strings tie together). Without this feature, the model matched the real data perfectly.

The Conclusion: A Good Start, But Needs Tuning

The paper concludes that the Pythia8 + Angantyr model is a great tool for simulating heavy-ion collisions, but it's not perfect yet.

• For small crashes (Peripheral): The model is excellent. It predicts how particles behave very accurately.
• For big crashes (Central): The model struggles. It needs a "dedicated tuning" (like a mechanic adjusting the engine specifically for a race car) to handle the extreme conditions of a central collision.

In Summary:
The scientists used a computer simulation to watch how particle couples separate after a nuclear crash. They found that light particles (pions) get squeezed together in big crashes, while heavy or "strange" particles (kaons and protons) stay spread out. The computer model is a great "practice run" for small crashes, but to perfectly simulate the biggest, most violent crashes, the physicists need to tweak the code to better understand the complex rules of the quantum dance floor.

Technical Summary

1. Problem Statement

The study addresses the challenge of modeling Balance Functions (BF) in heavy-ion collisions. The BF is a critical observable used to measure charge-anticharge correlations, providing insights into:

• Particle production time: Early production leads to broader BFs due to diffusion and collective expansion, while late production results in narrower BFs.
• Quark-Gluon Plasma (QGP) lifetime: The width of the BF serves as a "clock" for hadronization.
• Collective effects: Specifically, radial flow, which tightens correlations and narrows the BF.

While experimental data (e.g., from ALICE) shows distinct centrality dependencies for different particle species (pions, kaons, protons), standard event generators often struggle to reproduce these features quantitatively, particularly in central collisions. The authors aim to evaluate the capability of the Pythia8.3 + Angantyr model (which lacks a dense medium description) to reproduce these experimental trends and to disentangle the effects of resonance decays and Bose-Einstein correlations (BEC).

2. Methodology

• Event Generator: The study utilizes Pythia8.3 with the Angantyr model, designed for heavy-ion collisions (Pb–Pb) at $\sqrt{s_{NN}} = 2.76$ TeV. This model relies on the Lund string fragmentation mechanism and includes soft/hard scattering, Initial/Final State Radiation (ISR/FSR), and Multi-Parton Interactions (MPI).
• Configurations:
  • Tuning: The Monash 2013 tune was used.
  • Variables: Simulations were run with and without Color Reconnection (CR) and Bose-Einstein Correlations (BEC) to isolate their effects.
  • Resonance Handling: Specific decay channels were isolated or removed to study resonance contributions:
    • Pions: $\rho^0 \to \pi^+\pi^-$ and $\omega \to \pi^+\pi^-$.
    • Kaons: $\phi \to K^+K^-$.
  • Data Volume: 300 million events were generated.
• Analysis:
  • Balance functions were calculated for $\pi\pi$, $KK$, and $pp$ pairs.
  • Correlations were measured in terms of pseudo-rapidity difference ($\Delta\eta$ or $\Delta y$) and azimuthal angle difference ($\Delta\phi$).
  • Results were compared against ALICE experimental data across various centrality classes (0–5% to 70–90%).
  • Kinematic cuts were applied to match ALICE acceptance (e.g., $0.2 < p_T < 2.0$ GeV/c for pions/kaons).

3. Key Contributions

• Systematic Evaluation of Angantyr: A comprehensive assessment of the Pythia8.3 + Angantyr model's ability to describe identified particle balance functions in Pb–Pb collisions without a hydrodynamic medium description.
• Disentanglement of Physical Mechanisms: The study explicitly separates the effects of resonance decays and BEC on the shape and width of the balance functions for different particle species.
• Centrality Dependence Analysis: Detailed comparison of how the model handles the transition from peripheral to central collisions for pions, kaons, and protons.

4. Key Results

A. Pions ($\pi\pi$)

• Centrality Trend: The experimental data shows pion BFs become narrower from peripheral to central collisions. The model qualitatively reproduces this trend when Color Reconnection (CR) is included, but fails quantitatively in central collisions (0–5%).
• Resonance Effects: A "dip" observed at $\Delta y = 0$ and $\Delta \phi = 0$ in the data is attributed to resonance decays ($\rho, \omega$). When these contributions are removed in the simulation, the dip disappears, and the BF becomes narrower.
• BEC Effects: Including BEC allows the model to reproduce the dip at $\Delta y = 0$ for peripheral collisions. However, the model fails to reproduce this dip for central collisions, likely because the Angantyr model cannot simulate the large system size ($\sim 5$ fm) required for strong BEC effects in central events (limited to $\sim 2$ fm in the model).
• Conclusion: The Monash 2013 tune describes peripheral collisions well but requires a dedicated heavy-ion tuning to accurately describe central Pb–Pb data.

B. Kaons ($KK$)

• Centrality Trend: Unlike pions, kaon BFs show little to no dependence on centrality (remains broad). This is consistent with the theory that strange quarks are produced early in the collision, making their correlations less sensitive to the subsequent evolution of the medium.
• Model Performance: The model with CR agrees well with ALICE data across all centrality classes.
• Resonance Effects: The $\phi$ meson decay dominates the kaon BF. Removing $\phi$ contributions significantly reduces the yield and widens the BF.
• BEC Effects: Similar to pions, including BEC helps reproduce the dip at $\Delta y = 0$ for mid-to-peripheral centralities (30–90%), but the model struggles with the centralmost collisions due to system size limitations.

C. Protons ($pp$)

• Centrality Trend: Proton BFs are independent of centrality, suggesting baryon-antibaryon pairs are produced early.
• Model Performance: Interestingly, the model without Color Reconnection provides a better description of the proton data than the model with CR.
• Dip Feature: Experimental data shows a dip at $\Delta y = 0, \Delta \phi = 0$ attributed to baryon-antibaryon interactions (annihilation/regeneration). Since Angantyr lacks hadronic rescattering and medium effects, it cannot reproduce this dip.

D. Widths and Integrals

• Pions: The width ($\sigma_{\Delta y}$) decreases from peripheral to central collisions in data. The model shows a weak dependence with CR but almost none without it.
• Kaons & Protons: The widths remain nearly constant across centrality in both data and model, consistent with early production.
• Integrals: The integrated BF for pions decreases from central to peripheral collisions in data; the model captures this trend only weakly with CR.

5. Significance and Conclusion

The paper concludes that while Pythia8.3 + Angantyr successfully captures the qualitative features of balance functions (such as the narrowing of pion BFs and the centrality independence of kaons/protons), it has quantitative limitations in central heavy-ion collisions.

• Limitations: The lack of a dense medium description and the inability to scale the system size beyond 2 fm in the model prevent accurate reproduction of BEC-induced dips and the precise width of BFs in central collisions.
• Implication: To fully describe central Pb–Pb collisions, a dedicated heavy-ion tuning of the Angantyr framework is necessary.
• Physical Insight: The study reinforces the understanding that particle species with different quark contents (light vs. strange) and production times exhibit distinct centrality dependencies, and that resonance decays and quantum statistics (BEC) are crucial for explaining the fine structure (dips) of the balance functions.