Relativistic Hydrodynamic Simulations Research Articles

ALMA and RATIR observations of GRB 131030A

Abstract We report on the first open-use based Atacama Large Millimeter/submm Array (ALMA) 345 GHz observation for the late afterglow phase of GRB 131030A. The ALMA observation constrained a deep limit at 17.1 d for the afterglow and host galaxy. We also identified a faint submillimeter source (ALMA J2300−0522) near the GRB 131030A position. The deep limit at 345 GHz and multifrequency observations obtained using Swift and RATIR yielded forward-shock modeling with a two-dimensional relativistic hydrodynamic jet simulation and described X-ray excess in the afterglow. The excess was inconsistent with the synchrotron self-inverse Compton radiation from the forward shock. The host galaxy of GRB 131030A and optical counterpart of ALMA J2300−0522 were also identified in the Subaru image. Based on the deep ALMA limit for the host galaxy, the 3σ upper limits of IR luminosity and the star formation rate (SFR) are estimated as LIR < 1.11 × 1011 L⊙ and SFR <18.7 (M⊙ yr−1), respectively. Although the separation angle from the burst location (3${^{\prime\prime}_{.}}$5) was rather large, ALMA J2300−0522 may be one component of the GRB 131030A host galaxy, according to previous host galaxy cases.

Supernova ejecta with a relativistic wind from a central compact object: a unified picture for extraordinary supernovae

The hydrodynamical interaction between freely expanding supernova ejecta and a relativistic wind injected from the central region is studied in analytic and numerical ways. As a result of the collision between the ejecta and the wind, a geometrically thin shell surrounding a hot bubble forms and expands in the ejecta. We use a self-similar solution to describe the early dynamical evolution of the shell and carry out a two-dimensional special relativistic hydrodynamic simulation to follow further evolution. The Rayleigh–Taylor instability inevitably develops at the contact surface separating the shocked wind and ejecta, leading to the complete destruction of the shell and the leakage of hot gas from the hot bubble. The leaking hot materials immediately catch up with the outermost layer of the supernova ejecta and thus different layers of the ejecta are mixed. We present the spatial profiles of hydrodynamical variables and the kinetic energy distributions of the ejecta. We stop the energy injection when a total energy of 1052 erg, which is 10 times larger than the initial kinetic energy of the supernova ejecta, is deposited into the ejecta and follow the subsequent evolution. From the results of our simulations, we consider expected emission from supernova ejecta powered by the energy injection at the centre and discuss the possibility that superluminous supernovae and broad-lined Ic supernovae could be produced by similar mechanisms.

Impact of the gravity of a Schwarzschild black hole upon the Rossby wave instability

In an early work the Rossby Wave Instability (RWI) has been proposed to explain variability thought to originate in the close vicinity of black-holes but this was done in the pseudo-Newtonian approach. Here we present the first general relativistic hydrodynamics simulations of this instability not only proving its theorized existence in a full general relativistic (GR) environment but also studying the effect of the strong gravity on the instability. To that end we performed a set of simulations increasingly closer to the black hole with our new GR version of the MPI-AMRVAC code.This allows us to study the minute changes in the behaviour of the instability.We found that the pseudo-Newtonian approach gives adequate results provided that time shifting induced by the black hole gravity is taken into account. Hence, to view the disc as a distant observer would a full GR ray-tracing post treatment of the simulations is a must.

Solving 3D relativistic hydrodynamical problems with weighted essentially nonoscillatory discontinuous Galerkin methods

Discontinuous Galerkin (DG) methods coupled to WENO algorithms allow high order convergence for smooth problems and for the simulation of discontinuities and shocks. In this work, we investigate WENO-DG algorithms in the context of numerical general relativity, in particular for general relativistic hydrodynamics. We implement the standard WENO method at different orders, a compact (simple) WENO scheme, as well as an alternative subcell evolution algorithm. To evaluate the performance of the different numerical schemes, we study non-relativistic, special relativistic, and general relativistic testbeds. We present the first three-dimensional simulations of general relativistic hydrodynamics, albeit for a fixed spacetime background, within the framework of WENO-DG methods. The most important testbed is a single TOV-star in three dimensions, showing that long term stable simulations of single isolated neutron stars can be obtained with WENO-DG methods.

The influence of circumnuclear environment on the radio emission from TDE jets

Dozens of stellar tidal disruption events (TDEs) have been identified at optical, UV and X-ray wavelengths. A small fraction of these, most notably Swift J1644+57, produce radio synchrotron emission, consistent with a powerful, relativistic jet shocking the surrounding circumnuclear gas. The dearth of similar non-thermal radio emission in the majority of TDEs may imply that powerful jet formation is intrinsically rare, or that the conditions in galactic nuclei are typically unfavorable for producing a detectable signal. Here we explore the latter possibility by constraining the radial profile of the gas density encountered by a TDE jet using a one-dimensional model for the circumnuclear medium which includes mass and energy input from a stellar population. Near the jet Sedov radius of 10$^{18}$ cm, we find gas densities in the range of $n_{18} \sim$ 0.1$-$1000 cm$^{-3}$ across a wide range of plausible star formation histories. Using one- and two-dimensional relativistic hydrodynamical simulations, we calculate the synchrotron radio light curves of TDE jets (as viewed both on and off-axis) across the allowed range of density profiles. We find that bright radio emission would be produced across the plausible range of nuclear gas densities by jets as powerful as Swift J1644+57, and we quantify the relationship between the radio luminosity and jet energy. We use existing radio detections and upper limits to constrain the energy distribution of TDE jets. Radio follow up observations several months to several years after the TDE candidate will strongly constrain the energetics of any relativistic flow.

MONTE CARLO RADIATION TRANSFER SIMULATIONS OF PHOTOSPHERIC EMISSION IN LONG-DURATION GAMMA-RAY BURSTS

ABSTRACT We present MCRaT, a Monte Carlo Radiation Transfer code for self-consistently computing the light curves and spectra of the photospheric emission from relativistic, unmagnetized jets. We apply MCRaT to a relativistic hydrodynamic simulation of a long-duration gamma-ray burst jet, and present the resulting light curves and time-dependent spectra for observers at various angles from the jet axis. We compare our results to observational results and find that photospheric emission is a viable model to explain the prompt phase of long-duration gamma-ray bursts at the peak frequency and above, but faces challenges when reproducing the flat spectrum below the peak frequency. We finally discuss possible limitations of these results both in terms of the hydrodynamics and the radiation transfer and how these limitations could affect the conclusions that we present.

GR-AMRVAC code applications: accretion onto compact objects, boson stars versus black holes

In the close vicinity of a compact object strong gravity imprints its signature onto matter. Systems that contain at least one compact object are observed to exhibit extreme physical properties and typically emit highly energetic radiation. The nature of the compact objects that produce the strongest gravitational fields is to date not settled. General relativistic numerical simulations of fluid dynamics around black holes, neutron stars, and other compact objects such as boson stars (BSs) may give invaluable insights into this fundamental question. In order to study the behavior of fluid in the strong gravity regime of an arbitrary compact object we develop a new general relativistic hydrodynamics code. To this end we extend the existing versatile adaptive mesh refinement code MPI-AMRVAC into a general relativistic hydrodynamics framework and adapt it for the use of numerically given spacetime metrics. In the present article we study accretion flows in the vicinity of various types of BSs whose numerical metrics are calculated by the KADATH spectral solver library. We design specific tests to check the reliability of any code intending to study BSs and compare the solutions with those obtained in the context of Schwarzschild black holes. We perform the first ever general relativistic hydrodynamical simulations of gas accretion by a BS. The behavior of matter at small distances from the center of a BS differs notably from the black hole case. In particular we demonstrate that in the context of Bondi spherical accretion the mass accretion rate onto non-rotating BSs remains constant whereas it increases for Schwarzschild black holes. We also address the scenario of non-spherical accretion onto BSs and show that this may trigger mass ejection from the interior of the BS. This striking feature opens the door to forthcoming investigations regarding accretion-ejection flows around such types of compact objects.

Relativistic jet feedback in high-redshift galaxies – I. Dynamics

We present the results of three dimensional relativistic hydrodynamic simulations of interaction of AGN jets with a dense turbulent two-phase interstellar medium, which would be typical of high redshift galaxies. We describe the effect of the jet on the evolution of the density of the turbulent ISM. The jet driven energy bubble affects the gas to distances up to several kiloparsecs from the injection region. The shocks resulting from such interactions create a multi-phase ISM and radial outflows. One of the striking result of this work is that low power jets (P_jet < 10^{43} erg/s) although less efficient in accelerating clouds, are trapped in the ISM for a longer time and hence affect the ISM over a larger volume. Jets of higher power drill through with relative ease. Although the relativistic jets launch strong outflows, there is little net mass ejection to very large distances, supporting a galactic fountain scenario for local feedback.

Coupling hydrodynamics and radiation calculations for star-jet interactions in active galactic nuclei

Stars and their winds can contribute to the non-thermal (NT) emission in extragalactic jets. Given the complexity of jet-star interactions, the properties of the resulting emission are strongly linked to those of the emitting flows. We simulate the interaction between a stellar wind and a relativistic extragalactic jet and use the hydrodynamic results to compute the NT emission under different conditions. We perform relativistic axisymmetric hydrodynamical simulations of a relativistic jet interacting with a supersonic, non-relativistic stellar wind. We compute the corresponding streamlines out of the simulation results, and calculate the injection, evolution, and emission of NT particles accelerated in the jet shock, focusing on electrons or $e^\pm$-pairs. Several cases are explored, considering different jet-star interaction locations, magnetic fields and observer lines of sight. The jet luminosity and star properties are fixed, but the results are easily scalable under changes of these parameters. Individual jet-star interactions produce synchrotron and inverse Compton (IC) emission that peaks from X-rays to MeV energies (depending on the magnetic field), and at $\sim 100-1000$ GeV (depending on the stellar type), respectively. The radiation spectrum is hard in the scenarios explored here due to non-radiative cooling dominance, as low-energy electrons are efficiently advected even under relatively high magnetic fields. Interactions of jets with cold stars lead to an even harder IC spectrum because of the Klein-Nishina effect in the cross-section. Doppler boosting has a strong impact on the observer luminosity. The emission levels for individual interactions found here are in the line of previous, more approximate, estimates, strengthening the hypothesis that collective jet-star interactions could significantly contribute at high energies under efficient particle acceleration.

Parametric initial conditions for core-collapse supernova simulations

We investigate a method to construct parametrized progenitor models for core-collapse supernova simulations. Different from all modern core-collapse supernova studies, which rely on progenitor models from stellar evolution calculations, we follow the methodology of Baron & Cooperstein (1990) to construct initial models. Choosing parametrized spatial distributions of entropy and electron fraction as a function of mass coordinate and solving the equation of hydrostatic equilibrium, we obtain the initial density structures of our progenitor models. First, we calculate structures with parameters fitting broadly the evolutionary model s11.2 of Woosley et al. (2002). We then demonstrate the reliability of our method by performing general relativistic hydrodynamic simulations in spherical symmetry with the isotropic diffusion source approximation to solve the neutrino transport. Our comprehensive parameter study shows that initial models with a small central entropy ($\lesssim 0.4\,k_B$ nucleon$^{-1}$) can explode even in spherically symmetric simulations. Models with a large entropy ($\gtrsim 6\,k_B$ nucleon$^{-1}$) in the Si/O layer have a rather large explosion energy ($\sim 4\times 10^{50}$ erg) at the end of the simulations, which is still rapidly increasing.

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.

Spectral evolution of flaring blazars from numerical simulations

High resolution Very Long Baseline Interferometry (VLBI) observations of Active Galactic Nuclei (AGN) revealed traveling and stationary or quasi-stationary radio-components in several blazar jets. The traveling ones are in general interpreted as shock waves generated by pressure perturbations injected at the jet nozzle. The stationary features can be interpreted as recollimation shocks in non-pressure matched jets if they show a quasi-symmetric bump in the spectral index distribution. In some jets there may be interactions between the two kinds of shocks. These shock--shock interactions are observable with VLBI techniques, and their signature should also be imprinted on the single--dish light curves. We performed relativistic hydrodynamic (RHD) simulations of over-pressured and pressure-matched jets. To simulate the shock interaction we injected a perturbation at the jet nozzle once a steady-state was reached. We computed the non-thermal emission (including adiabatic and synchotron losses) resulting from the simulation. We show that the injection of perturbations in a jet can produce a bump in emission at GHz frequencies previous to the main flare, which is produced when the perturbation fills the jet in the observer's frame. The detailed analysis of our simulations and the non-thermal emission calculations show that interaction between a recollimation shock and traveling shock produce a typical and clear signature in both the single--dish light curves and in the VLBI observations: the flaring peaks are higher and delayed with respect to the evolution of a perturbation through a conical jet. This fact can allow to detect such interactions for stationary components lying outside of the region in where the losses are dominated by inverse Compton scattering.

Jet–intracluster medium interaction in Hydra A – II. The effect of jet precession

We present three-dimensional relativistic hydrodynamical simulations of a precessing jet interacting with the intracluster medium and compare the simulated jet structure with the observed structure of the Hydra A northern jet. For the simulations, we use jet parameters obtained in the parameter space study of the first paper in this series and probe different values for the precession period and precession angle. We find that for a precession period P ≈ 1 Myr and a precession angle ψ ≈ 20°, the model reproduces (i) the curvature of the jet, (ii) the correct number of bright knots within 20 kpc at approximately correct locations and (iii) the turbulent transition of the jet to a plume. The Mach number of the advancing bow shock ≈1.85 is indicative of gentle cluster atmosphere heating during the early stages of the AGN's activity.

Jet-induced modifications of the characteristic of the bulk nuclear matter

We present our studies on jet induced modifications of the characteristics of bulk nuclear matter. To describe such matter, we use efficient relativistic hydrodynamic simulations in (3+1)-dimension, employing the Graphics Processing Unit (GPU) in the parallel programming framework. We use Cartesian coordinates in the calculations to ensure a high spatial resolution that is constant throughout the evolution of the system. We show our results on how jets modify the hydrodynamics fields and discuss the implications.

PHOTOSPHERIC EMISSION FROM COLLAPSAR JETS IN 3D RELATIVISTIC HYDRODYNAMICS

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.

Grid-based Methods in Relativistic Hydrodynamics and Magnetohydrodynamics

An overview of grid-based numerical methods used in relativistic hydrodynamics (RHD) and magnetohydrodynamics (RMHD) is presented. Special emphasis is put on a comprehensive review of the application of high-resolution shock-capturing methods. Results of a set of demanding test bench simulations obtained with different numerical methods are compared in an attempt to assess the present capabilities and limits of the various numerical strategies. Applications to three astrophysical phenomena are briefly discussed to motivate the need for and to demonstrate the success of RHD and RMHD simulations in their understanding. The review further provides FORTRAN programs to compute the exact solution of the Riemann problem in RMHD, and to simulate 1D RMHD flows in Cartesian coordinates.Electronic Supplementary MaterialSupplementary material is available for this article at 10.1007/lrca-2015-3.

GENERAL RELATIVISTIC HYDRODYNAMIC SIMULATION OF ACCRETION FLOW FROM A STELLAR TIDAL DISRUPTION

We study how the matter dispersed when a supermassive black hole tidally disrupts a star joins an accretion flow. Combining a relativistic hydrodynamic simulation of the stellar disruption with a relativistic hydrodynamics simulation of the tidal debris motion, we track such a system until ~80% of the stellar mass bound to the black hole has settled into an accretion flow. Shocks near the stellar pericenter and also near the apocenter of the most tightly-bound debris dissipate orbital energy, but only enough to make the characteristic radius comparable to the semi-major axis of the most-bound material, not the tidal radius as previously thought. The outer shocks are caused by post-Newtonian effects, both on the stellar orbit during its disruption and on the tidal forces. Accumulation of mass into the accretion flow is non-monotonic and slow, requiring ~3--10x the orbital period of the most tightly-bound tidal streams, while the inflow time for most of the mass may be comparable to or longer than the mass accumulation time. Deflection by shocks does, however, remove enough angular momentum and energy from some mass for it to move inward even before most of the mass is accumulated into the accretion flow. Although the accretion rate rises sharply and then decays roughly as a power-law, its maximum is ~0.1x the previous expectation, and the duration of the peak is ~5x longer than previously predicted. The geometric mean of the black hole mass and stellar mass inferred from a measured event timescale is therefore ~0.2x the value given by classical theory.

Non-thermal emission from standing relativistic shocks: an application to red giant winds interacting with AGN jets

Galactic and extragalactic relativistic jets have rich environments that are full of moving objects, such as stars and dense clumps. These objects can enter into the jets and generate shocks and non-thermal emission. We characterize the emitting properties of the downstream region of a standing shock formed due to the interaction of a relativistic jet with an obstacle. We focus on the case of red giants interacting with an extragalactic jet. We perform relativistic axisymmetric hydrodynamical simulations of a relativistic jet meeting an obstacle of very large inertia. The results are interpreted in the framework of a red giant whose dense and slow wind interacts with the jet of an active galactic nucleus. Assuming that particles are accelerated in the standing shock generated in the jet as it impacts the red giant wind, we compute the non-thermal particle distribution, the Doppler boosting enhancement, and the non-thermal luminosity in gamma rays. The available non-thermal energy from jet-obstacle interactions is potentially enhanced by a factor of $\sim 100$ when accounting for the whole surface of the shock induced by the obstacle, instead of just the obstacle section. The observer gamma-ray luminosity, including the flow velocity and Doppler boosting effects, can be ~300(g/10)^2 times higher than when the emitting flow is assumed at rest and only the obstacle section is considered, where g is the jet Lorentz factor. For a whole population of red giants inside the jet of an AGN, the predicted persistent gamma-ray luminosities may be potentially detectable for a jet pointing to the observer. Obstacles interacting with relativistic outflows, for instance clouds and populations of stars for extragalactic jets, or stellar wind inhomogeneities in microquasar jets and in winds of pulsars in binaries, should be taken into account when investigating the non-thermal emission from these sources.

The SS433 jet from subparsec to parsec scales

Context. Relativistic jets associated with compact objects, as in the X-ray binary SS433, are known to be multiscale because they spawn over many orders of magnitude in distance. Here we model the precessing SS433 jet and study its dynamics from 𝒪 (0.01) to 𝒪(1) parsec scales.Aims. We aim to solve the discrepancy between the observations on a 0.1 pc scale of SS433, where the jet is clearly precessing with an angle of 20°, and the larger scale observations where the jet of SS433 interacts with the associated supernova remnant W50, requiring a precessing angle of 10°. Methods. We use 3D special relativistic hydrodynamical simulations on a domain of a scale of 1 pc. We use the finite volume code MPI-AMRVAC, solving the relativistic variant of the Euler equations. To cover lengthscale variations from 𝒪(0.001) pc as the jet beam width up to the domain size, we take full advantage of code parallelization and its adaptive mesh refinement scheme.Results. We found that by means of a simple hydrodynamical process, the jet of SS433 can transit from a precessing jet with an angle of 20°, to a continuous hollow non-precessing jet with a smaller opening angle of about 10°. Successive windings of the precessing jet helix undergo gradual deceleration by ISM interaction, to ultimately merge in a hollow straight jet at distances where the ram pressure of individual jet elements match the ISM pressure at about 0.068 pc from the source.Conclusions. We solve the discrepancy with an elegant and simple model that does not require the jet of SS433 to undergo any temporal changes in jet injection dynamics, but does so as a consequence of a hydrodynamically enforced spatial recollimation. Our simulation thus serves to validate simpler model prescriptions for SS433 on large scales, where a continuous jet profile suffices.

Simulations of an inhomogeneous stellar wind interacting with a pulsar wind in a binary system

Binary systems containing a massive star and a non-accreting pulsar present strong interaction between the stellar and the pulsar winds. The properties of this interaction, which largely determine the non-thermal radiation in these systems, strongly depend on the structure of the stellar wind, which can be clumpy or strongly anisotropic, as in Be stars. We study numerically the influence of inhomogeneities in the stellar wind on the structure of the two-wind interaction region. We carried out for the first time axisymmetric, relativistic hydrodynamical simulations, with Lorentz factors of ~6 and accounting for the impact of instabilities, to study the impact in the two-wind interaction structure of an over-dense region of the stellar wind. We also followed the evolution of this over-dense region or clump as it faces the impact of the pulsar wind. For typical system parameters, and adopting a stellar wind inhomogeneity with a density contrast >~10, clumps with radii of a few percent of the binary size can significantly perturb the two-wind interaction region, pushing the two-wind interface to <~40% of the initial distance to the pulsar. After it is shocked, the inhomogeneity quickly expands and is disrupted when it reaches the smallest distance to the pulsar. It eventually fragments, being advected away from the binary system. The whole interaction region is quite unstable, and the shocked pulsar wind can strongly change under small perturbations. We confirm the sensitive nature of the two-wind interaction structure to perturbations, in particular when the stellar wind is inhomogeneous. For realistic over-dense regions of the stellar wind, the interaction region may shrink by a factor of a few, with the shocked flow presenting a complex spatial and temporal pattern. This can lead to strong variations in the non-thermal radiation.