We assemble the equations of general relativistic magnetohydrodynamics (MHD) in 3+1 form. These consist of the complete coupled set of Maxwell equations for the electromagnetic field, Einstein's equations for the gravitational field, and the equations of relativistic MHD for a perfectly conducting ideal gas. The adopted form of the equations is suitable for evolving numerically a relativistic MHD fluid in a dynamical spacetime characterized by a strong gravitational field.

Thomas W. Baumgarte Stuart L. Shapiro

Solving dynamical problems in general relativity requires the full machinery of numerical relativity. Wilson has proposed a simpler but approximate scheme for systems near equilibrium, like binary neutron stars. We test the scheme on isolated, rapidly rotating, relativistic stars. Since these objects are in equilibrium, it is crucial that the approximation work well if we are to believe its predictions for more complicated systems like binaries. Our results are very encouraging.

Gregory B. Cook Stuart L. Shapiro Saul A. Teukolsky

We present a method for constructing equilibrium disks with net angular momentum in general relativity. The method solves the relativistic Vlasov equation coupled to Einstein's equations for the gravitational field. We apply the method to construct disks that are relativistic versions of Newtonian Kalnajs disks. In Newtonian gravity these disks are analytic, and are stable against ring formation for certain ranges of their velocity dispersion. We investigate the existence of fully general relativistic equilibrium sequences for differing values of the velocity dispersion. These models are the first rotating, relativistic disk solutions of the collisionless Boltzman equation.

A. Katrin Schenk Stuart L. Shapiro Saul A. Teukolsky

We experiment with modifications of the BSSN form of the Einstein field equations (a reformulation of the ADM equations) and demonstrate how these modifications affect the stability of numerical black hole evolution calculations. We use excision to evolve both non-rotating and rotating Kerr-Schild black holes in octant and equatorial symmetry, and without any symmetry assumptions, and obtain accurate and stable simulations for specific angular momenta J/M of up to about 0.9M.

Hwei-Jang Yo Thomas W. Baumgarte Stuart L. Shapiro

The merger of binary neutron stars is likely to lead to differentially rotating remnants. In this paper we numerically construct models of differentially rotating neutron stars in general relativity and determine their maximum allowed mass. We model the stars adopting a polytropic equation of state and tabulate maximum allowed masses as a function of differential rotation and stiffness of the equation of state. We also provide a crude argument that yields a qualitative estimate of the effect of stiffness and differential rotation on the maximum allowed mass.

Nicholas D. Lyford Thomas W. Baumgarte Stuart L. Shapiro

Primordial black holes (PBHs), if captured by neutron stars (NSs), would emit a characteristic gravitational wave (GW) signal as they orbit inside the host star. We identify a specific and qualitatively new feature of these signals, namely quasi-periodic beats caused by the precession of noncircular PBH orbits. We demonstrate numerically and analytically that the beat frequency depends rather sensitively on the NS structure, so that hypothetical future observations with next-generation GW detectors would provide valuable constraints on the nuclear equation of state.

Thomas W. Baumgarte Stuart L. Shapiro

Long-period radio transients have unusual properties that challenge their interpretation as pulsars or magnetars. We examine whether they might instead be powered by primordial black holes (PBHs) making repeated passages through a host star, thereby providing a signature of elusive dark-matter candidates. We demonstrate that constraints derived from the transients' period and period derivative alone already rule out this scenario for most potential host stars. While white dwarfs may satisfy these constraints, they are unlikely to capture PBHs in the required mass range.

Thomas W. Baumgarte Stuart L. Shapiro

Diffferentially rotating stars can support significantly more mass in equilibrium than nonrotating or uniformly rotating stars, according to general relativity. The remnant of a binary neutron star merger may give rise to such a ``hypermassive'' object. While such a star may be dynamically stable against gravitational collapse and bar formation, the radial stabilization due to differential rotation is likely to be temporary. Magnetic braking and viscosity combine to drive the star to uniform rotation, even if the seed magnetic field and the viscosity are small. This process inevitably leads to delayed collapse, which will be accompanied by a delayed gravitational wave burst and, possibly, a gamma-ray burst. We provide a simple, Newtonian, MHD calculation of the braking of differential rotation by magnetic fields and viscosity. The star is idealized as a differentially rotating, infinite cylinder consisting of a homogeneous, incompressible conducting gas. We solve analytically the simplest case in which the gas has no viscosity and the star resides in an exterior vacuum. We treat numerically cases in which the gas has internal viscosity and the star is embedded in an exterior, low-density, conducting medium. Our evolution calculations are presented to stimulate more realistic MHD simulations in full 3+1 general relativity. They serve to identify some of the key physical and numerical parameters, scaling behavior and competing timescales that characterize this important process.

Stuart L. Shapiro

The collapse of a uniformaly rotating, supermassive star (SMS) to a supermassive black hole (SMBH) has been followed recently by means of hydrodynamic simulations in full general relativity. The initial SMS of arbitrary mass M in these simulations rotates uniformly at the mass--shedding limit and is marginally unstable to radial collapse. The final black hole has mass M_h/M = 0.9 and and spin J_h/M_h^2 = 0.75, approximately. The remaining mass goes into a disk of mass M_disk/M = 0.1, also approximately. Here we show that these black hole and disk parameters can be calculated analytically from the initial stellar density and angular momentum distribution. The analytic calculation thereby corroborates and provides a simple physical explanation for the computational discovery that SMS collapse inevitably terminates in the simultaneous formation of a SMBH and a rather substantial ambient disk. This disk arises even though the total spin of the progenitor star, J/M^2 = 0.97, is safely below the Kerr limit. The calculation performed here applies to any marginally unstable n = 3 polytrope uniformly rotating at the break--up speed, independent of stellar mass or the source of internal pressure. It illustrates how the black hole and disk parameters can be determined for the collapse of other types of stars with different initial density and rotation profiles.

Stuart L. Shapiro Masaru Shibata

Simulations in general relativity show that the outcome of collapse of a marginally unstable, uniformly rotating star spinning at the mass-shedding limit depends critically on the equation of state. For a very stiff equation of state, which is likely to characterize a neutron star, essentially all of the mass and angular momentum of the progenitor are swallowed by the Kerr black hole formed during the collapse, leaving nearly no residual gas to form a disk. For a soft equation of state with an adiabatic index \Gamma - 4/3 << 1, which characterizes a very massive or supermassive star supported predominantly by thermal radiation pressure, as much as 10% of the mass of the progenitor avoids capture and goes into a disk about the central hole. We present a semi-analytic calculation that corroborates these numerical findings and shows how the final outcome of such a collapse may be determined from simple physical considerations. In particular, we employ a simple energy variational principle with an approximate, post-Newtonian energy functional to determine the structure of a uniformly rotating, polytropic star at the onset of collapse as a function of polytropic index n, where \Gamma = 1+1/n. We then use this data to calculate the mass and spin of the final black hole and ambient disk. We show that the fraction of the total mass that remains in the disk falls off sharply as 3-n (equivalently, \Gamma - 4/3) increases.

Stuart L. Shapiro