The standard model for GRB afterglow emission treats the accelerated electron population as a simple power law, for . However, in standard Fermi shock acceleration, a substantial fraction of the ...swept-up particles do not enter the acceleration process at all. Additionally, if acceleration is efficient, then the nonlinear back-reaction of accelerated particles on the shock structure modifies the shape of the nonthermal tail of the particle spectra. Both of these modifications to the standard synchrotron afterglow impact the luminosity, spectra, and temporal variation of the afterglow. To examine the effects of including thermal particles and nonlinear particle acceleration on afterglow emission, we follow a hydrodynamical model for an afterglow jet and simulate acceleration at numerous points during the evolution. When thermal particles are included, we find that the electron population is at no time well fitted by a single power law, though the highest-energy electrons are; if the acceleration is efficient, then the power-law region is even smaller. Our model predicts hard-soft-hard spectral evolution at X-ray energies, as well as an uncoupled X-ray and optical light curve. Additionally, we show that including emission from thermal particles has drastic effects (increases by factors of 100 and 30, respectively) on the observed flux at optical and GeV energies. This enhancement of GeV emission makes afterglow detections by future γ-ray observatories, such as CTA, very likely.
We include a general form for the scattering mean free path, λmfp(p), in a nonlinear Monte Carlo model of relativistic shock formation and Fermi acceleration. Particle-in-cell simulations, as well as ...analytic work, suggest that relativistic shocks tend to produce short-scale, self-generated magnetic turbulence that leads to a scattering mean free path with a stronger momentum dependence than the λmfp ∝ p dependence for Bohm diffusion. In unmagnetized shocks, this turbulence is strong enough to dominate the background magnetic field so the shock can be treated as parallel regardless of the initial magnetic field orientation, making application to γ-ray bursts, pulsar winds, type Ibc supernovae, and extragalactic radio sources more straightforward and realistic. In addition to changing the scale of the shock precursor, we show that, when nonlinear effects from efficient Fermi acceleration are taken into account, the momentum dependence of λmfp(p) has an important influence on the efficiency of cosmic ray production as well as the accelerated particle spectral shape. These effects are absent in non-relativistic shocks and do not appear in relativistic shock models unless nonlinear effects are self-consistently described. We show, for limited examples, how the changes in Fermi acceleration translate to changes in the intensity and spectral shape of γ-ray emission from proton–proton interactions and pion-decay radiation.
We present results from a Monte Carlo simulation of a parallel collisionless shock undergoing particle acceleration. Our simulation, which contains parameterized scattering and a particular thermal ...leakage injection model, calculates the feedback between accelerated particles ahead of the shock, which influence the shock precursor and "smooth" the shock, and thermal particle injection. We show that there is a transition between nonrelativistic shocks, where the acceleration efficiency can be extremely high and the nonlinear compression ratio can be substantially greater than the Rankine-Hugoniot value, and fully relativistic shocks, where diffusive shock acceleration is less efficient and the compression ratio remains at the Rankine-Hugoniot value. This transition occurs in the trans-relativistic regime and, for the particular parameters we use, occurs around a shock Lorentz factor gamma sub(0) = 1.5. We also find that nonlinear shock smoothing dramatically reduces the acceleration efficiency presumed to occur with large-angle scattering in ultra-relativistic shocks. Our ability to seamlessly treat the transition from ultra-relativistic to trans-relativistic to nonrelativistic shocks may be important for evolving relativistic systems, such as gamma-ray bursts and Type Ibc supernovae. We expect a substantial evolution of shock accelerated spectra during this transition from soft early on to much harder when the blast-wave shock becomes nonrelativistic.
We model the broadband emission from supernova remnant (SNR) RX J1713.7-3946 including, for the first time, a consistent calculation of thermal X-ray emission together with non-thermal emission in a ...nonlinear diffusive shock acceleration model. Our model tracks the evolution of the SNR including the plasma ionization state between the forward shock and the contact discontinuity. We use a plasma emissivity code to predict the thermal X-ray emission spectrum assuming the initially cold electrons are heated either by Coulomb collisions with the shock-heated protons (the slowest possible heating), or come into instant equilibration with the protons. For either electron heating model, electrons reach {approx}>10{sup 7} K rapidly and the X-ray line emission near 1 keV is more than 10 times as luminous as the underlying thermal bremsstrahlung continuum. Since recent Suzaku observations show no detectable line emission, this places strong constraints on the unshocked ambient medium density and on the relativistic electron-to-proton ratio. For the uniform circumstellar medium (CSM) models that we consider, the low densities and high relativistic electron-to-proton ratios required to match the Suzaku X-ray observations definitively rule out pion decay as the emission process producing GeV-TeV photons. We show that leptonic models, where inverse-Compton scattering against the cosmic background radiation dominates the GeV-TeV emission, produce better fits to the broadband thermal and non-thermal observations in a uniform CSM.
The giant molecular clouds (MCs) found in the Milky Way and similar galaxies play a crucial role in the evolution of these systems. The supernova explosions that mark the death of massive stars in ...these regions often lead to interactions between the supernova remnants (SNRs) and the clouds. These interactions have a profound effect on our understanding of SNRs. Shocks in SNRs should be capable of accelerating particles to cosmic ray (CR) energies with efficiencies high enough to power Galactic CRs. X-ray and
γ
-ray studies have established the presence of relativistic electrons and protons in some SNRs and provided strong evidence for diffusive shock acceleration as the primary acceleration mechanism, including strongly amplified magnetic fields, temperature and ionization effects on the shock-heated plasmas, and modifications to the dynamical evolution of some systems. Because protons dominate the overall energetics of the CRs, it is crucial to understand this hadronic component even though electrons are much more efficient radiators and it can be difficult to identify the hadronic component. However, near MCs the densities are sufficiently high to allow the
γ
-ray emission to be dominated by protons. Thus, these interaction sites provide some of our best opportunities to constrain the overall energetics of these particle accelerators. Here we summarize some key properties of interactions between SNRs and MCs, with an emphasis on recent X-ray and
γ
-ray studies that are providing important constraints on our understanding of cosmic rays in our Galaxy.
Observations of young supernova remnants (SNRs) in X-rays and γ-rays have provided conclusive evidence for particle acceleration to at least TeV energies. Analysis of high-spatial-resolution X-ray ...maps of young SNRs has indicated that the particle acceleration process is accompanied by strong nonadiabatic amplification of magnetic fields. If Fermi acceleration is the mechanism producing the energetic cosmic rays (CRs), the amplified magnetic field must be turbulent, and CR-driven instabilities are among the most probable mechanisms for converting the shock ram pressure into magnetic turbulence. The development and evolution of strong magnetic turbulence in collisionless plasmas forming SNR shells are complicated phenomena which include the amplification of magnetic modes, anisotropic mode transformations at shocks, as well as the nonlinear physics of turbulent cascades. Polarized X-ray synchrotron radiation from ultrarelativistic electrons accelerated in the SNR shock is produced in a thin layer immediately behind the shock and is not subject to the Faraday depolarization effect. These factors open up possibilities to study some properties of magnetic turbulence, and here we present polarized X-ray synchrotron maps of SNR shells assuming different models of magnetic turbulence cascades. It is shown that different models of anisotropic turbulence can be distinguished by measuring the predominant polarization angle direction. We discuss the detection of these features in Tycho's SNR with the coming generation of X-ray polarimeters such as the Imaging X-ray Polarimetry Explorer.
We have modelled the simultaneous first-order Fermi shock acceleration of protons, electrons, and helium nuclei by relativistic shocks. By parametrizing the particle diffusion, our steady-state Monte ...Carlo simulation allows us to follow particles from particle injection at non-relativistic thermal energies to above PeV energies, including the non-linear smoothing of the shock structure due to cosmic ray (CR) backpressure. We observe the mass-to-charge (A/Z) enhancement effect believed to occur in efficient Fermi acceleration in non-relativistic shocks and we parametrize the transfer of ion energy to electrons seen in particle-in-cell (PIC) simulations. For a given set of environmental and model parameters, the Monte Carlo simulation determines the absolute normalization of the particle distributions and the resulting synchrotron, inverse Compton, and pion-decay emission in a largely self-consistent manner. The simulation is flexible and can be readily used with a wide range of parameters typical of γ-ray burst (GRB) afterglows. We describe some preliminary results for photon emission from shocks of different Lorentz factors and outline how the Monte Carlo simulation can be generalized and coupled to hydrodynamic simulations of GRB blast waves. We assume Bohm diffusion for simplicity but emphasize that the non-linear effects we describe stem mainly from an extended shock precursor where higher energy particles diffuse further upstream. Quantitative differences will occur with different diffusion models, particularly for the maximum CR energy and photon emission, but these non-linear effects should be qualitatively similar as long as the scattering mean-free path is an increasing function of momentum.
We present a nonlinear Monte Carlo model of efficient diffusive shock acceleration where the magnetic turbulence responsible for particle diffusion is calculated self-consistently from the resonant ...cosmic-ray (CR) streaming instability, together with non-resonant short- and long-wavelength CR-current-driven instabilities. We include the backpressure from CRs interacting with the strongly amplified magnetic turbulence which decelerates and heats the super-Alfvenic flow in the extended shock precursor. Uniquely, in our plane-parallel, steady-state, multi-scale model, the full range of particles, from thermal (~eV) injected at the viscous subshock to the escape of the highest energy CRs (~PeV) from the shock precursor, are calculated consistently with the shock structure, precursor heating, magnetic field amplification, and scattering center drift relative to the background plasma. In addition, we show how the cascade of turbulence to shorter wavelengths influences the total shock compression, the downstream proton temperature, the magnetic fluctuation spectra, and accelerated particle spectra. A parameter survey is included where we vary shock parameters, the mode of magnetic turbulence generation, and turbulence cascading. From our survey results, we obtain scaling relations for the maximum particle momentum and amplified magnetic field as functions of shock speed, ambient density, and shock size.
To better model the efficient production of cosmic rays (CRs) in supernova remnants (SNRs) with the associated coupling between CR production and SNR dynamics, we have generalized an existing ...cr-hydro-NEI code to include the following processes: (1) an explicit calculation of the upstream precursor structure including the position-dependent flow speed, density, temperature, and magnetic field strength; (2) a momentum- and space-dependent CR diffusion coefficient; (3) an explicit calculation of magnetic field amplification; (4) calculation of the maximum CR momentum using the amplified magnetic field; (5) a finite Alfven speed for the particle scattering centers; and (6) the ability to accelerate a superthermal seed population of CRs, as well as the ambient thermal plasma. While a great deal of work has been done modeling SNRs, most work has concentrated on either the continuum emission from relativistic electrons or ions or the thermal emission from the shock heated plasma. Our generalized code combines these elements and describes the interplay between CR production and SNR evolution, including the nonlinear coupling of efficient diffusive shock acceleration, based mainly on the work of R Blasi and coworkers, and a non-equilibrium ionization (NEI) calculation of thermal X-ray line emission. We believe that our generalized model will provide a consistent modeling platform for SNRs, including those interacting with molecular clouds, and improve the interpretation of current and future observations, including the high-quality spectra expected from Astro-H. SNR RX J 1713.7-3946 is modeled as an example.
There is now substantial evidence that the progenitors of some core-collapse supernovae undergo enhanced or extreme mass loss prior to explosion. The imprint of this mass loss is observed in the ...spectra and dynamics of the expanding blast wave on timescales of days to years after core collapse, and the effects on the spectral and dynamical evolution may linger long after the supernova has evolved into the remnant stage. In this paper, we present, for the first time, largely self-consistent end-to-end simulations for the evolution of a massive star from the pre-main sequence, up to and through core collapse, and into the remnant phase. We present three models and compare and contrast how the progenitor mass-loss history impacts the dynamics and spectral evolution of the supernovae and supernova remnants. We study a model that only includes steady mass loss, a model with enhanced mass loss over a period of ∼5000 yr prior to core collapse, and a model with extreme mass loss over a period of ∼500 yr prior to core collapse. The models are not meant to address any particular supernova or supernova remnant, but rather to highlight the important role that the progenitor evolution plays in the observable qualities of supernovae and supernova remnants. Through comparisons of these three different progenitor evolution scenarios, we find that the mass loss in late stages (during and after core carbon burning) can have a profound impact on the dynamics and spectral evolution of the supernova remnant centuries after core collapse.