Relativistic Stellar Pulsations With Near-Zone Boundary Conditions

A new method is presented here for evaluating approximately the pulsation modes of relativistic stellar models. This approximation relies on the fact that gravitational radiation influences these modes only on timescales that are much longer than the basic hydrodynamic timescale of the system. This makes it possible to impose the boundary conditions on the gravitational potentials at the surface of the star rather than in the asymptotic wave zone of the gravitational field. This approximation is tested here by predicting the frequencies of the outgoing non-radial hydrodynamic modes of non-rotating stars. The real parts of the frequencies are determined with an accuracy that is better than our knowledge of the exact frequencies (about 0.01%) except in the most relativistic models where it decreases to about 0.1%. The imaginary parts of the frequencies are determined with an accuracy of approximately M/R, where M is the mass and R is the radius of the star in question.Comment: 10 pages (REVTeX 3.1), 5 figs., 1 table, fixed minor typos, published in Phys. Rev. D 56, 2118 (1997

Solving Einstein's Equations With Dual Coordinate Frames

A method is introduced for solving Einstein's equations using two distinct coordinate systems. The coordinate basis vectors associated with one system are used to project out components of the metric and other fields, in analogy with the way fields are projected onto an orthonormal tetrad basis. These field components are then determined as functions of a second independent coordinate system. The transformation to the second coordinate system can be thought of as a mapping from the original ``inertial'' coordinate system to the computational domain. This dual-coordinate method is used to perform stable numerical evolutions of a black-hole spacetime using the generalized harmonic form of Einstein's equations in coordinates that rotate with respect to the inertial frame at infinity; such evolutions are found to be generically unstable using a single rotating coordinate frame. The dual-coordinate method is also used here to evolve binary black-hole spacetimes for several orbits. The great flexibility of this method allows comoving coordinates to be adjusted with a feedback control system that keeps the excision boundaries of the holes within their respective apparent horizons.Comment: Updated to agree with published versio

On unbounded bodies with finite mass: asymptotic behaviour

There is introduced a class of barotropic equations of state (EOS) which become polytropic of index $n = 5$ at low pressure. One then studies asymptotically flat solutions of the static Einstein equations coupled to perfect fluids having such an EOS. It is shown that such solutions, in the same manner as the vacuum ones, are conformally smooth or analytic at infinity, when the EOS is smooth or analytic, respectively.Comment: 6 page

Gravitational Radiation Instability in Hot Young Neutron Stars

We show that gravitational radiation drives an instability in hot young rapidly rotating neutron stars. This instability occurs primarily in the l=2 r-mode and will carry away most of the angular momentum of a rapidly rotating star by gravitational radiation. On the timescale needed to cool a young neutron star to about T=10^9 K (about one year) this instability can reduce the rotation rate of a rapidly rotating star to about 0.076\Omega_K, where \Omega_K is the Keplerian angular velocity where mass shedding occurs. In older colder neutron stars this instability is suppressed by viscous effects, allowing older stars to be spun up by accretion to larger angular velocities.Comment: 4 Pages, 2 Figure

R-Modes in Superfluid Neutron Stars

The analogs of r-modes in superfluid neutron stars are studied here. These modes, which are governed primarily by the Coriolis force, are identical to their ordinary-fluid counterparts at the lowest order in the small angular-velocity expansion used here. The equations that determine the next order terms are derived and solved numerically for fairly realistic superfluid neutron-star models. The damping of these modes by superfluid ``mutual friction'' (which vanishes at the lowest order in this expansion) is found to have a characteristic time-scale of about 10^4 s for the m=2 r-mode in a ``typical'' superfluid neutron-star model. This time-scale is far too long to allow mutual friction to suppress the recently discovered gravitational radiation driven instability in the r-modes. However, the strength of the mutual friction damping depends very sensitively on the details of the neutron-star core superfluid. A small fraction of the presently acceptable range of superfluid models have characteristic mutual friction damping times that are short enough (i.e. shorter than about 5 s) to suppress the gravitational radiation driven instability completely.Comment: 15 pages, 8 figure

Second-order rotational effects on the r-modes of neutron stars

Techniques are developed here for evaluating the r-modes of rotating neutron stars through second order in the angular velocity of the star. Second-order corrections to the frequencies and eigenfunctions for these modes are evaluated for neutron star models. The second-order eigenfunctions for these modes are determined here by solving an unusual inhomogeneous hyperbolic boundary-value problem. The numerical techniques developed to solve this unusual problem are somewhat non-standard and may well be of interest beyond the particular application here. The bulk-viscosity coupling to the r-modes, which appears first at second order, is evaluated. The bulk-viscosity timescales are found here to be longer than previous estimates for normal neutron stars, but shorter than previous estimates for strange stars. These new timescales do not substantially affect the current picture of the gravitational radiation driven instability of the r-modes either for neutron stars or for strange stars.Comment: 13 pages, 5 figures, revte

Generalized r-Modes of the Maclaurin Spheroids

Analytical solutions are presented for a class of generalized r-modes of rigidly rotating uniform density stars---the Maclaurin spheroids---with arbitrary values of the angular velocity. Our analysis is based on the work of Bryan; however, we derive the solutions using slightly different coordinates that give purely real representations of the r-modes. The class of generalized r-modes is much larger than the previously studied `classical' r-modes. In particular, for each l and m we find l-m (or l-1 for the m=0 case) distinct r-modes. Many of these previously unstudied r-modes (about 30% of those examined) are subject to a secular instability driven by gravitational radiation. The eigenfunctions of the `classical' r-modes, the l=m+1 case here, are found to have particularly simple analytical representations. These r-modes provide an interesting mathematical example of solutions to a hyperbolic eigenvalue problem.Comment: 12 pages, 3 figures; minor changes and additions as will appear in the version to be published in Physical Review D, January 199

Time-Independent Gravitational Fields

This article reviews, from a global point of view, rigorous results on time independent spacetimes. Throughout attention is confined to isolated bodies at rest or in uniform rotation in an otherwise empty universe. The discussion starts from first principles and is, as much as possible, self-contained.Comment: 47 pages, LaTeX, uses Springer cl2emult styl

Reducing orbital eccentricity in binary black hole simulations

Binary black hole simulations starting from quasi-circular (i.e., zero radial velocity) initial data have orbits with small but non-zero orbital eccentricities. In this paper the quasi-equilibrium initial-data method is extended to allow non-zero radial velocities to be specified in binary black hole initial data. New low-eccentricity initial data are obtained by adjusting the orbital frequency and radial velocities to minimize the orbital eccentricity, and the resulting ($\sim 5$ orbit) evolutions are compared with those of quasi-circular initial data. Evolutions of the quasi-circular data clearly show eccentric orbits, with eccentricity that decays over time. The precise decay rate depends on the definition of eccentricity; if defined in terms of variations in the orbital frequency, the decay rate agrees well with the prediction of Peters (1964). The gravitational waveforms, which contain $\sim 8$ cycles in the dominant l=m=2 mode, are largely unaffected by the eccentricity of the quasi-circular initial data. The overlap between the dominant mode in the quasi-circular evolution and the same mode in the low-eccentricity evolution is about 0.99.Comment: 27 pages, 9 figures; various minor clarifications; accepted to the "New Frontiers" special issue of CQ