Numerical relativity of compact binaries in the 21st century.

We review the dramatic progress in the simulations of compact objects and compact-object binaries that has taken place in the first two decades of the twenty-first century. This includes simulations of the inspirals and violent mergers of binaries containing black holes and neutron stars, as well as simulations of black-hole formation through failed supernovae and high-mass neutron star-neutron star mergers. Modeling such events requires numerical integration of the field equations of general relativity in three spatial dimensions, coupled, in the case of neutron-star containing binaries, with increasingly sophisticated treatment of fluids, electromagnetic fields, and neutrino radiation. However, it was not until 2005 that accurate long-term evolutions of binaries containing black holes were even possible (Pretorius 2005 Phys. Rev. Lett. 95 121101, Campanelli et al 2006 Phys. Rev. Lett. 96 111101, Baker et al 2006 Phys. Rev. Lett. 96 111102). Since then, there has been an explosion of new results and insights into the physics of strongly-gravitating system. Particular emphasis has been placed on understanding the gravitational wave and electromagnetic signatures from these extreme events. And with the recent dramatic discoveries of gravitational waves from merging black holes by the Laser Interferometric Gravitational Wave Observatory and Virgo, and the subsequent discovery of both electromagnetic and gravitational wave signals from a merging neutron star-neutron star binary, numerical relativity became an indispensable tool for the new field of multimessenger astronomy.

Effects of Neutron-Star Dynamic Tides on Gravitational Waveforms within the Effective-One-Body Approach.

Extracting the unique information on ultradense nuclear matter from the gravitational waves emitted by merging neutron-star binaries requires robust theoretical models of the signal. We develop a novel effective-one-body waveform model that includes, for the first time, dynamic (instead of only adiabatic) tides of the neutron star as well as the merger signal for neutron-star-black-hole binaries. We demonstrate the importance of the dynamic tides by comparing our model against new numerical-relativity simulations of nonspinning neutron-star-black-hole binaries spanning more than 24 gravitational-wave cycles, and to other existing numerical simulations for double neutron-star systems. Furthermore, we derive an effective description that makes explicit the dependence of matter effects on two key parameters: tidal deformability and fundamental oscillation frequency.

Collapse of magnetized hypermassive neutron stars in general relativity.

Hypermassive neutron stars (HMNSs)--equilibrium configurations supported against collapse by rapid differential rotation--are possible transient remnants of binary neutron-star mergers. Using newly developed codes for magnetohydrodynamic simulations in dynamical spacetimes, we are able to track the evolution of a magnetized HMNS in full general relativity for the first time. We find that secular angular momentum transport due to magnetic braking and the magnetorotational instability results in the collapse of an HMNS to a rotating black hole, accompanied by a gravitational wave burst. The nascent black hole is surrounded by a hot, massive torus undergoing quasistationary accretion and a collimated magnetic field. This scenario suggests that HMNS collapse is a possible candidate for the central engine of short gamma-ray bursts.

Magnetized hypermassive neutron-star collapse: a central engine for short gamma-ray bursts.

A hypermassive neutron star (HMNS) is a possible transient formed after the merger of a neutron-star binary. In the latest axisymmetric magnetohydrodynamic simulations in full general relativity, we find that a magnetized HMNS undergoes "delayed" collapse to a rotating black hole (BH) as a result of angular momentum transport via magnetic braking and the magnetorotational instability. The outcome is a BH surrounded by a massive, hot torus with a collimated magnetic field. The torus accretes onto the BH at a quasisteady accretion rate [FORMULA: SEE TEXT]; the lifetime of the torus is approximately 10 ms. The torus has a temperature [FORMULA: SEE TEXT], leading to copious ([FORMULA: SEE TEXT]) thermal radiation that could trigger a fireball. Therefore, the collapse of a HMNS is a promising scenario for generating short-duration gamma-ray bursts and an accompanying burst of gravitational waves and neutrinos.

Dynamical determination of the innermost stable circular orbit of binary neutron stars.

We determine the innermost stable circular orbit (ISCO) of binary neutron stars (BNSs) by performing dynamical simulations in full general relativity. Evolving quasiequilibrium (QE) binaries that begin at different separations, we bracket the location of the ISCO by distinguishing stable circular orbits from unstable plunges. We study Gamma=2 polytropes of varying compactions in both corotational and irrotational equal-mass binaries. For corotational binaries, we find an ISCO orbital angular frequency somewhat smaller than that determined by applying turning-point methods to QE initial data. For the irrotational binaries, the initial data sequences terminate before reaching a turning point, but we find that the ISCO frequency is reached prior to the termination point. Our findings suggest that the ISCO frequency varies with compaction but does not depend strongly on the stellar spin.