Improved outer boundary conditions for Einstein's field equations

In a recent article, we constructed a hierarchy B_L of outer boundary conditions for Einstein's field equations with the property that, for a spherical outer boundary, it is perfectly absorbing for linearized gravitational radiation up to a given angular momentum number L. In this article, we generalize B_2 so that it can be applied to fairly general foliations of spacetime by space-like hypersurfaces and general outer boundary shapes and further, we improve B_2 in two steps: (i) we give a local boundary condition C_2 which is perfectly absorbing including first order contributions in 2M/R of curvature corrections for quadrupolar waves (where M is the mass of the spacetime and R is a typical radius of the outer boundary) and which significantly reduces spurious reflections due to backscatter, and (ii) we give a non-local boundary condition D_2 which is exact when first order corrections in 2M/R for both curvature and backscatter are considered, for quadrupolar radiation.Comment: accepted Class. Quant. Grav. numerical relativity special issue; 17 pages and 1 figur

Testing outer boundary treatments for the Einstein equations

Various methods of treating outer boundaries in numerical relativity are compared using a simple test problem: a Schwarzschild black hole with an outgoing gravitational wave perturbation. Numerical solutions computed using different boundary treatments are compared to a `reference' numerical solution obtained by placing the outer boundary at a very large radius. For each boundary treatment, the full solutions including constraint violations and extracted gravitational waves are compared to those of the reference solution, thereby assessing the reflections caused by the artificial boundary. These tests use a first-order generalized harmonic formulation of the Einstein equations. Constraint-preserving boundary conditions for this system are reviewed, and an improved boundary condition on the gauge degrees of freedom is presented. Alternate boundary conditions evaluated here include freezing the incoming characteristic fields, Sommerfeld boundary conditions, and the constraint-preserving boundary conditions of Kreiss and Winicour. Rather different approaches to boundary treatments, such as sponge layers and spatial compactification, are also tested. Overall the best treatment found here combines boundary conditions that preserve the constraints, freeze the Newman-Penrose scalar Psi_0, and control gauge reflections.Comment: Modified to agree with version accepted for publication in Class. Quantum Gra

Stable radiation-controlling boundary conditions for the generalized harmonic Einstein equations

This paper is concerned with the initial-boundary value problem for the Einstein equations in a first-order generalized harmonic formulation. We impose boundary conditions that preserve the constraints and control the incoming gravitational radiation by prescribing data for the incoming fields of the Weyl tensor. High-frequency perturbations about any given spacetime (including a shift vector with subluminal normal component) are analyzed using the Fourier-Laplace technique. We show that the system is boundary-stable. In addition, we develop a criterion that can be used to detect weak instabilities with polynomial time dependence, and we show that our system does not suffer from such instabilities. A numerical robust stability test supports our claim that the initial-boundary value problem is most likely to be well-posed even if nonzero initial and source data are included.Comment: 27 pages, 4 figures; more numerical results and references added, several minor amendments; version accepted for publication in Class. Quantum Gra

Nonradial oscillations of quark stars

Recently, it has been reported that a candidate for a quark star may have been observed. In this article, we pay attention to quark stars with radiation radii in the reported range. We calculate nonradial oscillations of $f$-, $w$- and $w_{\rm II}$-modes. Then, we find that the dependence of the $f$-mode quasi-normal frequency on the bag constant and stellar radiation radius is very strong and different from that of the lowest $w_{\rm II}$-mode quasi-normal frequency. Furthermore we deduce a new empirical formula between the $f$-mode frequency of gravitational waves and the parameter of the equation of state for quark stars. The observation of gravitational waves both of the $f$-mode and of the lowest $w_{\rm II}$-mode would provide a powerful probe for the equation of state of quark matter and the properties of quark stars.Comment: 13 pages, 6 figures, accepted for publication in Phys.Rev.

The rotational modes of relativistic stars: Numerical results

We study the inertial modes of slowly rotating, fully relativistic compact stars. The equations that govern perturbations of both barotropic and non-barotropic models are discussed, but we present numerical results only for the barotropic case. For barotropic stars all inertial modes are a hybrid mixture of axial and polar perturbations. We use a spectral method to solve for such modes of various polytropic models. Our main attention is on modes that can be driven unstable by the emission of gravitational waves. Hence, we calculate the gravitational-wave growth timescale for these unstable modes and compare the results to previous estimates obtained in Newtonian gravity (i.e. using post-Newtonian radiation formulas). We find that the inertial modes are slightly stabilized by relativistic effects, but that previous conclusions concerning eg. the unstable r-modes remain essentially unaltered when the problem is studied in full general relativity.Comment: RevTeX, 29 pages, 31 eps figure

All static spherically symmetric perfect fluid solutions of Einstein's Equations

An algorithm based on the choice of a single monotone function (subject to boundary conditions) is presented which generates all regular static spherically symmetric perfect fluid solutions of Einstein's equations. For physically relevant solutions the generating functions must be restricted by non-trivial integral-differential inequalities. Nonetheless, the algorithm is demonstrated here by the construction of an infinite number of previously unknown physically interesting exact solutions.Comment: Final form to appear in Phys Rev D. Includes a number of clarification

Quasi-Normal Mode Expansion for Linearized Waves in Gravitational Systems

The quasinormal modes (QNM's) of gravitational systems modeled by the Klein-Gordon equation with effective potentials are studied in analogy to the QNM's of optical cavities. Conditions are given for the QNM's to form a complete set, i.e., for the Green's function to be expressible as a sum over QNM's, answering a conjecture by Price and Husain [Phys. Rev. Lett. {\bf 68}, 1973 (1992)]. In the cases where the QNM sum is divergent, procedures for regularization are given. The crucial condition for completeness is the existence of spatial discontinuities in the system, e.g., the discontinuity at the stellar surface in the model of Price and Husain.Comment: 12 pages, WUGRAV-94-

Nonlinear mode coupling in rotating stars and the r-mode instability in neutron stars

We develop the formalism required to study the nonlinear interaction of modes in rotating Newtonian stars in the weakly nonlinear regime. The formalism simplifies and extends previous treatments. At linear order, we elucidate and extend slightly a formalism due to Schutz, show how to decompose a general motion of a rotating star into a sum over modes, and obtain uncoupled equations of motion for the mode amplitudes under the influence of an external force. Nonlinear effects are added perturbatively via three-mode couplings. We describe a new, efficient way to compute the coupling coefficients, to zeroth order in the stellar rotation rate, using spin-weighted spherical harmonics. We apply this formalism to derive some properties of the coupling coefficients relevant to the nonlinear interactions of unstable r-modes in neutron stars, postponing numerical integrations of the coupled equations of motion to a later paper. From an astrophysical viewpoint, the most interesting result of this paper is that many couplings of r-modes to other rotational modes (modes with zero frequencies in the non-rotating limit) are small: either they vanish altogether because of various selection rules, or they vanish to lowest order in the angular velocity. In zero-buoyancy stars, the coupling of three r-modes is forbidden entirely and the coupling of two r-modes to one hybrid rotational mode vanishes to zeroth order in rotation frequency. In incompressible stars, the coupling of any three rotational modes vanishes to zeroth order in rotation frequency.Comment: 62 pages, no figures. Corrected error in computation of coupling coefficients, added new selection rule and an appendix on energy and angular momentum of mode

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

Unequal Mass Binary Black Hole Plunges and Gravitational Recoil

We present results from fully nonlinear simulations of unequal mass binary black holes plunging from close separations well inside the innermost stable circular orbit with mass ratios q = M_1/M_2 = {1,0.85,0.78,0.55,0.32}, or equivalently, with reduced mass parameters $\eta=M_1M_2/(M_1+M_2)^2 = {0.25, 0.248, 0.246, 0.229, 0.183}$. For each case, the initial binary orbital parameters are chosen from the Cook-Baumgarte equal-mass ISCO configuration. We show waveforms of the dominant l=2,3 modes and compute estimates of energy and angular momentum radiated. For the plunges from the close separations considered, we measure kick velocities from gravitational radiation recoil in the range 25-82 km/s. Due to the initial close separations our kick velocity estimates should be understood as a lower bound. The close configurations considered are also likely to contain significant eccentricities influencing the recoil velocity.Comment: 12 pages, 5 figures, to appear in "New Frontiers" special issue of CQ