Source author record

K. Belczynski

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

42works
6topics
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

42 published item(s)

preprint2022arXiv

Black hole - black hole total merger mass and the origin of LIGO/Virgo sources

The LIGO-Virgo-KAGRA (LVK) Collaboration has reported nearly 100 BH-BH mergers. The LVK provides estimates of rates, masses, effective spins, and redshifts for these mergers. Yet, the formation channel(s) of the mergers remains uncertain. One way to search for a formation site is to contrast properties of detected BH-BH mergers with different models of BH-BH merger formation. Our study is designed to investigate the usefulness of the total BH-BH merger mass and its evolution with redshift in establishing the origin of gravitational-wave sources. We find that the average intrinsic BH-BH total merger mass shows exceptionally different behavior for the models that we adopt for our analysis. In the local universe (z=0) the average merger mass changes from Mtot,int=25Msun for CE binary evolution and open clusters formation channels, to Mtot,int=30Msun for the stable-RLOF binary channel, to Mtot,int=45Msun for the globular cluster channel. These differences are even more pronounced at larger redshifts. However, these differences are diminished when considering LVK O3 detector sensitivity. Comparison with LVK O3 data shows that none of our adopted models can match the data despite large errors on BH-BH masses and redshifts. We emphasize that our conclusions are derived from a small set of 6 models that are subject to numerous known uncertainties. We also note that BH-BH mergers may originate from a mix of several channels and that other (than those adopted here) BH-BH formation channels may exist.

preprint2022arXiv

Black hole - neutron star mergers: The first mass gap and kilonovae

Observations of X-ray binaries indicate a dearth of compact objects in the mass range from $\sim 2-5 M_{\odot}$. The existence of this (first mass) gap has been used to discriminate between proposed engines behind core-collapse supernovae. From LIGO/Virgo observations of binary compact remnant masses, several candidate first mass gap objects (either neutron stars (NSs) or black holes (BHs)) were identified during the O3 science run. Motivated by these new observations, we study the formation of BH-NS mergers in the framework of isolated classical binary evolution, using population synthesis methods to evolve large populations of binary stars (Population I and II) across cosmic time. We present results on the NS to BH mass ratios ($q=M_{\rm NS}/M_{\rm BH}$) in merging systems, showing that although systems with a mass ratio as low as $q=0.02$ can exist, typically BH-NS systems form with moderate mass ratios $q=0.1-0.2$. If we adopt a delayed supernova engine, we conclude that $\sim 30\%$ of BH-NS mergers may host at least one compact object in the first mass gap (FMG$^\circ$). Even allowing for uncertainties in the processes behind compact object formation, we expect the fraction of BH-NS systems ejecting mass during the merger to be small (from $\sim 0.6-9\%$). In our reference model, we assume: (i) the formation of compact objects within the FMG, (ii) natal NS/BH kicks decreased by fallback, (iii) low BH spins due to Tayler-Spruit angular momentum transport in massive stars. We find that $\lesssim 1\%$ of BH-NS mergers will have any mass ejection and about the same percentage will produce kilonova bright enough to have a chance of being detected with a large (Subaru-class) $8$m telescope. Interestingly, all these mergers will have both a BH and an NS in the FMG.

preprint2021arXiv

The Uncertain Future of Massive Binaries Obscures the Origin of LIGO/Virgo Sources

The LIGO/Virgo gravitational--wave observatories have detected 50 BH-BH coalescences. This sample is large enough to have allowed several recent studies to draw conclusions about the branching ratios between isolated binaries versus dense stellar clusters as the origin of double BHs. It has also led to the exciting suggestion that the population is highly likely to contain primordial black holes. Here we demonstrate that such conclusions cannot yet be robust, because of the large current uncertainties in several key aspects of binary stellar evolution. These include the development and survival of a common envelope, the mass and angular momentum loss during binary interactions, mixing in stellar interiors, pair-instability mass loss and supernova outbursts. Using standard tools such as the population synthesis codes StarTrack and COMPAS and the detailed stellar evolution code MESA, we examine as a case study the possible future evolution of Melnick 34, the most massive known binary star system. We show that, despite its well-known orbital architecture, various assumptions regarding stellar and binary physics predict a wide variety of outcomes: from a close BH-BH binary (which would lead to a potentially detectable coalescence), through a wide BH-BH binary (which might be seen in microlensing observations), or a Thorne-Zytkow object, to a complete disruption of both objects by pair-instability supernovae. Thus since the future of massive binaries is inherently uncertain, sound predictions about the properties of BH-BH systems are highly challenging at this time. Consequently, drawing conclusions about the formation channels for the LIGO/Virgo BH-BH merger population is premature.

preprint2020arXiv

Orbital period dependent spatial scatter of Galactic Black Hole Transients

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.

preprint2020arXiv

Synthetic catalog of black holes in the Milky Way

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

preprint2020arXiv

The evolutionary roads leading to low effective spins, high black hole masses, and O1/O2 rates of LIGO/Virgo binary black holes

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.

preprint2020arXiv

The Formation of a 70 Msun Black Hole at High Metallicity

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.

preprint2016arXiv

P13 an ULX that is a potential progenitor of merging BH-NS system

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 ÷0.6$yr$^{-1}$. This is the fourth empirical estimate of BH-NS merger rate.

preprint2016arXiv

Stability of mass transfer from massive giants: double black-hole binary formation and ultra-luminous X-ray sources

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.

preprint2016arXiv

The Effect of Pair-Instability Mass Loss on Black Hole Mergers

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.

preprint2015arXiv

A Model of the Stochastic Gravitational-Wave Background due to Core Collapse to Black Holes

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.

preprint2015arXiv

Cosmic Neutron Star Merger Rate and Gravitational Waves constrained by the R Process Nucleosynthesis

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.

preprint2015arXiv

Double Compact Objects III: Gravitational Wave Detection Rates

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.

preprint2015arXiv

Merger rates of double neutron stars and stellar origin black holes: The Impact of Initial Conditions on Binary Evolution Predictions

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.

preprint2015arXiv

On the formation and evolution of the first Be star in a black hole binary MWC 656

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.

preprint2014arXiv

Effect of metallicity on the gravitational-wave signal from the cosmological population of compact binary coalescences

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.

preprint2013arXiv

Modeling the Redshift Evolution of the Normal Galaxy X-ray Luminosity Function

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.

preprint2012arXiv

Gravitational wave background from population III binaries

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.

preprint2012arXiv

Impact of star formation inhomogeneities on merger rates and interpretation of LIGO results

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)

preprint2012arXiv

Missing Black Holes Unveil The Supernova Explosion Mechanism

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.

preprint2012arXiv

On the brightness distribution of Type Ia supernovae from violent white dwarf mergers

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

preprint2012arXiv

Probing the X-Ray Binary Populations of the Ring Galaxy NGC 1291

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

preprint2012arXiv

The Population of Helium-Merger Progenitors: Observational Predictions

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.

preprint2011arXiv

Compact Remnant Mass Function: Dependence on the Explosion Mechanism and Metallicity

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.

preprint2011arXiv

Delay Times and Rates for Type Ia Supernovae and Thermonuclear Explosions from Double-detonation Sub-Chandrasekhar Mass Models

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.

preprint2011arXiv

The eccentricity distribution of compact binaries

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

preprint2010arXiv

Binary compact object coalescence rates: The role of elliptical galaxies

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]

preprint2010arXiv

Black Hole Spin-Orbit Misalignment in Galactic X-ray Binaries

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.

preprint2010arXiv

Double and single recycled pulsars: an evolutionary puzzle?

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.

preprint2010arXiv

Predictions for the Rates of Compact Binary Coalescences Observable by Ground-based Gravitational-wave Detectors

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.

preprint2010arXiv

Sensitivity to Gravitational Waves from Compact Binary Coalescences Achieved during LIGO's Fifth and Virgo's First Science Run

We summarize the sensitivity achieved by the LIGO and Virgo gravitational wave detectors for compact binary coalescence (CBC) searches during LIGO's fifth science run and Virgo's first science run. We present noise spectral density curves for each of the four detectors that operated during these science runs which are representative of the typical performance achieved by the detectors for CBC searches. These spectra are intended for release to the public as a summary of detector performance for CBC searches during these science runs.

preprint2010arXiv

The effect of metallicity on the detection prospects for gravitational waves

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.

preprint2010arXiv

The LISA Gravitational Wave Foreground: A Study of Double White Dwarfs

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.

preprint2010arXiv

Type Ia Supernovae and Accretion Induced Collapse

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.

preprint2009arXiv

Probing Electron-Capture Supernovae: X-Ray Binaries in Starbursts

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.

preprint2007arXiv

Constraining population synthesis models via empirical binary compact object merger and supernovae rates

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)

preprint2007arXiv

Short Gamma-Ray Bursts and Binary Mergers in Spiral and Elliptical Galaxies: Redshift Distribution and Hosts

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)

preprint2006arXiv

Gravitational lensing as a probe of compact object population in the Galaxy

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.

preprint2005arXiv

Bounds on Expected Black Hole Spins in Inspiraling Binaries

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.

preprint2003arXiv

X-Ray Binary Populations: The Luminosity Function of NGC1569

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.

preprint2001arXiv

Binary Population Synthesis: Methods, Normalization, and Surprises

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.