Does the Compact Radio Jet in PG 1700+518 Drive a Molecular Outflow?

Radio jets play an important role in quasar feedback, but direct observations showing how the jets interact with the multi-phase interstellar medium of galaxy disks are few and far between. In this work, we provide new millimeter interferometric observations of PG 1700+518 in order to investigate the effect of its radio jet on the surrounding molecular gas. PG 1700 is a radio-quiet, low-ionization broad absorption line quasar whose host galaxy has a nearby interacting companion. On sub-kiloparsec scales, the ionized gas is driven to high velocities by a compact radio jet that is identified by radio interferometry. We present observations from the NOrthern Extended Millimeter Array (NOEMA) interferometer with a 3.8 arcsec (16 kpc) synthesized beam where we detect the CO (1-0) emission line at $30\sigma$ significance with a total flux of $3.12\pm0.02$ Jy km s$^{-1}$ and a typical velocity dispersion of $125\pm5$ km s$^{-1}$. Despite the outflow in ionized gas, we find no concrete evidence that the CO gas is being affected by the radio jet on size scales of a kiloparsec or more. However, a $\sim\!1$ arcsec drift in the spatial centroid of the CO emission as a function of velocity across the emission line and the compact nature of the jet hint that higher spatial resolution observations may reveal a signal of interaction between the jet and molecular gas.Comment: Updated version with minor changes after last round of review. Accepted to Ap

Orientation and quasar black hole mass estimation

We have constructed a sample of 386 radio-loud quasars with z < 0.75 from the Sloan Digital Sky Survey in order to investigate orientation effects on black hole mass estimates. Orientation is estimated using radio core dominance measurements based on FIRST survey maps. Black hole masses are estimated from virial-based scaling relationships using H-beta, and compared to the stellar velocity dispersion (sigma_*), predicted using the Full Width at Half Maximum (FWHM) of [O III] 5007, which tracks mass via the M-sigma_* relation. We find that the FWHM of Hbeta correlates significantly with radio core dominance and biases black hole mass determinations that use it, but that this is not the case for sigma_* based on [O III] 5007. The ratio of black hole masses predicted using orientation-biased and unbiased estimates, which can be determined for radio-quiet as well as radio-loud quasars, is significantly correlated with radio core dominance. Although there is significant scatter, this mass ratio calculated in this way may in fact serve as an orientation estimator. We additionally note the existence of a small population radio core-dominated quasars with extremely broad H-beta emission lines that we hypothesise may represent recent black hole mergers.Comment: 9 pages, 6 figures, accepted by MNRA

Updating quasar bolometric luminosity corrections. III. [O III] bolometric corrections

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<z<1.0002$). We adopted two functional forms to calculate $L_{\textrm{iso}}$, the bolometric luminosity determined under the assumption of isotropy: $L_{\textrm{iso}}=A\,L_{[O\,III]}$ for comparison with the literature and log$(L_{iso})=B+C\,$log$(L_{[O\,III]})$, which better characterizes the data. We also explored whether "Eigenvector 1", which describes the range of quasar spectral properties and quantifies their diversity, introduces scatter into the $L_{[O\,III]}-L_{iso}$ relationship. We found that the [O III] bolometric correction can be significantly improved by adding a term including the equivalent width ratio $R_{Fe\,II}\equiv EW_{Fe\,II}/EW_{H\beta}$, which is an Eigenvector 1 indicator. Inclusion of $R_{Fe\,II}$ in predicting $L_{iso}$ is significant at nearly the $3\sigma$ level and reduces the scatter and systematic offset of the luminosity residuals. Typically, [O III] bolometric corrections are adopted for Type 2 sources where the quasar continuum is not observed and in these cases, $R_{Fe\,II}$ cannot be measured. We searched for an alternative measure of Eigenvector 1 that could be measured in the optical spectra of Type 2 sources but were unable to identify one. Thus, the main contribution of this work is to present an improved [O III] bolometric correction based on measured bolometric luminosities and highlight the Eigenvector 1 dependence of the correction in Type 1 sources.Comment: Accepted for publication in MNRA

Orientation and accretion in a representative sample of active galactic nuclei

We highlight a representative sample of active galactic nuclei selected independent of orientation. The defining characteristic of the selection is sophisticated matching between the $0.1<z<0.6$ Sloan Digital Sky Survey quasars from the Seventh Data Release to the Westerbork Northern Sky Survey at 325 MHz and the subsequent application of a total radio luminosity cut. The resulting sample is complete down to the limiting luminosity and unbiased by orientation. Compared to orientation samples in the literature this approach yields less bias with redshift, relatively more lobe-dominated sources including those with radio lobes and no visible core, and a distribution of radio core dominance that is consistent with expectations from a uniform distribution of inclinations with solid angle. We measure properties of the optical spectra, and use the sample to investigate the orientation dependence of the velocity width of the broad H$\beta$ emission line. We recover the known orientation dependence, but the sharp envelope of previous studies where only edge-on sources display the broadest lines, is absent. Scatter in this diagram is not attributable solely to black hole mass, Eddington ratio, or contamination in the sample from compact steep spectrum sources. A physical framework for quasar beaming and a disk-like broad-line region can describe the representative sample when it is expanded to include additional parameters, in particular jet properties and the broad-line region velocity field. These points serve to illustrate the critical role of sample selection in the interpretation of observable properties as indicators of physical parameters of quasar central engines.Comment: 28 pages, 10 figures, 6 tables, submitted to Ap

The orientation dependence of quasar single-epoch black hole mass scaling relationships

Black hole masses are estimated for radio-loud quasars using several self-consistent scaling relationships based on emission-line widths and continuum luminosities. The emission lines used, H-beta, Mg II, and C IV, have different dependencies on orientation as estimated by radio core dominance. We compare differences in the log of black hole masses estimated from different emission lines and show that they depend on radio core dominance in the sense that core-dominated, jet-on objects have systematically smaller H-beta and Mg II determined masses compared to those from C IV, while lobe-dominated edge-on objects have systematically larger H-beta and Mg II determined masses compared to those from C IV. The effect is consistent with the H-beta line width, and to a lesser extent that of Mg II, being dependent upon orientation in the sense of a axisymmetric velocity field plus a projection effect. The size of the effect is nearly an order of magnitude in black hole mass going from one extreme orientation to the other. We find that radio spectral index is a good proxy for radio core dominance and repeating this analysis with radio spectral index yields similar results. Accounting for orientation could in principle significantly reduce the scatter in black hole mass scaling relationships, and we quantify and offer a correction for this effect cast in terms of radio core dominance and radio spectral index.Comment: 16 pages, 16 figures, accepted to MNRA

Does Size Matter? The Underlying Intrinsic Size Distribution of Radio Sources and Implications for Unification by Orientation

Unification by orientation is a ubiquitous concept in the study of active galactic nuclei. A gold standard of the orientation paradigm is the hypothesis that radio galaxies and radio-loud quasars are intrinsically the same, but are observed over different ranges of viewing angles. Historically, strong support for this model was provided by the projected sizes of radio structure in luminous radio galaxies, which were found to be significantly larger than those of quasars, as predicted due to simple geometric projection. Recently, this test of the simplest prediction of orientation-based models has been revisited with larger samples that cover wider ranges of fundamental properties---and no clear difference in projected sizes of radio structure is found. Cast solely in terms of viewing angle effects, these results provide convincing evidence that unification of these objects solely through orientation fails. However, it is possible that conflicting results regarding the role orientation plays in our view of radio sources simply result from insufficient sampling of their intrinsic size distribution. We test this possibility using Monte-Carlo simulations constrained by real sample sizes and properties. We develop models for the real intrinsic size distribution of radio sources, simulate observations by randomly sampling intrinsic sizes and viewing angles, and analyze how likely each sample is to support or dispute unification by orientation. We find that, while it is possible to reconcile conflicting results purely within a simple, orientation-based framework, it is very unlikely. We analyze the effects that sample size, relative numbers of radio galaxies and quasars, the critical angle that separates the two subclasses, and the shape of the intrinsic size distribution have on this type of test.Comment: 13 pages, 8 figues. Accepted for publication in Ap

The behavior of quasar C IV emission-line properties with orientation

With a quasar sample designed for studying orientation effects, we investigate the orientation dependence of characteristics of the C IV 1549 broad emission line in approximately 50 Type 1 quasars with z=0.1-1.4. Orientation is measured for the sample via radio core dominance. In our analysis we include measurements of the continuum luminosity and the optical-to-X-ray spectral slope, spectral properties commonly included in the suite known as "Eigenvector 1", and the full-width at half maximum, full-width at one-quarter-maximum, shape, blueshift, and equivalent width of the C IV broad emission line. We also investigate a new prescription that we recently developed for predicting the velocity line width of the H-beta broad emission line based on the velocity line width of the C IV line and the ratio of continuum subtracted peak fluxes of Si IV + O IV] at 1400 A to C IV. In addition to a correlation analysis of the ultraviolet spectral properties and radio core dominance, we provide composite spectra of edge-on and face-on sources for this sample. In particular, we highlight the orientation dependence of the velocity line width predicted for H-beta. We find that this predicted line width depends on orientation in a manner similar to the true velocity line width of H-beta, where no such dependence is observed for C IV. This is an indication that orientation information concerning the line emitting regions can be extracted from ultraviolet spectra. [abridged]Comment: 12 pages, 4 figures, 2 tables, accepted to MNRA

A Multi-wavelength Analysis of Binary-AGN Candidate PSO J334.2028+01.4075

We present analysis of the first Chandra observation of PSO J334.2028+01.4075 (PSO J334), targeted as a binary-AGN candidate based on periodic variations of the optical flux. With no prior targeted X-ray coverage for PSO J334, our new 40 ksec Chandra observation allows for the opportunity to differentiate between a single or binary-AGN system, and if a binary, can characterize the mode of accretion. Simulations show that the two expected accretion disk morphologies for binary-AGN systems are (i) a "cavity", where the inner region of the accretion disk is mostly empty and emission is truncated blueward of the wavelength associated with the temperature of the innermost ring, or (ii) "minidisks", where there is substantial accretion from the cirumbinary disk onto one or both of the members of the binary, each with their own shock-heated thin-disk accretion system. We find the X-ray emission to be well-fit with an absorbed power-law, incompatible with the simple cavity scenario. Further, we construct an SED of PSO J334 by combining radio through X-ray observations and find that the SED agrees well with that of a normal AGN, most likely incompatible with the minidisk scenario. Other analyses, such as locating the quasar on IR color-color diagrams and analyzing the quasar mass predicted by the fundamental plane of black hole activity, further highlight the similarity of PSO J334 with respect to normal AGN. On the multi-wavelength fronts we investigated, we find no evidence supporting PSO J334 as a binary-AGN system, though our analysis remains insensitive to some binary configurations.Comment: 9 pages, 2 figures, accepted to Ap

Likelihood for Detection of Sub-parsec Supermassive Black Hole Binaries in Spectroscopic Surveys

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 comment

AGN Triality of Triple Mergers: Detection of Faint X-ray Point Sources

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 Ap