Muon content of ultrahigh-energy air showers: Yakutsk data versus simulations

A sample of 33 extensive air showers (EASs) with estimated primary energies above 2 × 10 19 eV and high-quality muon data recorded by the Yakutsk EAS array is analyzed. The observed muon density is compared event-by-event 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 good agreement with the expected composition in which the light component corresponds to protons and the heavy component to iron-like nuclei. With SIBYLL, simulated muon densities for iron primaries are a factor of ∼ 1.5 less than those observed for the heavy component for the same electromagnetic signal. Assuming a two-component proton-iron composition and the EPOS model, the fraction of protons with energies E > 10 19 eV is 0.52 −0.20 +0.19 at the 95% C.L.