-
Time-dependent models of AGN disks with radiation from embedded stellar-mass black holes
Authors:
Marguerite Epstein-Martin,
Hiromichi Tagawa,
Zoltan Haiman,
Rosalba Perna
Abstract:
The brightest steady sources of radiation in the universe, active galactic nuclei (AGN), are powered by gas accretion onto a central supermassive black hole (SMBH). The large sizes and accretion rates implicated in AGN accretion disks are expected to lead to gravitational instability and fragmentation, effectively cutting off mass inflow to the SMBH. Radiative feedback from disk-embedded stars has…
▽ More
The brightest steady sources of radiation in the universe, active galactic nuclei (AGN), are powered by gas accretion onto a central supermassive black hole (SMBH). The large sizes and accretion rates implicated in AGN accretion disks are expected to lead to gravitational instability and fragmentation, effectively cutting off mass inflow to the SMBH. Radiative feedback from disk-embedded stars has been invoked to yield marginally stable, steady-state solutions in the outer disks. Here, we examine the consequences of this star formation with a semi-analytical model in which stellar-mass black hole (sBH) remnants in the disk provide an additional source of stabilizing radiative feedback. Assuming star formation seeds the embedded sBH population, we model the time-evolving feedback from both stars and the growing population of accreting sBHs. We find that in the outer disk, the luminosity of the sBHs quickly dominates that of their parent stars. However, because sBHs consume less gas than stars to stabilize the disk, the presence of the sBHs enhances the mass flux to the inner disk. As a result, star formation persists over the lifetime of the AGN, damped in the outer disk, but amplified in a narrow ring in the inner disk. Heating from the embedded sBHs significantly modifies the disk's temperature profile and hardens its spectral energy distribution, and direct emission from the sBHs adds a new hard X ray component, resembling a Compton reflection "hump".
△ Less
Submitted 15 May, 2024;
originally announced May 2024.
-
Indication for a compact object next to a LIGO-Virgo binary black hole merger
Authors:
Wen-Biao Han,
Shu-Cheng Yang,
Hiromichi Tagawa,
Ye Jiang,
Ping Shen,
Qianyun Yun,
Chen Zhang,
Xing-Yu Zhong
Abstract:
The astrophysical origin of binary black hole (BBH) mergers remains uncertain though many events have been observed by the LIGO-Virgo-KAGRA network. Such mergers are predicted to originate in the vicinity of massive black holes (MBHs). Especially, GW190814, due to its secondary mass and mass ratio being beyond the expectations of isolated stellar evolution theories, is a promising event that has h…
▽ More
The astrophysical origin of binary black hole (BBH) mergers remains uncertain though many events have been observed by the LIGO-Virgo-KAGRA network. Such mergers are predicted to originate in the vicinity of massive black holes (MBHs). Especially, GW190814, due to its secondary mass and mass ratio being beyond the expectations of isolated stellar evolution theories, is a promising event that has happened in an active galactic nucleus(AGN) disk. In this model, a compact object resides in the vicinity of a merging BBH. Here we report multiple pieces of evidence pointing to the fact that GW190814 is a BBH merging near a compact object. The orbital motion of BBHs around the third body produces a line-of-sight acceleration (LSA) and induces a varying Doppler shift. Using a waveform template that considers LSA, we perform Bayesian inference on a few BBH events with a high signal-to-noise ratio in the gravitational-wave transient catalog (GWTC). Compared to the model for isolated BBH mergers, we obtain significantly higher network signal-to-noise ratios for GW190814 by that with the LSA and constrain the LSA to $a = 0.0014^{+0.0014}_{-0.0022} ~c~\mathrm{s}^{-1}$. In addition, the logarithmic Bayes factor for the LSA case over the isolated case is $16.6$, which means the LSA model is significantly preferred by the GW data. We conclude that this is the first indication showing merging BBHs are located near a compact object.
△ Less
Submitted 3 January, 2024;
originally announced January 2024.
-
Shock cooling and breakout emission for optical flares associated with gravitational wave events
Authors:
Hiromichi Tagawa,
Shigeo S. Kimura,
Zoltán Haiman,
Rosalba Perna,
Imre Bartos
Abstract:
The astrophysical origin of stellar-mass black hole (BH) mergers discovered through gravitational waves (GWs) is widely debated. Mergers in the disks of active galactic nuclei (AGN) represent promising environments for at least a fraction of these events, with possible observational clues in the GW data. An additional clue to unveil AGN merger environments is provided by possible electromagnetic e…
▽ More
The astrophysical origin of stellar-mass black hole (BH) mergers discovered through gravitational waves (GWs) is widely debated. Mergers in the disks of active galactic nuclei (AGN) represent promising environments for at least a fraction of these events, with possible observational clues in the GW data. An additional clue to unveil AGN merger environments is provided by possible electromagnetic emission from post-merger accreting BHs. Associated with BH mergers in AGN disks, emission from shocks emerging around jets launched by accreting merger remnants is expected. In this paper we compute the properties of the emission produced during breakout and the subsequent adiabatic expansion phase of the shocks, and we then apply this model to optical flares suggested to be possibly associated with GW events. We find that the majority of the reported flares can be explained by the breakout and the shock cooling emission. If these events are real, then the merging locations of binaries are constrained depending on the emission processes. If the optical flares are produced by shock cooling emission, they would display moderate color evolution, possibly color variations among different events, a positive correlation between the delay time and the duration of flares, and accompanying breakout emission in X-ray bands before the optical flares. If the breakout emission dominates the observed lightcurve, it is expected that the color is distributed in a narrow range in the optical band, and the delay time from GW to electromagnetic emission is longer than $\sim 2$ days. Hence, further explorations of the distributions of delay times, color evolution of the flares, and associated X-ray emission will be useful to test the proposed emission model for the observed flares.
△ Less
Submitted 27 October, 2023;
originally announced October 2023.
-
High-energy electromagnetic, neutrino, and cosmic-ray emission by stellar-mass black holes in disks of active galactic nuclei
Authors:
Hiromichi Tagawa,
Shigeo S. Kimura,
Zoltán Haiman
Abstract:
Some Seyfert galaxies are detected in high-energy gamma rays, but the mechanism and site of gamma-ray emission are unknown. Also, the origins of the cosmic high-energy neutrino and MeV gamma-ray backgrounds have been veiled in mystery since their discoveries. We propose emission from stellar-mass BHs (sBHs) embedded in disks of active galactic nuclei (AGN) as their possible sources. These sBHs are…
▽ More
Some Seyfert galaxies are detected in high-energy gamma rays, but the mechanism and site of gamma-ray emission are unknown. Also, the origins of the cosmic high-energy neutrino and MeV gamma-ray backgrounds have been veiled in mystery since their discoveries. We propose emission from stellar-mass BHs (sBHs) embedded in disks of active galactic nuclei (AGN) as their possible sources. These sBHs are predicted to launch jets due to the Blandford-Znajek mechanism, which can produce intense electromagnetic, neutrino, and cosmic-ray emissions. We investigate whether these emissions can be the sources of cosmic high-energy particles. We find that emission from internal shocks in the jets can explain gamma rays from nearby radio-quiet Seyfert galaxies including NGC1068, if the Lorentz factor of the jets ($Γ_{\rm j}$) is high. On the other hand, for moderate $Γ_{\rm j}$, the emission can significantly contribute to the background gamma-ray and neutrino intensities in the $\sim {\rm MeV}$ and $\lesssim {\rm PeV}$ bands, respectively. Furthermore, for moderate $Γ_{\rm j}$ with efficient amplification of the magnetic field and cosmic-ray acceleration, the neutrino emission from NGC1068 and the ultrahigh-energy cosmic rays can be explained. These results suggest that the neutrino flux from NGC1068 as well as the background intensities of ${\rm MeV}$ gamma rays, neutrinos, and the ultrahigh-energy cosmic rays can be explained by a unified model. Future MeV gamma-ray satellites will test our scenario for neutrino emission.
△ Less
Submitted 12 July, 2023;
originally announced July 2023.
-
Flares from stars crossing active galactic nuclei disks on low-inclination orbits
Authors:
Hiromichi Tagawa,
Zoltán Haiman
Abstract:
The origin of the recently discovered new class of transients, X-ray quasi-periodic eruptions (QPEs), remains a puzzle. Due to their periodicity and association with active galactic nuclei (AGN), it is natural to relate these eruptions to stars or compact objects in tight orbits around supermassive black holes (SMBHs). In this paper, we predict the properties of emission from bow shocks produced b…
▽ More
The origin of the recently discovered new class of transients, X-ray quasi-periodic eruptions (QPEs), remains a puzzle. Due to their periodicity and association with active galactic nuclei (AGN), it is natural to relate these eruptions to stars or compact objects in tight orbits around supermassive black holes (SMBHs). In this paper, we predict the properties of emission from bow shocks produced by stars crossing AGN disks, and compare them to the observed properties of QPEs. We find that when a star's orbit is retrograde and has a low inclination ($\lesssim 40^\circ$) with respect to the AGN disk and the star is massive ($\gtrsim 10$ $M_{\odot}$), the breakout emission from the bow shock can explain the observed duration ($\sim$ hours) and X-ray luminosity ($\sim$few$\times10^{42}~{\rm erg~s^{-1}}$) of QPEs. This model can further explain various observed features of QPEs, such as their complex luminosity evolution, the gradual decline of luminosity of the flares over several years, the evolution of the hardness ratio, the modulation of the luminosity during quiescent phases, and the preference of the central SMBHs to have low masses.
△ Less
Submitted 29 August, 2023; v1 submitted 7 April, 2023;
originally announced April 2023.
-
Observable signatures of stellar-mass black holes in active galactic nuclei
Authors:
Hiromichi Tagawa,
Shigeo S. Kimura,
Zoltán Haiman,
Rosalba Perna,
Imre Bartos
Abstract:
Stellar-mass black holes (BHs) are predicted to be embedded in the disks of active galactic nuclei (AGN) due to gravitational drag and in-situ star formation. However, clear evidence for AGN disk-embedded BHs is currently lacking. Here, as possible electromagnetic signatures of these BHs, we investigate breakout emission from shocks emerging around Blandford-Znajek jets launched from accreting BHs…
▽ More
Stellar-mass black holes (BHs) are predicted to be embedded in the disks of active galactic nuclei (AGN) due to gravitational drag and in-situ star formation. However, clear evidence for AGN disk-embedded BHs is currently lacking. Here, as possible electromagnetic signatures of these BHs, we investigate breakout emission from shocks emerging around Blandford-Znajek jets launched from accreting BHs in AGN disks. We assume that the majority of the highly super-Eddington flow reaches the BH, produces a strong jet, and the jet produces feedback that shuts off accretion and thus leads to episodic flaring. While these assumptions are highly uncertain at present, they predict a breakout emission characterized by luminous thermal emission in the X-ray bands, and bright, broadband non-thermal emission from the infrared to the gamma-ray bands. The flare duration depends on the BH's distance $r$ from the central supermassive BH, varying between $10^3-10^6$ s for $r \sim 0.01-1$ pc. This emission can be discovered by current and future infrared, optical, and X-ray wide-field surveys and monitoring campaigns of nearby AGNs.
△ Less
Submitted 3 March, 2023;
originally announced March 2023.
-
Observable signature of merging stellar-mass black holes in active galactic nuclei
Authors:
Hiromichi Tagawa,
Shigeo S. Kimura,
Zoltán Haiman,
Rosalba Perna,
Imre Bartos
Abstract:
The origin of stellar-mass black hole mergers discovered through gravitational waves is being widely debated. Mergers in the disks of active galactic nuclei (AGN) represent a promising source of origin, with possible observational clues in the gravitational wave data. Beyond gravitational waves, a unique signature of AGN-assisted mergers is electromagnetic emission from the accreting black holes.…
▽ More
The origin of stellar-mass black hole mergers discovered through gravitational waves is being widely debated. Mergers in the disks of active galactic nuclei (AGN) represent a promising source of origin, with possible observational clues in the gravitational wave data. Beyond gravitational waves, a unique signature of AGN-assisted mergers is electromagnetic emission from the accreting black holes. Here we show that jets launched by accreting black holes merging in an AGN disk can be detected as peculiar transients by infrared, optical, and X-ray observatories We further show that this emission mechanism can explain the possible associations between gravitational wave events and the optical transient ZTF19abanrhr and the proposed gamma-ray counterparts GW150914-GBM and LVT151012-GBM. We demonstrate how these associations, if genuine, can be used to reconstruct the properties of these events' environments. Searching for infrared and X-ray counterparts to similar electromagnetic transients in the future, once host galaxies are localized by optical observations, could provide a smoking gun signature of the mergers' AGN origin.
△ Less
Submitted 18 March, 2023; v1 submitted 17 January, 2023;
originally announced January 2023.
-
Gravitational wave source populations: Disentangling an AGN component
Authors:
V. Gayathri,
Daniel Wysocki,
Y. Yang,
R. O Shaughnessy,
Z. Haiman,
H. Tagawa,
I. Bartos
Abstract:
The astrophysical origin of the over 90 compact binary mergers discovered by the LIGO and Virgo gravitational wave observatories is an open question. While the unusual mass and spin of some of the discovered objects constrain progenitor scenarios, the observed mergers are consistent with multiple interpretations. A promising approach to solve this question is to consider the observed distributions…
▽ More
The astrophysical origin of the over 90 compact binary mergers discovered by the LIGO and Virgo gravitational wave observatories is an open question. While the unusual mass and spin of some of the discovered objects constrain progenitor scenarios, the observed mergers are consistent with multiple interpretations. A promising approach to solve this question is to consider the observed distributions of binary properties and compare them to expectations from different origin scenarios. Here we describe a new hierarchical population analysis framework to assess the relative contribution of different formation channels simultaneously. For this study we considered binary formation in AGN disks along with phenomenological models, but the same framework can be extended to other models. We find that high-mass and high-mass-ratio binaries appear more likely to have an AGN origin compared to the same origin as lower-mass events. Future observations of high-mass black hole mergers could further disentangle the AGN component from other channels.
△ Less
Submitted 10 January, 2023;
originally announced January 2023.
-
Can stellar-mass black hole growth disrupt disks of active galactic nuclei? The role of mechanical feedback
Authors:
Hiromichi Tagawa,
Shigeo S. Kimura,
Zoltán Haiman,
Rosalba Perna,
Hidekazu Tanaka,
Imre Bartos
Abstract:
Stellar-mass BHs (sBHs) are predicted to be embedded in active galactic nuclei (AGN) disks due to gravitational drag and in-situ star formation. However, we find that due to a high gas density in an AGN disk environment, compact objects may rapidly grow to intermediate-mass BHs and deplete matter from the AGN disk unless accretion is suppressed by some feedback process(es). These consequences are…
▽ More
Stellar-mass BHs (sBHs) are predicted to be embedded in active galactic nuclei (AGN) disks due to gravitational drag and in-situ star formation. However, we find that due to a high gas density in an AGN disk environment, compact objects may rapidly grow to intermediate-mass BHs and deplete matter from the AGN disk unless accretion is suppressed by some feedback process(es). These consequences are inconsistent with AGN observations and the dynamics of the Galactic center. Here we consider mechanical feedback mechanisms for the reduction of gas accretion. Rapidly accreting sBHs launch winds and/or jets via the Blandford-Znajek mechanism, which produce high-pressure shocks and cocoons. Such a shock and cocoon can spread laterally in the plane of the disk, eject the outer regions of a circum-sBH disk (CsBD) and puncture a hole in the AGN disk with horizontal size comparable to the disk scale-height. Since the depletion timescale of the bound CsBD is much shorter than the resupply timescale of gas to the sBH, the time-averaged accretion rate onto sBHs is reduced by this process by a factor of $\sim 10$--$100$. This feedback mechanism can therefore help alleviate the sBH over-growth and AGN-disk depletion problems. On the other hand, we find that cocoons of jets can unbind a large fraction of the gas accreting in the disks of less massive SMBHs, which may help explain the dearth of high-Eddington ratio AGNs with SMBH mass $\lesssim10^5{\rm M_\odot}$.
△ Less
Submitted 2 December, 2021;
originally announced December 2021.
-
Micro Tidal Disruption Events in Active Galactic Nuclei
Authors:
Y. Yang,
I. Bartos,
G. Fragione,
Z. Haiman,
M. Kowalski,
S. Marka,
R. Perna,
H. Tagawa
Abstract:
Active galactic nuclei (AGNs) can act as black hole assembly lines, funneling some of the stellar-mass black holes from the vicinity of the galactic center into the inner plane of the AGN disk where the black holes can merge through dynamical friction and gravitational wave emission. Here, we show that stars near the galactic center are also brought into the AGN disk, where they can be tidally dis…
▽ More
Active galactic nuclei (AGNs) can act as black hole assembly lines, funneling some of the stellar-mass black holes from the vicinity of the galactic center into the inner plane of the AGN disk where the black holes can merge through dynamical friction and gravitational wave emission. Here, we show that stars near the galactic center are also brought into the AGN disk, where they can be tidally disrupted by the stellar-mass black holes in the disk. Such micro-tidal disruption events (micro-TDEs) could be useful probe of stellar interaction with the AGN disk. We find that micro-TDEs in AGNs occur at a rate of $\sim170$ Gpc$^{-3}$yr$^{-1}$. Their cleanest observational probe may be the detection of tidal disruption in AGNs by heavy supermassive black holes ($M_{\bullet}\gtrsim10^{8}$ M$_{\odot}$) so that cannot tidally disrupt solar-type stars. We discuss two such TDE candidates observed to date (ASASSN-15lh and ZTF19aailpwl).
△ Less
Submitted 5 May, 2021;
originally announced May 2021.
-
Signatures of Hierarchical Mergers in Black Hole Spin and Mass distribution
Authors:
Hiromichi Tagawa,
Zoltán Haiman,
Imre Bartos,
Bence Kocsis,
Kazuyuki Omukai
Abstract:
Recent gravitational wave (GW) observations by LIGO/Virgo show evidence for hierarchical mergers, where the merging BHs are the remnants of previous BH merger events. These events may carry important clues about the astrophysical host environments of the GW sources. In this paper, we present the distributions of the effective spin parameter ($χ_\mathrm{eff}$), the precession spin parameter (…
▽ More
Recent gravitational wave (GW) observations by LIGO/Virgo show evidence for hierarchical mergers, where the merging BHs are the remnants of previous BH merger events. These events may carry important clues about the astrophysical host environments of the GW sources. In this paper, we present the distributions of the effective spin parameter ($χ_\mathrm{eff}$), the precession spin parameter ($χ_\mathrm{p}$), and the chirp mass ($m_\mathrm{chirp}$) expected in hierarchical mergers. Under a wide range of assumptions, hierarchical mergers produce (i) a monotonic increase of the average of the typical total spin for merging binaries, which we characterize with ${\bar χ}_\mathrm{typ}\equiv \overline{(χ_\mathrm{eff}^2+χ_\mathrm{p}^2)^{1/2}}$, up to roughly the maximum $m_\mathrm{chirp}$ among first-generation (1g) BHs, and (ii) a plateau at ${\bar χ}_\mathrm{typ}\sim 0.6$ at higher $m_\mathrm{chirp}$. We suggest that the maximum mass and typical spin magnitudes for 1g BHs can be estimated from ${\bar χ}_\mathrm{typ}$ as a function of $m_\mathrm{chirp}$. The GW data observed in LIGO/Virgo O1--O3a prefers an increase in ${\bar χ}_\mathrm{typ}$ at low $m_\mathrm{chirp}$, which is consistent with the growth of the BH spin magnitude by hierarchical mergers, at $\sim 2 σ$ confidence. A Bayesian analysis suggests that 1g BHs have the maximum mass of $\sim 15$--$30\,M_\odot$ if the majority of mergers are of high-generation BHs (not among 1g-1g BHs), which is consistent with mergers in active galactic nucleus disks and/or nuclear star clusters, while if mergers mainly originate from globular clusters, 1g BHs are favored to have non-zero spin magnitudes of $\sim 0.3$. We also forecast that signatures for hierarchical mergers in the ${\bar χ}_\mathrm{typ}$ distribution can be confidently recovered once the number of GW events increases to $\gtrsim O(100)$.
△ Less
Submitted 10 August, 2021; v1 submitted 19 April, 2021;
originally announced April 2021.
-
Accretion-Induced Collapse of Neutron Stars in the Disks of Active Galactic Nuclei
Authors:
Rosalba Perna,
Hiromichi Tagawa,
Zoltan Haiman,
Imre Bartos
Abstract:
The disks of active galactic nuclei (AGNs) have emerged as a rich environment for the evolution of stars and their compact remnants. The very dense medium favors rapid accretion, while torques and migration traps enhance binary formation and mergers. Both long and short gamma-ray bursts (GRBs) are hence expected. We show that AGN disks constitute an ideal environment for another interesting phenom…
▽ More
The disks of active galactic nuclei (AGNs) have emerged as a rich environment for the evolution of stars and their compact remnants. The very dense medium favors rapid accretion, while torques and migration traps enhance binary formation and mergers. Both long and short gamma-ray bursts (GRBs) are hence expected. We show that AGN disks constitute an ideal environment for another interesting phenomenon: the accretion induced collapse (AIC) of neutron stars (NSs) to black holes (BHs). Rapid accretion in the dense disks can cause NSs to grow to the point of exceeding the maximum mass allowed by their equation of state. General relativistic magnetohydrodynamical simulations have shown that electromagnetic signatures are expected if the NS is surrounded by a mini-disk prior to collapse, which then rapidly accretes onto the BH, and/or if the NS is highly magnetized, from reconnection of the magnetosphere during collapse. Here we compute the rates of AICs and their locations within the disks for both isolated NSs, and for (initially stable) NSs formed from NS-NS mergers. We find that the global AIC rates are $\sim 0.07-20$~Gpc$^{-3}$~yr$^{-1}$, and we discuss their observable prospects and signatures as they emerge from the dense disk environments.
△ Less
Submitted 19 March, 2021;
originally announced March 2021.
-
Mass-gap Mergers in Active Galactic Nuclei
Authors:
Hiromichi Tagawa,
Bence Kocsis,
Zoltan Haiman,
Imre Bartos,
Kazuyuki Omukai,
Johan Samsing
Abstract:
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 efficien…
▽ More
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.
△ Less
Submitted 20 December, 2020; v1 submitted 30 November, 2020;
originally announced December 2020.
-
Eccentric Black Hole Mergers in Active Galactic Nuclei
Authors:
Hiromichi Tagawa,
Bence Kocsis,
Zoltan Haiman,
Imre Bartos,
Kazuyuki Omukai,
Johan Samsing
Abstract:
The astrophysical origin of gravitational wave (GW) transients is a timely open question in the wake of discoveries by LIGO/Virgo. In active galactic nuclei (AGNs), binaries form and evolve efficiently by interaction with a dense population of stars and the gaseous AGN disk. Previous studies have shown that stellar-mass black hole (BH) mergers in such environments can explain the merger rate and t…
▽ More
The astrophysical origin of gravitational wave (GW) transients is a timely open question in the wake of discoveries by LIGO/Virgo. In active galactic nuclei (AGNs), binaries form and evolve efficiently by interaction with a dense population of stars and the gaseous AGN disk. Previous studies have shown that stellar-mass black hole (BH) mergers in such environments can explain the merger rate and the number of suspected hierarchical mergers observed by LIGO/Virgo. The binary eccentricity distribution can provide further information to distinguish between astrophysical models. Here we derive the eccentricity distribution of BH mergers in AGN disks. We find that eccentricity is mainly due to binary-single (BS) interactions, which lead to most BH mergers in AGN disks having a significant eccentricity at $0.01\,\mathrm{Hz}$, detectable by LISA. If BS interactions occur in isotropic-3D directions, then $8$--$30\%$ of the mergers in AGN disks will have eccentricities at $10\,\mathrm{Hz}$ above $e_{10\,\rm Hz}\gtrsim 0.03$, detectable by LIGO/Virgo/KAGRA, while $5$--$17\%$ of mergers have $e_{10\,\rm Hz}\geq 0.3$. On the other hand, if BS interactions are confined to the AGN-disk plane due to torques from the disk, with 1-20 intermediate binary states during each interaction, or if BHs can migrate to $\lesssim10^{-3}\,\mathrm{pc}$ from the central supermassive black hole, then $10$--$70\%$ of the mergers will be highly eccentric ($e_{10\,\rm Hz} \geq 0.3$), consistent with the possible high eccentricity in GW190521.
△ Less
Submitted 18 December, 2020; v1 submitted 20 October, 2020;
originally announced October 2020.
-
AGN as Potential Factories for Eccentric Black Hole Mergers
Authors:
J. Samsing,
I. Bartos,
D. J. D'Orazio,
Z. Haiman,
B. Kocsis,
N. W. C. Leigh,
B. Liu,
M. E. Pessah,
H. Tagawa
Abstract:
There is some weak evidence that the black hole merger named GW190521 had a non-zero eccentricity. In addition, the component black holes' masses exceeded the limit predicted by stellar evolution. The large masses can be explained by successive mergers, which may be efficient in gas disks surrounding active galactic nuclei (AGN), but it is difficult to maintain an eccentric orbit all the way to th…
▽ More
There is some weak evidence that the black hole merger named GW190521 had a non-zero eccentricity. In addition, the component black holes' masses exceeded the limit predicted by stellar evolution. The large masses can be explained by successive mergers, which may be efficient in gas disks surrounding active galactic nuclei (AGN), but it is difficult to maintain an eccentric orbit all the way to the merger, as basic physics would argue for circularization. Here we show that AGN-disk environments can lead to an excess of eccentric mergers, if the interactions between single and binary black holes are frequent, and occur with mutual inclinations of less than a few degrees. We further illustrate that this eccentric population has a different distribution of the inclination between the spin vectors of the black holes and their orbital angular momentum at merger, referred to as the spin-orbit tilt, compared to the remaining circular mergers.
△ Less
Submitted 21 April, 2022; v1 submitted 19 October, 2020;
originally announced October 2020.
-
Black Hole Formation in the Lower Mass Gap through Mergers and Accretion in AGN Disks
Authors:
Y. Yang,
V. Gayathri,
I. Bartos,
Z. Haiman,
M. Safarzadeh,
H. Tagawa
Abstract:
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 gravita…
▽ More
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.
△ Less
Submitted 9 July, 2020;
originally announced July 2020.
-
Spin Evolution of Stellar-mass Black Hole Binaries in Active Galactic Nuclei
Authors:
Hiromichi Tagawa,
Zoltan Haiman,
Imre Bartos,
Bence Kocsis
Abstract:
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 i…
▽ More
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 $χ_\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 $χ_\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, $χ_\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 $χ_\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.
△ Less
Submitted 4 July, 2020; v1 submitted 24 April, 2020;
originally announced April 2020.
-
Cosmic Evolution of Stellar-mass Black Hole Merger Rate in Active Galactic Nuclei
Authors:
Y. Yang,
I. Bartos,
Z. Haiman,
B. Kocsis,
S. Márka,
H. Tagawa
Abstract:
Binary black hole mergers encode information about their environment and the astrophysical processes that led to their formation. Measuring the redshift dependence of their merger rate will help probe the formation and evolution of galaxies and the evolution of the star formation rate. Here we compute the cosmic evolution of the merger rate for stellar-mass binaries in the disks of Active Galactic…
▽ More
Binary black hole mergers encode information about their environment and the astrophysical processes that led to their formation. Measuring the redshift dependence of their merger rate will help probe the formation and evolution of galaxies and the evolution of the star formation rate. Here we compute the cosmic evolution of the merger rate for stellar-mass binaries in the disks of Active Galactic Nuclei (AGNs). We focus on recent evolution out to redshift $z=2$, covering the accessible range of current Earth-based gravitational-wave observatories. On this scale, the AGN population density is the main contributor to redshift-dependence. We find that the AGN-assisted merger rate does not meaningfully evolve with redshift, differentiating this channel from field binaries and some other dynamical formation scenarios.
△ Less
Submitted 19 March, 2020;
originally announced March 2020.
-
Formation and Evolution of Compact Object Binaries in AGN Disks
Authors:
Hiromichi Tagawa,
Zoltan Haiman,
Bence Kocsis
Abstract:
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…
▽ More
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.
△ Less
Submitted 29 May, 2020; v1 submitted 17 December, 2019;
originally announced December 2019.
-
Making a supermassive star by stellar bombardment
Authors:
Hiromichi Tagawa,
Zoltan Haiman,
Bence Kocsis
Abstract:
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 s…
▽ More
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.
△ Less
Submitted 26 March, 2020; v1 submitted 23 September, 2019;
originally announced September 2019.
-
Hypervelocity Stars from a Supermassive Black Hole-Intermediate-mass Black Hole binary
Authors:
Alexander Rasskazov,
Giacomo Fragione,
Nathan W. C. Leigh,
Hiromichi Tagawa,
Alberto Sesana,
Adrian Price-Whelan,
Elena Maria Rossi
Abstract:
In this paper we consider a scenario where the currently observed hypervelocity stars in our Galaxy have been ejected from the Galactic center as a result of dynamical interactions with an intermediate-mass black hole (IMBH) orbiting the central supermassive black hole (SMBH). By performing 3-body scattering experiments, we calculate the distribution of the ejected stars' velocities given various…
▽ More
In this paper we consider a scenario where the currently observed hypervelocity stars in our Galaxy have been ejected from the Galactic center as a result of dynamical interactions with an intermediate-mass black hole (IMBH) orbiting the central supermassive black hole (SMBH). By performing 3-body scattering experiments, we calculate the distribution of the ejected stars' velocities given various parameters of the IMBH-SMBH binary: IMBH mass, semimajor axis and eccentricity. We also calculate the rates of change of the BH binary orbital elements due to those stellar ejections. One of our new findings is that the ejection rate depends (although mildly) on the rotation of the stellar nucleus (its total angular momentum). We also compare the ejection velocity distribution with that produced by the Hills mechanism (stellar binary disruption) and find that the latter produces faster stars on average. Also, the IMBH mechanism produces an ejection velocity distribution which is flattened towards the BH binary plane while the Hills mechanism produces a spherically symmetric one. The results of this paper will allow us in the future to model the ejection of stars by an evolving BH binary and compare both models with \textit{Gaia} observations, for a wide variety of environments (galactic nuclei, globular clusters, the Large Magellanic Clouds, etc.).
△ Less
Submitted 23 May, 2019; v1 submitted 29 October, 2018;
originally announced October 2018.
-
Metallicity dependence of the Hercules stream in Gaia/RAVE data -- explanation by non-closed orbits
Authors:
Kohei Hattori,
Naoteru Gouda,
Taihei Yano,
Nobuyuki Sakai,
Hiromichi Tagawa,
Junichi Baba,
Jun Kumamoto
Abstract:
The origin of the Hercules stream, the most prominent velocity substructure in the Solar neighbour disc stars, is still under debate. Recent accurate measurements of position, velocity, and metallicity provided by Tycho Gaia Astrometric Solution (TGAS) and RAdial Velocity Experiments (RAVE) have revealed that the Hercules stream is most clearly seen in the metal-rich region ([Fe/H] > 0), while it…
▽ More
The origin of the Hercules stream, the most prominent velocity substructure in the Solar neighbour disc stars, is still under debate. Recent accurate measurements of position, velocity, and metallicity provided by Tycho Gaia Astrometric Solution (TGAS) and RAdial Velocity Experiments (RAVE) have revealed that the Hercules stream is most clearly seen in the metal-rich region ([Fe/H] > 0), while it is not clearly seen in lower metallicity region ([Fe/H] < -0.25). By using a large number of chemo-dynamical 2D test-particle simulations with a rotating bar and/or spiral arms, we find that the observed [Fe/H] dependence of the Hercules stream is a natural consequence of the inside-out formation of the stellar disc and the existence of highly non-closed orbits in the rotating frame of the bar or spiral arms. Our successful models that reproduce the observed properties of the Hercules stream include not only fast-bar-only and fast-bar+spiral models, but also slow-bar+spiral models. This indicates that it is very difficult to estimate the pattern speed of the bar or spiral arms based only on the observations of the Hercules stream in the Solar neighbourhood. As a by-product of our simulations, we make some predictions about the locations across the Galactic plane where we can observe velocity bimodality that is not associated with the Hercules stream. These predictions can be tested by the Gaia Data Release 2, and such a test will improve our understanding of the evolution of the Milky Way stellar disc.
△ Less
Submitted 22 January, 2019; v1 submitted 5 April, 2018;
originally announced April 2018.
-
Merger of multiple accreting black holes concordant with gravitational wave events
Authors:
Hiromichi Tagawa,
Masayuki Umemura
Abstract:
Recently, advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO) has detected black hole (BH) merger events, most of which are sourced by BHs more massive than $30~M_\odot$. Especially, the observation of GW170104 suggests dynamically assembled binaries favoring a distribution of misaligned spins. It has been argued that mergers of unassociated BHs can be engendered through a "chance…
▽ More
Recently, advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO) has detected black hole (BH) merger events, most of which are sourced by BHs more massive than $30~M_\odot$. Especially, the observation of GW170104 suggests dynamically assembled binaries favoring a distribution of misaligned spins. It has been argued that mergers of unassociated BHs can be engendered through a "chance meeting" in a multiple BH system under gas-rich environments. In this paper, we consider the merger of unassociated BHs, concordant with the massive BH merger events. To that end, we simulate a multiple BH system with a post-Newtonian $N$-body code incorporating gas accretion and general relativistic effects. As a result, we find that gas dynamical friction effectively promotes three-body interaction of BHs in dense gas of $n_\mathrm{gas}\gtrsim 10^6\,\mathrm{cm}^{-3}$, so that BH mergers can take place within $30$ Myr. This scenario predicts an isotropic distribution of spin tilts. In the concordant models with GW150914, the masses of seed BHs are required to be $\gtrsim 25~M_\odot$. The potential sites of such "chance meeting" BH mergers are active galactic nucleus (AGN) disks and dense interstellar clouds. Assuming the LIGO O1, we roughly estimate the event rates for PopI BHs and PopIII BHs in AGN disks to be respectively $\simeq 1-2\,\mathrm{yr}^{-1}$ and $\simeq 1\,\mathrm{yr}^{-1}$. Multiple episodes of AGNs may enhance the rates by roughly an order of magnitude. For massive PopI BHs in dense interstellar clouds, the rate is $\simeq 0.02\,\mathrm{yr}^{-1}$. Hence, high-density AGN disks are a more plausible site for mergers of chance meeting BHs.
△ Less
Submitted 21 February, 2018;
originally announced February 2018.
-
Compact object mergers driven by gas fallback
Authors:
Hiromichi Tagawa,
Takayuki R. Saitoh,
Bence Kocsis
Abstract:
Recently several gravitational wave detections have shown evidence for compact object mergers. However, the astrophysical origin of merging binaries is not well understood. Stellar binaries are typically at much larger separations than what is needed for the binaries to merge due to gravitational wave emission, which leads to the so-called final AU problem. In this letter we propose a new channel…
▽ More
Recently several gravitational wave detections have shown evidence for compact object mergers. However, the astrophysical origin of merging binaries is not well understood. Stellar binaries are typically at much larger separations than what is needed for the binaries to merge due to gravitational wave emission, which leads to the so-called final AU problem. In this letter we propose a new channel for mergers of compact object binaries which solves the final AU problem. We examine the binary evolution following gas expansion due to a weak failed supernova explosion, neutrino mass loss, core disturbance, or envelope instability. In such situations the binary is possibly hardened by ambient gas. We investigate the evolution of the binary system after a shock has propagated by performing smoothed particle hydrodynamics simulations. We find that significant binary hardening occurs when the gas mass bound to the binary exceeds that of the compact objects. This mechanism represents a new possibility for the pathway to mergers for gravitational wave events.
△ Less
Submitted 26 June, 2018; v1 submitted 1 February, 2018;
originally announced February 2018.
-
The validation of made-to-measure method for reconstruction of phase space distribution functions
Authors:
Hiromichi Tagawa,
Naoteru Gouda,
Taihei Yano,
Takuji Hara
Abstract:
We investigate how accurately phase space distribution functions (DFs) in galactic models can be reconstructed by a made-to-measure (M2M) method, which constructs $N$-particle models of stellar systems from photometric and various kinematic data. The advantage of the M2M method is that this method can be applied to various galactic models without assumption of the spatial symmetries of gravitation…
▽ More
We investigate how accurately phase space distribution functions (DFs) in galactic models can be reconstructed by a made-to-measure (M2M) method, which constructs $N$-particle models of stellar systems from photometric and various kinematic data. The advantage of the M2M method is that this method can be applied to various galactic models without assumption of the spatial symmetries of gravitational potentials adopted in galactic models, and furthermore, numerical calculations of the orbits of the stars cannot be severely constrained by the capacities of computer memories. The M2M method has been applied to various galactic models. However, the degree of accuracy for the recovery of DFs derived by the M2M method in galactic models has never been investigated carefully. Therefore, we show the degree of accuracy for the recovery of the DFs for the anisotropic Plummer model and the axisymmetric Stäckel model, which have analytic solutions of the DFs. Furthermore, this study provides the dependence of the degree of accuracy for the recovery of the DFs on various parameters and a procedure adopted in this paper. As a result, we find that the degree of accuracy for the recovery of the DFs derived by the M2M method for the spherical target model is a few percent, and more than ten percent for the axisymmetric target model.
△ Less
Submitted 15 August, 2016;
originally announced August 2016.
-
Mergers of accreting stellar-mass black holes
Authors:
Hiromichi Tagawa,
Masayuki Umemura,
Naoteru Gouda
Abstract:
We present post-Newtonian $N$-body simulations on mergers of accreting stellar-mass black holes (BHs), where such general relativistic effects as the pericenter shift and gravitational wave (GW) emission are taken into consideration. The attention is concentrated on the effects of the dynamical friction and the Hoyle-Lyttleton mass accretion by ambient gas. We consider a system composed of ten BHs…
▽ More
We present post-Newtonian $N$-body simulations on mergers of accreting stellar-mass black holes (BHs), where such general relativistic effects as the pericenter shift and gravitational wave (GW) emission are taken into consideration. The attention is concentrated on the effects of the dynamical friction and the Hoyle-Lyttleton mass accretion by ambient gas. We consider a system composed of ten BHs with initial mass of $30~M_\odot$. As a result, we show that mergers of accreting stellar-mass BHs are classified into four types: a gas drag-driven, an interplay-driven, a three body-driven, or an accretion-driven merger. We find that BH mergers proceed before significant mass accretion, even if the accretion rate is $\sim10$ Eddington accretion rate, and then all BHs can merge into one heavy BH. Using the simulation results for a wide range of parameters, we derive a critical accretion rate ($\dot{m}_{\rm c}$), below which the BH growth is promoted faster by mergers. Also, it is found that the effect of the recoil by the GW emission can reduce $\dot{m}_{\rm c}$ especially in gas number density higher than $10^8~{\rm cm}^{-3}$, and enhance the escape probability of merged BHs. Very recently, a gravitational wave event, GW150914, as a result of the merger of a $\sim 30~M_\odot$ BH binary has been detected (Abbott et al. 2016). Based on the present simulations, the BH merger in GW150914 is likely to be driven by three-body encounters accompanied by a few $M_\odot$ of gas accretion, in high-density environments like dense interstellar clouds or galactic nuclei.
△ Less
Submitted 28 July, 2016; v1 submitted 28 February, 2016;
originally announced February 2016.
-
Early Cosmic Merger of Multiple Black Holes
Authors:
Hiromichi Tagawa,
Masayuki Umemura,
Naoteru Gouda,
Taihei Yano,
Yuki Yamai
Abstract:
We perform numerical simulations on the merger of multiple black holes (BHs) in primordial gas at early cosmic epochs. We consider two cases of BH mass: $M_{BH} = 30 M_{\odot}$ and $M_{BH} = 10^4 M_{\odot}$. Attention is concentrated on the effect of the dynamical friction by gas in a host object. The simulations incorporate such general relativistic effects as the pericentre shift and gravitation…
▽ More
We perform numerical simulations on the merger of multiple black holes (BHs) in primordial gas at early cosmic epochs. We consider two cases of BH mass: $M_{BH} = 30 M_{\odot}$ and $M_{BH} = 10^4 M_{\odot}$. Attention is concentrated on the effect of the dynamical friction by gas in a host object. The simulations incorporate such general relativistic effects as the pericentre shift and gravitational wave emission. As a result, we find that multiple BHs are able to merge into one BH within 100 Myr in a wide range of BH density. The merger mechanism is revealed to be categorized into three types: gas-drag-driven merger (type A), interplay-driven merger (type B), and three-body-driven merger (type C). We find the relation between the merger mechanism and the ratio of the gas mass within the initial BH orbit ($M_{gas}$) to the total BH mass ($ΣM_{BH}$). Type A merger occurs if $M_{gas} \gtrsim 10^5 ΣM_{BH}$, type B if $M_{gas} \lesssim 10^5 ΣM_{BH}$, and type C if $M_{gas} \ll 10^5 ΣM_{BH}$. Supposing the gas and BH density based on the recent numerical simulations on first stars, all the BH remnants from first stars are likely to merge into one BH through the type B or C mechanism. Also, we find that multiple massive BHs ($M_{BH} = 10^4 M_{\odot}$) distributed over several parsec can merge into one BH through the type B mechanism, if the gas density is higher than $5\times 10^6$ cm$^{-3}$. The present results imply that the BH merger may contribute significantly to the formation of supermassive BHs at high redshift epochs.
△ Less
Submitted 12 May, 2015;
originally announced May 2015.