Abstract The broad emission lines (BELs) emitted by active galactic nuclei respond to variations in the ionizing continuum emission from the accretion disk surrounding the central supermassive black hole (SMBH). This reverberation response provides insights into the structure and dynamics of the broad-line region (BLR). In 2024, we introduced a new forward-modeling tool, the Broad Emission Line Mapping Code (BELMAC), which simulates the velocity-resolved reverberation response of the BLR to an input light curve. In this work, we describe a new version of BELMAC, which uses photoionization models to calculate the cloud luminosities for selected BELs. We investigated the reverberation responses of Hα, Hβ, MgIIλ2800, and CIVλ1550 for models representing a disk-like BLR with Keplerian rotation, radiatively driven outflows, and inflows. The line responses generally provide a good indication of the respective luminosity-weighted radii. However, there are situations where the BLR exhibits a negative response to the driving continuum, causing overestimates of the luminosity-weighted radius. The virial mass derived from the models can differ dramatically from the actual SMBH mass, depending mainly on the disk inclination and velocity field. In single-zone models, the BELs exhibit similar responses and profile shapes; two-zone models, such as a Keplerian disk and a biconical outflow, can reproduce observed differences between high- and low-ionization lines. Radial flows produce asymmetric line profile shapes due to both anisotropic cloud emission and electron scattering in an intercloud medium. These competing attenuation effects complicate the interpretation of profile asymmetries.
more »
« less
A New Diagnostic Method for Distinguishing Binary from Single Supermassive Black Holes Using Broad Emission Lines
During cosmic timescales, supermassive binary black holes (SMBBHs) form by galaxy merg- ers, where each galaxy hosts a supermassive black hole (SMBH) at its center. By studying SMBBHs, we can gain insights into galaxy evolution and black hole growth. However, the typical separation between black holes in SMBBHs is usually below 1 pc, making them dif- ficult to resolve using direct imaging or photometry. To be able to distinguish binary black holes (BBHs) from typical AGN powered by single black holes (SBHs), we conducted this research project to develop a new diagnostic method to identify a unique feature of SMBBHs, which can be used to distinguish AGN powered by BBHs from those powered by SBHs. The basic idea of this method is that BBHs have different configurations compared with those of SBHs, such as the circumbinary disk enveloping the whole binary system, the mini- disks around each black hole, and the streams between the circumbinary disk and minidisks. It is these different configurations of accretion disks of black holes that lead to the difference in the spectral energy distributions (SEDs), which in turn will differently photoionize the broad line region (BLR) and produce different strengths of broad emission lines. Thus, it is these differences in strengths of broad emission lines that we expect to see between two different configurations of black holes, either binary or single black holes. By identifying these differences in line strengths, we can distinguish between binary and single black holes. In order to achieve this goal of distinguishing BBHs from SBHs, we organized our project in the steps below: (1) obtained BBH SEDs from Gutiérrez et al., 2022, (2) generated SBH SEDs using XSPEC modeling code OPTXAGNF, (3) produced single-cloud models, which represent the BLR, input BBH and SBH SEDs into the models, and simulated the single0cloud response with a photoionization code Cloudy, (4) built cloud-ensemble models, which represent a more realistic BLR, input BBH and SBH SEDs into the models, and simulated the emission-line response within those clouds using a broad emission line mapping code BELMAC. In step (3), we have explored the differences in line ratios between BBH and SBH for a representative single-cloud photoionization model. The emission lines we used here are: Si IV λ1400Å, C III] λ1909Å, C IV λ1549Å, Mg II λ2798Å, and Lyα λ1216Å. It turned i out that differences do exist, but they are too small to be identified in observational data. Furthermore, we have investigated the line equivalent widths predicted by SBH and BBH models respectively. By doing so, we found some apparent differences between BBHs and SBHs in some specific emission lines: Lyα, CIV, and Hα. However, these differences vanish at the highest mass of black holes (109M⊙). In step (4), we continued the investigation of the equivalent width between SBH and BBH BLR cloud-ensemble models and found that some emission lines show the difference between BBHs and SBHs, such as CIV in the case of BH mass 107M⊙, U ∝ r−2, and log n/1 cm−3 = 10.5, 11.0. For the highest mass in the case of U ∝ r−2, the results are similar to the one in single-cloud models: no difference is shown between BBHs and SBHs across all emission lines. Most importantly, in this step, we found that for CIV λ1549 Å in the case of a black hole mass of 107M⊙, s = −2, and log n/1 cm−3 = 10.5, only the SBH EW falls inside the range of the observed range in SDSS DR7 Quasar Catalog while the BBH EW falls outside the range and becomes an outlier. This is what we want to find to distinguish BBHs from SBHs in the observational data.
more »
« less
- Award ID(s):
- 2009508
- PAR ID:
- 10644402
- Publisher / Repository:
- Rochester Institute of Technology
- Date Published:
- Format(s):
- Medium: X
- Institution:
- Rochester Institute of Technology
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Fast empirical models of the broad emission line region (BLR) are a powerful tool to interpret velocity-resolved reverberation mapping (RM) data, estimate the mass of the supermassive black holes, and gain insight into its geometry and kinematics. Much of the effort so far has been devoted to describing the emissivity of one emission line at a time. We present here an alternative approach aimed at describing the underlying BLR gas distribution, by exploiting simple numerical recipes to connect it with emissivity. This approach is a step toward describing multiple emission lines originating from the same gas and allows us to clarify some issues related to the interpretation of RM data. We illustrate this approach—implemented in the codeCARAMEL-gas—using three data sets covering the Hβemission line (Mrk 50, Mrk 1511, Arp 151) that have been modeled using the emissivity-based version of the code. As expected, we find differences in the parameters describing the BLR gas and emissivity distribution, but the emissivity-weighted lag measurements and all other model parameters including black hole mass and overall BLR morphology and kinematics are consistent with the previous measurements. We also model the Hαemission line for Arp 151 using both the gas- and emissivity-based BLR models. We find ionization stratification in the BLR with Hαarising at larger radii than Hβ, while all other model parameters are consistent within the uncertainties.more » « less
-
Abstract There are few observed high-mass X-ray binaries (HMXBs) that harbor massive black holes (BHs), and none are likely to result in a binary black hole (BBH) that merges within a Hubble time; however, we know that massive merging BBHs exist from gravitational-wave (GW) observations. We investigate the role that X-ray and GW observational selection effects play in determining the properties of their respective detected binary populations. We find that, as a result of selection effects, detectable HMXBs and detectable BBHs form at different redshifts and metallicities, with detectable HMXBs forming at much lower redshifts and higher metallicities than detectable BBHs. We also find disparities in the mass distributions of these populations, with detectable merging BBH progenitors pulling to higher component masses relative to the full detectable HMXB population. Fewer than 3% of detectable HMXBs host BHs >35M⊙in our simulated populations. Furthermore, we find the probability that a detectable HMXB will merge as a BBH system within a Hubble time is ≃0.6%. Thus, it is unsurprising that no currently observed HMXB is predicted to form a merging BBH with high probability.more » « less
-
ABSTRACT Galactic nuclei are promising sites for stellar origin black hole (BH) mergers, as part of merger hierarchies in deep potential wells. We show that binary black hole (BBH) merger rates in active galactic nuclei (AGNs) should always exceed merger rates in quiescent galactic nuclei (nuclear star clusters, NSCs) around supermassive black holes (SMBHs) without accretion discs. This is primarily due to average binary lifetimes in AGNs that are significantly shorter than those in NSCs. The lifetime difference comes from rapid hardening of BBHs in AGNs, such that their semimajor axes are smaller than the hard–soft boundary of their parent NSC; this contrasts with the large average lifetime to merger for BBHs in NSCs around SMBHs, due to binary ionization mechanisms. Secondarily, merger rates in AGNs are enhanced by gas-driven binary formation mechanisms. Formation of new BHs in AGN discs is a minor contributor to the rate differences. With the gravitational wave detection of several BBHs with at least one progenitor in the upper mass gap, and signatures of dynamical formation channels in the χeff distribution, we argue that AGNs could contribute $$\sim 25{\!-\!}80{{\ \rm per\ cent}}$$ of the LIGO–Virgo measured rate of $$\sim 24\, \rm {Gpc}^{-3} \rm {yr}^{-1}$$.more » « less
-
Abstract Numerical relativity (NR) simulations of binary black hole (BBH) systems provide the most accurate gravitational wave predictions, but at a high computational cost—especially when the black holes have nearly extremal spins (i.e. spins near the theoretical upper limit) or very unequal masses. Recently, the technique of reduced order modeling has enabled the construction of ‘surrogate models’ trained on an existing set of NR waveforms. Surrogate models enable the rapid computation of the gravitational waves emitted by BBHs. Typically these models are used for interpolation to compute gravitational waveforms for BBHs with mass ratios and spins within the bounds of the training set. Because simulations with nearly extremal spins are so technically challenging, surrogate models almost always rely on training sets with only moderate spins. In this paper, we explore how well surrogate models can extrapolate to nearly extremal spins when the training set only includes moderate spins. For simplicity, we focus on one-dimensional surrogate models trained on NR simulations of BBHs with equal masses and equal, aligned spins. We assess the performance of the surrogate models at higher spin magnitudes by calculating the mismatches between extrapolated surrogate model waveforms and NR waveforms, by calculating the differences between extrapolated and NR measurements of the remnant black-hole mass, and by testing how the surrogate model improves as the training set extends to higher spins. We find that while extrapolation in this one-dimensional case is viable for current detector sensitivities, surrogate models for next-generation detectors should use training sets that extend to nearly extremal spins.more » « less
An official website of the United States government

