Cosmic Ray Acceleration Research Articles

Spatial Coincidence between Ultra-high Energy Cosmic Rays and TeV Gamma Rays in the Direction of GRB 980425/SN 1998bw

Gamma-ray bursts (GRBs) have long been suspected as possible ultra-high energy cosmic ray (UHECR) accelerators. In this brief note, I report that GRB 980425/SN 1998bw falls within the region of interest (ROI) with the highest significance in an all-sky blind search for magnetically induced effects in the arrival directions of UHECRs conducted by the Pierre Auger Collaboration with events detected up to to 2018 August 31. There is also report in the literature of delayed TeV emission in archival Fermi LAT observations from the direction of GRB 980425/SN 1998bw. The combined probability that two distinct cosmic ray acceleration signatures in two different multimessenger experiments may appear at the same spatial location by chance is estimated to be between 1.62 × 10−3 and 0.0157.

Further evidence for ultrahigh-energy cosmic ray acceleration in starburst-driven superwinds

Very recently, the Pierre Auger and Telescope Array collaborations reported strong evidence for a correlation between the highest energy cosmic rays and nearby starburst galaxies, with a global significance post-trial of $4.7\ensuremath{\sigma}$. It is well-known that the collective effect of supernovae and winds from massive stars in the central region of these galaxies drives a galactic-scale superwind that can shock heat and accelerate ambient interstellar or circumgalactic gas. In previous work we showed that, for reasonable source parameters, starburst-driven superwinds can be the carriers of ultrahigh-energy cosmic ray acceleration. In this paper we assess the extent to which one can approach the archaeological ``inverse'' problem of deciphering properties of superwind evolution from present-day IR emission of their host galaxies. We show that the Outer Limits galaxy NGC 891 could provide ``smoking-gun evidence'' for the starburst-driven superwind model of ultrahigh-energy cosmic rays.

Solar Cosmic Ray Acceleration by a Shock Wave in the Lower Solar Corona on May 7, 1978

Based on the theory of diffusive shock acceleration of charged particles, we have theoretically studied the spectra of protons produced by a shock wave driven by a coronal mass ejection (CME) in the lower solar corona with known solar plasma parameters for the solar cosmic ray (SCR) event recorded near the Earth’s orbit on May 7, 1978 (ground level enhancement no. 31, GLE31) using numerical methods. The proton flux data from the CPME instrument installed on the IMP-8 spacecraft and the worldwide network of ground-based neutron monitors combined with the particle measurements by the telescopes on the IMP-7 and IMP-8 satellites have been used to analyze the event. In our calculations, to satisfactorily reproduce the proton spectrum measured at the Earth’s orbit, the CME velocity was assumed to be constant and equal to 600 km s–1. The shock acceleration of SCRs proceeded up to 1.4 R⊙ solar radii for 256 s.

Study of the effect of turbulent interstellar medium on the morphology of young supernova remnants

Context. Supernova remnants (SNRs) are one of the main sources of galactic cosmic-ray acceleration. This acceleration, believed to happen at the blast wave front, leads to energy loss at the shock front. This results in the apparent proximity between the blast wave and the contact discontinuity. Aims. In this article, we study the effect that turbulent-like density perturbations of the interstellar medium (ISM) have on the evolution of young SNRs. We focus on the impact these fluctuations have on the structure of SNRs and more precisely on the resulting distance between blast wave and contact discontinuity. As cosmic-ray acceleration is necessary to explain this distance, this study indirectly puts into question the cosmic-ray acceleration at the blast wave front. Methods. We performed a set of purely hydrodynamic three-dimensional simulations without cosmic-ray acceleration in a co-expanding frame. We randomly initialised the density variation of the ISM following a Kolmogorov power law. The resulting ratios of radii between blast wave, contact discontinuity, and reverse shock are then compared to astronomical observations. Results. The addition of density perturbation does not significantly change the average ratio of radii. However, the simulations show a higher growth of interfacial instabilities in the presence of a turbulent ISM. The resulting deformation of the contact discontinuity could explain the proximity between contact discontinuity and blast wave. The deformations also explain the plateau in the radial distribution of the line-of-sight velocity associated with the observations of Tycho. Conclusions. Density perturbations of the ISM should not be neglected in the simulation of young SNRs as they have an impact on the structure of these latter objects that is comparable to that of cosmic-ray acceleration.

Solar Cosmic Ray Acceleration by a Shock Wave in the Lower Solar Corona on May 7, 1978

Solar Cosmic Ray Acceleration by a Shock Wave in the Lower Solar Corona on May 7, 1978

Multiwavelength analysis of Galactic Supernova Remnants

The origin of Galactic Cosmic Rays (CRs) and the possibility of Supernova Remnants (SNRs) being potential CR accelerators is still an open debate. The charged CRs can be detected indirectly by the γ-ray observatories through the π 0 production and consequent decay, leading to the generation of high-energy γ-rays. The goal of the study is to identify qualitative and quantitative trends in favour of hadronic scenario and search for SNRs which could be potential accelerators up to PeV energies (PeVatrons).We have performed a Multiwavelength (MWL) study using different radiative models to evaluate the hadronic contribution. The spectral energy distributions (SEDs) of selected SNRs are modeled using the Naima [1] package. Two different radiative scenarios are considered, pure leptonic and lepto-hadronic scenarios and different methods are used to evaluate their importance.This study shows that the lepto-hadronic scenario is favored for most SNRs. Two particular indicators of hadronic contribution come from the data around the π 0 production threshold and the data above a few TeV. The hard rise at the π 0 production threshold cannot be explained by leptonic processes. More data in this region would be valuable for these studies. For some SNRs, an important hadronic contribution is observed up to a few TeV, thus making them promising PeVatron candidates. In this high-energy region where the leptonic processes are expected to be suppressed, more data is required to help distinguish between the leptonic and hadronic origin of γ-ray emission. In the future, we intend to use the obtained model parameters to simulate data for CTA and assess its capability to identify PeVatrons.

Nonlinear evolution of the Weibel instability with relativistic laser pulses

The Weibel instability is investigated using relativistic intense short laser pulses. A relativistic short laser pulse can generate a sub-relativistic high-density collisionless plasma. By irradiating double parallel planar targets with two relativistic laser pulses, sub-relativistic collisionless counterstreaming plasmas are created. Since the growth rate of the Weibel instability is proportional to the plasma density and velocity, the spatial and temporal scales of the Weibel instability can be much smaller than that from nanosecond large laser facilities. Recent theoretical and numerical studies have revealed that astrophysical collisionless shocks in sub-relativistic regimes in the absence and presence of an ambient magnetic field play essential roles in cosmic ray acceleration. With experimental verification in mind, we discuss the possible experimental models on the Weibel instability with intense short laser pulses. In order to show the experimental feasibility, we perform 2D particle-in-cell simulations in the absence of an external magnetic field as the first step and discuss the optimum conditions to realize the nonlinear evolutions of the Weibel instability in laboratories.

Searches for Neutrinos from Large High Altitude Air Shower Observatory Ultra-high-energy γ-Ray Sources Using the IceCube Neutrino Observatory

Galactic PeV cosmic-ray accelerators (PeVatrons) are Galactic sources theorized to accelerate cosmic rays up to PeV in energy. The accelerated cosmic rays are expected to interact hadronically with nearby ambient gas or the interstellar medium, resulting in γ-rays and neutrinos. Recently, the Large High Altitude Air Shower Observatory (LHAASO) identified 12 γ-ray sources with emissions above 100 TeV, making them candidates for PeVatrons. While at these high energies the Klein–Nishina effect exponentially suppresses leptonic emission from Galactic sources, evidence for neutrino emission would unequivocally confirm hadronic acceleration. Here, we present the results of a search for neutrinos from these γ-ray sources and stacking searches testing for excess neutrino emission from all 12 sources as well as their subcatalogs of supernova remnants and pulsar wind nebulae with 11 yr of track events from the IceCube Neutrino Observatory. No significant emissions were found. Based on the resulting limits, we place constraints on the fraction of γ-ray flux originating from the hadronic processes in the Crab Nebula and LHAASO J2226+6057.

Probing LHAASO galactic PeVatrons through gamma-ray and neutrino correspondence

ABSTRACT Recently, Large High Altitude Air Shower Observatory (LHAASO) has detected several Galactic point sources of ultra high energy (UHE; Eγ > 100 TeV) gamma-rays. These gamma-rays are possibly created in leptonic or hadronic interactions of cosmic rays (CRs) of PeV energies. In the hadronic channel (p − p interaction), the gamma-rays are accompanied by neutrinos. The detection of neutrinos is therefore crucial in understanding CR acceleration in such objects. To estimate the neutrino flux, we adopt the two LHAASO sources (J2226+6057, J1908+0621) found to be spatially associated with the supernova remnants (SNR G106.3+2.7, SNR G40.5–0.5). For these two sources, the detected TeV-PeV gamma-ray spectra are found to be unusually hard (with spectral index ∼1.8). We develop a model of gamma-ray and neutrino emission based on the above two prototypes. The neutrino fluxes from these two sources are found to be below the IceCube sensitivity, but are detectable in upcoming IceCube-Gen2 and KM3NeT experiments. We further estimate the neutrino fluxes from similar other 10 LHAASO PeVatron sources and most of them are found to be detectable in IceCube-Gen2 and KM3NeT. Finally, we explore our model parameters, in particular the spectral power-law index and estimate the future potential of the neutrino detectors to probe CR acceleration in such Galactic sources.

Exploring the role of composition and mass loading on the properties of hadronic jets

ABSTRACTAstrophysical jets are relativistic outflows that remain collimated for remarkably many orders of magnitude. Despite decades of research, the origin of cosmic rays (CRs) remains unclear, but jets launched by both supermassive black holes in the centre of galaxies and stellar-mass black holes harboured in X-ray binaries (BHXBs) are among the candidate sources for CR acceleration. When CRs accelerate in astrophysical jets, they initiate particle cascades that form γ-rays and neutrinos. In the so-called hadronic scenario, the population of accelerated CRs requires a significant amount of energy to properly explain the spectral constraints, similarly to a purely leptonic scenario. The amount of energy required often exceeds the Eddington limit or even the total energy available within the jets. The exact energy source for the accelerated protons is unclear, but due to energy conservation along the jets, it is believed to come from the jet itself via transfer of energy from the magnetic fields or kinetic energy from the outflow. To address this hadronic energy issue and to self-consistently evolve the energy flux along the flows, we explore a novel treatment for including hadronic content, in which instabilities along the jet/wind border play a critical role. We discuss the impact of the different jet compositions on the jet dynamics for a pair dominated and an electron-proton jet and, consequently, the emitted spectrum, accounting for both leptonic and hadronic processes. Finally, we discuss the implications of this mass-loading scenario to address the proton energy issue.

The cosmic ray ionization and γ-ray budgets of star-forming galaxies.

Cosmic rays in star-forming galaxies are a dominant source of both diffuse γ-ray emission and ionization in gas too deeply shielded for photons to penetrate. Though the cosmic rays responsible for γ-rays and ionization are of different energies, they are produced by the same star formation-driven sources, and thus galaxies' star formation rates, γ-ray luminosities, and ionization rates should all be linked. In this paper, we use up-to-date cross-sectiondata to determine this relationship, finding that cosmic rays in a galaxy of star formation rate [Formula: see text] and gas depletion time t dep produce a maximum primary ionization rate ζ ≈ 1×10-16(t dep/Gyr)-1 s-1 and a maximum γ-ray luminosity [Formula: see text] erg s-1 in the 0.1-100 GeV band. These budgets imply either that the ionization rates measured in Milky Way molecular clouds include a significant contribution from local sources that elevate them above the Galactic mean, or that CR-driven ionization in the Milky Way is enhanced by sources not linked directly to star formation. Our results also imply that ionization rates in starburst systems are only moderately enhanced compared to those in the Milky Way. Finally, we point out that measurements of γ-ray luminosities can be used to place constraints on galactic ionization budgets in starburst galaxies that are nearly free of systematic uncertainties on the details of cosmic ray acceleration.

Geminga SNR: Possible Candidate of Local Cosmic-Ray Factory (II)

Accurate measurements of the energy spectrum and anisotropy can help us discover local cosmic-ray accelerators. Our recent works have shown that spectral hardening above 200 GeV in the energy spectra and transition of large-scale anisotropy at ∼100 TeV are of an unifying origin. Less than 100 TeV, both spectral hardening and anisotropy explicitly indicate the dominant contribution from nearby sources. Recent observations of CR anisotropy suggest that this phase is consistent with the locally regular magnetic field (LRMF) of the interstellar boundary explorer (IBEX) below 100 TeV. In this work, we further investigate the parameter space of sources allowed by the observational energy spectra and amplitude and phase of dipole anisotropy. To obtain the best-fit source parameters, a numerical algorithm is to compute the parameter posterior distributions based on Bayesian inference. We found that by combining the observations of the energy spectrum and anisotropy, the parameters of the model can be well constrained. The LRMF and the effect of the corresponding anisotropic diffusion are considered in this work. Finally, the phase results’ right ascension (R.A.)=3.2 h below 100 TeV was obtained by fitting, which is in general agreement with the experimental observations. Since the Geminga SNR is very close to the mean of the fitted parameters, it could be a candidate for a local cosmic-ray accelerator.

Where are Milky Way’s hadronic PeVatrons?

Observations of the Milky Way at TeV--PeV energies reveal a bright diffuse flux of hadronic cosmic rays and also bright point sources of gamma rays. If the gamma-ray sources are hadronic cosmic-ray accelerators, then they must also be neutrino sources. However, no neutrino sources have been detected. Where are they? We introduce a new population-based approach to probe Milky Way hadronic PeVatrons, demanding consistency between diffuse and point-source PeV-range data on cosmic rays, gamma rays, and neutrinos. For the PeVatrons, two extreme scenarios are allowed: (1) the hadronic cosmic-ray accelerators and the gamma-ray sources are the same objects, so that bright neutrino sources exist and improved telescopes can detect them, versus (2) the hadronic cosmic-ray accelerators and the gamma-ray sources are distinct, so that there are no detectable neutrino sources. The latter case is possible if hadronic accelerators have sufficiently thin column densities. We quantify present constraints and future prospects, showing how to reveal the nature of the hadronic PeVatrons.

Revealing ultra-high-energy cosmic ray acceleration with multi-messenger observations of the nearby GRB 980425/SN 1998bw

The origin of ultra-high energy cosmic rays (UHECRs) is one of the most mystifying issues in astroparticle physics. It has been suggested that gamma-ray bursts (GRBs) are excellent acceleration sites for cosmic rays. The propagation of UHECRs from the GRB host galaxy to the Earth should generate delayed secondary photons and neutrinos. Here we present a dedicated search for delayed UHECR and neutrino emission centered around the position of nearby GRB 980425/SN 1998bw. Located at a distance of 36.9 Mpc, GRB 980425/SN 1998bw is well within the Greisen-Zatsepin-Kuzmin (GZK) distance horizon. We find no evidence for UHECR or neutrino clustering around the GRB 980425/SN 1998bw position between 2004 and 2020. Under ideal propagation conditions, we propose that it might be possible to detect an excess from delayed UHECRs around GRB 980425/SN 1998bw within the next 100 years if the intergalactic magnetic field (IGMF) strength is B ≤ 3 × 10-13 G.

Cosmic-Ray Acceleration and Nonthermal Radiation at Accretion Shocks in the Outer Regions of Galaxy Clusters

Cosmology models predict that external accretion shocks form in the outer region of galaxy clusters owing to supersonic gas infall from filaments and voids in the cosmic web. They are characterized by high sonic and Alfvénic Mach numbers, M s ∼ 10–102 and M A ∼ 102–103, and propagate into weakly magnetized plasmas of β ≡ P g /P B ≳ 102. Although strong accretion shocks are expected to be efficient accelerators of cosmic rays (CRs), nonthermal signatures of shock-accelerated CRs around clusters have not been confirmed, and detailed acceleration physics at such shocks has yet to be understood. In this study, we first establish through two-dimensional particle-in-cell simulations that at strong high-β shocks electrons can be pre-energized via stochastic Fermi acceleration owing to the ion Weibel instability in the shock transition region, possibly followed by injection into diffusive shock acceleration. Hence, we propose that the models derived from conventional thermal leakage injection may be employed for the acceleration of electrons and ions at accretion shocks as well. Applying these analytic models to numerical shock zones identified in structure formation simulations, we estimate nonthermal radiation, such as synchrotron and inverse Compton (IC) emission due to CR electrons and π 0-decay γ-rays due to CR protons, around simulated clusters. Our models with injection parameter Q ≈ 3.5–3.8 predict synthetic synchrotron maps, which seem consistent with recent radio observations of the Coma Cluster. However, the detection of nonthermal IC X-rays and γ-rays from accretion shocks would be quite challenging. We suggest that the proposed analytic models may be adopted as generic recipes for CR production at cosmological shocks.

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.

The AMS-02 detector on the ISS - Status and highlights after 11 years on orbit

The Alpha Magnetic Spectrometer, AMS-02, is a magnetic spectrometer detector operating on the International Space Station (ISS) since May the 19th, 2011. More than 200 billion events have been collected by the instrument in the first 11 years of data taking, providing detailed and novel insights on the composition and energy spectra of cosmic rays up to TeV energies. We review the most recent AMS-02 measurements and the advances in the understanding of cosmic ray origin, acceleration and propagation physics.

Ultrahigh-energy cosmic-ray signature in GRB 221009A

The brightest long gamma-ray burst (GRB) detected so far by the Swift-BAT and Fermi-GBM telescopes, GRB 221009A, provides an unprecedented opportunity for understanding the high-energy processes in extreme transient phenomena. We find that the conventional leptonic models for the afterglow emission from this source, synchrotron and synchrotron-self-Compton, have difficulties explaining the observation of ≳10 TeV γ rays (as high as 18 TeV) by the LHAASO detector. We modeled the γ-ray spectrum estimated in the energy range 0.1–1 GeV by the Fermi-LAT detector. The flux predicted by our leptonic models is severely attenuated at > 1 TeV due to γγ pair production with extragalactic background light, and hence an additional component is required at ≳10 TeV. Ultrahigh-energy cosmic rays can be accelerated in the GRB blast wave, and their propagation induces an electromagnetic cascade in the extragalactic medium. The line-of-sight component of this flux can explain the emission at ≳10 TeV detected by LHAASO, which requires a fraction of the GRB blast wave energy to be in ultrahigh-energy cosmic rays. This could be an indication of ultrahigh-energy cosmic-ray acceleration in GRBs.

Cosmic ray acceleration by multiple shocks in the jets of Active Galactic Nuclei

Active galactic nuclei are one of the most promising sources for accelerating particles up to the highest energies. In this contribution, we present a scenario in which cosmic rays are accelerated in multiple shocks created by the interaction of relativistic AGN jets with the winds of embedded massive stars. We solve the Fokker-Planck equation considering escape and radiative losses as well as the collective effect of the shocks and the reacceleration of the particles. Finally, we calculate the maximum energies that the particles can achieve and discuss the possibility of producing ultra-high energy cosmic rays in this astrophysical situation.

Investigating Galactic cosmic rays with <i>γ</i>-ray astronomy

Cosmic rays (CRs) are one of the most important components in the interstellar medium (ISM), and the origin of CRs remains a mystery. The diffusion of CRs in turbulent magnetic fields erases the information on the distribution of CR accelerators to a large extent. The energy dependent diffusion of CRs also significantly modifies the initial (acceleration) spectra of CRs. In this regard, <i>γ</i>-rays, the secondary products of interactions of CRs with gas and photons in the ISM, provide us with more information about the origin of CRs. More specifically, the <i>γ</i>-ray emissions associated with gas, can be used to study the distribution of CRs throughout the Galaxy; discrete <i>γ</i>-ray sources can elucidate the locations of individual CR accelerators. Here, the current status and prospects in these fields are reviewed.