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.

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.

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.