The recent release of the second Gravitational-Wave Transient Catalog (GWTC-2) has increased significantly the number of known GW events, enabling unprecedented constraints on formation models of ...compact binaries. One pressing question is to understand the fraction of binaries originating from different formation channels, such as isolated field formation versus dynamical formation in dense stellar clusters. In this paper, we combine the cosmic binary population synthesis suite and the cmc code for globular cluster evolution to create a mixture model for black hole binary formation under both formation scenarios. For the first time, these code bodies are combined self-consistently, with cmc itself employing cosmic to track stellar evolution. We then use a deep-learning enhanced hierarchical Bayesian analysis to continuously sample over and constrain the common envelope efficiency α assumed in cosmic, the initial cluster virial radius r v adopted in cmc, and the intrinsic mixture fraction f between each channel. Under specific assumptions about other uncertain aspects of isolated binary and globular cluster evolution, we report the median and 90% confidence interval of three physical parameters for the intrinsic population ... This simultaneous constraint agrees with observed properties of globular clusters in the Milky Way and is an important first step in the pathway toward learning the astrophysics of compact binary formation through GW observations. (ProQuest: … denotes formulae omitted.)
Using state-of-the-art dynamical simulations of globular clusters, including radiation reaction during black hole encounters and a cosmological model of star cluster formation, we create a realistic ...population of dynamically formed binary black hole mergers across cosmic space and time. We show that in the local universe, 10% of these binaries form as the result of gravitational-wave emission between unbound black holes during chaotic resonant encounters, with roughly half of those events having eccentricities detectable by current ground-based gravitational-wave detectors. The mergers that occur inside clusters typically have lower masses than binaries that were ejected from the cluster many Gyrs ago. Gravitational-wave captures from globular clusters contribute 1 – 2 Gpc−3 yr−1 to the binary merger rate in the local universe, increasing to ≳ 10 Gpc−3 yr−1 at z ∼ 3 . Finally, we discuss some of the technical difficulties associated with post-Newtonian scattering encounters, and how care must be taken when measuring the binary parameters during a dynamical capture.
We explore the long-term evolution of mass-transferring white dwarf (WD) binaries undergoing both direct-impact and disk accretion and explore implications of such systems to gravitational-wave (GW) ...astronomy. We cover a broad range of initial component masses and show that these systems, the majority of which lie within the Laser Interferometer Space Antenna (LISA) sensitivity range, exhibit prominent negative orbital frequency evolution (chirp) for a significant fraction of their lifetimes. Using a galactic population synthesis, we predict ∼2700 of these systems will be observable with a negative chirp of 0.1 yr−2 by a space-based GW detector like LISA. We also show that detections of mass-transferring double WD systems by LISA may provide astronomers with unique ways of probing the physics governing close compact object binaries.
Abstract
LIGO’s third observing run (O3) has reported several neutron star–black hole (NSBH) merger candidates. From a theoretical point of view, NSBH mergers have received less attention in the ...community than either binary black holes, or binary neutron stars. Here we examine single–single (sin–sin) gravitational wave (GW) captures in different types of star clusters—galactic nuclei, globular clusters, and young stellar clusters—and compare the merger rates from this channel to other proposed merger channels in the literature. There are currently large uncertainties associated with every merger channel, making a definitive conclusion about the origin of NSBH mergers impossible. However, keeping these uncertainties in mind, we find that sin–sin GW capture is unlikely to significantly contribute to the overall NSBH merger rate. In general, it appears that isolated binary evolution in the field or in clusters, and dynamically interacting binaries in triple configurations, may result in a higher merger rate.
The formation and evolution of binary stars are critical components of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling or merging compact ...binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observations of electromagnetic and gravitational wave sources. Binary population synthesis is a tool that balances physical modeling with simulation speed to produce large binary populations on timescales of days. We present a community-developed binary population synthesis suite, COSMIC, which is designed to simulate compact-object binary populations and their progenitors. As a proof of concept, we simulate the Galactic population of compact binaries and their gravitational wave signals observable by the Laser Interferometer Space Antenna.
Abstract
Compact-object binaries including a white dwarf component are unique among gravitational-wave sources because their evolution is governed not just by general relativity and tides, but also ...by mass transfer. While the black hole and neutron star binaries observed with ground-based gravitational-wave detectors are driven to inspiral due to the emission of gravitational radiation—manifesting as a “chirp-like” gravitational-wave signal—the astrophysical processes at work in double white dwarf (DWD) systems can cause the inspiral to stall and even reverse into an outspiral. The dynamics of the DWD outspiral thus encode information about tides, which tell us about the behavior of electron-degenerate matter. We carry out a population study to determine the effect of the strength of tides on the distributions of the DWD binary parameters that the Laser Interferometer Space Antenna (LISA) will be able to constrain. We find that the strength of tidal coupling parameterized via the tidal synchronization timescale at the onset of mass transfer affects the distribution of gravitational-wave frequencies and frequency derivatives for detectably mass-transferring DWD systems. Using a hierarchical Bayesian framework informed by binary population synthesis simulations, we demonstrate how this parameter can be inferred using LISA observations. By measuring the population properties of DWDs, LISA will be able to probe the behavior of electron-degenerate matter.
Abstract
Theoretical modeling of massive stars predicts a gap in the black hole (BH) mass function above ∼40–50
M
⊙
for BHs formed through single star evolution, arising from (pulsational) ...pair-instability supernovae (PISNe). However, in dense star clusters, dynamical channels may exist that allow construction of BHs with masses in excess of those allowed from single star evolution. The detection of BHs in this so-called “upper-mass gap” would provide strong evidence for the dynamical processing of BHs prior to their eventual merger. Here, we explore in detail the formation of BHs with masses within or above the pair-instability gap through collisions of young massive stars in dense star clusters. We run a suite of 68 independent cluster simulations, exploring a variety of physical assumptions pertaining to growth through stellar collisions, including primordial cluster mass segregation and the efficiency of envelope stripping during collisions. We find that as many as ∼20% of all BH progenitors undergo one or more collisions prior to stellar collapse and up to ∼1% of all BHs reside within or above the pair-instability gap through the effects of these collisions. We show that these BHs readily go on to merge with other BHs in the cluster, creating a population of massive BH mergers at a rate that may compete with the “multiple-generation” merger channel described in other analyses. This has clear relevance for the formation of very massive BH binaries as recently detected by the Laser Interferometer Gravitational-Wave Observatory/Virgo in GW190521. Finally, we describe how stellar collisions in clusters may provide a unique pathway to PISNe and briefly discuss the expected rate of these events and other electromagnetic transients.
Globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the so-called core-collapsed and non-core-collapsed clusters. Here, we use our ...Hénon-type Monte Carlo code, CMC, to explore initial cluster parameters that map into this bimodality. Remarkably, we find that by varying the initial size of clusters (specified in our initial conditions in terms of the initial virial radius, rv) within a relatively narrow range consistent with the measured radii of young star clusters in the local universe (rv 0.5-5 pc), our models reproduce the variety of present-day cluster properties. Furthermore, we show that stellar-mass black holes (BHs) play an intimate role in this mapping from initial conditions to the present-day structural features of GCs. We identify "best-fit" models for three GCs with known observed BH candidates, NGC 3201, M22, and M10, and show that these clusters harbor populations of ∼50-100 stellar-mass BHs at present. As an alternative case, we also compare our models to the core-collapsed cluster NGC 6752 and show that this cluster likely contains few BHs at present. Additionally, we explore the formation of BH binaries in GCs and demonstrate that these systems form naturally in our models in both detached and mass-transferring configurations with a variety of companion stellar types, including low-mass main-sequence stars, white dwarfs, and sub-subgiants.
Over 100 millisecond radio pulsars (MSPs) have been observed in globular clusters (GCs), motivating theoretical studies of the formation and evolution of these sources through stellar evolution ...coupled to stellar dynamics. Here we study MSPs in GCs using realistic N-body simulations with our Cluster Monte Carlo code. We show that neutron stars (NSs) formed in electron-capture supernovae (including both accretion-induced and merger-induced collapse of white dwarfs) can be spun up through mass transfer to form MSPs. Both NS formation and spin-up through accretion are greatly enhanced through dynamical interaction processes. We find that our models for average GCs at the present day with masses 2 × 105 M can produce up to 10-20 MSPs, while a very massive GC model with mass 106 M can produce close to 100. We show that the number of MSPs is anti-correlated with the total number of stellar-mass black holes (BHs) retained in the host cluster. The radial distributions are also affected: MSPs are more concentrated toward the center in a host cluster with a smaller number of retained BHs. As a result, the number of MSPs in a GC could be used to place constraints on its BH population. Some intrinsic properties of MSP systems in our models (such as the magnetic fields and spin periods) are in good overall agreement with observations, while others (such as the distribution of binary companion types) are less so, and we discuss the possible reasons for such discrepancies. Interestingly, our models also demonstrate the possibility of dynamically forming NS-NS and NS-BH binaries in GCs, although the predicted numbers are very small.
We explore the formation of mass-transferring binary systems containing black holes (BHs) within globular clusters (GC). We show that it is possible to form mass-transferring BH binaries with main ...sequence, giant, and white dwarf companions with a variety of orbital parameters in GCs spanning a large range in present-day properties. All mass-transferring BH binaries found in our models at late times are dynamically created. The BHs in these systems experienced a median of ∼30 dynamical encounters within the cluster before and after acquiring the donor. Furthermore, we show that the presence of mass-transferring BH systems has little correlation with the total number of BHs within the cluster at any time. This is because the net rate of formation of BH-non-BH binaries in a cluster is largely independent of the total number of retained BHs. Our results suggest that the detection of a mass-transferring BH binary in a GC does not necessarily indicate that the host cluster contains a large BH population.