Interacting Binaries with Eccentric Orbits. III. Orbital Evolution due to Direct Impact and Self-Accretion

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 Ap

Mass Transfer in Binary Stars using SPH. I. Numerical Method

Close interactions and mass transfer in binary stars can lead to the formation of many different exotic stellar populations, but detailed modeling of mass transfer is a computationally challenging problem. Here, we present an alternate Smoothed Particle Hydrodynamics approach to the modeling of mass transfer in binary systems that allows a better resolution of the flow of matter between main-sequence stars. Our approach consists of modeling only the outermost layers of the stars using appropriate boundary conditions and ghost particles. We arbitrarily set the radius of the boundary and find that our boundary treatment behaves physically and conserves energy well. In particular, when used with our binary relaxation procedure, our treatment of boundary conditions is also shown to evolve circular binaries properly for many orbits. The results of our first simulation of mass transfer are also discussed and used to assess the strengths and limitations of our method. We conclude that it is well suited for the modeling of interacting binary stars. The method presented here represents a convenient alternative to previous hydrodynamical techniques aimed at modeling mass transfer in binary systems since it can be used to model both the donor and the accretor while maintaining the density profiles taken from realistic stellar models.Comment: 10 pages, 10 figures. Accepted for publication in The Astrophysical Journa

Are Supernova Kicks Responsible for X-ray Binary Ejection from Young Clusters?

Recent Chandra observations of interacting and starburst galaxies have led us to investigate the apparent correlation between the positions of young star clusters and Chandra point sources. Assumed to be X-ray binaries (XRBs), these point sources do not seem to coincide with the massive (~1e5 Msun), young (1-50 Myr) stellar clusters that can easily form systems capable of such emission. We use a sophisticated binary evolution and population synthesis code (StarTrack) and a simplified cluster model to track both the X-ray luminosity and position of XRBs as a function of time. These binaries are born within the cluster potential with self-consistent positions and velocities and we show that a large fraction (~70%) can be ejected from the parent due to supernova explosions and associated systemic velocities. For brighter sources and cluster masses below ~1e6 Msun, we find that the average number of bright XRBs per cluster remains near or below unity, consistent with current observations.Comment: 5 pages, 1 figure. Accepted for publication in Astrophysical Journal Letter

Angular Momentum Changes Due to Direct Impact Accretion in a Simplified Binary System

We model a circular mass-transferring binary system to calculate the exchange of angular momentum between stellar spins and the orbit due to direct impact of the mass transfer stream onto the surface of the accretor. We simulate mass transfer by calculating the ballistic motion of a point mass ejected from the $L_1$ point of the donor star, conserving the total linear and angular momentum of the system, and treating the stars as uniform density spheres with main sequence radii determined by their masses. We show that, contrary to previous assumptions in the literature, direct impact does not always act as a sink of orbital angular momentum and may in fact increase it by facilitating the transfer of angular momentum from the spin of the donor to the orbit. Here, we show an example of the exchange of angular momentum, as well as a measure of the orbital angular momentum changes for a variety of binary star systems with main sequence components.Comment: 2 pages, 2 figures, Conference Proceedings for the International Conference on Binaries, Mykonos, Greece. Updated Version of Fig. 1b, correcting a scaling error. Results remain unchanged, but the numerical scaling factors have been decrease

Probing Electron-Capture Supernovae: X-Ray Binaries in Starbursts

Presenting population models of high-mass X-ray binaries (HMXBs) formed after bursts of star formation, we investigate the effect of electron-capture supernovae (ECS) of massive ONeMg white dwarfs and the hypothesis that ECS events are associated with typically low supernova kicks imparted to the nascent neutron stars. We identify an interesting ECS bump in the time evolution of HMXB numbers; this bump is caused by significantly increased production of wind-fed HMXBs 20-60 Myr post starburst. The amplitude and age extent of the ECS bump depend on the strength of ECS kicks and the mass range of ECS progenitors. We also find that ECS-HMXBs form through a specific evolutionary channel that is expected to lead to binaries with Be donors in wide orbits. These characteristics, along with their sensitivity to ECS properties, provide us with an intriguing opportunity to probe ECS physics and progenitors through studies of starbursts of different ages. Specifically, the case of the Small Magellanic Cloud, with a significant observed population of Be HMXBs and starburst activity 30-60 Myr ago, arises as a promising laboratory for understanding the role of electron-capture supernovae in neutron star formation.Comment: 5 pages, 3 figures, Published by ApJ in 07/0

The Red Rectangle: Its Shaping Mechanism and its Source of Ultraviolet Photons

The proto-planetary Red Rectangle nebula is powered by HD 44179, a spectroscopic binary (P = 318 d), in which a luminous post-AGB component is the primary source of both luminosity and current mass loss. Here, we present the results of a seven-year, eight-orbit spectroscopic monitoring program of HD 44179, designed to uncover new information about the source of the Lyman/far-ultraviolet continuum in the system as well as the driving mechanism for the bipolar outflow producing the current nebula. Our observations of the H-alpha line profile around the orbital phase of superior conjunction reveal the secondary component to be the origin of the fast (max. v~560$km s$^{-1}$) bipolar outflow in the Red Rectangle. The variation of total H-alpha flux from the central H II region with orbital phase also identifies the secondary or its surroundings as the source of the far-ultraviolet ionizing radiation in the system. The estimated mass of the secondary (~0.94 M$\sun$) and the speed of the outflow suggest that this component is a main sequence star and not a white dwarf, as previously suggested. We identify the source of the Lyman/far-ultraviolet continuum in the system as the hot, inner region (T$_{max} \ge 17,000$K) of an accretion disk surrounding the secondary, fed by Roche lobe overflow from the post-AGB primary at a rate of about$2 - 5\times10^{-5}$M$\sun$yr$^{-1}$. The total luminosity of the accretion disk around the secondary is estimated to be at least 300 L$\sun$, about 5% of the luminosity of the entire system. (abridged)Comment: Accepted for publication in Ap

Stellar Collisions and the Interior Structure of Blue Stragglers

Collisions of main sequence stars occur frequently in dense star clusters. In open and globular clusters, these collisions produce merger remnants that may be observed as blue stragglers. Detailed theoretical models of this process require lengthy hydrodynamic computations in three dimensions. However, a less computationally expensive approach, which we present here, is to approximate the merger process (including shock heating, hydrodynamic mixing, mass ejection, and angular momentum transfer) with simple algorithms based on conservation laws and a basic qualitative understanding of the hydrodynamics. These algorithms have been fine tuned through comparisons with the results of our previous hydrodynamic simulations. We find that the thermodynamic and chemical composition profiles of our simple models agree very well with those from recent SPH (smoothed particle hydrodynamics) calculations of stellar collisions, and the subsequent stellar evolution of our simple models also matches closely that of the more accurate hydrodynamic models. Our algorithms have been implemented in an easy to use software package, which we are making publicly available (see http://vassun.vassar.edu/~lombardi/mmas/). This software could be used in combination with realistic dynamical simulations of star clusters that must take into account stellar collisions.Comment: This revised version has 37 pages, 13 figures, 4 tables; submitted to ApJ; for associated software package, see http://vassun.vassar.edu/~lombardi/mmas/ This revised version presents additional comparisons with SPH results and slightly improved merger recipe