Shock Acceleration Mechanism Research Articles

IceCube’s astrophysical neutrino energy spectrum from CPT violation

The 6-year dataset of high-energy starting events (HESE) at IceCube indicates a spectrum of astrophysical neutrinos much softer than expected from the Fermi shock acceleration mechanism. On the other hand, IceCube's up-going muon neutrino dataset and Fermi-LAT's gamma-ray spectrum point to an $E^{-2}$ neutrino spectrum. If the HESE data above 200 TeV are fit with the latter flux, an excess at lower energies ensues, which then suggests a multicomponent spectrum. We show that the HESE dataset can be explained by a single $E^{-2}$ power-law neutrino flux from a muon-damped $p\gamma$ source if neutrino interactions are modified by CPT violation. The low-energy excess is naturally explained by the pileup of events from superluminal neutrino decay, and there is no cutoff at high energies due to the contribution of subluminal antineutrinos. The best-fit scenario with CPT violation also predicts the observation of Glashow resonance events in the near future.

Modeling a Single SEP Event from Multiple Vantage Points Using the iPATH Model

Abstract Using the recently extended 2D improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, we model an example gradual solar energetic particle event as observed at multiple locations. Protons and ions that are energized via the diffusive shock acceleration mechanism are followed at a 2D coronal mass ejection-driven shock where the shock geometry varies across the shock front. The subsequent transport of energetic particles, including cross-field diffusion, is modeled by a Monte Carlo code that is based on a stochastic differential equation method. Time intensity profiles and particle spectra at multiple locations and different radial distances, separated in longitudes, are presented. The results shown here are relevant to the upcoming Parker Solar Probe mission.

Cosmic Ray Production in Supernovae

We give a brief review of the origin and acceleration of cosmic rays (CRs), emphasizing the production of CRs at different stages of supernova evolution by the first-order Fermi shock acceleration mechanism. We suggest that supernovae with trans-relativistic outflows, despite being rather rare, may accelerate CRs to energies above $10^{18}\mbox{ eV}$ over the first year of their evolution. Supernovae in young compact clusters of massive stars, and interaction powered superluminous supernovae, may accelerate CRs well above the PeV regime. We discuss the acceleration of the bulk of the galactic CRs in isolated supernova remnants and re-acceleration of escaped CRs by the multiple shocks present in superbubbles produced by associations of OB stars. The effects of magnetic field amplification by CR driven instabilities, as well as superdiffusive CR transport, are discussed for nonthermal radiation produced by nonlinear shocks of all speeds including trans-relativistic ones.

Analysis and interpretation of inner-heliospheric SEP events with the ESA Standard Radiation Environment Monitor (SREM) onboard the INTEGRAL and Rosetta Missions

Using two heliospheric vantage points, we study 22 solar energetic particle (SEP) events, 14 of which were detected at both locations. SEP proton events were detected during the declining phase of solar cycle 23 (November 2003–December 2006) by means of two nearly identical Standard Radiation Environment Monitor (SREM) units in energies ranging between 12.6 MeV and 166.3 MeV. In this work we combine SREM data with diverse solar and interplanetary measurements, aiming to backtrace solar eruptions from their impact in geospace (i.e., from L1 Lagrangian point to Earth’s magnetosphere) to their parent eruptions at the Sun’s low atmosphere. Our SREM SEP data support and complement a consistent inner-heliospheric description of solar eruptions (solar flares and coronal mass ejections [CMEs]) and their magnetospheric impact. In addition, they provide useful information on the understanding of the origin, acceleration, and propagation of SEP events at multi-spacecraft settings. All SEP events in our sample originate from major eruptions consisting of major (>M-class) solar flares and fast (>1800 km/s, on average), overwhelmingly (>78%) halo, CMEs. All but one SEP event studied are unambiguously associated with shock-fronted CMEs, suggesting a CME-driven shock acceleration mechanism. Moreover, a significant correlation is found between the SEP event peak and the onset of the storm sudden commencement, that might help improve prediction of magnetospheric disturbances. In general, SEP events correlate better with interplanetary (i.e., in-situ; L1-based) than with solar eruption features. Our findings support (a) the routine use of cost-effective SREM units, or future improvements thereof, for the detection of SEP events and (b) their implementation in multi-spacecraft settings as a means to improve both the physical understanding of SEP events and their forecasting.

Acceleration of collimated 45\u2009MeV protons by collisionless shocks driven in low-density, large-scale gradient plasmas by a 1020\u2009W/cm2, 1\u2009\xb5m laser

A new type of proton acceleration stemming from large-scale gradients, low-density targets, irradiated by an intense near-infrared laser is observed. The produced protons are characterized by high-energies (with a broad spectrum), are emitted in a very directional manner, and the process is associated to relaxed laser (no need for high-contrast) and target (no need for ultra-thin or expensive targets) constraints. As such, this process appears quite effective compared to the standard and commonly used Target Normal Sheath Acceleration technique (TNSA), or more exploratory mechanisms like Radiation Pressure Acceleration (RPA). The data are underpinned by 3D numerical simulations which suggest that in these conditions a Low Density Collisionless Shock Acceleration (LDCSA) mechanism is at play, which combines an initial Collisionless Shock Acceleration (CSA) to a boost procured by a TNSA-like sheath field in the downward density ramp of the target, leading to an overall broad spectrum. Experiments performed at a laser intensity of 1020 W/cm2 show that LDCSA can accelerate, from ~1% critical density, mm-scale targets, up to 5 × 109 protons/MeV/sr/J with energies up to 45(±5) MeV in a collimated (~6° half-angle) manner.

Modeling Particle Acceleration and Transport at a 2‐D CME‐Driven Shock

AbstractWe extend our earlier Particle Acceleration and Transport in the Heliosphere (PATH) model to study particle acceleration and transport at a coronal mass ejection (CME)‐driven shock. We model the propagation of a CME‐driven shock in the ecliptic plane using the ZEUS‐3D code from 20 solar radii to 2 AU. As in the previous PATH model, the initiation of the CME‐driven shock is simplified and modeled as a disturbance at the inner boundary. Different from the earlier PATH model, the disturbance is now longitudinally dependent. Particles are accelerated at the 2‐D shock via the diffusive shock acceleration mechanism. The acceleration depends on both the parallel and perpendicular diffusion coefficients κ|| and κ⊥ and is therefore shock‐obliquity dependent. Following the procedure used in Li, Shalchi, et al. (), we obtain the particle injection energy, the maximum energy, and the accelerated particle spectra at the shock front. Once accelerated, particles diffuse and convect in the shock complex. The diffusion and convection of these particles are treated using a refined 2‐D shell model in an approach similar to Zank et al. (). When particles escape from the shock, they propagate along and across the interplanetary magnetic field. The propagation is modeled using a focused transport equation with the addition of perpendicular diffusion. We solve the transport equation using a backward stochastic differential equation method where adiabatic cooling, focusing, pitch angle scattering, and cross‐field diffusion effects are all included. Time intensity profiles and instantaneous particle spectra as well as particle pitch angle distributions are shown for two example CME shocks.

Collimated protons accelerated from an overdense gas jet irradiated by a 1\u2009\xb5m wavelength high-intensity short-pulse laser

We have investigated proton acceleration in the forward direction from a near-critical density hydrogen gas jet target irradiated by a high intensity (1018 W/cm2), short-pulse (5 ps) laser with wavelength of 1.054 μm. We observed the signature of the Collisionless Shock Acceleration mechanism, namely quasi-monoenergetic proton beams with small divergence in addition to the more commonly observed electron-sheath driven proton acceleration. The proton energies we obtained were modest (~MeV), but prospects for improvement are offered through further tailoring the gas jet density profile. Also, we observed that this mechanism is very robust in producing those beams and thus can be considered as a future candidate in laser-driven ion sources driven by the upcoming next generation of multi-PW near-infrared lasers.

A possible explanation of the knee of cosmic light component spectrum from 100 TeV to 3 PeV**Supported by the National Natural Science Foundation of China (11433004, 11363006, 11103016, 11173020), Top Talents Program of Yunnan Province (2015HA030) and the Natural Science Foundation of Yunnan Province(2015FB103)

A mixed hydrogen and helium (H + He) spectrum with a clear steepening at ∼ 700 TeV has been detected by the ARGO-YBJ experiments. In this paper, we demonstrate that the observed H + He spectrum can be reproduced well with a model of cosmic rays escaping from the supernova remnants (SNRs) in our Galaxy. In this model, particles are accelerated in a SNR through a non-linear diffusive shock acceleration mechanism. Three components of high energy light nuclei escaped from the SNR are considered. It should be noted that the proton spectrum observed by KASCADE can be also explained by this model given a higher acceleration efficiency.

Suzakuobservations of the merging galaxy cluster Abell 2255: The northeast radio relic

We present the results of deep 140 ks Suzaku X-ray observations of the north-east (NE) radio relic of the merging galaxy cluster Abell2255. The temperature structure of Abell2255 is measured out to 0.9 times the virial radius (1.9 Mpc) in the NE direction for the first time. The Suzaku temperature map of the central region suggests a complex temperature distribution, which agrees with previous work. Additionally, on a larger-scale, we confirm that the temperature drops from 6 keV around the cluster center to 3 keV at the outskirts, with two discontinuities at {\it r}$\sim$5\arcmin~(450 kpc) and $\sim$12\arcmin~(1100 kpc) from the cluster center. Their locations coincide with surface brightness discontinuities marginally detected in the XMM-Newton image, which indicates the presence of shock structures. From the temperature drop, we estimate the Mach numbers to be ${\cal M}_{\rm inner}\sim$1.2 and, ${\cal M}_{\rm outer}\sim$1.4. The first structure is most likely related to the large cluster core region ($\sim$350--430 kpc), and its Mach number is consistent with the XMM-Newton observation (${\cal M}\sim$1.24: Sakelliou & Ponman 2006). Our detection of the second temperature jump, based on the Suzaku key project observation, shows the presence of a shock structure across the NE radio relic. This indicates a connection between the shock structure and the relativistic electrons that generate radio emission. Across the NE radio relic, however, we find a significantly lower temperature ratio ($T_1/T_2\sim1.44\pm0.16$ corresponds to~${\cal M}_{\rm X-ray}\sim1.4$) than the value expected from radio wavelengths, based on the standard diffusive shock acceleration mechanism ($T_1/T_2>$ 3.2 or ${\cal M}_{\rm Radio}>$ 2.8).

Positive charge prevalence in cosmic rays: Room for dark matter in the positron spectrum

The unexpected energy spectrum of the positron/electron ratio is interpreted astrophysically, with a possible exception of the 100-300 GeV range. The data indicate that this ratio, after a decline between $0.5-8$ GeV, rises steadily with a trend towards saturation at 200-400GeV. These observations (except for the trend) appear to be in conflict with the diffusive shock acceleration (DSA) mechanism, operating in a \emph{single} supernova remnant (SNR) shock. We argue that $e^{+}/e^{-}$ ratio can still be explained by the DSA if positrons are accelerated in a \emph{subset} of SNR shocks which: (i) propagate in clumpy gas media, and (ii) are modified by accelerated CR \emph{protons}. The protons penetrate into the dense gas clumps upstream to produce positrons and, \emph{charge the clumps positively}. The induced electric field expels positrons into the upstream plasma where they are shock-accelerated. Since the shock is modified, these positrons develop a harder spectrum than that of the CR electrons accelerated in other SNRs. Mixing these populations explains the increase in the $e^{+}/e^{-}$ ratio at $E>8$ GeV. It decreases at $E<8$ GeV because of a subshock weakening which also results from the shock modification. Contrary to the expelled positrons, most of the antiprotons, electrons, and heavier nuclei, are left unaccelerated inside the clumps. Scenarios for the 100-300 GeV AMS-02 fraction exceeding the model prediction, including, but not limited to, possible dark matter contribution, are also discussed.

The γ $\gamma$ -ray sources in classical novae

The interaction between novae ejecta and accretion disks or circumbinary (CB) disks surrounding nova systems based on the diffusive shock acceleration (DSA) mechanism are investigated. Results suggest that interaction between the novae ejecta and CB disks most probably produce $\gamma$ -rays. Both leptonic and hadronic scenarios can feasibly produce $\gamma$ -rays in CB disks. A grid is calculated in order to discuss the $\gamma$ -ray sources in classical novae (CNe), and results indicated that the mass of a white dwarf and the orbit period in CNe can greatly affect the production of $\gamma$ -rays. According to the proposed grid and synthesis population method, it is estimated that the percentage of CNe that can produce $\gamma$ -rays ranges from approximately 0.17 % to 51 % for leptonic scenarios when the magnetic field strength in the shock region ranges from 0.005 G to 0.02 G. The occurrence rate of $\gamma$ -ray CNe ranges from approximately $1~\mbox{yr}^{-1} \mbox{ to } 27~\mbox{yr}^{-1}$ in the galaxy. The corresponding percentage ranges from approximately 64 % to 97 % in hadronic scenarios, while the occurrence rate of $\gamma$ -ray CNe ranges from approximately $35~\mbox{yr}^{-1} \mbox{ to } 52~\mbox{yr}^{-1}$ in the galaxy.

PLASMA COMPRESSION IN MAGNETIC RECONNECTION REGIONS IN THE SOLAR CORONA

ABSTRACT It has been proposed that particles bouncing between magnetized flows converging in a reconnection region can be accelerated by the first-order Fermi mechanism. Analytical considerations of this mechanism have shown that the spectral index of accelerated particles is related to the total plasma compression within the reconnection region, similarly to the case of the diffusive shock acceleration mechanism. As a first step to investigate the efficiency of Fermi acceleration in reconnection regions in producing hard energy spectra of particles in the solar corona, we explore the degree of plasma compression that can be achieved at reconnection sites. In particular, we aim to determine the conditions for the strong compressions to form. Using a two-dimensional resistive MHD numerical model, we consider a set of magnetic field configurations where magnetic reconnection can occur, including a Harris current sheet, a force-free current sheet, and two merging flux ropes. Plasma parameters are taken to be characteristic of the solar corona. Numerical simulations show that strong plasma compressions (≥4) in the reconnection regions can form when the plasma heating due to reconnection is efficiently removed by fast thermal conduction or the radiative cooling process. The radiative cooling process that is negligible in the typical 1 MK corona can play an important role in the low corona/transition region. It is found that plasma compression is expected to be strongest in low-beta plasma β ∼ 0.01–0.07 at reconnection magnetic nulls.

Protostars: Forges of cosmic rays?

Galactic cosmic rays (CR) are particles presumably accelerated in supernova remnant shocks that propagate in the interstellar medium up to the densest parts of molecular clouds, losing energy and their ionisation efficiency because of the presence of magnetic fields and collisions with molecular hydrogen. Recent observations hint at high levels of ionisation and at the presence of synchrotron emission in protostellar systems, which leads to an apparent contradiction. We want to explain the origin of these CRs accelerated within young protostars as suggested by observations. Our modelling consists of a set of conditions that has to be satisfied in order to have an efficient CR acceleration through diffusive shock acceleration. We analyse three main acceleration sites, then we follow the propagation of these particles through the protostellar system up to the hot spot region. We find that jet shocks can be strong accelerators of CR protons, which can be boosted up to relativistic energies. Other promising acceleration sites are protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate CR protons. In contrast, accretion flow shocks are too weak to efficiently accelerate CRs. Though CR electrons are weakly accelerated, they can gain a strong boost to relativistic energies through re-acceleration in successive shocks. We suggest a mechanism able to accelerate both CR protons and electrons through the diffusive shock acceleration mechanism, which can be used to explain the high ionisation rate and the synchrotron emission observed towards protostellar sources. The existence of an internal source of energetic particles can have a strong and unforeseen impact on the ionisation of the protostellar disc, on the star and planet formation processes, and on the formation of pre-biotic molecules.

Turbulent cosmic ray reacceleration and the curved radio spectrum of the radio relic in the Sausage Cluster

Abstract It has often been thought that the northern radio relic in the galaxy cluster CIZA J2242.8+5301 (the “Sausage” Cluster) is associated with cosmic ray (CR) electrons that are accelerated at a shock through the diffusive shock acceleration (DSA) mechanism. However, recent radio observations have shown that the radio spectrum is curved, which is inconsistent with the prediction of a simple DSA model. Moreover, the CR electron spectrum before being affected by radiative cooling seems to be too hard for DSA. In this study, we show that these facts are natural consequences if the electrons are reaccelerated in turbulence downstream of the shock. In this model, DSA is not the main mechanism for generating high-energy electrons. We find that the mean free path of the electrons should be much shorter than the Coulomb mean free path for efficient reacceleration. The scale of the turbulent eddies must be smaller than the width of the relic. We also predict hard X-ray spectra of inverse Compton scattering of photons.

TURBULENT COSMIC-RAY REACCELERATION AT RADIO RELICS AND HALOS IN CLUSTERS OF GALAXIES

Radio relics are synchrotron emission found on the periphery of galaxy clusters. From the position and the morphology, it is often believed that the relics are generated by cosmic-ray (CR) electrons accelerated at shocks through a diffusive shock acceleration (DSA) mechanism. However, some radio relics have harder spectra than the prediction of the standard DSA model. One example is observed in the cluster 1RXS J0603.3+4214, which is often called the ``Toothbrush Cluster''. Interestingly, the position of the relic is shifted from that of a possible shock. In this study, we show that these discrepancies in the spectrum and the position can be solved if turbulent (re)acceleration is very effective behind the shock. This means that for some relics turbulent reacceleration may be the main mechanism to produce high-energy electrons, contrary to the common belief that it is the DSA. Moreover, we show that for efficient reacceleration, the effective mean free path of the electrons has to be much smaller than their Coulomb mean free path. We also study the merging cluster 1E 0657-56 or the ``Bullet Cluster'', in which a radio relic has not been found at the position of the prominent shock ahead of the bullet. We indicate that a possible relic at the shock is obscured by the observed large radio halo that is generated by strong turbulence behind the shock. We propose a simple explanation if the morphological differences of radio emission among the Toothbrush, the Bullet, and the Sausage (CIZA J2242.8+5301) Clusters.

The widest frequency radio relic spectra: observations from 150 MHz to 30 GHz

Radio relics are patches of diffuse synchrotron radio emission that trace shock waves. Relics are thought to form when intra-cluster medium electrons are accelerated by cluster merger induced shock waves through the diffusive shock acceleration mechanism. In this paper, we present observations spanning 150 MHz to 30 GHz of the `Sausage' and `Toothbrush' relics from the Giant Metrewave and Westerbork telescopes, the Karl G. Jansky Very Large Array, the Effelsberg telescope, the Arcminute Microkelvin Imager and Combined Array for Research in Millimeter-wave Astronomy. We detect both relics at 30 GHz, where the previous highest frequency detection was at 16 GHz. The integrated radio spectra of both sources clearly steepen above 2 GHz, at the >6$\sigma$ significance level, supports the spectral steepening previously found in the `Sausage' and the Abell 2256 relic. Our results challenge the widely adopted simple formation mechanism of radio relics and suggest more complicated models have to be developed that, for example, involve re-acceleration of aged seed electrons.

Cosmic-ray acceleration in young protostars

The main signature of the interaction between cosmic rays and molecular clouds is the high ionisation degree. This decreases towards the densest parts of a cloud, where star formation is expected, because of energy losses and magnetic effects. However recent observations hint to high levels of ionisation in protostellar systems, therefore leading to an apparent contradiction that could be explained by the presence of energetic particles accelerated within young protostars. Our modelling consists of a set of conditions that has to be satisfied in order to have an efficient particle acceleration through the diffusive shock acceleration mechanism. We find that jet shocks can be strong accelerators of protons which can be boosted up to relativistic energies. Another possibly efficient acceleration site is located at protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate protons. Our results demonstrate the possibility of accelerating particles during the early phase of a proto-Solar-like system and can be used as an argument to support available observations. The existence of an internal source of energetic particles can have a strong and unforeseen impact on the star and planet formation process as well as on the formation of pre-biotic molecules.

Cosmic-ray acceleration at collisionless astrophysical shocks using Monte-Carlo simulations

Context. The di usive shock acceleration mechanism has been widely accepted as the acceleration mechanism for galactic cosmic rays. While self-consistent hybrid simulations have shown how power-law spectra are produced, detailed information on the interplay of di usive particle motion and the turbulent electromagnetic fields responsible for repeated shock crossings are still elusive. Aims. The framework of test-particle theory is applied to investigate the e ect of di usive shock acceleration by inspecting the obtained cosmic-ray energy spectra. The resulting energy spectra can be obtained this way from the particle motion and, depending on the prescribed turbulence model, the influence of stochastic acceleration through plasma waves can be studied. Methods. A numerical Monte-Carlo simulation code is extended to include collisionless shock waves. This allows one to trace the trajectories of test particle while they are being accelerated. In addition, the di usion coe cients can be obtained directly from the particle motion, which allows for a detailed understanding of the acceleration process. Results. The classic result of an energy spectrum with E 2 is only reproduced for parallel shocks, while, for all other cases, the energy spectral index is reduced depending on the shock obliqueness. Qualitatively, this can be explained in terms of the di usion coe cients in the directions that are parallel and perpendicular to the shock front.

Active Galactic Nuclei: Jets as the Source of Hadrons and Neutrinos

Active galactic nuclei are extragalactic sources, and their relativistic hot-plasma jets are believed to be the main candidates of the cosmic-ray origin, above the so-called knee region of the cosmic-ray spectrum. Relativistic shocks, either single or multiple, have been observed or been theorized to be forming within relativistic jet channels in almost all active galactic nuclei sources. The acceleration of non-thermal particles (e.g. electrons, protons) via the shock Fermi acceleration mechanism, is believed to be mainly responsible for the power-law energy distribution of the observed cosmic-rays, which in very high energies can consequently radiate high energy gamma-rays and neutrinos, through related radiation channels. Here, we will focus on the primary particle (hadronic) shock acceleration mechanism, and we will present a comparative simulation study of the properties of single and multiple relativistic shocks, which occur in AGN jets. We will show that the role of relativistic (quasi-parallel either quasi-perpendicular) shocks, is quite important since it can dramatically alter the primary CR spectral indices and acceleration eciencies. These properties being carried onto gamma-ray and neutrino radiation characteristics, makes the combination of them a quite appealing theme for relativistic plasma and shock acceleration physics, as well as observational cosmic-ray, gamma-ray and neutrino astronomy.

Major solar eruptions and high-energy particle events during solar cycle 24

We report on a study of all major solar eruptions that occurred on the front-side of the Sun during the rise to peak phase of cycle 24 (first 62 months) in order to understand the key factors affecting the occurrence of large solar energetic particle events (SEPs) and the ground levels enhancement (GLE) events. The eruptions involve major flares with soft X-ray peak flux >/= 5.0 x10-5 Wm-2 (i.e., flare size >/= M5.0) and accompanying coronal mass ejections (CMEs). The selection criterion was based on the fact that the only front-side GLE in cycle 24 (GLE 71) had a flare size of M5.1. Only ~37% of the major eruptions from the western hemisphere resulted in large SEP events. Almost the same number of large SEP events was produced in weaker eruptions (flare size <M5.0), suggesting that the soft X-ray flare is not a good indicator of SEP or GLE events. On the other hand, the CME speed is a better indicator of SEP and GLE events because it is consistently high supporting the shock acceleration mechanism for SEPs and GLEs. We found the CME speed, magnetic connectivity to Earth, and ambient conditions as the main factors that contribute to the lack of high energy particle events during cycle 24. Several eruptions poorly connected to Earth (eastern-hemisphere or behind-the-west-limb events) resulted in very large SEP events detected by the STEREO spacecraft. Some very fast CMEs, likely to have accelerated particles to GeV energies, did not result in a GLE event because of poor latitudinal connectivity. The stringent latitudinal requirement suggests that the highest energy particles are likely accelerated in the nose part of shocks. There were also well-connected fast CMEs, which did not seem to have accelerated high energy particles due to possible unfavorable ambient conditions (high Alfven speed, overall reduction in acceleration efficiency in cycle 24).