On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of 1.0×10(-21). It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ringdown of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203,000 years, equivalent to a significance greater than 5.1σ. The source lies at a luminosity distance of 410(-180)(+160)  Mpc corresponding to a redshift z=0.09(-0.04)(+0.03). In the source frame, the initial black hole masses are 36(-4)(+5)M⊙ and 29(-4)(+4)M⊙, and the final black hole mass is 62(-4)(+4)M⊙, with 3.0(-0.5)(+0.5)M⊙c(2) radiated in gravitational waves. All uncertainties define 90% credible intervals. These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.

/pdf/the-observation-of-gravitational-waves-from-a-binary-black-58srpupgg9.pdf

The Observation of Gravitational Waves from a Binary Black Hole Merger

We report on gravitational wave discoveries from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 15:00 UTC and 1 October 2019 15:00. By imposing a false-alarm-rate threshold of two per year in each of the four search pipelines that constitute our search, we present 39 candidate gravitational wave events. At this threshold, we expect a contamination fraction of less than 10%. Of these, 26 candidate events were reported previously in near real-time through GCN Notices and Circulars; 13 are reported here for the first time. The catalog contains events whose sources are black hole binary mergers up to a redshift of ~0.8, as well as events whose components could not be unambiguously identified as black holes or neutron stars. For the latter group, we are unable to determine the nature based on estimates of the component masses and spins from gravitational wave data alone. The range of candidate events which are unambiguously identified as binary black holes (both objects $\geq 3~M_\odot$) is increased compared to GWTC-1, with total masses from $\sim 14~M_\odot$ for GW190924_021846 to $\sim 150~M_\odot$ for GW190521. For the first time, this catalog includes binary systems with significantly asymmetric mass ratios, which had not been observed in data taken before April 2019. We also find that 11 of the 39 events detected since April 2019 have positive effective inspiral spins under our default prior (at 90% credibility), while none exhibit negative effective inspiral spin. Given the increased sensitivity of Advanced LIGO and Advanced Virgo, the detection of 39 candidate events in ~26 weeks of data (~1.5 per week) is consistent with GWTC-1.

GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run

We report on the population of 47 compact binary mergers detected with a false-alarm rate of 0.01 are dynamically assembled. Third, we estimate merger rates, finding RBBH = 23.9-+8.614.3 Gpc-3 yr-1 for BBHs and RBNS = 320-+240490 Gpc-3 yr-1 for binary neutron stars. We find that the BBH rate likely increases with redshift (85% credibility) but not faster than the star formation rate (86% credibility). Additionally, we examine recent exceptional events in the context of our population models, finding that the asymmetric masses of GW190412 and the high component masses of GW190521 are consistent with our models, but the low secondary mass of GW190814 makes it an outlier.

/pdf/population-properties-of-compact-objects-from-the-second-34o8o76ct0.pdf

Population Properties of Compact Objects from the Second LIGO-Virgo Gravitational-Wave Transient Catalog

Many questions in physical cosmology regarding the thermal history of the intergalactic medium, chemical enrichment, reionization, etc., are thought to be intimately related to the nature and evolution of pregalactic structure. In particular, the efficiency of primordial star formation and the primordial initial mass function are of special interest. We present results from high-resolution three-dimensional adaptive mesh re—nement simulations that follow the collapse of primordial molecular clouds and their subsequent fragmentation within a cosmologically representative volume. Comoving scales from 128 kpc down to 1 pc are followed accurately. Dark matter dynamics, hydrodynamics, and all relevant chemical and radiative processes (cooling) are followed self-consistently for a cluster-normalized cold dark matter (CDM) structure formation model. Primordial molecular clouds with D105 solar masses are assembled by mergers of multiple objects that have formed hydrogen molecules in the gas phase with a fractional abundance of As the subclumps merge, cooling lowers the temperature to D200 K in a ii cold (10~4. pocket ˇˇ at the center of the halo. Within this cold pocket, a quasi-hydrostatically contracting core with mass D200 and number densities cm~3 are found. We —nd that less than 1% of the primor- M _ Z105 dial gas in such small-scale structures cools and collapses to sufficiently high densities to be available for primordial star formation. Furthermore, it is worthwhile to note that this study achieved the highest dynamic range covered by structured adaptive mesh techniques in cosmological hydrodynamics to date. Subject headings: cosmology: theorygalaxies: formationmethods: numerical

The formation and fragmentation of primordial molecular clouds

(abridged) We use large cosmological simulations to study the origin of primordial star-forming clouds in a Lambda CDM universe, by following the formation of dark matter halos and the cooling of gas within them. To model the physics of chemically pristine gas, we employ a non-equilibrium treatment of the chemistry of 9 species and include cooling by molecular hydrogen. We explore the hierarchical growth of bound structures forming at redshifts z = 25 - 30 with total masses in the range 10^5 - 10^6 Msun. The complex interplay between the gravitational formation of dark halos and the thermodynamic and chemical evolution of the gas clouds compromises analytic estimates of the critical H2 fraction. Dynamical heating from mass accretion and mergers opposes relatively inefficient cooling by molecular hydrogen, delaying the production of star-forming clouds in rapidly growing halos. We also investigate the impact of photo-dissociating ultra-violet (UV) radiation on the formation of primordial gas clouds. We consider two extreme cases by first including a uniform radiation field in the optically thin limit and secondly by accounting for the maximum effect of gas self-shielding in virialized regions. In both the cases we consider, the overall impact can be described by computing an equilibrium H2 abundance for the radiation flux and defining an effective shielding factor. 
Based on our numerical results, we develop a semi-analytic model of the formation of the first stars, and demonstrate how it can be coupled with large N-body simulations to predict the star formation rate in the early universe.

/pdf/simulations-of-early-structure-formation-primordial-gas-34visb3q35.pdf

Simulations of Early Structure Formation: Primordial Gas Clouds

The astrophysical origin of gravitational wave (GW) events discovered by LIGO/VIRGO remains an outstanding puzzle. In active galactic nuclei (AGN), compact-object binaries form, evolve, and interact with a dense star cluster and a gas disk. An important question is whether and how binaries merge in these environments. To address this question, we have performed one-dimensional $N$-body simulations combined with a semi-analytical model which includes the formation, disruption, and evolution of binaries self-consistently. We point out that binaries can form in single-single interactions by the dissipation of kinetic energy in a gaseous medium. This ``gas capture'' binary formation channel contributes up to $97\,\%$ of gas-driven mergers and leads to a high merger rate in AGN disks even without pre-existing binaries. We find the merger rate to be in the range $\sim 0.02-60\,\mathrm{Gpc^{-3}yr^{-1}}$. The results are insensitive to the assumptions on gaseous hardening processes: we find that once they are formed, binaries merge efficiently via binary-single interactions even if these gaseous processes are neglected. We find that the average number of mergers per BH is $0.4$, and the probability for repeated mergers in 30 Myr is $\sim 0.21-0.45$. High BH masses due to repeated mergers, high eccentricities, and a significant Doppler drift of GWs are promising signatures which distinguish this merger channel from others. Furthermore, we find that gas-capture binaries reproduce the distribution of LMXBs in the Galactic center, including an outer cutoff at $\sim1$ pc due to the competition between migration and hardening by gas torques.

Formation and Evolution of Compact Object Binaries in AGN Disks

The recently discovered gravitational wave sources GW190521 and GW190814 have shown evidence of BH mergers with masses and spins that could be outside of the range expected from isolated stellar evolution. These merging objects could have undergone previous mergers. Such hierarchical mergers are predicted to be frequent in active galactic nuclei (AGN) disks, where binaries form and evolve efficiently by dynamical interactions and gaseous dissipation. Here we compare the properties of these observed events to the theoretical models of mergers in AGN disks, which are obtained by performing one-dimensional $N$-body simulations combined with semi-analytical prescriptions. The high BH masses in GW190521 are consistent with mergers of high-generation (high-g) BHs where the initial progenitor stars had high metallicity, 2g BHs if the original progenitors were metal-poor, or 1g BHs that had gained mass via super-Eddington accretion. Other measured properties related to spin parameters in GW190521 are also consistent with mergers in AGN disks. Furthermore, mergers in the lower mass gap or those with low mass ratio as found in GW190814 and GW190412 are also reproduced by mergers of 2g-1g or 1g-1g objects with significant accretion in AGN disks. Finally, due to gas accretion, the massive neutron star merger reported in GW190425 can be produced in an AGN disk.

Mass-gap Mergers in Active Galactic Nuclei

The heaviest neutron stars and lightest black holes expected to be produced by stellar evolution leave the mass-range $2.2$ M$_{\odot}\lesssim m \lesssim 5$ M$_\odot$ largely unpopulated. Objects found in this so-called lower mass gap likely originate from a distinct astrophysical process. Such an object, with mass $2.6$ M$_\odot$ was recently detected in the binary merger GW190814 through gravitational waves by LIGO/Virgo. Here we show that black holes in the mass gap are naturally assembled through mergers and accretion in AGN disks, and can subsequently participate in additional mergers. We compute the properties of AGN-assisted mergers involving neutron stars and black holes, accounting for accretion. We find that mergers in which one of the objects is in the lower mass gap represent up to $4$% of AGN-assisted mergers detectable by LIGO/Virgo. The lighter object of GW190814, with mass $2.6$ M$_\odot$, could have grown in an AGN disk through accretion. We find that the unexpectedly high total mass of 3.4 M$_\odot$ observed in the neutron star merger GW190425 may also be due to accretion in an AGN disk.

Black Hole Formation in the Lower Mass Gap through Mergers and Accretion in AGN Disks

The astrophysical origin of gravitational wave (GW) events is one of the most timely problems in the wake of the LIGO/Virgo discoveries. In active galactic nuclei (AGN), binaries form and evolve efficiently by dynamical interactions and gaseous dissipation. Previous studies have suggested that binary black hole (BBH) mergers in AGN disks can contribute significantly to BBH mergers observed by GW interferometers. Here we examine the distribution of the effective spin parameter $\chi_\mathrm{eff}$ of this GW source population. We extend our semi-analytical model of binary formation and evolution in AGN disks by following the evolution of the binary orbital angular momenta and black hole (BH) spins. BH spins change due to gas accretion and BH mergers, while the binary orbital angular momenta evolve due to gas accretion and binary-single interactions. We find that the distribution of $\chi_\mathrm{eff}$ predicted by our AGN model is similar to the distribution observed during LIGO/Virgo O1 and O2. On the other hand, if radial migration of BHs is inefficient, $\chi_\mathrm{eff}$ is skewed toward higher values compared with the observed distribution, because of the paucity of scattering events that would randomize spin directions relative to the orbital plane. We suggest that high binary masses and the positive correlation between binary mass and the standard deviation of $\chi_\mathrm{eff}$ for chirp masses up to $\approx 20$ $\mathrm{M}_\odot$, can be possible signatures for mergers originating in AGN disks. Finally, hierarchical mergers in AGN disks naturally produce properties of the recent GW event GW190412, including a low mass ratio, a high primary BH spin, and a significant spin component in the orbital plane.

Spin Evolution of Stellar-mass Black Hole Binaries in Active Galactic Nuclei

Approximately two hundred supermassive black holes (SMBHs) have been discovered within the first $\sim$Gyr after the Big Bang. One pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs). A possible obstacle to this scenario is that the collapsing gas fragments and forms a cluster of main-sequence stars. Here we raise the possibility that stellar collisions may be sufficiently frequent and energetic to inhibit the contraction of the massive protostar, avoiding strong UV radiation driven outflows, and allowing it to continue growing into an SMS. We investigate this scenario with semianalytic models incorporating star formation, gas accretion, dynamical friction from stars and gas, stellar collisions, and gas ejection. We find that when the collapsing gas fragments at a density of $\lesssim 3\times 10^{10}\,\mathrm{cm^{-3}}$, the central protostar contracts due to infrequent stellar mergers, and in turn photoevaporates the remaining collapsing gas, resulting in the formation of a $\lesssim 10^4~{\rm M_\odot}$ object. On the other hand, when the collapsing gas fragments at higher densities (expected for a metal-poor cloud with $Z\lesssim10^{-5}\,{\rm Z_\odot}$ with suppressed ${\rm H_2}$ abundance) the central protostar avoids contraction and keeps growing via frequent stellar mergers, reaching masses as high as $\sim 10^5-10^6\,{\rm M_\odot}$. We conclude that frequent stellar mergers represent a possible pathway to form massive BHs in the early universe.

Hiromichi Tagawa

Papers

Formation and Evolution of Compact Object Binaries in AGN Disks

Mass-gap Mergers in Active Galactic Nuclei

Black Hole Formation in the Lower Mass Gap through Mergers and Accretion in AGN Disks

Spin Evolution of Stellar-mass Black Hole Binaries in Active Galactic Nuclei

Making a Supermassive Star by Stellar Bombardment