
The goal of this report is to give a comprehensive overview of the rich field
of forward physics, with a special attention to the topics that can be studied
at the LHC. The report starts presenting a selection of the Monte Carlo
simulation tools currently available, chapter 2, then enters the rich
phenomenology of QCD at low, chapter 3, and high, chapter 4, momentum transfer,
while the unique scattering conditions of central exclusive production are
analyzed in chapter 5. The last two experimental topics, Cosmic Ray and Heavy
Ion physics are presented in the chapter 6 and 7 respectively. Chapter 8 is
dedicated to the BFKL dynamics, multiparton interactions, and saturation. The
report ends with an overview of the forward detectors at LHC. Each chapter is
correlated with a comprehensive bibliography, attempting to provide to the
interested reader with a wide opportunity for further studies.

A Regge pole model for PomeronPomeron total cross section in the resonance
region $\sqrt{M^2}\le$ 5 GeV is presented. The cross section is saturated by
directchannel contributions from the Pomeron as well as from two different $f$
trajectories, accompanied by the isolated f$_0(500)$ resonance which dominates
the $\sqrt{M^{2}}\lesssim 1$ GeV region. A slowly varying background is taken
into account. The calculated PomeronPomeron total cross section cannot be
measured directly, but is an essential part of central diffractive processes.
In preparation of future calculations of central resonance production at the
hadron level, and corresponding measurements at the LHC, we normalize the
PomeronPomeron cross section at large masses $\sigma_{t}^{PP}
(\sqrt{M^2}\rightarrow \infty) \approx$ 1 mb as suggested by QCDmotivated
estimates.

We analyze a modification of the BFKL kernel for the adjoint representation
of the colour group in the maximally supersymmetric (N=4) YangMills theory in
the limit of a large number of colours, related to the modification of the
eigenvalues of the kernel suggested by S. Bondarenko and A. Prygarin in order
to reach the Hermitian separability of the eigenvalues. We restore the modified
kernel in the momentum space. It turns out that the modification is related
only to the real part of the kernel and that the correction to the kernel can
not be presented by a single analytic function in the entire momentum region,
which contradicts the known properties of the kernel.

We perform an examination of discontinuities of multiple production
amplitudes, which are required for further development of the BFKL approach. It
turns out that the discontinuities of 2 $\to$ 2 + n amplitudes obtained in the
BFKL approach contradict to the BDS ansatz for amplitudes with maximal helicity
violation in N = 4 supersymmetric YangMills theory with large number of
colours starting with n = 2. Explicit expressions for the discontinuities of
the 2 $\to$ 3 and 2 $\to$ 4 amplitudes in the invariant mass of pairs of
produced gluons are obtained in the planar N=4 SYM in the nexttoleading
logarithmic approximation. These expressions can be used for checking the
conjectured duality between the lightlike Wilson loops and the MHV amplitudes.

By using a Reggepole model for vector meson production (VMP), that
successfully describes the HERA data, we analyze the connection of VMP cross
sections in photoninduced reactions at HERA with those in ultraperipheral
collisions at the Large Hadron Collider (LHC). The role of the lowenergy
behaviour of VMP cross sections in {\gamma}p collisions is scrutinized.

We show that a growth of the protonproton total cross section with energy
can be entirely attributed to the purely perturbative mechanism. The infrared
regularization at rather short distances, $R_c\simeq 0.3$ fm, allows to extend
the BFKL technique from deep inelastic to hadronhadron scattering. With the
account of the absorption corrections our results are in agreement with the LHC
data on $\sigma_{tot}^{pp}$.

A Pomeron model applicable to both "`soft" and "hard" processes is suggested
and tested against the highenergy data from virtual photoninduced reactions.
The Pomeron is universal, containing two terms, a "soft" and a "hard" one,
whose relative weight varies with $\widetilde {Q^2}=Q^2+M_V^2$, where $Q^2$ is
the virtuality of the incoming photon and $M_V$ is the mass of the produced
vector particle. With a small number of adjustable parameters, the model fits
all available data on vector meson production and deeply virtual Compton
scattering from HERA. Furthermore, we attempt to apply the model to
hadroninduced reactions, by using highenergy data from protonproton
scattering.

A Reggeometric (Regge+Geometry) model, based on the observed proportionality
between the forward slope of the differential cross section and the interaction
radius, the latter depending on virtuality Q2 of the incoming virtual photon
and on the mass M 2 of the produced particle, is constructed. The objective of
this study is the dependence of the Reggepole amplitude on the virtuality Q2
and masses of the external particles, which remains an open problem for the
theory. The present analysis is based on the HERA data on Deeply Virtual
Compton Scattering (DVCS) and exclusive diffractive Vector Meson Production
(VMP). We treat each class of reactions separately, anticipating a further
study that will include both a soft and a hard component of the unique pomeron.

It is shown that in the nexttoleading approximation of N=4 SUSY the BFKL
equation for twogluon composite states in the adjoint representation of the
gauge group can be reduced to a form which is invariant under Moebius
transformation in the momentum space. The corresponding similarity
transformation of its integral kernel is constructed in an explicit way.

We investigate, both analytically and numerically, the phase diagram of
threedimensional Z(N) lattice gauge theories at finite temperature for N > 4.
These models, in the strong coupling limit, are equivalent to a generalized
version of vector Potts models in two dimension, with Polyakov loops playing
the role of Z(N) spins. It is argued that the effective spin models have two
phase transitions of infinite order (i.e. BKT). Using a cluster algorithm we
confirm this conjecture, locate the position of the critical points and extract
various critical indices.

We perform an analytical and numerical study of the phase transitions in
threedimensional Z(N) lattice gauge theories at finite temperature for N>4
exploiting equivalence of these models with a generalized version of the
twodimensional vector Potts models in the limit of vanishing spatial coupling.
In this limit the Polyakov loops play the role of Z(N) spins. The effective
couplings of these twodimensional spin models are calculated explicitly. It is
argued that the effective spin models have two phase transitions of BKT type.
This is confirmed by largescale Monte Carlo simulations. Using a cluster
algorithm we locate the position of the critical points and study the critical
behavior across both phase transitions in details. In particular, we determine
various critical indices, compute the helicity modulus, the average action and
the specific heat. A scaling formula for the critical points with N is
proposed.

A feasible mechanism of unitarization of amplitudes of deep inelastic
scattering at small values of Bjorken $x$ is the gluon fusion. However, its
efficiency depends crucially on the vacuum color screening effect which
accompanies the multiplication and the diffusion of BFKL gluons from small to
large distances. From the fits to lattice data on field strength correlators
the propagation length of perturbative gluons is $R_c\simeq 0.20.3$ fermi. The
probability to find a perturbative gluon with short propagation length at large
distances is suppressed exponentially. It changes the pattern of (dif)fusion
dramatically. The magnitude of the fusion effect appears to be controlled by
the new dimensionless parameter $\sim R_c^2/8B$, with the diffraction cone
slope $B$ standing for the characteristic size of the interaction region. It
should slowly $\propto 1/\ln Q^2$ decrease at large $Q^2$. Smallness of the
ratio $R_c^2/8B$ makes the nonlinear effects rather weak even at lowest
Bjorken $x$ available at HERA. We report the results of our studies of the
nonlinear BFKL equation which has been generalized to incorporate the running
coupling and the screening radius $R_c$ as the infrared regulator.

The perturbative nonlinear (NL) effects in the small$x$ evolution of the
gluon densities depend crucially on the infrared (IR) regularization. The IR
regulator, $R_c$, is determined by the scale of the nonperturbative
fluctuations of QCD vacuum. From the instanton models and from the lattice
$R_c\sim 0.3$ fm. For perturbative gluons with the propagation length $R_c=
0.26$ fm the linear BFKL gives a good description of the proton structure
function $F_2(x,Q^2)$ in a wide range of $x$ and $Q^2$. The NL effects turn out
to be rather weak and amount to the 10% correction to $F_2(x,Q^2)$ for $x\lsim
10^{5}$. Much more pronounced NL effects were found in the nonlinear model,
described in the literature, with a very soft IR regularization corresponding
to the IR cutoff at $\simeq \Lambda^{1}_{QCD}$. The latter issue is also
commented below.

As it was recently shown, the colour singlet BFKL kernel, taken in Moebius
representation in the space of impact parameters, can be written in
quasiconformal shape, which is unbelievably simple compared with the
conventional form of the BFKL kernel in momentum space. It was also proved that
the total kernel is completely defined by its Moebius representation. In this
paper we calculated the difference between standard and quasiconformal BFKL
kernels in momentum space and discovered that it is rather simple. Therefore we
come to the conclusion that the simplicity of the quasiconformal kernel is
caused mainly by using the impact parameter space.

We investigate both analytically and numerically the renormalization group
equations in 2D Z(N) vector models. The position of the critical points of the
two phase transitions for N>4 is established and the critical index \nu\ is
computed. For N=7, 17 the critical points are located by Monte Carlo
simulations and some of the corresponding critical indices are determined. The
behavior of the helicity modulus is studied for N=5, 7, 17. Using these and
other available Monte Carlo data we discuss the scaling of the critical points
with N and some other open theoretical problems.

We discuss heavy quark contributions to the neutrinonucleon total cross
section at very high energies, well above the real top production threshold.
The topbottom weak current is found to generate strong leftright asymmetry of
neutrinonucleon interactions. We separate contributions of different helicity
states and make use of the $\bkappa$factorization to derive simple and
practically useful formulas for the lefthanded ($F_L$) and righthanded
($F_R$) components of the conventional structure function $2xF_3=F_LF_R$ in
terms of the integrated gluon density. We show that $F_L\gg F_R$ and,
consequently, $xF_3\approx F_T$, where $F_T$ is the transverse structure
function. The conventional structure function $F_2=F_S+F_T$ at $Q^2\ll m_t^2$
appears to be dominated by its scalar (also known as longitudinal) component
$F_S$ and the hierarchy $F_S\gg F_L\gg F_R$ arises naturally. We evaluate the
total neutrinonucleon cross section at ultrahigh energies within the color
dipole BFKL formalism.

We study the connection between complete representations of gauge invariant
operators and their Moebius representations acting in a limited space of
functions. The possibility to restore the complete representations from Moebius
forms in the coordinate space is proven and a method of restoration is worked
out. The operators for transition from the standard BFKL kernel to the
quasiconformal one are found both in Moebius and total representations.

Currently available estimates of the gluonfusion effect in ultrahigh energy
neutrinonucleon interactions as well as in DIS on protons suffer from
uncertainty in defining the scattering profile function $\Gamma(b)$. Indeed,
the area, $S$, in the impact parameter space populated with interacting gluons
varies by a factor of $4  5$ from one analysis to another. To get rid of
uncertainties we specify the dipolenucleon partialwave amplitude $\Gamma(b)$
which meets the restrictions imposed by both the total dipolenucleon cross
section and the small angle elastic scattering amplitude. The area $S$ becomes
a well defined quantity proportional to the diffraction cone slope. We solve
numerically the nonlinear color dipole BFKL equation and evaluate the UHE
neutrinonucleon total cross section. Our finding is that the saturation is a
rather weak effect, $\lsim 25%$, up to $E_{\nu}\sim 10^{12}$ GeV.

Exclusive J/\Psi photo and electroproduction is studied in the framework of
the analytic Smatrix theory. The differential and integrated elastic cross
sections are calculated using the Modified Dual Amplitude with Mandelstam
Analyticity (MDAMA) model. The model is applied to the description of the
available experimantal data and proves to be valid in a wide region of the
kinematical variables s, t, and Q^2. Our amplitude can be used also as a
universal background parametrization for the extraction of tiny resonance
signals.

It is shown that in the ultrahigh energy neutrino interactions the higher
twist corrections brought about by the nonconservation of the topbottom
current dramatically change the longitudinal structure function, $F_L$. To the
Double Leading Log Approximation simple and numerically accurate formulas for
$F_L$ and $\sigma^{\nu N}$ are derived.

The overall hardness scale of the ultrahigh energy neutrinonucleon
interactions is usually estimated as $Q^2\sim m_W^2$. The effect of
nonconservation of weak currents pushes this scale up to the top quark mass
squared and changes dynamics of the scattering process. The Double Leading Log
Approximation provides simple and numerically accurate formula for the
topbottom contribution to the total cross section $\sigma^{\nu N}$.
Corresponding correction to $\sigma^{\nu N}$ appears to be numerically large.
It is comparable with the leading contribution evaluated in the massless quark
approximation.

An eikonalized elastic protonproton and protonantiproton scattering
amplitude F(s,t), calculated from QCD as a finite sum of gluon ladders, is
compared with the existing experimental data on the total cross section and the
ratio \rho(s,0)=Re F(s,0)/Im F(s,0) of the real part to the imaginary part of
the forward amplitude. Predictions for the expected LHC energies are given.

We demonstrate that the ambiguity of the lowx evolution kernels in the
nexttoleading order (NLO) permits one to match the Mobius form of the BFKL
kernel and the kernel of the colour dipole model and to construct the Mobius
invariant NLO BFKL kernel in N=4 supersymmetric YangMills theory.

The nonconservation of charmedstrange current in the neutrino deep
inelastic scattering ($\nu$DIS) strongly affects the longitudinal structure
function, $F_L$, at small values of Bjorken $x$. The corresponding correction
to $F_L$ is a higher twist effect enhanced at small$x$ by the rapidly growing
gluon density factor. As a result, the component of $F_L$ induced by the
charmedstrange current prevails over the lightquark component and dominates
$F_L=F_L^{cs}+F_L^{ud}$ at $x\lsim 0.01$ and $Q^2\sim m_c^2$. The color dipole
analysis clarifies the physics behind the phenomenon and provides a
quantitative estimate of the effect.

Exclusive J/Psi electroproduction is studied in the framework of the analytic
Smatrix theory. The differential and integrated elastic cross sections are
calculated using the Modified Dual Amplitude with Mandelstam Analyticity
(MDAMA) model. The model is applied to the description of the available
experimantal data and proves to be valid in a wide region of the kinematical
variables s, t and Q^2. Our amplitude can be used also as a universal
background parametrization for the extraction of tiny resonance signals.