Search

Tidal disruption events (TDEs) are an important way to probe the properties of stellar populations surrounding supermassive black holes. Observed spectra of several TDEs, such as ASASSN-14li, show high nitrogen to carbon abundance ratios, leading to questions about their progenitors. Disrupting an intermediate- or high-mass star that has undergone CNO processing, increasing the nitrogen in its core, could lead to an enhanced nitrogen TDE. Galactic nuclei present a conducive environment for high-velocity stellar collisions that can lead to high mass loss, stripping the carbon- and hydrogen-rich envelopes of the stars and leaving behind the enhanced nitrogen cores. TDEs of these stripped stars may therefore exhibit even more extreme nitrogen enhancement. Using the smoothed particle hydrodynamics code StarSmasher, we provide a parameter space study of high-velocity stellar collisions involving intermediate-mass stars, analyzing the composition of the collision products. We conclude that high-velocity stellar collisions can form products that have abundance ratios similar to those observed in the motivating TDEs. Furthermore, we show that stars that have not experienced high CNO processing can yield low-mass collision products that retain even higher nitrogen to carbon abundance ratios. We analytically estimate the mass fallback for a typical TDE of several collision products to demonstrate consistency between our models and TDE observations. Lastly, we discuss how the extended collision products, with high central to average density ratios, can be related to repeated partial TDEs like ASASSN-14ko and G objects in the Galactic Center., Comment: 26 pages, 10 figures, 3 tables, submitted to ApJ

The MeerKAT telescope represents an outstanding opportunity for radio pulsar timing science with its unique combination of a large collecting area and aperture efficiency (effective area $\sim$7500 m$^2$), system temperature ($T<20$K), high slew speeds (1-2 deg/s), large bandwidths (770 MHz at 20cm wavelengths), southern hemisphere location (latitude $\sim -30^\circ$) and ability to form up to four sub-arrays. The MeerTime project is a five-year program on the MeerKAT array by an international consortium that will regularly time over 1000 radio pulsars to perform tests of relativistic gravity, search for the gravitational wave signature induced by supermassive black hole binaries in the timing residuals of millisecond pulsars, explore the interiors of neutron stars through a pulsar glitch monitoring programme, explore the origin and evolution of binary pulsars, monitor the swarms of pulsars that inhabit globular clusters and monitor radio magnetars. In addition to these primary programmes, over 1000 pulsars will have their arrival times monitored and the data made immediately public. The MeerTime pulsar backend comprises two server-class machines each of which possess four Graphics Processing Units. Up to four pulsars can be coherently dedispersed simultaneously up to dispersion measures of over 1000 pc cm$^{-3}$. All data will be provided in psrfits format. The MeerTime backend will be capable of producing coherently dedispersed filterbank data for timing multiple pulsars in the cores of globular clusters that is useful for pulsar searches of tied array beams. All MeerTime data will ultimately be made available for public use, and any published results will include the arrival times and profiles used in the results., Comment: 15 pages, MeerKAT Science: On the Pathway to the SKA, 25-27 May, 2016, Stellenbosch, South Africa, available at: https://pos.sissa.it/277/011/pdf

We present theoretical models for stellar black hole (BH) properties in young, massive star clusters. Using a Monte Carlo code for stellar dynamics, we model realistic star clusters with $N\simeq 5\times10^5$ stars and significant binary fractions (up to 50%) with self-consistent treatments of stellar dynamics and stellar evolution. We compute the formation rates and characteristic properties of single and binary BHs for various representative ages, cluster parameters, and metallicities. Because of dynamical interactions and supernova (SN) kicks, more single BHs end up retained in clusters compared to BHs in binaries. We also find that the ejection of BHs from a cluster is a strong function of initial density. In low-density clusters (where dynamical effects are negligible), it is mainly SN kicks that eject BHs from the cluster, whereas in high-density clusters (initial central density $\rho_c(0) \sim 10^5 \, M_\odot\, {\rm pc}^{-3} $ in our models) the BH ejection rate is enhanced significantly by dynamics. Dynamical interactions of binary systems in dense clusters also modify the orbital period and eccentricity distributions while also increasing the probability of a BH having a more massive companion.

The rapid circularization and synchronization of the stellar components in an eccentric binary system at the onset of Roche lobe overflow (RLO) is a fundamental assumption common to all binary stellar evolution and population synthesis codes, even though the validity of this assumption is questionable both theoretically and observationally. Here we calculate the evolution of the orbital elements of an eccentric binary through the direct three-body integration of a massive particle ejected through the inner Lagrangian point of the donor star at periastron. The trajectory of this particle leads to three possible outcomes: direct accretion (DA) onto the companion star within a single orbit, self-accretion (SA) back onto the donor star within a single orbit, or a quasi-periodic orbit around the companion star. We calculate the secular evolution of the binary orbit in the first two cases and conclude that DA can increase or decrease the orbital semi-major axis and eccentricity, while SA always decreases the orbital both orbital elements. In cases where mass overflow contributes to circularizing the orbit, circularization can set in on timescales as short as a few per cent of the mass transfer timescale. In cases where mass overflow increases the eccentricity, the orbital evolution is governed by competition between mass overflow and tidal torques. In the absence of tidal torques, mass overflow resulting in DI can lead to substantially subsynchronously rotating donor stars. Contrary to common assumptions, DI furthermore does not always provide a strong sink of orbital angular momentum in close mass-transferring binaries; in fact we instead find that a significant part can be returned to the orbit during the particle orbit. The formulation presented here can be combined with stellar and binary evolution codes to generate a better picture of the evolution of eccentric, RLO binary star systems., Comment: 15 pages, 10 figures, Accepted for publication in ApJ

Using our recently improved Monte Carlo evolution code, we study the evolution of the binary fraction in globular clusters. In agreement with previous N-body simulations, we find generally that the hard binary fraction in the core tends to increase with time over a range of initial cluster central densities for initial binary fractions <~ 90%. The dominant processes driving the evolution of the core binary fraction are mass segregation of binaries into the cluster core and preferential destruction of binaries there. On a global scale, these effects and the preferential tidal stripping of single stars tend to roughly balance, leading to overall cluster binary fractions that are roughly constant with time. Our findings suggest that the current hard binary fraction near the half-mass radius is a good indicator of the hard primordial binary fraction. However, the relationship between the true binary fraction and the fraction of main-sequence stars in binaries (which is typically what observers measure) is non-linear and rather complicated. We also consider the importance of soft binaries, which not only modify the evolution of the binary fraction, but can drastically change the evolution of the cluster as a whole. Finally, we describe in some detail the recent addition of single and binary stellar evolution to our cluster evolution code., Comment: 8 pages, 7 figures in emulateapj format. Submitted to ApJ

We investigate the secular evolution of the orbital semi-major axis and eccentricity due to mass transfer in eccentric binaries, allowing for both mass and angular momentum loss from the system. Adopting a delta function mass transfer rate at the periastron of the binary orbit, we find that, depending on the initial binary properties at the onset of mass transfer, the orbital semi-major axis and eccentricity can either increase or decrease at a rate linearly proportional to the magnitude of the mass transfer rate at periastron. The range of initial binary mass ratios and eccentricities that leads to increasing orbital semi-major axes and eccentricities broadens with increasing degrees of mass loss from the system and narrows with increasing orbital angular momentum loss from the binary. Comparison with tidal evolution timescales shows that the usual assumption of rapid circularization at the onset of mass transfer in eccentric binaries is not justified, irrespective of the degree of systemic mass and angular momentum loss. This work extends our previous results for conservative mass transfer in eccentric binaries and can be incorporated into binary evolution and population synthesis codes to model non-conservative mass transfer in eccentric binaries., Comment: 6 pages, 3 figures, Accepted for publication in ApJ. To appear September 2009

As the closest example of a galactic nucleus, the Galactic center presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environs. Detailed studies of stellar dynamics deep in the potential well of a galaxy, with exisiting and future large ground-based telescopes, offer several exciting directions in the coming decade. First, it will be possible to obtain precision measurements of the Galaxy's central potential, providing both a unique test of General Relativity (GR) and a detection of the extended dark matter distribution that is predicted to exist around the SMBH. Tests of gravity have not previously been possible on the mass scale of a SMBH. Similarly, only upper limits on the extended matter distribution on small scales currently exist; detection of dark matter on these scales is an important test of Lambda-CDM and the detection of stellar remnants would reveal a population that may dominate the stellar dynamics on the smallest scales. Second, our detailed view of the SMBH and its local gas and stellar environment provides insight into how SMBHs at the centers of galaxies form, grow and interact with their environs as well as on the exotic processes at work in the densest stellar clusters in the Universe. The key questions, still unanswered, of when and how SMBHs formed in the early universe, and the myriad ways in which feedback from SMBHs can affect structure formation, can be informed by directly observing the physical processes operating at the SMBH., Comment: An ASTRO2010 White Paper. Animations and high resolution images can be found at http://www.astro.ucla.edu/~ghezgroup/gc/pictures/Future_GCorbits.shtml

Mergers of double neutron stars are considered the most likely progenitors for short gamma-ray bursts. Indeed such a merger can produce a black hole with a transient accreting torus of nuclear matter (Lee & Ramirez-Ruiz 2007, Oechslin & Janka 2006), and the conversion of a fraction of the torus mass-energy to radiation can power a gamma-ray burst (Nakar 2006). Using available binary pulsar observations supported by our extensive evolutionary calculations of double neutron star formation, we demonstrate that the fraction of mergers that can form a black hole -- torus system depends very sensitively on the (largely unknown) maximum neutron star mass. We show that the available observations and models put a very stringent constraint on this maximum mass under the assumption that a black hole formation is required to produce a short gamma-ray burst in a double neutron star merger. Specifically, we find that the maximum neutron star mass must be within 2 - 2.5 Msun. Moreover, a single unambiguous measurement of a neutron star mass above 2.5 Msun would exclude a black hole -- torus central engine model of short gamma-ray bursts in double neutron star mergers. Such an observation would also indicate that if in fact short gamma-ray bursts are connected to neutron star mergers, the gamma-ray burst engine is best explained by the lesser known model invoking a highly magnetized massive neutron star (e.g., Usov 1992; Kluzniak & Ruderman 1998; Dai et al. 2006; Metzger, Quataert & Thompson 2007)., Comment: 3 pages of text + 4 figures: ApJ accepted (some revisions)

We investigate the secular evolution of the orbital semi-major axis and eccentricity due to mass transfer in eccentric binaries, assuming conservation of total system mass and orbital angular momentum. Assuming a delta function mass transfer rate centered at periastron, we find rates of secular change of the orbital semi-major axis and eccentricity which are linearly proportional to the magnitude of the mass transfer rate at periastron. The rates can be positive as well as negative, so that the semi-major axis and eccentricity can increase as well as decrease in time. Adopting a delta-function mass-transfer rate of $10^{-9} M_\sun {\rm yr}^{-1}$ at periastron yields orbital evolution timescales ranging from a few Myr to a Hubble time or more, depending on the binary mass ratio and orbital eccentricity. Comparison with orbital evolution timescales due to dissipative tides furthermore shows that tides cannot, in all cases, circularize the orbit rapidly enough to justify the often adopted assumption of instantaneous circularization at the onset of mass transfer. The formalism presented can be incorporated in binary evolution and population synthesis codes to create a self-consistent treatment of mass transfer in eccentric binaries., Comment: 16 pages, 8 figures, Accepted by The Astrophysical Journal

In this paper, the second of a series, we study the stellar dynamical and evolutionary processes leading to the formation of compact binaries containing neutron stars (NSs) in dense globular clusters (GCs). For this study, 70 dense clusters were simulated independently, with a total stellar mass ~2x10^7Msun, exceeding the total mass of all dense GCs in our Galaxy. We find that, in order to reproduce the empirically derived formation rate of low-mass X-ray binaries (LMXBs), we must assume that NSs can be formed via electron-capture supernovae (ECS) with typical natal kicks smaller than in core-collapse supernovae. Our results explain the observed dependence of the number of LMXBs on ``collision number'' as well as the large scatter observed between different GCs. We predict that the number of quiescent LMXBs in different GCs should not have a strong metallicity dependence. In our cluster model the following mass-gaining events create populations of MSPs that do not match the observations: (i) accretion during a common envelope event with a NS formed through ECS, and (ii) mass transfer (MT) from a WD donor. Some processes lead only to a mild recycling. In addition, for MSPs, we distinguish low-magnetic-field (long-lived) and high-magnetic-field (short-lived) populations. With this distinction and by considering only those mass-gaining events that appear to lead to NS recycling, we obtain good agreement of our models with the numbers and characteristics of observed MSPs in 47 Tuc and Terzan 5, as well as with the cumulative statistics for MSPs detected in GCs of different dynamical properties. We find that significant production of merging double NSs potentially detectable as short gamma-ray bursts occurs only in very dense, most likely core-collapsed GCs. (abridged), Comment: 25 pages, 7 figures, 12 tables, MNRAS accepted

We consider the formation of double white dwarfs (DWDs) through dynamical interactions in globular clusters. Such interactions can give rise to eccentric DWDs, in contrast to the exclusively circular population expected to form in the Galactic disk. We show that for a 5-year Laser Interferometer Space Antenna (LISA) mission and distances as far as the Large Magellanic Cloud, multiple harmonics from eccentric DWDs can be detected at a signal-to-noise ratio higher than 8 for at least a handful of eccentric DWDs, given their formation rate and typical lifetimes estimated from current cluster simulations. Consequently the association of eccentricity with stellar-mass LISA sources does not uniquely involve neutron stars, as is usually assumed. Due to the difficulty of detecting (eccentric) DWDs with present and planned electromagnetic observatories, LISA could provide unique dynamical identifications of these systems in globular clusters., Comment: Published in ApJ 665, L59

This article was co-authored by all invited speakers at the Joint Discussion on `Neutron Stars and Black Holes in Star Clusters,' which took place during the IAU General Assembly in Prague, Czech Republic, on August 17 and 18, 2006. Each section presents a short summary of recent developments in a key area of research, incorporating the main ideas expressed during the corresponding panel discussion at the meeting. Our meeting, which had close to 300 registered participants, was broadly aimed at the large community of astronomers around the world working on the formation and evolution of compact objects and interacting binary systems in dense star clusters, such as globular clusters and galactic nuclei. Great advances have occurred in this field during the past few years, including the introduction of fundamentally new theoretical paradigms for the formation and evolution of compact objects in binaries as well as countless new discoveries by astronomers that have challenged many accepted models. Some of the highlights include: a nearly complete census of all the millisecond pulsars in 47 Tuc; first detections of many new radio pulsars in other clusters, particularly Terzan 5; detailed studies of X-ray binary populations and their luminosity functions in many galaxies and extragalactic globular clusters; increasing evidence for intermediate-mass black holes in clusters and greatly improved theoretical understanding of their possible formation processes., Comment: 29 pages, proceedings of IAU JD06, to appear in Highlights of Astronomy, Volume 14, ed. K. A. van der Hucht

In this paper, the first of a series, we study the stellar dynamical and evolutionary processes leading to the formation of compact binaries containing white dwarfs in dense globular clusters. We examine the processes leading to the creation of X-ray binaries such as cataclysmic variables and AM CVn systems. Using numerical simulations, we identify the dominant formation channels and we predict the expected numbers and characteristics of detectable systems, emphasizing how the cluster sources differ from the field population. We explore the dependence of formation rates on cluster properties and we explain in particular why the distribution of cataclysmic variables has only a weak dependence on cluster density. We also discuss the frequency of dwarf nova outbursts in globular clusters and their connection with moderately strong white dwarf magnetic fields. We examine the rate of Type Ia supernovae via both single and double degenerate channels in clusters and we argue that those rates may contribute to the total SN Ia rate in elliptical galaxies. Considering coalescing white dwarf binaries we discuss possible constraints on the common envelope evolution of their progenitors and we derive theoretical expectations for gravitational wave detection by LISA., Comment: 17 pages (MNRAS style), 18 figures, submitted to MNRAS (the version with full resolution figures is available at http://www.cita.utoronto.ca/~nata/lit.html)

We present a comprehensive description of the population synthesis code StarTrack. The original code has been significantly modified and updated. Special emphasis is placed here on processes leading to the formation and further evolution of compact objects (white dwarfs, neutron stars, and black holes). Both single and binary star populations are considered. The code now incorporates detailed calculations of all mass-transfer phases, a full implementation of orbital evolution due to tides, as well as the most recent estimates of magnetic braking. This updated version of StarTrack can be used for a wide variety of problems, with relevance to many current and planned observatories, e.g., studies of X-ray binaries (Chandra, XMM-Newton), gravitational radiation sources (LIGO, LISA), and gamma-ray burst progenitors (HETE-II, Swift). The code has already been used in studies of Galactic and extra-galactic X-ray binary populations, black holes in young star clusters, Type Ia supernova progenitors, and double compact object populations. Here we describe in detail the input physics, we present the code calibration and tests, and we outline our current studies in the context of X-ray binary populations., Comment: 40 pages, 14 figures, accepted to ApJS (major extensions/revisions)

(abridged) We report the results of SPH calculations of parabolic collisions between a subgiant or slightly evolved red-giant star and a neutron star (NS). Such collisions are likely to form ultracompact X-ray binaries (UCXBs) observed today in old globular clusters. In particular, we compute collisions of a 1.4 Msun NS with realistically modelled parent stars of initial masses 0.8 and 0.9 Msun, each at three different evolutionary stages (corresponding to three different radii R). The distance of closest approach for the initial orbit varies from 0.04 R (nearly head-on) to 1.3 R (grazing). These collisions lead to the formation of a tight binary, composed of the NS and the subgiant or red-giant core, embedded in an extremely diffuse common envelope (CE) typically of mass ~0.1 to 0.3 Msun. Our calculations follow the binary for many hundreds of orbits, ensuring that the orbital parameters we determine at the end of the calculations are close to final. Some of the fluid initially in the envelope of the (sub)giant, from 0.003 to 0.023 Msun in the cases we considered, is left bound to the NS. The eccentricities of the resulting binaries range from about 0.2 for our most grazing collision to about 0.9 for the nearly head-on cases. In almost all the cases we consider, gravitational radiation alone will cause sufficiently fast orbital decay to form a UCXB within a Hubble time, and often on a much shorter timescale. Our hydrodynamics code implements the recent SPH equations of motion derived with a variational approach by Springel & Hernquist and by Monaghan. Numerical noise is reduced by enforcing an analytic constraint equation that relates the smoothing lengths and densities of SPH particles. We present tests of these new methods to help demonstrate their improved accuracy., Comment: 41 pages, 17 figures, accepted by ApJ

While many observed ultra-luminous X-ray sources (ULXs, Lx > 10^39 erg s^-1) could be extragalactic X-ray binaries (XRBs) emitting close to the Eddington limit, the highest-luminosity ULXs (Lx > 3x10^39 erg s^-1) exceed the isotropic Eddington luminosity for even high-stellar-mass accreting black hole XRBs. It has been suggested that these highest-luminosity ULXs may contain accreting intermediate-mass black hole (IMBH) binaries. We consider this hypothesis for dense, young (about 100 Myr) stellar clusters where we assume that a 50-500 solar mass central IMBH has formed through runaway growth of a massive star. Using numerical simulations of the dynamics and evolution of the central black hole's captured companions, we obtain estimates of the incidence of mass transfer phases and possible ULX activity throughout the IMBH's evolutionary history. We find that, although it is common for the central black hole to acquire binary companions, there is a very low probability that these interacting binaries will become observable ULX sources., Comment: 11 pages, 8 figures, submitted to ApJ

Using an updated population synthesis code we study the formation and evolution of black holes (BHs) in young star clusters following a massive starburst. This study continues and improves on the initial work described by Belczynski, Sadowski & Rasio (2004). In our new calculations we account for the possible ejections of BHs and their progenitors from clusters because of natal kicks imparted by supernovae and recoil following binary disruptions. The results indicate that the properties of both retained BHs in clusters and ejected BHs (forming a field population) depend sensitively on the depth of the cluster potential. In particular, most BHs ejected from binaries are also ejected from clusters with central escape speeds Vesc < 100 km/s. Conversely, most BHs remaining in binaries are retained by clusters with Vesc > 50 km/s. BHs from single star evolution are also affected significantly: about half of the BHs originating from primordial single stars are ejected from clusters with Vesc < 50 km/s. Our results lay a foundation for theoretical studies of the formation of BH X-ray binaries in and around star clusters, including possible ultra-luminous sources, as well as merging BH--BH binaries detectable with future gravitational-wave observatories., Comment: 35 pages, 8 tables, 17 figures; resubmitted to ApJ (revised version)

One of the most remarkable properties of extrasolar planets is their high orbital eccentricities. Observations have shown that at least 20% of these planets, including some with particularly high eccentricities, are orbiting a component of a wide binary star system. The presence of a distant binary companion can cause significant secular perturbations to the orbit of a planet. In particular, at high relative inclinations, a planet can undergo a large-amplitude eccentricity oscillation. This so-called "Kozai mechanism" is effective at a very long range, and its amplitude is purely dependent on the relative orbital inclination. In this paper, we address the following simple question: assuming that every host star with a detected giant planet also has a (possibly unseen, e.g., substellar) distant companion, with reasonable distributions of orbital parameters and masses, how well could secular perturbations reproduce the observed eccentricity distribution of planets? Our calculations show that the Kozai mechanism consistently produces an excess of planets with very high (e >0.6) and very low (e < 0.1) eccentricities. The paucity of near-circular orbits in the observed sample cannot be explained solely by the Kozai mechanism, because, even with high enough inclinations, the Kozai mechanism often fails to produce significant eccentricity perturbations when there are other competing sources of orbital perturbations on secular timescales, such as general relativity. On the other hand, the Kozai mechanism can produce many highly eccentric orbits. Indeed the overproduction of high eccentricities observed in our models could be combined with plausible circularizing mechanisms (e.g., friction from residual gas) to create more intermediate eccentricities (e=0.1-0.6)., Comment: 24 pages, 6 figures, ApJ, in press, minor changes to reflect the accepted version

Bright, ultracompact X-ray binaries observed in dense star clusters, such as Galactic globular clusters, must have formed relatively recently, since their lifetimes as persistent bright sources are short (e.g., ~10^8 yr above 10^36 erg/s for a 1.4 Msun neutron star accreting from a degenerate helium companion with an initial mass of ~0.2 Msun). Therefore, we can use the present conditions in a cluster core to study possible dynamical formation processes for these sources. Here we show that direct physical collisions between neutron stars and red giants can provide a sufficient formation rate to explain the observed numbers of bright sources. These collisions produce tight, eccentric neutron star -- white dwarf binaries that decay to contact by gravitational radiation on timescales ~10^6-10^10 yr, usually shorter and often much shorter than the cluster age., Comment: 4 pages, 1 figure, in emulateapj style. ApJ Letters, in press, includes corrections in proofs

We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code (StarTrack), and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing 3-body and 4-body interactions (Fewbody). We find that the combination of stellar evolution and dynamical interactions (binary-single and binary-binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster like 47 Tuc, assuming an initial binary fraction of 100%, is only about 5-10%. We show that this is in good agreement with recent HST observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modeling of globular clusters, suggesting that ``realistic models'' should incorporate much larger initial binary fractions than has usually been done in the past., Comment: 14 pages, MNRAS style, 9 figures, MNRAS accepted

The population of compact binaries in dense stellar systems is affected strongly by frequent dynamical interactions between stars and their interplay with the stellar evolution. In this contribution, we consider these effects on binaries with a white dwarf accretor, in particular cataclysmic variables and AM CVns. We examine which processes can successfully lead to the creation of such X-ray binaries. Using numerical simulations, we identify predominant formation channels and predict the expected numbers of detectable systems. We discuss also why the distribution of cataclysmic variables has a weaker dependence upon the cluster density than the distribution of quiescent low-mass X-ray binaries and why dwarf nova outbursts may not occur among globular cluster cataclysmic variables., Comment: 8 pages, 4 figures (needs AIP style file). To appear in the conference proceedings of `Interacting Binaries: Accretion, Evolution & Outcomes' (Cefalu, July 4-10 2004)

In dense stellar systems the frequent dynamical interactions between stars play a crucial role in the formation and evolution of compact binaries. We study these processes using a novel approach combining a state-of-the-art binary population synthesis code with a simple treatment of dynamical interactions in dense star cluster cores. Here we focus on the dynamical and evolutionary processes leading to the formation of compact binaries containing white dwarfs in dense globular clusters. We demonstrate that dynamics can increase by factors ~ 2 - 100 the production rates of interesting binaries such as cataclysmic variables, ``nonflickerers'' (He white dwarfs with a heavier dark companion), merging white dwarf binaries with total masses above the Chandrasekhar limit, and white dwarf binaries emitting gravitational waves in the LISA band., Comment: to appear in Proceedings of IAU Colloquium 194 "Compact binaries in the Galaxy and beyond" (Rev. Mex. A&A Conf. Series), eds. G. Tovmassian and E. Sion

Physical collisions between stars occur frequently in dense star clusters, either via close encounters between two single stars, or during strong dynamical interactions involving binary stars. Here we study stellar collisions that occur during binary-single and binary-binary interactions, by performing numerical scattering experiments. Our results include cross sections, branching ratios, and sample distributions of parameters for various outcomes. For interactions of hard binaries containing main-sequence stars, we find that the normalized cross section for at least one collision to occur (between any two of the four stars involved) is essentially unity, and that the probability of collisions involving more than two stars is significant. Hydrodynamic calculations have shown that the effective radius of a collision product can be 2-30 times larger than the normal main-sequence radius for a star of the same total mass. We study the effect of this expansion, and find that it increases the probability of further collisions considerably. We discuss these results in the context of recent observations of blue stragglers in globular clusters with masses exceeding twice the main-sequence turnoff mass. We also present Fewbody, a new, freely available numerical toolkit for simulating small-N gravitational dynamics that is particularly suited to performing scattering experiments., Comment: Accepted for publication in MNRAS; 29 pages, 22 figures, 7 tables

We study the evolution of binary stars in globular clusters using a novel approach combining a state-of-the-art population synthesis code with a simple treatment of dynamical interactions in the dense cluster core. We find that the combination of stellar evolution and dynamical interactions (binary-single and binary-binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster like 47 Tuc, assuming an initial binary fraction of 100%, is only about 5%. We show that this is in good agreement with recent HST observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modeling of globular clusters, suggesting that ``realistic models'' should incorporate much larger initial binary fractions than has usually been done in the past., Comment: 4 pages, ApJ style, 2 figures, 2 tables, ApJL submitted

We study the dynamical evolution of globular clusters using our 2D Monte Carlo code with the inclusion of primordial binary interactions for equal-mass stars. We use approximate analytical cross sections for energy generation from binary-binary and binary-single interactions. After a brief period of slight contraction or expansion of the core over the first few relaxation times, all clusters enter a much longer phase of stable "binary burning" lasting many tens of relaxation times. The structural parameters of our models during this phase match well those of most observed globular clusters. At the end of this phase, clusters that have survived tidal disruption undergo deep core collapse, followed by gravothermal oscillations. Our results clearly show that the presence of even a small fraction of binaries in a cluster is sufficient to support the core against collapse significantly beyond the normal core collapse time predicted without the presence of binaries. For tidally truncated systems, collapse is easily delayed sufficiently that the cluster will undergo complete tidal disruption before core collapse. As a first step toward the eventual goal of computing all interactions exactly using dynamical three- and four-body integration, we have incorporated an exact treatment of binary-single interactions in our code. We show that results using analytical cross sections are in good agreement with those using exact three-body integration, even for small binary fractions where binary-single interactions are energetically most important., Comment: Accepted for publication in ApJ. Minor changes to reflect accepted version. 28 pages, 17 figures; some figures low resolution. Full resolution paper available at http://www.mit.edu/~fregeau/paper3.pdf

We have calculated the evolution of 60 model binary systems consisting of helium stars in the mass range of M_He= 2.5-6Msun with a 1.4Msun neutron star companion to investigate the formation of double neutron star systems.Orbital periods ranging from 0.09 to 2 days are considered, corresponding to Roche lobe overflow starting from the helium main sequence to after the ignition of carbon burning in the core. We have also examined the evolution into a common envelope phase via secular instability, delayed dynamical instability, and the consequence of matter filling the neutron star's Roche lobe. The survival of some close He-star neutron-star binaries through the last mass transfer episode (either dynamically stable or unstable mass transfer phase) leads to the formation of extremely short-period double neutron star systems (with P<~0.1days). In addition, we find that systems throughout the entire calculated mass range can evolve into a common envelope phase, depending on the orbital period at the onset of mass transfer. The critical orbital period below which common envelope evolution occurs generally increases with M_He. In addition, a common envelope phase may occur during a short time for systems characterized by orbital periods of 0.1-0.5 days at low He-star masses (~ 2.6-3.3Msun). The existence of a short-period population of double neutron stars increases the predicted detection rate of inspiral events by ground-based gravitational-wave detectors and impacts their merger location in host galaxies and their possible role as gamma-ray burst progenitors. We use a set of population synthesis calculations and investigate the implications of the mass-transfer results for the orbital properties of DNS populations., Comment: 30 pages, Latex (AASTeX), 1 table, 8 figures. To appear in ApJ, v592 n1 July 20, 2003

We present the results of a new study of mass segregation in two-component star clusters, based on a large number of numerical N-body simulations using our recently developed dynamical Monte Carlo code. Specifically, we follow the dynamical evolution of clusters containing stars with individual masses m_1 as well as a tracer population of objects with individual masses m_2=\mu m_1, using N=10^5 total stars. For heavy tracers, which could represent stellar remnants such as neutron stars or black holes in a globular cluster, we characterize in a variety of ways the tendency for these objects to concentrate in or near the cluster core. In agreement with simple theoretical arguments, we find that the characteristic time for this mass segregation process varies as 1/\mu. For models with very light tracers (\mu <~ 10^-2), which could represent free-floating planets or brown dwarfs, we find the expected depletion of light objects in the cluster core, but also sometimes a significant enhancement in the halo. Using these results we estimate the optical depth to gravitational microlensing by planetary mass objects or brown dwarfs in typical globular clusters. For some initial conditions, the optical depth in the halo due to very low-mass objects could be much greater than that of luminous stars. If we apply our results to M22, using the recent null detection of Sahu, Anderson, & King (2001), we find an upper limit of ~25% at the 63% confidence level for the current mass fraction of M22 in the form of very low-mass objects., Comment: Accepted for publication in ApJ. Minor revisions reflecting the new results of Sahu et al. on M22. 13 pages in emulateapj style, including 9 figures and 3 tables

The radio millisecond pulsar PSR B1620-26 is part of an extraordinary triple star system in the globular cluster M4. The inner companion to the neutron star is thought to be a white dwarf of mass m1 ~0.3 Msun in an orbit of period ~0.5 year. The nature and orbital characteristics of the second, more distant companion, have remained a mystery for many years. A theoretical analysis of the latest available radio pulsar timing data is presented here, allowing us to determine approximately the mass and orbital parameters of the second companion. Remarkably, the current best-fit parameters correspond to a second companion of **planetary mass**, with m2 ~0.007 Msun, in an orbit of eccentricity e2 ~0.45 and with a large semimajor axis a2 ~60 AU. The short dynamical lifetime of this very wide triple in M4 suggests that large numbers of such planets must be present in globular clusters. We also address the question of the anomalously high eccentricity of the inner binary pulsar. While this eccentricity could have been induced during the same dynamical interaction that created the triple, we find that it could also naturally arise from long-term secular perturbation effects in the triple, combining the general relativistic precession of the inner orbit with the Newtonian gravitational perturbation by the outer planet., Comment: 13 pages, to appear in Evolution of Binary and Multiple Star Systems, a Meeting in Celebration of Peter Eggleton's 60th Birthday, Bormio, Italy, ASP Conference Series, eds. P. Podsiadlowski et al

Globular clusters house a population of compact binaries that will be interesting gravitational wave sources for LISA. We provide estimates for the numbers of sources of several categories and discuss the sensitivity of LISA to detecting these sources. The estimated total number of detectable sources ranges from about 10 to about 1000 with gravitational wave frequencies above 1 mHz. These sources are typically undetectable by any other means and thus offer an opportunity for doing true gravitational-wave astronomy. The detection of these sources would provide information about both binary star evolution and the dynamics of globular clusters., Comment: Contribution to Proceedings of 3rd LISA Symposium 7 pages, added references

We have developed a new parallel supercomputer code based on Henon's Monte Carlo method for simulating the dynamical evolution of globular clusters. This new code allows us to calculate the evolution of a cluster containing a realistic number of stars (N ~ 10^5 - 10^6) in about a day of computing time. The discrete, star-by-star representation of the cluster in the simulation allows us to treat naturally a number of important processes, including single and binary star evolution, all dynamical interactions of single stars and primordial binaries, and tidal interactions with the Galaxy., Comment: 7 pages, with 2 postscript figures, to appear in Dynamics of Star Clusters and the Milky Way, eds. R. Spurzem et al. (ASP Conference Series)

Coalescing compact binaries are thought to be involved in a wide variety of astrophysical phenomena. In particular, they are important sources of gravitational radiation for both ground-based and space-based laser-interferometer detectors, and they may be sources of supernova explosions or gamma-ray bursts. Mergers of two white dwarfs may produce neutron stars with peculiar properties, including perhaps millisecond radio pulsars sometimes accompanied by planets (as observed in PSR 1257+12). According to a widely held belief, the coalescence of two neutron stars should produce a rapidly rotating black hole surrounded by an accretion disk or torus, but this is by no means certain. This review paper focuses on the final hydrodynamic coalescence and merger of double neutron stars and double white dwarfs, and addresses the question of the nature of the final merger products., Comment: 14 pages, 1 figure, lecture presented at the NATO ASI on ''The Neutron Star - Black Hole Connection'' in Elounda, Crete, June 1999

We present the results of Monte Carlo simulations for the dynamical evolution of star clusters containing two stellar populations with individual masses m1 and m2 > m1, and total masses M1 and M2 < M1. We use both King and Plummer model initial conditions and we perform simulations for a wide range of individual and total mass ratios, m2/m1 and M2/M1. We ignore the effects of binaries, stellar evolution, and the galactic tidal field. The simulations use N = 10^5 stars and follow the evolution of the clusters until core collapse. We find that the departure from energy equipartition in the core follows approximately the theoretical predictions of Spitzer (1969) and Lightman & Fall (1978), and we suggest a more exact condition that is based on our results. We find good agreement with previous results obtained by other methods regarding several important features of the evolution, including the pre-collapse distribution of heavier stars, the time scale on which equipartition is approached, and the extent to which core collapse is accelerated by a small subpopulation of heavier stars. We briefly discuss the possible implications of our results for the dynamical evolution of primordial black holes and neutron stars in globular clusters., Comment: 31 pages, including 13 figures, to appear in ApJ

We present a new theoretical analysis of the PSR B1620-26 triple system in the globular cluster M4, based on the latest radio pulsar timing data, which now include measurements of five time derivatives of the pulse frequency. These data allow us to determine the mass and orbital parameters of the second companion completely (up to the usual unknown orbital inclination angle $i_2$). The current best-fit parameters correspond to a second companion of planetary mass, $m_2 \sin i_2 \simeq 7\times10^{-3} M_\odot$, in an orbit of eccentricity $e_2\simeq 0.45$ and semimajor axis $a_2\simeq 60 $AU. Using numerical scattering experiments, we study a possible formation scenario for the triple system, which involves a dynamical exchange interaction between the binary pulsar and a primordial star-planet system. The current orbital parameters of the triple are consistent with such a dynamical origin, and suggest that the separation of the parent star-planet system was very large, $\go 50 $AU. We also examine the possible origin of the anomalously high eccentricity of the inner binary pulsar. While this eccentricity could have been induced during the same dynamical interaction that created the triple, we find that it could equally well arise from long-term secular perturbation effects in the triple, combining the general relativistic precession of the inner orbit with the Newtonian gravitational perturbation of the planet. The detection of a planet in this system may be taken as evidence that large numbers of extrasolar planetary systems, not unlike those discovered recently in the solar neighborhood, also exist in old star clusters., Comment: 34 pages, including 11 figures, to appear in ApJ, Jan 1, 2000. Revised version, including some significant additions

We present a general method for determining the masses and orbital parameters of binary millisecond pulsars with long orbital periods (P_orb >> 1 yr), using timing data in the form of pulse frequency derivatives. We apply our method to analyze the properties of the second companion in the PSR B1620-26 triple system. We use the latest timing data for this system to constrain the mass and orbital parameters of the second companion. We find that all possible solutions have a mass m_2 in the range 2.4 10^-4 M_sun <= m_2 sin i_2 <= 1.2 10^-2 M_sun, i.e., almost certainly excluding a second companion of stellar mass and suggesting instead that the system contains a planet or a brown dwarf. Using Monte-Carlo realizations of the triple configuration in three dimensions we find the most probable value of m_2 to be 0.010(5) M_sun, corresponding to a distance of 38(6) AU from the center of mass of the inner binary (the errors indicate 80% confidence intervals). We also apply our method to analyze the planetary system around PSR B1257+12, where a distant, giant planet may be present in addition to the three well-established Earth-mass planets. We find that the simplest interpretation of the frequency derivatives implies the presence of a fourth planet with a mass of ~100 M_earth in a circular orbit of radius ~40 AU., Comment: 30 pages, Latex, 10 Postscript figures, uses aaspp4.sty. ApJ submitted. Also available at http://ensor.mit.edu/~rasio

Previous timing data for PSR~B1620$-$26 were consistent with a second companion mass $m_2$ anywhere in the range $\sim10^{-3}-1\,M_\odot$, i.e., from a Jupiter-type planet to a star. We present the latest timing parameters for the system, including a significant change in the projected semi-major axis of the inner binary, a marginal detection of the fourth time derivative of the pulse frequency, and the pulsar proper motion (which is in agreement with published values for the proper motion of M4), and use them to further constrain the mass $m_2$ and the orbital parameters. Using the observed value of $\stackrel{\ldots.}{f}$, we obtain a one-parameter family of solutions, all with $m_2 \lo 10^{-2}\,M_\odot$, i.e., excluding stellar masses. Varying $\stackrel{\ldots.}{f}$ within its formal $1\sigma$ error bar does not affect the mass range significantly. However, if we vary $\stackrel{\ldots.}{f}$ within a $4\sigma$ error bar, we find that stellar-mass solutions are still possible. We also calculate the predicted rate of change of the projected semi-major axis of the inner binary and show that it agrees with the measured value., Comment: uses PASP LaTeX macros, 6 pages, to appear in IAU Colloquium 160, Pulsars: Problems and Progress, also available from http://ensor.mit.edu/~rasio/

The 4.2-day orbit of the newly discovered planet around 51~Pegasi is formally unstable to tidal dissipation. However, the orbital decay time in this system is longer than the main-sequence lifetime of the central star. Given our best current understanding of tidal interactions, a planet of Jupiter's mass around a solar-like star could have dynamically survived in an orbit with a period as short as $\sim10\,$hr. Since radial velocities increase with decreasing period, we would expect to find those planets close to the tidal limit first and, unless this is a very unusual system, we would expect to find many more. We also consider the tidal stability of planets around more evolved stars and we re-examine in particular the question of whether the Earth can dynamically survive the red-giant phase in the evolution of the Sun., Comment: AAS LaTeX macros v.4, 14 pages, 2 postscript figures, also available from http://ensor.mit.edu/~rasio/, to appear in ApJ

We report the results of new SPH calculations of parabolic collisions between two main-sequence stars in a globular cluster. Such collisions are directly relevant to the formation of blue stragglers. In particular, we consider parent stars of mass $M/M_{TO}=0.2,0.5,0.75$, and $1$, where $M_{TO}$ is the cluster turnoff mass (typically about $0.8\,M_\odot$). Low-mass stars (with $M=0.2 M_{TO}$ or $0.5 M_{TO}$) are assumed to be fully convective and are therefore modeled as $n=1.5$ polytropes. Stars at the turnoff (with $M= M_{TO}$) are assumed to be mostly radiative and are modeled as $n=3$ polytropes. Intermediate-mass stars (with $M=0.75 M_{TO}$) are modeled as composite polytropes consisting of a radiative core with polytropic index $n=3$ and a convective envelope with $n=1.5$. We focus our study on the question of hydrodynamic mixing of helium and hydrogen, which plays a crucial role in determining the observable characteristics of blue stragglers. In all cases we find that there is negligible hydrodynamic mixing of helium into the outer envelope of the merger remnant. The amount of hydrogen mixed into the core of the merger depends strongly on the entropy profiles of the two colliding stars. In contrast to what has been done in previous studies, blue stragglers formed by direct stellar collisions should not necessarily be assumed to have initially homogeneous composition profiles., Comment: uses AAS latex macros v3, 66 pages with tables and figures, submitted to ApJ, complete postscript also available at http://ensor.mit.edu/~rasio/papers

In addition to their possible relevance to gamma-ray bursts, coalescing binary neutron stars have long been recognized as important sources of gravitational radiation that should become detectable with the new generation of laser interferometers such as LIGO. Hydrodynamics plays an essential role near the end of the coalescence when the two stars finally merge together into a single object. The shape of the corresponding burst of gravitational waves provides a direct probe into the interior structure of a neutron star and the nuclear equation of state. The interpretation of the gravitational waveform data will require detailed theoretical models of the complicated three-dimensional hydrodynamic processes involved. Here we review the results of our recent work on this problem, using both approximate quasi-analytic methods and large-scale numerical hydrodynamics calculations on supercomputers. We also discuss briefly the coalescence of white-dwarf binaries, which are also associated with a variety of interesting astrophysical phenomena., Comment: to appear in Compact Stars in Binaries, proceedings of IAU Symposium 165, ed. E. P. J. van den Heuvel, 12 pages, uuencoded compressed postscript, IAS preprint AST 94/55

The hierarchical triple system containing the millisecond pulsar PSR B1620-26 in M4 is the first triple star system ever detected in a globular cluster. Such systems should form in globular clusters as a result of dynamical interactions between binaries. We propose that the triple system containing PSR B1620-26 formed through an exchange interaction between a wide primordial binary and a {\it pre-existing\/} binary millisecond pulsar. This scenario would have the advantage of reconciling the $\sim10^9\,$yr timing age of the pulsar with the much shorter lifetime of the triple system in the core of M4., Comment: 13 pages, uuencoded compressed postscript with figures, IASSNS-AST 94/47

We have searched deep optical images of the globular cluster M4 for a possible optical counterpart for the third body in the PSR B1620-26 system. We identify as a possible optical counterpart a star located within $0.3''$ of the nominal pulsar position with $V\approx 20$. This magnitude is consistent with a main sequence cluster member with $M\approx 0.45M_\odot$. However, this star may be blended with fainter stars, and the probability of a chance superposition is non-negligible., Comment: 8 pages in uuencoded compressed postscript, to appear in ApJ Letters, Yale Preprint 25-07-94

We present a new numerical study of the equilibrium and stability properties of close binary systems using the smoothed-particle hydrodynamics (SPH) technique. We adopt a simple polytropic equation of state $p=K\rho^\gam$ with $\gam=5/3$ and $K=\,$constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio $q\lo0.4$ we find that the (synchronized) Roche limit configuration is secularly unstable. Degenerate binary configurations remain hydrodynamically stable all the way to the Roche limit for all mass ratios $q\ne1$. Dynamically unstable mass transfer can also lead to the rapid coalescence of a binary system, but the details of the hydrodynamic evolution are quite different. We discuss the implications of our results for double white-dwarf and W Ursae Majoris systems., Comment: 26 pages, LaTeX with AAS macros, 14 figures available upon request from best@sns.ias.edu, IAS-AST-94-24

A simple analytic calculation is presented for the secular evolution of the eccentricities in a hierarchical triple system such as the one containing the millisecond pulsar PSR 1620-26 in M4. If the second companion of PSR 1620-26 is of stellar mass ($m_2\ga0.1\,M_\odot$), an eccentricity as large as that observed today for the inner binary ($e_1\simeq0.03$) could very well have been induced by secular perturbations. In contrast, such a large eccentricity cannot be induced by a second companion of planetary mass., Comment: to appear in Millisecond Pulsars: the Decade of Surprises, ed. A. S. Fruchter, M. Tavani, and D. Backer (ASP Conf. Series), uuencoded compressed postscript, 5 pages, IAS-AST-94-25

Recent observations by the ROSAT X-ray satellite of the Vela supernova remnant have revealed, in addition to the previously identified compact nebula, a nearly circular emitting region with a radius of about 4 degrees. The Vela pulsar is slightly off the center of this circular region, consistent with its measured proper motion of about 100km/s and an age of about 10^4 yr. The emitting region is bounded by the main supernova shock. Just outside the shock, the X-ray image reveals several well-defined V-shaped features extending radially outwards. These features are most likely wakes produced by objects moving supersonically through the outside medium. The shapes and orientations of the wakes suggest that these objects have been ejected from the center of the supernova explosion. Their present positions indicate that they have been moving with a mean velocity of a few thousand km/s. We show that pre-existing objects, such as planets in orbit around the progenitor star, could not have been accelerated to sufficiently high velocities or would have been destroyed. Instead, we propose that the observed objects are fragments ejected during the formation of the neutron star. Fragmentation during gravitational collapse is a natural consequence of both convective and rotational instabilities., Comment: 8 pages, TEX, IAS preprint

The unusually large eccentricity ($e_1=0.025$) of the low-mass binary millisecond pulsar PSR B1620-26 can be explained naturally as arising from the secular perturbation of a second, more distant companion. Such a triple configuration has been proposed recently as the most likely cause of the anomalous second period derivative of the pulsar. The current timing data are consistent with a second companion mass $m_2$ as low as $\sim10^{-3}\,M_\odot$, i.e., comparable to that of Jupiter. However, {\em if\/} the eccentricity is indeed produced by secular perturbations, then the second companion must be another star, most likely of mass $m_2\lo1M_\odot$ and in a very eccentric ($e_2\go0.5$) orbit of period $P_2\sim10^2$--$10^3\,$yr. A second companion of planetary mass cannot induce the observed eccentricity. Independent of the mass of the second companion, small changes in the binary pulsar's orbit should become detectable with just a few additional years of timing data. This detection would provide direct confirmation of the triple nature of the system, and an accurate measurement of the effects would place important new constraints on the orbital parameters., Comment: 11 pages, uuencoded compressed postscript includes figures, IAS-AST-94-2098

We have performed a series of three-dimensional hydrodynamic calculations of binary coalescence using the smoothed particle hydrodynamics (SPH) method. The initial conditions are exact polytropic equilibrium configurations with $\gam > 5/3$, on the verge of dynamical instability. We calculate the emission of gravitational radiation in the quadrupole approximation. The fully nonlinear development of the instability is followed until a new equilibrium configuration is reached. We find that the properties of this final configuration depend sensitively on both the compressibility and mass ratio. An {\em axisymmetric} merged configuration is always produced when $\gam\lo2.3$. As a consequence, the emission of gravitational radiation shuts off abruptly right after the onset of dynamical instability. In contrast, {\em triaxial\/} merged configurations are obtained when $\gam\go2.3$, and the system continues to emit gravitational waves after the final coalescence. Systems with mass ratios $q\ne1$ typically become dynamically unstable before the onset of mass transfer. Stable mass transfer from one neutron star to another in a close binary is therefore probably ruled out. The maximum amplitude $h_{max}$ and peak luminosity $L_{max}$ of the gravitational waves emitted during the final coalescence are nearly independent of $\gam$, but depend very sensitively on the mass ratio $q$., Comment: 27 pages, uuencoded compressed postscript, 16 figures upon request from BEST@sns.ias.edu, IAS-AST-94-4

The formation of quasar black holes during the hydrodynamic collapse of protogalactic gas clouds is discussed. The dissipational collapse and long- term dynamical evolution of these systems is analysed using three-dimensional numerical simulations. The calculations focus on the final collapse stages of the inner baryonic component and therefore ignore the presence of dark matter. Two types of initial conditions are considered: uniformly rotating spherical clouds, and irrotational ellipsoidal clouds. In both cases the clouds are initially cold, homogeneous, and not far from rotational support ($T/|W| \approx 0.1$). Although the details of the dynamical evolution depend sensitively on the initial conditions, the qualitative features of the final configurations do not. Most of the gas is found to fragment into small dense clumps, that eventually make up a spheroidal component resembling a galactic bulge. About 5\% of the initial mass remains in the form of a smooth disk of gas supported by rotation in the gravitational potential well of the outer spheroid. If a central seed black hole of mass $\gsim10^6M_\odot$ forms, it can grow by steady accretion from the disk and reach a typical quasar black hole mass $\sim10^8M_\odot$ in less than $5\times 10^8\,$yr. In the absence of a sufficiently massive seed, dynamical instabilities in a strongly self- gravitating inner region of the disk will inhibit steady accretion of gas and may prevent the immediate formation of a quasar., Comment: 20 pages, uuencoded compressed postscript, 4 figures available upon request from BEST@sns.ias.edu, IAS-AST-94-3

We present a new analytic study of the equilibrium and stability properties of close binary systems containing polytropic components. Our method is based on the use of ellipsoidal trial functions in an energy variational principle. We consider both synchronized and nonsynchronized systems, constructing the compressible generalizations of the classical Darwin and Darwin-Riemann configurations. Our method can be applied to a wide variety of binary models where the stellar masses, radii, spins, entropies, and polytropic indices are all allowed to vary over wide ranges and independently for each component. We find that both secular and dynamical instabilities can develop before a Roche limit or contact is reached along a sequence of models with decreasing binary separation. High incompressibility always makes a given binary system more susceptible to these instabilities, but the dependence on the mass ratio is more complicated. As simple applications, we construct models of double degenerate systems and of low-mass main-sequence-star binaries. We also discuss the orbital evolution of close binary systems under the combined influence of fluid viscosity and secular angular momentum losses from processes like gravitational radiation. We show that the existence of global fluid instabilities can have a profound effect on the terminal evolution of coalescing binaries. The validity of our analytic solutions is examined by means of detailed comparisons with the results of recent numerical fluid calculations in three dimensions., Comment: 37 pages, plain TeX, postscript figures and tables in separate uufile, IAS AST 93/44

We study the importance of hydrodynamic effects on the evolution of coalescing binary neutron stars. Using an approximate energy functional constructed from equilibrium solutions for polytropic binary configurations, we incorporate hydrodynamic effects into the calculation of the orbital decay driven by gravitational wave emission. In particular, we follow the transition between the quasi-static, secular decay of the orbit at large separation and the rapid dynamical evolution of configurations approaching contact. We show that a purely Newtonian hydrodynamic instability can significantly accelerate the coalescence at small separation. Such an instability occurs in all close binary configurations containing sufficiently incompressible stars. Calculations are performed for various neutron star masses, radii, spins, and effective polytropic indices. Typically, we find that the radial infall velocity just prior to contact is about 10\% of the tangential orbital velocity. Post-Newtonian effects can move the stability limit to a larger binary separation, and may induce an even larger radial velocity. We also consider the possibility of mass transfer from one neutron star to the other. We show that stable mass transfer is impossible except when the mass of one of the components is very small (less than about 0.4 solar mass) and the viscosity is high enough to maintain corotation., Comment: Latex with aaspp style (AAS latex macros v. 3.0). Figures and tables: latex tables and postscript figures in a separate file. IASSNS number AST93/22

A model based on the binary configuration of the PSR1829-10 pulsar (Bailes et al., 1991) is used to show that the formation of a binary pulsar with a planet-size companion, large original separation, and small eccentricity could result from the rapid evaporation of a much more massive binary companion by the pulsar's radiation. Such an evaporation process is known to be taking place in at least two other binary pulsars: PSR1957 + 20 (Fruchter et al., 1990; Ryba and Taylor, 1991) and PSR1744 - 24A (Lyne et al., 1990). It is shown here that, about one million years ago, the companion mass and binary separation could have been comparable to those currently observed in the eclipsing binary pulsar PSR1957 + 20.