Search for: All records

Abstract Binary supermassive black holes (SMBHs) are consequences of galaxy mergers and dominate the low-frequency gravitational-wave background. Finding binary SMBHs in existing time-domain observations has proven difficult, as their periodic, electromagnetic signals can be confused with the natural variability of single quasars. In this work, we investigate the effects of host-galaxy contamination and survey design (cadence and duration) on the detectability of binary SMBHs with the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We simulate millions of LSST light curves of single and binary quasars, with a distribution of quasar and host-galaxy properties motivated by empirical observations and the anticipated LSST detection space. We then apply simple sinusoidal curve fits as a potential computationally inexpensive detection method. We find that host-galaxy contamination will increase false-positive rates and decrease binary parameter recovery rates. Lower-mass, lower-luminosity binary systems are most likely to be negatively affected by host-galaxy contamination. We also find that monitoring duration affects binary detection more than survey effective cadence for this detection method. As the light-curve duration increases, false-positive rates are suppressed and binary parameter recovery rates, especially for binary periods, are improved. Increasing the light-curve duration from 5 to 10 yr shows the most dramatic improvement for successful binary detection and false-positive rejection, with additional improvement from extending the light-curve duration to 20 yr. The observation duration increase is especially critical for recovering binary periods that are longer than a decade. 

Abstract Supermassive black holes (SMBHs) are found in the centers of massive galaxies, and galaxy mergers should eventually lead to SMBH mergers. Quasar activity has long been associated with galaxy mergers, so here we investigate if supermassive black hole binaries (SMBHBs) are preferentially found in quasars. Our multimessenger investigation folds together a gravitational-wave background signal from NANOGrav, a sample of periodic active galactic nucleus candidates from the Catalina Real-Time Transient Survey, and a quasar mass function, to estimate an upper limit on the fraction of quasars which could host an SMBHB. We find at 95% confidence that quasars are at most 5 times as likely to host an SMBHB as a random galaxy. Pulsar timing arrays may therefore be more likely to find SMBHBs by prioritizing quasars over a random selection of galaxies in their searches. 

Abstract Multiyear observations from the Sloan Digital Sky Survey (SDSS) Reverberation Mapping ​(RM) project have significantly increased the number of quasars with reliable RM lag measurements. We statistically analyze target properties, light-curve characteristics, and survey design choices to identify factors crucial for successful and efficient RM surveys. Analyzing 172 high-confidence (“gold”) lag measurements from SDSS-RM for the Hβ, Mgii, and Civemission lines, we find that the Durbin–Watson statistic (a statistical test for residual correlation) is the most significant predictor of light curves suitable for lag detection. The variability signal-to-noise ratio and emission-line placement on the detector also correlate with successful lag measurements. We further investigate the impact of the observing cadence on the survey design by analyzing the effect of reducing observations in the first year of SDSS-RM. Our results demonstrate that a modest reduction in the observing cadence to ∼1.5 weeks between observations can retain approximately 90% of the lag measurements compared to twice-weekly observations in the initial year. Provided similar and uniform sampling in subsequent years, this adjustment has a minimal effect on the overall recovery of lags across all emission lines. These results provide valuable inputs for optimizing future RM surveys. 

We use JWST Near-Infrared Spectrograph observations from the Cosmic Evolution Early Release survey, GLASS-JWST ERS (GLASS), and JWST Advanced Deep Extragalactic Survey to measure rest-frame optical emission-line ratios of 89 galaxies atz > 4. The stacked spectra of galaxies with and without a broad-line feature reveal a difference in the [Oiii]λ4364 and Hγratios. This motivated our investigation of the [Oiii]λ4364/Hγversus [Neiii]/[Oii] diagram. We define two active galactic nucleus (AGN)/star formation (SF) classification lines based on 21,048 Sloan Digital Sky Survey galaxies atz ∼ 0. After applying a redshift correction to the AGN/SF lines, we find 69.2% of broad-line active galactic nuclei (BLAGN) continue to land in the AGN region of the diagnostic, largely due to the [Neiii]/[Oii] ratio. However, 33.0% of non-BLAGN land is in the AGN region as well. The [Oiii]λ4364/Hγversus [Neiii]/[Oii] diagram does not robustly separate BLAGN from non-broad-line galaxies atz> 4. This could be due to star-forming galaxies having harder ionization, or these galaxies contain a narrow line AGN, which are not accounted for. We further inspected galaxies without broad emission lines in each region of [Oiii]λ4364/Hγversus [Neiii]/[Oii] diagram and found that they have slightly stronger Ciii]λ1908 fluxes and equivalent width when landing in the BLAGN region. However, the cause of this higher ionization is unclear and may be revealed by observing UV lines. 

Abstract We present dynamical modeling of the broad-line region (BLR) of the highly variable active galactic nucleus (AGN) SDSS J141041.25+531849.0 (z= 0.359) using photometric and spectroscopic monitoring data from the Sloan Digital Sky Survey (SDSS) Reverberation Mapping project and the current fifth-generation SDSS Black Hole Mapper program, spanning from early 2013 to early 2023. We model the geometry and kinematics of the BLR in the Hβ, Hα, and Mgiiemission lines for three different time periods to measure the potential change of structure within the BLR across time and line species. We find a moderately face-on $(i_{\mathrm{full}-\mathrm{state}}=29.{68}_{-3.62}^{+4.74}\mathrm{deg})$thick-disk $({\theta }_{\mathrm{opn},\mathrm{full}-\mathrm{state}}=42.{04}_{-3.96}^{+4.32}\mathrm{deg})$geometry for most BLRs, with a joint estimate for the mass of the supermassive black hole for each of three time periods, yielding ${\log }_{10}(M_{\mathrm{BH}}/M_{\odot })=8.10_{-0.03}^{+0.03}$when using the full data set. The inferred individual virial factorf∼ 1.6 is moderately smaller than the average factor for a local sample of dynamically modeled AGNs. There is strong evidence for nonvirial motion, with over 70% of clouds on inflowing/outflowing orbits. We analyze the change in model parameters across emission lines, finding the radii of BLRs for the emission lines are consistent with the following relative sizesR_Hβ ≲ R_MgII ≲ R_Hα. Comparing results across time, we findR_low-state ≲ R_high-state, with the change in BLR size for Hβbeing more significant than for the other two lines. The data also reveal complex, time-evolving, and potentially transient dynamics of the BLR gas over a decade-long timescale, encouraging for future dynamical modeling of fine-scale BLR kinematics. 

JWST spectroscopy has discovered a population ofz ≳ 3.5 galaxies with broad Balmer emission lines and narrow forbidden lines that are consistent with hosting active galactic nuclei (AGN). Many of these systems, now known as “little red dots,” are compact and have unique colors that are very red in the optical/near-infrared and blue in the ultraviolet. The relative contribution of galaxy starlight and AGN to these systems remains uncertain, especially for the galaxies with unusual blue+red spectral energy distributions. In this work, we use Balmer decrements to measure the independent dust attenuation of the broad and narrow emission-line components of a sample of 29 broad-line AGN identified from three public JWST spectroscopy surveys: CEERS, JADES, and RUBIES. Stacking the narrow components from the spectra of 25 sources with broad Hαand no broad Hβresults in a median narrow Hα/Hβ= $2.47_{-0.05}^{+0.05}$(consistent withA_v = 0) and broad Hα/Hβ>8.85 (A_v > 3.63). The narrow and broad Balmer decrements imply little to no attenuation of the narrow emission lines, which are consistent with being powered by star formation and located on larger physical scales. Meanwhile, the lower limit in the broad Hα/Hβdecrement, with broad Hβundetected in the stacked spectrum of 25 broad HαAGN, implies significant dust attenuation of the broad-line emitting region that is presumably associated with the central AGN. Our results indicate that these systems, on average, are consistent with heavily dust-attenuated AGN powering the red parts of their SED, while their blue UV emission is powered by unattenuated star formation in the host galaxy. 

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. 

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 massM_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 (∼10⁹–10¹¹M_⊙) and black hole mass (∼10⁷–10⁹M_⊙) 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 $\log \left(M_{\mathrm{BH}}/M_{\odot }\right)={7.01}_{-0.33}^{+0.23}+{1.74}_{-0.64}^{+0.64}\log \left(M_{*,\mathrm{host}}/{10}^{10}{M}_{\odot }\right)$, with an intrinsic scatter of ${0.47}_{-0.17}^{+0.24}$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. 

We present a velocity-resolved reverberation mapping analysis of the hypervariable quasar RM160 (SDSS J141041.25+531849.0) atz= 0.359 with 153 spectroscopic epochs of data representing a 10 yr baseline (2013–2023). We split the baseline into two regimes based on the 3× flux increase in the light curve: a “low state” phase during the years 2013–2019 and a “high state” phase during the years 2022–2023. The velocity-resolved lag profiles (VRLPs) indicate that gas with different kinematics dominates the line emission in different states. The HβVRLP begins with a signature of inflow onto the broad-line region (BLR) in the low state, while in the high state it is flatter with less signature of inflow. The HαVRLP begins consistent with a virialized BLR in the low state, while in the high state shows a signature of inflow. The differences in the kinematics between the Balmer lines and between the low state and the high state suggests complex BLR dynamics. We find that the BLR radius and velocity (both FWHM andσ) do not obey a constant virial product throughout the monitoring period. We find that the BLR lags and continuum luminosity are correlated, consistent with rapid response of the BLR gas to the illuminating continuum. The BLR kinematic profile changes in unpredictable ways that are not related to continuum changes and reverberation lag. Our observations indicate that nonvirial kinematics can significantly contribute to observed line profiles, suggesting caution for black hole mass estimation in luminous and highly varying quasars like RM160.