The study of multiparticle correlations and collectivity in the hot and dense matter created in collisions of heavy-ions as well as those of smaller systems such as proton--heavy-ion collisions has progressed to the point where detailed knowledge of the three-dimensional structure of the proton is needed to confront experimental data. We discuss results from a simple proof-of-principle initial state parton model which reproduces many features of the data on proton--heavy-ion collisions that are often ascribed to hydrodynamic flow. We outline how this model can be further improved with particular emphasis on missing elements in our understanding of the dynamical spatial and momentum structure of the proton.

Predictions for cold nuclear matter effects on charged hadrons, identified light hadrons, quarkonium and heavy flavor hadrons, Drell-Yan dileptons, jets, photons, gauge bosons and top quarks produced in $p+$Pb collisions at $\sqrt{s_{_{NN}}} = 8.16$ TeV are compiled and, where possible, compared to each other. Predictions of the normalized ratios of $p+$Pb to $p+p$ cross sections are also presented for most of the observables, providing new insights into the expected role of cold nuclear matter effects. In particular, the role of nuclear parton distribution functions on particle production can now be probed over a wider range of phase space than ever before.

We provide an assessment of the energy dependence of key measurements within the scope of the machine parameters for a U.S. based Electron-Ion Collider (EIC) outlined in the EIC White Paper. We first examine the importance of the physics underlying these measurements in the context of the outstanding questions in nuclear science. We then demonstrate, through detailed simulations of the measurements, that the likelihood of transformational scientific insights is greatly enhanced by making the energy range and reach of the EIC as large as practically feasible.

Recent classical-statistical numerical simulations have established the "bottom-up" thermalization scenario of Baier et al. as the correct weak coupling effective theory for thermalization in ultrarelativistic heavy-ion collisions. We perform a parametric study of photon production in the various stages of this bottom-up framework to ascertain the relative contribution of the off-equilibrium "Glasma" relative to that of a thermalized Quark-Gluon Plasma. Taking into account the constraints imposed by the measured charged hadron multiplicities at RHIC and the LHC, we find that Glasma contributions are important especially for large values of the saturation scale at both energies. These non-equilibrium effects should therefore be taken into account in studies where weak coupling methods are employed to compute photon yields.

We perform first classical-statistical real time lattice simulations of topological transitions in the non-equilibrium Glasma of weakly coupled but highly occupied gauge fields created immediately after the collision of ultra-relativistic nuclei. Simplifying our description by employing SU(2) gauge fields, and neglecting their longitudinal expansion, we find that the rate of topological transitions is initially strongly enhanced relative to the thermal sphaleron transition rate and decays with time during the thermalization process. Qualitative features of the time dependence of this non-equilibrium transition rate can be understood when expressed in terms of the magnetic screening length, which we also extract non-perturbatively. A detailed investigation of auto-correlation functions of the Chern-Simons number ($N_{CS}$) reveals non-Markovian features of the evolution distinct from previous simulations of non-Abelian plasmas in thermal equilibrium.

We examine the origins of azimuthal correlations observed in high energy proton-nucleus collisions by considering the simple example of the scattering of uncorrelated partons off color fields in a large nucleus. We demonstrate how the physics of fluctuating color fields in the color glass condensate (CGC) effective theory generates these azimuthal multiparticle correlations and compute the corresponding Fourier coefficients v_n within different CGC approximation schemes. We discuss in detail the qualitative and quantitative differences between the different schemes. We will show how a recently introduced color field domain model that captures key features of the observed azimuthal correlations can be understood in the CGC effective theory as a model of non-Gaussian correlations in the target nucleus.

In Phys. Rev. Lett. 114 (2015) 6, 061601, we reported on a new universality class for longitudinally expanding systems, encompassing strongly correlated non-Abelian plasmas and $N$-component self-interacting scalar field theories. Using classical-statistical methods, we showed that these systems share the same self-similar scaling properties for a wide range of momenta in a limit where particles are weakly coupled but their occupancy is high. Here we significantly expand on our previous work and delineate two further self-similar regimes. One of these occurs in the deep infrared (IR) regime of very high occupancies, where the nonequilibrium dynamics leads to the formation of a Bose-Einstein Condensate. The universal IR scaling exponents and the spectral index characterizing the isotropic IR distributions are described by an effective theory derived from a systematic large-$N$ expansion at next-to-leading order. Remarkably, this effective theory can be cast as a vertex-resummed kinetic theory. The other novel self-similar regime occurs close to the hard physical scale of the theory, and sets in only at later times. We argue that the important role of the infrared dynamics ensures that key features of our results for scalar and gauge theories cannot be reproduced consistently in conventional kinetic theory frameworks.

This document provides a summary of the discussions during the recent joint QCD Town Meeting at Temple University of the status of and future plans for the research program of the relativistic heavy-ion community. A list of compelling questions is formulated, and a number of recommendations outlining the greatest research opportunities and detailing the research priorities of the heavy-ion community, voted on and unanimously approved at the Town Meeting, are presented. They are supported by a broad discussion of the underlying physics and its relation to other subfields. Areas of overlapping interests with the "QCD and Hadron Structure" ("cold QCD") subcommunity, in particular the recommendation for the future construction of an Electron-Ion Collider, are emphasized. The agenda of activities of the "hot QCD" subcommunity at the Town Meeting is attached.

We argue that the proton multiplicities measured in Roman pot detectors at an electron ion collider can be used to determine centrality classes in incoherent diffractive scattering. Incoherent diffraction probes the fluctuations in the interaction strengths of multi-parton Fock states in the nuclear wavefunctions. In particular, the saturation scale that characterizes this multi-parton dynamics is significantly larger in central events relative to minimum bias events. As an application, we study the centrality dependence of incoherent diffractive vector meson production. We identify an observable which is simultaneously very sensitive to centrality triggered parton fluctuations and insensitive to details of the model.

Isolated quantum systems in extreme conditions can exhibit unusually large occupancies per mode. This over-population gives rise to new universality classes of many-body systems far from equilibrium. We present theoretical evidence that important aspects of non-Abelian plasmas in the ultra-relativistic limit admit a dual description in terms of a Bose condensed scalar field theory.

This White Paper presents the science case of an Electron-Ion Collider (EIC), focused on the structure and interactions of gluon-dominated matter, with the intent to articulate it to the broader nuclear science community. It was commissioned by the managements of Brookhaven National Laboratory (BNL) and Thomas Jefferson National Accelerator Facility (JLab) with the objective of presenting a summary of scientific opportunities and goals of the EIC as a follow-up to the 2007 NSAC Long Range plan. This document is a culmination of a community-wide effort in nuclear science following a series of workshops on EIC physics and, in particular, the focused ten-week program on "Gluons and quark sea at high energies" at the Institute for Nuclear Theory in Fall 2010. It contains a brief description of a few golden physics measurements along with accelerator and detector concepts required to achieve them, and it benefited from inputs from the users' communities of BNL and JLab. This White Paper offers the promise to propel the QCD science program in the U.S., established with the CEBAF accelerator at JLab and the RHIC collider at BNL, to the next QCD frontier.

Different thermalization scenarios for systems with large fields have been proposed in the literature based on classical-statistical lattice simulations approximating the underlying quantum dynamics. We investigate the range of validity of these simulations for condensate driven as well as fluctuation dominated initial conditions for the example of a single component scalar field theory. We show that they lead to the same phenomenon of turbulent thermalization for the whole range of (weak) couplings where the classical-statistical approach is valid. In the turbulent regime we establish the existence of a dual cascade characterized by universal scaling exponents and scaling functions. This complements previous investigations where only the direct energy cascade has been studied for the single component theory. A proposed alternative thermalization scenario for stronger couplings is shown to be beyond the range of validity of classical-statistical simulations.

Predictions for charged hadron, identified light hadron, quarkonium, photon, jet and gauge bosons in p+Pb collisions at sqrt s_NN = 5 TeV are compiled and compared. When test run data are available, they are compared to the model predictions.

We report recent results on the role of instabilities in the isotropization and thermalization of a longitudinally expanding system of quantum fields.

This report is based on a ten-week program on "Gluons and the quark sea at high-energies", which took place at the Institute for Nuclear Theory in Seattle in Fall 2010. The principal aim of the program was to develop and sharpen the science case for an Electron-Ion Collider (EIC), a facility that will be able to collide electrons and positrons with polarized protons and with light to heavy nuclei at high energies, offering unprecedented possibilities for in-depth studies of quantum chromodynamics. This report is organized around four major themes: i) the spin and flavor structure of the proton, ii) three-dimensional structure of nucleons and nuclei in momentum and configuration space, iii) QCD matter in nuclei, and iv) Electroweak physics and the search for physics beyond the Standard Model. Beginning with an executive summary, the report contains tables of key measurements, chapter overviews for each of the major scientific themes, and detailed individual contributions on various aspects of the scientific opportunities presented by an EIC.

Quantum fluctuations are believed to play an important role in the thermalization of classical fields in inflationary cosmology but their relevance for isotropization/thermalization of the classical fields produced in heavy ion collisions is not completely understood. We consider a scalar $\phi^4$ toy model coupled to a strong external source, like in the Color Glass Condensate description of the early time dynamics of ultrarelativistic heavy ion collisions. The leading order classical evolution of the scalar fields is significantly modified by the rapid growth of time-dependent quantum fluctuations, necessitating an all order resummation of such "secular" terms. We show that the resummed expressions cause the system to evolve in accordance with ideal hydrodynamics. We comment briefly on the thermalization of our quantum system and the extension of our results to a gauge theory.

The near-side ridge observed in A+A collisions at RHIC has been described as arising from the radial flow of Glasma flux tubes formed at very early times in the collisions. We investigate the viability of this scenario by performing a non-perturbative numerical computation of double inclusive gluon production in the Glasma. Our results support the conjecture that the range of transverse color screening of correlations determining the size of the flux tubes is a semi-hard scale, albeit with non-trivial structure. We discuss our results in the context of ridge correlations in the RHIC heavy ion experiments.

We provide a broad overview of the theoretical status and phenomenological applications of the Color Glass Condensate effective field theory describing universal properties of saturated gluons in hadron wavefunctions that are extracted from deeply inelastic scattering and hadron-hadron collision experiments at high energies.

Long range rapidity correlations in A+A collisions are sensitive to strong color field dynamics at early times after the collision. These can be computed in a factorization formalism \cite{GelisLV5} which expresses the $n$-gluon inclusive spectrum at arbitrary rapidity separations in terms of the multi-parton correlations in the nuclear wavefunctions. This formalism includes all radiative and rescattering contributions, to leading accuracy in $\alpha_s\Delta Y$, where $\Delta Y$ is the rapidity separation between either one of the measured gluons and a projectile, or between the measured gluons themselves. In this paper, we use a mean field approximation for the evolution of the nuclear wavefunctions to obtain a compact result for inclusive two gluon correlations in terms of the unintegrated gluon distributions in the nuclear projectiles. The unintegrated gluon distributions satisfy the Balitsky-Kovchegov equation, which we solve with running coupling and with initial conditions constrained by existing data on electron-nucleus collisions. Our results are valid for arbitrary rapidity separations between measured gluons having transverse momenta $p_\perp,q_\perp\gtrsim \qs$, where $\qs$ is the saturation scale in the nuclear wavefunctions. We compare our results to data on long range rapidity correlations observed in the near-side ridge at RHIC and make predictions for similar long range rapidity correlations at the LHC.

We calculate proton and nuclear diffractive structure functions in the IPsat (Kowalski-Teaney) dipole model. This parametrization has previously been shown to provide good agreement with inclusive F_2 measurements and exclusive vector meson measurements at HERA. We discuss how the impact parameter dependence crucially affects our analysis, in particular for small beta.

We obtain a novel result in QCD for long range rapidity correlations between gluons produced in the collision of saturated high energy hadrons or nuclei. This result, obtained in a high energy factorization framework, provides strong justification for the Glasma flux tube picture of coherent strong color fields. Our formalism can be applied to "near side ridge" events at RHIC and in future studies of long range rapidity correlations at the LHC.

We review recent work on understanding the next to leading order corrections to the classical fields that dominate the initial stages of a heavy ion collision. We have recently shown that the leading ln(1/x) divergences of these corrections to gluon multiplicities can be factorized into the JIMWLK evolution of the color charge density distributions.

We compute diffractive structure functions for both protons and nuclei in the framework of Color Glass Condensate models with impact parameter dependence. These models have previously been shown to provide good agreement with inclusive F_2 measurements and exclusive vector meson measurements at HERA. For nuclei, they provide good (parameter free) agreement with the inclusive F_2 data. We demonstrate good agreement of our computations with HERA measurements on inclusive diffraction. We extend our analysis to nuclei and predict the pattern of enhancement and suppression of the diffractive structures functions that can be measured at an Electron Ion Collider. We discuss how the impact parameter dependence crucially affects our analysis, in particular for large invariant masses at fixed Q^2.

We extend previous results (arXiv:0804.2630 [hep-ph]) on factorization in high-energy nucleus-nucleus collisions by computing the inclusive multigluon spectrum to next-to-leading order. The factorization formula is strictly valid for multigluon emission in a slice of rapidity of width \Delta Y< 1/\alpha_s. Our results show that often neglected disconnected graphs dominate the inclusive multigluon spectrum and are essential to prove factorization for this quantity. These results provide a dynamical framework for the Glasma flux tube picture of the striking "ridge"-like correlation seen in heavy ion collisions.

A future high-energy electron-ion collider would explore the non-linear weakly-coupled regime of QCD, and test the Color Glass Condensate (CGC) approach to high-energy scattering. Hard diffraction in deep inelastic scattering off nuclei will provide many fundamental measurements. In this work, the nuclear diffractive structure function F_{2,A}^D is predicted in the CGC framework, and the features of nuclear enhancement and suppression are discussed.