-
Numerical simulation of photospheric emission in long gamma-ray bursts: prompt correlations, spectral shapes, and polarizations
Authors:
Hirotaka Ito,
Jin Matsumoto,
Shigehiro Nagataki,
Donald C. Warren,
Maxim V. Barkov,
Daisuke Yonetoku
Abstract:
We explore the properties of photospheric emission in the context of long gamma-ray bursts (LGRBs) using three numerical models that combine relativistic hydrodynamical simulations and Monte Carlo radiation transfer calculations in three dimensions. Our simulations confirm that the photospheric emission gives rise to correlations between the spectral peak energy and luminosity that agree with the…
▽ More
We explore the properties of photospheric emission in the context of long gamma-ray bursts (LGRBs) using three numerical models that combine relativistic hydrodynamical simulations and Monte Carlo radiation transfer calculations in three dimensions. Our simulations confirm that the photospheric emission gives rise to correlations between the spectral peak energy and luminosity that agree with the observed Yonetoku, Amati, and Golenetskii correlations. It is also shown that the spectral peak energy and luminosity correlate with the bulk Lorentz factor, as indicated in the literature. On the other hand, synthetic spectral shapes tend to be narrower than those of the observations. The result indicates that an additional physical process that can provide non-thermal broadening is needed to reproduce the spectral features. Furthermore, the polarization analysis finds that, while the degree of polarization is low for the emission from the jet core ($Π< 4~\%$), it tends to increase with the viewing angle outside the core and can be as high as $Π\sim 20-40~\%$ in an extreme case. This suggests that the typical GRBs show systematically low polarization compared to softer, dimmer counterparts (X-ray-rich GRBs and X-ray flashes). Interestingly, our simulations indicate that photospheric emission exhibits large temporal variation in the polarization position angle ($Δψ\sim 90^{\circ}$), which may be compatible with those inferred in observations. A notable energy dependence of the polarization property is another characteristic feature found in the current study. Particularly, the difference in the position angle among different energy bands can be as large as $\sim 90^{\circ}$.
△ Less
Submitted 19 July, 2023;
originally announced July 2023.
-
Synchrotron Self-Compton Afterglow Closure Relations and Fermi-LAT Detected Gamma-Ray Bursts
Authors:
N. Fraija,
M. G. Dainotti,
S. Ugale,
D. Jyoti,
D. C. Warren
Abstract:
The Fermi Large Area Telescope (Fermi-LAT) Collaboration reported the Second Gamma-ray Burst Catalog (2FLGC), which comprises a subset of 29 bursts with photon energies above 10 GeV. Although the standard synchrotron forward-shock model has successfully explained the Gamma-ray burst (GRB) afterglow observations, energetic photons higher than 10 GeV from these transient events can hardly be describ…
▽ More
The Fermi Large Area Telescope (Fermi-LAT) Collaboration reported the Second Gamma-ray Burst Catalog (2FLGC), which comprises a subset of 29 bursts with photon energies above 10 GeV. Although the standard synchrotron forward-shock model has successfully explained the Gamma-ray burst (GRB) afterglow observations, energetic photons higher than 10 GeV from these transient events can hardly be described in this scenario. We present the closure relations (CRs) of synchrotron self-Compton (SSC) afterglow model in the adiabatic and radiative scenario and when the central engine injects continuous energy into the blastwave to study the evolution of the spectral and temporal indexes of those bursts reported in 2FLGC. We consider the SSC afterglow model evolving in stellar-wind and interstellar medium, and the CRs as a function of the radiative parameter, the energy injection index, and the electron spectral index for $1<p<2$ and $ 2\leq p$. We select all GRBs that have been modeled with both a simple or a broken power law in the 2FLGC. We found that the CRs of the SSC model can satisfy a significant fraction of burst that cannot be interpreted in the synchrotron scenario, even though those that require an intermediate density profile (e.g., GRB 130427A) or an atypical fraction of total energy given to amplify the magnetic field ($\varepsilon_B$). The value of this parameter in the SSC model ranges ($\varepsilon_B\approx 10^{-5} - 10^{-4}$) when the cooling spectral break corresponds to the Fermi-LAT band for typical values of GRB afterglow. The analysis shows that ISM is preferred for the scenario without energy injection and the stellar wind medium for an energy injection scenario.
△ Less
Submitted 23 June, 2022;
originally announced June 2022.
-
The double detonation of a double degenerate system, from Type Ia supernova explosion to its supernova remnant
Authors:
Gilles Ferrand,
Ataru Tanikawa,
Donald C. Warren,
Shigehiro Nagataki,
Samar Safi-Harb,
Anne Decourchelle
Abstract:
Type Ia supernovae (SNe) are believed to be caused by the thermonuclear explosion of a white dwarf (WD), but the nature of the progenitor system(s) is still unclear. Recent theoretical and observational developments have led to renewed interest in double degenerate models, in particular the "helium-ignited violent merger" or "dynamically-driven double-degenerate double-detonation" (D$^6$). In this…
▽ More
Type Ia supernovae (SNe) are believed to be caused by the thermonuclear explosion of a white dwarf (WD), but the nature of the progenitor system(s) is still unclear. Recent theoretical and observational developments have led to renewed interest in double degenerate models, in particular the "helium-ignited violent merger" or "dynamically-driven double-degenerate double-detonation" (D$^6$). In this paper we take the output of an existing D$^6$ SN model and carry it into the supernova remnant (SNR) phase up to 4000 years after the explosion, past the time when all the ejecta have been shocked. Assuming a uniform ambient medium, we reveal specific signatures of the explosion mechanism and spatial variations intrinsic to the ejecta. The first detonation produces an ejecta tail visible at early times, while the second detonation leaves a central density peak in the ejecta that is visible at late times. The SNR shell is off-centre at all times, because of an initial velocity shift due to binary motion. The companion WD produces a large conical shadow in the ejecta, visible in projection as a dark patch surrounded by a bright ring. This is a clear and long-lasting feature that is localized, and its impact on the observed morphology is dependent on the viewing angle of the SNR. These results offer a new way to diagnose the explosion mechanism and progenitor system using observations of a Type Ia SNR.
△ Less
Submitted 8 February, 2022;
originally announced February 2022.
-
Examining Two-Dimensional Luminosity-time Correlations for Gamma Ray Burst Radio Afterglows with VLA and ALMA
Authors:
Delina Levine,
Maria Dainotti,
Kevin J. Zvonarek,
Nissim Fraija,
Donald C. Warren,
Poonam Chandra,
Nicole Lloyd-Ronning
Abstract:
Gamma-ray burst (GRB) afterglow emission can be observed from sub-TeV to radio wavelengths, though only 6.6\% of observed GRBs present radio afterglows. We examine GRB radio light curves (LCs) to look for the presence of radio plateaus, resembling the plateaus observed in X-ray and optical. We analyze 404 GRBs from the literature with observed radio afterglow and fit 82 GRBs with at least 5 data p…
▽ More
Gamma-ray burst (GRB) afterglow emission can be observed from sub-TeV to radio wavelengths, though only 6.6\% of observed GRBs present radio afterglows. We examine GRB radio light curves (LCs) to look for the presence of radio plateaus, resembling the plateaus observed in X-ray and optical. We analyze 404 GRBs from the literature with observed radio afterglow and fit 82 GRBs with at least 5 data points with a broken power law (BPL) model, requiring 4 parameters. From these, we find 18 GRBs that present a break feature resembling a plateau. We conduct the first multi-wavelength study of the Dainotti correlation between the luminosity $L_a$ and the rest-frame time of break $T_a^*$ for those 18 GRBs, concluding that the correlation exists and resembles the corresponding correlation in X-ray and optical wavelengths after correction for evolutionary effects. We compare the $T_a^*$ for the radio sample with $T_a^*$ values in X-ray and optical data (Dainotti et al. 2013, 2020b), finding significantly later break times in radio. We propose that this late break time and compatibility in slope suggests either a long-lasting plateau or the passage of a spectral break in the radio band. We also correct the distribution of the isotropic energy $E_{\rm iso}$ vs. the rest-frame burst duration $T^*_{90}$ for evolutionary effects and conclude that there is no significant difference between the $T^*_{90}$ distribution for the radio LCs with a break and those without.
△ Less
Submitted 19 November, 2021;
originally announced November 2021.
-
A semi-analytic afterglow with thermal electrons and synchrotron self-Compton emission
Authors:
Donald C. Warren,
Maria Dainotti,
Maxim V. Barkov,
Bjorn Ahlgren,
Hirotaka Ito,
Shigehiro Nagataki
Abstract:
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 emission 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) e…
▽ More
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 emission 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.
△ Less
Submitted 15 September, 2021;
originally announced September 2021.
-
Closure relations during the plateau emission of Swift GRBs and the fundamental plane
Authors:
M. G. Dainotti,
A. Ł. Lenart,
N. Fraija,
S. Nagataki,
D. C. Warren,
B. De Simone,
G. Srinivasaragavan,
A. Mata
Abstract:
The Neil Gehrels Swift observatory observe Gamma-Ray bursts (GRBs) plateaus in X-rays. We test the reliability of the closure relations through the fireball model when dealing with the GRB plateau emission. We analyze 455 X-ray lightcurves (LCs) collected by \emph{Swift} from 2005 (January) until 2019 (August) for which the redshift is both known and unknown using the phenomenological Willingale 2…
▽ More
The Neil Gehrels Swift observatory observe Gamma-Ray bursts (GRBs) plateaus in X-rays. We test the reliability of the closure relations through the fireball model when dealing with the GRB plateau emission. We analyze 455 X-ray lightcurves (LCs) collected by \emph{Swift} from 2005 (January) until 2019 (August) for which the redshift is both known and unknown using the phenomenological Willingale 2007 model. Using these fits, we analyze the emission mechanisms and astrophysical environments of these GRBs through the closure relations within the time interval of the plateau emission. Finally, we test the 3D fundamental plane relation (Dainotti relation) which connects the prompt peak luminosity, the time at the end of the plateau (rest-frame), and the luminosity at that time, on the GRBs with redshift, concerning groups determined by the closure relations. This allows us to check if the intrinsic scatter σ_{int} of any of these groups is reduced compared to previous literature. The most fulfilled environments for the electron spectral distribution, p>2, are Wind Slow Cooling (SC) and ISM Slow Cooling for cases in which the parameter q, which indicates the flatness of the plateau emission and accounts for the energy injection, is =0 and =0.5, respectively, both in the cases with known and unknown redshifts.
We also find that for the sGRBs All ISM Environments with $q=0$ have the smallest σ_{int}=0.04 \pm 0.15 in terms of the fundamental plane relation holding a probability of occurring by chance of p=0.005. We have shown that the majority of GRBs presenting the plateau emission fulfil the closure relations, including the energy injection, with a particular preference for the Wind SC environment. The subsample of GRBs that fulfil given relations can be used as possible standard candles and can suggest a way to reduce the intrinsic scatter of these studied relationships.
△ Less
Submitted 22 May, 2021;
originally announced May 2021.
-
From supernova to supernova remnant: comparison of thermonuclear explosion models
Authors:
Gilles Ferrand,
Donald C. Warren,
Masaomi Ono,
Shigehiro Nagataki,
Friedrich K. Roepke,
Ivo R. Seitenzahl,
Florian Lach,
Hiroyoshi Iwasaki,
Toshiki Sato
Abstract:
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 mo…
▽ More
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.
△ Less
Submitted 9 November, 2020;
originally announced November 2020.
-
The maximum energy of shock-accelerated electrons in a microturbulent magnetic field
Authors:
Donald C. Warren,
Catherine A. A. Beauchemin,
Maxim V. Barkov,
Shigehiro Nagataki
Abstract:
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 particular that of electrons. In this work we perform Monte Carlo (MC) simulations of electrons encountering shocks with microturbulent fields. The simulati…
▽ More
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 particular 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.
△ Less
Submitted 30 October, 2020; v1 submitted 13 October, 2020;
originally announced October 2020.
-
3D modeling from the onset of the SN to the full-fledged SNR: Role of an initial ejecta anisotropy on matter mixing
Authors:
Antonio Tutone,
Salvatore Orlando,
Marco Miceli,
Sabina Ustamujic,
Masaomi Ono,
Shigehiro Nagataki,
Gilles Ferrand,
Emanuele Greco,
Giovanni Peres,
Donald C. Warren,
Fabrizio Bocchino
Abstract:
The aim of this work is to bridge the gap between CC SNe and their remnants by investigating how post-explosion anisotropies in the ejecta influence the structure and chemical properties of the remnant at later times. We performed three-dimensional magneto-hydrodynamical simulations starting soon after the SN event and following the evolution of the system in the circumstellar medium (consisting o…
▽ More
The aim of this work is to bridge the gap between CC SNe and their remnants by investigating how post-explosion anisotropies in the ejecta influence the structure and chemical properties of the remnant at later times. We performed three-dimensional magneto-hydrodynamical simulations starting soon after the SN event and following the evolution of the system in the circumstellar medium (consisting of the wind of the stellar progenitor), for 5000 years, obtaining the physical scenario of a SNR. Here we focused the analysis on the case of a progenitor red supergiant of 19.8 M_sun. We also investigated how a post-explosion large-scale anisotropy in the SN affects the ejecta distribution and the matter mixing of heavy elements in the remnant, during the first 5000 years of evolution. In the case of a spherically symmetric SN explosion without large-scale anisotropies, the remnant roughly keeps memory of the original onion-like layering of ejecta soon after the SN event. Nevertheless, as the reverse shock hits the ejecta, the element distribution departs from a homologous expansion, because of the slowing down of the outermost ejecta layers due to interaction with the reverse shock. In the case of a large-scale anisotropy developed after the SN, we found that the chemical stratification in the ejecta can be strongly modified and the original onion-like layering is not preserved. The anisotropy may cause spatial inversion of ejecta layers, for instance leading to Fe/Si-rich ejecta outside the O shell, and may determine the formation of Fe/Si-rich jet-like features that may protrude the remnant outline. The level of matter mixing and the properties of the jet-like feature are sensitive to the initial physical (density and velocity) and geometrical (size and position) initial characteristics of the anisotropy.
△ Less
Submitted 2 September, 2020;
originally announced September 2020.
-
From the supernova to the supernova remnant: the three-dimensional imprint of a thermonuclear explosion
Authors:
Gilles Ferrand,
Donald C. Warren,
Masaomi Ono,
Shigehiro Nagataki,
Friedrich K. Roepke,
Ivo R. Seitenzahl
Abstract:
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…
▽ More
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 Tycho's 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.
△ Less
Submitted 16 April, 2019;
originally announced April 2019.
-
The photospheric origin of the Yonetoku relation in gamma-ray bursts
Authors:
Hirotaka Ito,
Jin Matsumoto,
Shigehiro Nagataki,
Donald C. Warren,
Maxim V. Barkov,
Daisuke Yonetoku
Abstract:
Long duration gamma-ray bursts (GRBs), the brightest events since the Big Bang itself, are believed to originate in an ultra-relativistic jet breaking out from a massive stellar envelope. Despite decades of study, there is still no consensus on their emission mechanism. One unresolved question is the origin of the tight correlation between the spectral peak energy Ep and peak luminosity Lp discove…
▽ More
Long duration gamma-ray bursts (GRBs), the brightest events since the Big Bang itself, are believed to originate in an ultra-relativistic jet breaking out from a massive stellar envelope. Despite decades of study, there is still no consensus on their emission mechanism. One unresolved question is the origin of the tight correlation between the spectral peak energy Ep and peak luminosity Lp discovered in observations. This Yonetoku relation is the tightest correlation found in the properties of the prompt phase of GRB emission, providing the best diagnostic for the radiation mechanism. Here we present 3D hydrodynamical simulations, and post-process radiation transfer calculations, of photospheric emission from a relativistic jet. Our simulations reproduce the Yonetoku relation as a natural consequence of viewing angle. Although jet dynamics depend sensitively on luminosity, the Ep-Lp correlation holds regardless. This result strongly suggests that photospheric emission is the dominant component in the prompt phase of GRBs.
△ Less
Submitted 3 April, 2019; v1 submitted 2 June, 2018;
originally announced June 2018.
-
Synchrotron self-absorption in GRB afterglows: the effects of a thermal electron population
Authors:
Donald C. Warren,
Maxim V. Barkov,
Hirotaka Ito,
Shigehiro Nagataki,
Tanmoy Laskar
Abstract:
In the standard synchrotron afterglow model, a power law of electrons is responsible for all aspects of photon production and absorption. Recent numerical work has shown that the vast majority of particles in the downstream medium are actually "thermal" particles, which were shock-heated but did not enter the Fermi acceleration process (the name stands in contrast to the nonthermal high-energy tai…
▽ More
In the standard synchrotron afterglow model, a power law of electrons is responsible for all aspects of photon production and absorption. Recent numerical work has shown that the vast majority of particles in the downstream medium are actually "thermal" particles, which were shock-heated but did not enter the Fermi acceleration process (the name stands in contrast to the nonthermal high-energy tail, rather than connoting a Maxwellian distribution). There are substantial differences at optical and higher energies when these thermal electrons participate in the afterglow, but early work along these lines ignored the radio end of the electromagnetic spectrum. We report here on an extension of previous Monte Carlo simulations of gamma-ray burst afterglows. The model now includes the synchrotron self-absorption (SSA) process and so can simulate afterglows across the entire EM spectrum, and several orders of magnitude in time. In keeping with earlier work, inclusion of the thermal electrons increases the SSA frequency by a factor of 30, and the radio intensity by a factor of 100. Furthermore, these changes happen with no modification to the late optical or X-ray afterglow. Our results provide very strong evidence that thermal electrons must be considered in any multiwavelength model for afterglows.
△ Less
Submitted 16 April, 2018;
originally announced April 2018.
-
Nonlinear particle acceleration and thermal particles in GRB afterglows
Authors:
Donald C. Warren,
Donald C. Ellison,
Maxim V. Barkov,
Shigehiro Nagataki
Abstract:
The standard model for GRB afterglow emission treats the accelerated electron population as a simple power law, $N(E) \propto E^{-p}$ for $p \gtrsim 2$. 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 backreaction of accelerated particles on the…
▽ More
The standard model for GRB afterglow emission treats the accelerated electron population as a simple power law, $N(E) \propto E^{-p}$ for $p \gtrsim 2$. 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 backreaction of accelerated particles on the shock structure modifies the shape of the non-thermal 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 (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.
△ Less
Submitted 16 January, 2017;
originally announced January 2017.
-
Particle spectra and efficiency in nonlinear relativistic shock acceleration: survey of scattering models
Authors:
Donald C. Ellison,
Donald C. Warren,
Andrei M. Bykov
Abstract:
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 to produce short-scale, self-generated magnetic turbulence that leads to a scattering mean free path (mfp) with a stronger momentum dependence than the…
▽ More
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 to 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.
△ Less
Submitted 9 December, 2015;
originally announced December 2015.
-
Photospheric Emission from Collapsar Jets in 3D Relativistic Hydrodynamics
Authors:
Hirotaka Ito,
Jin Matsumoto,
Shigehiro Nagataki,
Donald C. Warren,
Maxim V. Barkov
Abstract:
We explore the photospheric emission from a relativistic jet breaking out from a massive stellar envelope based on relativistic hydrodynamical simulations and post-process radiation transfer calculations in three dimensions. To investigate the impact of three-dimensional (3D) dynamics on the emission, two models of injection conditions are considered for the jet at the center of the progenitor sta…
▽ More
We explore the photospheric emission from a relativistic jet breaking out from a massive stellar envelope based on relativistic hydrodynamical simulations and post-process radiation transfer calculations in three dimensions. To investigate the impact of three-dimensional (3D) dynamics on the emission, two models of injection conditions are considered for the jet at the center of the progenitor star: one with periodic precession and another without precession. We show that structures developed within the jet due to the interaction with the stellar envelope, as well as due to the precession, have a significant imprint on the resulting emission. Particularly, we find that the signature of precession activity by the central engine is not smeared out and can be directly observed in the light curve as a periodic signal. We also show non-thermal features, which can account for observations of gamma-ray bursts, are produced in the resulting spectra, even though only thermal photons are injected initially and the effect of non-thermal particles is not considered.
△ Less
Submitted 27 November, 2015; v1 submitted 11 November, 2015;
originally announced November 2015.
-
Electron and Ion Acceleration in Relativistic Shocks with Applications to GRB Afterglows
Authors:
Donald C. Warren,
Donald C. Ellison,
Andrei M. Bykov,
Shiu-Hang Lee
Abstract:
We have modeled the simultaneous first-order Fermi shock acceleration of protons, electrons, and helium nuclei by relativistic shocks. By parameterizing the particle diffusion, our steady-state Monte Carlo simulation allows us to follow particles from particle injection at nonthermal thermal energies to above PeV energies, including the nonlinear smoothing of the shock structure due to cosmic-ray…
▽ More
We have modeled the simultaneous first-order Fermi shock acceleration of protons, electrons, and helium nuclei by relativistic shocks. By parameterizing the particle diffusion, our steady-state Monte Carlo simulation allows us to follow particles from particle injection at nonthermal thermal energies to above PeV energies, including the nonlinear 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 parameterize 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 gamma-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 nonlinear 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 nonlinear effects should be qualitatively similar as long as the scattering mean free path is an increasing function of momentum.
△ Less
Submitted 9 June, 2015;
originally announced June 2015.
-
Monte Carlo Simulations of Nonlinear Particle Acceleration in Parallel Trans-relativistic Shocks
Authors:
Donald C. Ellison,
Donald C. Warren,
Andrei M. Bykov
Abstract:
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 tha…
▽ More
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 = 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.
△ Less
Submitted 23 August, 2013;
originally announced August 2013.
-
Three-dimensional numerical investigations of the morphology of type Ia SNRs
Authors:
D. C. Warren,
J. M. Blondin
Abstract:
We explore the morphology of Type Ia supernova remnants (SNRs) using three-dimensional hydrodynamics modeling and an exponential density profile. Our model distinguishes ejecta from the interstellar medium (ISM), and tracks the ionization age of shocked ejecta, both of which allow for additional analysis of the simulated remnants. We also include the adiabatic index as a free parameter, which affe…
▽ More
We explore the morphology of Type Ia supernova remnants (SNRs) using three-dimensional hydrodynamics modeling and an exponential density profile. Our model distinguishes ejecta from the interstellar medium (ISM), and tracks the ionization age of shocked ejecta, both of which allow for additional analysis of the simulated remnants. We also include the adiabatic index as a free parameter, which affects the compressibility of the fluid and emulates the efficiency of cosmic ray acceleration by shock fronts. In addition to generating 3-D images of the simulations, we compute line-of-sight projections through the remnants for comparison against observations of Tycho's SNR and SN 1006. We find that several features observed in these two remnants, such as the separation between the fluid discontinuities and the presence of ejecta knots ahead of the forward shock, can be generated by smooth ejecta without any initial clumpiness. Our results are consistent with SN 1006 being dynamically younger than Tycho's SNR, and more efficiently accelerating cosmic rays at its forward shock. We conclude that clumpiness is not a necessary condition to reproduce many observed features of Type Ia supernova remnants, particularly the radial profiles and the fleecy emission from ejecta at the central region of both remnants.
△ Less
Submitted 29 October, 2012;
originally announced October 2012.