Effect of hyperon bulk viscosity on neutron-star r-modes

Neutron stars are expected to contain a significant number of hyperons in addition to protons and neutrons in the highest density portions of their cores. Following the work of Jones, we calculate the coefficient of bulk viscosity due to nonleptonic weak interactions involving hyperons in neutron-star cores, including new relativistic and superfluid effects. We evaluate the influence of this new bulk viscosity on the gravitational radiation driven instability in the r-modes. We find that the instability is completely suppressed in stars with cores cooler than a few times 10^9 K, but that stars rotating more rapidly than 10-30% of maximum are unstable for temperatures around 10^10 K. Since neutron-star cores are expected to cool to a few times 10^9 K within seconds (much shorter than the r-mode instability growth time) due to direct Urca processes, we conclude that the gravitational radiation instability will be suppressed in young neutron stars before it can significantly change the angular momentum of the star.Comment: final PRD version, minor typos etc correcte

Gauge drivers for the generalized harmonic Einstein equations

The generalized harmonic representation of Einstein's equations is manifestly hyperbolic for a large class of gauge conditions. Unfortunately most of the useful gauges developed over the past several decades by the numerical relativity community are incompatible with the hyperbolicity of the equations in this form. This paper presents a new method of imposing gauge conditions that preserves hyperbolicity for a much wider class of conditions, including as special cases many of the standard ones used in numerical relativity: e.g., K freezing, Gamma freezing, Bona-Massó slicing, conformal Gamma drivers, etc. Analytical and numerical results are presented which test the stability and the effectiveness of this new gauge-driver evolution system

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

Shear viscosity of neutron matter from realistic nucleon-nucleon interactions

The calculation of transport properties of Fermi liquids, based on the formalism developed by Abrikosov and Khalatnikov, requires the knowledge of the probability of collisions between quasiparticles in the vicinity of the Fermi surface. We have carried out a numerical study of the shear viscosity of pure neutron matter, whose value plays a pivotal role in determining the stability of rotating neutron stars, in which these processes are described using a state-of-the-art nucleon-nucleon potential model. Within our approach medium modifications of the scattering cross section are consistently taken into account, through an effective interaction obtained from the matrix elements of the bare interaction between correlated states. Inclusion of medium effects lead to a large increase of the viscosity at densities larger than $\sim 0.1$ fm^{-3}.Comment: 4 pages, 4 figures. Corrected typo

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

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

Leptonic contribution to the bulk viscosity of nuclear matter

For beta-equilibrated nuclear matter we estimate the contribution to the bulk viscosity from purely leptonic processes, namely the conversion of electrons to and from muons. For oscillation frequencies in the kiloHertz range, we find that this process provides the dominant contribution to the bulk viscosity when the temperature is well below the critical temperature for superconductivity or superfluidity of the nuclear matter.Comment: 15 pages, LaTeX, new appendix and general clarifications in response to referee comment

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

Data analysis of continuous gravitational wave: Fourier transform-II

In this paper we obtain the Fourier Transform of a continuous gravitational wave. We have analysed the data set for (i) one year observation time and (ii) arbitrary observation time, for arbitrary location of detector and source taking into account the effects arising due to rotational as well as orbital motion of the earth. As an application of the transform we considered spin down and N-component signal analysis.Comment: Accepted in MNRAS, 14 pages, 4 figure