TUS (Tracking Ultraviolet Set-up) is the world's first orbital detector of ultra-high-energy cosmic rays (UHECRs). It was launched into orbit on 28th April 2016 as a part of the scientific payload of the Lomonosov satellite. The main aim of the mission was to test the technique of measuring the ultraviolet fluorescence and Cherenkov radiation of extensive air showers generated by primary cosmic rays with energies above ~100 EeV in the Earth atmosphere from space. During its operation for 1.5 years, TUS registered almost 80,000 events with a few of them satisfying conditions anticipated for extensive air showers (EASs) initiated by UHECRs. Here we discuss an event registered on 3rd October 2016. The event was measured in perfect observation conditions as an ultraviolet track in the nocturnal atmosphere of the Earth, with the kinematics and the light curve similar to those expected from an EAS. A reconstruction of parameters of a primary particle gave the zenith angle around 44$^\circ$ but an extreme energy not compatible with the cosmic ray energy spectrum obtained with ground-based experiments. We discuss in details all conditions of registering the event, explain the reconstruction procedure and its limitations and comment on possible sources of the signal, both of anthropogenic and astrophysical origin. We believe this detection represents a significant milestone in the space-based observation of UHECRs because it proves the capability of an orbital telescope to detect light signals with the apparent motion and light shape similar to what are expected from EASs. This is important for the on-going development of the future missions KLYPVE-EUSO and POEMMA, aimed for studying UHECRs from space.

The SPHERE-2 balloon-borne detector designed for extensive air shower (EAS) observations using EAS optical Vavilov-Cherenkov radiation (``Cherenkov light''), reflected from the snow-covered surface of Lake Baikal is described. We briefly discuss the concept behind the reflected Cherenkov light method, characterize the conditions at the experimental site and overview the construction of the tethered balloon used to lift the SPHERE-2 telescope above the surface. This paper is mainly dedicated to a detailed technical description of the detector, including its optical system, sensitive elements, electronics, and data acquisition system (DAQ). The results of some laboratory and field tests of the optical system are presented.

TUS (Tracking Ultraviolet Set-up), the first orbital detector of extreme energy cosmic rays (EECRs), those with energies above 50 EeV, was launched into orbit on April 28, 2016, as a part of the Lomonosov satellite scientific payload. The main aim of the mission is to test a technique of registering fluorescent and Cherenkov radiation of extensive air showers generated by EECRs in the atmosphere with a space telescope. We present preliminary results of its operation in a mode dedicated to registering extensive air showers in the period from August 16, 2016, to November 4, 2016. No EECRs have been conclusively identified in the data yet, but the diversity of ultraviolet emission in the atmosphere was found to be unexpectedly rich. We discuss typical examples of data obtained with TUS and their possible origin. The data is important for obtaining more accurate estimates of the nocturnal ultraviolet glow of the atmosphere, necessary for successful development of more advanced orbital EECR detectors including those of the KLYPVE (K-EUSO) and JEM-EUSO missions.

The origin and nature of extreme energy cosmic rays (EECRs), which have energies above the 50 EeV, the Greisen-Zatsepin-Kuzmin (GZK) energy limit, is one of the most interesting and complicated problems in modern cosmic-ray physics. Existing ground-based detectors have helped to obtain remarkable results in studying cosmic rays before and after the GZK limit, but have also produced some contradictions in our understanding of cosmic ray mass composition. Moreover, each of these detectors covers only a part of the celestial sphere, which poses problems for studying the arrival directions of EECRs and identifying their sources. As a new generation of EECR space detectors, TUS (Tracking Ultraviolet Set-up), KLYPVE and JEM-EUSO, are intended to study the most energetic cosmic-ray particles, providing larger, uniform exposures of the entire celestial sphere. The TUS detector, launched on board the Lomonosov satellite on April 28, 2016, from Vostochny Cosmodrome in Russia, is the first of these. It employs a single-mirror optical system and a photomultiplier tube matrix as a photo-detector and will test the fluorescent method of measuring EECRs from space. Utilizing the Earth's atmosphere as a huge calorimeter, it is expected to detect EECRs with energies above 100 EeV. It will also be able to register slower atmospheric transient events: atmospheric fluorescence in electrical discharges of various types including precipitating electrons escaping the magnetosphere and from the radiation of meteors passing through the atmosphere. We describe the design of the TUS detector and present results of different ground-based tests and simulations.

TUS (Track Ultraviolet Setup) is the first space experiment aimed to check the possibility of registering extreme energy cosmic rays (EECRs) at E>50 EeV by measuring the fluorescence signal of extensive air showers in the atmosphere. The detector operates as a part of the scientific payload of the Lomonosov satellite for more than a year. We describe an algorithm of searching for EECR events in the TUS data and briefly discuss a number of candidates selected by formal criteria.

On April 28, 2016, the Lomonosov satellite, equipped with a number of scientific instruments, was launched into orbit. Here we present briefly some of the results obtained with the first orbital telescope of extreme energy cosmic rays TUS and by a group of detectors aimed at multi-messenger observations of gamma-ray bursts.

We report theoretical and experimental results of on-going feasibility studies to detect cosmic neutrinos acoustically in Lake Baikal. In order to examine ambient noise conditions and to develop respective pulse detection techniques a prototype device was created. The device is operating at a depth of 150 m at the site of the Baikal Neutrino Telescope and is capable to detect and classify acoustic signals with different shapes, as well as signals from neutrino-induced showers.

Pairwise forces between particles in cosmological N-body simulations are generally softened to avoid hard collisions. Physically, this softening corresponds to treating the particles as diffuse clouds rather than point masses. For particles of unequal mass (and hence unequal softening length), computing the softened force involves a nontrivial double integral over the volumes of the two particles. We show that Plummer force softening is consistent with this interpretation of softening while spline softening is not. We provide closed-form expressions and numerical implementation for pairwise gravitational force laws for pairs of particles of general softening scales $\epsilon_1$ and $\epsilon_2$ assuming the commonly used cloud profiles: NGP, CIC, TSC, and PQS. Similarly, we generalize Plummer force law into pairs of particles of general softenings. We relate our expressions to the gaussian, Plummer and spline force softenings known from literature. Our expressions allow possible inclusions of pointlike particles such as stars or supermassive black holes.

We study the effect of primordial nongaussianity on large-scale structure, focusing upon the most massive virialized objects. Using analytic arguments and N-body simulations, we calculate the mass function and clustering of dark matter halos across a range of redshifts and levels of nongaussianity. We propose a simple fitting function for the mass function valid across the entire range of our simulations. We find pronounced effects of nongaussianity on the clustering of dark matter halos, leading to strongly scale-dependent bias. This suggests that the large-scale clustering of rare objects may provide a sensitive probe of primordial nongaussianity. We very roughly estimate that upcoming surveys can constrain nongaussianity at the level |fNL| <~ 10, competitive with forecasted constraints from the microwave background.

We investigate the origin of halo assembly bias, the dependence of halo clustering on assembly history. We relate halo assembly to peak properties measured in the Lagrangian space of the initial linear Gaussian random density field, and show how these same Lagrangian properties determine large-scale bias. We focus on the two regimes where assembly bias has been observed to be significant: at masses very large and very small compared to the nonlinear mass scale. At high masses, we show that assembly bias is expected from the statistics of the peaks of Gaussian random fluctuations, and we show that the extent of assembly bias found in N-body simulations of rare halos is in excellent agreement with our theoretical prediction. At low masses, we argue that assembly bias largely arises from a sub-population of low mass halos whose mass accretion has ceased. Due to their arrested development, these halos naturally become unbiased, in contrast to their anti-biased peers. We show that a simple toy model incorporating these effects can roughly reproduce the bias trends found in N-body simulations.

This is to announce a recent release of version 1.0.1a9 of parallel load balanced adaptive P3M cosmological N-body code named GRACOS. The code has been under development over years and is available for the download at http://www.gracos.org. GRACOS features: embedded script environment, cosmological initial conditions generator including non-gaussian initial conditions, particle data imager, mass density power spectrum estimator, particle position and velocity integrator with a choice between KDK and DKD integration schemes with Plummer force softening, file input and output with four supported serial and one dynamic distributed data format, standard installation procedure, and version specific documentation. GRACOS is released under the GNU General Public License (GPL). High efficiency is achieved with a number of implemented techniques such as timer-based load balancing using Hilbert space filling curve, adaptive P3M method for short range force computation, run-level compression in interprocessor communication, dynamic allocation for irregular domains and various sorting algorithms.

The project of an EAS Cherenkov array in the Tunka valley/Siberia with an area of about 1 km2 is presented. The new array will have a ten times bigger area than the existing Tunka-25 array and will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range from 10^15 to 10^18 eV.

We present a parallel implementation of the particle-particle/particle-mesh (P3M) algorithm for distributed memory clusters. The GRACOS (GRAvitational COSmology) code uses a hybrid method for both computation and domain decomposition. Long-range forces are computed using a Fourier transform gravity solver on a regular mesh; the mesh is distributed across parallel processes using a static one-dimensional slab domain decomposition. Short-range forces are computed by direct summation of close pairs; particles are distributed using a dynamic domain decomposition based on a space-filling Hilbert curve. A nearly-optimal method was devised to dynamically repartition the particle distribution so as to maintain load balance even for extremely inhomogeneous mass distributions. Tests using $800^3$ simulations on a 40-processor beowulf cluster showed good load balance and scalability up to 80 processes. We discuss the limits on scalability imposed by communication and extreme clustering and suggest how they may be removed by extending our algorithm to include adaptive mesh refinement.

The project of an EAS Cherenkov array in the Tunka valley/Siberia with an area of about 1 km^2 is presented. The new array will have a ten times bigger area than the existing Tunka-25 array and will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range from 10^15 to 10^18 eV.