We calculate the long-term evolution of angular momentum in double white dwarf binaries undergoing direct impact accretion over a broad range of parameter space. We allow the rotation rate of both components to vary, and account for the exchange of angular momentum between the spins of the white dwarfs and the orbit, while conserving the total angular momentum. We include gravitational, tidal, and mass transfer effects in the orbital evolution, and allow the Roche radius of the donor star to vary with both the stellar mass and the rotation rate. We examine the long-term stability of these systems, focusing in particular on those systems that may be progenitors of AM CVn or Type Ia Supernovae. We find that our analysis yields an increase in the predicted number of stable systems compared to that in previous studies. Additionally, we find that by properly accounting for the effects of asynchronism between the donor and the orbit on the Roche-lobe size, we eliminate oscillations in the orbital parameters which are found in previous studies. Removing these oscillations can reduce the peak mass transfer rate in some systems, keeping them from entering an unstable mass transfer phase.

We examine the exchange of angular momentum between the component spins and the orbit in semi-detached double white dwarf binaries undergoing mass transfer through direct impact of the transfer stream. We approximate the stream as a series of discrete massive particles ejected in the ballistic limit at the inner Lagrangian point of the donor toward the accretor. This work improves upon similar earlier studies in a number of ways. First, we self-consistently calculate the total angular momentum of the orbit at all times. This includes changes in the orbital angular momentum during the ballistic trajectory of the ejected mass, as well as changes during the ejection/accretion due to the radial component of the particle's velocity. Second, we calculate the particle's ballistic trajectory for each system, which allows us to determine the precise position and velocity of the particle upon accretion. We can then include specific information about the radius of the accretor as well as the angle of impact. Finally, we ensure that the total angular momentum is conserved, which requires the donor star spin to vary self-consistently. We calculate the angular momentum change of the orbit and each binary component across the entire parameter space of direct impact double white dwarf binaries. We find a significant decrease in the amount of angular momentum removed from the orbit during mass transfer, as well as cases where this process increases the angular momentum of the orbit at the expense of the spin angular momentum of the donor. We conclude that, unlike earlier claims in the literature, mass transfer through direct impact need not destabilize the binary and that the quantity and sign of the orbital angular momentum transfer depends on the binary properties. This stabilization may significantly impact the population synthesis calculations of double white dwarf progenitors.

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.

We investigate the characteristics of young (< 20 Myr) and bright (Lx > 1e36 erg/s) High-Mass X-ray Binaries (HMXBs) and find the population to be strongly metallicity-dependent. We separate the model populations among two distinct formation pathways: (1) systems undergoing active Roche Lobe Overflow (RLO), and (2) wind accretion systems with donors in the (super)giant (SG) stage, which we find to dominate the HMXB population. We find metallicity to primarily affect the number of systems which move through each formation pathway, rather than the observable parameters of systems which move through each individual pathway. We discuss the most important model parameters affecting the HMXB population at both low and high metallicities. Using these results, we show that (1) the population of ultra-luminous X-Ray sources can be consistently described by very bright HMXBs which undergo stable Roche Lobe overflow with mild super-Eddington accretion and (2) the HMXB population of the bright starburst galaxy NGC~1569 is likely dominated by one extremely metal-poor starburst cluster.

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.

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.

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.

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.

We investigate the existence and properties of equipotential surfaces and Lagrangian points in non-synchronous, eccentric binary star and planetary systems under the assumption of quasi-static equilibrium. We adopt a binary potential that accounts for non-synchronous rotation and eccentric orbits, and calculate the positions of the Lagrangian points as functions of the mass ratio, the degree of asynchronism, the orbital eccentricity, and the position of the stars or planets in their relative orbit. We find that the geometry of the equipotential surfaces may facilitate non-conservative mass transfer in non-synchronous, eccentric binary star and planetary systems, especially if the component stars or planets are rotating super-synchronously at the periastron of their relative orbit. We also calculate the volume-equivalent radius of the Roche lobe as a function of the four parameters mentioned above. Contrary to common practice, we find that replacing the radius of a circular orbit in the fitting formula of Eggleton (1983) with the instantaneous distance between the components of eccentric binary or planetary systems does not always lead to a good approximation to the volume-equivalent radius of the Roche-lobe. We therefore provide generalized analytic fitting formulae for the volume-equivalent Roche lobe radius appropriate for non-synchronous, eccentric binary star and planetary systems. These formulae are accurate to better than 1% throughout the relevant 2-dimensional parameter space that covers a dynamic range of 16 and 6 orders of magnitude in the two dimensions.

We present a spectroscopic analysis of HST/STIS and FOS low- and intermediate-resolution spectroscopy of 55 stars (turn-off stars, horizontal branch stars and blue stragglers) in four globular clusters (47 Tucanae, M3, NGC6752, and NGC6397). Stars were analyzed with non-Local Thermodynamic Equilibrium model atmospheres, and values for their effective temperatures and gravities and some rotation rates were obtained. Using photometric fluxes, we also obtained radii, luminosities and spectroscopic masses.

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.