We use a post-Newtonian diagnostic tool to examine numerically generated quasiequilibrium initial data sets for non-spinning double neutron star and neutron star-black hole binary systems. The PN equations include the effects of tidal interactions, parametrized by the compactness of the neutron stars and by suitable values of ``apsidal'' constants, which measure the degree of distortion of stars subjected to tidal forces. We find that the post-Newtonian diagnostic agrees well with the double neutron star initial data, typically to better than half a percent except where tidal distortions are becoming extreme. We show that the differences could be interpreted as representing small residual eccentricity in the initial orbits. In comparing the diagnostic with preliminary numerical data on neutron star-black hole binaries, we find less agreement.

The bending angle of light is a central quantity in the theory of gravitational lensing. We develop an analytical perturbation framework for calculating the bending angle of light rays lensed by a Schwarzschild black hole. Using a perturbation parameter given in terms of the gravitational radius of the black hole and the light ray's impact parameter, we determine an invariant series for the strong-deflection bending angle that extends beyond the standard logarithmic deflection term used in the literature. In the process, we discovered an improvement to the standard logarithmic deflection term. Our perturbation framework is also used to derive as a consistency check, the recently found weak deflection bending angle series. We also reformulate the latter series in terms of a more natural invariant perturbation parameter, one that smoothly transitions between the weak and strong deflection series. We then compare our invariant strong deflection bending-angle series with the numerically integrated exact formal bending angle expression, and find less than 1% discrepancy for light rays as far out as twice the critical impact parameter. The paper concludes by showing that the strong and weak deflection bending angle series together provide an approximation that is within 1% of the exact bending angle value for light rays traversing anywhere between the photon sphere and infinity.

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 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 propose a novel approach for studying $\nu_\mu \to \nu_\tau$ oscillations with extragalactic neutrinos. Active Galactic Nuclei and Gamma Ray Bursts are believed to be sources of ultrahigh energy muon neutrinos. With distances of 100 Mpc or more, they provide an unusually long baseline for possible detection of $\nu_\mu \to \nu_\tau$ with mixing parameters $\Delta m^2$ down to $10^{-17}$eV$^2$, many orders of magnitude below the current accelerator experiments. By solving the coupled transport equations, we show that high-energy $\nu_\tau$'s, as they propagate through the earth, cascade down in energy, producing the enhancement of the incoming $\nu_\tau$ flux in the low energy region, in contrast to the high-energy $\nu_\mu$'s, which get absorbed. For an AGN quasar model we find the $\nu_\tau$ flux to be a factor of 2 to 2.5 larger than the incoming flux in the energy range between $10^2$ GeV and $10^4$ GeV, while for a GRB fireball model, the enhancement is 10%-27% in the same energy range and for zero nadir angle. This enhancement decreases with larger nadir angle, thus providing a novel way to search for $\nu_\tau$ appearance by measuring the angular dependence of the muons. To illustrate how the cascade effect and the $\nu_\tau$ final flux depend on the steepness of the incoming $\nu_\tau$, we show the energy and angular distributions for several generic cases of the incoming tau neutrino flux, $F_\nu^0 \sim E^{-n}$ for n=1,2 and 3.6. We show that for the incoming flux that is not too steep, the signal for the appearance of high-energy $\nu_\tau$ is the enhanced production of lower energy $\mu$ and their distinctive angular dependence, due to the contribution from the $\tau$ decay into $\mu$ just below the detector.