
We point out that current calculations of inclusive twoparticle correlations
in pA collisions based on the Color Glass Condensate approach exhibit a
contribution from HanburyBrownTwiss correlations. These HBT correlations are
quite distinct from the standard ones, in that they are apparent for particles
widely separated in rapidity. The transverse size of the emitter which is
reflected in these correlations is the gluonic size of the proton. This raises
an interesting possibility of measuring the proton size directly by the HBT
effect of particle pairs produced in pA collisions.

We show explicitly how the high energy QCD evolution generated by the KLWMIJ
Hamiltonian can be cast in the form of the QCD Reggeon Field Theory. We show
how to reduce the KLWMIJ Hamitonian to physical color singlet degrees of
freedom. We suggest a natural way of defining the Pomeron and other Reggeons in
the framework of the KLWMIJ evolution and derive the QCD Reggeon Field Theory
Hamiltonian which includes several lowest Reggeon operators. This Hamiltonian
generates evolution equations for all Reggeons in the case of dilutedense
scattering, including the nonlinear BalitskyKovchegov equation for the
Pomeron. We also find explicit expressions for the Reggeon conjugate operators
in terms of QCD operators, and derive their evolution equations. This provides
a natural and unambiguous framework for reggeization procedure introduced in
\cite{BW, BE}. The Bartels triple Pomeron vertex is inherited directly from the
RFT Hamiltonian. For simplicity in the bulk of the paper we work in the large
$N_c$ limit.

We point out that the scenario for UV completion by "classicalization",
proposed recently is in fact Wilsonian in the classical Wilsonian sense. It
corresponds to the situation when a field theory has a nontrivial UV fixed
point governed by a higher dimensional operator. Provided the kinetic term is a
relevant operator around this point the theory will flow in the IR to the free
scalar theory. Physically, "classicalization", if it can be realized, would
correspond to a situation when the fluctuations of the field operator in the UV
are smaller than in the IR. As a result there exists a clear tension between
the "classicalization" scenario and constraints imposed by unitarity on a
quantum field theory, making the existence of classicalizing unitary theories
questionable.

Time evolution of a "little bang" created in heavy ion collisions can be
divided into two phases, the preequilibrium and hydrodynamic. At what moment
the evolution becomes hydrodynamic and is there any universality in the
hydrodynamic flow? To answer these questions we briefly discuss various
versions of hydrodynamics and their applicability conditions. In particular, we
elaborate on the idea of "universal" (allorder resumed) hydrodynamics and
propose a simple new model for it. The model is motivated by results obtained
recently via the AdS/CFT correspondence. Finally, charged hadron multiplicities
in heavy ion collisions at the RHIC and LHC are discussed. At the freezout, the
multiplicities can be related to total entropy produced in the collision.
Assuming the universal hydrodynamics to hold, we calculate the entropy
production in the hydro stage of the collision. We end up speculating about a
connection between the multiplicity growth and the temperature dependence of
the QGP viscosity.

The HERA electronproton collider has collected 100 pb$^{1}$ of data since
its startup in 1992, and recently moved into a highluminosity operation mode,
with upgraded detectors, aiming to increase the total integrated luminosity per
experiment to more than 500 pb$^{1}$. HERA has been a machine of excellence
for the study of QCD and the structure of the proton. The Large Hadron Collider
(LHC), which will collide protons with a centreofmass energy of 14 TeV, will
be completed at CERN in 2007. The main mission of the LHC is to discover and
study the mechanisms of electroweak symmetry breaking, possibly via the
discovery of the Higgs particle, and search for new physics in the TeV energy
scale, such as supersymmetry or extra dimensions. Besides these goals, the LHC
will also make a substantial number of precision measurements and will offer a
new regime to study the strong force via perturbative QCD processes and
diffraction. For the full LHC physics programme a good understanding of QCD
phenomena and the structure function of the proton is essential. Therefore, in
March 2004, a oneyearlong workshop started to study the implications of HERA
on LHC physics. This included proposing new measurements to be made at HERA,
extracting the maximum information from the available data, and
developing/improving the theoretical and experimental tools. This report
summarizes the results achieved during this workshop.

We develop the formalizm for calculating semi inclusive observables at high
energy in the JIMWLK/KLWMIJ approach. We consider several examples including
diffractive processes, elastic and inclusive over the target degrees of
freedom, scattering with fixed total transverse momentum transfer and inclusive
gluon production. We discuss evolution of these observables with respect to
various rapidity variables involved in their definitions (total rapidity,
rapidity gap, width of diffractive interval etc.). We also discuss the dipole
model limit of these observables and derive closed simple (as opposed to
functional) differential equations in this approximation. We point out that
there are some interesting differences between the full JIMWLK/KLWMIJ evolution
and the dipole model evolution of diffractive cross section. In particular we
show that in the dipole approximation the target does not diffract beyond the
valence rapidity interval, consistently with the intuition suggested by the
Pomeron fan diagramms. On the other hand such diffractive processes do exist in
the full JIMWLK/KLWMIJ approach, albeit suppressed by the factor 1/N_c^2.

The HERA electronproton collider has collected 100 pb$^{1}$ of data since
its startup in 1992, and recently moved into a highluminosity operation mode,
with upgraded detectors, aiming to increase the total integrated luminosity per
experiment to more than 500 pb$^{1}$. HERA has been a machine of excellence
for the study of QCD and the structure of the proton. The Large Hadron Collider
(LHC), which will collide protons with a centreofmass energy of 14 TeV, will
be completed at CERN in 2007. The main mission of the LHC is to discover and
study the mechanisms of electroweak symmetry breaking, possibly via the
discovery of the Higgs particle, and search for new physics in the TeV energy
scale, such as supersymmetry or extra dimensions. Besides these goals, the LHC
will also make a substantial number of precision measurements and will offer a
new regime to study the strong force via perturbative QCD processes and
diffraction. For the full LHC physics programme a good understanding of QCD
phenomena and the structure function of the proton is essential. Therefore, in
March 2004, a oneyearlong workshop started to study the implications of HERA
on LHC physics. This included proposing new measurements to be made at HERA,
extracting the maximum information from the available data, and
developing/improving the theoretical and experimental tools. This report
summarizes the results achieved during this workshop.

A second workshop on small x physics, within the Small x Collaboration, was
held in Lund in June 2002 with the aim of overviewing recent theoretical
progress in this area and summarizing the experimental status.
This paper is dedicated to the memory of Jan Kwiecinski, who died
unexpectedly on August 29, 2003.

We calculate the high energy behavior of quarkantiquark exchange in
$\gamma^* \gamma^*$ elastic scattering by summing, to all orders, the leading
double logarithmic contributions of the QCD ladder diagrams. Motivation comes
from the LEP data for $\sigma_{tot}^{\gamma^* \gamma^*}$ which indicate the
need for secondary reggeon exchange. We show that, for large photon
virtualities, this exchange is calculable in pQCD. This applies, in particular,
to parts of the LEP kinematic region.

We calculate inclusive photoproduction crosssections for both proton and
nuclei targets, using a saturation model derived from an approximate solution
to the BalitskyKovchegov nonlinear evolution equation. This paper extends our
hypothesis, that one can successfully use the dipole picture and saturation
models to describe soft processes, instead of the soft Pomeron Regge
parametrization. Our fit is compatible with our previous fit to DIS data, and
utilizes the same phenomenological parameters as in our paper devoted to soft
hadronic interactions. Using the Glauber formalism, and no additional
parameters, we calculate the crosssections for photoproduction on various
nuclei, and compare our results with the relevant data.

We pose the question: how far can we push the QCD dipole model and the idea
of QCD saturation to the kinematical regions traditionally associated with soft
nonperturbative physics? The answer is that it works sufficiently well allowing
a smooth transition from hard (BFKL) pomeron to the soft one.

We apply our solution of the nonlinear evolution equation to the case of
$J/\psi$ photo and DIS production on nucleons and nuclei targets. The
uncertainty in the $J/\psi$ wave function normalization due to Fermi motion is
treated as a free parameter. We obtain good reproduction of the HERA
experimental data on a proton target. Calculations of $J/\psi$ mesons coherent
production on nuclei targets are presented and discussed. Our analysis supports
the conclusions reached in our previous studies, stressing the importance of
nonlinear evolution in the kinematical domain of highdensity QCD.

We attempt to describe soft hadron interactions in the framework of
saturation models, one based upon the BalitskyKovchegov nonlinear equation
and another one due to GolecBiernat and W\"{u}sthoff. For $pp$, $Kp$, and $\pi
p$ scattering the relevant hadronic wave functions are formulated, and total,
elastic crosssections, and the forward elastic slope are calculated and
compared to experimental data. The saturation mechanism leads to reasonable
reproduction of the data for the quantities analyzed, except for the forward
elastic slope, where the predicted increase with energy is too moderate.

The process of single diffractive dissociation off nuclei is considered on a
basis of solutions to the nonlinear evolution equation. The relevant saturation
scales $Q_{s A}^D(x,x_0)$ are determined and their dependences on Bjorken $x$,
atomic number $A$, and minimal rapidity gap $Y_0\equiv \ln 1/x_0$ are
investigated. The solutions are shown to posses the geometrical scaling in a
broad kinematic region for $x$ well below $x_0$. A detailed theoretical
discussion of saturation scale and scaling is presented.
The ratio $\sigma^{diff}/\sigma^{tot}$ is computed for several nuclei. We
predict that at $x\simeq 10^{4}$ this ratio is of the order 25% which is much
larger to the one of proton. This result indicates a possibility to observe a
very strong nuclear shadowing.

A new approach to global QCD analysis is developed. The main ingredients are
two QCDbased evolution equations. The first one is the BalitskyKovchegov
nonlinear equation, which sums higher twists while preserving unitarity. The
second equation is linear and it is responsible for the correct short distance
behavior of the theory, namely it includes the DGLAP kernel. Our approach
allows extrapolation of the parton distributions to very high energies
available at the LHC as well as very low photon virtualities, $Q^2\ll 1 {\rm
GeV^2}$.
All existing low $x$ data on the $F_2$ structure function is reproduced using
one fitting parameter. The resulting $\chi^2/df=1$.
Analyzing the parameter $\lambda\equiv \partial\ln F_2/\partial(\ln 1/x)$ at
very low $x$ and $Q^2$ well below $1 {\rm GeV^2}$ we find $\lambda \simeq 0.08
 0.1$. A result which agrees with the "soft pomeron" intercept without
involving soft physics.

We discuss the possibility that mass effects, beyond phase space, may have
substantial influence on the direct determination of $V_{tq}$ in decays of top
quarks. In principle, our results are valid both for singly produced top and
for $t {\bar t}$ pair production. The mass effect is practically irrelevant for
the extraction of $V_{tb}$ from Tevatron data, but it may have implications for
higher energy multijet top decay.

The ratio of the diffractive production to the total cross section in DIS is
computed as a function of the produced mass. The analysis is based on the
solution to the nonlinear evolution equation for the diffraction dissociation
in DIS.
The obtained ratios almost do not depend on the central mass energy in
agreement with the HERA experimental data. This independence is argued to be a
consequence of the scaling phenomena displayed by the cross sections.
As a weakness point a significant discrepancy between the data and the
obtained results is found in the absolute values of the ratios. Several
explanatory reasons are discussed.

This paper presents the first numerical solution to the nonlinear evolution
equation for diffractive dissociation processes in deep inelastic scattering.
It is shown that the solution depends on one scaling variable $\tau = Q^2/Q^{D
2}_s(x,x_0)$, where $Q^D_s(x,x_0)$ is the saturation scale for the diffraction
processes. The dependence of the saturation scale $Q^D_s(x,x_0)$ on both $x$
and $x_0$ is investigated, ($Y_0 = \ln(1/x_0)$ is a minimal rapidity gap for
the diffraction process). The $x$  dependence of $Q^D_s$ turns out to be the
same as of the saturation scale in the total inclusive DIS cross section. In
our calculations $Q^D_s(x,x_0)$ reveals only mild dependence on $x_0$. The
scaling is shown to hold for $x \ll x_0$ but is violated at $ x \sim x_0$.

The numerical solutions of the nonlinear evolution equation are shown to
display the ``geometric'' scaling recently discovered in the experimental data.
The phenomena hold both for proton and nucleus targets for all $x$ below
$10^{2}$ and $0.25 {\rm GeV^{2}}\le Q^2 \le 2.5\times10^3 {\rm GeV^{2}}$. The
scaling is practically exact (few percent error) in the saturation region. In
addition, an approximate scaling is found in the validity domain of the linear
evolution where it holds with about 10% accuracy.
Basing on the scaling phenomena we determine the saturation scale $Q_s(x)$
and study both its $x$dependence and the atomic number dependence for the
nuclei.

We suggest a new procedure for extrapolating the parton distributions from
HERA to much higher energies. The procedure suggested consists of two steps.
First, we solve the nonlinear evolution equation. Second, we introduce a
correcting function for which we write a DGLAPtype linear evolution equation.
The nonlinear equation is solved numerically and estimates for the saturation
scale, as well as for the gluon density at THERA and LHC energies are made.

We present the numerical solution of the nonlinear evolution equation for
DIS on nuclei for $x = 10^{2} \div 10^{7}$. We demonstrate that the solution
to the nonlinear evolution equation is quite different from the Glauber 
Mueller formula which was used as the initial condition for the equation.
We illustrate the energy profit for performing DIS experiments on nuclei.
However, it turns out that the gain is quite modest: $x_{Au} \simeq 5 x_{\rm
proton} $ for the same parton density.
We find that the saturation scale $Q^2_s \propto A^{{1/3}}$. For gold the
saturation scale $Q_{s,Au} \simeq 1.5 GeV$ at $x= 10^{3}$. Such a large value
leads to considerable contribution of the high density QCD phase to RHIC data
and reveals itself in essential damping for both $xG_A$ and $F_{2A}$.

We suggest a new procedure for extrapolating the parton distributions from
HERA energies to higher energies at THERA and LHC. The procedure suggested
consists of two steps: first, we solve the nonlinear evolution equation which
includes the higher twists contributions, however this equation is deficient
due to the low ($log(1/x)$) accuracy of our calculations. Second, we introduce
a correcting function for which we write a DGLAP type linear evolution
equation. We show that this correcting function is small in the whole kinematic
region and decreases at low $x$. The nonlinear evolution equation is solved
numerically and first estimates for the saturation scale, as well as for the
value of the gluon density at THERA and LHC energies are made. We show that
nonlinear effects lead to damping of the gluon density by a factor of
$2 \div 3$ at $x \approx 10^{7}$.

These notes are a summary of our efforts to answer the question in the title.
Our answer is in the affirmative as: (i) HERA data indicate a large value of
the gluon structure function; (ii) no contradictions with the asymptotic
predictions of high density QCD have been observed; and (iii) the numerical
estimates of our model give a natural description of the size of deviation from
the routine DGLAP explanation. We discuss the alternative approaches and
possible new experiments.

We gereralize the KochegovMcLerran formula for the ratio
$\sigma^{DD}/\sigma_{tot}$ in perturbative QCD, using MuellerGlauber approach
for shadowing corrections and AGK cutting rules. We investigate several
phenomenological approaches with the goal of obtaining results consisent with
experimental data. We fail to reproduce the observed weak energy dependence of
the ratio, and conclude that the soft nonperturbative contribution present at
short distances must also be included.

We calculate, in the GlauberMueller approach, the ratio of the diffractive
dissociation cross section to the total cross section on nuclei. We observe
that shadowing corrections (SC) are significant, in the calculated $\ddt$ ratio
but they do not lead to a smooth energy dependence of the ratio for a proton
target as was observed at HERA.