Showing 1–50 of 354 results for author: Haiman, Z

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.

Comments: 21 pages, 15 figures

Relativistic Binary Precession: Impact on Eccentric Binary Accretion and Multi-Messenger Astronomy

Authors: Stanislav DeLaurentiis, Zoltan Haiman, John Ryan Westernacher-Schneider, Luke Major Krauth, Jordy Davelaar, Jonathan Zrake, Andrew MacFadyen

Abstract: Recent hydrodynamical simulations have shown that circumbinary gas disks drive the orbits of binary black holes to become eccentric, even when general relativistic corrections to the orbit are significant. Here, we study the general relativistic (GR) apsidal precession of eccentric equal-mass binary black holes in circumbinary disks (CBDs) via two-dimensional hydrodynamical simulations. We perform… ▽ More Recent hydrodynamical simulations have shown that circumbinary gas disks drive the orbits of binary black holes to become eccentric, even when general relativistic corrections to the orbit are significant. Here, we study the general relativistic (GR) apsidal precession of eccentric equal-mass binary black holes in circumbinary disks (CBDs) via two-dimensional hydrodynamical simulations. We perform a suite of simulations comparing precessing and non-precessing binaries across a range of eccentricities, semi-major axes, and precession rates. We find that the GR precession of the binary's semi-major axis can introduce a dominant modulation in the binary's accretion rate and the corresponding high-energy electromagnetic light-curves. We discuss the conditions under which this occurs and its detailed characteristics and mechanism. Finally, we discuss the potential to observe these precession signatures in electromagnetic and gravitational wave (GW) observations, as well as the precession signal's unique importance as a potential tool to constrain the mass, eccentricity, and semi-major axis of binary merger events. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 19 pages, 17 figures. Submitted to ApJ

Detecting environmental effects in gravitational waves from binaries perturbed by periodic forces

Authors: Lorenz Zwick, Christopher Tiede, Alessandro A. Trani, Andrea Derdzinski, Zoltan Haiman, Daniel J. D'Orazio, Johan Samsing

Abstract: We study the gravitational wave (GW) emission of sources perturbed by periodic dynamical forces which do not cause secular evolution in the orbital elements. We construct a corresponding post-Newtonian waveform model and provide estimates for the detectability of the resulting GW phase perturbations, for both space-based and future ground-based detectors. We validate our results by performing a se… ▽ More We study the gravitational wave (GW) emission of sources perturbed by periodic dynamical forces which do not cause secular evolution in the orbital elements. We construct a corresponding post-Newtonian waveform model and provide estimates for the detectability of the resulting GW phase perturbations, for both space-based and future ground-based detectors. We validate our results by performing a set of Bayesian parameter recovery experiments with post-Newtonian waveforms. We find that, in stark contrast to the more commonly studied secular dephasing, periodic phase perturbations do not suffer from degeneracies with any of the tested vacuum binary parameters. We discuss the applications of our findings to a range of possible astrophysical scenarios, finding that such periodic perturbations may be detectable for massive black hole binaries embedded in circum-binary discs, extreme mass-ratio inspirals in accretion discs, as well as stellar-mass compact objects perturbed by tidal fields. We argue that modelling conservative sub-orbital dynamics opens up a promising new avenue to detect environmental effects in binary sources of GWs that should be included in state-of-the-art waveform templates. △ Less

Submitted 15 May, 2024; v1 submitted 9 May, 2024; originally announced May 2024.

Comments: Submitted to PRD. Comments welcome

Identifying the electromagnetic counterparts of LISA massive black hole binaries in archival LSST data

Authors: Chengcheng Xin, Zoltan Haiman

Abstract: The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will catalogue the light-curves of up to 100 million quasars. Among these there can be up to approximately 100 ultra-compact massive black hole (MBH) binaries, which 5-15 years later can be detected in gravitational waves (GWs) by the Laser Interferometer Space Antenna (LISA). Here we assume that GWs from a MBH binary have been… ▽ More The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will catalogue the light-curves of up to 100 million quasars. Among these there can be up to approximately 100 ultra-compact massive black hole (MBH) binaries, which 5-15 years later can be detected in gravitational waves (GWs) by the Laser Interferometer Space Antenna (LISA). Here we assume that GWs from a MBH binary have been detected by LISA, and we assess whether or not its electromagnetic (EM) counterpart can be uniquely identified in archival LSST data as a periodic quasar. We use the binary's properties derived from the LISA waveform, such as the past evolution of its orbital frequency, its total mass, distance and sky localization, to predict the redshift, magnitude and historical periodicity of the quasar expected in the archival LSST data. We then use Monte Carlo simulations to compute the false alarm probability, i.e. the number of quasars in the LSST catalogue matching these properties by chance, based on the (extrapolated) quasar luminosity function, the sampling cadence of LSST, and intrinsic ``damped random walk (DRW)" quasar variability. We perform our analysis on four fiducial LISA binaries, with total masses and redshifts of $(M_{\rm bin}/{\rm M_{\odot}},z) = (3\times10^5,0.3)$, $(3\times10^6,0.3)$, $(10^7,0.3)$ and $(10^7,1)$. While DRW noise and aliasing due to LSST's cadence can produce false periodicities by chance, we find that the frequency chirp of the LISA source during the LSST observations washes out these noise peaks and allows the genuine source to stand out in Lomb-Scargle periodograms. We find that all four fiducial binaries yield excellent chances to be uniquely identified, with false alarm probabilities below $10^{-5}$, a week or more before their merger. This then enables deep follow-up EM observations targeting the individual EM counterparts during their inspiral stage. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 9 pages, 7 figures and 1 table

The Santa Barbara Binary-Disk Code Comparison

Authors: Paul C. Duffell, Alexander J. Dittmann, Daniel J. D'Orazio, Alessia Franchini, Kaitlin M. Kratter, Anna B. T. Penzlin, Enrico Ragusa, Magdalena Siwek, Christopher Tiede, Haiyang Wang, Jonathan Zrake, Adam M. Dempsey, Zoltan Haiman, Alessandro Lupi, Michal Pirog, Geoffrey Ryan

Abstract: We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed… ▽ More We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, $0.2 (a^2 Ω_B /ν)$ binary orbits, where $ν$ is the kinematic viscosity. We did not investigate dependence on binary mass ratio, eccentricity, disk temperature, or disk viscosity; therefore, these benchmarks may act as guides towards expanding converged solutions to the wider parameter space but might need to be updated in a future study that investigates dependence on system parameters. We find the most major discrepancies between codes resulted from the dimensionality of the setup (3D vs 2D disks). Beyond this, we find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D. △ Less

Submitted 20 June, 2024; v1 submitted 20 February, 2024; originally announced February 2024.

Comments: ApJ Accepted

LISA Definition Study Report

Authors: Monica Colpi, Karsten Danzmann, Martin Hewitson, Kelly Holley-Bockelmann, Philippe Jetzer, Gijs Nelemans, Antoine Petiteau, David Shoemaker, Carlos Sopuerta, Robin Stebbins, Nial Tanvir, Henry Ward, William Joseph Weber, Ira Thorpe, Anna Daurskikh, Atul Deep, Ignacio Fernández Núñez, César García Marirrodriga, Martin Gehler, Jean-Philippe Halain, Oliver Jennrich, Uwe Lammers, Jonan Larrañaga, Maike Lieser, Nora Lützgendorf , et al. (86 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the e… ▽ More The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 155 pages, with executive summary and table of contents

Self-lensing flares from black hole binaries III: general-relativistic ray tracing of circumbinary accretion simulations

Authors: Luke Major Krauth, Jordy Davelaar, Zoltán Haiman, John Ryan Westernacher-Schneider, Jonathan Zrake, Andrew MacFadyen

Abstract: Self-lensing flares (SLFs) are expected to be produced once or twice per orbit by an accreting massive black hole binary (MBHB), if the eclipsing MBHBs are observed close to edge-on. SLFs can provide valuable electromagnetic (EM) signatures to accompany the gravitational waves (GWs) detectable by the upcoming Laser Interferometer Space Antenna (LISA). EM follow-ups are crucial for, e.g., sky-local… ▽ More Self-lensing flares (SLFs) are expected to be produced once or twice per orbit by an accreting massive black hole binary (MBHB), if the eclipsing MBHBs are observed close to edge-on. SLFs can provide valuable electromagnetic (EM) signatures to accompany the gravitational waves (GWs) detectable by the upcoming Laser Interferometer Space Antenna (LISA). EM follow-ups are crucial for, e.g., sky-localization, and constraining the Hubble constant and the graviton mass. We use high-resolution two-dimensional viscous hydrodynamical simulations of a circumbinary disk (CBD) embedding a MBHB. We then use very high-cadence output of these hydrodynamical simulation inputs for a general-relativistic ray-tracing code to produce synthetic spectra and phase-folded light curves. Our main results show a significant periodic amplification of the flux with the characteristic shape of a sharp flare with a central dip, as the foreground black hole (BH) transits across the minidisk and shadow of the background BH, respectively. These corroborate previous conclusions based on the microlensing approximation and analytical toy models of the emission geometry. We also find that at lower inclinations, without some occlusion of the minidisk emission by the CBD, shocks from quasi-periodic mass-trading between the minidisks can produce bright flares which can mimic SLFs and could hinder their identification. △ Less

Submitted 13 November, 2023; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: 14 pages, 11 figures, submitted to journal, split Fig. 1 by frequency, fixed some typos

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.

Comments: 15 pages, 3 figures

Diagnosing the massive-seed pathway to high-redshift black holes: statistics of the evolving black hole to host galaxy mass ratio

Authors: Matthew T. Scoggins, Zoltán Haiman

Abstract: Supermassive black holes (SMBHs) with masses of $\sim 10^9 {\rm M_\odot}$ within the first billion year of the universe challenge our conventional understanding of black hole formation and growth. One pathway to these SMBHs proposes that supermassive stars (SMSs) born in pristine atomic cooling haloes (ACHs) yield massive seed BHs evolving to these early SMBHs. This scenario leads to an overly mas… ▽ More Supermassive black holes (SMBHs) with masses of $\sim 10^9 {\rm M_\odot}$ within the first billion year of the universe challenge our conventional understanding of black hole formation and growth. One pathway to these SMBHs proposes that supermassive stars (SMSs) born in pristine atomic cooling haloes (ACHs) yield massive seed BHs evolving to these early SMBHs. This scenario leads to an overly massive BH galaxy (OMBG), in which the BH to stellar mass ratio is initially $M_{\rm bh}/M_* \geq 1$, well in excess of the typical values of $\sim 10^{-3}$ at low redshifts. Previously, we have investigated two massive seed BH candidates from the \texttt{Renaissance} simulation and found that they remain outliers on the $M_{\rm bh}-M_{*}$ relation until the OMBG merges with a much more massive halo at $z{=}8$. In this work, we use Monte-Carlo merger trees to investigate the evolution of the $M_{\rm bh}-M_{*}$ relation for $50,000$ protogalaxies hosting massive BH seeds, across $10,000$ trees that merge into a $10^{12} {\rm M_\odot}$ halo at $z{=}6$. We find that up to $60\%$ (depending on growth parameters) of these OMBGs remain strong outliers for several 100 Myr, down to redshifts detectable with {\it JWST} and with sensitive X-ray telescopes. This represents a way to diagnose the massive-seed formation pathway for early SMBHs. We expect to find ${\sim} 0.1{-}1$ of these objects per {\it JWST} NIRCam field per unit redshift at $z\gtrsim 6$. Recently detected SMBHs with masses of $\sim 10^7~{\rm M_\odot}$ and low inferred stellar-mass hosts may be examples of this population. △ Less

Submitted 7 June, 2024; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: 14 pages, 9 figures, accepted for publication in MNRAS

Black Hole Binaries in AGN Accretion Discs II: Gas Effects on Black Hole Satellite Scatterings

Authors: Connar Rowan, Henry Whitehead, Tjarda Boekholt, Bence Kocsis, Zoltán Haiman

Abstract: The black hole (BH) binaries in active galactic nuclei (AGN) are expected to form mainly through scattering encounters in the ambient gaseous medium. Recent simulations, including our own, have confirmed this formation pathway is highly efficient. We perform 3D smoothed particle hydrodynamics (SPH) simulations of BH scattering encounters in AGN disks. Using a range of impact parameters, we probe t… ▽ More The black hole (BH) binaries in active galactic nuclei (AGN) are expected to form mainly through scattering encounters in the ambient gaseous medium. Recent simulations, including our own, have confirmed this formation pathway is highly efficient. We perform 3D smoothed particle hydrodynamics (SPH) simulations of BH scattering encounters in AGN disks. Using a range of impact parameters, we probe the necessary conditions for binary capture and how different orbital trajectories affect the dissipative effects from the gas. We identify a single range of impact parameters, typically of width $\sim0.86-1.59$ binary Hill radii depending on AGN disk density, that reliably leads to binary formation. The periapsis of the first encounter is the primary variable that determines the outcome of the initial scattering. We find an associated power-law between the energy dissipated and the periapsis depth to be $ΔE\propto r^{-b}$ with $b=0.42\pm0.16$, where deeper encounters dissipate more energy. Excluding accretion physics does not significantly alter these results. We identify the region of parameter space in initial energy vs impact parameter where a scattering leads to binary formation. Based on our findings, we provide a ready-to-use analytic criterion that utilises these two pre-encounter parameters to determine the outcome of an encounter, with a reliability rate of >90\%. As the criterion is based directly on our simulations, it provides a reliable and highly physically motivated criterion for predicting binary scattering outcomes which can be used in population studies of BH binaries and mergers around AGN. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: 20 pages, 19 figures

The Enhanced Population of Extreme Mass-Ratio Inspirals in the LISA Band from Supermassive Black Hole Binaries

Authors: Smadar Naoz, Zoltan Haiman

Abstract: Extreme mass ratio inspirals (EMRIs) take place when a stellar-mass black hole (BH) merges with a supermassive black hole (SMBH). The gravitational wave emission from such an event is expected to be detectable by the future Laser Interferometer Space Antenna (LISA) and other mHz detectors. It was recently suggested that the EMRI rate in SMBH binary systems is orders of magnitude higher than the EM… ▽ More Extreme mass ratio inspirals (EMRIs) take place when a stellar-mass black hole (BH) merges with a supermassive black hole (SMBH). The gravitational wave emission from such an event is expected to be detectable by the future Laser Interferometer Space Antenna (LISA) and other mHz detectors. It was recently suggested that the EMRI rate in SMBH binary systems is orders of magnitude higher than the EMRI rate around a single SMBH with the same total mass. Here we show that this high rate can produce thousands of SMBH-BH sources at a redshift of unity. We predict that LISA may detect a few hundred of these EMRIs with signal-to-noise ratio above SNR>=8 within a four-year mission lifetime. The remaining sub-threshold sources will contribute to a large confusion noise, which is approximately an order of magnitude above LISA's sensitivity level. Finally, we suggest that the individually detectable systems, as well as the background noise from the sub-threshold EMRIs, can be used to constrain the SMBH binary fraction in the low-redshift Universe. △ Less

Submitted 18 September, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Comments: 9 pages, 5 figures, accepted to ApJ letter

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.

Comments: 14 pages, 8 figures, accepted in ApJ

Eccentric Minidisks in Accreting Binaries

Authors: John Ryan Westernacher-Schneider, Jonathan Zrake, Andrew MacFadyen, Zoltán Haiman

Abstract: We show that gas disks around the components of an orbiting binary system (so-called minidisks) may be susceptible to a resonant instability which causes the minidisks to become significantly eccentric. Eccentricity is injected by, and also induces, regular impacts between the minidisks at roughly the orbital period of the binary. Such eccentric minidisks are seen in vertically integrated, two-dim… ▽ More We show that gas disks around the components of an orbiting binary system (so-called minidisks) may be susceptible to a resonant instability which causes the minidisks to become significantly eccentric. Eccentricity is injected by, and also induces, regular impacts between the minidisks at roughly the orbital period of the binary. Such eccentric minidisks are seen in vertically integrated, two-dimensional simulations of a circular, equal-mass binary accreting from a circumbinary gas disk with a $Γ$-law equation of state. Minidisk eccentricity is suppressed by the use of an isothermal equation of state. However, the instability still operates, and can be revealed in a minimal disk-binary simulation by removing the circumbinary disk, and feeding the minidisks from the component positions. Minidisk eccentricity is also suppressed when the gravitational softening length is large ($\gtrsim 4\%$ of the binary semi-major axis), suggesting that its absence could be an artifact of widely adopted numerical approximations; a follow-up study in three dimensions with well-resolved, geometrically thin minidisks (aspect ratios $\lesssim 0.02$) may be needed to assess whether eccentric minidisks can occur in real astrophysical environments. If they can, the electromagnetic signature may be important for discriminating between binary and single black hole scenarios for quasi-periodic oscillations in active galactic nuclei; in turn, this might aid in targeted searches with pulsar timing arrays for individual supermassive black hole binary sources of low-frequency gravitational waves. △ Less

Submitted 13 February, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 5 Mb, 22 pages, video content URLs provided in the captions to Figures 5, 6, 9, 11, and 15. Version 2 is close to the published, open-access version (see the DOI to the journal version below)

Journal ref: The Astrophysical Journal, 962:76 (17pp), 2024 February 10

Massive Black Hole Binaries as LISA Precursors in the Roman High Latitude Time Domain Survey

Authors: Zoltán Haiman, Chengcheng Xin, Tamara Bogdanović, Pau Amaro Seoane, Matteo Bonetti, J. Andrew Casey-Clyde, Maria Charisi, Monica Colpi, Jordy Davelaar, Alessandra De Rosa, Daniel J. D'Orazio, Kate Futrowsky, Poshak Gandhi, Alister W. Graham, Jenny E. Greene, Melanie Habouzit, Daryl Haggard, Kelly Holley-Bockelmann, Xin Liu, Alberto Mangiagli, Alessandra Mastrobuono-Battisti, Sean McGee, Chiara M. F. Mingarelli, Rodrigo Nemmen, Antonella Palmese , et al. (5 additional authors not shown)

Abstract: With its capacity to observe $\sim 10^{5-6}$ faint active galactic nuclei (AGN) out to redshift $z\approx 6$, Roman is poised to reveal a population of $10^{4-6}\, {\rm M_\odot}$ black holes during an epoch of vigorous galaxy assembly. By measuring the light curves of a subset of these AGN and looking for periodicity, Roman can identify several hundred massive black hole binaries (MBHBs) with 5-12… ▽ More With its capacity to observe $\sim 10^{5-6}$ faint active galactic nuclei (AGN) out to redshift $z\approx 6$, Roman is poised to reveal a population of $10^{4-6}\, {\rm M_\odot}$ black holes during an epoch of vigorous galaxy assembly. By measuring the light curves of a subset of these AGN and looking for periodicity, Roman can identify several hundred massive black hole binaries (MBHBs) with 5-12 day orbital periods, which emit copious gravitational radiation and will inevitably merge on timescales of $10^{3-5}$ years. During the last few months of their merger, such binaries are observable with the Laser Interferometer Space Antenna (LISA), a joint ESA/NASA gravitational wave mission set to launch in the mid-2030s. Roman can thus find LISA precursors, provide uniquely robust constraints on the LISA source population, help identify the host galaxies of LISA mergers, and unlock the potential of multi-messenger astrophysics with massive black hole binaries. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: White Paper for the Nancy Grace Roman Space Telescope's Core Community Surveys (https://roman.gsfc.nasa.gov/science/ccs_white_papers.html)

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): JWST Reveals a Filamentary Structure around a z=6.61 Quasar

Authors: Feige Wang, Jinyi Yang, Joseph F. Hennawi, Xiaohui Fan, Fengwu Sun, Jaclyn B. Champagne, Tiago Costa, Melanie Habouzit, Ryan Endsley, Zihao Li, Xiaojing Lin, Romain A. Meyer, Jan-Torge Schindler, Yunjing Wu, Eduardo Bañados, Aaron J. Barth, Aklant K. Bhowmick, Rebekka Bieri, Laura Blecha, Sarah Bosman, Zheng Cai, Luis Colina, Thomas Connor, Frederick B. Davies, Roberto Decarli , et al. (34 additional authors not shown)

Abstract: We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) gal… ▽ More We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) galaxy redshift survey at 3-4 $μ$m among JWST Cycle-1 programs and provide extensive legacy values for studying the formation of the earliest supermassive black holes (SMBHs), the assembly of galaxies, early metal enrichment, and cosmic reionization. In this first ASPIRE paper, we report the discovery of a filamentary structure traced by the luminous quasar J0305-3150 and ten [OIII] emitters at $z=6.6$. This structure has a 3D galaxy overdensity of $δ_{\rm gal}=12.6$ over 637 cMpc$^3$, one of the most overdense structures known in the early universe, and could eventually evolve into a massive galaxy cluster. Together with existing VLT/MUSE and ALMA observations of this field, our JWST observations reveal that J0305-3150 traces a complex environment where both UV-bright and dusty galaxies are present, and indicate that the early evolution of galaxies around the quasar is not simultaneous. In addition, we discovered 31 [OIII] emitters in this field at other redshifts, $5.3<z<6.7$, with half of them situated at $z\sim5.4$ and $z\sim6.2$. This indicates that star-forming galaxies, such as [OIII] emitters, are generally clustered at high redshifts. These discoveries demonstrate the unparalleled redshift survey capabilities of NIRCam WFSS and the potential of the full ASPIRE survey dataset. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: accepted for publication in ApJL

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): A First Look at the Rest-frame Optical Spectra of $z > 6.5$ Quasars Using JWST

Authors: Jinyi Yang, Feige Wang, Xiaohui Fan, Joseph F. Hennawi, Aaron J. Barth, Eduardo Bañados, Fengwu Sun, Weizhe Liu, Zheng Cai, Linhua Jiang, Zihao Li, Masafusa Onoue, Jan-Torge Schindler, Yue Shen, Yunjing Wu, Aklant K. Bhowmick, Rebekka Bieri, Laura Blecha, Sarah Bosman, Jaclyn B. Champagne, Luis Colina, Thomas Connor, Tiago Costa, Frederick B. Davies, Roberto Decarli , et al. (31 additional authors not shown)

Abstract: Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''… ▽ More Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE)" program. Our JWST spectra cover the quasars' emission between rest frame $\sim$ 4100 and 5100 Å. The profiles of these quasars' broad H$β$ emission lines span a FWHM from 3000 to 6000 $\rm{km~s^{-1}}$. The H$β$-based virial black hole (BH) masses, ranging from 0.6 to 2.1 billion solar masses, are generally consistent with their MgII-based BH masses. The new measurements based on the more reliable H$β$ tracer thus confirm the existence of billion solar-mass BHs in the reionization epoch. In the observed [OIII] $λλ$4960,5008 doublets of these luminous quasars, broad components are more common than narrow core components ($\le~1200~\rm{km~s^{-1}}$), and only one quasar shows stronger narrow components than broad. Two quasars exhibit significantly broad and blueshifted [OIII] emission, thought to trace galactic-scale outflows, with median velocities of $-610~\rm{km~s^{-1}}$ and $-1430~\rm{km~s^{-1}}$ relative to the [CII] $158\,μ$m line. All eight quasars show strong optical FeII emission, and follow the Eigenvector 1 relations defined by low-redshift quasars. The entire ASPIRE program will eventually cover 25 quasars and provide a statistical sample for the studies of the BHs and quasar spectral properties. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 13 pages, 4 figures, accepted for publication in ApJL

NuSTAR Observations of Candidate Subparsec Binary Supermassive Black Holes

Authors: M. Lynne Saade, Murray Brightman, Daniel Stern, Thomas Connor, S. G. Djorgovski, Daniel J. D'Orazio, K. E. S. Ford, Matthew J. Graham, Zoltan Haiman, Hyunsung D. Jun, Elias Kammoun, Ralph P. Kraft, Barry McKernan, Alexei Vikhlinin, Dominic J. Walton

Abstract: We present analysis of NuSTAR X-ray observations of three AGN that were identified as candidate subparsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey based on apparent periodicity in their optical light curves. Simulations predict that close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs. We previously o… ▽ More We present analysis of NuSTAR X-ray observations of three AGN that were identified as candidate subparsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey based on apparent periodicity in their optical light curves. Simulations predict that close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs. We previously observed these AGN with Chandra and found no differences between their low energy X-ray properties and the larger AGN population. However some models predict differences to be more prominent at energies higher than probed by Chandra. We find that even at the higher energies probed by NuSTAR, the spectra of these AGN are indistinguishable from the larger AGN population. This could rule out models predicting large differences in the X-ray spectra in the NuSTAR bands. Alternatively, it might mean that these three AGN are not binary SMBHs. △ Less

Submitted 27 March, 2024; v1 submitted 12 April, 2023; originally announced April 2023.

The importance of live binary evolution in numerical simulations of binaries embedded in circumbinary discs

Authors: Alessia Franchini, Alessandro Lupi, Alberto Sesana, Zoltan Haiman

Abstract: The shrinking of a binary orbit driven by the interaction with a gaseous circumbinary disc, initially advocated as a potential way to catalyze the binary merger, has been recently debated in the case of geometrically thick (i.e. with $H/R\gtrsim 0.1$) discs. However, a clear consensus is still missing mainly owing to numerical limitations, such as fixed orbit binaries or lack of resolution inside… ▽ More The shrinking of a binary orbit driven by the interaction with a gaseous circumbinary disc, initially advocated as a potential way to catalyze the binary merger, has been recently debated in the case of geometrically thick (i.e. with $H/R\gtrsim 0.1$) discs. However, a clear consensus is still missing mainly owing to numerical limitations, such as fixed orbit binaries or lack of resolution inside the cavity carved by the binary in its circumbinary disc. In this work, we asses the importance of evolving the binary orbit by means of hydrodynamic simulations performed with the code {\sc gizmo} in meshless-finite-mass mode. In order to model the interaction between equal mass circular binaries and their locally isothermal circumbinary discs, we enforce hyper-Lagrangian resolution inside the cavity. We find that fixing the binary orbit ultimately leads to an overestimate of the gravitational torque that the gas exerts on the binary, and in an underestimate of the torque due to the accretion of material onto the binary components. Furthermore, we find that the modulation of the accretion rate on the binary orbital period is strongly suppressed in the fixed orbit simulation, while it is clearly present in the live binary simulations. This has potential implications for the prediction of the observable periodicities in massive black hole binary candidates. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: 7 pages, 8 figures, accepted for publication in MNRAS

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.

Comments: 12 pages, 8 figures, accepted in MNRAS

Disappearing thermal X-ray emission as a tell-tale signature of merging massive black hole binaries

Authors: Luke Major Krauth, Jordy Davelaar, Zoltán Haiman, John Ryan Westernacher-Schneider, Jonathan Zrake, Andrew MacFadyen

Abstract: The upcoming Laser Interferometer Space Antenna (LISA) is expected to detect gravitational waves (GWs) from massive black hole binaries (MBHB). Finding the electromagnetic (EM) counterparts for these GW events will be crucial for understanding how and where MBHBs merge, measuring their redshifts, constraining the Hubble constant and the graviton mass, and for other novel science applications. Howe… ▽ More The upcoming Laser Interferometer Space Antenna (LISA) is expected to detect gravitational waves (GWs) from massive black hole binaries (MBHB). Finding the electromagnetic (EM) counterparts for these GW events will be crucial for understanding how and where MBHBs merge, measuring their redshifts, constraining the Hubble constant and the graviton mass, and for other novel science applications. However, due to poor GW sky localisation, multi-wavelength, time-dependent electromagnetic (EM) models are needed to identify the right host galaxy among many candidates. We studied merging MBHBs embedded in a circumbinary disc using high-resolution two-dimensional simulations, with a $Γ$-law equation of state, incorporating viscous heating, shock heating, and radiative cooling. We simulate the binary from large separation until after merger, allowing us to model the decoupling of the binary from the circumbinary disc (CBD). We compute the EM signatures and identify distinct features before, during, and after the merger. Our main result is a multi-band EM signature: we find that the MBHB produces strong thermal X-ray emission until 1-2 days prior to the merger. However, as the binary decouples from the CBD, the X-ray-bright minidiscs rapidly shrink in size, become disrupted, and the accretion rate drops precipitously. As a result, the thermal X-ray luminosity drops by orders of magnitude, and the source remains X-ray dark for several days after the merger, regardless of any post-merger effects such as GW recoil or mass loss. Looking for the abrupt spectral change where the thermal X-ray disappears is a tell-tale EM signature of LISA mergers that does not require extensive pre-merger monitoring. △ Less

Submitted 6 October, 2023; v1 submitted 5 April, 2023; originally announced April 2023.

Comments: 14 pages, 16 figures, 1 table, accepted for publication in MNRAS

Tidal Peeling Events: low-eccentricity tidal disruption of a star by a stellar-mass black hole

Authors: Chengcheng Xin, Zoltan Haiman, Rosalba Perna, Yihan Wang, Taeho Ryu

Abstract: Close encounters between stellar-mass black holes (BHs) and stars occur frequently in dense star clusters and in the disks of active galactic nuclei (AGNs). Recent studies have shown that in highly eccentric close encounters, the star can be tidally disrupted by the BH (micro-tidal disruption event, or micro-TDE), resulting in rapid mass accretion and possibly bright electromagnetic signatures. He… ▽ More Close encounters between stellar-mass black holes (BHs) and stars occur frequently in dense star clusters and in the disks of active galactic nuclei (AGNs). Recent studies have shown that in highly eccentric close encounters, the star can be tidally disrupted by the BH (micro-tidal disruption event, or micro-TDE), resulting in rapid mass accretion and possibly bright electromagnetic signatures. Here we consider a scenario in which the star might approach the stellar-mass BH in a gradual, nearly circular inspiral, under the influence of dynamical friction on a circum-binary gas disk or three-body interactions in a star cluster. We perform hydro-dynamical simulations of this scenario using the smoothed particle hydrodynamics code PHANTOM. We find that the mass of the star is slowly stripped away by the BH. We call this gradual tidal disruption a "tidal-peeling event", or a TPE. Depending on the initial distance and eccentricity of the encounter, TPEs might exhibit significant accretion rates and orbital evolution distinct from those of a typical (eccentric) micro-TDE. △ Less

Submitted 22 March, 2023; originally announced March 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.

Comments: 10 pages, 4 figures, accepted in ApJL

Absence of the predicted 2022 October outburst of OJ 287 and implications for binary SMBH scenarios

Authors: S. Komossa, D. Grupe, A. Kraus, M. A. Gurwell, Z. Haiman, F. K. Liu, A. Tchekhovskoy, L. C. Gallo, M. Berton, R. Blandford, J. L. Gomez, A. G. Gonzalez

Abstract: The project MOMO (Multiwavelength Observations and Modelling of OJ 287) was set up to test predictions of binary supermassive black hole (SMBH) scenarios and to understand disk-jet physics of the blazar OJ 287. After a correction, the precessing binary (PB) SMBH model predicted the next main outburst of OJ 287 in 2022 October, no longer in July, making the outburst well observable and the model te… ▽ More The project MOMO (Multiwavelength Observations and Modelling of OJ 287) was set up to test predictions of binary supermassive black hole (SMBH) scenarios and to understand disk-jet physics of the blazar OJ 287. After a correction, the precessing binary (PB) SMBH model predicted the next main outburst of OJ 287 in 2022 October, no longer in July, making the outburst well observable and the model testable. We have densely covered this period in our ongoing multi-frequency radio, optical, UV, and X-ray monitoring. The predicted outburst was not detected. Instead, OJ 287 was at low optical-UV emission levels, declining further into November. The predicted thermal bremsstrahlung spectrum was not observed either, at any epoch. Further, applying scaling relations, we estimate a SMBH mass of OJ 287 of 10^8 M_sun. The latest in a sequence of deep low-states that recur every 1-2 yrs is used to determine an upper limit on the Eddington ratio and on the accretion-disk luminosity. This limit is at least a factor of 10 lower than required by the PB model with its massive primary SMBH of >10^{10} M_sun. All these results favor alternative binary SMBH models of OJ 287 that neither require strong orbital precession nor a very large mass of the primary SMBH. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: MNRAS Letter, in press

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.

Comments: 21 pages, 6 figures, accepted in ApJ

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.

Comments: 7 pages, 4 figures, and 1 table

Cosmological constraints from HSC survey first-year data using deep learning

Authors: Tianhuan Lu, Zoltán Haiman, Xiangchong Li

Abstract: We present cosmological constraints from the Subaru Hyper Suprime-Cam (HSC) first-year weak lensing shear catalogue using convolutional neural networks (CNNs) and conventional summary statistics. We crop 19 $3\times3\,\mathrm{{deg}^2}$ sub-fields from the first-year area, divide the galaxies with redshift $0.3\le z\le1.5$ into four equally-spaced redshift bins, and perform tomographic analyses. We… ▽ More We present cosmological constraints from the Subaru Hyper Suprime-Cam (HSC) first-year weak lensing shear catalogue using convolutional neural networks (CNNs) and conventional summary statistics. We crop 19 $3\times3\,\mathrm{{deg}^2}$ sub-fields from the first-year area, divide the galaxies with redshift $0.3\le z\le1.5$ into four equally-spaced redshift bins, and perform tomographic analyses. We develop a pipeline to generate simulated convergence maps from cosmological $N$-body simulations, where we account for effects such as intrinsic alignments (IAs), baryons, photometric redshift errors, and point spread function errors, to match characteristics of the real catalogue. We train CNNs that can predict the underlying parameters from the simulated maps, and we use them to construct likelihood functions for Bayesian analyses. In the $Λ$ cold dark matter model with two free cosmological parameters $Ω_\mathrm{m}$ and $σ_8$, we find $Ω_\mathrm{m}=0.278_{-0.035}^{+0.037}$, $S_8\equiv(Ω_\mathrm{m}/0.3)^{0.5}σ_8=0.793_{-0.018}^{+0.017}$, and the IA amplitude $A_\mathrm{IA}=0.20_{-0.58}^{+0.55}$. In a model with four additional free baryonic parameters, we find $Ω_\mathrm{m}=0.268_{-0.036}^{+0.040}$, $S_8=0.819_{-0.024}^{+0.034}$, and $A_\mathrm{IA}=-0.16_{-0.58}^{+0.59}$, with the baryonic parameters not being well-constrained. We also find that statistical uncertainties of the parameters by the CNNs are smaller than those from the power spectrum (5--24 percent smaller for $S_8$ and a factor of 2.5--3.0 smaller for $Ω_\mathrm{m}$), showing the effectiveness of CNNs for uncovering additional cosmological information from the HSC data. With baryons, the $S_8$ discrepancy between HSC first-year data and Planck 2018 is reduced from $\sim2.2\,σ$ to $0.3\text{--}0.5\,σ$. △ Less

Submitted 1 March, 2023; v1 submitted 3 January, 2023; originally announced January 2023.

Comments: 22 pages, 14 figures

Black Hole Binary Formation in AGN Discs: From Isolation to Merger

Authors: Connar Rowan, Tjarda Boekholt, Bence Kocsis, Zoltán Haiman

Abstract: Motivated by the increasing number of detections of merging black holes by LIGO-VIRGO-KAGRA, black hole (BH) binary mergers in the discs of active galactic nuclei (AGN) is investigated as a possible merger channel. In this pathway, BH encounters in the gas disc form mutually bound black hole binary systems through interaction with the gas in the disc and subsequently inspiral through gravitational… ▽ More Motivated by the increasing number of detections of merging black holes by LIGO-VIRGO-KAGRA, black hole (BH) binary mergers in the discs of active galactic nuclei (AGN) is investigated as a possible merger channel. In this pathway, BH encounters in the gas disc form mutually bound black hole binary systems through interaction with the gas in the disc and subsequently inspiral through gravitational torques induced by the local gas. To determine the feasibility of this merger pathway, we present the first 3D global hydrodynamic simulations of the formation and evolution of a stellar-mass BH binaries AGN discs with three different AGN disc masses and five different initial radial separations. These 15 simulations show binary capture of prograde and retrograde binaries can be successful in a range of disc densities including cases well below that of a standard radiatively efficient alpha disc, identifying that the majority of these captured binaries are then subsequently hardened by the surrounding gas. The eccentricity evolution depends strongly on the orbital rotation where prograde binaries are governed by gravitational torques form their circumbinary mini-disc, with eccentricities being damped, while for retrograde binaries the eccentricities are excited to > 0.9 by accretion torques. In two cases, retrograde binaries ultimately undergo a close periapsis passage which results in a merger via gravitational waves after only a few thousand binary orbits. Thus, the merger timescale can be far shorter than the AGN disc lifetime. These simulations support an efficient AGN disc merger pathway for BHs. △ Less

Submitted 23 June, 2023; v1 submitted 12 December, 2022; originally announced December 2022.

Comments: 27 pages, 28 Figures, accepted for publication in MNRAS

The X-ray view of optically selected dual AGN

Authors: Alessandra De Rosa, Cristian Vignali, Paola Severgnini, Stefano Bianchi, Tamara Bogdanović, Maria Charisi, Matteo Guainazzi, Zoltan Haiman, S. Komossa, Zsolt Paragi, Miguel Perez-Torres, Enrico Piconcelli, Lorenzo Ducci, Manali Parvatikar, Roberto Serafinelli

Abstract: We present a study of optically selected dual AGN with projected separations of 3--97~kpc. Using multi-wavelength (MWL) information (optical, X-rays, mid-IR), we characterized the intrinsic nuclear properties of this sample and compared them with those of isolated systems. Among the 124 X-ray detected AGN candidates, 52 appear in pairs and 72 as single X-ray sources. Through MWL analysis, we confi… ▽ More We present a study of optically selected dual AGN with projected separations of 3--97~kpc. Using multi-wavelength (MWL) information (optical, X-rays, mid-IR), we characterized the intrinsic nuclear properties of this sample and compared them with those of isolated systems. Among the 124 X-ray detected AGN candidates, 52 appear in pairs and 72 as single X-ray sources. Through MWL analysis, we confirmed the presence of the AGN in a fraction >80\% of the detected targets in pairs (42 over 52). X-ray spectral analysis confirms the trend of increasing AGN luminosity with decreasing separation, suggesting that mergers may have contributed in triggering more luminous AGN. Through X/mid-IR ratio $vs$ X-ray colors, we estimated a fraction of Compton-thin AGN (with 10$^{22}$ cm$^{-2}$ $<$ N$_{\rm H} <$10$^{24}$ cm$^{-2}$) of about 80\%, while about 16\% are Compton thick (CT, with N$_{\rm H}>$10$^{24}$ cm$^{-2}$) sources. These fractions of obscured sources are larger than those found in samples of isolated AGN, confirming that pairs of AGN show higher obscuration. This trend is further confirmed by comparing the de-reddened [O\ III] emission with the observed X-ray luminosity. However, the derived fraction of Compton-thick sources in this sample at early stage of merging is lower than reported for late-merging dual-AGN samples. Comparing N$_{\rm H}$ from X-rays with that derived from E(B-V) from Narrow Line Regions, we find that the absorbing material is likely associated with the torus or the Broad Line Regions. We also explored the X-ray detection efficiency of dual-AGN candidates, finding that, when observed properly (at on-axis positions and with long exposures), X-ray data represent a powerful way to confirm and investigate dual-AGN systems. △ Less

Submitted 12 December, 2022; originally announced December 2022.

Comments: 12 pages, 4 figures. Accepted for publication in MNRAS

Gas dynamical friction as a binary formation mechanism in AGN discs

Authors: Stanislav DeLaurentiis, Marguerite Epstein-Martin, Zoltan Haiman

Abstract: In this paper, we study how gaseous dynamical friction (DF) affects the motion of fly-by stellar-mass black holes (sBHs) embedded in active galactic nucleus (AGN) discs. We perform 3-body integrations of the interaction of two co-planar sBHs in nearby, initially circular orbits around the supermassive black hole (SMBH). We find that DF can facilitate the formation of gravitationally bound near-Kep… ▽ More In this paper, we study how gaseous dynamical friction (DF) affects the motion of fly-by stellar-mass black holes (sBHs) embedded in active galactic nucleus (AGN) discs. We perform 3-body integrations of the interaction of two co-planar sBHs in nearby, initially circular orbits around the supermassive black hole (SMBH). We find that DF can facilitate the formation of gravitationally bound near-Keplerian binaries in AGN discs, and we delineate the discrete ranges of impact parameters and AGN disc parameters for which such captures occur. We also report trends in the bound binaries' eccentricity and sense of rotation (prograde or retrograde with respect to the background AGN disc) as a function of the impact parameter of the initial encounter. While based on an approximate description of gaseous friction, our results suggest that binary formation in AGN discs should be common and may produce both prograde and retrograde, as well as both circular and eccentric binaries. △ Less

Submitted 14 July, 2023; v1 submitted 5 December, 2022; originally announced December 2022.

Comments: 16 pages, 13 figures, Published in MNRAS

Journal ref: Volume 523, Issue 1, July 2023, Pages 1126:1139

Detection of stellar light from quasar host galaxies at redshifts above 6

Authors: Xuheng Ding, Masafusa Onoue, John D. Silverman, Yoshiki Matsuoka, Takuma Izumi, Michael A. Strauss, Knud Jahnke, Camryn L. Phillips, Junyao Li, Marta Volonteri, Zoltan Haiman, Irham Taufik Andika, Kentaro Aoki, Shunsuke Baba, Rebekka Bieri, Sarah E. I. Bosman, Connor Bottrell, Anna-Christina Eilers, Seiji Fujimoto, Melanie Habouzit, Masatoshi Imanishi, Kohei Inayoshi, Kazushi Iwasawa, Nobunari Kashikawa, Toshihiro Kawaguchi , et al. (19 additional authors not shown)

Abstract: The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting the… ▽ More The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting their underlying, previously-undetected host galaxies. Here we report rest-frame optical images and spectroscopy of two HSC-SSP quasars at $z>6$ with JWST. Using NIRCam imaging at 3.6$μ$m and 1.5$μ$m and subtracting the light from the unresolved quasars, we find that the host galaxies are massive (stellar masses of $13\times$ and $3.4\times$ $10^{10}$ M$_{\odot}$, respectively), compact, and disk-like. NIRSpec medium-resolution spectroscopy shows stellar absorption lines in the more massive quasar, confirming the detection of the host. Velocity-broadened gas in the vicinity of these quasars enables measurements of their black hole masses ($1.4\times 10^9$ and $2.0\times$ $10^{8}$ M$_{\odot}$, respectively). Their location in the black hole mass - stellar mass plane is consistent with the distribution at low redshift, suggesting that the relation between black holes and their host galaxies was already in place less than a billion years after the Big Bang. △ Less

Submitted 23 June, 2023; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: Matched to the published Nature version of the article. 27 pages, 4 main figures, 1 table, 6 supplementary figures, 2 supplementary table

Report of the Topical Group on Cosmic Probes of Fundamental Physics for for Snowmass 2021

Authors: Rana X. Adhikari, Luis A. Anchordoqui, Ke Fang, B. S. Sathyaprakash, Kirsten Tollefson, Tiffany R. Lewis, Kristi Engel, Amin Aboubrahim, Ozgur Akarsu, Yashar Akrami, Roberto Aloisio, Rafael Alves Batista, Mario Ballardini, Stefan W. Ballmer, Ellen Bechtol, David Benisty, Emanuele Berti, Simon Birrer, Alexander Bonilla, Richard Brito, Mauricio Bustamante, Robert Caldwell, Vitor Cardoso, Sukanya Chakrabarti, Thomas Y. Chen , et al. (96 additional authors not shown)

Abstract: Cosmic Probes of Fundamental Physics take two primary forms: Very high energy particles (cosmic rays, neutrinos, and gamma rays) and gravitational waves. Already today, these probes give access to fundamental physics not available by any other means, helping elucidate the underlying theory that completes the Standard Model. The last decade has witnessed a revolution of exciting discoveries such as… ▽ More Cosmic Probes of Fundamental Physics take two primary forms: Very high energy particles (cosmic rays, neutrinos, and gamma rays) and gravitational waves. Already today, these probes give access to fundamental physics not available by any other means, helping elucidate the underlying theory that completes the Standard Model. The last decade has witnessed a revolution of exciting discoveries such as the detection of high-energy neutrinos and gravitational waves. The scope for major developments in the next decades is dramatic, as we detail in this report. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: Report of theTopical Group on Cosmic Probes of Fundamental Physics, for the U.S. decadal Particle Physics Planning Exercise (Snowmass 2021)

Self-regulation of black hole accretion via jets in early protogalaxies

Authors: Kung-Yi Su, Greg L. Bryan, Zoltán Haiman, Rachel S. Somerville, Christopher C. Hayward, Claude-André Faucher-Giguère

Abstract: The early growth of black holes (BHs) in high-redshift galaxies is likely regulated by their feedback on the surrounding gas. While radiative feedback has been extensively studied, the role of mechanical feedback has received comparatively less scrutiny to date. Here we use high-resolution parsec-scale hydrodynamical simulations to study jet propagation and its effect on BH accretion onto 100… ▽ More The early growth of black holes (BHs) in high-redshift galaxies is likely regulated by their feedback on the surrounding gas. While radiative feedback has been extensively studied, the role of mechanical feedback has received comparatively less scrutiny to date. Here we use high-resolution parsec-scale hydrodynamical simulations to study jet propagation and its effect on BH accretion onto 100 ${\rm M_\odot}$ BHs in the dense, low-metallicity gas expected in early protogalaxies. As the jet propagates, it shocks the surrounding gas and forms a jet cocoon. The cocoon consists of a rapidly-cooling cold phase at the interface with the background gas and an over-pressured subsonic phase of reverse shock-heated gas filling the cocoon interior. We systematically vary the background gas density and temperature, BH feedback efficiency, and the jet model. We found that the jet cocoon width roughly follows a scaling derived by assuming momentum conservation in the jet propagation direction, and energy conservation in the lateral directions. Depending on the assumed gas and jet properties, the cocoon either stays elongated out to a large radius or isotropizes before reaching the Bondi radius, forming a nearly spherical bubble. Lower jet velocities and higher background gas densities result in self-regulation to higher momentum fluxes and elongated cocoons. In all cases, the outward momentum flux of the cocoon balances the inward momentum flux of the inflowing gas near the Bondi radius, which ultimately regulates BH accretion. The larger the distance the jet cocoon reaches, the longer the variability timescale of the BH accretion rate. Overall, the average accretion rate always remains below the Bondi rate, and exceeds the Eddington rate only if the ambient medium is dense and cold, and/or the jet is weak. We derive the combination of jet and ambient gas parameters yielding super-Eddington growth. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: 16 pages, 11 figures

Radiative feedback on supermassive star formation: the massive end of the Population III initial mass function

Authors: Daisuke Toyouchi, Kohei Inayoshi, Wenxiu Li, Zoltán Haiman, Rolf Kuiper

Abstract: Supermassive stars (SMSs) with masses of $M_\ast \simeq 10^4$--$10^5~{\rm M_\odot}$ are invoked as possible seeds of high-redshift supermassive black holes, but it remains under debate whether their protostar indeed acquires sufficient mass via gas accretion overcoming radiative feedback. We investigate protostellar growth in dynamically heated atomic-cooling haloes (ACHs) found in recent cosmolog… ▽ More Supermassive stars (SMSs) with masses of $M_\ast \simeq 10^4$--$10^5~{\rm M_\odot}$ are invoked as possible seeds of high-redshift supermassive black holes, but it remains under debate whether their protostar indeed acquires sufficient mass via gas accretion overcoming radiative feedback. We investigate protostellar growth in dynamically heated atomic-cooling haloes (ACHs) found in recent cosmological simulations, performing three-dimensional radiation hydrodynamical (RHD) simulations that consider stellar evolution under variable mass accretion. We find that one of the ACHs feeds the central protostar at rates exceeding a critical value, above which the star evolves in a cool bloating phase and hardly produces ionizing photons. Consequently, the stellar mass reaches $M_\ast \gtrsim 10^4~{\rm M_\odot}$ unimpeded by radiative feedback. In the other ACH, where the mass supply rate is lower, the star spends most of its life as a hot main-sequence star, emitting intense ionizing radiation. Then, the stellar mass growth is terminated around $500~{\rm M_\odot}$ by photoevaporation of the circumstellar disk. A series of our RHD simulations provide a formula of the final stellar mass determined either by stellar feedback or their lifetime as a function of the mass supply rate from the parent cloud in the absence of stellar radiation. Combining the results with the statistical properties of SMS-forming clouds in high-redshift quasar progenitor haloes, we construct a top-heavy mass distribution of primordial stars over $M_\ast \simeq 100$--$10^5~{\rm M_\odot}$, approximately following a power-law spectrum of $\propto M_\ast^{-1.3}$ with a steeper decline at $M_\ast \gtrsim 2 \times 10^4~{\rm M_\odot}$. Their massive BH remnants would be further fed via the dense debris disk, powering "milli-quasars" with a bolometric luminosity of $L_{\rm bol}~\gtrsim~10^{43}~{\rm erg~s^{-1}}$. △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: 16 pages, 12 figures, and 3 tables, submitted to MNRAS

How long do high-redshift massive black hole seeds remain outliers in black hole vs. host galaxy relations?

Authors: Matthew T. Scoggins, Zoltán Haiman, John H. Wise

Abstract: The existence of $10^9\ {\rm M_\odot}$ supermassive black holes (SMBHs) within the first billion years of the universe remains a puzzle in our conventional understanding of black hole formation and growth. Several suggested formation pathways for these SMBHs lead to a heavy seed, with an initial black hole mass of $10^4-10^6~{\rm M_\odot}$. This can lead to an overly massive BH galaxy (OMBG), whos… ▽ More The existence of $10^9\ {\rm M_\odot}$ supermassive black holes (SMBHs) within the first billion years of the universe remains a puzzle in our conventional understanding of black hole formation and growth. Several suggested formation pathways for these SMBHs lead to a heavy seed, with an initial black hole mass of $10^4-10^6~{\rm M_\odot}$. This can lead to an overly massive BH galaxy (OMBG), whose nuclear black hole's mass is comparable to or even greater than the surrounding stellar mass: the black hole to stellar mass ratio is $M_{\rm bh}/M_* \gg 10^{-3}$, well in excess of the typical values at lower redshift. We investigate how long these newborn BHs remain outliers in the $M_{\rm bh}-M_{*}$ relation, by exploring the subsequent evolution of two OMBGs previously identified in the \texttt{Renaissance} simulations. We find that both OMBGs have $M_{\rm bh}/M_* > 1$ during their entire life, from their birth at $z\approx 15$ until they merge with much more massive haloes at $z\approx 8$. We find that the OMBGs are spatially resolvable from their more massive, $10^{11}~{\rm M_\odot}$, neighboring haloes until their mergers are complete at $z\approx 8$. This affords a window for future observations with {\it JWST} and sensitive X-ray telescopes to diagnose the heavy-seed scenario, by detecting similar OMBGs and establishing their uniquely high black hole-to-stellar mass ratio. △ Less

Submitted 13 December, 2022; v1 submitted 19 May, 2022; originally announced May 2022.

Comments: 15 pages, 10 figures. Accepted in MNRAS

Supercritical growth pathway to overmassive black holes at cosmic dawn: coevolution with massive quasar hosts

Authors: Haojie Hu, Kohei Inayoshi, Zoltán Haiman, Wenxiu Li, Eliot Quataert, Rolf Kuiper

Abstract: Observations of the most luminous quasars at high redshifts ($z > 6$) have revealed that the largest supermassive black holes (SMBHs) at those epochs tend to be substantially overmassive relative to their host galaxies compared to the local relations, suggesting they experienced rapid early growth phases. We propose an assembly model for the SMBHs that end up in rare massive… ▽ More Observations of the most luminous quasars at high redshifts ($z > 6$) have revealed that the largest supermassive black holes (SMBHs) at those epochs tend to be substantially overmassive relative to their host galaxies compared to the local relations, suggesting they experienced rapid early growth phases. We propose an assembly model for the SMBHs that end up in rare massive $\sim10^{12}~M_{\odot}$ host halos at $z \sim 6-7$, applying a kinetic feedback prescription for BHs accreting above the Eddington rate, provided by radiation hydrodynamic simulations for the long-term evolution of the accretion-flow structure. The large inflow rates into these halos during their assembly enable the formation of $>10^9~M_{\odot}$ SMBHs by $z \sim 6$, even starting from stellar-mass seeds at $z \sim 30$, and even in the presence of outflows that reduce the BH feeding rate, especially at early times. This mechanism also naturally yields a high BH-to-galaxy mass ratio of $> 0.01$ before the SMBH mass reaches $M_{\rm BH} > 10^9~M_{\odot}$ by $z \sim 6$. These fast-growing SMBH progenitors are bright enough to be detected by upcoming observations with the James Webb Space Telescope over a wide range of redshift ($7 < z < 15$), regardless of how they were seeded. △ Less

Submitted 30 June, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: 4 figures, accepted for publication in ApJ

The science case and challenges of space-borne sub-millimeter interferometry

Authors: Leonid I. Gurvits, Zsolt Paragi, Ricardo I. Amils, Ilse van Bemmel, Paul Boven, Viviana Casasola, John Conway, Jordy Davelaar, M. Carmen Díez-González, Heino Falcke, Rob Fender, Sándor Frey, Christian M. Fromm, Juan D. Gallego-Puyol, Cristina García-Miró, Michael A. Garrett, Marcello Giroletti, Ciriaco Goddi, José L. Gómez, Jeffrey van der Gucht, José Carlos Guirado, Zoltán Haiman, Frank Helmich, Ben Hudson, Elizabeth Humphreys , et al. (29 additional authors not shown)

Abstract: Ultra-high angular resolution in astronomy has always been an important vehicle for making fundamental discoveries. Recent results in direct imaging of the vicinity of the supermassive black hole in the nucleus of the radio galaxy M87 by the millimeter VLBI system Event Horizon Telescope and various pioneering results of the Space VLBI mission RadioAstron provided new momentum in high angular reso… ▽ More Ultra-high angular resolution in astronomy has always been an important vehicle for making fundamental discoveries. Recent results in direct imaging of the vicinity of the supermassive black hole in the nucleus of the radio galaxy M87 by the millimeter VLBI system Event Horizon Telescope and various pioneering results of the Space VLBI mission RadioAstron provided new momentum in high angular resolution astrophysics. In both mentioned cases, the angular resolution reached the values of about 10-20 microrcseconds. Further developments toward at least an order of magnitude "sharper" values are dictated by the needs of astrophysical studies and can only be achieved by placing millimeter and submillimeter wavelength interferometric systems in space. A concept of such the system, called Terahertz Exploration and Zooming-in for Astrophysics (THEZA), has been proposed in the framework of the ESA Call for White Papers for the Voayage 2050 long term plan in 2019. In the current paper we discuss several approaches for addressing technological challenges of the THEZA concept. In particular, we consider a novel configuration of a space-borne millimeter/sub-millimeter antenna which might resolve several bottlenecks in creating large precise mechanical structures. The paper also presents an overview of prospective space-qualified technologies of low-noise analogue front-end instrumentation for millimeter/sub-millimeter telescopes, data handling and processing. The paper briefly discusses approaches to the interferometric baseline state vector determination and synchronisation and heterodyning system. In combination with the original ESA Voyage 2050 White Paper, the current work sharpens the case for the next generation microarcsceond-level imaging instruments and provides starting points for further in-depth technology trade-off studies. △ Less

Submitted 27 April, 2022; v1 submitted 19 April, 2022; originally announced April 2022.

Comments: The paper on 29 pages contains 12 figures, accepted for publication in Acta Astronautica, based on the presentation at the 72nd IAC, Dubai, UAE, October 2021

Journal ref: Acta Astronautica 196 (2022) 314

Long-term evolution of supercritical black hole accretion with outflows: a subgrid feedback model for cosmological simulations

Authors: Haojie Hu, Kohei Inayoshi, Zoltán Haiman, Eliot Quataert, Rolf Kuiper

Abstract: We study the long-term evolution of the global structure of axisymmetric accretion flows onto a black hole (BH) at rates substantially higher than the Eddington value ($\dot{M}_{\rm Edd}$), performing two-dimensional hydrodynamical simulations with and without radiative diffusion. In the high-accretion optically-thick limit, where the radiation energy is efficiently trapped within the inflow, the… ▽ More We study the long-term evolution of the global structure of axisymmetric accretion flows onto a black hole (BH) at rates substantially higher than the Eddington value ($\dot{M}_{\rm Edd}$), performing two-dimensional hydrodynamical simulations with and without radiative diffusion. In the high-accretion optically-thick limit, where the radiation energy is efficiently trapped within the inflow, the accretion flow becomes adiabatic and comprises of turbulent gas in the equatorial region and strong bipolar outflows. As a result, the mass inflow rate decreases toward the center as $\dot{M}_{\rm in}\propto r^{p}$ with $p\sim 0.5-0.7$ and a small fraction of the inflowing gas feeds the nuclear BH. Thus, super-Eddington accretion is sustained only when a larger amount of gas is supplied from larger radii at $> 100-1000~\dot{M}_{\rm Edd}$. The global structure of the flow settles down to a quasi-steady state in millions of the orbital timescale at the BH event horizon, which is $> 10-100$ times longer than that addressed in previous (magneto-)RHD simulation studies. Energy transport via radiative diffusion accelerates the outflow near the poles in the inner region but does not change the overall properties of the accretion flow compared to the cases without diffusion. Based on our simulation results, we provide a mechanical feedback model for super-Eddington accreting BHs. This can be applied as a sub-grid model in large-scale cosmological simulations that do not sufficiently resolve galactic nuclei, and to the formation of the heaviest gravitational-wave sources via accretion in dense environments. △ Less

Submitted 1 June, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 17 pages, 13 figures, 1 table (re-submitted to APJ, after revision)

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Comparing weak lensing peak counts in baryonic correction models to hydrodynamical simulations

Authors: Max E. Lee, Tianhuan Lu, Zoltán Haiman, Jia Liu, Ken Osato

Abstract: Next-generation weak lensing (WL) surveys, such as by the Vera Rubin Observatory's LSST, the $\textit{Roman}$ Space Telescope, and the $\textit{Euclid}$ space mission, will supply vast amounts of data probing small, highly nonlinear scales. Extracting information from these scales requires higher-order statistics and the controlling of related systematics such as baryonic effects. To account for b… ▽ More Next-generation weak lensing (WL) surveys, such as by the Vera Rubin Observatory's LSST, the $\textit{Roman}$ Space Telescope, and the $\textit{Euclid}$ space mission, will supply vast amounts of data probing small, highly nonlinear scales. Extracting information from these scales requires higher-order statistics and the controlling of related systematics such as baryonic effects. To account for baryonic effects in cosmological analyses at reduced computational cost, semi-analytic baryonic correction models (BCMs) have been proposed. Here, we study the accuracy of BCMs for WL peak counts, a well studied, simple, and effective higher-order statistic. We compare WL peak counts generated from the full hydrodynamical simulation IllustrisTNG and a baryon-corrected version of the corresponding dark matter-only simulation IllustrisTNG-Dark. We apply galaxy shape noise expected at the depths reached by DES, KiDS, HSC, LSST, $\textit{Roman}$, and $\textit{Euclid}$. We find that peak counts in BCMs are (i) accurate at the percent level for peaks with $\mathrm{S/N}<4$, (ii) statistically indistinguishable from IllustrisTNG in most current and ongoing surveys, but (iii) insufficient for deep future surveys covering the largest solid angles, such as LSST and $\textit{Euclid}$. We find that BCMs match individual peaks accurately, but underpredict the amplitude of the highest peaks. We conclude that existing BCMs are a viable substitute for full hydrodynamical simulations in cosmological parameter estimation from beyond-Gaussian statistics for ongoing and future surveys with modest solid angles. For the largest surveys, BCMs need to be refined to provide a more accurate match, especially to the highest peaks. △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: 12 pages, 10 figures

Self-lensing flares from black hole binaries II: observing black hole shadows via light-curve tomography

Authors: Jordy Davelaar, Zoltán Haiman

Abstract: Supermassive black hole (BH) binaries are thought to produce self-lensing flares (SLF) when the two BHs are aligned with the line-of-sight. If the binary orbit is observed nearly edge-on, we find a distinct feature in the light curve imprinted by the BH shadow from the lensed BH. We study this feature by ray-tracing in a binary model and predict that 1\% of the current binary candidates could show… ▽ More Supermassive black hole (BH) binaries are thought to produce self-lensing flares (SLF) when the two BHs are aligned with the line-of-sight. If the binary orbit is observed nearly edge-on, we find a distinct feature in the light curve imprinted by the BH shadow from the lensed BH. We study this feature by ray-tracing in a binary model and predict that 1\% of the current binary candidates could show this feature. Our BH tomography method proposed here could make it possible to extract BH shadows that are spatially unresolvable by high-resolution VLBI. △ Less

Submitted 27 December, 2021; v1 submitted 10 December, 2021; originally announced December 2021.

Comments: 6 pages, 3 figures, 1 table. Submitted to journal

Self-lensing flares from black hole binaries I: general-relativistic ray tracing of black hole binaries

Authors: Jordy Davelaar, Zoltán Haiman

Abstract: The self-lensing of a massive black hole binary (MBHB), which occurs when the two BHs are aligned close to the line of sight, is expected to produce periodic, short-duration flares. Here we study the shapes of self-lensing flares (SLFs) via general-relativistic ray tracing in a superimposed binary BH metric, in which the emission is generated by geometrically thin accretion flows around each compo… ▽ More The self-lensing of a massive black hole binary (MBHB), which occurs when the two BHs are aligned close to the line of sight, is expected to produce periodic, short-duration flares. Here we study the shapes of self-lensing flares (SLFs) via general-relativistic ray tracing in a superimposed binary BH metric, in which the emission is generated by geometrically thin accretion flows around each component. The suite of models covers eccentric binary orbits, black hole spins, unequal mass binaries, and different emission model geometries. We explore the above parameter space, and report how the light curves change as a function of, e.g., binary separation, inclination, and eccentricity. We also compare our light curves to those in the microlensing approximation, and show how strong deflections, as well as time-delay effects, change the size and shape of the SLF. If gravitational waves (GWs) from the inspiraling MBHB are observed by LISA, SLFs can help securely identify the source and localizing it on the sky, and to constrain the graviton mass by comparing the phasing of the SLFs and the GWs. Additionally, when these systems are viewed edge-on the SLF shows a distinct dip that can be directly correlated with the BH shadow size. This opens a new way to measure BH shadow sizes in systems that are unresolvable by current VLBI facilities. △ Less

Submitted 27 December, 2021; v1 submitted 10 December, 2021; originally announced December 2021.

Comments: 16 pages, 15 figures, 1 table. Submitted to journal

Multi-Messenger Constraints on Magnetic Fields in Merging Black Hole-Neutron Star Binaries

Authors: Daniel J. D'Orazio, Zoltan Haiman, Janna Levin, Johan Samsing, Alejandro Vigna-Gomez

Abstract: The LIGO-Virgo-KAGRA Collaboration recently detected gravitational waves (GWs) from the merger of black-hole-neutron-star (BHNS) binary systems GW200105 and GW200115. No coincident electromagnetic (EM) counterparts were detected. While the mass ratio and BH spin in both systems were not sufficient to tidally disrupt the NS outside of the BH event horizon, other, magnetospheric mechanisms for EM em… ▽ More The LIGO-Virgo-KAGRA Collaboration recently detected gravitational waves (GWs) from the merger of black-hole-neutron-star (BHNS) binary systems GW200105 and GW200115. No coincident electromagnetic (EM) counterparts were detected. While the mass ratio and BH spin in both systems were not sufficient to tidally disrupt the NS outside of the BH event horizon, other, magnetospheric mechanisms for EM emission exist in this regime and depend sensitively on the NS magnetic field strength. Combining GW measurements with EM flux upper limits, we place upper limits on the NS surface magnetic field strength above which magnetospheric emission models would have generated an observable EM counterpart. We consider fireball models powered by the black-hole battery mechanism, where energy is output in gamma-rays over $\lesssim1$~second. Consistency with no detection by Fermi-GBM or INTEGRAL SPI-ACS constrains the NS surface magnetic field to $\lesssim10^{15}$~G. Hence, joint GW detection and EM upper limits rule out the theoretical possibility that the NSs in GW200105 and GW200115, and the putative NS in GW190814, retain $\gtrsim10^{15}$~G dipolar magnetic fields until merger. They also rule out formation scenarios where strongly magnetized magnetars quickly merge with BHs. We alternatively rule out operation of the BH-battery powered fireball mechanism in these systems. This is the first multi-messenger constraint on NS magnetic fields in BHNS systems and a novel approach to probe fields at this point in NS evolution. This demonstrates the constraining power that multi-messenger analyses of BHNS mergers have on BHNS formation scenarios, the magnetic-field evolution in NSs, and the physics of BHNS magnetospheric interactions. △ Less

Submitted 3 December, 2021; originally announced December 2021.

Comments: Comments welcome

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.

Comments: 19 pages, 9 figures, revised version after the first referee report

Multi-band light curves from eccentric accreting supermassive black hole binaries

Authors: John Ryan Westernacher-Schneider, Jonathan Zrake, Andrew MacFadyen, Zoltán Haiman

Abstract: We use long-run, high-resolution hydrodynamics simulations to compute the multi-wavelength light curves (LCs) from thermal disk emission around accreting equal-mass supermassive black hole (BH) binaries, with a focus on revealing binary eccentricity. LCs are obtained by modeling the disk thermodynamics with an adiabatic equation of state, a local blackbody cooling prescription, and corrections to… ▽ More We use long-run, high-resolution hydrodynamics simulations to compute the multi-wavelength light curves (LCs) from thermal disk emission around accreting equal-mass supermassive black hole (BH) binaries, with a focus on revealing binary eccentricity. LCs are obtained by modeling the disk thermodynamics with an adiabatic equation of state, a local blackbody cooling prescription, and corrections to approximate the effects of radiation pressure. We find that modulation of multi-band LCs on the orbital time scale are generally in-phase (to within $\sim\,$2% of a binary orbital period), but they contain pulse substructure in the time domain that is not necessarily reflected in BH accretion rates $\dot M$. We thus predict that binary-hosting AGN will exhibit highly correlated, in-phase, periodic brightness modulations in their low-energy disk emission. However, detectability of these modulations in multi-wavelength observing campaigns could be seriously compromised because observed stochastic variability in AGNs typically has a higher amplitude than our proposed signal. It is possible that observations over temporal baselines of many binary periods may make the signal more prominent, but this would need to be analyzed carefully. If jet emission is predicted by $\dot{M}$, then we predict a weaker correlation with low-energy disk emission due to the differing sub-peak structure. For the binary parameters we explore, we show that LC variability due to hydrodynamics likely dominates Doppler brightening for all equal-mass binaries with disk Mach numbers $\lesssim 20$. A promising signature of eccentricity is weak or absent "lump" periodicity. We find hints that a significant lag exists between $\dot{M}$ and low-energy disk emission for circular binaries, but they are in-phase for eccentric binaries, which might explain some "orphan" blazar flares with no $γ$-ray counterpart. △ Less

Submitted 17 October, 2022; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: v6: accepted version (PRD)

How Binaries Accrete: Hydrodynamic Simulations with Passive Tracer Particles

Authors: Christopher Tiede, Jonathan Zrake, Andrew MacFadyen, Zoltan Haiman

Abstract: Linear analysis of gas flows around orbiting binaries suggests that a centrifugal barrier ought to clear a low-density cavity around the binary and inhibit mass transfer onto it. Modern hydrodynamics simulations have confirmed the low-density cavity, but show that any mass flowing from large scales into the circumbinary disk is eventually transferred onto the binary components. Even though many nu… ▽ More Linear analysis of gas flows around orbiting binaries suggests that a centrifugal barrier ought to clear a low-density cavity around the binary and inhibit mass transfer onto it. Modern hydrodynamics simulations have confirmed the low-density cavity, but show that any mass flowing from large scales into the circumbinary disk is eventually transferred onto the binary components. Even though many numerical studies confirm this picture, it is still not understood precisely how gas parcels overcome the centrifugal barrier and ultimately accrete. We present a detailed analysis of the binary accretion process, using an accurate prescription for evolving grid-based hydrodynamics with Lagrangian tracer particles that track the trajectories of individual gas parcels. We find that binary accretion can be described in four phases: (1) gas is viscously transported through the circumbinary disk up to the centrifugal barrier at the cavity wall, (2) the cavity wall is tidally distorted into accretion streams consisting of near-ballistic gas parcels on eccentric orbits, (3) the portion of each stream moving inwards of an ``accretion horizon'' radius $\bar r \simeq a$ -- the radius beyond which no material is returned to the cavity wall -- becomes bound to a minidisk orbiting an individual binary component, and (4) the minidisk gas accretes onto the binary component through the combined effect of viscous and tidal stresses. △ Less

Submitted 13 June, 2022; v1 submitted 8 November, 2021; originally announced November 2021.

Comments: Published in ApJ; 18 pages, 16 figures

Journal ref: The Astrophysical Journal, Volume 932, Number 1, 2022 June 13

Simultaneously constraining cosmology and baryonic physics via deep learning from weak lensing

Authors: Tianhuan Lu, Zoltán Haiman, José Manuel Zorrilla Matilla

Abstract: Ongoing and planned weak lensing (WL) surveys are becoming deep enough to contain information on angular scales down to a few arcmin. To fully extract information from these small scales, we must capture non-Gaussian features in the cosmological WL signal while accurately accounting for baryonic effects. In this work, we account for baryonic physics via a baryonic correction model that modifies th… ▽ More Ongoing and planned weak lensing (WL) surveys are becoming deep enough to contain information on angular scales down to a few arcmin. To fully extract information from these small scales, we must capture non-Gaussian features in the cosmological WL signal while accurately accounting for baryonic effects. In this work, we account for baryonic physics via a baryonic correction model that modifies the matter distribution in dark matter-only $N$-body simulations, mimicking the effects of galaxy formation and feedback. We implement this model in a large suite of ray-tracing simulations, spanning a grid of cosmological models in $Ω_\mathrm{m}-σ_8$ space. We then develop a convolutional neural network (CNN) architecture to learn and constrain cosmological and baryonic parameters simultaneously from the simulated WL convergence maps. We find that in a Hyper-Suprime Cam (HSC)-like survey, our CNN achieves a 1.7$\times$ tighter constraint in $Ω_\mathrm{m}-σ_8$ space ($1σ$ area) than the power spectrum and 2.1$\times$ tighter than the peak counts, showing that the CNN can efficiently extract non-Gaussian cosmological information even while marginalizing over baryonic effects. When we combine our CNN with the power spectrum, the baryonic effects degrade the constraint in $Ω_\mathrm{m}-σ_8$ space by a factor of 2.4, compared to the much worse degradation by a factor of 4.7 or 3.7 from either method alone. △ Less

Submitted 17 January, 2022; v1 submitted 22 September, 2021; originally announced September 2021.

Comments: 12 pages, 6 figures

Nucleosynthetic signatures of primordial origin around supermassive black holes

Authors: Phoebe Upton Sanderbeck, Simeon Bird, Zoltan Haiman

Abstract: If primordial black holes (PBHs) seeded the supermassive black holes (SMBHs) at the centers of high-redshift quasars, then the gas surrounding these black holes may reveal nucleosynthetic clues to their primordial origins. We present predictions of altered primordial abundances around PBHs massive enough to seed SMBHs at z~6-7.5. We find that if PBHs with initial masses of ~10^5 M$_{\odot}$ are re… ▽ More If primordial black holes (PBHs) seeded the supermassive black holes (SMBHs) at the centers of high-redshift quasars, then the gas surrounding these black holes may reveal nucleosynthetic clues to their primordial origins. We present predictions of altered primordial abundances around PBHs massive enough to seed SMBHs at z~6-7.5. We find that if PBHs with initial masses of ~10^5 M$_{\odot}$ are responsible for such SMBHs, they may produce primordial Deuterium and Helium fractions enhanced by >~ 10%, and Lithium abundance depleted by >~10%, at distances of up to ~ a comoving kiloparsec away from the black hole after decoupling. We estimate that ~ 10^8 M$_{\odot}$ of gas is enhanced (or depleted) by at least one percent. Evidence of these modified primordial Deuterium, Helium, and Lithium abundances could still be present if this circum-PBH gas remains unaccreted by the SMBH and in or near the host galaxies of high-redshift quasars. Measuring the abundance anomalies will be challenging, but could offer a novel way to reveal the primordial origin of such SMBH seeds. △ Less

Submitted 17 November, 2021; v1 submitted 10 September, 2021; originally announced September 2021.

Comments: 7 pages, 4 figures

Journal ref: Phys. Rev. D 104, 103022 (2021)

Cosmological Constraints from Weak Lensing Peaks: Can Halo Models Accurately Predict Peak Counts?

Authors: Alina Sabyr, Zoltán Haiman, José Manuel Zorrilla Matilla, Tianhuan Lu

Abstract: In order to extract full cosmological information from next-generation large and high-precision weak lensing (WL) surveys (e.g. Euclid, Roman, LSST), higher-order statistics that probe the small-scale, non-linear regime of large scale structure (LSS) need to be utilized. WL peak counts, which trace overdensities in the cosmic web, are one promising and simple statistic for constraining cosmologica… ▽ More In order to extract full cosmological information from next-generation large and high-precision weak lensing (WL) surveys (e.g. Euclid, Roman, LSST), higher-order statistics that probe the small-scale, non-linear regime of large scale structure (LSS) need to be utilized. WL peak counts, which trace overdensities in the cosmic web, are one promising and simple statistic for constraining cosmological parameters. The physical origin of WL peaks have previously been linked to dark matter halos along the line of sight and this peak-halo connection has been used to develop various semi-analytic halo-based models for predicting peak counts. Here, we study the origin of WL peaks and the effectiveness of halo-based models for WL peak counts using a suite of ray-tracing N-body simulations. We compare WL peaks in convergence maps from the full simulations to those in maps created from only particles associated with halos -- the latter playing the role of a "perfect" halo model. We find that while halo-only contributions are able to replicate peak counts qualitatively well, halos do not explain all WL peaks. Halos particularly underpredict negative peaks, which are associated with local overdensities in large-scale underdense regions along the line of sight. In addition, neglecting non-halo contributions to peaks counts leads to a significant bias on the parameters ($Ω_{\rm m}$, $σ_{8}$) for surveys larger than $\geq$ 100 deg$^{2}$. We conclude that other elements of the cosmic web, outside and far away from dark matter halos, need to be incorporated into models of WL peaks in order to infer unbiased cosmological constraints. △ Less

Submitted 28 February, 2022; v1 submitted 1 September, 2021; originally announced September 2021.

Comments: 15 pages, 15 figures, updated to match PRD version

Journal ref: Phys. Rev. D 105, 023505 (2022)

The Effect of Mission Duration on LISA Science Objectives

Authors: Pau Amaro Seoane, Manuel Arca Sedda, Stanislav Babak, Christopher P. L. Berry, Emanuele Berti, Gianfranco Bertone, Diego Blas, Tamara Bogdanović, Matteo Bonetti, Katelyn Breivik, Richard Brito, Robert Caldwell, Pedro R. Capelo, Chiara Caprini, Vitor Cardoso, Zack Carson, Hsin-Yu Chen, Alvin J. K. Chua, Irina Dvorkin, Zoltan Haiman, Lavinia Heisenberg, Maximiliano Isi, Nikolaos Karnesis, Bradley J. Kavanagh, Tyson B. Littenberg , et al. (16 additional authors not shown)

Abstract: The science objectives of the LISA mission have been defined under the implicit assumption of a 4 yr continuous data stream. Based on the performance of LISA Pathfinder, it is now expected that LISA will have a duty cycle of $\approx 0.75$, which would reduce the effective span of usable data to 3 yr. This paper reports the results of a study by the LISA Science Group, which was charged with asses… ▽ More The science objectives of the LISA mission have been defined under the implicit assumption of a 4 yr continuous data stream. Based on the performance of LISA Pathfinder, it is now expected that LISA will have a duty cycle of $\approx 0.75$, which would reduce the effective span of usable data to 3 yr. This paper reports the results of a study by the LISA Science Group, which was charged with assessing the additional science return of increasing the mission lifetime. We explore various observational scenarios to assess the impact of mission duration on the main science objectives of the mission. We find that the science investigations most affected by mission duration concern the search for seed black holes at cosmic dawn, as well as the study of stellar-origin black holes and of their formation channels via multi-band and multi-messenger observations. We conclude that an extension to 6 yr of mission operations is recommended. △ Less

Submitted 12 January, 2022; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 50 pages, 19 figures, 5 tables. Matches version published in GERG

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.

Comments: 7 pages, 1 figure