Quasars are the most luminous of active galactic nuclei (AGN), and are perhaps responsible for quenching star formation in their hosts. The Stripe 82X catalog covers 31.3 deg$^2$ of the Stripe 82 field, of which the 15.6 deg$^2$ covered with XMM-Newton is also covered by Herschel/SPIRE. We have 2500 X-ray detected sources with multi-wavelength counterparts, and 30% of these are unobscured quasars, with $L_X > 10^{44}\,$erg/s and $M_B < -23$. We define a new population of quasars which are unobscured, have X-ray luminosities in excess of $10^{44}\,$erg/s, have broad emission lines, and yet are also bright in the far-infrared, with a 250$\mu$m flux density of $S_{\rm 250}>30$mJy. We refer to these Herschel-detected, unobscured quasars as "Cold Quasars". A mere 4% (21) of the X-ray- and optically-selected unobscured quasars in Stripe 82X are detected at 250$\mu$m. These Cold Quasars lie at $z\sim1-3$, have $L_{\rm IR}>10^{12}\,L_\odot$, and have star formation rates of $\sim200-1400\,M_\odot$/yr. Cold Quasars are bluer in the mid-IR than the full quasar population, and 72% of our Cold Quasars have WISE W3 $<$ 11.5 [Vega], while only 19% of the full quasar sample meets this criteria. Crucially, Cold Quasars have on average $\sim9\times$ as much star formation as the main sequence of star forming galaxies at similar redshifts. Although dust-rich, unobscured quasars have occasionally been noted in the literature before, we argue that they should be considered as a separate class of quasars due to their high star formation rates. This phase is likely short-lived, as the central engine and immense star formation consume the gas reservoir. Cold Quasars are type-1 blue quasars that reside in starburst galaxies.

We present the results of a systematic search for quasars in the Catalina Real-time Transient Survey exhibiting both strong photometric and spectroscopic variability over a decadal baseline. We identify 73 sources with specific patterns of optical and mid-IR photometric behavior and a defined spectroscopic change. These "Changing-State" quasars (CSQs) form a higher luminosity sample to complement existing sets of "Changing-Look" AGN and quasars in the literature. The CSQs (by selection) exhibit larger photometric variability than the CLQs. The spectroscopic variability is marginally stronger in the CSQs than CLQs as defined by the change in H$\beta$/[OIII] ratio. We find 36 sources with declining H$\beta$ flux, 37 sources with increasing H$\beta$ flux and discover seven sources with $z > 0.8$, further extending the redshift arm. Our CSQ sample compares to the literature CLQ objects in similar distributions of H$\beta$ flux ratios and differential Eddington ratios between high (bright) and low (dim) states. Taken as a whole, we find that this population of extreme varying quasars is associated with changes in the Eddington ratio and the timescales imply cooling/heating fronts propagating through the disk.

We present the discovery of a gravitationally lensed dust-reddened QSO at $z=2.517$ discovered in a survey for red QSOs by infrared selection. $Hubble~Space~Telescope$ imaging in the WFC3/IR F160W and F125W filters reveals a quadruply lensed system in a cusp configuration. We find that compared to the central image of the cusp, the nearby, brightest image is anomalous by a factor of $\sim7-11$. Although the source is extremely bright in the mid-infrared, a magnification by a factor of $\sim50-120$ places an upper limit of 1.35 mJy on the intrinsic mid-infrared brightness, well below the $WISE~W4$ detection limit of 6 mJy. We find that this QSO is moderately reddened, with $E(B-V)=0.7$ and that $\sim1\%$ of the intrinsic spectrum is leaked back into the line of sight resulting in an upturn in its UV spectrum. We conclude that the QSO's reddening is intrinsic and not due to the lens. Consistent with previous red quasar samples, this source exhibits outflows in its spectrum as well as morphological properties suggestive of it being in a merger-driven transitional phase. Depending on how $L_{\rm bol}$ is computed, the quasar's accretion rate may be as high as $0.26~L_{\rm Edd}$. We detect two Lyman limit systems, at $z=2.102$ and $z=2.431$, with absorption by metal lines likely at small impact parameter to the QSO, and a putative lens redshift of $z=0.599$. Given the rarity of quad lenses, the discovery of this source allows detailed study of a less luminous, more typical infrared-selected quasar at high redshift.

We present a spectroscopically complete sample of 147 infrared-color-selected AGN down to a 22 $\mu$m flux limit of 20 mJy over the $\sim$270 deg$^2$ of the SDSS Stripe 82 region. Most of these sources are in the QSO luminosity regime ($L_{\rm bol} \gtrsim 10^{12} L_\odot$) and are found out to $z\simeq3$. We classify the AGN into three types, finding: 57 blue, unobscured Type-1 (broad-lined) sources; 69 obscured, Type-2 (narrow-lined) sources; and 21 moderately-reddened Type-1 sources (broad-lined and $E(B-V) > 0.25$). We study a subset of this sample in X-rays and analyze their obscuration to find that our spectroscopic classifications are in broad agreement with low, moderate, and large amounts of absorption for Type-1, red Type-1 and Type-2 AGN, respectively. We also investigate how their X-ray luminosities correlate with other known bolometric luminosity indicators such as [O III] line luminosity ($L_{\rm [OIII]}$) and infrared luminosity ($L_{6 \mu{\rm m}}$). While the X-ray correlation with $L_{\rm [OIII]}$ is consistent with previous findings, the most infrared-luminous sources appear to deviate from established relations such that they are either under-luminous in X-rays or over-luminous in the infrared. Finally, we examine the luminosity function (LF) evolution of our sample, and by AGN type, in combination with the complementary, infrared-selected, AGN sample of Lacy et al. (2013), spanning over two orders of magnitude in luminosity. We find that the two obscured populations evolve differently, with reddened Type-1 AGN dominating the obscured AGN fraction ($\sim$30%) for $L_{5 \mu{\rm m}} > 10^{45}$ erg s$^{-1}$, while the fraction of Type-2 AGN with $L_{5 \mu{\rm m}} < 10^{45}$ erg s$^{-1}$ rises sharply from 40% to 80% of the overall AGN population.

Multi-wavelength surveys covering large sky volumes are necessary to obtain an accurate census of rare objects such as high luminosity and/or high redshift active galactic nuclei (AGN). Stripe 82X is a 31.3 deg$^2$ X-ray survey with $Chandra$ and $XMM$-Newton observations overlapping the legacy Sloan Digital Sky Survey (SDSS) Stripe 82 field, which has a rich investment of multi-wavelength coverage from the ultraviolet to the radio. The wide-area nature of this survey presents new challenges for photometric redshifts for AGN compared to previous work on narrow-deep fields because it probes different populations of objects that need to be identified and represented in the library of templates. Here we present an updated X-ray plus multi-wavelength matched catalog, including $Spitzer$ counterparts, and estimated photometric redshifts for 5961 (96% of a total of 6181) X-ray sources, which have a normalized median absolute deviation, $\sigma_{\rm nmad}$ = 0.06 and an outlier fraction, $\eta$ = 13.7%. The populations found in this survey, and the template libraries used for photometric redshifts, provide important guiding principles for upcoming large-area surveys such as $eROSITA$ and 3$XMM$ (in X-ray) and the Large Synoptic Survey Telescope (LSST; optical).

We present results of a ground-based near-infrared campaign with Palomar TripleSpec, Keck NIRSPEC, and Gemini GNIRS to target two samples of reddened active galactic nucleus (AGN) candidates from the 31 deg$^2$ Stripe 82 X-ray survey. One sample, which is $\sim$89\% complete to $K<16$ (Vega), consists of eight confirmed AGNs, four of which were identified with our follow-up program, and is selected to have red $R-K$ colors ($>4$, Vega). The fainter sample ($K>17$, Vega) represents a pilot program to follow-up four sources from a parent sample of 34 that are not detected in the single-epoch SDSS catalog and have {\it WISE} quasar colors. All twelve sources are broad-line AGNs (at least one permitted emission line has a FWHM exceeding 1300 km s$^{-1}$) and span a redshift range $0.59 < z < 2.5$. Half the ($R-K$)-selected AGNs have features in their spectra suggestive of outflows. When comparing these sources to a matched sample of blue Type 1 AGNs, we find the reddened AGNs are more distant ($z > 0.5$) and a greater percentage have high X-ray luminosities ($L_{\rm X,full} > 10^{44}$ erg s$^{-1}$). Such outflows and high luminosities may be consistent with the paradigm that reddened broad-line AGNs represent a transitory phase in AGN evolution as described by the major merger model for black hole growth. Results from our pilot program demonstrate proof-of-concept that our selection technique is successful in discovering reddened quasars at $z > 1$ missed by optical surveys.

We obtained $XMM-Newton$ observations of two highly luminous dust-reddened quasars, F2M1113+1244 and F2M1656+3821, that appear to be in the early, transitional phase predicted by merger-driven models of quasar/galaxy co-evolution. These sources have been well-studied at optical through mid-infrared wavelengths and are growing relatively rapidly, with Eddington ratios $>30\%$. Their black hole masses are relatively small compared to their host galaxies placing them below the $M_{\rm BH} - L_{\rm bulge}$ relation. We find that for both sources, an absorbed power-law model with $1-3\%$ of the intrinsic continuum scattered or leaked back into the line-of-sight best fits their X-ray spectra. We measure the absorbing column density ($N_H$) and constrain the dust-to-gas ratios in these systems, finding that they lie well below the Galactic value. This, combined with the presence of broad emission lines in their optical and near-infrared spectra, suggests that the dust absorption occurs far from the nucleus, in the host galaxy, while the X-rays are mostly absorbed in the nuclear, dust-free region within the sublimation radius. We also compare the quasars' absorption-corrected, rest-frame X-ray luminosities ($2-10$ keV) to their rest-frame infrared luminosities (6$\mu$m) and find that red quasars, similar to other populations of luminous obscured quasars, are either underluminous in X-rays or overluminous in the infrared.

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 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 days 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. (2016) 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 disk: superluminous supernovae, tidal disruption events, and mergers of stellar mass black holes.

CRTS J084133.15+200525.8 is an optically bright quasar at z=2.345 that has shown extreme spectral variability over the past decade. Photometrically, the source had a visual magnitude of V~17.3 between 2002 and 2008. Then, over the following five years, the source slowly brightened by approximately one magnitude, to V~16.2. Only ~1 in 10,000 quasars show such extreme variability, as quantified by the extreme parameters derived for this quasar assuming a damped random walk model. A combination of archival and newly acquired spectra reveal the source to be an iron low-ionization broad absorption line (FeLoBAL) quasar with extreme changes in its absorption spectrum. Some absorption features completely disappear over the 9 years of optical spectra, while other features remain essentially unchanged. We report the first definitive redshift for this source, based on the detection of broad H-alpha in a Keck/MOSFIRE spectrum. Absorption systems separated by several 1000 km/s in velocity show coordinated weakening in the depths of their troughs as the continuum flux increases. We interpret the broad absorption line variability to be due to changes in photoionization, rather than due to motion of material along our line of sight. This source highlights one sort of rare transition object that astronomy will now be finding through dedicated time-domain surveys.

We describe the first data release from the Spitzer-IRAC Equatorial Survey (SpIES); a large-area survey of 115 deg^2 in the Equatorial SDSS Stripe 82 field using Spitzer during its 'warm' mission phase. SpIES was designed to probe sufficient volume to perform measurements of quasar clustering and the luminosity function at z > 3 to test various models for "feedback" from active galactic nuclei (AGN). Additionally, the wide range of available multi-wavelength, multi-epoch ancillary data enables SpIES to identify both high-redshift (z > 5) quasars as well as obscured quasars missed by optical surveys. SpIES achieves 5{\sigma} depths of 6.13 {\mu}Jy (21.93 AB magnitude) and 5.75 {\mu}Jy (22.0 AB magnitude) at 3.6 and 4.5 microns, respectively - depths significantly fainter than WISE. We show that the SpIES survey recovers a much larger fraction of spectroscopically-confirmed quasars (98%) in Stripe 82 than are recovered by WISE (55%). This depth is especially powerful at high-redshift (z > 3.5), where SpIES recovers 94% of confirmed quasars, whereas WISE only recovers 25%. Here we define the SpIES survey parameters and describe the image processing, source extraction, and catalog production methods used to analyze the SpIES data. In addition to this survey paper, we release 234 images created by the SpIES team and three detection catalogs: a 3.6 {\mu}m-only detection catalog containing 6.1 million sources, a 4.5 {\mu}m-only detection catalog containing 6.5 million sources, and a dual-band detection catalog containing 5.4 million sources.

Some reddened quasars appear to be transitional objects in the merger-induced black hole growth/galaxy evolution paradigm, where a heavily obscured nucleus starts to be unveiled by powerful quasar winds evacuating the surrounding cocoon of dust and gas. Hard X-ray observations are able to peer through this gas and dust, revealing the properties of circumnuclear obscuration. Here, we present NuSTAR and XMM-Newton/Chandra observations of FIRST-2MASS selected red quasars F2M 0830+3759 and F2M 1227+3214. We find that though F2M 0830+3759 is moderately obscured ($N_{\rm H,Z} = 2.1\pm0.2 \times10^{22}$ cm$^{-2}$) and F2M 1227+3214 is mildly absorbed ($N_{\rm H,Z} = 3.4^{+0.8}_{-0.7}\times10^{21}$ cm$^{-2}$) along the line-of-sight, heavier global obscuration may be present in both sources, with $N_{\rm H,S} = 3.7^{+4.1}_{-2.6} \times 10^{23}$ cm$^{-2}$ and $< 5.5\times10^{23}$ cm$^{-2}$, for F2M 0830+3759 and F2M 1227+3214, respectively. F2M 0830+3759 also has an excess of soft X-ray emission below 1 keV which is well accommodated by a model where 7% of the intrinsic AGN X-ray emission is scattered into the line-of-sight. While F2M 1227+3214 has a dust-to-gas ratio ($E(B-V)$/$N_{\rm H}$) consistent with the Galactic value, the $E(B-V)$/$N_{\rm H}$ value for F2M 0830+3759 is lower than the Galactic standard, consistent with the paradigm that the dust resides on galactic scales while the X-ray reprocessing gas originates within the dust-sublimation zone of the broad-line-region. The X-ray and 6.1$\mu$m luminosities of these red quasars are consistent with the empirical relations derived for high-luminosity, unobscured quasars, extending the parameter space of obscured AGN previously observed by NuSTAR to higher luminosities.

Capitalizing on the all-sky coverage of {\it WISE}, and the 35\% and 50\% sky coverage from SDSS and Pan-STARRS, respectively, we explore the efficacy of $m_{R}$ (optical) - $m_{3.4 \mu m}$ (mid-infrared), hereafter $R-W1$, as a color diagnostic to identify obscured supermassive black hole accretion in wide-area X-ray surveys. We use the $\sim$16.5 deg$^2$ Stripe 82 X-ray survey data as a test-bed to compare $R-W1$ with $R-K$, an oft-used obscured AGN selection criterion, and examine where different classes of objects lie in this parameter space. Most stars follow a well-defined path in $R-K$ vs. $R-W1$ space. We demonstrate that optically normal galaxies hosting X-ray AGN at redshifts $0.5<z<1$ can be recovered with an $R-W1>4$ color-cut, while they typically are not selected as AGN based on their $W1-W2$ colors. Additionally, different observed X-ray luminosity bins favor different regions in $R-W1$ parameter space: moderate luminosity AGN ($10^{43}$ erg s$^{-1} < L_{\rm 0.5-10 keV} < 10^{44}$ erg s$^{-1}$) tend to have red colors while the highest luminosity AGN ($L_{\rm 0.5-10 keV} > 10^{45}$ erg s$^{-1}$) have bluer colors; higher spectroscopic completeness of the Stripe 82X sample is needed to determine whether this is a selection effect or an intrinsic property. Finally, we parameterize X-ray obscuration of Stripe 82X AGN by calculating their hardness ratios (HRs) and find no clear trends between HR and optical reddening. Our results will help inform best-effort practices in following-up obscured AGN candidates in current and future wide-area, shallow X-ray surveys, including the all-sky {\it eROSITA} mission.

We release the next installment of the Stripe 82 X-ray survey point-source catalog, which currently covers 31.3 deg$^2$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 Legacy field. In total, 6181 unique X-ray sources are significantly detected with {\it XMM-Newton} ($>5\sigma$) and {\it Chandra} ($>4.5\sigma$). This catalog release includes data from {\it XMM-Newton} cycle AO 13, which approximately doubled the Stripe 82X survey area. The flux limits of the Stripe 82X survey are $8.7\times10^{-16}$ erg s$^{-1}$ cm$^{-2}$, $4.7\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$, and $2.1\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the soft (0.5-2 keV), hard (2-10 keV), and full bands (0.5-10 keV), respectively, with approximate half-area survey flux limits of $5.4\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$, $2.9\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$, and $1.7\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$. We matched the X-ray source lists to available multi-wavelength catalogs, including updated matches to the previous release of the Stripe 82X survey; 88\% of the sample is matched to a multi-wavelength counterpart. Due to the wide area of Stripe 82X and rich ancillary multi-wavelength data, including coadded SDSS photometry, mid-infrared {\it WISE} coverage, near-infrared coverage from UKIDSS and VHS, ultraviolet coverage from {\it GALEX}, radio coverage from FIRST, and far-infrared coverage from {\it Herschel}, as well as existing $\sim$30\% optical spectroscopic completeness, we are beginning to uncover rare objects, such as obscured high-luminosity AGN at high-redshift. The Stripe 82X point source catalog is a valuable dataset for constraining how this population grows and evolves, as well as for studying how they interact with the galaxies in which they live.

Red quasars are thought to be an intermediate population between merger-driven star-forming galaxies in dust-enshrouded phase and normal quasars. If so, they are expected to have high accretion ratios, but their intrinsic dust extinction hampers reliable determination of Eddington ratios. Here, we compare the accretion rates of 16 red quasars at $z \sim 0.7$ to those of normal type 1 quasars at the same redshift range. The red quasars are selected by their red colors in optical through near-infrared (NIR) and radio detection. The accretion rates of the red quasars are derived from the P$\beta$ line in NIR spectra, which is obtained by the SpeX on the Infrared Telescope Facility (IRTF) in order to avoid the effects of dust extinction. We find that the measured Eddington ratios ($L_{\rm bol}$/$L_{\rm Edd} \simeq 0.69$) of red quasars are significantly higher than those of normal type 1 quasars, which is consistent with a scenario in which red quasars are the intermediate population and the black holes of red quasars grow very rapidly during such a stage.

Hierarchical assembly models predict a population of supermassive black hole (SMBH) binaries. These are not resolvable by direct imaging but may be detectable via periodic variability (or nanohertz frequency gravitational waves). Following our detection of a 5.2 year periodic signal in the quasar PG 1302-102 (Graham et al. 2015), we present a novel analysis of the optical variability of 243,500 known spectroscopically confirmed quasars using data from the Catalina Real-time Transient Survey (CRTS) to look for close (< 0.1 pc) SMBH systems. Looking for a strong Keplerian periodic signal with at least 1.5 cycles over a baseline of nine years, we find a sample of 111 candidate objects. This is in conservative agreement with theoretical predictions from models of binary SMBH populations. Simulated data sets, assuming stochastic variability, also produce no equivalent candidates implying a low likelihood of spurious detections. The periodicity seen is likely attributable to either jet precession, warped accretion disks or periodic accretion associated with a close SMBH binary system. We also consider how other SMBH binary candidates in the literature appear in CRTS data and show that none of these are equivalent to the identified objects. Finally, the distribution of objects found is consistent with that expected from a gravitational wave-driven population. This implies that circumbinary gas is present at small orbital radii and is being perturbed by the black holes. None of the sources is expected to merge within at least the next century. This study opens a new unique window to study a population of close SMBH binaries that must exist according to our current understanding of galaxy and SMBH evolution.

We used the Hubble Space Telescope WFC3 near-infrared camera to image the host galaxies of a sample of eleven luminous, dust-reddened quasars at z ~ 2 -- the peak epoch of black hole growth and star formation in the Universe -- to test the merger-driven picture for the co-evolution of galaxies and their nuclear black holes. The red quasars come from the FIRST+2MASS red quasar survey and a newer, deeper, UKIDSS+FIRST sample. These dust-reddened quasars are the most intrinsically luminous quasars in the Universe at all redshifts, and may represent the dust-clearing transitional phase in the merger-driven black hole growth scenario. Probing the host galaxies in rest-frame visible light, the HST images reveal that 8/10 of these quasars have actively merging hosts, while one source is reddened by an intervening lower redshift galaxy along the line-of-sight. We study the morphological properties of the quasar hosts using parametric Sersic fits as well as the non-parametric estimators (Gini coefficient, M_{20} and asymmetry). Their properties are heterogeneous but broadly consistent with the most extreme morphologies of local merging systems such as Ultraluminous Infrared galaxies. The red quasars have a luminosity range of log(L_bol) = 47.8 - 48.3 (erg/s) and the merger fraction of their AGN hosts is consistent with merger-driven models of luminous AGN activity at z=2, which supports the picture in which luminous quasars and galaxies co-evolve through major mergers that trigger both star formation and black hole growth.

Quasars have long been known to be variable sources at all wavelengths. Their optical variability is stochastic, can be due to a variety of physical mechanisms, and is well-described statistically in terms of a damped random walk model. The recent availability of large collections of astronomical time series of flux measurements (light curves) offers new data sets for a systematic exploration of quasar variability. Here we report on the detection of a strong, smooth periodic signal in the optical variability of the quasar PG 1302-102 with a mean observed period of 1,884 $\pm$ 88 days. It was identified in a search for periodic variability in a data set of light curves for 247,000 known, spectroscopically confirmed quasars with a temporal baseline of $\sim9$ years. While the interpretation of this phenomenon is still uncertain, the most plausible mechanisms involve a binary system of two supermassive black holes with a subparsec separation. Such systems are an expected consequence of galaxy mergers and can provide important constraints on models of galaxy formation and evolution.

We compare quasar selection techniques based on their optical variability using data from the Catalina Real-time Transient Survey (CRTS). We introduce a new technique based on Slepian wavelet variance (SWV) that shows comparable or better performance to structure functions and damped random walk models but with fewer assumptions. Combining these methods with WISE mid-IR colors produces a highly efficient quasar selection technique which we have validated spectroscopically. The SWV technique also identifies characteristic timescales in a time series and we find a characteristic rest frame timescale of ~54 days, confirmed in the light curves of ~18000 quasars from CRTS, SDSS and MACHO data, and anticorrelated with absolute magnitude. This indicates a transition between a damped random walk and $P(f) \propto f^{-1/3}$ behaviours and is the first strong indication that a damped random walk model may be too simplistic to describe optical quasar variability.

We have analyzed the {\it XMM-Newton} and {\it Chandra} data overlapping $\sim$16.5 deg$^2$ of Sloan Digital Sky Survey Stripe 82, including $\sim$4.6 deg$^2$ of proprietary {\it XMM-Newton} data that we present here. In total, 3362 unique X-ray sources are detected at high significance. We derive the {\it XMM-Newton} number counts and compare them with our previously reported {\it Chandra} Log$N$-Log$S$ relations and other X-ray surveys. The Stripe 82 X-ray source lists have been matched to multi-wavelength catalogs using a maximum likelihood estimator algorithm. We discovered the highest redshift ($z=5.86$) quasar yet identified in an X-ray survey. We find 2.5 times more high luminosity (L$_x \geq 10^{45}$ erg s$^{-1}$) AGN than the smaller area {\it Chandra} and {\it XMM-Newton} survey of COSMOS and 1.3 times as many identified by XBo\"otes. Comparing the high luminosity AGN we have identified with those predicted by population synthesis models, our results suggest that this AGN population is a more important component of cosmic black hole growth than previously appreciated. Approximately a third of the X-ray sources not detected in the optical are identified in the infrared, making them candidates for the elusive population of obscured high luminosity AGN in the early universe.

We present the results of a pilot survey to find dust-reddened quasars by matching the FIRST radio catalog to the UKIDSS near-infrared survey, and using optical data from SDSS to select objects with very red colors. The deep K-band limit provided by UKIDSS allows for finding more heavily-reddened quasars at higher redshifts as compared with previous work using FIRST and 2MASS. We selected 87 candidates with K<=17.0 from the UKIDSS Large Area Survey (LAS) First Data Release (DR1) which covers 190 deg2. These candidates reach up to ~1.5 magnitudes below the 2MASS limit and obey the color criteria developed to identify dust-reddened quasars. We have obtained 61 spectroscopic observations in the optical and/or near-infrared as well as classifications in the literature and have identified 14 reddened quasars with E(B-V)>0.1, including three at z>2. We study the infrared properties of the sample using photometry from the WISE Observatory and find that infrared colors improve the efficiency of red quasar selection, removing many contaminants in an infrared-to-optical color-selected sample alone. The highest-redshift quasars (z > 2) are only moderately reddened, with E(B-V) ~ 0.2-0.3. We find that the surface density of red quasars rises sharply with faintness, comprising up to 17% of blue quasars at the same apparent K-band flux limit. We estimate that to reach more heavily reddened quasars (i.e., E(B-V) > 0.5) at z>2 and a depth of K=17 we would need to survey at least ~2.5 times more area.

We present the discovery of a plausible super-luminous supernova (SLSN), found in the archival data of Sloan Digital Sky Survey (SDSS) Stripe 82, called PSN 000123+000504. The supernova peaked at M_g<-21.3 mag in the second half of September 2005, but was missed by the real-time supernova hunt. The observed part of the light curve (17 epochs) showed that the rise to the maximum took over 30 days, while the decline time lasted at least 70 days (observed frame), closely resembling other SLSNe of SN2007bi type. Spectrum of the host galaxy reveals a redshift of z=0.281 and the distance modulus of \mu=40.77 mag. Combining this information with the SDSS photometry, we found the host galaxy to be an LMC-like irregular dwarf galaxy with the absolute magnitude of M_B=-18.2+/-0.2 mag and the oxygen abundance of 12+log[O/H]=8.3+/-0.2. Our SLSN follows the relation for the most energetic/super-luminous SNe exploding in low-metallicity environments, but we found no clear evidence for SLSNe to explode in low-luminosity (dwarf) galaxies only. The available information on the PSN 000123+000504 light curve suggests the magnetar-powered model as a likely scenario of this event. This SLSN is a new addition to a quickly growing family of super-luminous SNe.

We present a measurement of the Type I quasar luminosity function at z=5 using a large sample of spectroscopically confirmed quasars selected from optical imaging data. We measure the bright end (M_1450<-26) with Sloan Digital Sky Survey (SDSS) data covering ~6000 deg^2, then extend to lower luminosities (M_1450<-24) with newly discovered, faint z~5 quasars selected from 235 deg^2 of deep, coadded imaging in the SDSS Stripe 82 region (the celestial equator in the Southern Galactic Cap). The faint sample includes 14 quasars with spectra obtained as ancillary science targets in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), and 59 quasars observed at the MMT and Magellan telescopes. We construct a well-defined sample of 4.7<z<5.1 quasars that is highly complete, with 73 spectroscopic identifications out of 92 candidates. Our color selection method is also highly efficient: of the 73 spectra obtained, 71 are high redshift quasars. These observations reach below the break in the luminosity function (M_1450* ~ -27). The bright end slope is steep (beta <~ -4), with a constraint of beta < -3.1 at 95% confidence. The break luminosity appears to evolve strongly at high redshift, providing an explanation for the flattening of the bright end slope reported previously. We find a factor of ~2 greater decrease in the number density of luminous quasars (M_1450<-26) from z=5 to z=6 than from z=4 to z=5, suggesting a more rapid decline in quasar activity at high redshift than found in previous surveys. Our model for the quasar luminosity function predicts that quasars generate ~30% of the ionizing photons required to keep the universe ionized at z=5.

We present the first results of a wide area X-ray survey within the Sloan Digital Sky Survey (SDSS) Stripe 82, a 300 deg$^2$ region of the sky with a substantial investment in multi-wavelength coverage. We analyzed archival {\it Chandra} observations that cover 7.5 deg$^2$ within Stripe 82 ("Stripe 82 ACX"), reaching 4.5$\sigma$ flux limits of 7.9$\times10^{-16}$, 3.4$\times10^{-15}$ and 1.8$\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the soft (0.5-2 keV), hard (2-7 keV) and full (0.5-7 keV) bands, to find 774, 239 and 1118 X-ray sources, respectively. Three hundred twenty-one sources are detected only in the full band and 9 sources are detected solely in the soft band. Utilizing data products from the {\it Chandra} Source Catalog, we construct independent Log$N$-Log$S$ relationships, detailing the number density of X-ray sources as a function of flux, which show general agreement with previous {\it Chandra} surveys. We compare the luminosity distribution of Stripe 82 ACX with the smaller, deeper CDF-S + E-CDFS surveys and with {\it Chandra}-COSMOS, illustrating the benefit of wide-area surveys in locating high luminosity AGN. We also investigate the differences and similarities of X-ray and optical selection to uncover obscured AGN in the local Universe. Finally, we estimate the population of AGN we expect to find with increased coverage of 100 deg$^2$ or 300 deg$^2$, which will provide unprecedented insight into the high redshift, high luminosity regime of black hole growth currently under-represented in X-ray surveys.

We compute black hole masses and bolometric luminosities for 57 active galactic nuclei (AGN) in the redshift range 1.25 < z < 2.67, selected from the GOODS-South deep multi-wavelength survey field via their X-ray emission. We determine host galaxy morphological parameters by separating the galaxies from their central point sources in deep HST images, and host stellar masses and colors by multi-wavelength SED fitting. 90% of GOODS AGN at these redshifts have detected rest-frame optical nuclear point sources; bolometric luminosities range from 2e43 - 2e46 erg/s. The black holes are growing at a range of accretion rates, with at least 50% of the sample having L/L_Edd < 0.1. 70% of host galaxies have stellar masses M* > 1e10 M_sun, with a range of colors suggesting a complex star formation history. We find no evolution of AGN bolometric luminosity within the sample, and no correlation between AGN bolometric luminosity and host stellar mass, color or morphology. Fully half the sample of host galaxies is disk-dominated, with another 25% having strong disk components. Fewer than 15% of the systems appear to be at some stage of a major merger. These moderate-luminosity AGN hosts are therefore inconsistent with a dynamical history dominated by mergers strong enough to destroy disks, indicating minor mergers or secular processes dominate the co-evolution of galaxies and their central black holes at z ~ 2.

We present mid-infrared spectra and photometry of thirteen redshift 0.4<z<1 dust-reddened quasars obtained with Spitzer IRS and MIPS. We compare properties derived from their infrared spectral energy distributions (intrinsic AGN luminosity and far-infrared luminosity from star formation) to the host luminosities and morphologies from HST imaging, and black hole masses estimated from optical and/or near-infrared spectroscopy. Our results are broadly consistent with models in which most dust reddened quasars are an intermediate phase between a merger-driven starburst triggering a completely obscured AGN, and a normal, unreddened quasar. We find that many of our objects have high accretion rates, close to the Eddington limit. These objects tend to fall below the black hole mass -- bulge luminosity relation as defined by local galaxies, whereas most of our low accretion rate objects are slightly above the local relation, as typical for normal quasars at these redshifts. Our observations are therefore most readily interpreted in a scenario in which galaxy stellar mass growth occurs first by about a factor of three in each merger/starburst event, followed sometime later by black hole growth by a similar amount. We do not, however, see any direct evidence for quasar feedback affecting star formation in our objects, for example in the form of a relationship between accretion rate and star formation. Five of our objects, however, do show evidence for outflows in the OIII 5007 Angstrom emission line profile, suggesting that the quasar activity is driving thermal winds in at least some members of our sample.

We present a sample of 120 dust-reddened quasars identified by matching radio sources detected at 1.4 GHz in the FIRST survey with the near-infrared 2MASS catalog and color-selecting red sources. Optical and/or near-infrared spectroscopy provide broad wavelength sampling of their spectral energy distributions that we use to determine their reddening, characterized by E(B-V). We demonstrate that the reddening in these quasars is best-described by SMC-like dust. This sample spans a wide range in redshift and reddening (0.1 < z < 3, 0.1 < E(B-V) < 1.5), which we use to investigate the possible correlation of luminosity with reddening. At every redshift, dust-reddened quasars are intrinsically the most luminous quasars. We interpret this result in the context of merger-driven quasar/galaxy co-evolution where these reddened quasars are revealing an emergent phase during which the heavily obscured quasar is shedding its cocoon of dust prior to becoming a "normal" blue quasar. When correcting for extinction, we find that, depending on how the parent population is defined, these red quasars make up < 15-20% of the luminous quasar population. We estimate, based on the fraction of objects in this phase, that its duration is 15-20% as long as the unobscured, blue quasar phase.

We explore the nature of heavily obscured quasar host galaxies at z~2 using deep Hubble Space Telescope WFC3/IR imaging of 28 Dust Obscured Galaxies (DOGs) to investigate the role of major mergers in driving black hole growth. The high levels of obscuration of the quasars selected for this study act as a natural coronagraph, blocking the quasar light and allowing a clear view of the underlying host galaxy. The sample of heavily obscured quasars represents a significant fraction of the cosmic mass accretion on supermassive black holes as the quasars have inferred bolometric luminosities around the break of the quasar luminosity function. We find that only a small fraction (4%, at most 11-25%) of the quasar host galaxies are major mergers. Fits to their surface brightness profiles indicate that 90% of the host galaxies are either disk dominated, or have a significant disk. This disk-like host morphology, and the corresponding weakness of bulges, is evidence against major mergers and suggests that secular processes are the predominant driver of massive black hole growth. Finally, we suggest that the co-incidence of mergers and AGN activity is luminosity dependent, with only the most luminous quasars being triggered mostly by major mergers.

One of the key open questions in cosmology today pertains to understanding when, where and how super massive black holes form, while it is clear that mergers likely play a significant role in the growth cycles of black holes, how supermassive black holes form, and how galaxies grow around them. Here, we present Hubble Space Telescope WFC3/IR grism observations of a clumpy galaxy at z=1.35, with evidence for 10^6 - 10^7 Msun rapidly growing black holes in separate sub-components of the host galaxy. These black holes could have been brought into close proximity as a consequence of a rare multiple galaxy merger or they could have formed in situ. Such holes would eventually merge into a central black hole as the stellar clumps/components presumably coalesce to form a galaxy bulge. If we are witnessing the in-situ formation of multiple black holes, their properties can inform seed formation models and raise the possibility that massive black holes can continue to emerge in star-forming galaxies as late as z=1.35 (4.8 Gyr after the Big Bang).

[Abridged] We investigate the average physical properties and star formation histories of the most UV-luminous star-forming galaxies at z~3.7. Our results are derived from analyses of the average spectral energy distributions (SEDs), constructed from stacked optical to infrared photometry, of a sample of the 1,902 most UV-luminous star-forming galaxies found in 5.3 square degrees of the NOAO Deep Wide-Field Survey. We bin the sample according to UV luminosity, and find that the shape of the average SED in the rest-frame optical and infrared is fairly constant with UV luminosity: i.e., more UV luminous galaxies are, on average, also more luminous at longer wavelengths. In the rest-UV, however, the spectral slope (measured at 0.13-0.28 um) rises steeply with the median UV luminosity from -1.8 at L L* to -1.2 in the brightest bin (L~4-5L*). We use population synthesis analyses to derive the average physical properties of these galaxies and find that: (1) L_UV, and thus star formation rates (SFRs), scale closely with stellar mass such that more UV-luminous galaxies are also more massive; (2) The median ages indicate that the stellar populations are relatively young (200-400 Myr) and show little correlation with UV luminosity; and (3) More UV-luminous galaxies are dustier than their less-luminous counterparts, such that L~4-5L* galaxies are extincted up to A(1600)=2 mag while L L* galaxies have A(1600)=0.7-1.5 mag. Based on these observations, we argue that the average star formation histories of UV-luminous galaxies are better described by models in which SFR increases with time in order to simultaneously reproduce the tight correlation between the observed SFR and stellar mass, and the universally young ages of these galaxies. We demonstrate the potential of measurements of the SFR-M* relation at multiple redshifts to discriminate between simple models of star formation histories.

We present an updated determination of the z ~ 4 QSO luminosity function (QLF), improving the quality of the determination of the faint end of the QLF presented in Glikman et al. (2010). We have observed an additional 43 candidates from our survey sample, yielding one additional QSO at z = 4.23 and increasing the completeness of our spectroscopic follow-up to 48% for candidates brighter than R = 24 over our survey area of 3.76 deg2. We study the effect of using K-corrections to compute the rest-frame absolute magnitude at 1450A compared with measuring M1450 directly from the object spectra. We find a luminosity-dependent bias: template-based K-corrections overestimate the luminosity of low-luminosity QSOs, likely due to their reliance on templates derived from higher luminosity QSOs. Combining our sample with bright quasars from the Sloan Digital Sky Survey and using spectrum-based M1450 for all the quasars, we fit a double-power-law to the binned QLF. Our best fit has a bright-end slope, {\alpha} = 3.3\pm0.2, and faint-end slope, {\beta} = 1.6(+0.8/-0.6). Our new data revise the faint-end slope of the QLF down to flatter values similar to those measured at z ~ 3. The break luminosity, though poorly constrained, is at M* = -24.1(+0.7/-1.9), approximately 1 - 1.5 mag fainter than at z ~ 3. This QLF implies that QSOs account for about half the radiation needed to ionize the IGM at these redshifts.

We have identified a sample of cool field brown dwarf candidates using IRAC data from the Spitzer Deep, Wide-Field Survey (SDWFS). The candidates were selected from 400,000 SDWFS sources with [4.5] <= 18.5 mag and required to have [3.6]-[4.5] >= 1.5 and [4.5] - [8.0] <= 2.0 on the Vega system. The first color requirement selects objects redder than all but a handful of presently known brown dwarfs with spectral classes later than T7, while the second eliminates 14 probable reddened AGN. Optical detection of 4 of the remaining 18 sources implies they are likely also AGN, leaving 14 brown dwarf candidates. For two of the brightest candidates (SDWFS J143524.44+335334.6 and SDWFS J143222.82+323746.5), the spectral energy distributions including near-infrared detections suggest a spectral class of ~ T8. The proper motion is < 0.25 "/yr, consistent with expectations for a luminosity inferred distance of >70 pc. The reddest brown dwarf candidate (SDWFS J143356.62+351849.2) has [3.6] - [4.5]=2.24 and H - [4.5] > 5.7, redder than any published brown dwarf in these colors, and may be the first example of the elusive Y-dwarf spectral class. Models from Burrows et al. (2003) predict larger numbers of cool brown dwarfs should be found for a Chabrier (2003) mass function. Suppressing the model [4.5] flux by a factor of two, as indicated by previous work, brings the Burrows models and observations into reasonable agreement. The recently launched Wide-field Infrared Survey Explorer (WISE) will probe a volume ~40x larger and should find hundreds of brown dwarfs cooler than T7.

We have conducted a spectroscopic survey to find faint quasars (-26.0 < M_{1450} < -22.0) at redshifts z=3.8-5.2 in order to measure the faint end of the quasar luminosity function at these early times. Using available optical imaging data from portions of the NOAO Deep Wide-Field Survey and the Deep Lens Survey, we have color-selected quasar candidates in a total area of 3.76 deg^2. Thirty candidates have R <= 23 mags. We conducted spectroscopic followup for 28 of our candidates and found 23 QSOs, 21 of which are reported here for the first time, in the 3.74 < z <5.06 redshift range. We estimate our survey completeness through detailed Monte Carlo simulations and derive the first measurement of the density of quasars in this magnitude and redshift interval. We find that the binned luminosity function is somewhat affected by the K-correction used to compute the rest-frame absolute magnitude at 1450A. Considering only our R <= 23 sample, the best-fit single power-law (Phi \propto L^beta) gives a faint-end slope beta = -1.6+/-0.2. If we consider our larger, but highly incomplete sample going one magnitude fainter, we measure a steeper faint-end slope -2 < beta < -2.5. In all cases, we consistently find faint-end slopes that are steeper than expected based on measurements at z ~ 3. We combine our sample with bright quasars from the Sloan Digital Sky Survey to derive parameters for a double-power-law luminosity function. Our best fit finds a bright-end slope, alpha = -2.4+/-0.2, and faint-end slope, beta = -2.3+/-0.2, without a well-constrained break luminosity. This is effectively a single power-law, with beta = -2.7+/-0.1. We use these results to place limits on the amount of ultraviolet radiation produced by quasars and find that quasars are able to ionize the intergalactic medium at these redshifts.

We identify 3,113 highly variable objects in 7,200 square degrees of the Palomar-QUEST Survey, which each varied by more than 0.4 magnitudes simultaneously in two broadband optical filters on timescales from hours to roughly 3.5 years. The primary goal of the selection is to find blazars by their well-known violent optical variability. Because most known blazars have been found in radio and/or X-ray wavelengths, a sample discovered through optical variability may have very different selection effects, elucidating the range of behavior possible in these systems. A set of blazars selected in this unusual manner will improve our understanding of the physics behind this extremely variable and diverse class of AGN. The object positions, variability statistics, and color information are available using the Palomar-QUEST CasJobs server. The time domain is just beginning to be explored over large sky areas; we do not know exactly what a violently variable sample will hold. About 20% of the sample has been classified in the literature; over 70% of those objects are known or likely AGN. The remainder largely consists of a variety of variable stars, including a number of RR Lyrae and cataclysmic variables.

The clustering of quasars on small scales yields fundamental constraints on models of quasar evolution and the buildup of supermassive black holes. This paper describes the first systematic survey to discover high redshift binary quasars. Using color-selection and photometric redshift techniques, we searched 8142 deg^2 of SDSS imaging data for binary quasar candidates, and confirmed them with follow-up spectroscopy. Our sample of 27 high redshift binaries (24 of them new discoveries) at redshifts 2.9 < z < 4.3 with proper transverse separations 10 kpc < R_{\perp} < 650 kpc increases the number of such objects known by an order of magnitude. Eight members of this sample are very close pairs with R_{\perp} < 100 kpc, and of these close systems four are at z > 3.5. The completeness and efficiency of our well-defined selection algorithm are quantified using simulated photometry and we find that our sample is ~ 50% complete. Our companion paper uses this knowledge to make the first measurement of the small scale clustering (R < 1 Mpc/h comoving) of high-redshift quasars. High redshift binaries constitute exponentially rare coincidences of two extreme (M >~ 10^9 Msun) supermassive black holes. At z ~ 4 there is about one close binary per 10 Gpc^3, thus these could be the highest sigma peaks, the analogs of superclusters, in the early Universe.

Digital synoptic sky surveys pose several new object classification challenges. In surveys where real-time detection and classification of transient events is a science driver, there is a need for an effective elimination of instrument-related artifacts which can masquerade as transient sources in the detection pipeline, e.g., unremoved large cosmic rays, saturation trails, reflections, crosstalk artifacts, etc. We have implemented such an Artifact Filter, using a supervised neural network, for the real-time processing pipeline in the Palomar-Quest (PQ) survey. After the training phase, for each object it takes as input a set of measured morphological parameters and returns the probability of it being a real object. Despite the relatively low number of training cases for many kinds of artifacts, the overall artifact classification rate is around 90%, with no genuine transients misclassified during our real-time scans. Another question is how to assign an optimal star-galaxy classification in a multi-pass survey, where seeing and other conditions change between different epochs, potentially producing inconsistent classifications for the same object. We have implemented a star/galaxy multipass classifier that makes use of external and a priori knowledge to find the optimal classification from the individually derived ones. Both these techniques can be applied to other, similar surveys and data sets.

We present results on a survey to find extremely dust-reddened Type-1 Quasars. Combining the FIRST radio survey, the 2MASS Infrared Survey and the Sloan Digital Sky Survey, we have selected a candidate list of 122 potential red quasars. With more than 80% spectroscopically identified objects, well over 50% are classified as dust-reddened Type 1 quasars, whose reddenings (E(B-V)) range from approximately 0.1 to 1.5 magnitudes. They lie well off the color selection windows usually used to detect quasars and many fall within the stellar locus, which would have made it impossible to find these objects with traditional color selection techniques. The reddenings found are much more consistent with obscuration happening in the host galaxy rather than stemming from the dust torus. We find an unusually high fraction of Broad Absorption Line (BAL) quasars at high redshift, all but one of them belonging to the Low Ionization BAL (LoBAL) class and many also showing absorption the metastable FeII line (FeLoBAL). The discovery of further examples of dust-reddened LoBAL quasars provides more support for the hypothesis that BAL quasars (at least LoBAL quasars) represent an early stage in the lifetime of the quasar. The fact that we see such a high fraction of BALs could indicate that the quasar is in a young phase in which quasar feedback from the BAL winds is suppressing star formation in the host galaxy.

We conducted an exploratory search for quasars at z~ 6 - 8, using the Early Data Release from United Kingdom Infrared Deep Sky survey (UKIDSS) cross-matched to panoramic optical imagery. High redshift quasar candidates are chosen using multi-color selection in i,z,Y,J,H and K bands. After removal of apparent instrumental artifacts, our candidate list consisted of 34 objects. We further refined this list with deeper imaging in the optical for ten of our candidates. Twenty-five candidates were followed up spectroscopically in the near-infrared and in the optical. We confirmed twenty-five of our spectra as very low-mass main-sequence stars or brown dwarfs, which were indeed expected as the main contaminants of this exploratory search. The lack of quasar detection is not surprising: the estimated probability of finding a single z>6 quasar down to the limit of UKIDSS in the 27.3 square degrees of the EDR is <5%. We find that the most important limiting factor in this work is the depth of the available optical data. Experience gained in this pilot project can help refine high-redshift quasar selection criteria for subsequent UKIDSS data releases.

There is an increasing number of large, digital, synoptic sky surveys, in which repeated observations are obtained over large areas of the sky in multiple epochs. Likewise, there is a growth in the number of (often automated or robotic) follow-up facilities with varied capabilities in terms of instruments, depth, cadence, wavelengths, etc., most of which are geared toward some specific astrophysical phenomenon. As the number of detected transient events grows, an automated, probabilistic classification of the detected variables and transients becomes increasingly important, so that an optimal use can be made of follow-up facilities, without unnecessary duplication of effort. We describe a methodology now under development for a prototype event classification system; it involves Bayesian and Machine Learning classifiers, automated incorporation of feedback from follow-up observations, and discriminated or directed follow-up requests. This type of methodology may be essential for the massive synoptic sky surveys in the future.

We describe briefly the Palomar-Quest (PQ) digital synoptic sky survey, including its parameters, data processing, status, and plans. Exploration of the time domain is now the central scientific and technological focus of the survey. To this end, we have developed a real-time pipeline for detection of transient sources. We describe some of the early results, and lessons learned which may be useful for other, similar projects, and time-domain astronomy in general. Finally, we discuss some issues and challenges posed by the real-time analysis and scientific exploitation of massive data streams from modern synoptic sky surveys.

Combining radio observations with optical and infrared color selection -- demonstrated in our pilot study to be an efficient selection algorithm for finding red quasars -- we have obtained optical and infrared spectroscopy for 120 objects in a complete sample of 156 candidates from a sky area of 2716 square degrees. Consistent with our initial results, we find our selection criteria -- J-K>1.7, R-K>4.0 -- yield a ~50% success rate for discovering quasars substantially redder than those found in optical surveys. Comparison with UVX- and optical color-selected samples shows that >~ 10% of the quasars are missed in a magnitude-limited survey. Simultaneous two-frequency radio observations for part of the sample indicate that a synchrotron continuum component is ruled out as a significant contributor to reddening the quasars' spectra. We go on to estimate extinctions for our objects assuming their red colors are caused by dust. Continuum fits and Balmer decrements suggest E(B-V) values ranging from near zero to 2.5 magnitudes. Correcting the K-band magnitudes for these extinctions, we find that for K <= 14.0, red quasars make up between 25% and 60% of the underlying quasar population; owing to the incompleteness of the 2MASS survey at fainter K-band magnitudes, we can only set a lower limit to the radio-detected red quasar population of >20-30%.

We report the discovery of two low-luminosity quasars at z~4, both of which show prominent N IV] 1486A emission. This line is extremely rare in quasar spectra at any redshift; detecting it in two out of a sample of 23 objects (i.e., ~ 9% of the sample) is intriguing and is likely due to the low-luminosity, high-redshift quasar sample we are studying. This is still a poorly explored regime, where contributions from associated, early starbursts may be significant. One interpretation of this line posits photoionization by very massive young stars. Seeing N IV] 1486A emission in a high-redshift quasar may thus be understood in the context of co-formation and early co-evolution of galaxies and their supermassive black holes. Alternatively, we may be seeing a phenomenon related to the early evolution of quasar broad emission line regions. The non-detection (and possibly even broad absorption) of N V 1240A line in the spectrum of one of these quasars may support that interpretation. These two objects may signal a new faint quasar population or an early AGN evolutionary stage at high redshifts.

We report the discovery of the first known probable case of a physical triple quasar (not a gravitational lens). A previously known double system, QQ 1429-008 at z = 2.076, is shown to contain a third, fainter QSO component at the same redshift within the measurement errors. Deep optical and IR imaging at the Keck and VLT telescopes has failed to reveal a plausible lensing galaxy group or a cluster, and moreover, we are unable to construct any viable lensing model which could lead to the observed distribution of source positions and relative intensities of the three QSO image components. Furthermore, there are hints of differences in broad-band spectral energy distributions of different components, which are more naturally understood if they are physically distinct AGN. Therefore, we conclude that this system is most likely a physical triple quasar, the first such close QSO grouping known at any redshift. The projected component separations in the restframe are ~ 30 - 50 kpc for the standard concordance cosmology, typical of interacting galaxy systems. The existence of this highly unusual system supports the standard picture in which galaxy interactions lead to the onset of QSO activity.

We present a technique to explore the radio sky into the nanoJansky regime by employing image stacking using the FIRST survey. We first discuss the non-intuitive relationship between the mean and median values of a distribution that is dominated by noise, followed by an analysis of the systematic effects present in FIRST's 20cm VLA snapshot images. Image stacking allows us to recover the properties of source populations with fluxes a factor of 30 or more below the rms noise level. Mean estimates of radio flux density, luminosity, etc., are derivable for any source class having arcsecond positional accuracy. We use this technique to compute the mean radio properties for 41,295 quasars from the SDSS DR3 catalog. There is a tight correlation between optical and radio luminosity, with the radio luminosity increasing as the 0.85 power of optical luminosity. This implies declining radio-loudness with optical luminosity: the most luminous objects (M=-28.5) have average radio-to-optical ratios 3 times lower than the least luminous objects (M=-20). There is also a striking correlation between optical color and radio loudness: quasars that are either redder or bluer than the norm are brighter radio sources, with objects 0.8 magnitudes redder than the SDSS composite spectrum having radio-loudness ratios that are higher by a factor of 10. We explore the longstanding question of whether a radio-loud/radio-quiet dichotomy exists in quasars, finding that optical selection effects probably dominate the distribution function of radio loudness, which has at most a modest (~20%) inflection between the radio-loud and radio-quiet ends of the distribution. We also find, surprisingly, that broad absorption line quasars have higher mean radio flux densities, with the greatest disparity arising in the rare low-ionization BAL subclass.

We review some of the recent developments and challenges posed by the data analysis in modern digital sky surveys, which are representative of the information-rich astronomy in the context of Virtual Observatory. Illustrative examples include the problems of an automated star-galaxy classification in complex and heterogeneous panoramic imaging data sets, and an automated, iterative, dynamical classification of transient events detected in synoptic sky surveys. These problems offer good opportunities for productive collaborations between astronomers and applied computer scientists and statisticians, and are representative of the kind of challenges now present in all data-intensive fields. We discuss briefly some emergent types of scalable scientific data analysis systems with a broad applicability.

We present a near-infrared quasar composite spectrum spanning the wavelength range 0.58 - 3.5 um. The spectrum has been constructed from observations of 27 quasars obtained at the NASA IRTF telescope and satisfying the criteria Ks < 14.5 and Mi < -23; the redshift range is 0.118 < z < 0.418. The signal-to-noise is moderate, reaching a maximum of 150 between 1.6 and 1.9 um. While a power-law fit to the continuum of the composite spectrum requires two breaks, a single power-law slope of alpha=-0.92 plus a 1260 K blackbody provides an excellent description of the spectrum from H-alpha to 3.5 um, strongly suggesting the presence of significant quantities of hot dust in this blue-selected quasar sample. We measure intensities and line widths for ten lines, finding that the Paschen line ratios rule out Case B recombination. We compute K-corrections for the J, H, K, and Spitzer 3.6 um bands which will be useful in analyzing observations of quasars up to z=10.

We have constructed a sample of bright near-infrared sources which are detected at radio wavelengths but undetected on the POSS I plates in order to search for a population of dust-obscured quasars. Optical and infrared spectroscopic followup of the sample has led to the discovery of seventeen heavily reddened quasars (B-K>6.5), fourteen of which are reported here for the first time. This has allowed us to define a region in the R-K, J-K color plane in which 50% of the radio-selected objects are highly reddened quasars. We compare the surface density of this previously overlooked population to that of UVX- radio-selected quasars, finding that they make up ~20% of the total quasar population for K < 15.5.

We report the discovery of a z=2.65 low-ionization iron broad absorption line quasar, FIRST J100424.9+122922, which is gravitationally-lensed by a galaxy at z~0.95. The object was discovered as part of a program to find very red quasars by matching the FIRST radio survey with the 2-MASS near-infrared survey. J100424.9+122922 is the second lensed system to be found in this program, suggesting that many gravitational lenses are probably missed from conventional optical quasar surveys. We have made a simple lens model and a rough estimate of the reddening in the immediate environment of the quasar which suggests that the quasar is intrinsically very luminous and is accreting at close to the Eddington limit of its ~10^9 M_sun black hole. The lensing galaxy has a small amount of dust which is responsible for some excess reddening observed in the fainter image of the quasar, but is otherwise a fairly typical massive elliptical galaxy. We model the selection effects working against the detection of red quasars in both lensed and unlensed samples. We show that these selection effects are very effective at removing even lightly-reddened high redshift quasars from magnitude-limited samples, whether they are lensed or not. This suggests that the red quasar population in general could be very large, and in particular the class of iron broad absorption line quasars of which J100424.9+122922 is a member may be much larger than their rarity in magnitude-limited samples would suggest.

In a survey of quasar candidates selected by matching the FIRST and 2MASS catalogs, we have found two extraordinarily red quasars. FIRST J013435.7-093102 is a 1 Jy source at z=2.216 and has B-K > 10, while FIRST J073820.1+275045 is a 2.5 mJy source at z=1.985 with B-K = 8.4. FIRST J073820.1+275045 has strong absorption lines of MgII and CIV in the rest frame of the quasar and is highly polarized in the rest frame ultraviolet, strongly favoring the interpretation that its red spectral energy distribution is caused by dust reddening local to the quasar. FIRST J073820.1+275045 is thus one of the few low radio-luminosity, highly dust-reddened quasars known. The available observational evidence for FIRST J013435.7-093102 leads us to conclude that it too is reddened by dust. We show that FIRST J013435.7-093102 is gravitationally lensed, increasing the number of known lensed, extremely dust-reddened quasars to at least three, including MG0414-0534 and PKS1830-211. We discuss the implications of whether these objects are reddened by dust in the host or lensing galaxies. If reddened by their local environment, then we estimate that between 10 and 20% of the radio-loud quasar population is reddened by dust in the host galaxy. The discovery of FIRST J073820.1+275045 and objects now emerging from X-ray surveys suggests the existence of an analogous radio-quiet red quasar population. Such objects will be entirely missed by standard radio or optical quasar surveys. If dust in the lensing galaxies is primarily responsible for the extreme redness of the lensed quasars, then an untold number of gravitationally lensed quasars are being overlooked.