Sources Of Cosmic Rays Research Articles

General constraints on sources of high-energy cosmic rays from interaction losses

Sources of high-energy cosmic rays are presently unknown, but can be constrained in various ways. Some of these constraints can be graphically presented on the so-called Hillas diagram. Previous versions of this diagram determined the range of geometrical sizes and magnetic fields of potential astrophysical accelerators, taking into account geometrical criteria and radiation losses. In this work, we update the Hillas diagram for protons, taking into account the losses associated with the $p\ensuremath{\gamma}$ and photopair interactions, relating the allowed regions to the source electromagnetic luminosity. The strongest constraints are obtained for bright compact sources such as the central regions of active galactic nuclei.

Cosmic rays and random magnetic traps.

The spatial distribution of cosmic ray (CR) particles in the interstellar medium (ISM) is of major importance in radio astronomy, where its knowledge is essential for the interpretation of observations, and in theoretical astrophysics, where CRs contribute to the structure and dynamics of the ISM. Local inhomogeneities in interstellar magnetic field strength and structure can affect the local diffusivity and ensemble dynamics of the CR particles. Magnetic traps (regions between magnetic mirrors located on the same magnetic line) can lead to especially strong and persistent features in the CR spatial distribution. Using test particle simulations, we study the spatial distribution of an ensemble of CR particles (both protons and electrons) in various magnetic field configurations, from an idealized axisymmetric trap to those that emerge in intermittent (dynamo-generated) random magnetic fields. We demonstrate that both the inhomogeneity in the CR sources and the energy losses by the CR particles can lead to persistent local inhomogeneities in the CR distribution and that the protons and electrons have different spatial distributions. Our results can have profound implications for the interpretation of the synchrotron emission from astronomical objects, and in particular its random fluctuations.

Effects of the Galactic Magnetic Field on the UHECR Correlation Studies with Starburst Galaxies

We estimate the biases caused by the coherent deflection of cosmic rays due to the Galactic magnetic field (GMF) in maximum likelihood analyses for searches of ultrahigh-energy cosmic-ray sources in the literature. We simulate mock event data sets with a set of assumptions for the starburst galaxy source model, coherent deflection by a GMF model, and mixed-mass composition. We then conduct a maximum likelihood analysis without accounting for the GMF in the same manner as previous studies. We find that the anisotropic fraction f ani is estimated to be systematically lower than the true value. We estimate the true parameters that are compatible with the best-fit parameters that were reported, and find that except for a narrow region with a large anisotropic fraction and a small separation angular scale, a wide parameter space is still compatible with the experimental results. We also develop a maximum likelihood method that takes the GMF model into account and confirm in Monte Carlo simulations that we can estimate the true parameters within a 1σ contour under the ideal condition that we know the event-by-event mass and the GMF.

Curious case of the maximum rigidity distribution of cosmic-ray accelerators

Models of sources of ultra high energy cosmic rays (UHECR)s usually assume that these sources, despite significant variability between them, accelerate particles to the same maximum energy. Here, the authors improve on this assumption by assigning their array of sources a broken-power-law distribution of the maximum energies. They observe that a power law distribution necessitates that the source to source variation be small.

Effects of large-scale magnetic fields on the observed composition of ultrahigh-energy cosmic rays

Ultra high-energy (UHE) cosmic rays (CRs) from distant sources interact with intergalactic radiation fields, leading to their spallation and attenuation. They are also deflected in intergalactic magnetic fields (IGMFs), particularly those associated with Mpc-scale structures. These deflections extend the propagation times of CR particles, forming a magnetic horizon for each CR species. The cumulative cooling and interactions of a CR ensemble also modifies their spectral shape and composition observed on Earth. We construct a transport formulation to calculate the observed UHE CR spectral composition for 4 classes of source population. The effects on CR propagation brought about by IGMFs are modeled as scattering processes during transport, by centers associated with cosmic filaments. Our calculations demonstrate that IGMFs can have a marked effect on observed UHE CRs, and that source population models are degenerate with IGMF properties. Interpretation of observations, including the endorsement or rejection of any particular source classes, thus needs careful consideration of the structural properties and evolution of IGMFs. Future observations providing tighter constraints on IGMF properties will significantly improve confidence in assessing UHE CR sources and their intrinsic CR production properties.

Multiwavelength study of the galactic PeVatron candidate LHAASO J2108+5157

Context. Several new ultrahigh-energy (UHE) γ-ray sources have recently been discovered by the Large High Altitude Air Shower Observatory (LHAASO) collaboration. These represent a step forward in the search for the so-called Galactic PeVatrons, the enigmatic sources of the Galactic cosmic rays up to PeV energies. However, it has been shown that multi-TeV γ-ray emission does not necessarily prove the existence of a hadronic accelerator in the source; indeed this emission could also be explained as inverse Compton scattering from electrons in a radiation-dominated environment. A clear distinction between the two major emission mechanisms would only be made possible by taking into account multi-wavelength data and detailed morphology of the source. Aims. We aim to understand the nature of the unidentified source LHAASO J2108+5157, which is one of the few known UHE sources with no very high-energy (VHE) counterpart. Methods. We observed LHAASO J2108+5157 in the X-ray band with XMM-Newton in 2021 for a total of 3.8 hours and at TeV energies with the Large-Sized Telescope prototype (LST-1), yielding 49 hours of good-quality data. In addition, we analyzed 12 years of Fermi-LAT data, to better constrain emission of its high-energy (HE) counterpart 4FGL J2108.0+5155. We used naima and jetset software packages to examine the leptonic and hadronic scenario of the multi-wavelength emission of the source. Results. We found an excess (3.7σ) in the LST-1 data at energies E > 3 TeV. Further analysis of the whole LST-1 energy range, assuming a point-like source, resulted in a hint (2.2σ) of hard emission, which can be described with a single power law with a photon index of Γ = 1.6 ± 0.2 the range of 0.3 − 100 TeV. We did not find any significant extended emission that could be related to a supernova remnant (SNR) or pulsar wind nebula (PWN) in the XMM-Newton data, which puts strong constraints on possible synchrotron emission of relativistic electrons. We revealed a new potential hard source in Fermi-LAT data with a significance of 4σ and a photon index of Γ = 1.9 ± 0.2, which is not spatially correlated with LHAASO J2108+5157, but including it in the source model we were able to improve spectral representation of the HE counterpart 4FGL J2108.0+5155. Conclusions. The LST-1 and LHAASO observations can be explained as inverse Compton-dominated leptonic emission of relativistic electrons with a cutoff energy of 100−30+70 TeV. The low magnetic field in the source imposed by the X-ray upper limits on synchrotron emission is compatible with a hypothesis of a PWN or a TeV halo. Furthermore, the spectral properties of the HE counterpart are consistent with a Geminga-like pulsar, which would be able to power the VHE-UHE emission. Nevertheless, the lack of a pulsar in the neighborhood of the UHE source is a challenge to the PWN/TeV-halo scenario. The UHE γ rays can also be explained as π0 decay-dominated hadronic emission due to interaction of relativistic protons with one of the two known molecular clouds in the direction of the source. Indeed, the hard spectrum in the LST-1 band is compatible with protons escaping a shock around a middle-aged SNR because of their high low-energy cut-off, but the origin of the HE γ-ray emission remains an open question.

Testing hadronic and photohadronic interactions as responsible for ultrahigh energy cosmic rays and neutrino fluxes from starburst galaxies

We test the hypothesis of starburst galaxies as 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 {\it SimProp}, where hadronic processes in the environment of sources are implemented for the first time. Taking into account a star-formation-rate 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 exceed theoutcome in neutrinos from photo-hadronic interactions in the source environment and in the extra-galactic space.

Calculation of rescaling factors and nuclear multiplication of muons in extensive air showers

Recent results obtained from leading cosmic ray experiments indicate that simulations using LHC-tuned hadronic interaction models underestimate the number of muons in extensive air showers compared to experimental data. This is the so-called muon deficit problem. Determination of the muon component in the air shower is crucial for inferring the mass of the primary particle, which is a key ingredient in the efforts to pinpoint the sources of ultra-high energy cosmic rays. In this paper, we present a new method to derive the muon signal in detectors, which uses the difference between the total reconstructed (‘data’) and simulated signals, and is in turn related to the muon signal which is roughly independent of the zenith angle, but depends on the mass of the primary cosmic ray. Such a method offers an opportunity not only to test/calibrate the hadronic interaction models, but also to derive the β exponent, which describes an increase of the number of muons in a shower as a function of the energy and mass of the primary cosmic ray. Detailed simulations show a dependence of the β exponent on hadronic interaction properties, thus the determination of this parameter is important for understanding the muon deficit problem. We validate the method by using Monte-Carlo simulations for the EPOS-LHC and QGSJetII-04 hadronic interaction models, and show that this method allows us to recover the ratio of the muon signal between EPOS-LHC and QGSJetII-04 and the average β exponent for the studied system, within less than a few percent. This is a consequence of the good recovery of the muon signal for each primary included in the analysis.

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.

Efficient charged particle propagation methods

In astrophysics, the search for sources of the highest-energy cosmic rays continues. For further progress, not only ever better observatories but also ever more realistic numerical simulations are needed. We compare different approaches for numerical test simulations of ultra-high-energy cosmic rays in the intergalactic magnetic field and show that all methods provide correct statistical propagation characteristics of the particles in means of their diffusive behaviour. Through convergence tests, we show that the necessary requirements for the methods differ and ultimately reveal significant differences in the required simulation time.

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.

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.

Revisiting GeV-scale annihilating dark matter with the AMS-02 positron fraction

Antimatter cosmic-rays are used to probe new phenomena in physics, including dark matter annihilation. We use the cosmic-ray positron fraction spectrum by the Alpha Magnetic Spectrometer, to search for such an annihilation signal in the Galaxy. We focus on dark matter with mass between 5 and 120 GeV, producing high-energy electrons and positrons. In these cosmic-ray energies the interplay of multiple astrophysical sources and phenomena, makes this search highly sensitive to the underlying astrophysical background assumptions. We use a vast public library of astrophysical models for the cosmic-ray positron fraction background, to derive robust upper limits on the dark matter's annihilation cross section for a number of annihilation channels. This library accounts for different types of cosmic-ray sources and uncertainties on their distribution in space and time. Also, it accounts for uncertainties on those sources' output, their injected into the interstellar medium cosmic-ray spectra and for uncertainties on cosmic-ray propagation. For any given dark matter particle mass and annihilation channel, upper limits on the annihilation cross section are given by bands that stretch a full order of magnitude in its value. Our work provides weaker limits compared to earlier results, that are however robust to all the relevant astrophysical uncertainties. Between 5 and 15 GeV, we find indications for a possible excess flux of cosmic-ray electrons and positrons. That excess is found for most, but not all of our astrophysical background parameter space, and its significance can vary appreciably. Further scrutiny is necessary to improve the understanding of these lower energy cosmic rays. Finally, we note that even if an excess signal is found in these energies, the current background uncertainties do not allow us to accurately deduce its underlying particle properties.

The depth of the shower maximum of air showers measured with AERA

The Auger Engineering Radio Array (AERA), as part of the Pierre Auger Observatory, is an array of radio antennas probing the nature of ultra-high energy cosmic rays at energies around the transition from Galactic to extra-galactic origin. It measures the MHz radio emission of extensive air showers produced by cosmic rays hitting our atmosphere. The elemental composition of cosmic rays is a crucial piece of information in determining what the sources of cosmic rays are and how cosmic rays are accelerated. This composition can be obtained from the mass-sensitive parameterXmax, the depth of the shower maximum. We reconstructXmaxwith a likelihood analysis, by comparing the measured radio footprint on the ground to an ensemble of footprints from Monte-Carlo CORSIKA/CoREAS air shower simulations. We compare ourXmaxreconstruction with fluorescenceXmaxmeasurements on a per-event basis, a setup unique to the Pierre Auger Observatory, and show the methods to be fully compatible. Furthermore, we extensively validate our reconstruction by identifying and correcting for systematic uncertainties. We determine the resolution of our method as a function of energy and reach a precision better than 15 g cm−2at the highest energies. With a bias-free set of around 600 showers, we find a light to light-mixed composition at energies between 1017.5to 1018.8eV, also in agreement with the Auger fluorescence measurements.

Locating sources of a high energy cosmic ray extensive air shower using HiSPARC data

In the physics community, there is continued interest in high energy cosmic rays and their co-location to stellar sources, including to supernovae. Observational history and corresponding theories link the high energies associated with these types of cosmic rays to those produced in various latter evolution stages of stellar bodies. We hypothesize that high-energy cosmic ray air showers have a source region associated with specific stellar bodies. Here, we present our analysis of the locations of possible high-energy cosmic ray sources by using data from University of Twente cosmic ray detectors that are part of the High School Project on Astrophysics Research with Cosmics (HiSPARC) network. As an example of the potential of the HiSPARC network to aid cosmic ray source location, we analyzed event count data from a high-energy air shower recorded during the first week of January 2014. The data was then processed such that the direction of the air shower could be ascertained to provide the likely area of air shower origin. Angle and time data were then inputted in to Stellarium™, a free open-source planetarium application from which sources could be located. Our results suggest two sources are candidates for the origins of the high-energy air shower: the Cepheid variable Polaris and the binary system, Iota Ursae Majoris A. In sum, we have shown the potential of the HiSPARC network to aid location of stellar sources that may be responsible for detected cosmic air showers.

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.

The transition from Galactic to extragalactic cosmic rays: The high–energy end of the Galactic spectrum

Understanding the transition from Galactic to extragalactic cosmic rays (CRs) is essential to make sense of the Local cosmic ray spectrum. Several models have been proposed to account for this transition in the 0.1 - 10 × 1018 eV range. For instance: ankle models, where the change from a steep Galactic component to a hard extragalactic spectrum occurs in the 4 − 10 × 1018 eV region, dip models, where the interactions of CR protons with the CMB producing electron-positron pairs shape the ankle, or mixed composition models, in which extragalactic CRs are composed of nuclei of various types. In all these scenarios, the low-energy part of the transition involves the high-energy part of the Galactic component. Therefore, any information on the Galactic component, such as maximum energy, chemical composition, and spectrum after propagation, is crucial to understanding the Galactic-extragalactic transition. We briefly review the high-energy part of the CR spectrum expected from the best potential sources of Galactic CRs.

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.

Measurements of Cosmic Ray Mass Composition with the IceCube Neutrino Observatory

The IceCube Neutrino Observatory is a multi-component detector at the South Pole. Besides studying high-energy neutrinos, it is capable of measuring high-energy cosmic rays from PeV to EeV. This energy region is thought to cover the transition from galactic to extragalactic sources of cosmic rays. The observatory consists of the deep in-ice IceCube array, which measures the high-energy (≥500 GeV) muonic component, and the IceTop surface array, which is sensitive to the electromagnetic and low-energy muonic part of an air shower. The primary energy and the mass composition can be measured simultaneously by applying statistical methods including modern machine-learning techniques to reconstruct cosmic ray air showers. In this contribution, we will discuss recent improvements to the reconstruction techniques, the mass composition sensitivity, and an outlook on future improved measurements with the full surface scintillator/radio array to mitigate snow accumulation and measure the air shower maximum Xmax using imaging air-Cherenkov telescopes IceAct.