
Quarkgluon plasma produced at the early stage of ultrarelativistic heavy ion
collisions is unstable, if weakly coupled, due to the anisotropy of its
momentum distribution. Chromomagnetic fields are spontaneously generated and
can reach magnitudes much exceeding typical values of the fields in
equilibrated plasma. We consider a high energy test parton traversing an
unstable plasma that is populated with strong fields. We study the momentum
broadening parameter $\hat q$ which determines the radiative energy loss of the
test parton. We develop a formalism which gives $\hat q$ as the solution of an
initial value problem, and we focus on extremely oblate plasmas which are
physically relevant for relativistic heavy ion collisions. The parameter $\hat
q$ is found to be strongly dependent on time. For short times it is of the
order of the equilibrium value, but at later times $\hat q$ grows exponentially
due to the interaction of the test parton with unstable modes and becomes much
bigger than the value in equilibrium. The momentum broadening is also strongly
directionally dependent and is largest when the test parton velocity is
transverse to the beam axis. Consequences of our findings for the phenomenology
of jet quenching in relativistic heavy ion collisions are briefly discussed.

Production of $d$, $t$, and $^3$He nuclei in central Pb+Pb interactions was
studied at five collision energies ($\sqrt{s_{NN}}=$ 6.3, 7.6, 8.8, 12.3, and
17.3 GeV) with the NA49 detector at the CERN SPS. Transverse momentum spectra,
rapidity distributions, and particle ratios were measured. Yields are compared
to predictions of statistical models. Phasespace distributions of light nuclei
are discussed and compared to those of protons in the context of a coalescence
approach. The coalescence parameters $B_2$ and $B_3$, as well as coalescence
radii for $d$ and $^3$He were determined as a function of transverse mass at
all energies.

We look for fluctuations expected for the QCD critical point using an
intermittency analysis in the transverse momentum phase space of protons
produced around midrapidity in the 12.5\% most central C+C, Si+Si and Pb+Pb
collisions at the maximum SPS energy of 158$A$~GeV. We find evidence of
powerlaw fluctuations for the Si+Si data. The fitted powerlaw exponent
$\phi_{2} = 0.96^{+0.38}_{0.25}\text{ (stat.)}$ $\pm 0.16\text{ (syst.)}$ is
consistent with the va\lue expected for critical fluctuations. Powerlaw
fluctuations had previously also been observed in lowmass $\pi^+ \pi^$ pairs
in the same Si+Si collisions.

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.

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.

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.

Measurements of charged pion and kaon production are presented in centrality
selected Pb+Pb collisions at 40A GeV and 158A GeV beam energy as well as in
semicentral C+C and Si+Si interactions at 40A GeV. Transverse mass spectra,
rapidity spectra and total yields are determined as a function of centrality.
The systemsize and centrality dependence of relative strangeness production in
nucleusnucleus collisions at 40A GeV and 158A GeV beam energy are derived from
the data presented here and published data for C+C and Si+Si collisions at 158A
GeV beam energy. At both energies a steep increase with centrality is observed
for small systems followed by a weak rise or even saturation for higher
centralities. This behavior is compared to calculations using transport models
(UrQMD and HSD), a percolation model and the corecorona approach.

The energy loss of a fast parton scattering elastically in a weakly coupled
quarkgluon plasma is formulated as an initial value problem. The approach is
designed to study an unstable plasma, but it also reproduces the well known
result of energy loss in an equilibrium plasma. A twostream system, which is
unstable due to longitudinal chromoelectric modes, is discussed here some
detail. In particular, a strong time and directional dependence of the energy
loss is demonstrated.

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.

Interaction cross sections and charged pion spectra in p+C interactions at 31
GeV/c were measured with the large acceptance NA61/SHINE spectrometer at the
CERN SPS. These data are required to improve predictions of the neutrino flux
for the T2K long baseline neutrino oscillation experiment in Japan. A set of
data collected during the first NA61/SHINE run in 2007 with an isotropic
graphite target with a thickness of 4% of a nuclear interaction length was used
for the analysis. The measured p+C inelastic and production cross sections are
257.2 + 1.9 + 8.9 mb and 229.3 + 1.9 + 9.0 mb, respectively. Inclusive
production cross sections for negatively and positively charged pions are
presented as a function of laboratory momentum in 10 intervals of the
laboratory polar angle covering the range from 0 up to 420 mrad. The spectra
are compared with predictions of several hadron production models.

Kaons and protons carry large parts of two conserved quantities, strangeness
and baryon number. It is argued that their correlation and thus also
fluctuations are sensitive to conditions prevailing at the anticipated
partonhadron phase boundary. Fluctuations of the $(\mathrm{K}^+ +
\mathrm{K}^)/(\mathrm{p}+\bar{\mathrm{p}})$ and $\mathrm{K}^+/\mathrm{p}$
ratios have been measured for the first time by NA49 in central Pb+Pb
collisions at 5 SPS energies between $\sqrt{s_{NN}}$= 6.3 GeV and 17.3 GeV.
Both ratios exhibit a change of sign in $\sigma_{\mathrm{dyn}}$, a measure of
nonstatistical fluctuations, around $\sqrt{s_{NN}}$ = 8 GeV. Below this
energy, $\sigma_{\mathrm{dyn}}$ is positive, indicating higher fluctuation
compared to a mixed event background sample, while for higher energies,
$\sigma_{\mathrm{dyn}}$ is negative, indicating correlated emission of kaons
and protons. The results are compared to UrQMD calculations which which give a
good description at the higher SPS energies, but fail to reproduce the
transition to positive values.

The momentum correlation between protons and lambda particles emitted from
central Pb+Pb collisions at sqrt(s_{NN}) = 17.3 GeV was studied by the NA49
experiment at the CERN SPS. A clear enhancement is observed for small relative
momenta (q_{inv} < 0.2 GeV). By fitting a theoretical model, which uses the
strong interaction between the proton and the lambda in a given pair, to the
measured data a value for the effective source size is deduced. Assuming a
static Gaussian source distribution we derive an effective radius parameter of
R_G = 3.02 \pm 0.20$(stat.)^{+0.44}_{0.16}(syst.) fm.

Eventbyevent fluctuations of the chemical composition of the hadronic final
state of relativistic heavyion collisions carry valuable information on the
properties of strongly interacting matter produced in the collisions. However,
in experiments incomplete particle identification distorts the observed
fluctuation signals. The effect is quantitatively studied and a new technique
for measuring chemical fluctuations, the identity method, is proposed. The
method fully eliminates the effect of incomplete particle identification. The
application of the identity method to experimental data is explained.

The average transverse momentum transfer per unit path length to a fast
parton scattering elastically in a perturbative quarkgluon plasma is related
to the radiative energy loss of the parton. We first calculate the momentum
transfer coefficient $\hat q$ in terms of a classical Langevin problem and then
define it quantummechanically through scattering matrix element. After
treating the well known case of a quarkgluon plasma in equilibrium we consider
an offequilibrium unstable plasma. As a specific example, we treat the
twostream plasma with unstable modes of longitudinal chromoelectric field. In
the presence of the instabilities, $\hat q$ is shown to exponentially grow in
time.

Phi meson production is studied by the NA49 Collaboration in central Pb+Pb
collisions at 20A, 30A, 40A, 80A and 158A GeV beam energy. The data are
compared with measurements at lower and higher energies and to microscopic and
thermal models. The energy dependence of yields and spectral distributions is
compatible with the assumption that partonic degrees of freedom set in at low
SPS energies.

The latest NA49 results on eventbyevent transverse momentum fluctuations
are presented for central Pb+Pb interactions over the whole SPS energy range
(20A  158A GeV). Two different methods are applied: evaluating the
$\Phi_{p_{T}}$ fluctuation measure and studying twoparticle transverse
momentum correlations. The obtained results are compared to predictions of the
UrQMD model. The results on the energy dependence are compared to the NA49 data
on the system size dependence. The NA61 (SHINE, NA49future) strategy of
searching of the QCD critical endpoint is also discussed.

We propose a method to measure the subdiffusion parameter $\alpha$ and
subdiffusion coefficient $D_{\alpha}$ which are defined by means of the
relation $<x^2> =\frac{2D_\alpha} {\Gamma(1+\alpha)} t^\alpha$ where $<x^2>$
denotes a mean square displacement of a random walker starting from $x=0$ at
the initial time $t=0$. The method exploits a membrane system where a substance
of interest is transported in a solvent from one vessel to another across a
thin membrane which plays here only an auxiliary role. We experimentally study
a diffusion of glucose and sucrose in a gel solvent, and we precisely determine
the parameters $\alpha$ and $D_{\alpha}$, using a fully analytic solution of
the fractional subdiffusion equation.

We propose a method to extract from experimental data the subdiffusion
parameter $\alpha$ and subdiffusion coefficient $D_\alpha$ which are defined by
means of the relation $<x^2> =2D_\alpha/\Gamma(1+\alpha) t^\alpha$ where
$<x^2>$ denotes a mean square displacement of a random walker starting from
$x=0$ at the initial time $t=0$. The method exploits a membrane system where a
substance of interest is transported in a solvent from one vessel to another
across a thin membrane which plays here only an auxiliary role. Using such a
system, we experimentally study a diffusion of glucose and sucrose in a gel
solvent. We find a fully analytic solution of the fractional subdiffusion
equation with the initial and boundary conditions representing the system under
study. Confronting the experimental data with the derived formulas, we show a
subdiffusive character of the sugar transport in gel solvent. We precisely
determine the parameter $\alpha$, which is smaller than 1, and the subdiffusion
coefficient $D_\alpha$.

We discuss recently measured eventbyevent fluctuations of transverse
momentum and of multiplicity in relativistic heavyion collisions. It is shown
that the nonmonotonic behavior of the p_Tfluctuations as a function of
collision centrality can be fully explained by the observed nonmonotonic
multiplicity fluctuations. A possible mechanism responsible for the
multiplicity fluctuations is also considered.

The eventbyevent p_Tfluctuations due to the temperature variation are
considered. The socalled Phimeasure is computed and shown to be a linear
function of the temperature variance. The fluctuations of temperature naturally
explain the data on Phi(p_T) in protonproton and central PbPb collisions but
independent measurements of the temperature fluctuations are needed to confirm
the explanation. Feasibility of such an eventbyevent measurement is
discussed.