Source author record

Raju Venugopalan

Raju Venugopalan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

82works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

82 published item(s)

preprint2022arXiv

Dijet impact factor in DIS at next-to-leading order in the Color Glass Condensate

We compute the next-to-leading order impact factor for inclusive dijet production in deeply inelastic electron-nucleus scattering at small $x_{\rm Bj}$. Our computation, performed in the framework of the Color Glass Condensate effective field theory, includes all real and virtual contributions in the gluon shock wave background of all-twist lightlike Wilson line correlators. We demonstrate explicitly that the rapidity evolution of these correlators, to leading logarithmic accuracy, is described by the JIMWLK Hamiltonian. When combined with the next-to-leading order JIMWLK Hamiltonian, our results for the impact factor improve the accuracy of the inclusive dijet cross section to $\mathcal{O}(α_s^2 \ln(x_f/x_{\rm Bj}))$, where $x_f$ is a rapidity factorization scale. These results are an essential ingredient in assessing the discovery potential of inclusive dijets to uncover the physics of gluon saturation at the Electron-Ion Collider.

preprint2022arXiv

Measurement of Bell-type inequalities and quantum entanglement from $Λ$-hyperon spin correlations at high energy colliders

Spin correlations of $Λ$-hyperons embedded in the QCD strings formed in high energy collider experiments provide unique insight into their locality and entanglement features. We show from general considerations that while the Clauser-Horne-Shimony-Holt inequality is less stringent for such states, they provide a benchmark for quantum-to-classical transitions induced by varying i) the associated hadron multiplicity, ii) the spin of nucleons, iii) the separation in rapidity between pairs, and iv) the kinematic regimes accessed. These studies also enable the extraction of quantitative measures of quantum entanglement. We first explore such questions within a simple model of a QCD string composed of singlets of two partial distinguishable fermion flavors and compare analytical results to those obtained on quantum hardware. We further discuss a class of spin Hamiltonians that model the dynamics of $Λ$ spin correlations. Prospects for extracting quantum features of QCD strings from hyperon measurements at current and future colliders are outlined.

preprint2022arXiv

The role of the chiral anomaly in polarized deeply inelastic scattering II: Topological screening and transitions from emergent axion-like dynamics

In [1], we demonstrated that the structure function $g_1(x_B,Q^2)$ measured in polarized deeply inelastic scattering (DIS) is dominated by the triangle anomaly in both the Bjorken limit of large $Q^2$ and the Regge limit of small $x_B$. In the worldline formulation of quantum field theory, the triangle anomaly arises from the imaginary part of the worldline effective action. We show explicitly how a Wess-Zumino-Witten term coupling the topological charge density to a primordial isosinglet ${\bar η}$ arises in this framework. We demonstrate the fundamental role played by this contribution both in topological mass generation of the $η^\prime$ and in the cancellation of the off-forward pole arising from the triangle anomaly in the proton's helicity $Σ(Q^2)$. We recover the striking result by Shore and Veneziano that $Σ\propto \sqrt{χ'(0)}$, where $χ'$ is the slope of the QCD topological susceptibility in the forward limit. We construct an axion-like effective action for $g_1$ at small $x_B$ that describes the interplay between gluon saturation and the topology of the QCD vacuum. In particular, we outline the role of "over-the-barrier" sphaleron-like transitions in spin diffusion at small $x_B$. Such topological transitions can be measured in polarized DIS at a future Electron-Ion Collider.

preprint2021arXiv

Classicalization and unitarization of wee partons in QCD and Gravity: The CGC-Black Hole correspondence

We discuss a remarkable correspondence between the description of Black Holes as highly occupied condensates of $N$ weakly interacting gravitons and that of Color Glass Condensates (CGCs) as highly occupied gluon states. In both cases, the dynamics of "wee partons" in Regge asymptotics is controlled by emergent semi-hard scales that lead to perturbative unitarization and classicalization of $2\rightarrow N$ particle amplitudes at weak coupling. In particular, they attain a maximal entropy permitted by unitarity, bounded by the inverse coupling $α$ of the respective constituents. Strikingly, this entropy is equal to the area measured in units of the Goldstone constant corresponding to the spontaneous breaking of Poincar{é} symmetry by the corresponding graviton or gluon condensate. In gravity, the Goldstone constant is the Planck scale, and gives rise to the Bekenstein-Hawking entropy. Likewise, in the CGC, the corresponding Goldstone scale is determined by the onset of gluon screening. We point to further similarities in Black Hole formation, thermalization and decay, to that of the Glasma matter formed from colliding CGCs in ultrarelativistic nuclear collisions, which decays into a Quark-Gluon Plasma.

preprint2020arXiv

Computing real time correlation functions on a hybrid classical/quantum computer

Quantum devices may overcome limitations of classical computers in studies of nuclear structure functions and parton Wigner distributions of protons and nuclei. In this talk, we discuss a worldline approach to compute nuclear structure functions in the high energy Regge limit of QCD using a hybrid quantum computer, by expressing the fermion determinant in the QCD path integral as a quantum mechanical path integral over $0+1$-dimensional fermionic and bosonic world-lines in background gauge fields. Our simplest example of computing the well-known dipole model result for the structure function $F_2$ in the high energy Regge limit is feasible with NISQ era technology using few qubits and shallow circuits. This example can be scaled up in complexity and extended in scope to compute structure functions, scattering amplitudes and other real-time correlation functions in QCD, relevant for example to describe non-equilibrium transport of quarks and gluons in a Quark-Gluon-Plasma.

preprint2020arXiv

Deeply inelastic scattering structure functions on a hybrid quantum computer

We outline a strategy to compute deeply inelastic scattering structure functions using a hybrid quantum computer. Our approach takes advantage of the representation of the fermion determinant in the QCD path integral as a quantum mechanical path integral over 0+1-dimensional fermionic and bosonic worldlines. The proper time evolution of these worldlines can be determined on a quantum computer. While extremely challenging in general, the problem simplifies in the Regge limit of QCD, where the interaction of the worldlines with gauge fields is strongly localized in proper time and the corresponding quantum circuits can be written down. As a first application, we employ the Color Glass Condensate effective theory to construct the quantum algorithm for a simple dipole model of the $F_2$ structure function. We outline further how this computation scales up in complexity and extends in scope to other real-time correlation functions.

preprint2020arXiv

Inclusive prompt photon-jet correlations as a probe of gluon saturation in electron-nucleus scattering at small $x$

We compute the differential cross-section for inclusive prompt photon$+$quark production in deeply inelastic scattering of electrons off nuclei at small $x$ ($e+A$ DIS) in the framework of the Color Glass Condensate effective field theory. The result is expressed as a convolution of the leading order (in the strong coupling $α_{\mathrm{s}}$) impact factor for the process and universal dipole matrix elements, in the limit of hard photon transverse momentum relative to the nuclear saturation scale $Q_{s,A}(x)$. We perform a numerical study of this process for the kinematics of the Electron-Ion Collider (EIC), exploring in particular the azimuthal angle correlations between the final state photon and quark. We observe a systematic suppression and broadening pattern of the back-to-back peak in the relative azimuthal angle distribution, as the saturation scale is increased by replacing proton targets with gold nuclei. Our results suggest that photon+jet final states in inclusive $e+A$ DIS at high energies are in general a promising channel for exploring gluon saturation that is complementary to inclusive and diffractive dijet production. They also provide a sensitive empirical test of the universality of dipole matrix elements when compared to identical measurements in proton-nucleus collisions. However because photon+jet correlations at small $x$ in EIC kinematics require jet reconstruction at small $k_\perp$, it will be important to study their feasibility relative to photon-hadron correlations.

preprint2020arXiv

The role of the chiral anomaly in polarized deeply inelastic scattering I: Finding the triangle graph inside the box diagram in Bjorken and Regge asymptotics

We revisit the role of the chiral "triangle" anomaly in deeply inelastic scattering (DIS) of electrons off polarized protons employing a powerful worldline formalism. We demonstrate how the triangle anomaly appears at high energies in the DIS box diagram for the polarized proton structure function $g_1(x_B,Q^2)$ in both the Bjorken limit of large $Q^2$ and in the Regge limit of small $x_B$. We show that the operator product expansion is not required to extract the anomaly in either asymptotics though it is sufficient in the Bjorken limit. Likewise, the infrared pole in the anomaly arises in both limits. The leading contribution to $g_1$, in both Bjorken and Regge asymptotics, is therefore given by the expectation value of the topological charge density, generalizing a result previously argued by Jaffe and Manohar to hold for the first moment of $g_1$. In follow-up work, we will show how our results motivate the derivation of a helicity-dependent effective action incorporating the physics of the anomaly at small $x_B$ and shall discuss the QCD evolution of $g_1(x_B,Q^2)$ in this framework.

preprint2019arXiv

Extracting many-body correlators of saturated gluons with precision from inclusive photon+dijet final states in deeply inelastic scattering

We highlight the principal results of a computation in the Color Glass Condensate effective field theory (CGC EFT) of the next-to-leading order (NLO) impact factor for inclusive photon+dijet production at Bjorken $x_{\rm Bj} \ll 1$ in deeply inelastic electron-nucleus (e+A DIS) collisions. When combined with extant results for next-to-leading log $x_{\rm Bj}$ JIMWLK renormalization group (RG) evolution of gauge invariant two-point ("dipole") and four-point ("quadrupole") correlators of light-like Wilson lines, the inclusive photon+dijet e+A DIS cross-section can be determined to $\sim 10$\% accuracy. Our computation simultaneously provides the ingredients to compute fully inclusive DIS, inclusive photon, inclusive dijet and inclusive photon+jet channels to the same accuracy. This makes feasible quantitative extraction of many-body correlators of saturated gluons and precise determination of the saturation scale $Q_{S,A}(x_{\rm Bj})$ at a future Electron-Ion Collider. An interesting feature of our NLO result is the structure of the violation of the soft gluon theorem in the Regge limit. Another is the appearance in gluon emission of time-like non-global logs which also satisfy JIMWLK RG evolution.

preprint2019arXiv

NLO impact factor for inclusive photon$+$dijet production in $e+A$ DIS at small $x$

We compute the next-to-leading order (NLO) impact factor for inclusive photon $+$dijet production in electron-nucleus (e+A) deeply inelastic scattering (DIS) at small $x$. An important ingredient in our computation is the simple structure of ``shock wave" fermion and gluon propagators. This allows one to employ standard momentum space Feynman diagram techniques for higher order computations in the Regge limit of fixed $Q^2\gg Λ_{\rm QCD}^2$ and $x\rightarrow 0$. Our computations in the Color Glass Condensate (CGC) effective field theory include the resummation of all-twist power corrections $Q_s^2/Q^2$, where $Q_s$ is the saturation scale in the nucleus. We discuss the structure of ultraviolet, collinear and soft divergences in the CGC, and extract the leading logs in $x$; the structure of the corresponding rapidity divergences gives a nontrivial first principles derivation of the JIMWLK renormalization group evolution equation for multiparton lightlike Wilson line correlators. Explicit expressions are given for the $x$-independent $O(α_s)$ contributions that constitute the NLO impact factor. These results, combined with extant results on NLO JIMWLK evolution, provide the ingredients to compute the inclusive photon $+$ dijet cross-section at small $x$ to $O(α_s^3 \ln(x))$. First results for the NLO impact factor in inclusive dijet production are recovered in the soft photon limit. A byproduct of our computation is the LO photon+ 3 jet (quark-antiquark-gluon) cross-section.

preprint2016arXiv

Energy dependence of the ridge in high multiplicity proton-proton collisions

We demonstrate that the recent measurement of azimuthally collimated, long range rapidity ("ridge") correlations in $\sqrt{s}=13$ TeV proton-proton (p+p) collisions by the ATLAS collaboration at the LHC are in agreement with expectations from the color glass condensate effective theory of high energy QCD. The observation that the integrated near side yield as a function of multiplicity is independent of collision energy is a natural consequence of the fact that multiparticle production is driven by a single semi-hard saturation scale in the color glass condensate framework. We argue further that the azimuthal structure of these recent ATLAS ridge measurements strongly constrain hydrodynamic interpretations of such correlations in high multiplicity p+p collisions.

preprint2016arXiv

Mass ordering of spectra from fragmentation of saturated gluon states in high multiplicity proton-proton collisions

The mass ordering of mean transverse momentum $\left<p_T\right>$ and of the Fourier harmonic coefficient $v_2 (p_T)$ of azimuthally anisotropic particle distributions in high energy hadron collisions is often interpreted as evidence for the hydrodynamic flow of the matter produced. We investigate an alternative initial state interpretation of this pattern in high multiplicity proton-proton collisions at the LHC. The QCD Yang-Mills equations describing the dynamics of saturated gluons are solved numerically with initial conditions obtained from the Color Glass Condensate based IP-Glasma model. The gluons are subsequently fragmented into various hadron species employing the well established Lund string fragmentation algorithm of the PYTHIA event generator. We find that this ab initio initial state approach reproduces characteristic features of bulk spectra, in particular the particle mass dependence of $\left<p_T\right>$ and $v_2 (p_T)$.

preprint2016arXiv

Universal off-equilibrium scaling of critical cumulants in the QCD phase diagram

Exploiting the universality between the QCD critical point and the three dimensional Ising model, closed form expressions derived (arXiv:1506.00645 ) for non-equilibrium critical cumulants on the crossover side of the critical point reveal that they can differ both in magnitude and sign from equilibrium expectations. We demonstrate here that key elements of the Kibble-Zurek framework of non-equilibrium phase transitions can be employed to describe the dynamics of these critical cumulants. Our results suggest that observables sensitive to critical dynamics in heavy-ion collisions should be expressible as universal scaling functions, thereby providing powerful model independent guidance in searches for the QCD critical point.

preprint2015arXiv

$J/ψ$ production and suppression in high energy proton-nucleus collisions

We apply a Color Glass Condensate+Non-Relativistic QCD (CGC+NRQCD) framework to compute $J/ψ$ production in deuteron-nucleus collisions at RHIC and proton-nucleus collisions at the LHC. Our results match smoothly at high $p_\perp$ to a next-to-leading order perturbative QCD + NRQCD computation. Excellent agreement is obtained for $p_\perp$ spectra at RHIC and LHC for central and forward rapidities, as well as for the normalized ratio $R_{pA}$ of these results to spectra in proton-proton collisions. In particular, we observe that the $R_{pA}$ data is strongly bounded by our computations of the same for each of the individual NRQCD channels; this result provides strong evidence that our description is robust against uncertainties in initial conditions and hadronization mechanisms.

preprint2015arXiv

Azimuthal anisotropies in p+Pb collisions from classical Yang-Mills dynamics

We compute single and double inclusive gluon distributions in classical Yang-Mills simulations of proton-lead collisions and extract the associated transverse momentum dependent Fourier harmonics $v_2(p_T)$ and $v_3(p_T)$. Gluons have a large $v_2$ in the initial state, while odd harmonics such as $v_3$ vanish identically at the initial time $τ=0^{+}$. By the time $τ\lesssim 0.4\,{\rm fm/c}$ final state effects in the classical Yang-Mills evolution generate a non-zero $v_{3}$ and only mildly modify the gluon $v_{2}$. Unlike hydrodynamic flow, these momentum space anisotropies are uncorrelated with the global spatial anisotropy of the collision. A principal ingredient for the generation of $v_2$ and $v_3$ in this framework is the event-by-event breaking of rotational invariance in domains the size of the inverse of the saturation scale $Q_s$. In contrast to our findings in p+Pb collisions Yang-Mills simulations of lead-lead collisions generate much smaller values of $v_{2,3} (p_T)$ and additional collective flow effects are needed to explain experimental data. This is because the locally generated anisotropy due to the breaking of rotational invariance is depleted with the increase in the number of uncorrelated domains.

preprint2015arXiv

Probing short-range nucleon-nucleon interactions with an Electron-Ion Collider

We derive the cross-section for exclusive vector meson production in high energy deeply inelastic scattering off a deuteron target that disintegrates into a proton and a neutron carrying large relative momentum in the final state. This cross-section can be expressed in terms of a novel gluon Transition Generalized Parton Distribution (T-GPD); the hard scale in the final state makes the T-GPD sensitive to the short distance nucleon-nucleon interaction. We perform a toy model computation of this process in a perturbative framework and discuss the time scales that allow the separation of initial and final state dynamics in the T-GPD. We outline the more general computation based on the factorization suggested by the toy computation: in particular, we discuss the relative role of "point-like" and "geometric" Fock configurations that control the parton dynamics of short range nucleon-nucleon scattering. With the aid of exclusive $J/Ψ$ production data at HERA, as well as elastic nucleon-nucleon cross-sections, we estimate rates for exclusive deuteron photo-disintegration at a future Electron-Ion Collider (EIC). Our results, obtained using conservative estimates of EIC integrated luminosities, suggest that center-of-mass energies $s_{NN}\sim 12$ GeV$^2$ of the neutron-proton subsystem can be accessed. We argue that the high energies of the EIC can address outstanding dynamical questions regarding the short-range quark-gluon structure of nuclear forces by providing clean gluon probes of such "knockout" exclusive reactions in light and heavy nuclei.

preprint2015arXiv

Real time evolution of non-Gaussian cumulants in the QCD critical regime

We derive a coupled set of equations that describe the non-equilibrium evolution of cumulants of critical fluctuations for space-time trajectories on the cross-over side of the QCD phase diagram. In particular, novel expressions are obtained for the non-equilibrium evolution of non-Gaussian Skewness and Kurtosis cumulants. Utilizing a simple model of the space-time evolution of a heavy- ion collision, we demonstrate that, depending on the relaxation rate of critical fluctuations, Skewness and Kurtosis can differ significantly in magnitude as well as in sign from equilibrium expectations. Memory effects are important and shown to persist even for trajectories that skirt the edge of the critical regime. We use phenomenologically motivated parameterizations of freeze-out curves, and of the beam energy dependence of the net baryon chemical potential, to explore the implications of our model study for the critical point search in heavy-ion collisions.

preprint2014arXiv

Collective effects in light-heavy ion collisions

We present results for the azimuthal anisotropy of charged hadron distributions in A+A, p+A, d+A, and He3+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 root-s=2.76 TeV are well described by the model, the same quantities in root-s=5.02 TeV 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 He3+Au collisions at root-s=200 GeV. For d+Au and 3He+Au collisions we expect the detailed substructure of the nucleon to become less important.

preprint2014arXiv

Comprehensive description of $J/ψ$ production in proton-proton collisions at collider energies

We employ a small $x$ Color Glass Condensate (CGC)+ Non-Relativistic QCD (NRQCD) formalism to compute $J/ψ$ production at low $p_\perp$ in proton-proton collisions at collider energies. Very good agreement is obtained for total cross-sections, rapidity distributions and low momentum $p_\perp$ distributions. Similar agreement is obtained for $ψ^\prime$ production. We observe an overlap region in $p_\perp$ where our results match smoothly to those obtained in a next-to-leading order (NLO) collinearly factorized NRQCD formalism. The relative contribution of color singlet and color octet contributions can be quantified in the CGC+NRQCD framework, with the former contributing approximately $10\%$ of the total cross-section.

preprint2014arXiv

Eccentric protons? Sensitivity of flow to system size and shape in p+p, p+Pb and Pb+Pb collisions

We determine the transverse system size of the initial non-equilibrium 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 are consistent with recent measurements of Hanbury-Brown-Twiss (HBT) radii by the ALICE collaboration. Azimuthal anisotropy coefficients v_n 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 v_n 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 v_2 and v_3 coefficients. We argue that the measured v_n 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.

preprint2014arXiv

Initial state geometry and fluctuations in Au+Au, Cu+Au and U+U collisions at RHIC

We study within the IP-Glasma and two-component MC-Glauber models the effects of initial state geometry and fluctuations on multiplicities and eccentricities for several collision species at the Relativistic Heavy Ion Collider (RHIC). These include copper-gold (Cu+Au), gold-gold (Au+Au), and uranium-uranium (U+U) collisions. The multiplicity densities per participant pair are very similar in all systems studied. Ellipticities vary strongly between collision systems, most significantly for central collisions, while fluctuation driven odd moments vary little between systems. Event-by-event distributions of eccentricities in mid-central collisions are wider in Cu+Au relative to Au+Au and U+U systems. An anti-correlation between multiplicity and eccentricity is observed in ultra central U+U collisions which is weaker in the IP-Glasma model than the two-component MC-Glauber model. In ultra central Au+Au collisions the two models predict opposite signs for the slope of this correlation. Measurements of elliptic flow as a function of multiplicity in such central events can therefore be used to discriminate between models with qualitatively different particle production mechanisms.

preprint2014arXiv

Quarkonium production in high energy proton-nucleus collisions: CGC meets NRQCD

We study the production of heavy quarkonium states in high energy proton-nucleus collisions. Following earlier work of Blaizot, Fujii, Gelis, and Venugopalan, we systematically include both small $x$ evolution and multiple scattering effects on heavy quark pair production within the Color Glass Condensate (CGC) framework. We obtain for the first time expressions in the Non-Relativistic QCD (NRQCD) factorization formalism for heavy quarkonium differential cross sections as a function of transverse momentum and rapidity. We observe that the production of color singlet heavy quark pairs is sensitive to both "quadrupole" and "dipole" Wilson line correlators, whose energy evolution is described by the Balitsky-JIMWLK equations. In contrast, the color octet channel is sensitive to dipole correlators alone. In a quasi-classical approximation, our results for the color singlet channel reduce to those of Dominguez et. al. [1]. We compare our results to those obtained combining the CGC with the color evaporation model and point to qualitative differences in the two approaches.

preprint2014arXiv

Thermalization of the world's smallest fluids: recent developments

The late Gerry Brown was not shy to tackle complex scientific problems that took time to play out but yielded in the end a deeper understanding of many-body phenomena. In this note, prepared for a memorial volume in his honor, we provide a perspective on a couple of outstanding scientific puzzles that have their origin in our understanding of the thermalization of matter in ultrarelativistic heavy ion collisions, and possibly, in high multiplicity proton-proton and proton-nucleus collisions.

preprint2014arXiv

Turbulent thermalization process in high-energy heavy-ion collisions

We discuss the onset of the thermalization process in high-energy heavy-ion collisions from a weak coupling perspective, using classical-statistical real-time lattice simulations as a first principles tool to study the pre-equilibrium dynamics. Most remarkably, we find that the thermalization process is governed by a universal attractor, where the space-time evolution of the plasma becomes independent of the initial conditions and exhibits the self-similar dynamics characteristic of wave turbulence. We discuss the consequences of our weak coupling results for the thermalization process in heavy-ion experiments and briefly comment on the use of weak coupling techniques at larger values of the coupling.

preprint2013arXiv

Comparison of the Color Glass Condensate to di-hadron correlations in proton-proton and proton-nucleus collisions

We perform a detailed comparison of long range rapidity correlations in the Color Glass Condensate framework to high multiplicity di-hadron data in proton-proton and proton-lead collisions from the CMS, ALICE and ATLAS experiments at the LHC. The overall good agreement thus far of the non-trivial systematics of theory with data is strongly suggestive of gluon saturation and the presence of subtle quantum interference effects between rapidity separated gluons. In particular, the yield of pairs collimated in their relative azimuthal angle $Δϕ\sim 0$, is sensitive to the shape of unintegrated gluon distributions in the hadrons that are renormalization group evolved in rapidity from the beam rapidities to those of the measured hadrons. We present estimates for the collimated di-hadron yield expected in central deuteron-gold collisions at RHIC.

preprint2013arXiv

Initial state geometry and the role of hydrodynamics in proton-proton, proton-nucleus and deuteron-nucleus collisions

We apply the successful Monte Carlo Glauber and IP-Glasma initial state models of heavy ion collisions to the much smaller size systems produced in proton-proton, proton-nucleus and deuteron- nucleus collisions. We observe a significantly greater sensitivity of the initial state geometry to details of multi-particle production in these models compared to nucleus-nucleus collisions. In particular, we find that the size of the system produced in p+A collisions is very similar to the one produced in p+p collisions, and predict comparable Hanbury-Brown-Twiss radii in the absence of flow in both systems. Differences in the eccentricities computed in the models are large, while differences amongst the generated flow coefficients v_2 and v_3 are smaller. For a large number of participants in proton-lead collisions, the v_2 generated in the IP-Glasma model is comparable to the value obtained in proton-proton collisions. Viscous corrections to flow are large over characteristic lifetimes in the smaller size systems. In contrast, viscous contributions are significantly diminished over the longer space-time evolution of a heavy ion collision.

preprint2013arXiv

IP-Sat: Impact-Parameter dependent Saturation model; revised

In this talk, we present a global analysis of available small-x data on inclusive DIS and exclusive diffractive processes, including the latest data from the combined HERA analysis on reduced cross sections within the Impact-Parameter dependent Saturation (IP-Sat) Model. The impact-parameter dependence of dipole amplitude is crucial in order to have a unified description of both inclusive and exclusive diffractive processes. With the parameters of model fixed via a fit to the high-precision reduced cross-section, we compare model predictions to data for the structure functions, the longitudinal structure function, the charm structure function, exclusive vector mesons production and Deeply Virtual Compton Scattering (DVCS). Excellent agreement is obtained for the processes considered at small x in a wide range of Q^2.

preprint2013arXiv

Multiplicity distributions in p+p, p+A and A+A collisions from Yang-Mills dynamics

We compute transverse momentum and momentum integrated multiplicity distributions consistently in the IP-Glasma model for proton-proton and proton-lead collisions at the LHC, in deuteron-gold collisions at RHIC, and in heavy ion collisions at both RHIC and LHC energies. Several sources of sub-nucleon scale contributions to the multiplicity distributions are identified. Our results, which are constrained by inclusive and diffractive deeply inelastic scattering data from HERA, are compared to measured distributions for a range of collision energies. These results are an essential first step in quantifying the relative role of initial and final state effects on multiparticle correlations in light and heavy ion collisions.

preprint2013arXiv

Universal attractor in a highly occupied non-Abelian plasma

We study the thermalization process in highly occupied non-Abelian plasmas at weak coupling. The non-equilibrium dynamics of such systems is classical in nature and can be simulated with real-time lattice gauge theory techniques. We provide a detailed discussion of this framework and elaborate on the results reported in~\cite{Berges:2013eia} along with novel findings. We demonstrate the emergence of universal attractor solutions, which govern the non-equilibrium evolution on large time scales both for non-expanding and expanding non-Abelian plasmas. The turbulent attractor for a non-expanding plasma drives the system close to thermal equilibrium on a time scale $t\sim Q^{-1} α_s^{-7/4}$. The attractor solution for an expanding non-Abelian plasma leads to a strongly interacting albeit highly anisotropic system at the transition to the low-occupancy or quantum regime. This evolution in the classical regime is, within the uncertainties of our simulations, consistent with the ``bottom up'' thermalization scenario~\cite{Baier:2000sb}. While the focus of this paper is to understand the non-equilibrium dynamics in weak coupling asymptotics, we also discuss the relevance of our results for larger couplings in the early time dynamics of heavy ion collision experiments.

preprint2012arXiv

Analysis of combined HERA data in the Impact-Parameter dependent Saturation model

The Impact-Parameter dependent Saturation Model (IP-Sat) is a simple dipole model that incorporates key features of the physics of gluon saturation and matches smoothly to the perturbative QCD dipole expression at large Q^2 for a given x. It was previously shown that the model gives a good description of HERA data suggesting evidence for gluon saturation effects at small x. The model has also been applied to proton-proton and proton-nucleus collisions and provides the basis for the IP-Glasma model of initial conditions in heavy ion collisions. Here we present a reanalysis of available data in electron-proton collisions at small Bjorken-x, including the recently released combined data from the ZEUS and H1 collaborations. We first confront the model to the high precision combined data for the reduced cross-section and obtain its parameters. With these parameters fixed, we compare model results to data for the structure function F_2, the longitudinal structure function F_L, the charm structure function F_2^{c\bar{c}}, exclusive vector meson (J/ψ, ϕand ρ) production and Deeply Virtual Compton Scattering (DVCS). Excellent agreement is obtained for the processes considered at small x in a wide range of Q^2. Our results strongly hint at universality of the IP-Sat dipole amplitude and the extracted impact-parameter distribution of the proton. They also provide a benchmark for further refinements in studies of QCD saturation at colliders.

preprint2012arXiv

Azimuthal collimation of long range rapidity correlations by strong color fields in high multiplicity hadron-hadron collisions

The azimuthal collimation of di-hadrons with large rapidity separations in high multiplicity p+p collisions at the LHC is described in the Color Glass Condensate (CGC) effective theory [1] by N_c^2 suppressed multi-ladder QCD diagrams that are enhanced α_S^(-8) due to gluon saturation in hadron wavefunctions. We show that quantitative computations in the CGC framework are in good agreement with data from the CMS experiment on per trigger di-hadron yields and predict further systematics of these yields with varying trigger pT and charged hadron multiplicity. Radial flow generated by re-scattering is strongly limited by the structure of the p+p di-hadron correlations. In contrast, radial flow explains the systematics of identical measurements in heavy ion collisions.

preprint2012arXiv

Event-by-event anisotropic flow in heavy-ion collisions from combined Yang-Mills and viscous fluid dynamics

Anisotropic flow coefficients v_1-v_5 in heavy ion collisions are computed by combining a classical Yang-Mills description of the early time glasma flow with the subsequent relativistic viscous hydrodynamic evolution of matter through the quark-gluon plasma and hadron gas phases. The glasma dynamics, as realized in the IP-Glasma model, takes into account event-by-event geometric fluctuations in nucleon positions and intrinsic sub-nucleon scale color charge fluctuations; the pre-equilibrium flow of matter is then matched to the MUSIC algorithm describing viscous hydrodynamic flow and particle production at freeze-out. The IP-Glasma+MUSIC model describes well both transverse momentum dependent and integrated v_n data measured at the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC). The model also reproduces the event-by-event distributions of v_2, v_3 and v_4 measured by the ATLAS collaboration. The implications of our results for better understanding of the dynamics of the glasma as well as for the extraction of transport properties of the quark-gluon plasma are outlined.

preprint2012arXiv

Event-by-event gluon multiplicity, energy density and eccentricities at RHIC and LHC

The event-by-event multiplicity distribution, the energy densities and energy density weighted eccentricity moments epsilon_n (up to n=6) at early times in heavy-ion collisions at both RHIC (root-s=200 GeV) and LHC (root-s=2.76 TeV) are computed in the IP-Glasma model. This framework combines the impact parameter dependent saturation model (IP-Sat) for nucleon parton distributions (constrained by HERA deeply inelastic scattering data) with an event-by-event classical Yang-Mills description of early-time gluon fields in heavy-ion collisions. The model produces multiplicity distributions that are convolutions of negative binomial distributions without further assumptions or parameters. The eccentricity moments are compared to the MC-KLN model; a noteworthy feature is that fluctuation dominated odd moments are consistently larger than in the MC-KLN model.

preprint2012arXiv

Explanation of systematics of CMS p+Pb high multiplicity di-hadron data at $\sqrt{s}_{\rm NN} = 5.02$ TeV

In a recent article (arXiv:1210.3890), we showed that high multiplicity di-hadron proton- proton (p+p) data from the CMS experiment are in excellent agreement with computations in the Color Glass Condensate (CGC) Effective Field Theory (EFT). This agreement of the theory with several hundred data points provides a non-trivial description of both nearside ("ridge") and away-side azimuthal collimations of long range rapidity correlations in p+p collisions. Our prediction in arXiv:1210.3890 for proton-lead (p+Pb) collisions is consistent with results from the recent CMS p+Pb run at $\sqrt{s}_{\rm NN} = 5.02$ TeV for the largest track multiplicity $N_{\rm track}\sim 40$ we considered. The CMS p+Pb data shows the following striking features: i) a strong dependence of the ridge yield on $N_{\rm track}$, with a significantly larger signal than in p+p for the same $N_{\rm track}$, ii) a stronger $p_T$ dependence than in p+p for large $N_{\rm track}$, and iii) a nearside collimation for large $N_{\rm track}$ comparable to the awayside for the lower $p_T = p_{T}^{\rm trig.}=p_{T}^{\rm assoc.}$ di-hadron windows. We show here that these systematic features of the CMS p+Pb di-hadron data are all described by the CGC (with parameters fixed by the p+p data) when we extend our prediction in arXiv:1210.3890 to rarer high multiplicity events. We also predict the azimuthally collimated yield for yet unpublished windows in the $p_{T}^{\rm trig.}$ and $p_{T}^{\rm assoc.}$ matrix.

preprint2012arXiv

Fluctuating Glasma initial conditions and flow in heavy ion collisions

We compute initial conditions in heavy-ion collisions within the Color Glass Condensate (CGC) framework by combining the impact parameter dependent saturation model (IP-Sat) with the classical Yang-Mills description of initial Glasma fields. In addition to fluctuations of nucleon positions, this IP-Glasma description includes quantum fluctuations of color charges on the length-scale determined by the inverse nuclear saturation scale Q_s. The model naturally produces initial energy fluctuations that are described by a negative binomial distribution. The ratio of triangularity to eccentricity is close to that in a model tuned to reproduce experimental flow data. We compare transverse momentum spectra and v_(2,3,4)(p_T) of pions from different models of initial conditions using relativistic viscous hydrodynamic evolution.

preprint2012arXiv

Initial state fluctuations and higher harmonic flow in heavy-ion collisions

A framework combining Yang-Mills dynamics of the pre-equilibrium glasma with relativistic viscous hydrodynamic evolution of the quark-gluon plasma and hadron gas phases is presented. Event-by-event fluctuations of nucleon positions and color charges are taken into account, leading to negative binomial fluctuations of gluon multiplicities. Experimental anisotropic flow coefficients v2-v5 of charged hadron distributions in heavy-ion collisions at the Large Hadron Collider are well described. Furthermore, event-by-event distributions of v2, v3 and v4 measured by the ATLAS collaboration are reproduced.

preprint2012arXiv

Instability induced pressure isotropization in a longitudinally expanding system

In two previous works [arXiv:1009.4363,arXiv:1107.0668], we studied the time evolution of a system of real scalar fields with quartic coupling which shares important features with the Color Glass Condensate description of heavy ion collisions. Our primary objective was to understand how such a system, when initialized with a non-perturbatively large classical field configuration, reaches thermal equilibrium. An essential goal of these works was to highlight the role played by the quantum fluctuations. However, these studies considered only a system confined within a box of fixed volume. In the present paper, we extend this work to a system that expands in the longitudinal direction thereby more closely mimicking a heavy ion collision. We conclude that the microscopic processes that drive the system towards equilibrium are able to keep up with the expansion of the system; the pressure tensor becomes isotropic despite the anisotropic expansion.

preprint2012arXiv

QCD saturation at the LHC: comparisons of models to p+p and A+A data and predictions for p+Pb collisions

In a previous paper (arXiv:1011.1895), we showed that saturation models, constrained by e+p HERA data on inclusive and diffractive cross-sections, are in good agreement with p+p data at LHC in the soft sector. Particularly impressive was the agreement of saturation models with the multiplicity distribution as a function of $n_{\rm ch.}$. In this paper, we extend these studies further and consider the agreement of these models with data on bulk distributions in A+A collisions. We compare our results to data on central and forward particle production in d+Au collisions at RHIC and make predictions for inclusive distributions in p+Pb collisions at the LHC.

preprint2012arXiv

The dynamics of strongly correlated gluons at high energies

We describe some of the recent progress in our understanding of the dynamics of strongly correlated gluons at high parton densities. Computations in the Color Glass Condensate effective field theory provide a good description of inclusive and semi-inclusive final states in DIS, p+p and p+A collisions at small x. In nucleus-nucleus collisions, they provide an ab initio description of entropy generation, decoherence, isotropization and the onset of hydrodynamic flow. The successful description of a wide range of phenomena from RHIC to LHC is outlined and possibilities for more stringent tests noted.

preprint2012arXiv

Thermalization and Bose-Einstein Condensation in Overpopulated Glasma

We report recent progress on understanding the thermalization of the quark-gluon plasma during the early stage in a heavy ion collision. The initially high overpopulation in the far-from-equilibrium gluonic matter ("Glasma") is shown to play a crucial role. The strongly interacting nature (and thus fast evolution) naturally arises as an {\em emergent property} of this pre-equilibrium matter where the intrinsic coupling is weak but the highly occupied gluon states coherently amplify the scattering. A possible transient Bose-Einstein Condensate is argued to form dynamically on a rather general ground. We develop a kinetic approach for describing its evolution toward thermalization, and based on that we find approximate scaling solutions as well as numerically study the onset of condensation.

preprint2011arXiv

Bose--Einstein Condensation and Thermalization of the Quark Gluon Plasma

In ultra-relativistic heavy ion collisions, the matter formed shortly after the collision is a dense, out of equilibrium, system of gluons characterized by a semi-hard momentum scale $Q_{\rm s}$. Simple power counting arguments indicate that this system is over-occupied: the gluon occupation number is parametrically large when compared to a system in thermal equilibrium with the same energy density. On short time scales, soft elastic scatterings tend to drive the system towards the formation of a Bose--Einstein condensate that contains a large fraction of the gluons while contributing little to the energy density. The lifetime and existence of this condensate depends on whether inelastic processes, that occur on the same time scale as the elastic ones, preferably increase or decrease the number of gluons. During this overpopulated stage, and all the way to thermalization, the system behaves as a strongly interacting fluid, even though the elementary coupling constant is small. We argue that while complete isotropization may never be reached, the system may yet evolve for a long time with a fixed anisotropy between average longitudinal and transverse momenta.

preprint2011arXiv

Ekpyrosis and inflationary dynamics in heavy ion collisions: the role of quantum fluctuations

We summarize recent significant progress in the development of a first-principles formalism to describe the formation and evolution of matter in very high energy heavy ion collisions. The key role of quantum fluctuations both before and after a collision is emphasized. Systematic computations are now feasible to address early time dynamics essential to quantifying properties of strongly interacting quark-gluon matter.

preprint2011arXiv

Inclusive hadron distributions in p+p collisions from saturation models of HERA DIS data

Dipole models based on various saturation scenarios provide reasonable fits to small-x DIS inclusive, diffractive and exclusive data from HERA. Proton un-integrated gluon distributions extracted from such fits are employed in a $k_\bot$-factorization framework to calculate inclusive gluon distributions at various energies. The n-particle multiplicity distribution predicted in the Glasma flux tube approach shows good agreement with data over a wide range of energies. Hadron inclusive transverse momentum distributions expressed in terms of the saturation scale demonstrate universal behavior over a wider kinematic range systematically with increasing center of mass energies.

preprint2011arXiv

Renormalization group evolution of multi-gluon correlators in high energy QCD

Many-body QCD in leading high energy Regge asymptotics is described by the Balitsky-JIMWLK hierarchy of renormalization group equations for the x evolution of multi-point Wilson line correlators. These correlators are universal and ubiquitous in final states in deeply inelastic scattering and hadronic collisions. For instance, recently measured di-hadron correlations at forward rapidity in deuteron-gold collisions at the Relativistic Heavy Ion Collider (RHIC) are sensitive to four and six point correlators of Wilson lines in the small x color fields of the dense nuclear target. We evaluate these correlators numerically by solving the functional Langevin equation that describes the Balitsky-JIMWLK hierarchy. We compare the results to mean-field Gaussian and large N_c approximations used in previous phenomenological studies. We comment on the implications of our results for quantitative studies of multi-gluon final states in high energy QCD.

preprint2011arXiv

Saturation models of HERA DIS data and inclusive hadron distributions in p+p collisions at the LHC

Unintegrated gluon distributions sensitive to the transverse spatial distribution of gluons in the proton are extracted from data on exclusive and diffractive final states at HERA in the dipole approach. These unintegrated gluon distributions can be used to compute inclusive hadron production in p+p collisions at the LHC. In this paper, we consider a number of saturation models with differering dynamical assumptions that give good fits to the available HERA data. We apply these models to study the rapidity and transverse momentum dependence of the LHC data up to $\sqrt{s} = 7$ TeV. We examine the sensitivity of these results to parameters that are not constrained by the HERA data and comment on similarities and differences with previous work. We compute the n-particle inclusive multiplicity distribution and show that the LHC p+p results are in agreement with predictions for multi-particle production in the Color Glass Condensate approach. This result has significant ramifications for the interpretation of multi-particle correlations in high multiplicity events at the LHC.

preprint2011arXiv

The initial spectrum of fluctuations in the little bang

High parton densities in ultra-relativistic nuclear collisions suggest a description of these collisions wherein the high energy nuclear wavefunctions and the initial stages of the nuclear collision are dominated by classical fields. This underlying paradigm can be significantly improved by including quantum fluctuations around the classical background fields. One class of these contributes to the energy evolution of multi-parton correlators in the nuclear wavefunctions. Another dominant class of unstable quantum fluctuations grow rapidly with proper time $τ$ after the collision. These secular terms appear at each loop order; the leading contributions can be resummed to all loop orders to obtain expressions for final state observables. The all-order result can be expressed in terms of the spectrum of fluctuations on the initial proper time surface. We compute, in $A^τ=0$ gauge, the essential elements in this fluctuation spectrum--the small quantum fluctuation modes in the classical background field. With our derivation in QCD, we have all the ingredients to compute inclusive quantities in heavy ion collisions at early times including i) all--order leading logs in Bjorken $x_{1,2}$ of the two nuclei, ii) all strong multiple scattering contributions, and iii) all-order leading secular terms. In the simpler analogous formalism for a scalar $ϕ^4$ theory, numerical analysis of the behavior of the energy-momentum tensor is strongly suggestive of early hydrodynamic flow in the system. In QCD, in addition to studying the possible early onset of hydrodynamic behavior, additional important applications of our results include a) the computation of sphaleron transitions off-equilibrium, and b) "jet quenching", or medium modification of parton spectra, in strong color fields at early times.

preprint2010arXiv

From many body wee parton dynamics to perfect fluid: a standard model for heavy ion collisions

We discuss a standard model of heavy ion collisions that has emerged both from the experimental results of the RHIC program and associated theoretical developments. We comment briefly on the impact of early results of the LHC program on this picture. We consider how this standard model of heavy ion collisions could be solidified or falsified in future experiments at RHIC, the LHC and a future Electron-Ion Collider.

preprint2010arXiv

The ridge in proton-proton collisions at the LHC

We show that the key features of the CMS result on the ridge correlation seen for high multiplicity events in sqrt(s)=7TeV proton-proton collisions at the LHC can be understood in the Color Glass Condensate framework of high energy QCD. The same formalism underlies the explanation of the ridge events seen in A+A collisions at RHIC, albeit it is likely that flow effects may enhance the magnitude of the signal in the latter.

preprint2008arXiv

Glasma flux tubes and the near side ridge phenomenon at RHIC

We investigate the consequences of long range rapidity correlations in the Glasma. Particles produced locally in the transverse plane are correlated by approximately boost invariant flux tubes of longitudinal color electric and magnetic fields that are formed when two sheets of Colored Glass Condensate pass through one another, each acquiring a modified color charge density in the collision. We argue that such long range rapidity correlations persist during the evolution of the Quark Gluon Plasma formed later in the collision. When combined with transverse flow, these correlations reproduce many of the features of the recently observed ridge events in heavy ion collisions at RHIC.

preprint2007arXiv

High energy scattering in Quantum Chromodynamics

In this series of three lectures, we discuss several aspects of high energy scattering among hadrons in Quantum Chromodynamics. The first lecture is devoted to a description of the parton model, Bjorken scaling and the scaling violations due to the evolution of parton distributions with the transverse resolution scale. The second lecture describes parton evolution at small momentum fraction x, the phenomenon of gluon saturation and the Color Glass Condensate (CGC). In the third lecture, we present the application of the CGC to the study of high energy hadronic collisions, with emphasis on nucleus-nucleus collisions. In particular, we provide the outline of a proof of high energy factorization for inclusive gluon production.

preprint2006arXiv

Collective non-Abelian instabilities in a melting Color Glass Condensate

We present first results for 3+1-D simulations of SU(2) Yang-Mills equations for matter expanding into the vacuum after a heavy ion collision. Violations of boost invariance cause a Weibel instability leading soft modes to grow with proper time $τ$ as $\exp(Γ\sqrt{g^2μτ})$, where $g^2μ$ is a scale arising from the saturation of gluons in the nuclear wavefunction. The scale for the growth rate $Γ$ is set by a plasmon mass, defined as $ω_{\rm pl}= κ_0 \sqrt{\frac{g^2μ}τ}$, generated dynamically in the collision. We compare the numerical ratio $Γ/κ_0$ to the corresponding value predicted by the Hard Thermal Loop formalism for anisotropic plasmas.

preprint2006arXiv

The Unstable Glasma

We discuss results from 3+1-D numerical simulations of SU(2) Yang--Mills equations for an unstable Glasma expanding into the vacuum after a high energy heavy ion collision. We expand on our earlier work on a non-Abelian Weibel instability in such a system and study the behavior of the instability in greater detail on significantly larger lattices than previously. We establish the time scale for the onset of the instability and demonstrate that the growth rate is robust as one approaches the continuum limit. For large violations of boost invariance, non-Abelian effects cause the growth of soft modes to saturate. At late times, we observe significant creation of longitudinal pressure and a systematic trend towards isotropy. These time scales however are significantly larger than those required for early thermalization in heavy ion collisions. We discuss additional effects in the produced Glasma that may speed up thermalization.

preprint2005arXiv

Quantitative study of the violation of kt-factorization in hadroproduction of quarks at collider energies

We demonstrate the violation of kT-factorization for quark production in high energy hadronic collisions. This violation is quantified in the Color Glass Condensate framework and studied as a function of the quark mass, the quark transverse momentum, and the saturation scale Qsat, which is a measure of large parton densities. At x values where parton densities are large but leading twist shadowing effects are still small, violations of kT-factorization can be significant - especially for lighter quarks. At very small x, where leading twist shadowing is large, we show that violations of kT-factorization are relatively weaker.

preprint2003arXiv

Classical Gluodynamics of High Energy Nuclear Collisions: an Erratum and an Update

We comment on the relation of our previous work on the classical gluodynamics of high energy nuclear collisions to recent work by Lappi (hep-ph/0303076). While our results for the non-perturbative number liberation coefficient agree, those for the energy disagree by a factor of 2. This discrepancy can be traced to an overall normalization error in our non-perturbative formula for the energy. When corrected for, all previous results are in excellent agreement with those of Lappi. The implications of the results of these two independent computations for RHIC phenomenology are noted.

preprint2003arXiv

Finite-Size Effects on Nucleation in a First-Order Phase Transition

We discuss finite-size effects on homogeneous nucleation in first-order phase transitions. We study their implications for cosmological phase transitions and to the hadronization of a quark-gluon plasma generated in high-energy heavy ion collisions. Very general arguments allow us to show that the finite size of the early universe has virtually no relevance in the process of nucleation and in the growth of cosmological bubbles during the primordial quark-hadron and the electroweak phase transitions. In the case of high-energy heavy ion collisions, finite-size effects play an important role in the late-stage growth of hadronic bubbles.

preprint2003arXiv

The Color Glass Condensate and High Energy Scattering in QCD

At very high energies or small values of Bjorken x, the density of partons, per unit transverse area, in hadronic wavefunctions becomes very large leading to a saturation of partonic distributions. When the scale corresponding to the density per unit transverse area, the saturation scale Q_s, becomes large (Q_s\gg Λ_{QCD}), the coupling constant becomes weak (α_S(Q_s)\ll 1) which suggests that the high energy limit of QCD may be studied using weak coupling techniques. This simple idea can be formalized in an effective theory, the Color Glass Condensate (CGC), which describes the behavior of the small x components of the hadronic wavefunction in QCD. The Green functions of the theory satisfy Wilsonian renormalization group equations which reduce to the standard linear QCD evolution equations in the limit of low parton densities. The effective theory has a rich structure that has been explored using analytical and numerical techniques. The CGC can be applied to study a wide range of high energy scattering experiments from Deep Inelastic Scattering at HERA and the proposed Electron Ion Collider (EIC) to proton/deuterium-nucleus and nucleus-nucleus experiments at the RHIC and LHC colliders.

preprint2003arXiv

The Cronin Effect, Quantum Evolution and the Color Glass Condensate

We show that the numerical solution of the classical SU(3) Yang-Mills equations of motion in the McLerran-Venugopalan model for gluon production in central heavy ion collisions leads to a suppresion at low $p_t$ and an enhancement at the intermediate $p_t$ region as compared to peripheral heavy ion and pp collisions at the same energy. Our results are compared to previous, Color Glass Condensate inspired calculations of gluon production in heavy ion collisions. We revisit the predictions of the Color Glass Condensate model for $pA$ ($dA$) collisions in Leading Order and show that quantum evolution--in particular the phenomenon of geometric scaling and change of anomalous dimensions--preserves the Cronin enhancement of $pA$ cross section (when normalized to the leading twist term) in the Leading Order approximation even though the $p_t$ spectrum can change. We comment on the case when gluon radiation is included.

preprint2002arXiv

Elliptic flow of colored glass in high energy heavy ion collisions

We compute the elliptic flow generated by classical gluon fields in a high energy nuclear collision. The classical gluon fields are described by a typical momentum scale, the saturation scale Lambda_s, which is, for RHIC energies, of the order of 1-2 GeV. A significant elliptic flow is generated only over time scales on the order of the system size R. The flow is dominated by soft modes $p_{T} \sim Lambda_s / 4$ which linearize at very late times $tau \sim R \gg 1/Lambda_s$. We discuss the implications of our result for the theoretical interpretation of the RHIC data.

preprint2002arXiv

Generalized scaling of the transverse mass spectrum at the Relativistic Heavy-Ion Collider

We argue that the transverse mass spectra of identified hadrons as measured in gold-gold collisions at BNL's Relativistic Heavy-Ion Collider (RHIC) follows a generalized scaling law. Such a scaling behavior is motivated by the idea of a Color Glass Condensate, or more generally, saturation of the gluon density. In particular, we describe the shapes of transverse mass spectra as a function of centrality. This scaling of the transverse mass spectrum is shown to be consistent with previously observed scaling of multiplicity with centrality.

preprint2001arXiv

Small x physics and the initial conditions in heavy ion collisions

At very high energies, the high parton densities (characterized by a semi-hard saturation scale Λ_s) ensure that parton distributions can be described by a classical effective field theory with remarkable properties analogous to those of spin glass systems. This color glass condensate (CGC) of gluons also provides the initial conditions for multi-particle production in high energy nuclear collisions. In this talk, we briefly summarize recent theoretical and phenomenological progress in the CGC approach to small x physics. In particular, we discuss recent numerical work on the real time gluodynamics of partons after a nuclear collision. The implications of this work for the theoretical study of thermalization in nuclear collisions and on the phenomenological interpretation of results of the recent RHIC experiments are also discussed.

preprint2000arXiv

From colored glass condensate to gluon plasma: equilibration in high energy heavy ion collisions

The initial distribution of gluons at the very early times after a high energy heavy ion collision is described by the bulk scale $Q_s$ of gluon saturation in the nuclear wavefunction. The subsequent evolution of the system towards kinetic equilibrium is described by a non-linear Landau equation for the single particle distributions \cite{Mueller1,Mueller2}. In this paper, we solve this equation numerically for the idealized initial conditions proposed by Mueller, and study the evolution of the system to equilibrium. We discuss the sensitivity of our results on the dynamical screening of collinear divergences. In a particular model of dynamical screening, the convergence to the hydrodynamic limit is seen to be rapid relative to hydrodynamic time scales. The equilibration time, the initial temperature, and the chemical potential are shown to have a strong functional dependence on the initial gluon saturation scale $Q_s$.

preprint2000arXiv

Wong's equations and the small x effective action in QCD

We propose a new form for the small x effective action in QCD. This form of the effective action is motivated by Wong's equations for classical, colored particles in non-Abelian background fields. We show that the BFKL equation, which sums leading logarithms in x, is efficiently reproduced with this form of the action. We argue that this form of the action may be particularly useful in computing next-to- leading-order results in QCD at small x.

preprint1999arXiv

Classical methods in DIS and nuclear scattering at small x

In hadrons and nuclei at very small x, parton distributions saturate at a scale Q_s(x). Since the occupation number is large, and $Q_s(x)>>Λ_{QCD}$, classical weak coupling methods may be used to study this novel regime of non-linear classical fields in QCD. In these lectures, we apply these methods to compute structure functions in deeply inelastic scattering (DIS) and the energy density of gluons produced in high energy nuclear collisions.

preprint1997arXiv

Boost covariant gluon distributions in large nuclei

It has been shown recently that there exist analytical solutions of the Yang-Mills equations for non-Abelian Weizsäcker-Williams fields which describe the distribution of gluons in large nuclei at small x. These solutions however depend on the color charge distribution at large rapidities. We here construct a model of the color charge distribution of partons in the fragmentation region and use it to compute the boost covariant momentum distributions of wee gluons. The phenomenological applications of our results are discussed.

preprint1995arXiv

Wee partons in large nuclei: from virtual dream to hard reality

We construct a weak coupling, many body theory to compute parton distributions in large nuclei for $x\ll A^{-1/3}$. The wee partons are highly coherent, non--Abelian Weizsäcker--Williams fields. Radiative corrections to the classical results are discussed. The parton distributions for a single nucleus provide the initial conditions for the dynamical evolution of matter formed in ultrarelativistic nuclear collisions.

preprint1993arXiv

Computing Quark and Gluon Distribution Functions for Very Large Nuclei

We argue that the distribution functions for quarks and gluons are computable at small x for sufficiently large nuclei, perhaps larger than can be physically realized. For such nuclei, we argue that weak coupling methods may be used. We show that the computation of the distribution functions can be recast as a many body problem with a modified propagator, a coupling constant which depends on the multiplicity of particles per unit rapidity per unit area, and for non-abelian gauge theories, some extra media dependent vertices. We explicitly compute the distribution function for gluons to lowest order, and argue how they may be computed in higher order.