Computation of gravitational waves from inspiraling binary neutron stars in quasiequilibrium circular orbits : Formulation and calibration

After a short review of prominent properties of gravitational waves and the newly born gravitational astronomy, we focus on theoretical aspects. Analytic approximation methods in general relativity have played a crucial role in the recent discoveries of gravitational waves. They are used to build theoretical template banks for searching and analyzing the signals in the ground-based detectors LIGO and Virgo, and, further ahead, space-based LISA-like detectors. In particular, t...

To be observed and analyzed by the network of gravitational wave detectors on ground (LIGO, VIRGO, etc.) and by the future detectors in space (LISA, etc.), inspiralling compact binaries --- binary star systems composed of neutron stars and/or black holes in their late stage of evolution --- require high-accuracy templates predicted by general relativity theory. The gravitational waves emitted by these very relativistic systems can be accurately modelled using a high-order pos...

We investigate close binary neutron stars in quasiequilibrium states in a general relativistic framework. We assume conformal flatness for the spatial metric and irrotational velocity field for the neutron stars. We adopt the polytropic equation of state. The computation is performed for the polytropic index n(=0.5, 0.66667, 0.8, 1, 1.25), and compactness of neutron stars M/R(=0.03 - 0.3). Results of this paper are as follows. (i) The sequences of the irrotational binary are ...

Extending our previous studies, we perform high-resolution simulations of inspiraling binary neutron stars in numerical relativity. We thoroughly carry through a convergence study in our currently available computational resources with the smallest grid spacing of $\approx 63$--86~meter for the neutron-star radius 10.9--13.7\,km. The estimated total error in the gravitational-wave phase is of order 0.1~rad for the total phase of $\gtrsim 210$\,rad in the last $\sim 15$--16 in...

We develop a model for frequency-domain gravitational waveforms from inspiraling binary neutron stars. Our waveform model is calibrated by comparison with hybrid waveforms constructed from our latest high-precision numerical-relativity waveforms and the SEOBNRv2T waveforms in the frequency range of $10$--$1000\,{\rm Hz}$. We show that the phase difference between our waveform model and the hybrid waveforms is always smaller than $0.1\, {\rm rad}$ for the binary tidal deformab...

Despite the recent rapid progress in numerical relativity, a convergence order less than the second has so far plagued codes solving the Einstein-Euler system of equations. We report simulations of the inspiral of binary neutron stars in quasi-circular orbits computed with a new code employing high-order, high-resolution shock-capturing, finite-differencing schemes that, for the first time, go beyond the second-order barrier. In particular, without any tuning or alignment, we...

We construct closed-form gravitational waveforms (GWs) with tidal effects for the coalescence and merger of binary neutron stars. The method relies on a new set of eccentricity-reduced and high-resolution numerical relativity (NR) simulations and is composed of three steps. First, tidal contributions to the GW phase are extracted from the time-domain NR data. Second, those contributions are employed to fix high-order coefficients in an effective and resummed post-Newtonian ex...

The gravitational wave signal from an inspiralling binary neutron star system will contain detailed information about tidal coupling in the system, and thus, about the internal physics of the neutron stars. To extract this information will require highly accurate models for the gravitational waveform. We present here a calculation of the gravitational wave signal from a binary with quadrupolar tidal interactions which includes all post-1-Newtonian-order effects in both the co...

Two relations, the virial relation $M_{\rm ADM}=M_{\rm K}$ and the first law in the form $\delta M_{\rm ADM}=\Omega \delta J$, should be satisfied by a solution and a sequence of solutions describing binary compact objects in quasiequilibrium circular orbits. Here, $M_{\rm ADM}$, $M_{\rm K}$, $J$, and $\Omega$ are the ADM mass, Komar mass, angular momentum, and orbital angular velocity, respectively. $\delta$ denotes an Eulerian variation. These two conditions restrict the al...

We report on numerical results from an independent formalism to describe the quasi-equilibrium structure of nonsynchronous binary neutron stars in general relativity. This is an important independent test of controversial numerical hydrodynamic simulations which suggested that nonsynchronous neutron stars in a close binary can experience compression and even collapse prior to the last stable circular orbit. We show that the interior density indeed increases as irrotational bi...