Search

We present the pc-scale radio spectra of a representative sample of 13 Palomar-Green radio-quiet quasars, based on our new Very Long Baseline Array (VLBA) observations at 8.4 and 23.6 GHz and our earlier VLBA studies at 1.5 and 5.0 GHz. The radio core emission generally exhibits a flat spectrum at 1.5-5.0 GHz, which indicates a compact optically thick synchrotron source on a scale smaller than the broad-line region (BLR) radius R_BLR ~ 0.01-0.1 pc. The 8.4-23.6 GHz spectral slope remains flat in four objects indicating the inner radius of the radio source R_in < 0.1 R_BLR, and becomes steep in four other objects indicating R_in ~ 0.5 R_BLR. The flat 8.4-23.6 GHz slope sources may be associated with a continuous ejection starting at the accretion disk corona. The steep 8.4-23.6 GHz slope sources may be produced by an interaction of an AGN-driven wind with the BLR gas or a low-power jet extending to the BLR scale. Seven of these eight objects, which have a flat or steep 8.4-23.6 GHz slope, reside at the Eddington ratios L/L_Edd < 0.3, and four of the remaining five objects, where the 8.4-23.6 GHz fluxes are too faint to significantly constrain the slope, reside at L/L_Edd > 0.3. The 8.4-23.6 GHz radio emission in the high L/L_Edd objects may be weak due to more extended emission from a radiation pressure driven wind. Future sub-millimeter observations can further constrain the inward radial extent of the radio emission down to the coronal scale., Comment: Accepted for publication in ApJ

Does a broad-line region (BLR) wind in radio-quiet (RQ) active galactic nuclei (AGN) extend to pc scales and produce radio emission? We explore the correlations between a pc-scale radio wind and the BLR wind in a sample of 19 RQ Palomar-Green (PG) quasars. The radio wind is defined based on the spectral slope and the compactness of the emission at 1.5-5 GHz, and the BLR wind is defined by the excess blue wing in the C IV emission line profile. The five objects with both radio and BLR wind indicators are found at high Eddington ratios L/L_Edd (> 0.66), and eight of the nine objects with neither radio nor BLR winds reside at low L/L_Edd (< 0.28). This suggests that the BLR wind and the radio wind in RQ AGN are related to a radiation pressure driven wind. Evidence for free-free absorption by AGN photoionized gas, which flattens the spectral slope, is found in two objects. Radio outflows in three low L/L_Edd (0.05-0.12) objects are likely from a low-power jet, as suggested by additional evidence. The presence of a mild equatorial BLR wind in four intermediate L/L_Edd (0.2-0.4) objects can be tested with future spectropolarimetry., Comment: Accepted for publication in ApJ

The supermassive black holes ($M_{\rm BH} \sim 10^{6}$$-$$10^{10}~M_\odot$) that power luminous active galactic nuclei (AGNs), i.e., quasars, generally show a correlation between thermal disk emission in the ultraviolet (UV) and coronal emission in hard X-rays. In contrast, some "massive" black holes (mBHs; $M_{\rm BH} \sim 10^{5}$$-$$10^{6}~M_\odot$) in low-mass galaxies present curious X-ray properties with coronal radiative output up to 100$\times$ weaker than expected. To examine this issue, we present a pilot study incorporating Very Large Array radio observations of a sample of 18 high-accretion-rate (Eddington ratios $L_{\rm bol}/L_{\rm Edd} > 0.1$), mBH-powered AGNs ($M_{\rm BH} \sim 10^{6}~M_\odot$) with Chandra X-ray coverage. Empirical correlations previously revealed in samples of radio-quiet, high-Eddington AGNs indicate that the radio$-$X-ray luminosity ratio, $L_{\rm R}/L_{\rm X}$, is approximately constant. Through multiwavelength analysis, we instead find that the X-ray-weaker mBHs in our sample tend toward larger values of $L_{\rm R}/L_{\rm X}$ even though they remain radio-quiet per their optical$-$UV properties. This trend results in a tentative but highly intriguing correlation between $L_{\rm R}/L_{\rm X}$ and X-ray weakness, which we argue is consistent with a scenario in which X-rays may be preferentially obscured from our line of sight by a "slim" accretion disk. We compare this observation to weak emission-line quasars (AGNs with exceptionally weak broad-line emission and a significant X-ray-weak fraction) and conclude by suggesting that our results may offer a new observational signature for finding high-accretion-rate AGNs., Comment: Accepted for publication in ApJ. 26 pages, 8 figures, 6 tables

We present analysis of a homogeneous, optically selected, volume-limited ($0.2

We present radio spectra spanning $0.1 - 10$ GHz for the sample of heavily obscured luminous quasars with extremely red mid-infrared-optical colors and compact radio emission. The spectra are constructed from targeted 10 GHz observations and archival radio survey data, which together yield $6-11$ flux density measurements for each object. Our suite of Python tools for modeling the radio spectra is publicly available on Github. Our primary result is that most (61%) of the sample have peaked or curved radio spectra and many (36%) could be classified as Gigahertz Peaked Spectrum (GPS) sources. This indicates compact emission regions likely arising from recently triggered radio jets. Assuming synchrotron self-absorption (SSA) generates the peaks, we infer compact source sizes ($3 - 100$ pc) with strong magnetic fields ($6 - 100$ mG) and young ages ($30 - 10^4$ years). Conversely, free-free absorption (FFA) could also create peaks due to the high column densities associated with the deeply embedded nature of the sample. However, we find no correlations between the existence or frequency of the peaks and any parameters of the MIR emission. The high-frequency spectral indices are steep ($\alpha \approx -1$) and correlate, weakly, with the ratio of MIR photon energy density to magnetic energy density, suggesting that the spectral steepening could arise from inverse Compton scattering off the intense MIR photon field. This study provides a foundation for combining multi-frequency and mixed-resolution radio survey data for understanding the impact of young radio jets on the ISM and star formation rates of their host galaxies., Comment: 48 pages, 17 figures, published in Astrophysical Journal

We have performed a search over 3440 deg$^2$ of Epoch 1 (2017-2019) of the Very Large Array Sky Survey to identify unobscured quasars in the optical ($0.2 < z < 3.2$) and obscured active galactic nuclei (AGN) in the infrared that have brightened dramatically in the radio over the past one to two decades. These sources would have been previously classified as "radio-quiet" quasars based on upper limits from the Faint Images of the Radio Sky at Twenty Centimeters survey (1993-2011), but they are now consistent with "radio-loud" quasars ($L_{\rm 3\,GHz} = 10^{40 - 42} \,\, {\rm erg} \,{\rm s}^{-1}$). A quasi-simultaneous, multiband ($\sim1-18$ GHz) follow-up study of 14 sources with the VLA has revealed compact sources ($<0.1^{\prime \prime}$ or $<1$ kpc) with peaked radio spectral shapes. The high-amplitude variability over decadal timescales at 1.5 GHz (100% to $>$2500%), but roughly steady fluxes over a few months at 3 GHz, are inconsistent with extrinsic variability due to propagation effects, thus favoring an intrinsic origin. We conclude that our sources are powerful quasars hosting compact/young jets. This challenges the generally accepted idea that "radio-loudness" is a property of the quasar/AGN population that remains fixed on human timescales. Our study suggests that frequent episodes of short-lived AGN jets that do not necessarily grow to large scales may be common at high redshift. We speculate that intermittent but powerful jets on subgalactic scales could interact with the interstellar medium, possibly driving feedback capable of influencing galaxy evolution., Comment: 28 pages, 9 figures, accepted to ApJ

We present new sub-arcsecond-resolution Karl G. Jansky Very Large Array (VLA) imaging at 10 GHz of 155 ultra-luminous ($L_{\rm bol}\sim10^{11.7-14.2} L_\odot$) and heavily obscured quasars with redshifts $z \sim0.4-3$. The sample was selected to have extremely red mid-infrared (MIR)-optical color ratios based on data from Wide-Field Infrared Survey Explorer (WISE) along with a detection of bright, unresolved radio emission from the NRAO VLA Sky Survey (NVSS) or Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) Survey. Our high-resolution VLA observations have revealed that the majority of the sources in our sample (93 out of 155) are compact on angular scales $<0.2^{\prime \prime}$ ($\leq 1.7$ kpc at $z \sim2$). The radio luminosities, linear extents, and lobe pressures of our sources are similar to young radio active galactic nuclei (AGN; e.g., Gigahertz Peaked Spectrum, GPS, and Compact Steep Spectrum, CSS, sources), but their space density is considerably lower. Application of a simple adiabatic lobe expansion model suggests relatively young dynamical ages ($\sim10^{4-7}$ years), relatively high ambient ISM densities ($\sim1-10^4$ cm$^{-3}$), and modest lobe expansion speeds ($\sim30-10,000$ km s$^{-1}$). Thus, we find our sources to be consistent with a population of newly triggered, young jets caught in a unique evolutionary stage in which they still reside within the dense gas reservoirs of their hosts. Based on their radio luminosity function and dynamical ages, we estimate only $\sim20\%$ of classical large scale FRI/II radio galaxies could have evolved directly from these objects. We speculate that the WISE-NVSS sources might first become GPS or CSS sources, of which some might ultimately evolve into larger radio galaxies., Comment: Accepted for publication in ApJ, 38 pages, 14 figures, 6 tables

Most of the stars today reside in galactic spheroids, whose properties are tightly tied to the supermassive black holes (MBHs) at their centers, implying that the accretion activity onto MBHs leaves a lasting imprint on the evolution of their host galaxies. Despite the importance of this so-called MBH-galaxy co-evolution, the physical mechanisms responsible for driving this relationship - such as the dominant mode of energetic feedback from active galactic nuclei (AGN) - remain a poorly understand aspect of galaxy assembly. A key challenge for identifying and characterizing AGN during the peak epoch of galaxy assembly and beyond is the presence of large columns of gas and dust, which fuels the growth of their MBHs but effectively obscures them from view in optical and X-ray studies. The high sensitivity of the ngVLA will capture emission from AGN in an extinction-free manner out to $z \sim 6$ and beyond. At lower-redshifts ($z \sim 2$), the high angular resolution of the ngVLA will enable spatially-resolved studies capable of localizing the sites of actively growing MBHs within their host galaxies during the peak epoch of cosmic assembly., Comment: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)

We have detected in ALMA observations CO J = 6 - 5 emission from the nucleus of the Seyfert galaxy NGC 1068. The low-velocity (up to +/- 70 km/s relative to systemic) CO emission resolves into a 12x7 pc structure, roughly aligned with the nuclear radio source. Higher-velocity emission (up to +/- 400 km/s) is consistent with a bipolar outflow in a direction nearly perpendicular (roughly 80 degrees) to the nuclear disk. The position-velocity diagram shows that in addition to the outflow, the velocity field may also contain rotation about the disk axis. These observations provide compelling evidence in support of the disk-wind scenario for the AGN obscuring torus., Comment: Accepted for publication in the Astrophysical Journal Letters. 6 pages, 3 figures

We present ALMA 870 micron (345 GHz) data for 49 high redshift (0.47

We present observations and analysis of an unusual [C II] emission line in the very luminous QSO SDSS J155426.16+193703.0 at z~4.6. The line is extremely broad (FWHM 735 km/s) and seems to have a flat-topped or double-peaked line profile. A velocity map of the line shows a gradient across the source that indicates large-scale rotation of star-forming gas. Together, the velocity map and line profile suggest the presence of a massive rotating disc with a dynamical mass M_dyn > 5x10^10 M_sun. Using the assumption of a rotating disc origin, we employ an empirical relation between galaxy disc circular velocity and bulge velocity dispersion (sigma) to estimate that sigma > 310 km/s, subject to a correction for the unknown disc inclination. This result implies that this source is consistent with the local M--sigma relation, or offset at most by an order of magnitude in black hole mass. In contrast, the assumption of a bulge origin for the [C II] emission line would lead to a conclusion that the black hole is nearly two orders of magnitude more massive than predicted by the M--sigma relation, similar to previous findings for other high-redshift QSOs. As disc rotation may be a common origin for [C II] emission at high redshifts, these results stress that careful consideration of dynamical origins is required when using observations of this line to derive properties of high-redshift galaxies., Comment: 12 pages, 5 figures, accepted for publications in MNRAS

We present a sample of 1317 Infrared-Faint Radio Sources (IFRSs) that, for the first time, are reliably detected in the infrared, generated by cross-correlating the Wide-Field Infrared Survey Explorer (WISE) all-sky survey with major radio surveys. Our IFRSs are brighter in both radio and infrared than the first generation IFRSs that were undetected in the infrared by the Spitzer Space Telescope. We present the first spectroscopic redshifts of IFRSs, and find that all but one of the IFRSs with spectroscopy has z > 2. We also report the first X-ray counterparts of IFRSs, and present an analysis of radio spectra and polarization, and show that they include Gigahertz-Peaked Spectrum, Compact Steep Spectrum, and Ultra-Steep Spectrum sources. These results, together with their WISE infrared colours and radio morphologies, imply that our sample of IFRSs represents a population of radio-loud Active Galactic Nuclei at z > 2. We conclude that our sample consists of lower-redshift counterparts of the extreme first generation IFRSs, suggesting that the fainter IFRSs are at even higher redshift., Comment: 23 pages, 17 figures. Submitted to MNRAS

We present near-infrared spectra of young radio quasars [P(1.4GHz) ~ 26-27 W/Hz] selected from the Wide-Field Infrared Survey Explorer. The detected objects have typical redshifts of z ~ 1.6-2.5 and bolometric luminosities ~ 10^47 erg/s. Based on the intensity ratios of narrow emission lines, we find that these objects are mainly powered by active galactic nuclei (AGNs), although star formation contribution cannot be completely ruled out. The host galaxies experience moderate levels of extinction, A(V) ~ 0-1.3 mag. The observed [O III] luminosities and rest-frame J-band magnitudes constrain the black hole masses to lie in the range ~ 10^8.9-10^9.7 solar mass. From the empirical correlation between black hole mass and host galaxy mass, we infer stellar masses of ~ 10^11.3-10^12.2 solar mass. The [O III] line is exceptionally broad, with full width at half maximum ~1300 to 2100 km/s, significantly larger than that of ordinary distant quasars. We argue that these large line widths can be explained by jet-induced outflows, as predicted by theoretical models of AGN feedback., Comment: To appear in The Astrophysical Journal Letters

We used the 1.4 GHz NVSS to study radio sources in two color-selected QSO samples: a volume-limited sample of 1313 QSOs defined by M_i < -23 in the redshift range 0.2 < z < 0.45 and a magnitude-limited sample of 2471 QSOs with m_r < 18.5 and 1.8 < z < 2.5. About 10% were detected above the 2.4 mJy NVSS catalog limit and are powered primarily by AGNs. The space density of the low-redshift QSOs evolves as rho proportional to (1+z)^6. In both redshift ranges the flux-density distributions and luminosity functions of QSOs stronger than 2.4 mJy are power laws, with no features to suggest more than one kind of radio source. Extrapolating the power laws to lower luminosities predicts the remaining QSOs should be extremely radio quiet, but they are not. Most were detected statistically on the NVSS images with median peak flux densities S_p(mJy/beam) ~ 0.3 and 0.05 in the low- and high-redshift samples, corresponding to 1.4 GHz spectral luminosities log[L(W/Hz)] ~ 22.7$ and 24.1, respectively. We suggest that the faint radio sources are powered by star formation at rates ~20 M_sun per year in the moderate luminosity (median M_i ~ -23.4) low-redshift QSOs and ~500 M_sun per year in the very luminous (M_i} ~ -27.5) high-redshift QSOs. Such luminous starbursts [ log(L / L_sun) ~ 11.2 and 12.6, respectively] are consistent with "quasar mode" accretion in which cold gas flows fuel both AGN and starburst., Comment: 17 pages, 12 figures Accepted for publication in the ApJ

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $\Lambda$CDM cosmology with H0 = 70 km/s/Mpc, $\Omega_{\rm M}$ = 0.3, and $\Omega_{\Lambda}$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys., Comment: Accepted for publication in A&A, Catalog available online at http://www.sdss3.org/dr9/algorithms/qso_catalog.php

Despite decades of study, it remains unclear whether there are distinct radio-loud and radio-quiet populations of quasi-stellar objects (QSOs). Early studies were limited by inhomogeneous QSO samples, inadequate sensitivity to probe the radio-quiet population, and degeneracy between redshift and luminosity for flux-density-limited samples. Our new 6 GHz EVLA observations allow us for the first time to obtain nearly complete (97%) radio detections in a volume-limited color-selected sample of 179 QSOs more luminous than M_i = -23 from the Sloan Digital Sky Survey (SDSS) Data Release Seven in the narrow redshift range 0.2 < z < 0.3. The dramatic improvement in radio continuum sensitivity made possible with the new EVLA allows us, in 35 minutes of integration, to detect sources as faint as 20 microJy, or log[L_6 (W/Hz)] ~ 21.5 at z = 0.25, well below the radio luminosity, log[L_6 (W/Hz)] ~ 22.5, that separates star-forming galaxies from radio-loud active galactic nuclei (AGNs) driven by accretion onto a super-massive black hole. We calculate the radio luminosity function (RLF) for these QSOs using three constraints: (a) EVLA 6 GHz observations for log[L_6 (W/Hz)] < 23.5, (b) NRAO-VLA Sky Survey (NVSS) observations for log[L_6 (W/Hz)] > 23.5, and (c) the total number of SDSS QSOs in our volume-limited sample. We show that the RLF can be explained as a superposition of two populations, dominated by AGNs at the bright end and star formation in the QSO host galaxies at the faint end., Comment: 11 pages, 2 figures. Accepted for publication. This Letter will appear in the ApJL EVLA special issue

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU., Comment: Accepted by PASA; see http://www.atnf.csiro.au/people/rnorris/papers/n252.pdf for higher-resolution figures

(Abridged) With the largest homogeneous quasar sample with optical spectra and robust radio morphology classifications assembled to date, we investigate quasar radio and optical properties with unprecedented statistical power. The data consist of 4714 radio quasars from FIRST with S_{20}>2mJy and SDSS spectra. Radio morphology classes include core-only (core), core-lobe (lobe), core-jet (jet), lobe-core-lobe (triple), and double-lobe. We examine the optical colors of radio-morphology subsamples and find that radio quasars with core emission unresolved by FIRST (on 5" scale) have a redder color distribution than radio-quiet quasars (S_{20}<1mJy); other classes of radio quasars have optical color distributions similar to the radio-quiet quasars. This analysis also suggests that optical colors of z<2.7 SDSS quasars are not strongly (<0.1mag) biased blue. We show that the radio core-to-lobe flux density ratio (R) and the radio-to-optical (i-band) ratio of the quasar core (RI) are correlated, suggesting that both parameters are indicative of line-of-sight orientation. We investigate spectral line equivalent widths as a function of R and RI, including the [OIII] narrow line doublet and the CIV {\lambda}1549 and MgII {\lambda}2799 broad lines. We find that the rest equivalent widths (EWs) of the broad lines correlate positively with RI at the 4-8{\sigma} level. But we find no strong dependence of EW on R, in contrast to previous results. One interpretation is that EWs increase as the line-of-sight angle to the radio-jet axis decreases. These results are in stark contrast to commonly accepted orientation-based theories, which suggest that continuum emission should increase as the angle to the radio-jet axis decreases, resulting in smaller EW of emission lines (assumed isotropic). Finally, we find that the Baldwin effect in our sample does not depend on quasar radio morphology., Comment: 22 pages, 16 figures, 9 tables. Accepted for publication in the Astronomical Journal

We apply Principal Component Analysis (PCA) to ~100,000 stellar spectra obtained by the Sloan Digital Sky Survey (SDSS). In order to avoid strong non-linear variation of spectra with effective temperature, the sample is binned into 0.02 mag wide intervals of the g-r color (-0.20

We conduct a search for radio stars by combining radio and optical data from the Faint Images of the Radio Sky at Twenty cm survey (FIRST) and the Sloan Digital Sky Survey (SDSS). The faint limit of SDSS makes possible a homogeneous search for radio emission from stars of low optical luminosity. We select a sample of 112 candidate radio stars in the magnitude range $15

(Abridged) We construct a catalog of radio sources detected by the GB6 (6 cm), FIRST and NVSS (20 cm), and WENSS (92 cm) radio surveys, and the SDSS optical survey. The 2.7 million entries in the publicly-available master catalog are comprised of the closest three FIRST to NVSS matches (within 30 arcsec) and vice-versa, and unmatched sources from each survey. Entries are supplemented by data from the other radio and optical surveys, where available. We perform data analysis a ~3000 deg^2 region of sky where the surveys overlap, which contains 140,000 NVSS-FIRST sources, of which 64,000 are detected by WENSS and 12,000 by GB6. About one third of each sample is detected by SDSS. An automated classification method based on 20 cm fluxes defines three radio morphology classes: complex, resolved, and compact. Radio color-magnitude- morphology diagrams for these classes show structure suggestive of strong underlying physical correlations. Complex and resolved sources tend to have a steep spectral slope (alpha ~ -0.8) that is nearly constant from 6 to 92 cm, while the compact class contains a significant number of flat-spectrum (alpha ~ 0) sources. In the optically-detected sample, quasars dominate the flat-spectrum compact sources while steep-spectrum and resolved objects contain substantial numbers of both quasars and galaxies. Differential radio counts of quasars and galaxies are similar at bright flux levels (>100 mJy at 20 cm), while at fainter levels the quasar counts are significantly reduced below galaxy counts. The optically-undetected sample is strongly biased toward steep-spectrum sources. In samples of quasars and galaxies with SDSS spectra, we find that radio properties such as spectral slope, morphology, and radio loudness are correlated with optical color and luminosity., Comment: 42 pages, 22 figures, 8 tables; a version with high resolution figures is available at http://www.astro.washington.edu/users/akimball/radiocat/ . Published in the Astronomical Journal. Replacement comment: data URL has changed

The Sloan Digital Sky Survey (SDSS) and Two Micron All Sky Survey (2MASS) are rich resources for studying stellar astrophysics and the structure and formation history of the Galaxy. As new surveys and instruments adopt similar filter sets, it is increasingly important to understand the properties of the ugrizJHKs stellar locus, both to inform studies of `normal' main sequence stars as well as for robust searches for point sources with unusual colors. Using a sample of ~600,000 point sources detected by SDSS and 2MASS, we tabulate the position and width of the ugrizJHKs stellar locus as a function of g-i color, and provide accurate polynomial fits. We map the Morgan-Keenan spectral type sequence to the median stellar locus by using synthetic photometry of spectral standards and by analyzing 3000 SDSS stellar spectra with a custom spectral typing pipeline. We develop an algorithm to calculate a point source's minimum separation from the stellar locus in a seven-dimensional color space, and use it to robustly identify objects with unusual colors, as well as spurious SDSS/2MASS matches. Analysis of a final catalog of 2117 color outliers identifies 370 white-dwarf/M dwarf (WDMD) pairs, 93 QSOs, and 90 M giant/carbon star candidates, and demonstrates that WDMD pairs and QSOs can be distinguished on the basis of their J-Ks and r-z colors. We also identify a group of objects with correlated offsets in the u-g vs. g-r and g-r vs. r-i color-color spaces, but subsequent follow-up is required to reveal the nature of these objects. Future applications of this algorithm to a matched SDSS-UKIDSS catalog may well identify additional classes of objects with unusual colors by probing new areas of color-magnitude space., Comment: 23 pages in emulateapj format, 17 figures, 7 tables. Accepted for publication in the Astronomical Journal. To access a high-resolution version of this paper, as well as machine readable tables and an archive of 'The Hammer' spectral typing suite, see http://www.cfa.harvard.edu/~kcovey v2 -- fixed typos in Table 7 (mainly affecting lines for M8-M10 III stars)