Source author record

Krzysztof Belczynski

Krzysztof Belczynski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2022arXiv

A Nuclear Equation of State Inferred from Stellar r-process Abundances

Binary neutron star mergers (NSMs) have been confirmed as one source of the heaviest observable elements made by the rapid neutron-capture (r-) process. However, modeling NSM outflows -- from the total ejecta masses to their elemental yields -- depends on the unknown nuclear equation of state (EOS) that governs neutron-star structure. In this work, we derive a phenomenological EOS by assuming that NSMs are the dominant sources of the heavy element material in metal-poor stars with r-process abundance patterns. We start with a population synthesis model to obtain a population of merging neutron star binaries and calculate their EOS-dependent elemental yields. Under the assumption that these mergers were responsible for the majority of r-process elements in the metal-poor stars, we find parameters representing the EOS for which the theoretical NSM yields reproduce the derived abundances from observations of metal-poor stars. For our proof-of-concept assumptions, we find an EOS that is slightly softer than, but still in agreement with, current constraints, e.g., by the Neutron Star Interior Composition Explorer, with $R_{1.4}=12.25\pm 0.03$~km and $M_{\textrm TOV}$ of $2.17\pm 0.03$~M$_\odot$(statistical uncertainties, neglecting modeling systematics).

preprint2022arXiv

Modelling the formation of the first two neutron star-black hole mergers, GW200105 and GW200115: metallicity, chirp masses and merger remnant spins

The two neutron star-black hole mergers (GW200105 and GW200115) observed in gravitational waves by advanced LIGO and Virgo, mark the first ever discovery of such binaries in nature. We study these two neutron star-black hole systems through isolated binary evolution, using a grid of population synthesis models. Using both mass and spin observations (chirp mass, effective spin and remnant spin) of the binaries, we probe their different possible formation channels in different metallicity environments. Our models only support LIGO data when assuming the black hole is non spinning. Our results show a strong preference that GW200105 and GW200115 formed from stars with sub-solar metallicities $Z\lesssim 0.005$. Only two metal-rich ($Z=0.02$) models are in agreement with GW200115. We also find that chirp mass and remnant spins jointly aid in constraining the models, whilst the effective spin parameter does not add any further information. We also present the observable (i.e. post selection effects) median values of spin and mass distribution from all our models, which maybe used as a reference for future mergers. Further, we show that the remnant spin parameter distribution exhibits distinguishable features in different neutron star-black hole sub-populations. We find that non-spinning, first born black holes dominate significantly the merging neutron star-black hole population, ensuring electromagnetic counterparts to such mergers a rare affair.

preprint2022arXiv

The Effect of Supernova Convection On Neutron Star and Black Hole Masses

Our understanding of the convective-engine paradigm driving core-collapse supernovae has been used for 2 decades to predict the remnant mass distribution from stellar collapse. These predictions improve as our understanding of this engine increases. In this paper, we review our current understanding of convection (in particular, the growth rate of convection) in stellar collapse and study its effect on the remnant mass distribution. We show how the depth of the mass gap between neutron stars and black holes can help probe this convective growth. We include a study of the effects of stochasticity in both the stellar structure and the convective seeds caused by stellar burning. We study the role of rotation and its effect on the pair-instability mass gap. Under the paradigm limiting stellar rotation to those stars in tight binaries, we determine the effect of rotation on the remnant mass distribution.

preprint2020arXiv

BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

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.

preprint2020arXiv

Formation of low-spinning 100 Msun black holes

It is speculated that a merger of two massive stellar-origin BHs in a dense stellar environment may lead to the formation of a massive BH in the pair-instability mass gap (50-135 Msun). Such a merger-formed BH is expected to typically have a high spin (a=0.7). If such a massive BH acquires another BH it may lead to another merger detectable by LIGO/Virgo in gravitational waves. Acquiring a companion may be hindered by gravitational-wave kick/recoil, which accompanies the first merger and may quickly remove the massive BH from its parent globular or nuclear cluster. We test whether it is possible for a massive merger-formed BH in the pair-instability gap to be retained in its parent cluster and have low spin. Such a BH would be indistinguishable from a primordial BH. We employed results from numerical relativity calculations of black hole mergers to explore the range of gravitational-wave recoil velocities for various combinations of merging BH masses and spins. We compared merger-formed massive BH speeds with typical escape velocities from globular and nuclear clusters. We show that a globular cluster is highly unlikely to form and retain a 100 Msun BH if the spin of the BH is low (a<0.3) as such BHs acquire high recoil speeds (>200 km/s) that exceed typical escape speeds from globular clusters (50 km/s). However, a very low-spinning (a=0.1) and massive (100 Msun) BH could be formed and retained in a galactic nuclear star cluster. Even though such massive merger-formed BHs 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). A future detection of a massive BH in the pair-instability mass gap with low spin would therefore not be proof of the existence of primordial BHs, which are sometimes claimed to have low spins and arbitrarily high masses.

preprint2020arXiv

Non-interacting black hole binaries with Gaia and LAMOST

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.

preprint2020arXiv

Phase-dependent study of near-infrared disk emission lines in LB-1

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$β$ and Pa$γ$ (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$.

preprint2016arXiv

Compact Binary Merger Rates: Comparison with LIGO/Virgo Upper Limits

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.

preprint2016arXiv

Spectroscopy of Kerr black holes with Earth- and space-based interferometers

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.

preprint2016arXiv

The first gravitational-wave source from the isolated evolution of two 40-100 Msun stars

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.

preprint2015arXiv

Distinguishing types of compact-object binaries using the gravitational-wave signatures of their mergers

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.

preprint2015arXiv

Nature of the Extreme Ultraluminous X-ray Sources

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.

preprint2014arXiv

Black Hole X-ray Transients: The Formation Puzzle

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.

preprint2014arXiv

Comparison of LIGO/Virgo upper limits with predicted compact binary merger rates

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.

preprint2014arXiv

Prospects for Detection of Extragalactic Stellar Black Hole Binaries in the Nearby Universe

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.

preprint2014arXiv

Stellar Mergers Are Common

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.

preprint2014arXiv

The effect of helium accretion efficiency on rates of Type Ia supernovae: double-detonations in accreting binaries

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.

preprint2014arXiv

The Formation and Gravitational-Wave Detection of Massive Stellar Black-Hole Binaries

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.

preprint2013arXiv

Double Compact Objects II: Cosmological Merger Rates

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.

preprint2012arXiv

Common Proper Motion Wide White Dwarf Binaries Selected From The Sloan Digital Sky Survey

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.

preprint2012arXiv

Cyg X-3: a Galactic double black hole or black hole-neutron star progenitor

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.

preprint2012arXiv

Cyg X-3: a low-mass black hole or a neutron star

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.

preprint2012arXiv

Double Compact Objects I: The Significance Of The Common Envelope On Merger Rates

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.

preprint2012arXiv

High Mass X-ray Binaries: Future Evolution and Fate

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

preprint2012arXiv

The X-ray Star Formation Story as Told by Lyman Break Galaxies in the 4 Ms CDF-S

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.

preprint2012arXiv

X-ray Binary Evolution Across Cosmic Time

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.

preprint2011arXiv

On the apparent lack of be x-ray binaries with black holes in the galaxy and in the Magellanic clouds

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.

preprint2011arXiv

The coalescence rates of double black holes

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.

preprint2011arXiv

The fate of Cyg X-1: an empirical lower limit on BH-NS merger rate

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.

preprint2010arXiv

Double Compact Objects as Low-frequency Gravitational Wave Sources

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

preprint2010arXiv

IC10~X-1/NGC300~X-1: the very immediate progenitors of BH-BH binaries

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.

preprint2010arXiv

On The Maximum Mass of Stellar Black Holes

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.

preprint2010arXiv

Spectra of Type Ia Supernovae from Double Degenerate Mergers

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.

preprint2009arXiv

On the Apparent Lack of Be X-ray Binaries with Black Holes

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.

preprint2009arXiv

On The Origin Of The Highest Redshift Gamma-Ray Bursts

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.

preprint2009arXiv

Rates and Delay Times of Type Ia Supernovae

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

preprint2007arXiv

On the rarity of double black hole binaries: consequences for gravitational-wave detection

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.

preprint2006arXiv

The Lack of Gamma-Ray Bursts from Population III Binaries

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.

preprint2004arXiv

On the Chandra X-ray Sources in the Galactic Center

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.

preprint2002arXiv

A Comprehensive Study of Binary Compact Objects as Gravitational Wave Sources: Evolutionary Channels, Rates, and Physical Properties

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.

preprint1999arXiv

Distribution of black hole binaries around galaxies

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.