This is the second in a series of papers on the construction and validation of a three-dimensional code for the solution of the coupled system of the Einstein equations and of the general relativistic hydrodynamic equations, and on the application of this code to problems in general relativistic astrophysics. In particular, we report on the accuracy of our code in the long-term dynamical evolution of relativistic stars and on some new physics results obtained in the process of code testing. The tests involve single non-rotating stars in stable equilibrium, non-rotating stars undergoing radial and quadrupolar oscillations, non-rotating stars on the unstable branch of the equilibrium configurations migrating to the stable branch, non-rotating stars undergoing gravitational collapse to a black hole, and rapidly rotating stars in stable equilibrium and undergoing quasi-radial oscillations. The numerical evolutions have been carried out in full general relativity using different types of polytropic equations of state using either the rest-mass density only, or the rest-mass density and the internal energy as independent variables. New variants of the spacetime evolution and new high resolution shock capturing (HRSC) treatments based on Riemann solvers and slope limiters have been implemented and the results compared with those obtained from previous methods. Finally, we have obtained the first eigenfrequencies of rotating stars in full general relativity and rapid rotation. A long standing problem, such frequencies have not been obtained by other methods. Overall, and to the best of our knowledge, the results presented in this paper represent the most accurate long-term three-dimensional evolutions of relativistic stars available to date.

Shapiro put forth a conjecture stating that neutron stars in head-on collisions (infalling from infinity) will not collapse to black holes before neutrino cooling, independent of the mass of the neutron stars. In a previous paper we carried out a numerical simulation showing a counter example based on 1.4 $M_{\odot}$ neutron stars, and provided an analysis explaining why Shapiro's argument was not applicable for this case. A recent paper by Shapiro put forth an argument suggesting that numerical simulations of the 1.4 $M_{\odot}$ collisions could not disprove the conjecture with the accuracy that is presently attainable. We show in this paper that this argument is not applicable for the same reason that the Shapiro conjecture is not applicable to the 1.4 $M_{\odot}$ neutron star collision, namely, the collision is too dynamical to be treated by quasi-equilibrium arguments.

We study the question of prompt vs. delayed collapse in the head-on collision of two neutron stars. We show that the prompt formation of a black hole is possible, contrary to a conjecture of Shapiro which claims that collapse is delayed until after neutrino cooling. We discuss the insight provided by Shapiro's conjecture and its limitation. An understanding of the limitation of the conjecture is provided in terms of the many time scales involved in the problem. General relativistic simulations in the Einstein theory with the full set of Einstein equations coupled to the general relativistic hydrodynamic equations are carried out in our study.

This is the first in a series of papers on the construction and validation of a three-dimensional code for general relativistic hydrodynamics, and its application to general relativistic astrophysics. This paper studies the consistency and convergence of our general relativistic hydrodynamic treatment and its coupling to the spacetime evolutions described by the full set of Einstein equations with a perfect fluid source. The numerical treatment of the general relativistic hydrodynamic equations is based on high resolution shock capturing schemes. These schemes rely on the characteristic information of the system. A spectral decomposition for general relativistic hydrodynamics suitable for a general spacetime metric is presented. Evolutions based on three different approximate Riemann solvers coupled to four different discretizations of the Einstein equations are studied and compared. The coupling between the hydrodynamics and the spacetime (the right and left hand side of the Einstein equations) is carried out in a treatment which is second order accurate in {\it both} space and time. Convergence tests for all twelve combinations with a variety of test beds are studied, showing consistency with the differential equations and correct convergence properties. The test-beds examined include shocktubes, Friedmann-Robertson-Walker cosmology tests, evolutions of self-gravitating compact (TOV) stars, and evolutions of relativistically boosted TOV stars. Special attention is paid to the numerical evolution of strongly gravitating objects, e.g., neutron stars, in the full theory of general relativity, including a simple, yet effective treatment for the surface region of the star (where the rest mass density is abruptly dropping to zero).