Until recently, black holes (BHs) could be discovered only through accretion from other stars in X-ray binaries, or in merging double compact objects. Improvements in astrometric and spectroscopic measurements have made it possible to detect BHs also in non-interacting BH binaries (nBHB) through a precise analysis of the companion's motion. In this study, using an updated version of the \startrack\ binary-star population modelling code and a detailed model of the Milky Way (MW) galaxy we calculate the expected number of detections for \gaia\ and \lamost\ surveys. We develop a formalism to convolve the binary population synthesis output with a realistic stellar density distribution, star-formation history (SFH), and chemical evolution for the MW, which produces a probability distribution function of the predicted compact-binary population over the MW. This avoids the additional statistical uncertainty which is introduced by methods which Monte Carlo sample from binary population synthesis output to produce one potential specific realisation of the MW compact-binary distribution, and our method is also comparatively fast to such Monte Carlo realisations. Specifically, we predict $\sim41$--$340$ nBHBs to be observed by \gaia, although the numbers may drop to $\sim10$--$70$ if the recent ($\lesssim100\;$ Myr) star formation is low ($\sim1\;M_\odot$/yr ). For \lamost\ we predict $\lesssim14$ detectable nBHBs, which is lower partially because its field-of-view covers just $\sim6\%$ of the Galaxy.

The mass, origin and evolutionary stage of the binary system LB-1 has been the subject of intense debate, following the claim that it hosts an $\sim$70$M_{\odot}$ black hole, in stark contrast with the expectations for stellar remnants in the Milky Way. We conducted a high-resolution, phase-resolved spectroscopic study of the near-infrared Paschen lines in this system, using the 3.5-m telescope at Calar Alto Observatory. We find that Pa$\beta$ and Pa$\gamma$ (after proper subtraction of the stellar absorption component) are well fitted with a standard double-peaked model, typical of disk emission. We measured the velocity shifts of the red and blue peaks at 28 orbital phases: the line center has an orbital motion in perfect antiphase with the stellar motion, and the radial velocity amplitude ranges from 8 to 13 km/s for different choices of lines and profile modelling. We interpret this curve as proof that the disk is tracing the orbital motion of the primary, ruling out the circumbinary disk and the hierarchical triple scenarios. The phase-averaged peak-to-peak half-separation (proxy for the projected rotational velocity of the outer disk) is $\sim$70 km s$^{-1}$, larger than the stellar orbital velocity and also inconsistent with a circumbinary disk. From those results, we infer a primary mass 4--8 times higher than the secondary mass. Moreover, we show that the ratio of the blue and red peaks (V/R intensity ratio) has a sinusoidal behaviour in phase with the secondary star, which can be interpreted as the effect of external irradiation by the secondary star on the outer disk. Finally, we briefly discuss our findings in the context of alternative scenarios recently proposed for LB-1. Definitive tests between alternative solutions will require further astrometric data from $Gaia$.

There remain significant uncertainties in the origin and evolution of black holes in binary systems, in particular regarding their birth sites and the influence of natal kicks. These are long-standing issues, but their debate has been reinvigorated in the era of gravitational wave detections and the improving precision of astrometric measurements. Using recent and archival characterisation of Galactic black hole X-ray binaries (BHXBs), we report here an apparent anticorrelation between P{orb} (system orbital periods) and scatter in z (elevation above the Galactic plane). The absence of long period sources at high z is not an obvious observational bias, and two possible explanatory scenarios are qualitatively explored: (1) a disc origin for BHXBs followed by natal kicks producing the scatter in z, with only the tightest binaries preferentially surviving strong kicks; (2) a halo origin, with P{orb} shortening through dynamical interactions in globular clusters (GCs). For the latter case, we show a correspondence in z-scatter between BHXBs and the GCs with most compact core radii of <0.1pc. However, the known absence of outbursting BHXB transients within Galactic GCs remains puzzling in this case, in contrast to the multitude of known GC neutron star XRBs. These results provide an interesting observational constraint for any black hole binary evolutionary model to satisfy.

The LIGO/Virgo collaboration has reported the detection of GW190412, a black hole-black hole (BH-BH) merger with the most unequal masses to date: m_1=24.4--34.7 Msun and m_2=7.4--10.1 Msun, corresponding to a mass ratio of q=0.21--0.41 (90 per cent probability range). Additionally, GW190412's effective spin was estimated to be X_eff=0.14--0.34, with the spin of the primary BH in the range a_spin=0.17--0.59. Based on this and prior detections, >10 per cent of BH-BH mergers have q<0.4. Major BH-BH formation channels (i.e., dynamics in dense stellar systems, classical isolated binary evolution, or chemically homogeneous evolution) tend to produce BH-BH mergers with comparable masses (typically with q>0.5). Here we test whether the classical isolated binary evolution channel can produce mergers resembling GW190412. We show that our standard binary evolution scenario, with the typical assumptions on input physics we have used in the past, produces such mergers. We provide an explicit example of an unequal mass BH-BH merger, which forms at low metallicity (Z=0.002) from two massive stars (77 and 35 Msun) and results in a BH-BH merger with m_1=27.0 Msun, m_2=9.9 Msun (q=0.37), with primary black hole spin a_spin=0.19 and an effective spin parameter X_eff=0.335. For this particular model of the input physics the overall BH-BH merger rate density in the local Universe (z=0) is: 73.5 Gpc^-3 yr^-1, while for systems with q<0.21, 0.28, 0.41, and 0.59 the rate density is: 0.01, 0.12, 6.8, and 22.2 Gpc^-3 yr^-1, respectively. The results from our standard model are consistent with the masses and spins of the black holes in GW190412, as well as with the LIGO/Virgo estimate of the fraction of unequal-mass BH-BH mergers.

The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code's BSE sector, such as the 'rapid' and the 'delayed' supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster.

We present an open-access database which includes a synthetic catalog of black holes in the Milky Way. To calculate evolution of single and binary stars we used updated population synthesis code StarTrack. We applied a new model of star formation history and chemical evolution of Galactic disk, bulge and halo synthesized from observational and theoretical data. We find that at the current moment Milky Way (disk+bulge+halo) contains about 1.2 x 10^8 single black holes with average mass of about 14 Msun and 9.3 x 10^6 BHs in binary systems with average mass of 19 Msun. We present basic statistical properties of BH populations such as distributions of single and binary BH masses, velocities, orbital parameters or numbers of BH binary systems in different evolutionary configurations. We find that the most massive BHs are formed in mergers of binary systems, such as BH-MS, BH+He, BH-BH. The metallicity of stellar population has a significant impact on the final BH mass due to the stellar winds. Therefore the most massive single BH in our simulation, 113 Msun, originates from a merger of a helium star and a black hole in a low metallicity stellar environment in Galactic halo. The most massive BH in binary system is 60 Msun and was also formed in Galactic halo. We constrain that only 0.006% of total Galactic halo mass (including dark matter) could be hidden in the form of stellar origin BHs which are not detectable by current observational surveys. Galactic binary BHs are minority (10% of all Galactic BHs) and most of them are in BH-BH systems. The current Galactic merger rates for two considered common envelope models which are: 3-81 Myr^-1 for BH-BH, 1-9 Myr^-1, for BH-NS and 14-59 Myr^-1 for NS-NS systems. Data files are available at https://bhc.syntheticuniverse.org/.

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

It is speculated that a merger of two massive astrophysical black holes in dense stellar environment may lead to the formation of a massive black hole in the pair-instability mass gap (70/80-125 Msun). Such a merger-formed black hole is expected to typically have a high spin (a=0.7), however, with a broad spectrum of spins being allowed (0<a<1). If such massive black hole acquires another black hole it may lead to another merger detectable by LIGO/Virgo in gravitational waves. We show that it is highly unlikely to form in this way and retain a 100 Msun black hole in a globular cluster if the black hole's spin is low (a<0.3). Massive merger-formed black holes with low spins acquire high recoil speeds (> 200 km/s) from gravitational-wave kick during formation that exceed typical escape speeds from globular clusters (50 km/s). However, a very low-spinning (a=0.1) and massive (100 Msun) black hole could be formed this way and retained in a galactic nuclear star cluster. Even though massive, merger-formed black holes with such low spins acquire high speeds during formation (400 km/s), they may avoid ejection since massive nuclear clusters have high escape velocities (300-500 km/s). This adds to already existing astrophysical scenarios of the formation of massive black holes with low spins. A future detection of a massive black hole in the pair-instability mass gap with low spin is therefore not a proof of the existence of primordial black holes, which are sometimes claimed to have low spins and arbitrarily high masses (not affected by pair-instability).

A 70Msun BH was discovered in Milky Way disk in a long period and almost circular detached binary system (LB-1) with a high metallicity 8Msun B star companion. Current consensus on the formation of BHs from high metallicity stars limits the black hole mass to be below 20Msun. Using simple evolutionary model, we show that the formation of a 70Msun BH in high metallicity environment is possible if stellar wind mass loss rates are reduced by factor of 5. As observations indicate, a fraction of massive stars have surface magnetic fields which may quench the wind mass-loss, independently of stellar mass and metallicity. We also computed detailed stellar evolution models and we confirm such a scenario. A non-rotating 85Msun model at Z=0.014 with decreased winds ends up as a 71Msun star prior core-collapse with a 32Msun helium core and a 28Msun CO core. Such star avoids pair-instability pulsation supernova mass loss and may form a 70Msun BH in the direct collapse. Stars that can form such BHs expand to significant size with radius of R>600Rsun, exceeding the size of LB-1 orbit. Therefore, we can explain the formation of BHs upto 70Msun at high metallicity and this result is independent from LB-1. However, if LB-1 hosts a massive BH we are unable to explain how such a binary star system could have formed without invoking some exotic scenarios.

If confirmed, a wide binary system of 70 $M_{\odot}$ black hole (BH) and an 8 $M_{\odot}$ main sequence star (LB-1) is observed to reside in the Milky Way (MW). We show that long term evolution of an 8 $M_{\odot}$ star around a BH with mass between 5-70 $M_{\odot}$ makes them visible as ultra-luminous X-ray (ULX) sources in the sky. Given the expected ULX phase lifetime ($\approx0.1$ Myr) and their lack of detection in the MW, we conclude that the frequency of an 8-20 $M_{\odot}$ star to be in binary around a stellar mass BH should be less than ($f<2\times10^{-3}$). This is in tension with Liu et al. (2019) claimed detection frequency of LB-1 like system around 8-20$M_{\odot}$ stars ($f\approx3\times10^{-2}$). Moreover, the 8 $M_{\odot}$ star is likely to end as a neutron star (NS) born with a very small kick from an electron capture supernova (ECSN), leaving behind a wide NS-BH binary. So far less than 1\% of all the detectable pulsars in the MW are mapped and there has been no detection of any pulsars in binary systems around BHs which sets an upper bound of about 100 possible pulsar-BH systems in the MW. We show if the NS is born from ECSN, a frequency upper limit of ($f=5\times10^{-4}$) for stars with masses $\approx 8-20~M_{\odot}$ in the MW to have a BH companion. The rate discrepancy will further increase as more pulsars are mapped in the MW, yet these searches would not be able to rule out the Liu et al. detection frequency if NSs are instead born in core collapse SNe with the commonly inferred high kick velocities.

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

Modeling the evolution of the elements in the Milky Way is a multidisciplinary and challenging task. In addition to simulating the 13 billion years evolution of our Galaxy, chemical evolution simulations must keep track of the elements synthesized and ejected from every astrophysical site of interest (e.g., supernova, compact binary merger). The elemental abundances of such ejecta, which are a fundamental input for chemical evolution codes, are usually taken from theoretical nucleosynthesis calculations performed by the nuclear astrophysics community. Therefore, almost all chemical evolution predictions rely on the nuclear physics behind those calculations. In this proceedings, we highlight the impact of nuclear physics uncertainties on galactic chemical evolution predictions. We demonstrate that nuclear physics and galactic evolution uncertainties both have a significant impact on interpreting the origin of neutron-capture elements in our Solar System. Those results serve as a motivation to create and maintain collaborations between the fields of nuclear astrophysics and galaxy evolution.

In large and complicated stellar systems like galaxies it is difficult to predict the number and characteristics of a black hole population. Such populations may be modelled as an aggregation of homogeneous (i.e. having uniform star formation history and the same initial chemical composition) stellar populations. Using realistic evolutionary models we predict the abundances and properties of black holes formed from binaries in these environments. We show that the black hole population will be dominated by single black holes originating from binary disruptions and stellar mergers. Furthermore, we discuss how black hole populations are influenced by such factors as initial parameters, metallicity, initial mass function, and natal kick models. As an example application of our results, we estimate that about 26 microlensing events to happen every year in the direction of the Galactic Bulge due to black holes in a survey like OGLE-IV. Our results may be used to perform in-depth studies related to realistic black hole populations, e.g. observational predictions for space survey missions like Gaia, or Einstein Probe. We prepared a publicly available database with the raw data from our simulations to be used for more in-depth studies.

Joint gravitational-wave detections of stellar-mass black-hole binaries by ground- and space-based observatories will provide unprecedented opportunities for fundamental physics and astronomy. We present a semianalytic method to estimate multiband event rates by combining selection effects of ground-based interferometers (like LIGO/Virgo) and space missions (like LISA). We forecast the expected number of multiband detections first by using information from current LIGO/Virgo data, and then through population synthesis simulations of binary stars. We estimate that few to tens of LISA detections can be used to predict mergers detectable on the ground. Conversely, hundreds of events could potentially be extracted from the LISA data stream using prior information from ground detections. In general, the merger signal of binaries observable by LISA is strong enough to be unambiguously identified by both current and future ground-based detectors. Therefore third-generation detectors will not increase the number of multiband detections compared to LIGO/Virgo. We use population synthesis simulations of isolated binary stars to explore some of the stellar physics that could be constrained with multiband events, and we show that specific formation pathways might be overrepresented in multiband events compared to ground-only detections.

We have analyzed how anisotropic emission of radiation affects the observed sample of ultraluminous X-ray sources (ULXs) by performing simulations of the evolution of stellar populations, employing recent developments in stellar and binary physics, and by utilizing a geometrical beaming model motivated by theory and observation. Whilst ULXs harboring black hole accretors (BH ULXs) are typically emitting isotropically, the majority of ULXs with neutron star accretors (NS ULXs) are found to be beamed. These findings confirm previous assertions that a significant fraction of ULXs are hidden from view due to a substantial misalignment of the emission beam and the line-of-sight. We find the total number of NS ULXs in regions with constant star formation, solar metallicity, and ages above ~1 Gyr to be higher than the BH ULXs, although observationally both populations are comparable. For lower metallicities BH ULX dominate both the total and observed ULX populations. As far as burst star-formation is concerned, young ULX populations are dominated by BH ULXs, but this changes as the population ages and, post star-formation, NS ULXs dominate both the observed and total population of ULXs. We also compare our simulation output to a previous analytical prediction for the relative ratio of BH to NS ULXs in idealized flux-limited observations and find broad agreement for all but the lowest metallicities. In so doing we find that in such surveys the observed ULX population should be heavily dominated by black-hole systems rather than by systems containing neutron stars.

Probing the origin of r-process elements in the universe represents a multi-disciplinary challenge. We review the observational evidence that probe the properties of r-process sites, and address them using galactic chemical evolution simulations, binary population synthesis models, and nucleosynthesis calculations. Our motivation is to define which astrophysical sites have significantly contributed to the total mass of r-process elements present in our Galaxy. We found discrepancies with the neutron star (NS-NS) merger scenario. Assuming they are the only site, the decreasing trend of [Eu/Fe] at [Fe/H]\,$>-1$ in the disk of the Milky Way cannot be reproduced while accounting for the delay-time distribution (DTD) of coalescence times ($\propto~t^{-1}$) derived from short gamma-ray bursts and population synthesis models. Steeper DTD functions ($\propto~t^{-1.5}$) or power laws combined with a strong burst of mergers before the onset of Type~Ia supernovae can reproduce the [Eu/Fe] trend, but this scenario is inconsistent with the similar fraction of short gamma-ray bursts and Type~Ia supernovae occurring in early-type galaxies, and reduces the probability of detecting GW170817 in an early-type galaxy. One solution is to assume an extra production site of Eu that would be active in the early universe, but would fade away with increasing metallicity. If this is correct, this extra site could be responsible for roughly 50% of the Eu production in the early universe, before the onset of Type~Ia supernovae. Rare classes of supernovae could be this additional r-process source, but hydrodynamic simulations still need to ensure the conditions for a robust r-process pattern.

Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario for neutron star binaries, the system might pass through a phase where a neutron star spirals into the envelope of its giant star companion. These phases lead to mass accretion onto the neutron star. Accretion onto these common-envelope-phase neutron stars can eject matter that has undergone burning near to the neutron star surface. This paper presents nucleosynthetic yields of this ejected matter, using population synthesis models to study the importance of these nucleosynthetic yields in a galactic chemical evolution context. Depending on the extreme conditions in temperature and density found in the accreted material, both proton-rich and neutron-rich nucleosynthesis can be obtained, with efficient production of neutron rich isotopes of low Z material at the most extreme conditions, and proton rich isotopes, again at low Z, in lower density models. Final yields are found to be extremely sensitive to the physical modeling of the accretion phase. We show that neutron stars accreting in binary common envelopes might be a new relevant site for galactic chemical evolution, and therefore more comprehensive studies are needed to better constrain nucleosynthesis in these objects.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

We investigate evolutionary pathways leading to neutron star formation through the collapse of oxygen-neon white dwarf (ONe WD) stars in interacting binaries. We consider (1) non-dynamical mass transfer where an ONe WD approaches the Chandrasekhar mass leading to accretion-induced collapse (AIC) and (2) dynamical timescale merger-induced collapse (MIC) between an ONe WD and another WD. We present rates, delay times, and progenitor properties for two different treatments of common envelope evolution. We show that AIC neutron stars are formed via many different channels and the most dominant channel depends on the adopted common envelope physics. Most AIC and MIC neutron stars are born shortly after star formation, though some have delay times >10 Gyr. The shortest delay time (25-50 Myr) AIC neutron stars have stripped-envelope, compact, helium-burning star donors, though many prompt AIC neutron stars form via wind-accretion from an asymptotic giant branch star. The longest delay time AIC neutron stars, which may be observed as young milli-second pulsars among globular clusters, have a red giant or main sequence donor at the time of NS formation and will eventually evolve into NS + helium WD binaries. We discuss AIC & MIC binaries as potential gravitational wave sources for LISA. Neutron stars created via AIC undergo a LMXB phase, offering an electromagnetic counterpart for those shortest orbital period sources that LISA could identify. The formation of neutron stars from interacting WDs in binaries is likely to be a key mechanism for the production of LIGO/Virgo gravitational wave sources (NS-NS and BH-NS mergers) in globular clusters.

(Abridged) The explosion mechanism of electron-capture supernovae (ECSNe) remains equivocal. We attempt to constrain the explosion mechanism (neutron-star-forming implosion or thermonuclear explosion) and the frequency of occurrence of ECSNe using nucleosynthesis simulations of the latter scenario, population synthesis, the solar abundance distribution, pre-solar meteoritic oxide grain isotopic ratio measurements and the white dwarf mass-radius relation. Tracer particles from 3d hydrodynamic simulations were post-processed with a large nuclear reaction network in order to determine the complete compositional state of the bound ONeFe remnant and the ejecta, and population synthesis simulations were performed in order to estimate the ECSN rate with respect to the CCSN rate. The 3d deflagration simulations drastically overproduce the neutron-rich isotopes $^{48}$Ca, $^{50}$Ti, $^{54}$Cr, $^{60}$Fe and several of the Zn isotopes relative to their solar abundances. Using the solar abundance distribution as our constraint, we place an upper limit on the frequency of thermonuclear ECSNe as 1$-$3~\% the frequency at which core-collapse supernovae (FeCCSNe) occur. This is on par with or 1~dex lower than the estimates for ECSNe from single stars. The upper limit from the yields is also in relatively good agreement with the predictions from our population synthesis simulations. The $^{54}$Cr/$^{52}$Cr and $^{50}$Ti/$^{48}$Ti isotopic ratios in the ejecta are a near-perfect match with recent measurements of extreme pre-solar meteoritc oxide grains, and $^{53}$Cr/$^{52}$Cr can also be matched if the ejecta condenses before mixing with the interstellar medium. Theoretical mass-radius relations for the bound ONeFe WD remnants of these explosions are apparently consistent with several observational WD candidates.

The first neutron star-neutron star merger was detected in by LIGO/Virgo in a galaxy in which the majority of star formation was taking place a long time ago (11 Gyr). LIGO/Virgo estimated that local cosmic NS-NS merger rate is 110-3840 Gpc^-3 yr^-1. Only some extreme evolutionary models can generate NS-NS merger rates in old host galaxies consistent with the LIGO/Virgo estimate. However, we show that these models generate rates exceeding empirical Galactic NS-NS merger rates based on the large population of Milky Way NS-NS binaries. Typically, current evolutionary models produce NS-NS merger rates that are consistent with the Milky Way empirical rates. However, these models generate local cosmic NS-NS merger rate in old host galaxies that are below the LIGO/Virgo estimate. The reason behind this tension is the predicted delay time distribution between star formation and NS-NS mergers that favors short delays. Evolutionary models produce a generic steep power-law NS-NS delay time distribution. This limits NS-NS merger rates in old host galaxies. However, we show that such distribution is consistent with observations of Galactic NS-NS binaries; 50% of which show very long merger times (much longer than Hubble time). Once model distributions are convolved with continuous prolonged (10 Gyr) star formation in the Galactic disk, then 20-70% (depending on a model) of the predicted NS-NS population has very long current Galactic merger times. Although NS-NS binaries are formed predominantly with short delay times, many of short delay time systems merge and do not make it to the present, while long delay time systems survive and contribute to the current Galactic NS-NS population. This study highlights the tension between the current evolutionary predictions and the observation of the first NS-NS merger in an old host galaxy.

Symbiotic stars are interacting binaries with one of the longest orbital periods. Since they can contain a massive white dwarf with a high accretion rate they are considered a possible type Ia supernovae (SNe Ia) progenitors. Among symbiotic binaries there are systems containing a Mira donor, which can have orbital periods of a few tens of years and more. This subclass of symbiotic stars due to their very large separation usually was not considered promising SNe Ia progenitors. We analysed evolution of one of the well studied symbiotic star with a Mira donor, V407 Cyg. We showed that the standard evolution model predicts that the system will not become a SN Ia. However, by simply adding a Wind Roche Lobe Overflow as one of the mass transfer modes we predict that the white dwarf in V407 Cyg will reach the Chandrasekhar limit in 40-200 Myr.

The distributions of the initial main-sequence binary parameters are one of the key ingredients in obtaining evolutionary predictions for compact binary (BH-BH / BH-NS / NS-NS) merger rates. Until now, such calculations were done under the assumption that initial binary parameter distributions were independent. Here, we implement empirically derived inter-correlated distributions of initial binary parameters primary mass (M1), mass ratio (q), orbital period (P), and eccentricity (e). Unexpectedly, the introduction of inter-correlated initial binary parameters leads to only a small decrease in the predicted merger rates by a factor of 2 $-$ 3 relative to the previously used non-correlated initial distributions. The formation of compact object mergers in the isolated classical binary evolution favors initial binaries with stars of comparable masses (q = 0.5 $-$ 1) at intermediate orbital periods (log P (days) = 2 $-$ 4). New distributions slightly shift the mass ratios towards smaller values with respect to the previously used flat q distribution, which is the dominant effect decreasing the rates. New orbital periods only negligibly increase the number of progenitors. Additionally, we discuss the uncertainty of merger rate predictions associated with possible variations of the massive-star initial mass function (IMF). We argue that evolutionary calculations should be normalized to a star formation rate (SFR) that is obtained from the observed amount of UV light at wavelength 1500{\AA} (SFR indicator). In this case, contrary to recent reports, the uncertainty of the IMF does not affect the rates by more than a factor of 2. Any change to the IMF slope for massive stars requires a change of SFR in a way that counteracts the impact of IMF variations on the merger rates. In contrast, we suggest that the uncertainty in cosmic SFR at low metallicity can be a significant factor at play.

Binaries that merge within the local Universe originate from progenitor systems that formed at different times and in various environments. The efficiency of formation of double compact objects is highly sensitive to metallicity of the star formation. Therefore, to confront the theoretical predictions with observational limits resulting from gravitational waves observations one has to account for the formation and evolution of progenitor stars in the chemically evolving Universe. In particular, this requires knowledge of the distribution of cosmic star formation rate at different metallicities and times, probed by redshift (SFR(Z,z)). We investigate the effect of the assumed SFR(Z,z) on the properties of merging double compact objects, in particular on their merger rate densities. Using a set of binary evolution models from Chruslinska et al. (2018) we demonstrate that the reported tension between the merger rates of different types of double compact objects and current observational limits in some cases can be resolved if a SFR(Z,z) closer to that expected based on observations of local star-forming galaxies is used, without the need for changing the assumptions about the evolution of progenitor stars of different masses. This highlights the importance of finding tighter constraints on SFR(Z,z) and understanding the associated uncertainties.

We study the expected spin misalignments of merging binary black holes (BHs) formed in isolation by combining state-of-the-art population-synthesis models with efficient post-Newtonian evolutions, thus tracking sources from stellar formation to gravitational-wave detection. We present extensive predictions of the properties of sources detectable by both current and future interferometers. We account for the fact that detectors are more sensitive to spinning BH binaries with suitable spin orientations and find that this significantly impacts the population of sources detectable by LIGO, while this is not the case for 3rd-generation detectors. We find that three formation pathways, differentiated by the order of core collapse and common-envelope phases, dominate the observed population, and that their relative importance critically depends on the recoils imparted to BHs at birth. Our models suggest that measurements of the "effective-spin" parameter $\chi_{\rm eff}$ will allow for powerful constraints. For instance, we find that the role of spin magnitudes and spin directions in $\chi_{\rm eff}$ can be largely disentangled, and that the symmetry of the effective-spin distribution is a robust indicator of the binary's formation history. Our predictions for individual spin directions and their precessional morphologies confirm and extend early toy models, while exploring substantially more realistic and broader sets of initial conditions. Our main conclusion is that specific subpopulations of BH binaries will exhibit distinctive precessional dynamics: these classes include (but are not limited to) sources where stellar tidal interactions act on sufficiently short timescales, and massive binaries produced in pulsational pair-instability supernovae. Measurements of BH spin orientations have enormous potential to constrain specific evolutionary processes in the lives of massive binary stars.

Anisotropy of gravitational wave (GW) emission results in a net momentum gained by the black hole (BH) merger product, leading to a recoil velocity up to $\sim10^3\text{ km s}^{-1}$, which may kick it out of a globular cluster (GC). We estimate GW kick retention fractions of merger products assuming different models for BH spin magnitude and orientation (MS0 - random, MS1 - spin as a function of mass and metalicity, MS2 - constant value of $0.5$). We check how they depend on BH-BH merger time and properties of the cluster. We analyze the implications of GW kick retention fractions on intermediate massive BH (IMBH) formation by repeated mergers in a GC. We also calculate final spin of the merger product, and investigate how it correlates with effective spin of the binary. We used data from MOCCA (MOnte Carlo Cluster simulAtor) GC simulations to get a realistic sample of BH-BH mergers, assigned each BH spin value according to a studied model, and calculated recoil velocity and final spin based on most recent theoretical formulas. We discovered that for physically motivated models, GW kick retention fractions are about $30\%$ and display small dependence on assumptions about spin, but are much more prone to cluster properties. In particular, we discovered a strong dependence of GW kick retention fractions on cluster density. We also show that GW kick retention fractions are high in final life stages of the cluster, but low at the beginning. Finally, we derive formulas connecting final spin with effective spin for primordial binaries, and with maximal effective spin for dynamical binaries.

Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physical processes such as mass transfer stability criteria, supernova kick distributions, remnant mass distributions and common-envelope energetics. We compute the likelihood of observing the orbital period-eccentricity distribution of the Galactic DNS population under each of our population synthesis models, allowing us to quantitatively compare the models. We find that mass transfer from a stripped post-helium-burning secondary (case BB) onto a neutron star is most likely dynamically stable. We also find that a natal kick distribution composed of both low (Maxwellian $\sigma=30\rm~km~s^{-1}$) and high ($\sigma=265\rm~km~s^{-1}$) components is preferred over a single high-kick component. We find that the observed DNS mass distribution can place strong constraints on model assumptions.

We compare binary evolution models with different assumptions about black-hole natal kicks to the first gravitational-wave observations performed by the LIGO detectors. Our comparisons attempt to reconcile merger rate, masses, spins, and spin-orbit misalignments of all current observations with state-of-the-art formation scenarios of binary black holes formed in isolation. We estimate that black holes (BHs) should receive natal kicks at birth of the order of $\sigma\simeq 200$ (50) km/s if tidal processes do (not) realign stellar spins. Our estimate is driven by two simple factors. The natal kick dispersion $\sigma$ is bounded from above because large kicks disrupt too many binaries (reducing the merger rate below the observed value). Conversely, the natal kick distribution is bounded from below because modest kicks are needed to produce a range of spin-orbit misalignments. A distribution of misalignments increases our models' compatibility with LIGO's observations, if all BHs are likely to have natal spins. Unlike related work which adopts a concrete BH natal spin prescription, we explore a range of possible BH natal spin distributions. Within the context of our models, for all of the choices of $\sigma$ used here and within the context of one simple fiducial parameterized spin distribution, observations favor low BH natal spin.

Binary black holes (BBHs) are one of the endpoints of isolated binary evolution, and their mergers a leading channel for gravitational wave events. Here, using the evolutionary code \textsc{StarTrack}, we study the statistical properties of the BBH population from isolated binary evolution for a range of progenitor star metallicities and BH natal kicks. We compute the mass function and the distribution of the primary BH spin $a$ as a result of mass accretion during the binary evolution, and find that this is not an efficient process to spin up BHs, producing an increase by at most $a\sim$~0.2--0.3 for very low natal BH spins. We further compute the distribution of merger sites within the host galaxy, after tracking the motion of the binaries in the potentials of a massive spiral, a massive elliptical, and a dwarf galaxy. We find that a fraction of 70-90\% of mergers in massive galaxies and of 40-60\% in dwarfs (range mostly sensitive to the natal kicks) is expected to occur inside of their hosts. The number density distribution at the merger sites further allows us to estimate the broadband luminosity distribution that BBH mergers would produce, \textit{if} associated with a kinetic energy release in an outflow, {which, as a reference, we assume at the level inferred for the \textit{Fermi} GBM counterpart to GW150914, with the understanding that current limits from the O1 and O2 runs would require such emission to be produced within a jet of angular size within $\lesssim 50^\circ$.}

The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms of double compact objects: classical isolated binary star evolution, dynamical evolution in globular clusters and nuclear cluster formation to test whether they are likely to produce NS-NS mergers in old host galaxies. Our simulations with optimistic assumptions show current NS-NS merger rates at the level of 10^-2 yr^-1 from binary stars, 5 x 10^-5 yr^-1 from globular clusters and 10^-5 yr^-1 from nuclear clusters for all local elliptical galaxies (within 100 Mpc^3). These models are thus in tension with the detection of GW170817 with an observed rate 1.5 yr^-1 (per 100 Mpc^3; LIGO/Virgo estimate). Our results imply that either (i) the detection of GW170817 by LIGO/Virgo at their current sensitivity in an elliptical galaxy is a statistical coincidence; or that (ii) physics in at least one of our three models is incomplete in the context of the evolution of stars that can form NS-NS mergers; or that (iii) another very efficient (unknown) formation channel with a long delay time between star formation and merger is at play.

Rapidly accreting white dwarfs (RAWDs) have been proposed as contributors to the chemical evolution of heavy elements in the Galaxy. Here, we test this scenario for the first time and determine the contribution of RAWDs to the solar composition of first-peak neutron-capture elements. We add the metallicity-dependent contribution of RAWDs to the one-zone galactic chemical evolution code OMEGA according to RAWD rates from binary stellar population models combined with metallicity-dependent $i$-process stellar yields calculated following the models of Denissenkov et al. (2017). With this approach we find that the contribution of RAWDs to the evolution of heavy elements in the Galaxy could be responsible for a significant fraction of the solar composition of Kr, Rb, Sr, Y, Zr, Nb, and Mo ranging from $2$ to $45\%$ depending on the element, the enrichment history of the Galactic gas, and the total mass ejected per RAWD. This contribution could explain the missing solar Lighter Element Primary Process for some elements (e.g., Sr, Y, and Zr). We do not overproduce any isotope relative to the solar composition, but $^{96}$Zr is produced in a similar amount. The $i$ process produces efficiently the Mo stable isotopes $^{95}$Mo and $^{97}$Mo. When nuclear reaction rate uncertainties are combined with our GCE uncertainties, the upper limits for the predicted RAWD contribution increase by a factor of $1.5-2$ for Rb, Sr, Y, and Zr, and by 3.8 and 2.4 for Nb and Mo, respectively. We discuss the implication of the RAWD stellar evolution properties on the single degenerate Type Ia supernova scenario.

Some of the heavy elements, such as gold and europium (Eu), are almost exclusively formed by the rapid neutron capture process (r-process). However, it is still unclear which astrophysical site between core-collapse supernovae and neutron star - neutron star (NS-NS) mergers produced most of the r-process elements in the universe. Galactic chemical evolution (GCE) models can test these scenarios by quantifying the frequency and yields required to reproduce the amount of europium (Eu) observed in galaxies. Although NS-NS mergers have become popular candidates, their required frequency (or rate) needs to be consistent with that obtained from gravitational wave measurements. Here we address the first NS-NS merger detected by LIGO/Virgo (GW170817) and its associated Gamma-ray burst and analyze their implication on the origin of r-process elements. The range of NS-NS merger rate densities of 320-4740 Gpc$^{-3}$ yr$^{-1}$ provided by LIGO/Virgo is remarkably consistent with the range required by GCE to explain the Eu abundances in the Milky Way with NS-NS mergers, assuming the solar r-process abundance pattern for the ejecta. Under the same assumption, this event has produced about 1-5 Earth masses of Eu, and 3-13 Earth masses of gold. When using theoretical calculations to derive Eu yields, constraining the role of NS-NS mergers becomes more challenging because of nuclear astrophysics uncertainties. This is the first study that directly combines nuclear physics uncertainties with GCE calculations. If GW170817 is a representative event, NS-NS mergers can produce Eu in sufficient amounts and are likely to be the main r-process site.

We revisit double neutron star (DNS) formation in the classical binary evolution scenario in light of the recent LIGO/Virgo DNS detection (GW170817). The observationally estimated Galactic DNS merger rate of $R_{\rm MW}=21^{+28}_{-14}$ Myr$^{-1}$, based on 3 Galactic DNS systems, fully supports our standard input physics model with $R_{\rm MW} =24$ Myr$^{-1}$. This estimate for the Galaxy translates in a non-trivial way (due to cosmological evolution of progenitor stars in chemically evolving Universe) into a local ($z\approx0$) DNS merger rate density of $R_{\rm local}=48$ Gpc$^{-3}$yr$^{-1}$, which {\em is not} consistent with the current LIGO/Virgo DNS merger rate estimate ($1540^{+3200}_{-1220}$ Gpc$^{-3}$yr$^{-1}$). Within our study of the parameter space we find solutions that allow for DNS merger rates as high as $R_{\rm local} \approx 600^{+600}_{-300}$ Gpc$^{-3}$yr$^{-1}$ which are thus consistent with the LIGO/Virgo estimate. However, our corresponding BH-BH merger rates for the models with high DNS merger rates exceed the current LIGO/Virgo estimate of local BH-BH merger rate ($12$-$213$ Gpc$^{-3}$yr$^{-1}$). Apart from being particularly sensitive to the common envelope treatment, DNS merger rates are rather robust against variations of several of the key factors probed in our study (e.g. mass transfer, angular momentum loss, natal kicks). This might suggest that either common envelope development/survival works differently for DNS ($\sim$ 10-20 Msun stars) than for BH-BH ($\sim$ 40-100 Msun stars) progenitors, or high BH natal kicks are needed to meet observational constraints for both types of binaries. Note that our conclusion is based on a limited number of (21) evolutionary models and is valid only within this particular DNS and BH-BH isolated binary formation scenario.

Recently, several ultraluminous X-ray (ULX) sources were shown to host a neutron star (NS) accretor. We perform a suite of evolutionary calculations which show that, in fact, NSs are the dominant type of ULX accretor. Although black holes (BH) dominate early epochs after the star-formation burst, NSs outweigh them after a few 100 Myr and may appear as late as a few Gyr after the end of the star formation episode. If star formation is a prolonged and continuous event (i.e., not a relatively short burst), NS accretors dominate ULX population at any time in solar metallicity environment, whereas BH accretors dominate when the metallicity is sub-solar. Our results show a very clear (and testable) relation between the companion/donor evolutionary stage and the age of the system. A typical NS ULX consists of a $\sim1.3\,M_\odot$ NS and $\sim1.0\,M_\odot$ Red Giant. A typical BH ULX consist of a $\sim8\,M_\odot$ BH and $\sim6\,M_\odot$ main-sequence star. Additionally, we find that the very luminous ULXs ($L_X\gtrsim10^{41}$ erg/s) are predominantly BH systems ($\sim9\,M_\odot$) with Hertzsprung gap donors ($\sim2\,M_\odot$). Nevertheless, some NS ULX systems may also reach extremely high X-ray luminosities ($\gtrsim10^{41}$ erg/s).

Very wide binaries (> 500 AU) are subject to numerous encounters with flying-by stars in the Galactic field and can be perturbated into highly eccentric orbits (e ~ 0.99). For such systems tidal interactions at close pericenter passages can lead to orbit circularization and possibly mass transfer, consequently producing X-Ray binaries without the need for common envelope. We test this scenario for the case of Black Hole Low-Mass X-Ray Binaries (BH LMXBs) by performing a population synthesis from primordial binaries with numerical treatment of random stellar encounters. We test various models for the threshold pericenter distance under which tidal forces cause circularization. We estimate that fly-by interactions can produce a current population of ~ 60$-$220 BH LMXBs in the Galactic field. The results are sensitive to the assumption on tidal circularization efficiency and zero to very small BH natal kicks of a few km/s are required. We show that the most likely donors are low-mass stars (< 1 Msun; at the onset of mass transfer) as observed in the population of known sources (~ 20). However, the low number of systems formed along this route is in tension with most recent observational estimate of the number of dormant BH LMXBs in the Galaxy 10$^4-$10$^8$ (Tetarenko et al. 2016). If indeed the numbers are so high, alternative formation channels of BHs with low-mass donors need to be identified.

The role of compact binary mergers as the main production site of r-process elements is investigated by combining stellar abundances of Eu observed in the Milky Way, galactic chemical evolution (GCE) simulations, binary population synthesis models, and Advanced LIGO gravitational wave measurements. We compiled and reviewed seven recent GCE studies to extract the frequency of neutron star - neutron star (NS-NS) mergers that is needed in order to reproduce the observed [Eu/Fe] vs [Fe/H] relationship. We used our simple chemical evolution code to explore the impact of different analytical delay-time distribution (DTD) functions for NS-NS mergers. We then combined our metallicity-dependent population synthesis models with our chemical evolution code to bring their predictions, for both NS-NS mergers and black hole - neutron star mergers, into a GCE context. Finally, we convolved our results with the cosmic star formation history to provide a direct comparison with current and upcoming Advanced LIGO measurements. When assuming that NS-NS mergers are the exclusive r-process sites, and that the ejected r-process mass per merger event is 0.01 Msun, the number of NS-NS mergers needed in GCE studies is about 10 times larger than what is predicted by standard population synthesis models. These two distinct fields can only be consistent with each other when assuming optimistic rates, massive NS-NS merger ejecta, and low Fe yields for massive stars. Population synthesis models and GCE simulations are in agreement with the current upper limit (O1) established by Advanced LIGO during their first run of observations. Upcoming measurements will provide important constraints on the actual local NS-NS merger rate, will provide insights on the plausibility of the GCE requirement, and will help to define whether or not compact binary mergers can be the dominant source of r-process elements in the Universe.

Single star evolution does not allow extremely low-mass stars to cross the classical instability strip (IS) during the Hubble time. However, within binary evolution framework low-mass stars can appear inside the IS once the mass transfer (MT) is taken into account. Triggered by a discovery of low-mass 0.26 Msun RR Lyrae-like variable in a binary system, OGLE-BLG-RRLYR-02792, we investigate the occurrence of similar binary components in the IS, which set up a new class of low-mass pulsators. They are referred to as Binary Evolution Pulsators (BEPs) to underline the interaction between components, which is crucial for substantial mass loss prior to the IS entrance. We simulate a population of 500 000 metal-rich binaries and report that 28 143 components of binary systems experience severe MT (loosing up to 90% of mass), followed by at least one IS crossing in luminosity range of RR Lyrae (RRL) or Cepheid variables. A half of these systems enter the IS before the age of 4 Gyr. BEPs display a variety of physical and orbital parameters, with the most important being the BEP mass in range 0.2-0.8 Msun, and the orbital period in range 10-2500 d. Based on the light curve only, BEPs can be misclassified as genuine classical pulsators, and as such they would contaminate genuine RRL and classical Cepheid variables at levels of 0.8% and 5%, respectively. We state that the majority of BEPs will remain undetected and we discuss relevant detection limitations.

We study the contribution of binary black hole (BH-BH) mergers from the first, metal-free stars in the Universe (Pop III) to gravitational wave detection rates. Our study combines initial conditions for the formation of Pop III stars based on N-body simulations of binary formation (including rates, binary fraction, initial mass function, orbital separation and eccentricity distributions) with an updated model of stellar evolution specific for Pop III stars. We find that the merger rate of these Pop III BH-BH systems is relatively small (< 0.1 Gpc^-3 yr^-1) at low redshifts (z<2), where it can be compared with the LIGO empirical estimate of 9-240 Gpc^-3 yr^-1 (Abbott et al. 2016). The predicted rates are even smaller for Pop III double neutron star and black hole neutron star mergers. Our rates are compatible with those of Hartwig et al. (2016), but significantly smaller than those found in previous work (Bond & Carr 1984; Belczynski et al. 2004; Kinugawa et al. 2014, 2016). We explain the reasons for this discrepancy by means of detailed model comparisons and point out that (i) identification of Pop III BH-BH mergers may not be possible by advanced LIGO, and (ii) the level of stochastic gravitational wave background from Pop III mergers may be lower than recently estimated (Kowalska et al. 2012; Inayoshi et al. 2016; Dvorkin et al. 2016). We further estimate gravitational wave detection rates for third-generation interferometric detectors. Our calculations are relevant for low to moderately rotating Pop III stars. We can now exclude significant (> 1 per cent) contribution of these stars to low-redshift BH-BH mergers. However, it remains to be tested whether (and at what level) rapidly spinning Pop III stars (homogeneous evolution) can contribute to BH-BH mergers in the local Universe.

We employ population synthesis method to model the double neutron star (DNS) population and test various possibilities on natal kick velocities gained by neutron stars after their formation. We first choose natal kicks after standard core collapse SN from a Maxwellian distribution with velocity dispersion of sigma=265 km/s as proposed by Hobbs et al. (2005) and then modify this distribution by changing the velocity dispersion towards smaller and larger kick values. We also take into account the possibility of NS formation through electron capture supernova. In this case we test two scenarios: zero natal kick or small natal kick, drawn from Maxwellian distribution with sigma = 26.5 km/s. We calculate the present-day orbital parameters of binaries and compare the resulting eccentricities with those known for observed DNSs. As an additional test we calculate Galactic merger rates for our model populations and confront them with observational limits. We do not find any model unequivocally consistent with both observational constraints simultaneously.The models with low kicks after CCSN for binaries with the second NS forming through core collapse SN are marginally consistent with the observations. This means that either 14 observed DNSs are not representative of the intrinsic Galactic population, or that our modeling of DNS formation needs revision.

We estimate the potential of present and future interferometric gravitational-wave detectors to test the Kerr nature of black holes through "gravitational spectroscopy," i.e. the measurement of multiple quasinormal mode frequencies from the remnant of a black hole merger. Using population synthesis models of the formation and evolution of stellar-mass black hole binaries, we find that Voyager-class interferometers will be necessary to perform these tests. Gravitational spectroscopy in the local Universe may become routine with the Einstein Telescope, but a 40-km facility like Cosmic Explorer is necessary to go beyond $z\sim 3$. In contrast, eLISA-like detectors should carry out a few - or even hundreds - of these tests every year, depending on uncertainties in massive black hole formation models. Many space-based spectroscopical measurements will occur at high redshift, testing the strong gravity dynamics of Kerr black holes in domains where cosmological corrections to general relativity (if they occur in nature) must be significant.

Mergers of two stellar origin black holes are a prime source of gravitational waves and are under intensive investigations. One crucial ingredient in their modeling has so far been neglected. Pair-instability pulsation supernovae with associated severe mass loss may suppress formation of massive black holes, decreasing black hole merger rates for the highest black hole masses. The mass loss associated with pair-instability pulsation supernovae limits the Population I/II stellar-origin black hole mass to 50 Msun, in tension with earlier predictions that the maximum black hole mass could be as high as 100 Msun. Suppression of double black hole merger rates by pair-instability pulsation supernovae is negligible for our evolutionary channel. Our standard evolutionary model with inclusion of pair-instability pulsation supernovae and with pair-instability supernovae is fully consistent with the LIGO detections of black hole mergers: GW150914, GW151226 and LVT151012. The LIGO observations seem to exclude high (>400 km/s) BH natal kicks. We predict the detection of several up to 60 BH-BH mergers with total mass 10--150 Msun (most likely range: 20--80 Msun) in the forthcoming 60 effective days of LIGO O2 observations.

The merger of two massive 30 Msun black holes has been detected in gravitational waves (1,GW150914). This discovery validates recent predictions (2-4) that massive binary black holes would constitute the first detection. However, previous calculations have not sampled the relevant binary black hole progenitors---massive, low-metallicity binary stars---with sufficient accuracy and input physics to enable robust predictions to better than several orders of magnitude (5-10). Here we report a suite of high-precision numerical simulations of binary black hole formation via the evolution of isolated binary stars, providing a framework to interpret GW150914 and predict the properties of subsequent binary black hole gravitational-wave events. Our models imply that these events form in an environment where the metallicity is less than 10 percent of solar; have initial masses of 40-100 Msun; and interact through mass transfer and a common envelope phase. Their progenitors likely form either at 2 Gyr, or somewhat less likely, at 11 Gyr after the Big Bang. Most binary black holes form without supernova explosions, and their spins are nearly unchanged since birth, but do not have to be parallel. The classical field formation of binary black holes proposed in this study, with low natal kicks and restricted common envelope evolution, produces 40 times more binary black holes than dynamical formation channels involving globular clusters (11) and is comparable to the rate from homogeneous evolution channels (12-15). Our calculations predict detections of about 1,000 black hole mergers per year with total mass of 20-80 Msun once second generation ground-based gravitational wave observatories reach full sensitivity.

The mass transfer in binaries with massive donors and compact companions, when the donors rapidly evolve after their main sequence, is one of the dominant formation channels of merging double stellar-mass black hole binaries. This mass transfer was previously postulated to be unstable and was expected to lead to a common envelope event. The common envelope event then would end with either double black hole formation, or with the merger of the two stars. We re-visit the stability of this mass transfer, and find that for a large range of the binary orbital separations this mass transfer is stable. This newly found stability allows us to reconcile the theoretical rate for double black hole binary mergers predicted by population synthesis studies, and the empirical rate obtained by LIGO. Futhermore, the stability of the mass transfer leads to the formation of ultra-luminous X-ray sources. The theoretically predicted formation rates of ultra-luminous X-ray sources powered by a stellar-mass BH, as well as the range of produced X-ray luminosity, can explain the observed bright ultra-luminous X-ray sources.

We compare evolutionary predictions of double compact object merger rate densities with initial and forthcoming LIGO/Virgo upper limits. We find that: (i) Due to the cosmological reach of advanced detectors, current conversion methods of population synthesis predictions into merger rate densities are insufficient. (ii) Our optimistic models are a factor of 18 below the initial LIGO/Virgo upper limits for BH-BH systems, indicating that a modest increase in observational sensitivity (by a factor of 2.5) may bring the first detections or first gravitational wave constraints on binary evolution. (iii) Stellar-origin massive BH-BH mergers should dominate event rates in advanced LIGO/Virgo and can be detected out to redshift z=2 with templates including inspiral, merger, and ringdown. Normal stars (<150 Msun) can produce such mergers with total redshifted mass up to 400 Msun. (iv) High black hole natal kicks can severely limit the formation of massive BH-BH systems (both in isolated binary and in dynamical dense cluster evolution), and thus would eliminate detection of these systems even at full advanced LIGO/Virgo sensitivity. We find that low and high black hole natal kicks are allowed by current observational electromagnetic constraints. (v) The majority of our models yield detections of all types of mergers (NS-NS, BH-NS, BH-BH) with advanced detectors. Numerous massive BH-BH merger detections will indicate small (if any) natal kicks for massive BHs.

We have studied the future evolution of a recently discovered ULX source P13 in NGC7793. This source was shown to contain a 5-15$M_{\odot}$ black hole and a massive 18-23$M_{\odot}$ B9Ia companion on a 64 day orbit. For low black hole mass (5-10$M_{\odot}$) and high companion mass $\gtrsim 20 M_{\odot}$ the binary is predicted to initiate a common envelope evolution in near future and significantly decrease the orbital separation (6 hr orbit). This leads to a high probability ($\sim70\%$) of the system surviving a supernova explosion that will form a neutron star out of the companion. About one third of the surviving BH-NS systems will merge within Hubble time and be a source of high frequency gravitational radiation. We estimate that the chances of detection of BH-NS systems with advanced LIGO/Virgo that form via P13-evolutionary channel are at the level of 0.1 yr$^{-1}$ with wide range of allowed probability $0 \div 0.6$yr$^{-1}$. This is the fourth empirical estimate of BH-NS merger rate.

SMC 3 is one of the most interesting symbiotic stars. This binary contains a bright K-type giant transferring mass to a massive white dwarf comanion, which makes it is a very promising SN Ia candidate. We discuss the evolutionary status of the system using results of population synthesis code.

The initial mass function (IMF), binary fraction and distributions of binary parameters (mass ratios, separations and eccentricities) are indispensable input for simulations of stellar populations. It is often claimed that these are poorly constrained significantly affecting evolutionary predictions. Recently, dedicated observing campaigns provided new constraints on the initial conditions for massive stars. Findings include a larger close binary fraction and a stronger preference for very tight systems. We investigate the impact on the predicted merger rates of neutron stars and black holes. Despite the changes with previous assumptions, we only find an increase of less than a factor 2 (insignificant compared with evolutionary uncertainties of typically a factor 10-100). We further show that the uncertainties in the new initial binary properties do not significantly affect (within a factor of 2) our predictions of double compact object merger rates. An exception is the uncertainty in IMF (variations by a factor of 6 up and down). No significant changes in the distributions of final component masses, mass ratios, chirp masses and delay times are found. We conclude that the predictions are, for practical purposes, robust against uncertainties in the initial conditions concerning binary parameters with exception of the IMF. This eliminates an important layer of the many uncertain assumptions affecting the predictions of merger detection rates with the gravitational wave detectors aLIGO/aVirgo.

The cosmic evolution of the neutron star merger (NSM) rate can be deduced from the observed cosmic star formation rate. This allows to estimate the rate expected in the horizon of the gravitational wave detectors advanced Virgo and ad LIGO and to compare those rates with independent predictions. In this context, the rapid neutron capture process, or r process, can be used as a constraint assuming NSM is the main astrophysical site for this nucleosynthetic process. We compute the early cosmic evolution of a typical r process element, Europium. Eu yields from NSM are taken from recent nucleosynthesis calculations. The same approach allows to compute the cosmic rate of Core Collapse SuperNovae (CCSN) and the associated evolution of Eu. We find that the bulk of Eu observations at high iron abundance can be rather well fitted by either CCSN or NSM scenarios. However, at lower metallicity, the early Eu cosmic evolution favors NSM as the main astrophysical site for the r process. A comparison between our calculations and spectroscopic observations at very low metallicities allows to constrain the coalescence timescale in the NSM scenario to about 0.1 to 0.2 Gyr. These values are in agreement with the coalescence timescales of some observed binary pulsars. Finally, the cosmic evolution of Eu is used to put constraints on the NSM rate, the merger rate in the horizon of the gravitational wave detectors advanced Virgo/ad LIGO, as well as the expected rate of electromagnetic counterparts to mergers (kilonovae) in large near-infrared surveys.

In this proof-of-concept study we demonstrate that in a binary system mass can be transferred toward an accreting compact object at extremely high rate. If the transferred mass is efficiently converted to X-ray luminosity (with disregard of the classical Eddington limit) or if the X-rays are focused into a narrow beam then binaries can form extreme ULX sources with the X-ray luminosity of Lx>10^42 erg/s. For example, Lasota & King argued that the brightest known ULX (HLX-1) is a regular binary system with a rather low-mass compact object (a stellar-origin black hole or a neutron star). The predicted formation efficiencies and lifetimes of binaries with the very high mass transfer rates are large enough to explain all observed systems with extreme X-ray luminosities. These systems are not only limited to binaries with stellar-origin black hole accretors. Noteworthy, we have also identified such objects with neutron stars. Typically, a 10 Msun black hole is fed by a massive (10 Msun) Hertzsprung gap donor with Roche lobe overflow rate of 10^-3 Msun/yr (2600 MEdd). For neutron star systems the typical donors are evolved low-mass (2 Msun) helium stars with Roche lobe overflow rate of 10^-2 Msun/yr. Our study does not prove that any particular extreme ULX is a regular binary system, but it demonstrates that any ULX, including the most luminous ones, may potentially be a short-lived phase in the life of a binary star.

The unprecedented range of second-generation gravitational-wave (GW) observatories calls for refining the predictions of potential sources and detection rates. The coalescence of double compact objects (DCOs)---i.e., neutron star-neutron star (NS-NS), black hole-neutron star (BH-NS), and black hole-black hole (BH-BH) binary systems---is the most promising source of GWs for these detectors. We compute detection rates of coalescing DCOs in second-generation GW detectors using the latest models for their cosmological evolution, and implementing inspiral-merger-ringdown (IMR) gravitational waveform models in our signal-to-noise ratio calculations. We find that: (1) the inclusion of the merger/ringdown portion of the signal does not significantly affect rates for NS-NS and BH-NS systems, but it boosts rates by a factor $\sim 1.5$ for BH-BH systems; (2) in almost all of our models BH-BH systems yield by far the largest rates, followed by NS-NS and BH-NS systems, respectively, and (3) a majority of the detectable BH-BH systems were formed in the early Universe in low-metallicity environments. We make predictions for the distributions of detected binaries and discuss what the first GW detections will teach us about the astrophysics underlying binary formation and evolution.

We find that the formation of MWC 656 (the first Be binary containing a black hole) involves a common envelope phase and a supernova explosion. This result supports the idea that a rapidly rotating Be star can emerge out of a common envelope phase, which is very intriguing because this evolutionary stage is thought to be too fast to lead to significant accretion and spin up of the B star. We predict $\sim 10-100$ of B BH binaries to currently reside in the Galactic disk, among which around $1/3$ contain a Be star, but there is only a small chance to observe a system with parameters resembling MWC 656. If MWC 656 is representative of intrinsic Galactic Be BH binary population, it may indicate that standard evolutionary theory needs to be revised. This would pose another evolutionary problem in understanding BH binaries, with BH X-ray Novae formation issue being the prime example. The future evolution of MWC 656 with a $\sim 5$ M$_{\odot}$ black hole and with a $\sim 13$ M$_{\odot}$ main sequence companion on a $\sim 60$ day orbit may lead to the formation of a coalescing BH-NS system. The estimated Advanced LIGO/Virgo detection rate of such systems is up to $\sim 0.2$ yr$^{-1}$. This empirical estimate is a lower limit as it is obtained with only one particular evolutionary scenario, the MWC 656 binary. This is only a third such estimate available (after Cyg X-1 and Cyg X-3), and it lends additional support to the existence of so far undetected BH--NS binaries.

Superposition of gravitational waves generated by astrophysical sources is expected to give rise to the stochastic gravitational-wave background. We focus on the background generated by the ring-down of black holes produced in the stellar core collapse events across the universe. We systematically study the parameter space in this model, including the most recent information about the star formation rate and about the population of black holes as a function of redshift and of metallicity. We investigate the accessibility of this gravitational wave background to the upcoming gravitational-wave detectors, such as Advanced LIGO and Einstein Telescope.

We analyze the distinguishability of populations of coalescing binary neutron stars, neutron-star black-hole binaries, and binary black holes, whose gravitational-wave signatures are expected to be observed by the advanced network of ground-based interferometers LIGO and Virgo. We consider population-synthesis predictions for plausible merging binary distributions in mass space, along with measurement accuracy estimates from the main gravitational-wave parameter-estimation pipeline. We find that for our model compact-object binary mass distribution, we can always distinguish binary neutron stars and black-hole--neutron-star binaries, but not necessarily black-hole--neutron-star binaries and binary black holes; however, with a few tens of detections, we can accurately identify the three subpopulations and measure their respective rates.

Neutron star (binary neutron star and neutron star - black hole) mergers are believed to produce short-duration gamma-ray bursts. They are also believed to be the dominant source of gravitational waves to be detected by the advanced LIGO and the dominant source of the heavy r-process elements in the universe. Whether or not these mergers produce short-duration GRBs depends sensitively on the fate of the core of the remnant (whether, and how quickly, it forms a black hole). In this paper, we combine the results of merger calculations and equation of state studies to determine the fate of the cores of neutron star mergers. Using population studies, we can determine the distribution of these fates to compare to observations. We find that black hole cores form quickly only for equations of state that predict maximum non-rotating neutron star masses below 2.3-2.4 solar masses. If quick black hole formation is essential in producing gamma-ray bursts, LIGO observed rates compared to GRB rates could be used to constrain the equation of state for dense nuclear matter.

Recent studies on stellar evolution have shown that the properties of compact objects strongly depend on the metallicity of the environment in which they were formed. Using some very simple assumptions on the metallicity of the stellar populations, we explore how this property affects the unresolved gravitational-wave background from extragalactic compact binaries. We obtained a suit of models using population synthesis code, estimated the gravitational-wave background they produce, and discuss its detectability with second- (advanced LIGO, advanced Virgo) and third- (Einstein Telescope) generation detectors. Our results show that the background is dominated by binary black holes for all considered models in the frequency range of terrestrial detectors, and that it could be detected in most cases by advanced LIGO/Virgo, and with Einstein Telescope with a very high signal-to-noise ratio. The observed peak in a gravitational wave spectrum depends on the metallicity of the stellar population.

Stellar mass black hole binaries have individual masses between 10-80 solar masses. These systems may emit gravitational waves at frequencies detectable at Megaparsec distances by space-based gravitational wave observatories. In a previous study, we determined the selection effects of observing these systems with detectors similar to the Laser Interferometer Space Antenna by using a generated population of binary black holes that covered a reasonable parameter space and calculating their signal-to-noise ratio. We further our study by populating the galaxies in our nearby (less than 30 Mpc) universe with binary black hole systems drawn from a distribution found in the Synthetic Universe to ultimately investigate the likely event rate of detectable binaries from galaxies in the nearby universe.

If binaries consisting of two 100 Msun black holes exist they would serve as extraordinarily powerful gravitational-wave sources, detectable to redshifts of z=2 with the advanced LIGO/Virgo ground-based detectors. Large uncertainties about the evolution of massive stars preclude definitive rate predictions for mergers of these massive black holes. We show that rates as high as hundreds of detections per year, or as low as no detections whatsoever, are both possible. It was thought that the only way to produce these massive binaries was via dynamical interactions in dense stellar systems. This view has been challenged by the recent discovery of several stars with mass above 150 Msun in the R136 region of the Large Magellanic Cloud. Current models predict that when stars of this mass leave the main sequence, their expansion is insufficient to allow common envelope evolution to efficiently reduce the orbital separation. The resulting black-hole--black-hole binary remains too wide to be able to coalesce within a Hubble time. If this assessment is correct, isolated very massive binaries do not evolve to be gravitational-wave sources. However, other formation channels exist. For example, the high multiplicity of massive stars, and their common formation in relatively dense stellar associations, opens up dynamical channels for massive black hole mergers (e.g., via Kozai cycles or repeated binary-single interactions). We identify key physical factors that shape the population of very massive black-hole--black-hole binaries. Advanced gravitational-wave detectors will provide important constraints on the formation and evolution of very massive stars.

There are 19 confirmed BH binaries in the Galaxy. 16 of them are X-ray transients hosting a ~5-15 Msun BH and a Roche-lobe overflowing low-mass companion. Companion masses are found mostly in 0.1-1 Msun mass range with peak at 0.6 Msun. The formation of these systems is believed to involve a common envelope phase, initiated by a BH progenitor, expected to be a massive star >20 Msun. It was realized that it may be very problematic for a low-mass companion to eject a massive envelope of the black hole progenitor. It invoked suggestions that an intermediate-mass companion ejects the envelope, and then is shredded by the Roche-lobe overflow to its current low-mass. But this creates another issue; a temperature mismatch between hot models and the observed cool low-mass donors. Finally, the main driver of Roche-lobe overflow that is believed to be magnetic braking does not seem to follow any theoretically calculated models. Number of ideas were put forward to explain various parts of this conundrum; pre-main sequence donor nature, alternative approach to magnetic braking and common envelope energy was revisited. We test various proposals and models to show that no overall solution exists so far. We argue that common envelope physics is not crucial in the understanding BH transient physical properties, but may affect significantly their formation rates. Our failure most likely indicates that either the current evolutionary models for low-mass stars and magnetic braking are not realistic or that the intrinsic population of BH transients is quite different from the observed one.

The observed Galactic rate of stellar mergers or the initiation of common envelope phases brighter than M_V=-3 (M_I=-4) is of order 0.5 (0.3)/year with 90% confidence statistical uncertainties of 0.24-1.1 (0.14-0.65) and factor of 2 systematic uncertainties. The (peak) luminosity function is roughly dN/dL L^(-1.4+/-0.3), so the rates for events more luminous than V1309 Sco (M_V=-7 mag) or V838Mon (M_V=-10 mag) are lower at r~0.1/year and 0.03/year, respectively. The peak luminosity is a steep function of progenitor mass, L M^(2-3). This very roughly parallels the scaling of luminosity with mass on the main sequence, but the transients are ~2000-4000 times more luminous at peak. Combining these, the mass function of the progenitors, dN/dM M^(-2.0+/-0.8), is consistent with the initial mass function, albeit with broad uncertainties. These observational results are also broadly consistent with the estimates of binary population synthesis models. While extragalactic variability surveys can better define the rates and properties of the high luminosity events, systematic, moderate depth (I>16 mag) surveys of the Galactic plane are needed to characterize the low luminosity events. The existing Galactic samples are only ~20% complete and Galactic surveys are (at best) reaching a typical magnitude limit of <13 mag.

The double-detonation explosion scenario of Type Ia supernovae has gained increased support from the SN Ia community as a viable progenitor model, making it a promising candidate alongside the well-known single degenerate and double degenerate scenarios. We present delay times of double-detonation SNe, in which a sub-Chandrasekhar mass carbon-oxygen white dwarf accretes non-dynamically from a helium-rich companion. One of the main uncertainties in quantifying SN rates from double-detonations is the (assumed) retention efficiency of He-rich matter. Therefore, we implement a new prescription for the treatment of accretion/accumulation of He-rich matter on white dwarfs. In addition, we test how the results change depending on which criteria are assumed to lead to a detonation in the helium shell. In comparing the results to our standard case (Ruiter et al. 2011), we find that regardless of the adopted He accretion prescription, the SN rates are reduced by only 25% if low-mass He shells (< 0.05 Msun) are sufficient to trigger the detonations. If more massive (0.1 Msun) shells are needed, the rates decrease by 85% and the delay time distribution is significantly changed in the new accretion model - only SNe with prompt (< 500 Myr) delay times are produced. Since theoretical arguments favour low-mass He shells for normal double-detonation SNe, we conclude that the rates from double-detonations are likely to be high, and should not critically depend on the adopted prescription for accretion of He.

We compare the current LIGO/Virgo upper limits on double compact object volumetric merger rates with our theoretical predictions. Our optimistic models are a factor of 3 below the existing upper limits for massive BH-BH systems with total mass 50-70 Msun, suggesting that a small increase in observational sensitivity may bring the first detections. The LIGO/Virgo gravitational wave observatories are currently being upgraded to advanced design sensitivity. If a sizeable population of BH-BH binaries is detected, the maximum total binary mass of this population will discriminate between two general families of common envelope models. If no binaries are detected, the new upper limits will provide astrophysically useful information about the environment and physical processes (e.g., metallicity of host galaxies or BH natal kicks) crucial to the formation of binaries containing black holes. For NS-NS systems, our predicted rates are 3 orders of magnitude below the current upper limits; even if advanced instruments reach their design sensitivities (factors of 10 times in distance, and 1,000 times in volumetric rate) the detection of NS-NS systems is not assured. However, we note that although our predicted NS-NS merger rates are consistent with estimates derived from Galactic NS-NS binaries and short GRBs, they are on the low side of these empirical estimates.

The development of advanced gravitational wave (GW) observatories, such as Advanced LIGO and Advanced Virgo, provides impetus to refine theoretical predictions for what these instruments might detect. In particular, with the range increasing by an order of magnitude, the search for GW sources is extending beyond the "local" Universe and out to cosmological distances. Double compact objects (neutron star-neutron star (NS-NS), black hole-neutron star (BH-NS) and black hole-black hole (BH-BH) systems) are considered to be the most promising gravitational wave sources. In addition, NS-NS and/or BH-NS systems are thought to be the progenitors of gamma ray bursts (GRBs), and may also be associated with kilonovae. In this paper we present the merger event rates of these objects as a function of cosmological redshift. We provide the results for four cases, each one investigating a different important evolution parameter of binary stars. Each case is also presented for two metallicity evolution scenarios. We find that (i) in most cases NS-NS systems dominate the merger rates in the local Universe, while BH-BH mergers dominate at high redshift; (ii) BH-NS mergers are less frequent than other sources per unit volume, for all time; and (iii) natal kicks may alter the observable properties of populations in a significant way, allowing the underlying models of binary evolution and compact object formation to be easily distinguished. This is the second paper in a series of three. The third paper will focus on calculating the detection rates of mergers by gravitational wave telescopes.

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

Emission from X-ray binaries (XRBs) is a major component of the total X-ray luminosity of normal galaxies, so X-ray studies of high redshift galaxies allow us to probe the formation and evolution of X-ray binaries on very long timescales. In this paper, we present results from large-scale population synthesis models of binary populations in galaxies from z = 0 to 20. We use as input into our modeling the Millennium II Cosmological Simulation and the updated semi-analytic galaxy catalog by Guo et al. (2011) to self-consistently account for the star formation history (SFH) and metallicity evolution of each galaxy. We run a grid of 192 models, varying all the parameters known from previous studies to affect the evolution of XRBs. We use our models and observationally derived prescriptions for hot gas emission to create theoretical galaxy X-ray luminosity functions (XLFs) for several redshift bins. Models with low CE efficiencies, a 50% twins mass ratio distribution, a steeper IMF exponent, and high stellar wind mass loss rates best match observational results from Tzanavaris & Georgantopoulos (2008), though they significantly underproduce bright early-type and very bright (Lx > 10d41) late-type galaxies. These discrepancies are likely caused by uncertainties in hot gas emission and SFHs, AGN contamination, and a lack of dynamically formed Low-mass XRBs. In our highest likelihood models, we find that hot gas emission dominates the emission for most bright galaxies. We also find that the evolution of the normal galaxy X-ray luminosity density out to z = 4 is driven largely by XRBs in galaxies with X-ray luminosities between 10d40 and 10d41 erg/s.

High redshift galaxies permit the study of the formation and evolution of X-ray binary populations on cosmological timescales, probing a wide range of metallicities and star-formation rates. In this paper, we present results from a large scale population synthesis study that models the X-ray binary populations from the first galaxies of the universe until today. We use as input to our modeling the Millennium II Cosmological Simulation and the updated semi-analytic galaxy catalog by Guo et al. (2011) to self-consistently account for the star formation history and metallicity evolution of the universe. Our modeling, which is constrained by the observed X-ray properties of local galaxies, gives predictions about the global scaling of emission from X-ray binary populations with properties such as star-formation rate and stellar mass, and the evolution of these relations with redshift. Our simulations show that the X-ray luminosity density (X-ray luminosity per unit volume) from X-ray binaries in our Universe today is dominated by low-mass X-ray binaries, and it is only at z>2.5 that high-mass X-ray binaries become dominant. We also find that there is a delay of ~1.1 Gyr between the peak of X-ray emissivity from low-mass Xray binaries (at z~2.1) and the peak of star-formation rate density (at z~3.1). The peak of the X-ray luminosity from high-mass X-ray binaries (at z~3.9), happens ~0.8 Gyr before the peak of the star-formation rate density, which is due to the metallicity evolution of the Universe.

We investigate the brightness distribution expected for thermonuclear explosions that might result from the ignition of a detonation during the violent merger of white dwarf (WD) binaries. Determining their brightness distribution is critical for evaluating whether such an explosion model could be responsible for a significant fraction of the observed population of SNe Ia. We argue that the brightness of an explosion realized via the violent merger model is determined by the mass of nickel 56 produced in the detonation of the primary WD. We use a set of sub-Chandrasekhar mass WD detonation models to derive a relationship between primary WD mass and expected peak bolometric brightness. We use this relationship to convert the masses of merging primary WDs from binary population models to a predicted distribution of explosion brightness. We find a striking similarity between the shape of our theoretical peak-magnitude distribution and that observed for SNe Ia: our model produces a brightness distribution that roughly covers the range and matches the shape of the one observed for SNe Ia. However, this agreement hinges on a particular phase of mass accretion during binary evolution: the primary WD gains ~0.15-0.35 Msun from a slightly-evolved helium star companion. We also find that with high probability, violent WD mergers involving the most massive primaries (which should produce bright SNe) have delay times <500 Myr (Abridged).

Cyg X-3 is a highly interesting accreting X-ray binary, emitting from the radio to high-energy gamma-rays. It consists of a compact object wind-fed by a Wolf-Rayet (WR) star, but the masses of the components and the mass-loss rate have been a subject of controversies. Here, we determine its masses, inclination, and the mass-loss rate using our derived relationship between the mass-loss rate and the mass for WR stars of the WN type, published infrared and X-ray data, and a relation between the mass-loss rate and the binary period derivative (observed to be >0 in Cyg X-3). Our obtained mass-loss rate is almost identical to that from two independent estimates and consistent with other ones, which strongly supports the validity of this solution. The found WR and compact object masses are 10.3_{-2.8}^{+3.9}, 2.4_{-1.1}^{+2.1} solar masses, respectively. Thus, our solution still allows for the presence of either a neutron star or a black hole, but the latter only with a low mass. However, the radio, infrared and X-ray properties of the system suggest that the compact object is a black hole. Such a low-mass black-hole could be formed via accretion-induced collapse or directly from a supernova.

The helium-merger gamma-ray burst progenitor is produced by the rapid accretion onto a compact remnant (neutron star or black hole) when it undergoes a common envelope inspiral with its companion's helium core. This merger phase produces a very distinct environment around these outbursts and recent observations suggest that, in some cases, we are detecting the signatures of the past merger in the GRB afterglow. These observations allow us, for the first time, to study the specific features of the helium merger progenitor. In this paper, we couple population synthesis calculations to our current understanding of gamma-ray burst engines and common envelope evolution to make observational predictions for the helium-merger gamma-ray burst population. Many mergers do not produce GRB outbursts and we discuss the implications of these mergers with the broader population of astrophysical transients.

We present results from deep X-ray stacking of >4000 high redshift galaxies from z~1 to 8 using the 4 Ms Chandra Deep Field South (CDF-S) data, the deepest X-ray survey of the extragalactic sky to date. The galaxy samples were selected using the Lyman break technique based primarily on recent HST ACS and WFC3 observations. Based on such high specific star formation rates (sSFRs): log SFR/M* > -8.7, we expect that the observed properties of these LBGs are dominated by young stellar populations. The X-ray emission in LBGs, eliminating individually detected X-ray sources (potential AGN), is expected to be powered by X-ray binaries and hot gas. We find, for the first time, evidence of evolution in the X-ray/SFR relation. Based on X-ray stacking analyses for z<4 LBGs (covering ~90% of the Universe's history), we find that the 2-10 keV X-ray luminosity evolves weakly with redshift (z) and SFR as log LX = 0.93 log (1+z) + 0.65 log SFR + 39.80. By comparing our observations with sophisticated X-ray binary population synthesis models, we interpret that the redshift evolution of LX/SFR is driven by metallicity evolution in HMXBs, likely the dominant population in these high sSFR galaxies. We also compare these models with our observations of X-ray luminosity density (total 2-10 keV luminosity per Mpc^3) and find excellent agreement. While there are no significant stacked detections at z>5, we use our upper limits from 5<z<8 LBGs to constrain the SMBH accretion history of the Universe around the epoch of reionization.

There are no known double black hole (BH-BH) or black hole-neutron star (BH-NS) systems. We argue that Cyg X-3 is a very likely BH-BH or BH-NS progenitor. This Galactic X-ray binary consists of a compact object, wind-fed by a Wolf-Rayet (WR) type companion. Based on a comprehensive analysis of observational data, it was recently argued that Cyg X-3 harbors a 2-4.5 Msun BH and a 7.5-14.2 Msun WR companion. We find that the fate of such a binary leads to the prompt (<1 Myr) formation of a close BH-BH system for the high end of the allowed WR mass (M_WR>13 Msun). For the low- to mid-mass range of the WR star (M_WR=7-10 Msun) Cyg X-3 is most likely (probability 70%) disrupted when WR ends up as a supernova. However, with smaller probability, it may form a wide (15%) or a close (15%) BH-NS system. The advanced LIGO/VIRGO detection rate for mergers of BH-BH systems from the Cyg X-3 formation channel is 10 per year, while it drops down to 0.1 per year for BH-NS systems. If Cyg X-3 in fact hosts a low mass BH and massive WR star, it lends additional support for the existence of BH-BH/BH-NS systems.

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 \msun$), 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. {\em (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.

Wide binaries made up of two white dwarfs (WDs) receive far less attention than their tight counterparts. However, our tests using the binary population synthesis code {\tt StarTrack} indicate that, for any set of reasonable initial conditions, there exists a significant observable population of double white dwarfs (WDWDs) with orbital separations of 10$^2$ to 10$^5$ AU. We adapt the technique of Dhital et al.\ to search for candidate common proper motion WD companions separated by $<10\amin$ around the $>$12,000 spectroscopically confirmed hydrogen-atmosphere WDs recently identified in the Sloan Digital Sky Survey. Using two techniques to separate random alignments from high-confidence pairs, we find nine new high-probability wide WDWDs and confirm three previously identified candidate wide WDWDs. This brings the number of known wide WDWDs to 45; our new pairs are a significant addition to the sample, especially at small proper motions ($<$200 mas/yr) and large angular separations ($>$10\asec). Spectroscopic follow-up and an extension of this method to a larger, photometrically selected set of SDSS WDs may eventually produce a large enough dataset for WDWDs to realize their full potential as testbeds for theories of stellar evolution.

BH-NS and BH-BH systems are among the most promising gravitational wave sources detectable by advanced LIGO/VIRGO and the Einstein Telescope. Although the rates of these systems may be above those of NS-NS mergers, BH-NS and BH-BH systems are difficult to detect, and thusfar none have been observed. But the progenitors of BH-NS and BH-BH binary systems may have been observed, in the form of High-Mass X-ray Binaries (HMXBs). In this paper, we continue work studying these potential progenitors of these important gravitational wave sources. In the first two papers of the series, we have demonstrated that IC10 X-1 and NGC300 X-1 are direct progenitors of BH-BH systems and that Cyg X-1 may form, alas with a very low probability, a BH-NS system. Here, we analyze the Galactic binaries GX 301-2, Vela X-1, XTEJ1855-026, 4U1907+09, Cir X-1 and extra-galactic LMC X-1, LMC X-3, M33 X-1. In each case, we find that the future evolution will not allow the formation of a BH-NS nor a BH-BH system. Most of these binaries will soon merge in the common envelope phase, with a compact object sinking into a helium-rich core of a stellar companion. This "helium-merger" may be a progenitor for long duration gamma-ray bursts (GRBs). Based on the observed HMXB population, the rate of helium-mergers may make up a sizable fraction of long-duration GRBs. Due to this high number of potential GRB progenitors, a chance that a Galactic HMXB has caused one of the recent major mass extinction events is significant (10-20%).

It is firmly established that the stellar mass distribution is smooth, covering the range 0.1-100 Msun. It is to be expected that the masses of the ensuing compact remnants correlate with the masses of their progenitor stars, and thus it is generally thought that the remnant masses should be smoothly distributed from the lightest white dwarfs to the heaviest black holes. However, this intuitive prediction is not borne out by observed data. In the rapidly growing population of remnants with observationally determined masses, a striking mass gap has emerged at the boundary between neutron stars and black holes. The heaviest neutron stars reach a maximum of two solar masses, while the lightest black holes are at least five solar masses. Over a decade after the discovery, the gap has become a significant challenge to our understanding of compact object formation. We offer new insights into the physical processes that bifurcate the formation of remnants into lower mass neutron stars and heavier black holes. Combining the results of stellar modeling with hydrodynamic simulations of supernovae, we both explain the existence of the gap, and also put stringent constraints on the inner workings of the supernova explosion mechanism. In particular, we show that core-collapse supernovae are launched within 100-200 milliseconds of the initial stellar collapse, implying that the explosions are driven by instabilities with a rapid (10-20 ms) growth time. Alternatively, if future observations fill in the gap, this will be an indication that these instabilities develop over a longer (>200 milliseconds) timescale.

Within the next decade, ground based gravitational wave detectors are in principle capable of determining the compact object merger rate per unit volume of the local universe to better than 20% with more than 30 detections. We argue that the stellar models are sensitive to heterogeneities (in age and metallicity at least) in such a way that the predicted merger rates are subject to an additional 30-50% systematic errors unless these heterogeneities are taken into account. Without adding new electromagnetic constraints on massive binary evolution or relying on more information from each merger (e.g., binary masses and spins), as few as the $\simeq 5$ merger detections could exhaust the information available in a naive comparison to merger rate predictions. As a concrete example, we use a nearby-universe catalog to demonstrate that no one tracer of stellar content can constrain merger rates without introducing a systematic error of order $O(30%)$ at 90% confidence. More generally, we argue that theoretical binary evolution can depend sufficiently sensitively on star-forming conditions -- even assuming no uncertainty in binary evolution model -- that the \emph{distribution} of star forming conditions must be incorporated to reduce the systematic error in merger rate predictions below roughly 40%. (Abridged)

The current star formation models imply that the binary fraction of population III stars is non zero. The evolution of such binaries must have led to formation of compact object binaries. In this paper we estimate the gravitational wave background originating in such binaries and discuss its observability. The properties of the population III binaries are investigated using a binary population synthesis code. We numerically model the background and we take into account the evolution of eccentric binaries. The gravitational wave background from population III binaries dominates the spectrum below 100 Hz. If the binary fraction is larger than 0.01 the background will be detectable by LISA and DECIGO. Gravitational wave background from population III binaries will dominate the spectrum below 100 Hz. LISA, ET and DECIGO should either see it easily or, in case of non detection, provide very strong constraints on the properties of the population III stars.

We present Chandra studies of the X-ray binary (XRB) populations in the bulge and ring regions of the ring galaxy NGC 1291. We detect 169 X-ray point sources in the galaxy, 75 in the bulge and 71 in the ring, utilizing the four available Chandra observations totaling an effective exposure of 179 ks. We report photometric properties of these sources in a point-source catalog. There are ~40% of the bulge sources and ~25% of the ring sources showing >3\sigma long-term variability in their X-ray count rate. The X-ray colors suggest that a significant fraction of the bulge (~75%) and ring (~65%) sources are likely low-mass X-ray binaries (LMXBs). The spectra of the nuclear source indicate that it is a low-luminosity AGN with moderate obscuration; spectral variability is observed between individual observations. We construct 0.3-8.0 keV X-ray luminosity functions (XLFs) for the bulge and ring XRB populations, taking into account the detection incompleteness and background AGN contamination. We reach 90% completeness limits of ~1.5\times10^{37} and ~2.2\times10^{37} erg/s for the bulge and ring populations, respectively. Both XLFs can be fit with a broken power-law model, and the shapes are consistent with those expected for populations dominated by LMXBs. We perform detailed population synthesis modeling of the XRB populations in NGC 1291, which suggests that the observed combined XLF is dominated by an old LMXB population. We compare the bulge and ring XRB populations, and argue that the ring XRBs are associated with a younger stellar population than the bulge sources, based on the relative overdensity of X-ray sources in the ring, the generally harder X-ray color of the ring sources, the overabundance of luminous sources in the combined XLF, and the flatter shape of the ring XLF.

In the Galaxy and in the Magellanic Clouds there are 170 Be X-ray binaries known to-date. Out of those, 111 host a neutron star, and for the reminder the nature of a companion is not known. None, so far, is known to host a black hole. This disparity is referred to as a missing Be -- black hole X-ray binary problem. The stellar population synthesis calculations following the formation of Be X-ray binaries in the Galaxy (Belczynski and Ziolkowski 2009) demonstrate that there is no problem of the missing Be+BH X-ray binaries for the Galaxy (the expected number of Be -- black hole X-ray binaries is 0 to 2, which is entirely consistent with the observed Galactic sample). However, the preliminary calculations for Magellanic Clouds indicate that there is a problem of the missing Be+BH X-ray binaries for the Clouds (the expected number is about 6, while none is observed). We believe, that to remove the discrepancy, one has to take into account a different history of the star formation rate in the Magellanic Clouds, with the respect to the Galaxy. New stellar population synthesis calculations are currently being carried out. An updated (as of November 2011) list of all 170 Be X-ray binaries known presently in the Galaxy and in the Magellanic Clouds is included.

The recent distance determination allowed precise estimation of the orbital parameters of Cyg X-1, which contains a massive 14.8 Msun BH with a 19.2 Msun O star companion. This system appears to be the clearest example of a potential progenitor of a BH-NS system. We follow the future evolution of Cyg X-1, and show that it will soon encounter a Roche lobe overflow episode, followed shortly by a Type Ib/c supernova and the formation of a NS. It is demonstrated that in majority of cases (70%) the supernova and associated natal kick disrupts the binary due to the fact that the orbit expanded significantly in the Roche lobe overflow episode. In the reminder of cases (30%) the newly formed BH-NS system is too wide to coalesce in the Hubble time. Only sporadically (1%) a Cyg X-1 like binary may form a coalescing BH-NS system given a favorable direction and magnitude of the natal kick. If Cyg X-1 like channel (comparable mass BH-O star bright X-ray binary) is the only or dominant way to form BH-NS binaries in the Galaxy we can estimate the empirical BH-NS merger rate in the Galaxy at the level of 0.001 per Myr. This rate is so low that the detection of BH-NS systems in gravitational radiation is highly unlikely, generating Advanced LIGO/VIRGO detection rates at the level of only 1 per century. If BH-NS inspirals are in fact detected, it will indicate that the formation of these systems proceeds via some alternative and yet unobserved channels.

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 solar masses 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.

Currently three isolated radio pulsars and one binary radio pulsar with no evidence of any previous recycling are known in 97 surveyed Galactic globular clusters. As pointed out by Lyne et al., the presence of these pulsars cannot be explained by core-collapse supernovae, as is commonly assumed for their counterparts in the Galactic disk. We apply a Bayesian analysis to the results from surveys for radio pulsars in globular clusters and find the number of potentially observable non-recycled radio pulsars present in all clusters to be < 3600. Accounting for beaming and retention considerations, the implied birth rate for any formation scenario for all 97 clusters is < 0.25 pulsars per century assuming a Maxwellian distribution of velocities with a dispersion of 10 km s^{-1}. The implied birth rates for higher velocity dispersions are substantially higher than inferred for such pulsars in the Galactic disk. This suggests that the velocity dispersion of young pulsars in globular clusters is significantly lower than those of disk pulsars. These numbers may be substantial overestimates due to the fact that the currently known sample of young pulsars is observed only in metal-rich clusters. We propose that young pulsars may only be formed in globular clusters with metallicities with log[Fe/H] > -0.6. In this case, the potentially observable population of such young pulsars is 447^{+1420}_{-399} (the error bars give the 95% confidence interval) and their birth rate is 0.012^{+0.037}_{-0.010} pulsars per century. The mostly likely creation scenario to explain these pulsars is the electron capture supernova of a OMgNe white dwarf.

We present theoretical delay times and rates of thermonuclear explosions that are thought to produce Type Ia supernovae, including the double-detonation sub-Chandrasekhar mass model, using the population synthesis binary evolution code StarTrack. If detonations of sub-Chandrasekhar mass carbon-oxygen white dwarfs following a detonation in an accumulated layer of helium on the white dwarf's surface ("double-detonation" models) are able to produce thermonuclear explosions which are characteristically similar to those of SNe Ia, then these sub-Chandrasekhar mass explosions may account for at least some substantial fraction of the observed SN Ia rate. Regardless of whether all double-detonations look like 'normal' SNe Ia, in any case the explosions are expected to be bright and thus potentially detectable. Additionally, we find that the delay time distribution of double-detonation sub-Chandrasekhar mass SNe Ia can be divided into two distinct formation channels: the 'prompt' helium-star channel with delay times <500 Myr (~10% of all sub-Chandras), and the 'delayed' double white dwarf channel, with delay times >800 Myr spanning up to a Hubble time (~90%). These findings coincide with recent observationally-derived delay time distributions which have revealed that a large number of SNe Ia are prompt with delay times <500 Myr, while a significant fraction also have delay times spanning ~1 Gyr to a Hubble time.

We present the summary of the recent investigations of double black hole binaries in context of their formation and merger rates. In particular we discuss the spectrum of black hole masses, the formation scenarios in the local Universe and the estimates of detection rates for gravitational radiation detectors like LIGO and VIRGO. Our study is based on observed properties of known Galactic and extra-galactic stellar mass black holes and evolutionary predictions. We argue that the binary black holes are the most promising source of gravitational radiation.

Using the StarTrack binary population synthesis code we model the population of double neutron stars in the Galaxy. We include a detailed treatment of the spin evolution of each pulsar due to processes such as spin-down and spin-up during accretion events as well as magnetic field decay. We also model the spatial distribution of double neutron stars by including their natal kicks and subsequent propagation in the Galactic gravitational potential. This synthetic pulsar population is compared to the observed sample of double neutron stars taking into account the selection effects of detection in the radio band, to determine the most likely evolutionary parameters. With these parameters we determine the properties of the double neutron star binaries detectable in gravitational waves by the high frequency interferometers LIGO and VIRGO. In particular, we discuss the distributions of chirp masses and mass ratios in samples selected by their radio or gravitational wave emission.

The current gravitational wave detectors have reached their operational sensitivity and are nearing detection of compact object binaries. In the coming years, we expect that the Advanced LIGO/VIRGO will start taking data. At the same time, there are plans for third generation ground-based detectors such as the Einstein Telescope, and space detectors such as DECIGO. We discuss the eccentricity distribution of inspiral compact object binaries during they inspiral phase. We analyze the expected distributions of eccentricities at three frequencies that are characteristic of three future detectors: Advanced LIGO/VIRGO (30 Hz), Einstein Telescope (3 Hz), and DECIGO (0.3 Hz). We use the StarTrack binary population code to investigate the properties of the population of compact binaries in formation. We evolve their orbits until the point that they enter a given detector sensitivity window and analyze the eccentricity distribution at that time. We find that the eccentricities of BH-BH and BH-NS binaries are quite small when entering the Advanced LIGO/VIRGO detector window for all considered models of binary evolution. Even in the case of the DECIGO detector, the typical eccentricities of BH-BH binaries are below 10^{-4}, and the BH-NS eccentricities are smaller than 10^{-3}. Some fraction of NS-NS binaries may have significant eccentricities. Within the range of considered models, we found that a fraction of between 0.2% and 2% NS-NS binaries will have an eccentricity above 0.01 for the Advanced LIGO/VIRGO detectors. For the ET detector, this fraction is between 0.4% and 4%, and for the DECIGO detector it lies between 2% and 27%.

We investigate the future evolution of two extragalactic X-ray binaries: IC10 X-1 and NGC300 X-1. Each of them consists of a high mass BH ($\sim 20-30 \msun$) accreting from a massive WR star companion ($\gtrsim 20 \msun$), and both are located in low metallicity galaxies. We analyze the current state of the systems and demonstrate that both systems will very quickly ($\lesssim 0.3$ Myr) form close BH-BH binaries with the short coalescence time ($\sim 3$ Gyr) and large chirp mass ($\sim 15 \msun$). The formation of BH-BH system seems unavoidable, as {\em (i)} WR companions are well within their Roche lobes and they do not expand so no Roche lobe overflow is expected, {\em (ii)} even intense WR wind mass loss does not remove sufficient mass to prohibit the formation of the second BH, {\em (ii)} even if BH receives the large natal kick, the systems are very closely bound and are almost impossible to disrupt. As there are two such immediate BH-BH progenitor systems within 2 Mpc and as the current gravitational wave instruments LIGO/VIRGO (initial stage) can detect such massive BH-BH mergers out to $\sim 200$ Mpc, the empirically estimated detection rate of such inspirals is $R=3.36^{+8.29}_{-2.92}$ at the 99% confidence level. If there is no detection in the current LIGO/VIRGO data (unreleased year of $s6$ run), the existence of these two massive BH systems poses an interesting challenge. Either the gravitational radiation search is not sensitive to massive inspirals or there is some fundamental misunderstanding of stellar evolution physics leading directly to the formation of BH-BH binaries.

We study the Galactic field population of double compact objects (NS-NS, BH-NS, BH-BH binaries) to investigate the number (if any) of these systems that can potentially be detected with LISA at low gravitational-wave frequencies. We calculate the Galactic numbers and physical properties of these binaries and show their relative contribution from the disk, bulge and halo. Although the Galaxy hosts 10^5 double compact object binaries emitting low-frequency gravitational waves, only a handful of these objects in the disk will be detectable with LISA, but none from the halo or bulge. This is because the bulk of these binaries are NS-NS systems with high eccentricities and long orbital periods (weeks/months) causing inefficient signal accumulation (small number of signal bursts at periastron passage in 1 yr of LISA observations) rendering them undetectable in the majority of these cases. We adopt two evolutionary models that differ in their treatment of the common envelope phase that is a major (and still mostly unknown) process in the formation of close double compact objects. Depending on the adopted evolutionary model, our calculations indicate the likely detection of about 4 NS-NS binaries and 2 BH-BH systems (model A; likely survival of progenitors through CE) or only a couple of NS-NS binaries (model B; suppression of the double compact object formation due to CE mergers).

Using the population synthesis binary evolution code StarTrack, we present theoretical rates and delay times of Type Ia supernovae arising from various formation channels. These channels include binaries in which the exploding white dwarf reaches the Chandrasekhar mass limit (DDS, SDS, and helium-rich donor scenario) as well as the sub-Chandrasekhar mass scenario, in which a white dwarf accretes from a helium-rich companion and explodes as a SN Ia before reaching the Chandrasekhar mass limit. We find that using a common envelope parameterization employing energy balance with alpha=1 and lambda=1, the supernova rates per unit mass (born in stars) of sub-Chandrasekhar mass SNe Ia exceed those of all other progenitor channels at epochs t=0.7 - 4 Gyr for a burst of star formation at t=0. Additionally, the delay time distribution of the sub-Chandrasekhar model can be divided in to two distinct evolutionary channels: the `prompt' helium-star channel with delay times < 500 Myr, and the `delayed' double white dwarf channel with delay times > 800 Myr spanning up to a Hubble time. These findings are in agreement with recent observationally-derived delay time distributions which predict that a large number of SNe Ia have delay times < 1 Gyr, with a significant fraction having delay times < 500 Myr. We find that the DDS channel is also able to account for the observed rates of SNe Ia. However, detailed simulations of white dwarf mergers have shown that most of these mergers will not lead to SNe Ia but rather to the formation of a neutron star via accretion-induced collapse. If this is true, our standard population synthesis model predicts that the only progenitor channel which can account for the rates of SNe Ia is the sub-Chandrasekhar mass scenario, and none of the other progenitors considered can fully account for the observed rates.

The merger of two white dwarfs (a.k.a. double degenerate merger) has often been cited as a potential progenitor of type Ia supernovae. Here we combine population synthesis, merger and explosion models with radiation-hydrodynamics light-curve models to study the implications of such a progenitor scenario on the observed type Ia supernova population. Our standard model, assuming double degenerate mergers do produce thermonuclear explosions, produces supernova light-curves that are broader than the observed type Ia sample. In addition, we discuss how the shock breakout and spectral features of these double degenerate progenitors will differ from the canonical bare Chandrasekhar-massed explosion models. We conclude with a discussion of how one might reconcile these differences with current observations.

Double white dwarfs are expected to be a source of confusion-limited noise for the future gravitational wave observatory LISA. In a specific frequency range, this 'foreground noise' is predicted to rise above the instrumental noise and hinder the detection of other types of signals, e.g., gravitational waves arising from stellar mass objects inspiraling into massive black holes. In many previous studies only detached populations of compact object binaries have been considered in estimating the LISA gravitational wave foreground signal. Here, we investigate the influence of compact object detached and Roche-Lobe Overflow Galactic binaries on the shape and strength of the LISA signal. Since >99% of remnant binaries which have orbital periods within the LISA sensitivity range are white dwarf binaries, we consider only these binaries when calculating the LISA signal. We find that the contribution of RLOF binaries to the foreground noise is negligible at low frequencies, but becomes significant at higher frequencies, pushing the frequency at which the foreground noise drops below the instrumental noise to >6 mHz. We find that it is important to consider the population of mass transferring binaries in order to obtain an accurate assessment of the foreground noise on the LISA data stream. However, we estimate that there still exists a sizeable number (~11300) of Galactic double white dwarf binaries which will have a signal-to-noise ratio >5, and thus will be potentially resolvable with LISA. We present the LISA gravitational wave signal from the Galactic population of white dwarf binaries, show the most important formation channels contributing to the LISA disc and bulge populations and discuss the implications of these new findings.

We investigate the statistics of isolated recycled pulsars and double neutron star binaries in the Galactic disk. Since recycled pulsars are believed to form through accretion and spinup in close binaries, the isolated objects presumably originate from disrupted progenitors of double neutron stars. There are a comparable number of double neutron star systems compared to isolated recycled pulsars. We find that standard evolutionary models cannot explain this fact, predicting several times the number of isolated recycled pulsars than those in double neutron star systems. We demonstrate, through population synthesis calculations, that the velocity distribution of isolated recycled pulsars is broader than for binary systems. When this is accounted for in a model for radio pulsar survey selection effects, which include the effects of Doppler smearing for the double neutron star binaries, we find that there is a small (25%) bias towards the detection of double neutron star systems. This bias, however, is not significant enough to explain the observational discrepancy if standard (sigma = 265 km/s) neutron star natal kick velocities are invoked in binary population syntheses. Population syntheses in which the 1D Maxwellian velocity dispersion of the natal kick is sigma=170 km/s are consistent with the observations. These conclusions further support earlier findings the neutron stars formed in close interacting binaries receive significantly smaller natal kicks than the velocities of Galactic single pulsars would seem to indicate.

Data from the SDSS (300,000 galaxies) indicates that recent star formation (within the last 1 billion years) is bimodal: half the stars form from gas with high amounts of metals (solar metallicity), and the other half form with small contribution of elements heavier than Helium (10-30% solar). Theoretical studies of mass loss from the brightest stars derive significantly higher stellar-origin BH masses (30-80 Msun) than previously estimated for sub-solar compositions. We combine these findings to estimate the probability of detecting gravitational waves (GWs) arising from the inspiral of double compact objects. Our results show that a low metallicity environment significantly boosts the formation of double compact object binaries with at least one BH. In particular, we find the GW detection rate is increased by a factor of 20 if the metallicity is decreased from solar (as in all previous estimates) to a 50-50 mixture of solar and 10% solar metallicity. The current sensitivity of the two largest instruments to NS-NS binary inspirals (VIRGO: 9 Mpc; LIGO: 18) is not high enough to ensure a first detection. However, our results indicate that if a future instrument increased the sensitivity to 50-100 Mpc, a detection of GWs would be expected within the first year of observation. It was previously thought that NS-NS inspirals were the most likely source for GW detection. Our results indicate that BH-BH binaries are 25-times more likely sources than NS-NS systems and that we are on the cusp of GW detection.

We estimate binary compact object merger detection rates for LIGO, including the binaries formed in ellipticals long ago. Specifically, we convolve hundreds of model realizations of elliptical- and spiral-galaxy population syntheses with a model for elliptical- and spiral-galaxy star formation history as a function of redshift. Our results favor local merger rate densities of 4\times 10^{-3} {Mpc}^{-3}{Myr}^{-1} for binary black holes (BH), 3\times 10^{-2} {Mpc}^{-3}{Myr}^{-1} for binary neutron stars (NS), and 10^{-2} {Mpc}^{-3}{Myr}^{-1} for BH-NS binaries. Mergers in elliptical galaxies are a significant fraction of our total estimate for BH-BH and BH-NS detection rates; NS-NS detection rates are dominated by the contribution from spiral galaxies. Using only models that reproduce current observations of Galactic NS-NS binaries, we find slightly higher rates for NS-NS and largely similar ranges for BH-NS and BH-BH binaries. Assuming a detection signal-to-noise ratio threshold of 8 for a single detector (as part of a network), corresponding to radii \Cv of the effective volume inside of which a single LIGO detector could observe the inspiral of two 1.4 M_\sun neutron stars of 14 Mpc and 197 Mpc, for initial and advanced LIGO, we find event rates of any merger type of 2.9* 10^{-2} -- 0.46 and 25-400 per year (at 90% confidence level), respectively. We also find that the probability P_{detect} of detecting one or more mergers with this single detector can be approximated by (i) P_{detect}\simeq 0.4+0.5\log (T/0.01{yr}), assuming \Cv=197 {Mpc} and it operates for T years, for T between 2 days and 0.1 {yr}); or by (ii) P_{detect}\simeq 0.5 + 1.5 \log \Cv/32{Mpc}, for one year of operation and for $\Cv$ between 20 and 70 Mpc. [ABRIDGED]

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 per Myr per Milky Way Equivalent Galaxy (MWEG), although the rate could plausibly range from 1 per Myr per MWEG to 1000 per Myr per MWEG. 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 0.0002 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.

We present the spectrum of compact object masses: neutron stars and black holes that originate from single stars in different environments. In particular, we calculate the dependence of maximum black hole mass on metallicity and on some specific wind mass loss rates (e.g., Hurley et al. and Vink et al.). Our calculations show that the highest mass black holes observed in the Galaxy M_bh = 15 Msun in the high metallicity environment (Z=Zsun=0.02) can be explained with stellar models and the wind mass loss rates adopted here. To reach this result we had to set Luminous Blue Variable mass loss rates at the level of about 0.0001 Msun/yr and to employ metallicity dependent Wolf-Rayet winds. With such winds, calibrated on Galactic black hole mass measurements, the maximum black hole mass obtained for moderate metallicity (Z=0.3 Zsun=0.006) is M_bh,max = 30 Msun. This is a rather striking finding as the mass of the most massive known stellar black hole is M_bh = 23-34 Msun and, in fact, it is located in a small star forming galaxy with moderate metallicity. We find that in the very low (globular cluster-like) metallicity environment the maximum black hole mass can be as high as M_bh,max = 80 Msun (Z=0.01 Zsun=0.0002). It is interesting to note that X-ray luminosity from Eddington limited accretion onto an 80 Msun black hole is of the order of about 10^40 erg/s and is comparable to luminosities of some known ULXs. We emphasize that our results were obtained for single stars only and that binary interactions may alter these maximum black hole masses (e.g., accretion from a close companion). This is strictly a proof-of-principle study which demonstrates that stellar models can naturally explain even the most massive known stellar black holes.

In black hole X-ray binaries, a misalignment between the spin axis of the black hole and the orbital angular momentum can occur during the supernova explosion that forms the compact object. In this letter we present population synthesis models of Galactic black hole X-ray binaries, and study the probability density function of the misalignment angle, and its dependence on our model parameters. In our modeling, we also take into account the evolution of misalignment angle due to accretion of material onto the black hole during the X-ray binary phase. The major factor that sets the misalignment angle for X-ray binaries is the natal kick that the black hole may receive at its formation. However, large kicks tend to disrupt binaries, while small kicks allow the formation of XRBs and naturally select systems with small misalignment angles. Our calculations predict that the majority (>67%) of Galactic field BH XRBs have rather small (>10 degrees) misalignment angles, while some systems may reach misalignment angles as high as ~90 degrees and even higher. This results is robust among all population synthesis models. The assumption of small small misalignment angles is extensively used to observationally estimate black hole spin magnitudes, and for the first time we are able to confirm this assumption using detailed population synthesis calculations.

In the Galaxy there are 64 Be X-ray binaries known to-date. Out of those, 42 host a neutron star, and for the reminder the nature of a companion is not known. None, so far, is known to host a black hole. There seems to be no apparent mechanism that would prevent formation or detection of Be stars with black holes. This disparity is referred to as a missing Be -- black hole X-ray binary problem. We point out that current evolutionary scenarios that lead to the formation of Be X-ray binaries predict that the ratio of these binaries with neutron stars to the ones with black holes is rather high F_NStoBH=10-50, with the more likely formation models providing the values at the high end. The ratio is a natural outcome of (i) the stellar initial mass function that produces more neutron stars than black holes and (ii) common envelope evolution (i.e. a major mechanism involved in the formation of interacting binaries) that naturally selects progenitors of Be X-ray binaries with neutron stars (binaries with comparable mass components have more likely survival probabilities) over ones with black holes (which are much more likely to be common envelope mergers). A comparison of this ratio (i.e. F_NStoBH=30) with the number of confirmed Be -- neutron star X-ray binaries (42) indicates that the expected number of Be -- black hole X-ray binaries is of the order of only 0-2. This is entirely consistent with the observed Galactic sample.

GRB 080913 and GRB 090423 are the most distant gamma-ray bursts (GRBs) known to-date, with spectroscopically determined redshifts of z=6.7 and z=8.1, respectively. The detection of bursts at this early epoch of the Universe significantly constrains the nature of GRBs and their progenitors. We perform population synthesis studies of the formation and evolution of early stars, and calculate the resulting formation rates of short and long-duration GRBs at high redshift. The peak of the GRB rate from Population II stars occurs at z=7 for a model with efficient/fast mixing of metals, while it is found at z=3 for an inefficient/slow metallicity evolution model. We show that in the redshift range 6<z<10 essentially all GRBs originate from Population II stars, regardless of metallicity evolution model. These stars (having small, but non-zero metallicity) are the most likely progenitors for both long GRBs (collapsars) and short GRBs (NS-NS or BH-NS mergers) at this epoch. Although the predicted intrinsic rates of long and short GRBs are similar at these high redshifts, observational selection effects lead to higher (factor of 10) observed rates for long GRBs. We conclude that the two recently observed high-z GRB events are most likely long GRBs originating from Population II collapsars.

We review the formation scenarios for binary black holes, and show that their coalescence rate depends very strongly on the outcome of the second mass transfer. However, the observations of IC10 X-1, an binary with a massive black hole accreting from a Wolf-Rayet star proves that this mass transfer can be stable. We analyze the future evolution of IC10 X-1 and show that it is very likely to form a binary black hole system merging in a few Gyrs. We estimate the coalescence rate density of such systems to be $ 0.06 {\rm Mpc}^{-3} {\rm Myr}^{-1}$, and the detection rate for the current LIGO/VIRGO of $ 0.69 {\rm yr}^{-1} $, a much higher value than the expected double neutron star rate. Thus the first detection of a gravitational wave source is likely to be a coalescence of a binary black hole.

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 analyze the evolution of binary stars to calculate synthetic rates and delay times of the most promising Type Ia Supernovae progenitors. We present and discuss evolutionary scenarios in which a white dwarf reaches the Chandrasekhar-mass and potentially explodes in a Type Ia supernova. We consider: Double Degenerate (DDS), Single Degenerate (SDS), and AM Canum Venaticorum scenarios. The results are presented for two different star formation histories; burst (elliptical-like galaxies) and continuous (spiral-like galaxies). It is found that delay times for the DDS in our standard model (with common envelope efficiency alpha = 1) follow a power-law distribution. For the SDS we note a wide range of delay times, while AM CVn progenitors produce a short burst of SNe Ia at early times. We point out that only the rates for two merging carbon-oxygen white dwarfs, the only systems found in the DDS, are consistent with the observed rates for typical Milky Way-like spirals. We also note that DDS progenitors are the dominant population in elliptical galaxies. The fact that the delay time distribution for the DDS follows a power-law implies more Type Ia supernovae (per unit mass) in young rather than in aged populations. Our results do not exclude other scenarios, but strongly indicate that the DDS is the dominant channel generating SNe Ia in spiral galaxies, at least in the framework of our adopted evolutionary models. Since it is believed that white dwarf mergers cannot produce a thermonuclear explosion given the current understanding of accreting white dwarfs, either the evolutionary calculations along with accretion physics are incorrect, or the explosion calculations are inaccurate and need to be revisited (Abridged).

This white paper briefly describes the astrophysics of ultra-compact binaries, with emphasis of the challenges and opportunities in the next decade.

We propose a physically motivated and self-consistent prescription for the modeling of transient neutron star (NS) low-mass X-ray binary (LMXB) properties, such as duty cycle (DC), outburst duration and recurrence time. We apply this prescription to the population synthesis (PS) models of field LMXBs presented by Fragos et al. (2008), and compare the transient LMXB population to the Chandra X-ray survey of the two elliptical galaxies NGC 3379 and NGC 4278, which revealed several transient sources (Brassington et al., 2008, 2009). We are able to exclude models with a constant DC for all transient systems, while models with a variable DC based on the properties of each system are consistent with the observed transient populations. We predict that the majority of the observed transient sources in these two galaxies are LMXBs with red giant donors. Our comparison suggests that LMXBs formed through evolution of primordial field binaries are dominant in globular cluster (GC) poor elliptical galaxies, while they still have a significant contribution in GC rich ones.

Compact object binaries, mostly double with dwarfs, are believed to be a potential source of confusion-limited noise for the Laser Interferometer Space Antenna (LISA). In a specific frequency range, this noise may rise above the instrumental noise sources and hence hinder detection of other types of signals, e.g. extreme mass ratio inspirals of stellar mass objects into supermassive black holes. In most previous studies only detached populations of compact object binaries have been considered. Here, we investigate the influence of Galactic populations of compact mass transferring binaries on the shape and strength of the LISA signal and compare our results with the signal for detached binaries only. Our population synthesis includes all binary systems containing two compact remnants. It is found that 99.5% of these systems are double white dwarfs, and therefore we consider only binaries consisting of two white dwarfs when calculating the LISA signal. We find that the mass transferring binaries contributing to the confusion limited noise dominate the population of LISA systems by number (60%), due to their long lifetimes (few Gyrs). However, due to the very specific physical properties of mass transferring systems, it is found that their contribution to the noise is negligible at low frequencies and detached systems will set the magnitude of the confusion-limited signal in this range. However, the mass transferring systems begin to dominate the signal at higher frequencies (f>6 mHz). The transition frequency at which the confusion limited regime becomes important is found at 4 mHz. The addition of mass transferring systems tends to slightly increase the transition frequency (from 3.5 to 4 mHz).

Galactic double white dwarfs were postulated as a source of confusion limited noise for LISA, the future space-based gravitational wave observatory. Until very recently, the Galactic population consisted of a relatively well studied disk population, a somewhat studied smaller bulge population and a mostly unknown, but potentially large halo population. It has been argued that the halo population may produce a signal that is much stronger (factor of ~5 in spectral amplitude) than the disk population. However, this surprising result was not based on an actual calculation of a halo white dwarf population but was derived on (i) the assumption that one can extrapolate the halo population properties from those of the disk population and (ii) the postulated (unrealistically) high number of white dwarfs in the halo. We perform the first calculation of a halo white dwarf population using population synthesis models. Our comparison with the signal arising from double white dwarfs in the Galactic disk+bulge clearly shows that it is impossible for the double white dwarf halo signal to exceed that of the rest of the Galaxy. Using microlensing results to give an upper limit on the content of white dwarfs in the halo (~30 % baryonic mass in white dwarfs), our predicted halo signal is a factor of 10 lower than the disk+bulge signal. Even in the implausible case where all of the baryonic halo mass is found in white dwarfs, the halo signal does not become comparable to that of the disk+bulge, and thus would still have a negligible effect on the detection of other LISA sources.

The evolution leading to the formation of a neutron star in the very young Westerlund 1 star cluster is investigated. The turnoff mass has been estimated to be 35 Msun, indicating a cluster age ~ 3-5 Myr. The brightest X-ray source in the cluster, CXO J164710.2-455216, is a slowly spinning (10 s) single neutron star and potentially a magnetar. Since this source was argued to be a member of the cluster, the neutron star progenitor must have been very massive (M_zams > 40 Msun) as noted by Muno et al. (2006). Since such massive stars are generally believed to form black holes (rather than neutron stars), the existence of this object poses a challenge for understanding massive star evolution. We point out while single star progenitors below M_zams < 20 Msun form neutron stars, binary evolution completely changes the progenitor mass range. In particular, we demonstrate that mass loss in Roche lobe overflow enables stars as massive as 50-80 Msun, under favorable conditions, to form neutron stars. If the very high observed binary fraction of massive stars in Westerlund 1 (> 70 percent) is considered, it is natural that CXO J164710.2-455216 was formed in a binary which was disrupted in a supernova explosion such that it is now found as a single neutron star. Hence, the existence of a neutron star in a given stellar population does not necessarily place stringent constraints on progenitor mass when binary interactions are considered. It is concluded that the existence of a neutron star in Westerlund 1 cluster is fully consistent with the generally accepted framework of stellar evolution.

The evolution of the spin and tilt of black holes in compact black hole - neutron star and black hole - black hole binary systems is investigated within the framework of the coalescing compact star binary model for short gamma ray bursts via the population synthesis method. Based on recent results on accretion at super critical rates in slim disk models, estimates of natal kicks, and the results regarding fallback in supernova models, we obtain the black hole spin and misalignment. It is found that the spin parameter, a_spin}, is less than 0.5 for initially non rotating black holes and the tilt angle, i_tilt, is less than 45 deg for 50% of the systems in black hole - neutron star binaries. Upon comparison with the results of black hole - neutron star merger calculations we estimate that only a small fraction (~ 0.01) of these systems can lead to the formation of a torus surrounding the coalesced binary potentially producing a short-hard gamma ray burst. On the other hand, for high initial black hole spin parameters (a_spin>0.6) this fraction can be significant (~ 0.4). It is found that the predicted gravitational radiation signal for our simulated population does not significantly differ from that for non rotating black holes. Due to the (i) insensitivity of signal detection techniques to the black hole spin and the (ii) predicted overall low contribution of black hole binaries to the signal we find that the detection of gravitational waves are not greatly inhibited by current searches with non spinning templates. It is pointed out that the detection of a black hole - black hole binary inspiral system with LIGO or VIRGO may provide a direct measurement of the initial spin of a black hole.

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).

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)

We compare the X-ray spectra and luminosities, in the 2-8 keV band, of known and suspected cataclysmic variables (CVs) in different environments, assessing the nature of these source populations. These objects include nearby CVs observed with ASCA; the Galactic Center X-ray source population identified by Muno et al.; and likely CVs identified in globular clusters. Both of the latter have been suggested to be dominated by magnetic CVs. We find that the brighter objects in both categories are likely to be magnetic CVs, but that the fainter objects are likely to include a substantial contribution from normal CVs. The strangely hard spectra observed from the Galactic Center sources reflect the high and variable extinction, which is significantly greater than the canonical 6e22 /cm2 over much of the region, and the magnetic nature of many of the brightest CVs. The total numbers of faint Galactic Center sources are compatible with expectations of the numbers of CVs in this field.

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.

We present theoretical models for the formation and evolution of populations of low-mass X-ray binaries (LMXB) in the two elliptical galaxies NGC 3379 and NGC 4278. The models are calculated with the recently updated StarTrack code (Belczynski et al., 2006), assuming only a primordial galactic field LMXB population. StarTrack is an advanced population synthesis code that has been tested and calibrated using detailed binary star calculations and incorporates all the important physical processes of binary evolution. The simulations are targeted to modeling and understanding the origin of the X-ray luminosity functions (XLF) of point sources in these galaxies. For the first time we explore the population XLF down to luminosities of 3X10^36 erg/s, as probed by the most recent observational results (Kim et al., 2006). We consider models for the formation and evolution of LMXBs in galactic fields with different CE efficiencies, stellar wind prescriptions, magnetic braking laws and initial mass functions. We identify models that produce an XLF in excellent agreement with the observations both in shape and absolute normalization. We also find that the treatment of the outburst luminosity of transient systems remains a crucial factor for the determination of the XLF since the modeled populations are dominated by transient X-ray systems.

Using a population synthesis approach, we compute the total merger rate in the local Universe for double neutron stars, double black holes, and black hole -- neutron star binaries. These compact binaries are the prime source candidates for gravitational-wave detection by LIGO and VIRGO. We account for mergers originating from field populations and from dense stellar clusters. For both populations we use the same treatment of stellar evolution. Our results indicate that the merger rates of double neutron stars and black hole -- neutron star binaries are strongly dominated by field populations, while merging black hole binaries are formed much more effectively in dense stellar clusters. The overall merger rate of double compact objects depends sensitively on the (largely unknown) initial mass fraction contained in dense clusters f_cl. For f_cl < 0.0001, the Advanced LIGO detection rate will be dominated by field populations of double neutron star mergers, with a small but significant number of detections ~20 yr^-1. However for a higher mass fraction in clusters, f_cl > 0.001, the detection rate will be dominated by numerous mergers of double black holes originating from dense clusters, and it will be considerably higher, ~25 - 300 yr^-1. In addition, we show that, once mergers of double black holes are detected, it is easy to differentiate between systems formed in the field and in dense clusters, since the chirp mass distributions are strikingly different. Finally, we point out that there may exist a population of merging black hole binaries in intergalactic space.

We investigate the effect of including a significant ``binary twin'' population (binaries with almost equal mass stars, q = M2/M1 > 0.95) for the production of double compact objects and some resulting consequences, including LIGO inspiral rate and some properties of short-hard gamma-ray bursts. We employ very optimistic assumptions on the twin fraction (50%) among all binaries, and therefore our calculations place an upper limits on the influence of twins on double compact object populations. We show that for LIGO the effect of including twins is relatively minor: although the merger rates does indeed increase when twins are considered, the rate increase is fairly small (1.5). Also, chirp mass distribution for double compact objects formed with or without twins are almost indistinguishable. If double compact object are short-hard GRB progenitors, including twins in population synthesis calculations does not alter significantly the earlier rate predictions for the event rate. However, for one channel of binary evolution, introducing twins more than doubles the rate of ``very prompt'' NS-NS mergers (time to merger less than 1 Myr) compared to models with the ``flat'' q distribution. In that case, 70% of all NS-NS binaries merge within 100 Myr after their formation, indicating a possibility of a very significant population of ``prompt'' short-hard gamma-ray bursts, associated with star forming galaxies. We also point out that, independent of assumptions, fraction of such prompt neutron star mergers is always high, 35--70%. We note that recent observations (e.g., Berger et al.) indicate that fraction of short-hard GRBs found in young hosts is at least 40% and possibly even 80%.

To test whether the short GRB rates, redshift distribution and host galaxies are consistent with current theoretical predictions, we use avery large database of population synthesis calculations to examine BH-NS and NS-NS merger rates in the universe, factoring in (i) the star formation history of the universe, (ii) a heterogeneous population of star-forming galaxies, including spirals and ellipticals, and (iii) a simple flux-limited selection model for short GRB detection. When we require our models reproduce the known short GRB rates and redshift measurements (and, for NS-NS, the merger rates extrapolated from binary pulsars in the Galaxy), a small fraction of models reproduce all observations, both when we assume a NS-NS and a BH-NS origin for bursts. Most commonly models produce mergers preferentially in spiral galaxies if short GRBs arise from NS-NS mergers alone. Model universes where present-day binary mergers occur preferentially in elliptical galaxies necessarily include a significant fraction of binaries with long delay times between birth and merger (often $O(10{\rm Gyr})$). Though long delays occur, almost all of our models predict that a higher proportion of short GRBs should occur at moderate to high redshift (e.g., $z>1$) than has presently been observed, in agreement with recent observations which suggest a selection bias towards successful follow-up of low-redshift short GRBs. Finally, if only a fraction of BH-NS mergers have the right combination of masses and spins to make GRBs, then at best only a small fraction of BH-NS models could be consistent with all {\em current} available data. (Abridged)

We present results of a population synthesis study aimed at examining the role of spin-kick alignment in producing a correlation between the spin period of the first-born neutron star and the orbital eccentricity of observed double neutron star binaries in the Galactic disk. We find spin-kick alignment to be compatible with the observed correlation, but not to alleviate the requirements for low kick velocities suggested in previous population synthesis studies. Our results furthermore suggest low- and high-eccentricity systems may form through two distinct formation channels distinguished by the presence or absence of a stable mass transfer phase before the formation of the second neutron star. The presence of highly eccentric systems in the observed sample of double neutron stars may furthermore support the notion that neutron stars accrete matter when moving through the envelope of a giant companion.

The observed samples of supernovae (SN) and double compact objects (DCOs) provide several critical constraints on population-synthesis models: the parameters of these models must be carefully chosen to reproduce, among other factors, (i) the formation rates of double neutron star (NS-NS) binaries and of white dwarf-neutron star (WD-NS) binaries, estimated from binary samples, and (ii) the type II and Ib/c supernova rates. Even allowing for extremely conservative accounting of the uncertainties in observational and theoretical predictions, we find only a few plausible population synthesis models (roughly 9%) are consistent with DCO and SN rates empirically determined from observations. As a proof of concept, we describe the information that can be extracted about population synthesis models given such stringent observational tests, including surprisingly good agreement with the neutron star kick distributions inferred from pulsar proper-motion measurements. In the present study, we find that the current observational constraints favor: kicks described by a single Maxwellian with a typical velocity of about 300km/s; mass-loss fractions during non-conservative, but stable, mass transfer episodes of about 90%; and common envelope parameters of about 0.2-0.5. Finally, we use the subset of astrophysically consistent models to predict the rates at which black hole-neutron star (BH-NS) and NS-NS binaries merge in the Milky Way and the nearby Universe, assuming Milky-Way-like galaxies dominate. (Abridged)

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.

Binary population synthesis calculations and associated predictions, especially event rates, are known to depend on a significant number of input model parameters with different degrees of sensitivity. At the same time, for systems with relatively low formation rates, such simulations are heavily computationally demanding and therefore the needed explorations of the high-dimensional parameter space require major -- often prohibitive -- computational resources. In the present study, to better understand several key event rates involving binary evolution and binaries with two compact objects in Milky Way-like galaxies and to provide ways of reducing the computational costs of complete parameter space explorations: (i) we perform a methodical parameter study of the \emph{StarTrack} population synthesis code ; and (ii) we develop a formalism and methodology for the derivation of {\em multi-dimensional fits} for event rates. We significantly generalize our earlier study, and we focus on ways of thoroughly assessing the accuracy of the fits. We anticipate that the efficient tools developed here can be applied in lieu of large-scale population calculations and will facilitate the exploration of the dependence of rate predictions on a wide range binary evolution parameters. Such explorations can then allow the derivation of constraints on these parameters, given empirical rate constraints and accounting for fitting errors. Here we describe in detail the principles and practice behind constructing these fits, estimating their accuracy, and comparing them with observations in a manner that accounts for their errors.

Double black hole binaries are among the most important sources of gravitational radiation for ground-based 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), but are still quite high for double neutron stars (~20/yr). 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.

Current observations of double neutron stars provide us with a wealth of information that we can use to investigate their evolutionary history and the physical conditions of neutron star formation. Understanding this history and formation conditions further allow us to make theoretical predictions for the formation of other double compact objects with one or two black hole components and assess the detectability of such systems by ground-based gravitational-wave interferometers. In this paper we summarize our group's body of work in the past few years and we place our conclusions and current understanding in the framework of other work in this area of astrophysical research.

The population of solitary compact objects in the Galaxy is very diffcult to investigate. In this paper we analyze the possibility of using microlensing searches to detect and to analyze the properties of the solitary black holes and neutron stars. Evolution of single and binary stars is considered using the StarTrack population synthesis code. We investigate the properties of the Galactic population of compact objects numerically. We find that the compact object lensing events are concentrated in a region with the radius of $\approx 5$ degrees around the Galactic center. The distribution of masses of the lenses for the models we consider differs but only slightly from the underlying massdistribution. The expected detection rates are of the order of a few per year.

We study the evolution of first star (Population III) binaries. Under specific conditions, these stars may produce high redshift gamma-ray bursts (GRBs). We demonstrate that the occurrence rate of GRBs does not depend sensitively on evolutionary parameters in the population synthesis models. We show that the first binaries may form a very small group (< 1%) of fast rotating stars through binary tidal interactions that make GRBs. This finding is contrary to the Bromm & Loeb assumption that all stars in close Population III binaries will be spun up by tides and produce a GRB. We find that there is simply not enough fast rotating stars in Population III binaries to expect detection with SWIFT. Predicted detection rates, even with very optimistic assumptions on binary fraction, evolutionary parameters and GRB detection, are very small: 0.1-0.01 per year.

We revisit the formation of PSR J0737-3039, taking into account the most recent observational constraints. We show that the most likely kick velocity and progenitor parameters depend strongly on the consideration of the full five-dimensional PDF for the magnitude and direction of the kick velocity imparted to pulsar B at birth, the mass of pulsar B's pre-supernova helium star progenitor, and the pre-supernova orbital separation, and on the adopted prior assumptions. The priors consist of the transverse systemic velocity, the age of the system, and the treatment of the unknown radial velocity. Since the latter cannot be determined from observation, we adopt a statistical approach and use theoretical radial-velocity distributions obtained from population synthesis calculations for coalescing double neutron stars. We find that the prior assumptions about the pre-supernova helium star mass affect the derived most likely parameters significantly: when the minimum helium star mass required for neutron star formation is assumed to be 2.1Msun, the most likely kick velocity ranges from 70-180km/s; when masses lower than 2.1Msun are assumed to allow neutron star formation, the most likely kick velocity can be as low as a few km/s, although the majority of the considered models still yield most likely kick velocities of 50-170km/s. We also show that the proximity of the double pulsar to the Galactic plane and the small proper motion do not pose stringent constraints on the kick velocity and progenitor mass of pulsar B. Instead, the constraints imposed by the orbital dynamics of asymmetric supernova explosions turn out to be much more restrictive. We conclude that the currently available observational constraints cannot be used to favor a specific core-collapse and neutron star formation mechanism. (abridged)

Recent pulsar surveys have increased the number of observed double neutron stars (DNS) in our galaxy enough so that observable trends in their properties are starting to emerge. In particular, it has been noted that the majority of DNS have eccentricities less than 0.3, which are surprisingly low for binaries that survive a supernova explosion that we believe imparts a significant kick to the neutron star. To investigate this trend, we generate many different theoretical distributions of DNS eccentricities using Monte Carlo population synthesis methods. We determine which eccentricity distributions are most consistent with the observed sample of DNS binaries. In agreement with Chaurasia & Bailes (2005), assuming all double neutron stars are equally as probable to be discovered as binary pulsars, we find that highly eccentric, coalescing DNS are less likely to be observed because of their accelerated orbital evolution due to gravitational wave emission and possible early mergers. Based on our results for coalescing DNS, we also find that models with vanishingly or moderately small kicks (sigma < about 50 km/s) are inconsistent with the current observed sample of such DNS. We discuss the implications of our conclusions for DNS merger rate estimates of interest to ground-based gravitational-wave interferometers. We find that, although orbital evolution due to gravitational radiation affects the eccentricity distribution of the observed sample, the associated upwards correction factor to merger rate estimates is rather small (typically 10-40%).

We present a theoretical study of double compact objects as potential short/hard gamma-ray burst (GRB) progenitors. An updated population synthesis code StarTrack is used to calculate properties of double neutron stars and black-hole neutron star binaries. We obtain their formation rates, estimate merger times and finally predict their most likely merger locations and afterglow properties for different types of host galaxies. Our results serve for a direct comparison with the recent HETE-II and SWIFT observations of several short bursts, for which afterglows and host galaxies were detected. We also discuss the possible constraints these observations put on the evolutionary models of double compact object formation. We emphasize that our double compact object models can successfully reproduce at the same time short GRBs within both young, star-forming galaxies (e.g., GRB 050709 and GRB 051221A), as well as within old, elliptical hosts (e.g., GRB 050724 and probably GRB 050509B)

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.

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.

We analyze the properties of VIRGO detector with the aim of studying its ability to search for coalescing black hole binaries. We focus on the remnants of the Population III stars, which currently should be massive black holes ($\sim 100-1000 M_\odot$), some of them bound in binary systems. The coalescence of such binaries due to emission of gravitational waves may be currently observable. We use a binary population synthesis to model the evolution of Population III binaries. We calculate the signal to noise ratios of gravitational waves emitted by the system in each of the coalescence phase: inspiral, merger and ringdown, and provide simple formulae for the signal to noise ratio as a function of masses of the binaries. We estimate the detection rates for the VIRGO interferometer and also compare them with the estimates for the current LIGO. We show that these expected rates are similar to, or larger than the expected rates from coalescences of Population I and II compact object binaries.

Recent Chandra observations have revealed a large population of faint X-ray point sources in the Galactic Centre. The observed population consists of about 2000 faint sources in the luminosity range ~10^31-10^33 erg/s. The majority of these sources (70%) are described by hard spectra, while the rest are rather soft. The nature of these sources still remains unknown. Belczynski & Taam (2004) demonstrated that X-ray binaries with neutron star or black hole accretors may account for most of the soft sources, but are not numerous enough to account for the observed number and X-ray properties of the faint hard sources. A population synthesis calculation of the Galactic Centre region has been carried out. Our results indicate that the numbers and X-ray luminosities of intermediate polars are consistent with the observed faint hard Galactic Centre population.

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.

As a first step towards understanding the angular momentum evolution history of black holes in merging black-hole/neutron-star binaries, we perform population synthesis calculations to track the distribution of accretion histories of compact objects in such binaries. We find that there are three distinct processes which can possibly contribute to the black-hole spin magnitude: a birth spin for the black hole, imparted at either (i) the collapse of a massive progenitor star to a black hole or (ii) the accretion-induced collapse of a neutron star to a black hole; and (iii) an accretion spin-up when the already formed black hole [via (i) or (ii)] goes through an accretion episode (through an accretion disk or a common-envelope phase). Our results show that, with regard to accretion-induced spinup in merging BH-NS binaries [method (iii) above], only {\em accretion episodes associated with common-envelope phases and hypercritical accretion rates} occur in the formation history of merging black hole/neutron star binaries. Lacking unambiguous experimental information about BH birth spins [i.e., regarding the results of processes (i) and (ii)], we choose two fiducial values for the BH birth angular momentum parameter a=J/M^2, consistent with observations of (i) NS birth spins (a roughly 0) and (ii) X-ray binaries (a=0.5). Using these two fiducial values and a conservative upper bound on the specific angular momentum of accreted matter, we discuss the expected range of black hole spins in the binaries of interest. We conclude with comments on the significance of these results for ground-based gravitational-wave searches of inspiral signals from black hole binaries.

The observed sample of double neutron-star (NS-NS) binaries presents a challenge to population-synthesis models of compact object formation: the parameters entering into these models must be carefully chosen so as to match (i) the observed star formation rate and (ii) the formation rate of NS-NS binaries, which can be estimated from the observed sample and the selection effects related to the discoveries with radio-pulsar surveys. In this paper, we select from an extremely broad family of possible population synthesis models those few (2%) which are consistent with the observed sample of NS-NS binaries. To further sharpen the constraints the observed NS-NS population places upon our understanding of compact-object formation processes, we separate the observed NS-NS population into two channels: (i) merging NS-NS binaries, which will inspiral and merge through the action of gravitational waves within $10 $ Gyr, and (ii) wide NS-NS binaries, consisting of all the rest. With the subset of astrophysically consistent models, we explore the implications for the rates at which double black hole (BH-BH), black hole-neutron star (BH-NS), and NS-NS binaries will merge through the emission of gravitational waves.

We use the StarTrack population synthesis code to discuss potential progenitor models of SN Ia: single degenerate scenario, semi-detached double white dwarf binary and double degenerate scenario. Using the most recent calculations of accretion onto white dwarfs, we consider SN Ia explosions of various white dwarfs in binary systems. We calculate the evolutionary delay times from the zero age main sequence to the explosion, and verify which scenarios can satisfy the constraints found by Strolger et al. (2004), i.e. a delay time of 2-4 Gyrs. It is found that the semi-detached double white dwarf binary model is inconsistent with observations. Only SN Ia progenitors in the double degenerate scenario are compatible with observations in all tested evolutionary models, with a characteristic delay time of about 3 Gyrs. If double degenerate mergers are excluded as SN Ia progenitors, we find that the single degenerate scenario may explain the observed delay times for models with low common envelope efficiency. It is also noted that the delay time distributions within the various tested SN Ia scenarios vary significantly and potentially can be used as an additional constraint on progenitor models.

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.

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.

Recent deep Chandra surveys of the Galactic center region have revealed the existence of a faint, hard X-ray source population. While the nature of this population is unknown, it is likely that several types of stellar objects contribute. For sources involving binary systems, accreting white dwarfs and accreting neutron stars with main sequence companions have been proposed. Among the accreting neutron star systems, previous studies have focused on stellar wind-fed sources. In this paper, we point out that binary systems in which mass transfer occurs via Roche lobe overflow (RLOF) can also contribute to this X-ray source population. A binary population synthesis study of the Galactic center region has been carried out, and it is found that evolutionary channels for neutron star formation involving the accretion induced collapse of a massive ONeMg white dwarf, in addition to the core collapse of massive stars, can contribute to this population. The RLOF systems would appear as transients with quiescent luminosities, above 2 keV, in the range from 10^31-10^32 ergs/s. The results reveal that RLOF systems primarily contribute to the faint X-ray source population in the Muno et al. (2003) survey and wind-fed systems can contribute to the less sensitive Wang et al. (2002) survey. However, our results suggest that accreting neutron star systems are not likely to be the major contributor to the faint X-ray source population in the Galactic center.

Formation rates of compact-object binaries are often derived from population synthesis calculations. However, such calculations depend sensitively on a relatively large number of model input parameters. Given considerable uncertainty in those model parameters, the predicted inspiral rates for double compact objects relevant to gravitational-wave interferometric detectors have been shown to be are uncertain by several orders of magnitude. Typically, inspiral rates are estimated for only a small set of models with a remarkably poor coverage of the highly multi-dimensional parameter space (primarily because of limited computer resources). Here, using as an example seven population-synthesis model parameters, we show that it is possible to derive fits of double-compact-object inspiral rates dependent simultaneously on all seven parameters. We find these fits to be accurate to 50% for binary black holes and to 40% for binary neutron stars. The availability of such fits implies that (i) depending on the problem of interest, it is not necessary to complete large numbers of computationally demanding population synthesis calculations; and (ii) for the first time, the sufficient exploration of the relevant phase space and the assessment of the uncertainties involved is not limited by computational resources and becomes feasible. Finally, we have also produced a histogram of the (a priori likely) binary black hole inpsiral rate, assuming our population synthesis models are equally likely. This histogram, effectively an a priori probability distribution for the BH-BH inspiral rate, suggests that merger rate is conservatively bounded below by 10^(-8)/yr/Milky-Way-galaxy.

In this paper we consider the population of eccentric binaries with a neutron star and a white dwarf that has been revealed in our galaxy in recent years through binary pulsar observations. We apply our statistical analysis method (Kim, Kalogera, & Lorimer 2003)and calculate the Galactic formation rate of these binaries empirically. We then compare our results with rate predictions based on binary population synthesis from various research groups and for various ranges of model input parameters. For our reference moel, we find the Galactic formation rate of these eccentric systems to be ~7 per Myr, about an order of magnitude smaller than results from binary evolution estimations. However, the empirical estimates are calculated with no correction for pulsar beaming, and therefore they should be taken as lower limits. Despite uncertainties that exceed an order of magnitude, there is significant overlap of the various rate calculations. This consistency lends confidence that our current understanding of the formation of these eccentric NS-WD binaries is reasonable.

We summarize our results on the Galactic merger rate of double neutron stars (DNS) in view of the recent discovery of PSR J0737-3039. We also present previously unpublished results for the 'global' probability distribution of merger rate values that incorporate the presently known systematics from the radio pulsar luminosity function. The most likely value obtained from the global distribution is only ~15 per Myr, but a re-analysis of the current pulsar sample and radio luminosities is needed for a reliable assessment of the best fitting distribution. Finally, we use our theoretical understanding of DNS formation to calculate a possible upper limit on the DNS merger rate from current Type Ib/c supernova rate estimates.

We examine the selection effects that determine how the population of inspiraling binary compact objects (BCOs) is reflected by those potentially observed with ground-based interferometers like LIGO. We lay the ground-work for the interpretation of future observations in terms of constraints on the real population and, correspondingly, binary star evolution models. To determine the extra-galactic population of inspiraling binaries we combine data on distance and blue luminosity from galaxy catalogs with current models of the galactic BCO mass distribution to simulate the physical distribution of binaries in the nearby universe. We use Monte Carlo methods to determine the fraction of binaries observable by the LIGO detectors from each galaxy as a function of the BCO chirp mass. We examine separately the role of source distance, sky position, time of detection, and binary system chirp mass on detection efficiency and selection effects relevant to the three LIGO detectors. Finally, we discuss the implications of the nearby geography of space on anticipated GW detection and compare our results to previous studies, which have assumed uniform galaxy volume density and fixed chirp mass for binary compact objects. From these considerations, actual BCO inspiral observations or significant upper limits on the coalescence rate anticipated in the near future by ground-based interferometers can be used to improve our knowledge of the galactic binary inspiral rate and to constrain models of radio pulsar characteristics and binary star evolution channels leading to neutron star or black hole binaries.

We use the well tested StarTrack binary population synthesis code to examine the properties of the population of compact object binaries. We calculate the distribution of masses and mass ratios taking into account weights introduced by observability in gravitational waves during inspiral. We find that in the observability weighted distribution of double neutron star binaries there are two peaks: one for nearly equal mass systems, and one for the systems consisting of a low and a high mass neutron star, $q=0.6-0.7$. The observability weighted distribution of black hole neutron star binaries is concentrated on the systems with the mass ratio $q=0.3-0.5$, while for the double black hole binaries the observability weighted distribution is dominated by the massive, nearly equal mass binaries with $q>0.7$.

Relativistic spin-orbit and spin-spin couplings has been shown to modify the gravitational waveforms expected from inspiraling binaries with a black hole and a neutron star. As a result inspiral signals may be missed due to significant losses in signal-to-noise ratio, if precession effects are ignored in gravitational-wave searches. We examine the sensitivity of the anticipated loss of signal-to-noise ratio on two factors: the accuracy of the precessing waveforms adopted as the true signals and the expected distributions of spin-orbit tilt angles, given the current understanding of their physical origin. We find that the results obtained using signals generated by approximate techniques are in good agreement with the ones obtained by integrating the 2PN equations. This shows that a complete account of all high-order post-Newtonian effects is usually not necessary for the determination of detection efficiencies. Based on our current astrophysical expectations, large tilt angles are not favored and as a result the decrease in detection rate varies rather slowly with respect to the black hole spin magnitude and is within 20--30% of the maximum possible values.

We present theoretical models of black hole (BH) populations in young stellar environments, such as starbursts and young star clusters. Using a population synthesis approach we compute the formation rates and characteristic properties of single and binary BHs for various representative ages and choices of parameters. We find that most of the BHs (typically 80% for an initial 50% binary fraction) are single, but with many originating from primordial binaries (which either merged into a single massive star or were disrupted following a supernova explosion). A smaller but significant fraction (typically 20%) of the BHs remain in binary systems. Main-sequence stars are the most frequent BH companions, but massive BH--BH binaries are the next most numerous group. The most massive BHs found in our simulations reach 80 Msun, and are formed through mergers of massive binary components. If formed in a dense star cluster such a massive stellar BH may become the seed for growth to a more massive (``intermediate-mass'') BH. Although we do not include dynamical interactions, our results provide realistic initial conditions for N-body simulations of dense star clusters (e.g., globular clusters) including primordial BHs.

Evolution of first population of massive metal-free binary stars is followed. Due to the low metallicity, the stars are allowed to form with large initial masses and to evolve without significant mass loss. Evolution at zero metallicity, therefore, may lead to the formation of massive remnants. In particular, black holes of intermediate-mass (100-500 Msun) are expected to have formed in early Universe, in contrast to the much lower mass stellar black holes (10 Msun) being formed at present. Following a natural assumption, that some of these Population III stars have formed in binaries, the physical properties of first stellar binary black holes are presented. We find that a significant fraction of such binary black holes coalesces within the Hubble time. We point out that burst of gravitational waves from the final coalescences and the following ringdown of these binary black hole mergers can be observed in the interferometric detectors. We estimate that advanced LIGO detection rate of such mergers is at least several events per year with high signal to noise ratio (>10).

Recent RXTE and Chandra discoveries of low mass X-ray binaries with ultra-short orbital periods have initiated theoretical work on the origins of these peculiar systems. Using the StarTrack population synthesis code the formation and evolution of X-ray ultracompact binaries (UCBs) in the Galactic field are analyzed. The relative number of UCBs with a neutron star or a black hole accretor populating our Galaxy is predicted. Our results demonstrate that standard evolutionary scenarios involving primordial binaries can be sufficient to produce the UCBs in the Galactic field without requiring additional processes associated with the dense stellar environments in the cores of globular clusters. In contrast to previous studies we find that the majority of the immediate progenitors of these systems consist of a hydrogen exhausted donor with an ONeMg white dwarf. The evolution of these systems leads to the accretion induced collapse of the white dwarf to a neutron star, which can play an important role in the formation of a majority of Galactic UCBs. We predict that with an increase in the number of X-ray active UCBs hosting neutron stars by an order of magnitude, a system with a black hole accretor may be found.

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.

Using the population synthesis code StarTrack we construct the first synthetic X-ray binary populations for direct comparison with the X-ray luminosity function (XLF) of NGC 1569 observed with Chandra. Our main goal is to examine whether it is possible to reproduce the XLF shape with our models, given the current knowledge for the star-formation history of this starburst galaxy. We thus produce hybrid models meant to represent the two stellar populations: one old, metal-poor with continuous star-formation for 1.5 Gyr and another recent and metal-rich population. To examine the validity of the models we compare XLFs calculated for varying ages of the populations and varying relative weights for the star-formation rates in the two populations. We find that, for typical binary evolution parameters, it is indeed possible to quite closely match the observed XLF shape. The robust match is achieved for an age of the young population and a ratio of star formation rates in the two populations that are within factors of 1.5 and 2, respectively, of those inferred from HST observations of NGC 1569. In view of this encouraging first step, we discuss the implications of our X-ray binary models and their potential as tools to study binary populations in galaxies.

Using the StarTrack binary population synthesis code we investigate the properties of population of compact object binaries. Taking into account the selection effects we calculate the expected properties of the observed binaries.We analyze possible constraints on the stellar evolution models and find that an observed sample of about one hundred mergers will yield strong constraints on the binary evolution scenarios.

Compact object mergers are possible progenitors of short bursts. We analyze properties of compact object mergers using the StarTrack population synthesis code, and find that the double neutron star population is dominated by short lived systems, thus they merge within host galaxies, while black hole neutron star binaries merge outside of the host galaxies. We summarize the differences between these two types of mergers.

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.

We use the StarTrack stellar population synthesis code to analyze properties of double compact object binaries as sources of gravitational waves. Since the distribution of lifetimes of these objects extends up to the Hubble time we conclude that a proper calculation of the expected rate must include a full cosmological model. We present such model, calculate the expected coalescence rates, and analyze the intrinsic sensitivity of these rates to the model assumptions. We find that the rate alone is a very poor estimator of the underlying stellar evolution model. However we show that the distribution of observed chirp masses is very sensitive to the underlying stellar evolution model, while it is very insensitive to the underlying cosmology, star formation rate history and variation of detector sensitivity.

We estimate the observed distribution of chirp masses of compact object binaries for the gravitational wave detectors. The stellar binary evolution is modeled using the {\em StarTrack} population synthesis code. The distribution of the predicted "observed" chirp masses vary with variation of different parameters describing stellar binary evolution. We estimate the sensitivity of the observed distribution to variation of these parameters and show which of the parameters can be constrained after observing 20, 100, and 500 compact object mergers. As a general feature of all our models we find that the population of observed binaries is dominated by the double black hole mergers.

We present an analysis of optical and near infrared spectra of the eclipsing symbiotic system FN Sgr. In particular, we have determined for the first time spectroscopic orbits based of the radial velocity curves for both components. A set of blue absorption lines resembling an A-F type star is present in all our spectra. They seem to be associated with the hot component but the best orbital fitting corresponds to a larger eccentricity. We have also studied spectral changes and photometric variations in function of both the hot component activity and the orbital motion.

We present an analysis of optical and near infrared spectra of the symbiotic system AE Ara, composed of an M5.5 III giant similar to the Galactic Bulge M giants and a hot luminous companion. In particular, we have determined for the first time spectroscopic orbits based of the radial velocity curves for the red and the hot component from the M-giant absorption lines and from the wings of He II 4686 emission profiles, respectively. We have also studied spectral changes and photometric variations in function of both the orbital phase and activity. The resulting physical parameters of the binary components and the nature of the hot component activity are briefly discussed.

Using the StarTrack population synthesis code we compute the distribution of masses of merging compact object (black hole or neutron star) binaries. The shape of the mass distribution is sensitive to some of the parameters governing the stellar binary evolution. We discuss the possibility of constraining stellar evolution models using mass measurements obtained from the detection of compact object inspiral with the upcoming gravitational-wave observatories.

Using the StarTrack population synthesis code we analyze the possible formation channels of the Galactic microquasar GRS 1915+105 harboring the most massive known stellar black hole. We find that we are able to reproduce the evolution and the current state of GRS 1915+105 in a model with lowered wind mass loss rates in massive stars. We conclude that the existence of GRS 1915+105 implies that stellar winds in massive stars are a factor of two weaker, but still within observational bounds, than used in stellar models. We discuss the mass transfer initiated by the massive primary in binaries which are the possible precursors of systems harboring massive black holes. We show that such mass transfer evolution is neither clear case B nor C, but instead is initiated during the core helium burning phase of the massive donor star, the progenitor of the black hole. We also argue that massive black holes do not receive any kicks at their formation and are formed through direct collapse of progenitor star without accompanying supernova explosion.

Using the StarTrack population synthesis code we analyze the formation channels possibly available to double neutron star binaries and find that they can be richer than previously thought. We identify a group of short lived, tight binaries, which do not live long enough to escape their host galaxies, despite their large center-of-mass velocities. We present our most recent results on all possible evolutionary paths leading to the formation of double neutron stars, calculate their coalescence rates, and also revisit the question of the distribution of merger sites around host galaxies. For a wide variety of binary evolution models and galaxy potentials, we find that most of neutron star mergers take place within galaxies. Our results stem from allowing for radial and common envelope evolution of helium-rich stars (testable in the future with detailed stellar-structure and hydrodynamic calculations) and indicate that double neutron star binaries may not be excluded as Gamma-Ray Burst (GRB) progenitors solely on the basis of their spatial distribution around host galaxies. We also find, in contrast to Bethe & Brown (1998), that in a significant fraction of common envelope (CE) phases neutron stars do not accrete enough material to become black holes, and thus the channels involving CEs are still open for the formation of double neutron stars.

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). 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: double neutron stars, black hole-neutron star systems and double black holes. We also examine their formation channels. 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 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.

We compute and present the distribution in mass of single and binary neutron stars, strange stars, and black holes. The calculations were performed using a stellar population synthesis code. We follow evolution of massive single stars as well as binaries with high mass primaries. The final product of the latter evolution can be either a binary composed of a white dwarf and a compact object (neutron star, black hole, strange star), two compact objects in a binary, or two single stars if the system was disrupted. We find in binaries a population of black holes which are more massive than single black holes which are a product of either binary or single evolution. We also find that if quark stars exist at all, their population can be as large as the population of black holes.

GRB data accumulated over the years have shown that the distribution of their time duration is bimodal. While there is some evidence that long bursts are associated with star-forming regions, nothing is known regarding the class of short bursts. Their very short timescales are hard to explain with the collapse of a massive star, but would be naturally produced by the merger of two compact objects, such as two neutron stars (NS), or a neutron star and a black hole (BH). As for the case of long bursts, afterglow obervations for short bursts should help reveal their origin. By using updated population synthesis code calculations, we simulate a cosmological population of merging NS-NS and NS-BH, and compute the distribution of their galactic off-sets, the density distribution of their environment, and, if indeed associated with GRBs, their expected afterglow characteristics.

Recently much work in studying Gamma-Ray Burst has been devoted to revealing the nature of outburst mechanism and studies of GRB afterglows. These issues have also been closely followed by the quest for identifying GRB progenitors. In this paper we consider the proposed binary star progenitors of GRBs: white dwarf neutron star binaries, white dwarf black hole binaries, helium core neutron star mergers, helium core black hole mergers, double neutron stars and neutron star black hole binaries. Using population synthesis methods we calculate merger rates of these binary progenitors and we compare them to the observed BATSE GRB rate. We also calculate the distribution of merger sites around host galaxies and compare them to the observed locations of GRB afterglows with respect to their hosts. We find that the rates of binary GRB progenitors in our standard model are lower than the observed GRB rates if GRBs are highly collimated. However, the uncertainty in the population synthesis results is too large to make this a firm conclusion. Although some observational signatures seem to point to collapsars as progenitors of long GRBs, we find that mergers of WD-NS, He-NS, He-BH, and NS-NS systems also trace the star formation regions of their host galaxies, as it is observed for long GRBs. We also speculate about possible progenitors of short-duration GRBs. For these, the most likely candidates are still mergers of compact objects. We find that the locations NS-NS and NS-BH mergers with respect to their hosts are significantly different. This may allow to distinguish between these two progenitor models, once current and near future missions, such as HETE-II or SWIFT, measure the locations of short GRBs.

In this paper we present a brief overview of population synthesis methods with a discussion of their main advantages and disadvantages. In the second part, we present some recent results from synthesis models of close binary compact objects with emphasis on the predicted rates, their uncertainties, and the model input parameters the rates are most sensitive to. We also report on a new evolutionary path leading to the formation of close double neutron stars (NS), with the unique characteristic that none of the two NS ever had the chance to be recycled by accretion. Their formation rates turn out to be comparable to or maybe even higher than those of recycled NS-NS binaries (like the ones observed), but their detection probability as binary pulsars is much smaller because of their short lifetimes. We discuss the implications of such a population for gravitational-wave detection of NS-NS inspiral events, and possibly for gamma-ray bursts and their host galaxies.

We report on a new evolutionary path leading to the formation of close double neutron stars (NS), with the unique characteristic that none of the two NS ever had the chance to be recycled by accretion. The existence of this channel stems from the evolution of helium-rich stars (cores of massive NS progenitors), which has been neglected in most previous studies of double compact object formation. We find that these non-recycled NS-NS binaries are formed from bare carbon-oxygen cores in tight orbits, with formation rates comparable to or maybe even higher than those of recycled NS-NS binaries. On the other hand, their detection probability as binary pulsars is greatly reduced (by about 1000) relative to recycled pulsars, because of their short lifetimes. We conclude that, in the context of gravitational-wave detection of NS-NS inspiral events, this new type of binaries calls for an increase of the rate estimates derived from the observed NS-NS with recycled pulsars, typically by factors of 1.5-3 or even higher.

Double compact objects (neutron stars and black holes) found in binaries with small orbital separations are known to spiral in and are expected to coalesce eventually because of the emission of gravitational waves. Such inspiral and merger events are thought to be primary sources for ground based gravitational-wave interferometric detectors (such as LIGO). Here, we present a brief review of estimates of coalescence rates and we examine the origin and relative importance of uncertainties associated with the rate estimates. For the case of double neutron star systems, we compare the most recent rate estimates to upper limits derived in a number of different ways. We also discuss the implications of the formation of close binaries with two non-recycled pulsars.

We compute and present the distribution in mass of single and binary neutron stars, strange stars, and black holes. The calculations were performed using a stellar population synthesis code. We follow all phases of single and binary evolution, starting from a ZAMS binary and ending in the creation of one compact object (neutron star, black hole, strange star) and a white dwarf, or two compact objects (single or binary). We assume that neutron stars are formed in the collapse of iron/nickel cores in the mass range M0 < M < M1, quark stars in the range M1 < M < M2, and black holes for core masses M > M2 and find that the population of quark stars can easily be as large as the population of black holes, even if there is only a small mass window for their formation.

Recently much work in Gamma-Ray Burst (GRB) studies was devoted to revealing the nature of outburst mechanism and to looking for GRB progenitors. Several types of progenitors were proposed for GRBs. Most promising objects are collapsars, compact object binaries, Helium star mergers and recently discussed supernovae. In this paper we consider four proposed binary star progenitors of GRBs: double neutron star (NS-NS), black hole neutron star (BH-NS), black hole white dwarf (BH-WD) mergers and Helium star mergers (He-BH). Helium star mergers are a possible outcome of common envelope evolution of a compact object entering the envelope of a giant with a helium core. Using population synthesis we calculate number of the binary progenitors and show that BH-WD and Helium star mergers dominate population of the proposed binary progenitors. Comparison of distribution of different binary mergers around galaxies they are born in, with localization of GRB afterglows in their host galaxies shows that only Helium mergers may be responsible for GRBs with observed afterglows while it excludes NS-NS, BH-NS and BH-WD systems as GRB progenitors. Assuming that all GRBs come from Helium star mergers and comparing numbers of Helium star mergers with observed BATSE GRB rate let us derive upper limit on GRB collimation to be 4 degrees.

We present a new catalogue of symbiotic stars. In our list we include 188 symbiotic stars as well as 28 objects suspected of being symbiotic. For each star, we present basic observational material: coordinates, V and K magnitudes, ultraviolet (UV), infrared (IR), X-ray and radio observations. We also list the spectral type of the cool component, the maximum ionization potential observed, references to finding charts, spectra, classifications and recent papers discussing the physical parameters and nature of each object. Moreover, we present the orbital photometric ephemerides and orbital elements of known symbiotic binaries, pulsational periods for symbiotic Miras, Hipparcos parallaxes and information about outbursts and flickering.

We use a stellar binary population synthesis code to find the lifetimes and velocities of several types of possible GRB progenitors: double neutron stars, black hole neutron stars, black hole white dwarfs, helium star mergers. Assuming that they are born in different types of galaxies we compute their spatial distribution and compare it with the observed locations of GRB afterglows within their hosts. We discuss constraints on the compact object merger model of GRBs imposed by this comparison and find that the observations of afterglows and their host galaxies appear inconsistent with the GRB compact object merger model.

Using stellar binary population synthesis code we calculate the production rates and lifetimes of several types of possible GRB progenitors. We consider mergers of double neutron stars, black hole neutron stars, black hole white dwarfs and helium star mergers. We calibrate the results with the measured star formation rate history. We discuss the viability of each GRB model, and alternatively assuming that all bursts are connected with one model we constrain the required collimation of GRBs. We also show the importance of widely used evolutionary parameters on the merger rates of calculated binary populations.

Compact object mergers are one of the favorite models of GRBs. It has been noted that in opposition to the collapsars, compact object mergers do not necessarily take place in the host galaxies, and may travel outside of them. With the discovery of afterglows and identification of host galaxies one can measure the distribution of GRBs with respect to their host galaxies. This distribution has been calculated using different population synthesis codes, and for different galactic potentials In this paper we compare the distributions of different types of compact object binaries: double neutron star systems (NS-NS), black hole -neutron star systems (BH-NS) and double black holes (BH-BH). We calculate the orbits and distributions of the projected distances on the sky for two extreme cases: a massive galaxy like the Milky Way, and empty space (corresponding to e.g. a globular cluster), and consider a wide range of possible kick velocity distributions. We find that BH-NS are more likely gamma-ray burst counterparts, since they lie close to the host galaxies, contrary to the NS-NS binaries.

Compact object mergers are one of the currently favored models for the origin of GRBs. The discovery of optical afterglows and identification of the nearest, presumably host, galaxies allows the analysis of the distribution of burst sites with respect to these galaxies. Using a model of stellar binary evolution we synthesize a population of compact binary systems which merge within the Hubble time. We include the kicks in the supernovae explosions and calculate orbits of these binaries in galactic gravitational potentials. We present the resulting distribution of merger sites and discuss the results in the framework of the observed GRB afterglows.

Mergers of compact objects may lead to different astrophysical phenomena: they may provide sources of observable gravitational radiation, and also may be connected with gamma-ray bursts. Estimate of the rate with which such mergers take place are based on assumptions about various parameters describing the binary evolution. The distribution of one of these parameters - the kick velocity a neutron star receives at birth will strongly influence the number and orbital parameters of compact objects binaries. We calculate these effect using population synthesis of binary stars and show that the expected compact object merger rate changes by a factor of 30 when the kick velocity varies within its current observational bounds.

The amplitude of the ellipsoidal variability, the mass function, and the evolutionary limits on the component masses have been used to constrain the binary system parameters. Contrary to all previous studies, our analysis shows that the mass ratio of T CrB q(M_giant/M_hot)=0.6, which implies a low mass binary system, with the stellar masses M_giant=0.7 solar mass, and M_hot=1.2 solar mass. This result strongly supports the thermonuclear runaway model for this recurrent nova, and solves all controversies about the nature of the hot component and the physical causes of its eruptions.