Many codes have been developed to study highly relativistic, magnetized flows around and inside compact objects. Depending on the adopted formalism, some of these codes evolve the vector potential $\mathbf{A}$, and others evolve the magnetic field $\mathbf{B}=\nabla\times\mathbf{A}$ directly. Given that these codes possess unique strengths, sometimes it is desirable to start a simulation using a code that evolves $\mathbf{B}$ and complete it using a code that evolves $\mathbf{A}$. Thus transferring the data from one code to another would require an inverse curl algorithm. This paper describes two new inverse curl techniques in the context of Cartesian numerical grids: a cell-by-cell method, which scales approximately linearly with the numerical grid, and a global linear algebra approach, which has worse scaling properties but is generally more robust (e.g., in the context of a magnetic field possessing some nonzero divergence). We demonstrate these algorithms successfully generate smooth vector potential configurations in challenging special and general relativistic contexts.

During the merger of a black hole and a neutron star, baryonic mass can become unbound from the system. Because the ejected material is extremely neutron-rich, the r-process rapidly synthesizes heavy nuclides as the material expands and cools. In this work, we map general relativistic models of black hole-neutron star (BHNS) mergers into a Newtonian smoothed particle hydrodynamics (SPH) code and follow the evolution of the thermodynamics and morphology of the ejecta until the outflows become homologous. We investigate how the subsequent evolution depends on our mapping procedure and find that the results are robust. Using thermodynamic histories from the SPH particles, we then calculate the expected nucleosynthesis in these outflows while varying the level of neutrino irradiation coming from the postmerger accretion disk. We find that the ejected material robustly produces r-process nucleosynthesis even for unrealistically high neutrino luminosities, due to the rapid velocities of the outflow. Nonetheless, we find that neutrinos can have an impact on the detailed pat- tern of the r-process nucleosynthesis. Electron neutrinos are captured by neutrons to produce protons while neutron capture is occurring. The produced protons rapidly form low mass seed nuclei for the r-process. These low mass seeds are eventually incorporated into the first r-process peak at A~78, producing mainly Ge and Se. We consider the mechanism of this process in detail and discuss if it can impact galactic chemical evolution of the first peak r-process nuclei.

To make relevant predictions about observable emission, hydrodynamical simulation codes must employ schemes that account for radiative losses, but the large dimensionality of accurate radiative transfer schemes is often prohibitive. Stamatellos and collaborators introduced a scheme for smoothed particle hydrodynamics (SPH) simulations based on the notion of polytropic pseudo-clouds that uses only local quantities to estimate cooling rates. The computational approach is extremely efficient and works well in cases close to spherical symmetry, such as in star formation problems. Unfortunately, the method, which takes the local gravitational potential as an input, can be inaccurate when applied to non-spherical configurations, limiting its usefulness when studying disks or stellar collisions, among other situations of interest. Here, we introduce the "pressure scale height method," which incorporates the fluid pressure scale height into the determination of column densities and cooling rates, and show that it produces more accurate results across a wide range of physical scenarios while retaining the computational efficiency of the original method. The tested models include spherical polytropes as well as disks with specified density and temperature profiles. We focus on applying our techniques within an SPH code, although our method can be implemented within any particle-based Lagrangian or grid-based Eulerian hydrodynamic scheme. Our new method may be applied in a broad range of situations, including within the realm of stellar interactions, collisions, and mergers.

We review the current status of studies of the coalescence of binary neutron star systems. We begin with a discussion of the formation channels of merging binaries and we discuss the most recent theoretical predictions for merger rates. Next, we turn to the quasi-equilibrium formalisms that are used to study binaries prior to the merger phase and to generate initial data for fully dynamical simulations. The quasi-equilibrium approximation has played a key role in developing our understanding of the physics of binary coalescence and, in particular, of the orbital instability processes that can drive binaries to merger at the end of their lifetimes. We then turn to the numerical techniques used in dynamical simulations, including relativistic formalisms, (magneto-)hydrodynamics, gravitational-wave extraction techniques, and nuclear microphysics treatments. This is followed by a summary of the simulations performed across the field to date, including the most recent results from both fully relativistic and microphysically detailed simulations. Finally, we discuss the likely directions for the field as we transition from the first to the second generation of gravitational-wave interferometers and while supercomputers reach the petascale frontier.

Numerical calculations of merging black hole binaries indicate that asymmetric emission of gravitational radiation can kick the merged black hole at up to thousands of km/s, and a number of systems have been observed recently whose properties are consistent with an active galactic nucleus containing a supermassive black hole moving with substantial velocity with respect to its broader accretion disk. We study here the effect of an impulsive kick delivered to a black hole on the dynamical evolution of its accretion disk using a smoothed particle hydrodynamics code, focusing attention on the role played by the kick angle with respect to the orbital angular momentum vector of the pre-kicked disk. We find that for more vertical kicks, for which the angle between the kick and the normal vector to the disk $\theta\lesssim 30^\circ$, a gap remains present in the inner disk, in accordance with the prediction from an analytic collisionless Keplerian disk model, while for more oblique kicks with $\theta\gtrsim 45^\circ$, matter rapidly accretes toward the black hole. There is a systematic trend for higher potential luminosities for more oblique kick angles for a given black hole mass, disk mass and kick velocity, and we find large amplitude oscillations in time in the case of a kick oriented $60^\circ$ from the vertical.

We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The Toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, compact objects, and collapsing stars, as well as a full suite of analysis tools. The Einstein Toolkit is currently based on the Cactus Framework for high-performance computing and the Carpet adaptive mesh refinement driver. It implements spacetime evolution via the BSSN evolution system and general-relativistic hydrodynamics in a finite-volume discretization. The toolkit is under continuous development and contains many new code components that have been publicly released for the first time and are described in this article. We discuss the motivation behind the release of the toolkit, the philosophy underlying its development, and the goals of the project. A summary of the implemented numerical techniques is included, as are results of numerical test covering a variety of sample astrophysical problems.

The tidal breakup of binary star systems by the supermassive black hole (SMBH) in the center of the galaxy has been suggested as the source of both the observed sample of hypervelocity stars (HVSs) in the halo of the Galaxy and the S-stars that remain in tight orbits around Sgr A*. Here, we use a post-Newtonian N-body code to study the dynamics of main-sequence binaries on highly elliptical bound orbits whose periapses lie close to the SMBH, determining the properties of ejected and bound stars as well as collision products. Unlike previous studies, we follow binaries that remain bound for several revolutions around the SMBH, finding that in the case of relatively large periapses and highly inclined binaries the Kozai resonance can lead to large periodic oscillations in the internal binary eccentricity and inclination. Collisions and mergers of the binary elements are found to increase significantly for multiple orbits around the SMBH, while HVSs are primarily produced during a binary's first passage. This process can lead to stellar coalescence and eventually serve as an important source of young stars at the galactic center.

Black hole-neutron star (BHNS) binaries are expected to be among the leading sources of gravitational waves observable by ground-based detectors, and may be the progenitors of short-hard gamma ray bursts (SGRBs) as well. Here, we discuss our new fully general relativistic calculations of merging BHNS binaries, which use high-accuracy, low-eccentricity, conformal thin-sandwich configurations as initial data. Our evolutions are performed using the moving puncture method and include a fully relativistic, high-resolution shock-capturing hydrodynamics treatment. Focusing on systems in which the neutron star is irrotational and the black hole is nonspinning with a 3:1 mass ratio, we investigate the inspiral, merger, and disk formation in the system. We find that the vast majority of material is promptly accreted and no more than 3% of the neutron star's rest mass is ejected into a tenuous, gravitationally bound disk. We find similar results for mass ratios of 2:1 and 1:1, even when we reduce the NS compaction in the 2:1 mass ratio case. These ambient disks reach temperatures suitable for triggering SGRBs, but their masses may be too small to produce the required total energy output. We measure gravitational waveforms and compute the effective strain in frequency space, finding measurable differences between our waveforms and those produced by binary black hole mergers within the advanced LIGO band. These differences appear at frequencies corresponding to the emission that occurs when the NS is tidally disrupted and accreted by the black hole. The resulting information about the radius of the neutron star may be used to constrain the neutron star equation of state.

We construct new models of black hole-neutron star binaries in quasiequilibrium circular orbits by solving Einstein's constraint equations in the conformal thin-sandwich decomposition together with the relativistic equations of hydrostationary equilibrium. We adopt maximal slicing, assume spatial conformal flatness, and impose equilibrium boundary conditions on an excision surface (i.e., the apparent horizon) to model the black hole. In our previous treatment we adopted a "leading-order" approximation for a parameter related to the black-hole spin in these boundary conditions to construct approximately nonspinning black holes. Here we improve on the models by computing the black hole's quasilocal spin angular momentum and setting it to zero. As before, we adopt a polytropic equation of state with adiabatic index Gamma=2 and assume the neutron star to be irrotational. In addition to recomputing several sequences for comparison with our earlier results, we study a wider range of neutron star masses and binary mass ratios. To locate the innermost stable circular orbit we search for turning points along both the binding energy and total angular momentum curves for these sequences. Unlike for our previous approximate boundary condition, these two minima now coincide. We also identify the formation of cusps on the neutron star surface, indicating the onset of tidal disruption. Comparing these two critical binary separations for different mass ratios and neutron star compactions we distinguish those regions that will lead to a tidal disruption of the neutron star from those that will result in the plunge into the black hole of a neutron star more or less intact, albeit distorted by tidal forces.

We construct quasiequilibrium sequences of black hole-neutron star binaries in general relativity. We solve Einstein's constraint equations in the conformal thin-sandwich formalism, subject to black hole boundary conditions imposed on the surface of an excised sphere, together with the relativistic equations of hydrostatic equilibrium. In contrast to our previous calculations we adopt a flat spatial background geometry and do not assume extreme mass ratios. We adopt a Gamma=2 polytropic equation of state and focus on irrotational neutron star configurations as well as approximately nonspinning black holes. We present numerical results for ratios of the black hole's irreducible mass to the neutron star's ADM mass in isolation of M_{irr}^{BH}/M_{ADM,0}^{NS} = 1, 2, 3, 5, and 10. We consider neutron stars of baryon rest mass M_B^{NS}/M_B^{max} = 83% and 56%, where M_B^{max} is the maximum allowed rest mass of a spherical star in isolation for our equation of state. For these sequences, we locate the onset of tidal disruption and, in cases with sufficiently large mass ratios and neutron star compactions, the innermost stable circular orbit. We compare with previous results for black hole-neutron star binaries and find excellent agreement with third-order post-Newtonian results, especially for large binary separations. We also use our results to estimate the energy spectrum of the outgoing gravitational radiation emitted during the inspiral phase for these binaries.

We present new sequences of general relativistic, quasiequilibrium black hole-neutron star binaries. We solve for the gravitational field in the conformal thin-sandwich decomposition of Einstein's field equations, coupled to the equations of relativistic hydrostatic equilibrium for a perfect fluid. We account for the black hole by solving these equations in the background metric of a Schwarzschild black hole whose mass is much larger than that of the neutron star. The background metric is treated in Kerr-Schild as well as isotropic coordinates. For the neutron star, we assume a polytropic equation of state with adiabatic index Gamma=2, and solve for both irrotational and corotational configurations. By comparing the results of irrotational and synchronized configurations with the same background metric, we conclude that the effect of the rotation on the location of tidal break-up is only on the order of a few percent. The different choices in the background also lead to differences of order a few percent, which may be an indication of the level to which these configurations approximate quasiequilibrium.

Many of the recent numerical simulations of binary black holes in vacuum adopt the moving puncture approach. This successful approach avoids the need to impose numerical excision of the black hole interior and is easy to implement. Here we wish to explore how well the same approach can be applied to moving black hole punctures in the presence of relativistic hydrodynamic matter. First, we evolve single black hole punctures in vacuum to calibrate our BSSN (Baumgarte-Shapiro-Shibata-Nakamura) implementation and to confirm that the numerical solution for the exterior spacetime is invariant to any ``junk'' (i.e., constraint-violating) initial data employed in the black hole interior. Then we focus on relativistic Bondi accretion onto a moving puncture Schwarzschild black hole as a numerical testbed for our high-resolution shock-capturing relativistic hydrodynamics scheme. We find that the hydrodynamical equations can be evolved successfully in the interior without imposing numerical excision. These results help motivate the adoption of the moving puncture approach to treat the binary black hole-neutron star problem using conformal thin-sandwich initial data.

We follow the inspiral and merger of equal-mass black holes (BHs) by the moving puncture technique and demonstrate that both the exterior solution and the asymptotic gravitational waveforms are unchanged when the initial interior solution is replaced by constraint-violating ``junk'' initial data. We apply this result to evolve conformal thin-sandwich (CTS) binary BH initial data by filling their excised interiors with arbitrary, but smooth, initial data and evolving with standard puncture gauge choices. The waveforms generated for both puncture and filled-CTS initial data are remarkably similar, and there are only minor differences between irrotational and corotational CTS BH binaries. Even the interior solutions appear to evolve to the same constraint-satisfying solution at late times, independent of the initial data.

We construct quasiequilibrium sequences of black hole-neutron star binaries for arbitrary mass ratios by solving the constraint equations of general relativity in the conformal thin-sandwich decomposition. We model the neutron star as a stationary polytrope satisfying the relativistic equations of hydrodynamics, and account for the black hole by imposing equilibrium boundary conditions on the surface of an excised sphere (the apparent horizon). In this paper we focus on irrotational configurations, meaning that both the neutron star and the black hole are approximately nonspinning in an inertial frame. We present results for a binary with polytropic index n=1, mass ratio M_{irr}^{BH}/M_{B}^{NS}=5 and neutron star compaction M_{ADM,0}^{NS}/R_0=0.0879, where M_{irr}^{BH} is the irreducible mass of the black hole, M_{B}^{NS} the neutron star baryon rest-mass, and M_{ADM,0}^{NS} and R_0 the neutron star Arnowitt-Deser-Misner mass and areal radius in isolation, respectively. Our models represent valid solutions to Einstein's constraint equations and may therefore be employed as initial data for dynamical simulations of black hole-neutron star binaries.

We have calculated the first dynamical evolutions of merging black hole-neutron star binaries that treat the combined spacetime in a nonperturbative general relativistic framework. Using the conformal flatness approximation, we have studied how the location of the tidal disruption radius with respect to the the black hole horizon and innermost stable circular orbit (ISCO) affects the qualitative evolution of the system. Based on simple arguments, we show that for a binary mass ratio q>~0.24, tidal disruption occurs outside the ISCO, while the opposite is true for q<~0.24. When tidal disruption occurs sufficiently far outside the ISCO, mass is transferred unstably from the neutron star to the black hole, resulting in the complete disruption of the neutron star. When tidal disruption occurs slightly within the ISCO, we find that some of the mass forms an extremely hot disk around the black hole. The resulting configurations in this case are excellent candidates for the progenitors of short-hard gamma ray bursts.

The recent localization of some short-hard gamma ray bursts (GRBs) in galaxies with low star formation rates has lent support to the suggestion that these events result from compact object binary mergers. We discuss how new simulations in general relativity are helping to identify the central engine of short-hard GRBs. Motivated by our latest relativistic black hole-neutron star merger calculations, we discuss a scenario in which these events may trigger short-hard GRBs, and compare this model to competing relativistic models involving binary neutron star mergers and the delayed collapse of hypermassive neutron stars. Distinguishing features of these models may help guide future GRB and gravitational wave observations to identify the nature of the sources.

We calculate the first dynamical evolutions of merging black hole-neutron star binaries that construct the combined black hole-neutron star spacetime in a general relativistic framework. We treat the metric in the conformal flatness approximation, and assume that the black hole mass is sufficiently large compared to that of the neutron star so that the black hole remains fixed in space. Using a spheroidal spectral methods solver, we solve the resulting field equations for a neutron star orbiting a Schwarzschild black hole. The matter is evolved using a relativistic, Lagrangian, smoothed particle hydrodynamics (SPH) treatment. We take as our initial data recent quasiequilibrium models for synchronized neutron star polytropes generated as solutions of the conformal thin-sandwich (CTS) decomposition of the Einstein field equations. We are able to construct from these models relaxed SPH configurations whose profiles show good agreement with CTS solutions. Our adiabatic evolution calculations for neutron stars with low compactness show that mass transfer, when it begins while the neutron star orbit is still outside the innermost stable circular orbit, is more unstable than is typically predicted by analytical formalisms. This dynamical mass loss is found to be the driving force in determining the subsequent evolution of the binary orbit and the neutron star, which typically disrupts completely within a few orbital periods. The majority of the mass transferred onto the black hole is accreted promptly; a significant fraction (~30%) of the mass is shed outward as well, some of which will become gravitationally unbound and ejected completely from the system. The remaining portion forms an accretion disk around the black hole, and could provide the energy source for short-duration gamma ray bursts.

We review recent progress in understanding the hydrodynamics of compact binary mergers using relativistic smoothed particle hydrodynamics (SPH) codes. Recent results are discussed for both double neutron stars and black hole - neutron star binaries. The gravitational wave signals from these mergers contain a lot of information about the masses and radii of the stars and the equation of state of the neutron star fluid. For black hole - neutron star mergers, the final fate of the fluid and the gravitational wave emission depend strongly on the mass ratio and the black hole spin.

We calculate the evolution of planets undergoing a strong tidal encounter using smoothed particle hydrodynamics (SPH), for a range of periastron separations. We find that outside the Roche limit, the evolution of the planet is well-described by the standard model of linear, non-radial, adiabatic oscillations. If the planet passes within the Roche limit at periastron, however, mass can be stripped from it, but in no case do we find enough energy transferred to the planet to lead to complete disruption. In light of the three new extrasolar planets discovered with periods shorter than two days, we argue that the shortest-period cases observed in the period-mass relation may be explained by a model whereby planets undergo strong tidal encounters with stars, after either being scattered by dynamical interactions into highly eccentric orbits, or tidally captured from nearly parabolic orbits. Although this scenario does provide a natural explanation for the edge found for planets at twice the Roche limit, it does not explain how such planets will survive the inevitable expansion that results from energy injection during tidal circularization.

We present the first results from our new general relativistic, Lagrangian hydrodynamics code, which treats gravity in the conformally flat (CF) limit. The evolution of fluid configurations is described using smoothed particle hydrodynamics (SPH), and the elliptic field equations of the CF formalism are solved using spectral methodsin spherical coordinates. The code was tested on models for which the CF limit is exact, finding good agreement with the classical Oppenheimer-Volkov solution for a relativistic static spherical star as well as the exact semi-analytic solution for a collapsing spherical dust cloud. By computing the evolution of quasi-equilibrium neutron star binary configurations in the absence of gravitational radiation backreaction, we have confirmed that these configurations can remain dynamically stable all the way to the development of a cusp. With an approximate treatment of radiation reaction, we have calculated the complete merger of an irrotational binary configuration from the innermost point on an equilibrium sequence through merger and remnant formation and ringdown, finding good agreement withprevious relativistic calculations. In particular, we find that mass loss is highly suppressed by relativistic effects, but that, for a reasonably stiff neutron star equation of state, the remnant is initially stable against gravitational collapse because of its strong differential rotation. The gravity wave signal derived from our numerical calculation has an energy spectrum which matches extremely well with estimates based solely on quasi-equilibrium results, deviating from the Newtonian power-law form at frequencies below 1 kHz, i.e., within the reach of advanced interferometric detectors.

We have designed and tested a new relativistic Lagrangian hydrodynamics code, which treats gravity in the conformally flat approximation to general relativity. We have tested the resulting code extensively, finding that it performs well for calculations of equilibrium single-star models, collapsing relativistic dust clouds, and quasi-circular orbits of equilibrium solutions. By adding in a radiation reaction treatment, we compute the full evolution of a coalescing binary neutron star system. We find that the amount of mass ejected from the system, much less than a percent, is greatly reduced by the inclusion of relativistic gravitation. The gravity wave energy spectrum shows a clear divergence away from the Newtonian point-mass form, consistent with the form derived from relativistic quasi-equilibrium fluid sequences.

The recent passage of Jupiter by the quasar QSO J0842+1835 at a separation of 3.7 arcminutes on September 8, 2002, combined with recent advances in interferometric radio timing, has allowed for the first measurement of higher-order post-Newtonian terms in the Shapiro time delay which depend linearly on the velocity of the gravitating body. Claims have been made that these measurements also allow for the measurement of the propagation speed of the gravitational force. This conclusion disagrees with recent calculations done in the parameterized post-Newtonian (PPN) model, which find no dependence of the velocity-dependent terms in the time delay on the speed of gravity to the stated order. Here, to test out these claims and counterclaims, we calculate the time delay in the limit of an instantaneous gravitational force, and find that the velocity-dependent terms are in complete agreement with previous PPN calculations, with no dependence on the speed of gravity. We conclude that the speed of gravity cannot be determined by measuring these terms in the Shapiro time delay, and suggest a reason why other groups mistakenly came to the opposite conclusion.

Coalescing binary neutron stars (NS) are expected to be an important source of gravitational waves (GW) detectable by laser interferometers. We present here a simple method for determining the compactness ratio M/R of NS based on the observed deviation of the GW energy spectrum from point-mass behavior at the end of an inspiral event. Our method is based on the properties of quasi-equilibrium binary NS sequences and does not require the computation of the full GW signal h(t). Combined with the measurement of the NS masses from the GW signal during inspiral, the determination of M/R will allow very strong constraints to be placed on the equation of state of nuclear matter at high densities.

Using our new post-Newtonian (PN) smoothed particle hydrodynamics (SPH) code, we have studied numerically the mergers of neutron star binaries with irrotational initial configurations. Here we describe a new method for constructing numerically accurate initial conditions for irrotational binary systems with circular orbits in PN gravity. We then compute the 3D hydrodynamic evolution of these systems until the two stars have completely merged, and we determine the corresponding GW signals. We present results for systems with different binary mass ratios, and for neutron stars represented by polytropes with $\Gamma=2$ or $\Gamma=3$. Compared to mergers of corotating binaries, we find that irrotational binary mergers produce similar peak GW luminosities, but they shed almost no mass at all to large distances. The dependence of the GW signal on numerical resolution for calculations performed with N>10^5 SPH particles is extremely weak, and we find excellent agreement between runs utilizing N=10^5 and N=10^6 SPH particles (the largest SPH calculation ever performed to study such irrotational binary mergers). We also compute GW energy spectra based on all calculations reported here and in our previous works. We find that PN effects lead to clearly identifiable features in the GW energy spectrum of binary neutron star mergers, which may yield important information about the nuclear equation of state at extreme densities.

The final burst of gravitational radiation emitted by coalescing binary neutron stars carries direct information about the neutron star fluid, and, in particular, about the equation of state of nuclear matter at extreme densities. The final merger may also be accompanied by a detectable electromagnetic signal, such as a gamma-ray burst. In this paper, we summarize the results of theoretical work done over the past decade that has led to a detailed understanding of this hydrodynamic merger process for two neutron stars, and we discuss the prospects for the detection and physical interpretation of the gravity wave signals by ground-based interferometers such as LIGO. We also present results from our latest post-Newtonian SPH calculations of binary neutron star coalescence, using up to 10^6 SPH particles to compute with higher spatial resolution than ever before the merger of an initially irrotational system. We discuss the detectability of our calculated gravity wave signals based on power spectra.

Using our new Post-Newtonian SPH (smoothed particle hydrodynamics) code, we study the final coalescence and merging of neutron star (NS) binaries. We vary the stiffness of the equation of state (EOS) as well as the initial binary mass ratio and stellar spins. Results are compared to those of Newtonian calculations, with and without the inclusion of the gravitational radiation reaction. We find a much steeper decrease in the gravity wave peak strain and luminosity with decreasing mass ratio than would be predicted by simple point-mass formulae. For NS with softer EOS (which we model as simple $\Gamma=2$ polytropes) we find a stronger gravity wave emission, with a different morphology than for stiffer EOS (modeled as $\Gamma=3$ polytropes as in our previous work). We also calculate the coalescence of NS binaries with an irrotational initial condition, and find that the gravity wave signal is relatively suppressed compared to the synchronized case, but shows a very significant second peak of emission. Mass shedding is also greatly reduced, and occurs via a different mechanism than in the synchronized case. We discuss the implications of our results for gravity wave astronomy with laser interferometers such as LIGO, and for theoretical models of gamma-ray bursts (GRBs) based on NS mergers.

We continue our exploration of collisionally merged stars in the blue straggler region of the color-magnitude diagram. We report the results of new SPH calculations of parabolic collisions between two main-sequence stars, with the initial structure and composition profiles of the parent stars having been determined from stellar evolution calculations. Parallelization of the SPH code has permitted much higher numerical resolution of the hydrodynamics. We also present evolutionary tracks for the resulting collision products, which emerge as rapidly rotating blue stragglers. The rotating collision products are brighter, bluer and remain on the main sequence longer than their non-rotating counterparts. In addition, they retain their rapid rotation rates throughout their main sequence lifetime. Rotationally-induced mixing strongly affects the evolution of the collision products, although it is not sufficient to mix the entire star. We discuss the implications of these results for studies of blue straggler populations in clusters. This work shows that off-axis collision products cannot become blue stragglers unless they lose a large fraction of their initial angular momentum. The mechanism for this loss is not apparent, although some possibilities are discussed.

We present the first results from our Post-Newtonian (PN) Smoothed Particle Hydrodynamics (SPH) code, which has been used to study the coalescence of binary neutron star (NS) systems. The Lagrangian particle-based code incorporates consistently all lowest-order (1PN) relativistic effects, as well as gravitational radiation reaction, the lowest-order dissipative term in general relativity. We test our code on sequences of single NS models of varying compactness, and we discuss ways to make PN simulations more relevant to realistic NS models. We also present a PN SPH relaxation procedure for constructing equilibrium models of synchronized binaries, and we use these equilibrium models as initial conditions for our dynamical calculations of binary coalescence. Though unphysical, since tidal synchronization is not expected in NS binaries, these initial conditions allow us to compare our PN work with previous Newtonian results. We compare calculations with and without 1PN effects, for NS with stiff equations of state, modeled as polytropes with $\Gamma=3$. We find that 1PN effects can play a major role in the coalescence, accelerating the final inspiral and causing a significant misalignment in the binary just prior to final merging. In addition, the character of the gravitational wave signal is altered dramatically, showing strong modulation of the exponentially decaying waveform near the end of the merger. We also discuss briefly the implications of our results for models of gamma-ray bursts at cosmological distances.