Showing 1–50 of 98 results for author: Winchen, T

CRPropa 3.2: a public framework for high-energy astroparticle simulations

Authors: Sophie Aerdker, Rafael Alves Batista, Julia Becker Tjus, Julien Dörner, Andrej Dundovic, Björn Eichmann, Antonius Frie, Christopher Heiter, Mario Hoerbe, Karl-Heinz Kampert, Lukas Merten, Gero Müller, Patrick Reichherzer, Simone Rossoni, Andrey Saveliev, Leander Schlegel, Günter Sigl, Arjen van Vliet, Tobias Winchen

Abstract: CRPropa is a Monte Carlo framework for simulating the propagation of (ultra-) high-energy particles in the Universe, including cosmic rays, gamma rays, electrons, and neutrinos. It covers energies from ZeV down to GeV for gamma rays and electrons, and TeV for cosmic rays and neutrinos, supporting various astrophysical environments such as the surroundings of astrophysical sources, galactic, and ex… ▽ More CRPropa is a Monte Carlo framework for simulating the propagation of (ultra-) high-energy particles in the Universe, including cosmic rays, gamma rays, electrons, and neutrinos. It covers energies from ZeV down to GeV for gamma rays and electrons, and TeV for cosmic rays and neutrinos, supporting various astrophysical environments such as the surroundings of astrophysical sources, galactic, and extragalactic environments. The newest version, CRPropa 3.2, represents a significant leap forward towards a universal multi-messenger framework, opening up the possibility for many more astrophysical applications. This includes extensions to simulate cosmic-ray acceleration and particle interactions within astrophysical source environments, a full Monte Carlo treatment of electromagnetic cascades, improved ensemble-averaged Galactic propagation, significant performance improvements for cosmic-ray tracking through magnetic fields, and a user-friendly implementation of custom photon fields, among many more enhancements. This contribution will give an overview of the new features and present several applications to cosmic-ray and gamma-ray propagation. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: PoS 444 (38th ICRC) (2023) 1471

CRPropa 3.2 -- an advanced framework for high-energy particle propagation in extragalactic and galactic spaces

Authors: Rafael Alves Batista, Julia Becker Tjus, Julien Dörner, Andrej Dundovic, Björn Eichmann, Antonius Frie, Christopher Heiter, Mario R. Hoerbe, Karl-Heinz Kampert, Lukas Merten, Gero Müller, Patrick Reichherzer, Andrey Saveliev, Leander Schlegel, Günter Sigl, Arjen van Vliet, Tobias Winchen

Abstract: The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, largely due to the inflow of data collected by large-scale cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 is a Monte Carlo c… ▽ More The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, largely due to the inflow of data collected by large-scale cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 is a Monte Carlo code for simulating the propagation of high-energy particles in the Universe. This version represents a major leap forward, significantly expanding the simulation framework and opening up the possibility for many more astrophysical applications. This includes, among others: efficient simulation of high-energy particles in diffusion-dominated domains, self-consistent and fast modelling of electromagnetic cascades with an extended set of channels for photon production, and studies of cosmic-ray diffusion tensors based on updated coherent and turbulent magnetic-field models. Furthermore, several technical updates and improvements are introduced with the new version, such as: enhanced interpolation, targeted emission of sources, and a new propagation algorithm (Boris push). The detailed description of all novel features is accompanied by a discussion and a selected number of example applications. △ Less

Submitted 29 July, 2022; originally announced August 2022.

Comments: 32 pages, 9 figures

Journal ref: JCAP 09 (2022) 035

A burst storm from the repeating FRB 20200120E in an M81 globular cluster

Authors: K. Nimmo, J. W. T. Hessels, M. P. Snelders, R. Karuppusamy, D. M. Hewitt, F. Kirsten, B. Marcote, U. Bach, A. Bansod, E. D. Barr, J. Behrend, V. Bezrukovs, S. Buttaccio, R. Feiler, M. P. Gawroński, M. Lindqvist, A. Orbidans, W. Puchalska, N. Wang, T. Winchen, P. Wolak, J. Wu, J. Yuan

Abstract: The repeating fast radio burst (FRB) source FRB 20200120E is exceptional because of its proximity and association with a globular cluster. Here we report $60$ bursts detected with the Effelsberg telescope at 1.4 GHz. We observe large variations in the burst rate, and report the first FRB 20200120E `burst storm', where the source suddenly became active and 53 bursts (fluence $\geq 0.04$ Jy ms) occu… ▽ More The repeating fast radio burst (FRB) source FRB 20200120E is exceptional because of its proximity and association with a globular cluster. Here we report $60$ bursts detected with the Effelsberg telescope at 1.4 GHz. We observe large variations in the burst rate, and report the first FRB 20200120E `burst storm', where the source suddenly became active and 53 bursts (fluence $\geq 0.04$ Jy ms) occurred within only 40 minutes. We find no strict periodicity in the burst arrival times, nor any evidence for periodicity in the source's activity between observations. The burst storm shows a steep energy distribution (power-law index $α= 2.39\pm0.12$) and a bi-modal wait-time distribution, with log-normal means of 0.94$^{+0.07}_{-0.06}$ s and 23.61$^{+3.06}_{-2.71}$ s. We attribute these wait-time distribution peaks to a characteristic event timescale and pseudo-Poisson burst rate, respectively. The secondary wait-time peak at $\sim1$ s is $\sim50\times$ longer than the $\sim24$ ms timescale seen for both FRB 20121102A and FRB 20201124A -- potentially indicating a larger emission region, or slower burst propagation. FRB 20200120E shows order-of-magnitude lower burst durations and luminosities compared with FRB 20121102A and FRB 20201124A. Lastly, in contrast to FRB 20121102A, which has observed dispersion measure (DM) variations of $Δ{\rm DM} >1$ pc cm$^{-3}$ on month-to-year timescales, we determine that FRB 20200120E's DM has remained stable ($Δ{\rm DM} <0.15$ pc cm$^{-3}$) over $>10$ months. Overall, the observational characteristics of FRB 20200120E deviate quantitatively from other active repeaters, but it is unclear whether it is qualitatively a different type of source. △ Less

Submitted 23 January, 2023; v1 submitted 8 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS

CRPropa 3.2: a framework for high-energy astroparticle propagation

Authors: Rafael Alves Batista, Julia Becker Tjus, Julien Dörner, Andrej Dundovic, Björn Eichmann, Antonius Frie, Christopher Heiter, Mario R. Hoerbe, Karl-Heinz Kampert, Lukas Merten, Gero Müller, Patrick Reichherzer, Andrey Saveliev, Leander Schlegel, Günter Sigl, Arjen van Vliet, Tobias Winchen

Abstract: The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, in large part owing to the inflow of high-quality data collected by present cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 i… ▽ More The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, in large part owing to the inflow of high-quality data collected by present cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 is a Monte Carlo code for simulating the propagation of high-energy particles in the Universe. This new version represents a step further towards a more complete simulation framework for multimessenger studies. Some of the new developments include: cosmic-ray acceleration, support for particle interactions within astrophysical sources, full Monte Carlo treatment of electromagnetic cascades, improved ensemble-averaged Galactic propagation, and a number of technical enhancements. Here we present some of these novel features and some applications to gamma- and cosmic-ray propagation. △ Less

Submitted 4 July, 2021; originally announced July 2021.

Comments: 8 pages, 4 figures; Proceedings of the 37th International Cosmic Ray Conference (ICRC2021)

Depth of shower maximum and mass composition of cosmic rays from 50 PeV to 2 EeV measured with the LOFAR radio telescope

Authors: A. Corstanje, S. Buitink, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, G. K. Krampah, P. Mitra, K. Mulrey, A. Nelles, H. Pandya, J. P. Rachen, O. Scholten, S. ter Veen, S. Thoudam, G. Trinh, T. Winchen

Abstract: We present an updated cosmic-ray mass composition analysis in the energy range $10^{16.8}$ to $10^{18.3}$ eV from 334 air showers measured with the LOFAR radio telescope, and selected for minimal bias. In this energy range, the origin of cosmic rays is expected to shift from galactic to extragalactic sources. The analysis is based on an improved method to infer the depth of maximum $X_{\rm max}$ o… ▽ More We present an updated cosmic-ray mass composition analysis in the energy range $10^{16.8}$ to $10^{18.3}$ eV from 334 air showers measured with the LOFAR radio telescope, and selected for minimal bias. In this energy range, the origin of cosmic rays is expected to shift from galactic to extragalactic sources. The analysis is based on an improved method to infer the depth of maximum $X_{\rm max}$ of extensive air showers from radio measurements and air shower simulations. We show results of the average and standard deviation of $X_{\rm max}$ versus primary energy, and analyze the $X_{\rm max}$-dataset at distribution level to estimate the cosmic ray mass composition. Our approach uses an unbinned maximum likelihood analysis, making use of existing parametrizations of $X_{\rm max}$-distributions per element. The analysis has been repeated for three main models of hadronic interactions. Results are consistent with a significant light-mass fraction, at best fit $23$ to $39$ $\%$ protons plus helium, depending on the choice of hadronic interaction model. The fraction of intermediate-mass nuclei dominates. This confirms earlier results from LOFAR, with systematic uncertainties on $X_{\rm max}$ now lowered to 7 to $9$ $\mathrm{g/cm^2}$. We find agreement in mass composition compared to results from Pierre Auger Observatory, within statistical and systematic uncertainties. However, in line with earlier LOFAR results, we find a slightly lower average $X_{\rm max}$. The values are in tension with those found at Pierre Auger Observatory, but agree with results from other cosmic ray observatories based in the Northern hemisphere. △ Less

Submitted 2 May, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: 24 pages, 14 figures. Accepted for publication in Phys. Rev. D

Journal ref: Phys. Rev. D 103, 102006 (2021)

The FRAM robotic telescope for atmospheric monitoring at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker , et al. (355 additional authors not shown)

Abstract: FRAM (F/Photometric Robotic Atmospheric Monitor) is a robotic telescope operated at the Pierre Auger Observatory in Argentina for the purposes of atmospheric monitoring using stellar photometry. As a passive system which does not produce any light that could interfere with the observations of the fluorescence telescopes of the observatory, it complements the active monitoring systems that use lase… ▽ More FRAM (F/Photometric Robotic Atmospheric Monitor) is a robotic telescope operated at the Pierre Auger Observatory in Argentina for the purposes of atmospheric monitoring using stellar photometry. As a passive system which does not produce any light that could interfere with the observations of the fluorescence telescopes of the observatory, it complements the active monitoring systems that use lasers. We discuss the applications of stellar photometry for atmospheric monitoring at optical observatories in general and the particular modes of operation employed by the Auger FRAM. We describe in detail the technical aspects of FRAM, the hardware and software requirements for a successful operation of a robotic telescope for such a purpose and their implementation within the FRAM system. △ Less

Submitted 26 July, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: Published version, 27 pages, 8 figures

Journal ref: JINST 16 P06027 (2021)

Deep-Learning based Reconstruction of the Shower Maximum $X_{\mathrm{max}}$ using the Water-Cherenkov Detectors of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker , et al. (348 additional authors not shown)

Abstract: The atmospheric depth of the air shower maximum $X_{\mathrm{max}}$ is an observable commonly used for the determination of the nuclear mass composition of ultra-high energy cosmic rays. Direct measurements of $X_{\mathrm{max}}$ are performed using observations of the longitudinal shower development with fluorescence telescopes. At the same time, several methods have been proposed for an indirect e… ▽ More The atmospheric depth of the air shower maximum $X_{\mathrm{max}}$ is an observable commonly used for the determination of the nuclear mass composition of ultra-high energy cosmic rays. Direct measurements of $X_{\mathrm{max}}$ are performed using observations of the longitudinal shower development with fluorescence telescopes. At the same time, several methods have been proposed for an indirect estimation of $X_{\mathrm{max}}$ from the characteristics of the shower particles registered with surface detector arrays. In this paper, we present a deep neural network (DNN) for the estimation of $X_{\mathrm{max}}$. The reconstruction relies on the signals induced by shower particles in the ground based water-Cherenkov detectors of the Pierre Auger Observatory. The network architecture features recurrent long short-term memory layers to process the temporal structure of signals and hexagonal convolutions to exploit the symmetry of the surface detector array. We evaluate the performance of the network using air showers simulated with three different hadronic interaction models. Thereafter, we account for long-term detector effects and calibrate the reconstructed $X_{\mathrm{max}}$ using fluorescence measurements. Finally, we show that the event-by-event resolution in the reconstruction of the shower maximum improves with increasing shower energy and reaches less than $25~\mathrm{g/cm^{2}}$ at energies above $2\times 10^{19}~\mathrm{eV}$. △ Less

Submitted 27 July, 2021; v1 submitted 8 January, 2021; originally announced January 2021.

Comments: Published version, 29 pages, 12 figures

Report number: FERMILAB-PUB-21-084-AD-AE-SCD-TD

Journal ref: JINST 16 P07019 (2021)

Calibration of the underground muon detector of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker , et al. (348 additional authors not shown)

Abstract: To obtain direct measurements of the muon content of extensive air showers with energy above $10^{16.5}$ eV, the Pierre Auger Observatory is currently being equipped with an underground muon detector (UMD), consisting of 219 10 $\mathrm{m^2}$-modules, each segmented into 64 scintillators coupled to silicon photomultipliers (SiPMs). Direct access to the shower muon content allows for the study of b… ▽ More To obtain direct measurements of the muon content of extensive air showers with energy above $10^{16.5}$ eV, the Pierre Auger Observatory is currently being equipped with an underground muon detector (UMD), consisting of 219 10 $\mathrm{m^2}$-modules, each segmented into 64 scintillators coupled to silicon photomultipliers (SiPMs). Direct access to the shower muon content allows for the study of both of the composition of primary cosmic rays and of high-energy hadronic interactions in the forward direction. As the muon density can vary between tens of muons per m$^2$ close to the intersection of the shower axis with the ground to much less than one per m$^2$ when far away, the necessary broad dynamic range is achieved by the simultaneous implementation of two acquisition modes in the read-out electronics: the binary mode, tuned to count single muons, and the ADC mode, suited to measure a high number of them. In this work, we present the end-to-end calibration of the muon detector modules: first, the SiPMs are calibrated by means of the binary channel, and then, the ADC channel is calibrated using atmospheric muons, detected in parallel to the shower data acquisition. The laboratory and field measurements performed to develop the implementation of the full calibration chain of both binary and ADC channels are presented and discussed. The calibration procedure is reliable to work with the high amount of channels in the UMD, which will be operated continuously, in changing environmental conditions, for several years. △ Less

Submitted 14 April, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: Published version, 25 pages, 12 figures, 0 tables

Report number: FERMILAB-PUB-20-678-AD-E-TD

Journal ref: JINST (2021) 16 P04003

Design, upgrade and characterization of the silicon photomultiplier front-end for the AMIGA detector at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker, J. A. Bellido , et al. (335 additional authors not shown)

Abstract: AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to complement the study of ultra-high-energy cosmic rays (UHECR) by measuring the muon content of extensive air showers (EAS). It consists of an array of 61 water Cherenkov detectors on a denser spacing in combination with underground scintillation detectors used for muon density measurement. Each det… ▽ More AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to complement the study of ultra-high-energy cosmic rays (UHECR) by measuring the muon content of extensive air showers (EAS). It consists of an array of 61 water Cherenkov detectors on a denser spacing in combination with underground scintillation detectors used for muon density measurement. Each detector is composed of three scintillation modules, with 10 m$^2$ detection area per module, buried at 2.3 m depth, resulting in a total detection area of 30 m$^2$. Silicon photomultiplier sensors (SiPM) measure the amount of scintillation light generated by charged particles traversing the modules. In this paper, the design of the front-end electronics to process the signals of those SiPMs and test results from the laboratory and from the Pierre Auger Observatory are described. Compared to our previous prototype, the new electronics shows a higher performance, higher efficiency and lower power consumption, and it has a new acquisition system with increased dynamic range that allows measurements closer to the shower core. The new acquisition system is based on the measurement of the total charge signal that the muonic component of the cosmic ray shower generates in the detector. △ Less

Submitted 25 January, 2021; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 40 pages, 33 figures

Report number: FERMILAB-PUB-20-605-AD-E-TD

Journal ref: JINST 16 (2021) P01026

A search for ultra high energy neutrinos from TXS 0506+056 using the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker, J. A. Bellido , et al. (342 additional authors not shown)

Abstract: Results of a search for ultra-high-energy neutrinos with the Pierre Auger Observatory from the direction of the blazar TXS 0506+056 are presented. They were obtained as part of the follow-up that stemmed from the detection of high-energy neutrinos and gamma rays with IceCube, \textit{Fermi}-LAT, MAGIC, and other detectors of electromagnetic radiation in several bands. The Pierre Auger Observatory… ▽ More Results of a search for ultra-high-energy neutrinos with the Pierre Auger Observatory from the direction of the blazar TXS 0506+056 are presented. They were obtained as part of the follow-up that stemmed from the detection of high-energy neutrinos and gamma rays with IceCube, \textit{Fermi}-LAT, MAGIC, and other detectors of electromagnetic radiation in several bands. The Pierre Auger Observatory is sensitive to neutrinos in the energy range from 100 PeV to 100 EeV and in the zenith angle range from $θ=60^\circ$ to $θ=95^\circ$, where the zenith angle is measured from the vertical direction. No neutrinos from the direction of TXS 0506+056 have been found. The results were analyzed in three periods: One of 6 months around the detection of IceCube-170922A, coinciding with a flare period of TXS 0506+056, a second one of 110 days during which the IceCube collaboration found an excess of 13 neutrinos from a direction compatible with TXS 0506+056, and a third one from 2004 January 1 up to 2018 August 31, over which the Pierre Auger Observatory has been taking data. The sensitivity of the Observatory is addressed for different spectral indices by considering the fluxes that would induce a single expected event during the observation period. For indices compatible with those measured by the IceCube collaboration the expected number of neutrinos at the Observatory is well-below one. Spectral indices as hard as 1.5 would have to apply in this energy range to expect a single event to have been detected. △ Less

Submitted 21 October, 2020; originally announced October 2020.

Comments: Published version, 4 figures

Journal ref: Published in The Astrophysical Journal, 902:105, 2020

Features of the energy spectrum of cosmic rays above $2.5{\times} 10^{18}$ eV using the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker , et al. (358 additional authors not shown)

Abstract: We report a measurement of the energy spectrum of cosmic rays above $2.5{\times} 10^{18}$ eV based on $215,030$ events. New results are presented: at about $1.3{\times} 10^{19}$ eV, the spectral index changes from $2.51 \pm 0.03 \textrm{ (stat.)} \pm 0.05 \textrm{ (sys.)}$ to $3.05 \pm 0.05 \textrm{ (stat.)}\pm 0.10\textrm{ (sys.)}$, evolving to… ▽ More We report a measurement of the energy spectrum of cosmic rays above $2.5{\times} 10^{18}$ eV based on $215,030$ events. New results are presented: at about $1.3{\times} 10^{19}$ eV, the spectral index changes from $2.51 \pm 0.03 \textrm{ (stat.)} \pm 0.05 \textrm{ (sys.)}$ to $3.05 \pm 0.05 \textrm{ (stat.)}\pm 0.10\textrm{ (sys.)}$, evolving to $5.1\pm0.3\textrm{ (stat.)} \pm 0.1\textrm{ (sys.)}$ beyond $5{\times} 10^{19}$ eV, while no significant dependence of spectral features on the declination is seen in the accessible range. These features of the spectrum can be reproduced in models with energy-dependent mass composition. The energy density in cosmic rays above $5{\times} 10^{18}$ eV is $(5.66 \pm 0.03 \textrm{ (stat.)} \pm 1.40 \textrm{ (sys.)} ) {\times} 10^{53}~$erg Mpc$^{-3}$. △ Less

Submitted 6 October, 2020; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: Published version; 2 figures, 1 table

Report number: FERMILAB-PUB-20-431-E-TD

Journal ref: Phys. Rev. Lett. 125, 121106 (2020)

Measurement of the cosmic-ray energy spectrum above $2.5{\times} 10^{18}$ eV using the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker , et al. (358 additional authors not shown)

Abstract: We report a measurement of the energy spectrum of cosmic rays for energies above $2.5 {\times} 10^{18}~$eV based on 215,030 events recorded with zenith angles below $60^\circ$. A key feature of the work is that the estimates of the energies are independent of assumptions about the unknown hadronic physics or of the primary mass composition. The measurement is the most precise made hitherto with th… ▽ More We report a measurement of the energy spectrum of cosmic rays for energies above $2.5 {\times} 10^{18}~$eV based on 215,030 events recorded with zenith angles below $60^\circ$. A key feature of the work is that the estimates of the energies are independent of assumptions about the unknown hadronic physics or of the primary mass composition. The measurement is the most precise made hitherto with the accumulated exposure being so large that the measurements of the flux are dominated by systematic uncertainties except at energies above $5 {\times} 10^{19}~$eV. The principal conclusions are: (1) The flattening of the spectrum near $5 {\times} 10^{18}~$eV, the so-called "ankle", is confirmed. (2) The steepening of the spectrum at around $5 {\times} 10^{19}~$eV is confirmed. (3) A new feature has been identified in the spectrum: in the region above the ankle the spectral index $γ$ of the particle flux ($\propto E^{-γ}$) changes from $2.51 \pm 0.03~{\rm (stat.)} \pm 0.05~{\rm (sys.)}$ to $3.05 \pm 0.05~{\rm (stat.)} \pm 0.10~{\rm (sys.)}$ before changing sharply to $5.1 \pm 0.3~{\rm (stat.)} \pm 0.1~{\rm (sys.)}$ above $5 {\times} 10^{19}~$eV. (4) No evidence for any dependence of the spectrum on declination has been found other than a mild excess from the Southern Hemisphere that is consistent with the anisotropy observed above $8 {\times} 10^{18}~$eV. △ Less

Submitted 6 October, 2020; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: Published version; 16 figures, 6 table

Report number: FERMILAB-PUB-20-432-E-TD

Journal ref: Phys. Rev. D 102, 062005 (2020)

Reconstruction of Events Recorded with the Surface Detector of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker , et al. (356 additional authors not shown)

Abstract: Cosmic rays arriving at Earth collide with the upper parts of the atmosphere, thereby inducing extensive air showers. When secondary particles from the cascade arrive at the ground, they are measured by surface detector arrays. We describe the methods applied to the measurements of the surface detector of the Pierre Auger Observatory to reconstruct events with zenith angles less than $60^\circ$ us… ▽ More Cosmic rays arriving at Earth collide with the upper parts of the atmosphere, thereby inducing extensive air showers. When secondary particles from the cascade arrive at the ground, they are measured by surface detector arrays. We describe the methods applied to the measurements of the surface detector of the Pierre Auger Observatory to reconstruct events with zenith angles less than $60^\circ$ using the timing and signal information recorded using the water-Cherenkov detector stations. In addition, we assess the accuracy of these methods in reconstructing the arrival directions of the primary cosmic ray particles and the sizes of the induced showers. △ Less

Submitted 5 November, 2020; v1 submitted 17 July, 2020; originally announced July 2020.

Comments: Published version, 36 pages, 14 figures

Report number: FERMILAB-PUB-20-487-AD-E-TD

Journal ref: JINST 15 (2020) P10021

Studies on the response of a water-Cherenkov detector of the Pierre Auger Observatory to atmospheric muons using an RPC hodoscope

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker , et al. (353 additional authors not shown)

Abstract: 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 knowled… ▽ More 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. △ Less

Submitted 9 September, 2020; v1 submitted 8 July, 2020; originally announced July 2020.

Comments: Published version, 25 pages, 9 figures

Report number: FERMILAB-PUB-20-285-PPD-TD

Journal ref: JINST 15 (2020) P09002

The initial stage of cloud lightning imaged in high-resolution

Authors: O. Scholten, B. M. Hare, J. Dwyer, C. Sterpka, I. Kolmašová, O. Santolík, R. Lán, L. Uhlíř, S. Buitink, A. Corstanje, H. Falcke, T. Huege, J. R. Hörandel, G. K. Krampah, P. Mitra, K. Mulrey, A. Nelles, H. Pandya, A. Pel, J. P. Rachen, T. N. G. Trinh, S. ter Veen, S. Thoudam, T. Winchen

Abstract: With LOFAR we have been able to image the development of lightning flashes with meter-scale accuracy and unprecedented detail. We discuss the primary steps behind our most recent lightning mapping method. To demonstrate the capabilities of our technique we show and interpret images of the first few milliseconds of two intra-cloud flashes. In all our flashes the negative leaders propagate in the ch… ▽ More With LOFAR we have been able to image the development of lightning flashes with meter-scale accuracy and unprecedented detail. We discuss the primary steps behind our most recent lightning mapping method. To demonstrate the capabilities of our technique we show and interpret images of the first few milliseconds of two intra-cloud flashes. In all our flashes the negative leaders propagate in the charge layer below the main negative charge. Among several interesting features we show that in about 2~ms after initiation the Primary Initial Leader triggers the formation of a multitude (more than ten) negative leaders in a rather confined area of the atmosphere. From these only one or two continue to propagate after about 30~ms to extend over kilometers horizontally while another may propagate back to the initiation point. We also show that normal negative leaders can transition into an initial-leader like state, potentially in the presence of strong electric fields. In addition, we show some initial breakdown pulses that occurred during the primary initial leader, and even during two "secondary" initial leaders that developed out of stepped leaders. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: Submitted to Journal of geophysics research: Atmospheres

Radio emission from negative lightning leader steps reveals inner meter-scale structure

Authors: B. M. Hare, O. Scholten, J. Dwyer, U. Ebert, S. Nijdam, A. Bonardi, S. Buitink, A. Corstanje, H. Falcke, T. Huege, J. R. Hörandel, G. K. Krampah, P. Mitra, K. Mulrey, B. Neijzen, A. Nelles, H. Pandya, J. P. Rachen, L. Rossetto, T. N. G. Trinh, S. ter Veen, T. Winchen

Abstract: We use the Low Frequency ARray (LOFAR) to probe the dynamics of the stepping process of negatively-charged plasma channels (negative leaders) in a lightning discharge. We observe that at each step of a leader, multiple pulses of VHF (30~--~80 MHz) radiation are emitted in short-duration bursts ($<10\ μ$s). This is evidence for streamer formation during corona flashes that occur with each leader st… ▽ More We use the Low Frequency ARray (LOFAR) to probe the dynamics of the stepping process of negatively-charged plasma channels (negative leaders) in a lightning discharge. We observe that at each step of a leader, multiple pulses of VHF (30~--~80 MHz) radiation are emitted in short-duration bursts ($<10\ μ$s). This is evidence for streamer formation during corona flashes that occur with each leader step, which has not been observed before in natural lightning and it could help explain X-ray emission from lightning leaders, as X-rays from laboratory leaders tend to be associated with corona flashes. Surprisingly we find that the stepping length is very similar to what was observed near the ground, however with a stepping time that is considerably larger, which as yet is not understood. These results will help to improve lightning propagation models, and eventually lightning protection models. △ Less

Submitted 7 July, 2020; originally announced July 2020.

Comments: 5 pages

Journal ref: Phys. Rev. Lett. 124, 105101 (2020)

Reconstructing air shower parameters with LOFAR using event specific GDAS atmospheres

Authors: P. Mitra, A. Bonardi, A. Corstanje, S. Buitink, G. K Krampah, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, K. Mulrey, A. Nelles, H. Pandya, J. P. Rachen, L. Rossetto, O. Scholten, S. ter Veen, T. N. G. Trinh, T. Winchen

Abstract: The limited knowledge of atmospheric parameters like humidity, pressure, temperature, and the index of refraction has been one of the important systematic uncertainties in reconstructing the depth of the shower maximum from the radio emission of air showers. Current air shower Monte Carlo simulation codes like CORSIKA and the radio plug-in CoREAS use various averaged parameterized atmospheres. How… ▽ More The limited knowledge of atmospheric parameters like humidity, pressure, temperature, and the index of refraction has been one of the important systematic uncertainties in reconstructing the depth of the shower maximum from the radio emission of air showers. Current air shower Monte Carlo simulation codes like CORSIKA and the radio plug-in CoREAS use various averaged parameterized atmospheres. However, time-dependent and location-specific atmospheric models are needed for the cosmic ray analysis method used for LOFAR data. There, dedicated simulation sets are used for each detected cosmic ray, to take into account the actual atmospheric conditions at the time of the measurement. Using the Global Data Assimilation System (GDAS), a global atmospheric model, we have implemented time-dependent, realistic atmospheric profiles in CORSIKA and CoREAS. We have produced realistic event-specific atmospheres for all air showers measured with LOFAR, an event set spanning several years and many different weather conditions. A complete re-analysis of our data set shows that for the majority of data, our previous correction factor performed rather well; we found only a small systematic shift of 2 g/cm$^2$ in the reconstructed $X_{\rm max}$. However, under extreme weather conditions, for example, very low air pressure, the shift can be up to 15 g/cm$^2$. We provide a correction formula to determine the shift in $X_{\rm max}$ resulting from a comparison of simulations done using the US-Std atmosphere and the GDAS-based atmosphere. △ Less

Submitted 2 June, 2020; originally announced June 2020.

Comments: Accepted for publication in Astroparticle Physics. arXiv admin note: text overlap with arXiv:1911.02859

On the cosmic-ray energy scale of the LOFAR radio telescope

Authors: K. Mulrey, S. Buitink, A. Corstanje, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, G. K. Krampah, P. Mitra, A. Nelles, H. Pandya, J. P. Rachen, O. Scholten, S. ter Veen, S. Thoudam, T. N. G. Trinh, T. Winchen

Abstract: Cosmic rays are routinely measured at LOFAR, both with a dense array of antennas and with the LOFAR Radboud air shower Array (LORA) which is an array of plastic scintillators. In this paper, we present two results relating to the cosmic-ray energy scale of LOFAR. First, we present the reconstruction of cosmic-ray energy using radio and particle techniques along with a discussion of the event-by-ev… ▽ More Cosmic rays are routinely measured at LOFAR, both with a dense array of antennas and with the LOFAR Radboud air shower Array (LORA) which is an array of plastic scintillators. In this paper, we present two results relating to the cosmic-ray energy scale of LOFAR. First, we present the reconstruction of cosmic-ray energy using radio and particle techniques along with a discussion of the event-by-event and absolute scale uncertainties. The resulting energies reconstructed with each method are shown to be in good agreement, and because the radio-based reconstructed energy has smaller uncertainty on an event-to-event basis, LOFAR analyses will use that technique in the future. Second, we present the radiation energy of air showers measured at LOFAR and demonstrate how radiation energy can be used to compare the energy scales of different experiments. The radiation energy scales quadratically with the electromagnetic energy in an air shower, which can in turn be related to the energy of the primary particle. Once the local magnetic field is accounted for, the radiation energy allows for a direct comparison between the LORA particle-based energy scale and that of the Pierre Auger Observatory. They are shown to agree to within (6$\pm$20)% for a radiation energy of 1 MeV, where the uncertainty on the comparison is dominated by the antenna calibrations of each experiment. This study motivates the development of a portable radio array which will be used to cross-calibrate the energy scales of different experiments using radiation energy and the same antennas, thereby significantly reducing the uncertainty on the comparison. △ Less

Submitted 29 September, 2020; v1 submitted 27 May, 2020; originally announced May 2020.

Search for magnetically-induced signatures in the arrival directions of ultra-high-energy cosmic rays measured at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker , et al. (350 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 27 July, 2020; v1 submitted 22 April, 2020; originally announced April 2020.

Comments: Published version; correction of a minor typo

Report number: FERMILAB-PUB-20-177-E-TD

Journal ref: Published in JCAP06(2020)017

Cosmic-ray anisotropies in right ascension measured by the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker , et al. (351 additional authors not shown)

Abstract: 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 e… ▽ More 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$σ$ 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 $α_{\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. △ Less

Submitted 14 February, 2020; originally announced February 2020.

Comments: Accepted by ApJ

Report number: FERMILAB-PUB-20-076

Journal ref: Published in Astrophys. J 891 (2020)142

Galactic Magnetic Field Bias on Inferences from UHECR Data

Authors: Björn Eichmann, Tobias Winchen

Abstract: A consequence of Liouville's theorem indicates that the recently observed large scale anisotropy in the arrival direction of Ultra-High-Energy Cosmic Rays (UHECRs) cannot be produced by the Galactic magnetic field, thus this anisotropy already needs to be present outside our Galaxy. But in this case, the observed energy spectrum and composition of UHECRs differs from the one outside of the Milky W… ▽ More A consequence of Liouville's theorem indicates that the recently observed large scale anisotropy in the arrival direction of Ultra-High-Energy Cosmic Rays (UHECRs) cannot be produced by the Galactic magnetic field, thus this anisotropy already needs to be present outside our Galaxy. But in this case, the observed energy spectrum and composition of UHECRs differs from the one outside of the Milky Way, due to the suppression or the amplification of the UHECR flux from certain directions by the Galactic magnetic field. In this work, we investigate this effect for the case of a dipole and a quadrupole anisotropy, respectively, for the widely-used JF12 magnetic field model. We investigate boundaries on the maximal amplitude of the observed anisotropy and the maximal charge number of UHECRs. Furthermore, the flux modification is discussed in the light of the Auger data on the recent dipole and also the chemical composition. We find that this modification effect yields a modification of the observed flux of up to $\sim 10\%$ for the investigated magnetic field model and the observed dipole, in particular for a heavy chemical composition of UHECRs as suggested by the 'EPOS-LHC' model. △ Less

Submitted 27 January, 2020; v1 submitted 6 January, 2020; originally announced January 2020.

Contributions of the LOFAR Cosmic Ray Key Science Project to the 36th International Cosmic Ray Conference (ICRC 2019)

Authors: A. Bonardi, S. Buitink, A. Corstanje, H. Falcke, B. M. Hare, J. R. Horandel, T. Huege, G. Krampah, P. Mitra, K. Mulrey, A. Nelles, H. Pandya, J. P. Rachen, L. Rossetto, O. Scholten, S. ter Veen, T. N. G. Trinh, T. Winchen

Abstract: This is a collection of papers that have been contributed by the LOFAR Cosmic Ray Key Science Project (CRKSP) to the 36th International Cosmic Ray Conference held in Madison, Wisconsin, on July 24th to August 1st, 2019 (ICRC 2019). All papers contained here have been individually published in PoS(ICRC2019) with paper numbers 205, 362, 352, 416, and 363, in the order they appear in this collection.… ▽ More This is a collection of papers that have been contributed by the LOFAR Cosmic Ray Key Science Project (CRKSP) to the 36th International Cosmic Ray Conference held in Madison, Wisconsin, on July 24th to August 1st, 2019 (ICRC 2019). All papers contained here have been individually published in PoS(ICRC2019) with paper numbers 205, 362, 352, 416, and 363, in the order they appear in this collection. Minor modifications to the PoS versions have been applied where appropriate. △ Less

Submitted 7 November, 2019; originally announced November 2019.

Comments: 42 pages, 17 figures, Collection of papers presented at the ICRC2019

The Pierre Auger Observatory: Contributions to the 36th International Cosmic Ray Conference (ICRC 2019)

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (361 additional authors not shown)

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

Submitted 19 September, 2019; originally announced September 2019.

Comments: compendium of 29 proceedings for the ICRC 2019

Probing the origin of ultra-high-energy cosmic rays with neutrinos in the EeV energy range using the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (367 additional authors not shown)

Abstract: 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 $τ$ neutrinos with nearly tangential trajectories relative to the earth. No neutrino candidates were found in… ▽ More 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 $τ$ 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_ν^{-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_ν/{\rm d}E_ν< 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. △ Less

Submitted 16 October, 2019; v1 submitted 18 June, 2019; originally announced June 2019.

Comments: Published version, 27 pages, 7 figures, 2 tables

Report number: FERMILAB-PUB-19-280

Journal ref: JCAP 10 (2019) 022

Limits on point-like sources of ultra-high-energy neutrinos with the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (367 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 8 November, 2019; v1 submitted 18 June, 2019; originally announced June 2019.

Comments: 25 pages, 9 figures

Report number: FERMILAB-PUB-19-279

Journal ref: JCAP 11 (2019) 004

NuRadioMC: Simulating the radio emission of neutrinos from interaction to detector

Authors: Christian Glaser, Daniel García-Fernández, Anna Nelles, Jaime Alvarez-Muñiz, Steven W. Barwick, Dave Z. Besson, Brian A. Clark, Amy Connolly, Cosmin Deaconu, Krijn de Vries, Jordan C. Hanson, Ben Hokanson-Fasig, R. Lahmann, Uzair Latif, Stuart A. Kleinfelder, Christopher Persichilli, Yue Pan, Carl Pfender, Ilse Plaisier, Dave Seckel, Jorge Torres, Simona Toscano, Nick van Eijndhoven, Abigail Vieregg, Christoph Welling , et al. (2 additional authors not shown)

Abstract: NuRadioMC is a Monte Carlo framework designed to simulate ultra-high energy neutrino detectors that rely on the radio detection method. This method exploits the radio emission generated in the electromagnetic component of a particle shower following a neutrino interaction. NuRadioMC simulates everything from the neutrino interaction in a medium, the subsequent Askaryan radio emission, the propagat… ▽ More NuRadioMC is a Monte Carlo framework designed to simulate ultra-high energy neutrino detectors that rely on the radio detection method. This method exploits the radio emission generated in the electromagnetic component of a particle shower following a neutrino interaction. NuRadioMC simulates everything from the neutrino interaction in a medium, the subsequent Askaryan radio emission, the propagation of the radio signal to the detector and finally the detector response. NuRadioMC is designed as a modern, modular Python-based framework, combining flexibility in detector design with user-friendliness. It includes a state-of-the-art event generator, an improved modelling of the radio emission, a revisited approach to signal propagation and increased flexibility and precision in the detector simulation. This paper focuses on the implemented physics processes and their implications for detector design. A variety of models and parameterizations for the radio emission of neutrino-induced showers are compared and reviewed. Comprehensive examples are used to discuss the capabilities of the code and different aspects of instrumental design decisions. △ Less

Submitted 3 February, 2020; v1 submitted 4 June, 2019; originally announced June 2019.

Comments: replaced with published version

Journal ref: Eur. Phys. J. C 80, 77 (2020)

Multi-Messenger Physics with the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (368 additional authors not shown)

Abstract: 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 co… ▽ More 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. △ Less

Submitted 26 April, 2019; originally announced April 2019.

Comments: 28 pages, 7 figures, Review published in Frontiers in Astronomy and Space Sciences

Report number: FERMILAB-PUB-19-135

Journal ref: April 2019 | Volume 6 | Article 24

Cosmic Ray Physics with the LOFAR Radio Telescope

Authors: T Winchen, A Bonardi, S Buitink, A Corstanje, H Falcke, B M Hare, J R Hörandel, P Mitra, K Mulrey, A Nelles, J P Rachen, L Rossetto, P Schellart, O Scholten, S ter Veen, S Thoudam, T N G Trinh

Abstract: The LOFAR radio telescope is able to measure the radio emission from cosmic ray induced air showers with hundreds of individual antennas. This allows for precision testing of the emission mechanisms for the radio signal as well as determination of the depth of shower maximum $X_{\max}$, the shower observable most sensitive to the mass of the primary cosmic ray, to better than 20 g/cm$^2$. With a d… ▽ More The LOFAR radio telescope is able to measure the radio emission from cosmic ray induced air showers with hundreds of individual antennas. This allows for precision testing of the emission mechanisms for the radio signal as well as determination of the depth of shower maximum $X_{\max}$, the shower observable most sensitive to the mass of the primary cosmic ray, to better than 20 g/cm$^2$. With a densely instrumented circular area of roughly 320 m$^2$, LOFAR is targeting for cosmic ray astrophysics in the energy range $10^{16}$ - $10^{18}$ eV. In this contribution we give an overview of the status, recent results, and future plans of cosmic ray detection with the LOFAR radio telescope. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: Proceedings of the 26th Extended European Cosmic Ray Symposium (ECRS), Barnaul/Belokurikha, 2018

CRPropa - A Toolbox for Cosmic Ray Simulations

Authors: R. Alves Batista, J. Becker Tjus, A. Dundovic, M. Erdmann, C. Heiter, K. -H. Kampert, D. Kuempel, L. Merten, G. Müller, G. Sigl, A. v. Vliet, D. Walz, T. Winchen, M. Wirtz

Abstract: The astrophysical interpretation of recent experimental observations of cosmic rays relies increasingly on Monte Carlo simulations of cosmic ray propagation and acceleration. Depending on the energy range of interest, several different propagation effects inside the Milky Way as well as in extragalactic space have to be taken into account when interpreting the data. With the CRPropa framework we a… ▽ More The astrophysical interpretation of recent experimental observations of cosmic rays relies increasingly on Monte Carlo simulations of cosmic ray propagation and acceleration. Depending on the energy range of interest, several different propagation effects inside the Milky Way as well as in extragalactic space have to be taken into account when interpreting the data. With the CRPropa framework we aim to provide a toolbox for according simulations. In recent versions of CRPropa, the ballistic single particle propagation mode aiming primarily at extragalactic cosmic rays has been complemented by a solver for the differential transport equation to address propagation of galactic cosmic rays. Additionally, modules have been developed to address cosmic ray acceleration and many improvements have been added for simulations of electromagnetic secondaries. In this contribution we will give an overview of the CRPropa simulation framework with a focus on the latest improvements and highlight selected features by example applications. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: Proceedings of the 26th Extended European Cosmic Ray Symposium (ECRS), Barnaul/Belokurikha, 2018

Status of the Lunar Detection Mode for Cosmic Particles of LOFAR

Authors: T. Winchen, A. Bonardi, S. Buitink, A. Corstanje, H. Falcke, B. M. Hare, J. R. Hörandel, P. Mitra, K. Mulrey, A. Nelles, J. P. Rachen, L. Rossetto, P. Schellart, O. Scholten, S. ter Veen, S. Thoudam, T. N. G. Trinh

Abstract: Cosmic particles hitting Earth's moon produce radio emission via the Askaryan effect. If the resulting radio ns-pulse can be detected by radio telescopes, this technique potentially increases the available collective area for ZeV scale particles by several orders of magnitude compared to current experiments. The LOw Frequency ARray (LOFAR) is the largest radio telescope operating in the optimum fr… ▽ More Cosmic particles hitting Earth's moon produce radio emission via the Askaryan effect. If the resulting radio ns-pulse can be detected by radio telescopes, this technique potentially increases the available collective area for ZeV scale particles by several orders of magnitude compared to current experiments. The LOw Frequency ARray (LOFAR) is the largest radio telescope operating in the optimum frequency regime for this technique. In this contribution, we report on the status of the implementation of the lunar detection mode at LOFAR. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: Proceedings of the 26th Extended European Cosmic Ray Symposium (ECRS), Barnaul/Belokurikha, 2018

Calibration of the LOFAR low-band antennas using the Galaxy and a model of the signal chain

Authors: K. Mulrey, A. Bonardi, S. Buitink, A. Corstanje, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, P. Mitra, A. Nelles, J. P. Rachen, L. Rossetto, P. Schellart, O. Scholten, S. ter Veen, S. Thoudam, T. N. G. Trinh, T. Winchen

Abstract: The LOw-Frequency ARray (LOFAR) is used to make precise measurements of radio emission from extensive air showers, yielding information about the primary cosmic ray. Interpreting the measured data requires an absolute and frequency-dependent calibration of the LOFAR system response. This is particularly important for spectral analyses, because the shape of the detected signal holds information abo… ▽ More The LOw-Frequency ARray (LOFAR) is used to make precise measurements of radio emission from extensive air showers, yielding information about the primary cosmic ray. Interpreting the measured data requires an absolute and frequency-dependent calibration of the LOFAR system response. This is particularly important for spectral analyses, because the shape of the detected signal holds information about the shower development. We revisit the calibration of the LOFAR antennas in the range of 30 - 80 MHz. Using the Galactic emission and a detailed model of the LOFAR signal chain, we find an improved calibration that provides an absolute energy scale and allows for the study of frequency-dependent features in measured signals. With the new calibration, systematic uncertainties of 13% are reached, and comparisons of the spectral shape of calibrated data with simulations show promising agreement. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Comments: 23 pages, 10 figures

Data-driven estimation of the invisible energy of cosmic ray showers with the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (367 additional authors not shown)

Abstract: 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 producti… ▽ More 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. △ Less

Submitted 6 November, 2019; v1 submitted 23 January, 2019; originally announced January 2019.

Comments: Published version, 18 pages, 10 figures, 4 tables

Report number: FERMILAB-PUB-19-040

Journal ref: Phys. Rev. D 100, 082003 (2019)

Measurement of the average shape of longitudinal profiles of cosmic-ray air showers at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvar ez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (363 additional authors not shown)

Abstract: 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 Gais… ▽ More 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. △ Less

Submitted 16 May, 2019; v1 submitted 12 November, 2018; originally announced November 2018.

Comments: Published version, 20 pages, 6 figures, 2 tables

Report number: FERMILAB-PUB-18-633

Journal ref: JCAP03(2019)018

RadioPropa - A Modular Raytracer for In-Matter Radio Propagation

Authors: Tobias Winchen

Abstract: Experiments for radio detection of UHE particles such as e.g. ARA/ARIANNA or NuMoon require detailed understanding of the propagation of radio waves in the surrounding matter. The index of refraction in e.g. polar ice or lunar rock may have a complex spatial structure that makes detailed simulations of the radio propagation necessary to design the respective experiments and analyze their data. Her… ▽ More Experiments for radio detection of UHE particles such as e.g. ARA/ARIANNA or NuMoon require detailed understanding of the propagation of radio waves in the surrounding matter. The index of refraction in e.g. polar ice or lunar rock may have a complex spatial structure that makes detailed simulations of the radio propagation necessary to design the respective experiments and analyze their data. Here, we present RadioPropa as a new modular ray tracing code that solves the eikonal equation with a Runge-Kutta method in arbitrary refractivity fields. RadioPropa is based on the cosmic ray propagation code CRPropa, which has been forked to allow efficient incorporation of the required data structures for ray tracing while retaining its modular design. This allows for the setup of versatile simulation geometries as well as the easy inclusion of additional physical effects such as e.g. partial reflection on boundary layers in the simulations. We discuss the principal design of the code as well as its performance in example applications. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: Proceedings of the ARENA Conference 2018

Contributions of the LOFAR Cosmic Ray Key Science Project to the 35th International Cosmic Ray Conference (ICRC 2017)

Authors: A. Bonardi, S. Buitink, A. Corstanje, H. Falcke, B. M. Hare, J. R. Hörandel, P. Mitra, K. Mulrey, A. Nelles, J. P. Rachen, L. Rossetto, P. Schellart, O. Scholten, S. ter Veen, S. Thoudam, T. N. G. Trinh, T. Winchen

Abstract: Contributions of the LOFAR Cosmic Ray Key Science Project to the 35th International Cosmic Ray Conference (ICRC 2017) Contributions of the LOFAR Cosmic Ray Key Science Project to the 35th International Cosmic Ray Conference (ICRC 2017) △ Less

Submitted 21 November, 2017; originally announced November 2017.

Comments: compendium of 8 proceedings for the ICRC 2017

The Pierre Auger Observatory: Contributions to the 35th International Cosmic Ray Conference (ICRC 2017)

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. F. M. Albuquerque, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, N. Arsene, H. Asorey, P. Assis, J. Aublin, G. Avila, A. M. Badescu, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker, J. A. Bellido , et al. (373 additional authors not shown)

Abstract: Contributions of the Pierre Auger Collaboration to the 35th International Cosmic Ray Conference (ICRC 2017), 12-20 July 2017, Bexco, Busan, Korea. Contributions of the Pierre Auger Collaboration to the 35th International Cosmic Ray Conference (ICRC 2017), 12-20 July 2017, Bexco, Busan, Korea. △ Less

Submitted 2 October, 2017; v1 submitted 22 August, 2017; originally announced August 2017.

Comments: compendium of 26 proceedings for the ICRC 2017

Overview of lunar detection of ultra-high energy particles and new plans for the SKA

Authors: Clancy W. James, Jaime Alvarez-Muñiz, Justin D. Bray, Stijn Buitink, Rustam D. Dagkesamanskii, Ronald D. Ekers, Heino Falcke, Ken Gayley, Tim Huege, Maaijke Mevius, Rob Mutel, Olaf Scholten, Ralph Spencer, Sander ter Veen, Tobias Winchen

Abstract: The lunar technique is a method for maximising the collection area for ultra-high-energy (UHE) cosmic ray and neutrino searches. The method uses either ground-based radio telescopes or lunar orbiters to search for Askaryan emission from particles cascading near the lunar surface. While experiments using the technique have made important advances in the detection of nanosecond-scale pulses, only at… ▽ More The lunar technique is a method for maximising the collection area for ultra-high-energy (UHE) cosmic ray and neutrino searches. The method uses either ground-based radio telescopes or lunar orbiters to search for Askaryan emission from particles cascading near the lunar surface. While experiments using the technique have made important advances in the detection of nanosecond-scale pulses, only at the very highest energies has the lunar technique achieved competitive limits. This is expected to change with the advent of the Square Kilometre Array (SKA), the low-frequency component of which (SKA-low) is predicted to be able to detect an unprecedented number of UHE cosmic rays. In this contribution, the status of lunar particle detection is reviewed, with particular attention paid to outstanding theoretical questions, and the technical challenges of using a giant radio array to search for nanosecond pulses. The activities of SKA's High Energy Cosmic Particles Focus Group are described, as is a roadmap by which this group plans to incorporate this detection mode into SKA-low observations. Estimates for the sensitivity of SKA-low phases 1 and 2 to UHE particles are given, along with the achievable science goals with each stage. Prospects for near-future observations with other instruments are also described. △ Less

Submitted 18 April, 2017; originally announced April 2017.

Comments: 6 pages, 3 figures, presented at the 7th International Conference on Acoustic and Radio EeV Neutrino Detection Activities (ARENA 2016), Groningen, The Netherlands

Journal ref: EPJ Web Conf., Volume 135, 04001 (2017)

Muon Counting using Silicon Photomultipliers in the AMIGA detector of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, E. J. Ahn, I. Al Samarai, I. F. M. Albuquerque, I. Allekotte, P. Allison, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, M. Ambrosio, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, F. Arqueros, N. Arsene, H. Asorey, P. Assis, J. Aublin, G. Avila, A. M. Badescu , et al. (400 additional authors not shown)

Abstract: AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory designed to extend its energy range of detection and to directly measure the muon content of the cosmic ray primary particle showers. The array will be formed by an infill of surface water-Cherenkov detectors associated with buried scintillation counters employed for muon counting. Each counter is com… ▽ More AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory designed to extend its energy range of detection and to directly measure the muon content of the cosmic ray primary particle showers. The array will be formed by an infill of surface water-Cherenkov detectors associated with buried scintillation counters employed for muon counting. Each counter is composed of three scintillation modules, with a 10 m$^2$ detection area per module. In this paper, a new generation of detectors, replacing the current multi-pixel photomultiplier tube (PMT) with silicon photo sensors (aka. SiPMs), is proposed. The selection of the new device and its front-end electronics is explained. A method to calibrate the counting system that ensures the performance of the detector is detailed. This method has the advantage of being able to be carried out in a remote place such as the one where the detectors are deployed. High efficiency results, i.e. 98 % efficiency for the highest tested overvoltage, combined with a low probability of accidental counting ($\sim$2 %), show a promising performance for this new system. △ Less

Submitted 4 October, 2017; v1 submitted 17 March, 2017; originally announced March 2017.

Comments: 24 pages, 12 figures

Report number: FERMILAB-PUB-16-656-AD-CD-TD

Journal ref: JINST 12 (2017) P03002

Thunderstorm electric fields probed by extensive air showers through their polarized radio emission

Authors: T. N. G. Trinh, O. Scholten, A. Bonardi, S. Buitink, A. Corstanje, U. Ebert, J. E. Enriquez, H. Falcke, J. R. Horandel, B. M. Hare, P. Mitra, K. Mulrey, A. Nelles, J. P. Rachen, L. Rossetto, C. Rutjes, P. Schellart, S. Thoudam, S. ter Veen, T. Winchen

Abstract: We observe a large fraction of circular polarization in radio emission from extensive air showers recorded during thunderstorms, much higher than in the emission from air showers measured during fair-weather circumstances. We show that the circular polarization of the air showers measured during thunderstorms can be explained by the change in the direction of the transverse current as a function o… ▽ More We observe a large fraction of circular polarization in radio emission from extensive air showers recorded during thunderstorms, much higher than in the emission from air showers measured during fair-weather circumstances. We show that the circular polarization of the air showers measured during thunderstorms can be explained by the change in the direction of the transverse current as a function of altitude induced by atmospheric electric fields. Thus by using the full set of Stokes parameters for these events, we obtain a good characterization of the electric fields in thunderclouds. We also measure a large horizontal component of the electric fields in the two events that we have analysed. △ Less

Submitted 20 March, 2017; v1 submitted 14 March, 2017; originally announced March 2017.

Calibration of the Logarithmic-Periodic Dipole Antenna (LPDA) Radio Stations at the Pierre Auger Observatory using an Octocopter

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. Al Samarai, I. F. M. Albuquerque, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, F. Arqueros, N. Arsene, H. Asorey, P. Assis, J. Aublin, G. Avila, A. M. Badescu, A. Balaceanu, F. Barbato, R. J. Barreira Luz , et al. (380 additional authors not shown)

Abstract: 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 broad… ▽ More 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°. △ Less

Submitted 13 June, 2018; v1 submitted 5 February, 2017; originally announced February 2017.

Comments: Published version. Updated online abstract only. Manuscript is unchanged with respect to v2. 39 pages, 15 figures, 2 tables

Report number: FERMILAB-PUB-17-046

Journal ref: JINST 12 (2017) T10005

The effect of the atmospheric refractive index on the radio signal of extensive air showers

Authors: A. Corstanje, A. Bonardi, S. Buitink, H. Falcke, J. R. Hörandel, P. Mitra, K. Mulrey, A. Nelles, J. P. Rachen, L. Rossetto, P. Schellart, O. Scholten, S. ter Veen, S. Thoudam, G. Trinh, T. Winchen

Abstract: For the interpretation of measurements of radio emission from extensive air showers, an important systematic uncertainty arises from natural variations of the atmospheric refractive index $n$. At a given altitude, the refractivity $N=10^6\, (n-1)$ can have relative variations on the order of $10 \%$ depending on temperature, humidity, and air pressure. Typical corrections to be applied to $N$ are… ▽ More For the interpretation of measurements of radio emission from extensive air showers, an important systematic uncertainty arises from natural variations of the atmospheric refractive index $n$. At a given altitude, the refractivity $N=10^6\, (n-1)$ can have relative variations on the order of $10 \%$ depending on temperature, humidity, and air pressure. Typical corrections to be applied to $N$ are about $4\%$. Using CoREAS simulations of radio emission from air showers, we have evaluated the effect of varying $N$ on measurements of the depth of shower maximum $X_{\rm max}$. For an observation band of 30 to 80 MHz, a difference of $4 \%$ in refractivity gives rise to a systematic error in the inferred $X_{\rm max}$ between 3.5 and 11 $\mathrm{g/cm^2}$, for proton showers with zenith angles ranging from 15 to 50 degrees. At higher frequencies, from 120 to 250 MHz, the offset ranges from 10 to 22 $\mathrm{g/cm^2}$. These offsets were found to be proportional to the geometric distance to $X_{\rm max}$. We have compared the results to a simple model based on the Cherenkov angle. For the 120 to 250 MHz band, the model is in qualitative agreement with the simulations. In typical circumstances, we find a slight decrease in $X_{\rm max}$ compared to the default refractivity treatment in CoREAS. While this is within commonly treated systematic uncertainties, accounting for it explicitly improves the accuracy of $X_{\rm max}$ measurements. △ Less

Submitted 25 January, 2017; originally announced January 2017.

Comments: 13 pages, 5 figures. Accepted for publication in Astroparticle Physics

Combined fit of spectrum and composition data as measured by the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. Al Samarai, I. F. M. Albuquerque, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, F. Arqueros, N. Arsene, H. Asorey, P. Assis, J. Aublin, G. Avila, A. M. Badescu, A. Balaceanu, R. J. Barreira Luz, J. J. Beatty , et al. (375 additional authors not shown)

Abstract: 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 dist… ▽ More 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. △ Less

Submitted 26 February, 2018; v1 submitted 21 December, 2016; originally announced December 2016.

Comments: Corrected typos in values of SD vertical spectrum energy bias (B.5); Added expression for SD inclined spectrum energy resolution (B.6), erroneously claimed to have been published before

Report number: FERMILAB-PUB-16-618

Journal ref: JCAP04(2017)038

Realtime processing of LOFAR data for the detection of nano-second pulses from the Moon

Authors: T. Winchen, A. Bonardi, S. Buitink, A. Corstanje, J. E. Enriquez, H. Falcke, J. R. Hörandel, P. Mitra, K. Mulrey, A. Nelles, J. P. Rachen, L. Rossetto, P. Schellart, O. Scholten, S. Thoudam, T. N. G. Trinh, S. ter Veen

Abstract: The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. Even lower fluxes of neutrinos with energies above $10^{22}$ eV are predicted in certain Grand-Unifying-Theories (GUTs) and e.g.\ models for super-heavy dark matter (SHDM). The significant increase in detector volume required to… ▽ More The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. Even lower fluxes of neutrinos with energies above $10^{22}$ eV are predicted in certain Grand-Unifying-Theories (GUTs) and e.g.\ models for super-heavy dark matter (SHDM). The significant increase in detector volume required to detect these particles can be achieved by searching for the nano-second radio pulses that are emitted when a particle interacts in Earth's moon with current and future radio telescopes. In this contribution we present the design of an online analysis and trigger pipeline for the detection of nano-second pulses with the LOFAR radio telescope. The most important steps of the processing pipeline are digital focusing of the antennas towards the Moon, correction of the signal for ionospheric dispersion, and synthesis of the time-domain signal from the polyphased-filtered signal in frequency domain. The implementation of the pipeline on a GPU/CPU cluster will be discussed together with the computing performance of the prototype. △ Less

Submitted 20 December, 2016; originally announced December 2016.

Comments: Proceedings of the 22nd International Conference on Computing in High Energy and Nuclear Physics (CHEP2016), USA

Journal ref: Journal of Physics: Conference Series 898(3), 032004, 2017

A targeted search for point sources of EeV photons with the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. Al Samarai, I. F. M. Albuquerque, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, F. Arqueros, N. Arsene, H. Asorey, P. Assis, J. Aublin, G. Avila, A. M. Badescu, A. Balaceanu, R. J. Barreira Luz, J. J. Beatty , et al. (375 additional authors not shown)

Abstract: Simultaneous measurements of air showers with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for EeV photon point sources. Several Galactic and extragalactic candidate objects are grouped in classes to reduce the statistical penalty of many trials from that of a blind search and are analyzed for a significant excess above the background expectation.… ▽ More Simultaneous measurements of air showers with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for EeV photon point sources. Several Galactic and extragalactic candidate objects are grouped in classes to reduce the statistical penalty of many trials from that of a blind search and are analyzed for a significant excess above the background expectation. The presented search does not find any evidence for photon emission at candidate sources, and combined $p$-values for every class are reported. Particle and energy flux upper limits are given for selected candidate sources. These limits significantly constrain predictions of EeV proton emission models from non-transient Galactic and nearby extragalactic sources, as illustrated for the particular case of the Galactic center region. △ Less

Submitted 21 March, 2017; v1 submitted 13 December, 2016; originally announced December 2016.

Comments: Published version

Report number: FERMILAB-PUB-16-655

Journal ref: The Astrophysical Journal Letters, 837:L25 (7pp), 2017 March 10

Energy Spectrum of Fast Second Order Fermi Accelerators as Sources of Ultra-High-Energy Cosmic Rays

Authors: Tobias Winchen, Stijn Buitink

Abstract: Stochastic acceleration of cosmic rays in second order Fermi processes is usually considered too slow to reach ultra-high energies, except in specific cases. In this paper we present the energy spectrum obtained from second order Fermi acceleration in highly turbulent magnetic fields as e.g.\ found in the outskirts of AGN jets in situations where it can be sufficiently fast to accelerate particles… ▽ More Stochastic acceleration of cosmic rays in second order Fermi processes is usually considered too slow to reach ultra-high energies, except in specific cases. In this paper we present the energy spectrum obtained from second order Fermi acceleration in highly turbulent magnetic fields as e.g.\ found in the outskirts of AGN jets in situations where it can be sufficiently fast to accelerate particles to the highest observed energies. We parametrize the resulting non-power-law spectra and show that these can describe the cosmic ray energy spectrum and mass-composition data at the highest energies if propagation effects are taken into account. △ Less

Submitted 18 April, 2018; v1 submitted 12 December, 2016; originally announced December 2016.

Journal ref: Astroparticle Physics 102 (2018) 25 - 31

Search for photons with energies above 10$^{18}$ eV using the hybrid detector of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. Al Samarai, I. F. M. Albuquerque, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, F. Arqueros, N. Arsene, H. Asorey, P. Assis, J. Aublin, G. Avila, A. M. Badescu, A. Balaceanu, R. J. Barreira Luz, J. J. Beatty , et al. (375 additional authors not shown)

Abstract: A search for ultra-high energy photons with energies above 1 EeV is performed using nine years of data collected by the Pierre Auger Observatory in hybrid operation mode. An unprecedented separation power between photon and hadron primaries is achieved by combining measurements of the longitudinal air-shower development with the particle content at ground measured by the fluorescence and surface d… ▽ More A search for ultra-high energy photons with energies above 1 EeV is performed using nine years of data collected by the Pierre Auger Observatory in hybrid operation mode. An unprecedented separation power between photon and hadron primaries is achieved by combining measurements of the longitudinal air-shower development with the particle content at ground measured by the fluorescence and surface detectors, respectively. Only three photon candidates at energies 1 - 2 EeV are found, which is compatible with the expected hadron-induced background. Upper limits on the integral flux of ultra-high energy photons of 0.038, 0.010, 0.009, 0.008 and 0.007 km$^{-2}$ sr$^{-1}$ yr$^{-1}$ are derived at 95% C.L. for energy thresholds of 1, 2, 3, 5 and 10 EeV. These limits bound the fractions of photons in the all-particle integral flux below 0.14%, 0.17%, 0.42%, 0.86% and 2.9%. For the first time the photon fraction at EeV energies is constrained at the sub-percent level. The improved limits are below the flux of diffuse photons predicted by some astrophysical scenarios for cosmogenic photon production. The new results rule-out the early top-down models $-$ in which ultra-high energy cosmic rays are produced by, e.g., the decay of super-massive particles $-$ and challenge the most recent super-heavy dark matter models. △ Less

Submitted 28 September, 2020; v1 submitted 5 December, 2016; originally announced December 2016.

Comments: Corrected version after erratum published in JCAP

Report number: FERMILAB-PUB-16-657

Journal ref: JCAP04(2017)009; JCAP09(2020)E02

Multi-resolution anisotropy studies of ultrahigh-energy cosmic rays detected at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I. Al Samarai, I. F. M. Albuquerque, I. Allekotte, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, F. Arqueros, N. Arsene, H. Asorey, P. Assis, J. Aublin, G. Avila, A. M. Badescu, A. Balaceanu, R. J. Barreira Luz, C. Baus , et al. (378 additional authors not shown)

Abstract: We report a multi-resolution search for anisotropies in the arrival directions of cosmic rays detected at the Pierre Auger Observatory with local zenith angles up to $80^\circ$ and energies in excess of 4 EeV ($4 \times 10^{18}$ eV). This search is conducted by measuring the angular power spectrum and performing a needlet wavelet analysis in two independent energy ranges. Both analyses are complem… ▽ More We report a multi-resolution search for anisotropies in the arrival directions of cosmic rays detected at the Pierre Auger Observatory with local zenith angles up to $80^\circ$ and energies in excess of 4 EeV ($4 \times 10^{18}$ eV). This search is conducted by measuring the angular power spectrum and performing a needlet wavelet analysis in two independent energy ranges. Both analyses are complementary since the angular power spectrum achieves a better performance in identifying large-scale patterns while the needlet wavelet analysis, considering the parameters used in this work, presents a higher efficiency in detecting smaller-scale anisotropies, potentially providing directional information on any observed anisotropies. No deviation from isotropy is observed on any angular scale in the energy range between 4 and 8 EeV. Above 8 EeV, an indication for a dipole moment is captured; while no other deviation from isotropy is observed for moments beyond the dipole one. The corresponding $p$-values obtained after accounting for searches blindly performed at several angular scales, are $1.3 \times 10^{-5}$ in the case of the angular power spectrum, and $2.5 \times 10^{-3}$ in the case of the needlet analysis. While these results are consistent with previous reports making use of the same data set, they provide extensions of the previous works through the thorough scans of the angular scales. △ Less

Submitted 20 June, 2017; v1 submitted 21 November, 2016; originally announced November 2016.

Comments: Published version. Added journal reference and DOI. Added Report Number

Report number: FERMILAB-PUB-16-556-AD-AE-CD-TD

Journal ref: JCAP 06 (2017) 026

Measurement of the circular polarization in radio emission from extensive air showers confirms emission mechanisms

Authors: O. Scholten, T. N. G. Trinh, A. Bonardi, S. Buitink, P. Correa, A. Corstanje, Q. Dorosti Hasankiadeh, H. Falcke, J. R. Hörandel, P. Mitra, K. Mulrey, A. Nelles, J. P. Rachen, L. Rossetto, P. Schellart, S. Thoudam, S. ter Veen, K. D. de Vries, T. Winchen

Abstract: We report here on a novel analysis of the complete set of four Stokes parameters that uniquely determine the linear and/or circular polarization of the radio signal for an extensive air shower. The observed dependency of the circular polarization on azimuth angle and distance to the shower axis is a clear signature of the interfering contributions from two different radiation mechanisms, a main co… ▽ More We report here on a novel analysis of the complete set of four Stokes parameters that uniquely determine the linear and/or circular polarization of the radio signal for an extensive air shower. The observed dependency of the circular polarization on azimuth angle and distance to the shower axis is a clear signature of the interfering contributions from two different radiation mechanisms, a main contribution due to a geomagnetically-induced transverse current and a secondary component due to the build-up of excess charge at the shower front. The data, as measured at LOFAR, agree very well with a calculation from first principles. This opens the possibility to use circular polarization as an investigative tool in the analysis of air shower structure, such as for the determination of atmospheric electric fields. △ Less

Submitted 19 November, 2016; v1 submitted 2 November, 2016; originally announced November 2016.

Comments: Accepted for publication in Phys. Rev. D

Testing Hadronic Interactions at Ultrahigh Energies with Air Showers Measured by the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, E. J. Ahn, I. Al Samarai, I. F. M. Albuquerque, I. Allekotte, J. Allen, P. Allison, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, M. Ambrosio, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, F. Arqueros, N. Arsene, H. Asorey, P. Assis, J. Aublin, G. Avila , et al. (413 additional authors not shown)

Abstract: Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre… ▽ More Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons. △ Less

Submitted 31 October, 2016; v1 submitted 26 October, 2016; originally announced October 2016.

Comments: Published version. Added journal reference and DOI. Added Report Number

Report number: FERMILAB-PUB-16-504-AD-AE-CD-TD

Journal ref: Phys. Rev. Lett. 117, 192001 (2016)

Evidence for a mixed mass composition at the `ankle' in the cosmic-ray spectrum

Authors: The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, E. J. Ahn, I. Al Samarai, I. F. M. Albuquerque, I. Allekotte, P. Allison, A. Almela, J. Alvarez Castillo, J. Alvarez-Muñiz, M. Ambrosio, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, F. Arqueros, N. Arsene, H. Asorey, P. Assis, J. Aublin, G. Avila, A. M. Badescu , et al. (401 additional authors not shown)

Abstract: We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the… ▽ More We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the electromagnetic and muonic components. It allows an accurate determination of the spread of primary masses in the cosmic-ray flux. Up till now, constraints on the spread of primary masses have been dominated by systematic uncertainties. The present correlation measurement is not affected by systematics in the measurement of the depth of shower maximum or the signal in the water Cherenkov stations. The analysis relies on general characteristics of air showers and is thus robust also with respect to uncertainties in hadronic event generators. The observed correlation in the energy range around the `ankle' at $\lg(E/{\rm eV})=18.5-19.0$ differs significantly from expectations for pure primary cosmic-ray compositions. A light composition made up of proton and helium only is equally inconsistent with observations. The data are explained well by a mixed composition including nuclei with mass $A > 4$. Scenarios such as the proton dip model, with almost pure compositions, are thus disfavoured as the sole explanation of the ultrahigh-energy cosmic-ray flux at Earth. △ Less

Submitted 22 November, 2016; v1 submitted 27 September, 2016; originally announced September 2016.

Comments: Published version. Added journal reference and DOI. Added Report Number

Report number: FERMILAB-PUB-16-412-AD-AE-CD-TD

Journal ref: Phys. Lett. B762 (2016) 288-295