Merger of binary neutron stars with realistic equations of state in full general relativity
- Graduate School of Arts and Sciences, University of Tokyo, Komaba, Meguro, Tokyo 153-8902 (Japan)
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801 (United States)
- Astrophysical Sector, SISSA, via Beirut 2/4, Trieste 34013 (Italy)
We present numerical results of three-dimensional simulations for the merger of binary neutron stars in full general relativity. Hybrid equations of state are adopted to mimic realistic nuclear equations of state. In this approach, we divide the equations of state into two parts as P=P{sub cold}+P{sub th}. P{sub cold} is the cold part for which we assign a fitting formula for realistic equations of state of cold nuclear matter slightly modifying the formula developed by Haensel and Potekhin. We adopt the SLy and FPS equations of state for which the maximum allowed Arnowitt-Deser-Misner (ADM) mass of cold and spherical neutron stars is {approx_equal}2.04M{sub {center_dot}} and 1.80M{sub {center_dot}}, respectively. P{sub th} denotes the thermal part which is written as P{sub th}=({gamma}{sub th}-1){rho}({epsilon}-{epsilon}{sub cold}), where {rho}, {epsilon}, {epsilon}{sub cold}, and {gamma}{sub th} are the baryon rest-mass density, total specific internal energy, specific internal energy of the cold part, and the adiabatic constant, respectively. Simulations are performed for binary neutron stars of the total ADM mass in the range between 2.4M{sub {center_dot}} and 2.8M{sub {center_dot}} with the rest-mass ratio Q{sub M} to be in the range 0.9 < or approx. Q{sub M}{<=}1. It is found that if the total ADM mass of the system is larger than a threshold M{sub thr}, a black hole is promptly formed in the merger irrespective of the mass ratios. In the other case, the outcome is a hypermassive neutron star of a large ellipticity, which results from the large adiabatic index of the realistic equations of state adopted. The value of M{sub thr} depends on the equation of state: M{sub thr}{approx}2.7M{sub {center_dot}} and {approx}2.5M{sub {center_dot}} for the SLy and FPS equations of state, respectively. Gravitational waves are computed in terms of a gauge-invariant wave extraction technique. In the formation of the hypermassive neutron star, quasiperiodic gravitational waves of a large amplitude and of frequency between 3 and 4 kHz are emitted. The estimated emission time scale is < or approx.100 ms, after which the hypermassive neutron star collapses to a black hole. Because of the long emission time, the effective amplitude may be large enough to be detected by advanced laser interferometric gravitational wave detectors if the distance to the source is smaller than {approx}100 Mpc. Thermal properties of the outcome formed after the merger are also analyzed to approximately estimate the neutrino emission energy.
- OSTI ID:
- 20709092
- Journal Information:
- Physical Review. D, Particles Fields, Vol. 71, Issue 8; Other Information: DOI: 10.1103/PhysRevD.71.084021; (c) 2005 The American Physical Society; Country of input: International Atomic Energy Agency (IAEA); ISSN 0556-2821
- Country of Publication:
- United States
- Language:
- English
Similar Records
Long-term general relativistic simulation of binary neutron stars collapsing to a black hole
Merger of black hole and neutron star in general relativity: Tidal disruption, torus mass, and gravitational waves
Related Subjects
AMPLITUDES
BARYONS
BINARY STARS
BLACK HOLES
COMPUTERIZED SIMULATION
EQUATIONS OF STATE
GAUGE INVARIANCE
GENERAL RELATIVITY THEORY
GRAVITATIONAL WAVE DETECTORS
GRAVITATIONAL WAVES
KHZ RANGE
NEUTRINOS
NEUTRON STARS
NUCLEAR MATTER
REST MASS
STAR EVOLUTION
THERMODYNAMIC PROPERTIES
THREE-DIMENSIONAL CALCULATIONS