We present a fully three-dimensional model providing initial conditions for energy and net-baryon density distributions in heavy ion collisions at arbitrary collision energy. The model includes the ...dynamical deceleration of participating nucleons or valence quarks, depending on the implementation. The duration of the deceleration continues until the string spanned between colliding participants is assumed to thermalize, which is either after a fixed proper time, or a uctuating time depending on sampled final rapidities. Energy is deposited in space-time along the string, which in general will span a range of space-time rapidities and proper times. We study various observables obtained directly from the initial state model, including net-baryon rapidity distributions, 2-particle rapidity correlations, as well as the rapidity decorrelation of the transverse geometry. Their dependence on the model implementation and parameter values is investigated. Here, we also present the implementation of the model with 3+1 dimensional hydrodynamics, which involves the addition of source terms that deposit energy and net-baryon densities produced by the initial state model at proper times greater than the initial time for the hydrodynamic simulation.
We show within the saturation framework that measurements of exclusive vector meson production at high energy provide evidence for strong geometric fluctuations of the proton. In comparison, the ...effect of saturation scale and color charge fluctuations is weak. This knowledge will allow detailed future measurements of the incoherent cross section to tightly constrain the fluctuating geometry of the proton as a function of the parton momentum fraction x.
We present calculations of bulk properties and multiparticle correlations in a large variety of collision systems within a hybrid formalism consisting of IP-Glasma initial conditions, Music viscous ...relativistic hydrodynamics, and UrQMD microscopic hadronic transport. In particular, we study heavy ion collisions at the Large Hadron Collider (LHC), including Pb + Pb, Xe + Xe, and O + O collisions, and Au + Au, U + U, Ru + Ru, Zr + Zr, and O + O collisions at the Relativistic Heavy Ion Collider (RHIC). We further study asymmetric systems, including p + Au, d + Au, 3He + Au, and p + Pb collisions at various energies as well as p+p collisions at 0.5 and 13 TeV. Here, we describe experimental observables in all heavy ion systems well with one fixed set of parameters, validating the energy and system dependence of the framework. As a result, many observables in the smaller systems are also well described, although they test the limits of the model.
We review recent developments in the theoretical description and understanding of multi-particle correlation measurements in collisions of small projectiles (p/d/3He) with heavy nuclei (Au, Pb) as ...well as proton+proton collisions. We focus on whether the physical processes responsible for the observed long range rapidity correlations and their azimuthal structure are the same in small systems as in heavy ion collisions. In the latter they are interpreted as generated by the initial spatial geometry being transformed into momentum correlations by strong final state interactions. However, explicit calculations show that also initial state momentum correlations are present and could contribute to observables in small systems. If strong final state interactions are present in small systems, recent developments show that results are sensitive to the shape of the proton and its fluctuations.
We determine the transverse system size of the initial nonequilibrium Glasma state and of the hydrodynamically evolving fireball as a function of produced charged particles in p+p, p+Pb, and Pb+Pb ...collisions at the Large Hadron Collider. Our results show features similar to those of recent measurements of Hanbury Brown-Twiss (HBT) radii by the ALICE Collaboration. Azimuthal anisotropy coefficients vn generated by combining the early time Glasma dynamics with viscous fluid dynamics in Pb+Pb collisions are in excellent agreement with experimental data for a wide range of centralities. In particular, event-by-event distributions of the vn values agree with the experimental data out to fairly peripheral centrality bins. In striking contrast, our results for p+Pb collisions significantly underestimate the magnitude and do not reproduce the centrality dependence of data for v2 and v3 coefficients. We argue that the measured vn data and HBT radii strongly constrain the shapes of initial parton distributions across system sizes that would be compatible with a flow interpretation in p+Pb collisions. Alternately, additional sources of correlations may be required to describe the systematics of long-range rapidity correlations in p+p and p+Pb collisions.
Multi-particle correlation observables in the Relativistic Heavy Ion Collider small system scan are computed in a framework that contains both initial state momentum anisotropies from the Color Glass ...Condensate effective theory and final state hydrodynamic evolution. The initial state is computed using the IP-Glasma model and coupled to viscous relativistic hydrodynamic simulations, which are followed by microscopic hadronic transport. All parameters of the calculation were previously constrained using experimental data on Au+Au collisions at the same center of mass energy. We find that the qualitative features of the experimental data, such as the system and centrality dependence of the charged hadron momentum anisotropy, can only be reproduced when final state interactions are present. On the other hand, we also demonstrate that the details of the initial state are crucially important for the quantitative description of observables in the studied small systems, as neglecting the initial transverse flow profile or the initial shear stress tensor, which contain information on the momentum anisotropy from the Color Glass Condensate, has dramatic effects on the produced final state anisotropy. We further show that the initial state momentum anisotropy is correlated with the observed elliptic flow in all small systems, with the effect increasing with decreasing multiplicity. We identify the precise measurement of v2 in d+Au and Au+Au collisions at RHIC energy at the same multiplicity as a means to reveal effects of the initial state momentum anisotropy.
Experimental processes that are sensitive to parton Wigner distributions provide a powerful tool to advance our understanding of proton structure. In this work, we compute gluon Wigner and Husimi ...distributions of protons within the color glass condensate framework, which includes a spatially dependent McLerran-Venugopalan initial configuration and the explicit numerical solution of the Jalilian-Marian–Iancu–McLerran–Weigert–Leonidov–Kovner equations. We determine the leading anisotropy of the Wigner and Husimi distributions as a function of the angle between the impact parameter and transverse momentum. We study experimental signatures of these angular correlations at a proposed electron-ion collider by computing coherent diffractive dijet production cross sections in e+p collisions within the same framework. Specifically, we predict the elliptic modulation of the cross section as a function of the relative angle between the nucleon recoil and dijet transverse momentum for a wide kinematical range. We further predict its dependence on the collision energy, which is dominated by the growth of the proton with decreasing x.
We study coherent diffractive dijet production in electron-hadron and electron-nucleus collisions within the dipole picture. We provide semianalytic results for the differential cross section and ...elliptic anisotropy in the angle between dijet transverse momentum and hadron recoil momentum. We demonstrate the direct relation between angular moments of the dipole amplitude in coordinate space and angular moments of the diffractive dijet cross sections. To perform explicit calculations, we employ two different saturation models, extended to include the target geometry. In the limit of large photon virtuality or quark masses, we find fully analytic results that allow direct insight into how the differential cross section and elliptic anisotropy depend on the saturation scale, target geometry, and kinematic variables. We further provide numerical results for more general kinematics in collisions at a future electron-ion collider, and study the effects of approaching the saturated regime on diffractive dijet observables.
This work presents a three-dimensional dynamical initialization model for relativistic heavy-ion collisions, implementing local energy-momentum conservation and baryon charge fluctuations at string ...junctions. Constraining parameters using experimental data from p + p collisions at various collision energies, the model provides a very good description of the charged hadron and net proton rapidity distributions in Au +Au collisions from 7.7 to 200 GeV and Pb + Pb collisions at 8.77 and 17.3 GeV. Here we demonstrate the importance of fluctuations of baryon densities to string junctions for describing net-proton distributions at collision energies of 62.4 and 200 GeV. Including this improved baryon stopping description along with the requirement of strangeness neutrality also yields a good description of identified particle yields as functions of the collision energy above 7.7 GeV. We further study asymmetric p + Al and (p,d, 3He) + Au collisions at the top RHIC energy and p + Pb, Xe + Xe, and Pb + Pb collisions at LHC energies. We identify the produced particle rapidity distributions in asymmetric collision systems as particularly useful for constraining models of the early-time longitudinal dynamics.