Core-Collapse Supernovae

MHD Simulations of Non-rotating Stellar Core Collapse with Neutrinos

Calibration of advanced Virgo and reconstruction of the gravitational wave signal h(t) during the observing run O2

In August 2017, Advanced Virgo joined Advanced LIGO for the end of the O2 run, leading to the first gravitational waves detections with the three-detector network. This paper describes the Advanced Virgo calibration and the gravitational wave strain h(t) reconstruction during O2. The methods are the same as the ones developed for the initial Virgo detector and have already been described in previous publications, this paper summarizes the differences and emphasis is put on estimating systematic uncertainties. Three versions of the h(t) signal have been computed for the Virgo O2 run, an online version and two post-run reprocessed versions with improved detector calibration and reconstruction…

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.

The CoCoNuT code: from neutron star oscillations to supernova explosions

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…

Termination of the MRI via parasitic instabilities in core-collapse supernovae: influence of numerical methods

We study the influence of numerical methods and grid resolution on the termination of the magnetorotational instability (MRI) by means of parasitic instabilities in three-dimensional shearing-disc simulations reproducing typical conditions found in core-collapse supernovae. Whether or not the MRI is able to amplify weak magnetic fields in this context strongly depends, among other factors, on the amplitude at which its growth terminates. The qualitative results of our study do not depend on the numerical scheme. In all our models, MRI termination is caused by Kelvin-Helmholtz instabilities, consistent with theoretical predictions. Quantitatively, however, there are differences, but numerica…

Axisymmetric simulations of magnetorotational core collapse: approximate inclusion of general relativistic effects

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 TO…

Magnetic field amplification in collapsing, non-rotating stellar cores

Nucleosynthesis in jet-driven and jet-associated supernovae

In contrast to regular core-collapse supernovae, explosions of rapidly rotating massive stars can develop jets, fast collimated outflows directed along the rotational axis. Depending on the rate of rotation and the magnetic field strength before collapse as well as on possible mechanisms amplifying the magnetic field, such a core can explode magnetorotationally rather than via the standard supernova mechanism based on neutrino heating. This scenario can explain the highest kinetic energies observed in the class of hypernovae. On longer time scales, rotation and magnetic fields can play an important role in the engine of long gamma-ray burst powered by proto-magnetars or hyperaccreting black…

Search for intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network

Gravitational wave astronomy has been firmly established with the detection of gravitational waves from the merger of ten stellar mass binary black holes and a neutron star binary. This paper reports on the all-sky search for gravitational waves from intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network. The search uses three independent algorithms: two based on matched filtering of the data with waveform templates of gravitational wave signals from compact binaries, and a third, model-independent algorithm that employs no signal model for the incoming signal. No intermediate mass black hole binary event was detected in this sear…

Interpreting Observations of GRBs with Numerical Simulations

Hydromagnetic instabilities and magnetic field amplification in core collapse supernovae

Some of the most violent events in the universe, the gamma ray burst, could be related to the gravitational collapse of massive stellar cores. The recent association of long GRBs to some class of type Ic supernova seems to support this view. In such scenario fast rotation, strong magnetic fields and general relativistic effects are key ingredients. It is thus important to understand the mechanism that amplifies the magnetic field under that conditions. I present global simulations of the magneto-rotational collapse of stellar cores in general relativity and semi-global simulations of hydromagnetic instabilities under core collapse conditions. I discuss effect of the magneto-rotational insta…

Magnetorotational core collapse of possible GRB progenitors – I. Explosion mechanisms

We investigate the explosion of stars with zero-age main-sequence masses between 20 and 35 solar masses and varying degrees of rotation and magnetic fields including ones commonly considered progenitors of gamma-ray bursts (GRBs). The simulations, combining special relativistic magnetohydrodynamics, a general relativistic approximate gravitational potential, and two-moment neutrino transport, demonstrate the viability of different scenarios for the post-bounce evolution. Having formed a highly massive proto-neutron star (PNS), several models launch successful explosions, either by the standard supernova mechanism based on neutrino heating and hydrodynamic instabilities or by magnetorotation…

Local simulations of the magnetized Kelvin-Helmholtz instability in neutron-star mergers

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 …

Simulations of core collapse using two-moment neutrino transport

Spectral multi-flavour neutrino transport for sipernova simulations

Semi-global simulations of the magneto-rotational instability in core collapse supernovae

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. W…

Powering Short GRBs by Mergers of Moderately Magnetized Neutron Stars

Protomagnetar and black hole formation in high-mass stars

Using axisymmetric simulations coupling special relativistic MHD, an approximate post-Newtonian gravitational potential and two-moment neutrino transport, we show different paths for the formation of either protomagnetars or stellar mass black holes. The fraction of prototypical stellar cores which should result in collapsars depends on a combination of several factors, among which the structure of the progenitor star and the profile of specific angular momentum are probably the foremost. Along with the implosion of the stellar core, we also obtain supernova-like explosions driven by neutrino heating and hydrodynamic instabilities or by magneto-rotational effects in cores of high-mass stars…

Numerical viscosity in simulations of the two-dimensional Kelvin-Helmholtz instability

The Kelvin-Helmholtz instability serves as a simple, well-defined setup for assessing the accuracy of different numerical methods for solving the equations of hydrodynamics. We use it to extend our previous analysis of the convergence and the numerical dissipation in models of the propagation of waves and in the tearing-mode instability in magnetohydrodynamic models. To this end, we perform two-dimensional simulations with and without explicit physical viscosity at different resolutions. A comparison of the growth of the modes excited by our initial perturbations allows us to estimate the effective numerical viscosity of two spatial reconstruction schemes (fifth-order monotonicity preservin…

Search for Gravitational-wave Signals Associated with Gamma-Ray Bursts during the Second Observing Run of Advanced LIGO and Advanced Virgo

We present the results of targeted searches for gravitational-wave transients associated with gamma-ray bursts during the second observing run of Advanced LIGO and Advanced Virgo, which took place from 2016 November to 2017 August. We have analyzed 98 gamma-ray bursts using an unmodeled search method that searches for generic transient gravitational waves and 42 with a modeled search method that targets compact-binary mergers as progenitors of short gamma-ray bursts. Both methods clearly detect the previously reported binary merger signal GW170817, with p-values of <9.38 10-6 (modeled) and 3.1 10-4 (unmodeled). We do not find any significant evidence for gravitational-wave signals associate…

A Fermi Gamma-Ray Burst Monitor Search for Electromagnetic Signals Coincident with Gravitational-wave Candidates in Advanced LIGO's First Observing Run

We present a search for prompt gamma-ray counterparts to compact binary coalescence gravitational wave (GW) candidates from Advanced LIGO's first observing run (O1). As demonstrated by the multimessenger observations of GW170817/GRB 170817A, electromagnetic and GW observations provide complementary information about the astrophysical source and, in the case of weaker candidates, may strengthen the case for an astrophysical origin. Here we investigate low-significance GW candidates from the O1 compact-binary coalescence searches using the Fermi Gamma-ray Burst Monitor (GBM), leveraging its all-sky and broad energy coverage. Candidates are ranked and compared to background to measure signific…

High-order methods for the simulation of hydromagnetic instabilities in core-collapse supernovae

AbstractWe present an assessment of the accuracy of a recently developed MHD code used to study hydromagnetic flows in supernovae and related events. The code, based on the constrained transport formulation, incorporates unprecedented ultra-high-order methods (up to 9th order) for the reconstruction and the most accurate approximate Riemann solvers. We estimate the numerical resistivity of these schemes in tearing instability simulations.

Search for GW signals associated with GRBs

We present the results of targeted searches for gravitational-wave transients associated with gamma-ray bursts during the second observing run of Advanced LIGO and Advanced Virgo, which took place from 2016 November to 2017 August. We have analyzed 98 gamma-ray bursts using an unmodeled search method that searches for generic transient gravitational waves and 42 with a modeled search method that targets compact-binary mergers as progenitors of short gamma-ray bursts. Both methods clearly detect the previously reported binary merger signal GW170817, with p-values of <9.38x10^-6^ (modeled) and 3.1x10^-4^ (unmodeled). We do not find any significant evidence for gravitational-wave signals assoc…