Engineering PapersSearch

Engineering topics

Teukolsky, Saul A.

Publications and source records attributed to Teukolsky, Saul A..

Computing supernova collapse to neutron stars and black holes

We present a new numerical code for spherical hydrodynamics in general relativity. The code can handle gravitational collapse to a neutron star or to a black hole without the appearance of singularities. Moreover, the variables and equations in the code are very similar to those appearing in traditional Lagrangian supernova codes. Any such existing code can thus be easily adapted to treat collapse where the final fate is uncertain and may be either a neutron star or a black hole. The code is based on the formulation of Hernandez & Misner, in which retarded time is used as coordinate. This prevents the computational grid from penetrating inside any black hole that may form. We present the equations and a complete finite difference scheme for the adiabatic evolution of a fluid that obeys a gamma-law equation of state. We summarize the results of several testbed calculations performed to check our code. We also give the transformation of the analytic Oppenheimer-Snyder solution for homogeneous dust collapse to our coordinate system.

Baumgarte, Thomas W.

Rapidly rotating neutron stars in general relativity: Realistic equations of state

We construct equilibrium sequences of rotating neutron stars in general relativity. We compare results for 14 nuclear matter equations of state. We determine a number of important physical parameters for such stars, including the maximum mass and maximum spin rate. The stability of the configurations to quasi-radial perturbations is assessed. We employ a numerical scheme particularly well suited to handle rapid rotation and large departures from spherical symmetry. We provide an extensive tabulation of models for future reference. Two classes of evolutionary sequences of fixed baryon rest mass and entropy are explored: normal sequences, which behave very much like Newtonian sequences, and supramassive sequences, which exist for neutron stars solely because of general relativistic effects. Adiabatic dissipation of energy and angular momentum causes a star to evolve in quasi-stationary fashion along an evolutionary sequence. Supramassive sequences have masses exceeding the maximum mass of a nonrotating neutron star. A supramassive star evolves toward eventual catastrophic collapse to a black hole. Prior to collapse, the star actually spins up as it loses angular momentum, an effect that may provide an observable precursor to gravitational collapse to a black hole.

Cook, Gregory B.

Rapidly rotating polytropes in general relativity

We construct an extensive set of equilibrium sequences of rotating polytropes in general relativity. We determine a number of important physical parameters of such stars, including maximum mass and maximum spin rate. The stability of the configurations against quasi-radial perturbations is diagnosed. Two classes of evolutionary sequences of fixed rest mass and entropy are explored: normal sequences which behave very much like Newtonian evolutionary sequences, and supramassive sequences which exist solely because of relativistic effects. Dissipation leading to loss of angular momentum causes a star to evolve in a quasi-stationary fashion along an evolutionary sequence. Supramassive sequences evolve towards eventual catastrophic collapse to a black hole. Prior to collapse, the star must spin up as it loses angular momentum, an effect which may provide an observational precursor to gravitational collapse to a black hole.

Cook, Gregory B.

A hydrodynamical model for the explosion of a neutron star just below the minimum mass

The instability of a neutron star at the minimum mass is investigated using a hydrodynamical model to follow the evolution of the unstable star. A detailed analysis of the linear stability of the equilibrium configurations near the minimum mass is performed by solving the radial eigenvalue problem for the fundamental mode. The minimum mass configuration M(mmc) of 0.196 solar mass is found to be stable to small perturbations. Stability to radial perturbations is lost only at a lower critical mass of 0.16 solar mass, corresponding to about 0.8 M(mmc). It is shown that a simple dynamical model constructed using a 3-polytrope equation of state for hot dense matter reproduces the key dynamical features of the instability in the explosion phase.

Colpi, Monica

Orbital evolution of the PSR1257+12 planetary system

The detection of orbital perturbation effects in the PSR1257+12 timing data would provide irrefutable confirmation that planets are indeed orbiting the pulsar. Here we give an overview of how perturbation effects are expected to affect the orbital elements of the two planets over the next few years. In particular, we give a simple calculation of resonant perturbations, including the nonlinear effects which could be important if sin i less than about 0.1. We also present a new analysis of the effects of close encounters between the two planets and we discuss their detectability.

Rasio, Frederic A.

Spin-up of a rapidly rotating star by angular momentum loss - Effects of general relativity

It has recently been shown that a rapidly rotating Newtonian star can spin up by radiating angular momentum. Extremely fast pulsars losing energy and angular momentum by magnetic dipole radiation or gravitational radiation may exhibit this behavior. Here, we show that this phenomenon is more widespread for rapidly rotating stars in general relativity. We construct and tabulate polytropic sequences of fully relativistic rotating stars of constant rest mass and entropy. We find that the range of adiabatic indices allowing spin-up extends somewhat above 4/3 because of the nonlinear effects of relativistic gravity. In addition, there is a new class of 'supramassive' stars which will inevitably spin up by losing angular momentum regardless of their equation of state. A supramassive star, spinning up via angular momentum loss, will ultimately evolve until it becomes unstable to catastrophic collapse to a black hole. Spin-up in a rapidly rotating star may thus be an observational precursor to such collapse.

Cook, Gregory B.

Equilibrium stellar systems with spindle singularities

Equilibrium sequences of axisymmetric Newtonian clusters that tend toward singular states are constructed. The distribution functions are chosen to be of the form f = f(E, Jz). The numerical method then determines the density and gravitational potential self-consistently to satisfy Poisson's equation. For the prolate models, spindle singularities arise from the depletion of angular momentum near the symmetry axis. While the resulting density enhancement is confined to the region near the axis, the influence of the spindle extends much further out through its tidal gravitational field. Centrally condensed prolate clusters may contain strong-field regions even though the spindle mass is small and the mean cluster eccentricity is not extreme. While the calculations performed here are entirely Newtonian, the issue of singularities is an important topic in general relativity. Equilibrium solutions for relativistic star clusters can provide a testing ground for exploring this issue. The methods used in this paper for building nonspherical clusters can be extended to relativistic systems.

Shapiro, Stuart L.

Testing nuclear theory using the 0.5 ms pulsar

It is argued that a uniformly rotating neutron star with period of 0.5 ms or less may be unable to spin down to become a slowly rotating neutron star if any of the currently proposed equations of state are correct. Asssuming that the neutron star in SN1987A is 'typical' and does not collapse to a black hole as it spins down to a rotational period much larger than 0.5 ms, then its existence may actually invalidate all previously viable nuclear equations of state. An approximate but analytical test that can be used to identify untenable equations of state in the future is proposed.

Shapiro, Stuart L.