We study the dynamical evolution of a phase-transition-induced collapse neutron star to a hybrid star, which consists of a mixture of hadronic matter and strange quark matter. The collapse is triggered by a sudden change of equation of state, which result in a large amplitude stellar oscillation. The evolution of the system is simulated by using a 3D Newtonian hydrodynamic code with a high resolution shock capture scheme. We find that both the temperature and the density at the neutrinosphere are oscillating with acoustic frequency. However, they are nearly 180$^{\circ}$ out of phase. Consequently, extremely intense, pulsating neutrino/antineutrino fluxes will be emitted periodically. Since the energy and density of neutrinos at the peaks of the pulsating fluxes are much higher than the non-oscillating case, the electron/positron pair creation rate can be enhanced dramatically. Some mass layers on the stellar surface can be ejected by absorbing energy of neutrinos and pairs. These mass ejecta can be further accelerated to relativistic speeds by absorbing electron/positron pairs, created by the neutrino and antineutrino annihilation outside the stellar surface. The possible connection between this process and the cosmological Gamma-ray Bursts is discussed.

The standard model for Type II supernovae explosion, confirmed by the detection of the neutrinos emitted during the supernova explosion, predicts the formation of a compact object, usually assumed to be a neutron star. However, the lack of the detection of a neutron star or pulsar formed in the SN 1987A still remains an unsolved mystery. In this paper we suggest that the newly formed neutron star at the center of SN1987A may undergo a phase transition after the neutrino trapping time scale (~10 s). Consequently the compact remnant of SN 1987A may be a strange quark star, which has a softer equation of state than that of neutron star matter. Such a phase transition can induce the stellar collapse and result in a large amplitude stellar oscillations. We use a three dimensional Newtonian hydrodynamic code to study the time evolution of the temperature and density at the neutrinosphere. Extremely intense pulsating neutrino fluxes, with submillisecond period and with neutrino energy (> 30 MeV) can be emitted because the oscillations of the temperature and density are out of phase almost 180 degree. If this is true we predict that the current X-ray emission from the compact remnant of SN 1987A will be lower than 10^34 erg s-1, and it should be a thermal bremsstrahlung spectrum for a bare strange star with surface temperature of around ~10^7 K.

Previously we found that large amplitude $r$-modes could decay catastrophically due to nonlinear hydrodynamic effects. In this paper we found the particular coupling mechanism responsible for this catastrophic decay, and identified the fluid modes involved. We find that for a neutron star described by a polytropic equation of state with polytropic index $\Gamma=2$, the coupling strength of the particular three-mode interaction causing the decay is strong enough that the usual picture of the $r$-mode instability with a flow pattern dominated by that of an $r$-mode can only be valid for the dimensionless $r$-mode amplitude less than $O(10^{-2})$.

(Abridged) We study the gravitational radiation from gravitational collapses of rapidly rotating neutron stars induced by a phase-transition from normal nuclear matter to a mixed phase of quark and nuclear matter in the core of the stars. The study is based on self-consistent three dimensional hydrodynamic simulations with Newtonian gravity and a high resolution shock capturing scheme, and the quadrupole formula of gravitational radiation. The quark matter of the mixed phase is described by the MIT bag model and the normal nuclear matter is described by an ideal fluid equation of state (EOS). 1. We determined the magnitudes of the emitted gravitational waves for several collapse scenarios. 2. We determined the types and frequencies of the fluid oscillation modes excited by the process. In particular, we find that the gravitational wave signals produced by the collapses are dominated by the fundamental quadrupole and quasi-radial modes of the final equilibrium configurations. In some collapse scenarios, we find that the oscillations are damped out within a few dynamical timescales due to the growth of differential rotations and the formation of strong shock waves. 3. We showed that the spectrum of the gravitational wave signals is sensitive to the EOS, implying that the detection of such gravitational waves could provide useful constraints on the EOS of newly born quark stars. 4. For the range of rotation periods we have studied, we found no sign of the development of nonaxisymmetric dynamical instabilities in the collapse process.

We have carried out numerical simulations of strongly gravitating systems based on the Einstein equations coupled to the relativistic hydrodynamic equations using adaptive mesh refinement (AMR) techniques. We show AMR simulations of NS binary inspiral and coalescence carried out on a workstation having an accuracy equivalent to that of a $1025^3$ regular unigrid simulation, which is, to the best of our knowledge, larger than all previous simulations of similar NS systems on supercomputers. We believe the capability opens new possibilities in general relativistic simulations.

This paper reports on our effort in modeling realistic astrophysical neutron star binaries in general relativity. We analyze under what conditions the conformally flat quasiequilibrium (CFQE) approach can generate ``astrophysically relevant'' initial data, by developing an analysis that determines the violation of the CFQE approximation in the evolution of the binary described by the full Einstein theory. We show that the CFQE assumptions significantly violate the Einstein field equations for corotating neutron stars at orbital separations nearly double that of the innermost stable circular orbit (ISCO) separation, thus calling into question the astrophysical relevance of the ISCO determined in the CFQE approach. With the need to start numerical simulations at large orbital separation in mind, we push for stable and long term integrations of the full Einstein equations for the binary neutron star system. We demonstrate the stability of our numerical treatment and analyze the stringent requirements on resolution and size of the computational domain for an accurate simulation of the system.

We study the gravitational wave signals emitted from phase-transition induced collapses of rapidly rotating neutron stars to strange stars by performing 3D numerical simulations. Our preliminary results suggest that the complete conversion of neutron stars to strange stars would occur within a fraction of millisecond. We also find that the gravitational waves generated from the collapse process may be detectable by the advanced LIGO for reasonable source distance. In addition, the study such gravitational wave signals would put useful constraint on the parameters of QCD.

An approach to general relativity based on conformal flatness and quasiequilibrium (CFQE) assumptions has played an important role in the study of the inspiral dynamics and in providing initial data for fully general relativistic numerical simulations of coalescing compact binaries. However, the regime of validity of the approach has never been established. To this end, we develop an analysis that determines the violation of the CFQE approximation in the evolution of the binary described by the full Einstein theory. With this analysis, we show that the CFQE assumption is significantly violated even at relatively large orbital separations in the case of corotational neutron star binaries. We also demonstrate that the innermost stable circular orbit (ISCO) determined in the CFQE approach for corotating neutron star binaries may have no astrophysical significance.

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.

We perform fully relativistic hydrodynamic simulations of the deceleration and lateral expansion of a relativistic jet as it expands into an ambient medium. The hydrodynamic calculations use a 2D adaptive mesh refinement (AMR) code, which provides adequate resolution of the thin shell of matter behind the shock. We find that the sideways propagation is different than predicted by simple analytic models. The physical conditions at the sides of the jet are found to be significantly different than at the front of the jet, and most of the emission occurs within the initial opening angle of the jet. The light curves, as seen by observers at different viewing angles with respect to the jet axis, are then calculated assuming synchrotron emission. For an observer along the jet axis, we find a sharp achromatic `jet break' in the light curve at frequencies above the typical synchrotron frequency, at $t_{jet}\approx 5.8(E_{52}/n_1)^{1/3}(\theta_0/0.2)^{8/3}$ days, while the temporal decay index $\alpha$ ($F_{\nu}\propto t^{\alpha}$) after the break is steeper than $-p$ ($\alpha=-2.85$ for $p=2.5$). At larger viewing angles $t_{jet}$ increases and the jet break becomes smoother.

We describe fully relativistic three dimensional calculations of the slowing down and spreading of a relativistic jet by an external medium like the ISM. We calculate the synchrotron spectra and light curves using the conditions determined by the hydrodynamic calculations. Preliminary results with a moderate resolution are presented here. Higher resolution calculations are in progress.

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.

We study the thermal evolution of a pulsar after a glitch in which the energy is released from a relative compact region. A set of relativistic thermal transport and energy balance equations is used to study the thermal evolution, without making the assumption of spherical symmetry. We use an exact cooling model to solve this set of differential equtions. Our results differ significantly from those obtained under the assumption of spherical symmetry. Even for young pulsars with a hot core like the Vela pulsar, we find that a detectable hot spot can be observed after a glitch. The results suggest that the intensity variation and the relative phases of hard X-ray emissions in different epoches can provide important information on the equation of state.

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).