Efficiency of internal shocks in magnetized relativistic jets

We study the dynamic and radiative efficiency of conversion of kinetic-to-thermal/magnetic energy by internal shocks in relativistic magnetized outflows. A parameter study of a large number of collisions of cylindrical shells is performed. We explore how, while keeping the total flow luminosity constant, the variable fluid magnetization influences the efficiency and find that the interaction of shells in a mildly magnetized jet yields higher dynamic, but lower radiative efficiency than in a non-magnetized flow. A multi-wavelength radiative signature of different shell magnetization is computed assuming that relativistic particles are accelerated at internal shocks.Comment: 4 pages, 2 figures, proceedings of the meeting "HEPRO III: High Energy Phenomena in Relativistic Outflows" (Barcelona, June 2011), fixed the bibliography error

The radio afterglow of Swift J1644+57 reveals a powerful jet with fast core and slow sheath

We model the non-thermal transient Swift J1644+57 as resulting from a relativistic jet powered by the accretion of a tidally-disrupted star onto a super-massive black hole. Accompanying synchrotron radio emission is produced by the shock interaction between the jet and the dense circumnuclear medium, similar to a gamma-ray burst afterglow. An open mystery, however, is the origin of the late-time radio rebrightening, which occurred well after the peak of the jetted X-ray emission. Here, we systematically explore several proposed explanations for this behavior by means of multi-dimensional hydrodynamic simulations coupled to a self-consistent radiative transfer calculation of the synchrotron emission. Our main conclusion is that the radio afterglow of Swift J1644+57 is not naturally explained by a jet with a one-dimensional top-hat angular structure. However, a more complex angular structure comprised of an ultra-relativistic core (Lorentz factor $\Gamma \sim 10$) surrounded by a slower ($\Gamma \sim$ 2) sheath provides a reasonable fit to the data. Such a geometry could result from the radial structure of the super-Eddington accretion flow or as the result of jet precession. The total kinetic energy of the ejecta that we infer of $\sim$ few $10^{53}\,$erg requires a highly efficient jet launching mechanism. Our jet model providing the best fit to the light curve of the on-axis event Swift J1644+57 is used to predict the radio light curves for off-axis viewing angles. Implications for the presence of relativistic jets from TDEs detected via their thermal disk emission, as well as the prospects for detecting orphan TDE afterglows with upcoming wide-field radio surveys and resolving the jet structure with long baseline interferometry, are discussed.Comment: Accepted for publication in MNRA

A method for computing synchrotron and inverse-Compton emission from hydrodynamic simulations of supernova remnants

The observational signature of supernova remnants (SNRs) is very complex, in terms of both their geometrical shape and their spectral properties, dominated by non-thermal synchrotron and inverse-Compton scattering. We propose a post-processing method to analyse the broad-band emission of SNRs based on three-dimensional hydrodynamical simulations. From the hydrodynamical data, we estimate the distribution of non-thermal electrons accelerated at the shock wave and follow the subsequent evolution as they lose or gain energy by adiabatic expansion or compression and emit energy by radiation. As a first test case, we use a simulation of a bipolar supernova expanding into a cloudy medium. We find that our method qualitatively reproduces the main observational features of typical SNRs and produces fluxes that agree with observations to within a factor of a few. allowing for further use in more extended sets of models.Comment: 15 pages, 3 figures; accepted, HEDLA 2014 special issue of High Energy Density Physic

Multiwavelength afterglow light curves from magnetized GRB flows

We use high-resolution relativistic MHD simulations coupled with a radiative transfer code to compute multiwavelength afterglow light curves of magnetized ejecta of gamma-ray bursts interacting with a uniform circumburst medium. The aim of our study is to determine how the magnetization of the ejecta at large distance from the central engine influences the afterglow emission, and to assess whether observations can be reliably used to infer the strength of the magnetic field. We find that, for typical parameters of the ejecta, the emission from the reverse shock peaks for magnetization $\sigma_0 \sim 0.01 - 0.1$ of the flow, and that it is greatly suppressed for higher $\sigma_0$. The emission from the forward shock shows an achromatic break shortly after the end of the burst marking the onset of the self-similar evolution of the blast wave. Fitting the early afterglow of GRB 990123 and 090102 with our numerical models we infer respective magnetizations of $\sigma_0 \sim 0.01$ and $\sigma_0 \sim 0.1$ for these bursts. We argue that the lack of observed reverse shock emission from the majority of the bursts can be understood if \sigma_0 \simmore 0.1, since we obtain that the luminosity of the reverse shock decreases significantly for $\sigma_0 \sim 1$. For ejecta with \sigma_0 \simmore 0.1 our models predict that there is sufficient energy left in the magnetic field, at least during an interval of ~10 times the burst duration, to produce a substantial emission if the magnetic energy can be dissipated (for instance, due to resistive effects) and radiated away.Comment: 9 pages, 9 figures. Submitted to MNRAS

Internal shocks in relativistic outflows: collisions of magnetized shells

(Abridged): We study the collision of magnetized irregularities (shells) in relativistic outflows in order to explain the origin of the generic phenomenology observed in the non-thermal emission of both blazars and gamma-ray bursts. We focus on the influence of the magnetic field on the collision dynamics, and we further investigate how the properties of the observed radiation depend on the strength of the initial magnetic field and on the initial internal energy density of the flow. The collisions of magnetized shells and the radiation resulting from these collisions are calculated using the 1D relativistic magnetohydrodynamics code MRGENESIS. The interaction of the shells with the external medium prior to their collision is also determined using an exact solver for the corresponding 1D relativistic magnetohydrodynamic Riemann problem. Our simulations show that two magnetization parameters - the ratio of magnetic energy density and thermal energy density, \alpha_B, and the ratio of magnetic energy density and mass-energy density, \sigma - play an important role in the pre-collision phase, while the dynamics of the collision and the properties of the light curves depend mostly on the magnetization parameter \sigma. The interaction of the shells with the external medium changes the flow properties at their edges prior to the collision. For sufficiently dense shells moving at large Lorentz factors (\simgt 25) these properties depend only on the magnetization parameter \sigma. Internal shocks in GRBs may reach maximum efficiencies of conversion of kinetic into thermal energy between 6% and 10%, while in case of blazars, the maximum efficiencies are \sim 2%.Comment: 17 pages, 18 figures. 2 new references have been added. Accepted for publication in Astronomy and Astrophysic

Simulations of the Magneto-rotational Instability in Core-Collapse Supernovae

We assess the importance of the magneto-rotational instability in core-collapse supernovae by an analysis of the growth rates of unstable modes in typical post-collapse systems and by numerical simulations of simplified models. The interplay of differential rotation and thermal stratification defines different instability regimes which we confirm in our simulations. We investigate the termination of the growth of the MRI by parasitic instabilities, establish scaling laws characterising the termination amplitude, and study the long-term evolution of the saturated turbulent state.Comment: 6 pages, 1 figure. To appear in Proceedings of 4th International Conference on Numerical Modeling of Space Plasma Flows (Chamonix 2009

General Relativistic Simulations of the Collapsar Scenario

We are exploring the viability of the collapsar model for long-soft gamma-ray bursts. For this we perform state-of-the-art general relativistic hydrodynamic simulations in a dynamically evolving space-time with the CoCoNuT code. We start from massive low metallicity stellar models evolved up to core gravitational instability, and then follow the subsequent evolution until the system collapses forming a compact remnant. A preliminary study of the collapse outcome is performed by varying the typical parameters of the scenario, such as the initial stellar mass, metallicity, and rotational profile of the stellar progenitor. 1D models (without rotation) have been used to test our newly developed neutrino leakage scheme. This is a fundamental piece of our approach as it allows the central remnant (in all cases considered, a metastable high-mass neutron star) to cool down, eventually collapsing to a black hole. In two dimensions, we show that sufficiently fast rotating cores lead to the formation of Kerr black holes, due to the fall-back of matter surrounding the compact remnant, which has not been successfully unbounded by a precedent supernova shock.Comment: 4 pages, 1 figure, Proceedings of the Spanish Relativity Meeting 201