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 and 2 matched in redshift and [O III] luminosity, we produce composite rest-frame 0.2-15 mu m SEDs based on SDSS, UKIDSS, and Wide-field Infrared Survey Explorer 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 active galactic nucleus (AGN) relative to the host galaxy increases with AGN luminosity (L-[O III]) 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 MIR and optical photometry.

This is the second paper in a series on a new luminous z similar to 5 quasar survey using optical and near-infrared colors. Here we present a new determination of the bright end of the quasar luminosity function (QLF) at z similar to 5. Combining. our 45 new quasars with previously known quasars that satisfy our selections, we construct the largest uniform luminous z similar to 5 quasar sample to date, with 99 quasars in the range of 4.7 <= z < 5.4 and -29 < M-1450 <= -26.8, within the Sloan Digital Sky Survey (SDSS) footprint. We use a modified 1/V-a method including flux limit correction to derive a binned QLF, and we model the parametric QLF using maximum likelihood estimation. With the faint-end slope of the QLF fixed as alpha = -2.03 from previous deeper samples, the best fit of our QLF gives a flatter bright end slope beta = -3.58 +/- 0.24 and a fainter break magnitude M-1450(*) = -26.98 +/- 0.23 than previous studies at similar redshift. Combined with previous work at lower and higher redshifts, our result is consistent with a luminosity evolution and density evolution model. Using the best-fit QLF, the contribution of quasars to the ionizing background at z similar to 5 is found to be 18%-45% with a clumping factor C of 2-5. Our sample suggests an evolution of radio loud fraction with optical luminosity but no obvious evolution with redshift.

We present initial results from the first systematic survey of luminous z similar to 5.5 quasars. Quasars at z similar to 5.5, the post-reionization epoch, are crucial tools to explore the evolution of intergalactic medium, quasar evolution, and the early super-massive black hole growth. However, it has been very challenging to select quasars at redshifts 5.3 <= z <= 5.7 using conventional color selections, due to their similar optical colors to late-type stars, especially M dwarfs, resulting in a glaring redshift gap in quasar redshift distributions. We develop a new selection technique for z similar to 5.5 quasars based on optical, near-IR, and mid-IR photometric data from Sloan Digital Sky Survey (SDSS), UKIRT InfraRed Deep Sky Surveys-Large Area Survey (ULAS), VISTA Hemisphere Survey (VHS), and Wide Field Infrared Survey Explorer. From our pilot observations in the SDSS-ULAS/VHS area, we have discovered 15 new quasars at 5.3. z. 5.7 and 6 new lower redshift quasars, with SDSS z band magnitude brighter than 20.5. Including other two z similar to 5.5 quasars already published in our previous work, we now construct a uniform quasar sample at 5.3 <= z <= 5.7, with 17 quasars in a similar to 4800 square degree survey area. For further application in a larger survey area, we apply our selection pipeline to do a test selection by using the new wide field J-band photometric data from a preliminary version of the UKIRT Hemisphere Survey (UHS). We successfully discover the first UHS selected z similar to 5.5 quasar.

We present the discovery of nine quasars at z similar to 6 identified in the Sloan Digital Sky Survey (SDSS) imaging data. This completes our survey of z similar to 6 quasars in the SDSS footprint. Our final sample consists of 52 quasars at 5.7 < z << 6.4, including 29 quasars with z(AB) <= 20 mag selected from 11,240 deg(2) of the SDSS single-epoch imaging survey (the main survey), 10 quasars with 20 <= z(AB) <= 20.5 selected from 4223 deg2 of the SDSS overlap regions (regions with two or more imaging scans), and 13 quasars down to z(AB) approximate to 22 mag from the 277 deg2 in Stripe 82. They span a wide luminosity range of -29.0 <= M-1450 <= -24.5. This well-defined sample is used to derive the quasar luminosity function (QLF) at z similar to 6. After combining our SDSS sample with two faint (M-1450 >= -23 mag) quasars from the literature, we obtain the parameters for a double power-law fit to the QLF. The bright-end slope beta of the QLF is well constrained to be beta = -2.8 +/- 0.2. Due to the small number of low-luminosity quasars, the faint-end slope a and the characteristic magnitude M*(1450) are less well constrained, with alpha = -1.90(-0.44)(+0.58) and M* = -25.2(-3.8)(+1.2) mag. The spatial density of luminous quasars, parametrized as rho(M-1450 < -26, z) = rho(z = 6)10(k(z-6)), drops rapidly from z similar to 5 to 6, with k = -0.72 +/- 0.11. Based on our fitted QLF and assuming an intergalactic medium (IGM) clumping factor of C = 3, we find that the observed quasar population cannot provide enough photons to ionize the z similar to 6 IGM at similar to 90% confidence. Quasars may still provide a significant fraction of the required photons, although much larger samples of faint quasars are needed for more stringent constraints on the quasar contribution to reionization.

Luminous high-redshift quasars can be used to probe of the intergalactic medium in the early universe because their UV light is absorbed by the neutral hydrogen along the line of sight. They help us to measure the neutral hydrogen fraction of the high-z universe, shedding light on the end of reionization epoch. In this paper, we present a discovery of a new quasar (PSO J006.1240+39.2219) at redshift z = 6.61 +/- 0.02 from Panoramic Survey Telescope & Rapid Response System 1.Including this quasar, there are nine quasars above z > 6.5 up to date. The estimated continuum brightness is M-1450 = -25.96 +/- 0.08. PSO J006.1240+39.2219 has a strong Ly alpha emission compared with typical low-redshift quasars, but the measured near-zone region size is R-NZ = 3.2 +/- 1.1 proper megaparsecs, which is consistent with other quasars at z similar to 6.

We report the results of a faint quasar survey in a one-square-degree field. The aim is to test the Y - K/g - z and J - K/i - Y color selection criteria for quasars at faint magnitudes to obtain a complete sample of quasars based on deep optical and near-infrared color-color selection and to measure the faint end of the quasar luminosity function (QLF) over a wide redshift range. We carried out a quasar survey based on the Y - K/g - z and J - K/i - Y quasar selection criteria, using the deep Y-band data obtained from our CFHT/WIRCam Y-band images in a two-degree field within the F22 field of the VIMOS VLT deep survey, optical co-added data from Sloan Digital Sky Survey Stripe 82 and deep near-infrared data from the UKIDSS Deep Extragalactic Survey in the same field. We discovered 25 new quasars at 0.5 < z < 4.5 and i < 22.5 mag within one-square-degree field. The survey significantly increases the number of faint quasars in this field, especially at z similar to 2-3. It confirms that our color selections are highly complete in a wide redshift range (z < 4.5), especially over the quasar number density peak at z similar to 2-3, even for faint quasars. Combining all previous known quasars and new discoveries, we construct a sample with 109 quasars and measure the binned QLF and parametric QLF. Although the sample is small, our results agree with a pure luminosity evolution at lower redshift and luminosity evolution and density evolution model at redshift z > 2.5.

We present a new algorithm to estimate quasar photometric redshifts (photo-zs), by considering the asymmetries in the relative flux distributions of quasars. The relative flux models are built with multivariate Skew-t distributions in the multidimensional space of relative fluxes as a function of redshift and magnitude. For 151,392 quasars in the SDSS, we achieve a photo-z accuracy, defined as the fraction of quasars with the difference between the photo-z z(p) and the spectroscopic redshift z(s), vertical bar Delta z vertical bar=vertical bar z(s)-z(p)vertical bar/(1 + z(s)) within 0.1, of 74%. Combining the WISE W1 and W2 infrared data with the SDSS data, the photo-z accuracy is enhanced to 87%. Using the Pan-STARRS1 or DECaLS photometry with WISE W1 and W2 data, the photo-z accuracies are 79% and 72%, respectively. The prior probabilities as a function of magnitude for quasars, stars, and galaxies are calculated, respectively, based on (1) the quasar luminosity function, (2) the Milky Way synthetic simulation with the Besancon model, and (3) the Bayesian Galaxy Photometric Redshift estimation. The relative fluxes of stars are obtained with the Padova isochrones, and the relative fluxes of galaxies are modeled through galaxy templates. We test our classification method to select quasars using the DECaLS g, r, z, and WISE W1 and W2 photometry. The quasar selection completeness is higher than 70% for a wide redshift range 0.5 < z < 4.5, and a wide magnitude range 18 < r < 21.5 mag. Our photo-z regression and classification method has the potential to extend to future surveys. The photo-z code will be publicly available.

We present a survey of the [C II] 158 mu m line and underlying far-infrared (FIR) dust continuum emission in a sample of 27 greater than or similar to 6 quasars using the Atacama Large Millimeter Array (ALMA) at similar to 1 '' resolution. The [C II] line was significantly detected (at > 5-sigma) in 23 sources (85%). We find typical line luminosities of L-[C (II]) = 10(9-10) L-circle dot, and an average line width of similar to 385 km s(-1). The [C II]-to-far-infrared luminosity ratios ([C II]/FIR) in our sources span one order of magnitude, highlighting a variety of conditions in the star-forming medium. Four quasar host galaxies are clearly resolved in their [C II] emission on a few kpc scales. Basic estimates of the dynamical masses of the host galaxies give masses between 2 x 10(10) and 2 x 10(11) M-circle dot, i.e., more than an order of magnitude below what is expected from local scaling relations, given the available limits on the masses of the central black holes (> 3 x 10(8) M-circle dot, assuming Eddington-limited accretion). In stacked ALMA [C II] spectra of individual sources in our sample, we find no evidence of a deviation from a single Gaussian profile. The quasar luminosity does not strongly correlate with either the [C II] luminosity or equivalent width. This survey (with typical on-source integration times of 8 minutes) showcases the unparalleled sensitivity of ALMA at millimeter wavelengths, and offers a unique reference sample for the study of the first massive galaxies in the universe.

We present reverberation mapping results from the first year of combined spectroscopic and photometric observations of the Sloan Digital Sky Survey Reverberation Mapping Project. We successfully recover reverberation time delays between the g+i band emission and the broad H beta emission line for a total of 44 quasars, and for the broad Ha emission line in 18 quasars. Time delays are computed using the JAVELIN and CREAM software and the traditional interpolated cross-correlation function (ICCF): using well-defined criteria, we report measurements of 32 H beta and 13 Ha lags with JAVELIN, 42 H beta and 17 Ha lags with CREAM, and 16 H beta and eight Ha lags with the ICCF. Lag values are generally consistent among the three methods, though we typically measure smaller uncertainties with JAVELIN and CREAM than with the ICCF, given the more physically motivated light curve interpolation and more robust statistical modeling of the former two methods. The median redshift of our H beta-detected sample of quasars is 0.53, significantly higher than that of the previous reverberation mapping sample. We find that in most objects, the time delay of the Ha emission is consistent with or slightly longer than that of H beta. We measure black hole masses using our measured time delays and line widths for these quasars. These black hole mass measurements are mostly consistent with expectations based on the local M-BH-sigma* relationship, and are also consistent with single-epoch black hole mass measurements. This work increases the current sample size of reverberation-mapped active galaxies by about two-thirds and represents the first large sample of reverberation mapping observations beyond the local universe (z < 0.3).

There is a large degree of variety in the optical variability of quasars and it is unclear whether this is all attributable to a single (set of) physical mechanism(s). We present the results of a systematic search for major flares in active galactic nucleus (AGN) in the Catalina Real-time Transient Survey as part of a broader study into extreme quasar variability. Such flares are defined in a quantitative manner as being atop of the normal, stochastic variability of quasars. We have identified 51 events from over 900 000 known quasars and high-probability quasar candidates, typically lasting 900 d and with a median peak amplitude of Delta m = 1.25 mag. Characterizing the flare profile with a Weibull distribution, we find that nine of the sources are well described by a single-point single-lens model. This supports the proposal by Lawrence et al. that microlensing is a plausible physical mechanism for extreme variability. However, we attribute the majority of our events to explosive stellar-related activity in the accretion disc: superluminous supernovae, tidal disruption events and mergers of stellar mass black holes.