The Advanced LIGO and Advanced Virgo gravitational-wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies. We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star, a neutron star–black hole binary and a binary black hole, where we show a cross comparison of results obtained using three independent sampling algorithms. These systems were analyzed with nonspinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations. We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analyzing a set of 100 signals drawn from the prior. We discuss the computational cost of these algorithms, and describe the general and problem-specific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence parameter space.

Parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library

In a single process, the merger of binary neutron star systems combines extreme gravity, the copious emission of gravitational waves, complex microphysics and electromagnetic processes, which can lead to astrophysical signatures observable at the largest redshifts. We review here the recent progress in understanding what could be considered Einstein's richest laboratory, highlighting in particular the numerous significant advances of the last decade. Although special attention is paid to the status of models, techniques and results for fully general-relativistic dynamical simulations, a review is also offered on the initial data and advanced simulations with approximate treatments of gravity. Finally, we review the considerable amount of work carried out on the post-merger phase, including black-hole formation, torus accretion onto the merged compact object, the connection with gamma-ray burst engines, ejected material, and its nucleosynthesis.

Binary Neutron Star Mergers: A Review of Einstein's Richest Laboratory

The grand challenges of contemporary fundamental physics-dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem-all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress. This write-up is an initiative taken within the framework of the European Action on 'Black holes, Gravitational waves and Fundamental Physics'. © 2019 IOP Publishing Ltd.

/pdf/black-holes-gravitational-waves-and-fundamental-physics-a-2xpderzol2.pdf

Black holes, gravitational waves and fundamental physics: a roadmap

Merging binaries of neutron-stars are not only strong sources of gravitational waves, but also have the potential of revealing states of matter at densities and temperatures not accessible in laboratories. A crucial and long-standing question in this context is whether quarks are deconfined as a result of the dramatic increase in density and temperature following the merger. We present the first fully general-relativistic simulations of merging neutron-stars including quarks at finite temperatures that can be switched off consistently in the equation of state. Within our approach, we can determine clearly what signatures a quark-hadron phase transition would leave in the gravitational-wave signal. We show that if after the merger the conditions are met for a phase transition to take place at several times nuclear saturation density, they would lead to a postmerger signal considerably different from the one expected from the inspiral, that can only probe the hadronic part of the equations of state, and to an anticipated collapse of the merged object. We also show that the phase transition leads to a very hot and dense quark core that, when it collapses to a black hole, produces a ringdown signal different from the hadronic one. Finally, in analogy with what is done in heavy-ion collisions, we use the evolution of the temperature and density in the merger remnant to illustrate the properties of the phase transition in a QCD phase diagram.

Signatures of Quark-Hadron Phase Transitions in General-Relativistic Neutron-Star Mergers.

Magnetohydrodynamic turbulence is important in many high-energy astrophysical systems, where instabilities can amplify the local magnetic field over very short timescales. Specifically, the magnetorotational instability and dynamo action have been suggested as a mechanism for the growth of magnetar-strength magnetic fields (of 10^(15) gauss and above) and for powering the explosion of a rotating massive star. Such stars are candidate progenitors of type Ic-bl hypernovae, which make up all supernovae that are connected to long γ-ray bursts. The magnetorotational instability has been studied with local high-resolution shearing-box simulations in three dimensions, and with global two-dimensional simulations, but it is not known whether turbulence driven by this instability can result in the creation of a large-scale, ordered and dynamically relevant field. Here we report results from global, three-dimensional, general-relativistic magnetohydrodynamic turbulence simulations. We show that hydromagnetic turbulence in rapidly rotating protoneutron stars produces an inverse cascade of energy. We find a large-scale, ordered toroidal field that is consistent with the formation of bipolar magnetorotationally driven outflows. Our results demonstrate that rapidly rotating massive stars are plausible progenitors for both type Ic-bl supernovae and long γ-ray bursts, and provide a viable mechanism for the formation of magnetars. Moreover, our findings suggest that rapidly rotating massive stars might lie behind potentially magnetar-powered superluminous supernovae.

A large-scale dynamo and magnetoturbulence in rapidly rotating core-collapse supernovae

We present the new general-relativistic magnetohydrodynamics (GRMHD) capabilities of the Einstein Toolkit, an open-source community-driven numerical relativity and computational relativistic astrophysics code. The GRMHD extension of the Toolkit builds upon previous releases and implements the evolution of relativistic magnetised fluids in the ideal MHD limit in fully dynamical spacetimes using the same shock-capturing techniques previously applied to hydrodynamical evolution. In order to maintain the divergence-free character of the magnetic field, the code implements both hyperbolic divergence cleaning and constrained transport schemes. We present test results for a number of MHD tests in Minkowski and curved spacetimes. Minkowski tests include aligned and oblique planar shocks, cylindrical explosions, magnetic rotors, Alfv\'en waves and advected loops, as well as a set of tests designed to study the response of the divergence cleaning scheme to numerically generated monopoles. We study the code's performance in curved spacetimes with spherical accretion onto a black hole on a fixed background spacetime and in fully dynamical spacetimes by evolutions of a magnetised polytropic neutron star and of the collapse of a magnetised stellar core. Our results agree well with exact solutions where these are available and we demonstrate convergence. All code and input files used to generate the results are available on this http URL This makes our work fully reproducible and provides new users with an introduction to applications of the code.

/pdf/grhydro-a-new-open-source-general-relativistic-2zvgqrr1fq.pdf

GRHydro: A new open source general-relativistic magnetohydrodynamics code for the Einstein Toolkit

In the extreme violence of merger and mass accretion, compact objects like black holes and neutron stars are thought to launch some of the most luminous outbursts of electromagnetic and gravitational wave energy in the Universe. Modeling these systems realistically is a central problem in theoretical astrophysics, but has proven extremely challenging, requiring the development of numerical relativity codes that solve Einstein's equations for the spacetime, coupled to the equations of general relativistic (ideal) magnetohydrodynamics (GRMHD) for the magnetized fluids. Over the past decade, the Illinois Numerical Relativity (ILNR) Group's dynamical spacetime GRMHD code has proven itself as a robust and reliable tool for theoretical modeling of such GRMHD phenomena. However, the code was written "by experts and for experts" of the code, with a steep learning curve that would severely hinder community adoption if it were open-sourced. Here we present IllinoisGRMHD, which is an open-source, highly-extensible rewrite of the original closed-source GRMHD code of the ILNR Group. Reducing the learning curve was the primary focus of this rewrite, with the goal of facilitating community involvement in the code's use and development, as well as the minimization of human effort in generating new science. IllinoisGRMHD also saves computer time, generating roundoff-precision identical output to the original code on adaptive-mesh grids, but nearly twice as fast at scales of hundreds to thousands of cores.

/pdf/illinoisgrmhd-an-open-source-user-friendly-grmhd-code-for-2t7xfrzcmh.pdf

IllinoisGRMHD: An Open-Source, User-Friendly GRMHD Code for Dynamical Spacetimes

We provide additional information on our recent study of the electromagnetic emission produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. As anticipated in a recent letter, our results show that although a dual-jet structure is present, the associated luminosity is ~100 times smaller than the total one, which is predominantly quadrupolar. Here we discuss the details of our implementation of the equations in which the force-free condition is not implemented at a discrete level, but rather obtained via a damping scheme which drives the solution to satisfy the correct condition. We show that this is important for a correct and accurate description of the current sheets that can develop in the course of the simulation. We also study in greater detail the three-dimensional charge distribution produced as a consequence of the inspiral and show that during the inspiral it possesses a complex but ordered structure which traces the motion of the two black holes. Finally, we provide quantitative estimates of the scaling of the electromagnetic emission with frequency, with the diffused part having a dependence that is the same as the gravitational-wave one and that scales as L^(non-coll)_(EM) ≈ Ω^((10/3)–(8/3)), while the collimated one scales as L^(coll)_(EM) ≈ Ω^((5/3)–(6/3)), thus with a steeper dependence than previously estimated. We discuss the impact of these results on the potential detectability of dual jets from supermassive black holes and the steps necessary for more accurate estimates.

/pdf/accurate-simulations-of-binary-black-hole-mergers-in-force-14tvcuqn13.pdf

Accurate Simulations of Binary Black Hole Mergers in Force-free Electrodynamics

We revisit the suggestion that dual jets can be produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. By performing independent calculations of the late inspiral and merger, and by computing the electromagnetic (EM) emission in a way which is consistent with estimates using the Poynting flux, we show that a dual-jet structure is present but energetically subdominant with respect to a non-collimated and predominantly quadrupolar emission, which is similar to the one computed when the binary is in electrovacuum. While our findings set serious restrictions on the detectability of dual jets from coalescing binaries, they also increase the chances of detecting an EM counterpart from these systems.

/pdf/on-the-detectability-of-dual-jets-from-binary-black-holes-3cnogs6fqg.pdf

On the detectability of dual jets from binary black holes

The understanding of strong-field dynamics near black-hole horizons is a long-standing and challenging problem in general relativity. Recent advances in numerical relativity and in the geometric characterization of black-hole horizons open new avenues into the problem. In this first paper in a series of two, we focus on the analysis of the recoil occurring in the merger of binary black holes, extending the analysis initiated in [L. Rezzolla, R. P. Macedo, and J. L. Jaramillo, Phys. Rev. Lett. 104, 221101 (2010).] with Robinson-Trautman spacetimes. More specifically, we probe spacetime dynamics through the correlation of quantities defined at the black-hole horizon and at null infinity. The geometry of these hypersurfaces responds to bulk gravitational fields acting as test screens in a scattering perspective of spacetime dynamics. Within a $3+1$ approach we build an effective-curvature vector from the intrinsic geometry of dynamical-horizon sections and correlate its evolution with the flux of Bondi linear momentum at large distances. We employ this setup to study numerically the head-on collision of nonspinning black holes and demonstrate its validity to track the qualitative aspects of recoil dynamics at infinity. We also make contact with the suggestion that the antikick can be described in terms of a ``slowness parameter'' and how this can be computed from the local properties of the horizon. In a companion paper [J. L. Jaramillo, R. P. Macedo, P. Moesta, and L. Rezzolla, following article, Phys. Rev. D 85, 084031 (2012).] we will further elaborate on the geometric aspects of this approach and on its relation with other approaches to characterize dynamical properties of black-hole horizons.

/pdf/black-hole-horizons-as-probes-of-black-hole-dynamics-i-post-4llqp86vm4.pdf

Philipp Moesta

Papers

GRHydro: A new open source general-relativistic magnetohydrodynamics code for the Einstein Toolkit

IllinoisGRMHD: An Open-Source, User-Friendly GRMHD Code for Dynamical Spacetimes

Accurate Simulations of Binary Black Hole Mergers in Force-free Electrodynamics

On the detectability of dual jets from binary black holes

Black-hole horizons as probes of black-hole dynamics I: post-merger recoil in head-on collisions