Abstract
Motivated by a recent paper by the Potsdam numerical relativity group, we have constructed a new numerical code for hydrodynamic simulation of axisym- metric systems in full general relativity. In this code, we solve the Einstein field equation using Cartesian coordinates with appropriate boundary conditions. On the other hand, the hydrodynamic equations are solved in cylindrical coordinates. Using this code, we perform simulations to study axisymmetric collapse of rotat- ing stars, which thereby become black holes or new compact stars, in full general relativity. To investigate the effects of rotation on the criterion for prompt col- lapse to black holes, we first adopt a polytropic equation of state, P = Kρ , where P, ρ, and K are the pressure, rest mass density, and polytropic constant, with Γ = 2. In this case, the collapse is adiabatic (i.e., no change in entropy), and we can focus on the bare effect of rotation. As the initial conditions, we prepare rigidly and differentially rotating stars in equilibrium and then decrease the pres- sure to induce collapse. In this paper, we consider cases in which q ≡ J/M 2 g < 1, where J and Mg are the angular momentum and the gravitational mass. It is found that the criterion of black hole formation is strongly dependent on the an- gular momentum parameter q. For q < 0.5, the criterion is not strongly sensitive to q; more precisely, if the rest mass is slightly larger than the maximum allowed value of spherical stars, a black hole is formed. However, for q <∼ 1, it changes significantly: For q ≃ 0.9, the maximum allowed rest mass becomes ∼ 70-80% larger than that for spherical stars. These findings depend only weakly on the rotational profiles given initially. We then report the results for simulations em- ploying a Γ-law equation of state P = (Γ − 1)ρε, where ε is the specific internal energy, to study effects of shock heating. We find that the effects of shock heating are particularly important for preventing prompt collapse to black holes in the case of large q (i.e., q = O(1)).
All Related Versions