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 (PTAs) 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 pulsar timing arrays, in providing an independent statistical test of the nature of SMBHB candidates discovered in electromagnetic searches.

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.

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.

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

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.

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.

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.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

We analyze a set of optical-to-near-infrared long-slit nuclear spectra of 16 infrared-luminous spiral galaxies. All of the studied sources present H$_2$ emission, which reflects the star-forming nature of our sample, and they clearly display H I emission lines in the optical. Their continua contain many strong stellar absorption lines, with the most common features due to Ca I, Ca II, Fe I, Na I, Mg I, in addition to prominent absorption bands of TiO, VO, ZrO, CN and CO. We report a homogeneous set of equivalent width (EW) measurements for 45 indices, from optical to NIR species for the 16 star-forming galaxies as well as for 19 early type galaxies where we collected the data from the literature. This selected set of emission and absorption-feature measurements can be used to test predictions of the forthcoming generations of stellar population models. We find correlations among the different absorption features and propose here correlations between optical and NIR indices, as well as among different NIR indices, and compare them with model predictions. While for the optical absorption features the models consistently agree with the observations,the NIR indices are much harder to interpret. For early-type spirals the measurements agree roughly with the models, while for star-forming objects they fail to predict the strengths of these indices.

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.

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.

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.

Active galactic nuclei (AGN) that show strong rest-frame optical/UV variability in their blue continuum and broad line emission are classified as "changing-look" AGN, or at higher luminosities changing look quasars (CLQs). These surprisingly large and sometimes rapid transitions challenge accepted models of quasar physics and duty cycles, offer several new avenues for study of quasar host galaxies, and open a wider interpretation of the cause of differences between broad and narrow line AGN. To better characterize extreme quasar variability, we present follow-up spectroscopy as part of a comprehensive search for CLQs across the full SDSS footprint using spectroscopically confirmed quasars from the SDSS DR7 catalog. Our primary selection requires large-amplitude (|\Delta g|>1 mag, |\Delta r|>0.5 mag) variability over any of the available time baselines probed by the SDSS and Pan-STARRS 1 surveys. We employ photometry from the Catalina Sky Survey to verify variability behavior in CLQ candidates where available, and confirm CLQs using optical spectroscopy from the William Herschel, MMT, Magellan, and Palomar telescopes. For our adopted S/N threshold on variability of broad H\beta emission, we find 17 new CLQs, yielding a confirmation rate of >~ 20%. These candidates are at lower Eddington ratio relative to the overall quasar population which supports a disk-wind model for the broad line region. Based on our sample, the CLQ fraction increases from 10% to roughly half as the continuum flux ratio between repeat spectra at 3420 Angstroms increases from 1.5 to 6. We release a catalog of over 200 highly variable candidates to facilitate future CLQ searches.

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.

We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007.

The majority of low-ionization nuclear emission-line regions (LINERs) harbor supermassive black holes with very low accretion rates. However, the accretion flows do not produce enough ionizing photons to power the emission lines emitted on scales of ~100 pc, and therefore additional sources of power are required. We present and analyze Hubble Space Telescope spectra of three nearby luminous LINERs that are spatially resolved on scales of < 9 pc. The targets have multiple indicators of an accreting black hole, as well as a deficient ionizing photon budget. We measure diagnostic emission line ratios as a function of distance from the nucleus and compare them to models for different excitation mechanisms: shocks, photoionization by the accreting black hole, and photoionization by young or old hot stars. We also consider the kinematics of the line-emitting gas, as revealed by the widths and shifts of the emission lines. We conclude that, in LINERs with low-luminosity active nuclei, shocks by jets or other outflows are crucial in exciting the gas in and around the nucleus, as suggested by other authors. The physical model that best describes our targets comprises a low-luminosity, accretion-powered active nucleus that photoionizes the gas within ~20 pc of the galaxy center, and shock excitation of the gas at larger distances.

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.

We explore the kinematics of ionized gas via the [O III] $\lambda$5007 emission lines in active galactic nuclei (AGN) selected on the basis of their mid-infrared (IR) emission, and split into obscured and unobscured populations based on their optical-IR colors. After correcting for differences in redshift distributions, we provide composite spectra of spectroscopically and photometrically defined obscured/Type 2 and unobscured/Type 1 AGN from 3500 to 7000 \AA. The IR-selected obscured sources contain a mixture of narrow-lined Type 2 AGN and intermediate sources that have broad H$\alpha$ emission and significantly narrower H$\beta$. Using both [OIII] luminosities and AGN luminosities derived from optical-IR spectral energy distribution fitting, we find evidence for enhanced large-scale obscuration in the obscured sources. In matched bins of luminosity we find that the obscured population typically has broader, more blueshifted \OIII\ emission than in the unobscured sample, suggestive of more powerful AGN-driven outflows. This trend is not seen in spectroscopically classified samples, and is unlikely to be entirely explained by orientation effects. In addition, outflow velocities increase from small to moderate AGN $E(B-V)$ values, before flattening out (as traced by FWHM) and even decreasing (as traced by blueshift). While difficult to fully interpret in a single physical model, due to both the averaging over populations and the spatially-averaged spectra, these results agree with previous findings that simple geometric unification models are insufficient for the IR-selected AGN population, and may fit into an evolutionary model for obscured and unobscured AGN.

We present optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC\,2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find lags consistent with geometrically thin accretion-disk models that predict a lag-wavelength relation of $\tau \propto \lambda^{4/3}$. However, the observed lags are larger than predictions based on standard thin-disk theory by factors of 3.3 for MCG+08-11-011 and 2.3 for NGC\,2617. These differences can be explained if the mass accretion rates are larger than inferred from the optical luminosity by a factor of 4.3 in MCG+08-11-011 and a factor of 1.3 in NGC\,2617, although uncertainty in the SMBH masses determines the significance of this result. While the X-ray variability in NGC\,2617 precedes the UV/optical variability, the long 2.6 day lag is problematic for coronal reprocessing models.

We present the first results from an optical reverberation mapping campaign executed in 2014, targeting the active galactic nuclei (AGN) MCG+08-11-011, NGC 2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting observational properties, including a "changing look" AGN and a broad-line radio galaxy. Based on continuum-H$\beta$ lags, we measure black hole masses for all five targets. We also obtain H$\gamma$ and He{\sc ii}\,$\lambda 4686$ lags for all objects except 3C 382. The He{\sc ii}\,$\lambda 4686$ lags indicate radial stratification of the BLR, and the masses derived from different emission lines are in general agreement. The relative responsivities of these lines are also in qualitative agreement with photoionization models. These spectra have extremely high signal-to-noise ratios (100--300 per pixel) and there are excellent prospects for obtaining velocity-resolved reverberation signatures.

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.

As astronomers increasingly exploit the information available in the time domain, spectroscopic variability in particular opens broad new channels of investigation. Here we describe the selection algorithms for all targets intended for repeat spectroscopy in the Time Domain Spectroscopic Survey (TDSS), part of the extended Baryon Oscillation Spectroscopic Survey within the Sloan Digital Sky Survey-IV. Also discussed are the scientific rationale and technical constraints leading to these target selections. The TDSS includes a large "Repeat Quasar Spectroscopy" (RQS) program delivering ~13,000 repeat spectra of confirmed SDSS quasars, and several smaller "Few-Epoch Spectroscopy" (FES) programs targeting specific classes of quasars as well as stars. The RQS program aims to provide a large and diverse quasar data set for studying variations in quasar spectra on timescales of years, a comparison sample for the FES quasar programs, and opportunity for discovering rare, serendipitous events. The FES programs cover a wide variety of phenomena in both quasars and stars. Quasar FES programs target broad absorption line quasars, high signal-to-noise ratio normal broad line quasars, quasars with double-peaked or very asymmetric broad emission line profiles, binary supermassive black hole candidates, and the most photometrically variable quasars. Strongly variable stars are also targeted for repeat spectroscopy, encompassing many types of eclipsing binary systems, and classical pulsators like RR Lyrae. Other stellar FES programs allow spectroscopic variability studies of active ultracool dwarf stars, dwarf carbon stars, and white dwarf/M dwarf spectroscopic binaries. We present example TDSS spectra and describe anticipated sample sizes and results.

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.

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.

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.

We report the detection of extended Lyman-$\alpha$ emission from the host galaxy of SDSS~J2222+2745, a strongly lensed quasar at $z = 2.8$. Spectroscopic follow-up clearly reveals extended Lyman-$\alpha$ in emission between two images of the central active galactic nucleus (AGN). We reconstruct the lensed quasar host galaxy in the source plane by applying a strong lens model to HST imaging, and resolve spatial scales as small as $\sim$200 parsecs. In the source plane we recover the host galaxy morphology to within a few hundred parsecs of the central AGN, and map the extended Lyman-$\alpha$ emission to its physical origin on one side of the host galaxy at radii $\sim$0.5-2 kpc from the central AGN. There are clear morphological differences between the Lyman-$\alpha$ and rest-frame ultraviolet stellar continuum emission from the quasar host galaxy. Furthermore, the relative velocity profiles of quasar Lyman-$\alpha$, host galaxy Lyman-$\alpha$, and metal lines in outflowing gas reveal differences in the absorbing material affecting the AGN and host galaxy. These data indicate the presence of patchy local intervening gas in front of the central quasar and its host galaxy. This interpretation is consistent with the central luminous quasar being obscured across a substantial fraction of its surrounding solid angle, resulting in strong anisotropy in the exposure of the host galaxy to ionizing radiation from the AGN. This work demonstrates the power of strong lensing-assisted studies to probe spatial scales that are currently inaccessible by other means.

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.

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.

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.

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.

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$.

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.

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.

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.

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.

The masses of the black holes powering quasars represent a fundamental parameter of active galaxies. Estimates of quasar black hole masses using single-epoch spectra are quite uncertain, and require quantitative improvement. We recently identified a correction for C IV $\lambda$1549-based scaling relationships used to estimate quasar black hole masses that relies on the continuum-subtracted peak flux ratio of the ultraviolet emission-line blend Si IV + OIV] (the $\lambda$1400 feature) to that of C IV. This parameter correlates with the suite of associated quasar spectral properties collectively known as "Eigenvector 1" (EV1). Here we use a sample of 85 quasars with quasi-simultaneous optical-ultraviolet spectrophotometry to demonstrate how biases in the average EV1 properties can create systematic biases in C IV-based black hole mass scaling relationships. This effect results in nearly an order of magnitude moving from objects with small $<$ peak $\lambda$1400/C IV $>$, which have overestimated black hole masses, to objects with large $<$ peak $\lambda1400$/C IV $>$, which have underestimated values. We show that existing reverberation-mapped samples of quasars with ultraviolet spectra -- used to calibrate C IV-based scaling relationships -- have significant EV1 biases that result in predictions of black hole masses nearly 50\% too high for the average quasar. We offer corrections and suggestions to account for this bias.

Studies of radio-loud (RL) broad absorption line (BAL) quasars indicate that popular orientation-based BAL models fail to account for all observations. Are these results extendable to radio-quiet (RQ) BAL quasars? Comparisons of RL and RQ BAL quasars show that many of their properties are quite similar. Here we extend these analyses to the rest-frame ultraviolet (UV) spectral properties, using a sample of 73 RL and 473 RQ BAL quasars selected from the Sloan Digital Sky Survey (SDSS). Each RQ quasar is individually matched to a RL quasar in both redshift (over the range $1.5 < z < 3.5$) and continuum luminosity. We compare several continuum, emission line, and absorption line properties, as well as physical properties derived from these measurements. Most properties in the samples are statistically identical, though we find slight differences in the velocity structure of the BALs that cause apparent differences in CIV emission line properties. Differences in the velocities may indicate an interaction between the radio jets and the absorbing material. We also find that UV FeII emission is marginally stronger in RL BAL quasars. All of these differences are subtle, so in general we conclude that RL and RQ BAL QSOs are not fundamentally different objects, except in their radio properties. They are therefore likely to be driven by similar physical phenomena, suggesting that results from samples of RL BAL quasars can be extended to their RQ counterparts.

We analyze orbital solutions for 48 massive multiple-star systems in the Cygnus OB2 Association, 23 of which are newly presented here, to find that the observed distribution of orbital periods is approximately uniform in log P for P<45 d, but it is not scale-free. Inflections in the cumulative distribution near 6 d, 14, d, and 45 d, suggest key physical scales of ~0.2, ~0.4, and ~1 A.U. where yet-to-be-identified phenomena create distinct features. No single power law provides a statistically compelling prescription, but if features are ignored, a power law with exponent beta = -0.22 provides a crude approximation over P=1.4 -- 2000 d, as does a piece-wise linear function with a break near 45 d. The cumulative period distribution flattens at P > 45 d, even after correction for completeness, indicating either a lower binary fraction or a shift toward low-mass companions. A high degree of similarity (91% likelihood) between the Cyg OB2 period distribution and that of other surveys suggests that the binary properties at P<25 d are determined by local physics of disk/clump fragmentation and are relatively insensitive to environmental and evolutionary factors. Fully 30% of the unbiased parent sample is a binary with period P < 45 d. Completeness corrections imply a binary fraction near 55% for P < 5000 d. The observed distribution of mass ratios 0.2<q<1 is consistent with uniform, while the observed distribution of eccentricities 0.1<e<0.6 is consistent with uniform plus an excess of e ~ 0 systems. We identify six stars, all supergiants, that exhibit aperiodic velocity variations of ~30 km/s attributed to atmospheric fluctuations.

Quasar luminosity functions are a fundamental probe of the growth and evolution of supermassive black holes. Measuring the intrinsic luminosity function is difficult in practice, due to a multitude of observational and systematic effects. As sample sizes increase and measurement errors drop, characterizing the systematic effects is becoming more important. It is well known that the continuum emission from the accretion disk of quasars is anisotropic --- in part due to its disk-like structure --- but current luminosity function calculations effectively assume isotropy over the range of unobscured lines of sight. Here, we provide the first steps in characterizing the effect of random quasar orientations and simple models of anisotropy on observed luminosity functions. We find that the effect of orientation is not insignificant and exceeds other potential corrections such as those from gravitational lensing of foreground structures. We argue that current observational constraints may overestimate the intrinsic luminosity function by as much as a factor of ~2 on the bright end. This has implications for models of quasars and their role in the Universe, such as quasars' contribution to cosmological backgrounds.

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]

We investigate the orientation dependence of the spectral energy distributions in a sample of radio-loud quasars. Selected specifically to study orientation issues, the sample contains 52 sources with redshifts in the range 0.16<z<1.4 and measurements of radio core dominance, a radio orientation indicator. Measured properties include monochromatic luminosities at a range of wavelengths between the infrared and X-rays, integrated infrared luminosity, spectral slopes, and the covering fraction of the obscuring circumnuclear dust. We estimate dust covering fraction assuming that the accretion disk emits anisotropically and discuss the shortcomings and technical difficulties of this calculation. Luminosities are found to depend on orientation, with face-on sources factors of a 2-3 brighter than more edge-on sources, depending on wavelength. The degree of anisotropy varies very little with wavelength such that the overall shape of the spectral energy distribution does not vary significantly with orientation. In the infrared, we do not observe a decrease in anisotropy with increasing wavelength. The spectral slopes and estimates of covering fraction are not significantly orientation dependent. We construct composite spectral energy distributions as a function of radio core dominance and find that these illustrate the results determined from the measured properties.

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.

Currently, the ability to produce black hole mass estimates using the CIV 1549 line that are consistent with Hb mass estimates is uncertain, due in large part to disagreement between velocity line width measurements for the two lines. This discrepancy between Hb and CIV arises from the fact that both line profiles are treated the same way in single-epoch scaling relationships based on the assumption that the broad-line region is virialized, even though a non-virialized emission component is often significant in the CIV line and absent or weak in the Hb line. Using quasi-simultaneous optical and ultra-violet spectra for a sample of 85 optically bright quasars with redshifts in the range z=0.03-1.4, we present a significant step along the path to rehabilitating the CIV line for black hole mass estimates. We show that the residuals of velocity line width between CIV and Hb are significantly correlated with the peak flux ratio of SiIV+OIV] 1400 to CIV. Using this relationship, we develop a prescription for estimating black hole masses from the ultra-violet spectrum that agree better with Hb-based masses than the traditional CIV masses. The scatter between Hb and CIV masses is initially 0.43 dex in our sample and is reduced to 0.33 dex when using our prescription. The peak flux ratio of SiIV+OIV] 1400 to CIV is an ultraviolet indicator of the suite of spectral properties commonly known as "Eigenvector 1", thus the reduction in scatter between CIV and Hb black hole masses is essentially due to removing an Eigenvector 1 bias in CIV-based masses.

We have obtained IRTF/SpeX spectra of eight moderate-redshift (z=0.7-2.4), radio-selected (log R*~0.4-1.9) broad absorption line (BAL) quasars. The spectra cover the rest-frame optical band. We compare the optical properties of these quasars to those of canonically radio-quiet (log R*<1) BAL quasars at similar redshifts and to low-redshift quasars from the Palomar-Green catalog. As with previous studies of BAL quasars, we find that [OIII] 5007 is weak, and optical FeII emission is strong, a rare combination in canonically radio-loud (log R*>1) quasars. With our measurements of the optical properties, particularly the Balmer emission line widths and the continuum luminosity, we have used empirical scaling relations to estimate black hole masses and Eddington ratios. These lie in the range (0.4-2.6)x10^9 M_sun and 0.1-0.9, respectively. Despite their comparatively extreme radio properties relative to most BAL quasars, their optical properties are quite consistent with those of radio-quiet BAL quasars and dissimilar to those of radio-loud non-BAL quasars. While BAL quasars generally appear to have low values of [OIII]/FeII, an extreme of "Eigenvector 1", the Balmer line widths and Eddington ratios do not appear to significantly differ from those of unabsorbed quasars at similar redshifts and luminosities.

If broad absorption line (BAL) quasars represent a high covering fraction evolutionary state (even if this is not the sole factor governing the presence of BALs), it is expected that they should show an excess of mid-infrared radiation compared to normal quasars. Some previous studies have suggested that this is not the case. We perform the first analysis of the IR properties of radio-loud BAL quasars, using IR data from WISE and optical (rest-frame ultraviolet) data from SDSS, and compare the BAL quasar sample with a well-matched sample of unabsorbed quasars. We find a statistically significant excess in the mid- to near-infrared luminosities of BAL quasars, particularly at rest-frame wavelengths of 1.5 and 4 microns. Our sample was previously used to show that BALs are observed along many lines of sight towards quasars, but with an overabundance of more edge-on sources, suggesting that orientation factors into the appearance of BALs. The evidence here---of a difference in IR luminosities between BAL quasars and unabsorbed quasars---may be ascribed to evolution. This suggests that a merging of the current BAL paradigms is needed to fully describe the class.

We present optical spectroscopy of a sample of 38 post-starburst quasars (PSQs) at z ~ 0.3, 29 of which have morphological classifications based on Hubble Space Telescope imaging. These broad-lined active galactic nuclei (AGNs) possess the spectral signatures of massive intermediate-aged stellar populations making them potentially useful for studying connections between nuclear activity and host galaxy evolution. We model the spectra in order to determine the ages and masses of the host stellar populations, and the black hole masses and Eddington fractions of the AGNs. Our model components include an instantaneous starburst, a power-law, and emission lines. We find the PSQs have MBH ~ 10^8 Msun accreting at a few percent of Eddington luminosity and host ~ 10^10.5 Msun stellar populations which are several hundred Myr to a few Gyr old. We investigate relationships among these derived properties, spectral properties, and morphologies. We find that PSQs hosted in spiral galaxies have significantly weaker AGN luminosities, older starburst ages, and narrow emission-line ratios diagnostic of ongoing star-formation when compared to their early-type counterparts. We conclude that the early-type PSQs are likely the result of major mergers and were likely luminous infrared galaxies in the past, while spiral PSQs with more complex star-formation histories are triggered by less dramatic events (e.g., harassment, bars). We provide diagnostics to distinguish the early-type and spiral hosts when high spatial resolution imaging is not available.

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.

We recently presented radio observations of a large sample of radio-loud broad absorption line (BAL) quasars from the SDSS and FIRST surveys, as well as a well matched sample of unabsorbed quasars, primarily to measure their radio spectral indices and estimate ensemble orientations. Here, we analyze the SDSS spectra of these samples and compare the rest-frame ultraviolet properties of radio-loud BAL and non-BAL quasars. Ultraviolet properties include the continuum shape, emission-line measurements of C IV, Al III, C III], Fe II, and Mg II, and BAL properties including the balnicity index (BI), absorption index (AI), and minimum and maximum outflow velocities. We find that radio-loud BAL quasars have similar ultraviolet properties compared to radio-loud non-BAL sources, though they do appear to have redder continua and stronger Fe II emission, which is consistent with what is found for radio-quiet BAL sources. No correlations exist between outflow properties and orientation (radio spectral index), suggesting that BAL winds along any line of sight are driven by the same mechanisms. There are also few correlations between spectral index and other properties. We conclude that BAL outflows occur along all lines of sight with similar strengths and velocities.

We present measurements and statistical properties of the optical and ultraviolet emission lines present in the spectra of 85 bright quasars which have detailed spectral energy distributions. This heterogeneous sample has redshifts up to z=1.5 and is comprised of three subsamples that may be of particular utility: ultraviolet excess Palomar-Green quasars, quasars with far-ultraviolet coverage from FUSE, and radio-loud quasars selected to have similar extended radio luminosity originally selected for orientation studies. Most of the objects have quasi-simultaneous optical-ultraviolet spectra, with significant coverage in the radio-to-X-ray wavebands. The parameters of all strong emission lines are measured by detailed spectral fitting. Many significant correlations previously found among quasar emission-line properties are also present in this sample, e.g., the Baldwin effect, the optical correlations collectively known as eigenvector 1, and others. Finally, we use our measurements plus scaling relationships to estimate black hole masses and Eddington fractions. We show the mass estimates from different emission lines are usually in agreement within a factor of 2, but nearly a third show larger differences. We suggest using multiple mass scaling relationships to estimate black hole masses when possible, and adopting a median of the estimates as the black hole mass for individual objects. Line measurements and derived AGN properties will be used for future studies examining the relationships among quasar emission lines and their spectral energy distributions.

We present infrared bolometric luminosity corrections derived from the detailed spectral energy distributions of 62 bright quasars of low- to moderate-redshift (z=0.03-1.4). At 1.5, 2, 3, 7, 12, 15, and 24 microns we provide bolometric corrections of the mathematical forms L_iso=\zeta \lambda L_\lambda and log(L_iso)=A+B log(\lambda L_\lambda). Bolometric corrections for radio-loud and radio-quiet objects are consistent within 95% confidence intervals, so we do not separate them. Bolometric luminosities estimated using these corrections are typically smaller than those derived from some commonly used in the literature. We investigate the possibility of a luminosity dependent bolometric correction and find that, while the data are consistent with such a correction, the dispersion is too large and the luminosity range too small to warrant such a detailed interpretation. Bolometric corrections at 1.5 $\mu$m are appropriate for objects with properties that fall in the range log(L_bol)=45.4-47.3 and bolometric corrections at all other wavelengths are appropriate for objects with properties that fall in the range log(L_bol)=45.1-47.0.

Massive binary stars may constitute a substantial fraction of progenitors to supernovae and gamma-ray bursts, and the distribution of their orbital characteristics holds clues to the formation process of massive stars. As a contribution to securing statistics on OB-type binaries, we report the discovery and orbital parameters for five new systems as part of the Cygnus OB2 Radial Velocity Survey. Four of the new systems (MT070, MT174, MT267, and MT734 (a.k.a. VI Cygni #11) are single-lined spectroscopic binaries while one (MT103) is a double-lined system (B1V+B2V). MT070 is noteworthy as the longest period system yet measured in Cyg OB2, with P=6.2 yr. The other four systems have periods ranging between 4 and 73 days. MT174 is noteworthy for having a probable mass ratio q<0.1, making it a candidate progenitor to a low-mass X-ray binary. These measurements bring the total number of massive binaries in Cyg OB2 to 25, the most currently known in any single cluster or association.

Bolometric corrections are used in quasar studies to quantify total energy output based on a measurement of a monochromatic luminosity. First, we enumerate and discuss the practical difficulties of determining such corrections, then we present bolometric luminosities between 1 \mu m and 8 keV rest frame and corrections derived from the detailed spectral energy distributions of 63 bright quasars of low to moderate redshift (z = 0.03-1.4). Exploring several mathematical fittings, we provide practical bolometric corrections of the forms L_iso=\zeta \lambda L_{\lambda} and log(L_iso)=A+B log(\lambda L_{\lambda}) for \lambda= 1450, 3000, and 5100 \AA, where L_iso is the bolometric luminosity calculated under the assumption of isotropy. The significant scatter in the 5100 \AA\ bolometric correction can be reduced by adding a first order correction using the optical slope, \alpha_\lambda,opt. We recommend an adjustment to the bolometric correction to account for viewing angle and the anisotropic emission expected from accretion discs. For optical/UV monochromatic luminosities, radio-loud and radio-quiet bolometric corrections are consistent within 95% confidence intervals so we do not make separate radio-loud and radio-quiet corrections. In addition, we provide several bolometric corrections to the 2-10 keV X-ray luminosity, which are shown to have very large scatter. Separate radio-loud and radio-quiet corrections are warranted by the X-ray data.

We have produced the next generation of quasar spectral energy distributions (SEDs), essentially updating the work of Elvis et al. (1994) by using high-quality data obtained with several space and ground-based telescopes, including NASA's Great Observatories. We present an atlas of SEDs of 85 optically bright, non-blazar quasars over the electromagnetic spectrum from radio to X-rays. The heterogeneous sample includes 27 radio-quiet and 58 radio-loud quasars. Most objects have quasi-simultaneous ultraviolet-optical spectroscopic data, supplemented with some far-ultraviolet spectra, and more than half also have Spitzer mid-infrared IRS spectra. The X-ray spectral parameters are collected from the literature where available. The radio, far-infrared, and near-infrared photometric data are also obtained from either the literature or new observations. We construct composite spectral energy distributions for radio-loud and radio-quiet objects and compare these to those of Elvis et al., finding that ours have similar overall shapes, but our improved spectral resolution reveals more detailed features, especially in the mid and near-infrared.