
Basis Lightfront Quantization (BLFQ) has recently been developed as a
promising nonperturbative technique. Using BLFQ, we investigate the Generalized
Parton Distributions (GPDs) in a nonperturbative framework for a dressed
electron in QED. We evaluate lightfront wave functions and carry out overlap
calculations to obtain GPDs. We also perform perturbative calculations in the
corresponding basis spaces to demonstrate that they compare reasonably with the
BLFQ results.

The presence of nuclear medium and collective phenomena which involve several
nucleons modify the parton distribution functions of nuclei (nPDFs) compared to
those of a free nucleon. These modifications have been investigated by
different groups using global analyses of high energy nuclear reaction world
data resulting in modern nPDF parametrizations with error estimates, such as
EPS09(s), HKN07 and nDS. These phenomenological nPDF sets roughly agree within
their uncertainty bands, but have antiquarks for large$x$ and gluons for the
whole $x$range poorly constrained by the available data. In the kinematics
accessible at the LHC this has negative impact on the interpretation of the
heavyion collision data, especially for the $p + A$ benchmarking runs. The EMC
region is also sensitive to the proper definition of $x$, where the nuclear
binding effects have to be taken into account, and for heavy nuclei one also
needs to take into account that a fraction of the nucleus momentum is carried
by the equivalent photons which modifies the momentum sum rule. We study how
these effects affect the predictions for the nuclear modification ratios at the
LHC kinematics using a model where we combine theoretical input for the leading
twist nuclear shadowing (the FGS model) and the EKS98s/EPS09s nPDF set where
the spatial dependence is formulated as a power series of the nuclear thickness
functions $T_A$.

We employ Hamiltonian lightfront quantum field theory in a basis function
approach to solve the nonperturbative problem of an electron in a strong
scalar transverse confining potential. We evaluate both the invariant mass
spectra and the anomalous magnetic moment of the lowest state for this
twoscale system. The weak external field limit of the anomalous magnetic
moment agrees with the result of QED perturbation theory within the anticipated
accuracy.

Nonperturbative solutions to the quantumfield theory is a topic of current
and broad interest, especially for the heavy ion and laser physics communities,
since they investigate particle production in the presence of strong external
fields. We have solved a nonperturbative QED + external field problem of an
electron in a strong transverse confining potential using Hamiltonian
lightfront quantum field theory in a basis function representation. The
invariant mass spectra and the anomalous magnetic moment of the lowest state
for this twoscale system are also evaluated. With this method the perturbative
QED results are also reproduced with a good accuracy. We also discuss the
extension of the method to other problems as well.

Hamiltonian lightfront quantum field theory constitutes a framework for the
nonperturbative solution of invariant masses and correlated parton amplitudes
of selfbound systems. By choosing the lightfront gauge and adopting a basis
function representation, we obtain a large, sparse, Hamiltonian matrix for mass
eigenstates of gauge theories that is solvable by adapting the ab initio
nocore methods of nuclear manybody theory. Full covariance is recovered in
the continuum limit, the infinite matrix limit. There is considerable freedom
in the choice of the orthonormal and complete set of basis functions with
convenience and convergence rates providing key considerations. Here, we use a
twodimensional harmonic oscillator basis for transverse modes that corresponds
with eigensolutions of the softwall AdS/QCD model obtained from lightfront
holography. We outline our approach and present illustrative features of some
noninteracting systems in a cavity. We illustrate the first steps towards
solving QED by obtaining the mass eigenstates of an electron in a cavity in
small basis spaces and discuss the computational challenges.

Nonperturbative Hamiltonian lightfront quantum field theory presents
opportunities and challenges that bridge particle physics and nuclear physics.
Fundamental theories, such as Quantum Chromodynmamics (QCD) and Quantum
Electrodynamics (QED) offer the promise of great predictive power spanning
phenomena on all scales from the microscopic to cosmic scales, but new tools
that do not rely exclusively on perturbation theory are required to make
connection from one scale to the next. We outline recent theoretical and
computational progress to build these bridges and provide illustrative results
for nuclear structure and quantum field theory. As our framework we choose
lightfront gauge and a basis function representation with twodimensional
harmonic oscillator basis for transverse modes that corresponds with
eigensolutions of the softwall AdS/QCD model obtained from lightfront
holography.

Hamiltonian lightfront quantum field theory constitutes a framework for the
nonperturbative solution of invariant masses and correlated parton amplitudes
of selfbound systems. By choosing lightfront gauge and adopting a basis
function representation, we obtain a large, sparse, Hamiltonian matrix for mass
eigenstates of gauge theories that is solvable by adapting the ab initio
nocore methods of nuclear manybody theory. Full covariance is recovered in
the continuum limit, the infinite matrix limit. We outline our approach and
discuss the computational challenges.

We propose a Parton Distribution Function (PDF) fitting technique which is
based on an interactive neural network algorithm using SelfOrganizing Maps
(SOMs). SOMs are visualization algorithms based on competitive learning among
spatiallyordered neurons. Our SOMs are trained with stochastically generated
PDF samples. On every optimization iteration the PDFs are clustered on the SOM
according to a userdefined feature and the most promising candidates are
selected as a seed for the subsequent iteration. Our main goal is thus to
provide a fitting procedure that, at variance with the global analyses and
standard neural network approaches, allows for an increased control of the
systematic bias by enabling user interaction in the various stages of the
fitting process.

Neural network algorithms have been recently applied to construct Parton
Distribution Function (PDF) parametrizations which provide an alternative to
standard global fitting procedures. We propose a technique based on an
interactive neural network algorithm using SelfOrganizing Maps (SOMs). SOMs
are a class of clustering algorithms based on competitive learning among
spatiallyordered neurons. Our SOMs are trained on selections of stochastically
generated PDF samples. The selection criterion for every optimization iteration
is based on the features of the clustered PDFs. Our main goal is to provide a
fitting procedure that, at variance with the standard neural network
approaches, allows for an increased control of the systematic bias by enabling
user interaction in the various stages of the process.

Following a previous detailed study of unpolarized generalized parton
distribution functions in the nonsinglet sector, and at zero values of the
skewness variable, $\zeta$, we propose a physically motivated parametrization
that is valid at $\zeta \neq 0$. Our method makes use of information from the
nucleon form factor data, from deep inelastuc scattering parton distribution
functions, and from lattice results on the Mellin moments of generalized parton
distributions. It provides, therefore, a step towards a model independent
extraction of generalized distributions from the data, alternative to the
mathematical ansatz of double distributions. Comparisons with recent
experimental data on the proton are shown.

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.

We propose a physically motivated parametrization for the unpolarized
generalized parton distributions, H and E, valid at both zero and nonzero
values of the skewness variable, \zeta. At \zeta=0, H and E are determined
using constraints from simultaneous fits of experimental data on both the
nucleon elastic form factors and the deep inelastic structure functions.
Lattice calculations of the higher moments constrain the parametrization at
\zeta > 0. Our method provides a step towards a model independent extraction of
generalized distributions from the data that is alternative to the mathematical
ansatz of double distributions.

We compute the initial energy density and net baryon number density in 5%
most central Pb+Pb collisions at $\sqrt s=5.5$ TeV from pQCD + (final state)
saturation, and describe the evolution of the produced system with
boostinvariant transversely expanding hydrodynamics. In addition to the total
multiplicity at midrapidity, we give predictions for the multiplicity of
charged hadrons, pions, kaons and (anti)protons, for the total transverse
energy and netbaryon number, as well as for the $p_T$spectrum of charged
hadrons, pions and kaons. We also predict the region of applicability of
hydrodynamics by comparing these results with high$p_T$ hadron spectra
computed from pQCD and energy losses.

We present an alternative algorithm to global fitting procedures to construct
Parton Distribution Functions (PDFs) parametrizations. The proposed algorithm
uses SelfOrganizing Maps (SOMs) which at variance with the standard Neural
Networks, are based on competitivelearning. SOMs generate a nonuniform
projection from a high dimensional data s pace onto a low dimensional one
(usually 1 or 2 dimensions) by clustering similar PDF representations together.
The SOMs are trained on progressively narrower selections of data samples. The
selection criterion is that of convergence towards a neighborhood of the
experimental data. All available data sets on deep inelastic scattering in the
kinematical region of 0.001 < x < 0.75, and 1 <Q^2 < 100 GeV^2, with a cut on
the final state invariant mass, W^2 > 10 GeV^2 were implemented. The proposed
fitting procedure, at variance with standard neural network approaches, allows
for an increased control of the systematic bias by enabling the user to
directly control the data selection procedure at various stages of the process.

We study the transverse spectra of hadrons in nearly central $AA$ collisions
at RHIC and LHC in a broad transverse momentum range Low$p_T$ spectra are
calculated by using boostinvariant hydrodynamics with initial energy and
netbaryon densities from the EKRT pQCD+saturation model. High$p_T$ spectra
are obtained from pQCD jet calculation including the energy loss of the parton
in the matter prior to its fragmentation to final hadrons.

We study the hadron spectra in nearly central $A$+$A$ collisions at RHIC and
LHC in a broad transverse momentum range. We cover the low$p_T$ spectra using
longitudinally boostinvariant hydrodynamics with initial energy and netbaryon
number densities from the perturbative QCD (pQCD)+saturation model. Buildup of
the transverse flow and sensitivity of the spectra to a single decoupling
temperature $\Tdec$ are studied. Comparison with RHIC data at $\ssNN=130$ and
200 GeV suggests a rather high value $\Tdec=150$ MeV. The high$p_T$ spectra
are computed using factorized pQCD cross sections, nuclear parton
distributions, fragmentation functions, and describing partonic energy loss in
the quarkgluon plasma by quenching weights. Overall normalization is fixed on
the basis of p+$\bar{\rm p}$(p) data and the strength of energy loss is
determined from RHIC Au+Au data. Uncertainties are discussed. With constraints
from RHIC data, we predict the $p_T$ spectra of hadrons in 5 % most central
Pb+Pb collisions at the LHC energy $\ssNN=5500$ GeV. Due to the closed
framework for primary production, we can also predict the netbaryon number at
midrapidity, as well as the strength of partonic energy losses at the LHC. Both
at the LHC and RHIC, we recognize a rather narrow crossover region in the $p_T$
spectra, where the hydrodynamic and pQCD fragmentation components become of
equal size. We argue that in this crossover region the two contributions are to
a good approximation mutually independent. In particular, our results suggest a
wider $p_T$region of applicability for hydrodynamical models at the LHC than
at RHIC.

We study the suppression of highpT hadron spectra in nuclear collisions,
supplementing the perturbative QCD factorized formalism with radiative parton
energy loss. We find that the nuclear modification factor which quantifies the
degree of suppression, is almost pTindependent both for RHIC (in agreement
with data) and for the LHC. This is a consequence of the shape of the partonic
pTspectrum in elementary collisions which implies that for the same value of
the nuclear modification factor at higher pT, an increasingly smaller fraction
of parton energy loss is needed. When the values of the timeaveraged transport
coefficient exceed 5 GeV^2/fm, the nuclear modification factor gradually loses
its sensitivity to the corresponding produced energy density. This is due to
particle production in the outer corona of the medium, which remains almost
unsuppressed even for extreme densities. Thus, even for the highest
experimentally accessible transverse momentum at the LHC and in contrast to
jets, the measurement of leading partons via leading hadrons is not a
penetrating probe of the dense matter, but a rather fragile probe which
fragments for the opacities reached below the skin of the medium. Relating the
transport coefficient to the energy density produced in the collision region,
we find and discuss a phenomenon reminiscent of the opacity problem of elliptic
flow: namely, the interaction of the hard parent parton with the medium appears
to be much stronger than that expected for perturbative interactions of the
hard parton with an ideal quark gluon plasma.

This manuscript is the outcome of the subgroup ``PDFs, shadowing and $pA$
collisions'' from the CERN workshop ``Hard Probes in Heavy Ion Collisions at
the LHC''. In addition to the experimental parameters for $pA$ collisions at
the LHC, the issues discussed are factorization in nuclear collisions, nuclear
parton distributions (nPDFs), hard probes as the benchmark tests of
factorization in $pA$ collisions at the LHC, and semihard probes as
observables with potentially large nuclear effects. Also, novel QCD phenomena
in $pA$ collisions at the LHC are considered. The importance of the $pA$
program at the LHC is emphasized.

The effects of the first nonlinear corrections to the DGLAP evolution
equations are studied by using the recent HERA data for the structure function
$F_2(x,Q^2)$ of the free proton and the parton distributions from CTEQ5L and
CTEQ6L as a baseline. By requiring a good fit to the H1 data, we determine
initial parton distributions at $Q_0^2=1.4$ GeV$^2$ for the nonlinear scale
evolution. We show that the nonlinear corrections improve the agreement with
the $F_2(x,Q^2)$ data in the region of $x\sim 3\cdot 10^{5}$ and $Q^2\sim 1.5$
GeV$^2$ without paying the price of obtaining a worse agreement at larger
values of $x$ and $Q^2$. For the gluon distribution the nonlinear effects are
found to play an increasingly important role at $x\lsim 10^{3}$ and
$Q^2\lsim10$ GeV$^2$, but rapidly vanish at larger values of $x$ and $Q^2$.
Consequently, contrary to CTEQ6L, the obtained gluon distribution at $Q^2=1.4$
GeV$^2$ shows a powerlike growth at small $x$. Relative to the CTEQ6L gluons,
an enhancement up to a factor $\sim6$ at $x=10^{5}$, $Q_0^2=1.4$ GeV$^2$
reduces to a negligible difference at $Q^2\gsim 10$ GeV$^2$.

The effects of the first nonlinear corrections to the DGLAP equations are
studied in light of the HERA data. Saturation limits are determined in the
DGLAP+GLRMQ approach for the free proton and for the Pb nucleus.

The DGLAP analyses of the nuclear parton distribution functions (nPDF) based
on the global fits to the data are reviewed, and the results from EKS98 and HKM
are compared. The usefulness of measuring hard probes in $pA$ collisions, at
the LHC in particular, is demonstrated.

We explain how the constraints from present experimental data can be used to
obtain the nPDF in the framework of LO DGLAP evolution. We will also compare
the only two available sets of this type and comment on the important
information that neutrino factories could provide.

The $Q^2$ dependence of the ratios of nuclear structure functions $F_2^A$ is
studied by performing QCD evolution of nuclear parton distribution functions.
The log $Q^2$ slope of these ratios is very sensitive to the nuclear gluon
distribution function. Taking different parametrizations, we show that the NMC
data on the $Q^2$ dependence of $F_2^{\rm Sn}/F_2^{\rm C}$ rule out the case
where nuclear shadowing (suppression) of gluons at $x\sim 0.01$ is much larger
than the shadowing observed in the ratio $F_2^A/F_2^{\rm D}$. We also take into
account modifications to the DGLAP evolution by including gluon fusion terms
and see that the effect is small at present energies, and, in any case, a
strong gluon shadowing is not favored. The region studied ($x \sim 0.01$) is
the most relevant for RHIC multiplicities.

We compute the distributions of charged particles at large transverse momenta
in $p\bar p(p)$, $pA$ and $AA$ collisions in the framework of perturbative QCD,
by using collinear factorization and the modern PDFs and fragmentation
functions. At the highest cmsenergies the shape of the spectra measured in
$p\bar p(p)$ collisions at large $q_T$ can be well explained. The difference
between the data and the lowestorder computation is quantified in terms of a
constant $K$factor for each energy. The $K$factor is found to systematically
decrease with growing $\sqrt s$. Also a lower limit for the partonic transverse
momentum, $p_0$, is extracted for each $\sqrt s$ based on the comparison with
the measurements. A systematic increase of $p_0$ as a function of $\sqrt s$ is
found. Nuclear effects in the chargedparticle spectra in $pA$ and $AA$
collisions at RHIC and LHC are studied in the framework of collinear
factorization by applying the EKS98 nuclear corrections to the parton
distributions. The nuclear effects are shown to mostly enhance the computed
spectra. A comparison with the recent PHENIX data from central and peripheral
Au+Au collisions at RHIC is done.

The $Q^2$ dependence of the ratios of the cross sections of deep inelastic
leptonnucleus scattering is studied in the framework of leading twist, lowest
order perturbative QCD. The $\log Q^2$ slope of the ratio $F_2^{\rm
Sn}/F_2^{\rm C}$ is computed by using the DGLAP evolution equations, and shown
to be sensitive to the nuclear gluon distribution functions. Four different
parametrizations for the nuclear effects of parton distributions are studied.
We show that the NMC data on the $Q^2$ dependence of $F_2^{\rm Sn}/F_2^{\rm C}$
rule out the case where nuclear shadowing (suppression) of gluons at $x\sim
0.01$ is much larger than the shadowing observed in the ratio $F_2^A/F_2^{\rm
D}$. We also show that the possible nonlinear correction terms due to gluon
fusion in the evolution equations do not change this conclusion. Some
consequences for computation of RHIC multiplicities, which probe the region
$x\gsim0.01$, are also discussed.