
We present our current best estimate of the plausible observing scenarios for
the Advanced LIGO, Advanced Virgo and KAGRA gravitationalwave detectors over
the next several years, with the intention of providing information to
facilitate planning for multimessenger astronomy with gravitational waves. We
estimate the sensitivity of the network to transient gravitationalwave 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
gravitationalwave 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 skyarea 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.

The observational signature of supernova remnants (SNRs) is very complex, in
terms of both their geometrical shape and their spectral properties, dominated
by nonthermal synchrotron and inverseCompton scattering. We propose a
postprocessing method to analyse the broadband emission of SNRs based on
threedimensional hydrodynamical simulations. From the hydrodynamical data, we
estimate the distribution of nonthermal 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.

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.04622). 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.

Context. Global MHD simulations show KelvinHelmholtz (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
backreaction of the magnetic field onto the flow. Aims. We investigate the
amplification of initially weak fields in KH unstable shear flows, and the
backreaction of the field onto the flow. Methods. We use a highresolution
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 smallscale
features. The magnetic field is at most in equipartition with the decaying
shear flow. (abridged)

We assess the importance of the magnetorotational instability in
corecollapse supernovae by an analysis of the growth rates of unstable modes
in typical postcollapse 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 longterm evolution of the saturated turbulent state.

Possible effects of magnetic fields in core collapse supernovae rely on an
efficient amplification of the weak precollapse fields. It has been suggested
that the magnetorotational 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 semiglobal
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
subKeplerian 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]

We have performed 2.5D and 3D simulations of conical jets driven by the
rotation of an ordered, largescale 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.

The possibility that gammaray 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 quasyisotropic 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 nonlinear 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).

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 jetlike outflow.
[Abstract abbreviated]

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 relaxingTVD
method for the solution of the ideal MHD equations, and apply the
constrainttransport method to guarantee a divergencefree evolution of the
magnetic field. We neglect effects due to neutrino transport and employ a
simplified equation of state. The precollapse 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)