We present the next-to-leading-order event-by-event 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 2-particle cumulants, in $\sqrt{s_{NN}}=5.44$ TeV Xe+Xe collisions at the CERN Large Hadron Collider. Our prediction for the 0-5 \% 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 shear-viscosity-to-entropy ratio. As we discuss here, the new nuclear mass-number systematics especially in the flow harmonics serves as a welcome further constraint for describing the space-time evolution of a heavy-ion 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, Drell-Yan 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 heavy-ion collisions from NLO perturbative QCD using a saturation conjecture to control soft particle production, and describe the subsequent space-time 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 viscosity-to-entropy 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 low-transverse-momentum bulk observables in $\sqrt{s_{NN}} = 5.023$ TeV Pb+Pb collisions at the LHC from the event-by-event next-to-leading-order perturbative-QCD + saturation + viscous hydrodynamics ("EKRT") model. In particular, we consider the centrality dependence of charged hadron multiplicity, flow coefficients of the azimuth-angle asymmetries and correlations of event-plane 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 event-plane 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 0-5\% centrality class, obtained using the two temperature-dependent shear-viscosity-to-entropy 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 power-law 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 event-by-event pQCD + saturation + hydro ("EKRT") framework for high-energy heavy-ion collisions, where we compute the produced fluctuating QCD-matter energy densities from next-to-leading order (NLO) perturbative QCD (pQCD) using saturation to control soft particle production, and describe the space-time 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 event-plane angles, with LHC data. Our systematic multi-energy 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 viscosity-to-entropy 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 event-by-event perturbative-QCD + saturation + hydro ("EKRT") framework for ultrarelativistic heavy-ion collisions, where we compute the produced fluctuating QCD-matter energy densities from next-to-leading order perturbative QCD using a saturation conjecture to control soft particle production, and describe the space-time 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 azimuth-angle asymmetries, against the LHC and RHIC measurements. We compare also the computed event-by-event probability distributions of relative fluctuations of elliptic flow, and event-plane angle correlations, with the experimental data from Pb+Pb collisions at the LHC. We show how such a systematic multi-energy and multi-observable analysis tests the initial state calculation and the applicability region of hydrodynamics, and in particular how it constrains the temperature dependence of the shear viscosity-to-entropy ratio of QCD matter in its different phases in a remarkably consistent manner.

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

We calculate initial conditions for the hydrodynamical evolution in ultrarelativistic heavy-ion collisions at the LHC and RHIC in an improved next-to-leading 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 next-to-leading order perturbative QCD and a conjecture of saturation to suppress the production of low-energy partons, we calculate the initial energy densities and formation times for the dissipative fluid dynamical evolution of the quark-gluon plasma produced in ultrarelativistic heavy-ion 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 electron-ion 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 next-to-leading 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 collinear-safe 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 final-state 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 QCD-matter initial conditions for longitudinally boost-invariant ideal hydrodynamics. We compute the charged-particle 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 four-momentum 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 electron-deuteron and electron-ion scattering in a $(Q^2, 1/x)$ range extended by four orders of magnitude as compared to previous lepton-nucleus DIS experiments for novel investigations of neutron's and nuclear structure, the initial conditions of Quark-Gluon Plasma formation and further quantum chromodynamic phenomena. The LHeC may be realised either as a ring-ring or as a linac-ring 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 high-luminosity phase. It thus represents a major opportunity for progress in particle physics exploiting the investment made in the LHC.

In hydrodynamicalmodeling of heavy-ion 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 nucleon-nucleon interactions between the colliding nucleon pairs and nucleon-nucleon correlations inside the colliding nuclei. We compare the black disk model and a probabilistic profile function approach for the inelastic nucleon-nucleon interactions, and study the effects of initial state correlations using state-of-theart modeling of these.

In hydrodynamical modeling of heavy-ion collisions, the initial-state spatial anisotropies are translated into momentum anisotropies of the final-state particle distributions. Thus, understanding the origin of the initial-state 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 nucleon-nucleon interactions between the colliding nucleon pairs and nucleon-nucleon correlations inside the colliding nuclei. We discuss the differences between the black disk model and a probabilistic profile function approach for the inelastic nucleon-nucleon interactions, and investigate the influence of initial-state correlations using state-of-the-art modeling of these.

We present a next-to-leading 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, Drell-Yan 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 EMC-effect 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 user-friendly way to calculate how the nPDF uncertainties propagate to any factorizable nuclear cross-section. 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 next-to-leading order (NLO) global analysis of the nuclear parton distribution functions (nPDFs) called EPS09 - a higher order successor to the well-known leading-order (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 Drell-Yan 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 parameter-space vicinity of the best-fit to a collection of nPDF error sets. These error sets provide the end-user a way to compute how the PDF-uncertainties 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 leading-order global DGLAP analysis of nuclear parton distribution functions (nPDFs), supplementing the traditionally used data from deep inelastic lepton-nucleus scattering and Drell-Yan dilepton production in proton-nucleus 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$ forward-rapidity 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 perfect-fluid 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 EKS98-analysis 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 DGLAP-evolving nPDFs are constrained by nuclear hard process data from deep inelastic $l+A$ scattering (DIS) and the Drell-Yan (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 16-parameter 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 Paschos-Wolfenstein (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 20-30 % nuclear modifications to the $u_V$ and $d_V$ distributions account for the change induced in the PW ratio by the NuTeV-suggested 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 Drell-Yan 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.