Cosmic Ray Acceleration Research Articles

A Cosmic Ray Acceleration Mechanism Based on Background Flow Velocity Inhomogeneities Yielding Power-law Spectra

In astrophysics, one significant challenge lies in understanding the acceleration of cosmic rays, which leads to the occurrence of a power law. In this article, momentum transport generated by the combined effects of pitch-angle diffusion and background flow velocity inhomogeneities is proposed to obtain a cosmic rays acceleration mechanism, starting from the well-known focused transport equation describing particle diffusion and acceleration. The inhomogeneities of background flow velocity are ubiquitous in the astrophysical environment. The equation for the isotropic part of the distribution function of charged energetic particles is derived, and its solution is obtained, demonstrating the form of momentum power laws of cosmic rays. In addition, if it is assumed that cosmic rays penetrate compressive MHD waves or turbulence, for quasi-steady states, the spectral index δ of the momentum power law spectrum of cosmic rays is found to be in the range [−5, −3], which includes the observed power law indices of galactic cosmic rays. The results obtained in this article demonstrate that the mechanism proposed in this article, along with shock acceleration, may also contribute to the acceleration of galactic cosmic rays. Furthermore, when momentum convection effect and higher-order momentum derivative terms are considered, the indices of power laws should be smaller than −5. This may explain the power laws of solar energetic particle events.

Solar Flares and the Aether Dynamic

This study investigates the dynamics of Solar flares using the "closed fluid dynamic principle" proposed by Muhammad Aslam Musakhail, which reinforces the pre-Einsteinian aether theory. By defining the "aether force" as the difference between relativistic total mass and rest mass, the principle provides a novel perspective on the relativistic transitions in Solar flare phenomena. The analysis builds on Parker's force-free model and Melrose’s resistive slab theory, extending them to describe the dual-phase dynamics of Solar flares: the Alfvénic phase, where massive fermions (v<c) propagate as Alfvén waves, and the heat-dissipation phase, where fermions transition to massless states (v=c), driven by current sources. Key results reveal that energy dissipation scales with resistivity, and the temporal evolution of Poynting flux and magnetic fields highlights distinct transitions between the phases. Numerical simulations demonstrate the exponential decay of axial magnetic fields and the helical organization of flux tubes. These findings validate theoretical predictions and provide insights into particle acceleration, magnetic reconnection, and energy transfer mechanisms during Solar flares. This work not only advances the understanding of Solar flare energetics but also establishes a mathematical framework that can be extended to other astrophysical phenomena, such as cosmic ray acceleration and stellar flares. Future studies should focus on numerical validation and observational testing to refine the dual-phase model and its broader applicability.

Numerical simulations of laser-driven experiments of ion acceleration in stochastic magnetic fields

We present numerical simulations used to interpret laser-driven plasma experiments at the GSI Helmholtz Centre for Heavy Ion Research. The mechanisms by which non-thermal particles are accelerated in astrophysical environments, e.g., the solar wind, supernova remnants, and gamma ray bursts, is a topic of intense study. When shocks are present, the primary acceleration mechanism is believed to be first-order Fermi, which accelerates particles as they cross a shock. Second-order Fermi acceleration can also contribute, utilizing magnetic mirrors for particle energization. Despite this mechanism being less efficient, the ubiquity of magnetized turbulence in the universe necessitates its consideration. Another acceleration mechanism is the lower-hybrid drift instability, arising from gradients of both density and magnetic field, which produce lower-hybrid waves with an electric field that energizes particles as they cross these waves. With the combination of high-powered laser systems and particle accelerators, it is possible to study the mechanisms behind cosmic-ray acceleration in the laboratory. In this work, we combine experimental results and high-fidelity three-dimensional simulations to estimate the efficiency of ion acceleration in a weakly magnetized interaction region. We validate the FLASH magneto-hydrodynamic code with experimental results and use OSIRIS particle-in-cell code to verify the initial formation of the interaction region, showing good agreement between codes and experimental results. We find that the plasma conditions in the experiment are conducive to the lower-hybrid drift instability, yielding an increase in energy ΔE of ∼ 264 keV for 242 MeV calcium ions.

High-Statistics Measurement of the Cosmic-Ray Electron Spectrum with H.E.S.S.

Owing to their rapid cooling rate and hence loss-limited propagation distance, cosmic-ray electrons and positrons (CRe) at very high energies probe local cosmic-ray accelerators and provide constraints on exotic production mechanisms such as annihilation of dark matter particles. We present a high-statistics measurement of the spectrum of CRe candidate events from 0.3 to 40TeV with the High Energy Stereoscopic System, covering 2 orders of magnitude in energy and reaching a proton rejection power of better than 10^{4}. The measured spectrum is well described by a broken power law, with a break around 1TeV, where the spectral index increases from Γ_{1}=3.25±0.02(stat)±0.2(sys) to Γ_{2}=4.49±0.04(stat)±0.2(sys). Apart from the break, the spectrum is featureless. The absence of distinct signatures at multi-TeV energies imposes constraints on the presence of nearby CRe accelerators and the local CRe propagation mechanisms.

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.

The remarkable microquasar S26: A super-Eddington PeVatron

Context. S26 is an extragalactic microquasar with the most powerful jets ever discovered. They have a kinetic luminosity of Lj ∼ 5 × 1040 erg s−1. This implies that the accretion power to the black hole should be super-Eddington, of the order of Lacc ∼ Lj. However, the observed X-ray flux of this system indicates an apparent very sub-Eddington accretion luminosity of LX ≈ 1037 erg s−1. Aims. We aim to characterize the nature of S26, explain the system emission, and study the feasibility of super-Eddington microquasars as potential PeVatron sources. Methods. 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. Results. 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 nonthermal 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 ∼1036 erg s−1. Conclusions. 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.

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 yr of observation data from Fermi Large Area Telescope, we detected an extended gamma-ray source in this region with a significance of $\sim 13 \sigma$. 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 of 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.

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.