
The study of multiparticle correlations and collectivity in the hot and dense
matter created in collisions of heavyions as well as those of smaller systems
such as protonheavyion collisions has progressed to the point where detailed
knowledge of the threedimensional structure of the proton is needed to
confront experimental data. We discuss results from a simple proofofprinciple
initial state parton model which reproduces many features of the data on
protonheavyion 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, DrellYan 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 ElectronIon 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 classicalstatistical numerical simulations have established the
"bottomup" thermalization scenario of Baier et al. as the correct weak
coupling effective theory for thermalization in ultrarelativistic heavyion
collisions. We perform a parametric study of photon production in the various
stages of this bottomup framework to ascertain the relative contribution of
the offequilibrium "Glasma" relative to that of a thermalized QuarkGluon
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 nonequilibrium effects should therefore be taken into account
in studies where weak coupling methods are employed to compute photon yields.

We perform first classicalstatistical real time lattice simulations of
topological transitions in the nonequilibrium Glasma of weakly coupled but
highly occupied gauge fields created immediately after the collision of
ultrarelativistic 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 nonequilibrium
transition rate can be understood when expressed in terms of the magnetic
screening length, which we also extract nonperturbatively. A detailed
investigation of autocorrelation functions of the ChernSimons number
($N_{CS}$) reveals nonMarkovian features of the evolution distinct from
previous simulations of nonAbelian plasmas in thermal equilibrium.

We examine the origins of azimuthal correlations observed in high energy
protonnucleus 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 nonGaussian 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
nonAbelian plasmas and $N$component selfinteracting scalar field theories.
Using classicalstatistical methods, we showed that these systems share the
same selfsimilar 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
selfsimilar 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 BoseEinstein 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
nexttoleading order. Remarkably, this effective theory can be cast as a
vertexresummed kinetic theory. The other novel selfsimilar 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 heavyion 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 heavyion
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 ElectronIon 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 multiparton Fock states in the nuclear wavefunctions.
In particular, the saturation scale that characterizes this multiparton
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 overpopulation gives rise to new universality
classes of manybody systems far from equilibrium. We present theoretical
evidence that important aspects of nonAbelian plasmas in the
ultrarelativistic limit admit a dual description in terms of a Bose condensed
scalar field theory.

This White Paper presents the science case of an ElectronIon Collider (EIC),
focused on the structure and interactions of gluondominated 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
followup to the 2007 NSAC Long Range plan. This document is a culmination of a
communitywide effort in nuclear science following a series of workshops on EIC
physics and, in particular, the focused tenweek 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 classicalstatistical 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
classicalstatistical 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 classicalstatistical 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 tenweek program on "Gluons and the quark sea at
highenergies", 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 ElectronIon 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
indepth studies of quantum chromodynamics. This report is organized around
four major themes: i) the spin and flavor structure of the proton, ii)
threedimensional 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 timedependent 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 nearside 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
nonperturbative 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
semihard scale, albeit with nontrivial 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 hadronhadron 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
multiparton 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
BalitskyKovchegov equation, which we solve with running coupling and with
initial conditions constrained by existing data on electronnucleus 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 nearside 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
(KowalskiTeaney) 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 [hepph]) on factorization in
highenergy nucleusnucleus collisions by computing the inclusive multigluon
spectrum to nexttoleading 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 highenergy electronion collider would explore the nonlinear
weaklycoupled regime of QCD, and test the Color Glass Condensate (CGC)
approach to highenergy 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.