The mass distribution of neutron stars and stellar-mass black holes provides vital clues into the nature of stellar core collapse and the physical engine responsible for supernova explosions. Using recent advances in our understanding of supernova engines, we derive mass distributions of stellar compact remnants. We provide analytical prescriptions for compact object masses for major population synthesis codes. In an accompanying paper, Belczynski et al., we demonstrate that these qualitatively new results for compact objects can explain the observed gap in the remnant mass distribution between ∼ 2 − 5 M ⊙ and that they place strong constraints on the nature of the supernova engine. Here, we show that advanced gravitational radiation detectors (like LIGO/VIRGO or the Einstein Telescope) will be able to further test the supernova explosion engine models once double black hole inspirals are detected.

We present an up-to-date, comprehensive summary of the rates for all types of compact binary coalescence sources detectable by the initial and advanced versions of the ground-based gravitational-wave detectors LIGO and Virgo. Astrophysical estimates for compact-binary coalescence rates depend on a number of assumptions and unknown model parameters and are still uncertain. The most confident among these estimates are the rate predictions for coalescing binary neutron stars which are based on extrapolations from observed binary pulsars in our galaxy. These yield a likely coalescence rate of 100 Myr−1 per Milky Way Equivalent Galaxy (MWEG), although the rate could plausibly range from 1 Myr−1 MWEG−1 to 1000 Myr−1 MWEG−1 (Kalogera et al 2004 Astrophys. J. 601 L179; Kalogera et al 2004 Astrophys. J. 614 L137 (erratum)). We convert coalescence rates into detection rates based on data from the LIGO S5 and Virgo VSR2 science runs and projected sensitivities for our advanced detectors. Using the detector sensitivities derived from these data, we find a likely detection rate of 0.02 per year for Initial LIGO–Virgo interferometers, with a plausible range between 2 × 10−4 and 0.2 per year. The likely binary neutron–star detection rate for the Advanced LIGO–Virgo network increases to 40 events per year, with a range between 0.4 and 400 per year.

A new generation of ground-based interferometric detectors for gravitational waves is currently under construction or has entered the commissioning phase (LIGO, VIRGO, GEO600, TAMA). The purpose of these detectors is to observe gravitational waves from astrophysical sources and help improve our understanding of the source origin and physical properties. In this paper we study the most promising candidate sources for these detectors: inspiraling double compact objects. We use population synthesis methods to calculate the properties and coalescence rates of compact object binaries: double neutron stars, black hole-neutron star systems and double black holes. We also examine the formation channels available to double compact object binaries. We explicitly account for the evolution of low-mass helium stars and investigate the possibility of common-envelope evolution involving helium stars as well as two evolved stars. As a result we identify a significant number of new formation channels for double neutron stars, in particular, leading to populations with very distinct properties. We discuss the theoretical and observational implications of such populations, but we also note the need for hydrodynamical calculations to settle the question of whether such common-envelope evolution is possible. We also present and discuss the physical properties of compact object binaries and identify a number of robust, qualitative features as well as their origin. Using the calculated coalescence rates we compare our results to earlier studies and derive expected detection rates for the Laser Interferometer Gravitational-wave Observatory (LIGO). We find that our most optimistic estimate for the first LIGO detectors reach a couple of events per year and our most pessimistic estimate for advanced LIGO detectors exceed ≃ 10 events per year.

The last decade of observational and theoretical developments in stellar and binary evolution provides an opportunity to incorporate major improvements to the predictions from populations synthesis models. We compute the Galactic merger rates for NS-NS, BH-NS, and BH-BH mergers with the StarTrack code. The most important revisions include: updated wind mass loss rates (allowing for stellar mass black holes up to 80 M ⊙ ), a realistic treatment of the common envelope phase (a process that can affect merger rates by 2-3 orders of magnitude), and a qualitatively new neutron star/black hole mass distribution (consistent with the observed "mass gap"). Our findings include: (i) The binding energy of the envelope plays a pivotal role in determining whether a binary merges within a Hubble time. (ii) Our description of natal kicks from supernovae plays an important role, especially for the formation of BH-BH systems. (iii) The masses of BH-BH systems can be substantially increased in the case of low metallicities or weak winds. (iv) Certain combinations of parameters underpredict the Galactic NS-NS merger rate, and can be ruled out. (v) Models incorporating delayed supernovae do not agree with the observed NS/BH "mass gap", in accordance with our previous work. This is the first in a series of three papers. The second paper will study the merger rates of double compact objects as a function of redshift, star formation rate, and metallicity. In the third paper we will present the detection rates for gravitational wave observatories, using up-to-date signal waveforms and sensitivity curves.

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.