Comparing the signals measured by the surface and underground scintillator detectors of the Yakutsk Extensive Air Shower Array, we place upper limits on the integral flux and the fraction of primary cosmic-ray photons with energies E > 10^18 eV, E > 2*10^18 eV and E > 4*10^18 eV. The large collected statistics of the showers measured by large-area muon detectors provides a sensitivity to photon fractions < 10^-2, thus achieving precision previously unreachable at ultra-high energies.

The energies of the most energetic extensive air showers observed at the Yakutsk array have been estimated with help of the all detectors readings instead of using of the standard procedure with a parameter s(600). The energy of the most energetic extensive air shower observed at the Yakutsk array happened to be 200, 200, 180 and 165 EeV with the values of the Xi**2 function per one degree of freedom 0.9, 1., 0.9 and 1.1 for the primary protons and helium, oxygen and iron nuclei accordingly.

The fluxes of electrons, positrons, gammas, Cherenkov photons and muons in individual extensive air showers induced by the primary protons and helium, oxygen and iron nuclei at the level of observation have been estimated with help of the code CORSICA 6.616. The comparison show that the values of the function Xi**2 per one degree of freedom changes from 1.1 for iron nuclei to 0.9 for primary protons. As this difference is small all readings of detectors of the Vavilov-Cherenkov radiation have been used. At last, readings of underground detectors of muons with energies above 1 GeV have been exploited to make definite conclusion about chemical composition.

We analyse a sample of 33 extensive air showers (EAS) with estimated primary energies above 2\cdot 10^{19} eV and high-quality muon data recorded by the Yakutsk EAS array. We compare, event-by-event, the observed muon density to that expected from CORSIKA simulations for primary protons and iron, using SIBYLL and EPOS hadronic interaction models. The study suggests the presence of two distinct hadronic components, ``light'' and ``heavy''. Simulations with EPOS are in a good agreement with the expected composition in which the light component corresponds to protons and the heavy component to iron-like nuclei. With SYBILL, simulated muon densities for iron primaries are a factor of \sim 1.5 less than those observed for the heavy component, for the same electromagnetic signal. Assuming two-component proton-iron composition and the EPOS model, the fraction of protons with energies E>10^{19} eV is 0.52^{+0.19}_{-0.20} at 95% confidence level.

By making use of the data on the total signal and on the muon component of the air showers detected by the Yakutsk array, we analyze, in the frameworks of the recently suggested event-by-event approach, how large the fraction of primary gamma-rays at ultra-high energies can be. We derive upper limits on the photon fraction in the integral flux of primary cosmic rays. At the 95% confidence level (CL), these limits are 22% for primary energies E_0>4\cdot 10^{19}eV and 12% for E_0>2\cdot 10^{19}eV. Despite the presence of muonless events, the data are consistent with complete absence of photons at least at 95% CL. The sensitivity of the results to systematic uncertainties, in particular to those of the energy determination for non-photon primaries, is discussed.

We present the interpretation of the muon and scintillation signals of ultra-high-energy air showers observed by AGASA and Yakutsk extensive air shower array experiments. We consider case-by-case ten highest energy events with known muon content and conclude that at the 95% confidence level (C.L.) none of them was induced by a primary photon. Taking into account statistical fluctuations and differences in the energy estimation of proton and photon primaries, we derive an upper limit of 36% at 95% C.L. on the fraction of primary photons in the cosmic-ray flux above 10^20 eV. This result disfavors the Z-burst and superheavy dark-matter solutions to the GZK-cutoff problem.

It is proposed a new approach for estimating the composition of cosmic rays. It is found that the zenith angle distributions and muon components of EAS for energies E>10^19 eV and E>4.10^19 eV differ from each other. It is shown that the cosmic rays above E>4.10^19 eV is heavier than the cosmic rays at energy E~ 10^19 eV. According to our estimation the SUGAR array detected 8 showers above 10^20 eV. It is concluded that no sign of Greisen-Zatsepin-Kuzmin (GZK) cut off in the spectrum of cosmic rays and all cosmic rays are galactic.

The ratio of the muon flux density to charged particle flux density at distances of 300 and 600 m from the shower axis ($\rhom(300)/\rhos(300)$ and $\rhom(600)/\rhos(600)$) is measured. In addition, the energy dependence of $\rhom(1000)$ is analysed for showers with energies above $10^{18}$ eV. A comparison between the experimental data and calculations performed with the QGSJET model is given for the cases of primary proton, iron nucleus and gamma- ray. We conclude that the showers with $\E\ge3\times10^{18}$ eV can be formed by light nuclei with a pronounced fraction of protons and helium nuclei. It is not excluded however that a small part of showers with energies above $10^{19}$ eV could be initiated by primary gamma-rays.

The energy spectrum of primary cosmic rays with ultra high energies based on the Yakutsk EAS Array data is presented. For the largest events values of S600 and axis coordinates have been obtained using revised lateral distribution function. The affect of the arrival time distribution at several axis distance on estimated density for Yakutsk and AGASA is considered.

There are some discrepancies in the results on energy spectrum from Yakutsk, AGASA, and HiRes experiments. In this work differential energy spectrum of primary cosmic rays based on the Yakutsk EAS Array data is presented. For the largest events values of $S_{600}$ and axes coordinates have been obtained using revised lateral distribution function. Simulation of converter's response at large distances showed no considerable underestimation of particle density. Complex shape of spectrum in region of $\E > 10^{17}$ eV is confirmed. After adjustment of parameters and additional exposition at the Yakutsk array there are three events with energy $\E > 10^{20}$ eV.

In the work the data of the Yakutsk complex EAS array and their comparison with calculation in the case of primary nuclei of different chemical elements are presented. The calculation by QGSJET model have been used interpreting experimental data.

The energy spectrum of UHECRs is discussed on the basis of the Yakutsk array database analysis. In the region E=0.1 to 30 EeV the showers are detected under trigger-500, while at energies above 30 EeV the whole acceptance area for trigger-1000 is used in order to utilize all the data available in the region of GZK cutoff.

Arrival directions of extensive air showers by using world data are considered. It is shown that distributions of showers in zenith angle at E>10^{19}eV and E>4.10^{19}eV differ from each other. On this basis the conclusions are made: the estimation of the energy of inclined showers at E>4.10^{19}eV at Yakutsk and AGASA arrays is incorrect; the cosmic ray spectrum is not break up and most likely continued up to E=10^{20}eV and higher. It is found that ultrahigh energy cosmic rays are, probably, superheavy nuclei.