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

Toward stable 3D numerical evolutions of black-hole spacetimes

Three dimensional (3D) numerical evolutions of static black holes with excision are presented. These evolutions extend to about 8000M, where M is the mass of the black hole. This degree of stability is achieved by using growth-rate estimates to guide the fine tuning of the parameters in a multi-parameter family of symmetric hyperbolic representations of the Einstein evolution equations. These evolutions were performed using a fixed gauge in order to separate the intrinsic stability of the evolution equations from the effects of stability-enhancing gauge choices.Comment: 4 pages, 5 figures. To appear in Phys. Rev. D. Minor additions to text for clarification. Added short paragraph about inner boundary dependenc

Trapped gravitational wave modes in stars with R>3M

The possibility of trapped modes of gravitational waves appearing in stars with R>3M is considered. It is shown that the restriction to R<3M in previous studies of trapped modes, using uniform density models, is not essential. Scattering potentials are computed for another family of analytic stellar models showing the appearance of a deep potential well for one model with R>3M. However, the provided example is unstable, although it has a more realistic equation of state in the sense that the sound velocity is finite. On the other hand it is also shown that for some stable models belonging to the same family but having R<3M, the well is significantly deeper than that of the uniform density stars. Whether there are physically realistic equations of state which allow stable configurations with trapped modes therefore remains an open problem.Comment: 10 pages, 3 figures, LaTeX2

Gravitational waves from a test particle scattered by a neutron star: Axial mode case

Using a metric perturbation method, we study gravitational waves from a test particle scattered by a spherically symmetric relativistic star. We calculate the energy spectrum and the waveform of gravitational waves for axial modes. Since metric perturbations in axial modes do not couple to the matter fluid of the star, emitted waves for a normal neutron star show only one peak in the spectrum, which corresponds to the orbital frequency at the turning point, where the gravitational field is strongest. However, for an ultracompact star (the radius $R \lesssim 3M$), another type of resonant periodic peak appears in the spectrum. This is just because of an excitation by a scattered particle of axial quasinormal modes, which were found by Chandrasekhar and Ferrari. This excitation comes from the existence of the potential minimum inside of a star. We also find for an ultracompact star many small periodic peaks at the frequency region beyond the maximum of the potential, which would be due to a resonance of two waves reflected by two potential barriers (Regge-Wheeler type and one at the center of the star). Such resonant peaks appear neither for a normal neutron star nor for a Schwarzschild black hole. Consequently, even if we analyze the energy spectrum of gravitational waves only for axial modes, it would be possible to distinguish between an ultracompact star and a normal neutron star (or a Schwarzschild black hole).Comment: 21 pages, revtex, 11 figures are attached with eps files Accepted to Phys. Rev.

Nuclear symmetry energy and the r-mode instability of neutron stars

We analyze the role of the symmetry energy slope parameter $L$ on the {\it r}-mode instability of neutron stars. Our study is performed using both microscopic and phenomenological approaches of the nuclear equation of state. The microscopic ones include the Brueckner--Hartree--Fock approximation, the well known variational equation of state of Akmal, Pandharipande and Ravenhall, and a parametrization of recent Auxiliary Field Diffusion Monte Carlo calculations. For the phenomenological approaches, we use several Skyrme forces and relativisic mean field models. Our results show that the {\it r}-mode instability region is smaller for those models which give larger values of $L$. The reason is that both bulk ($\xi$) and shear ($\eta$) viscosities increase with $L$ and, therefore, the damping of the mode is more efficient for the models with larger $L$. We show also that the dependence of both viscosities on $L$ can be described at each density by simple power-laws of the type $\xi=A_{\xi}L^{B_\xi}$ and $\eta=A_{\eta}L^{B_\eta}$. Using the measured spin frequency and the estimated core temperature of the pulsar in the low-mass X-ray binary 4U 1608-52, we conclude that observational data seem to favor values of $L$ larger than $\sim 50$ MeV if this object is assumed to be outside the instability region, its radius is in the range $11.5-12$($11.5-13$) km, and its mass $1.4M_\odot$($2M_\odot$). Outside this range it is not possible to draw any conclusion on $L$ from this pulsar.Comment: 10 pages, 6 figures. Version published in Physical Review

General Relativistic Rossby-Haurwitz waves of a slowly and differentially rotating fluid shell

We show that, at first order in the angular velocity, the general relativistic description of Rossby-Haurwitz waves (the analogues of r-waves on a thin shell) can be obtained from the corresponding Newtonian one after a coordinate transformation. As an application, we show that the results recently obtained by Rezzolla and Yoshida (2001) in the analysis of Newtonian Rossby-Haurwitz waves of a slowly and differentially rotating, fluid shell apply also in General Relativity, at first order in the angular velocity.Comment: 4 pages. Comment to Class. Quantum Grav. 18(2001)L8

Equilibrium sequences of irrotational binary polytropic stars : The case of double polytropic stars

Solutions to equilibrium sequences of irrotational binary polytropic stars in Newtonian gravity are expanded in a power of $\epsilon=a_0/R$, where R and $a_0$ are the orbital separation of the binary system and the radius of each star for $R=\infty$. For each order of $\epsilon$, we should solve ordinary differential equations for arbitrary polytropic indices n. We show solutions for polytropic indices n= 0.5, 1, 1.5 and 2 up to $\epsilon^6$ orders. Our semi-analytic solutions can be used to check the validity of numerical solutions.Comment: 59 pages including 15 tables and 13 figures, revtex, accepted to Phys. Rev.

Topological and geometrical restrictions, free-boundary problems and self-gravitating fluids

Let (P1) be certain elliptic free-boundary problem on a Riemannian manifold (M,g). In this paper we study the restrictions on the topology and geometry of the fibres (the level sets) of the solutions f to (P1). We give a technique based on certain remarkable property of the fibres (the analytic representation property) for going from the initial PDE to a global analytical characterization of the fibres (the equilibrium partition condition). We study this analytical characterization and obtain several topological and geometrical properties that the fibres of the solutions must possess, depending on the topology of M and the metric tensor g. We apply these results to the classical problem in physics of classifying the equilibrium shapes of both Newtonian and relativistic static self-gravitating fluids. We also suggest a relationship with the isometries of a Riemannian manifold.Comment: 36 pages. In this new version the analytic representation hypothesis is proved. Please address all correspondence to D. Peralta-Sala

A relativistic formalism for computation of irrotational binary stars in quasi equilibrium states

We present relativistic hydrostatic equations for obtaining irrotational binary neutron stars in quasi equilibrium states in 3+1 formalism. Equations derived here are different from those previously given by Bonazzola, Gourgoulhon, and Marck, and have a simpler and more tractable form for computation in numerical relativity. We also present hydrostatic equations for computation of equilibrium irrotational binary stars in first post-Newtonian order.Comment: 5 pages, corrected eqs.(2.10), (2.11) and (3.1

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