Cosmic Ray Acceleration Research Articles

An X-Ray Shell Reveals the Supernova Explosion for Galactic Microquasar SS 433

How black holes are formed remains an open and fundamental question in astrophysics. Despite theoretical predictions, it lacks observations to understand whether the black hole formation experiences a supernova explosion. Here we report the discovery of an X-ray shell north of the Galactic microquasar SS 433 harboring a stellar-mass black hole spatially associated with radio continuum and polarization emissions and an H i cloud. Its spectrum can be reproduced by a 1 keV underionized plasma, from which the shell is inferred to have been created by a supernova explosion 20–30 kyr ago, and its properties constitute evidence for canonical supernova explosions to create some black holes. Our analysis precludes other possible origins including heated by jets or blown by disk winds. According to the lower mass limit of the compact object in SS 433, we roughly deduced that the progenitor should be more massive than 25 M ⊙. The existence of such a young remnant in SS 433 can also lead to new insights into the supercritical accretion in young microquasars and the γ-ray emission of this system. The fallback ejecta may provide accretion materials within tens of thousands of years, while the shock of the supernova remnant may play a crucial role in the cosmic-ray (re)acceleration.

Constraining the Diffusion Coefficient and Cosmic-Ray Acceleration Efficiency Using Gamma-Ray Emission from the Star-forming Region RCW 38

Stellar winds from massive stars may be significant sources of cosmic rays (CRs). To investigate this connection, we report a detailed study of gamma-ray emission near the young Milky Way star cluster (≈0.5 Myr old) in the star-forming region RCW 38 and compare this emission to its stellar wind properties and diffuse X-ray emission. Using 15 yr of Fermi-LAT data in the 0.2–300 GeV band, we find a significant (σ > 22) detection coincident with the star cluster, producing a total gamma-ray luminosity (extrapolated over 0.1–500 GeV) of L γ =(2.66 ± 0.92) × 1034 erg s−1 adopting a power-law spectral model (Γ = 2.34 ± 0.04). Using an empirical relationship and STARBURST99, we estimate the total wind power to be 8 × 1036 erg s−1, corresponding to a CR acceleration efficiency of η CR ≃ 0.4 for an assumed diffusion coefficient consistent with D = 1028 cm2 s−1. Alternatively, a lower acceleration efficiency of 0.1 can produce this L γ if the diffusion coefficient is smaller, D ≃ 2.5 × 1027 cm2 s−1. Additionally, we analyze Chandra X-ray data from the region and compare the hot-gas pressure to the CR pressure. We find the former is 4 orders of magnitude greater, suggesting that the CR pressure is not dynamically important relative to stellar winds. As RCW 38 is too young for supernovae to have occurred, the high CR acceleration efficiency in RCW 38 demonstrates that stellar winds may be an important source of Galactic CRs.

Search for Joint Multimessenger Signals from Potential Galactic Cosmic-Ray Accelerators with HAWC and IceCube

The origin of high-energy galactic cosmic rays is yet to be understood, but some galactic cosmic-ray accelerators can accelerate cosmic rays up to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding material or radiation, resulting in the production of gamma-rays and neutrinos. To optimize for the detection of such associated production of gamma-rays and neutrinos for a given source morphology and spectrum, a multimessenger analysis that combines gamma-rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework with IceCube Maximum Likelihood Analysis software and HAWC Accelerated Likelihood to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 yr of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that, for five of the sources, the gamma-ray emission observed by HAWC cannot be produced purely from hadronic interactions. We report the limit for the fraction of gamma-rays produced by hadronic interactions for these five sources.

Exploring the intermittency of magnetohydrodynamic turbulence by synchrotron polarization radiation

Context. Magnetohydrodynamic (MHD) turbulence plays a critical role in many key astrophysical processes, such as star formation, acceleration of cosmic rays, and heat conduction. However, its properties are still poorly understood. Aims. In this work, we explore how to extract the intermittency of compressible MHD turbulence from synthetic and real observations. Methods. We used three statistical methods, namely the probability distribution function, kurtosis, and scaling exponent of the multi-order structure function, to reveal the intermittency of MHD turbulence. Results. Our numerical results demonstrate that: (1) the synchrotron polarization intensity statistics can be used to probe the intermittency of magnetic turbulence, by which we can distinguish different turbulence regimes; (2) the intermittency of MHD turbulence is dominated by the slow mode in the sub-Alfvénic turbulence regime; and (3) the Galactic interstellar medium (ISM) in the low latitude region corresponds to the sub-Alfvénic and supersonic turbulence regime. Conclusions.We have successfully measured the intermittency of the Galactic ISM from synthetic and realistic observations.

Origin of the spectral features observed in the cosmic-ray spectrum

Recent measurements reveal the presence of several features in the cosmic-ray (CR) spectrum. In particular, the proton and helium spectra exhibit a spectral hardening at $ GV $ and a spectral steepening at $ TV $, followed by the well-known knee-like feature at $ PV $. The spectra of heavier nuclei also harden at $ GV $, while no claim can be currently made about the presence of the $ TV $ softening, due to low statistics. In addition, the B/C ratio also exhibits a hardening at $ GeV/n $ and seems to be rather shallow at $ TeV/n We propose a possible explanation of the observed spectral features in the framework of a composite diffusion scenario and considering different classes of sources. The proposed scenario is based on two assumptions. First, in the Galactic disk, where magnetic field lines are mainly oriented along the Galactic plane, particle scattering is assumed to be very inefficient. Therefore, the transport of CRs from the disk to the halo is set by the magnetic field line random walk induced by large-scale turbulence. Second, we propose that the spectral steepening at $ TV $ is related to the typical maximum rigidity reached in the acceleration of CRs by the majority of supernova remnants, while we assume that only a fraction of sources, contributing to $ 10-20<!PCT!> $ of the CR population, can accelerate particles up to sim PV rigidities. Within this framework we show that it is possible to reproduce the proton and helium spectra from GV to multi-PV; the $p$/He ratio; the spectra of CRs from lithium to iron; the $ p $ flux and the $ p /p$ ratio; and the abundance ratios B/C, B/O, C/O, Be/C, Be/O, and Be/B. We also discuss the Be Be ratio in view of the recent AMS02 preliminary measurements.

Detection of extended gamma-ray emission in the vicinity of Cl Danks 1 and 2

Abstract We report the detection of high-energy gamma-ray emission towards the G305 star-forming region. Using almost 15 years of observation data from Fermi Large Area Telescope, we detected an extended gamma-ray source in this region with a significance of ∼13σ. The gamma-ray radiation reveals a clear pion-bump feature and can be fitted with the power law parent proton spectrum with an index of −2.5. The total cosmic ray (CR) proton energy in the gamma-ray production region is estimated to be the order of $10^{49}\ \rm erg$. We further derived the CR radial distribution from both the gamma-ray emission and gas distribution and found it roughly obeys the 1/r type profile, though a constant profile is not ruled out. This is consistent with other similar systems and expected from the continuous injection of CRs by the central powerful young massive star cluster Danks 1 or Danks 2 in this region. Together with former detections of similar gamma-ray structures, such as Cygnus cocoon, Westerlund 1, Westerlund 2, NGC 3603, and W40, the detection supports the hypothesis that young massive star clusters are CR accelerators.

A Magnetized Strongly Turbulent Corona as the Source of Neutrinos from NGC 1068

The cores of active galactic nuclei are potential accelerators of 10–100 TeV cosmic rays, in turn producing high-energy neutrinos. This picture was confirmed by the compelling evidence of a TeV neutrino signal from the nearby active galaxy NGC 1068, leaving open the question of what is the site and mechanism of cosmic-ray acceleration. One candidate is the magnetized turbulence surrounding the central supermassive black hole. Recent particle-in-cell simulations of magnetized turbulence indicate that stochastic cosmic-ray acceleration is nonresonant, in contrast to the assumptions of previous studies. We show that this has important consequences on a self-consistent theory of neutrino production in the corona, leading to a more rapid cosmic-ray acceleration than previously considered. The turbulent magnetic-field fluctuations needed to explain the neutrino signal are consistent with a magnetically powered corona. We find that strong turbulence, with turbulent magnetic energy density higher than 1% of the rest-mass energy density, naturally explains the normalization of the IceCube neutrino flux, in addition to the neutrino spectral shape. Only a fraction of the protons in the corona, which can be directly inferred from the neutrino signal, are accelerated to high energies. Thus, in this framework, the neutrino signal from NGC 1068 provides a testbed for particle acceleration in magnetized turbulence.

The remarkable microquasar S26: A super-Eddington PeVatron

S26 is an extragalactic microquasar with the most powerful jets ever discovered. They have a kinetic luminosity of $L_ j $. This implies that the accretion power to the black hole should be super-Eddington, of the order of $L_ acc j $. However, the observed X-ray flux of this system indicates an apparent very sub-Eddington accretion luminosity of $L_ X We aim to characterize the nature of S26, explain the system emission, and study the feasibility of super-Eddington microquasars as potential PeVatron sources. We first analyze multi-epoch X-ray observations of S26 obtained with XMM-Newton and model the super-Eddington disk and its wind. We then develop a jet model and study the particle acceleration and radiative processes that occur in shocks generated near the base of the jet and in its terminal region. We find that the discrepancy between the jet and the apparent disk luminosities in S26 is caused by the complete absorption of the disk radiation by the wind ejected from the super-Eddington disk. The non-thermal X-rays are produced near the base of the jet, and the thermal X-rays are emitted in the terminal regions. The radio emission observed with the Australia Telescope Compact Array can be explained as synchrotron radiation produced at the reverse shock in the lobes. We also find that S26 can accelerate protons to PeV energies in both the inner jet and the lobes. The ultra-high energy protons accelerated in the lobes of S26 are injected into the interstellar medium with a total power of $ We conclude that S26 is a super-Eddington microquasar with a dense disk-driven wind that obscures the X-ray emission from the inner disk, and that the supercritical nature of the system allows the acceleration of cosmic rays to PeV energies.

On the Correlation between Young Massive Star Clusters and Gamma-Ray Unassociated Sources

Star clusters (SCs) are potential cosmic-ray accelerators and therefore are expected to emit high-energy radiation. However, a clear detection of gamma-ray emission from this source class has only been possible for a handful of cases. This could in principle result from two different reasons: either detectable SCs are limited to a small fraction of the total number of Galactic SCs, or gamma-ray-emitting SCs are not recognized as such and therefore are listed in the ensemble of unidentified sources. In this Letter we investigate this latter scenario by comparing available catalogs of SCs and H ii regions, obtained from Gaia and Wide-field Infrared Survey Explorer observations, to the gamma-ray GeV and TeV catalogs built from Fermi Large Area Telescope (LAT), H.E.S.S., and LHAASO data. The significance of the correlation between catalogs is evaluated by comparing the results with simulations of synthetic populations. A strong correlation emerges between Fermi-LAT-unidentified sources and H ii regions that trace massive SCs in the earliest (≲1–2 Myr) phase of their life, where no supernova explosions have happened yet, confirming that winds of massive stars can alone accelerate particles and produce gamma-ray emission at least up to GeV energies. The association with TeV energy sources is less evident. Similarly, no significant association is found between Gaia SCs and GeV nor TeV sources. We ascribe this fact to the larger extension of these objects but also to an intrinsic bias in the Gaia selection toward SCs surrounded by a lower target gas density, which would otherwise hinder the detection in the optical wave band.

Cosmic-ray acceleration in core-collapse supernova remnants with the wind termination shock

Cosmic-ray acceleration in core-collapse supernova remnants with the wind termination shock

Diagnostics of Magnetohydrodynamic Modes in the Interstellar Medium through Synchrotron Polarization Statistics

One of the biggest challenges in understanding magnetohydrodynamic (MHD) turbulence is identifying the plasma mode components from observational data. Previous studies on synchrotron polarization from the interstellar medium (ISM) suggest that the dominant MHD modes can be identified via statistics of Stokes parameters, which would be crucial for studying various ISM processes such as the scattering and acceleration of cosmic rays, star formation, and dynamo. In this paper, we present a numerical study of the synchrotron polarization analysis (SPA) method through systematic investigation of the statistical properties of the Stokes parameters. We derive the theoretical basis for our method from the fundamental statistics of MHD turbulence, recognizing that the projection of the MHD modes allows us to identify the modes dominating the energy fraction from synchrotron observations. Based on the discovery, we revise the SPA method using synthetic synchrotron polarization observations obtained from 3D ideal MHD simulations with a wide range of plasma parameters and driving mechanisms, and present a modified recipe for mode identification. We propose a classification criterion based on a new SPA+ fitting procedure, which allows us to distinguish between Alfvén mode and compressible/slow mode dominated turbulence. We further propose a new method to identify fast modes by analyzing the asymmetry of the SPA+ signature and establish a new asymmetry parameter to detect the presence of fast mode turbulence. Additionally, we confirm through numerical tests that the identification of the compressible and fast modes is not affected by Faraday rotation in both the emitting plasma and the foreground.

Probing particle acceleration in Abell 2256: From 16 MHz to gamma rays

Merging galaxy clusters often host spectacular diffuse radio synchrotron sources. These sources can be explained by a non-thermal pool of relativistic electrons that are accelerated by shocks and turbulence in the intracluster medium. The origin of the pool and details of the cosmic ray transport and acceleration mechanisms in clusters are still open questions. Due to the often extremely steep spectral indices of diffuse radio emission, it is best studied at low frequencies. However, the lowest frequency window available to ground-based telescopes (10−30 MHz) has remained largely unexplored as radio frequency interference and calibration problems related to the ionosphere become severe. Here, we present LOFAR observations from 16 to 168 MHz targeting the famous cluster Abell 2256. In the deepest-ever images at decametre wavelengths, we detected and resolved the radio halo, radio shock, and various steep spectrum sources. We measured standard single power-law behaviour for the radio halo and radio shock spectra, with spectral indices of α = −1.56 ± 0.02 from 24 to 1500 MHz and α = −1.00 ± 0.02 from 24 to 3000 MHz, respectively. Additionally, we found significant spectral index and curvature fluctuations across the radio halo, indicating an inhomogeneous emitting volume. In contrast to the straight power-law spectra of the large-scale diffuse sources, the various AGN-related sources showed extreme steepening towards higher frequencies and flattening towards low frequencies. We also discovered a new fossil plasma source with a steep spectrum between 23 and 144 MHz, with α = −1.9 ± 0.1. Finally, by comparing radio and gamma-ray observations, we ruled out purely hadronic models for the radio halo origin in Abell 2256, unless the magnetic field strength in the cluster is exceptionally high, which is unsupportable by energetic arguments and inconsistent with the knowledge of other cluster magnetic fields.

Interplay between the non-resonant streaming instability and self-generated pressure anisotropies

ABSTRACT The non-thermal particles escaping from collision-less shocks into the surrounding medium can trigger a non-resonant streaming instability that converts parts of their drift kinetic energy into large amplitude magnetic field perturbations, and promote the confinement and acceleration of high energy cosmic rays. We present simulations of the instability using an hybrid-Particle-in-Cell approach including Monte Carlo collisions, and demonstrate that the development of the non-resonant mode is associated with important ion pressure anisotropies in the background plasma. Depending on the initial conditions, the anisotropies may act on the instability by lowering its growth and trigger secondary micro-instabilities. Introducing collisions with neutrals yield a strong reduction of the magnetic field amplification as predicted by linear fluid theory. In contrast, Coulomb collisions in fully ionized plasmas are found to mitigate the self-generated pressure anisotropies and promote the growth of the magnetic field.

Different spectra of cosmic ray H, He, and heavier nuclei escaping compact star clusters

ABSTRACT Cosmic ray acceleration at the termination shock of compact star clusters has recently received much attention, mainly because of the detection of gamma-ray emission from some of such astrophysical sources. Here we focus on the acceleration of nuclei at the termination shock and we investigate the role played by proton energy losses and spallation reactions of nuclei, especially downstream of the shock. We show that for a reasonable choice of the mean gas density in the cavity excavated by the cluster wind, dominated by the presence of dense clouds, the spectrum of He nuclei escaping the bubble is systematically harder than the spectrum of hydrogen, in a manner that appears to be qualitatively consistent with the observed and yet unexplained phenomenon of discrepant hardening. We also find that, in this scenario, the spallation reactions of heavier nuclei are likely to be so severe that their spectra becomes very hard and with a low normalization, meaning that it is unlikely that heavy nuclei escaping star clusters can provide a sizeable contribution to the spectrum of cosmic rays at the Earth. Limitations and implications of this scenario are discussed.

Observational Evidence for Magnetic Field Amplification in SN 1006

We report the first observational evidence for magnetic field amplification in the northeast/southwest (NE/SW) shells of supernova remnant SN 1006, one of the most promising sites of cosmic ray acceleration. In previous studies, the strength of magnetic fields in these shells was estimated to be B SED ≃ 25 μG from the spectral energy distribution, where the synchrotron emission from relativistic electrons accounted for radio to X-rays, along with the inverse Compton emission extending from the GeV to TeV energy bands. However, the analysis of broadband radio data, ranging from 1.37 to 100 GHz, indicated that the radio spectrum steepened from α 1 = 0.52 ± 0.02 to α 2 = 1.34 ± 0.21 by Δα = 0.85 ± 0.21. This is naturally interpreted as a cooling break under a strong magnetic field of B brk ≥ 2 mG. Moreover, the high-resolution MeerKAT image indicated that the width of the radio NE/SW shells was broader than that of the X-ray shell by a factor of only 3−20, as measured by Chandra. Such narrow radio shells can be naturally explained if the magnetic field responsible for the radio emissions is B R ≥ 2 mG. Assuming that the magnetic field is locally enhanced by a factor of approximately a = 100 along the NE/SW shells, we argue that the filling factor, which is the volume ratio of such a magnetically enhanced region to that of the entire shell, must be as low as approximately k = 2.5 × 10−5.

Evidence for particle acceleration approaching PeV energies in the W51 complex

The γ-ray emission from the W51 complex is widely acknowledged to be attributed to the interaction between the cosmic rays (CRs) accelerated by the shock of supernova remnant (SNR) W51C and the dense molecular clouds in the adjacent star-forming region, W51B. However, the maximum acceleration capability of W51C for CRs remains elusive. Based on observations conducted with the Large High Altitude Air Shower Observatory (LHAASO), we report a significant detection of γ rays emanating from the W51 complex, with energies from 2 to 200 TeV. The LHAASO measurements, for the first time, extend the γ-ray emission from the W51 complex beyond 100 TeV and reveal a significant spectrum bending at tens of TeV. By combining the “π0-decay bump” featured data from Fermi-LAT, the broadband γ-ray spectrum of the W51 region can be well-characterized by a simple pp-collision model. The observed spectral bending feature suggests an exponential cutoff at ∼400 TeV or a power-law break at ∼200 TeV in the CR proton spectrum, most likely providing the first evidence of SNRs serving as CR accelerators approaching the PeV regime. Additionally, two young star clusters within W51B could also be theoretically viable to produce the most energetic γ rays observed by LHAASO. Our findings strongly support the presence of extreme CR accelerators within the W51 complex and provide new insights into the origin of Galactic CRs.

Very-high-energy γ-Ray Emission from Young Massive Star Clusters in the Large Magellanic Cloud

The Tarantula Nebula in the Large Magellanic Cloud is known for its high star formation activity. At its center lies the young massive star cluster R136, providing a significant amount of the energy that makes the nebula shine so brightly at many wavelengths. Recently, young massive star clusters have been suggested to also efficiently produce very high-energy cosmic rays, potentially beyond PeV energies. Here, we report the detection of very-high-energy γ-ray emission from the direction of R136 with the High Energy Stereoscopic System, achieved through a multicomponent, likelihood-based modeling of the data. This supports the hypothesis that R136 is indeed a very powerful cosmic-ray accelerator. Moreover, from the same analysis, we provide an updated measurement of the γ-ray emission from 30 Dor C, the only superbubble detected at TeV energies presently. The γ-ray luminosity above 0.5 TeV of both sources is (2–3) × 1035 erg s−1. This exceeds by more than a factor of 2 the luminosity of HESS J1646−458, which is associated with the most massive young star cluster in the Milky Way, Westerlund 1. Furthermore, the γ-ray emission from each source is extended with a significance of >3σ and a Gaussian width of about 30 pc. For 30 Dor C, a connection between the γ-ray emission and the nonthermal X-ray emission appears likely. Different interpretations of the γ-ray signal from R136 are discussed.

Exploring cosmic ray energy loss mechanisms: Insights from Bethe-Bloch, modified bethe-bloch, and inverse compton scattering equations

This study delves into the intricate dynamics of cosmic ray energy loss mechanisms, employing three distinct mathematical approaches: the Bethe-Bloch equation, the modified Bethe-Bloch equation, and the inverse compton scattering equation. The Bethe-Bloch equation elucidates the ionization and bremsstrahlung processes governing cosmic ray interactions with plasma mediums, providing foundational insights into energy loss phenomena. Expanding upon this, the modified Bethe-Bloch equation integrates additional factors such as medium density and magnetic fields, refining our understanding of cosmic ray propagation in diverse astrophysical environments. Moreover, the inverse Compton scattering equation uncovers the energy loss mechanisms associated with cosmic ray interactions with photons in the interstellar medium, contributing essential insights into high-energy astrophysical phenomena. Through mathematical derivations, numerical computations, and graphical analyses, our investigation reveals the intricate dependencies of final cosmic ray energy on various parameters, including initial energy, target photon density, and magnetic field strength. Furthermore, our findings have broader implications for astrophysical processes such as gamma-ray emission and cosmic-ray acceleration mechanisms. By elucidating the complex interplay between cosmic rays and their surrounding environments, this study advances our understanding of high-energy astrophysical phenomena and provides a framework for future research in this field.

Early-timeγ-ray constraints on cosmic-ray acceleration in the core-collapse SN 2023ixf with theFermiLarge Area Telescope

Context.While supernova remnants (SNRs) have been considered the most relevant Galactic cosmic ray (CR) accelerators for decades, core-collapse supernovae (CCSNe) could accelerate particles during the earliest stages of their evolution and hence contribute to the CR energy budget in the Galaxy. Some SNRs have indeed been associated with TeVγ-rays, yet proton acceleration efficiency during the early stages of an SN expansion remains mostly unconstrained.Aims.The multi-wavelength observation of SN 2023ixf, a Type II supernova (SN) in the nearby galaxy M 101 (at a distance of 6.85 Mpc), opens the possibility to constrain CR acceleration within a few days after the collapse of the red super-giant stellar progenitor. With this work, we intend to provide a phenomenological, quasi-model-independent constraint on the CR acceleration efficiency during this event at photon energies above 100 MeV.Methods.We performed a maximum-likelihood analysis ofγ-ray data from theFermiLarge Area Telescope up to one month after the SN explosion. We searched for high-energy, non-thermal emission from its expanding shock, and estimated the underlying hadronic CR energy reservoir assuming a power-law proton distribution consistent with standard diffusive shock acceleration.Results.We do not find significantγ-ray emission from SN 2023ixf. Nonetheless, our non-detection provides the first limit on the energy transferred to the population of hadronic CRs during the very early expansion of a CCSN.Conclusions.Under reasonable assumptions, our limits would imply a maximum efficiency on the CR acceleration of as low as 1%, which is inconsistent with the common estimate of 10% in generic SNe. However, this result is highly dependent on the assumed geometry of the circumstellar medium, and could be relaxed back to 10% by challenging spherical symmetry. Consequently, a more sophisticated, inhomogeneous characterisation of the shock and the progenitor’s environment is required before establishing whether or not Type II SNe are indeed efficient CR accelerators at early times.

On the Width of a Collisionless Shock and the Index of the Cosmic Rays It Accelerates

Despite being studied for many years, the structure of collisionless shocks is still not fully determined. Such shocks are known to be accelerators of cosmic rays (CRs), which, in turn, modify the shock structure. The shock width λ is known to be connected to the CR spectral index a. Here, we use an instability analysis to derive the shock width in the presence of CRs. We obtain an analytical expression connecting the shock width to the CR index and to the fraction of upstream particles that are accelerated. We find that when this fraction becomes larger than ∼30%, a new instability becomes dominant. The shock undergoes a transition where its width increases by a factor ∼8–10, and the CR acceleration effectively ends. Our analysis is valid for strong, nonrelativistic, and unmagnetized shocks. We discuss the implication of these results on the expected range of CR spectra and flux observed and on the structure of nonrelativistic collisionless shocks.