
Measurements of inclusive spectra and mean multiplicities of $\pi^\pm$,
K$^\pm$, p and $\bar{\textrm{p}}$ produced in inelastic p+p interactions at
incident projectile momenta of 20, 31, 40, 80 and 158 GeV/c ($\sqrt{s} = $ 6.3,
7.7, 8.8, 12.3 and 17.3 GeV, respectively) were performed at the CERN Super
Proton Synchrotron using the large acceptance NA61/SHINE hadron spectrometer.
Spectra are presented as function of rapidity and transverse momentum and are
compared to predictions of current models. The measurements serve as the
baseline in the NA61/SHINE study of the properties of the onset of
deconfinement and search for the critical point of strongly interacting matter.

We present measurements of $\rho^0$, $\omega$ and K$^{*0}$ spectra in
$\pi^{} + $C production interactions at 158 GeV/c and $\rho^0$ spectra at 350
GeV/c using the NA61/SHINE spectrometer at the CERN SPS. Spectra are presented
as a function of the Feynman's variable $x_\text{F}$ in the range $0 <
x_\text{F} < 1$ and $0 < x_\text{F} < 0.5$ for 158 GeV/c and 350 GeV/c
respectively. Furthermore, we show comparisons with previous measurements and
predictions of several hadronic interaction models. These measurements are
essential for a better understanding of hadronic shower development and for
improving the modeling of cosmic ray air showers.

A unified equation of state for quarkhadron matter is presented in the
generalized BethUhlenbeck form. It follows from a $\Phi$derivable approach to
the thermodynamic potential where the ansatz for the $\Phi$ functional contains
all 2PI diagrams at twoloop order formed with quark cluster Green's functions
for quark, diquark, meson and baryon propagators. We present numerical results
using an effective model for the generic behaviour of hadron masses and phase
shifts at finite temperature which shares basic features with recent
developments within the PNJL model for correlations in quark matter. We obtain
the transition between a hadron resonance gas phase and the quark gluon plasma
where the Mott dissociation of hadrons is encoded in the hadronic phase shifts.
The resulting thermodynamics is in very good agreement with recent lattice QCD
simulations.

Results on twoparticle $\Delta\eta\Delta\phi$ correlations in inelastic p+p
interactions at 20, 31, 40, 80, and 158~GeV/c are presented. The measurements
were performed using the large acceptance NA61/SHINE hadron spectrometer at the
CERN Super Proton Synchrotron. The data show structures which can be attributed
mainly to effects of resonance decays, momentum conservation, and quantum
statistics. The results are compared with the EPOS and UrQMD models.

A unified equation of state for quarkhadron matter is presented in the
generalized BethUhlenbeck form. It follows from a $\Phi$derivable approach to
the thermodynamic potential where the ansatz for the $\Phi$ functional contains
all 2PI diagrams at twoloop order formed with quark cluster Green's functions
for quark, diquark, meson and baryon propagators. We present numerical results
using an effective model for the generic behaviour of hadron masses and phase
shifts at finite temperature which shares basic features with recent
developments within the PNJL model for correlations in quark matter. We obtain
the transition between a hadron resonance gas phase and the quark gluon plasma
where the Mott dissociation of hadrons is encoded in the hadronic phase shifts.
The resulting thermodynamics is in very good agreement with recent lattice QCD
simulations.

Measurements of particle emission from a replica of the T2K 90 cmlong carbon
target were performed in the NA61/SHINE experiment at CERN SPS, using data
collected during a highstatistics run in 2009. An efficient use of the
longtarget measurements for neutrino flux predictions in T2K requires
dedicated reconstruction and analysis techniques. Fullycorrected differential
yields of $\pi^\pm$mesons from the surface of the T2K replica target for
incoming 31 GeV/c protons are presented. A possible strategy to implement these
results into the T2K neutrino beam predictions is discussed and the propagation
of the uncertainties of these results to the final neutrino flux is performed.

Measurements of multiplicity and transverse momentum fluctuations of charged
particles were performed in inelastic p+p interactions at 20, 31, 40, 80 and
158 GeV/c beam momentum. Results for the scaled variance of the multiplicity
distribution and for three strongly intensive measures of multiplicity and
transverse momentum fluctuations \$\Delta[P_{T},N]\$, \$\Sigma[P_{T},N]\$ and
\$\Phi_{p_T}\$ are presented. For the first time the results on fluctuations
are fully corrected for experimental biases. The results on multiplicity and
transverse momentum fluctuations significantly deviate from expectations for
the independent particle production. They also depend on charges of selected
hadrons. The stringresonance Monte Carlo models EPOS and UrQMD do not describe
the data. The scaled variance of multiplicity fluctuations is significantly
higher in inelastic p+p interactions than in central Pb+Pb collisions measured
by NA49 at the same energy per nucleon. This is in qualitative disagreement
with the predictions of the Wounded Nucleon Model. Within the statistical
framework the enhanced multiplicity fluctuations in inelastic p+p interactions
can be interpreted as due to eventbyevent fluctuations of the fireball energy
and/or volume.

Inclusive production of $\Lambda$hyperons was measured with the large
acceptance NA61/SHINE spectrometer at the CERN SPS in inelastic p+p
interactions at beam momentum of 158~\GeVc. Spectra of transverse momentum and
transverse mass as well as distributions of rapidity and x$_{_F}$ are
presented. The mean multiplicity was estimated to be
$0.120\,\pm0.006\;(stat.)\,\pm 0.010\;(sys.)$. The results are compared with
previous measurements and predictions of the EPOS, UrQMD and FRITIOF models.

Measurements of hadron production in p+C interactions at 31 GeV/c are
performed using the NA61/ SHINE spectrometer at the CERN SPS. The analysis is
based on the full set of data collected in 2009 using a graphite target with a
thickness of 4% of a nuclear interaction length. Inelastic and production cross
sections as well as spectra of $\pi^\pm$, $K^\pm$, p, $K^0_S$ and $\Lambda$ are
measured with high precision. These measurements are essential for improved
calculations of the initial neutrino fluxes in the T2K longbaseline neutrino
oscillation experiment in Japan. A comparison of the NA61/SHINE measurements
with predictions of several hadroproduction models is presented.

The NA61/SHINE experiment performs a detailed study of the onset of
deconfinement and search for critical point of hadronic matter by colliding
nuclei of different sizes at various beam momenta from 13A to 158A GeV/c.
Experimental setup and results on the theoretically expected signatures are
discussed.

We present an effective model for the generic behaviour of hadron masses and
phase shifts at finite temperature which shares basic features with recent
developments within the PNJL model for correlations in quark matter. On this
basis we obtain the transition between a hadron resonance gas phase and the
quark gluon plasma in the spirit of the generalized BethUhlenbeck approach
where the Mott dissociation of hadrons is encoded in the hadronic phase shifts.
Here we restrict ourselves to lowlying hadronic channels and perform a
discussion of recent lattice QCD thermodynamics results from this perspective.
We find agreement in the asymptotic regions while for the description of the
transition itself the inclusion of further hadronic channels as well as a
selfconsistent determination of the continuum thresholds is required.

We present experimental results on inclusive spectra and mean multiplicities
of negatively charged pions produced in inelastic p+p interactions at incident
projectile momenta of 20, 31, 40, 80 and 158 GeV/c ($\sqrt{s} = $ 6.3, 7.7,
8.8, 12.3 and 17.3 GeV, respectively). The measurements were performed using
the large acceptance NA61/SHINE hadron spectrometer at the CERN Super Proton
Synchrotron.
Twodimensional spectra are determined in terms of rapidity and transverse
momentum. Their properties such as the width of rapidity distributions and the
inverse slope parameter of transverse mass spectra are extracted and their
collision energy dependences are presented. The results on inelastic p+p
interactions are compared with the corresponding data on central Pb+Pb
collisions measured by the NA49 experiment at the CERN SPS.
The results presented in this paper are part of the NA61/SHINE ion program
devoted to the study of the properties of the onset of deconfinement and search
for the critical point of strongly interacting matter. They are required for
interpretation of results on nucleusnucleus and protonnucleus collisions.

An effective model reproducing the equation of state of hadronic matter as
obtained in recent lattice QCD simulations and from hadron resonance gas data
is presented. The hadronic phase is described by means of an extended
MottHagedorn resonance gas while the QGP phase is described by the extended
PNJL model. The dissociation of hadrons is obtained by including the state
dependent hadron resonance width. The model gives a quantitative estimate for
partial fractions of hadronic and partonic degrees of freedom above $T_c$.

Spectra of K0S mesons and Lambda hyperons were measured in p+C interactions
at 31 GeV/c with the large acceptance NA61/SHINE spectrometer at the CERN SPS.
The data were collected with an isotropic graphite target with a thickness of
4% of a nuclear interaction length. Interaction cross sections, charged pion
spectra, and charged kaon spectra were previously measured using the same data
set. Results on K0S and Lambda production in p+C interactions serve as
reference for the understanding of the enhancement of strangeness production in
nucleusnucleus collisions. Moreover, they provide important input for the
improvement of neutrino flux predictions for the T2K long baseline neutrino
oscillation experiment in Japan. Inclusive production cross sections for K0S
and Lambda are presented as a function of laboratory momentum in intervals of
the laboratory polar angle covering the range from 0 up to 240 mrad. The
results are compared with predictions of several hadron production models. The
K0S mean multiplicity in production processes <n_K0S> and the inclusive cross
section for K0S production were measured and amount to 0.127 + 0.005 (stat) +
0.022 (sys) and 29.0 + 1.6 (stat) + 5.0 (sys) mb, respectively.

A combined effective model reproducing the equation of state of hadronic
matter as obtained in recent lattice QCD simulations is presented. The model
reproduces basic physical characteristics encountered in dense hadronic matter
in the quarkgluon plasma (QGP) phase and the lower temperature hadron
resonance gas phase. The hadronic phase is described by means of an extended
MottHagedorn resonance gas while the QGP phase is described by the extended
PNJL model. The dissociation of hadrons is obtained by including the state
dependent hadron resonance width.

The T2K longbaseline neutrino oscillation experiment in Japan needs precise
predictions of the initial neutrino flux. The highest precision can be reached
based on detailed measurements of hadron emission from the same target as used
by T2K exposed to a proton beam of the same kinetic energy of 30 GeV. The
corresponding data were recorded in 20072010 by the NA61/SHINE experiment at
the CERN SPS using a replica of the T2K graphite target. In this paper details
of the experiment, data taking, data analysis method and results from the 2007
pilot run are presented. Furthermore, the application of the NA61/SHINE
measurements to the predictions of the T2K initial neutrino flux is described
and discussed.

A compound effective model reproducing the equation of state of hadronic
matter as obtained in recent lattice QCD simulations is presented.

Spectra of positively charged kaons in p+C interactions at 31 GeV/c were
measured with the NA61/SHINE spectrometer at the CERN SPS. The analysis is
based on the full set of data collected in 2007 with a graphite target with a
thickness of 4% of a nuclear interaction length. Interaction cross sections and
charged pion spectra were already measured using the same set of data. These
new measurements in combination with the published ones are required to improve
predictions of the neutrino flux for the T2K long baseline neutrino oscillation
experiment in Japan. In particular, the knowledge of kaon production is crucial
for precisely predicting the intrinsic electron neutrino component and the high
energy tail of the T2K beam. The results are presented as a function of
laboratory momentum in 2 intervals of the laboratory polar angle covering the
range from 20 up to 240 mrad. The kaon spectra are compared with predictions of
several hadron production models. Using the published pion results and the new
kaon data, the K+/\pi+ ratios are computed.

A kinetic equation approach is applied to model anomalous J/psi suppression
at RHIC and SPS by absorption in a hadron resonance gas which successfully
describes statistical hadron production in both experiments. The puzzling
rapidity dependence of the PHENIX data is reproduced as a geometric effect due
to a longer absorption path for J/psi production at forward rapidity.

The property of extended longitudinal scaling of rapidity distributions was
noticed recently over a broad range of beam energies. It is shown here that
this property is consistent with predictions of the statistical thermal model
up to the highest RHIC beam energies, however, we expect that at LHC energies
the rapidity distribution of produced particles will violate extended
longitudinal scaling.

The relativistic hydrodynamic approach is used to describe production of
strangeness and/or heavy quarks in ultrarelativistic heavy ion reactions.
Production processes are important ingredients of dissipative effects in the
hadronic liquid. Beyond viscosity also chemo and thermodiffusion processes
are considered. This also allows to specify chemical and thermal freezeout
conditions.

The properties of statistical ensembles with abelian charges close to the
thermodynamic limit are discussed. The finite volume corrections to the
probability distributions and particle density moments are calculated. Results
are obtained for statistical ensembles with both exact and average charge
conservation. A new class of variables (semiintensive variables) which differ
in the thermodynamic limit depending on how charge conservation is implemented
in the system is introduced. The thermodynamic limit behavior of these
variables is calculated through the next to leading order finite volume
corrections to the corresponding probability density distributions.

The probability distributions for charged particle numbers and their
densities are derived in statistical ensembles with conservation laws. It is
shown that if this limit is properly taken then the canonical and grand
canonical ensembles are equivalent. This equivalence is proven on the most
general, probability distribution level.

The entropy production occurring in relativistic hydrodynamical systems such
as the quarkgluon plasma (QGP) formed in highenergy nuclear collisions is
explored. We study mechanisms which change the composition of the fluid, i.e.
particle production and/or chemical reactions, along with chemo and
thermodiffusion. These effects complement the conventional dissipative effects
of shear viscosity, bulk viscosity, and heat conductivity.

Using the monomerdimer representation of the lattice Schwinger model, with
$N_f =1$ Wilson fermions in the strongcoupling regime ($\beta=0$), we
evaluate its partition function, $Z$, exactly on finite lattices. By studying
the zeroes of $Z(k)$ in the complex plane $(Re(k),Im(k))$ for a large number of
small lattices, we find the zeroes closest to the real axis for infinite
stripes in temporal direction and spatial extent $S=2$ and 3. We find evidence
for the existence of a critical value for the hopping parameter in the
thermodynamic limit $S\rightarrow \infty$ on the real axis at about $k_c \simeq
0.39$. By looking at the behaviour of quantities, such as the chiral
condensate, the chiral susceptibility and the third derivative of $Z$ with
respect to $1/2k$, close to the critical point $k_c$, we find some indications
for a continuous phase transition.