Recently the CERN ALICE experiment, in its dedicated cosmic ray run, observed muon bundles of very high multiplicities, thereby confirming similar findings from the LEP era at CERN (in the CosmoLEP project). Significant evidence for anisotropy of arrival directions of the observed high multiplicity muonic bundles is found. Estimated directionality suggests their possible extragalactic provenance. We argue that muonic bundles of highest multiplicity are produced by strangelets, hypothetical stable lumps of strange quark matter infiltrating our Universe.

Recently the CERN ALICE experiment, in its dedicated cosmic ray run, observed muon bundles of very high multiplicities, thereby confirming similar findings from the LEP era at CERN (in the CosmoLEP project). Originally it was argued that they apparently stem from the primary cosmic rays with a heavy masses. We propose an alternative possibility arguing that muonic bundles of highest multiplicity are produced by strangelets, hypothetical stable lumps of strange quark matter infiltrating our Universe. We also address the possibility of additionally deducing their directionality which could be of astrophysical interest. Significant evidence for anisotropy of arrival directions of the observed high multiplicity muonic bundles is found. Estimated directionality suggests their possible extragalactic provenance.

Results on two-particle $\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.

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 string-resonance 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 event-by-event fluctuations of the fireball energy and/or volume.

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 long-baseline neutrino oscillation experiment in Japan. A comparison of the NA61/SHINE measurements with predictions of several hadroproduction models is presented.

Results from the NA49 experiment at the CERN SPS are presented on event-by-event 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.

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 nucleus-nucleus 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.

We evaluate the exposure during nadir observations with JEM-EUSO, the Extreme Universe Space Observatory, on-board the Japanese Experiment Module of the International Space Station. Designed as a mission to explore the extreme energy Universe from space, JEM-EUSO will monitor the Earth's nighttime atmosphere to record the ultraviolet light from tracks generated by extensive air showers initiated by ultra-high energy cosmic rays. In the present work, we discuss the particularities of space-based observation and we compute the annual exposure in nadir observation. The results are based on studies of the expected trigger aperture and observational duty cycle, as well as, on the investigations of the effects of clouds and different types of background light. We show that the annual exposure is about one order of magnitude higher than those of the presently operating ground-based observatories.

It has been recently noticed that transverse momenta ($p_T$) distributions observed in high energy production processes exhibit remarkably universal scaling behavior. This is seen when they are in some suitable variable, replacing the usual $p_T$. On the other hand, it is also known that transverse momenta distributions in general follow a power-like Tsallis distribution, rather than an exponential Boltzmann-Gibbs, with a (generally energy dependent) nonextensivity parameter $q$. We now show that it is possible to choose a suitable variable such that $p_T$ distributions of particles produced in proton-proton interactions in a wide energy range can be fitted by the same Tsallis distribution (with the same, energy independent, value of the $q$-parameter). Similar scaling behavior in nucleus-nucleus collisions is also observed. The possible dynamical origin of the q parameter used in these fits will be discussed.

The T2K long-baseline 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 2007-2010 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 semi-central 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 system-size and centrality dependence of relative strangeness production in nucleus-nucleus 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 core-corona approach.

Contributions of the JEM-EUSO Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.

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.

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.

Energy spectrum of cosmic rays (CR) exhibits power-like behavior with very characteristic "knee" structure. We consider a generalized statistical model for the production process of cosmic rays which accounts for such behavior in a natural way either by assuming existence of temperature fluctuations in the source of CR, or by assuming specific temperature distribution of CR sources. Both possibilities yield the so called Tsallis statistics and lead to the power-like distribution. We argue that the "knee" structure arises as result of abrupt change of fluctuations in the source of CR. Its possible origin is briefly discussed.

Examples of joint probability distributions are studied in terms of Tsallis' nonextensive statistics both for correlated and uncorrelated variables, in particular it is explicitely shown how correlations in the system can make Tsallis entropy additive and that the effective nonextensivity parameter $q_N$ decreases towards unity when the number of variables $N$ increases. We demonstrate that Tsallis distribution of energies of particles in a system leads in natural way to the Negative Binomial multiplicity distribution in this system.

Stochastic networks represent very important subject of research because they have been found in almost all branches of modern science, including also sociology and economy. We provide a information theory point of view, mostly based on its nonextensive version, on their most characteristic properties illustrating it with some examples.

Recently new data from the Cosmo-LEP project appeared, this time from DELPHI detector. They essentially confirm the findings reported some time ago by ALEPH, namely the appearence of bundles of muons with unexpectedly high multiplicities, which so far cannot be accounted by present day models. We argue, using arguments presented by us some time ago, that this phenomenon could be regarded as one more candidate for the presence in the flux of cosmic rays entering the Earth's atmosphere from outer space nuggets od Strange Quark Matter (SQM) in form of so called strangelets.

Assuming that cosmic rays entering the Earth's atmosphere contain a small admixture of nuggets of strange quark matter in form of strangelets one can explain a number of apparently "strange" effects observed in different cosmic rays experiments. We shall demonstrate here that the mass spectrum of such strangelets filles the "nuclear desert" gap existing between the heaviest elements observed in Universe and the next "nuclear-like objects" represented by neutron and strange stars.

We discuss recently measured event-by-event fluctuations of transverse momentum and of multiplicity in relativistic heavy-ion collisions. It is shown that the non-monotonic behavior of the p_T-fluctuations as a function of collision centrality can be fully explained by the observed non-monotonic multiplicity fluctuations. A possible mechanism responsible for the multiplicity fluctuations is also considered.

We discuss recently measured event-by-event multiplicity fluctuations in relativistic heavy-ion collisions. It is shown that the non-monotonic behavior of the multiplicity fluctuations as a function of collision centrality can be fully explained by the correlations between produced particles.

We demonstrate that hadron production viewed as formation of specific stochastic network can explain in natural way the power-law distributions of transverse mass spectra of pions found recently, which seem to substitute the expected Boltzmann statistical factor.

We describe the CASTOR detector designed to probe the very forward, baryon-rich rapidity region in nucleus-nucleus collisions at the LHC. We present a phenomenological model describing the formation of a QGP fireball in a high baryochemical potential environment, and its subsequent decay into baryons and strangelets. The model explains Centauros and the long-penetrating component and makes predictions for the LHC. Simulations of Centauro-type events were done. To study the response of the apparatus to new effects different exotic species (DCC, Centauros, strangelets etc.) were passed through the deep calorimeter. The energy deposition pattern in the calorimeter appears to be a new clear signature of the QGP.

The event-by-event p_T-fluctuations due to the temperature variation are considered. The so-called Phi-measure 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 proton-proton and central Pb-Pb collisions but independent measurements of the temperature fluctuations are needed to confirm the explanation. Feasibility of such an event-by-event measurement is discussed.