Ultra-high Energy Cosmic Rays Research Articles

Analysing arrival directions of ultra-high-energy cosmic rays with convolutional neural networks

The problem of identification of ultra-high-energy cosmic ray (UHECR) sources is greatly complicated by the fact that even the highest energy cosmic rays may be deflected by tens of degrees in the galactic magnetic fields. We show that arrival directions of UHECRs from several nearest active galaxies form specific patterns in the sky, which can be effectively recognized by convolutional neural networks. We use one of the recently developed convnet implementations for images defined on the sphere to train the classifier that is able to detect patterns that can be present in the experimental data. We calculate the minimal detectable from-source event fractions for several realistic source candidates and discuss the method limitations.

A Simulation Study of Ultra-relativistic Jets. III. Particle Acceleration in FR-II Jets

We study the acceleration of ultra-high-energy cosmic rays (UHECRs) in Class II Fanaroff–Riley (FR-II) radio galaxies by performing Monte Carlo simulations for the transport, scattering, and energy change of CR particles injected into time-evolving jet flows that are realized through relativistic hydrodynamic simulations. Toward that end, we adopt physically motivated models for the magnetic field and particle scattering. By identifying the primary acceleration process among diffusive shock acceleration (DSA), turbulent shear acceleration (TSA), and relativistic shear acceleration (RSA), we find that CRs of E ≲ 1 EeV gain energy mainly through DSA in the jet‐spine flow and backflow containing many shocks and turbulence. After they attain E ≳ a few exaelectronvolts, CRs are energized mostly via RSA at the jet–backflow interface, reaching energies well above 1020 eV. TSA makes a relatively minor contribution. The time-asymptotic energy spectrum of escaping particles is primarily governed by the jet power, shifting to higher energies at more powerful jets. The UHECR spectrum fits well to a double power-law form, whose break energy, E break, corresponds to the size-limited maximum energy. It is close to below E break, while it follows above E break, decreasing more gradually than the exponential. The power-law slope of the high-energy end is determined by energy boosts via non-gradual shear acceleration across the jet–backflow interface and confinement by an elongated cocoon. We conclude that giant radio galaxies could be major contributors to the observed UHECRs.

Detection of ultra high energy cosmic rays and neutrinos with lunar orbital radio telescope

Particle cascades induced by ultra-high-energy (UHE) cosmic rays and neutrinos impacting on the lunar regolith usually radiate Cherenkov radio emissions due to the presence of excess negative charge, which is known as Askaryan effect. Several experiments have been carried out to detect the Cherenkov radio emissions in the lunar regolith. To prepare for future lunar Ultra-Long Wavelength (ULW, frequencies below 30 MHz) radio astronomy missions, we study the detection of the Cherenkov radio emissions with the ULW radio telescope that are operating at the lunar orbit. We have carried out instrument modelling and analytic calculations for the analysis of aperture, flux and event rate, and the analyses show the detectability of the Cherenkov radiation. Based on the properties of the Cherenkov radiation, we have demonstrated that the cosmic ray and neutrino events could be reconstructed with the three ULW vector antennas onboard the lunar satellites via measurements of the Askaryan radio pulse intensity, polarizations, etc. The results obtained by this study would be useful for future lunar radio explorer mission, where the detections of UHE cosmic rays and neutrinos could be successfully attempted.

Ultra High Energy Cosmic Rays an overview

We review the main experimental evidences on ultra high energy cosmic rays and their implications in the physics of these extremely energetic particles, also in connection with dark matter and cosmology. We discuss the basis of theoretical models aiming at explaining observations, highlighting the most relevant open questions in this fascinating field of research.

The science of ultra-high energy cosmic rays after more than 15 years of operation of the Pierre Auger Observatory

The Pierre Auger Observatory has been detecting ultra-high energy cosmic rays (UHECRs) for more than fifteen years. An essential feature of the Observatory is its hybrid design: cosmic rays above 100 PeV are detected through the observation of the associated air showers with different and complementary techniques, from surface detector arrays and fluorescence telescopes to radio antennas. The analyses of the multi-detector data have enabled high-statistics and high-precision studies of the energy spectrum, mass composition and distribution of arrival directions of UHECRs. The resulting picture is summarized in this contribution. While no discrete source of UHECRs has been identified so far, the extragalactic origin of the particles has been recently determined from the arrival directions above 8 EeV, and the ring is closing around nearby astrophysical sites at higher energies. Also, the established upper limits on fluxes of UHE neutrinos and photons have implications on dark matter and cosmological aspects that are also presented in this contribution.

Observational constraints on accelerators of ultra-high energy cosmic rays

We explore two generic hypotheses for tracing the sources of ultra-high energy cosmic rays (UHECRs) in the Universe: star formation rate density or stellar mass density. For each scenario, we infer a set of constraints for the emission mechanisms in the accelerators, for their energetics and for the abundances of elements at escape from their environments. From these constraints, we generate sky maps above 40 EeV expected from a catalog that comprises 410,761 galaxies out to 350 Mpc and provides a near-infrared flux-limited sample to map both stellar mass and star formation rate over the full sky. Considering a scenario of intermittent sources hosted in every galaxy, we show that the main features observed in arrival directions of UHECRs can in turn constrain the burst rate of the sources provided that magnetic-horizon effects are at play in clusters of galaxies.

Selected results from IceCube

Neutrino astronomy saw its birth with the discovery by IceCube of a difFuse flux at energies above 60 TeV with intensity comparable to a predicted upper limit to the flux from extra-galactic sources of ultra-high energy cosmic rays (UHECRs). While such an upper limit corresponds to the case of calorimetric sources, in which cosmic rays lose all their energy into photo-pion production, the first statistically significant coincident observation between neutrinos and gamma-rays was observed from a blazar of intriguing nature. A very-high-energy muon event, of most probable neutrino energy of 290 TeV for an E −2.13 spectrum, alerted other observatories triggering a large amount of investigations in many bands of the electromagnetic (EM) spectrum. A high gamma-ray state from the blazar was revealed soon after the event and in a follow up to about 40 days. A posteriori observations also in the optical and in the radio indicated a rise of the flux from the TXS 0506+056 blazar. A previous excess of events of duration of more than 100 d was observed by IceCube with higher significance than the alert itself. These observations triggered more complex modelling than simple one zone proton synchrotron models for proton acceleration in jets of active galactic nuclei (AGNs) and more observations across the EM spectrum.A second evidence was a steady excess of about 50 neutrino events with reconstructed soft spectrum in a sample of lower energy well reconstructed muon events than the alert event. A hot spot was identified in a catalogue of 110 gamma-ray intense emitters and starburst galaxies in a direction compatible to NGC 1068 with significance of 2.9σ. NGC 1068 hosts a mildly relativistic jet in a starburst galaxy, seen not from the jet direction but rather through the torus. This Seyfert II galaxy is at only 14.4 Mpc from the Earth. The source turned out to be also the hottest spot of an all-sky search. Analysed cumulatively, the catalogue excess was 3.3σ with the contribution of NGC 1068 and TXS 0506+056, as expected, and other 2 sources, PKS 1424+240, and GB6 J1542+6129, with similar features to TXS 0506+056, indicating that they might all be Flat Spectrum Radio Quasars (FSRQs).While all these observations and the directions of the measured events contributing to diffuse fluxes hint to their extra-galactic origin, a few percent level contribution might be the end of a lower energy ‘granted’ flux of neutrinos from interactions of cosmic rays in the Galactic Plane. This relevant observation is at the reach of IceCube and other neutrino telescopes. These aspects were discussed at the conference and are summarised in this write up.

High-energy neutrino astronomy — the neutrino connections to the cosmic-ray origin: present and future

High-energy neutrino astronomy has been blooming. In addition to the possible identification of the blazar and Seyfert II galaxies as neutrino emitters, the present neutrino data has indicated some hints to characterize or constrain the origin of cosmic rays. Being motivated by the observational fact that the astrophysical neutrino background energy flux is comparable to that of ultrahigh-energy (UHE) cosmic rays, we derive the generic requirements that a major fraction of UHE cosmic ray sources must meet, if they are also responsible for the ∼ 100 TeV-energy cosmic neutrino background radiation. The source parameters characterizing the cosmic ray – neutrino unified scheme, such as the photon radiation luminosity and the cosmic ray luminosity density, suggest that the yet-unidentified cosmic ray and neutrino origins can be transient objects visible in the optical/NIR wavelength band. We propose a viable scheme of multimessenger observations to identify the sources using the neutrino multiplet detections.

Multi-messenger implications of the Pierre Auger Observatory measurements

The Pierre Auger Observatory is not only the world’s largest cosmic-ray observatory, but also an efficient detector of cosmic particles other than nuclei. Using its enormous exposure, the Observatory is a competitive player in searches for ultra-high energy neutrinos, photons, and even particles that could be predicted by theories beyond the standard model of particle physics. We present our recent results of such searches, targeting both the diffuse fluxes of particles and specific transient events. Special attention is given to the potential of the Observatory to support or constrain the observations of the ANITA experiment interpreted as upward-going air showers. The search for correlations of arrival directions of ultra-high energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array with high-energy neutrinos detected by IceCube and ANTARES experiments is also discussed.

The Astroparticle Detectors Array—An Educational Project in Cosmic Ray Physics

ADA, short for Astroparticle Detectors Array, is an educational project aiming to detect cosmic radiation and possibly high-energy particles known as ultra-high-energy cosmic rays (UHECRs) or even to spot a supernova event. Its working process is the same as that used in professional cosmic ray observatories: it consists of simple detectors spread over the entire Italian territory and beyond. The detectors are hosted among high schools, associations, and private astronomical observatories. ADA has been operating since 2013 and was brought about with the intention of promoting astroparticle physics to any given level of outreach. Furthermore, ADA is becoming an interesting tool not only for teachers but also for independent and keen scientists. Over the years, we have discovered the importance of having a long series of data for studying the relation between and among cosmic rays, weather, and space weather and to investigate the main cause of oscillations in cosmic ray data. In this paper, we show what we find to be the most compelling results, such as the beautiful symmetry of the behavior between muons and the atmospheric temperature and, likewise, the evident anti-correlation between the intensity of the muons at ground level compared with solar activity.

UHECR Signatures and Sources

We discuss recent results on the clustering, composition and distribution of Ultra-High Energy Cosmic Rays (UHECR) in the sky; from the energy of several tens of EeV in the dipole anisotropy, up to the highest energy of a few narrow clusters, those of Hot Spots. Following the early UHECR composition records deviations from proton we noted that the UHECR events above 40 EeV can be made not just by any light or heavy nuclei, but mainly by the lightest ones as He,D, Li,Be. The remarkable Virgo absence and the few localized nearby extragalactic sources, such as CenA, NGC 253 and M82, are naturally understood: lightest UHECR nuclei cannot reach us from the Virgo distance of twenty Mpc, due to their nuclei fragility above a few Mpc distances. Their deflection and smearing in wide hot spots is better tuned to the lighter nuclei than to the preferred proton or heavy nuclei candidate courier. We note that these lightest nuclei still suffer of a partial photodistruction even from such close sources. Therefore, their distruption in fragments, within few tens EeV multiplet chain of events, have been expected and later on observed by Auger collaboration, nearly a decade ago. These multiplet presences, strongly correlate with the same CenA, NGC253 sources. The statistical weight of such correlation is reminded. We conclude that the same role of NGC 253 clustering at lower energies could also feed the Auger dipole anisotropy at lower energy ranges. Such lower energy anisotropy could be fed and integrated by nearest Vela, Crab, LMC and Cas A contributes. In our present UHECR model, based on lightest nuclei in local volumes of a few Mpcs, closest AGN, Star-Burst or very close SNR are superimposing their signals, frozen in different epochs, distances and directions, feeding small and wide anisotropy. Possible tests to confirm, or untangle the current model from alternative ones, are suggested and updated.

CRPropa high statistics simulations for UHECR anisotropy studies

Despite the great progress made by modern cosmic ray observatories, the origin and acceleration mechanism of ultra-high-energy cosmic rays (UHECRs) remains unsolved problems to this day. However, there is experimental evidence for an anisotropic component in the UHECR arrival direction at energies greater than few EeV. The search for UHECR sources is further complicated by two main factors: during extragalactic propagation UHECRs interact with background photon fields (like the CMB and the EBL) and, since they are electrically charged, they are deflected by extragalactic and galactic magnetic fields (EGMF and GMF). Moreover, the strength, structure and origin of the EGMF are still not well known, causing reconstructing the deflection of UHECRs to be a non-trivial task. In this work we consider several EGMF models obtained from constrained MagnetoHydroDynamics (MHD) simulations of our local Universe to study the propagation of UHECRs through such a structured environment. We simulate propagation, interactions and observation of UHECRs by using the Monte Carlo code CRPropa3. We also take into account the effect of the GMF by adopting a lensing procedure of the arrival UHECR sky map. We explore several combination of source distribution, EGMF structure and mass composition. As a reference, we also simulate scenarios without the EGMF and with a statistically homogeneous field. We use our simulation results to compute UHECR observables, such as the energy spectrum, the angular power spectrum and the arrival direction map (before and after the GMF) in order to obtain constraints on possible combinations of source distributions and EGMF models.

Examination of Xmax anisotropy for the next generation Ultra-high energy cosmic rays observations

The estimation of mass composition ultra-high energy cosmic rays is essential to understand their origin and generation mechanism. Recent experiments are expected to discover anisotropy of mass composition in ultra-high energy cosmic rays. Anisotropy of mass composition using Xmax is currently being performed, but the problem is that the statistics of cosmic rays with measured Xmax is very small. For the next generation experiments, it is important to know how much statistics can be accumulated to find Xmax anisotropy. Therefore, in this analysis, we examine Xmax anisotropic search in ultra-high energy cosmic rays by simulating a cosmic ray air showers under various conditions such as energy spectrum and mass composition.

The Snowmass UHECR White Paper on Ultra-High-Energy Cosmic Rays

This proceeding summarizes the talk given at the opening of the UHECR 2022 conference in L’Aquila on the whitepaper ‘Ultra-High-Energy Cosmic Rays: The Intersection of the Cosmic and Energy Frontiers’ [Astroparticle Physics 149 (2023) 102819 - arXiv:2205.05845] that has been prepared for the Snowmass survey in the USA. The whitepaper provides an overview of recent progress and open questions regarding the particle physics and astrophysics related to ultra-high-energy cosmic rays (UHECR) and outlines the connections between the particle and astrophysics aspects of cosmic rays. It also discusses what instrumentation is needed to address the major scientific questions in ultra-high-energy cosmic-ray physics. While the upgraded Pierre Auger Observatory and Telescope Array will remain the workhorses at the highest energies in the current decade, new experiments with significantly higher exposure are needed in the coming decade. Ground arrays featuring simultaneous detection of the position of the shower maximum and the size of the muonic component will enable particle astronomy by measuring the rigidity of individual events. They should be complemented by other detectors maximizing the total exposure. This can be achieved by a few next-generation experiments using the latest developments in detection and analysis techniques: GRAND as a ground-based radio array, and POEMMA as a space-borne stereo fluorescence telescope will feature complementary approaches to provide maximum exposure; IceCube-Gen2 with its surface array, and GCOS aim at increased statistics with high accuracy for particle physics and rigidity-based galactic and extra-galactic astrophysics. While designed to discover the astrophysical cosmic-ray sources at the highest energies, the same experiments also contribute to particle physics, e.g., by studying the muon puzzle in cosmic-ray air showers, and by their discovery potential for exciting new physics, such as certain Dark Matter candidates. With the full whitepaper available as a reference, this proceeding will briefly present the science cases of the experiments, highlighting their individual strengths and outlining how they complement each other.

Ultra High Energy Cosmic Ray Source Models: Successes, Challenges and General Predictions

Understanding the acceleration of Ultra High Energy Cosmic Rays is one of the great challenges of contemporary astrophysics. In this short review, we summarize the general observational constraints on their composition, spectrum and isotropy which indicate that nuclei heavier than single protons dominate their spectra above ˜ 5 EeV, that they are strongly suppressed above energies ˜ 50 EeV, and that the only significant departure from isotropy is a dipole. Constraints based upon photopion and photodisintegration losses allow their ranges and luminosity density to be estimated. Three general classes of source model are discussed - magnetospheric models (including neutron stars and black holes), jet models (including Gamma Ray Bursts, Active Galactic Nuclei and Tidal Disruption Events) and Diffusive Shock Acceleration models (involving large accretion shocks around rich clusters of galaxies). The value of constructing larger and more capable arrays to measure individual masses at the highest energies and probably identifying their sources is emphasized.

2022 report from the Auger-TA working group on UHECR arrival directions

After over 60 years, the powerful engines that accelerate ultra-high-energy cosmic rays (UHECRs) to the formidable energies at which we observe them from Earth remain mysterious. Assuming standard physics, we expect UHECR sources to lie within the local Universe (up to a few hundred Mpc). The distribution of matter in the local Universe is anisotropic, and we expect this anisotropy to be imprinted on the distribution of UHECR arrival directions. Even though intervening intergalactic and Galactic magnetic fields deflect charged UHECRs and can distort these anisotropies, some amount of information on the distribution of the sources is preserved. In this proceedings contribution, we present the results of the joint Pierre Auger Observatory and Telescope Array searches for (a) the largest-scale anisotropies (the harmonic dipole and quadrupole) and (b) correlations with a sample of nearby starburst galaxies and the 2MRS catalogue tracing stellar mass within 250 Mpc. This analysis updates our previous results with the most recent available data, notably with the addition of 3 years of new Telescope Array data. The main finding is a correlation between the arrival directions of 12.1%−3.1%+4.5% of UHECRs detected with E ≥ 38 EeV by Auger or with E ≳ 49 EeV by TA and the positions of nearby starburst galaxies on a 15.1°−3.0°+4.6° angular scale, with a 4.7σ post-trial significance, up from 4.2σ obtained in our previous study.

Testing hadronic and photo-hadronic interactions as responsible for UHECR and neutrino fluxes from Starburst Galaxies

We test the hypothesis that starburst galaxies are the sources of ultra-high energy cosmic rays and high-energy neutrinos. The computation of interactions of ultra-high energy cosmic rays in the starburst environment as well as in the propagation to the Earth is made using a modified version of the Monte Carlo code SimProp, where hadronic processes are implemented for the first time. Taking into account a star-formationrate distribution of sources, the fluxes of ultra-high energy cosmic rays and high-energy neutrinos are computed and compared with observations, and the explored parameter space for the source characteristics is discussed. We find that, depending on the density of the gas in the source environment, spallation reactions could hide the outcome in neutrinos from photo-hadronic interactions in the source environment and in extra-galactic space. We confirm that source-propagation models constitute a promising way to improve the discrimination power of models considering only ultra-high energy cosmic rays, on the way to unveiling the source class responsible for ultra-high energy cosmic rays and high-energy neutrinos.

High-energy Neutrino Emission from Espresso-reaccelerated Ions in Jets of Active Galactic Nuclei

We present a bottom-up calculation of the flux of ultrahigh-energy cosmic rays (UHECRs) and high-energy neutrinos produced by powerful jets of active galactic nuclei (AGNs). By propagating test particles in 3D relativistic magnetohydrodynamic jet simulations, including a Monte Carlo treatment of sub-grid pitch-angle scattering and attenuation losses due to realistic photon fields, we study the spectrum and composition of the accelerated UHECRs and estimate the amount of neutrinos produced in such sources. We find that UHECRs may not be significantly affected by photodisintegration in AGN jets where the espresso mechanism efficiently accelerates particles, consistent with Auger’s results that favor a heavy composition at the highest energies. Moreover, we present estimates and upper bounds for the flux of high-energy neutrinos expected from AGN jets. In particular, we find that (i) source neutrinos may account for a sizable fraction, or even dominate, the expected flux of cosmogenic neutrinos; (ii) neutrinos from the β-decay of secondary neutrons produced in nucleus photodisintegration end up in the teraelectronvolt to petaelectronvolt band observed by IceCube, but can hardly account for the observed flux; (iii) UHECRs accelerated via the espresso mechanism lead to nearly isotropic neutrino emission, which suggests that nearby radio galaxies may be more promising as potential sources. We discuss our results in light of multimessenger astronomy and current/future neutrino experiments.

Estimation and reduction of the biases by the galactic magnetic field on the UHECR correlation studies

The correlation studies between ultrahigh-energy cosmic ray (UHECR) anisotropy and source candidates are keys to understand the origin of UHECRs. Especially, studies with consideration of source models, magnetic fields, and mass compositions are necessary. We estimated the biases caused by the coherent deflection due to the galactic magnetic field (GMF) in the previous maximum-likelihood analysis for searching the UHECR sources (Aab et al. 2018, Abbasi et al. 2018). In our work with simulated mock datasets, we find that the anisotropic fraction fani is estimated systematically lower than the true value when we ignore the effect caused by the GMF. We also develop the maximum-likelihood method which includes the GMF model and confirm that the estimated parameters would be improved. We apply the method to the observational datasets obtained from the Telescope Array and Auger experiments.

Recent results from prototypes of the Fluorescence detector Array of Single-pixel Telescopes (FAST) in both hemispheres

The origin and nature of ultrahigh-energy cosmic rays (UHECRs) are of uppermost importance in astroparticle physics. Motivated by the need for an unprecedented aperture for further advancements, the Fluorescence detector Array of Single-pixel Telescopes (FAST) is a prospective next-generation, ground-based UHECR observatory that aims to cover an enormous area by deploying a large array of low-cost fluorescence telescopes. The full-scale FAST prototype consists of four 20 cm photomultiplier tubes at the focus of a segmented mirror 1.6 m in diameter. Three FAST prototypes have been installed at the Telescope Array Experiment in Utah, USA, and two prototypes at the Pierre Auger Observatory in Mendoza, Argentina, commencing remote observation of UHECRs in both hemispheres. We report on recent results of the full-scale FAST prototypes operated in both hemispheres, including telescope calibrations, atmospheric monitoring, ongoing electronics upgrades, development of sophisticated reconstruction methods and UHECR detections.