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 observations with 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 extragalactic 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.

No abstract

This work aims to present our current best physical understanding of common-envelope evolution (CEE). We highlight areas of consensus and disagree- ment, and stress ideas which should point the way forward for progress in this important but long-standing and largely unconquered problem. Unusually for CEE-related work, we mostly try to avoid relying on results from population synthesis or observations, in order to avoid potentially being misled by previous misunderstandings. As far as possible we debate all the relevant issues starting from physics alone, all the way from the evolution of the binary system immediately before CEE begins to the processes which might occur just after the ejection of the envelope. In particular, we include extensive discussion about the energy sources and sinks operating in CEE, and hence examine the foundations of the standard energy formalism. Special attention is also given to comparing the results of hydrodynamic simulations from different groups and to discussing the potential effect of initial conditions on the differences in the outcomes. We compare current numerical techniques for the problem of CEE and also whether more appropriate tools could and should be produced (including new formulations of computational hydrodynamics, and attempts to include 3D processes within 1D codes). Finally we explore new ways to link CEE with observations. We compare previous simulations of CEE to the recent outburst from V1309 Sco, and discuss to what extent post-common-envelope binaries and nebulae can provide information, e.g. from binary eccentricities, which is not currently being fully exploited.

Double black hole binaries are among the most important sources of gravitational radiation for groundbased detectors such as LIGO or VIRGO. Even if formed with lower efficiency than double neutron star binaries, they could dominate the predicted detection rates, since black holes are more massive than neutron stars and therefore could be detected at greater distances. Here we discuss an evolutionary process that can very significantly limit the formation of close double black hole binaries: the vast majority of their potential progenitors undergo a common envelope (CE) phase while the donor, one of the massive binary components, is evolving through the Hertzsprung gap. Our latest theoretical understanding of the CE process suggests that this will probably lead to a merger, inhibiting double black hole formation. Barring uncertainties in the physics of CE evolution, we use population synthesis calculations, and find that the corresponding reduction in the merger rate of double black holes formed in galactic fields is so great (by ∼ 500) that their contribution to inspiral detection rates for ground-based detectors could become relatively small (∼ 1 in 10) compared to double neutron star binaries. A similar process also reduces the merger rates for double neutron stars, by factor of ∼ 5, eliminating most of the previously predicted ultracompact NS-NS systems. Our predicted detection rates for Advanced LIGO are now much lower for double black holes (∼ 2 yr −1 ), but are still quite high for double neutron stars (∼ 20 yr −1 ). If double black holes were found to be dominant in the detected inspiral signals, this could indicate that they mainly originate from dense star clusters (not included here) or that our theoretical understanding of the CE phase requires significant revision.

The hydrodynamic evolution of the common envelope phase of a low mass binary composed of a 1.05M red giant and a 0.6M companion has been followed for five orbits of the system using a high resolution method in three spatial dimensions. During the rapid inspiral phase, the interaction of the companion with the red giant's extended atmosphere causes about 25% of the common envelope to be ejected from the system, with mass continuing to be lost at the end of the simulation at a rate ∼ 2M yr −1 . In the process the resulting loss of angular momentum and energy reduces the orbital separation by a factor of seven. After this inspiral phase the eccentricity of the orbit rapidly decreases with time. The gravitational drag dominates hydrodynamic drag at all times in the evolution, and the commonly-used Bondi-Hoyle-Lyttleton prescription for estimating the accretion rate onto the companion significantly overestimates the true rate. On scales comparable to the orbital separation, the gas flow in the orbital plane in the vicinity of the two cores is subsonic with the gas nearly corotating with the red giant core and circulating about the red giant companion. On larger scales, 90% of the outflow is contained within 30 degrees of the orbital plane, and the spiral shocks in this material leave an imprint on the density and velocity structure. Of the energy released by the inspiral of the cores, only about 25% goes toward ejection of the envelope.

No abstract

Type Ia supernovae are thought to be caused by thermonuclear explosions of a carbon-oxygen white dwarf in close binary systems. In the single-degenerate scenario (SDS), the companion star is non-degenerate and can be significantly affected by the explosion. We explore this interaction by means of multi-dimensional adaptive mesh refinement simulations using the FLASH code. We consider several different companion types, including main-sequence-like stars (MS), red giants (RG), and helium stars (He). In addition, we include the symmetry-breaking effects of orbital motion, rotation of the non-degenerate star, and Roche-lobe overflow. A detailed study of a sub-grid model for Type Ia supernovae is also presented. We find that the dependence of the unbound stellar mass and kick velocity on the initial binary separation can be fitted by power-law relations. By using the tracer particles in FLASH, the process leading to the unbinding of matter is dominated by ablation, which has usually been neglected in past analytical studies. The level of Ni/Fe contamination of the companion that results from the passage of the supernova ejecta is found to be a ∼ 10 −5 M ⊙ for the MS star, ∼ 10 −4 M ⊙ for the He star, and ∼ 10 −8 M ⊙ for the RG. The spinning MS companion star loses about half of its initial angular momentum during the impact, causing the rotational velocity to drop to a quarter of the original rotational velocity, suggesting that the Tycho G star is a promising progenitor candidate in the SDS.

Pair-instability and pulsational pair-instability supernovae (PPISN) have not been unambiguously observed so far. They are, however, promising candidates for the progenitors of the heaviest binary black hole (BBH) mergers detected. If these BBHs are the product of binary evolution, then PPISNe could occur in very close binaries. Motivated by this, we discuss the implications of a PPISN happening with a close binary companion, and what impact these explosions have on the formation of merging BBHs through binary evolution. For this, we have computed a set of models of metal-poor (Z /10) helium stars using the MESA software instrument. For PPISN progenitors with pre-explosion masses > 50M we find that, after a pulse, heat deposited throughout the layers of the star that remain bound cause it to expand to more than 100R for periods of 10 2 − 10 4 yrs depending on the mass of the progenitor. This results in long-lived phases of Roche-lobe overflow or even common-envelope events if there is a close binary companion, leading to additional electromagnetic transients associated to PPISN eruptions. If we ignore the effect of these interactions, we find that mass loss from PPISNe reduces the final black hole spin by ∼ 30%, induces eccentricities that can be detected by the LISA observatory, and can produce a double-peaked distribution of measured chirp masses in BBH mergers observed by ground-based detectors.