General relativistic hydrodynamics with viscosity: Contraction, catastrophic collapse, and disk formation in hypermassive neutron stars

Viscosity and magnetic fields drive differentially rotating stars toward uniform rotation, and this process has important consequences in many astrophysical contexts. For example, merging binary neutron stars can form a "hypermassive" remnant, i.e., a differentially rotating star with a mass greater than would be possible for a uniformly rotating star. The removal of the centrifugal support provided by differential rotation can lead to delayed collapse of the remnant to a black hole, accompanied by a delayed burst of gravitational radiation. Both magnetic fields and viscosity alter the structure of differentially rotating stars on secular time scales, and tracking this evolution presents a strenuous challenge to numerical hydrodynamic codes. Here, we present the first evolutions of rapidly rotating stars with shear viscosity in full general relativity. We self-consistently include viscosity in our relativistic hydrodynamic code by solving the fully relativistic Navier-Stokes equations. We perform these calculations both in axisymmetry and in full 3+1 dimensions. In axisymmetry, the resulting reduction in computational costs allows us to follow secular evolution with high resolution over dozens of rotation periods (thousands of M). We find that viscosity operating in a hypermassive star generically leads to the formation of a compact, uniformly rotating core surrounded by a low-density disk. These uniformly rotating cores are often unstable to gravitational collapse. We follow the collapse in such cases and determine the mass and the spin of the final black hole and ambient disk. However, viscous braking of differential rotation in hypermassive neutron stars does not always lead to catastrophic collapse, especially when viscous heating is substantial. The stabilizing influences of viscous heating, which generates enhanced thermal pressure, and centrifugal support prevent collapse in some cases, at least until the star cools. In all cases studied, the rest mass of the resulting disk is found to be 10-20 % of the original star, whether surrounding a uniformly rotating core or a rotating black hole. This study represents an important step toward understanding secular effects in relativistic stars and foreshadows more detailed, future simulations, including those involving magnetic fields.