Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

The study of the ultra-high energy cosmic rays, neutrinos and gamma rays is one of the most important challenges in astrophysics. The low fluxes of these particles do not allow one to detect them directly. The detection is performed by the measuring of the air-showers produced by the primary particles in the Earth's atmosphere. A radio detection of ultra-high energy air-showers is a cost-effective technique that provides a precise reconstruction of the parameters of primary particle and almost full duty cycle in comparison with other methods. The main challenge of the modern radio detectors is the development of efficient self-trigger technology, resistant to high-level background and radio frequency interference. Most of the modern radio detectors receive trigger generated by either particle or optical detectors. The development of the self trigger for the radio detector will significantly simplify the operation of existing instruments and allow one to access the main advantages of the radio method as well as open the way to the construction of the next generation of large-scale radio detectors. In the present work we discuss our progress in the solution of this problem, particularly the classification of broadband pulses.

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$\sigma$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $\alpha_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

With the Surface Detector array (SD) of the Pierre Auger Observatory we can detect neutrinos with energy between $10^{17}\,$eV and $10^{20}\,$eV from point-like sources across the sky, from close to the Southern Celestial Pole up to $60^\circ$ in declination, with peak sensitivities at declinations around $\sim -53^\circ$ and $\sim+55^\circ$, and an unmatched sensitivity for arrival directions in the Northern hemisphere. A search has been performed for highly-inclined air showers induced by neutrinos of all flavours with no candidate events found in data taken between 1 Jan 2004 and 31 Aug 2018. Upper limits on the neutrino flux from point-like steady sources have been derived as a function of source declination. An unrivaled sensitivity is achieved in searches for transient sources with emission lasting over an hour or less, if they occur within the field of view corresponding to the zenith angle range between $60^\circ$ and $~95^\circ$ where the SD of the Pierre Auger Observatory is most sensitive to neutrinos.

The determination of the primary energy of extensive air showers using the fluorescence detection technique requires an estimation of the energy carried away by particles that do not deposit all their energy in the atmosphere. This estimation is typically made using Monte Carlo simulations and thus depends on the assumed primary particle mass and on model predictions for neutrino and muon production. In this work we present a new method to obtain the invisible energy from events detected by the Pierre Auger Observatory. The method uses measurements of the muon number at ground level, and it allows us to reduce significantly the systematic uncertainties related to the mass composition and the high energy hadronic interaction models, and consequently to improve the estimation of the energy scale of the Observatory.

Neutrinos with energies above $10^{17}$ eV are detectable with the Surface Detector Array of the Pierre Auger Observatory. The identification is efficiently performed for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for Earth-skimming $\tau$ neutrinos with nearly tangential trajectories relative to the earth. No neutrino candidates were found in $\sim\,14.7$ years of data taken up to 31 August 2018. This leads to restrictive upper bounds on their flux. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high-energy neutrinos with an $E_\nu^{-2}$ spectrum in the energy range $1.0 \times 10^{17}~{\rm eV} - 2.5 \times 10^{19}~{\rm eV}$ is $E^2 {\rm d}N_\nu/{\rm d}E_\nu < 4.4 \times 10^{-9}~{\rm GeV~cm^{-2}~s^{-1}~sr^{-1}}$, placing strong constraints on several models of neutrino production at EeV energies and on the properties of the sources of ultra-high-energy cosmic rays.

Contributions of the Pierre Auger Collaboration to the 36th International Cosmic Ray Conference (ICRC 2019), 24 July - 1 August 2019, Madison, Wisconsin, USA.

Sparse digital antenna arrays constitute a promising detection technique for future large-scale cosmic-ray observatories. It has recently been shown that this kind of instrumentation can provide a resolution of the energy and of the shower maximum on the level of other cosmic-ray detection methods. Due to the dominant geomagnetic nature of the air-shower radio emission in the traditional frequency band of 30 to 80 MHz, the amplitude and polarization of the radio signal strongly depend on the azimuth and zenith angle of the arrival direction. Thus, the estimation of the efficiency and subsequently of the aperture of an antenna array is more complex than for particle or Cherenkov-light detectors. We have built a new efficiency model based on utilizing a lateral distribution function as a shower model, and a probabilistic treatment of the detection process. The model is compared to the data measured by the Tunka Radio Extension (Tunka-Rex), a digital antenna array with an area of about 1 km$^2$ located in Siberia at the Tunka Advanced Instrument for Cosmic rays and Gamma Ray Astronomy (TAIGA). Tunka-Rex detects radio emission of air showers using trigger from air-Cherenkov and particle detectors. The present study is an essential step towards the measurement of the cosmic-ray flux with Tunka-Rex, and is important for radio measurements of air showers in general.

The Tunka Radio Extension (Tunka-Rex) is a digital antenna array located in the Tunka Valley in Siberia, which measures the radio emission of cosmic-ray air-showers with energies up to EeV. Tunka-Rex is externally triggered by the Tunka-133 air-Cherenkov timing array (during nights) and by the Tunka-Grande array of particle detectors (remaining time). These three arrays comprise the cosmic-ray extension of the Tunka Advanced Instrument for cosmic rays and Gamma Astronomy (TAIGA). The configuration and analysis pipeline of Tunka-Rex have significantly changed over its runtime. Density of the antennas was tripled and the pipeline has become more developed forming now sophisticated piece of reconstruction software. During its lifecycle Tunka-Rex has demonstrated that a cost-effective and full duty-cycle radio detector can reconstruct the energy and shower maximum with a precision comparable to optical detectors. Moreover, it was shown that cosmic-ray instruments, that use different detection techniques and are placed in different locations, can be cross-calibrated via their radio extensions. These results show the prospects of application of the radio technique for future large-scale experiments for cosmic-ray and neutrino detection. For the time being Tunka-Rex has ceased active measurements and focuses on the data analysis and publication of corresponding software and data in an open-access data center with online analysis features. In this report we present the current status of the array and give an overview of the results achieved during these years as well as discuss upcoming improvements in instrumentation and data analysis, which can be applied for the future radio arrays.

The Tunka Radio Extension (Tunka-Rex) is a digital antenna array operating in the frequency band of 30-80 MHz, measuring the radio emission of air-showers induced by ultra-high energy cosmic rays. Tunka-Rex is co-located with the TAIGA experiment in Siberia and consists of 63 antennas, 57 of them in a densely instrumented area of about 1km2. The signals from the air showers are short pulses, which have a duration of tens of nanoseconds and are recorded in traces of about 5{\mu}s length. The Tunka-Rex analysis of cosmic-ray events is based on the reconstruction of these signals, in particular, their positions in the traces and amplitudes. This reconstruction suffers at low signal-to-noise ratios, i.e. when the recorded traces are dominated by background. To lower the threshold of the detection and increase the efficiency, we apply advanced methods of signal reconstruction, namely matched filtering and deep neural networks with autoencoder architecture. In the present work we show the comparison between the signal reconstructions obtained with these techniques, and give an example of the first reconstruction of the Tunka-Rex signals obtained with a deep neural networks.

The Tunka Radio Extension (Tunka-Rex) is a cosmic-ray detector operating since 2012. The detection principle of Tunka-Rex is based on the radio technique, which impacts data acquisition and storage. In this paper we give a first detailed overview of the concept of the Tunka-Rex Virtual Observatory (TRVO), a framework for open access to the Tunka-Rex data, which currently is under active development and testing. We describe the structure of the data, main features of the interface and possible applications of the TRVO.

The Tunka Radio Extension (Tunka-Rex) is a digital antenna array, which measures the radio emission of the cosmic-ray air-showers in the frequency band of 30-80 MHz. Tunka-Rex is co-located with TAIGA experiment in Siberia and consists of 63 antennas, 57 of them are in a densely instrumented area of about 1 km\textsuperscript{2}. In the present work we discuss the improvements of the signal reconstruction applied for the Tunka-Rex. At the first stage we implemented matched filtering using averaged signals as template. The simulation study has shown that matched filtering allows one to decrease the threshold of signal detection and increase its purity. However, the maximum performance of matched filtering is achievable only in case of white noise, while in reality the noise is not fully random due to different reasons. To recognize hidden features of the noise and treat them, we decided to use convolutional neural network with autoencoder architecture. Taking the recorded trace as an input, the autoencoder returns denoised trace, i.e. removes all signal-unrelated amplitudes. We present the comparison between standard method of signal reconstruction, matched filtering and autoencoder, and discuss the prospects of application of neural networks for lowering the threshold of digital antenna arrays for cosmic-ray detection.

The profile of the longitudinal development of showers produced by ultra-high energy cosmic rays carries information related to the interaction properties of the primary particles with atmospheric nuclei. In this work, we present the first measurement of the average shower profile in traversed atmospheric depth at the Pierre Auger Observatory. The shapes of profiles are well reproduced by the Gaisser-Hillas parametrization within the range studied, for E > 10^{17.8} eV. A detailed analysis of the systematic uncertainties is performed using 10 years of data and a full detector simulation. The average shape is quantified using two variables related to the width and asymmetry of the profile, and the results are compared with predictions of hadronic interaction models for different primary particles.

An overview of the multi-messenger capabilities of the Pierre Auger Observatory is presented. The techniques and performance of searching for Ultra-High Energy neutrinos, photons and neutrons are described. Some of the most relevant results are reviewed, such as stringent upper bounds that were placed to a flux of diffuse cosmogenic neutrinos and photons, bounds placed on neutrinos emitted from compact binary mergers that were detected by LIGO and Virgo during their first and second observing runs, as well as searches for high energy photons and neutrons from the Galactic center that constrain the properties of the putative Galactic PeVatron, observed by the H.E.S.S.\ collaboration. The observation of directional correlations between ultra-high energy cosmic rays and either high energy astrophysical neutrinos or specific source populations, weighted by their electromagnetic radiation, are also discussed. They constitute additional multi-messenger approaches aimed at identifying the sources of high energy cosmic rays.

The paper is a script of a lecture given at the ISAPP-Baikal summer school in 2018. The lecture gives an overview of the Tunka Advanced Instrument for cosmic rays and Gamma Astronomy (TAIGA) facility including historical introduction, description of existing and future setups, and outreach and open data activities.

The Tunka Radio Extension (Tunka-Rex) is an antenna array spread over an area of about 1~km$^2$. The array is placed at the Tunka Advanced Instrument for cosmic rays and Gamma Astronomy (TAIGA) and detects the radio emission of air showers in the band of 30 to 80~MHz. During the last years it was shown that a sparse array such as Tunka-Rex is capable of reconstructing the parameters of the primary particle as accurate as the modern instruments. Based on these results we continue developing our data analysis. Our next goal is the reconstruction of cosmic-ray energy spectrum observed only by a radio instrument. Taking a step towards it, we develop a model of aperture of our instrument and test it against hybrid TAIGA observations and Monte-Carlo simulations. In the present work we give an overview of the current status and results for the last five years of operation of Tunka-Rex and discuss prospects of the cosmic-ray energy estimation with sparse radio arrays.

The Tunka Radio Extension (Tunka-Rex) is a digital antenna array for the detection of radio emission from cosmic-ray air showers in the frequency band of 30 to 80 MHz and for primary energies above 100 PeV. The standard analysis of Tunka-Rex includes events with zenith angle of up to 50$^\circ$. This cut is determined by the efficiency of the external trigger. However, due to the air-shower footprint increasing with zenith angle and due to the more efficient generation of radio emission (the magnetic field in the Tunka valley is almost vertical), there are a number of ultra-high-energy inclined events detected by Tunka-Rex. In this work we present a first analysis of a subset of inclined events detected by Tunka-Rex. We estimate the energies of the selected events and test the efficiency of Tunka-Rex antennas for detection of inclined air showers.

The Tunka Radio Extension (Tunka-Rex) is a digital radio array operating in the frequency band of 30-80 MHz and detecting radio emission from air-showers produced by cosmic rays with energies above 100 PeV. The experiment is installed at the site of the TAIGA (Tunka Advanced Instrument for cosmic rays and Gamma Astronomy) observatory and performs joint measurements with the co-located particle and air-Cherenkov detectors in passive mode receiving a trigger from the latter. Tunka-Rex collects data since 2012, and during the last five years went through several upgrades. As a result the density of the antenna field was increased by three times since its commission. In this contribution we present the latest results of Tunka-Rex experiment, particularly an updated analysis and efficiency study, which have been applied to the measurement of the mean shower maximum as a function of energy for cosmic rays of energies up to EeV. The future plans are also discussed: investigations towards an energy spectrum of cosmic rays with Tunka-Rex and their mass composition using a combination of Tunka-Rex data with muon measurements by the particle detector Tunka-Grande.

We present a detailed study of the large-scale anisotropies of cosmic rays with energies above 4 EeV measured using the Pierre Auger Observatory. For the energy bins [4,8] EeV and $E\geq 8$ EeV, the most significant signal is a dipolar modulation in right ascension at energies above 8 EeV, as previously reported. In this paper we further scrutinize the highest-energy bin by splitting it into three energy ranges. We find that the amplitude of the dipole increases with energy above 4 EeV. The growth can be fitted with a power law with index $\beta=0.79\pm 0.19$. The directions of the dipoles are consistent with an extragalactic origin of these anisotropies at all the energies considered. Additionally we have estimated the quadrupolar components of the anisotropy: they are not statistically significant. We discuss the results in the context of the predictions from different models for the distribution of ultrahigh-energy sources and cosmic magnetic fields.

With the Auger Engineering Radio Array (AERA) of the Pierre Auger Observatory, we have observed the radio emission from 561 extensive air showers with zenith angles between 60$^\circ$ and 84$^\circ$. In contrast to air showers with more vertical incidence, these inclined air showers illuminate large ground areas of several km$^2$ with radio signals detectable in the 30 to 80\,MHz band. A comparison of the measured radio-signal amplitudes with Monte Carlo simulations of a subset of 50 events for which we reconstruct the energy using the Auger surface detector shows agreement within the uncertainties of the current analysis. As expected for forward-beamed radio emission undergoing no significant absorption or scattering in the atmosphere, the area illuminated by radio signals grows with the zenith angle of the air shower. Inclined air showers with EeV energies are thus measurable with sparse radio-antenna arrays with grid sizes of a km or more. This is particularly attractive as radio detection provides direct access to the energy in the electromagnetic cascade of an air shower, which in case of inclined air showers is not accessible by arrays of particle detectors on the ground.

We present an improved method for the precise reconstruction of cosmic ray air showers above $10^{17}$ eV with sparse radio arrays. The method is based on the comparison of predictions for radio pulse shapes by CoREAS simulations to measured pulses. We applied our method to the data of Tunka-Rex, a 1 km$^2$ radio array in Siberia operating in the frequency band of 30-80 MHz. Tunka-Rex is triggered by the air-Cherenkov detector Tunka-133 and by scintillators (Tunka-Grande). The instrument collects air-shower data since 2012. The present paper describes updated data and signal analyses of Tunka-Rex and details of a new method applied. After efficiency cuts, when Tunka-Rex reaches its full efficiency, the energy resolution of about 10% given by the new method has reached the limit of systematic uncertainties due to the calibration uncertainty and shower-to-shower fluctuations. At the same time the shower maximum reconstruction is significantly improved up to an accuracy of 35 g/cm$^2$ compared to the previous method based on the slope of the lateral distribution. We also define and now achieved conditions of the measurements, at which the shower maximum resolution of Tunka-Rex reaches a value of 25 g/cm$^2$ and becomes competitive to optical detectors. To check and validate our reconstruction and efficiency cuts we compare individual events to the reconstruction of Tunka-133. Furthermore, we compare the mean of shower maximum as a function of primary energy to the measurements of other experiments.

Tunka-Rex (Tunka Radio Extension) is an antenna array for cosmic-ray detection located in Siberia. Previous studies of cosmic rays with Tunka-Rex have shown high precision in determining the energy of the primary particle and the possibility to reconstruct the depth of the shower maximum. The next step is the reconstruction of the mass composition and the energy spectrum of cosmic rays. One of the main problems appearing within this task is to estimate the detection efficiency of the instrument, and the exposure of the observations. The detection efficiency depends on properties of the primary cosmic rays, such as energy and arrival direction, as well as on many parameters of the instrument: density of the array, efficiency of the receiving antennas, signal-detection threshold, data-acquisition acceptance, and trigger properties. More than that, the configuration of detector changes with time. During the measurements some parts of the detector can provide corrupted data or sometimes do not operate. All these features should be taken into account for an estimation of the detection efficiency. For each energy and arrival direction we estimate the detection probability and effective area of the instrument. To estimate the detection probability of a shower we use a simple Monte Carlo model, which predicts the size of the footprint of the radio emission as function of the primary energy and arrival direction (taking into account the geometry of Earth's magnetic field). Combining these approaches we calculate the event statistics and exposure for each run. This is the first accurate study of the exposure for irregular large-scale radio arrays taking into account most important features of detection, which will be used for the measurement of primary cosmic-ray spectra with Tunka-Rex.

The Tunka Radio Extension (Tunka-Rex) is an array of 63 antennas located in the Tunka Valley, Siberia. It detects radio pulses in the 30-80 MHz band produced during the air-shower development. As shown by Tunka-Rex, a sparse radio array with about 200 m spacing is able to reconstruct the energy and the depth of the shower maximum with satisfactory precision using simple methods based on parameters of the lateral distribution of amplitudes. The LOFAR experiment has shown that a sophisticated treatment of all individually measured amplitudes of a dense antenna array can make the precision comparable with the resolution of existing optical techniques. We develop these ideas further and present a method based on the treatment of time series of measured signals, i.e. each antenna station provides several points (trace) instead of a single one (amplitude or power). We use the measured shower axis and energy as input for CoREAS simulations: for each measured event we simulate a set of air-showers with proton, helium, nitrogen and iron as primary particle (each primary is simulated about ten times to cover fluctuations in the shower maximum due to the first interaction). Simulated radio pulses are processed with the Tunka-Rex detector response and convoluted with the measured signals. A likelihood fit determines how well the simulated event fits to the measured one. The positions of the shower maxima are defined from the distribution of chi-square values of these fits. When using this improved method instead of the standard one, firstly, the shower maximum of more events can be reconstructed, secondly, the resolution is increased. The performance of the method is demonstrated on the data acquired by the Tunka-Rex detector in 2012-2014.

Tunka-Rex is detector for radio emission produced by cosmic-ray air-showers located in Siberia, triggered by Tunka-133, a co-located air-Cherenkov detector during night, and by a scintillator array Tunka-Grande during day. Tunka-Rex demonstrates that the radio technique can provide a cost-effective extension of existing air-shower arrays. Operating in the frequency range of 30-80 MHz, Tunka-Rex is limited by the galactic background, and suffers from the local radio interferences. We investigate the possibilities of the improving of measured data using different approaches, particularly, the multivariate background suppression is considered, as well as improved likelihood fit of the lateral distribution of amplitudes.

We investigate features of the lateral distribution function (LDF) of the radio signal emitted by cosmic ray air-showers with primary energies $> 0.1$~EeV and its connection to air-shower parameters such as energy and shower maximum using CoREAS simulations made for the configuration of the Tunka-Rex antenna array. Taking into account all significant contributions to the total radio emission, such as by the geomagnetic effect, the charge excess, and the atmospheric refraction we parameterize the radio LDF. This parameterization is two-dimensional and has several free parameters. The large number of free parameters is not suitable for experiments of sparse arrays operating at low SNR (signal-to-noise ratios). Thus, exploiting symmetries, we decrease the number of free parameters and reduce the LDF to a simple one-dimensional function. The remaining parameters can be fit with a small number of points, i.e. as few as the signal from three antennas above detection threshold. Finally, we present a method for the reconstruction of air-shower parameters, in particular, energy and $X_{\mathrm{max}}$ (shower maximum), which can be reached with a theoretical accuracy of better than 15\% and 30~g/cm$^2$, respectively.