• We present our current best estimate of the plausible observing scenarios for the Advanced LIGO, Advanced Virgo and KAGRA gravitational-wave detectors over the next several years, with the intention of providing information to facilitate planning for multi-messenger astronomy with gravitational waves. We estimate the sensitivity of the network to transient gravitational-wave signals for the third (O3), fourth (O4) and fifth observing (O5) runs, including the planned upgrades of the Advanced LIGO and Advanced Virgo detectors. We study the capability of the network to determine the sky location of the source for gravitational-wave signals from the inspiral of binary systems of compact objects, that is BNS, NSBH, and BBH systems. The ability to localize the sources is given as a sky-area probability, luminosity distance, and comoving volume. The median sky localization area (90\% credible region) is expected to be a few hundreds of square degrees for all types of binary systems during O3 with the Advanced LIGO and Virgo (HLV) network. The median sky localization area will improve to a few tens of square degrees during O4 with the Advanced LIGO, Virgo, and KAGRA (HLVK) network. We evaluate sensitivity and localization expectations for unmodeled signal searches, including the search for intermediate mass black hole binary mergers.
• ### A method for computing synchrotron and inverse-Compton emission from hydrodynamic simulations of supernova remnants(1408.0896)

Nov. 13, 2014 astro-ph.SR, astro-ph.HE
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.
• ### Hydrodynamic simulations of the interaction of supernova shock waves with a clumpy environment: the case of the RX J0852.0-4622 (Vela Jr) supernova remnant(1408.0895)

Aug. 5, 2014 astro-ph.SR, astro-ph.HE
Observations in all electromagnetic bands show that many supernova remnants (SNRs) have a very aspherical shape. This can be the result of asymmetries in the supernova explosion or a clumpy circumstellar medium. We study the generation of inhomogeneities and the mixing of elements arising from these two sources in multidimensional hydrodynamic simulations of the propagation of a supernova blast wave into a cloudy environment. We model a specific SNR, Vela Jr (RX J0852.0-4622). By comparing our results with recent observations, we can constrain the properties of the explosion. We find that a very energetic explosion of several 10^{51} erg occurring roughly about 800 years ago is consistent with the shape and emission of the SNR, as well as a supernova with an energy closer to the canonical value of 10^{51} erg a few thousand years ago.
• ### Local simulations of the magnetized Kelvin-Helmholtz instability in neutron-star mergers(1003.6031)

March 31, 2010 astro-ph.SR
Context. Global MHD simulations show Kelvin-Helmholtz (KH) instabilities at the contact surface of two merging neutron stars. That region has been identified as the site of efficient amplification of magnetic fields. However, these global simulations, due to numerical limitations, were unable to determine the saturation level of the field strength, and thus the possible back-reaction of the magnetic field onto the flow. Aims. We investigate the amplification of initially weak fields in KH unstable shear flows, and the back-reaction of the field onto the flow. Methods. We use a high-resolution ideal MHD code to perform 2D and 3D local simulations of shear flows. Results. In 2D, the magnetic field is amplified in less than 0.01ms until it reaches locally equipartition with the kinetic energy. Subsequently, it saturates due to resistive instabilities that disrupt the KH vortex and decelerate the shear flow on a secular time scale. We determine scaling laws of the field amplification with the initial field strength and the grid resolution. In 3D, this hydromagnetic mechanism may be dominated by purely hydrodynamic instabilities limiting the amplification. We find maximum magnetic fields of 10^16 G locally, and r.m.s. maxima within the box of 10^15 G. However, such strong fields exist only for a short period. In the saturated state, the magnetic field is mainly oriented parallel to the shear flow for strong initial fields, while weaker initial fields tend to lead to a more balanced distribution of the field energy. In all models the flow shows small-scale features. The magnetic field is at most in equipartition with the decaying shear flow. (abridged)
• ### Simulations of the Magneto-rotational Instability in Core-Collapse Supernovae(0911.3222)

Nov. 17, 2009 astro-ph.SR
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.
• ### Semi-global simulations of the magneto-rotational instability in core collapse supernovae(0811.1652)

March 17, 2009 astro-ph
Possible effects of magnetic fields in core collapse supernovae rely on an efficient amplification of the weak pre-collapse fields. It has been suggested that the magneto-rotational instability (MRI) leads to rapid field growth. Although MRI studies exist for accretion discs, the application of their results to core collapse supernovae is inhibited as the physics of supernova cores is substantially different from that of accretion discs. We address the problem of growth and saturation of the MRI by means of semi-global simulations, which combine elements of global and local simulations by taking the presence of global background gradients into account and using a local computational grid. We analyze the dispersion relation of the MRI to identify different regimes of the instability. This analysis is complemented by simulations, where we consider a local computational box rotating at sub-Keplerian velocity, and where we allow for a radial entropy gradient. We identify six regimes of the MRI depending on the ratio of the entropy and angular velocity gradient. Our numerical models confirm the instability criteria and growth rates for all relevant regimes. The MRI grows exponentially within milliseconds the flow and magnetic field geometries being dominated by channel flows. The MRI growth ceases once the channels are disrupted by resistive instabilities (due to finite numerical conductivity), and MHD turbulence sets in. From an analysis of the growth rates of the resistive instabilities, we deduce scaling laws for the termination amplitude of the MRI which agree well with our numerical models. We determine the dependence of the development of coherent flow structures in the saturated state on the aspect ratio of the simulation boxes. [abridged]
• ### Kink instabilities in jets from rotating magnetic fields(0809.3165)

Sept. 18, 2008 astro-ph
We have performed 2.5D and 3D simulations of conical jets driven by the rotation of an ordered, large-scale magnetic field in a stratified atmosphere. The simulations cover about three orders of magnitude in distance to capture the centrifugal acceleration as well as the evolution past the Alfven surface. We find that the jets develop kink instabilities, the characteristics of which depend on the velocity profile imposed at the base of the flow. The instabilities are especially pronounced with a rigid rotation profile, which induces a shearless magnetic field. The jet's expansion appears to be limiting the growth of Alfven mode instabilities.
• ### Relativistic Outflows in Gamma-Ray Bursts(astro-ph/0701187)

Jan. 8, 2007 astro-ph
The possibility that gamma-ray bursts (GRBs) were not isotropic emissions was devised theoretically as a way to ameliorate the huge energetic budget implied by the standard fireball model for these powerful phenomena. However, the mechanism by which after the quasy-isotropic release of a few $10^{50}$erg yields a collimated ejection of plasma could not be satisfactory explained analytically. The reason being that the collimation of an outflow by its progenitor system depends on a very complex and non-linear dynamics. That has made necessary the use of numerical simulations in order to shed some light on the viability of some likely progenitors of GRBs. In this contribution I will review the most relevant features shown by these numerical simulations and how they have been used to validate the collapsar model (for long GRBs) and the model involving the merger of compact binaries (for short GRBs).
• ### Axisymmetric simulations of magnetorotational core collapse: Approximate inclusion of general relativistic effects(astro-ph/0602187)

Oct. 31, 2006 astro-ph
We continue our investigations of the magnetorotational collapse of stellar cores discussing simulations performed with a modified Newtonian gravitational potential that mimics general relativistic effects. The approximate TOV potential used in our simulations catches several features of fully relativistic simulations quite well. It is able to correctly reproduce the behavior of models which show a qualitative change both of the dynamics and the gravitational wave signal when switching from Newtonian to fully relativistic simulations. If this is not the case, the Newtonian and the approximate TOV models differ quantitatively. The collapse proceeds to higher densities with the approximate TOV potential allowing for a more efficient amplification of the magnetic field by differential rotation. Sufficiently strong magnetic fields brake down the core's rotation and trigger a contraction phase to higher densities. Several models exhibit two different kinds of shock generation. Due to magnetic braking, a first shock wave created during the initial centrifugal bounce does not suffice to eject any mass, and the core continues to collapse to supranuclear densities. Another stronger shock wave is generated during the second bounce as the core exceeds nuclear matter density. The gravitational wave signal of these models does not fit into the standard classification. Instead it belongs to the signal type IV introduced by us in the first paper of this series. This signal type is more frequent for the approximate relativistic potential than for the Newtonian one. Strongly magnetized models emit a substantial fraction of their GW power at very low frequencies. A flat spectrum between 10 Hz and > 100 kHz denotes the generation of a jet-like outflow. [Abstract abbreviated]
• ### Axisymmetric simulations of magneto--rotational core collapse: dynamics and gravitational wave signal(astro-ph/0510184)

Feb. 8, 2006 astro-ph
We have performed a comprehensive parameter study of the collapse of rotating, strongly magnetized stellar cores in axisymmetry to determine their gravitational wave signature based on the Einstein quadrupole formula. We use a Newtonian explicit magnetohydrodynamic Eulerian code based on the relaxing-TVD method for the solution of the ideal MHD equations, and apply the constraint-transport method to guarantee a divergence--free evolution of the magnetic field. We neglect effects due to neutrino transport and employ a simplified equation of state. The pre--collapse initial models are polytropes in rotational equilibrium with a prescribed degree of differential rotation and rotational energy (~ 1 % of the gravitational energy). The initial magnetic fields are purely poloidal the field strength ranging from 10^10 G to 10^13 G. The evolution of the core, whose collapse is initiated by reducing the gas pressure by a prescribed amount, is followed until a few ten milliseconds past core bounce. The initial magnetic fields are amplified mainly by the differential rotation of the core giving rise to a strong toroidal field component. The poloidal field component grows by compression during collapse, but does not change significantly after core bounce if (abbreviated)