
We evaluate the transverse momentum spectrum of J/Psi (up to semihard
momenta) in pA and AA collisions taking into account only the initial state
effects, but resumming them to all orders in \alpha_s^2 A^{1/3}. In our
previous papers we noticed that cold nuclear matter effects alone could not
explain the experimental data on rapidity and centrality dependencies of the
J/Psi yield in AA collisions indicating the existence of an additional
suppression mechanism. Our present calculations indicate that the discrepancy
persists and even increases at semihard transverse momenta, implying a
significant final state effect on J/Psi production in this kinematical domain.
The QCD dipole model we employ is only marginally applicable for J/Psi
production at midrapidity at RHIC energies but its use is justified in the
forward rapidity region. At LHC energies we can quantitatively evaluate the
magnitude of cold nuclear matter effects in the entire kinematical region of
interest. We present our calculations of J/Psi transverse momentum spectra in
pA and AA collisions at LHC and RHIC energies.

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.

We derive the formulae for the cross section of J/Psi production in high
energy pA and AA collisions taking into account the gluon saturation/color
glass condensate effects. We then perform the numerical calculations of the
corresponding nuclear modification factors and find a good agreement between
our calculations and the experimental data on J/Psi production in dA
collisions. We also observe that cold nuclear modification effects alone cannot
describe the data on J/Psi production in AA collisions. Additional final state
suppression (at RHIC) and enhancement (at LHC) mechanisms are required to
explain the experimental observations.

We extract the bulk viscosity of hot quarkgluon matter in the presence of
light quarks from the recent lattice data on the QCD equation of state. For
that purpose we extend the sum rule analysis by including the contribution of
light quarks. We also discuss the universal properties of bulk viscosity in the
vicinity of a second order phase transition, as it might occur in the chiral
limit of QCD at fixed strange quark mass and most likely does occur in
twoflavor QCD. We point out that a chiral transition in the O(4) universality
class at zero baryon density as well as the transition at the chiral critical
point which belongs to the Z(2) universality class both lead to the critical
behavior of bulk viscosity. In particular, the latter universality class
implies the divergence of the bulk viscosity, which may be used as a signature
of the critical point. We discuss the physical picture behind the dramatic
increase of bulk viscosity seen in our analysis, and devise possible
experimental tests of related phenomena.

Kubo's formula relates bulk viscosity to the retarded Green's function of the
trace of the energymomentum tensor. Using low energy theorems of QCD for the
latter we derive the formula which relates the bulk viscosity to the energy
density and pressure of hot matter. We then employ the available lattice QCD
data to extract the bulk viscosity as a function of temperature. We find that
close to the deconfinement temperature bulk viscosity becomes large, with
viscositytoentropy ratio zeta/s about 1.

We consider the production of open charm in heavy ion collisions in the
framework of the Color Glass Condensate. In the central rapidity region at
RHIC, for the charm quark yield we expect N(coll) (number of collisions)
scaling in the absence of finalstate effects. At higher energies, or forward
rapidities at RHIC, the saturation scale exceeds the charm quark mass; we find
that this results in the approximate N(part) (number of participants) scaling
of charm production in AA collisions and N(part)^A scaling in p(d)A collisions,
similarly to the production of high pT gluons discussed earlier. We also show
that the saturation phenomenon makes spectra harder as compared to the naive
parton model approach. We then discuss the energy loss of charm quarks in hot
and cold media and argue that the hardness of the spectrum implies very slow
dependence of the quenching factor on pT.

QCD at the classical level possesses scale invariance which is broken by
quantum effects. This ``dimensional transmutation'' phenomenon can be
mathematically described by formulating classical gluodynamics in a curved,
conformally flat, spacetime with nonvanishing cosmological constant. We study
QCD highenergy scattering in this theory. We find that the properties of the
scattering amplitude at small momentum transfer are determined by the energy
density of vacuum fluctuations. The approach gives rise to the power growth of
the total hadronhadron cross section with energy, i.e. the Pomeron. The
intercept of the Pomeron and the multiplicity of produced particles are
evaluated. We also speculate about a possible link between conformal anomaly
and TeV scale quantum gravity.

The nonlinear evolution equation for the scattering amplitude of colour
dipole off the heavy nucleus is solved in the double logarithmic approximation.
It is found that if the initial parton density in a nucleus is smaller then
some critical value, then the scattering amplitude is a function of one scaling
variable inside the saturation region, whereas if it is greater then the
critical value, then the scaling behaviour breaks down. Dependence of the
saturation scale on the number of nucleons is discussed as well.

These notes are a summary of our predictions for the new THERA project,
related to deep inelastic scattering in the region of ultra low $x$ ($x \to
10^{7}$). We collect here predictions that satisfy two criteria (i) they do
not depend on specific features of the model that we have to use to estimate a
possible effect; and (ii) they do not contradict the HERA data.

We develop a new approach for solving the nonlinear evolution equation in
the low $x_B$ region and show that the remarkable "geometric" scaling of its
solution holds not only in the saturation region, but in much wider kinematical
region. This is in a full agreement with experimental data (GolecBiernat,
Kwiecinski and Stasto).

We calculate the total cross section for two scalar particles scattering at
high energies in $\lambda\phi^3$ theory in five dimensions, four of which are
usual Minkowski ones and the fifth is compact. It is shown that the cross
section is dominated by exchange of Pomeron whose intercept is larger than in
usual fourdimensional case.