Merger and Postmerger of Binary Neutron Stars with a Quark-Hadron Crossover Equation of State.

Fully general-relativistic binary-neutron-star (BNS) merger simulations with quark-hadron crossover (QHC) equations of state (EOS) are studied for the first time. In contrast to EOS with purely hadronic matter or with a first-order quark-hadron phase transition (1PT), in the transition region QHC EOS show a peak in sound speed and thus a stiffening. We study the effects of such stiffening in the merger and postmerger gravitational (GW) signals. Through simulations in the binary-mass range 2.5

Binary neutron star mergers: a review of Einstein's richest laboratory.

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.

Gravitational wave content and stability of uniformly, rotating, triaxial neutron stars in general relativity.

Targets for ground-based gravitational wave interferometers include continuous, quasiperiodic sources of gravitational radiation, such as isolated, spinning neutron stars. In this work, we perform evolution simulations of uniformly rotating, triaxially deformed stars, the compressible analogs in general relativity of incompressible, Newtonian Jacobi ellipsoids. We investigate their stability and gravitational wave emission. We employ five models, both normal and supramassive, and track their evolution with different grid setups and resolutions, as well as with two different evolution codes. We find that all models are dynamically stable and produce a strain that is approximately one-tenth the average value of a merging binary system. We track their secular evolution and find that all our stars evolve toward axisymmetry, maintaining their uniform rotation, rotational kinetic energy, and angular momentum profiles while losing their triaxiality.

Constraining the equation of state of neutron stars from binary mergers.

Determining the equation of state of matter at nuclear density and hence the structure of neutron stars has been a riddle for decades. We show how the imminent detection of gravitational waves from merging neutron star binaries can be used to solve this riddle. Using a large number of accurate numerical-relativity simulations of binaries with nuclear equations of state, we find that the postmerger emission is characterized by two distinct and robust spectral features. While the high-frequency peak has already been associated with the oscillations of the hypermassive neutron star produced by the merger and depends on the equation of state, a new correlation emerges between the low-frequency peak, related to the merger process, and the total compactness of the stars in the binary. More importantly, such a correlation is essentially universal, thus providing a powerful tool to set tight constraints on the equation of state. If the mass of the binary is known from the inspiral signal, the combined use of the two frequency peaks sets four simultaneous constraints to be satisfied. Ideally, even a single detection would be sufficient to select one equation of state over the others. We test our approach with simulated data and verify it works well for all the equations of state considered.

Analytic modeling of tidal effects in the relativistic inspiral of binary neutron stars.

To detect the gravitational-wave (GW) signal from binary neutron stars and extract information about the equation of state of matter at nuclear density, it is necessary to match the signal with a bank of accurate templates. We present the two longest (to date) general-relativistic simulations of equal-mass binary neutron stars with different compactnesses, C=0.12 and C=0.14, and compare them with a tidal extension of the effective-one-body (EOB) model. The typical numerical phasing errors over the ≃22   GW cycles are Δϕ≃±0.24   rad. By calibrating only one parameter (representing a higher-order amplification of tidal effects), the EOB model can reproduce, within the numerical error, the two numerical waveforms essentially up to the merger. By contrast, the third post-Newtonian Taylor-T4 approximant with leading-order tidal corrections dephases with respect to the numerical waveforms by several radians.

Challenging the paradigm of singularity excision in gravitational collapse.

A paradigm deeply rooted in modern numerical relativity calculations prescribes the removal of those regions of the computational domain where a physical singularity may develop. We here challenge this paradigm by performing three-dimensional simulations of the collapse of uniformly rotating stars to black holes without excision. We show that this choice, combined with suitable gauge conditions and the use of minute numerical dissipation, improves dramatically the long-term stability of the evolutions. In turn, this allows for the calculation of the waveforms well beyond what was previously possible, providing information on the black-hole ringing and setting a new mark on the present knowledge of the gravitational-wave emission from the stellar collapse to a rotating black hole.