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

Axisymmetric general relativistic hydrodynamics: Long-term evolution of neutron stars and stellar collapse to neutron stars and black holes

We report a new implementation for axisymmetric simulation in full general relativity. In this implementation, the Einstein equations are solved using the Nakamura-Shibata formulation with the so-called cartoon method to impose an axisymmetric boundary condition, and the general relativistic hydrodynamic equations are solved using a high-resolution shock-capturing scheme based on an approximate Riemann solver. As tests, we performed the following simulations: (i) long-term evolution of non-rotating and rapidly rotating neutron stars, (ii) long-term evolution of neutron stars of a high-amplitude damping oscillation accompanied with shock formation, (iii) collapse of unstable neutron stars to black holes, and (iv) stellar collapses to neutron stars. The tests (i)--(iii) were carried out with the $\Gamma$-law equation of state, and the test (iv) with a more realistic parametric equation of state for high-density matter. We found that this new implementation works very well: It is possible to perform the simulations for stable neutron stars for more than 10 dynamical time scales, to capture strong shocks formed at stellar core collapses, and to accurately compute the mass of black holes formed after the collapse and subsequent accretion. In conclusion, this implementation is robust enough to apply to astrophysical problems such as stellar core collapse of massive stars to a neutron star and black hole, phase transition of a neutron star to a high-density star, and accretion-induced collapse of a neutron star to a black hole. The result for the first simulation of stellar core collapse to a neutron star started from a realistic initial condition is also presented.Comment: 28 pages, to appear in PRD 67, 0440XX (2003

Robustness of a high-resolution central scheme for hydrodynamic simulations in full general relativity

A recent paper by Lucas-Serrano et al. indicates that a high-resolution central (HRC) scheme is robust enough to yield accurate hydrodynamical simulations of special relativistic flows in the presence of ultrarelativistic speeds and strong shock waves. In this paper we apply this scheme in full general relativity (involving {\it dynamical} spacetimes), and assess its suitability by performing test simulations for oscillations of rapidly rotating neutron stars and merger of binary neutron stars. It is demonstrated that this HRC scheme can yield results as accurate as those by the so-called high-resolution shock-capturing (HRSC) schemes based upon Riemann solvers. Furthermore, the adopted HRC scheme has increased computational efficiency as it avoids the costly solution of Riemann problems and has practical advantages in the modeling of neutron star spacetimes. Namely, it allows simulations with stiff equations of state by successfully dealing with very low-density unphysical atmospheres. These facts not only suggest that such a HRC scheme may be a desirable tool for hydrodynamical simulations in general relativity, but also open the possibility to perform accurate magnetohydrodynamical simulations in curved dynamic spacetimes.Comment: 4 pages, to be published in Phys. Rev. D (brief report

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

Pulse Control of Decoherence in a Qubit Coupled with a Quantum Environment

We study the time evolution of a qubit linearly coupled with a quantum environment under a sequence of short pi pulses. Our attention is focused on the case where qubit-environment interactions induce the decoherence with population decay. We assume that the environment consists of a set of bosonic excitations. The time evolution of the reduced density matrix for the qubit is calculated in the presence of periodic short pi pulses. We confirm that the decoherence is suppressed if the pulse interval is shorter than the correlation time for qubit-environment interactions.Comment: 5 pages, 2figure

New criterion for direct black hole formation in rapidly rotating stellar collapse

We study gravitational collapse of rapidly rotating relativistic polytropes of the adiabatic index $\Gamma = 1.5$ and 2, in which the spin parameter $q \equiv J/M^{2} > 1$ where $J$ and $M$ are total angular momentum and gravitational mass, in full general relativity. First, analyzing initial distributions of the mass and the spin parameter inside stars, we predict the final outcome after the collapse. Then, we perform fully general relativistic simulations on assumption of axial and equatorial symmetries and confirm our predictions. As a result of simulations, we find that in contrast with the previous belief, even for stars with $q > 1$, the collapse proceeds to form a seed black hole at central region, and the seed black hole subsequently grows as the ambient fluids accrete onto it. We also find that growth of angular momentum and mass of the seed black hole can be approximately determined from the initial profiles of the density and the specific angular momentum. We define an effective spin parameter at the central region of the stars, $q_{c}$, and propose a new criterion for black hole formation as q_{c} \alt 1. Plausible reasons for the discrepancy between our and previous results are clarified.Comment: submitted to PR

Phonon-Metamorphosis in Ferromagnetic Manganite Films: Probing the Evolution of an Inhomogeneous State

The analysis of phonon anomalies provides valuable information about the cooperative dynamics of lattice, spin and charge degrees of freedom. Significant is the anomalous temperature dependence of the external modes observed in La$_{2/3}$Sr$_{1/3}$MnO$_{3}$ (LSMO) films. The two external modes merge close to the ferromagnetic to paramagnetic transition at $T_C$ and, moreover, two new modes evolve in this temperature range with strong resonances at slightly higher frequencies. We propose that this observed phonon metamorphosis probes the inhomogeneous Jahn-Teller distortion, manifest on the temperature scale $T_C$. The analysis is based on the first observation of all eight phonon modes in the metallic phase of LSMO and on susceptibility measurements which identify a Griffiths-like phase.Comment: 4 pages, 4 figure

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.

Simulation of merging binary neutron stars in full general relativity: Γ=2\Gamma=2 case

We have performed 3D numerical simulations for merger of equal mass binary neutron stars in full general relativity. 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 the adiabatic constant with $\Gamma=2$. As initial conditions, we adopt models of corotational and irrotational binary neutron stars in a quasi-equilibrium state which are obtained using the conformal flatness approximation for the three geometry as well as an assumption that a helicoidal Killing vector exists. In this paper, we pay particular attention to the final product of the coalescence. We find that the final product depends sensitively on the initial compactness parameter of the neutron stars : In a merger between sufficiently compact neutron stars, a black hole is formed in a dynamical timescale. As the compactness is decreased, the formation timescale becomes longer and longer. It is also found that a differentially rotating massive neutron star is formed instead of a black hole for less compact binary cases, in which the rest mass of each star is less than 70-80% of the maximum allowed mass of a spherical star. In the case of black hole formation, we roughly evaluate the mass of the disk around the black hole. For the merger of corotational binaries, a disk of mass $\sim 0.05-0.1M_*$ may be formed, where M_* is the total rest mass of the system. On the other hand, for the merger of irrotational binaries, the disk mass appears to be very small : < 0.01M_*.Comment: 27 pages, to appear in Phys. Rev.

Measuring black-hole parameters and testing general relativity using gravitational-wave data from space-based interferometers

Among the expected sources of gravitational waves for the Laser Interferometer Space Antenna (LISA) is the capture of solar-mass compact stars by massive black holes residing in galactic centers. We construct a simple model for such a capture, in which the compact star moves freely on a circular orbit in the equatorial plane of the massive black hole. We consider the gravitational waves emitted during the late stages of orbital evolution, shortly before the orbiting mass reaches the innermost stable circular orbit. We construct a simple model for the gravitational-wave signal, in which the phasing of the waves plays the dominant role. The signal's behavior depends on a number of parameters, including $\mu$, the mass of the orbiting star, $M$, the mass of the central black hole, and $J$, the black hole's angular momentum. We calculate, using our simplified model, and in the limit of large signal-to-noise ratio, the accuracy with which these quantities can be estimated during a gravitational-wave measurement. Our simplified model also suggests a method for experimentally testing the strong-field predictions of general relativity.Comment: ReVTeX, 16 pages, 5 postscript figure