We extend previous work on gamma-ray burst (GRB) afterglows involving hot thermal electrons at the base of a shock-accelerated tail. Using a physically-motivated electron distribution based on first-principles simulations, we compute broadband emissi
on from radio to TeV gamma-rays. For the first time, we present the effects of a thermal distribution of electrons on synchrotron self-Compton (SSC) emission. The presence of thermal electrons causes temporal and spectral structure across the entire observable afterglow, which is substantively different from models that assume a pure power-law distribution for the electrons. We show that early-time TeV emission is enhanced by more than an order of magnitude for our fiducial parameters, with a time-varying spectral index that does not occur for a pure power law of electrons. We further show that the X-ray closure relations take a very different, also time-dependent, form when thermal electrons are present; the shape traced out by the X-ray afterglows is a qualitative match to observations of the traditional decay phase.
Progress in the three-dimensional modeling of supernovae (SN) prompts us to revisit the supernova remnant (SNR) phase. We continue our study of the imprint of a thermonuclear explosion on the SNR it produces, that we started with a delayed-detonation
model of a Chandrasekhar-mass white dwarf. Here we compare two different types of explosion models, each with two variants: two delayed detonation models (N100ddt, N5ddt) and two pure deflagration models (N100def, N5def), where the N number parametrizes the ignition. The output of each SN simulation is used as input of a SNR simulation carried on until 500 yr after the explosion. While all SNR models become more spherical over time and overall display the theoretical structure expected for a young SNR, clear differences are visible amongst the models, depending on the geometry of the ignition and on the presence or not of detonation fronts. Compared to N100 models, N5 models have a strong dipole component, and produce asymmetric remnants. N5def produces a regular-looking, but offset remnant, while N5ddt produces a two-sided remnant. Pure deflagration models exhibit specific traits: a central over-density, because of the incomplete explosion, and a network of seam lines across the surface, boundaries between burning cells. Signatures from the SN dominate the morphology of the SNR up to 100 yr to 300 yr after the explosion, depending on the model, and are still measurable at 500 yr, which may provide a way of testing explosion models.
Relativistic shocks propagating into a medium with low magnetization are generated and sustained by small-scale but very strong magnetic field turbulence. This so-called microturbulence modifies the typical shock acceleration process, and in particul
ar that of electrons. In this work we perform Monte Carlo (MC) simulations of electrons encountering shocks with microturbulent fields. The simulations cover a three-dimensional parameter space in shock speed, acceleration efficiency, and peak magnetic field strength. From these, a Markov Chain Monte Carlo (MCMC) method was employed to estimate the maximum electron momentum from the MC-simulated electron spectra. Having estimated this quantity at many points well-distributed over an astrophysically relevant parameter space, an MCMC method was again used to estimate the parameters of an empirical formula that computes the maximum momentum of a Fermi-accelerated electron population anywhere in this parameter space. The maximum energy is well-approximated as a broken power-law in shock speed, with the break occurring when the shock decelerates to the point where electrons can begin to escape upstream from the shock.
Recent progress in the three-dimensional modeling of supernovae (SN) has shown the importance of asymmetries for the explosion. This calls for a reconsideration of the modeling of the subsequent phase, the supernova remnant (SNR), which has commonly
relied on simplified ejecta models. In this paper we bridge SN and SNR studies by using the output of a SN simulation as the input of a SNR simulation carried on until 500~yr. We consider the case of a thermonuclear explosion of a carbon-oxygen white dwarf star as a model for a Type Ia SN; specifically we use the N100 delayed detonation model of Seitenzahl et al 2013. In order to analyze the morphology of the SNR, we locate the three discontinuities that delineate the shell of shocked matter: the forward shock, the contact discontinuity, and the reverse shock, and we decompose their radial variations as a function of angular scale and time. Assuming a uniform ambient medium, we find that the impact of the SN on the SNR may still be visible after hundreds of years. Previous 3D simulations aiming at reproducing Tychos SNR, that started out from spherically symmetric initial conditions, failed to reproduce structures at the largest angular scales observed in X-rays. Our new simulations strongly suggest that the missing ingredient was the initial asymmetries from the SN itself. With this work we establish a way of assessing the viability of SN models based on the resulting morphology of the SNR.
We include a general form for the scattering mean free path in a nonlinear Monte Carlo model of relativistic shock formation and Fermi acceleration. Particle-in-cell (PIC) simulations, as well as analytic work, suggest that relativistic shocks tend t
o produce short-scale, self-generated magnetic turbulence that leads to a scattering mean free path (mfp) 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 gamma-ray bursts (GRBs), pulsar winds, Type Ibc supernovae, and extra-galactic 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 the mfp has an important influence on the efficiency of cosmic-ray production as well as the accelerated particle spectral shape. These effects are absent in nonrelativistic 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 gamma-ray emission from proton-proton interactions and pion-decay radiation.
