Some discrepancies have been reported between observed and simulated muon content of extensive air showers: the number of observed muons exceeded the expectations in HiRes-MIA, Yakutsk and Pierre Auger Observatory data. Here, we analyze the data of the Moscow State University Extensive Air Shower (EAS-MSU) array on E_mu>~10 GeV muons in showers caused by ~(10^17-10^18) eV primary particles and demonstrate that they agree with simulations (QGSJET-II-04 hadronic interaction model) once the primary composition inferred from the surface-detector data is assumed.

Results of the search for $\sim (10^{16} - 10^{17.5})$ eV primary cosmic-ray photons with the data of the Moscow State University (MSU) Extensive Air Shower (EAS) array are reported. The full-scale reanalysis of the data with modern simulations of the installation does not confirm previous indications of the excess of gamma-ray candidate events. Upper limits on the corresponding gamma-ray flux are presented. The limits are the most stringent published ones at energies $\sim 10^{17}$ eV.

The Moscow State University Extensive Air Shower (EAS-MSU) array studied high-energy cosmic rays with primary energies ~(1-500) PeV in the Northern hemisphere. The EAS-MSU data are being revisited following recently found indications to an excess of muonless showers, which may be interpreted as the first observation of cosmic gamma rays at ~100 PeV. In this paper, we present a complete Monte-Carlo model of the surface detector which results in a good agreement between data and simulations. The model allows us to study the performance of the detector and will be used to obtain physical results in further studies.

Results of a blind search for localized regions of an excessive flux of cosmic rays in the energy range from 50 TeV to 20 PeV with the data of the FIAN KLARA-Chronotron experiment, the EAS MSU array and the Prototype of the EAS-1000 array are presented. A number of regions with a significant excess of the registered flux over an expected isotropic background are found. Some of the regions are present in at least two of the data sets considered.

Archival EAS-MSU data are searched for anomalous muonless events which may be caused by primary gamma rays with energies between 10^15 eV and 10^18 eV. We consider a refined sample of high-quality data and confirm the previously reported detection of a non-zero gamma-ray flux at ~5x10^16 eV with a similar flux value but at somewhat lower statistical significance, corresponding to the depletion of the sample. We present upper limits on the flux below and above these energies, including the first constraints in the range (10^17-10^18) eV never studied by any other experiment.

We reanalyze archival EAS-MSU data in order to search for events with anomalously low content of muons with energies E_mu>10 GeV in extensive air showers with number of particles N_e>2*10^7. We confirm the first evidence for nonzero flux of primary cosmic gamma rays at energies E~10^17 eV. The estimated fraction of primary gamma rays in the flux of cosmic particles with energies E>5.4*10^16 eV is (0.43 +0.12-0.11)%, which corresponds to the intensity of (1.2 +0.4-0.3)*10^{-16} cm^{-2} s^{-1} sr^{-1}. The study of arrival directions does not favour any particular mechanism of the origin of the photon-like events.

The origin of ultra-high-energy cosmic rays (UHECRs) is one of the most intriguing problems of modern cosmic ray physics. We briefly review the main astrophysical models of their origin and the forthcoming orbital experiments TUS and JEM-EUSO, and discuss how the new data can help one solve the long-standing puzzle.

Recent results of Milagro, Tibet, ARGO-YBJ and IceCube experiments on the small-scale anisotropy of Galactic cosmic rays (CRs) with energies from units up to a few hundred TeV arise a question on a possible nature of the observed phenomenon, as well as on the anisotropy of CRs at higher energies. An analysis of a small-scale anisotropy of CRs with energies at around PeV registered with the EAS MSU array presented in the article, reveals a number of regions with an excessive flux. A typical size of the regions varies from 3 up to 12 degrees. We study correlation of these regions with positions of potential astrophysical sources of CRs and discuss a possible origin of the observed anisotropy.

Three regions of excessive flux of cosmic rays with energies of the order of PeV are found in the experimental data of the EAS MSU array at a confidence level greater than 4\sigma. For two of them, there are similar regions in the experimental data of the EAS-1000 Prototype array. One of the interesting features of the regions is the absence of supernova remnants in their vicinities, traditionally considered as the main sources of Galactic cosmic rays, but the presence of isolated pulsars, some of which are able to accelerate heavy nuclei up to energies close to PeV. In our opinion, this favors the assumption that isolated pulsars are able to contribute to the flux of Galactic cosmic rays more than is usually assumed.

An analysis of arrival directions of extensive air showers (EAS) registered with the EAS MSU and EAS-1000 prototype arrays has revealed a region of excessive flux of PeV cosmic rays in the direction toward pulsars PSR J1840+5640 and LAT PSR J1836+5925 at significance level up to 4.5sigma. The first of the pulsars was discovered almost 30 years ago and is a well-studied old radio pulsar located at the distance of 1.7pc from the Solar system. The second pulsar belongs to a new type of pulsars, discovered by the space gamma-ray observatory Fermi, pulsations of which are not observed in optical and radio wavelengths but only in the gamma-ray range of energies (gamma-ray-only pulsars). In our opinion, the existence of the region of excessive flux of cosmic rays registered with two different arrays provides a strong evidence that isolated pulsars can give a noticeable contribution to the flux of Galactic cosmic rays in the PeV energy range.

Angular distributions of extensive air showers with different number of charged particles in the range 2.5x10^5--4x10^7 are derived using the experimental data obtained with the EAS MSU array. Possible approximations of the obtained distributions with different empiric functions available in literature, are analysed. It is shown that the exponential function provides the best approximation of the angular distributions in the sense of the chi-squared criterion.

An analysis of arrival directions of extensive air showers generated by cosmic rays in the PeV energy range and registered with the EAS MSU and EAS-1000 Prototype arrays reveals a considerable number of regions of excessive flux of cosmic rays. We present results of comparative analysis of regions found in the two data sets, estimate probabilities of their appearance, and discuss correlation of their locations with coordinates of possible astrophysical sources of PeV cosmic rays.

New results of an analysis of arrival directions of extensive air showers registered with the EAS--1000 Prototype array from August 1997 till February 1999 are presented. The method of Alexandreas et al., which has been used for analysis of data registered with CYGNUS, Milagrito, HEGRA AIROBICC, KASCADE and a number of other experiments, is employed. The existence of zones of excessive flux of cosmic rays with energies in the region of the knee is confirmed, as well as closeness of the zones to coordinates of possible astrophysical cosmic ray sources.

We study arrival directions of 1.4x10^6 extensive air showers (EAS) registered with the EAS--1000 Prototype Array in the energy range 0.1--10 PeV. By applying an iterative algorithm that provides uniform distribution of the data with respect to sidereal time and azimuthal angles, we find a number of zones with excessive flux of cosmic rays (CRs) at >=3 sigma level. We compare locations of the zones with positions of galactic supernova remnants (SNRs), pulsars, open star clusters, and regions of ionized hydrogen and find remarkable coincidences, which may witness in favour of the hypothesis that certain objects of these types, including the SNRs Cassiopeia A, the Crab Nebula, the Monogem Ring and some other, provide a noticeable contribution to the flux of CRs in the PeV range of energies. In addition, we find certain signs of a contribution from the M33 galaxy and a number of comparatively nearby groups of active galactic nuclei and interacting galaxies, in particular those in the Virgo cluster of galaxies. The results also provide some hints for a search of possible sources of ultra-high energy (UHE) cosmic rays and support an earlier idea that a part of both UHE and PeV CRs may originate from the same astrophysical objects.

Bursts of the count rate of extensive air showers (EAS) lead to the appearance of clusters in time series that represent EAS arrival times. We apply methods of nonlinear time series analysis to twenty EAS cluster events found in the data set obtained with the EAS-1000 prototype array. In particular, we use the Grassberger-Procaccia algorithm to compute the correlation dimension of the time series in the vicinity of the clusters. We find that four cluster events produce signs of chaos in the corresponding time series. By applying a number of supplementary methods we assess that the nature of the observed behaviour of the correlation dimension is likely to be deterministic. We suggest a simple qualitative model that might explain an origin of clusters in general and "possibly chaotic" clusters in particular. Finally, we compare our conclusions with the results of similar investigations performed by the EAS-TOP and LAAS groups.

We perform cluster analysis of arrival times of extensive air showers (EAS) registered with the EAS-1000 Prototype array in the period from August, 1997 till February, 1999. We present twenty cluster events each consisting of one or several EAS clusters, study the dynamics of their development, and analyze the angular distribution of EAS in clusters. We find that there may be certain correlation between EAS clusters with mean electron number $\sim10^5$ on the one hand and gamma-ray bursts and ultrahigh energy cosmic rays on the other.

We apply certain methods of nonlinear time series analysis to the extensive air shower clusters found earlier in the data set obtained with the EAS-1000 Prototype array. In particular, we use the Grassberger-Procaccia algorithm to compute the correlation dimension of samples in the vicinity of the clusters. The validity of the results is checked by surrogate data tests and some additional quantities. We compare our conclusions with the results of similar investigations performed by the EAS-TOP and LAAS groups.

We have analysed arrival times of extensive air showers (EAS) registered with the EAS-1000 prototype array during the period from August, 1997 till February, 1999. Our analysis has revealed that though the vast majority of samples of consecutive time intervals between EAS arrival times obey the exponential distribution, there are sequences of showers that have another distribution and thus violate the homogeneity hypothesis. The search for correlation between such events and clusters of showers and events with big delays between arrival times was also carried out.