Search

The detection of a stochastic gravitational wave background by pulsar timing arrays suggests the presence of a supermassive black hole binary population. Although the observed spectrum generally aligns with predictions from orbital evolution driven by gravitational wave emission in circular orbits, there is a discernible preference for a turnover at the lowest observed frequencies. This turnover could indicate a significant hardening phase, transitioning from early environmental influences to later stages predominantly influenced by gravitational wave emission. In the vicinity of these binaries, the ejection of stars or dark matter particles through gravitational three-body slingshots efficiently extracts orbital energy, leading to a low-frequency turnover in the spectrum. By analyzing the NANOGrav 15-year data, we assess how the gravitational wave spectrum depends on the initial inner galactic profile prior to disruption by binary ejections, accounting for a range of initial binary eccentricities. Our findings suggest a parsec-scale galactic center density around $10^6\,M_\odot/\textrm{pc}^3$ across most of the parameter space, offering insights into the environmental effects on black hole evolution and combined matter density near galaxy centers., Comment: 15 pages, 5 figures

We present a velocity-resolved reverberation mapping analysis of the hypervariable quasar RM160 (SDSS J141041.25+531849.0) at z = 0.359 with 153 spectroscopic epochs of data representing a ten-year baseline (2013-2023). We split the baseline into two regimes based on the 3x 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 (VRLP) indicate that gas with different kinematics dominates the line emission in different states. The H\b{eta} VRLP begins with a signature of inflow onto the BLR in the 'low state', while in the 'high state' it is flatter with less signature of inflow. The H{\alpha} 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 {\sigma}) do not obey a constant virial product throughout the monitoring period. We find that 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 non-virial kinematics can significantly contribute to observed line profiles, suggesting caution for black-hole mass estimation in luminous and highly varying quasars like RM160., Comment: 23 pages, 17 figures

The physical processes that produce X-ray Quasi-Periodic Eruptions (QPEs) recently discovered from the nuclei of several low-redshift galaxies are mysterious. Several pieces of observational evidence strongly suggest a link between QPEs and Tidal Disruption Events (TDE). Previous studies also reveal that the morphologies of TDE host galaxies are highly concentrated, with high Sersic indicies, bulge-to-total light (B/T) ratios, and stellar surface mass densities relative to the broader galaxy population. We use these distinctive properties to test the link between QPEs and TDEs, by comparing these parameters of QPE host galaxies to TDE host galaxies. We employ archival Legacy Survey images of a sample of 9 QPE host galaxies and a sample of 13 TDE host galaxies, and model their surface brightness profiles. We show that QPE host galaxies have high Sersic indices of ~3, high B/T ratios of ~0.5, and high surface mass densities of ~10^10 Msun kpc^-2. These properties are similar to TDE host galaxies, but are in strong contrast to a mass- and redshift-matched control sample of galaxies. We also find tentative evidence that the central black holes in both QPE and TDE host galaxies are undermassive relative to their stellar mass. The morphological similarities between QPE and TDE host galaxies at the population level add to the mounting evidence of a physical link between these phenomena, and favor QPE models that also invoke TDEs., Comment: 12 pages, 4 figures, submitted to ApJ

The NANOGrav 15-year data provides compelling evidence for a stochastic gravitational-wave (GW) background at nanohertz frequencies. The simplest model-independent approach to characterizing the frequency spectrum of this signal consists in a simple power-law fit involving two parameters: an amplitude A and a spectral index \gamma. In this paper, we consider the next logical step beyond this minimal spectral model, allowing for a running (i.e., logarithmic frequency dependence) of the spectral index, \gamma_run(f) = \gamma + \beta \ln(f/f_ref). We fit this running-power-law (RPL) model to the NANOGrav 15-year data and perform a Bayesian model comparison with the minimal constant-power-law (CPL) model, which results in a 95% credible interval for the parameter \beta consistent with no running, \beta \in [-0.80,2.96], and an inconclusive Bayes factor, B(RPL vs. CPL) = 0.69 +- 0.01. We thus conclude that, at present, the minimal CPL model still suffices to adequately describe the NANOGrav signal; however, future data sets may well lead to a measurement of nonzero \beta. Finally, we interpret the RPL model as a description of primordial GWs generated during cosmic inflation, which allows us to combine our results with upper limits from big-bang nucleosynthesis, the cosmic microwave background, and LIGO-Virgo-KAGRA., Comment: 17 pages, 4 figures, 2 tables

Pulsar-timing-array experiments have reported evidence for a stochastic background of nanohertz gravitational waves consistent with the signal expected from a population of supermassive--black-hole binaries. Those analyses assume power-law spectra for intrinsic pulsar noise and for the background, as well as a Hellings--Downs cross-correlation pattern among the gravitational-wave--induced residuals across pulsars. These assumptions are idealizations that may not be realized in actuality. We test them in the NANOGrav 15 yr data set using Bayesian posterior predictive checks: after fitting our fiducial model to real data, we generate a population of simulated data-set replications, and use them to assess whether the optimal-statistic significance, inter-pulsar correlations, and spectral coefficients assume extreme values for the real data when compared to the replications. We confirm that the NANOGrav 15 yr data set is consistent with power-law and Hellings--Downs assumptions. We also evaluate the evidence for the stochastic background using posterior-predictive versions of the frequentist optimal statistic and of Bayesian model comparison, and find comparable significance (3.2\ $\sigma$ and 3\ $\sigma$ respectively) to what was previously reported for the standard statistics. We conclude with novel visualizations of the reconstructed gravitational waveforms that enter the residuals for each pulsar. Our analysis strengthens confidence in the identification and characterization of the gravitational-wave background as reported by NANOGrav., Comment: 20 pages, 14 Figures

The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 \sigma - 1.9 \sigma$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 \sigma - 2.1 \sigma$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26$~nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe., Comment: 10 pages, 8 figures, 1 appendix, submitted to ApJ

Periodic signatures in time-domain observations of quasars have been used to search for binary supermassive black holes. 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 supermassive black holes 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. 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: 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 supermassive black holes. However, pure electromagnetic detections of binaries using photometric periodicity with amplitude greater than 0.1 magnitude 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., Comment: 21 pages, 14 figures, 3 tables

The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations., Comment: Accepted by Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org. Edited to fix two equation typos (Eq.13 & 21), and minor text typos

Quasars at $z~{\gtrsim}~1$ most often have redshifts measured from rest-frame ultraviolet emission lines. One of the most common such lines, C IV ${\lambda}1549$, shows blueshifts up to ${\approx}~5000~\rm{km~s^{-1}}$, and in rare cases even higher. This blueshifting results in highly uncertain redshifts when compared to redshift determinations from rest-frame optical emission lines, e.g., from the narrow [O III] ${\lambda}5007$ feature. We present spectroscopic measurements for 260 sources at $1.55~{\lesssim}~z~{\lesssim}~3.50$ having $-28.0~{\lesssim}~M_i~{\lesssim}~-30.0$ mag from the Gemini Near Infrared Spectrograph - Distant Quasar Survey (GNIRS-DQS) catalog, augmenting the previous iteration which contained 226 of the 260 sources whose measurements are improved upon in this work. We obtain reliable systemic redshifts based on [O III] ${\lambda}5007$ for a subset of 121 sources which we use to calibrate prescriptions for correcting UV-based redshifts. These prescriptions are based on a regression analysis involving C IV full-width-at-half-maximum intensity and equivalent width, along with the UV continuum luminosity at a rest-frame wavelength of 1350 A. Applying these corrections can improve the accuracy and the precision in the C IV-based redshift by up to ${\sim}~850~\rm{km~s^{-1}}$ and ${\sim}~150~\rm{km~s^{-1}}$, respectively, which correspond to ${\sim}~8.5$ Mpc and ${\sim}~1.5$ Mpc in comoving distance at $z~=~2.5$. Our prescriptions also improve the accuracy of the best available multi-feature redshift determination algorithm by ${\sim}~100~\rm{km~s^{-1}}$, indicating that the spectroscopic properties of the C IV emission line can provide robust redshift estimates for high-redshift quasars. We discuss the prospects of our prescriptions for cosmological and quasar studies utilizing upcoming large spectroscopic surveys., Comment: 20 pages (AASTeX 6.3.1), 8 figures, accepted for publication in ApJ

Weak emission-line quasars (WLQs) are a subset of Type 1 quasars that exhibit extremely weak Ly$\alpha +$N V $\lambda$1240 and/or C IV $\lambda$1549 emission lines. We investigate the relationship between emission-line properties and accretion rate for a sample of 230 `ordinary' Type 1 quasars and 18 WLQs at $z < 0.5$ and $1.5 < z < 3.5$ that have rest-frame ultraviolet and optical spectral measurements. We apply a correction to the H$\beta$-based black-hole mass ($M_{\rm BH}$) estimates of these quasars using the strength of the optical Fe II emission. We confirm previous findings that WLQs' $M_{\rm BH}$ values are overestimated by up to an order of magnitude using the traditional broad emission-line region size-luminosity relation. With this $M_{\rm BH}$ correction, we find a significant correlation between H$\beta$-based Eddington luminosity ratios and a combination of the rest-frame C IV equivalent width and C IV blueshift with respect to the systemic redshift. This correlation holds for both ordinary quasars and WLQs, which suggests that the two-dimensional C IV parameter space can serve as an indicator of accretion rate in all Type 1 quasars across a wide range of spectral properties., Comment: 17 pages (AASTeX 6.3.1), 5 figures, accepted for publication in ApJ

We present a high-cadence multi-epoch analysis of dramatic variability of three broad emission lines (MgII, H$\beta$, and H$\alpha$) in the spectra of the luminous quasar ($\lambda L_{\lambda}$(5100\r{A}) = $4.7 \times 10^{44}$ erg s$^{-1}$) SDSS J141041.25+531849.0 at $z = 0.359$ with 127 spectroscopic epochs over 9 years of monitoring (2013-2022). We observe anti-correlations between the broad emission-line widths and flux in all three emission lines, indicating that all three broad emission lines "breathe" in response to stochastic continuum variations. We also observe dramatic radial velocity shifts in all three broad emission lines, ranging from $\Delta{v}$ $\sim$400 km s$^{-1}$ to $\sim$800 km s$^{-1}$, that vary over the course of the monitoring period. Our preferred explanation for the broad-line variability is complex kinematics in the broad-line region gas. We suggest a model for the broad-line variability that includes a combination of gas inflow with a radial gradient, an azimuthal asymmetry (e.g., a hot spot), superimposed on the stochastic flux-driven changes to the optimal emission region ("line breathing"). Similar instances of line-profile variability due to complex gas kinematics around quasars are likely to represent an important source of false positives in radial velocity searches for binary black holes, which typically lack the kind of high-cadence data we analyze here. The long-duration, wide-field, and many-epoch spectroscopic monitoring of SDSS-V BHM-RM provides an excellent opportunity for identifying and characterizing broad emission-line variability, and the inferred nature of the inner gas environment, of luminous quasars.

Feedback from active galactic nuclei (AGN) has proven to be a critical ingredient in the current picture of galaxy assembly and growth. However, observational constraints on AGN-driven outflows face technical challenges and as a result, the cold molecular gas outflow properties of type-1 AGN are not well known. We present new IRAM NOrthern Extended Milimeter Array (NOEMA) observations of CO (1-0) in F07599+6508, Z11598$-$0112, F13342+3932, and PG1440+356, all nearby type-1 AGN and ultraluminous infrared galaxies (ULIRGs). We achieve spatial resolution of 1-3 arcsec corresponding to physical scales of 2--8~kpc and spectral resolution of 15-60 km s$^{-1}$ which enables updated CO (1-0) redshifts and a detailed morphological view of the cold molecular gas in these sources. The CO (1-0) luminosities, $L_{CO}^{\prime}$, are in the range 2-12$\times 10^9$ K km s$^{-1}$ pc$^{2}$ and inferred molecular gas masses, $M$(H$_2$), are in the range 2-9$\times 10^9$ M$_{\odot}$. The velocity fields and gas distributions do not unambiguously identify any of these sources as having outflows. However, Z11598$-$0112 has signs of infalling material and after the subtraction of a rotating disk model PG 1440+356 shows complex kinematics in the residuals that may indicate an outflow or warped disk., Comment: 20 pages, 10 figures, accepted for publication in MNRAS

Recent claimed detections of tidal disruption events (TDEs) in multi-wavelength data have opened potential new windows into the evolution and properties of otherwise dormant supermassive black holes (SMBHs) in the centres of galaxies. At present, there are several dozen TDE candidates, which share some properties and differ in others. The range in properties is broad enough to overlap other transient types, such as active galactic nuclei (AGN) and supernovae (SNe), which can make TDE classification ambiguous. A further complication is that "TDE signatures" have not been uniformly observed to similar sensitivities or even targeted across all candidates. This chapter reviews those events that are unusual relative to other TDEs, including the possibility of TDEs in pre-existing AGN, and summarises those characteristics thought to best distinguish TDEs from continuously accreting AGN, strongly flaring AGN, SNe, and Gamma-Ray Bursts (GRBs), as well as other potential impostors like stellar collisions, "micro-TDEs," and circumbinary accretion flows. We conclude that multiple observables should be used to classify any one event as a TDE. We also consider the TDE candidate population as a whole, which, for certain host galaxy or SMBH characteristics, is distinguishable statistically from non-TDEs, suggesting that at least some TDE candidates do in fact arise from SMBH-disrupted stars., Comment: 57 pages, 17 figures. Accepted for publication in Springer Space Science Reviews. Chapter in ISSI review "The Tidal Disruption of Stars by Massive Black Holes" vol. 79. Corrected several typos from previous submission

"Changing-look" quasars are a new class of highly variable active galactic nuclei that have changed their spectral type over surprisingly short timescales of just a few years. The origin of this phenomenon is debated, but is likely to reflect some change in the accretion flow. To investigate the disk-corona systems in these objects, we measure optical/UV-X-ray spectral indices ($\alpha_{\rm OX}$) and Eddington ratios ($\lambda_{\rm Edd}$) of ten previously-discovered changing-look quasars at two or more epochs. By comparing these data with simulated results based on the behavior of X-ray binaries, we find possible similarities in spectral indices below 1% Eddington ratio. We further investigate the Eddington ratios of changing-look quasars before and after their spectral type changes, and find that changing-look quasars cross the 1% Eddington ratio boundary when their broad emission lines disappear/emerge. This is consistent with the disk-wind model as the origin of broad emission lines., Comment: 31 pages, 9 figures, accepted for publication in ApJ

We present results from a multi-wavelength analysis searching for multiple AGN systems in nearby (z<0.077) triple galaxy mergers. Combining archival Chandra, SDSS, WISE, and VLA observations, we quantify the rate of nearby triple AGN, as well as investigate possible connections between SMBH accretion and merger environments. Analyzing the multi-wavelength observations of 7 triple galaxy mergers, we find that 1 triple merger has a single AGN (NGC 3341); we discover, for the first time, 4 likely dual AGN (SDSS J1027+1749, SDSS J1631+2352, SDSS J1708+2153, and SDSS J2356-1016); we confirm one triple AGN system, SDSS J0849+1114; and 1 triple merger in our sample remains ambiguous (SDSS J0858+1822). Analyzing the WISE data, we find a trend of increasing N_H (associated with the primary AGN) as a function of increasing W1-W2 color, reflecting that the motions of gas and dust are coupled in merging environments, where large amount of both can be funneled into the active central region during mergers. Additionally, we find that the one triple AGN system in our sample has the highest levels of N_H and W1-W2 color, while the dual AGN candidates all have lower levels; these results are consistent with theoretical merger simulations that suggest higher levels of nuclear gas are more likely to activate AGN in mergers., Comment: 19 pages, 3 figures, accepted to ApJ

We present results from our X-ray analysis of the first systematic search for triple AGN in nearby (z<0.077) triple galaxy mergers. We analyze archival Chandra observations of 7 triple galaxy mergers with BAYMAX (Bayesian Analysis of Multiple AGN in X-rays), fitting each observation with single, dual, and triple X-ray point source models. In doing so, we conclude that 1 triple merger has one X-ray point source (SDSS J0858+1822, although it's unlikely to be an AGN); 5 triple mergers are likely composed of two X-ray point sources (NGC 3341, SDSS J1027+1749, SDSS J1631+2352, SDSS J1708+2153, and SDSS J2356$-$1016); and one system is composed of three X-ray point sources (SDSS J0849+1114). By fitting the individual X-ray spectra of each point source, we analyze the 2-7 keV luminosities as well as the levels of obscuration associated with each potential AGN. We find that 4/5 dual X-ray point source systems have primary and secondary point sources with bright X-ray luminosities (L_2-7 kev >10^40 erg s^-1), possibly associated with 4 new undetected dual AGN. The dual and triple point source systems are found to have physical separations between 3-9 kpc and flux ratios between 2x10^-3 - 0.84. A multi-wavelength analysis to determine the origin of the X-ray point sources discovered in this work is presented in our companion paper (Foord et al. 2020c)., Comment: 27 pages, 8 figures, accepted to ApJ

We present spectroscopic measurements for 226 sources from the Gemini Near Infrared Spectrograph - Distant Quasar Survey (GNIRS-DQS). Being the largest uniform, homogeneous survey of its kind, it represents a flux-limited sample ($m_{i}$ ${\lesssim}$ 19.0 mag, $H$ ${\lesssim}$ 16.5 mag) of Sloan Digital Sky Survey (SDSS) quasars at 1.5 ${\lesssim}$ $z$ ${\lesssim}$ 3.5 with a monochromatic luminosity (${\lambda}L_{\lambda}$) at 5100 ${\unicode{xC5}}$ in the range of $10^{44} - 10^{46}$ erg $\rm{s}^{-1}$. A combination of the GNIRS and SDSS spectra covers principal quasar diagnostic features, chiefly the C IV ${\lambda}$1549, Mg II ${\lambda}{\lambda}$2798, 2803, H${\beta}$ ${\lambda}$4861, and [O III] ${\lambda}{\lambda}$4959, 5007 emission lines, in each source. The spectral inventory will be utilized primarily to develop prescriptions for obtaining more accurate and precise redshifts, black hole masses, and accretion rates for all quasars. Additionally, the measurements will facilitate an understanding of the dependence of rest-frame ultraviolet-optical spectral properties of quasars on redshift, luminosity, and Eddington ratio, and test whether the physical properties of the quasar central engine evolve over cosmic time., Comment: 20 pages, 7 figures, accepted for publication in ApJS; catalog data can be obtained at http://physics.uwyo.edu/agn/GNIRS-DQS/spectra.html

Radio core dominance measurements, an indicator of jet orientation, sometimes rely on core flux density measurements from large-area surveys like Faint Images of the Radio Sky at Twenty cm (FIRST) that have an angular resolution of only 5''. Such low-resolution surveys often fail to resolve cores from the extended emission resulting in an erroneous measurement. We focus on investigating this resolution effect for a sample of 119 radio-loud quasars. We obtained continuum observations from NSF's Karl G. Jansky Very Large Array (VLA) at 10 GHz in A-array with a 0.2'' resolution. Our measurements show that at FIRST spatial resolution, core flux measurements are indeed systematically high even after considering the core-variability. For a handful of quasars, 10 GHz images reveal extended features, whereas the FIRST image shows a point source. We found that the resolution effect is more prominent for quasars with smaller angular sizes. We further computed two radio core dominance parameters R & R5100 for use in statistical orientation investigations with this sample. We also present the spectral energy distributions between 74 MHz and 1.4 GHz, which we used to measure the spectral index of the extended emission of these quasars. Our results empirically confirm that determination of radio core dominance requires high-spatial resolution data. We highlight the practical issues associated with the choice of frequency and resolution in the measurement of core and extended flux densities., Comment: Accepted by ApJ, 40 pages, 10 figures and 4 tables. All figures and machine readable tables are included in the source zip file (Download --> Other formats --> Source)

We estimate the merger timescale of spectroscopically-selected, subparsec supermassive black hole binary (SMBHB) candidates by comparing their expected contribution to the gravitational wave background (GWB) with the sensitivity of current pulsar timing array (PTA) experiments and in particular, with the latest upper limit placed by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). We find that the average timescale to coalescence of such SMBHBs is $\langle t_{\rm evol} \rangle > 6\times 10^4\,$yr, assuming that their orbital evolution in the PTA frequency band is driven by emission of gravitational waves. If some fraction of SMBHBs do not reside in spectroscopically detected active galaxies, and their incidence in active and inactive galaxies is similar, then the merger timescale could be $\sim 10$ times longer, $\langle t_{\rm evol} \rangle > 6\times 10^5\,$yr. These limits are consistent with the range of timescales predicted by theoretical models and imply that all the SMBHB candidates in our spectroscopic sample could be binaries without violating the observational constraints on the GWB. This result illustrates the power of the multi-messenger approach, facilitated by the PTAs, in providing an independent statistical test of the nature of SMBHB candidates discovered in electromagnetic searches., Comment: Published in ApJL (9 pages, 3 figures)

Velocity offsets in the broad Balmer lines of quasars and their temporal variations serve as indirect evidence for bound supermassive black hole binaries (SBHBs) at sub-parsec separations. In this work, we test the SBHB hypothesis for 14 quasars with double-peaked broad emission lines using their long-term (14--41 years) radial velocity curves. We improve on previous work by (a) using elliptical instead of circular orbits for the SBHBs, (b) adopting a statistical model for radial velocity jitter, (c) employing a Markov Chain Monte Carlo method to explore the orbital parameter space efficiently and build posterior distributions of physical parameters and (d) incorporating new observations. We determine empirically that jitter comprises approximately Gaussian distributed fluctuations about the smooth radial velocity curves that are larger than the measurement errors by factors of order a few. We initially treat jitter by enlarging the effective error bars and then verify this approach via a variety of Gaussian process models for it. We find lower mass limits for the hypothesized SBHBs in the range $10^8$--$10^{11}\;M_{\odot}$. For seven objects the SBHB scenario appears unlikely based on goodness-of-fit tests. For two additional objects the minimum SBHB masses are unreasonably large ($>10^{10}\;M_{\odot}$), strongly disfavoring the SBHB scenario. Using constraints on the orbital inclination angle (which requires some assumptions) makes the minimum masses of four more objects unreasonably large. We also cite physical and observational arguments against the SBHB hypothesis for nine objects. We conclude that the SBHB explanation is not the favoured explanation of double-peaked broad emission lines., Comment: 27 pages, accepted by MNRAS

Physical models of X-ray binary outbursts can aid in understanding the origin of 'changing-look' active galactic nuclei (AGN), if we can establish that these two black hole accretion phenomena are analogous. Previously, studies of the correlation between the UV-to-X-ray spectral index alpha_OX and Eddington ratio using single-epoch observations of changing-look AGN samples have revealed possible similarities to the spectral evolution of outbursting X-ray binaries. However, direct comparisons using multi-epoch UV/X-ray light curves of individual changing-look AGN undergoing dramatic changes in Eddington ratio have been scarce. Here, we use published Swift UV/X-ray light curves of two changing-look AGN (NGC 2617 and ZTF18aajupnt) to examine the evolution of their alpha_OX values during outburst. We show that the combination of these two changing-look AGN can trace out the predicted spectral evolution from X-ray binary outbursts, including the inversion in the evolution of alpha_OX as a function of Eddington ratio. We suggest that the spectral softening that is observed to occur below a critical Eddington ratio in both AGN and X-ray binaries is due to reprocessing of Comptonized X-ray emission by the accretion disk, based on the X-ray to UV reverberation lags previously observed in NGC 2617. Our results suggest that the physical processes causing the changing-look AGN phenomenon are similar to those in X-ray binary outbursts., Comment: 14 pages, 5 figures. ApJ submitted

We present a method for comparing the H$\beta$ emission-line profiles of observed supermassive black hole (SBHB) candidates and models of sub-parsec SBHBs in circumbinary disks. Using the approach based on principal component analysis we infer the values of the binary parameters for the spectroscopic SBHB candidates and evaluate the parameter degeneracies, representative of the uncertainties intrinsic to such measurements. We find that as a population, the SBHB candidates favor the average value of the semimajor axis corresponding to $\log(a/M) \approx 4.20\pm 0.42$ and comparable mass ratios, $q>0.5$. If the SBHB candidates considered are true binaries, this result would suggest that there is a physical process that allows initially unequal mass systems to evolve toward comparable mass ratios (e.g., accretion that occurs preferentially onto the smaller of the black holes) or point to some, yet unspecified, selection bias. Our method also indicates that the SBHB candidates equally favor configurations in which the mini-disks are coplanar or misaligned with the binary orbital plane. If confirmed for true SBHBs, this finding would indicate the presence of a physical mechanism that maintains misalignment of the mini-disks down to sub-parsec binary separations (e.g., precession driven by gravitational torques). The probability distributions of the SBHB parameters inferred for the observed SBHB candidates and our control group of AGNs are statistically indistinguishable, implying that this method can in principle be used to interpret the observed emission-line profiles once a sample of confirmed SBHBs is available but cannot be used as a conclusive test of binarity., Comment: Accepted to ApJ (31 pages, 19 figures), includes referee's comments

We present the first results from BAYMAX (Bayesian AnalYsis of Multiple AGN in X-rays), a tool that uses a Bayesian framework to quantitatively evaluate whether a given Chandra observation is more likely a single or dual point source. Although the most robust method of determining the presence of dual AGNs is to use X-ray observations, only sources that are widely separated relative to the instrument PSF are easy to identify. It becomes increasingly difficult to distinguish dual AGNs from single AGNs when the separation is on the order of Chandra's angular resolution (<1''). Using likelihood models for single and dual point sources, BAYMAX quantitatively evaluates the likelihood of an AGN for a given source. Specifically, we present results from BAYMAX analyzing the lowest-mass dual AGN candidate to date, SDSS J0914+0853, where archival Chandra data shows a possible secondary AGN ~0.3'' from the primary. Analyzing a new 50 ks Chandra observation, results from BAYMAX shows that SDSS J0914+0853 is most likely a single AGN with a Bayes factor of 13.5 in favor of a single point source model. Further, posterior distributions from the dual point source model are consistent with emission from a single AGN. We find the probability of SDSS J0914+0853 being a dual AGN system with a flux ratio f>0.3 and separation r>0.3'' to be very low. Overall, BAYMAX will be an important tool for correctly classifying candidate dual AGNs in the literature, and studying the dual AGN population where past spatial resolution limits have prevented systematic analyses., Comment: 16 pages, 5 figures, accepted to ApJ

Changing-look quasars are a newly-discovered class of luminous active galactic nuclei that undergo rapid ($\lesssim$10 year) transitions between Type 1 and Type 1.9/2, with an associated change in their continuum emission. We characterize the host galaxies of four faded changing-look quasars using broadband optical imaging. We use \textit{gri} images obtained with the Gemini Multi Object Spectrograph (GMOS) on Gemini North to characterize the surface brightness profiles of the quasar hosts and search for [O III] $\lambda4959,\lambda5007$ emission from spatially extended regions, or voorwerpjes, with the goal of using them to examine past luminosity history. Although we do not detect, voorwerpjes surrounding the four quasar host galaxies, we take advantage of the dim nuclear emission to characterize the colors and morphologies of the host galaxies. Three of the four galaxies show morphological evidence of merger activity or tidal features in their residuals. The three galaxies which are not highly distorted are fit with a single S\'ersic profile to characterize their overall surface brightness profiles. The single-S\'ersic fits give intermediate S\'ersic indices between the $n=1$ of disk galaxies and the $n=4$ of ellipticals. On a color-magnitude diagram, our changing-look quasar host galaxies reside in the blue cloud, with other AGN host galaxies and star-forming galaxies. On a color-S\'ersic index diagram the changing-look quasar hosts reside with other AGN hosts in the "green valley". Our analysis suggests that the hosts of changing-look quasars are predominantly disrupted or merging galaxies that resemble AGN hosts, rather than inactive galaxies., Comment: 20 pages, 5 figures

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts., Comment: Submitted to the Astro2020 decadal review. This is one of five core white papers written by members of the NANOGrav Collaboration, see also: J.Cordes et al., S.R.Taylor et al., X.Siemens et al., and E.Fonseca et al

Despite their factor of ~10^8 difference in black hole mass, several lines of evidence suggest possible similarities between black hole accretion flows in active galactic nuclei (AGN) and Galactic X-ray binaries. However, it is still unclear whether the geometry of the disk-corona system in X-ray binaries directly scale up to AGN, and whether this analogy still holds in different accretion states. We test this AGN/X-ray binary analogy, by comparing the observed correlations between the UV-to-X-ray spectral index (alpha_OX) and Eddington ratio in AGN to those predicted from observations of X-ray binary outbursts. This approach probes the geometry of their disk-corona systems as they transition between different accretion states. We use new Chandra X-ray and ground-based rest-UV observations of faded 'changing-look' quasars to extend this comparison to lower Eddington ratios of <10^-2, where observations of X-ray binaries predict a softening of alpha_OX in AGN. We find that the observed correlations between alpha_OX and Eddington ratio of AGN displays a remarkable similarity to accretion state transitions in prototypical X-ray binary outbursts, including an inversion of this correlation at a critical Eddington ratio of ~10^-2. Our results suggest that the structures of black hole accretion flows directly scale across a factor of ~10^8 in black hole mass and across different accretion states, enabling us to apply theoretical models of X-ray binaries to explain AGN phenomenology., Comment: 26 pages, 12 figures, published in ApJ, 883, 76

We present an improved semi-analytic model for calculation of the broad optical emission-line signatures from sub-parsec supermassive black hole binaries (SBHBs) in circumbinary disks. The second-generation model improves upon the treatment of radiative transfer by taking into account the effect of the radiation driven accretion disk wind on the properties of the emission-line profiles. Analysis of 42.5 million modeled emission-line profiles shows that correlations between the profile properties and SBHB parameters identified in the first-generation model are preserved, indicating that their diagnostic power is not diminished. The profile shapes are a more sensitive measure of the binary orbital separation and the degree of alignment of the black hole mini-disks, and are less sensitive to the SBHB mass ratio and orbital eccentricity. We also find that modeled profile shapes are more compatible with the observed sample of SBHB candidates than with our control sample of regular AGNs. Furthermore, if the observed sample of SBHBs is made up of genuine binaries, it must include compact systems with comparable masses, and misaligned mini-disks. We note that the model described in this paper can be used to interpret the observed emission-line profiles once a sample of confirmed SBHBs is available but cannot be used to prove that the observed SBHB candidates are true binaries., Comment: Accepted to ApJ (25 pages, 16 figures), includes referee's comments

Motivated by observational searches for sub-parsec supermassive black hole binaries (SBHBs) we develop a modular analytic model to determine the likelihood for detection of SBHBs by ongoing spectroscopic surveys. The model combines the parametrized rate of orbital evolution of SBHBs in circumbinary disks with the selection effects of spectroscopic surveys and returns a multivariate likelihood for SBHB detection. Based on this model we find that in order to evolve into the detection window of the spectroscopic searches from larger separations in less than a Hubble time, $10^8M_\odot$ SBHBs must, on average, experience angular momentum transport faster than that provided by a disk with accretion rate $0.06\,\dot{M}_E$. Spectroscopic searches with yearly cadence of observations are in principle sensitive to binaries with orbital separations $< {\rm few}\times 10^4\, r_g$ ($r_g = GM/c^2$ and $M$ is the binary mass), and for every one SBHB in this range there should be over 200 more gravitationally bound systems with similar properties, at larger separations. Furthermore, if spectra of all SBHBs in this separation range exhibit the AGN-like emission lines utilized by spectroscopic searches, the projection factors imply five undetected binaries for each observed $10^8M_\odot$ SBHB with mass ratio $0.3$ and orbital separation $10^4\,r_g$ (and more if some fraction of SBHBs is inactive). This model can be used to infer the most likely orbital parameters for observed SBHB candidates and to provide constraints on the rate of orbital evolution of SBHBs, if observed candidates are shown to be genuine binaries., Comment: Accepted to ApJ (16 pages, 11 figures). Revised version includes referee's comments

We present quasar bolometric corrections using the [O III] $\lambda5007$ narrow emission line luminosity based on the detailed spectral energy distributions of 53 bright quasars at low to moderate redshift ($0.0345

We have been spectroscopically monitoring 88 quasars selected to have broad H$\beta$ emission lines offset from their systemic redshift by thousands of km s$^{-1}$. By analogy with single-lined spectroscopic binary stars, we consider these quasars to be candidates for hosting supermassive black hole binaries (SBHBs). In this work we present new radial velocity measurements, typically 3-4 per object over a time period of up to 12 years in the observer's frame. In 29/88 of the SBHB candidates no variability of the shape of the broad H$\beta$ profile is observed, which allows us to make reliable measurements of radial velocity changes. Among these, we identify three objects that have displayed systematic and monotonic velocity changes by several hundred km s$^{-1}$ and are prime targets for further monitoring. Because the periods of the hypothetical binaries are expected to be long, we cannot hope to observe many orbital cycles during our lifetimes. Instead, we seek to evaluate the credentials of the SBHB candidates by attempting to rule out the SBHB hypothesis. In this spirit, we present a method for placing a lower limit on the period, and thus the mass, of the SBHBs under the assumption that the velocity changes we observe are due to orbital motion. Given the duration of our monitoring campaign and the uncertainties in the radial velocities, we were able to place a lower limit on the total mass in the range $4.7\times10^4-3.8\times10^8$ $M_{\scriptscriptstyle \odot}$, which does not yet allow us to rule out the SBHB hypothesis for any candidates., Comment: Accepted for publication in MNRAS, full version of Table 1 will be available in MNRAS online version of the paper

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV., Comment: Full information on DR13 available at http://www.sdss.org. Comments welcome to spokesperson@sdss.org. To be published in ApJS

Recent studies have found that obscured quasars cluster more strongly and are thus hosted by dark matter haloes of larger mass than their unobscured counterparts. These results pose a challenge for the simplest unification models, in which obscured objects are intrinsically the same as unobscured sources but seen through a dusty line of sight. There is general consensus that a structure like a "dusty torus" exists, meaning that this intrinsic similarity is likely the case for at least some subset of obscured quasars. However, the larger host halo masses of obscured quasars implies that there is a second obscured population that has an even higher clustering amplitude and typical halo mass. Here, we use simple assumptions about the host halo mass distributions of quasars, along with analytical methods and cosmological $N$-body simulations to isolate the signal from this population. We provide values for the bias and halo mass as a function of the fraction of the "non-torus obscured" population. Adopting a reasonable value for this fraction of $\sim$25% implies a non-torus obscured quasar bias that is much higher than the observed obscured quasar bias, because a large fraction of the obscured population shares the same clustering strength as the unobscured objects. For this non-torus obscured population, we derive a bias of $\sim$3, and typical halo masses of $\sim3\times10^{13}$ M$_{\odot}/h$ at $z=1$. These massive haloes are likely the descendants of high-mass unobscured quasars at high redshift, and will evolve into members of galaxy groups at $z=0$., Comment: 13 pages, 10 figures. Accepted for publication in MNRAS

The Time-Domain Spectroscopic Survey (TDSS) is an SDSS-IV eBOSS subproject primarily aimed at obtaining identification spectra of ~220,000 optically-variable objects systematically selected from SDSS/Pan-STARRS1 multi-epoch imaging. We present a preview of the science enabled by TDSS, based on TDSS spectra taken over ~320 deg^2 of sky as part of the SEQUELS survey in SDSS-III, which is in part a pilot survey for eBOSS in SDSS-IV. Using the 15,746 TDSS-selected single-epoch spectra of photometrically variable objects in SEQUELS, we determine the demographics of our variability-selected sample, and investigate the unique spectral characteristics inherent in samples selected by variability. We show that variability-based selection of quasars complements color-based selection by selecting additional redder quasars, and mitigates redshift biases to produce a smooth quasar redshift distribution over a wide range of redshifts. The resulting quasar sample contains systematically higher fractions of blazars and broad absorption line quasars than from color-selected samples. Similarly, we show that M-dwarfs in the TDSS-selected stellar sample have systematically higher chromospheric active fractions than the underlying M-dwarf population, based on their H-alpha emission. TDSS also contains a large number of RR Lyrae and eclipsing binary stars with main-sequence colors, including a few composite-spectrum binaries. Finally, our visual inspection of TDSS spectra uncovers a significant number of peculiar spectra, and we highlight a few cases of these interesting objects. With a factor of ~15 more spectra, the main TDSS survey in SDSS-IV will leverage the lessons learned from these early results for a variety of time-domain science applications., Comment: 17 pages, 14 figures, submitted to ApJ

We report the discovery of a new "changing-look" quasar, SDSS J101152.98+544206.4, through repeat spectroscopy from the Time Domain Spectroscopic Survey. This is an addition to a small but growing set of quasars whose blue continua and broad optical emission lines have been observed to decline by a large factor on a time scale of approximately a decade. The 5100 Angstrom monochromatic continuum luminosity of this quasar drops by a factor of > 9.8 in a rest-frame time interval of < 9.7 years, while the broad H-alpha luminosity drops by a factor of 55 in the same amount of time. The width of the broad H-alpha line increases in the dim state such that the black hole mass derived from the appropriate single-epoch scaling relation agrees between the two epochs within a factor of 3. The fluxes of the narrow emission lines do not appear to change between epochs. The light curve obtained by the Catalina Sky Survey suggests that the transition occurs within a rest-frame time interval of approximately 500 days. We examine three possible mechanisms for this transition suggested in the recent literature. An abrupt change in the reddening towards the central engine is disfavored by the substantial difference between the timescale to obscure the central engine and the observed timescale of the transition. A decaying tidal disruption flare is consistent with the decay rate of the light curve but not with the prolonged bright state preceding the decay, nor can this scenario provide the power required by the luminosities of the emission lines. An abrupt drop in the accretion rate onto the supermassive black hole appears to be the most plausible explanation for the rapid dimming., Comment: Submitted to MNRAS

The uncertain origin of the recently-discovered `changing-looking' quasar phenomenon -- in which a luminous quasar dims significantly to a quiescent state in repeat spectroscopy over ~10 year timescales -- may present unexpected challenges to our understanding of quasar accretion. To better understand this phenomenon, we take a first step to building a sample of changing-look quasars with a systematic but simple archival search for these objects in the Sloan Digital Sky Survey Data Release 12. By leveraging the >10 year baselines for objects with repeat spectroscopy, we uncover two new changing-look quasars, and a third discovered previously. Decomposition of the multi-epoch spectra and analysis of the broad emission lines suggest that the quasar accretion disk emission dims due to rapidly decreasing accretion rates (by factors of >2.5), while disfavoring changes in intrinsic dust extinction for the two objects where these analyses are possible. Broad emission line energetics also support intrinsic dimming of quasar emission as the origin for this phenomenon rather than transient tidal disruption events or supernovae. Although our search criteria included quasars at all redshifts and transitions from either quasar-like to galaxy-like states or the reverse, all of the clear cases of changing-look quasars discovered were at relatively low-redshift (z ~ 0.2 - 0.3) and only exhibit quasar-like to galaxy-like transitions., Comment: 15 pages, 8 figures. Updated to accepted version

We present new spectroscopic observations that are part of our continuing monitoring campaign of 88 quasars at z<0.7 whose broad H$\beta$ lines are offset from their systemic redshifts by a few thousand km/s. These quasars have been considered candidates for hosting supermassive black hole binaries (SBHBs) by analogy with single-lined spectroscopic binary stars. We present the data and describe our improved analysis techniques, which include an extensive evaluation of uncertainties. We also present a variety of measurements from the spectra that are of general interest and will be useful in later stages of our analysis. Additionally, we take this opportunity to study the variability of the optical continuum and integrated flux of the broad H$\beta$ line. We compare the variability properties of the SBHB candidates to those of a sample of typical quasars with similar redshifts and luminosities observed multiple times during the Sloan Digital Sky Survey. We find that the variability properties of the two samples are similar (variability amplitudes of 10-30% on time scales of approximately 1-7 years) and that their structure functions can be described by a common model with parameters characteristic of typical quasars. These results suggest that the broad-line regions of SBHB candidates have a similar extent as those of typical quasars. We discuss the implications of this result for the SBHB scenario and ensuing constraints on the orbital parameters., Comment: 19 pages in ApJ format. Accepted for publication in the Astrophysical Journal Supplements