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.

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.

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.