Compilation of papers presented by the JEM-EUSO Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin.

The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions.

EUSO-Balloon is a pathfinder mission for the Extreme Universe Space Observatory onboard the Japanese Experiment Module (JEM-EUSO). It was launched on the moonless night of the 25$^{th}$ of August 2014 from Timmins, Canada. The flight ended successfully after maintaining the target altitude of 38 km for five hours. One part of the mission was a 2.5 hour underflight using a helicopter equipped with three UV light sources (LED, xenon flasher and laser) to perform an inflight calibration and examine the detectors capability to measure tracks moving at the speed of light. We describe the helicopter laser system and details of the underflight as well as how the laser tracks were recorded and found in the data. These are the first recorded laser tracks measured from a fluorescence detector looking down on the atmosphere. Finally, we present a first reconstruction of the direction of the laser tracks relative to the detector.

We report the observation of a very high energy \gamma-ray source, whose position is coincident with HESS J1841-055. This source has been observed for 4.5 years by the ARGO-YBJ experiment from November 2007 to July 2012. Its emission is detected with a statistical significance of 5.3 standard deviations. Parameterizing the source shape with a two-dimensional Gaussian function we estimate an extension \sigma=(0.40(+0.32,-0.22}) degree, consistent with the HESS measurement. The observed energy spectrum is dN/dE =(9.0-+1.6) x 10^{-13}(E/5 TeV)^{-2.32-+0.23} photons cm^{-2} s^{-1} TeV^{-1}, in the energy range 0.9-50 TeV. The integral \gamma-ray flux above 1 TeV is 1.3-+0.4 Crab units, which is 3.2-+1.0 times the flux derived by HESS. The differences in the flux determination between HESS and ARGO-YBJ, and possible counterparts at other wavelengths are discussed.

As one of the brightest active blazars in both X-ray and very high energy $\gamma$-ray bands, Mrk 501 is very useful for physics associated with jets from AGNs. The ARGO-YBJ experiment is monitoring it for $\gamma$-rays above 0.3 TeV since November 2007. Starting from October 2011 the largest flare since 2005 is observed, which lasts to about April 2012. In this paper, a detailed analysis is reported. During the brightest $\gamma$-rays flaring episodes from October 17 to November 22, 2011, an excess of the event rate over 6 $\sigma$ is detected by ARGO-YBJ in the direction of Mrk 501, corresponding to an increase of the $\gamma$-ray flux above 1 TeV by a factor of 6.6$\pm$2.2 from its steady emission. In particular, the $\gamma$-ray flux above 8 TeV is detected with a significance better than 4 $\sigma$. Based on time-dependent synchrotron self-Compton (SSC) processes, the broad-band energy spectrum is interpreted as the emission from an electron energy distribution parameterized with a single power-law function with an exponential cutoff at its high energy end. The average spectral energy distribution for the steady emission is well described by this simple one-zone SSC model. However, the detection of $\gamma$-rays above 8 TeV during the flare challenges this model due to the hardness of the spectra. Correlations between X-rays and $\gamma$-rays are also investigated.

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

ARGO-YBJ is an air shower detector array with a fully covered layer of resistive plate chambers. It is operated with a high duty cycle and a large field of view. It continuously monitors the northern sky at energies above 0.3 TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period from 2007 November to 2010 February. This source was observed by the satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 421 was especially active in the first half of 2008. Many flares are observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and gamma-ray photons longer than 1 day is found. The evolution of the spectral energy distribution is investigated by measuring spectral indices at four different flux levels. Hardening of the spectra is observed in both X-ray and gamma-ray bands. The gamma-ray flux increases quadratically with the simultaneously measured X-ray flux. All these observational results strongly favor the synchrotron self-Compton process as the underlying radiative mechanism.

The sun blocks cosmic ray particles from outside the solar system, forming a detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ experiment in Tibet. Because the cosmic ray particles are positive charged, the magnetic field between the sun and the earth deflects them from straight trajectories and results in a shift of the shadow from the true location of the sun. Here we show that the shift measures the intensity of the field which is transported by the solar wind from the sun to the earth.

We use the results of a high resolution N-body simulation to investigate the role of the environment on the formation and evolution of galaxy-sized halos. Starting from a set of constrained initial conditions, we have produced a final configuration hosting a double cluster in one octant and a large void extending over two octants of the simulation box. We present results for two statistics: the relationship between 1-D velocity dispersion and mass and the probability distribution of the spin parameter $P(\lambda)$. The \svm relationship is well reproduced by the Truncated Isothermal Sphere (TIS) model introduced by Shapiro et al. (1999), although the slope is different from the original prediction. A series of \svm relationships for different values of the anisotropy parameter $\beta$, obtained using the theoretical predictions by Lokas and Mamon (2001) for NFW density profiles are found to be only marginally consistent with the data. Using some properties of the equilibrium TIS models, we construct subsamples of {\em fiducial} equilibrium TIS halos from each of the three subregions, and we study their properties. For these halos, we do find an environmental dependence of their properties, in particular of the spin parameter distribution $P(\lambda)$. We study in more detail the TIS model, and we find new relationships between the truncation radius and other structural parameters. No gravitationally bound halo is found having a radius larger than the critical value for gravithermal instability for TIS halos (\rt $\ge 34.2 r_{0}$, where $r_{0}$ is the core radius of the TIS solution). We do however find a dependence of this relationship on the environment, like for the $P(\lambda)$ statistics. These facts hint at a possible r\^{o}le of tidal fields at determining the statistical properties of halos.

We use Chandra data to map the gas temperature in the central region of the merging cluster A2142. The cluster is markedly nonisothermal; it appears that the central cooling flow has been disturbed but not destroyed by a merger. The X-ray image exhibits two sharp, bow-shaped, shock-like surface brightness edges or gas density discontinuities. However, temperature and pressure profiles across these edges indicate that these are not shock fronts. The pressure is reasonably continuous across these edges, while the entropy jumps in the opposite sense to that in a shock (i.e. the denser side of the edge has lower temperature, and hence lower entropy). Most plausibly, these edges delineate the dense subcluster cores that have survived a merger and ram pressure stripping by the surrounding shock-heated gas.

Using the results of a high mass resolution ($256^{3}$ particles) N-body simulation of a cluster-forming region, we study the statistic $\sigma_{v}-M$ for low mass halos (${\rm M} < 10^{12} {\rm M}_{\sun}$), and we compare it with three theoretical models. At the final redshift halos are well described by the Truncated Isothermal Sphere model recently introduced by Shapiro, Iliev & Raga (1999). We speculate that this is a consequence of the action of tidal fields on galaxy formation within clusters.

We present and discuss a method to identify substructures in combined angular-redshift samples of galaxies within Clusters. The method relies on the use of Discrete Wavelet Transform (hereafter DWT) and has already been applied to the analysis of the Coma cluster (Gambera et al. 1997). The main new ingredient of our method with respect to previous studies lies in the fact that we make use of a 3D data set rather than a 2D. We test the method on mock cluster catalogs with spatially localized substructures and on a N-body simulation. Our main conclusion is that our method is able to identify the existing substructures provided that: a) the subclumps are detached in part or all of the phase space, b) one has a statistically significant number of redshifts, increasing as the distance decreases due to redshift distortions; c) one knows {\it a priori} the scale on which substructures are to be expected. We have found that to allow an accurate recovery we must have both a significant number of galaxies ($\approx 200$ for clusters at z$\geq 0.4$ or about 800 at z$\leq$ 0.4) and a limiting magnitude for completeness $m_B=16$. The only true limitation to our method seems to be the necessity of knowing {\it a priori} the scale on which the substructure is to be found. This is an intrinsic drawback of the method and no improvement in numerical codes based on this technique could make up for it.

N-body algorithms for long-range unscreened interactions like gravity belong to a class of highly irregular problems whose optimal solution is a challenging task for present-day massively parallel computers. In this paper we describe a strategy for optimal memory and work distribution which we have applied to our parallel implementation of the Barnes & Hut (1986) recursive tree scheme on a Cray T3D using the CRAFT programming environment. We have performed a series of tests to find an " optimal data distribution " in the T3D memory, and to identify a strategy for the " Dynamic Load Balance " in order to obtain good performances when running large simulations (more than 10 million particles). The results of tests show that the step duration depends on two main factors: the data locality and the T3D network contention. Increasing data locality we are able to minimize the step duration if the closest bodies (direct interaction) tend to be located in the same PE local memory (contiguous block subdivison, high granularity), whereas the tree properties have a fine grain distribution. In a very large simulation, due to network contention, an unbalanced load arises. To remedy this we have devised an automatic work redistribution mechanism which provided a good Dynamic Load Balance at the price of an insignificant overhead.

By means of numerical methods, we solve the equations of motion for the collapse of a shell of baryonic matter, made of galaxies and substructure falling into the central regions of a cluster of galaxies, taking into account the effect of the dynamical friction. The parameters on which the dynamical friction mainly depends are: the peaks' height, the number of peaks inside a protocluster multiplied by the correlation function evaluated at the origin, the filtering radius and the nucleus radius of the protocluster of galaxies. We show how the collapse time (Tau) of the shell depends on these parameters. We give a formula that links the dynamical friction coefficient (Eta) o the parameters mentioned above and an analytic relation between the collapse time and (Eta). Finally, we obtain an analytical relation between (Tau) and the mean overdensity (mean Delta) within the shell. All the analytical relations that we find are in excellent agreement with the numerical integration.

Evidence for clustering within the Coma cluster is found by means of a multiscale analysis of the combined angular-redshift distribution. We have compiled a catalogue of 798 galaxy redshifts from published surveys from the region of the Coma cluster. We examine the presence of substructure and of voids at different scales ranging from $\sim 1 to \sim 16 h^{-1}$ Mpc, using subsamples of the catalogue, ranging from $cz=3000$ km/s to $cz=28000$ km/s. Our substructure detection method is based on the wavelet transform and on the segmentation analysis. The wavelet transform allows us to find out structures at different scales and the segmentation method allows us a quantitative statistical and morphological analysis of the sample. From the whole catalogue we select a subset of 320 galaxies, with redshifts between cz=5858 km/s and cz=8168 km/s that we identify as belonging to the central region of Coma and on which we have performed a deeper analysis, on scales ranging from $180 h^{-1}$ kpc to $1.44 h^{-1}$ Mpc. Our results are expressed in terms of the number of structures or voids and their sphericity for different values of the threshold detection and at all the scales investigated. According to our analysis, there is strong evidence for multiple hierarchical substructure, on scales ranging from a few hundreds of kpc to about $4 h^{-1}$ Mpc. The morphology of these substructures is rather spherical. On the scale of $720 h^{-1}$ kpc we find two main subclusters which where also found before, but our wavelet analysis shows even more substructures, whose redshift position is approximatively marked by these bright galaxies: NGC 4934 & 4840, 4889, 4898 & 4864, 4874 & 4839, 4927, 4875.

We solve numerically the equations of motion for the collapse of a shell of baryonic matter, made of galaxies and substructure of $ 10^{6} M_{\odot} \div 10^{9} M_{\odot}$, taking into account the dynamical friction and the parameters on which it depends: the peaks' height $\nu_{c}$, the number of peaks inside a protocluster $N_{tot}$ multiplied by the correlation function $\xi (r)$ evaluated at r=0, the filtering radius $R_{f}$ and the core radius of the cluster of galaxies, r_0. We show how the collapse time of a shell of baryonic matter depends on those parameters.

We describe a new parallel N-body code for cosmological simulations. The code is based on a work- and data sharing scheme, and is implemented within the Cray Research Corporation's CRAFT programming environment. Different data distribution schemes have been adopted for bodies' and tree's structures. Tests performed for two different types of initial distributions show that the performance scales almost ideally as a function of the size of the system and of the number of processors. We discuss the factors affecting the absolute speedup and how it can be increased with a better tree's data distribution scheme.