Arxiv:2012.00011V2 [Astro-Ph.HE] 20 Dec 2020 Mergers in Gas-Rich Environments (Mckernan Et Al
Total Page:16
File Type:pdf, Size:1020Kb
Draft version December 22, 2020 Preprint typeset using LATEX style emulateapj v. 12/16/11 MASS-GAP MERGERS IN ACTIVE GALACTIC NUCLEI Hiromichi Tagawa1, Bence Kocsis2, Zoltan´ Haiman3, Imre Bartos4, Kazuyuki Omukai1, Johan Samsing5 1Astronomical Institute, Graduate School of Science, Tohoku University, Aoba, Sendai 980-8578, Japan 2 Rudolf Peierls Centre for Theoretical Physics, Clarendon Laboratory, Parks Road, Oxford OX1 3PU, UK 3Department of Astronomy, Columbia University, 550 W. 120th St., New York, NY, 10027, USA 4Department of Physics, University of Florida, PO Box 118440, Gainesville, FL 32611, USA 5Niels Bohr International Academy, The Niels Bohr Institute, Blegdamsvej 17, 2100 Copenhagen, Denmark Draft version December 22, 2020 ABSTRACT The recently discovered gravitational wave sources GW190521 and GW190814 have shown evidence of BH mergers with masses and spins 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. Keywords: binaries: close { gravitational waves {galaxies: active { methods: numerical { stars: black holes 1. INTRODUCTION (Chatzopoulos & Wheeler 2012; Yoon et al. 2012). Like- Recently, several gravitational wave (GW) events were wise, the probability that the progenitor of the 85 M reported by LIGO and Virgo whose measured physi- BH avoided a pulsation pair-instability supernova due cal properties pose interesting constraints on their as- to uncertainties on the reaction rate of carbon burning trophysical origin. One of these is GW190521 (The is low (Farmer et al. 2019). An alternative proposal is LIGO Scientific Collaboration et al. 2020b; Abbott for Pop III stellar remnants to grow via accretion (Sa- et al. 2020d), in which the masses of the merging farzadeh & Haiman 2020); this requires the remnant BHs +21 +17 to spend a significant time in a dense gaseous environ- BHs (85−14 M and 66−18 M , but see Nitz & Ca- pano 2020) indicate that the merging BHs fall into the ment of their host protogalaxies. so-called \upper-mass gap", as they are heavier than Similarly, the neutron star (NS) merger GW190425 in- the maximum BH mass imposed by (pulsation) pair- dicates masses higher than in NS binaries observed in instability supernovae (Chatzopoulos & Wheeler 2012; the Milky Way (Abbott et al. 2020a). To explain this Farmer et al. 2019). The components' dimensionless event, a rapid merging pathway is proposed for the bi- +0:27 +0:24 nary evolution channel (Romero-Shaw et al. 2020; Sa- spins (a = 0:69−0:62 and 0:73−0:64) are also larger than the maximum natal spins predicted by stellar evolution- farzadeh et al. 2020b; Galaudage et al. 2020). Addi- ary models for angular momentum transfer (Fuller & Ma tionally, mergers with very unequal masses have been 2019). A possible electromagnetic counterpart has been reported { GW190412 (The LIGO Scientific Collabora- reported for GW190521, which, if real, would suggest tion et al. 2020a) and GW190814 (Abbott et al. 2020b) {, which are also atypical in stellar evolutionary models arXiv:2012.00011v2 [astro-ph.HE] 20 Dec 2020 mergers in gas-rich environments (McKernan et al. 2019; Graham et al. 2020). Further candidates for mergers in of isolated binaries (Gerosa et al. 2020b; Olejak et al. the upper-mass gap are reported in GW170729 (Abbott 2020). The mass of the secondary object of 2:6 M in et al. 2019), GW170817A (Zackay et al. 2019a), and six GW190814 falls in the so-called lower mass gap, pre- other candidates in O3 (Abbott et al. 2020c). A possible dicted by stellar-evolutionary models (Fryer et al. 2012; explanation proposed for the large masses in these events Zevin et al. 2020 but see Belczynski et al. 2012; Drozda et al. 2020) and obeyed by observations of X-ray bina- is isolated binary evolution, presuming that massive BHs ¨ can form from Population III (Pop III) or metal-poor ries (Ozel et al. 2012). The object in the lower mass gap in GW190814 and a non-zero spin for the primary BH stars (Tanikawa et al. 2020; Kinugawa et al. 2020; Farrell +0:16 et al. 2020; Liu & Bromm 2020; Liu et al. 2020; Belczyn- (0:43−0:26) in GW190412 are consistent with a scenario ski 2020). However, due to their expected rapid rotation, in which the merging compact objects (COs) had expe- stellar remnants are unlikely to retain their hydrogen en- rienced previous mergers or significant accretion. These velopes and produce BHs as massive as in GW190521 events suggest that growth by gas accretion or hierar- chical mergers may be common among COs (see e.g. E-mail: [email protected] O'Leary et al. 2016; Gerosa et al. 2020b; Abbott et al. 2 Tagawa et al. 2020d; Hamers & Safarzadeh 2020; Rodriguez et al. 2020; three-body encounters. Highly eccentric CO binaries also Safarzadeh et al. 2020a; Safarzadeh & Haiman 2020; form due to the GW-capture mechanism in single-single Yang et al. 2020a; Liu & Lai 2020). encounters (e.g. O'Leary et al. 2009) and BS exchange By analyzing the ensemble of events detected during interactions (Paper III). LIGO/Virgo's O1-O3a observing runs, Kimball et al. Regarding the interaction with gas, the velocities of all (2020) and Tiwari & Fairhurst(2020) found preference COs relative to the local AGN disk decrease due to ac- for at least one, but probably multiple hierarchical merg- cretion torque and gas dynamical friction. The binaries' ers in the detected sample. The results by Kimball semi-major axis decreases due to gas dynamical friction et al.(2020) strongly depend on the escape velocity at from the AGN disk. Their semi-major axis and eccen- the merger sites, with higher escape velocities favoring a tricity evolve also due to type I/II migration torque from larger number of hierarchical mergers. a circumbinary disk. COs migrate toward the SMBH due Active galactic nuclei (AGNs) represent environments to type I/II torques from the AGN disk. CO masses, BH where mass growth by gas accretion is possible (Levin spins, and the orbital angular momentum directions of 2007) and hierarchical mergers are frequent due to effi- binaries gradually evolve due gas accretion (Paper II). cient hardening by interaction with gas (e.g. Bartos et al. We also account for dynamical interactions with single 2017; Stone et al. 2017; McKernan et al. 2018, 2020b,a), stars/COs and CO binaries. The binaries' semi-major stars, and COs (Tagawa et al. 2020b, hereafter Paper I). axis, velocities, orbital angular momentum directions, Additionally, the high escape velocity in such systems and eccentricity evolve due to BS interactions as pre- helps retain recoiling merger remnants. Yang et al. scribed in Leigh et al.(2018), and the velocities of all (2020b) investigated the possibility of lower mass gap COs additionally evolve due to scattering. The soft bi- mergers by assuming that the objects are quickly deliv- naries, i.e. those with binding energy lower than that ered to migration traps of AGN disks where they undergo of scattering objects, are either softened or ultimately mergers. In Paper I, we performed one-dimensional (1D) disrupted during BS interactions. Hard binaries become N-body simulations combined with semi-analytical pre- harder during BS interactions, and exchange interactions scriptions for gaseous and stellar effects, to follow the may take place; the most massive pair is assumed to re- evolution of binary separation, along with the radius and main in the binary after the interaction (Paper III). inclination of the orbit around the SMBH. The results The binaries' separation and eccentricity decrease due showed that binaries form efficiently during dynamical to GW emission. When the pericenter distance of the encounters in gas, and mergers are facilitated by binary- binary becomes smaller than the innermost stable orbit, single (BS) interactions in the gap-forming inner regions. we assume that BHs merge, assign a kick velocity and In Tagawa et al.(2020a, hereafter Paper II), we extended mass loss due to GW emission, and prescribe BH spin this model to follow the evolution of BH spins, and bi- evolution at the merger as in Paper II. nary orbital angular momentum directions. In Tagawa et al.(2020c, hereafter Paper III), we also included the 2.2. Evolution of neutron stars effects of binary eccentricity, and found that mergers are While only stellar-mass BHs are followed in Pa- often highly eccentric in AGN disks if BS interactions are pers I/II/III, here we additionally follow the evolution confined to the AGN plane due to gaseous torques (see of NSs among the N-body particles. We assume that also Samsing et al. 2020) or if the objects migrate close stars with zero-age main-sequence masses of 8{20 M to the SMBH where they undergo GW capture events.