Kōji Uryū

New numerical method for constructing quasiequilibrium sequences of irrotational binary neutron stars in general relativity

We propose a new numerical method to compute quasi-equilibrium sequences of general relativistic irrotational binary neutron star systems. It is a good approximation to assume that (1) the binary star system is irrotational, i.e. the vorticity of the flow field inside component stars vanishes everywhere (irrotational flow), and (2) the binary star system is in quasi-equilibrium, for an inspiraling binary neutron star system just before the coalescence as a result of gravitational wave emission. We can introduce the velocity potential for such an irrotational flow field, which satisfies an elliptic partial differential equation (PDE) with a Neumann type boundary condition at the stellar surface. For a treatment of general relativistic gravity, we use the Wilson--Mathews formulation, which assumes conformal flatness for spatial components of metric. In this formulation, the basic equations are expressed by a system of elliptic PDEs. We have developed a method to solve these PDEs with appropriate boundary conditions. The method is based on the established prescription for computing equilibrium states of rapidly rotating axisymmetric neutron stars or Newtonian binary systems. We have checked the reliability of our new code by comparing our results with those of other computations available. We have also performed several convergence tests. By using this code, we have obtained quasi-equilibrium sequences of irrotational binary star systems with strong gravity as models for final states of real evolution of binary neutron star systems just before coalescence. Analysis of our quasi-equilibrium sequences of binary star systems shows that the systems may not suffer from dynamical instability of the orbital motion and that the maximum density does not increase as the binary separation decreases.

The influence of quark matter at high densities on binary neutron star mergers

We consider the influence of potential quark matter existing at high densities in neutron star interiors on gravitational waves (GW) emitted in a binary neutron star merger event. Two types of equations of state (EoS) at zero temperatures are used, one describing pure nuclear matter, the other nuclear matter with a phase transition to quark matter at very high densities. Binary equilibrium sequences close to the innermost stable circular orbit (ISCO) are calculated to determine the GW frequencies just before merger. It is found that EoS effects begin to play a role for gravitational masses larger than M∞ ≃ 1.5M⊙. The difference in the gravitational wave frequency at the ISCO grows to up to ≃ 10% for the maximal allowed mass given by the EoSs used. Then, we perform 3D hydrodynamic simulations for each EoS varying the initial mass and determine the characteristic GW frequencies of the merger remnants. The implications of quark matter show up mainly in a different collapse behaviour of the merger remnant. If the collapse does not take place immediately after merger, we find a phase difference between two EoS’s in the post-merger GW signal. We also compare the GW frequencies emitted by the merger remnant to values from simulations using a polytropic EoS and find an imprint of the non-constant adiabatic index of our EoSs. All calculations are based on the conformally flat (CF) approximation to general relativity and the GW signal from the merger simulation is extracted up to quadrupole order.

Properties of general relativistic, irrotational binary neutron stars in close quasiequilibrium orbits: Polytropic equations of state

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 always terminated at an innermost orbit where a cusp (inner Lagrange point)appears at the inner edges of the stellar surface. The binaries with cusps are found to be dynamically unstable for n=0.5 and stable for n > 0.8 irrespective of M/R 0.13 irrespective of n, which indicates that the realistic binary neutron stars satisfy a necessary condition (q<1) for formation of a black hole before the merger. (vi) The specific angular momentum of any mass element in irrotational binary neutron stars at the innermost orbit appears to be too small to form a disk around black holes formed after the merger.

Binary neutron stars: Equilibrium models beyond spatial conformal flatness.

Equilibria of binary neutron stars in close circular orbits are computed numerically in a waveless formulation: the full Einstein-relativistic-Euler system is solved on an initial hypersurface to obtain an asymptotically flat form of the 4-metric and an extrinsic curvature whose time derivative vanishes in a comoving frame. Two independent numerical codes are developed, and solution sequences that model inspiraling binary neutron stars during the final several orbits are successfully computed. The binding energy of the system near its final orbit deviates from earlier results of third post-Newtonian and of spatially conformally flat calculations. The new solutions may serve as initial data for merger simulations and as members of quasiequilibrium sequences to generate gravitational-wave templates, and may improve estimates of the gravitational-wave cutoff frequency set by the last inspiral orbit.

Numerical method for binary black hole/neutron star initial data: Code test

A new numerical method to construct binary black hole/neutron star initial data is presented. The method uses three spherical coordinate patches; two of these are centered at the binary compact objects and cover a neighborhood of each object; the third patch extends to the asymptotic region. As in the Komatsu-Eriguchi-Hachisu method, nonlinear elliptic field equations are decomposed into a flat space Laplacian and a remaining nonlinear expression that serves in each iteration as an effective source. The equations are solved iteratively, integrating a Greens function against the effective source at each iteration. Detailed convergence tests for the essential part of the code are performed for a few types of selected Greens functions to treat different boundary conditions. Numerical computation of the gravitational potential of a fluid source, and a toy model for a binary black hole field, are carefully calibrated with the analytic solutions to examine accuracy and convergence of the new code. As an example of the application of the code, an initial data set for binary black holes in the Isenberg-Wilson-Mathews formulation is presented, in which the apparent horizons are located using a method described in Appendix 1.

Models of helically symmetric binary systems

Results from helically symmetric scalar-field models and first results from a convergent helically symmetric binary neutron-star code are reported here; these are models stationary in the rotating frame of a source with constant angular velocity Q. In the scalar-field models and the neutron-star code, helical symmetry leads to a system of mixed elliptic-hyperbolic character. The scalar-field models involve nonlinear terms of the form ψ 3 , (∇ψ) 2 and ψ□ψ that mimic nonlinear terms of the Einstein equation. Convergence is strikingly different for different signs of each nonlinear term; it is typically insensitive to the iterative method used, and it improves with an outer boundary in the near zone. In the neutron-star code, one has no control on the sign of the source, and convergence has been achieved only for an outer boundary less than ∼1 wavelength from the source or for a code that imposes helical symmetry only inside a near zone of that size. The inaccuracy of helically symmetric solutions with appropriate boundary conditions should be comparable to the inaccuracy of a waveless formalism that neglects gravitational waves, and the (near zone) solutions we obtain for waveless and helically symmetric BNS codes with the same boundary conditions nearly coincide.

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

Gravitational waves from binary neutron stars in quasiequilibrium circular orbits are computed using an approximate method which we propose in this paper. In the first step of this method, we prepare general relativistic irrotational binary neutron stars in a quasiequilibrium circular orbit, neglecting gravitational waves. We adopt the so-called conformal flatness approximation for a three-metric to obtain the quasiequilibrium states in this paper. In the second step, we compute gravitational waves, solving linear perturbation equations in the background spacetime of the quasiequilibrium states. Comparing numerical results with post Newtonian waveforms and luminosity of gravitational waves from two point masses in circular orbits, we demonstrate that this method can produce accurate waveforms and luminosity of gravitational waves. It is shown that the effects of tidal deformation of neutron stars and strong general relativistic gravity modify the post Newtonian results for compact binary neutron stars in close orbits. We indicate that the magnitude of a systematic error in quasiequilibrium states associated with the conformal flatness approximation is fairly large for close and compact binary neutron stars. Several formulations for improving the accuracy of quasiequilibrium states are proposed.

Stationary Structures of Irrotational Binary Systems: Models for Close Binary Systems of Compact Stars

We propose a new numerical method to calculate irrotational binary systems composed of compressible gaseous stars in Newtonian gravity. Assuming irrotationality, i.e., vanishing of the vorticity vector everywhere in the star in the inertial frame, we can introduce the velocity potential for the flow field. Using this velocity potential we can derive a set of basic equations for stationary states that consist of (1) the generalized Bernoulli equation, (2) the Poisson equation for the Newtonian gravitational potential, and (3) the equation for the velocity potential with the Neumann-type boundary condition. We succeeded in developing a new code to compute numerically exact solutions to these equations for the first time. Such irrotational configurations of binary systems are appropriate models for realistic neutron star binaries composed of inviscid gases, just prior to coalescence of two stars caused by emission of gravitational waves. Accuracies of our numerical solutions are so high that we can compute reliable models for fully deformed final stationary configurations and hence determine the inner most stable circular orbit of binary neutron star systems under the approximations of weak gravity and inviscid limit.

Gravitational wave emission by cataclysmic variables: numerical models of semi-detached binaries

Gravitational wave emission is considered to be the driving force for the evolution of shortperiod cataclysmic binary stars, making them a potential test for the validity of General Relativity. In spite of continuous refinements of the physical description, a 10% mismatch exists between the theoretical minimum period (Pturn ’ 70 min) and the short-period cut-off (Pmin ’ 80 min) observed in the period distribution for cataclysmic variable binaries. A possible explanation for this mismatch was associated with the use of the Roche model. We here present a systematic comparison between self-consistent, numerically constructed sequences of hydrostatic models of binary stars and Roche models of semi-detached binaries. On the basis of our approach, we also derive a value for the minimum period of cataclysmic variable binaries. The results obtained through the comparison indicate that the Roche model is indeed very good, with deviations from the numerical solution which are of a few percent at most. Our results therefore suggest that additional sources of angular momentum loss or alternative explanations need to be considered in order to justify the mismatch.

Modeling differential rotations of compact stars in equilibriums

Outcomes of numerical relativity simulations of massive core collapses or binary neutron star mergers with moderate masses suggest formations of rapidly and differentially rotating neutron stars. Subsequent fall back accretion may also amplify the degree of differential rotations. We propose new formulations for modeling differential rotations of those compact stars, and present selected solutions of differentially rotating, stationary, and axisymmetric compact stars in equilibriums. For the cases when rotating stars reach break-up velocities, the maximum masses of such rotating models are obtained.