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.

During the past decade, semi-conductor photo sensors have replaced photo multiplier tubes in numerous applications featuring single-photon resolution, insensitivity to magnetic fields, higher robustness and enhanced photo detection efficiency at lower operation voltage and lower costs. This paper presents a measurement of gain, optical crosstalk and dark count rate of 127 Silicon Photo Multipliers (SiPMs) of type SensL MicroFJ-60035-TSV at nominal bias voltage. In SiPMs, optical crosstalk can fake higher multiplicities for single-photon hits. This paper demonstrates that optically coupled light guides reduce this probability significantly. The measured sensors offer a comparably small variation of their breakdown voltage with temperature of only $21.5\,\mathrm{mV/K}$. The manufacturer specifies a limited range of only $\pm 250\,\mathrm{mV/K}$ for their nominal breakdown voltage. This should allow for operation with small systematic errors without major efforts on temperature-based bias voltage compensation and without calibration of their individual breakdown voltages. This paper proves that the gain distribution of the measured sensors is consistent with the data sheet values.

Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $\nu_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In August 2016, 1ES~2344+514 was detected with the ground-based $\gamma$-ray telescope FACT during a high $\gamma$-ray state, triggering multi-wavelength (MWL) observations. We studied the MWL light curves of 1ES~2344+514 during the 2016 flaring state, using data from radio to VHE $\gamma$ rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, \textit{Swift}-UVOT, \textit{Swift}-XRT, \textit{Fermi}-LAT, FACT and MAGIC. With simultaneous observations of the flare, we built the broadband SED and studied it in the framework of a leptonic and an hadronic model. The VHE $\gamma$-ray observations show a flux level of 55\% of the Crab Nebula flux above 300\,GeV, similar to the historical maximum of 1995. The combination of MAGIC and \textit{Fermi}-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The $\Gamma$ index of the intrinsic spectrum in the VHE $\gamma$-ray band is $2.04\pm0.12_{\rm stat}\pm0.15_{\rm sys}$. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to $10^{18}$\,Hz

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.

Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars. The temporal evolution of the most prominent and variable segments of the SED is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. A narrow feature at ~3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux ($>0.3$ keV) measured during the entire Swift mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3$\sigma$. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4$\sigma$, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3$\sigma$. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: a) a pileup in the EED due to stochastic acceleration; b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and c) an emission from an IC pair cascade.

The IceCube Neutrino Observatory has revealed the existence of sources of high-energy astrophysical neutrinos. However, identification of the sources is challenging because astrophysical neutrinos are difficult to separate from the background of atmospheric neutrinos produced in cosmic-ray-induced particle cascades in the atmosphere. The efficient detection of air showers in coincidence with detected neutrinos can greatly reduce those backgrounds and increase the sensitivity of neutrino telescopes. Imaging Air Cherenkov Telescopes (IACTs) are sensitive to gamma-ray-induced (and cosmic-ray-induced) air showers in the 50 GeV to 50 TeV range, and can therefore be used as background-identifiers for neutrino observatories. This paper describes the feasibility of an array of small scale, wide field-of-view, cost-effective IACTs as an air shower veto for neutrino astronomy. A surface array of 250 to 750 telescopes would significantly improve the performance of a cubic kilometer-scale detector like IceCube, at a cost of a few percent of the original investment. The number of telescopes in the array can be optimized based on astronomical and geometrical considerations.

In this paper we describe the first results of a compact imaging air-Cherenkov telescope, IceAct, operating in coincidence with the IceCube Neutrino Observatory (IceCube) at the geographic South Pole. An array of IceAct telescopes (referred to as the IceAct project) is under consideration as part of the IceCube-Gen2 extension to IceCube. Surface detectors in general will be a powerful tool in IceCube-Gen2 for distinguishing astrophysical neutrinos from the dominant backgrounds of cosmic-ray induced atmospheric muons and neutrinos: the IceTop array is already in place as part of IceCube, but has a high energy threshold. Although the duty cycle will be lower for the IceAct telescopes than the present IceTop tanks, the IceAct telescopes may prove to be more effective at lowering the detection threshold for air showers. Additionally, small imaging air-Cherenkov telescopes in combination with IceTop, the deep IceCube detector or other future detector systems might improve measurements of the composition of the cosmic ray energy spectrum. In this paper we present measurements of a first 7-pixel imaging air Cherenkov telescope demonstrator, proving the capability of this technology to measure air showers at the South Pole in coincidence with IceTop and the deep IceCube detector.

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.

Mrk 421 is a high-synchrotron-peaked blazar featuring bright and persistent GeV and TeV emission. We use the longest and densest ongoing unbiased observing campaign obtained at TeV and GeV energies during 5.5 years with the FACT telescope and the Fermi-LAT detector. The contemporaneous multi-wavelength observations were used to characterize the variability of the source and to constrain the underlying physical mechanisms. We study and correlate light curves obtained by nine different instruments from radio to gamma rays and found two significant results. The TeV and X-ray light curves are very well correlated with lag, if any, shorter than a day. The GeV light curve varies independently and accurately leads the variations observed at long wavelengths, in particular in the radio band. We find that the observations match the predictions of leptonic models and suggest that the physical conditions vary along the jet, when the emitting region moves outwards.

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

Mrk 421 is a high-synchrotron-peaked blazar featuring bright and persistent GeV and TeV emission. We use multi-wavelength light curves of Mrk 421 spanning 5.5 years with FACT (TeV) and Fermi LAT (GeV) in the gamma rays, Swift BAT, Swift XRT and MAXI in the X-rays, together with optical and radio data and investigate the physical processes driving the emission and variability. Observations by FACT are continuous and not triggered, so the source was found in a wide range of flux states and more than 30 flares were identified from X-rays to TeV. The light curves in TeV and X-rays feature very similar flares with rise and decay times of a few days and zero lag, characteristic for electron processes. At least two parameters per flare, the amplitude and the cut-off energy, are required to explain the observed variability. In addition, the GeV light curve leads and is strongly correlated with the optical and radio light curves as expected from SSC emitting shock propagating in a conical jet.

We describe plans for the development of the Southern Wide-field Gamma-ray Observatory (SWGO), a next-generation instrument with sensitivity to the very-high-energy (VHE) band to be constructed in the Southern Hemisphere. SWGO will provide wide-field coverage of a large portion of the southern sky, effectively complementing current and future instruments in the global multi-messenger effort to understand extreme astrophysical phenomena throughout the universe. A detailed description of science topics addressed by SWGO is available in the science case white paper [1]. The development of SWGO will draw on extensive experience within the community in designing, constructing, and successfully operating wide-field instruments using observations of extensive air showers. The detector will consist of a compact inner array of particle detection units surrounded by a sparser outer array. A key advantage of the design of SWGO is that it can be constructed using current, already proven technology. We estimate a construction cost of 54M USD and a cost of 7.5M USD for 5 years of operation, with an anticipated US contribution of 20M USD ensuring that the US will be a driving force for the SWGO effort. The recently formed SWGO collaboration will conduct site selection and detector optimization studies prior to construction, with full operations foreseen to begin in 2026. Throughout this document, references to science white papers submitted to the Astro2020 Decadal Survey with particular relevance to the key science goals of SWGO, which include unveiling Galactic particle accelerators [2-10], exploring the dynamic universe [11-21], and probing physics beyond the Standard Model [22-25], are highlighted in red boldface.

The next generation magnetic spectrometer in space, AMS-100, is designed to have a geometrical acceptance of $100\,\mathrm{m}^2\,\mathrm{sr}$ and to be operated for at least ten years at the Sun-Earth Lagrange Point 2. Compared to existing experiments, it will improve the sensitivity for the observation of new phenomena in cosmic rays, and in particular in cosmic antimatter, by at least a factor of 1000. The magnet design is based on high temperature superconductor tapes, which allow the construction of a thin solenoid with a homogeneous magnetic field of 1 Tesla inside. The inner volume is instrumented with a silicon tracker reaching a maximum detectable rigidity of 100 TV and a calorimeter system that is 70 radiation lengths deep, equivalent to four nuclear interaction lengths, which extends the energy reach for cosmic-ray nuclei up to the PeV scale, i.e. beyond the cosmic-ray knee. Covering most of the sky continuously, AMS-100 will detect high-energy gamma rays in the calorimeter system and by pair conversion in the thin solenoid, reconstructed with excellent angular resolution in the silicon tracker.

Propagation of ultra-high energy photons in the solar magnetosphere gives rise to cascades comprising thousands of photons. We study the cascade development using Monte Carlo simulations and find that the photons in the cascades are spatially extended over hundreds of kilometers as they arrive at the top of the Earth's atmosphere. We compare results from simulations which use two models of the solar magnetic field, and show that although signatures of such cascades are different for the models used, for practical detection purpose in the ground-based detectors, they are similar.

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.

We outline the science motivation for SGSO, the Southern Gamma-Ray Survey Observatory. SGSO will be a next-generation wide field-of-view gamma-ray survey instrument, sensitive to gamma-rays in the energy range from 100 GeV to hundreds of TeV. Its science topics include unveiling galactic and extragalactic particle accelerators, monitoring the transient sky at very high energies, probing particle physics beyond the Standard Model, and the characterization of the cosmic ray flux. SGSO will consist of an air shower detector array, located in South America. Due to its location and large field of view, SGSO will be complementary to other current and planned gamma-ray observatories such as HAWC, LHAASO, and CTA.

The rate of extensive air-showers observed with imaging air-Cherenkov telescopes is zenith angle dependent. This effect originates from the increasing geometrical distance of the observed shower to the telescope with increasing zenith distance. This paper investigates how this alters the observed image and how this affects the trigger rate as a function of zenith angle. The discussed effects include the change of Cherenkov light yield, of absorption in the atmosphere, of photon density at the aperture and of the image size at the focal plane of the telescope. Based on a simple model for the atmosphere and well-known first principles on the development of extensive air-showers, the zenith angle dependence is expressed analytically. The assumption that most light is emitted from the shower core and mathematical approximations allow to derive an analytical expression describing the zenith angle dependence well with only three free parameters which are directly linked with the underlying physics. This suggests further investigations about how these fit parameters are linked to the properties of the atmosphere and the instrument. Using data published by the First G-APD Cherenkov Telescope, a good match of the fit functions with the data is obtained. For the trigger rate of cosmic rays, the obtained parameters are consistent with the naive expectation.

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.

Aims: To investigate the extension of the very-high-energy spectral tail of the Crab pulsar at energies above 400 GeV. Methods: We analyzed $\sim$320 hours of good quality data of Crab with the MAGIC telescope, obtained from February 2007 until April 2014. Results: We report the most energetic pulsed emission ever detected from the Crab pulsar reaching up to 1.5 TeV. The pulse profile shows two narrow peaks synchronized with the ones measured in the GeV energy range. The spectra of the two peaks follow two different power-law functions from 70 GeV up to 1.5 TeV and connect smoothly with the spectra measured above 10 GeV by the Large Area Telescope (LAT) on board of the Fermi satellite. When making a joint fit of the LAT and MAGIC data, above 10 GeV, the photon indices of the spectra differ by 0.5$\pm$0.1. Conclusions: We measured with the MAGIC telescopes the most energetic pulsed photons from a pulsar to date. Such TeV pulsed photons require a parent population of electrons with a Lorentz factor of at least $5\times 10^6$. These results strongly suggest IC scattering off low energy photons as the emission mechanism and a gamma-ray production region in the vicinity of the light cylinder.

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.

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

A general purpose instrument for imaging of Cherenkov light or fluorescence light emitted by extensive air showers is presented. Its refractive optics allows for a compact and light-weight design with a wide field-of-view of 12{\deg}. The optical system features a 0.5 m diameter Fresnel lens and a camera with 61 pixels composed of Winston cones and large-sized 6x6 mm photo sensors. As photo sensors, semi conductor light sensors (SiPMs) are utilized. The camera provides a high photon detection efficiency together with robust operation. The enclosed optics permit operation in regions of harsh environmental conditions. The low price of the telescope allows the production of a large number of telescopes and the application of the instrument in various projects, such as FAMOUS for the Pierre Auger Observatory, HAWC's Eye for HAWC or IceAct for IceCube. In this paper the novel design of this telescope and first measurements are presented.

An in-situ calibration of a logarithmic periodic dipole antenna with a frequency coverage of 30 MHz to 80 MHz is performed. Such antennas are part of a radio station system used for detection of cosmic ray induced air showers at the Engineering Radio Array of the Pierre Auger Observatory, the so-called Auger Engineering Radio Array (AERA). The directional and frequency characteristics of the broadband antenna are investigated using a remotely piloted aircraft (RPA) carrying a small transmitting antenna. The antenna sensitivity is described by the vector effective length relating the measured voltage with the electric-field components perpendicular to the incoming signal direction. The horizontal and meridional components are determined with an overall uncertainty of 7.4^{+0.9}_{-0.3} % and 10.3^{+2.8}_{-1.7} % respectively. The measurement is used to correct a simulated response of the frequency and directional response of the antenna. In addition, the influence of the ground conductivity and permittivity on the antenna response is simulated. Both have a negligible influence given the ground conditions measured at the detector site. The overall uncertainties of the vector effective length components result in an uncertainty of 8.8^{+2.1}_{-1.3} % in the square root of the energy fluence for incoming signal directions with zenith angles smaller than 60{\deg}.

Although the photon structure is most efficiently studied with the accelerator instruments, there is also a scientifically complementary potential in investigations on photons produced in the outer space. This potential is already being explored with gamma ray telescopes, ultra-high energy cosmic ray observatories and, since very recently, by the Cosmic-Ray Extremely Distributed Observatory (CREDO). Unlike the former instruments focused on detection of single photons, CREDO aims at the detection of cascades (ensembles) of photons originating even at astrophysical distances. If at least a part of such a cascade reaches Earth, it might produce a unique pattern composed of a number of air showers observable by an appropriately dense array of standard detectors. If the energies of air showers constituting the pattern are relatively low and if the typical distances between the neighbors are large, the ensemble character of the whole phenomenon might remain uncovered, unless the CREDO strategy is implemented.

We present a combined fit of a simple astrophysical model of UHECR sources to both the energy spectrum and mass composition data measured by the Pierre Auger Observatory. The fit has been performed for energies above $5 \cdot 10^{18}$ eV, i.e.~the region of the all-particle spectrum above the so-called "ankle" feature. The astrophysical model we adopted consists of identical sources uniformly distributed in a comoving volume, where nuclei are accelerated through a rigidity-dependent mechanism. The fit results suggest sources characterized by relatively low maximum injection energies, hard spectra and heavy chemical composition. We also show that uncertainties about physical quantities relevant to UHECR propagation and shower development have a non-negligible impact on the fit results.

A new analysis of the dataset from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above 20 EeV with zenith angles up to 80 deg recorded before 2017 April 30. Sky models have been created for two distinct populations of extragalactic gamma-ray emitters: active galactic nuclei from the second catalog of hard Fermi-LAT sources (2FHL) and starburst galaxies from a sample that was examined with Fermi-LAT. Flux-limited samples, which include all types of galaxies from the Swift-BAT and 2MASS surveys, have been investigated for comparison. The sky model of cosmic-ray density constructed using each catalog has two free parameters, the fraction of events correlating with astrophysical objects and an angular scale characterizing the clustering of cosmic rays around extragalactic sources. A maximum-likelihood ratio test is used to evaluate the best values of these parameters and to quantify the strength of each model by contrast with isotropy. It is found that the starburst model fits the data better than the hypothesis of isotropy with a statistical significance of 4.0 sigma, the highest value of the test statistic being for energies above 39 EeV. The three alternative models are favored against isotropy with 2.7-3.2 sigma significance. The origin of the indicated deviation from isotropy is examined and prospects for more sensitive future studies are discussed.

Any considerations on propagation of particles through the Universe must involve particle interactions: processes leading to production of particle cascades. While one expects existence of such cascades, the state of the art cosmic-ray research is oriented purely on a detection of single particles, gamma rays or associated extensive air showers. The natural extension of the cosmic-ray research with the studies on ensembles of particles and air showers is being proposed by the CREDO Collaboration. Within the CREDO strategy the focus is put on generalized super-preshowers (SPS): spatially and/or temporally extended cascades of particles originated above the Earth atmosphere, possibly even at astrophysical distances. With CREDO we want to find out whether SPS can be at least partially observed by a network of terrestrial and/or satellite detectors receiving primary or secondary cosmic-ray signal. This paper addresses electromagnetic SPS, e.g. initiated by VHE photons interacting with the cosmic microwave background, and the SPS signatures that can be seen by gamma-ray telescopes, exploring the exampleof Cherenkov Telescope Array. The energy spectrum of secondary electrons and photons in an electromagnetic super-preshower might be extended over awide range of energy, down to TeV or even lower, as it is evident from the simulation results. This means that electromagnetic showers induced by such particles in the Earth atmosphere could be observed by imaging atmospheric Cherenkov telescopes. We present preliminary results from the study of response of the Cherenkov Telescope Array to SPS events, including the analysis of the simulated shower images on the camera focal plane and implementedgeneric reconstruction chains based on the Hillas parameters.

The Cosmic-Ray Extremely Distributed Observatory (CREDO) is an infrastructure for global analysis of extremely extended cosmic-ray phenomena, so-called super-preshowers, beyond the capabilities of existing, discrete, detectors and observatories. To date cosmic-ray research has been focused on detecting single air showers, while the search for ensembles of cosmic-ray events induced by super-preshowers is a scientific terra incognita - CREDO explores this uncharted realm. Positive detection of super-preshowers would have an impact on ultra-high energy astrophysics, cosmology and the physics of fundamental particle interactions as they can theoretically be formed within both classical (photon-photon interactions) and exotic (Super Heavy Dark Matter particle decay and interaction) scenarios. Some super-preshowers are predicted to have a significant spatial extent - a unique signature only detectable with the existing cosmic-ray infrastructure taken as a global network. An obvious, although yet unprobed, super-preshower 'detection limit' would be located somewhere between an air shower, induced by a super-preshower composed of tightly collimated particles, and a super-preshower composed of particles spread so widely that only few of them can reach the Earth. CREDO will probe this detection limit, leading to either an observation of an as yet unseen physical phenomenon, or the setting upper limits to the existence of large extraterrestrial cascades which would constrain fundamental physics models. While CREDO's focus is on testing physics at energies close to the Grand Unified Theories range, the broader phenomena are expected to be composed of particles with energies ranging from GeV to ZeV. This motivates our advertising of this concept across the astroparticle physics community.