Showing 1–50 of 236 results for author: Becker, H

WisPerMed at "Discharge Me!": Advancing Text Generation in Healthcare with Large Language Models, Dynamic Expert Selection, and Priming Techniques on MIMIC-IV

Authors: Hendrik Damm, Tabea M. G. Pakull, Bahadır Eryılmaz, Helmut Becker, Ahmad Idrissi-Yaghir, Henning Schäfer, Sergej Schultenkämper, Christoph M. Friedrich

Abstract: This study aims to leverage state of the art language models to automate generating the "Brief Hospital Course" and "Discharge Instructions" sections of Discharge Summaries from the MIMIC-IV dataset, reducing clinicians' administrative workload. We investigate how automation can improve documentation accuracy, alleviate clinician burnout, and enhance operational efficacy in healthcare facilities.… ▽ More This study aims to leverage state of the art language models to automate generating the "Brief Hospital Course" and "Discharge Instructions" sections of Discharge Summaries from the MIMIC-IV dataset, reducing clinicians' administrative workload. We investigate how automation can improve documentation accuracy, alleviate clinician burnout, and enhance operational efficacy in healthcare facilities. This research was conducted within our participation in the Shared Task Discharge Me! at BioNLP @ ACL 2024. Various strategies were employed, including few-shot learning, instruction tuning, and Dynamic Expert Selection (DES), to develop models capable of generating the required text sections. Notably, utilizing an additional clinical domain-specific dataset demonstrated substantial potential to enhance clinical language processing. The DES method, which optimizes the selection of text outputs from multiple predictions, proved to be especially effective. It achieved the highest overall score of 0.332 in the competition, surpassing single-model outputs. This finding suggests that advanced deep learning methods in combination with DES can effectively automate parts of electronic health record documentation. These advancements could enhance patient care by freeing clinician time for patient interactions. The integration of text selection strategies represents a promising avenue for further research. △ Less

Submitted 18 May, 2024; originally announced May 2024.

Comments: 8 pages, 6 tables, 8 figures, submitted to: BioNLP 2024 and Shared Tasks @ ACL 2024

Comparative analysis of mathematical formulations for the two-dimensional guillotine cutting problem

Authors: Henrique Becker, Mateus Martin, Olinto Araujo, Luciana S. Buriol, Reinaldo Morabito

Abstract: About ten years ago, a paper proposed the first integer linear programming formulation for the constrained two-dimensional guillotine cutting problem (with unlimited cutting stages). Since, six other formulations followed, five of them in the last two years. This spike of interest gave no opportunity for a comprehensive comparison between the formulations. We review each formulation and compare th… ▽ More About ten years ago, a paper proposed the first integer linear programming formulation for the constrained two-dimensional guillotine cutting problem (with unlimited cutting stages). Since, six other formulations followed, five of them in the last two years. This spike of interest gave no opportunity for a comprehensive comparison between the formulations. We review each formulation and compare their empirical results over instance datasets of the literature. We adapt most formulations to allow for piece rotation. The possibility of adaptation was already predicted but not realized by the prior work. The results show the dominance of pseudo-polynomial formulations until the point instances become intractable by them, while more compact formulations keep achieving good primal solutions. Our study also reveals a small but consistent advantage of the Gurobi solver over the CPLEX solver in our context; that the choice of solver hardly benefits one formulation over another; and a mistake in the generation of the T instances, which should have the same optima with or without guillotine cuts. Our study also proposes hybridising the most recent formulation with a prior formulation for a restricted version of the problem. The hybridisations show a reduction of about 20% of the branch-and-bound time thanks to the symmetries broken by the hybridisation. △ Less

Submitted 9 August, 2023; originally announced August 2023.

Comments: 23 pages, 7 tables, 3 figures

MSC Class: 90-02 ACM Class: G.2.0

Exploring the Impact of Configurational Entropy on the Design and Development of CoNi-Based Superalloys for Sustainable Applications

Authors: Ahad Mohammadzadeha, Akbar Heidarzadeh, Hailey Becker, Jorge Valilla Robles, Alberto Meza, Manuel Avella, Miguel A. Monclus, Damien Tourret, Jose Manuel Torralba

Abstract: A comprehensive literature review on recently rediscovered Co- and/or CoNi-based superalloys, strengthened by the γ' phase, revealed a relationship between the configurational entropy of the system and the γ' solvus temperature. This study was conducted on a high Cr CoNi-based superalloy system with high configurational entropy to test our hypothesis based on the sustainable metallurgy framework.… ▽ More A comprehensive literature review on recently rediscovered Co- and/or CoNi-based superalloys, strengthened by the γ' phase, revealed a relationship between the configurational entropy of the system and the γ' solvus temperature. This study was conducted on a high Cr CoNi-based superalloy system with high configurational entropy to test our hypothesis based on the sustainable metallurgy framework. Thermodynamic calculations were performed to design the chemical compositions, followed by vacuum casting and heat treatments to produce the desired alloys. The microstructures were characterized using a scanning electron microscope, electron backscattered diffraction, transmission electron microscope, and differential thermal analysis. Microhardness and nanoindentation tests were employed to measure the mechanical properties. The results showed that both the configurational entropy and the type of alloying elements determine the final high-temperature performance of the alloys. We found that to enhance the higher γ' solvus temperature, the configurational entropy should be increased by adding γ' stabilizing elements. The microstructural and mechanical characteristics of the designed alloys before and after heat treatments are discussed in detail. The outcome of this study is beneficial for developing cobalt-based high-entropy superalloys with appropriate processing windows and freezing ranges for advanced sustainable manufacturing purposes, such as using powder bed fusion technologies. △ Less

Submitted 15 July, 2023; originally announced July 2023.

K-D Bonsai: ISA-Extensions to Compress K-D Trees for Autonomous Driving Tasks

Authors: Pedro H. E. Becker, José María Arnau, Antonio González

Abstract: Autonomous Driving (AD) systems extensively manipulate 3D point clouds for object detection and vehicle localization. Thereby, efficient processing of 3D point clouds is crucial in these systems. In this work we propose K-D Bonsai, a technique to cut down memory usage during radius search, a critical building block of point cloud processing. K-D Bonsai exploits value similarity in the data structu… ▽ More Autonomous Driving (AD) systems extensively manipulate 3D point clouds for object detection and vehicle localization. Thereby, efficient processing of 3D point clouds is crucial in these systems. In this work we propose K-D Bonsai, a technique to cut down memory usage during radius search, a critical building block of point cloud processing. K-D Bonsai exploits value similarity in the data structure that holds the point cloud (a k-d tree) to compress the data in memory. K-D Bonsai further compresses the data using a reduced floating-point representation, exploiting the physically limited range of point cloud values. For easy integration into nowadays systems, we implement K-D Bonsai through Bonsai-extensions, a small set of new CPU instructions to compress, decompress, and operate on points. To maintain baseline safety levels, we carefully craft the Bonsai-extensions to detect precision loss due to compression, allowing re-computation in full precision to take place if necessary. Therefore, K-D Bonsai reduces data movement, improving performance and energy efficiency, while guaranteeing baseline accuracy and programmability. We evaluate K-D Bonsai over the euclidean cluster task of Autoware.ai, a state-of-the-art software stack for AD. We achieve an average of 9.26% improvement in end-to-end latency, 12.19% in tail latency, and a reduction of 10.84% in energy consumption. Differently from expensive accelerators proposed in related work, K-D Bonsai improves radius search with minimal area increase (0.36%). △ Less

Submitted 30 August, 2023; v1 submitted 1 February, 2023; originally announced February 2023.

MSC Class: Article No. 18; 2018 Related DOI: https://doi.org/10.1145/3243176.3243184 Focus to learn more

Journal ref: ISCA'23 Proceedings of the 50th Annual International Symposium on Computer Architecture, Article No. 20, 2023

HR-EBSD analysis of in situ stable crack growth at the micron scale

Authors: Abdalrhaman Koko, Thorsten H. Becker, Elsiddig Elmukashfi, Nicola M. Pugno, Angus J. Wilkinson, T. James Marrow

Abstract: Understanding the local fracture resistance of microstructural features. such as brittle inclusions, coatings, and interfaces at the microscale under complex loading conditions is critical for microstructure-informed design of materials. In this study, a novel approach has been formulated to decompose the J-integral evaluation of the elastic energy release rate to the three-dimensional stress inte… ▽ More Understanding the local fracture resistance of microstructural features. such as brittle inclusions, coatings, and interfaces at the microscale under complex loading conditions is critical for microstructure-informed design of materials. In this study, a novel approach has been formulated to decompose the J-integral evaluation of the elastic energy release rate to the three-dimensional stress intensity factors directly from experimental measurements of the elastic deformation gradient tensors of the crack field by in situ high (angular) resolution electron backscatter diffraction (HR-EBSD). An exemplar study is presented of a quasi-static crack, inclined to the observed surface, propagating on low index {hkl} planes in a (001) single crystal silicon wafer. △ Less

Submitted 14 September, 2022; v1 submitted 21 June, 2022; originally announced June 2022.

SplitBFT: Improving Byzantine Fault Tolerance Safety Using Trusted Compartments

Authors: Ines Messadi, Markus Horst Becker, Kai Bleeke, Leander Jehl, Sonia Ben Mokhtar, Rüdiger Kapitza

Abstract: Byzantine fault-tolerant agreement (BFT) in a partially synchronous system usually requires 3f + 1 nodes to tolerate f faulty replicas. Due to their high throughput and finality property BFT algorithms build the core of recent permissioned blockchains. As a complex and resource-demanding infrastructure, multiple cloud providers have started offering Blockchain-as-a-Service. This eases the deployme… ▽ More Byzantine fault-tolerant agreement (BFT) in a partially synchronous system usually requires 3f + 1 nodes to tolerate f faulty replicas. Due to their high throughput and finality property BFT algorithms build the core of recent permissioned blockchains. As a complex and resource-demanding infrastructure, multiple cloud providers have started offering Blockchain-as-a-Service. This eases the deployment of permissioned blockchains but places the cloud provider in a central controlling position, thereby questioning blockchains' fault tolerance and decentralization properties and their underlying BFT algorithm. This paper presents SplitBFT, a new way to utilize trusted execution technology (TEEs), such as Intel SGX, to harden the safety and confidentiality guarantees of BFT systems thereby strengthening the trust in could-based deployments of permissioned blockchains. Deviating from standard assumptions, SplitBFT acknowledges that code protected by trusted execution may fail. We address this by splitting and isolating the core logic of BFT protocols into multiple compartments resulting in a more resilient architecture. We apply SplitBFT to the traditional practical byzantine fault tolerance algorithm (PBFT) and evaluate it using SGX. Our results show that SplitBFT adds only a reasonable overhead compared to the non-compartmentalized variant. △ Less

Submitted 24 May, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Dandelion: Certified Approximations of Elementary Functions

Authors: Heiko Becker, Mohit Tekriwal, Eva Darulova, Anastasia Volkova, Jean-Baptiste Jeannin

Abstract: Elementary function operations such as sin and exp cannot in general be computed exactly on today's digital computers, and thus have to be approximated. The standard approximations in library functions typically provide only a limited set of precisions, and are too inefficient for many applications. Polynomial approximations that are customized to a limited input domain and output accuracy can pro… ▽ More Elementary function operations such as sin and exp cannot in general be computed exactly on today's digital computers, and thus have to be approximated. The standard approximations in library functions typically provide only a limited set of precisions, and are too inefficient for many applications. Polynomial approximations that are customized to a limited input domain and output accuracy can provide superior performance. In fact, the Remez algorithm computes the best possible approximation for a given polynomial degree, but has so far not been formally verified. This paper presents Dandelion, an automated certificate checker for polynomial approximations of elementary functions computed with Remez-like algorithms that is fully verified in the HOL4 theorem prover. Dandelion checks whether the difference between a polynomial approximation and its target reference elementary function remains below a given error bound for all inputs in a given constraint. By extracting a verified binary with the CakeML compiler, Dandelion can validate certificates within a reasonable time, fully automating previous manually verified approximations. △ Less

Submitted 29 April, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Enhanced Formulation for Guillotine 2D Cutting Problems

Authors: Henrique Becker, Olinto Araujo, Luciana S. Buriol

Abstract: We advance the state of the art in Mixed-Integer Linear Programming (MILP) formulations for Guillotine 2D Cutting Problems by (i) adapting a previously known reduction to our preprocessing phase and by (ii) enhancing a previous formulation by cutting down its size and symmetries. Our focus is the Guillotine 2D Knapsack Problem with orthogonal and unrestricted cuts, constrained demand, unlimited st… ▽ More We advance the state of the art in Mixed-Integer Linear Programming (MILP) formulations for Guillotine 2D Cutting Problems by (i) adapting a previously known reduction to our preprocessing phase and by (ii) enhancing a previous formulation by cutting down its size and symmetries. Our focus is the Guillotine 2D Knapsack Problem with orthogonal and unrestricted cuts, constrained demand, unlimited stages, and no rotation -- however, the formulation may be adapted to many related problems. The code is available. Concerning the set of 59 instances used to benchmark the original formulation, and summing the statistics for all models generated, the enhanced formulation has only a small fraction of the variables and constraints of the original model (respectively, 3.07% and 8.35%). The enhanced formulation also takes about 4 hours to solve all instances while the original formulation takes 12 hours to solve 53 of them (the other six runs hit a three-hour time limit each). We integrate, to both formulations, a pricing framework proposed for the original formulation; the enhanced formulation keeps a significant advantage in this situation. Finally, in a recently proposed set of 80 harder instances, the enhanced formulation (with and without the pricing framework) found: 22 optimal solutions for the unrestricted problem (5 already known, 17 new); 22 optimal solutions for the restricted problem (all are new and they are not the same 22 of the optimal unrestricted solutions); better lower bounds for 25 instances; better upper bounds for 58 instances. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 22 pages, 6 tables

MSC Class: 68R05 (Primary); 68U99; 05D99; 52B99

The energy spectrum of cosmic rays beyond the turn-down around $10^{17}$ eV as measured with the surface detector of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, 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, J. A. Bellido , et al. (352 additional authors not shown)

Abstract: We present a measurement of the cosmic-ray spectrum above 100\,PeV using the part of the surface detector of the Pierre Auger Observatory that has a spacing of 750~m. An inflection of the spectrum is observed, confirming the presence of the so-called \emph{second-knee} feature. The spectrum is then combined with that of the 1500\,m array to produce a single measurement of the flux, linking this sp… ▽ More We present a measurement of the cosmic-ray spectrum above 100\,PeV using the part of the surface detector of the Pierre Auger Observatory that has a spacing of 750~m. An inflection of the spectrum is observed, confirming the presence of the so-called \emph{second-knee} feature. The spectrum is then combined with that of the 1500\,m array to produce a single measurement of the flux, linking this spectral feature with the three additional breaks at the highest energies. The combined spectrum, with an energy scale set calorimetrically via fluorescence telescopes and using a single detector type, results in the most statistically and systematically precise measurement of spectral breaks yet obtained. These measurements are critical for furthering our understanding of the highest energy cosmic rays. △ Less

Submitted 20 April, 2022; v1 submitted 27 September, 2021; originally announced September 2021.

Comments: Published version, 25 pages, 20 figures, 11 tables Updated April 2022 to fix typo in Table 10

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

Journal ref: Eur. Phys. J. C (2021) 81:966

High-Speed Tunable Microcavities Coupled to Rare-Earth Quantum Emitters

Authors: Kangwei Xia, Fiammetta Sardi, Colin Sauerzapf, Thomas Kornher, Hans-Werner Becker, Zsolt Kis, Laszlo Kovacs, Roman Kolesov, Jörg Wrachtrup

Abstract: Electro-optical control of on-chip photonic devices is an essential tool for efficient integrated photonics. Lithium niobate on insulator (LNOI) is an emerging platform for on-chip photonics due to its large electro-optic coefficient and high nonlinearity [1]. Integrating quantum emitters into LNOI would extend their versatile use in classic photonics to quantum computing and communication [2, 3].… ▽ More Electro-optical control of on-chip photonic devices is an essential tool for efficient integrated photonics. Lithium niobate on insulator (LNOI) is an emerging platform for on-chip photonics due to its large electro-optic coefficient and high nonlinearity [1]. Integrating quantum emitters into LNOI would extend their versatile use in classic photonics to quantum computing and communication [2, 3]. Here, we incorporate single rare-earth ions (REI) quantum emitters in electro-optical tunable lithium niobite (LN) thin films and demonstrate control of LN microcavities coupled to REI over a frequency range of 160 GHz with 5 \mus switching speed. Dynamical control of the cavities enables the modulation of the Purcell enhancement of the REIs with short time constants. Using the Purcell enhancement, we show evidence of detecting single Yb3+ ions in LN cavities. Coupling quantum emitters in fast tunable photonic devices is an efficient method to shape the waveform of the emitter [4]. It also offers a platform to encode quantum information in the integration of a spectral-temporal-spatial domain to achieve high levels of channel multiplexing, as well as an approach to generate deterministic single-photon sources [5, 6]. △ Less

Submitted 1 April, 2021; originally announced April 2021.

Extraction of the Muon Signals Recorded with the Surface Detector of the Pierre Auger Observatory Using Recurrent Neural Networks

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 Pierre Auger Observatory, at present the largest cosmic-ray observatory ever built, is instrumented with a ground array of 1600 water-Cherenkov detectors, known as the Surface Detector (SD). The SD samples the secondary particle content (mostly photons, electrons, positrons and muons) of extensive air showers initiated by cosmic rays with energies ranging from $10^{17}~$eV up to more than… ▽ More The Pierre Auger Observatory, at present the largest cosmic-ray observatory ever built, is instrumented with a ground array of 1600 water-Cherenkov detectors, known as the Surface Detector (SD). The SD samples the secondary particle content (mostly photons, electrons, positrons and muons) of extensive air showers initiated by cosmic rays with energies ranging from $10^{17}~$eV up to more than $10^{20}~$eV. Measuring the independent contribution of the muon component to the total registered signal is crucial to enhance the capability of the Observatory to estimate the mass of the cosmic rays on an event-by-event basis. However, with the current design of the SD, it is difficult to straightforwardly separate the contributions of muons to the SD time traces from those of photons, electrons and positrons. In this paper, we present a method aimed at extracting the muon component of the time traces registered with each individual detector of the SD using Recurrent Neural Networks. We derive the performances of the method by training the neural network on simulations, in which the muon and the electromagnetic components of the traces are known. We conclude this work showing the performance of this method on experimental data of the Pierre Auger Observatory. We find that our predictions agree with the parameterizations obtained by the AGASA collaboration to describe the lateral distributions of the electromagnetic and muonic components of extensive air showers. △ Less

Submitted 1 August, 2021; v1 submitted 22 March, 2021; originally announced March 2021.

Comments: Published version, 23 pages, 15 figures

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

Journal ref: JINST 16 P07016 (2021)

Measurement of the fluctuations in the number of muons in extensive air showers with 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. (343 additional authors not shown)

Abstract: We present the first measurement of the fluctuations in the number of muons in extensive air showers produced by ultrahigh energy cosmic rays. We find that the measured fluctuations are in good agreement with predictions from air shower simulations. This observation provides new insights into the origin of the previously reported deficit of muons in air shower simulations and constrains models of… ▽ More We present the first measurement of the fluctuations in the number of muons in extensive air showers produced by ultrahigh energy cosmic rays. We find that the measured fluctuations are in good agreement with predictions from air shower simulations. This observation provides new insights into the origin of the previously reported deficit of muons in air shower simulations and constrains models of hadronic interactions at ultrahigh energies. Our measurement is compatible with the muon deficit originating from small deviations in the predictions from hadronic interaction models of particle production that accumulate as the showers develop. △ Less

Submitted 27 April, 2021; v1 submitted 15 February, 2021; originally announced February 2021.

Comments: Published version, 11 pages, 3 figures, 1 table, Supplemental material appended

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

Journal ref: Phys. Rev. Lett. 126, 152002 (2021)

Lightning generation in moist convective clouds and constraints on the water abundance in Jupiter

Authors: Yury S. Aglyamov, Jonathan Lunine, Heidi N. Becker, Tristan Guillot, Seran G. Gibbard, Sushil Atreya, Scott J. Bolton, Steven Levin, Shannon T. Brown, Michael H. Wong

Abstract: Recent Juno observations have greatly extended the temporal and spatial coverage of lightning detection on Jupiter. We use these data to constrain a model of moist convection and lightning generation in Jupiter's atmosphere, and derive a roughly solar abundance of water at the base of the water cloud. Shallow lightning, observed by Juno (Becker et al., 2020, Nature, 584, 55-58) and defined as flas… ▽ More Recent Juno observations have greatly extended the temporal and spatial coverage of lightning detection on Jupiter. We use these data to constrain a model of moist convection and lightning generation in Jupiter's atmosphere, and derive a roughly solar abundance of water at the base of the water cloud. Shallow lightning, observed by Juno (Becker et al., 2020, Nature, 584, 55-58) and defined as flashes originating at altitudes corresponding to pressure less than 2 bars, is reproduced, as is lightning at a deeper range of pressures, including those below the water cloud base. It is found that the generation of lightning requires ammonia to stabilize liquid water at altitudes corresponding to sub-freezing temperatures. We find a range of local water abundances in which lightning is possible, including subsolar values of water--consistent with other determinations of deep oxygen abundance. △ Less

Submitted 28 January, 2021; originally announced January 2021.

Comments: 47 pages, 7 figures, to be published in Journal of Geophysical Research - Planets

Design and implementation of the AMIGA embedded system for data acquisition

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. (361 additional authors not shown)

Abstract: The Auger Muon Infill Ground Array (AMIGA) is part of the AugerPrime upgrade of the Pierre Auger Observatory. It consists of particle counters buried 2.3 m underground next to the water-Cherenkov stations that form the 23.5 km$^2$ large infilled array. The reduced distance between detectors in this denser area allows the lowering of the energy threshold for primary cosmic ray reconstruction down t… ▽ More The Auger Muon Infill Ground Array (AMIGA) is part of the AugerPrime upgrade of the Pierre Auger Observatory. It consists of particle counters buried 2.3 m underground next to the water-Cherenkov stations that form the 23.5 km$^2$ large infilled array. The reduced distance between detectors in this denser area allows the lowering of the energy threshold for primary cosmic ray reconstruction down to about 10$^{17}$ eV. At the depth of 2.3 m the electromagnetic component of cosmic ray showers is almost entirely absorbed so that the buried scintillators provide an independent and direct measurement of the air showers muon content. This work describes the design and implementation of the AMIGA embedded system, which provides centralized control, data acquisition and environment monitoring to its detectors. The presented system was firstly tested in the engineering array phase ended in 2017, and lately selected as the final design to be installed in all new detectors of the production phase. The system was proven to be robust and reliable and has worked in a stable manner since its first deployment. △ Less

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

Comments: Accepted for publication at JINST. Published version, 34 pages, 15 figures, 4 tables

Report number: FERMILAB-PUB-21-032-AD-AE-E-TD

Journal ref: JINST 16 T07008 (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)

Lassie: HOL4 Tactics by Example

Authors: Heiko Becker, Nathaniel Bos, Ivan Gavran, Eva Darulova, Rupak Majumdar

Abstract: Proof engineering efforts using interactive theorem proving have yielded several impressive projects in software systems and mathematics. A key obstacle to such efforts is the requirement that the domain expert is also an expert in the low-level details in constructing the proof in a theorem prover. In particular, the user needs to select a sequence of tactics that lead to a successful proof, a ta… ▽ More Proof engineering efforts using interactive theorem proving have yielded several impressive projects in software systems and mathematics. A key obstacle to such efforts is the requirement that the domain expert is also an expert in the low-level details in constructing the proof in a theorem prover. In particular, the user needs to select a sequence of tactics that lead to a successful proof, a task that in general requires knowledge of the exact names and use of a large set of tactics. We present Lassie, a tactic framework for the HOL4 theorem prover that allows individual users to define their own tactic language by example and give frequently used tactics or tactic combinations easier-to-remember names. The core of Lassie is an extensible semantic parser, which allows the user to interactively extend the tactic language through a process of definitional generalization. Defining tactics in Lassie thus does not require any knowledge in implementing custom tactics, while proofs written in Lassie retain the correctness guarantees provided by the HOL4 system. We show through case studies how Lassie can be used in small and larger proofs by novice and more experienced interactive theorem prover users, and how we envision it to ease the learning curve in a HOL4 tutorial. △ Less

Submitted 4 January, 2021; originally announced January 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

Reconstructing the late accretion history of the Moon

Authors: Meng-Hua Zhu, Natalia Artemieva, Alessandro Morbidelli, Qing-Zhu Yin, Harry Becker, Kai Wunnemann

Abstract: The importance of highly siderophile elements (HSEs) to track planetary late accretion has long been recognized. However, the precise nature of the Moon's accretional history remains enigmatic. There exists a significant mismatch of HSE budgets between the Earth and Moon, with the Earth disproportionally accreted far more HSEs than the Moon did. Several scenarios have been proposed to explain this… ▽ More The importance of highly siderophile elements (HSEs) to track planetary late accretion has long been recognized. However, the precise nature of the Moon's accretional history remains enigmatic. There exists a significant mismatch of HSE budgets between the Earth and Moon, with the Earth disproportionally accreted far more HSEs than the Moon did. Several scenarios have been proposed to explain this conundrum, including the delivery of HSEs to Earth by a few big impactors, the accretion of pebble-sized objects on dynamically cold orbits that enhanced the Earth's gravitational focusing factor, and the "sawtooth model" with much reduced impact flux before ~4.10 Gyr. However, most of these models assume a high impactor retention ratio f (fraction of impactor mass retained on the target) for the Moon. Here, we performed a series of impact simulations to quantify the f-value, followed by a Monte Carlo procedure enacting a monotonically decaying impact flux, to compute the mass accreted into lunar crust and mantle over their histories. We found that the average f-value for the Moon's entire impact history is about 3 times lower than previously estimated. Our results indicate that, to match the HSE budget of lunar crust and mantle, the retention of HSEs should have started ~ 4.35 Gyr ago, when most of lunar magma ocean was solidified. Mass accreted prior to 4.35 Gyr must have lost its HSE to the lunar core, presumably during the lunar mantle crystallization. The combination of a low impactor retention ratio and a late retention of HSEs in the lunar mantle provide a realistic explanation for the apparent deficit of Moon's late accreted mass relative to the Earth. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: 31 pages, 3 figures

Journal ref: Nature,2019

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

Frontiers in Planetary Rings Science

Authors: Shawn M. Brooks, Tracy M. Becker, K. Baillié, H. N. Becker, E. T. Bradley, J. E. Colwell, J. N. Cuzzi, I. de Pater, S. Eckert, M. El Moutamid, S. G. Edgington, P. R. Estrada, M. W. Evans, A. Flandes, R. G. French, Á. García, M. K. Gordon, M. M. Hedman, H. -W. Hsu, R. G. Jerousek, E. A. Marouf, B. K. Meinke, P. D. Nicholson, S. H. Pilorz, M. R. Showalter , et al. (3 additional authors not shown)

Abstract: We now know that the outer solar system is host to at least six diverse planetary ring systems, each of which is a scientifically compelling target with the potential to inform us about the evolution, history and even the internal structure of the body it adorns. These diverse ring systems represent a set of distinct local laboratories for understanding the physics and dynamics of planetary disks,… ▽ More We now know that the outer solar system is host to at least six diverse planetary ring systems, each of which is a scientifically compelling target with the potential to inform us about the evolution, history and even the internal structure of the body it adorns. These diverse ring systems represent a set of distinct local laboratories for understanding the physics and dynamics of planetary disks, with applications reaching beyond our Solar System. We highlight the current status of planetary rings science and the open questions before the community to promote continued Earth-based and spacecraft-based investigations into planetary rings. As future spacecraft missions are launched and more powerful telescopes come online in the decades to come, we urge NASA for continued support of investigations that advance our understanding of planetary rings, through research and analysis of data from existing facilities, more laboratory work and specific attention to strong rings science goals during future mission selections. △ Less

Submitted 27 August, 2020; originally announced August 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

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

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

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

An empirical analysis of exact algorithms for the unbounded knapsack problem

Authors: Henrique Becker, Luciana S. Buriol

Abstract: This work presents an empirical analysis of exact algorithms for the unbounded knapsack problem, which includes seven algorithms from the literature, two commercial solvers, and more than ten thousand instances. The terminating step-off, a dynamic programming algorithm from 1966, presented the lowest mean time to solve the most recent benchmark from the literature. The threshold and collective dom… ▽ More This work presents an empirical analysis of exact algorithms for the unbounded knapsack problem, which includes seven algorithms from the literature, two commercial solvers, and more than ten thousand instances. The terminating step-off, a dynamic programming algorithm from 1966, presented the lowest mean time to solve the most recent benchmark from the literature. The threshold and collective dominances are properties of the unbounded knapsack problem first discussed in 1998, and exploited by the current state-of-the-art algorithms. The terminating step-off algorithm did not exploit such dominances, but has an alternative mechanism for dealing with dominances which does not explicitly exploits collective and threshold dominances. Also, the pricing subproblems found when solving hard cutting stock problems with column generation can cause branch-and-bound algorithms to display worst-case times. The authors present a new class of instances which favors the branch-and-bound approach over the dynamic programming approach but do not have high amounts of simple, multiple and collective dominated items. This behaviour illustrates how the definition of hard instances changes among algorithm approachs. The codes used for solving the unbounded knapsack problem and for instance generation are all available online. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: 42 pages, 7 tables (one in conclusions, the rest in appendix), 2 diagrams, 1 algorithm, 8 plots

ACM Class: G.1.6; G.2.1

A New, Deep JVLA Radio Survey of M33

Authors: Richard L. White, Knox S. Long, Robert H. Becker, William P. Blair, David J. Helfand, P. Frank Winkler

Abstract: We have performed new 1.4 GHz and 5 GHz observations of the Local Group galaxy M33 with the Jansky Very Large Array. Our survey has a limiting sensitivity of 20 uJy (4-sigma) and a resolution of 5.9 arcsec (FWHM), corresponding to a spatial resolution of 24 pc at 817 kpc. Using a new multi-resolution algorithm, we have created a catalog of 2875 sources, including 675 with well-determined spectral… ▽ More We have performed new 1.4 GHz and 5 GHz observations of the Local Group galaxy M33 with the Jansky Very Large Array. Our survey has a limiting sensitivity of 20 uJy (4-sigma) and a resolution of 5.9 arcsec (FWHM), corresponding to a spatial resolution of 24 pc at 817 kpc. Using a new multi-resolution algorithm, we have created a catalog of 2875 sources, including 675 with well-determined spectral indices. We detect sources at the position of 319 of the X-ray sources in the Tuellmann et al. (2011) Chandra survey of M33, the majority of which are likely to be background galaxies. The radio source coincident with M33 X-8, the nuclear source, appears to be extended. Along with numerous H II regions or portions of H II region complexes, we detect 155 of the 217 optical supernova remnants included in the lists of Long et al. (2010) and Lee & Lee (2014), making this by far the largest sample of remnants at known distances with multiwavelength coverage. The remnants show a large dispersion in the ratio of radio to X-ray luminosity at a given diameter, a result that challenges the current generation of models for synchrotron radiation evolution in supernova remnants. See http://sundog.stsci.edu/m33 for access to catalogs and images. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 37 pages, 24 figures; Accepted for publication in the Astrophysical Journal; see http://sundog.stsci.edu/m33 for access to catalogs and images

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)

Possible Evidence of the Radio AGN Quenching of Neighboring Galaxies at z $\sim$ 1

Authors: Lu Shen, Adam R. Tomczak, Brian C. Lemaux, Debora Pelliccia, Lori M. Lubin, Neal A. Miller, Serena Perrotta, Christopher D. Fassnacht, Robert H. Becker, Roy R. Gal, Po-Feng Wu, Gordon Squires

Abstract: Using 57 Radio Active Galactic nuclei (RAGN) at 0.55 $\leq$ z $\leq$ 1.3 drawn from five fields of the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey, we study the effect of injection of energy from outbursts of RAGN on their spectroscopically-confirmed neighboring galaxies (SNGs). We observe an elevated fraction of quenched neighbors (fq) within 500 kpc projected r… ▽ More Using 57 Radio Active Galactic nuclei (RAGN) at 0.55 $\leq$ z $\leq$ 1.3 drawn from five fields of the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey, we study the effect of injection of energy from outbursts of RAGN on their spectroscopically-confirmed neighboring galaxies (SNGs). We observe an elevated fraction of quenched neighbors (fq) within 500 kpc projected radius of RAGN in the most dense local environments compared to those of non-RAGN control samples matched to the RAGN population in colour, stellar mass, and local environment at 2$σ$ significance. Further analyses show that there are offsets at similar significance between fqs of RAGN-SNGs and the appropriate control samples for galaxies specifically in cluster environments and those hosted by most massive cluster galaxies, which tentatively suggests that some negative feedback from the RAGN is occurring in these dense environments. In addition, we find that the median radio power of RAGN increases with increasing local overdensity, an effect which may lend itself to the quenching of neighboring galaxies. Furthermore, we find that, in the highest local overdensities, the fq of the sub-sample of lower stellar mass RAGN-SNGs is larger than that of the higher stellar mass RAGN-SNGs sub-sample, which indicates a more pronounced effect from RAGN on lower stellar mass galaxies. We propose a scenario in which RAGN residing within clusters might heat the intracluster medium (ICM) affecting both in situ star formation and any inflowing gas that remains in their neighboring galaxies. △ Less

Submitted 11 December, 2018; originally announced December 2018.

Comments: submitted to MNRAS, 15 pages, 2 tables and 6 figures, comments welcome

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

Large-scale cosmic-ray anisotropies above 4 EeV 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, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker, J. A. Bellido , et al. (352 additional authors not shown)

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

Submitted 5 December, 2018; v1 submitted 10 August, 2018; originally announced August 2018.

Comments: Published version, 16 pages, 7 figures, 8 tables

Report number: FERMILAB-PUB-18-393

Journal ref: ApJ 868 (2018) 4

Observation of inclined EeV air showers with the radio detector of 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, N. Arsene, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Balaceanu, F. Barbato, R. J. Barreira Luz, S. Baur, K. H. Becker , et al. (370 additional authors not shown)

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

Submitted 24 October, 2018; v1 submitted 14 June, 2018; originally announced June 2018.

Comments: Published version, 14 pages, 12 figures

Report number: FERMILAB-PUB-18-259

Journal ref: JCAP10(2018)026

Combining Tools for Optimization and Analysis of Floating-Point Computations

Authors: Heiko Becker, Pavel Pancheckha, Eva Darulova, Zachary Tatlock

Abstract: Recent renewed interest in optimizing and analyzing floating-point programs has lead to a diverse array of new tools for numerical programs. These tools are often complementary, each focusing on a distinct aspect of numerical programming. Building reliable floating point applications typically requires addressing several of these aspects, which makes easy composition essential. This paper describe… ▽ More Recent renewed interest in optimizing and analyzing floating-point programs has lead to a diverse array of new tools for numerical programs. These tools are often complementary, each focusing on a distinct aspect of numerical programming. Building reliable floating point applications typically requires addressing several of these aspects, which makes easy composition essential. This paper describes the composition of two recent floating-point tools: Herbie, which performs accuracy optimization, and Daisy, which performs accuracy verification. We find that the combination provides numerous benefits to users, such as being able to use Daisy to check whether Herbie's unsound optimizations improved the worst-case roundoff error, as well as benefits to tool authors, including uncovering a number of bugs in both tools. The combination also allowed us to compare the different program rewriting techniques implemented by these tools for the first time. The paper lays out a road map for combining other floating-point tools and for surmounting common challenges. △ Less

Submitted 7 May, 2018; originally announced May 2018.

Comments: 13 pages

A realization functor for abelian model categories

Authors: Hanno Becker

Abstract: We study liftings of abelian model structures to categories of chain complexes and construct a realization functor from the derived category of a Grothendieck abelian category equipped with a cofibrantly generated, hereditary abelian model structure to the homotopy category of that model structure. We study liftings of abelian model structures to categories of chain complexes and construct a realization functor from the derived category of a Grothendieck abelian category equipped with a cofibrantly generated, hereditary abelian model structure to the homotopy category of that model structure. △ Less

Submitted 8 March, 2018; originally announced March 2018.

Comments: 17 pages, comments welcome

Indication of anisotropy in arrival directions of ultra-high-energy cosmic rays through comparison to the flux pattern of extragalactic gamma-ray sources

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, G. Avila, A. M. Badescu, A. Balaceanu, F. Barbato, R. J. Barreira Luz, J. J. Beatty, K. H. Becker, J. A. Bellido , et al. (368 additional authors not shown)

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

Submitted 6 February, 2018; v1 submitted 18 January, 2018; originally announced January 2018.

Comments: Published version, 10 pages, 3 figures, 2 tables

Report number: FERMILAB-PUB-17-358

Journal ref: Aab. et al., The Astrophysical Journal Letters, 853:L29 (2018)

Inferences on Mass Composition and Tests of Hadronic Interactions from 0.3 to 100 EeV using the water-Cherenkov Detectors 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, F. Barbato, R. J. Barreira Luz , et al. (381 additional authors not shown)

Abstract: We present a new method for probing the hadronic interaction models at ultra-high energy and extracting details about mass composition. This is done using the time profiles of the signals recorded with the water-Cherenkov detectors of the Pierre Auger Observatory. The profiles arise from a mix of the muon and electromagnetic components of air-showers. Using the risetimes of the recorded signals we… ▽ More We present a new method for probing the hadronic interaction models at ultra-high energy and extracting details about mass composition. This is done using the time profiles of the signals recorded with the water-Cherenkov detectors of the Pierre Auger Observatory. The profiles arise from a mix of the muon and electromagnetic components of air-showers. Using the risetimes of the recorded signals we define a new parameter, which we use to compare our observations with predictions from simulations. We find, firstly, inconsistencies between our data and predictions over a greater energy range and with substantially more events than in previous studies. Secondly, by calibrating the new parameter with fluorescence measurements from observations made at the Auger Observatory, we can infer the depth of shower maximum for a sample of over 81,000 events extending from 0.3 EeV to over 100 EeV. Above 30 EeV, the sample is nearly fourteen times larger than currently available from fluorescence measurements and extending the covered energy range by half a decade. The energy dependence of the average depth of shower maximum is compared to simulations and interpreted in terms of the mean of the logarithmic mass. We find good agreement with previous work and extend the measurement of the mean depth of shower maximum to greater energies than before, reducing significantly the statistical uncertainty associated with the inferences about mass composition. △ Less

Submitted 19 October, 2017; originally announced October 2017.

Comments: 21 pages; 18 figures; version accepted for publication in Phys. Rev. D

Report number: FERMILAB-PUB-17-357

Search for High-energy Neutrinos from Binary Neutron Star Merger GW170817 with ANTARES, IceCube, and the Pierre Auger Observatory

Authors: A. Albert, M. Andre, M. Anghinolfi, M. Ardid, J. -J. Aubert, J. Aublin, T. Avgitas, B. Baret, J. Barrios-Marti, S. Basa, B. Belhorma, V. Bertin, S. Biagi, R. Bormuth, S. Bourret, M. C. Bouwhuis, H. Branzacs, R. Bruijn, J. Brunner, J. Busto, A. Capone, L. Caramete, J. Carr, S. Celli, R. Cherkaoui El Moursli , et al. (1916 additional authors not shown)

Abstract: The Advanced LIGO and Advanced Virgo observatories recently discovered gravitational waves from a binary neutron star inspiral. A short gamma-ray burst (GRB) that followed the merger of this binary was also recorded by the Fermi Gamma-ray Burst Monitor (Fermi-GBM), and the Anticoincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory (INTEGRAL), indicating par… ▽ More The Advanced LIGO and Advanced Virgo observatories recently discovered gravitational waves from a binary neutron star inspiral. A short gamma-ray burst (GRB) that followed the merger of this binary was also recorded by the Fermi Gamma-ray Burst Monitor (Fermi-GBM), and the Anticoincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory (INTEGRAL), indicating particle acceleration by the source. The precise location of the event was determined by optical detections of emission following the merger. We searched for high-energy neutrinos from the merger in the GeV--EeV energy range using the ANTARES, IceCube, and Pierre Auger Observatories. No neutrinos directionally coincident with the source were detected within $\pm500$ s around the merger time. Additionally, no MeV neutrino burst signal was detected coincident with the merger. We further carried out an extended search in the direction of the source for high-energy neutrinos within the 14-day period following the merger, but found no evidence of emission. We used these results to probe dissipation mechanisms in relativistic outflows driven by the binary neutron star merger. The non-detection is consistent with model predictions of short GRBs observed at a large off-axis angle. △ Less

Submitted 9 November, 2017; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: 22 pages, 2 figures

Report number: LIGO-P1700344

Observation of a Large-scale Anisotropy in the Arrival Directions of Cosmic Rays above $8 \times 10^{18}$ eV

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. (382 additional authors not shown)

Abstract: Cosmic rays are atomic nuclei arriving from outer space that reach the highest energies observed in nature. Clues to their origin come from studying the distribution of their arrival directions. Using $3 \times 10^4$ cosmic rays above $8 \times 10^{18}$ electron volts, recorded with the Pierre Auger Observatory from a total exposure of 76,800 square kilometers steradian year, we report an anisotro… ▽ More Cosmic rays are atomic nuclei arriving from outer space that reach the highest energies observed in nature. Clues to their origin come from studying the distribution of their arrival directions. Using $3 \times 10^4$ cosmic rays above $8 \times 10^{18}$ electron volts, recorded with the Pierre Auger Observatory from a total exposure of 76,800 square kilometers steradian year, we report an anisotropy in the arrival directions. The anisotropy, detected at more than the 5.2$σ$ level of significance, can be described by a dipole with an amplitude of $6.5_{-0.9}^{+1.3}$% towards right ascension $α_{d} = 100 \pm 10$ degrees and declination $δ_{d} = -24_{-13}^{+12}$ degrees. That direction indicates an extragalactic origin for these ultra-high energy particles. △ Less

Submitted 21 September, 2017; originally announced September 2017.

Comments: 19 pages (with supplementary material), 8 figures

Report number: FERMILAB-PUB-17-354

Journal ref: Science 357 (22 September 2017) 1266

Spectral Calibration of the Fluorescence Telescopes 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, F. Barbato, R. J. Barreira Luz , et al. (381 additional authors not shown)

Abstract: We present a novel method to measure precisely the relative spectral response of the fluorescence telescopes of the Pierre Auger Observatory. We used a portable light source based on a xenon flasher and a monochromator to measure the relative spectral efficiencies of eight telescopes in steps of 5 nm from 280 nm to 440 nm. Each point in a scan had approximately 2 nm FWHM out of the monochromator.… ▽ More We present a novel method to measure precisely the relative spectral response of the fluorescence telescopes of the Pierre Auger Observatory. We used a portable light source based on a xenon flasher and a monochromator to measure the relative spectral efficiencies of eight telescopes in steps of 5 nm from 280 nm to 440 nm. Each point in a scan had approximately 2 nm FWHM out of the monochromator. Different sets of telescopes in the observatory have different optical components, and the eight telescopes measured represent two each of the four combinations of components represented in the observatory. We made an end-to-end measurement of the response from different combinations of optical components, and the monochromator setup allowed for more precise and complete measurements than our previous multi-wavelength calibrations. We find an overall uncertainty in the calibration of the spectral response of most of the telescopes of 1.5% for all wavelengths; the six oldest telescopes have larger overall uncertainties of about 2.2%. We also report changes in physics measureables due to the change in calibration, which are generally small. △ Less

Submitted 2 October, 2017; v1 submitted 5 September, 2017; originally announced September 2017.

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

Report number: FERMILAB-PUB-17-355

Journal ref: Astroparticle Physics 95 (2017) 44-56

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

The Properties of Radio Galaxies and the Effect of Environment in Large Scale Structures at $z\sim1$

Authors: Lu Shen, Neal A. Miller, Brian C. Lemaux, Adam R. Tomczak, Lori M. Lubin, Nicholas Rumbaugh, Christopher D. Fassnacht, Robert H. Becker, Roy R. Gal, Po-Feng. Wu, Gordon Squires

Abstract: In this study we investigate 89 radio galaxies that are spectroscopically-confirmed to be members of five large scale structures in the redshift range of $0.65 \le z \le 0.96$. Based on a two-stage classification scheme, the radio galaxies are classified into three sub-classes: active galactic nucleus (AGN), hybrid, and star-forming galaxy (SFG). We study the properties of the three radio sub-clas… ▽ More In this study we investigate 89 radio galaxies that are spectroscopically-confirmed to be members of five large scale structures in the redshift range of $0.65 \le z \le 0.96$. Based on a two-stage classification scheme, the radio galaxies are classified into three sub-classes: active galactic nucleus (AGN), hybrid, and star-forming galaxy (SFG). We study the properties of the three radio sub-classes and their global and local environmental preferences. We find AGN hosts are the most massive population and exhibit quiescence in their star-formation activity. The SFG population has a comparable stellar mass to those hosting a radio AGN but are unequivocally powered by star formation. Hybrids, though selected as an intermediate population in our classification scheme, were found in almost all analyses to be a unique type of radio galaxies rather than a mixture of AGN and SFGs. They are dominated by a high-excitation radio galaxy (HERG) population. We discuss environmental effects and scenarios for each sub-class. AGN tend to be preferentially located in locally dense environments and in the cores of clusters/groups, with these preferences persisting when comparing to galaxies of similar colour and stellar mass, suggesting that their activity may be ignited in the cluster/group virialized core regions. Conversely, SFGs exhibit a strong preference for intermediate-density global environments, suggesting that dusty starbursting activity in LSSs is largely driven by galaxy-galaxy interactions and merging. △ Less

Submitted 7 August, 2017; originally announced August 2017.

Comments: 28 pages, 10 figures, accepted to MNRAS

A Verified Certificate Checker for Finite-Precision Error Bounds in Coq and HOL4

Authors: Heiko Becker, Nikita Zyuzin, Raphael Monat, Eva Darulova, Magnus O. Myreen, Anthony Fox

Abstract: Being able to soundly estimate roundoff errors of finite-precision computations is important for many applications in embedded systems and scientific computing. Due to the discrepancy between continuous reals and discrete finite-precision values, automated static analysis tools are highly valuable to estimate roundoff errors. The results, however, are only as correct as the implementations of the… ▽ More Being able to soundly estimate roundoff errors of finite-precision computations is important for many applications in embedded systems and scientific computing. Due to the discrepancy between continuous reals and discrete finite-precision values, automated static analysis tools are highly valuable to estimate roundoff errors. The results, however, are only as correct as the implementations of the static analysis tools. This paper presents a formally verified and modular tool which fully automatically checks the correctness of finite-precision roundoff error bounds encoded in a certificate. We present implementations of certificate generation and checking for both Coq and HOL4 and evaluate it on a number of examples from the literature. The experiments use both in-logic evaluation of Coq and HOL4, and execution of extracted code outside of the logics: we benchmark Coq extracted unverified OCaml code and a CakeML-generated verified binary. △ Less

Submitted 23 August, 2018; v1 submitted 7 July, 2017; originally announced July 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