Incompressible fluid binary systems with internal flows -- models of close binary neutron star systems with spin

In this paper we have examined numerically exact configurations of close binary systems composed of incompressible fluids with internal flows. Component stars of binary systems are assumed to be circularly orbiting each other but rotating nonsynchronously with the orbital motion, i.e. stars in binary systems have steady motions seen from a rotational frame of reference. We have computed several equilibrium sequences by taking fully into account the tidal effect of Newtonian gravity without approximation. We consider two binary systems consisting of either 1) a point mass and a fluid star or 2) a fluid star and a fluid star. Each of them corresponds to generalization of the Roche--Riemann or the Darwin--Riemann problem, respectively. Our code can be applied to various types of incompressible binary systems with various mass ratios and spin as long as the vorticity is constant. We compare these equilibrium sequences of binaries with approximate solutions which were studied extensively by Lai, Rasio and Shapiro (LRS) as models for neutron star--neutron star (NS--NS) binary systems or black hole--neutron star (BH--NS) binary systems. Our results coincide qualitatively with those of LRS but are different from theirs for configurations with small separations. For these binary systems, our sequences show that dynamical or secular instability sets in as the separation decreases. The quantitative errors of the ellipsoidal approximation amount to 2\sim 10% for configurations near the instability point. Compared to the results of LRS, the separation of the stars at the point where the instability sets in is larger and correspondingly the orbital frequency at the critical state is smaller for our models.Comment: 13 pages, 6 bitmapped ps files, to appear in MNRA

Circular solution of two unequal mass particles in post-Minkowski approximation

A Fokker action for post-Minkowski approximation with the first post-Newtonian correction is introduced in our previous paper, and a solution for the helically symmetric circular orbit is obtained. We present supplemental results for the circular solution of two unequal mass point-particles. Circular solutions for selected mass ratios are found numerically, and analytic formulas in the extreme mass ratio limit are derived. The leading terms of the analytic formulas agree with the first post-Newtonian formulas in this limit.Comment: 4 pages, 4 figures, 4/27/0

New code for equilibriums and quasiequilibrium initial data of compact objects. II. Convergence tests and comparisons of binary black hole initial data

COCAL is a code for computing equilibriums or quasiequilibrium initial data of single or binary compact objects based on finite difference methods. We present the results of supplementary convergence tests of COCAL code using time symmetric binary black hole data (Brill-Lindquist solution). Then, we compare the initial data of binary black holes on the conformally flat spatial slice obtained from COCAL and KADATH, where KADATH is a library for solving a wide class of problems in theoretical physics including relativistic compact objects with spectral methods. Data calculated from the two codes converge nicely towards each other, for close as well as largely separated circular orbits of binary black holes. Finally, as an example, a sequence of equal mass binary black hole initial data with corotating spins is calculated and compared with data in the literature.Comment: 9 pages, PRD in pres

Gravitational Waves from the Merger of Binary Neutron Stars in a Fully General Relativistic Simulation

We performed 3D numerical simulations of the merger of equal-mass binary neutron stars in full general relativity using a new large scale supercomputer. We take the typical grid size as (505,505,253) for (x,y,z) and the maximum grid size as (633,633,317). These grid numbers enable us to put the outer boundaries of the computational domain near the local wave zone and hence to calculate gravitational waveforms of good accuracy (within $\sim 10$ error) for the first time. To model neutron stars, we adopt a $\Gamma$-law equation of state in the form $P=(\Gamma-1)\rho\epsilon$, where P, $\rho$, \varep and $\Gamma$ are the pressure, rest mass density, specific internal energy, and adiabatic constant. It is found that gravitational waves in the merger stage have characteristic features that reflect the formed objects. In the case that a massive, transient neutron star is formed, its quasi-periodic oscillations are excited for a long duration, and this property is reflected clearly by the quasi-periodic nature of waveforms and the energy luminosity. In the case of black hole formation, the waveform and energy luminosity are likely damped after a short merger stage. However, a quasi-periodic oscillation can still be seen for a certain duration, because an oscillating transient massive object is formed during the merger. This duration depends strongly on the initial compactness of neutron stars and is reflected in the Fourier spectrum of gravitational waves. To confirm our results and to calibrate the accuracy of gravitational waveforms, we carried out a wide variety of test simulations, changing the resolution and size of the computational domain.Comment: 40 pages; pubslihed in Prog. Theor. Phys. 107 (2002), 26

Thermodynamics of magnetized binary compact objects

Binary systems of compact objects with electromagnetic field are modeled by helically symmetric Einstein-Maxwell spacetimes with charged and magnetized perfect fluids. Previously derived thermodynamic laws for helically-symmetric perfect-fluid spacetimes are extended to include the electromagnetic fields, and electric currents and charges; the first law is written as a relation between the change in the asymptotic Noether charge \dl Q and the changes in the area and electric charge of black holes, and in the vorticity, baryon rest mass, entropy, charge and magnetic flux of the magnetized fluid. Using the conservation laws of the circulation of magnetized flow found by Bekenstein and Oron for the ideal magnetohydrodynamic (MHD) fluid, and also for the flow with zero conducting current, we show that, for nearby equilibria that conserve the quantities mentioned above, the relation \dl Q=0 is satisfied. We also discuss a formulation for computing numerical solutions of magnetized binary compact objects in equilibrium with emphasis on a first integral of the ideal MHD-Euler equation.Comment: 21 pages, to appear in PR

Thermodynamics of binary black holes and neutron stars

We consider compact binary systems, modeled in general relativity as vacuum or perfect-fluid spacetimes with a helical Killing vector k^\alpha, heuristically, the generator of time-translations in a corotating frame. Systems that are stationary in this sense are not asymptotically flat, but have asymptotic behavior corresponding to equal amounts of ingoing and outgoing radiation. For black-hole binaries, a rigidity theorem implies that the Killing vector lies along the horizon's generators, and from this one can deduce the zeroth law (constant surface gravity of the horizon). Remarkably, although the mass and angular momentum of such a system are not defined, there is an exact first law, relating the change in the asymptotic Noether charge to the changes in the vorticity, baryon mass, and entropy of the fluid, and in the area of black holes. Binary systems with M\Omega small have an approximate asymptopia in which one can write the first law in terms of the asymptotic mass and angular momentum. Asymptotic flatness is precise in two classes of solutions used to model binary systems: spacetimes satisfying the post-Newtonian equations, and solutions to a modified set of field equations that have a spatially conformally flat metric. (The spatial conformal flatness formalism with helical symmetry, however, is consistent with maximal slicing only if replaces the extrinsic curvature in the field equations by an artificially tracefree expression in terms of the shift vector.) For these spacetimes, nearby equilibria whose stars have the same vorticity obey the relation \delta M = \Omega \delta J, from which one can obtain a turning point criterion that governs the stability of orbits.Comment: 26 pages, revised version, modified appendix A, to appear in PR

Darwin-Riemann Problems in Newtonian Gravity

In this paper, we have reviewed the present status of the theory of equilibrium configurations of compact binary star systems in Newtonian gravity. Evolutionary processes of compact binary star systems due to gravitational wave emission can be divided into three stages according to the time scales and configurations. The evolution is quasi-stationary until a merging process starts, since the time scale of the orbital change due to gravitational wave emission is longer than the orbital period. In this stage, equilibrium sequences can be applied to evolution of compact binary star systems. Along the equilibrium sequences, there appear several critical states where some instability sets in or configuration changes drastically. We have discussed relations among these critical points and have stressed the importance of the mass overflow as well as the dynamical instability of orbital motions. Concerning the equilibrium sequences of binary star systems, we have summarized classical results of incompressible ellipsoidal configurations. Recent results of compressible binary star systems obtained by the ellipsoidal approximation and by numerical computations have been shown and discussed. It is important to note that numerical computational solutions to {\it exact equations} show that compressibility may lead realistic neutron star binary systems to mass overflows instead of dynamical disruptions for a wide range of parameters.Comment: 17 pages, 10 figures, PTPTeX style files are include