Shock Acceleration Mechanism Research Articles

Studying the role of supernova remnants as cosmic ray PeVatron with LHAASO observations

It is commonly believed that galactic cosmic rays are produced in supernova remnants (SNRs) and accelerated via a diffusive shock acceleration (DSA) mechanism in supernova blast waves driven by expanding SNRs. The latest theoretical advancement of the diffusive shock acceleration hypothesis in SNRs shows that cosmic rays may be accelerated up to the knee energy of the observed cosmic ray spectrum under the amplified magnetic field scenario. There is, however, no empirical evidence to support SNRs as sources of hadrons with energies larger than a few tens of TeV. Very recently, LHAASO observatory reported the very high-energy gamma ray emission between TeV to PeV energy range from two SNRs, Cassiopeia A and IC 443. Above 25 TeV energies, non-detection of gamma-ray flux by LHAASO yields a strong upper limit. In this work, we investigate the implications of the acceleration of cosmic ray protons in the SNR on energetic gamma rays produced in the hadronic interaction of cosmic rays with ambient matter. Our findings imply that when we consider the highest attainable energy of cosmic ray protons in the SNRs to be about 100 TeV, the observed gamma-ray spectra from the two SNRs can be described consistently. Therefore, we conclude that the Cassiopia A and IC 443 SNRs are unlikely to be cosmic ray PeVatrons.

Supernova Remnants in Gamma Rays

In the 1960s, the remnants of supernova explosions (SNRs) were indicated as a possible source of galactic cosmic rays through the Diffusive Shock Acceleration (DSA) mechanism. Since then, the observation of gamma-ray emission from relativistic ions in these objects has been one of the main goals of high-energy astrophysics. A few dozen SNRs have been detected at GeV and TeV photon energies in the last two decades. However, these observations have shown a complex phenomenology that is not easy to reduce to the standard paradigm based on DSA acceleration. Although the understanding of these objects has greatly increased, and their nature as efficient electron and proton accelerators has been observed, it remains to be clarified whether these objects are the main contributors to galactic cosmic rays. Here, we review the observations of γ-ray emission from SNRs and the perspectives for the future.

Self-consistent modeling of the energetic storm particle event of November 10, 2012

Context. It is thought that solar energetic ions associated with coronal and interplanetary shock waves are accelerated to high energies by the diffusive shock acceleration mechanism. For this mechanism to be efficient, intense magnetic turbulence is needed in the vicinity of the shock. The enhanced turbulence upstream of the shock can be produced self-consistently by the accelerated particles themselves via streaming instability. Comparisons of quasi-linear-theory-based particle acceleration models that include this process with observations have not been fully successful so far, which has motivated the development of acceleration models of a different nature. Aims. Our aim is to test how well our self-consistent quasi-linear SOLar Particle Acceleration in Coronal Shocks (SOLPACS) simulation code, developed earlier to simulate proton acceleration in coronal shocks, models the particle foreshock region. Methods. We applied SOLPACS to model the energetic storm particle (ESP) event observed by the STEREO A spacecraft on November 10, 2012. Results. All but one main input parameter of SOLPACS are fixed by the in situ plasma measurements from the spacecraft. By comparing a simulated proton energy spectrum at the shock with the observed one, we were able to fix the last simulation input parameter related to the efficiency of particle injection to the acceleration process. A subsequent comparison of simulated proton time-intensity profiles in a number of energy channels with the observed ones shows a very good correspondence throughout the upstream region. Conclusions. Our results strongly support the quasi-linear description of the foreshock region.

An X-rays-to-radio investigation of the nuclear polarization from the radio-galaxy Centaurus A

ABSTRACT Centaurus A is one of the closest radio galaxies to Earth. Its proximity allowed us to extensively study its active galactic nucleus but the core emission mechanism remains elusive because of local strong dust and gas obscuration. The capability of polarimetry to shave-off contaminating emission has been exploited without success in the near-infrared by previous studies but the very recent measurement of the 2–8 keV polarization by the Imaging X-ray Polarimetry Explorer (IXPE) brought the question back to the fore. To determine what is the prevalent photon generation mechanism to the multiwavelength emission from the core of Centaurus A, we retrieved from the archives the panchromatic polarization measurements of the central compact component. We built the total and polarized flux spectral energy distributions of the core and demonstrated that synchrotron self-Compton models nicely fit the polarized flux from the radio to the X-ray band. The linear polarization of the synchrotron continuum is perpendicular to the jet radio axis from the optical to the radio band, and parallel to it at higher energies. The observed smooth rotation of the polarization angle in the ultraviolet band is attributed to synchrotron emission from regions that are getting closer to the particle acceleration site, where the orientation of the jet’s magnetic fields become perpendicular to the jet axis. This phenomenon support the shock acceleration mechanism for particle acceleration in Centaurus A, in line with IXPE observations of several high-synchrotron peak blazars.

Acceleration of suprathermal protons near an interplanetary shock

Context.Interplanetary collisionless shocks are known to be sources of energetic charged particles up to hundreds of MeV. However, the underlying acceleration mechanisms are still under debate.Aims.We determine the properties of suprathermal protons accelerated by the interplanetary shock on 2021 November 3 with the unprecedented high-resolution measurements by the SupraThermal Electron Proton sensor of the Energetic Particle Detector on board the Solar Orbiter spacecraft, in order to constrain the potential shock acceleration mechanisms.Methods.We first reconstructed the pitch-angle distributions (PADs) of suprathermal protons in the solar wind frame. Then, we studied the evolution of the PADs, the temporal flux profile, and the velocity distribution function of this proton population close to the shock and compared the observations to theoretical predictions.Results.We find that the suprathermal proton fluxes peak ∼12 to ∼24 s before the shock in the upstream region. The proton fluxes rapidly decrease by ∼50% in a thin layer (∼8000 km) adjacent to the shock in the downstream region and become constant farther downstream. Furthermore, the proton velocity distribution functions in the upstream (downstream) region fit a double power law,f(v)∼v−γ, at ∼1000 − 3600 km s−1, with aγof ∼3.4 ± 0.2 (∼4.3 ± 0.7) at velocities (v) below a break at ∼1800 ± 100 km s−1(∼1600 ± 200 km s−1) and aγof ∼5.8 ± 0.3 (∼5.8 ± 0.2) at velocities higher than this. These indices are all smaller than predicted by first-order Fermi acceleration. In addition, the proton PADs in the covered pitch-angle range show anisotropies in the direction away from the shock in the region close to the upstream region and become nearly isotropic farther upstream, while downstream of the shock, they show a tendency of anisotropies towards 90° PA.Conclusions.These results suggest that the acceleration of suprathermal protons at interplanetary shocks are dynamic on a timescale of ∼10 s, that is, few proton gyroperiods. Furthermore, shock-drift acceleration likely plays an important role in accelerating these suprathermal protons.

Dynamical Coupling between Anomalous Cosmic Rays and Solar Wind in Outer Heliosphere

Voyager 2 (V2) observed that pickup ions (PUIs) and anomalous cosmic rays (ACRs) have a significant influence on the solar wind structures in the outer heliosphere. In particular, the maxima in energetic particle intensities often lagged behind the shock front, while the flow velocity in some cases featured a precursor in front of the shock. These two effects are believed to be caused by the cease of ACR injection from PUIs at large heliocentric distances, and the backward diffusion of ACRs, respectively. This paper investigates the dynamical coupling between the ACRs and the solar wind in the outer heliosphere by means of a time-dependent numerical MHD simulation, in which the ACRs are treated as a massless fluid that only contributes its pressure to the system. Two types of inner boundary conditions were used, namely a synthetic short-term shock event and an extended period of data-driven solar wind variability based on the OMNI database. The lag of the ACR pressure maximum relative to the shock front, and the extended shock precursor were reproduced by the numerical results. The increase rate of the lag is related to the corresponding diffusion coefficient and the injection efficiency from PUIs to ACRs. The model was also applied to the termination shock, where simulations likewise showed that the peak in the ACR distribution can be located a short distance downstream of the shock front, indicating that the time-dependent diffusive shock acceleration mechanism is a candidate to interpret the lag between the ACR pressure peak and the shock front observed by V2.

Numerical Modeling of Spectral Hardening at a Finite-width Shock

Spectral hardening has been identified in solar flare hard X-ray observations for several decades and remains a puzzle. We examine spectral hardening under the diffusive shock acceleration mechanism using numerical simulations. The hardening is related to the finite width of the shock and is controlled by the shock Péclet number. We implement two different types of Monte Carlo simulations. The first is based on the backward stochastic differential equation method, where the Parker transport equation is solved by casting it to a set of stochastic different equations, and by following the trajectories of individual quasiparticles. In the second approach, we follow real particles and particles are assumed to move freely between scatterings from magnetic turbulence in the plasma. The scattering is modeled as either large-angle hard-sphere elastic collision, or small-angle pitch-angle scattering. We show that the results from these two approaches agree well with each other and agree with analytical results. We also use a Pan-spectrum form to fit the resulting spectra.

High-yield and high-angular-fluence neutron generation from deuterons accelerated by laser-driven collisionless shock

A bright collimated neutron source is an essential tool for global security missions and fundamental scientific research. In this paper, we study a compact high-yield and high-angular-fluence neutron source particularly suitable for high-energy neutron applications utilizing the breakup reaction of laser-driven deuterons in a 9Be converter. The neutron generation scaling from such a reaction is used to guide the choice and optimization of the acceleration process for bulk ions in a low density CD2 foam. In particular, the collisionless shock acceleration mechanism is exploited with proper choice in the laser and target parameter space to accelerate these ions toward energies above the temperature of the distribution. Particle-in-cell and Monte Carlo simulations are coupled here to investigate this concept and possible adverse effects as well as the contribution from the surface ions accelerated and the optimal converter design. The simulation results indicated that our design can be a practical approach to increase both the neutron yield and angular fluence of laser-driven neutron sources, reaching >1011 neutron/pulse (or >108 neutron/J) and >1011 neutron/sr (or >108 neutron/sr/J) with present-day kJ-class high-power lasers. Such developments will advance fundamental neutron science, high precision radiography, and other global security applications with laser-driven sources.

SN 2008iy circumstellar interaction: bright and lesser light effects

ABSTRACT Optical photometry and spectra of the luminous type IIn supernova SN 2008iy are analysed in detail with implications for cosmic ray acceleration and the radio emission. The light curve and expansion velocities indicate ejecta with the kinetic energy of 3 × 1051 erg to collide with the ∼10 M⊙ circumstellar envelope. The luminous H α is explained as originated primarily from circumstellar clouds interacting with the forward shock. For the first time the fluorescent O i 8446 Å emission is used to demonstrate that the cloud fragmentation cascade spans a scale range >2.3 dex. The narrow circumstellar H α permitted us to estimate the acceleration efficiency of cosmic rays. The found value is close to the efficiency inferred in the same way for other two SNe IIn, SN 1997eg and SN 2002ic. The efficiency of cosmic ray acceleration is utilized to reproduce the radio flux from SN 2008iy for the amplified magnetic field consistent with the saturated turbulent magnetic field in the diffusive shock acceleration mechanism.

Backstreaming Pickup Ions

Abstract Ions that are reflected at the shock front and escape back into the upstream region can play the role of ions that start to be accelerated by a diffusive shock acceleration mechanism. Backstreaming ions have been shown to be generated from a superthermal tail of the solar wind at sufficiently high upstream temperatures. The number of such ions was found to be low and they were not found at shock angles exceeding 50°. The mechanism of production is multiple reflection when an ion changes the direction of motion inside the ramp for the first time, due to the cross-shock potential. Since pickup ions (PUIs) constitute a strongly superthermal population of protons a substantially stronger production of backstreaming PUIs can be expected. We study the reflection of PUIs in a planar stationary shock front using test particle analysis. The used model is inspired by the observed profile of the termination shock. The influence of magnetic compression, the shock angle, and the overshoot are analyzed. It is found that generation of backstreaming PUIs in this shock is substantially more efficient than the generation of backstreaming protons from thermal solar wind. The fraction of backstreaming PUIs rapidly increases with the increase of magnetic compression and the decrease of the shock angle. Overshoot enhances production of backstreaming PUIs and allows it for larger shock angles. No backstreaming ions have been found for shock angles larger than 60°. The results of the test particle analysis are supported by full-particle simulations.

New high-frequency radio observations of the Cygnus Loop supernova remnant with the Italian radio telescopes

ABSTRACT Supernova remnants (SNRs) represent a powerful laboratory to study the cosmic ray acceleration processes at shocks, and their relation to the properties of the circumstellar medium. With the aim of studying the high-frequency radio emission and investigating the energy distribution of accelerated electrons and the magnetic field conditions, we performed single-dish observations of the large and complex Cygnus Loop SNR from 7.0–24.8 GHz with the Medicina and Sardinia Radio Telescopes, focusing on the northern filament (NGC 6992) and the southern shell. Both regions show a spectrum well fitted by a power-law function (S ∝ ν−α), with spectral index α = 0.45 ± 0.05 for NGC 6992 and α = 0.49 ± 0.01 for the southern shell and without any indication of a spectral break. The spectra are significantly flatter than the whole Cygnus Loop spectrum (α = 0.54 ± 0.01), suggesting a departure from the plain shock acceleration mechanisms, which for NGC 6992 could be related to the ongoing transition towards a radiative shock. We model the integrated spectrum of the whole SNR considering the evolution of the maximum energy and magnetic field amplification. Through the radio spectral parameters, we infer a magnetic field at the shock of 10 μG. This value is compatible with purely adiabatic compression of the interstellar magnetic field, suggesting that the amplification process is currently inefficient.

On the Challenges of Cosmic-Ray Proton Shock Acceleration in the Intracluster Medium

Galaxy clusters host the largest particle accelerators in the Universe: Shock waves in the intracluster medium (ICM), a hot and ionised plasma, that accelerate particles to high energies. Radio observations pick up synchrotron emission in the ICM, proving the existence of accelerated cosmic-ray electrons. However, a sign of cosmic-ray protons, in form of γ-rays. remains undetected. This is know as the missingγ-ray problem and it directly challenges the shock acceleration mechanism at work in the ICM.Over the last decade, theoretical and numerical studies focused on improving our knowledge on the microphysics that govern the shock acceleration process in the ICM. These new models are able to predict a γ-ray signal, produced by shock accelerated cosmic-ray protons, below the detection limits set modern γ-ray observatories. In this review, we summarise the latest advances in solving the missing γ-ray problem.

Fast Particle Acceleration at Perpendicular Shocks with Uniform Upstream Magnetic Field and Strong Downstream Turbulence

Abstract Shock waves produce relativistic particles via the diffusive shock acceleration (DSA) mechanism. Among various circumstances, fast acceleration has been expected for perpendicular shocks. We investigate the acceleration time and the energy spectrum of particles accelerated at a perpendicular shock. In our model, the upstream perpendicular magnetic field has no fluctuation, and the downstream region is highly turbulent. Then, the particle motion is the gyration in the upstream region and Bohm-like diffusion downstream. Under this situation, we derive an analytical form of the acceleration time. Using test particle simulations, the validity of our formula is verified. In addition, the energy spectrum of particles is the same as those predicted by standard DSA. Therefore, the presently proposed mechanism simultaneously realizes the rapid acceleration and the canonical spectrum, dN/dp ∝ p −2, even if there is no upstream magnetic amplification.

MMS Observations of Accelerated Interstellar Pickup He+ Ions at an Interplanetary Shock

Abstract Interstellar pickup ions (PUI) are interstellar neutrals that have been ionized in transit through the heliosphere via charge exchange or photoionization. These new PUIs then “freeze” into the solar wind (SW) and move with the bulk SW velocity (V SW). They also begin to gyrate around the local magnetic field with V SW, resulting in a maximum PUI velocity of 2 ∗ V SW. Understanding how He+ PUIs are accelerated at shocks in space provides valuable insights into shock dynamics and shock acceleration mechanisms. On 2018 January 8, while the Magnetospheric Multiscale (MMS) Mission was in the SW, it observed the forward shock of a corotating interaction region. In this work, we analyze the upstream and downstream velocity distributions of He+ for this quasiperpendicular, marginally supercritical (θ Bn = 67°, M A = 2.8) interplanetary shock. We derive average two-dimensional pitch-angle distributions in the field-aligned SW frame, as well as reduced one-dimensional velocity distributions for selected upstream and downstream intervals. We find that the shock accelerates He+ largely perpendicular to the magnetic field, consistent with shock reflection. Furthermore, we derive a measured ratio of accelerated He+ ions by estimating a maximum speed in the downstream frame, above which we assume that He+ ions have been accelerated by their interaction with the shock. The measured acceleration ratio, R Measured = 0.14, is then compared to a theoretical ratio derived from a simple model of the upstream distribution, R Theoretical = 0.06.

Extragalactic cosmic rays diffusing from two populations of sources

We consider the possibility of explaining the observed spectrum and composition of the cosmic rays with energies above ${10}^{17}\text{ }\text{ }\mathrm{eV}$ in terms of two different extragalactic populations of sources in the presence of a turbulent intergalactic magnetic field (including also a fading galactic cosmic-ray component). The populations are considered to be the superposition of different nuclear species having rigidity-dependent spectra. The first extragalactic population is dominant in the energy range ${10}^{17}\ensuremath{-}{10}^{18}\text{ }\text{ }\mathrm{eV}$ and consists of sources having a relatively large density ($>{10}^{\ensuremath{-}3}\text{ }\text{ }{\mathrm{Mpc}}^{\ensuremath{-}3}$) and a steep spectrum. The second extragalactic population dominates the cosmic-ray flux above a few EeV; it has a harder spectral slope and has a high-energy cutoff at a few $Z$ EeV (where $eZ$ is the associated cosmic-ray charge). This population has a lower density of sources ($<{10}^{\ensuremath{-}4}\text{ }\text{ }{\mathrm{Mpc}}^{\ensuremath{-}3}$), so that the typical intersource separation is larger than few tens of Mpc, being significantly affected by a magnetic horizon effect that strongly suppresses its flux for energies below $\ensuremath{\sim}Z$ EeV. We discuss how this scenario could be reconciled with the values of the cosmic-ray source spectral indices that are expected to result from the diffusive shock acceleration mechanism.

Shock-accelerated cosmic rays and streaming instability in the adaptive mesh refinement code Ramses

Cosmic rays (CRs) are thought to play a dynamically important role in several key aspects of galaxy evolution, including the structure of the interstellar medium, the formation of galactic winds, and the non-thermal pressure support of halos. We introduce a numerical model solving for the CR streaming instability and acceleration of CRs at shocks with a fluid approach in the adaptive mesh refinement code RAMSES. CR streaming is solved with a diffusion approach and its anisotropic nature is naturally captured. We introduce a shock finder for the RAMSES code that automatically detects shock discontinuities in the flow. Shocks are the loci for CR injection, and their efficiency of CR acceleration is made dependent on the upstream magnetic obliquity according to the diffuse shock acceleration mechanism. We show that the shock finder accurately captures shock locations and estimates the shock Mach number for several problems. The obliquity-dependent injection of CRs in the Sedov solution leads to situations where the supernova bubble exhibits large polar caps (homogeneous background magnetic field), or a patchy structure of the CR distribution (inhomogeneous background magnetic field). Finally, we combine both accelerated CRs with streaming in a simple turbulent interstellar medium box, and show that the presence of CRs significantly modifies the structure of the gas.

Particle acceleration in interstellar shocks

This review presents the fundamentals of the particle acceleration processes active in interstellar medium (ISM), which are essentially based on the so-called Fermi mechanism theory. More specifically, the review presents here in more details the first order Fermi acceleration process -- also known as diffusive shock acceleration (DSA) mechanism. In this case, acceleration is induced by the interstellar (IS) shock waves. These IS shocks are mainly associated with emission nebulae (HII regions, planetary nebulae and supernova remnants). Among all types of emission nebulae, the strongest shocks are associated with supernova remnants (SNRs). Due to this fact they also provide the most efficient manner to accelerate ISM particles to become {high energy particles}, i.e.~cosmic-rays (CRs). The review therefore focuses on the particle acceleration at the strong shock waves of supernova remnants.

Origin and impacts of the first cosmic rays

Nonthermal phenomena are ubiquitous in the Universe, and cosmic rays (CRs) play various roles in different environments. When, where, and how CRs are first generated since the Big Bang? We argue that blast waves from the first cosmic explosions at z~20 lead to Weibel mediated nonrelativistic shocks and CRs can be generated by the diffusive shock acceleration mechanism. We show that protons are accelerated at least up to sub-GeV energies, and the fast velocity component of supernova ejecta is likely to allow CRs to achieve a few GeV in energy. We discuss other possible accelerators of the first CRs, including accretion shocks due to the cosmological structure formation. These CRs can play various roles in the early universe, such as the ionization and heating of gas, the generation of magnetic fields, and feedbacks on the galaxy formation.

Cosmic-ray Spectrum Steepening in Supernova Remnants. I. Loss-free Self-similar Solution

Abstract The direct measurements of cosmic rays (CRs), after correction for the propagation effects in the interstellar medium, indicate that their source spectra are likely to be significantly steeper than the canonical E −2 spectrum predicted by the standard diffusive shock acceleration (DSA) mechanism. The DSA has long been held responsible for the production of galactic CRs in supernova remnant (SNR) shocks. The γ-ray “probes” of the acceleration spectra of CRs on the spot, inside the SNRs, lead to the same conclusion. We show that the steep acceleration spectrum can be attributed to the combination of (i) spherical expansion, (ii) tilting of the magnetic field along the shock surface, and (iii) shock deceleration. Because of (i) and (ii), the DSA is efficient only on two “polar caps” of a spherical shock where the local magnetic field is within ≃45° to its normal. The shock-produced spectrum observed edge-on steepens with the particle energy because the number of freshly accelerated particles with lower energies continually adds up to a growing acceleration region. We demonstrate the steepening effect by obtaining an exact self-similar solution for the particle acceleration at an expanding shock surface with an arbitrary energy dependence of particle diffusivity κ. We show that its increase toward higher energy steepens the spectrum, which deeply contrasts with the standard DSA spectrum where κ cancels out.

High-energy Emission from Interacting Supernovae: New Constraints on Cosmic-Ray Acceleration in Dense Circumstellar Environments

Abstract Supernovae (SNe) with strong interactions with circumstellar material (CSM) are promising candidate sources of high-energy neutrinos and gamma-rays and have been suggested as an important contributor to Galactic cosmic rays (CRs) beyond 1015 eV. Taking into account the shock dissipation by a fast velocity component of SN ejecta, we present comprehensive calculations of the nonthermal emission from SNe powered by shock interactions with a dense wind or CSM. Remarkably, we consider electromagnetic cascades in the radiation zone and subsequent attenuation in the pre-shock CSM. A new time-dependent phenomenological prescription provided by this work enables us to calculate gamma-ray, hard X-ray, radio, and neutrino signals, which originate from CRs accelerated by the diffusive shock acceleration (DSA) mechanism. We apply our results to Type IIn SN 2010jl and Type Ib/IIn SN 2014C, for which the model parameters can be determined from the multiwavelength data. For SN 2010jl, the more promising case, by using the the latest Fermi Large Area Telescope Pass 8 data release, we derive new constraints on the CR energy fraction, ϵ p ≲ 0.05–0.1. We also find that the late-time radio data of these interacting SNe are consistent with our model. Further multimessenger and multiwavelength observations of nearby interacting SNe should give us new insights into the DSA in dense environments, as well as pre-SN mass-loss mechanisms.