Collective effects in light-heavy ion collisions

Collective effects in light-heavy ion collisions

Björn Schenke and Raju Venugopalan Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA

We present results for the azimuthal anisotropy of charged hadron distributions in A+A, p+A, d+A, and He+A collisions within the IP-Glasma+music model. Obtained anisotropies are due to the fluid dynamic response of the system to the fluctuating initial geometry of the interaction region. While the elliptic and triangular anisotropies in peripheral Pb+Pb collisions at are well described by the model, the same quantities in p+Pb collisions underestimate the experimental data. This disagreement can be due to neglected initial state correlations or the lack of a detailed description of the fluctuating spatial structure of the proton, or both. We further present predictions for azimuthal anisotropies in p+Au, d+Au, and He+Au collisions at . For d+Au and He+Au collisions we expect the detailed substructure of the nucleon to become less important.

heavy ion collisions, fluctuations, fluid dynamics

00 \journalnameNuclear Physics A \runauthB. Schenke and R.Venugopalan \jidnupha \jnltitlelogoNuclear Physics A

1 Introduction

For the physical interpretation of observables in heavy ion, proton/deuteron-heavy ion and (high-multiplicity) proton+proton collisions, the existence of a sophisticated description of the multi-particle production mechanism and event-by-event fluctuations is essential. In the high energy limit the color glass condensate (CGC) framework Gelis et al. (2010) is the proper effective theory of quantum chromodynamics (QCD) that provides such description.

One particular implementation of the CGC is the IP-Glasma model Schenke et al. (2012a, b). It combines the IP-Sat dipole model Kowalski and Teaney (2003); Kowalski et al. (2008), which parametrizes the impact parameter and -dependence of the saturation scale, with the classical dynamics of produced gluon fields Krasnitz and Venugopalan (2000, 2001); Lappi (2003). With its parameters constrained by inclusive and diffractive deeply inelastic scattering (DIS) data from e+p scattering at HERA, the IP-Glasma model correctly describes the bulk features of various collision systems over a wide range of energies Schenke et al. (2014).

After briefly introducing the IP-Glasma model and the employed relativistic fluid dynamic simulation music, we present results for the azimuthal anisotropy of charged hadrons produced in peripheral Pb+Pb collisions, and p+Pb collisions of comparable multiplicity. We compare to experimental data from the CMS collaboration Chatrchyan et al. (2013). We then present predictions for the systematics of the transverse momentum dependent anisotropy coefficients in p+Au, d+Au, and He+Au collisions at top RHIC energies.

2 IP-Glasma + MUSIC

The IP-Glasma model Schenke et al. (2012a, b) relates the DIS constrained nuclear dipole cross-sections to the initial classical dynamics of highly occupied gluon fields produced in a nuclear collision. Given an initial distribution of color charges in the high energy nuclear wave-functions, the strong multiple scatterings of gluon fields are computed by event-by-event solutions of Yang-Mills equations. Both fluctuating distributions of nucleons in the nuclear wave-functions and intrinsic fluctuations of the color charge distributions are included. This results in “lumpy” transverse projections of the gluon field configurations that vary event to event. The scale of this lumpiness is given on average by the nuclear saturation scale which corresponds to distance scales smaller than the nucleon size Kowalski et al. (2008).

The IP-Glasma model provides the initial conditions for fluid dynamic calculations at a given time . The initial energy density and flow velocities , are extracted from the gluon fields’ energy-momentum tensor at every transverse position via the relation . In the results presented below, the viscous part of the energy momentum tensor is set to zero at the initial time of the fluid dynamic simulation. This is done because the gluon field strength tensor is very anisotropic (the longitudinal pressure is approximately zero). A full 3+1 dimensional simulation including quantum fluctuations could provide a mechanism for isotropization via instabilities Berges et al. (2014); Epelbaum and Gelis (2013).

We employ the viscous relativistic fluid dynamic simulation music Schenke et al. (2010, 2011a, 2011b), which is a 3+1 dimensional simulation. However, because the initial conditions from the IP-Glasma model are boost-invariant, it is used in its 2+1 dimensional mode.

3 p+Pb collisions at the LHC

Figure 1: Multiplicity dependence of the root-mean-square elliptic (left) and triangular (right) flow coefficients in Pb+Pb (open symbols) and p+Pb collisions (filled symbols) from the IP-Glasma+music model (connected triangles) compared to experimental data by the CMS collaboration Chatrchyan et al. (2013). Figure from Schenke and Venugopalan (2014).

It was shown in Schenke and Venugopalan (2014) that experimental data from heavy ion collisions is well described by the IP-Glasma+music model out to fairly peripheral centrality bins. The natural question that arises is how the model does in describing data from p+A collisions that produce similar multiplicities. Fig. 1 shows the multiplicity dependence of and in peripheral Pb+Pb collisions and p+Pb collisions with comparable multiplicity. While the agreement with experimental data from the CMS collaboration Chatrchyan et al. (2013) is fairly good in the Pb+Pb case, for p+Pb collisions is under-predicted by approximately a factor of 4. In p+Pb agrees for the lower multiplicities studied, but has a rather flat multiplicity dependence and underestimates the experimental data at the higher multiplicities.

One reason for the disagreement could be that all initial state correlations that lead to an elliptic anisotropy111When including fluctuations initial state correlations should also contribute to odd harmonics. Dusling and Venugopalan (2012, 2013a, 2013b, 2013c); Gyulassy et al. (2014); Dumitru and Giannini (2014); McLerran and Skokov (2014) are neglected. Additionally, the description of the proton in the IP-Glasma model is over-simplified. The proton shape is approximated by a sphere, and any deviations from that are due to the small scale color charge fluctuations. Since the interaction region in p+Pb collisions is dominated by the shape of the smaller projectile, initial geometries in p+Pb collisions have very small eccentricities. This inevitably leads to small coefficients. If there is indeed a large contribution to coefficients in p+Pb collisions from collective effects, our result indicates that the shape of the proton fluctuates significantly more than assumed in the IP-Glasma model. One could envision a description where the small gluon distributions are still concentrated around large valence quark positions, leading to much larger eccentricities and fluctuations. In this case p+A collisions could be used to determine the shape and fluctuations of gluon distributions in the proton at high energies.

4 p+Au, d+Au, and He+Au collisions at RHIC

To determine whether final state collective effects provide the dominant contribution to the measured azimuthal anisotropy, RHIC is now studying He+Au collisions that on average generate more triangular initial state configurations compared to p+Au or d+Au. If collectivity is the physical explanation for the observed anisotropies, we expect a larger in He+Au collisions compared to p+Au and d+Au collisions at the same energy. To make this expectation more quantitative, we present predictions from the IP-Glasma+music framework.

For deuteron-gold collisions (d+Au) we compute the nucleon distribution in the deuteron using the Hulthen form of its wave function Miller et al. (2007); Bzdak et al. (2013). For He, we use the same nucleon configurations as employed in Nagle et al. (2013). They are obtained from Green’s function Monte Carlo calculations using the AV18 + UIX model interaction Carlson and Schiavilla (1998).

For this comparative study we do not perform a detailed centrality selection, but instead sample the impact parameter between 0 and in all systems. We then compute the initial state distribution of the energy density and flow velocity at time and evolve the system using viscous fluid dynamics with until freeze-out at .

We present typical configurations of the initial energy density distribution in the transverse plane and final results for the transverse momentum dependent azimuthal anisotropy coefficients to in Fig. 2. While we find very small values for through in p+Au collisions, the additional nucleons and their position fluctuations generate larger in d+Au and He+Au collisions. The odd harmonics and are noticeably larger in He+Au collisions compared to d+Au collisions. This qualitative prediction can be compared to future measurements at RHIC.

Figure 2: Typical configurations of the initial energy density distribution for p+Au, d+Au and He+Au collisions (upper panel). Azimuthal anisotropy coefficients in the three collision systems (lower panel).

5 Conclusions

We have demonstrated that experimental results for and in proton-heavy ion collisions at the LHC are not well described by the IP-Glasma+music model. Reasons for this could be the neglected initial state correlations and/or the lack of a detailed description of the fluctuating subnucleonic structure of the proton. Our results for p+A collisions differ significantly from those in Bozek and Broniowski (2013); Werner et al. (2014); Kozlov et al. (2014); Bzdak and Ma (2014), suggesting that the details of the initial shape in small systems are of paramount importance.

We predict an increase of both and in He+Au collisions compared to d+Au collisions, while the even harmonics are comparable in both systems. The detailed substructure of the nucleon is expected to be less important for the initial state geometry in these collisions compared to p+A collisions.

Acknowledgments We thank Jamie Nagle for providing us with the He nucleon configurations. This work was supported under DOE Contract No. DE-AC02-98CH10886. BPS acknowledges support from a DOE Office of Science Early Career Award. This research used resources of the National Energy Research Scientific Computing Center, which is supported by the DOE Office of Science under Contract No. DE-AC02-05CH11231.


  • Gelis et al. (2010) F. Gelis, E. Iancu, J. Jalilian-Marian, R. Venugopalan, The Color Glass Condensate, Ann.Rev.Nucl.Part.Sci. 60 (2010) 463–489.
  • Schenke et al. (2012a) B. Schenke, P. Tribedy, R. Venugopalan, Fluctuating Glasma initial conditions and flow in heavy ion collisions, Phys. Rev. Lett. 108 (2012a) 252301.
  • Schenke et al. (2012b) B. Schenke, P. Tribedy, R. Venugopalan, Event-by-event gluon multiplicity, energy density and eccentricities at RHIC and LHC, Phys. Rev. C86 (2012b) 034908.
  • Kowalski and Teaney (2003) H. Kowalski, D. Teaney, An impact parameter dipole saturation model, Phys. Rev. D68 (2003) 114005.
  • Kowalski et al. (2008) H. Kowalski, T. Lappi, R. Venugopalan, Nuclear enhancement of universal dynamics of high parton densities, Phys. Rev. Lett. 100 (2008) 022303.
  • Krasnitz and Venugopalan (2000) A. Krasnitz, R. Venugopalan, The initial energy density of gluons produced in very high energy nuclear collisions, Phys. Rev. Lett. 84 (2000) 4309–4312.
  • Krasnitz and Venugopalan (2001) A. Krasnitz, R. Venugopalan, The initial gluon multiplicity in heavy ion collisions, Phys. Rev. Lett. 86 (2001) 1717–1720.
  • Lappi (2003) T. Lappi, Production of gluons in the classical field model for heavy ion collisions, Phys. Rev. C67 (2003) 054903.
  • Schenke et al. (2014) B. Schenke, P. Tribedy, R. Venugopalan, Multiplicity distributions in p+p, p+A and A+A collisions from Yang-Mills dynamics, Phys.Rev. C89 (2014) 024901.
  • Chatrchyan et al. (2013) S. Chatrchyan, et al., Multiplicity and transverse momentum dependence of two- and four-particle correlations in pPb and PbPb collisions, Phys.Lett. B724 (2013) 213–240.
  • Berges et al. (2014) J. Berges, K. Boguslavski, S. Schlichting, R. Venugopalan, Turbulent thermalization process in heavy-ion collisions at ultrarelativistic energies, Phys.Rev. D89 (2014) 074011.
  • Epelbaum and Gelis (2013) T. Epelbaum, F. Gelis, Pressure isotropization in high energy heavy ion collisions, Phys.Rev.Lett. 111 (2013) 232301.
  • Schenke et al. (2010) B. Schenke, S. Jeon, C. Gale, (3+1)D hydrodynamic simulation of relativistic heavy-ion collisions, Phys. Rev. C82 (2010) 014903.
  • Schenke et al. (2011a) B. Schenke, S. Jeon, C. Gale, Elliptic and triangular flow in event-by-event (3+1)D viscous hydrodynamics, Phys. Rev. Lett. 106 (2011a) 042301.
  • Schenke et al. (2011b) B. Schenke, S. Jeon, C. Gale, Higher flow harmonics from (3+1)D event-by-event viscous hydrodynamics, Phys. Rev. C85 (2011b) 024901.
  • Schenke and Venugopalan (2014) B. Schenke, R. Venugopalan, Eccentric protons? Sensitivity of flow to system size and shape in p+p, p+Pb and Pb+Pb collisions , arXiv:1405.3605.
  • Dusling and Venugopalan (2012) K. Dusling, R. Venugopalan, Azimuthal collimation of long range rapidity correlations by strong color fields in high multiplicity hadron-hadron collisions, Phys.Rev.Lett. 108 (2012) 262001.
  • Dusling and Venugopalan (2013a) K. Dusling, R. Venugopalan, Evidence for BFKL and saturation dynamics from dihadron spectra at the LHC, Phys.Rev. D87 (5) (2013a) 051502.
  • Dusling and Venugopalan (2013b) K. Dusling, R. Venugopalan, Explanation of systematics of CMS p+Pb high multiplicity di-hadron data at TeV, Phys.Rev. D87 (5) (2013b) 054014.
  • Dusling and Venugopalan (2013c) K. Dusling, R. Venugopalan, Comparison of the Color Glass Condensate to di-hadron correlations in proton-proton and proton-nucleus collisions, Phys.Rev. D87 (2013c) 094034.
  • Gyulassy et al. (2014) M. Gyulassy, P. Levai, I. Vitev, T. Biro, Non-Abelian Bremsstrahlung and Azimuthal Asymmetries in High Energy p+A Reactions , arXiv:1405.7825.
  • Dumitru and Giannini (2014) A. Dumitru, A. V. Giannini, Initial state angular asymmetries in high energy p+A collisions: spontaneous breaking of rotational symmetry by a color electric field and C-odd fluctuations , arXiv:1406.5781.
  • McLerran and Skokov (2014) L. McLerran, V. V. Skokov, The Eccentric Collective BFKL Pomeron , arXiv:1407.2651.
  • Miller et al. (2007) M. L. Miller, K. Reygers, S. J. Sanders, P. Steinberg, Glauber modeling in high energy nuclear collisions, Ann. Rev. Nucl. Part. Sci. 57 (2007) 205–243.
  • Bzdak et al. (2013) A. Bzdak, B. Schenke, P. Tribedy, R. Venugopalan, Initial state geometry and the role of hydrodynamics in proton-proton, proton-nucleus and deuteron-nucleus collisions, Phys.Rev. C87 (6) (2013) 064906.
  • Nagle et al. (2013) J. Nagle, A. Adare, S. Beckman, T. Koblesky, J. O. Koop, et al., Exploiting Intrinsic Triangular Geometry in Relativistic He3+Au Collisions to Disentangle Medium Properties , arXiv:1312.4565.
  • Carlson and Schiavilla (1998) J. Carlson, R. Schiavilla, Structure and dynamics of few nucleon systems, Rev.Mod.Phys. 70 (1998) 743–842.
  • Bozek and Broniowski (2013) P. Bozek, W. Broniowski, Correlations from hydrodynamic flow in p-Pb collisions, Phys.Lett. B718 (2013) 1557–1561, doi:10.1016/j.physletb.2012.12.051.
  • Werner et al. (2014) K. Werner, M. Bleicher, B. Guiot, I. Karpenko, T. Pierog, Evidence for flow in pPb collisions at 5 TeV from v2 mass splitting, Phys.Rev.Lett. 112 (2014) 232301, doi:10.1103/PhysRevLett.112.232301.
  • Kozlov et al. (2014) I. Kozlov, M. Luzum, G. Denicol, S. Jeon, C. Gale, Transverse momentum structure of pair correlations as a signature of collective behavior in small collision systems , arXiv:1405.3976.
  • Bzdak and Ma (2014) A. Bzdak, G.-L. Ma, Elliptic and triangular flow in p+Pb and peripheral Pb+Pb collisions from parton scatterings , arXiv:1406.2804.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description