Basis Light-front 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 light-front 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 heavy-ion 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 light-front quantum field theory in a basis function approach to solve the non-perturbative 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 two-scale system. The weak external field limit of the anomalous magnetic moment agrees with the result of QED perturbation theory within the anticipated accuracy.

Non-perturbative solutions to the quantum-field 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 non-perturbative QED + external field problem of an electron in a strong transverse confining potential using Hamiltonian light-front quantum field theory in a basis function representation. The invariant mass spectra and the anomalous magnetic moment of the lowest state for this two-scale 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 light-front quantum field theory constitutes a framework for the non-perturbative solution of invariant masses and correlated parton amplitudes of self-bound systems. By choosing the light-front 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 no-core methods of nuclear many-body 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 two-dimensional harmonic oscillator basis for transverse modes that corresponds with eigensolutions of the soft-wall AdS/QCD model obtained from light-front holography. We outline our approach and present illustrative features of some non-interacting 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.

Non-perturbative Hamiltonian light-front 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 light-front gauge and a basis function representation with two-dimensional harmonic oscillator basis for transverse modes that corresponds with eigensolutions of the soft-wall AdS/QCD model obtained from light-front holography.

Hamiltonian light-front quantum field theory constitutes a framework for the non-perturbative solution of invariant masses and correlated parton amplitudes of self-bound systems. By choosing light-front 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 no-core methods of nuclear many-body 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 Self-Organizing Maps (SOMs). SOMs are visualization algorithms based on competitive learning among spatially-ordered neurons. Our SOMs are trained with stochastically generated PDF samples. On every optimization iteration the PDFs are clustered on the SOM according to a user-defined 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 Self-Organizing Maps (SOMs). SOMs are a class of clustering algorithms based on competitive learning among spatially-ordered 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 non-singlet 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 non-zero 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 boost-invariant 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 net-baryon 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 Self-Organizing Maps (SOMs) which at variance with the standard Neural Networks, are based on competitive-learning. SOMs generate a non-uniform 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 boost-invariant hydrodynamics with initial energy and net-baryon 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 boost-invariant hydrodynamics with initial energy and net-baryon number densities from the perturbative QCD (pQCD)+saturation model. Build-up 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 quark-gluon 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 net-baryon 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 high-pT 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 pT-independent both for RHIC (in agreement with data) and for the LHC. This is a consequence of the shape of the partonic pT-spectrum 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 time-averaged 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 semi-hard 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 power-like 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 cms-energies 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 lowest-order 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 charged-particle 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 lepton--nucleus 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.