Periodic signatures in time-domain observations of quasars have been used to search for binary supermassive black holes (SMBHs). These searches, across existing time-domain surveys, have produced several hundred candidates. The general stochastic variability of quasars, however, can masquerade as a false-positive periodic signal, especially when monitoring cadence and duration are limited. In this work, we predict the detectability of binary SMBHs in the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We apply computationally inexpensive sinusoidal curve fits to millions of simulated LSST Deep Drilling Field light curves of both single, isolated quasars and binary quasars. The period and phase of simulated binary signals can generally be disentangled from quasar variability. Binary amplitude is overestimated and poorly recovered for two-thirds of potential binaries due to quasar accretion variability. Quasars with strong intrinsic variability can obscure a binary signal too much for recovery. We also find that the most luminous quasars mimic current binary candidate light curves and their properties: The false-positive rates are 60% for these quasars. The reliable recovery of binary period and phase for a wide range of input binary LSST light curves is promising for multi-messenger characterization of binary SMBHs. However, pure electromagnetic detections of binaries using photometric periodicity with amplitude greater than 0.1 mag will result in samples that are overwhelmed by false positives. This paper represents an important and computationally inexpensive way forward for understanding the true and false-positive rates for binary candidates identified by Rubin.
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
Abstract -
Abstract We explore reprocessing models for a sample of 17 hypervariable quasars, taken from the Sloan Digital Sky Survey Reverberation Mapping project, which all show coordinated optical luminosity hypervariability with amplitudes of factors ≳2 between 2014 and 2020. We develop and apply reprocessing models for quasar light curves in simple geometries that are likely to be representative of quasar inner environments. In addition to the commonly investigated thin-disk model, we include the thick-disk and hemisphere geometries. The thick-disk geometry could, for instance, represent a magnetically elevated disk, whereas the hemisphere model can be interpreted as a first-order approximation for any optically thick out-of-plane material caused by outflows/winds, warped/tilted disks, and so on. Of the 17 quasars in our sample, 11 are best-fitted by a hemisphere geometry, five are classified as thick disks, and both models fail for just one object. We highlight the successes and shortcomings of our thermal reprocessing models in case studies of four quasars that are representative of the sample. While reprocessing is unlikely to explain all of the variability that we observe in quasars, we present our classification scheme as a starting point for revealing the likely geometries of reprocessing for quasars in our sample and hypervariable quasars in general.
-
Abstract We measure the correlation between black hole mass
M BHand host stellar massM *for a sample of 38 broad-line quasars at 0.2 ≲z ≲ 0.8 (median redshiftz med= 0.5). The black hole masses are derived from a dedicated reverberation mapping program for distant quasars, and the stellar masses are derived from two-band optical+IR Hubble Space Telescope imaging. Most of these quasars are well centered within ≲1 kpc from the host galaxy centroid, with only a few cases in merging/disturbed systems showing larger spatial offsets. Our sample spans two orders of magnitude in stellar mass (∼109–1011M ⊙) and black hole mass (∼107–109M ⊙) and reveals a significant correlation between the two quantities. We find a best-fit intrinsic (i.e., selection effects corrected)M BH–M *,hostrelation of , with an intrinsic scatter of dex. Decomposing our quasar hosts into bulges and disks, there is a similarM BH–M *,bulgerelation with slightly larger scatter, likely caused by systematic uncertainties in the bulge–disk decomposition. TheM BH–M *,hostrelation atz med= 0.5 is similar to that in local quiescent galaxies, with negligible evolution over the redshift range probed by our sample. With direct black hole masses from reverberation mapping and the large dynamical range of the sample, selection biases do not appear to affect our conclusions significantly. Our results, along with other samples in the literature, suggest that the locally measured black hole mass–host stellar mass relation is already in place atz ∼ 1. -
Abstract This work studies the relationship between accretion-disk size and quasar properties, using a sample of 95 quasars from the Sloan Digital Sky Survey Reverberation Mapping Project with measured lags between the
g andi photometric bands. Our sample includes disk lags that are both longer and shorter than predicted by the Shakura and Sunyaev model, requiring explanations that satisfy both cases. Although our quasars each have one lag measurement, we explore the wavelength-dependent effects of diffuse broad-line region (BLR) contamination through our sample’s broad redshift range, 0.1 <z < 1.2. We do not find significant evidence of variable diffuse Feii and Balmer nebular emission in the rms spectra, nor from Anderson–Darling tests of quasars in redshift ranges with and without diffuse nebular emission falling in the observed-frame filters. Contrary to previous work, we do not detect a significant correlation between the measured continuum and BLR lags in our luminous quasar sample, similarly suggesting that our continuum lags are not dominated by diffuse nebular emission. Similar to other studies, we find that quasars with larger-than-expected continuum lags have lower 3000 Å luminosities, and we additionally find longer continuum lags with lower X-ray luminosities and black hole masses. Our lack of evidence for diffuse BLR contribution to the lags indicates that the anticorrelation between continuum lag and luminosity is not likely to be due to the Baldwin effect. Instead, these anticorrelations favor models in which the continuum lag increases in lower-luminosity active galactic nuclei, including scenarios featuring magnetic coupling between the accretion disk and X-ray corona, and/or ripples or rims in the disk. -
Abstract Spectroscopic studies of extreme-ionization galaxies (EIGs) are critical to our understanding of exotic systems throughout cosmic time. These EIGs exhibit spectral features requiring >54.42 eV photons: the energy needed to ionize helium into He2+fully and emit He
ii recombination lines. Spectroscopic studies of EIGs can probe exotic stellar populations or accretion onto intermediate-mass black holes (∼102–105M ⊙), which are the possibly key contributors to the reionization of the Universe. To facilitate the use of EIGs as probes of high-ionization systems, we focus on ratios constructed from several rest-frame UV/optical emission lines: [Oiii ]λ 5008, Hβ , [Neiii ]λ 3870, [Oii ]λ λ 3727, 3729, and [Nev ]λ 3427. These lines probe the relative intensity at energies of 35.12, 13.62, 40.96, 13.62, and 97.12 eV, respectively, covering a wider range of ionization than traced by other common rest-frame UV/optical techniques. We use the ratios of these lines ([Nev ]/[Neiii ] ≡ Ne53, [Oiii ]/Hβ , and [Neiii ]/[Oii ]), which are nearby in wavelength, mitigating the effects of dust attenuation and uncertainties in flux calibration. We make predictions from photoionization models constructed fromCloudy that use a broad range of stellar populations and black hole accretion models to explore the sensitivity of these line ratios to changes in the ionizing spectrum. We compare our models to observations from the Hubble Space Telescope and JWST of galaxies with strong high-ionization emission lines atz ∼ 0,z ∼ 2, and 5 <z < 8.5. We show that the Ne53 ratio can separate galaxies with ionization from “normal” stellar populations from those with active galactic nuclei and even “exotic” Population III models. We introduce new selection methods to identify galaxies with photoionization driven by Population III stars or intermediate-mass black hole accretion disks that could be identified in upcoming high-redshift spectroscopic surveys. -
Abstract We analyze a sample of 25 [Ne
v ] (λ 3426) emission-line galaxies at 1.4 <z < 2.3 using Hubble Space Telescope/Wide Field Camera 3 G102 and G141 grism observations from the CANDELS Lyα Emission at Reionization (CLEAR) survey. [Nev ] emission probes extremely energetic photoionization (creation potential of 97.11 eV) and is often attributed to energetic radiation from active galactic nuclei (AGNs), shocks from supernovae, or an otherwise very hard ionizing spectrum from the stellar continuum. In this work, we use [Nev ] in conjunction with other rest-frame UV/optical emission lines ([Oii ]λ λ 3726, 3729, [Neiii ]λ 3869, Hβ , [Oiii ]λ λ 4959, 5007, Hα +[Nii ]λ λ 6548, 6583, [Sii ]λ λ 6716, 6731), deep (2–7 Ms) X-ray observations (from Chandra), and mid-infrared imaging (from Spitzer) to study the origin of this emission and to place constraints on the nature of the ionizing engine. The majority of the [Nev ]-detected galaxies have properties consistent with ionization from AGNs. However, for our [Nev ]-selected sample, the X-ray luminosities are consistent with local (z ≲ 0.1) X-ray-selected Seyferts, but the [Nev ] luminosities are more consistent with those fromz ∼ 1 X-ray-selected QSOs. The excess [Nev ] emission requires either reduced hard X-rays or a ∼0.1 keV excess. We discuss possible origins of the apparent [Nev ] excess, which could be related to the “soft (X-ray) excess” observed in some QSOs and Seyferts and/or be a consequence of a complex/anisotropic geometry for the narrow-line region, combined with absorption from a warm, relativistic wind ejected from the accretion disk. We also consider implications for future studies of extreme high-ionization systems in the epoch of reionization (z ≳ 6) with the James Webb Space Telescope. -
Abstract We investigate spatially resolved emission-line ratios in a sample of 219 galaxies (0.6 <
z < 1.3) detected using the G102 grism on the Hubble Space Telescope Wide Field Camera 3 taken as part of the CANDELS Lyα Emission at Reionization survey to measure ionization profiles and search for low-luminosity active galactic nuclei (AGN). We analyze [Oiii ] and Hβ emission-line maps, enabling us to spatially resolve the [Oiii ]/Hβ emission-line ratio across the galaxies in the sample. We compare the [Oiii ]/Hβ ratio in galaxy centers and outer annular regions to measure ionization differences and investigate the potential of sources with nuclear ionization to host AGN. We investigate some of the individual galaxies that are candidates to host strong nuclear ionization and find that they often have low stellar mass and are undetected in X-rays, as expected for low-luminosity AGN in low-mass galaxies. We do not find evidence for a significant population of off-nuclear AGN or other clumps of off-nuclear ionization. We model the observed distribution of [Oiii ]/Hβ spatial profiles and find that most galaxies are consistent with a small or zero difference between their nuclear and off-nuclear line ratios, but 6%–16% of galaxies in the sample are likely to host nuclear [Oiii ]/Hβ that is ∼0.5 dex higher than in their outer regions. This study is limited by large uncertainties in most of the measured [Oiii ]/Hβ spatial profiles; therefore, deeper data, e.g., from deeper HST/WFC3 programs or from JWST/NIRISS, are needed to more reliably measure the spatially resolved emission-line conditions of individual high-redshift galaxies. -
ABSTRACT Supermassive black hole binaries (SMBHBs) are a natural outcome of galaxy mergers and should form frequently in galactic nuclei. Sub-parsec binaries can be identified from their bright electromagnetic emission, e.g. Active Galactic Nuclei (AGNs) with Doppler shifted broad emission lines or AGN with periodic variability, as well as from the emission of strong gravitational radiation. The most massive binaries (with total mass >108M⊙) emit in the nanohertz band and are targeted by Pulsar Timing Arrays (PTAs). Here we examine the synergy between electromagnetic and gravitational wave signatures of SMBHBs. We connect both signals to the orbital dynamics of the binary and examine the common link between them, laying the foundation for joint multimessenger observations. We find that periodic variability arising from relativistic Doppler boost is the most promising electromagnetic signature to connect with GWs. We delineate the parameter space (binary total mass/chirp mass versus binary period/GW frequency) for which joint observations are feasible. Currently multimessenger detections are possible only for the most massive and nearby galaxies, limited by the sensitivity of PTAs. However, we demonstrate that as PTAs collect more data in the upcoming years, the overlapping parameter space is expected to expand significantly.
-
Abstract We explore the possibility and practical limitations of using a three-line approach to measure both the slope and normalization of the dust attenuation law in individual galaxies. To do this, we focus on a sample of 11 galaxies with existing ground-based Balmer H
α and Hβ measurements from slit spectra, plus space-based grism constraints on Paschen-β . When accounting for observational uncertainties, we show that one galaxy has a well-constrained dust-law slope and normalization in the range expected from theoretical arguments; this galaxy therefore provides an example of what may be possible in the future. However, most of the galaxies are best fit by unusually steep or shallow slopes. We then explore whether additional astrophysical effects or observational biases could explain the elevated Paschen-β /Hα ratios driving these results. We find that galaxies with high Paschen-β /Hα ratios may be explained by slightly sub-unity covering fractions (>97%). Alternatively, differing slit losses for different lines can have a large impact on the results, emphasizing the importance of measuring all three lines with a consistent spectroscopic aperture. We conclude that, while the three-line approach to constraining the shape of the dust attenuation law in individual galaxies is promising, deep observations and a consistent observational strategy will be required to minimize observational biases and to disentangle the astrophysically interesting effect of differing covering fractions. The James Webb Space Telescope will provide more sensitive measurements of Balmer and Paschen lines for galaxies atz ≈ 0.3–2, enabling uniform constraints on the optical–infrared dust attenuation law and its intrinsic variation. -
Abstract We present an overview of the CANDELS Ly
α Emission At Reionization (CLEAR) survey. CLEAR is a 130 orbit program of the Hubble Space Telescope using the Wide Field Camera 3 (WFC3) IR G102 grism. CLEAR targets 12 pointings divided between the GOODS-N and GOODS-S fields of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). Combined with existing spectroscopic data from other programs, the full CLEAR data set includes spectroscopic imaging of these fields over 0.8–1.7μ m. In this paper, we describe the CLEAR survey, the survey strategy, the data acquisition, reduction, processing, and science products and catalogs released alongside this paper. The catalogs include emission line fluxes and redshifts derived from the combination of the photometry and grism spectroscopy for 6048 galaxies, primarily ranging from 0.2 ≲z ≲ 3. We also provide an overview of CLEAR’s science goals and results. In conjunction with this paper we provide links to electronic versions of the data products, including 1D+2D extracted spectra and emission line maps.