Showing 1–50 of 108 results for author: Day, J

Kitaev Interactions Through an Extended Superexchange Pathway in the jeff = 1/2 Ru3+ Honeycomb Magnet, RuP3SiO11

Authors: Aly H. Abdeldaim, Hlynur Gretarsson, Sarah J. Day, M. Duc Le, Gavin B. G. Stenning, Pascal Manuel, Robin S. Perry, Alexander A. Tsirlin, Gøran J. Nilsen, Lucy Clark

Abstract: Magnetic materials are composed of the simple building blocks of magnetic moments on a crystal lattice that interact via short-range magnetic exchange interactions. Yet from these simple building blocks emerges a remarkable diversity of magnetic states. Some of these, such as ferromagnetism, are familiar in our everyday lives, while others reveal the deep quantum mechanical origins of magnetism. A… ▽ More Magnetic materials are composed of the simple building blocks of magnetic moments on a crystal lattice that interact via short-range magnetic exchange interactions. Yet from these simple building blocks emerges a remarkable diversity of magnetic states. Some of these, such as ferromagnetism, are familiar in our everyday lives, while others reveal the deep quantum mechanical origins of magnetism. A prime example of the latter are quantum spin liquid (QSL) states in which -- unlike in a ferromagnet where magnetic moments are driven by their exchange interactions to adopt long-range order -- magnetic moments remain disordered at low temperatures but are simultaneously correlated over long length scales through quantum entanglement. A particularly promising theoretical model of a QSL is the Kitaev model, composed of unusual bond-dependent exchange interactions between magnetic moments on a honeycomb lattice. However, the Kitaev QSL is extremely challenging to realise experimentally as it is unstable to competing exchange interactions and crystal lattice perturbations that inevitably arise in real materials. This makes it essential to understand the relationship between the structure and interactions that may give rise to Kitaev interactions in new candidate materials. Here we show that the material requirements for the Kitaev QSL survive for an extended pseudo-edge-sharing superexchange pathway of Ru3+ 4d5 octahedra within the honeycomb layers of the inorganic framework solid, RuP3SiO11. Through materials synthesis and structural characterisation, resonant inelastic X-ray and neutron scattering experiments, we confirm the requisite jeff = 1/2 state of Ru3+ in RuP3SiO11 and resolve the hierarchy of exchange interactions that provide experimental access to an otherwise unexplored region of the extended Kitaev phase diagram. △ Less

Submitted 28 March, 2024; originally announced March 2024.

Comments: 38 Pages, 5 figures, 7 supporting information figures

Learning 3D object-centric representation through prediction

Authors: John Day, Tushar Arora, Jirui Liu, Li Erran Li, Ming Bo Cai

Abstract: As part of human core knowledge, the representation of objects is the building block of mental representation that supports high-level concepts and symbolic reasoning. While humans develop the ability of perceiving objects situated in 3D environments without supervision, models that learn the same set of abilities with similar constraints faced by human infants are lacking. Towards this end, we de… ▽ More As part of human core knowledge, the representation of objects is the building block of mental representation that supports high-level concepts and symbolic reasoning. While humans develop the ability of perceiving objects situated in 3D environments without supervision, models that learn the same set of abilities with similar constraints faced by human infants are lacking. Towards this end, we developed a novel network architecture that simultaneously learns to 1) segment objects from discrete images, 2) infer their 3D locations, and 3) perceive depth, all while using only information directly available to the brain as training data, namely: sequences of images and self-motion. The core idea is treating objects as latent causes of visual input which the brain uses to make efficient predictions of future scenes. This results in object representations being learned as an essential byproduct of learning to predict. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 21 pages, 11 figures. Project webpage can be found at https://jday54.github.io/opple_site/

ACM Class: I.2.10; I.4.8; I.4.6; I.4.10; I.2.6

Layered and Staged Monte Carlo Tree Search for SMT Strategy Synthesis

Authors: Zhengyang Lu, Stefan Siemer, Piyush Jha, Joel Day, Florin Manea, Vijay Ganesh

Abstract: Modern SMT solvers, such as Z3, offer user-controllable strategies, enabling users to tailor solving strategies for their unique set of instances, thus dramatically enhancing solver performance for their use case. However, this approach of strategy customization presents a significant challenge: handcrafting an optimized strategy for a class of SMT instances remains a complex and demanding task fo… ▽ More Modern SMT solvers, such as Z3, offer user-controllable strategies, enabling users to tailor solving strategies for their unique set of instances, thus dramatically enhancing solver performance for their use case. However, this approach of strategy customization presents a significant challenge: handcrafting an optimized strategy for a class of SMT instances remains a complex and demanding task for both solver developers and users alike. In this paper, we address this problem of automatic SMT strategy synthesis via a novel Monte Carlo Tree Search (MCTS) based method. Our method treats strategy synthesis as a sequential decision-making process, whose search tree corresponds to the strategy space, and employs MCTS to navigate this vast search space. The key innovations that enable our method to identify effective strategies, while keeping costs low, are the ideas of layered and staged MCTS search. These novel heuristics allow for a deeper and more efficient exploration of the strategy space, enabling us to synthesize more effective strategies than the default ones in state-of-the-art (SOTA) SMT solvers. We implement our method, dubbed Z3alpha, as part of the Z3 SMT solver. Through extensive evaluations across six important SMT logics, Z3alpha demonstrates superior performance compared to the SOTA synthesis tool FastSMT, the default Z3 solver, and the CVC5 solver on most benchmarks. Remarkably, on a challenging QF_BV benchmark set, Z3alpha solves 42.7% more instances than the default strategy in the Z3 SMT solver. △ Less

Submitted 30 April, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: Accepted at IJCAI 2024

Estimating Primary Substation Boundaries and the Value of Mapping the Electrical Network Infrastructure of Great Britain

Authors: Joseph Day, I. A. Grant Wilson, Daniel L. Donaldson, Edward Barbour, Bruno Cárdenas, Christopher R. Jones, Andrew J. Urquhart, Seamus D. Garvey

Abstract: Localised data aggregation in many countries including Great Britain (GB) is typically done to a geographical level with polygon boundaries that have a robust and trusted governance system in place. At a minimum this will mean there is confidence in a process to create a set of polygons that have unique identifiers coupled to geographical areas, and the ability to have these updated through a defi… ▽ More Localised data aggregation in many countries including Great Britain (GB) is typically done to a geographical level with polygon boundaries that have a robust and trusted governance system in place. At a minimum this will mean there is confidence in a process to create a set of polygons that have unique identifiers coupled to geographical areas, and the ability to have these updated through a defined code of practice. Examples found across many countries are in the delivery of post, such as postcodes and zip codes, and of the definition of census areas and municipal boundaries. The confidence in these boundaries allows different data to be aggregated by third parties, which itself provides greater levels of data over comparable geographical areas to enhance wider analysis and decision making. Here we combine publicly available datasets published from the six regional electricity Distribution Network Operators of GB to produce a new geospatial dataset with 4436 unique polygons defining the areas served by electrical primary substations. An example is also presented of the use of these polygons to link postcode level open government datasets on domestic energy consumption (2015-2020) from the Department of Energy Security and Net Zero (DESNZ). This results in another dataset with energy statistics aggregated to the geographical areas served by each primary substation across Great Britain. Therefore, we believe there is a compelling argument for countries to set up processes to create and update polygons that have a meaningful relationship to energy systems. This would allow more accurate energy systems analysis to be performed, ultimately leading to an accelerated or potentially lower cost transition to a net-zero world. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 58 pages, 26 figures

Knowledge Engineering for Wind Energy

Authors: Yuriy Marykovskiy, Thomas Clark, Justin Day, Marcus Wiens, Charles Henderson, Julian Quick, Imad Abdallah, Anna Maria Sempreviva, Jean-Paul Calbimonte, Eleni Chatzi, Sarah Barber

Abstract: With the rapid evolution of the wind energy sector, there is an ever-increasing need to create value from the vast amounts of data made available both from within the domain, as well as from other sectors. This article addresses the challenges faced by wind energy domain experts in converting data into domain knowledge, connecting and integrating it with other sources of knowledge, and making it a… ▽ More With the rapid evolution of the wind energy sector, there is an ever-increasing need to create value from the vast amounts of data made available both from within the domain, as well as from other sectors. This article addresses the challenges faced by wind energy domain experts in converting data into domain knowledge, connecting and integrating it with other sources of knowledge, and making it available for use in next generation artificially intelligent systems. To this end, this article highlights the role that knowledge engineering can play in the process of digital transformation of the wind energy sector. It presents the main concepts underpinning Knowledge-Based Systems and summarises previous work in the areas of knowledge engineering and knowledge representation in a manner that is relevant and accessible to domain experts. A systematic analysis of the current state-of-the-art on knowledge engineering in the wind energy domain is performed, with available tools put into perspective by establishing the main domain actors and their needs and identifying key problematic areas. Finally, guidelines for further development and improvement are provided. △ Less

Submitted 1 October, 2023; originally announced October 2023.

Journal ref: Wind Energ. Sci. 9 (2024) 883-917

Equilibrium measures for two-sided shift spaces via dimension theory

Authors: Vaughn Climenhaga, Jason Day

Abstract: Given a two-sided shift space on a finite alphabet and a continuous potential function, we give conditions under which an equilibrium measure can be described using a construction analogous to Hausdorff measure that goes back to the work of Bowen. This construction was previously applied to smooth uniformly and partially hyperbolic systems by the first author, Pesin, and Zelerowicz. Our results he… ▽ More Given a two-sided shift space on a finite alphabet and a continuous potential function, we give conditions under which an equilibrium measure can be described using a construction analogous to Hausdorff measure that goes back to the work of Bowen. This construction was previously applied to smooth uniformly and partially hyperbolic systems by the first author, Pesin, and Zelerowicz. Our results here apply to all subshifts of finite type and Hölder continuous potentials, but extend beyond this setting, and we also apply them to shift spaces with synchronizing words. △ Less

Submitted 22 May, 2024; v1 submitted 20 September, 2023; originally announced September 2023.

Comments: 41 pages

Quantifying the risk of workplace COVID-19 clusters in terms of commuter, workplace, and population characteristics

Authors: Christopher E. Overton, Rachel Abbey, Tarrion Baird, Rachel Christie, Owen Daniel, Julie Day, Matthew Gittins, Owen Jones, Robert Paton, Maria Tang, Tom Ward, Jack Wilkinson, Camilla Woodrow-Hill, Tim Aldridge, Yiqun Chen

Abstract: Objectives: To identify and quantify risk factors that contribute to clusters of COVID-19 in the workplace. Methods: We identified clusters of COVID-19 cases in the workplace and investigated the characteristics of the individuals, the workplaces, the areas they work, and the methods of commute to work, through data linkages based on Middle Layer Super Output Areas (MSOAs) in England between 20/… ▽ More Objectives: To identify and quantify risk factors that contribute to clusters of COVID-19 in the workplace. Methods: We identified clusters of COVID-19 cases in the workplace and investigated the characteristics of the individuals, the workplaces, the areas they work, and the methods of commute to work, through data linkages based on Middle Layer Super Output Areas (MSOAs) in England between 20/06/2021 and 20/02/2022. We estimated associations between potential risk factors and workplace clusters, adjusting for plausible confounders identified using a Directed Acyclic Graph (DAG). Results: For most industries, increased physical proximity in the workplace was associated with increased risk of COVID-19 clusters, while increased vaccination was associated with reduced risk. Commuter demographic risk factors varied across industry, but for the majority of industries, a higher proportion of black/african/caribbean ethnicities, and living in deprived areas, was associated with increased cluster risk. A higher proportion of commuters in the 60-64 age group was associated with reduced cluster risk. There were significant associations between gender, work commute methods, and staff contract type with cluster risk, but these were highly variable across industries. Conclusions: This study has used novel national data linkages to identify potential risk factors of workplace COVID-19 clusters, including possible protective effects of vaccination and increased physical distance at work. The same methodological approach can be applied to wider occupational and environmental health research. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Planetary Exploration Horizon 2061 Report Chapter 5: Enabling technologies for planetary exploration

Authors: Manuel Grande, Linli Guo, Michel Blanc, Advenit Makaya, Sami Asmar, David Atkinson, Anne Bourdon, Pascal Chabert, Steve Chien, John Day, Alberto G. Fairen, Anthony Freeman, Antonio Genova, Alain Herique, Wlodek Kofman, Joseph Lazio, Olivier Mousis, Gian Gabriele Ori, Victor Parro, Robert Preston, Jose A Rodriguez-Manfredi, Veerle Sterken, Keith Stephenson, Joshua Vander Hook, Hunter Waite , et al. (1 additional authors not shown)

Abstract: The main objective of this chapter is to present an overview of the different areas of key technologies that will be needed to fly the technically most challenging of the representative missions identified in chapter 4 (the Pillar 2 Horizon 2061 report). It starts with a description of the future scientific instruments which will address the key questions of Horizon 2061 described in chapter 3 (th… ▽ More The main objective of this chapter is to present an overview of the different areas of key technologies that will be needed to fly the technically most challenging of the representative missions identified in chapter 4 (the Pillar 2 Horizon 2061 report). It starts with a description of the future scientific instruments which will address the key questions of Horizon 2061 described in chapter 3 (the Pillar 1 Horizon 2061 report) and the new technologies that the next generations of space instruments will require (section 2). From there, the chapter follows the line of logical development and implementation of a planetary mission: section 3 describes some of the novel mission architectures that will be needed and how they will articulate interplanetary spacecraft and science platforms; section 4 summarizes the system-level technologies needed: power, propulsion, navigation, communication, advanced autonomy on board planetary spacecraft; section 5 describes the diversity of specialized science platforms that will be needed to survive, operate and return scientific data from the extreme environments that future missions will target; section 6 describes the new technology developments that will be needed for long-duration missions and semi-permanent settlements; finally, section 7 attempts to anticipate on the disruptive technologies that should emerge and progressively prevail in the decades to come to meet the long-term needs of future planetary missions. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 100 pages, 23 figures, Horizon 2061 is a science-driven, foresight exercise, for future scientific investigations

Nowcasting the 2022 mpox outbreak in England

Authors: Christopher E. Overton, Sam Abbott, Rachel Christie, Fergus Cumming, Julie Day, Owen Jones, Rob Paton, Charlie Turner, Thomas Ward

Abstract: In May 2022, a cluster of mpox cases were detected in the UK that could not be traced to recent travel history from an endemic region. Over the coming months, the outbreak grew, with over 3000 total cases reported in the UK, and similar outbreaks occurring worldwide. These outbreaks appeared linked to sexual contact networks between gay, bisexual and other men who have sex with men. Following the… ▽ More In May 2022, a cluster of mpox cases were detected in the UK that could not be traced to recent travel history from an endemic region. Over the coming months, the outbreak grew, with over 3000 total cases reported in the UK, and similar outbreaks occurring worldwide. These outbreaks appeared linked to sexual contact networks between gay, bisexual and other men who have sex with men. Following the COVID-19 pandemic, local health systems were strained, and therefore effective surveillance for mpox was essential for managing public health policy. However, the mpox outbreak in the UK was characterised by substantial delays in the reporting of the symptom onset date and specimen collection date for confirmed positive cases. These delays led to substantial backfilling in the epidemic curve, making it challenging to interpret the epidemic trajectory in real-time. Many nowcasting models exist to tackle this challenge in epidemiological data, but these lacked sufficient flexibility. We have developed a novel nowcasting model using generalised additive models to correct the mpox epidemic curve in England, and provide real-time characteristics of the state of the epidemic, including the real-time growth rate. This model benefited from close collaboration with individuals involved in collecting and processing the data, enabling temporal changes in the reporting structure to be built into the model, which improved the robustness of the nowcasts generated. △ Less

Submitted 17 February, 2023; originally announced February 2023.

Denoising neural networks for magnetic resonance spectroscopy

Authors: Natalie Klein, Amber J. Day, Harris Mason, Michael W. Malone, Sinead A. Williamson

Abstract: In many scientific applications, measured time series are corrupted by noise or distortions. Traditional denoising techniques often fail to recover the signal of interest, particularly when the signal-to-noise ratio is low or when certain assumptions on the signal and noise are violated. In this work, we demonstrate that deep learning-based denoising methods can outperform traditional techniques w… ▽ More In many scientific applications, measured time series are corrupted by noise or distortions. Traditional denoising techniques often fail to recover the signal of interest, particularly when the signal-to-noise ratio is low or when certain assumptions on the signal and noise are violated. In this work, we demonstrate that deep learning-based denoising methods can outperform traditional techniques while exhibiting greater robustness to variation in noise and signal characteristics. Our motivating example is magnetic resonance spectroscopy, in which a primary goal is to detect the presence of short-duration, low-amplitude radio frequency signals that are often obscured by strong interference that can be difficult to separate from the signal using traditional methods. We explore various deep learning architecture choices to capture the inherently complex-valued nature of magnetic resonance signals. On both synthetic and experimental data, we show that our deep learning-based approaches can exceed performance of traditional techniques, providing a powerful new class of methods for analysis of scientific time series data. △ Less

Submitted 31 October, 2022; originally announced November 2022.

Comments: 5 pages with appendix

A Generic Information Extraction System for String Constraints

Authors: Joel D. Day, Adrian Kröger, Mitja Kulczynski, Florin Manea, Dirk Nowotka, Danny Bøgsted Poulsen

Abstract: String constraint solving, and the underlying theory of word equations, are highly interesting research topics both for practitioners and theoreticians working in the wide area of satisfiability modulo theories. As string constraint solving algorithms, a.k.a. string solvers, gained a more prominent role in the formal analysis of string-heavy programs, especially in connection to symbolic code exec… ▽ More String constraint solving, and the underlying theory of word equations, are highly interesting research topics both for practitioners and theoreticians working in the wide area of satisfiability modulo theories. As string constraint solving algorithms, a.k.a. string solvers, gained a more prominent role in the formal analysis of string-heavy programs, especially in connection to symbolic code execution and security protocol verification, we can witness an ever-growing number of benchmarks collecting string solving instances from real-world applications as well as an ever-growing need for more efficient and reliable solvers, especially for the aforementioned real-world instances. Thus, it seems that the string solving area (and the developers, theoreticians, and end-users active in it) could greatly benefit from a better understanding and processing of the existing string solving benchmarks. In this context, we propose SMTQUERY: an SMT-LIB benchmark analysis tool for string constraints. SMTQUERY is implemented in Python 3, and offers a collection of analysis and information extraction tools for a comprehensive data base of string benchmarks (presented in SMT-LIB format), based on an SQL-centred language called QLANG. △ Less

Submitted 18 August, 2022; originally announced August 2022.

Cosmological implications of photon-flux upper limits at ultra-high energies in scenarios of Planckian-interacting massive particles for dark matter

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

Abstract: Using the data of the Pierre Auger Observatory, we report on a search for signatures that would be suggestive of super-heavy particles decaying in the Galactic halo. From the lack of signal, we present upper limits for different energy thresholds above ${\gtrsim}10^8$\,GeV on the secondary by-product fluxes expected from the decay of the particles. Assuming that the energy density of these super-h… ▽ More Using the data of the Pierre Auger Observatory, we report on a search for signatures that would be suggestive of super-heavy particles decaying in the Galactic halo. From the lack of signal, we present upper limits for different energy thresholds above ${\gtrsim}10^8$\,GeV on the secondary by-product fluxes expected from the decay of the particles. Assuming that the energy density of these super-heavy particles matches that of dark matter observed today, we translate the upper bounds on the particle fluxes into tight constraints on the couplings governing the decay process as a function of the particle mass. Instantons, which are non-perturbative solutions to Yang-Mills equations, can give rise to decay channels otherwise forbidden and transform stable particles into meta-stable ones. Assuming such instanton-induced decay processes, we derive a bound on the reduced coupling constant of gauge interactions in the dark sector: $α_X \lesssim 0.09$, for $10^{9} \lesssim M_X/\text{GeV} < 10^{19}$. Conversely, we obtain that, for instance, a reduced coupling constant $α_X = 0.09$ excludes masses $M_X \gtrsim 3\times 10^{13}~$GeV. In the context of dark matter production from gravitational interactions alone during the reheating epoch, we derive constraints on the parameter space that involves, in addition to $M_X$ and $α_X$, the Hubble rate at the end of inflation, the reheating efficiency, and the non-minimal coupling of the Higgs with curvature. △ Less

Submitted 15 December, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 15 pages, 8 figures, Accompanying paper of arXiv:2203.08854, accepted for publication in Phys. Rev. D

Journal ref: Phys.Rev.D 107 (2023) 4, 042002

Subsequences With Gap Constraints: Complexity Bounds for Matching and Analysis Problems

Authors: Joel D. Day, Maria Kosche, Florin Manea, Markus L. Schmid

Abstract: We consider subsequences with gap constraints, i.e., length-k subsequences p that can be embedded into a string w such that the induced gaps (i.e., the factors of w between the positions to which p is mapped to) satisfy given gap constraints $gc = (C_1, C_2, ..., C_{k-1})$; we call p a gc-subsequence of w. In the case where the gap constraints gc are defined by lower and upper length bounds… ▽ More We consider subsequences with gap constraints, i.e., length-k subsequences p that can be embedded into a string w such that the induced gaps (i.e., the factors of w between the positions to which p is mapped to) satisfy given gap constraints $gc = (C_1, C_2, ..., C_{k-1})$; we call p a gc-subsequence of w. In the case where the gap constraints gc are defined by lower and upper length bounds $C_i = (L^-_i, L^+_i) \in \mathbb{N}^2$ and/or regular languages $C_i \in REG$, we prove tight (conditional on the orthogonal vectors (OV) hypothesis) complexity bounds for checking whether a given p is a gc-subsequence of a string w. We also consider the whole set of all gc-subsequences of a string, and investigate the complexity of the universality, equivalence and containment problems for these sets of gc-subsequences. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Arrival Directions of Cosmic Rays above 32 EeV from Phase One of the Pierre Auger Observatory

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

Abstract: A promising energy range to look for angular correlation between cosmic rays of extragalactic origin and their sources is at the highest energies, above few tens of EeV ($1\:{\rm EeV}\equiv 10^{18}\:$eV). Despite the flux of these particles being extremely low, the area of ${\sim}\:3{,}000 \: \text{km}^2$ covered at the Pierre Auger Observatory, and the 17-year data-taking period of the Phase 1 of… ▽ More A promising energy range to look for angular correlation between cosmic rays of extragalactic origin and their sources is at the highest energies, above few tens of EeV ($1\:{\rm EeV}\equiv 10^{18}\:$eV). Despite the flux of these particles being extremely low, the area of ${\sim}\:3{,}000 \: \text{km}^2$ covered at the Pierre Auger Observatory, and the 17-year data-taking period of the Phase 1 of its operations, have enabled us to measure the arrival directions of more than 2,600 ultra-high energy cosmic rays above $32\:\text{EeV}$. We publish this data set, the largest available at such energies from an integrated exposure of $122{,}000 \: \text{km}^2\:\text{sr}\:\text{yr}$, and search it for anisotropies over the $3.4π$ steradians covered with the Observatory. Evidence for a deviation in excess of isotropy at intermediate angular scale, with ${\sim}\:15^\circ$ Gaussian spread or ${\sim}\:25^\circ$ top-hat radius, is obtained at the $4\:σ$ significance level for cosmic-ray energies above ${\sim}\:40\:\text{EeV}$. △ Less

Submitted 5 September, 2022; v1 submitted 27 June, 2022; originally announced June 2022.

Comments: Published version, 33 pages, 10 figures. Dataset and analysis software available at https://doi.org/10.5281/zenodo.6504276

Report number: FERMILAB-PUB-22-491-AD-PPD-SCD-TD

Journal ref: ApJ 935 170 (2022)

Investigating Hadronic Interactions at Ultra-High Energies with the Pierre Auger Observatory

Authors: Isabel Goos, :, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova , et al. (352 additional authors not shown)

Abstract: The development of an extensive air shower depends not only on the nature of the primary ultra-high-energy cosmic ray but also on the properties of the hadronic interactions. For energies above those achievable in human-made accelerators, hadronic interactions are only accessible through the studies of extensive air showers, which can be measured at the Pierre Auger Observatory. With its hybrid de… ▽ More The development of an extensive air shower depends not only on the nature of the primary ultra-high-energy cosmic ray but also on the properties of the hadronic interactions. For energies above those achievable in human-made accelerators, hadronic interactions are only accessible through the studies of extensive air showers, which can be measured at the Pierre Auger Observatory. With its hybrid detector design, the Pierre Auger Observatory measures both the longitudinal development of showers in the atmosphere and the lateral distribution of particles that arrive at the ground. This way, observables that are sensitive to hadronic interactions at ultra-high energies can be obtained. While the hadronic interaction cross-section can be assessed from the longitudinal profiles, the number of muons and their fluctuations measured with the ground detectors are linked to other physical properties. In addition to these direct studies, we discuss here how measurements of the atmospheric depth of the maximum of air-shower profiles and the characteristics of the muon signal at the ground can be used to test the self-consistency of the post-LHC hadronic models. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: Proceedings of VHEPU 56th Rencontres de Moriond 2022

A search for photons with energies above $2{\times}10^{17}$ eV using hybrid data from the low-energy extensions of the Pierre Auger Observatory

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

Abstract: Ultra-high-energy photons with energies exceeding $10^{17}$ eV offer a wealth of connections to different aspects of cosmic-ray astrophysics as well as to gamma-ray and neutrino astronomy. The recent observations of photons with energies in the $10^{15}$ eV range further motivate searches for even higher-energy photons. In this paper, we present a search for photons with energies exceeding… ▽ More Ultra-high-energy photons with energies exceeding $10^{17}$ eV offer a wealth of connections to different aspects of cosmic-ray astrophysics as well as to gamma-ray and neutrino astronomy. The recent observations of photons with energies in the $10^{15}$ eV range further motivate searches for even higher-energy photons. In this paper, we present a search for photons with energies exceeding $2{\times}10^{17}$ eV using about 5.5 years of hybrid data from the low-energy extensions of the Pierre Auger Observatory. The upper limits on the integral photon flux derived here are the most stringent ones to date in the energy region between $10^{17}$ and $10^{18}$ eV. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: Accepted for publication in ApJ

Report number: FERMILAB-PUB-22-441-AD-PPD-SCD-TD

Journal ref: ApJ 933 (2022) 125

No nematicity at the onset temperature of the pseudogap phase in the cuprate superconductor YBCO

Authors: G. Grissonnanche, O. Cyr-Choinière, J. Day, R. Liang, D. A. Bonn, W. N. Hardy, N. Doiron-Leyraud, L. Taillefer

Abstract: Electronic nematicity is the spontaneous loss of rotational symmetry in a metal, without breaking translational symmetry. In the cuprate superconductors, there is experimental evidence for nematicity, but its origin remains unclear. Here we investigate the onset of nematicity in the transport of charge by means of electric and thermoelectric measurements in underdoped YBa$_{\rm 2}$Cu$_{\rm 3}$O… ▽ More Electronic nematicity is the spontaneous loss of rotational symmetry in a metal, without breaking translational symmetry. In the cuprate superconductors, there is experimental evidence for nematicity, but its origin remains unclear. Here we investigate the onset of nematicity in the transport of charge by means of electric and thermoelectric measurements in underdoped YBa$_{\rm 2}$Cu$_{\rm 3}$O$_{\rm y}$, performed by passing the current (electrical or thermal) first along the $a$ axis then the $b$ axis of the orthorhombic structure in the same crystal, with a hole doping $p = 0.12$. Upon cooling, we observe no additional in-plane anisotropy -- beyond the background anisotropy due to the CuO chains -- in either the resistivity $ρ$ or the Seebeck coefficient $S$ as the temperature $T^{\star}$~for the onset of the pseudogap phase is crossed. We conclude that the pseudogap phase of cuprates is not nematic. However, at temperatures much lower than $T^{\star}$, a strong additional anisotropy is observed, most clearly in the Peltier coefficient $α= S / ρ$. We interpret it as nematicity associated with the development of charge order. △ Less

Submitted 14 February, 2023; v1 submitted 10 May, 2022; originally announced May 2022.

Comments: 7 pages, 4 figures

Journal ref: Physical Review X 13, 031010 (2023)

Formal Languages via Theories over Strings

Authors: Joel D. Day, Vijay Ganesh, Nathan Grewal, Florin Manea

Abstract: We investigate the properties of formal languages expressible in terms of formulas over quantifier-free theories of word equations, arithmetic over length constraints, and language membership predicates for the classes of regular, visibly pushdown, and deterministic context-free languages. In total, we consider 20 distinct theories and decidability questions for problems such as emptiness and univ… ▽ More We investigate the properties of formal languages expressible in terms of formulas over quantifier-free theories of word equations, arithmetic over length constraints, and language membership predicates for the classes of regular, visibly pushdown, and deterministic context-free languages. In total, we consider 20 distinct theories and decidability questions for problems such as emptiness and universality for formal languages over them. First, we discuss their relative expressive power and observe a rough division into two hierarchies based on whether or not word equations are present. Second, we consider the decidability status of several important decision problems, such as emptiness and universality. Note that the emptiness problem is equivalent to the satisfiability problem over the corresponding theory. Third, we consider the problem of whether a language in one theory is expressible in another and show several negative results in which this problem is undecidable. These results are particularly relevant in the context of normal forms in both practical and theoretical aspects of string solving. △ Less

Submitted 1 May, 2022; originally announced May 2022.

Non-local electrodynamics in ultra-pure PdCoO$_{2}$

Authors: Graham Baker, Timothy W. Branch, J. S. Bobowski, James Day, Davide Valentinis, Mohamed Oudah, Philippa McGuinness, Seunghyun Khim, Piotr Surówka, Yoshiteru Maeno, Roderich Moessner, Jörg Schmalian, Andrew P. Mackenzie, D. A. Bonn

Abstract: The motion of electrons in the vast majority of conductors is diffusive, obeying Ohm's law. However, the recent discovery and growth of high-purity materials with extremely long electronic mean free paths has sparked interest in non-ohmic alternatives, including viscous and ballistic flow. Although non-ohmic transport regimes have been discovered across a range of materials, including two-dimensio… ▽ More The motion of electrons in the vast majority of conductors is diffusive, obeying Ohm's law. However, the recent discovery and growth of high-purity materials with extremely long electronic mean free paths has sparked interest in non-ohmic alternatives, including viscous and ballistic flow. Although non-ohmic transport regimes have been discovered across a range of materials, including two-dimensional electron gases, graphene, topological semimetals, and the delafossite metals, determining their nature has proved to be challenging. Here, we report on a new approach to the problem, employing broadband microwave spectroscopy of the delafossite metal PdCoO$_{2}$ in three distinct sample geometries that would be identical for diffusive transport. The observed differences, which go as far as differing power laws, take advantage of the hexagonal symmetry of PdCoO$_{2}$. This permits a particularly elegant symmetry-based diagnostic for non-local electrodynamics, with the result favouring ballistic over strictly hydrodynamic flow. Furthermore, it uncovers a new effect for ballistic electron flow, owing to the highly facetted shape of the hexagonal Fermi surface. We combine our extensive dataset with an analysis of the Boltzmann equation to characterize the non-local regime in PdCoO$_{2}$. More broadly, our results highlight the potential of broadband microwave spectroscopy to play a central role in investigating exotic transport regimes in the new generation of ultra-high conductivity materials. △ Less

Submitted 1 March, 2023; v1 submitted 29 April, 2022; originally announced April 2022.

Limits to gauge coupling in the dark sector set by the non-observation of instanton-induced decay of Super-Heavy Dark Matter in the Pierre Auger Observatory data

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

Abstract: Instantons, which are non-perturbative solutions to Yang-Mills equations, provide a signal for the occurrence of quantum tunneling between distinct classes of vacua. They can give rise to decays of particles otherwise forbidden. Using data collected at the Pierre Auger Observatory, we search for signatures of such instanton-induced processes that would be suggestive of super-heavy particles decayi… ▽ More Instantons, which are non-perturbative solutions to Yang-Mills equations, provide a signal for the occurrence of quantum tunneling between distinct classes of vacua. They can give rise to decays of particles otherwise forbidden. Using data collected at the Pierre Auger Observatory, we search for signatures of such instanton-induced processes that would be suggestive of super-heavy particles decaying in the Galactic halo. These particles could have been produced during the post-inflationary epoch and match the relic abundance of dark matter inferred today. The non-observation of the signatures searched for allows us to derive a bound on the reduced coupling constant of gauge interactions in the dark sector: $α_X \lesssim 0.09$, for $10^{9} \lesssim M_X/{\rm GeV} < 10^{19}$. Conversely, we obtain that, for instance, a reduced coupling constant $α_X = 0.09$ excludes masses $M_X \gtrsim 3\times 10^{13}~$GeV. In the context of dark matter production from gravitational interactions alone, we illustrate how these bounds are complementary to those obtained on the Hubble rate at the end of inflation from the non-observation of tensor modes in the cosmological microwave background. △ Less

Submitted 15 December, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 9 pages, 3 figures, accepted for publication in Phys. Rev. Lett

Journal ref: Phys.Rev.Lett. 130 (2023) 6, 061001

Low-Metallicity Globular Clusters in the Low-Mass Isolated Spiral Galaxy NGC 2403

Authors: Duncan A. Forbes, Anna Ferré-Mateu, Jonah S. Gannon, Aaron J. Romanowsky, Jeffrey L. Carlin, Jean P. Brodie, Jacob Day

Abstract: The globular cluster (GC) systems of low-mass late-type galaxies, such as NGC 2403, have been poorly studied to date. As a low mass galaxy (M$_{\ast}$ = 7 $\times$ 10$^{9}$ M$_{\odot}$), cosmological simulations predict NGC 2403 to contain few, if any, accreted GCs. It is also isolated, with a remarkably undisturbed HI disk. Based on candidates from the literature, Sloan Digital Sky Survey and Hyp… ▽ More The globular cluster (GC) systems of low-mass late-type galaxies, such as NGC 2403, have been poorly studied to date. As a low mass galaxy (M$_{\ast}$ = 7 $\times$ 10$^{9}$ M$_{\odot}$), cosmological simulations predict NGC 2403 to contain few, if any, accreted GCs. It is also isolated, with a remarkably undisturbed HI disk. Based on candidates from the literature, Sloan Digital Sky Survey and Hyper Suprime-Cam imaging, we selected several GCs for follow-up spectroscopy using the Keck Cosmic Web Imager. From their radial velocities, and other properties, we identify 8 bona-fide GCs associated with either the inner halo or the disk of this bulgeless galaxy. A stellar population analysis suggests a wide range of GC ages from shortly after the Big Bang until the present day. We find all of the old GCs to be metal-poor with [Fe/H] $\le$ --1. The age--metallicity relation for the observed GCs suggests that they were formed over many Gyr from gas with a low effective yield, similar to that observed in the SMC. Outflows of enriched material may have contributed to the low yield. With a total system of $\sim$50 GCs expected, our study is the first step in fully mapping the star cluster history of NGC 2403 in both space and time. △ Less

Submitted 21 February, 2022; originally announced February 2022.

Comments: 9 pages, 7 figures, accepted for publication in MNRAS

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

Authors: The ANTARES collaboration, A. Albert, S. Alves, M. André, M. Anghinolfi, M. Ardid, S. Ardid, J. -J. Aubert, J. Aublin, B. Baret, S. Basa, B. Belhorma, M. Bendahman, V. Bertin, S. Biagi, M. Bissinger, J. Boumaaza, M. Bouta, M. C. Bouwhuis, H. Brânzaş, R. Bruijn, J. Brunner, J. Busto, B. Caiffi, D. Calvo , et al. (1025 additional authors not shown)

Abstract: For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for corre… ▽ More For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs. △ Less

Submitted 23 August, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 39 pages, 7 figures, 4 tables; updated source files including xml authorlist

Report number: FERMILAB-PUB-22-033-AD-PPD-SCD-TD

Journal ref: ApJ 934 164 (2022)

Testing effects of Lorentz invariance violation in the propagation of astroparticles with the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, K. Almeida Cheminant, 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, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, J. A. Bellido , et al. (352 additional authors not shown)

Abstract: Lorentz invariance violation (LIV) is often described by dispersion relations of the form $E_i^2=m_i^2+p_i^2+δ_{i,n} E^{2+n}$ with delta different based on particle type $i$, with energy $E$, momentum $p$ and rest mass $m$. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constrai… ▽ More Lorentz invariance violation (LIV) is often described by dispersion relations of the form $E_i^2=m_i^2+p_i^2+δ_{i,n} E^{2+n}$ with delta different based on particle type $i$, with energy $E$, momentum $p$ and rest mass $m$. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constraints on the LIV coefficients $δ_{i,n}$ tend to come from the highest energies. At sufficiently high energies, photons produced by cosmic ray interactions as they propagate through the Universe could be subluminal and unattenuated over cosmological distances. Cosmic ray interactions can also be modified and lead to detectable fingerprints in the energy spectrum and mass composition observed on Earth. The data collected at the Pierre Auger Observatory are therefore possibly sensitive to both the electromagnetic and hadronic sectors of LIV. In this article, we explore these two sectors by comparing the energy spectrum and the composition of cosmic rays and the upper limits on the photon flux from the Pierre Auger Observatory with simulations including LIV. Constraints on LIV parameters depend strongly on the mass composition of cosmic rays at the highest energies. For the electromagnetic sector, while no constraints can be obtained in the absence of protons beyond $10^{19}$ eV, we obtain $δ_{γ,0} > -10^{-21}$, $δ_{γ,1} > -10^{-40}$ eV$^{-1}$ and $δ_{γ,2} > -10^{-58}$ eV$^{-2}$ in the case of a subdominant proton component up to $10^{20}$ eV. For the hadronic sector, we study the best description of the data as a function of LIV coefficients and we derive constraints in the hadronic sector such as $δ_{\mathrm{had},0} < 10^{-19}$, $δ_{\mathrm{had},1} < 10^{-38}$ eV$^{-1}$ and $δ_{\mathrm{had},2}< 10^{-57}$ eV$^{-2}$ at 5$σ$ CL. △ Less

Submitted 19 January, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: Comments: Published version, 22 pages, 8 figures, 2 tables

Report number: FERMILAB-PUB-21-716-AD-PPD-SCD-TD

Journal ref: JCAP 01 (2022) 023

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

String Theories involving Regular Membership Predicates: From Practice to Theory and Back

Authors: Murphy Berzish, Joel D. Day, Vijay Ganesh, Mitja Kulczynski, Florin Manea, Federico Mora, Dirk Nowotka

Abstract: Widespread use of string solvers in formal analysis of string-heavy programs has led to a growing demand for more efficient and reliable techniques which can be applied in this context, especially for real-world cases. Designing an algorithm for the (generally undecidable) satisfiability problem for systems of string constraints requires a thorough understanding of the structure of constraints pre… ▽ More Widespread use of string solvers in formal analysis of string-heavy programs has led to a growing demand for more efficient and reliable techniques which can be applied in this context, especially for real-world cases. Designing an algorithm for the (generally undecidable) satisfiability problem for systems of string constraints requires a thorough understanding of the structure of constraints present in the targeted cases. In this paper, we investigate benchmarks presented in the literature containing regular expression membership predicates, extract different first order logic theories, and prove their decidability, resp. undecidability. Notably, the most common theories in real-world benchmarks are PSPACE-complete and directly lead to the implementation of a more efficient algorithm to solving string constraints. △ Less

Submitted 15 May, 2021; originally announced May 2021.

Comments: arXiv admin note: substantial text overlap with arXiv:2010.07253

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)

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)

Calibration of the underground muon detector of the Pierre Auger Observatory

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

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

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

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

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

Journal ref: JINST (2021) 16 P04003

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

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

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

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

Comments: 40 pages, 33 figures

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

Journal ref: JINST 16 (2021) P01026

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

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

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

Submitted 21 October, 2020; originally announced October 2020.

Comments: Published version, 4 figures

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

Evidence for transient atmospheres during eruptive outgassing on the Moon

Authors: Elishevah van Kooten, Frédéric Moynier, James Day

Abstract: Events following the giant impact formation of the Moon are thought to have led to volatile depletion and concurrent mass-dependent fractionation of the isotopes of moderately volatile elements (MVE). The detailed processes and conditions surrounding this episode remain obscured and are not unified by a single model for all volatile elements and compounds. Using available data, including new Zn is… ▽ More Events following the giant impact formation of the Moon are thought to have led to volatile depletion and concurrent mass-dependent fractionation of the isotopes of moderately volatile elements (MVE). The detailed processes and conditions surrounding this episode remain obscured and are not unified by a single model for all volatile elements and compounds. Using available data, including new Zn isotope data for eight lunar samples, we demonstrate that the isotopic fractionation of MVE in the Moon is best expressed by non-ideal Rayleigh distillation, approaching the fractionation factor α using the reduced masses of the evaporated isotopologues. With these calculations a best fit for the data is obtained when the lunar MVE isotope data is normalized to ordinary or enstatite chondrites (ΔMoon-OC,EC), rather than a bulk silicate Earth (BSE) composition. This analysis further indicates that the parent body from which the Moon formed cannot have partitioned S into its core based on S isotope compositions of lunar rocks. The best fit between ΔMoon-OC,EC and modelled non-ideal Rayleigh fractionation is defined by a slope that corresponds to a saturation index of 90 +/- 4 %. In contrast, the older Highland suite is defined by a saturation index of 75 +/- 2 %, suggesting the vapor phase pressure was higher during mare basalt eruptions. This provides the first tangible evidence that the Moon was veiled by a thin atmosphere during mare basalt eruption events spanning at least from 3.8 to 3 billion years ago and implies that MVE isotope fractionation dominantly occurred after the Moon had accreted. △ Less

Submitted 15 October, 2020; originally announced October 2020.

Comments: Accepted, preprint

Journal ref: Planetary Science Journal, 2020

An SMT Solver for Regular Expressions and Linear Arithmetic over String Length

Authors: Murphy Berzish, Mitja Kulczynski, Federico Mora, Florin Manea, Joel D. Day, Dirk Nowotka, Vijay Ganesh

Abstract: We present a novel length-aware solving algorithm for the quantifier-free first-order theory over regex membership predicate and linear arithmetic over string length. We implement and evaluate this algorithm and related heuristics in the Z3 theorem prover. A crucial insight that underpins our algorithm is that real-world instances contain a wealth of information about upper and lower bounds on len… ▽ More We present a novel length-aware solving algorithm for the quantifier-free first-order theory over regex membership predicate and linear arithmetic over string length. We implement and evaluate this algorithm and related heuristics in the Z3 theorem prover. A crucial insight that underpins our algorithm is that real-world instances contain a wealth of information about upper and lower bounds on lengths of strings under constraints, and such information can be used very effectively to simplify operations on automata representing regular expressions. Additionally, we present a number of novel general heuristics, such as the prefix/suffix method, that can be used in conjunction with a variety of regex solving algorithms, making them more efficient. We showcase the power of our algorithm and heuristics via an extensive empirical evaluation over a large and diverse benchmark of 57256 regex-heavy instances, almost 75% of which are derived from industrial applications or contributed by other solver developers. Our solver outperforms five other state-of-the-art string solvers, namely, CVC4, OSTRICH, Z3seq, Z3str3, and Z3-Trau, over this benchmark, in particular achieving a 2.4x speedup over CVC4, 4.4x speedup over Z3seq, 6.4x speedup over Z3-Trau, 9.1x speedup over Z3str3, and 13x speedup over OSTRICH. △ Less

Submitted 7 May, 2021; v1 submitted 14 October, 2020; originally announced October 2020.

Comments: 25 pages (main body 21 pages). 7 figures, 6 tables

Advancing the Scientific Frontier with Increasingly Autonomous Systems

Authors: Rashied Amini, Abigail Azari, Shyam Bhaskaran, Patricia Beauchamp, Julie Castillo-Rogez, Rebecca Castano, Seung Chung, John Day, Richard Doyle, Martin Feather, Lorraine Fesq, Jeremy Frank, P. Michael Furlong, Michel Ingham, Brian Kennedy, Ksenia Kolcio, Issa Nesnas, Robert Rasmussen, Glenn Reeves, Cristina Sorice, Bethany Theiling, Jay Wyatt

Abstract: A close partnership between people and partially autonomous machines has enabled decades of space exploration. But to further expand our horizons, our systems must become more capable. Increasing the nature and degree of autonomy - allowing our systems to make and act on their own decisions as directed by mission teams - enables new science capabilities and enhances science return. The 2011 Planet… ▽ More A close partnership between people and partially autonomous machines has enabled decades of space exploration. But to further expand our horizons, our systems must become more capable. Increasing the nature and degree of autonomy - allowing our systems to make and act on their own decisions as directed by mission teams - enables new science capabilities and enhances science return. The 2011 Planetary Science Decadal Survey (PSDS) and on-going pre-Decadal mission studies have identified increased autonomy as a core technology required for future missions. However, even as scientific discovery has necessitated the development of autonomous systems and past flight demonstrations have been successful, institutional barriers have limited its maturation and infusion on existing planetary missions. Consequently, the authors and endorsers of this paper recommend that new programmatic pathways be developed to infuse autonomy, infrastructure for support autonomous systems be invested in, new practices be adopted, and the cost-saving value of autonomy for operations be studied. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Comments: 10 pages (compared to 8 submitted to PSADS), 2 figures, submitted to National Academy of Sciences Planetary Science and Astrobiology Decadal Survey 2023-2032

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

Calculating Great Britains half-hourly electrical demand from publicly available data

Authors: IA Grant Wilson, Shivangi Sharma, Joseph Day, Noah Godfrey

Abstract: Here we present a method to combine half-hourly publicly available electrical generation and interconnector operational data for Great Britain to create a timeseries that approximates its electrical demand. We term the calculated electrical demand ESPENI that is an acronym for Elexon Sum Plus Embedded Net Imports. The method adds value to the original data by combining both transmission and distri… ▽ More Here we present a method to combine half-hourly publicly available electrical generation and interconnector operational data for Great Britain to create a timeseries that approximates its electrical demand. We term the calculated electrical demand ESPENI that is an acronym for Elexon Sum Plus Embedded Net Imports. The method adds value to the original data by combining both transmission and distribution generation data into a single dataset and adding ISO 8601 compatible datetimes to increase interoperability with other timeseries data. Data cleansing is undertaken by visually flagging data errors and then using simple linear interpolation to impute values to replace the flagged data. Publishing the method allows it to be further enhanced or adapted and to be considered and critiqued by a wider community. In addition, the published raw and cleaned data is a valuable resource that saves researchers considerable time in repeating the steps presented in the method to prepare the data for further analysis. The data is a public record of the decarbonisation of Great Britains electrical system since late 2008, widely seen as an example of rapid decarbonisation of an electrical system away from fossil fuel generation to lower carbon sources. △ Less

Submitted 15 September, 2021; v1 submitted 28 June, 2020; originally announced June 2020.

Comments: 33 pages, 3 Figures, 6 tables

On the microscopic origin of reversible and irreversible reactions of LiNixCoyMnx cathode materials: Ni-O hybrid bond formation vs. cationic and anionic redox

Authors: Karin Kleiner, Claire A. Murray, Cristina Grosu, Sarah J. Day, Martin Winter, Peter Nagel, Stefan Schuppler, Michael Merz

Abstract: Energy density limitations of layered oxides with different Ni contents, i.e., of the conventional cathode materials in Li-ion batteries, are investigated across the first discharge cycle using advanced spectroscopy and state-of-the-art diffraction. For the first time unambiguous experimental evidence is provided, that redox reactions in NCMs proceed via a reversible oxidation of Ni and a hybridiz… ▽ More Energy density limitations of layered oxides with different Ni contents, i.e., of the conventional cathode materials in Li-ion batteries, are investigated across the first discharge cycle using advanced spectroscopy and state-of-the-art diffraction. For the first time unambiguous experimental evidence is provided, that redox reactions in NCMs proceed via a reversible oxidation of Ni and a hybridization with O, and not, as widely assumed, via pure cationic or more recently discussed, pure anionic redox processes. Once Ni-O hybrid states are formed, the sites cannot be further oxidized. Instead, irreversible reactions set in which lead to a structural collapse and thus, the lack of ionic Ni limits the reversible capacity. Moreover, the degree of hybridization, which varies with the Ni content, triggers the electronic structure and the operation potential of the cathodes. With an increasing amount of Ni, the covalent character of the materials increases and the potential decreases. △ Less

Submitted 17 June, 2020; originally announced June 2020.

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

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

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

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

Comments: Published version; correction of a minor typo

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

Journal ref: Published in JCAP06(2020)017

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

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

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

Submitted 14 February, 2020; originally announced February 2020.

Comments: Accepted by ApJ

Report number: FERMILAB-PUB-20-076

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

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

On Solving Word Equations Using SAT

Authors: Joel D. Day, Thorsten Ehlers, Mitja Kulczynski, Florin Manea, Dirk Nowotka, Danny Bøgsted Poulsen

Abstract: We present Woorpje, a string solver for bounded word equations (i.e., equations where the length of each variable is upper bounded by a given integer). Our algorithm works by reformulating the satisfiability of bounded word equations as a reachability problem for nondeterministic finite automata, and then carefully encoding this as a propositional satisfiability problem, which we then solve using… ▽ More We present Woorpje, a string solver for bounded word equations (i.e., equations where the length of each variable is upper bounded by a given integer). Our algorithm works by reformulating the satisfiability of bounded word equations as a reachability problem for nondeterministic finite automata, and then carefully encoding this as a propositional satisfiability problem, which we then solve using the well-known Glucose SAT-solver. This approach has the advantage of allowing for the natural inclusion of additional linear length constraints. Our solver obtains reliable and competitive results and, remarkably, discovered several cases where state-of-the-art solvers exhibit a faulty behaviour. △ Less

Submitted 27 June, 2019; originally announced June 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

Mapping the structural transitions controlled by the trilinear coupling in Ca3-xSrxTi2O7

Authors: Marie Kratochvilova, Fei-Ting Huang, Maria-Teresa Fernandez Diaz, Milan Klicpera, Sarah J. Day, Stephen P. Thompson, Yoon-Seok Oh, Bin Gao, Sang-Wook Cheong, Je-Geun Park

Abstract: We present the results of the high-temperature neutron and x-ray diffraction experiments on the Ca3-xSrxTi2O7 (x = 0.5, 0.8, 0.85, 0.9) compounds. The ferro- to paraelectric transition in these hybrid improper ferroelectric materials arises from the so-called trilinear coupling. Depending on the Strontium content, various structures and phase transitions, different from theoretical predictions, em… ▽ More We present the results of the high-temperature neutron and x-ray diffraction experiments on the Ca3-xSrxTi2O7 (x = 0.5, 0.8, 0.85, 0.9) compounds. The ferro- to paraelectric transition in these hybrid improper ferroelectric materials arises from the so-called trilinear coupling. Depending on the Strontium content, various structures and phase transitions, different from theoretical predictions, emerge. The in-situ x-ray powder diffraction indicates a direct ferro- to paraelectric transition between the orthorhombic A21am and tetragonal undistorted I4/mmm phase for x < 0.6. We identified a reduction in the trilinear coupling robustness by increasing the Sr-doping level to lead to the emergence of the intermediate tetragonal P42/mnm phase and the gradual suppression of the orthorhombic phase. The observed character of the structure transitions and the Ca3-xSrxTi2O7 phase diagram are discussed in the framework of theoretical models of other related hybrid improper ferroelectric systems. △ Less

Submitted 4 June, 2019; originally announced June 2019.

Comments: 10 pages, 3 figures, 1 table

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