Clustering measurements of obscured and unobscured quasars show that obscured quasars reside in more massive dark matter halos than their unobscured counterparts. These results are inconsistent with simple unified (torus) scenarios, but might be explained by models in which the distribution of obscuring material depends on Eddington ratio or galaxy stellar mass. We test these possibilities by constructing simple physical models to compare to observed AGN populations. We find that previously observed relationships between obscuration and Eddington ratio or stellar mass are not sufficient reproduce the observed quasar clustering results ($\langle \log M_{\text{halo}}/M_{\odot} \rangle = 12.94 ^{+ 0.10}_{- 0.11}$ and $\langle \log M_{\text{halo}}/M_{\odot} \rangle = 12.49 ^{+ 0.08}_{- 0.08}$ for obscured and unobscured populations, respectively) while maintaining the observed fraction of obscured quasars (30-65$\%$). This work suggests that evolutionary models, in which obscuration evolves on the typical timescale for black hole growth, are necessary to understand the observed clustering of mid-IR selected quasars.

We present a simple semi-numerical model designed to explore black hole growth and galaxy evolution. This method builds on a previous model for black hole accretion that uses a semi-numerical galaxy formation model and universal Eddington ratio distribution to describe the full AGN population by independently connecting galaxy and AGN growth to the evolution of the host dark matter halos. We fit observed X-ray luminosity functions up to a redshift of z ~ 4, as well as investigate the evolution of the Eddington ratio distributions. We find that the Eddington ratio distribution evolves with redshift such that the slope of the low-Eddington accretion rate distribution increases with cosmic time, consistent with the behavior predicted in hydrodynamical simulations for galaxies with different gas fractions. We also find that the evolution of our average Eddington ratio is correlated with observed star formation histories, supporting a picture in which black holes and galaxies evolve together in a global sense. We further confirm the impact of luminosity limits on observed galaxy and halo properties by applying selection criteria to our fiducial model and comparing to surveys across a wide range of redshifts.

A primary aim of the Nuclear Spectroscopic Telescope Array (NuSTAR) mission is to find and characterize heavily obscured Active Galactic Nuclei (AGNs). Based on mid-infrared photometry from the Wide-Field Infrared Survey Explorer (WISE) and optical photometry from the Sloan Digital Sky Surveys, we have selected a large population of luminous obscured AGN (i.e., obscured quasars). Here we report NuSTAR observations of four WISE-selected heavily obscured quasars for which we have optical spectroscopy from the Southern African Large Telescope and W. M. Keck Observatory. Optical diagnostics confirm that all four targets are AGNs. With NuSTAR hard X-ray observations, three of the four objects are undetected, while the fourth has a marginal detection. We confirm that these objects have observed hard X-ray (10-40 keV) luminosities at or below ~ 10^43 erg s^-1. We compare X-ray and IR luminosities to obtain estimates of the hydrogen column densities (N_H) based on the suppression of the hard X-ray emission. We estimate N_H of these quasars to be at or larger than 10^25 cm^-2, confirming that WISE and optical selection can identify very heavily obscured quasars that may be missed in X-ray surveys, and do not contribute significantly to the cosmic X-ray background. From the optical Balmer decrements, we found that our three extreme obscured targets lie in highly reddened host environments. This galactic extinction is not adequate to explain the more obscured AGN, but it may imply a different scale of obscuration in the galaxy.

We report on the results of active galactic nuclei (AGNs) detection by NuSTAR performed in three extragalactic survey fields (COSMOS, UDS, ECDFS) in three hard bands, namely H1 (8-16 keV), H2 (16-24 keV) and VH (35-55 keV). The aggregated area of the surveys is $\sim 2.7$ deg$^2$. While a large number of sources is detected in the H1 band (72 at the $97\%$ level of reliability), the H2 band directly probing close to the peak of the Cosmic X-ray Background (CXB) returns four significant detections, and two tentative, although not significant, detections are found in the VH band. All the sources detected above 16 keV are also detected at lower energies. We compute the integral number counts for sources in such bands, which show broad consistency with population synthesis models of the CXB. We furthermore identify two Compton-thick AGNs, one in the COSMOS field, associated with a hard and faint Chandra source, and one in the UDS field, never detected in the X-ray band before. Both sources are at the same redshift $z \sim 1.25$, which shifts their Compton-hump into the H1 band, and were previously missed in the usually employed NuSTAR bands, confirming the potential of using the H1 band to discover obscured AGNs at $z > 1$ in deep surveys.

The rare close projection of two quasars on the sky provides the opportunity to study the host galaxy environment of a foreground quasar in absorption against the continuum emission of a background quasar. For over a decade the "Quasars probing quasars" series has utilized this technique to further the understanding of galaxy formation and evolution in the presence of a quasar at z>2, resolving scales as small as a galactic disc and from bound gas in the circumgalactic medium to the diffuse environs of intergalactic space. Presented here, is the public release of the quasar pair spectral database utilized in these studies. In addition to projected pairs at z>2, the database also includes quasar pair members at z<2, gravitational lens candidates and quasars closely separated in redshift that are useful for small-scale clustering studies. In total the database catalogs 5627 distinct objects, with 4083 lying within 5' of at least one other source. A spectral library contains 3582 optical and near-infrared spectra for 3028 of the cataloged sources. As well as reporting on 54 newly discovered quasar pairs, we outline the key contributions made by this series over the last ten years, summarize the imaging and spectroscopic data used for target selection, discuss the target selection methodologies, describe the database content and explore some avenues for future work. Full documentation for the spectral database, including download instructions, are supplied at the following URL http://specdb.readthedocs.io/en/latest/

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 model the projected angular two-point correlation function (2PCF) of obscured and unobscured quasars selected using the Wide-field Infrared Survey Explorer (WISE), at a median redshift of $z \sim 1$ using a five-parameter Halo Occupation Distribution (HOD) parameterization, derived from a cosmological hydrodynamic simulation by Chatterjee et al. The HOD parameterization was previously used to model the 2PCF of optically selected quasars and X-ray bright active galactic nuclei (AGN) at $z \sim 1$. The current work shows that a single HOD parameterization can be used to model the population of different kinds of AGN in dark matter halos suggesting the universality of the relationship between AGN and their host dark matter halos. Our results show that the median halo mass of central quasar hosts increases from optically selected ($4.1^{+0.3}_{-0.4} \times 10^{12} \; h^{-1} \; {M_{sun}}$) and infra-red (IR) bright unobscured populations ($6.3^{+6.2}_{-2.3} \times 10^{12} \; h^{-1} \; {M_{sun}}$) to obscured quasars ($10.0^{+2.6}_{-3.7} \times 10^{12} \; h^{-1} \; {M_{sun}}$), signifying an increase in the degree of clustering. The projected satellite fractions also increase from optically bright to obscured quasars and tend to disfavor a simple `orientation only' theory of active galactic nuclei unification. Our results also show that future measurements of the small-scale clustering of obscured quasars can constrain current theories of galaxy evolution where quasars evolve from an IR- bright obscured phase to the optically bright unobscured phase.

We present the results of an optical spectroscopic survey of 46 heavily obscured quasar candidates. Objects are selected using their mid-infrared (mid-IR) colours and magnitudes from the Wide-Field Infrared Survey Explorer (WISE) and their optical magnitudes from the Sloan Digital Sky Survey (SDSS). Candidate Active Galactic Nuclei (AGNs) are selected to have mid-IR colours indicative of quasar activity and lie in a region of mid-IR colour space outside previously published X-ray based selection regions. We obtain optical spectra for our sample using the Robert Stobie Spectrograph on the Southern African Large Telescope. Thirty objects (65%) have identifiable emission lines, allowing for the determination of spectroscopic redshifts. Other than one object at $z\sim2.6$, candidates have moderate redshifts ranging from $z=0.1$ to $0.8$ with a median of 0.3. Twenty-one (70%) of our objects with identified redshift (46% of the whole sample) are identified as AGNs through common optical diagnostics. We model the spectral energy distributions of our sample and found that all require a strong AGN component, with an average intrinsic AGN fraction at 8$\,\mu$m of 0.91. Additionally, the fits require large extinction coefficients with an average $E(B-V)_\textrm{AGN} = 17.8$ (average $A(V)_\textrm{AGN} = 53.4$). By focusing on the area outside traditional mid-IR photometric cuts, we are able to capture and characterise a population of deeply buried quasars that were previously unattainable through X-ray surveys alone.

We present observed mid-infrared and optical colors and composite spectral energy distributions (SEDs) of type 1 (broad-line) and 2 (narrow-line) quasars selected from Sloan Digital Sky Survey (SDSS) spectroscopy. A significant fraction of powerful quasars are obscured by dust, and are difficult to detect in optical photometric or spectroscopic surveys. However these may be more easily identified on the basis of mid-infrared (MIR) colors and SEDs. Using samples of SDSS type 1 type 2 matched in redshift and [OIII] luminosity, we produce composite rest-frame 0.2-15 micron SEDs based on SDSS, UKIDSS, and Wide-Field Infrared Survey Explorer (WISE) photometry and perform model fits using simple galaxy and quasar SED templates. The SEDs of type 1 and 2 quasars are remarkably similar, with the differences explained primarily by the extinction of the quasar component in the type 2 systems. For both types of quasar, the flux of the AGN relative to the host galaxy increases with AGN luminosity (L_[OIII]) and redder observed MIR color, but we find only weak dependencies of the composite SEDs on mechanical jet power as determined through radio luminosity. We conclude that luminous quasars can be effectively selected using simple MIR color criteria similar to those identified previously (W1-W2 > 0.7 [Vega]), although these criteria miss many heavily obscured objects. Obscured quasars can be further identified based on optical-IR colors (for example, (u-W3 [AB]) > 1.4(W1-W2 [Vega])+3.2). These results illustrate the power of large statistical studies of obscured quasars selected on the basis of mid-IR and optical photometry.

In studies of the connection between active galactic nuclei (AGN) and their host galaxies there is widespread disagreement on some key aspects stemming largely from a lack of understanding of the nature of the full underlying AGN population. Recent attempts to probe this connection utilize both observations and simulations to correct for a missed population, but presently are limited by intrinsic biases and complicated models. We take a simple simulation for galaxy evolution and add a new prescription for AGN activity to connect galaxy growth to dark matter halo properties and AGN activity to star formation. We explicitly model selection effects to produce an "observed" AGN population for comparison with observations and empirically motivated models of the local universe. This allows us to bypass the difficulties inherent in many models which attempt to infer the AGN population by inverting selection effects. We investigate the impact of selecting AGN based on thresholds in luminosity or Eddington ratio on the "observed" AGN population. By limiting our model AGN sample in luminosity, we are able to recreate the observed local AGN luminosity function and specific star formation-stellar mass distribution, and show that using an Eddington ratio threshold introduces less bias into the sample by selecting the full range of growing black holes, despite the challenge of selecting low mass black holes. We find that selecting AGN using these various thresholds yield samples with different AGN host galaxy properties.

Recent work has found evidence for a difference in the bias and dark matter halo masses of WISE-selected obscured and unobscured quasars, implying a distinction between these populations beyond random line-of-sight effects. However, the significance of this difference in the most up-to-date measurements is relatively weak, at $\sim$2$\sigma$ for individual measurements but bolstered by agreement from different techniques, including angular clustering and cross-correlations with cosmic microwave background (CMB) lensing maps. Here, we expand the footprint of previous work, aiming to improve the precision of both methods. In this larger area we correct for position dependent selection effects, in particular fluctuations of the WISE-selected quasar density as a function of Galactic latitude. We also measure the cross-correlation of the obscured and unobscured samples and confirm that they are well-matched in redshift, both centred at $z=1$. Combined with very similar detection fractions and magnitude distributions in the long-wavelength WISE bands, this redshift match strongly supports the fact that IR selection identifies obscured and unobscured quasars of similar bolometric luminosity. Finally, we perform cross-correlations with confirmed spectroscopic quasars, again confirming the results from other methods --- obscured quasars reside in haloes a factor of 3 times more massive than unobscured quasars. This difference is significant at the $\sim$5$\sigma$ level when the measurements are combined, strong support for the idea that obscuration in at least some quasars is tied to the larger environment, and may have an evolutionary component.

We present the results of a multi-wavelength study of a sample of high-redshift Radio Loud (RL) Broad Absorption Line (BAL) quasars. This way we extend to higher redshift previous studies on the radio properties, and broadband optical colors of these objects. We have se- lected a sample of 22 RL BAL quasars with 3.6 z 4.8 cross-correlating the FIRST radio survey with the SDSS. Flux densities between 1.25 and 9.5 GHz have been collected with the JVLA and Effelsberg-100m telescopes for 15 BAL and 14 non-BAL quasars used as compar- ison sample. We determine the synchrotron peak frequency, constraining their age. A large number of GigaHertz Peaked Spectrum (GPS) and High Frequency Peakers (HFP) sources has been found in both samples (80% for BAL and 71% for non-BAL QSOs), not suggesting a younger age for BAL quasars. The spectral index distribution provides information about the orientation of these sources, and we find statistically similar distributions for the BAL and non-BAL quasars in contrast to work done on lower redshift samples. Our sample may be too small to convincingly find the same effect, or might represent a real evolutionary effect based on the large fraction of young sources. We also study the properties of broadband colors in both optical (SDSS) and near- and mid-infrared (UKIDSS and WISE) bands, finding that also at high redshift BAL quasars tend to be optically redder than non-BAL quasars. However, these differences are no more evident at longer wavelength, when comparing colors of the two samples by mean of the WISE survey.

Recent measurements of the dark matter halo masses of infrared-selected obscured quasars are in tension --- some indicate that obscured quasars have higher halo mass compared to their unobscured counterparts, while others find no difference. The former result is inconsistent with the simplest models of quasar unification that rely solely on viewing angle, while the latter may support such models. Here, using empirical relationships between dark matter halo and supermassive black hole masses, we provide a simple evolutionary picture that naturally explains these findings and is motivated by more sophisticated merger-driven quasar fueling models. The model tracks the growth rate of haloes, with the black hole growing in spurts of quasar activity in order to "catch-up" with the $M_{\textrm{BH}}$ - $M_{\textrm{star}}$ - $M_{\textrm{halo}}$ relationship. The first part of the quasar phase is obscured and is followed by an unobscured phase. Depending on the luminosity limit of the sample, driven by observational selection effects, a difference in halo masses may or may not be significant. For high luminosity samples, the difference can be large (a few to 10 times higher masses in obscured quasars), while for lower luminosity samples the halo mass difference is very small, much smaller than current observational constraints. Such a simple model provides a qualitative explanation for the higher mass haloes of obscured quasars, as well as rough quantitative agreement with seemingly disparate results.

An important question in extragalactic astronomy concerns the distribution of black hole accretion rates of active galactic nuclei (AGN). Based on observations at X-ray wavelengths, the observed Eddington ratio distribution appears as a power law, while optical studies have often yielded a lognormal distribution. There is increasing evidence that these observed discrepancies may be due to contamination by star formation and other selection effects. Using a sample of galaxies from the Sloan Digital Sky Survey Data Release 7, we test if an intrinsic Eddington ratio distribution that takes the form of a Schechter function is consistent with previous work that suggests that young galaxies in optical surveys have an observed lognormal Eddington ratio distribution. We simulate the optical emission line properties of a population of galaxies and AGN using a broad instantaneous luminosity distribution described by a Schechter function near the Eddington limit. This simulated AGN population is then compared to observed galaxies via the positions on an emission line excitation diagram and Eddington ratio distributions. We present an improved method for extracting the AGN distribution using BPT diagnostics that allows us to probe over one order of magnitude lower in Eddington ratio counteracting the effects of dilution by star formation. We conclude that for optically selected AGN in young galaxies, the intrinsic Eddington ratio distribution is consistent with a possibly universal, broad power law with an exponential cutoff, as this distribution is observed in old optically selected galaxies and in X-rays.

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

Using the most recent releases of WISE and Planck data, we perform updated measurements of the bias and typical dark matter halo mass of infrared-selected obscured and unobscured quasars, using the angular autocorrelation function and cosmic microwave background (CMB) lensing cross-correlations. Since our recent work of this kind, the WISE Allwise catalogue was released with improved photometry, and the Planck mission was completed and released improved products. These new data provide a more reliable measurement of the quasar bias and provide an opportunity to explore the role of changing survey pipelines in results downstream. We present a comparison of IR color-selected quasars, split into obscured and unobscured populations based on optical-IR colors, selected from two versions of the WISE data. Which combination of data is used impacts the final results, particularly for obscured quasars, both because of mitigation of some systematics and because the newer catalogue provides a slightly different sample. We show that Allwise data is superior in several ways, though there may be some systematic trends with Moon contamination that were not present in the previous catalogue. We opt currently for the most conservative sample that meet our selection criteria in both the previous and new WISE catalogues. We measure a higher bias and halo mass for obscured quasars ($b_{\textrm{obsc}} \sim 2.1$, $b_{\textrm{unob}} \sim 1.8$) --- at odds with simple orientation models --- but at a reduced significance ($\sim$1.5$\sigma$) as compared to our work with previous survey data.

Extreme deconvolution (XD) of broad-band photometric data can both separate stars from quasars and generate probability density functions for quasar redshifts, while incorporating flux uncertainties and missing data. Mid-infrared photometric colors are now widely used to identify hot dust intrinsic to quasars, and the release of all-sky WISE data has led to a dramatic increase in the number of IR-selected quasars. Using forced-photometry on public WISE data at the locations of SDSS point sources, we incorporate this all-sky data into the training of the XDQSOz models originally developed to select quasars from optical photometry. The combination of WISE and SDSS information is far more powerful than SDSS alone, particularly at $z>2$. The use of SDSS$+$WISE photometry is comparable to the use of SDSS$+$ultraviolet$+$near-IR data. We release a new public catalogue of 5,537,436 (total; 3,874,639 weighted by probability) potential quasars with probability $P_{\textrm{QSO}} > 0.2$. The catalogue includes redshift probabilities for all objects. We also release an updated version of the publicly available set of codes to calculate quasar and redshift probabilities for various combinations of data. Finally, we demonstrate that this method of selecting quasars using WISE data is both more complete and efficient than simple WISE color-cuts, especially at high redshift. Our fits verify that above $z \sim 3$ WISE colors become bluer than the standard cuts applied to select quasars. Currently, the analysis is limited to quasars with optical counterparts, and thus cannot be used to find highly obscured quasars that WISE color-cuts identify in significant numbers.

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.

We cross-correlate a cosmic microwave background (CMB) lensing map with the projected space densities of quasars to measure the bias and halo masses of a quasar sample split into obscured and unobscured populations, the first application of this method to distinct quasar subclasses. Several recent studies of the angular clustering of obscured quasars have shown that these objects likely reside in higher-mass halos compared to their unobscured counterparts. This has important implications for models of the structure and geometry of quasars, their role in growing supermassive black holes, and mutual quasar/host galaxy evolution. However, the magnitude and significance of this difference has varied from study to study. Using data from \planck, \wise, and SDSS, we follow up on these results using the independent method of CMB lensing cross-correlations. The region and sample are identical to that used for recent angular clustering measurements, allowing for a direct comparison of the CMB-lensing and angular clustering methods. At $z \sim 1$, we find that the bias of obscured quasars is $b_q = 2.57 \pm 0.24$, while that of unobscured quasars is $b_q = 1.89 \pm 0.19$. This corresponds to halo masses of $\log (M_h / M_{\odot} h^{-1}) = 13.24_{-0.15}^{+0.14}$ (obscured) and $\log (M_h / M_{\odot} h^{-1}) = 12.71_{-0.13}^{+0.15}$ (unobscured). These results agree well with with those from angular clustering (well within $1\sigma$), and confirm that obscured quasars reside in host halos $\sim$3 times as massive as halos hosting unobscured quasars. This implies that quasars spend a significant portion of their lifetime in an obscured state, possibly more than one half of the entire active phase.

We present the results of an optical spectroscopic survey of a sample of 40 candidate obscured quasars identified on the basis of their mid-infrared emission detected by the Wide-Field Infrared Survey Explorer (WISE). Optical spectra for this survey were obtained using the Robert Stobie Spectrograph (RSS) on the Southern African Large Telescope (SALT). Our sample was selected with WISE colors characteristic of AGNs, as well as red optical to mid-IR colors indicating that the optical/UV AGN continuum is obscured by dust. We obtain secure redshifts for the majority of the objects that comprise our sample (35/40), and find that sources that are bright in the WISE W4 (22$\mu$m) band are typically at moderate redshift (<z> = 0.35$) while sources fainter in W4 are at higher redshifts (<z> = 0.73$). The majority of the sources have narrow emission lines, with optical colors and emission line ratios of our WISE-selected sources that are consistent with the locus of AGN on the rest-frame $g-z$ color vs. [NeIII]$\lambda$3869 / [OII]$\lambda\lambda$3726+3729 line ratio diagnostic diagram. We also use empirical AGN and galaxy templates to model the spectral energy distributions (SEDs) for the objects in our sample, and find that while there is significant variation in the observed SEDs for these objects, the majority require a strong AGN component. Finally, we use the results from our analysis of the optical spectra and the SEDs to compare our selection criteria to alternate criteria presented in the literature. These results verify the efficacy of selecting luminous obscured AGNs based on their WISE colors.

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.

Recent studies of luminous infrared-selected active galactic nuclei (AGN) suggest that the reddest, most obscured objects display a higher angular clustering amplitude, and thus reside in higher-mass dark matter halos. This is a direct contradiction to the prediction of the simplest unification-by-orientation models of AGN and quasars. However, clustering measurements depend strongly on the "mask" that removes low-quality data and describes the sky and selection function. We find that applying a robust, conservative mask to WISE-selected quasars yields a weaker but still significant difference in the bias between obscured and unobscured quasars. These findings are consistent with results from previous Spitzer surveys, and removes any scale dependence of the bias. For obscured quasars with $\langle z \rangle = 0.99$ we measure a bias of $b_q = 2.67 \pm 0.16$, corresponding to a halo mass of $\log (M_h / M_{\odot} h^{-1}) = 13.3 \pm 0.1$, while for unobscured sources with $\langle z \rangle = 1.04$ we find $b_q = 2.04 \pm 0.17$ with a halo mass $\log (M_h / M_{\odot} h^{-1} )= 12.8 \pm 0.1$. This improved measurement indicates that WISE-selected obscured quasars reside in halos only a few times more massive than the halos of their unobscured counterparts, a reduction in the factor of $\sim$10 larger halo mass as has been previously reported using WISE-selected samples. Additionally, an abundance matching analysis yields lifetimes for both obscured and unobscured quasar phases on the order of a few 100 Myr ($\sim$ 1\% of the Hubble time) --- however, the obscured phase lasts roughly twice as long, in tension with many model predictions.

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 report the results of a new 60 ks Chandra X-ray Observatory Advanced CCD Imaging Spectrometer S-array (ACIS-S) observation of the reddened, radio-selected, highly polarized `FeLoBAL' quasar FIRST J1556+3517. We investigated a number of models of varied sophistication to fit the 531-photon spectrum. These models ranged from simple power laws to power laws absorbed by hydrogen gas in differing ionization states and degrees of partial covering. Preferred fits indicate that the intrinsic X-ray flux is consistent with that expected for quasars of similarly high luminosity, i.e., an intrinsic, dereddened and unabsorbed optical to X-ray spectral index of -1.7. We cannot tightly constrain the intrinsic X-ray power-law slope, but find indications that it is flat (photon index Gamma = 1.7 or flatter at a >99% confidence for a neutral hydrogen absorber model). Absorption is present, with a column density a few times 10^23 cm^-2, with both partially ionized models and partially covering neutral hydrogen models providing good fits. We present several lines of argument that suggest the fraction of X-ray emissions associated with the radio jet is not large. We combine our Chandra data with observations from the literature to construct the spectral energy distribution of FIRST J1556+3517 from radio to X-ray energies. We make corrections for Doppler beaming for the pole-on radio jet, optical dust reddening, and X-ray absorption, in order to recover a probable intrinsic spectrum. The quasar FIRST J1556+3517 seems to be an intrinsically normal radio-quiet quasar with a reddened optical/UV spectrum, a Doppler-boosted but intrinsically weak radio jet, and an X-ray absorber not dissimilar from that of other broad absorption line quasars.

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 present new spectropolarimetric observations of 8 radio-loud broad absorption line (BAL) quasars, and combine these new data with our previous spectropolarimetric atlases (of both radio-loud and radio-quiet objects) in order to investigate the polarization properties of BAL quasars as a group. The total (radio-selected) sample includes 36 (26) high-ionization and 22 (15) low-ionization BAL quasars. On average, we confirm that broad emission lines are polarized at a level similar to or less than the continuum and broad absorption troughs are more highly polarized, but we note that these properties are not true for all individual objects. Of the whole sample, 18 (31%) have high (>2%) continuum polarization, including 45% of the LoBALs and 22% of HiBALs. We identify a few correlations between polarization and other quasar properties, as well as some interesting non-correlations. In particular, continuum polarization does not correlate with radio spectral index, which suggests that the polarization is not due to a standard geometry and preferred viewing angle to BAL quasars. The polarization also does not correlate with the amount of intrinsic dust reddening, indicating that the polarization is not solely due to direct light attenuation either. Polarization does appear to depend on the minimum BAL outflow velocity, confirming the results of previous studies, and it may correlate with the maximum outflow velocity. We also find that continuum polarization anti-correlates with the polarization in the C IV broad emission and broad absorption. These results suggest that the polarization of BAL quasars cannot be described by one simple model, and that the scatterer location and geometry can vary significantly from object to object.

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.

It was recently shown that there is a significant difference in the radio spectral index distributions of broad absorption line (BAL) quasars and unabsorbed quasars, with an overabundance of BAL quasars with steeper radio spectra. This result suggests that source orientation does play into the presence or absence of BAL features. In this paper we provide more quantitative analysis of this result based on Monte-Carlo simulations. While the relationship between viewing angle and spectral index does indeed contain a lot of scatter, the spectral index distributions are different enough to overcome that intrinsic variation. Utilizing two different models of the relationship between spectral index and viewing angle, the simulations indicate that the difference in spectral index distributions can be explained by allowing BAL quasar viewing angles to extend about 10 degrees farther from the radio jet axis than non-BAL sources, though both can be seen at small angles. These results show that orientation cannot be the only factor determining whether BAL features are present, but it does play a role.

New multi-frequency radio observations of a large sample of radio-selected BAL quasars, along with a very well matched sample of normal quasars, are presented. The observations were made one immediately after the other at 4.9 and 8.4 GHz with the goal of measuring the radio spectral index of each source. We have identified, for the first time, a significant difference in the spectral index distributions of BAL versus non-BAL quasars, with BAL sources showing an overabundance of steep-spectrum sources. This is the first direct observation suggesting that BAL quasars are more likely to be seen farther from the radio jet axis, although a range of orientations is needed to explain the width of the distribution. Utilizing a few different relationships between spectral index and viewing angle, we have also performed Monte-Carlo simulations to quantify the viewing angle to these sources. We find that the difference in the distributions of spectral index can be explained by allowing the BAL sources to have viewing angles extending about 10 degrees farther from the jet axis than non-BAL quasars.

We have built a sample of 74 radio-selected broad absorption line quasars from the Sloan Digital Sky Survey Data Release 5 (SDSS DR5) and Faint Images of the Radio Sky at Twenty Centimeters (FIRST), along with a well matched sample of 74 unabsorbed "normal" quasars. The sources have been observed with the NRAO Very Large Array/Expanded Very Large Array at 8.4 GHz (3.5 cm) and 4.9 GHz (6 cm). All sources have additional archival 1.4 GHz (21 cm) data. Here we present the measured radio fluxes, spectral indices, and our initial findings. The percentage of BAL quasars with extended structure (on the order of 10%) in our sample is similar to previous studies at similar resolutions, suggesting that BAL quasars are indeed generally compact, at least at arsecond resolutions. The majority of sources do not appear to be significantly variable at 1.4 GHz, but we find two previously unidentified BAL quasars that may fit into the "polar" BAL category. We also identify a significant favoring of steeper radio spectral index for BAL compared to non-BAL quasars. This difference is apparent for several different measures of the spectral index, and persists even when restricting the samples to only include compact objects. Because radio spectral index is a statistical indicator of viewing angle for large samples, these results suggest that BAL quasars do have a range of orientations but are more often observed farther from the jet axis compared to normal quasars.

We present spectropolarimetry of 19 confirmed and 4 possible bright, southern broad absorption line (BAL) quasars from the European Southern Observatory (ESO) Very Large Telescope (VLT). A wide range of redshifts is covered in the sample (from 0.9 to 3.4), and both low- and hi-ionization quasars are represented, as well as radio-loud and radio-quiet BALQSOs. We continue to confirm previously established spectropolarimetric properties of BALQSOs, including the generally rising continuum polarization with shorter wavelengths and comparatively large fraction with high broad-band polarization (6 of 19 with polarizations >2%). Emission lines are polarized less than or similar to the continuum, except in a few unusual cases, and absorption troughs tend to have higher polarizations. A search for correlations between polarization properties has been done, identifying 2 significant or marginally significant correlations. These are an increase in continuum polarization with decreasing optical luminosity (increasing absolute B magnitude) and decreasing \ion{C}{4} emission-line polarization with increased continuum polarization.

We report spectropolarimetry of 30 radio-selected broad absorption line (BAL) quasars with the Keck Observatory, 25 from the sample of Becker et al. (2000). Both high and low-ionization BAL quasars are represented, with redshifts ranging from 0.5 to 2.5. The spectropolarimetric properties of radio-selected BAL quasars are very similar to those of radio-quiet BAL quasars: a sizeable fraction (20%) show large continuum polarization (2-10%) usually rising toward short wavelengths, emission lines are typically less polarized than the continuum, and absorption line troughs often show large polarization jumps. There are no significant correlations between polarization properties and radio properties, including those indicative of system orientation, suggesting that BAL quasars are not simply normal quasars seen from an edge-on perspective.