
We present the nexttoleadingorder eventbyevent EKRT model predictions
for the centrality dependence of the charged hadron multiplicity in the
pseudorapidity interval $\eta\le 0.5$, and for the centrality dependence of
the charged hadron flow harmonics $v_n\{2\}$ obtained from 2particle
cumulants, in $\sqrt{s_{NN}}=5.44$ TeV Xe+Xe collisions at the CERN Large
Hadron Collider. Our prediction for the 05 \% central charged multiplicity is
$dN_{\rm ch}/d\eta =1218\pm 46$. We also predict $v_n\{2\}$ in Xe+Xe collisions
to increase more slowly from central towards peripheral collisions than those
in a Pb+Pb system. We find that at $10 \dots 50$\% centralities $v_2\{2\}$ is
smaller and $v_3\{2\}$ is larger than in the Pb+Pb system while $v_4\{2\}$ is
of the same magnitude in both systems. We also find that the ratio of flow
harmonics in Xe+Xe collisions and in Pb+Pb collisions shows a slight
sensitivity to the temperature dependence of the shearviscositytoentropy
ratio. As we discuss here, the new nuclear massnumber systematics especially
in the flow harmonics serves as a welcome further constraint for describing the
spacetime evolution of a heavyion system and for determining the shear
viscosity and other transport properties of strongly interacting matter.

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

Predictions made in Albacete {\it et al} prior to the LHC $p+$Pb run at
$\sqrt{s_{NN}} = 5$ TeV are compared to currently available data. Some
predictions shown here have been updated by including the same experimental
cuts as the data. Some additional predictions are also presented, especially
for quarkonia, that were provided to the experiments before the data were made
public but were too late for the original publication are also shown here.

We compute the initial energy densities produced in ultrarelativistic
heavyion collisions from NLO perturbative QCD using a saturation conjecture to
control soft particle production, and describe the subsequent spacetime
evolution of the system with hydrodynamics, event by event. The resulting
centrality dependence of the low$p_T$ observables from this pQCD + saturation
+ hydro ("EKRT") framework are then compared simultaneously to the LHC and RHIC
measurements. With such an analysis we can test the initial state calculation,
and constrain the temperature dependence of the shear viscositytoentropy
ratio $\eta/s$ of QCD matter. Using these constraints from the current RHIC and
LHC measurements we then predict the charged hadron multiplicities and flow
coefficients for the 5.023 TeV Pb+Pb collisions.

We compute predictions for various lowtransversemomentum bulk observables
in $\sqrt{s_{NN}} = 5.023$ TeV Pb+Pb collisions at the LHC from the
eventbyevent nexttoleadingorder perturbativeQCD + saturation + viscous
hydrodynamics ("EKRT") model. In particular, we consider the centrality
dependence of charged hadron multiplicity, flow coefficients of the
azimuthangle asymmetries and correlations of eventplane angles. The
centrality dependencies of the studied observables are predicted to be very
similar to those at 2.76 TeV, and the magnitudes of the flow coefficients and
eventplane angle correlations are predicted to be close to those at 2.76 TeV.
The flow coefficients may, however, offer slightly more discriminating power on
the temperature dependence of QCD matter viscosity than the 2.76 TeV
measurements. Our prediction for the multiplicity in the 05\% centrality
class, obtained using the two temperaturedependent shearviscositytoentropy
ratios that give the best overall fit to RHIC and LHC data is $dN_{\rm
ch}/d\eta\big_{\eta\le 0.5} =1876\dots2046$. We also predict a powerlaw
increase from 200 GeV Au+Au collisions at RHIC to 2.76 and 5.023 TeV Pb+Pb
collisions at the LHC, $dN_{\rm ch}/d\eta\big_{\eta\le 0.5} \propto
s^{0.164\dots0.174}$.

We introduce an eventbyevent pQCD + saturation + hydro ("EKRT") framework
for highenergy heavyion collisions, where we compute the produced fluctuating
QCDmatter energy densities from nexttoleading order (NLO) perturbative QCD
(pQCD) using saturation to control soft particle production, and describe the
spacetime evolution of the QCD matter with viscous hydrodynamics, event by
event (EbyE). We compare the computed centrality dependence of hadronic
multiplicities, p_T spectra and flow coefficients v_n against LHC and RHIC
data. We compare also the computed EbyE probability distributions of relative
fluctuations of v_n, as well as correlations of 2 and 3 eventplane angles,
with LHC data. Our systematic multienergy and observable analysis not only
tests the initial state calculation and applicability of hydrodynamics, but
also makes it possible to constrain the temperature dependence of the shear
viscositytoentropy ratio, eta/s(T), of QCD matter in its different phases.
Remarkably, we can describe all these different flow observables and
correlations consistently with eta/s(T) that is independent of the collision
energy.

We introduce an eventbyevent perturbativeQCD + saturation + hydro ("EKRT")
framework for ultrarelativistic heavyion collisions, where we compute the
produced fluctuating QCDmatter energy densities from nexttoleading order
perturbative QCD using a saturation conjecture to control soft particle
production, and describe the spacetime evolution of the QCD matter with
dissipative fluid dynamics, event by event. We perform a simultaneous
comparison of the centrality dependence of hadronic multiplicities, transverse
momentum spectra, and flow coefficients of the azimuthangle asymmetries,
against the LHC and RHIC measurements. We compare also the computed
eventbyevent probability distributions of relative fluctuations of elliptic
flow, and eventplane angle correlations, with the experimental data from Pb+Pb
collisions at the LHC. We show how such a systematic multienergy and
multiobservable analysis tests the initial state calculation and the
applicability region of hydrodynamics, and in particular how it constrains the
temperature dependence of the shear viscositytoentropy ratio of QCD matter in
its different phases in a remarkably consistent manner.

We review the recent results from the computation of saturated
nexttoleading order perturbative QCD minijet intial conditions combined with
viscous hydrodynamical evolution of ultrarelativistic heavyion collisions at
the LHC and RHIC. Comparison with experimental data is shown.

We calculate initial conditions for the hydrodynamical evolution in
ultrarelativistic heavyion collisions at the LHC and RHIC in an improved
nexttoleading order perturbative QCD + saturation framework. Using viscous
relativistic hydrodynamics, we show that we obtain a good simultaneous
description of the centrality dependence of charged particle multiplicities,
transverse momentum spectra and elliptic flow at the LHC and at RHIC. In
particular, we discuss how the temperature dependence of the shear viscosity is
constrained by these data.

Using nexttoleading order perturbative QCD and a conjecture of saturation
to suppress the production of lowenergy partons, we calculate the initial
energy densities and formation times for the dissipative fluid dynamical
evolution of the quarkgluon plasma produced in ultrarelativistic heavyion
collisions. We identify the framework uncertainties and demonstrate the
predictive power of the approach by a good global agreement with the measured
centrality dependence of charged particle multiplicities, transverse momentum
spectra and elliptic flow simultaneously for the Pb+Pb collisions at the LHC
and Au+Au at RHIC. In particular, the shear viscosity in the different phases
of QCD matter is constrained in this new framework simultaneously by all these
data.

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

The present note relies on the recently published conceptual design report of
the LHeC and extends the first contribution to the European strategy debate in
emphasising the role of the LHeC to complement and complete the high luminosity
LHC programme. The brief discussion therefore focuses on the importance of high
precision PDF and $\alpha_s$ determinations for the physics beyond the Standard
Model (GUTs, SUSY, Higgs). Emphasis is also given to the importance of high
parton density phenomena in nuclei and their relevance to the heavy ion physics
programme at the LHC.

This document provides a brief overview of the recently published report on
the design of the Large Hadron Electron Collider (LHeC), which comprises its
physics programme, accelerator physics, technology and main detector concepts.
The LHeC exploits and develops challenging, though principally existing,
accelerator and detector technologies. This summary is complemented by brief
illustrations of some of the highlights of the physics programme, which relies
on a vastly extended kinematic range, luminosity and unprecedented precision in
deep inelastic scattering. Illustrations are provided regarding high precision
QCD, new physics (Higgs, SUSY) and electronion physics. The LHeC is designed
to run synchronously with the LHC in the twenties and to achieve an integrated
luminosity of O(100) fb$^{1}$. It will become the cleanest high resolution
microscope of mankind and will substantially extend as well as complement the
investigation of the physics of the TeV energy scale, which has been enabled by
the LHC.

We bring the EKRT framework, which combines perturbative QCD (pQCD) minijet
production with gluon saturation and hydrodynamics, to nexttoleading order
(NLO) in pQCD as rigorously as possible. We chart the model uncertainties, and
study the viability and predictive power of the model in the light of the RHIC
and LHC measurements in central $A+A$ collisions. In particular, we introduce a
new set of measurement functions to define the infrared and collinearsafe
minijet transverse energy, $E_T$, in terms of which we formulate the
saturation. We update the framework with the EPS09 NLO nuclear parton
distributions (nPDFs), and study the propagation of the nPDF uncertainties into
the computed $E_T$, saturation scales and the finalstate multiplicities. The
key parameters, which need to be fixed using the measurements, are identified,
and their correlation is discussed. We convert the saturated minijet $E_T$ into
QCDmatter initial conditions for longitudinally boostinvariant ideal
hydrodynamics. We compute the chargedparticle multiplicities and identified
bulk hadron $p_T$ spectra in 5% most central Au+Au collisions at RHIC and Pb+Pb
at the LHC. We obtain an encouragingly good agreement with the experimental
data, simultaneously at RHIC and LHC, showing that the approach has a definite
predictive power.

The physics programme and the design are described of a new collider for
particle and nuclear physics, the Large Hadron Electron Collider (LHeC), in
which a newly built electron beam of 60 GeV, up to possibly 140 GeV, energy
collides with the intense hadron beams of the LHC. Compared to HERA, the
kinematic range covered is extended by a factor of twenty in the negative
fourmomentum squared, $Q^2$, and in the inverse Bjorken $x$, while with the
design luminosity of $10^{33}$ cm$^{2}$s$^{1}$ the LHeC is projected to
exceed the integrated HERA luminosity by two orders of magnitude. The physics
programme is devoted to an exploration of the energy frontier, complementing
the LHC and its discovery potential for physics beyond the Standard Model with
high precision deep inelastic scattering measurements. These are designed to
investigate a variety of fundamental questions in strong and electroweak
interactions. The physics programme also includes electrondeuteron and
electronion scattering in a $(Q^2, 1/x)$ range extended by four orders of
magnitude as compared to previous leptonnucleus DIS experiments for novel
investigations of neutron's and nuclear structure, the initial conditions of
QuarkGluon Plasma formation and further quantum chromodynamic phenomena. The
LHeC may be realised either as a ringring or as a linacring collider. Optics
and beam dynamics studies are presented for both versions, along with technical
design considerations on the interaction region, magnets and further
components, together with a design study for a high acceptance detector. Civil
engineering and installation studies are presented for the accelerator and the
detector. The LHeC can be built within a decade and thus be operated while the
LHC runs in its highluminosity phase. It thus represents a major opportunity
for progress in particle physics exploiting the investment made in the LHC.

In hydrodynamicalmodeling of heavyion collisions the initial state spatial
anisotropies translate into momentum anisotropies of the final state particle
distributions. Thus, understanding the origin of the initial anisotropies and
quantifying their uncertainties is important for the extraction of specific QCD
matter properties, such as viscosity, from the experimental data. In this work
we study the wounded nucleon approach in the Monte Carlo Glauber model
framework, focusing especially on the uncertainties which arise from the
modeling of the nucleonnucleon interactions between the colliding nucleon
pairs and nucleonnucleon correlations inside the colliding nuclei. We compare
the black disk model and a probabilistic profile function approach for the
inelastic nucleonnucleon interactions, and study the effects of initial state
correlations using stateoftheart modeling of these.

In hydrodynamical modeling of heavyion collisions, the initialstate spatial
anisotropies are translated into momentum anisotropies of the finalstate
particle distributions. Thus, understanding the origin of the initialstate
anisotropies and their uncertainties is important before extracting specific
QCD matter properties, such as viscosity, from the experimental data. In this
work we review the wounded nucleon approach based on the Monte Carlo Glauber
model, charting in particular the uncertainties arising from modeling of the
nucleonnucleon interactions between the colliding nucleon pairs and
nucleonnucleon correlations inside the colliding nuclei. We discuss the
differences between the black disk model and a probabilistic profile function
approach for the inelastic nucleonnucleon interactions, and investigate the
influence of initialstate correlations using stateoftheart modeling of
these.

We present a nexttoleading order (NLO) global DGLAP analysis of nuclear
parton distribution functions (nPDFs) and their uncertainties. Carrying out an
NLO nPDF analysis for the first time with three different types of experimental
input  deep inelastic $\ell$+A scattering, DrellYan dilepton production in
p+$A$ collisions, and inclusive pion production in d+Au and p+p collisions at
RHIC  we find that these data can well be described in a conventional
collinear factorization framework. Although the pion production has not been
traditionally included in the global analyses, we find that the shape of the
nuclear modification factor $R_{\rm dAu}$ of the pion $p_T$spectrum at
midrapidity retains sensitivity to the gluon distributions, providing evidence
for shadowing and EMCeffect in the nuclear gluons. We use the Hessian method
to quantify the nPDF uncertainties which originate from the uncertainties in
the data. In this method the sensitivity of $\chi^2$ to the variations of the
fitting parameters is mapped out to orthogonal error sets which provide a
userfriendly way to calculate how the nPDF uncertainties propagate to any
factorizable nuclear crosssection. The obtained NLO and LO nPDFs and the
corresponding error sets are collected in our new release called {\ttfamily
EPS09}. These results should find applications in precision analyses of the
signatures and properties of QCD matter at the LHC and RHIC.

In this talk, we introduce our recently completed nexttoleading order (NLO)
global analysis of the nuclear parton distribution functions (nPDFs) called
EPS09  a higher order successor to the wellknown leadingorder (LO) analysis
EKS98 and also to our previous LO work EPS08. As an extension to similar global
analyses carried out by other groups, we complement the data from deep
inelastic $l+A$ scattering and DrellYan dilepton measurements in p+$A$
collisions by inclusive midrapidity pion production data from d+Au collisions
at RHIC, which results in better constrained gluon distributions than before.
The most important new ingredient, however, is the detailed error analysis,
which employs the Hessian method and which allows us to map out the
parameterspace vicinity of the bestfit to a collection of nPDF error sets.
These error sets provide the enduser a way to compute how the
PDFuncertainties will propagate into the cross sections of his/her interest.
The EPS09 package to be released soon, will contain both the NLO and LO results
for the best fits and the uncertainty sets.

We present an improved leadingorder global DGLAP analysis of nuclear parton
distribution functions (nPDFs), supplementing the traditionally used data from
deep inelastic leptonnucleus scattering and DrellYan dilepton production in
protonnucleus collisions, with inclusive high$p_T$ hadron production data
measured at RHIC in d+Au collisions. With the help of an extended definition of
the $\chi^2$ function, we now can more efficiently exploit the constraints the
different data sets offer, for gluon shadowing in particular, and account for
the overall data normalization uncertainties during the automated $\chi^2$
minimization. The very good simultaneous fit to the nuclear hard process data
used demonstrates the feasibility of a universal set of nPDFs, but also
limitations become visible. The high$p_T$ forwardrapidity hadron data of
BRAHMS add a new crucial constraint into the analysis by offering a direct
probe for the nuclear gluon distributions  a sector in the nPDFs which has
traditionally been very badly constrained. We obtain a strikingly stronger
gluon shadowing than what has been estimated in previous global analyses. The
obtained nPDFs are released as a parametrization called EPS08.

We use perfectfluid hydrodynamical model to predict the elliptic flow
coefficients in Pb + Pb collisions at the Large Hadron Collider (LHC). The
initial state for the hydrodynamical calculation for central $A + A$ collisions
is obtained from the perturbative QCD + saturation (EKRT) model. The centrality
dependence of the initial state is modeled by the optical Glauber model. We
show that the baseline results obtained from the framework are in good
agreement with the data from the Relativistic Heavy Ion Collider (RHIC), and
show predictions for the $p_T$ spectra and elliptic flow of pions in Pb + Pb
collisions at the LHC. Also mass and multiplicity effects are discussed.

This writeup is a compilation of the predictions for the forthcoming Heavy
Ion Program at the Large Hadron Collider, as presented at the CERN Theory
Institute 'Heavy Ion Collisions at the LHC  Last Call for Predictions', held
from May 14th to June 10th 2007.

In this talk, we shortly report results from our recent global DGLAP analysis
of nuclear parton distributions. This is an extension of our former
EKS98analysis improved with an automated $\chi^2$ minimization procedure and
uncertainty estimates. Although our new analysis show no significant deviation
from EKS98, a sign of a significantly stronger gluon shadowing could be seen in
the RHIC BRAHMS data.

In this talk, we present the results from our recent global reanalysis of
nuclear parton distribution functions (nPDFs), where the DGLAPevolving nPDFs
are constrained by nuclear hard process data from deep inelastic $l+A$
scattering (DIS) and the DrellYan (DY) process in $p+A$ collisions, and by sum
rules. The main improvements over our earlier work {\em EKS98} are the
automated $\chi^2$ minimization, better controllable fit functions and
possibility for error estimates. The obtained 16parameter fit to N=514
datapoints is good, $\chi^2/{\rm d.o.f}=0.82$. Fit quality comparison and the
error estimates obtained show that the old {\em EKS98} parametrization is fully
consistent with the present automated reanalysis. Comparison with other global
nPDF analyses is presented as well. Within the DGLAP framework we also discuss
the possibility of incorporating a clearly stronger gluon shadowing, which is
suggested by the RHIC BRAHMS data from d+Au collisions.

By studying the PaschosWolfenstein (PW) ratio of deep inelastic $\nu{\rm
Fe}$ and $\bar \nu{\rm Fe}$ scattering cross sections, we show that it should
be possible to explain the NuTeV $\sin ^2 \theta_{\rm W}$ anomaly with quite
conventional physics, by introducing mutually different nuclear modifications
for the valence$u$ and valence$d$ quark distributions of the protons in iron.
Keeping the EKS98 nuclear modifications for $u_V+d_V$ as a baseline, we find
that some 2030 % nuclear modifications to the $u_V$ and $d_V$ distributions
account for the change induced in the PW ratio by the NuTeVsuggested increase
$\Delta \sin ^2 \theta_{\rm W}=0.005$. We show that introduction of such
nuclear modifications in $u_V$ and $d_V$ individually, does not lead into
contradiction with the present global DGLAP analyses of the nuclear parton
distributions, where deep inelastic $lA$ scattering data and DrellYan dilepton
data from $pA$ collisions are used as constraints. We thus suggest that the
NuTeV result serves as an important further constraint in pinning down the
nuclear effects of the bound nucleon PDFs. We also predict that if the NuTeV
anomaly is explained by this mechanism, the NOMAD experiment should see an
increase in the weak mixing angle quite close to the NuTeV result.