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

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

r-modes in Relativistic Superfluid Stars

We discuss the modal properties of the $r$-modes of relativistic superfluid neutron stars, taking account of the entrainment effects between superfluids. In this paper, the neutron stars are assumed to be filled with neutron and proton superfluids and the strength of the entrainment effects between the superfluids are represented by a single parameter $\eta$. We find that the basic properties of the $r$-modes in a relativistic superfluid star are very similar to those found for a Newtonian superfluid star. The $r$-modes of a relativistic superfluid star are split into two families, ordinary fluid-like $r$-modes ($r^o$-mode) and superfluid-like $r$-modes ($r^s$-mode). The two superfluids counter-move for the $r^s$-modes, while they co-move for the $r^o$-modes. For the $r^o$-modes, the quantity $\kappa\equiv\sigma/\Omega+m$ is almost independent of the entrainment parameter $\eta$, where $m$ and $\sigma$ are the azimuthal wave number and the oscillation frequency observed by an inertial observer at spatial infinity, respectively. For the $r^s$-modes, on the other hand, $\kappa$ almost linearly increases with increasing $\eta$. It is also found that the radiation driven instability due to the $r^s$-modes is much weaker than that of the $r^o$-modes because the matter current associated with the axial parity perturbations almost completely vanishes.Comment: 14 pages, 4 figures. To appear in Physical Review

R-mode oscillations and rocket effect in rotating superfluid neutron stars. I. Formalism

We derive the hydrodynamical equations of r-mode oscillations in neutron stars in presence of a novel damping mechanism related to particle number changing processes. The change in the number densities of the various species leads to new dissipative terms in the equations which are responsible of the {\it rocket effect}. We employ a two-fluid model, with one fluid consisting of the charged components, while the second fluid consists of superfluid neutrons. We consider two different kind of r-mode oscillations, one associated with comoving displacements, and the second one associated with countermoving, out of phase, displacements.Comment: 10 page

Dynamic modelling of nitrous oxide emissions from three Swedish sludge liquor treatment systems

The objective of this paper is to model the dynamics and validate the results of nitrous oxide (N2O) emissions from three Swedish nitrifying/denitrifying, nitritation and anammox systems treating real anaerobic digester sludge liquor. The Activated Sludge Model No. 1 is extended to describe N2O production by both heterotrophic and autotrophic denitrification. In addition, mass transfer equations are implemented to characterize the dynamics of N2O in the water and the gas phases. The biochemical model is simulated and validated for two hydraulic patterns: (1) a sequencing batch reactor; and (2) a moving-bed biofilm reactor. Results show that the calibrated model is partly capable of reproducing the behaviour of N2O as well as the nitritation/nitrification/denitrification dynamics. However, the results emphasize that additional work is required before N2O emissions from sludge liquor treatment plants can be generally predicted with high certainty by simulations. Continued efforts should focus on determining the switching conditions for different N2O formation pathways and, if full-scale data are used, more detailed modelling of the measurement devices might improve the conclusions that can be drawn.</jats:p

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

(In)finiteness of Spherically Symmetric Static Perfect Fluids

This work is concerned with the finiteness problem for static, spherically symmetric perfect fluids in both Newtonian Gravity and General Relativity. We derive criteria on the barotropic equation of state guaranteeing that the corresponding perfect fluid solutions possess finite/infinite extent. In the Newtonian case, for the large class of monotonic equations of state, and in General Relativity we improve earlier results

Towards absorbing outer boundaries in General Relativity

We construct exact solutions to the Bianchi equations on a flat spacetime background. When the constraints are satisfied, these solutions represent in- and outgoing linearized gravitational radiation. We then consider the Bianchi equations on a subset of flat spacetime of the form [0,T] x B_R, where B_R is a ball of radius R, and analyze different kinds of boundary conditions on \partial B_R. Our main results are: i) We give an explicit analytic example showing that boundary conditions obtained from freezing the incoming characteristic fields to their initial values are not compatible with the constraints. ii) With the help of the exact solutions constructed, we determine the amount of artificial reflection of gravitational radiation from constraint-preserving boundary conditions which freeze the Weyl scalar Psi_0 to its initial value. For monochromatic radiation with wave number k and arbitrary angular momentum number l >= 2, the amount of reflection decays as 1/(kR)^4 for large kR. iii) For each L >= 2, we construct new local constraint-preserving boundary conditions which perfectly absorb linearized radiation with l <= L. (iv) We generalize our analysis to a weakly curved background of mass M, and compute first order corrections in M/R to the reflection coefficients for quadrupolar odd-parity radiation. For our new boundary condition with L=2, the reflection coefficient is smaller than the one for the freezing Psi_0 boundary condition by a factor of M/R for kR > 1.04. Implications of these results for numerical simulations of binary black holes on finite domains are discussed.Comment: minor revisions, 30 pages, 6 figure

Physically Realistic Solutions to the Ernst Equation on Hyperelliptic Riemann Surfaces

We show that the class of hyperelliptic solutions to the Ernst equation (the stationary axisymmetric Einstein equations in vacuum) previously discovered by Korotkin and Neugebauer and Meinel can be derived via Riemann-Hilbert techniques. The present paper extends the discussion of the physical properties of these solutions that was begun in a Physical Review Letter, and supplies complete proofs. We identify a physically interesting subclass where the Ernst potential is everywhere regular except at a closed surface which might be identified with the surface of a body of revolution. The corresponding spacetimes are asymptotically flat and equatorially symmetric. This suggests that they could describe the exterior of an isolated body, for instance a relativistic star or a galaxy. Within this class, one has the freedom to specify a real function and a set of complex parameters which can possibly be used to solve certain boundary value problems for the Ernst equation. The solutions can have ergoregions, a Minkowskian limit and an ultrarelativistic limit where the metric approaches the extreme Kerr solution. We give explicit formulae for the potential on the axis and in the equatorial plane where the expressions simplify. Special attention is paid to the simplest non-static solutions (which are of genus two) to which the rigidly rotating dust disk belongs.Comment: 32 pages, 2 figures, uses pstricks.sty, updated version (October 7, 1998), to appear in Phys. Rev.

Exact relativistic treatment of stationary counter-rotating dust disks I: Boundary value problems and solutions

This is the first in a series of papers on the construction of explicit solutions to the stationary axisymmetric Einstein equations which describe counter-rotating disks of dust. These disks can serve as models for certain galaxies and accretion disks in astrophysics. We review the Newtonian theory for disks using Riemann-Hilbert methods which can be extended to some extent to the relativistic case where they lead to modular functions on Riemann surfaces. In the case of compact surfaces these are Korotkin's finite gap solutions which we will discuss in this paper. On the axis we establish for general genus relations between the metric functions and hence the multipoles which are enforced by the underlying hyperelliptic Riemann surface. Generalizing these results to the whole spacetime we are able in principle to study the classes of boundary value problems which can be solved on a given Riemann surface. We investigate the cases of genus 1 and 2 of the Riemann surface in detail and construct the explicit solution for a family of disks with constant angular velocity and constant relative energy density which was announced in a previous Physical Review Letter.Comment: 32 pages, 1 figure, to appear in Phys. Rev.