Elliptic flow of charged hadrons in d+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV using a multi-phase transport model

This study utilizes the AMPT model to demonstrate that the elliptic flow of charged hadrons in 200 GeV d+Au collisions is primarily driven by early-stage partonic interactions rather than hadronic rescattering, successfully reproducing experimental trends and highlighting the significance of parton scattering mechanisms in asymmetric collision systems.

The Gist

Imagine you are at a massive, chaotic concert where thousands of people are packed into a stadium. Now, imagine two different scenarios:

1. The Symmetric Clash: Two identical, perfectly round crowds crash into each other head-on.
2. The Asymmetric Crash: A small, dense group (a deuteron, which is like a tiny cluster of two people) smashes into a massive, round crowd (a gold nucleus, like a stadium full of people).

This paper is about the second scenario. Scientists are trying to understand what happens when these two very different groups collide at nearly the speed of light. Specifically, they are looking at how the "people" (particles) in the aftermath move in a specific, organized pattern called Elliptic Flow.

Here is a simple breakdown of what the researchers did and what they found, using everyday analogies.

The Big Question: Is it a "Soup" or a "Gas"?

When these particles collide, they create a super-hot, super-dense state of matter called Quark-Gluon Plasma (QGP). Think of this as a "perfect fluid" or a cosmic soup where the building blocks of matter (quarks and gluons) melt together and flow like water.

In a perfect, round collision (like two gold nuclei), this soup expands and squirts out in a specific oval shape (like squeezing a water balloon). This is "elliptic flow." But in an asymmetric collision (Deuteron + Gold), the shape is weird and lopsided. The big question was: Does this weird shape still create that organized "flow" pattern, or does it just scatter chaotically?

The Tool: The "Traffic Simulator"

To answer this, the researchers didn't just look at the data; they built a digital traffic simulator called the AMPT model.

Imagine you are a traffic engineer trying to predict how cars will move after a massive pile-up. You have to decide:

• Do the cars bounce off each other like billiard balls?
• Do they stick together and flow like a river?
• How long do they interact before they stop?

The researchers ran their simulation with two different "rulesets":

1. Default Mode: The particles act mostly like solid objects that bounce around.
2. String Melting Mode: The particles melt into a fluid-like soup first, then re-form. This is the "perfect fluid" scenario.

The Experiments: Changing the Rules

The team tweaked the settings in their simulator to see what changed the flow:

1. The "Stickiness" of the Particles (Cross-Section)
They changed how much the particles "bump" into each other.

• Analogy: Imagine the particles are people in a mosh pit. If they are very "sticky" (high cross-section), they bump into each other constantly, transferring energy and moving together in a wave. If they are "slippery" (low cross-section), they just slide past each other.
• Finding: When they made the particles "stickier," the organized flow (elliptic flow) got much stronger. This suggests that for the flow to happen, the particles need to interact heavily, like a fluid.

2. The "Hang Time" (Hadronic Cascading)
They changed how long the particles kept interacting after the initial crash.

• Analogy: After the concert ends, do the crowd members immediately leave, or do they linger, bumping into each other in the lobby for a while?
• Finding: Surprisingly, how long they lingered didn't change the main flow pattern much. The "dance" was mostly choreographed in the very first split second of the collision.

The Results: What the Data Said

The researchers compared their simulator's output to real data from two famous experiments (STAR and PHENIX) at the Relativistic Heavy Ion Collider (RHIC).

• The "Event Plane" vs. The "Participant Plane":
  This is the trickiest part. Imagine trying to guess the direction a crowd is moving.

  • Participant Plane: You look at where the people started (the geometry of the crash).
  • Event Plane: You look at where the people ended up (the final pattern).
  • In these weird, small collisions, the starting point and the ending point don't always match up perfectly because of random fluctuations (like a few people tripping early).
  • Finding: The simulator matched the real data best when they looked at the "Participant Plane" (the starting geometry) and assumed the particles were very "sticky" (high interaction).

• No "Baryon-Meson" Split:
  In big, symmetric collisions, heavy particles (like protons) and light particles (like pions) separate into different lanes based on their mass. In these small d+Au collisions, the researchers found no such separation.

  • Analogy: In a big river, heavy rocks and light leaves might sort themselves out. In this tiny, chaotic splash, everything is mixed up together. This suggests the "soup" might be smaller or shorter-lived than in big collisions.

The Bottom Line

This paper tells us that even in a small, asymmetric crash (Deuteron hitting Gold), the particles behave like a fluid. They don't just bounce randomly; they push against each other and flow in an organized oval shape.

However, this flow is very sensitive to how much the particles interact. If they don't bump into each other enough, the flow disappears. The study confirms that the "perfect fluid" behavior isn't just for giant collisions; it can happen in smaller systems too, provided the particles interact strongly enough in that first, fleeting moment.

In short: Even a small, lopsided crash can create a perfectly fluid "dance," but only if the dancers are close enough to hold hands and move together.

Technical Summary

Here is a detailed technical summary of the paper "Elliptic flow of charged hadrons in d+Au collisions at $\sqrt{s_{NN}} = 200$ GeV using a multi-phase transport model."

1. Problem Statement

The study addresses the ongoing debate regarding the nature of collective flow in asymmetric collision systems (specifically d+Au) at the Relativistic Heavy Ion Collider (RHIC). While collective behavior (elliptic flow, $v_2$) is well-established in symmetric heavy-ion collisions (Au+Au) as evidence of Quark-Gluon Plasma (QGP) formation, its origin in smaller, asymmetric systems like d+Au is less clear.

• Key Question: Does the observed $v_2$ in d+Au collisions arise from hydrodynamic expansion driven by initial spatial geometry and partonic interactions, or is it dominated by initial-state momentum correlations and final-state hadronic effects?
• Specific Challenge: Experimental data from STAR and PHENIX show discrepancies in $v_2$ measurements depending on the analysis method (Event Plane vs. Participant Plane) and the treatment of "non-flow" effects. Theoretical models need to clarify the roles of early-stage partonic scattering versus late-stage hadronic rescattering in generating these anisotropies.

2. Methodology

The authors utilized the AMPT (A Multi-Phase Transport) model to simulate 100 million minimum-bias d+Au collision events at $\sqrt{s_{NN}} = 200$ GeV.

• Model Variants:
  • AMPT-DEF (Default): Partons recombine with parent strings; strings fragment into hadrons via the Lund model.
  • AMPT-SM (String Melting): Excited strings are "melted" into partons, which then undergo scattering before hadronizing via quark coalescence. This mode explicitly includes a partonic phase.
• Simulation Parameters:
  • Parton Scattering Cross-section ($\sigma_{qq}$): Varied between 1.5, 3.0, and 6.0 mb to study the sensitivity of $v_2$ to partonic interaction strength.
  • Hadronic Cascading Time ($t_{had}$): Varied between 0.6 and 30 fm/c to assess the impact of final-state hadronic interactions.
  • Centrality: Events were categorized into 0-10% (central), 10-40% (mid-central), and 40-80% (peripheral) based on charged particle multiplicity ($|\eta| < 0.9$).
• Analysis Techniques:
  • $v_2$ Estimation: Calculated using two distinct reference planes:
    1. Event Plane ($\psi_2^{ep}$): Reconstructed using the $\eta$-subevent method (separating particles into forward and backward rapidity regions to suppress non-flow correlations).
    2. Participant Plane ($\psi_2^{pp}$): Calculated event-by-event based on the transverse positions of participating nucleons, representing the true initial geometric symmetry.
  • Particle Identification: Analyzed inclusive charged hadrons ($h^\pm$) and identified species ($\pi^\pm, K^\pm, p/\bar{p}$).
  • Scaling: Examined $v_2$ as a function of transverse kinetic energy ($K_{ET}$) scaled by the number of constituent quarks ($n_q$) to test for quark coalescence signatures.

3. Key Contributions

• Systematic Parameter Variation: The study provides a comprehensive scan of how $\sigma_{qq}$ and $t_{had}$ influence $v_2$ in asymmetric systems, isolating the contributions of the partonic phase versus the hadronic phase.
• Methodological Comparison: It explicitly contrasts $v_2$ results derived from the Event Plane method against the Participant Plane method within the same model framework, offering insight into why experimental results from different collaborations (STAR vs. PHENIX) might differ.
• Asymmetry Analysis: It investigates the unique centrality dependence of $v_2$ in d+Au compared to symmetric Au+Au collisions, specifically looking at the reversal of trends in peripheral collisions.

4. Key Results

A. Transverse Momentum ($p_T$) Dependence

• Event Plane ($v_2^{ep}$): Shows a monotonic increase with $p_T$, indicating growing collective motion.
• Participant Plane ($v_2^{pp}$): Increases up to $p_T \approx 2.0$ GeV/c and then saturates.
• Magnitude: $v_2^{ep}$ is consistently larger than $v_2^{pp}$ across all particle types, attributed to event-by-event fluctuations in the final-state particle distribution used to define the event plane.

B. Centrality Dependence

• Low $p_T$ (< 2.0 GeV/c): $v_2^{ep}$ shows negligible centrality dependence.
• High $p_T$ (> 2.0 GeV/c): $v_2^{ep}$ increases from central to peripheral collisions, a trend similar to symmetric heavy-ion collisions.
• Participant Plane Trend: $v_2^{pp}$ exhibits a centrality dependence across all $p_T$ ranges but follows an opposite trend to symmetric systems (decreasing from central to peripheral), highlighting the significant role of initial geometry fluctuations in asymmetric collisions.

C. Particle Type and Scaling

• Mass Ordering: Weak mass ordering is observed at low $p_T$, but no clear baryon-meson separation is seen at intermediate $p_T$ (2–5 GeV/c).
• NCQ Scaling: The study finds no adequate Number of Constituent Quark (NCQ) scaling when plotting $v_2/n_q$ vs. $K_{ET}/n_q$. This contrasts with symmetric Au+Au collisions and suggests that the medium in d+Au collisions at this energy may be dominated by hadronic interactions or lacks sufficient partonic collectivity to exhibit quark coalescence scaling.

D. Sensitivity to Model Parameters

• Partonic Cross-section ($\sigma_{qq}$): Increasing $\sigma_{qq}$ significantly enhances the magnitude of $v_2$. A larger cross-section implies stronger partonic collectivity.
• Hadronic Cascading Time ($t_{had}$): Varying $t_{had}$ from 0.6 to 30 fm/c has a negligible effect on $v_2^{ep}$ but affects $v_2^{pp}$. This suggests that the early-stage partonic phase is the primary driver of the observed elliptic flow, while final-state hadronic interactions play a minor role in shaping the anisotropy.

E. Comparison with Experimental Data

• STAR Data (0-10% centrality): The AMPT-SM model with a high cross-section ($\sigma_{qq} \approx 6.0$ mb) and Participant Plane angle best reproduces the STAR data.
• PHENIX Data (0-5% centrality): The AMPT-DEF model and AMPT-SM with a low cross-section ($\sigma_{qq} \approx 1.5$ mb) using the Event Plane angle align well with PHENIX data.
• Discrepancy: The ~10% difference between STAR and PHENIX experimental results is attributed to different non-flow subtraction techniques. The model confirms that the choice of reference plane (Event vs. Participant) and the treatment of non-flow are critical in interpreting small-system data.

5. Significance

This study underscores that early-stage partonic scattering is the dominant mechanism driving elliptic flow in d+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The lack of NCQ scaling and baryon-meson splitting suggests that while collective behavior exists, the system may not fully thermalize into a hydrodynamic QGP in the same way as central Au+Au collisions.

Crucially, the paper highlights the sensitivity of $v_2$ measurements to the analysis method. The divergence between Event Plane and Participant Plane results, and the resulting discrepancies between experimental collaborations, emphasizes that in asymmetric systems, "non-flow" effects and initial geometry fluctuations must be carefully disentangled to draw valid conclusions about the formation of a QGP. The AMPT model serves as a vital tool for interpreting these complex dependencies, suggesting that a single "true" $v_2$ value may not exist without specifying the reference frame and interaction timescales.