Investigating the elliptic anisotropy of identified particles in p–Pb collisions with a multi-phase transport model

NUCLEAR PHYSICS AND INTERDISCIPLINARY RESEARCH

Investigating the elliptic anisotropy of identified particles in p–Pb collisions with a multi-phase transport model

Si-Yu Tang，

Liang Zheng，

Xiao-Ming Zhang ，

Ren-Zhuo Wan

Nuclear Science and Techniques

Vol.35, No.2

Article number 32

Published in print Feb 2024

Available online 28 Mar 2024

DOI：10.1007/s41365-024-01387-4

52405

The elliptic azimuthal anisotropy coefficient (v₂) of the identified particles at midrapidity (|η|<0.8) was investigated in p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV using a multi-phase transport model (AMPT). The calculations of differential v₂ based on the advanced flow extraction method of light flavor hadrons (pions, kaons, protons, and Λ) in small collision systems were extended to a wider transverse momentum (p_T) range of up to 8 GeV/c for the first time. The string- melting version of the AMPT model provides a good description of the measured p_T-differential v₂ of the mesons but exhibits a slight deviation from the baryon v₂. In addition, we observed the features of mass ordering at low p_T and the approximate number-of-constituent-quark (NCQ) scaling at intermediate p_T. Moreover, we demonstrate that hadronic rescattering does not have a significant impact on v₂ in p–Pb collisions for different centrality selections, whereas partonic scattering dominates in generating the elliptic anisotropy of the final particles. This study provides further insight into the origin of collective-like behavior in small collision systems and has referential value for future measurements of azimuthal anisotropy.

Azimuthal anisotropySmall collision systemsTransport model

Introduction

The main goal of heavy-ion collisions at ultrarelativistic energies is to explore the deconfined state of strongly interacting matter created at a high energy density and temperature, known as the gluon plasma (QGP) [1, 2]. An important observation for investigating the transport properties of the QGP is anisotropic flow [3, 4], which is quantified by the flow harmonic coefficients vn obtained from the Fourier expansion of the azimuthal distribution of the produced particles [5, 6]: $\frac{d N}{d φ} \propto 1 + 2 \sum_{n = 1}^{\infty} v_{n} cos [n (φ - Ψ_{n})],$ (1) where φ is the azimuthal angle of the final-state particle angle and Ψn is the symmetry plane angle in the collision for the n-th harmonic [7, 8]. The second-order coefficient v₂, referred to as the elliptic flow, is derived from the initial state spatial anisotropy of the almond-shaped collision overlap region that is propagated to the final state momentum space. The magnitude of the elliptic flow is sensitive to the fundamental transport properties of the fireball, such as the temperature-dependent equation of state and the ratio of shear viscosity to entropy density (η/s) [9, 10].

Over the past few decades, various measurements of elliptic flow in heavy-ion collisions performed at the relativistic heavy-ion collider (RHIC) [11-14] and the Large Hadron Collider (LHC) [15-18] have helped build a full paradigm of the strongly coupled QGP. Comprehensive measurements of p_T-differential elliptic flow of the identified particles were conducted by the ALICE Collaboration [19, 20]. The observed mass-ordering effect (i.e., heavier particles have a smaller elliptic flow than lighter particles at the same p_T) at low p_T is well described by hydrodynamic calculations and is attributed to the radial expansion of the QGP [21]. At intermediate p_T, the grouping of v₂ of mesons and baryons was observed, with mesons exhibiting less v₂ than baryons. These behaviors can be explained by the hypothesis that baryons and mesons have different production mechanisms through quark coalescence, which has been further investigated using the number-of-constituent-quark (NCQ) scaling [22-25]. Interestingly, such flow-like phenomena have been observed in small-collision systems. Long-range double-ridge structures were first measured in high-multiplicity pp and p–Pb collisions by the ALICE, ATLAS, and CMS collaborations [26-28]. The measurement of elliptic and triangular azimuthal anisotropies in central ³He+Au, d+Au, and p+Au collisions performed by the STAR Collaboration [29] suggests that sub-nucleon fluctuations also play an important role in influencing the flow coefficients in these small collision systems. In addition, these measurements were extended to the identified particles associated with the discovery of a significant positive v₂ [30, 31]. The observed particle-mass dependence of v₂ is similar to that measured in heavy-ion collisions [30]; however, the origin of such collective-like behavior remains unclear. Several theoretical explanations relying on either the initial state or final state effects have been proposed to understand the origin of azimuthal anisotropies in small systems. Studies that extend hydrodynamics from large to small systems based on final-state effects can well describe v₂ of soft hadrons [32-36]; however, they are based on the strong assumption that there is sufficient scattering among constituents in small systems. Hydrodynamics combined with the linearized Boltzmann transport (LBT) model can also describe the identified particle v₂ in a high-multiplicity small-collision system at an intermediate p_T [37]. Color-Glass Condensate (CGC)-based models and IP-Glasma models that consider the effect of momentum correlations in the initial state can quantitatively describe some features of collectivity in p–Pb collisions [38, 39], but without clear conclusions, particularly regarding the dependence on collision systems and rapidity.

In addition, an approach called parton escape shows that few scatterings can also create sufficient azimuthal anisotropies, which have been investigated using multiphase transport (AMPT) [40, 41]. The v₂ values of light hadrons measured in p–Pb collisions are well described in AMPT, where the contribution of anisotropic parton escape rather than hydrodynamics plays an important role [41]. In this study, we extend the AMPT calculations of the p_T-differential v₂ for identified particles (π^±, K^±, $p (\bar{p})$ , $Λ (\bar{Λ})$ ) to higher p_T region in p–Pb collisions at 5.02 TeV, in order to systematically test the mass-ordering effect and baryon-meson grouping at low- and intermediate-p_T, respectively. We also investigate how the key mechanisms implemented in AMPT, such as the parton cascade and hadronic rescattering, affect elliptic anisotropy in small collision systems. In addition, various nonflow subtraction methods with different sensitivities to jet-like correlations were studied.

Model and Methodology

2.1

A multiphase transport model

The string-melting version of the AMPT model (v2.26t9b, available online) [40, 42] was employed in this study to calculate v₂ of the final-state particles in high-multiplicity p–Pb at 5.02 TeV. The AMPT model includes four main processes: initial conditions, partial scattering, hadronisation, and hadronic interactions. The initial conditions are generated from the heavy ion jet interaction generator (HIJING) model [43, 44], where minijet partons and soft-excited strings are produced and then converted to primordial hadrons based on Lund fragmentation. Under the string-melting mechanism, primordial hadrons are converted into partons, a process determined by their flavor and spin structures. Elastic scattering between the partons was simulated using Zhang’s parton cascade (ZPC) model [45], which includes two-body scattering with a cross-section described by the following simplified equation: $σ_{g g} \approx \frac{9 π α_{s}^{2}}{2 μ^{2}} .$ (2) In this study, the strong coupling constant α_s was set to 0.33, and the Debye screening mass μ = 2.2814 fm^-1, resulting in a total parton scattering cross section of σ= 3 mb. To isolate the effect of partonic scattering, σ is adjusted to be close to 0 by increasing μ (see set "w/o parton scat." in Tab. 1). Once the partonic interaction ceases, hadronization with a quark coalescence model is implemented to combine the nearest two (or three) quarks into mesons (or baryons) [40]. The formed hadrons enter the subsequent hadronic rescattering process using a relativistic transport (ART) model [46], in which both elastic and inelastic scattering are considered for baryon-baryon, baryon-meson, and meson-meson interactions. The hadronic interaction time was set by default to t_max = 30 fm/c. Alternatively, t_max is set to 0.4 fm/c to effectively turn off the hadron scattering process while still considering the resonance decay [47](see set "w/o hadron scat." in Tab. 1). In addition, the random orientation of the reaction plane was turned on and the shadowing effect was considered in this analysis.

Details of three configurations

Description	σ (mb)	t_max (fm/c)
w/ all	3	30
w/o parton scat.	~0	30
w/o hadron scat.	3	0.4

2.2

Two-particle correlation and nonflow subtraction

The two-particle correlation (2PC) method is widely used to extract the flow signal in small collision systems because it can suppress the non-flow contribution from long-range jet correlations [26-28, 30, 48]. Similar to Eq. 1, the azimuthal correlation between two emission particles can be represented by N^pairs pairs of emitted particles (labeled as C(Δφ)) as a function of the relative angle Δφ=φa-φb between particles a and b and expanded in the Fourier series as follows: $\begin{array}{l} C (Δ φ) = \frac{d N^{pair}}{d Δ φ} \propto 1 + 2 \sum_{n = 1}^{\infty} V_{n Δ} (p_{T}^{a}, p_{T}^{b}) cos [n (Δ φ)], \end{array}$ (3) where V_nΔ refers to the two-particle n-th order harmonic. In a pure hydrodynamic scenario, because particle emissions are independent, $V_{n Δ} (p_{T}^{a}, p_{T}^{b})$ can be factorized into the product of a single-particle flow $v_{n}^{a}$ and $v_{n}^{b}$ : $\begin{array}{l} V_{n Δ} (p_{T}^{a}, p_{T}^{b}) = v_{n} (p_{T}^{a}) v_{n} (p_{T}^{b}) . \end{array}$ (4) Based on the factorization assumption, vn of a single particle a can be obtained using the 3×2PC method, which was recently proposed by the PHENIX Collaboration [49]. This requires the formation of two-particle correlations between three groups of particles (labeled a, b and c) and the extraction of the flow coefficients for three combinations: $\begin{array}{l} v_{n} (p_{T}^{a}) = \sqrt{\frac{V_{n Δ} (p_{T}^{a}, p_{T}^{b}) V_{n Δ} (p_{T}^{a}, p_{T}^{c})}{V_{n Δ} (p_{T}^{b}, p_{T}^{c})}} . \end{array}$ (5) In small-collision systems, two main types of nonflow contributions to the flow signal are the near-side jet and away-side jet (recoil jet) correlations. The former can be effectively removed by introducing a large rapidity gap between the trigger and associated particles during the construction of the correlations. Several methods have been developed to subtract the latter [50]. A traditional approach is to directly subtract the correlation function distribution obtained from low-multiplicity events [27, 30] from that obtained from high-multiplicity events. This method assumes that the yield and shape of dijets are identical for both collision types as follows: $\begin{matrix} C^{HM} (Δ φ) - C^{LM} (Δ φ) \propto 1 + 2 \sum_{n = 1}^{\infty} V_{n Δ} cos [n (Δ φ)] \\ = a_{0} + 2 \sum_{n = 1}^{\infty} a_{n} cos [n (Δ φ)], \end{matrix}$ (6) where C^LM(Δφ) and C^HM(Δφ) represent the correlation function distributions obtained for low- and high-multiplicity events, respectively. This method relies on the "zero yield at minimum" (ZYAM) hypothesis [27, 30] that a flat combinatoric component should be subtracted from the correlation function in low-multiplicity events. Therefore, the fit parameter a₂ is the absolute modulation in the subtracted correlation function distribution and characterizes the modulation relative to a baseline, assuming that such a modulation is not present in the low-multiplicity class below the baseline. In this case, the flow coefficient V_nΔ is calculated as $\begin{array}{l} V_{n Δ} = a_{n} / (a_{0} + b), \end{array}$ (7) where b is the baseline, estimated using the minimum correlation function for low-multiplicity events. However, the measurement of jet-like correlations in p–Pb collisions indicates that the dependence of the dijet yield on the particle multiplicity cannot be ignored. In this case, a new template fit method was developed by the ATLAS collaboration [51], where the correlation function distribution obtained in high-multiplicity events is assumed to result from the superposition of the distribution obtained in low-multiplicity events scaled up by a multiplicative factor F and a constant modulated by cos(nΔφ) for n > 1, as shown in $\begin{array}{l} C (Δ φ) = F C^{LM} (Δ φ) + G (1 + 2 \sum_{n = 1}^{3} V_{n Δ} cos (n Δ φ)), \end{array}$ (8) where G denotes the normalization factor that maintains the integral of C(Δφ) equal to C^HM(Δφ). Furthermore, an improved template fitting method [52] developed in recent years was tested. It applies a correction procedure to the default template fit method by considering the multiplicity dependence of the remaining ridge in low-multiplicity events, as shown in $\begin{array}{l} V_{n Δ} = V_{n Δ} (tmp) - \frac{F G^{LM}}{G^{HM}} (V_{n Δ}^{2} (tmp) - V_{n Δ}^{2} (LM)), \end{array}$ (9) where V_nΔ(tmp) and $V_{n Δ}^{2} (LM)$ are obtained by using the default template method for high- and low-multiplicity events. All these nonflow subtraction methods are implemented in this study, and their different sensitivities to nonflow effect are also discussed.

Analysis procedures

To directly compare the AMPT calculations with the results from ALICE, we focused on the particles within the pseudorapidity range |η| < 0.8, aligning with the TPC acceptance in ALICE [53]. In the 3×2PC method, long-range correlations were constructed between the charged particles at mid-rapidity, forward rapidity (2.5 < η < 4), and backward rapidity (-4 < η < -2.5), that is, the central-forward correlation (-4.8 < Δη < -1.7), central-backward correlation (1.7 < Δη < 4.8), and backward-forward correlations (-8 < Δη < -5). In addition, the centrality classes are defined by counting the charged particles in the acceptance of the V0A detector [53], that is, 2.8<η<5.1.

The correlation function distribution C(Δφ) was obtained by correcting the number of particle pairs in the same events normalized to the number of trigger particles N_trig by using an event-mixing technique: $\begin{array}{l} C (Δ φ, Δ η) = \frac{1}{N_{trig}} \frac{d^{2} N_{pairs}}{d Δ η d Δ φ} = \frac{S (Δ φ, Δ η)}{B (Δ φ, Δ η)}, \end{array}$ (10) where $S (Δ φ, Δ η) = \frac{1}{N_{trig}} \frac{d^{2} N_{same}}{d Δ η d Δ φ}$ is the correlation function in same events and $B (Δ φ, Δ η) = α \frac{d^{2} N_{mixed}}{d Δ η d Δ φ}$ is the associated yield as a function of Δφ and Δφ in mixed events. Factor α is used to normalize B(Δφ, Δη) to unity in the Δη region of the maximal pair acceptance. The obtained 2-D correlation function C(Δφ, Δη) is projected onto Δφ axis, and we follow the nonflow subtraction procedures and factorization, as discussed in Eq. 5–Eq. 9, v₂ of the charged particles at |η| < 0.8 can be calculated.

Results and Discussions

We first investigated the p_T spectrum of the identified particles before performing the flow analysis. Figure 1 illustrates the p_T distribution of proton, pion, and kaon in 0–20% high-multiplicity p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV, which are obtained from AMPT with three different sets of configurations listed in Table 1 and ALICE experimental data [54]. The AMPT results, both with and without hadronic rescattering, are consistent. This behavior differs from previous findings in heavy-ion collisions, where the hadronic interaction significantly reduces the particle yield [55]. The spectrum obtained in the AMPT without considering the parton cascade process is enhanced compared to that obtained with partonic scattering, and this enhancement is more significant at a high p_T. This outcome is expected because partons experience energy loss during the parton cascade, which reduces the production of final-state particles. In addition, the ratios of the p_T spectra obtained from the AMPT calculations and data are shown. The AMPT model calculation reproduces the particle yields well at low and intermediate p_T values when both partonic and hadronic scattering are included; however, it overestimates the high p_T data because parton-parton inelastic collisions and, subsequently, hard parton fragmentation are absent in the model [42].

Fig. 1

(Color online) The p_T distribution of pions, kaons, and protons in 0–20% high-multiplicity p–Pb collisions at

\sqrt{s_{NN}}

= 5.02 TeV, obtained from AMPT model calculations, is compared to ALICE measurement [54]. The results in AMPT without hadronic scattering and partonic scattering are also presented

Figure 2 (left) shows the v₂ of pions, kaons, protons and Λ as a function of p_T in 0–20% high-multiplicity p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV, obtained in AMPT calculations with 3×2PC method. A comparison with the ALICE measurement for v₂ of charged hadrons, pions, kaons, and protons [30] and the CMS measurement for v₂ of $K_{s}^{0}$ and Λ [31] is also presented. The AMPT calculations applied the template fit method to suppress the away-side jet contribution and considered the ZYAM assumption to enable direct comparison with the observed data. The v₂ values of charged hadrons, pions, and kaons can be described well by AMPT calculations, but the v₂ values of baryons (protons and Λ) cannot be reproduced. In addition, the mass-ordering effect (i.e., the v₂ of baryons is lower than that of mesons) is reproduced for p_T <2 GeV/c. Owing to the advanced flow extraction method, the calculations of v₂ were extended to the high-p_T region, up to 8 GeV/c in the AMPT model for the first time. The v₂ values of protons and Λ are consistent, and both of them are observed to have a higher value of v₂ value than that of the mesons for 2<p_T< 7 GeV/c. The observed meson-baryon particle type grouping in heavy-ion collision flow measurements indicates collective behavior at the partonic level, leading to the coalescence of quarks into hadrons. The number-of-constituent-quarks (NCQ) scaling techniques described in [22] can be used for further studies of this grouping. v₂ and p_T in Fig. 2 (left) are replaced by v₂/n_q and p_T/n_q, where the n_q is the number of constituent quark in mesons (n_q = 2) and baryons (n_q = 3), as shown in Fig. 2 (right). v₂/n_q obtained from the data show approximate values at intermediate p_T; however, the results calculated in AMPT cannot reproduce the scaling in p_T/n_q>1 GeV/c. In order to consider the observed mass hierarchy of v₂, we also plot the v₂ of identified particle as a function of the transverse kinetic energy kE_T ( $k E_{T} = m_{T} - m_{0} = \sqrt{p_{T}^{2} + m_{0}^{2}} - m_{0}$ ), and its NCQ scaling in Fig. 3 (left), and Fig. 3 (right). All particle species showed a set of similar v₂ values after NCQ scaling in kE_T/n_q < 1 GeV, confirming that the quark degree of freedom in flowing matter can also be probed in the transport model. However, this NCQ scaling is violated for kE_T/n_q > 1 GeV. This may be attributed to the hadronization mechanism implemented in the AMPT model used in this study, where baryons are produced only after the formation of mesons by simply combining the three nearest partons, regardless of the relative momentum among the coalescing partons. This results in an underestimation of the baryon v₂ at intermediate p_T in this study. An improved coalescence model implemented in the newer AMPT [56] introduced a new coalescence parameter to control the relative probability of a quark forming a baryon instead of a meson precisely. This improvement could have different NCQ scaling on v₂ but requires more systematic studies. Further studies on v₂ calculations in small collision systems with other improved hadronization mechanisms, for example, considering the Wigner function [57] and hard parton fragmentation [58], should be performed in the future.

Fig. 2

(Color online) Left: the v₂ as a function of p_T in 0–20% high-multiplicity p–Pb collisions at

\sqrt{s_{NN}}

= 5.02 TeV, obtained from default AMPT model calculations with 3×2PC method, is compared to ALICE and CMS measurement [30, 31]. Right: the p_T-differential v₂ scaled by the number of constituent quark (n_q)

Fig. 3

(Color online) Left: the v₂ as a function of transverse kinetic energy (kE_T) in 0–20% high-multiplicity p–Pb collisions at

\sqrt{s_{NN}}

= 5.02 TeV, obtained from default AMPT model calculations with 3×2PC method, is compared to ALICE and CMS measurement [30, 31]. Right: the kE_T-differential v₂ scaled by the number of constituent quark (n_q)

We also extend our investigation to include a study of integrated v₂ within various centrality bins spanning the 0–60% range. We focus on the region where the NCQ scaling criterion is satisfied, that is, for transverse kinetic energies per constituent quark (kE_T/n_q) ranging from 0.4 to 1 GeV. The non-flow contribution was estimated and subtracted within the 60–100% centrality class by using the template fit method. As shown in Fig. 4, the v₂ values as a function of centrality exhibit a systematic decrease from central to peripheral collisions, reflecting the changing dynamic conditions and particle production mechanisms in different collision zones. Intriguingly, in the v₂ measurements, we observed a distinct mass-splitting phenomenon, with baryons and mesons exhibiting distinct elliptic flow patterns. Such a mass dependence in v₂ is similar to that in heavy-ion collisions at the LHC energies presented in a previous study [47]. This provides valuable insights into the collective behavior of different particle species within the evolving fireball created during these collisions.

Fig. 4

(Color online) The integrated v₂ in 0.4<kE_T/n_q<1 GeV for pion, kaon and proton varying with the centrality

Moreover, to gain a deeper understanding of the NCQ scaling properties, we explored the ratios of n_q-scaled integrated v₂ values for protons relative to pions and kaons relative to pions as functions of centrality. The results are shown in Fig. 5. A notable trend is observed in these ratios: they tend to approach unity as the collisions become more peripheral. It indicates that the collective flow of particles in low-multiplicity events may be approaching a behavior that is closer to the expected scaling behavior based on the number of constituent quark.

Fig. 5

(Color online) The ratio of integrated v₂ within 0.4<kE_T/n_q<1 GeV for proton over pion and kaon over pion varying with the centrality. The dash line represents the location of unity ratio

The effects of partonic and hadronic scattering on the elliptical anisotropy of the final-state particles were examined in this study. Figure 6 shows the calculated p_T-differential v₂ of pions, kaons, and protons in AMPT with and without considering hadronic rescattering process in 0—20% high-multiplicity p-Pb collisions. The results show that the ratio of the v₂ values with and without hadronic rescattering is consistent with unity for all particle species, indicating that the hadronic rescattering mechanism has almost no effect on v₂ in high-multiplicity p—Pb collisions. We also investigated the centrality dependence of the hadronic rescattering effects by calculating p_T-integrated v₂ in several centrality bins between 0 and 60%, as illustrated in Fig. 7. The results demonstrate that the influence of hadronic rescattering is independent of the centrality selection and has almost no impact on NCQ scaling in the range of 0.4<kE_T/n_q<1 GeV.

Fig. 6

(Color online) The p_T-differential v₂ of pions, kaons, and protons calculated in AMPT model with and without considering hadronic scattering. The ratios of the two sets are also presented

Fig. 7

(Color online) The integrated v₂ in 0.4<kE_T/n_q<1 GeV for pions, kaons, and protons calculated in AMPT model with and without considering hadronic scattering. The ratios of the two sets are also presented

On the other hand, when we set the parton scattering cross-section σ to zero but maintain the hadronic scatterings, the V_2Δ of charged particles for the central-forward (CF) and central-backward (CB) correlations is almost zero, as shown in Fig. 8. If both the partonic and hadronic scatterings are turned off, the results remain consistent with zero. This indicates that the elliptical anisotropy in high-multiplicity small-collision systems is mostly generated by parton scattering. Our conclusion is consistent with previous studies on the AMPT [41], which suggested that the majority of elliptic anisotropies comes from the anisotropic escape probability of partons.

Fig. 8

(Color online) The p_T-differential V_2Δ for central-forward (CF) and central-backward (CB) correlations calculated in AMPT model with and without considering partonic scattering

Finally, different non-flow subtraction methods were investigated in this study. Figure 9 (left) shows the p_T-differential v₂ of the charged particles calculated using the 3×2PC method in 0—20% high-multiplicity p—Pb collisions. Several nonflow subtraction methods are implemented. To demonstrate how the nonflow contribution is removed, v₂ obtained with a direct Fourier transform of the C(Δφ) correlation (as shown in Eq. 3). The results show significant suppression across all the subtraction methods, particularly at higher p_T values where jet correlations are dominant. The results obtained with peripheral subtraction and template fitting were consistent, indicating that the away-side jet contribution was automatically removed using the 3×2PC method, even though the dependence of the jet correlation on multiplicity was not considered in the peripheral subtraction method. The v₂ calculated using the improved template fit method was slightly lower than that from the template fit, and it was similar to the features observed in the ATLAS measurement [52]. The same conclusions were drawn for the extraction of the identified particles (pions, kaons, protons, and Λ) v₂.

Fig. 9

(Color online) The p_T-differential v₂ of charged hadrons calculated in AMPT with different nonflow subtraction methods

Summary

This study systematically investigated the elliptic anisotropy of identified particles (pions, kaons, protons, and Λ) in p—Pb collisions at 5.02 TeV using the AMPT model. We extended the calculation of v₂ to higher p_T regions, up to 8 GeV/c, using advanced nonflow subtraction techniques for the first time. We also examined the mass-ordering effect and baryon-meson grouping at low and intermediate p_T, respectively. We argue that, with the approximate NCQ scaling of baryons and mesons, v₂ can be reproduced well at kE_T/n_q<1 GeV for several centrality bins. Furthermore, we demonstrate that parton interactions can simultaneously decrease the yield of light hadrons and generate significant v₂. However, hadronic rescatterings had little influence on the elliptical anisotropy of the final-state particles. Thus, these findings indicate that the nonequilibrium anisotropic parton escape mechanism coupled with the quark coalescence model can also reproduce the hydro-like behavior of the identified particles observed in small collision systems. Overall, this study provides new insights into the existence of partonic collectivity in small collision systems.

References

E. V. Shuryak,

Quark-Gluon Plasma and Hadronic Production of Leptons, Photons and Psions

. Phys. Lett. B 78, 150 (1978). https://doi.org/10.1016/0370-2693(78)90370-2