No abstract

We employ a semianalytic technique, based on a modified WKB approach, to determine the complex normal-mode frequencies of Schwarzschild black holes. It yields a simple analytic formula that gives the real and imaginary parts of the frequency in terms of the parameters of the black hole and of the field whose perturbation is under study, and in terms of the quantity (n ++), where n =0, k 1 , k2, . . . , and labels the fundamental mode, first overtone mode, and so on. In the case of the fundamental gravitational normal modes of the Schwarzschild black hole, the WKB estimates agree with numerical results to better than 0.13% in the real part of the frequency and 0.22% in the imaginary part. The agreement for both the real and imaginary parts of the low overtones is better than 0.5%. The relative agreement improves with increasing angular harmonic.

This is the second in a series of papers on the construction and validation of a three-dimensional code for the solution of the coupled system of the Einstein equations and of the general relativistic hydrodynamic equations, and on the application of this code to problems in general relativistic astrophysics. In particular, we report on the accuracy of our code in the long-term dynamical evolution of relativistic stars and on some new physics results obtained in the process of code testing. The following aspects of our code have been validated: the generation of initial data representing perturbed general relativistic polytropic models ͑both rotating and nonrotating͒, the long-term evolution of relativistic stellar models, and the coupling of our evolution code to analysis modules providing, for instance, the detection of apparent horizons or the extraction of gravitational waveforms. The tests involve single nonrotating stars in stable equilibrium, nonrotating stars undergoing radial and quadrupolar oscillations, nonrotating stars on the unstable branch of the equilibrium configurations migrating to the stable branch, nonrotating stars undergoing gravitational collapse to a black hole, and rapidly rotating stars in stable equilibrium and undergoing quasiradial oscillations. We have carried out evolutions in full general relativity and compared the results to those obtained either with perturbation techniques, or with lower dimensional numerical codes, or in the Cowling approximation ͑in which all the perturbations of the spacetime are neglected͒. In all cases an excellent agreement has been found. The numerical evolutions have been carried out using different types of polytropic equations of state using either the rest-mass density only, or the rest-mass density and the internal energy as independent variables. New variants of the spacetime evolution and new high resolution shock capturing treatments based on Riemann solvers and slope limiters have been implemented and the results compared with those obtained from previous methods. In particular, we have found the ''monotonized central differencing'' limiter to be particularly effective in evolving the relativistic stellar models considered. Finally, we have obtained the first eigenfrequencies of rotating stars in full general relativity and rapid rotation. A long standing problem, such frequencies have not been obtained by other methods. Overall, and to the best of our knowledge, the results presented in this paper represent the most accurate long-term three-dimensional evolutions of relativistic stars available to date.

No abstract

We use a post-Newtonian diagnostic tool to examine numerically generated quasiequilibrium initial data sets for non-spinning double neutron star and neutron star-black hole binary systems. The PN equations include the effects of tidal interactions, parametrized by the compactness of the neutron stars and by suitable values of "apsidal" constants, which measure the degree of distortion of stars subjected to tidal forces. We find that the post-Newtonian diagnostic agrees well with the double neutron star initial data, typically to better than half a percent except where tidal distortions are becoming extreme. We show that the differences could be interpreted as representing small residual eccentricity in the initial orbits. In comparing the diagnostic with preliminary numerical data on neutron star-black hole binaries, we find less agreement.