Carpet is an air-shower array at Baksan, Russia, equipped with a large-area muon detector, which makes it possible to separate primary photons from hadrons. We report first results of the search for primary photons with energies E>100 TeV. The experiment's ongoing upgrade and future sensitivity are also discussed.

We present a summary of recent tests and measurements of hadronic interaction properties with air showers. This report has a special focus on muon density measurements. Several experiments reported deviations between simulated and recorded muon densities in extensive air showers, while others reported no discrepancies. We combine data from eight leading air shower experiments to cover shower energies from PeV to tens of EeV. Data are combined using the z-scale, a unified reference scale based on simulated air showers. Energy-scales of experiments are cross-calibrated. Above 10 PeV, we find a muon deficit in simulated air showers for each of the six considered hadronic interaction models. The deficit is increasing with shower energy. For the models EPOS-LHC and QGSJet-II.04, the slope is found significant at 8 sigma.

Early results of the search for E_\gamma>1 PeV cosmic photons from point sources with the data of Carpet-2, an air-shower array equipped with a 175 m^2 muon detector, are presented. They include 95% CL upper limits on PeV photon fluxes from stacked directions of high-energy IceCube neutrino events and from four predefined sources, Crab, Cyg X-3, Mrk 421 and Mrk 501. An insignificant excess of events from Mrk 421 will be further monitored. Prospects of the use of the upgraded installation, Carpet-3 (410 m^2 muon detector), scheduled to start data taking in 2019, for searches of E_\gamma>100 TeV photons, are briefly discussed.

Carpet-2 is an air-shower array at Baksan Valley, Russia, equipped with a large-area (175 m^2) muon detector, which makes it possible to separate primary photons from hadrons. We report the first results of the search for primary photons with energies E_\gamma>1 PeV, directionally associated with IceCube high-energy neutrino events, in the data obtained in 3080 days of Carpet-2 live time.

There are multiple techniques to determine the chemical composition of the ultra-high-energy cosmic rays. While most of the methods are primarily sensitive to the average atomic mass, it is challenging to discriminate between the two lightest elements: proton and helium. In this paper, the proton-to-helium ratio in the energy range $10^{18.0} \mbox{eV}$ to $10^{19.3} \mbox{eV}$ is estimated using the tail of the distribution of the depth of the shower maximum $X_{\mathrm{max}}$. Using the exponential decay scale $\Lambda$ measured by the Pierre Auger Observatory and the Telescope Array experiment we derive the 68\%\,CL constraints on the proton-to-helium ratio $\mathrm{p}/\mathrm{He} > 7.3 $ and $\mathrm{p}/\mathrm{He} > 0.43 $ for $10^{18.0} < E < 10^{18.5}$ eV and $10^{18.3} < E < 10^{19.3}$ eV correspondingly. It is shown that the result is conservative with respect to the admixture of heavier elements. We evaluate the impact of the hadronic interaction model uncertainty. The implications for the astrophysical models of the origin of cosmic rays and the safety of the future colliders are discussed.

The Sydney University Giant Air-shower Recorder (SUGAR) measured the energy spectrum of ultra-high-energy cosmic rays reconstructed from muon-detector readings, while the Pierre Auger Observatory, looking at the same Southern sky, uses the calorimetric fluorescence method for the same purpose. Comparison of their two spectra allows us to reconstruct the empirical dependence of the number of muons in the shower on the primary energy for energies between $10^{17}$ and $10^{18.5}$ eV. We compare this dependence with the predictions of hadronic interaction models \mbox{QGSJET-II-04} and \mbox{EPOS-LHC}. The empirically determined number of muons with energies above 0.75 GeV exceeds the simulated one by the factors $\sim$1.67 and $\sim$1.28 for $10^{17}$ eV proton and iron primaries, respectively. The muon excess grows moderately with the primary energy, increasing by an additional factor of $\sim 1.2$ for $10^{18.5}$ eV primaries.

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.

The luminosity of fading type Ia supernovae is governed by radioactive decays of 56Ni and 56Co. The decay rates are proportional to the Fermi coupling constant G_F and, therefore, are determined by the vacuum expectation value v of the Brout-Englert-Higgs field. We use publicly available sets of light curves of type Ia supernova at various redshifts to constrain possible spacetime variations of the 56Ni decay rate. The resulting constraint is not very tight; however, it is the only direct bound on the variation of the decay rate for redshifts up to z~1. We discuss potential applications of the result to searches for non-constancy of G_F and v.