A powerful hydrodynamic booster for relativistic jets

Velocities close to the speed of light are a robust observational property of the jets observed in microquasars and AGNs, and are expected to be behind much of the phenomenology of GRBs. Yet, the mechanism boosting relativistic jets to such large Lorentz factors is still essentially unknown. Building on recent general-relativistic, multidimensional simulations of progenitors of short GRBs, we discuss a new effect in relativistic hydrodynamics which can act as an efficient booster in jets. This effect is purely hydrodynamical and occurs when large velocities tangential to a discontinuity are present in the flow, yielding Lorentz factors $\Gamma \sim 10^2-10^3$ or larger in flows with moderate initial Lorentz factors. Although without a Newtonian counterpart, this effect can be explained easily through the most elementary hydrodynamical flow: i.e., a relativistic Riemann problem.Comment: 4 pages, 4 figures (1 in color). ApJ Letters accepte

Numerical study of broadband spectra caused by internal shocks in magnetized relativistic jets of blazars

The internal-shocks scenario in relativistic jets has been used to explain the variability of blazars' outflow emission. Recent simulations have shown that the magnetic field alters the dynamics of these shocks producing a whole zoo of spectral energy density patterns. However, the role played by magnetization in such high-energy emission is still not entirely understood. With the aid of \emph{Fermi}'s second LAT AGN catalog, a comparison with observations in the $\gamma$-ray band was performed, in order to identify the effects of the magnetic field.Comment: Proceedings of the meeting The Innermost Regions of Relativistic Jets and Their Magnetic Fields, June 10-14, 2013, Granada (Spain), 4 pages, 3 figure

Numerical simulations of the internal shock model in magnetized relativistic jets of blazars

The internal shocks scenario in relativistic jets is used to explain the variability of the blazar emission. Recent studies have shown that the magnetic field significantly alters the shell collision dynamics, producing a variety of spectral energy distributions and light-curves patterns. However, the role played by magnetization in such emission processes is still not entirely understood. In this work we numerically solve the magnetohydodynamic evolution of the magnetized shells collision, and determine the influence of the magnetization on the observed radiation. Our procedure consists in systematically varying the shell Lorentz factor, relative velocity, and viewing angle. The calculations needed to produce the whole broadband spectral energy distributions and light-curves are computationally expensive, and are achieved using a high-performance parallel code.Comment: 7 pages, 5 figures, proceeding of the "Swift: 10 Years of Discovery" conference (December 2014, Rome, Italy

R-process viable outflows are suppressed in global alpha-viscosity models of collapsar disks

Collapsar disks have been proposed to be rich factories of heavy elements, but the major question of whether their outflows are neutron-rich, and could therefore represent significant sites of the rapid neutron-capture (r-) process, or dominated by iron-group elements remains unresolved. We present the first global models of collapsars that start from a stellar progenitor and self-consistently describe the evolution of the disk, its composition, and its outflows in response to the imploding stellar mantle, using energy-dependent M1 neutrino transport and an alpha-viscosity to approximate turbulent angular-momentum transport. We find that a neutron-rich, neutrino-dominated accretion flow (NDAF) is established only marginally--either for short times or relatively low viscosities--because the disk tends to disintegrate into an advective disk (ADAF) already at relatively high mass-accretion rates, launching powerful outflows but preventing it from developing a hot, dense, and therefore neutron-rich core. Viscous outflows disrupt the star within ~100s with explosion energies close to that of hypernovae. If viscosity is neglected, a stable NDAF with disk mass of about 1Msun is formed but is unable to release neutron-rich ejecta, while it produces a relatively mild explosion powered by a neutrino-driven wind blown off its surface. With ejecta electron fractions close to 0.5, all models presumably produce large amounts of Ni56. Our results suggest that collapsar models based on the alpha-viscosity are inefficient r-process sites and that genuinely magnetohydrodynamic effects may be required to generate neutron-rich outflows. A relatively weak effective viscosity generated by magnetohydrodynamic turbulence would improve the prospects for obtaining neutron-rich ejecta.Comment: 7 pages, 4 figures, 1 table, slightly revised discussion, main results unchanged compared to v1, accepted to ApJ

Numerical simulations of the jetted tidal disruption event Swift J1644+57

In this work we focus on the technical details of the numerical simulations of the non-thermal transient Swift J1644+57, whose emission is probably produced by a two- component jet powered by a tidal disruption event. In this context we provide details of the coupling between the relativistic hydrodynamic simulations and the radiative transfer code. First, we consider the technical demands of one-dimensional simulations of a fast relativistic jet, and show to what extent (for the same physical parameters of the model) do the computed light curves depend on the numerical parameters of the different codes employed. In the second part we explain the difficulties of computing light curves from axisymmetric two dimensonal simulations and discuss a procedure that yields an acceptable tradeoff between the computational cost and the quality of the results

Magnetorotational supernovae: a nucleosynthetic analysis of sophisticated 3D models

Magnetorotational supernovae are a rare type of core-collapse supernovae where the magnetic field and rotation play a central role in the dynamics of the explosion. We present the post-processed nucleosynthesis of state-of-the-art neutrino-MHD supernova models that follow the post explosion evolution for few seconds. We find three different dynamical mechanisms to produce heavy r-process elements: (i) a prompt ejection of matter right after core bounce, (ii) neutron-rich matter that is ejected at late times due to a reconfiguration of the protoneutronstar shape, (iii) small amount of mass ejected with high entropies in the centre of the jet. We investigate total ejecta yields, including the ones of unstable nuclei such as 26Al, 44Ti, 56Ni, and 60Fe. The obtained 56Ni masses vary between 0.01−1M⊙. The latter maximum is compatible with hypernova observations. Furthermore, all of our models synthesize Zn masses in agreement with observations of old metal-poor stars. We calculate simplified light curves to investigate whether our models can be candidates for superluminous supernovae. The peak luminosities obtained from taking into account only nuclear heating reach up to a few ∼1043ergs−1. Under certain conditions, we find a significant impact of the 66Ni decay chain that can raise the peak luminosity up to ∼38 percent compared to models including only the 56Ni decay chain. This work reinforces the theoretical evidence on the critical role of magnetorotational supernovae to understand the occurrence of hypernovae, superluminous supernovae, and the synthesis of heavy elements

Relativistic Magnetohydrodynamics: Renormalized eigenvectors and full wave decomposition Riemann solver

We obtain renormalized sets of right and left eigenvectors of the flux vector Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. The renormalization procedure relies on the characterization of the degeneracy types in terms of the normal and tangential components of the magnetic field to the wavefront in the fluid rest frame. Proper expressions of the renormalized eigenvectors in conserved variables are obtained through the corresponding matrix transformations. Our work completes previous analysis that present different sets of right eigenvectors for non-degenerate and degenerate states, and can be seen as a relativistic generalization of earlier work performed in classical MHD. Based on the full wave decomposition (FWD) provided by the the renormalized set of eigenvectors in conserved variables, we have also developed a linearized (Roe-type) Riemann solver. Extensive testing against one- and two-dimensional standard numerical problems allows us to conclude that our solver is very robust. When compared with a family of simpler solvers that avoid the knowledge of the full characteristic structure of the equations in the computation of the numerical fluxes, our solver turns out to be less diffusive than HLL and HLLC, and comparable in accuracy to the HLLD solver. The amount of operations needed by the FWD solver makes it less efficient computationally than those of the HLL family in one-dimensional problems. However its relative efficiency increases in multidimensional simulations.Comment: 50 pages, 17 figures (2 in color). Submitted to ApJ Suppl. Se