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

A relativistic formalism for computation of irrotational binary stars in quasi equilibrium states

We present relativistic hydrostatic equations for obtaining irrotational binary neutron stars in quasi equilibrium states in 3+1 formalism. Equations derived here are different from those previously given by Bonazzola, Gourgoulhon, and Marck, and have a simpler and more tractable form for computation in numerical relativity. We also present hydrostatic equations for computation of equilibrium irrotational binary stars in first post-Newtonian order.Comment: 5 pages, corrected eqs.(2.10), (2.11) and (3.1

Gravitational waves from axisymmetrically oscillating neutron stars in general relativistic simulations

Gravitational waves from oscillating neutron stars in axial symmetry are studied performing numerical simulations in full general relativity. Neutron stars are modeled by a polytropic equation of state for simplicity. A gauge-invariant wave extraction method as well as a quadrupole formula are adopted for computation of gravitational waves. It is found that the gauge-invariant variables systematically contain numerical errors generated near the outer boundaries in the present axisymmetric computation. We clarify their origin, and illustrate it possible to eliminate the dominant part of the systematic errors. The best corrected waveforms for oscillating and rotating stars currently contain errors of magnitude $\sim 10^{-3}$ in the local wave zone. Comparing the waveforms obtained by the gauge-invariant technique with those by the quadrupole formula, it is shown that the quadrupole formula yields approximate gravitational waveforms besides a systematic underestimation of the amplitude of $O(M/R)$ where $M$ and $R$ denote the mass and the radius of neutron stars. However, the wave phase and modulation of the amplitude can be computed accurately. This indicates that the quadrupole formula is a useful tool for studying gravitational waves from rotating stellar core collapse to a neutron star in fully general relativistic simulations. Properties of the gravitational waveforms from the oscillating and rigidly rotating neutron stars are also addressed paying attention to the oscillation associated with fundamental modes

Possible explanation for star-crushing effect in binary neutron star simulations

A possible explanation is suggested for the controversial star-crushing effect seen in numerical simulations of inspiraling neutron star binaries by Wilson, Mathews and Marronetti (WMM). An apparently incorrect definition of momentum density in the momentum constraint equation used by WMM gives rise to a post-1-Newtonian error in the approximation scheme. We show by means of an analytic, post-1-Newtonian calculation that this error causes an increase of the stars' central densities which is of the order of several percent when the stars are separated by a few stellar radii, in agreement with what is seen in the simulations.Comment: 4 pages, 1 figure, uses revetx macros, minor revision

Various features of quasiequilibrium sequences of binary neutron stars in general relativity

Quasiequilibrium sequences of binary neutron stars are numerically calculated in the framework of the Isenberg-Wilson-Mathews (IWM) approximation of general relativity. The results are presented for both rotation states of synchronized spins and irrotational motion, the latter being considered as the realistic one for binary neutron stars just prior to the merger. We assume a polytropic equation of state and compute several evolutionary sequences of binary systems composed of different-mass stars as well as identical-mass stars with adiabatic indices gamma=2.5, 2.25, 2, and 1.8. From our results, we propose as a conjecture that if the turning point of binding energy (and total angular momentum) locating the innermost stable circular orbit (ISCO) is found in Newtonian gravity for some value of the adiabatic index gamma_0, that of the ADM mass (and total angular momentum) should exist in the IWM approximation of general relativity for the same value of the adiabatic index.Comment: Text improved, some figures changed or deleted, new table, 38 pages, 31 figures, accepted for publication in Phys. Rev.

Coupled charge and valley excitations in graphene quantum Hall ferromagnets

Graphene is a two-dimensional carbon material with a honeycomb lattice and Dirac-type low-energy spectrum. In a strong magnetic field, where Coulomb interactions dominate against disorder broadening, quantum Hall ferromagnetic states realize at integer fillings. Extending the quantum Hall ferromagnetism to the fractional filling case of massless Dirac fermions, we study the elementally charge excitations which couple with the valley degrees of freedom (so-called valley skyrmions). With the use of the density matrix renomalization group (DMRG) method, the excitation gaps are calculated and extrapolated to the thermodynamic limit. These results exhibit numerical evidences and criterions of the skyrmion excitations in graphene.Comment: 5 pages, 5 figure

Merger of binary neutron stars of unequal mass in full general relativity

We present results of three dimensional numerical simulations of the merger of unequal-mass binary neutron stars in full general relativity. A $\Gamma$-law equation of state $P=(\Gamma-1)\rho\epsilon$ is adopted, where $P$, $\rho$, \varep, and $\Gamma$ are the pressure, rest mass density, specific internal energy, and the adiabatic constant, respectively. We take $\Gamma=2$ and the baryon rest-mass ratio $Q_M$ to be in the range 0.85--1. The typical grid size is $(633,633,317)$ for $(x,y,z)$ . We improve several implementations since the latest work. In the present code, the radiation reaction of gravitational waves is taken into account with a good accuracy. This fact enables us to follow the coalescence all the way from the late inspiral phase through the merger phase for which the transition is triggered by the radiation reaction. It is found that if the total rest-mass of the system is more than $\sim 1.7$ times of the maximum allowed rest-mass of spherical neutron stars, a black hole is formed after the merger irrespective of the mass ratios. The gravitational waveforms and outcomes in the merger of unequal-mass binaries are compared with those in equal-mass binaries. It is found that the disk mass around the so formed black holes increases with decreasing rest-mass ratios and decreases with increasing compactness of neutron stars. The merger process and the gravitational waveforms also depend strongly on the rest-mass ratios even for the range $Q_M= 0.85$--1.Comment: 32 pages, PRD68 to be publishe

Temporal 1/f^\alpha Fluctuations from Fractal Magnetic Fields in Black Hole Accretion Flow

Rapid fluctuation with a frequency dependence of $1/f^{\alpha}$ (with $\alpha \simeq 1 - 2$) is characteristic of radiation from black-hole objects. Its origin remains poorly understood. We examine the three-dimensional magnetohydrodynamical (MHD) simulation data, finding that a magnetized accretion disk exhibits both $1/f^\alpha$ fluctuation (with $\alpha \simeq 2$) and a fractal magnetic structure (with the fractal dimension of $D \sim 1.9$). The fractal field configuration leads reconnection events with a variety of released energy and of duration, thereby producing $1/f^\alpha$ fluctuations.Comment: 5 pages, 4 figures. Accepted for publication in PASJ Letters, vol. 52 No.1 (Feb 2000