The realization that GRBs are narrowly beamed implied that the actual rate of GRBs is much larger than the observed one. There are 500 unobserved GRBs for each observed one. The lack of a clear trigger makes it hard to detect these unobserved GRBs as orphan afterglows. At late time, hundreds or thousands of years after a GRB, we expect to observe a GRB remnant (GRBR). These remnants could be distinguished from the more frequent SNRs using their different morphology. While SNRs are spherical, GRBRs that arise from a highly collimated flow, are expected to be initially nonspehrical. Using SPH simulations we follow the evolution of a GRBR and calculate the image of the remnant produced by bremsstrahlung and by synchrotron emission. We find that the GRBR becomes spherical after $ \sim 3000 {\rm yr}(E_{51}/n)^{1/3}$ at $R\sim 12{\rm pc} (E_{51}/n)^{1/3}$, where $E_{51}$ is the initial energy in units of $10^{51}{\rm erg}$ and $n$ is the surrounding ISM number density in ${\rm cm}^{-3}$. We expect $0.5 (E_{51}/n)^{1/3}$ non-spherical GRBs per galaxy. Namely, we expect $\sim ~ 20$ non spherical GRBRs with angular sizes $\sim 1 \mu$arcsec within a distance of 10Mpc. These results are modified if there is an underlying spherical supernova. In this case the GRBR will remain spherical only for $\sim150 {\rm yr} (E_{51}/n)^{1/3}$ and the number of non-spherical GRBRs is smaller by a factor of 10 and their size is smaller by a factor of 3.

We introduce an adaptation of the well known Tree+SPH numerical scheme to Post Newtonian (PN) hydrodynamics and gravity. Our code solves the (0+1+2.5)PN equations. These equations include Newtonian hydrodynamics and gravity (0PN), the first order relativistic corrections to those (1PN) and the lowest order gravitational radiation terms (2.5PN). We test various aspects of our code using analytically solvable test problems. We then proceed to study the 1PN effects on binary neutron star coalescence by comparing calculations with and without the 1PN terms. We find that the effect of the 1PN terms is rather small. The largest effect arises with a stiff equation of state for which the maximum rest mass density increases. This could induce black hole formation. The gravitational wave luminosity is also affected.

We model the intermediate time evolution of a "jetted" gamma-ray burst by two blobs of matter colliding with the interstellar medium. We follow the hydrodynamical evolution of this system numerically and calculate the bremsstrahlung and synchrotron images of the remnant. We find that for a burst energy of $10^{51}$ erg the remnant becomes spherical after $\sim 5000$ years when it collects $\sim 50M_\odot$ of interstellar mass. This result is independent of the exact details of the GRB, such as the opening angle. After this time a gamma-ray burst remnant has an expanding sphere morphology. The similarity to a supernova remnant makes it difficult distinguish between the two at this stage. The expected number of non-spherical gamma-ray burst remnants is $\sim0.05$ per galaxy for a beaming factor of 0.01 and a burst energy of $10^{51}$ erg. Our results suggest that that the double-shell object DEM L 316 is not a GRB remnant.

We study the long term evolution of a solar type star that is being disrupted by a super massive (10^6 solar mass) black hole. The evolution is followed from the disruption event, which turns the star into a long thin stream of gas, to the point where some of this gas returns to pericenter and begins its second orbit around the black hole. Following the evolution for this long allows us to determine the amount of mass that is accreted by the black hole. We find that approximately 75% of the returning mass is not accreted but instead becomes unbound, following the large compression characterizing the return to pericenter. The impact of a tidal disruption on the surrounding gas may therefore be like that of two consecutive supernove-type events.