Black Hole-black Hole Mergers Research Articles

Impact of common envelope development criteria on the formation of LIGO/Virgo sources

The treatment and criteria for development of unstable Roche lobe overflow (RLOF) that leads to the common envelope (CE) phase have hindered the area of evolutionary predictions for decades. In particular, the formation of black hole-black hole (BH-BH), black hole-neutron star (BH-NS), and neutron star-neutron star (NS-NS) merging binaries depends sensitively on the CE phase in classical isolated binary evolution model. All these mergers are now reported as LIGO/Virgo sources or source candidates. CE is even considered by some as a mandatory phase in the formation of BH-BH, BH-NS, or NS-NS mergers in binary evolution models. At the moment, there is no full first-principles model for the development of the CE. We employed the StarTrack population synthesis code to test the current advancements in studies on the stability of RLOF for massive donors to assess their effect on the LIGO/Virgo source population. In particular, we allowed for more restrictive CE development criteria for massive donors (M > 18 M⊙). We also tested a modified condition for switching between different types of stable mass transfer and between the thermal or nuclear timescale. The implemented modifications significantly influence the basic properties of merging double compact objects, sometimes in non-intuitive ways. For one of the tested models, with restricted CE development criteria, the local merger rate density for BH-BH systems increased by a factor of 2–3 due to the emergence of a new dominant formation scenario without any CE phase. We find that the changes in highly uncertain assumptions on RLOF physics may significantly affect: (i) the local merger rate density; (ii) shape of the mass and mass ratio distributions; and (iii) dominant evolutionary formation (with and without CE) scenarios of LIGO/Virgo sources. Our results demonstrate that without sufficiently strong constraints on RLOF physics, it is not possible to draw fully reliable conclusions about the population of double compact object systems based on population synthesis studies.

Gamma-ray bursts associating with gravitational waves from compact binary mergers

Neutron star-neutron star, neutron star-black hole, and black hole-black hole mergers are the most promising sources of high-frequency gravitational waves (GWs) in the universe. Simultaneously, these compact binary mergers can also produce various electromagnetic emission represented by short gamma-ray burst (SGRB) phenomena, which had been confirmed by the discoveries in the GW170817/GRB 170817A event. Therefore, these merger events are undoubtedly one of the most important targets of the multi-messenger astronomical researches. We provide a brief review on the state of art of the SGRB studies. First of all, the emission properties during the prompt, precursor, extended, and afterglow phases of SGRBs are described subsequently by following a theoretical understanding for relativistic jets. For completeness, we also pay many attentions to the associated kilonova/mergernova emission from a quasi-isotropic merger ejecta, the potential neutrino production, and the host galaxies of SGRBs, which can provide independent constraints on the origin of SGRBs. The observations in the GW170817/GRB 170817A event are specially introduced, in particular, focusing on the implications for the merger product and the jet structure. Finally, by considering of the viewing angle-dependence of SGRB emission, we discuss a new understanding of the luminosity function of SGRBs. Then, we can prospect the future observation by the GECAM satellites, which are specially designed to monitor the GW-associated SGRBs and planned to be launched in 2020.

The evolution of kicked stellar-mass black holes in star cluster environments

We consider how dynamical friction acts on black holes that receive a velocity kick while located at the center of a gravitational potential, analogous to a star cluster, due to either a natal kick or the anisotropic emission of gravitational waves during a black hole-black hole merger. Our investigation specifically focuses on how well various Chandrasekhar-based dynamical friction models can predict the orbital decay of kicked black holes with $m_{bh} \lesssim 100 M_\odot$ due to an inhomogeneous background stellar field. In general, the orbital evolution of a kicked black hole follows that of a damped oscillator where two-body encounters and dynamical friction serve as sources of damping. However, we find models for approximating the effects of dynamical friction do not accurately predict the amount of energy lost by the black hole if the initial kick velocity $v_{k}$ is greater than the stellar velocity dispersion $\sigma$. For all kick velocities, we also find that two-body encounters with nearby stars can cause the energy evolution of a kicked BH to stray significantly from standard dynamical friction theory as encounters can sometimes lead to an energy gain. For larger kick velocities, we find the orbital decay of a black hole departs from classical theory completely as the black hole's orbital amplitude decays linearly with time as opposed to exponentially. Therefore, we have developed a linear decay formalism which scales linearly with black hole mass and $\frac{v_{k}}{\sigma}$ in order to account for the variations in the local gravitational potential.

Measuring the distance-redshift relation with the cross-correlation of gravitational wave standard sirens and galaxies

Gravitational waves from inspiraling compact binaries are known to be an excellent absolute distance indicator, yet it is unclear whether electromagnetic counterparts of these events are securely identified for measuring their redshifts, especially in the case of black hole-black hole mergers such as the one recently observed with the Advanced LIGO. We propose to use the cross-correlation between spatial distributions of gravitational wave sources and galaxies with known redshifts as an alternative means of constraining the distance-redshift relation from gravitational waves. In our analysis, we explicitly include the modulation of the distribution of gravitational wave sources due to weak gravitational lensing. We show that the cross-correlation analysis in next-generation observations will be able to tightly constrain the relation between the absolute distance and the redshift, and therefore constrain the Hubble constant as well as dark energy parameters.

DISCOVERY OF THE MASSIVE OVERCONTACT BINARY VFTS 352: EVIDENCE FOR ENHANCED INTERNAL MIXING

The contact phase expected to precede the coalescence of two massive stars is poorly characterized due to the paucity of observational constraints. Here we report on the discovery of VFTS 352, an O-type binary in the 30 Doradus region, as the most massive and earliest spectral type overcontact system known to date. We derived the 3D geometry of the system, its orbital period $P_{\rm orb}=1.1241452(4)$ d, components' effective temperatures -- $T_1=42\,540\pm280$ K and $T_2=41\,120\pm290$ K -- and dynamical masses -- $M_1=28.63\pm0.30 M_{\odot}$ and $M_2 = 28.85\pm0.30 M_{\odot}$. Compared to single-star evolutionary models, the VFTS 352 components are too hot for their dynamical masses by about 2700 and 1100 K, respectively. These results can be explained naturally as a result of enhanced mixing, theoretically predicted to occur in very short-period tidally-locked systems. The VFTS 352 components are two of the best candidates identified so far to undergo this so-called chemically homogeneous evolution. The future of VFTS 352 is uncertain. If the two stars merge, a very rapidly rotating star will be produced. Instead, if the stars continue to evolve homogeneously and keep shrinking within their Roche Lobes, coalescence can be avoided. In this case, tides may counteract the spin down by winds such that the VFTS 352 components may, at the end of their life, fulfill the requirements for long gamma-ray burst (GRB) progenitors in the collapsar scenario. Independently of whether the VFTS 352 components become GRB progenitors, this scenario makes VFTS 352 interesting as a progenitor of a black hole binary, hence as a potential gravitational wave source through black hole-black hole merger.

What powers the most relativistic jets? – I. BL Lacs

The dramatic relativistic jets pointing directly at us in BL Lac objects can be well modelled by bulk motion beaming of synchrotron self-Compton emission powered by a low Eddington fraction accretion flow. Nearly 500 of these AGN are seen in the 2nd Fermi Large Area Telescope catalogue of AGN. We combine the jet models which describe individual spectra with the expected jet parameter scalings with mass and mass accretion rate to predict the expected number of Fermi detected sources given the number densities of AGN from cosmological simulations. We select only sources with Eddington scaled mass accretion rate $< 0.01$ (i.e. radiatively inefficient flows), and include cooling, orientation effects, and the effects of absorption from pair production on the extragalactic IR background. These models overpredict the number of Fermi detected BL Lacs by a factor of 1000! This clearly shows that one of the underlying assumptions is incorrect, almost certainly that jets do not scale simply with mass and accretion rate. The most plausible additional parameter which can affect the region producing the Fermi emission is black hole spin. We can reproduce the observed numbers of BL Lacs if such relativistic jets are only produced by the highest spin ($a_*>0.8$) black holes, in agreement with the longstanding spin-jet paradigm. This also requires that high spins are intrinsically rare, as predicted by the cosmological simulations for growing black hole mass via chaotic (randomly aligned) accretion episodes, where only the most massive black holes have high spin due to black hole-black hole mergers.

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.

Nonspinning black hole-neutron star mergers: A model for the amplitude of gravitational waveforms

Black hole-neutron star binary mergers display a much richer phenomenology than black hole-black hole mergers, even in the relatively simple case - considered in this paper - in which both the black hole and the neutron star are nonspinning. When the neutron star is tidally disrupted, the gravitational wave emission is radically different from the black hole-black hole case and it can be broadly classified in two groups, depending on the spatial extent of the disrupted material. We present a phenomenological model for the gravitational waveform amplitude in the frequency domain that encompasses the three possible outcomes of the merger: no tidal disruption, "mild", and "strong" tidal disruption. The model is calibrated to general relativistic numerical simulations using piecewise polytropic neutron star equations of state. It should prove useful to extract information on the nuclear equation of state from future gravitational-wave observations, and also to obtain more accurate estimates of black hole-neutron star merger event rates in second- and third-generation interferometric gravitational-wave detectors. We plan to extend and improve the model as longer and more accurate gravitational waveforms become available, and we will make it publicly available online as a Mathematica package. We also present in appendix analytical fits of the projected KAGRA noise spectral density, that should be useful in data analysis applications.

First direct comparison of nondisrupting neutron star-black hole and binary black hole merger simulations

We present the first direct comparison of numerical simulations of neutron star-black hole and black hole-black hole mergers in full general relativity. We focus on a configuration with nonspinning objects and within the most likely range of mass ratio for neutron star-black hole systems ($q=6$). In this region of the parameter space, the neutron star is not tidally disrupted prior to merger, and we show that the two types of mergers appear remarkably similar. The effect of the presence of a neutron star on the gravitational wave signal is not only undetectable by the next generation of gravitational wave detectors, but also too small to be measured in the numerical simulations: even the plunge, merger and ringdown signals appear in perfect agreement for both types of binaries. The characteristics of the post-merger remnants are equally similar, with the masses of the final black holes agreeing within $\ensuremath{\delta}{M}_{\mathrm{BH}}&lt;5\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}4}{M}_{\mathrm{BH}}$ and their dimensionless spins within $\ensuremath{\delta}{\ensuremath{\chi}}_{\mathrm{BH}}&lt;{10}^{\ensuremath{-}3}$. The rate of periastron advance in the mixed binary agrees with previously published binary black hole results, and we use the inspiral waveforms to place constraints on the accuracy of our numerical simulations independent of algorithmic choices made for each type of binary. Overall, our results indicate that nondisrupting neutron star-black hole mergers are exceptionally well modeled by black hole-black hole mergers, and that given the absence of mass ejection, accretion disk formation, or differences in the gravitational wave signals, only electromagnetic precursors could prove the presence of a neutron star in low-spin systems of total mass $\ensuremath{\sim}10{M}_{\ensuremath{\bigodot}}$, at least until the advent of gravitational wave detectors with a sensitivity comparable to that of the proposed Einstein Telescope.

Intrinsic selection biases of ground-based gravitational wave searches for high-mass black hole-black hole mergers

The next generation of ground-based gravitational wave detectors may detect a few mergers of comparable-mass $M\ensuremath{\simeq}100\ensuremath{-}1000{M}_{\ensuremath{\bigodot}}$ [``intermediate-mass'' (IMBH)] spinning black holes. Black hole spin is known to have a significant impact on the orbit, merger signal, and post-merger ringdown of any binary with non-negligible spin. In particular, the detection volume for spinning binaries depends significantly on the component black hole spins. We provide a fit to the single-detector and isotropic-network detection volume versus (total) mass and arbitrary spin for equal-mass binaries. Our analysis assumes matched filtering to all significant available waveform power (up to $l=6$ available for fitting, but only $l\ensuremath{\le}4$ significant) estimated by an array of 64 numerical simulations with component spins as large as ${S}_{1,2}/{M}^{2}\ensuremath{\le}0.8$. We provide a spin-dependent estimate of our uncertainty, up to ${S}_{1,2}/{M}^{2}\ensuremath{\le}1$. For the initial (advanced) LIGO detector, our fits are reliable for $M\ensuremath{\in}[100,500]{M}_{\ensuremath{\bigodot}}$ ($M\ensuremath{\in}[100,1600]{M}_{\ensuremath{\bigodot}}$). In the online version of this article, we also provide fits assuming incomplete information, such as the neglect of higher-order harmonics. We briefly discuss how a strong selection bias towards aligned spins influences the interpretation of future gravitational wave detections of IMBH-IMBH mergers.

Black Hole Spin Evolution

We consider a subset of the physical processes that determine the spin j = a/M of astrophysical black holes. These include: (1) Initial conditions. Recent models suggest that the collapse of supermassive stars are likely to produce black holes with j ~ 0.7. (2) Major mergers. The outcome of a nearly equal mass black hole-black hole merger is not yet known, but we review the current best guesses and analytic bounds. (3) Minor mergers. We recover the result of Blandford & Hughes that accretion of small companions with isotropically distributed orbital angular momenta results in spindown, with j ~ M^{-7/3}. (4) Accretion. We present new results from fully relativistic magnetohydrodynamic accretion simulations. These show that, at least for one sequence of flow models, spin equilibrium (dj/dt = 0) is reached for j ~ 0.9, far less than the canonical value 0.998 of Thorne that was derived in the absence of MHD effects. This equilibrium value may not apply to all accretion flows, particularly thin disks. Nevertheless, it opens the possibility that black holes that have grown primarily through accretion are not maximally rotating.