
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.

Results from the NA49 experiment at the CERN SPS are presented on
eventbyevent transverse momentum and multiplicity fluctuations of charged
particles, produced at forward rapidities in central Pb+Pb interactions at beam
momenta 20$A$, 30$A$, 40$A$, 80$A$, and 158$A$ GeV/c, as well as in systems of
different size ($p+p$, C+C, Si+Si, and Pb+Pb) at 158$A$ GeV/c. This publication
extends the previous NA49 measurements of the strongly intensive measure
$\Phi_{p_T}$ by a study of the recently proposed strongly intensive measures of
fluctuations $\Delta[P_T, N]$ and $\Sigma[P_T, N]$. In the explored kinematic
region transverse momentum and multiplicity fluctuations show no significant
energy dependence in the SPS energy range. However, a remarkable system size
dependence is observed for both $\Delta[P_T, N]$ and $\Sigma[P_T, N]$, with the
largest values measured in peripheral Pb+Pb interactions. The results are
compared with NA61/SHINE measurements in $p+p$ collisions, as well as with
predictions of the UrQMD and EPOS models.

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.

A novel approach, the identity method, was used for particle identification
and the study of fluctuations of particle yield ratios in Pb+Pb collisions at
the CERN Super Proton Synchrotron (SPS). This procedure allows to unfold the
moments of the unknown multiplicity distributions of protons (p), kaons (K),
pions ($\pi$) and electrons (e). Using these moments the excitation function of
the fluctuation measure $\nu_{\text{\text{dyn}}}$[A,B] was measured, with A and
B denoting different particle types. The obtained energy dependence of
$\nu_{\text{dyn}}$ agrees with previously published NA49 results on the related
measure $\sigma_{\text{dyn}}$. Moreover, $\nu_{\text{dyn}}$ was found to depend
on the phase space coverage for [K,p] and [K,$\pi$] pairs. This feature most
likely explains the reported differences between measurements of NA49 and those
of STAR in central Au+Au 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.

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.

The recently introduced multifractal network generator (MFNG), has been shown
to provide a simple and flexible tool for creating random graphs with very
diverse features. The MFNG is based on multifractal measures embedded in 2d,
leading also to isolated nodes, whose number is relatively low for realistic
cases, but may become dominant in the limiting case of infinitely large network
sizes. Here we discuss the relation between this effect and the information
dimension for the 1d projection of the link probability measure (LPM), and
argue that the node isolation can be avoided by a simple transformation of the
LPM based on rotation.

The yields of (anti)protons were measured by the NA49 Collaboration in
centrality selected Pb+Pb collisions at 40A GeV and 158A GeV. Particle
identification was obtained in the laboratory momentum range from 5 to 63 GeV/c
by the measurement of the energy loss dE/dx in the TPC detector gas. The
corresponding rapidity coverage extends 1.6 units from midrapidity into the
forward hemisphere. Transverse mass spectra, the rapidity dependences of the
average transverse mass, and rapidity density distributions were studied as a
function of collision centrality. The values of the average transverse mass as
well as the midrapidity yields of protons when normalized to the number of
wounded nucleons show only modest centrality dependences. In contrast, the
shape of the rapidity distribution changes significantly with collision
centrality, especially at 40A GeV. The experimental results are compared to
calculations of the HSD and UrQMD transport models.

We introduce a new approach to constructing networks with realistic features.
Our method, in spite of its conceptual simplicity (it has only two parameters)
is capable of generating a wide variety of network types with prescribed
statistical properties, e.g., with degree or clustering coefficient
distributions of various, very different forms. In turn, these graphs can be
used to test hypotheses, or, as models of actual data. The method is based on a
mapping between suitably chosen singular measures defined on the unit square
and sparse infinite networks. Such a mapping has the great potential of
allowing for graph theoretical results for a variety of network topologies. The
main idea of our approach is to go to the infinite limit of the singular
measure and the size of the corresponding graph simultaneously. A very unique
feature of this construction is that the complexity of the generated network is
increasing with the size. We present analytic expressions derived from the
parameters of the  to be iterated initial generating measure for such major
characteristics of graphs as their degree, clustering coefficient and
assortativity coefficient distributions. The optimal parameters of the
generating measure are determined from a simple simulated annealing process.
Thus, the present work provides a tool for researchers from a variety of fields
(such as biology, computer science, biology, or complex systems) enabling them
to create a versatile model of their network data.

In the last 30 years it was found that many combinatorial systems undergo
phase transitions. One of the most important examples of these can be found
among the random ksatisfiability problems (often referred to as kSAT), asking
whether there exists an assignment of Boolean values satisfying a Boolean
formula composed of clauses with k random variables each. The random 3SAT
problem is reported to show various phase transitions at different critical
values of the ratio of the number of clauses to the number of variables. The
most famous of these occurs when the probability of finding a satisfiable
instance suddenly drops from 1 to 0. This transition is associated with a rise
in the hardness of the problem, but until now the correlation between any of
the proposed phase transitions and the hardness is not totally clear. In this
paper we will first show numerically that the number of solutions universally
follows a lognormal distribution, thereby explaining the puzzling question of
why the number of solutions is still exponential at the critical point.
Moreover we provide evidence that the hardness of the closely related problem
of counting the total number of solutions does not show any phase
transitionlike behavior. This raises the question of whether the probability
of finding a satisfiable instance is really an order parameter of a phase
transition or whether it is more likely to just show a simple sharp threshold
phenomenon. More generally, this paper aims at starting a discussion where a
simple sharp threshold phenomenon turns into a genuine phase transition.

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 study the energy level structure of twodimensional charged particles in
inhomogeneous magnetic fields. In particular, for magnetic antidots the
magnetic field is zero inside the dot and constant outside. Such a device can
be fabricated with presentday technology. We present detailed semiclassical
studies of such magnetic antidot systems and provide a comparison with exact
quantum calculations. In the semiclassical approach we apply the BerryTabor
formula for the density of states and the BorhSommerfeld quantization rules.
In both cases we found good agreement with the exact spectrum in the weak
magnetic field limit. The energy spectrum for a given missing flux quantum is
classified in six possible classes of orbits and summarized in a socalled
phase diagram. We also investigate the current flow patterns of different
quantum states and show the clear correspondence with classical trajectories.

In this article we give an in depth overview of the recent advances in the
field of equilibrium networks. After outlining this topic, we provide a novel
way of defining equilibrium graph (network) ensembles. We illustrate this
concept on the classical random graph model and then survey a large variety of
recently studied network models. Next, we analyze the structural properties of
the graphs in these ensembles in terms of both local and global
characteristics, such as degrees, degreedegree correlations, component sizes,
and spectral properties. We conclude with topological phase transitions and
show examples for both continuous and discontinuous transitions.

In pseudo integrable systems diffractive scattering caused by wedges and
impurities can be described within the framework of Geometric Theory of
Diffraction (GDT) in a way similar to the one used in the Periodic Orbit Theory
of Diffraction (POTD). We derive formulas expressing the reflection and
transition matrix elements for one and many diffractive points and apply it for
impurity and wedge diffraction. Diffraction can cause backscattering in
situations, where usual semiclassical backscattering is absent causing an
erodation of ideal conductance steps. The length of diffractive periodic orbits
and diffractive loops can be detected in the power spectrum of the reflection
matrix elements. The tail of the power spectrum shows $\sim 1/l^{1/2}$ decay
due to impurity scattering and $\sim 1/l^{3/2}$ decay due to wedge scattering.
We think this is a universal sign of the presence of diffractive scattering in
pseudo integrable waveguides.