As part of an on-going study of radio transients in Epoch 1 (2017-2019) of the Very Large Array Sky Survey (VLASS), we have discovered a sample of 0.2<z<3.2 active galactic nuclei (AGN) selected in the optical/infrared that have recently brightened dramatically in the radio. These sources would have previously been classified as radio-quiet based on upper limits from the Faint Images of the Radio Sky at Twenty-centimeters (FIRST; 1993-2011) survey; however, they are now consistent with radio-loud quasars. We present a quasi-simultaneous, multi-band (1-18 GHz) VLA follow-up campaign of our sample of AGN with extreme radio variability. We conclude that the radio properties are most consistent with AGN that have recently launched jets within the past few decades, potentially making them among the youngest radio AGN known.

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.

Past studies of compact active galactic nuclei (AGNs), the dominant population at high radio frequencies, selected them using flat spectral index criteria. This biases the sample due to the steepening of AGN spectra at high radio frequencies. We improve upon this by selecting 3610 compact AGNs using their angular size information ($\sim$0.15 arcsec scale) from the Australia Telescope 20 GHz (AT20G) high-angular-resolution catalogue. We cross-match these against the Wide-field Infrared Survey Explorer All-WISE catalogue and present a catalogue with 3300 (91%) matches, 91 (3%) rejects and 219 (6%) nondetections that are excellent high redshift candidates. Of the matched compact AGNs, 92% exhibit QSO mid-infrared colours (W1-W2>0.5). Therefore, our sample of high frequency compact sources has a very high rate of identification with mid-infrared QSOs. We find counterparts for 88% of 387 compact steep-spectrum (CSS) sources in the AT20G survey, 82%$\pm$5% of which exhibit QSO mid-infrared colours and have moderate redshifts (median redshift = 0.82), while those dominated by host galaxy colours in mid-infrared have lower redshifts (median redshift = 0.13). The latter classified into late- and early-type galaxies using their mid-infrared colours shows a majority (68%$\pm$4%) have colours characteristic of late-type galaxies. Thus, we find that a larger fraction of these CSS sources are embedded in hosts with higher gas densities than average early-type galaxies. We compare mid-infrared colours of our AGNs against those reported for AGNs primarily selected using non-radio techniques. This shows that mid-infrared SED of high frequency selected compact radio AGN is comparatively less red, possibly due to contributions from their hosts.

We present an updated version of the Kimball & Ivezic (2008) radio catalog using updated radio and optical data, and new low-frequency radio data. The catalog, containing millions of radio sources, was created by consolidating large-area radio and optical surveys GB6 (6cm), FIRST (20cm), NVSS (20cm), WENSS (92cm), VLSSr (4m), and SDSS DR9 (optical). The region where all surveys overlap covers 3269 square degrees in the North Galactic Cap, and contains >160,000 20-cm sources, with about 12,000 detected in all five radio surveys and over one-third detected optically. Combining parameters from the sky surveys allows easy and efficient classification by radio and optical morphology and radio spectral index. The catalog is available at http://www.aoc.nrao.edu/~akimball/radiocat.shtml .

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey.

The nature and energetics of feedback from thermal winds in quasars can be constrained via observations of the Sunyaev-Zeldovich Effect (SZE) induced by the bubble of thermal plasma blown into the intergalactic medium by the quasar wind. In this letter, we present evidence that we have made the first detection of such a bubble, associated with the hyperluminous quasar HE0515-4414. The SZE detection is corroborated by the presence of extended emission line gas at the same position angle as the wind. Our detection appears on only one side of the quasar, consistent with the SZE signal arising from a combination of thermal and kinetic contributions. Estimates of the energy in the wind allow us to constrain the wind luminosity to the lower end of theoretical predictions, ~0.01% of the bolometric luminosity of the quasar. However, the age we estimate for the bubble, ~0.1 Gyr, and the long cooling time, ~0.6 Gyr, means that such bubbles may be effective at providing feedback between bursts of quasar activity.

The Sunyaev-Zeldovich Effect (SZE) can be used to detect the hot bubbles in the intergalactic medium blown by energetic winds from AGN and starbursts. By directly constraining the kinetic luminosity, age and total energy of the outflow, it offers the promise of greatly increasing our understanding of the effects of wind feedback on galaxy evolution. Detecting the SZE in these winds is very challenging, at the edge of what is possible using existing facilities. The scale of the signal (10-100 kpc) is, however, well matched to interferometers operating at mm wavelengths for objects at z~1. Thus this could become a major science area for the ngVLA, especially if the design of the core is optimized for sensitivity on angular scales of >1 arcsec in the 90 GHz band.

Telescope slew and settle time markedly reduces the efficiency of wide-field multi-epoch surveys for sensitive interferometers with small fields of view. The overheads can be mitigated through the use of On-the-Fly Mosaicing (OTFM), where the the antennas are driven at a non-sidereal rate and visibilities are recorded continuously. Here we introduce the OTFM technique for the VLA, and describe its implementation for the Caltech-NRAO Stripe 82 Survey (CNSS), a dedicated 5-epoch survey for slow transients at S band (2-4 GHz). We also describe the OTFSim tool for planning dynamically-scheduled OTFM observations on the VLA, the latest imaging capabilities for OTFM in CASA, and present a comparison of OTFM observations with pointed observations. Using the subset of our observations from the CNSS pilot and final surveys, we demonstrate that the wide-band and wide-field OTFM observations with the VLA can be imaged accurately, and that this technique offers a more efficient alternative to standard mosaicing for multi-epoch shallow surveys such as the CNSS and the VLA Sky Survey (VLASS). We envisage that the new OTFM mode will facilitate new synoptic surveys and high-frequency mapping experiments on the VLA.

We have observed the environments of a population of 33 heavily dust obscured, ultra-luminous, high-redshift galaxies, selected using WISE and NVSS at $z>$1.3 with the Infra-Red Array Camera on the $Spitzer$ Space Telescope over $\rm5.12\,'\times5.12\,'$ fields. Colour selections are used to quantify any potential overdensities of companion galaxies in these fields. We find no significant excess of galaxies with the standard colour selection for IRAC colours of $\rm[3.6]-[4.5]>-0.1$ consistent with galaxies at $z>$1.3 across the whole fields with respect to wide-area $Spitzer$ comparison fields, but there is a $\rm>2\sigma$ statistical excess within $\rm0.25\,'$ of the central radio-WISE galaxy. Using a colour selection of $\rm[3.6]-[4.5]>0.4$, 0.5 magnitudes redder than the standard method of selecting galaxies at $z>$1.3, we find a significant overdensity, in which $\rm76\%$ ($\rm33\%$) of the 33 fields have a surface density greater than the $\rm3\sigma$ ($\rm5\sigma$) level. There is a statistical excess of these redder galaxies within $\rm0.5\,'$, rising to a central peak $\rm\sim2$--4 times the average density. This implies that these galaxies are statistically linked to the radio-WISE selected galaxy, indicating similar structures to those traced by red galaxies around radio-loud AGN.

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.

Most massive galaxies are now thought to go through an Active Galactic Nucleus (AGN) phase one or more times. Yet, the cause of triggering and the variations in the intrinsic and observed properties of AGN population are still poorly understood. Young, compact radio sources associated with accreting supermassive black holes (SMBHs) represent an important phase in the life cycles of jetted AGN for understanding AGN triggering and duty cycles. The superb sensitivity and resolution of the ngVLA, coupled with its broad frequency coverage, will provide exciting new insights into our understanding of the life cycles of radio AGN and their impact on galaxy evolution. The high spatial resolution of the ngVLA will enable resolved mapping of young radio AGN on sub-kiloparsec scales over a wide range of redshifts. With broad continuum coverage from 1 to 116 GHz, the ngVLA will excel at estimating ages of sources as old as $30-40$ Myr at $z \sim 1$. In combination with lower-frequency ($\nu < 1$ GHz) instruments such as ngLOBO and the Square Kilometer Array, the ngVLA will robustly characterize the spectral energy distributions of young radio AGN.

Energetic feedback by Active Galactic Nuclei (AGNs) plays an important evolutionary role in the regulation of star formation (SF) on galactic scales. However, the effects of this feedback as a function of redshift and galaxy properties such as mass, environment and cold gas content remain poorly understood. The broad frequency coverage (1 to 116 GHz), high sensitivity (up to ten times higher than the Karl G. Jansky Very Large Array), and superb angular resolution (maximum baselines of at least a few hundred km) of the proposed next generation Very Large Array (ngVLA) are uniquely poised to revolutionize our understanding of AGNs and their role in galaxy evolution. Here, we provide an overview of the science related to AGN feedback that will be possible in the ngVLA era and present new continuum ngVLA imaging simulations of resolved radio jets spanning a wide range of intrinsic extents. We also consider key computational challenges and discuss exciting opportunities for multi-wavelength synergy with other next-generation instruments, such as the Square Kilometer Array and the James Webb Space Telescope. The unique combination of high-resolution, large collecting area, and wide frequency range will enable significant advancements in our understanding of the effects of jet-driven feedback on sub-galactic scales, particularly for sources with extents of a few pc to a few kpc such as young and/or lower-power radio AGNs, AGNs hosted by low-mass galaxies, radio jets that are interacting strongly with the interstellar medium of the host galaxy, and AGNs at high redshift.

We present the first results from our on-going Australia Telescope Compact Array survey of CO(1-0) in ALMA-identified submillimetre galaxies in the Extended Chandra Deep Field South. Strong detections of CO(1-0) emission from two submillimetre galaxies, ALESS 122.1 (z = 2.0232) and ALESS 67.1 (z = 2.1230), were obtained. We estimate gas masses of M_gas ~ 1.3 \times 10^{11} M_odot and M_gas ~ 1.0 \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively, adopting alpha_CO = 1.0. Dynamical mass estimates from the kinematics of the CO(1-0) line yields M_dyn (sin i)^2 = 2.1 +- 1.1 \times 10^{11} M_odot and (3.2 +- 0.9) \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively. This is consistent with the total baryonic mass estimates of these two systems. We examine star formation efficiency using the L_FIR versus L'_CO(1-0) relation for samples of local ULIRGs and LIRGs, and more distant star-forming galaxies, with CO(1-0) detections. We find some evidence of a shallower slope for ULIRGs and SMGs compared to less luminous systems, but a larger sample is required for definite conclusions. We determine gas-to-dust ratios of 170 +- 30 and 140 +- 30 for ALESS 122.1 and ALESS 67.1, respectively, showing ALESS 122.1 has an unusually large gas reservoir. By combining the 38.1 GHz continuum detection of ALESS 122.1 with 1.4 and 5.5 GHz data, we estimate that the free-free contribution to radio emission at 38.1 GHz is 34 +- 17 microJy, yielding a star formation rate (1400 +- 700 M_\odot yr^{-1}) consistent with that from the infrared luminosity.

We discuss 6 GHz JVLA observations covering a volume-limited sample of 178 low redshift ($0.2 < z < 0.3$) optically selected QSOs. Our 176 radio detections fall into two clear categories: (1) About $20$\% are radio-loud QSOs (RLQs) having spectral luminosities $L_6 \gtrsim 10^{\,23.2} \mathrm{~W~Hz}^{-1}$ primarily generated in the active galactic nucleus (AGN) responsible for the excess optical luminosity that defines a \emph{bona fide} QSO. (2) The radio-quiet QSOs (RQQs) have $10^{\,21} \lesssim L_6 \lesssim 10^{\,23.2} \mathrm{~W~Hz}^{-1}$ and radio sizes $\lesssim 10 \mathrm{~kpc}$, and we suggest that the bulk of their radio emission is powered by star formation in their host galaxies. "Radio silent" QSOs ($L_6 \lesssim 10^{\,21} \mathrm{~W~Hz}^{-1}$) are rare, so most RQQ host galaxies form stars faster than the Milky Way; they are not "red and dead" ellipticals. Earlier radio observations did not have the luminosity sensitivity $L_6 \lesssim 10^{\,21} \mathrm{~W~Hz}^{-1}$ needed to distinguish between such RLQs and RQQs. Strong, generally double-sided, radio emission spanning $\gg 10 \mathrm{~kpc}$ was found associated with 13 of the 18 RLQ cores having peak flux densities $S_\mathrm{p} > 5 \mathrm{~mJy~beam}^{-1}$ ($log(L) \gtrsim 24$). The radio luminosity function of optically selected QSOs and the extended radio emission associated with RLQs are both inconsistent with simple "unified" models that invoke relativistic beaming from randomly oriented QSOs to explain the difference between RLQs and RQQs. Some intrinsic property of the AGNs or their host galaxies must also determine whether or not a QSO appears radio loud.

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.

The Boolardy Engineering Test Array is a 6 x 12 m dish interferometer and the prototype of the Australian Square Kilometre Array Pathfinder (ASKAP), equipped with the first generation of ASKAP's phased array feed (PAF) receivers. These facilitate rapid wide-area imaging via the deployment of simultaneous multiple beams within a 30 square degree field of view. By cycling the array through 12 interleaved pointing positions and using 9 digitally formed beams we effectively mimic a traditional 1 hour x 108 pointing survey, covering 150 square degrees over 711 - 1015 MHz in 12 hours of observing time. Three such observations were executed over the course of a week. We verify the full bandwidth continuum imaging performance and stability of the system via self-consistency checks and comparisons to existing radio data. The combined three epoch image has arcminute resolution and a 1-sigma thermal noise level of 375 micro-Jy per beam, although the effective noise is a factor 3 higher due to residual sidelobe confusion. From this we derive a catalogue of 3,722 discrete radio components, using the 35 percent fractional bandwidth to measure in-band spectral indices for 1,037 of them. A search for transient events reveals one significantly variable source within the survey area. The survey covers approximately two-thirds of the Spitzer South Pole Telescope Deep Field. This pilot project demonstrates the viability and potential of using PAFs to rapidly and accurately survey the sky at radio wavelengths.

The Next-Generation Very Large Array (ngVLA) will be critical for understanding how galaxies are built and evolve at the earliest epochs. The sensitivity and frequency coverage will allow for the detection of cold gas and dust in `normal' distant galaxies, including the low-J transitions of molecular gas tracers such as CO, HNC, and HCO+; synchrotron and free-free continuum emission; and even the exciting possibility of thermal dust emission at the highest (z~7) redshifts. In particular, by enabling the total molecular gas reservoirs to be traced to unprecedented sensitivities across a huge range of epochs simultaneously -- something no other radio or submillimeter facility will be capable of -- the detection of the crucial low-J transitions of CO in a diverse body of galaxies will be the cornerstone of ngVLA's contribution to high-redshift galaxy evolution science. The ultra-wide bandwidths will allow a complete sampling of radio SEDs, as well as the detection of emission lines necessary for spectroscopic confirmation of elusive dusty starbursts. The ngVLA will also deliver unique contributions to our understanding of cosmic magnetism and to science accessible through microwave polarimetry. Finally, the superb angular resolution will move the field beyond detection experiments and allow detailed studies of the morphology and dynamics of these systems, including dynamical modeling of disks/mergers, determining the properties of outflows, measuring black hole masses from gas disks, and resolving multiple AGN nuclei. We explore the contribution of a ngVLA to these areas and more, as well as synergies with current and upcoming facilities including ALMA, SKA, large single-dish submillimeter observatories, GMT/TMT, and JWST.

We present ALMA 870 micron (345 GHz) data for 49 high redshift (0.47<z<2.85), luminous (11.7 < log L(bol) (Lsun) < 14.2) radio-powerful AGN, obtained to constrain cool dust emission from starbursts concurrent with highly obscured radiative-mode black hole (BH) accretion in massive galaxies which possess a small radio jet. The sample was selected from WISE with extremely steep (red) mid-infrared (MIR) colors and with compact radio emission from NVSS/FIRST. Twenty-six sources are detected at 870 microns, and we find that the sample has large mid- to far-infrared luminosity ratios consistent with a dominant and highly obscured quasar. The rest-frame 3 GHz radio powers are 24.7 < log P3.0 GHz (W/Hz) < 27.3, and all sources are radio-intermediate or radio-loud. BH mass estimates are 7.7 < log M(BH) (Msun) < 10.2. The rest frame 1-5 um SEDs are very similar to the "Hot DOGs" (Hot Dust Obscured Galaxies), and steeper (redder) than almost any other known extragalactic sources. ISM masses estimated for the ALMA detected sources are 9.9 < log M(ISM) (Msun) < 11.75 assuming a dust temperature of 30K. The cool dust emission is consistent with star formation rates (SFRs) reaching several thousand Msun/yr, depending on the assumed dust temperature, however we cannot rule out the alternative that the AGN powers all the emission in some cases. Our best constrained source has radiative transfer solutions with ~ equal contributions from an obscured AGN and a young (10-15 Myr) compact starburst.

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.

We present JCMT SCUBA-2 850microns submillimetre (submm) observations of 30 mid-infrared (mid-IR) luminous AGN, detected jointly by the WISE all-sky IR survey and the NVSS/FIRST radio survey. These rare sources are selected by their extremely red mid-infrared spectral energy distributions (SEDs) and compact radio counterparts. Further investigations show that they are highly obscured, have abundant warm AGN-heated dust and are thought to be experiencing intense AGN feedback. These galaxies appear to be consistent with an AGN-dominated galaxy, and could be a transient phase of merging galaxies. When comparing the number of submm galaxies (SMGs) detected serendipitously in the surrounding 1.5-arcmin to those in blank-field submm surveys, there is a very significant overdensity, of order 5, but no sign of radial clustering centred at our primary objects. The WISE/radio-selected AGN thus reside in 10-Mpc-scale overdense environments, that could be forming in pre-viralised clusters of galaxies. WISE/radio-selected AGNs appear to be the strongest signposts of high-density regions of active, luminous and dusty galaxies. SCUBA-2 850microns observations indicate that their submm fluxes are low compared to many popular AGN SED templates, hence the WISE/radio-selected AGNs have either less cold and/or more warm dust emission than normally assumed for typical AGN. Most of the targets are not detected, only four targets are detected at SCUBA-2 850microns, and have total IR luminosities >= 10^13 L_solar, if their redshifts are consistent with the subset of the 10 SCUBA-2 undetected targets with known redshifts, z ~ 0.44 - 2.86.

We present the detection of $^{12}$CO(2-1) in the $z = 4.44$ submillimetre galaxy ALESS65.1 using the Australia Telescope Compact Array. A previous ALMA study of submillimetre galaxies in the Extended Chandra Deep Field South determined the redshift of this optically and near-infrared undetected source through the measurement of [CII] 157.74 $\mu$m emission. Using the luminosity of the $^{12}$CO(2-1) emission we estimate the gas mass to be $M_{\rm gas} \sim 1.7 \times 10^{10}$ ${\rm M}_\odot$. The gas depletion timescale of ALESS65.1 is $\sim$ 25 Myr, similar to other high redshift submillimetre galaxies and consistent with $z > 4$ SMGs being the progenitors of massive "red-and-dead" galaxies at $z > 2$. The ratio of the [CII], $^{12}$CO and far-infrared luminosities implies a strong far-ultraviolet field of $G_0 \sim 10^{3.25}$, which is at the high end of the far-ultraviolet fields seen in local starbursts, but weaker than the far-ultraviolet fields of most nearby ULIRGs. The high ratio of $L_{\rm [CII]}/L_{\rm FIR} = 1.0 \times 10^{-3}$ observed in ALESS65.1, combined with $L_{\rm [CII]}/L_{\rm CO} \sim 2300$, is consistent with ALESS65.1 having more extended regions of intense star formation than local ULIRGs.

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.

A community meeting on the topic of "Radio Astronomy in the LSST Era" was hosted by the National Radio Astronomy Observatory in Charlottesville, VA (2013 May 6--8). The focus of the workshop was on time domain radio astronomy and sky surveys. For the time domain, the extent to which radio and visible wavelength observations are required to understand several classes of transients was stressed, but there are also classes of radio transients for which no visible wavelength counterpart is yet known, providing an opportunity for discovery. From the LSST perspective, the LSST is expected to generate as many as 1 million alerts nightly, which will require even more selective specification and identification of the classes and characteristics of transients that can warrant follow up, at radio or any wavelength. The LSST will also conduct a deep survey of the sky, producing a catalog expected to contain over 38 billion objects in it. Deep radio wavelength sky surveys will also be conducted on a comparable time scale, and radio and visible wavelength observations are part of the multi-wavelength approach needed to classify and understand these objects. Radio wavelengths are valuable because they are unaffected by dust obscuration and, for galaxies, contain contributions both from star formation and from active galactic nuclei. The workshop touched on several other topics, on which there was consensus including the placement of other LSST "Deep Drilling Fields," inter-operability of software tools, and the challenge of filtering and exploiting the LSST data stream. There were also topics for which there was insufficient time for full discussion or for which no consensus was reached, which included the procedures for following up on LSST observations and the nature for future support of researchers desiring to use LSST data products.

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.

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.

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.

We use SDSS photometry of 73 million stars to simultaneously obtain best-fit main-sequence stellar energy distribution (SED) and amount of dust extinction along the line of sight towards each star. Using a subsample of 23 million stars with 2MASS photometry, whose addition enables more robust results, we show that SDSS photometry alone is sufficient to break degeneracies between intrinsic stellar color and dust amount when the shape of extinction curve is fixed. When using both SDSS and 2MASS photometry, the ratio of the total to selective absorption, $R_V$, can be determined with an uncertainty of about 0.1 for most stars in high-extinction regions. These fits enable detailed studies of the dust properties and its spatial distribution, and of the stellar spatial distribution at low Galactic latitudes. Our results are in good agreement with the extinction normalization given by the Schlegel et al. (1998, SFD) dust maps at high northern Galactic latitudes, but indicate that the SFD extinction map appears to be consistently overestimated by about 20% in the southern sky, in agreement with Schlafly et al. (2010). The constraints on the shape of the dust extinction curve across the SDSS and 2MASS bandpasses support the models by Fitzpatrick (1999) and Cardelli et al. (1989). For the latter, we find an $R_V=3.0\pm0.1$(random) $\pm0.1$(systematic) over most of the high-latitude sky. At low Galactic latitudes (|b|<5), we demonstrate that the SFD map cannot be reliably used to correct for extinction as most stars are embedded in dust, rather than behind it. We introduce a method for efficient selection of candidate red giant stars in the disk, dubbed "dusty parallax relation", which utilizes a correlation between distance and the extinction along the line of sight. We make these best-fit parameters, as well as all the input SDSS and 2MASS data, publicly available in a user-friendly format.

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.

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.

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

We extend the Unified Radio Catalog, a catalog of sources detected by various (NVSS, FIRST, WENSS, GB6) radio surveys, and SDSS, to IR wavelengths by matching it to the IRAS Point and Faint Source catalogs. By fitting each NVSS-selected galaxy's NUV-NIR spectral energy distribution (SED) with stellar population synthesis models we add to the catalog star formation rates, stellar masses, and attenuations.We further add information about optical emission line properties for NVSS-selected galaxies with available SDSS spectroscopy. Using an NVSS 20cm (F_{1.4GHz} ge 2.5mJy) selected sample, matched to the SDSS spectroscopic ("main" galaxy and quasar) catalogs and IRAS data (0.04<z le 0.2) we perform an in depth analysis of the radio-FIR correlation for various types of galaxies, separated into i) quasars, ii) star forming, iii) composite, iv) Seyfert, v) LINER and vi) absorption line galaxies using the standard optical spectroscopic diagnostic tools. We utilize SED-based star formation rates to independently quantify the source of radio and FIR emission in our galaxies. Our results show that Seyfert galaxies have FIR/radio ratios lower than, but still within the scatter of, the canonical value due to an additional (likely AGN) contribution to their radio continuum emission. Furthermore, IR-detected absorption and LINER galaxies are on average strongly dominated by AGN activity in both their FIR and radio emission; however their average FIR/radio ratio is consistent with that expected for star forming galaxies. In summary, we find that most AGN-containing galaxies in our NVSS-IRAS-SDSS sample have FIR/radio flux ratios indistinguishable from those of the star-forming galaxies that define the radio-FIR correlation. Thus, attempts to separate AGNs from star-forming galaxies by their FIR/radio flux ratios alone can separate only a small fraction of the AGNs, such as the radio-loud quasars.

We model the time variability of ~9,000 spectroscopically confirmed quasars in SDSS Stripe 82 as a damped random walk. Using 2.7 million photometric measurements collected over 10 years, we confirm the results of Kelly et al. (2009) and Koz{\l}owski et al. (2010) that this model can explain quasar light curves at an impressive fidelity level (0.01-0.02 mag). The damped random walk model provides a simple, fast [O(N) for N data points], and powerful statistical description of quasar light curves by a characteristic time scale (tau) and an asymptotic rms variability on long time scales (SF_inf). We searched for correlations between these two variability parameters and physical parameters such as luminosity and black hole mass, and rest-frame wavelength. We find that tau increases with increasing wavelength with a power law index of 0.17, remains nearly constant with redshift and luminosity, and increases with increasing black hole mass with power law index of 0.21+/-0.07. The amplitude of variability is anti-correlated with the Eddington ratio, which suggests a scenario where optical fluctuations are tied to variations in the accretion rate. The radio-loudest quasars have systematically larger variability amplitudes by about 30%, when corrected for the other observed trends, while the distribution of their characteristic time scale is indistinguishable from that of the full sample. We do not detect any statistically robust differences in the characteristic time scale and variability amplitude between the full sample and the small subsample of quasars detected by ROSAT. Our results provide a simple quantitative framework for generating mock quasar light curves, such as currently used in LSST image simulations. (abridged)

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<g-r<0.90, roughly corresponding to MK spectral types A3 to K3), and PCA is applied independently for each bin. In each color bin, the first four eigenspectra are sufficient to describe the observed spectra within the measurement noise. We discuss correlations of eigencoefficients with metallicity and gravity estimated by the Sloan Extension for Galactic Understanding and Exploration (SEGUE) Stellar Parameters Pipeline. The resulting high signal-to-noise mean spectra and the other three eigenspectra are made publicly available. These data can be used to generate high quality spectra for an arbitrary combination of effective temperature, metallicity, and gravity within the parameter space probed by the SDSS. The SDSS stellar spectroscopic database and the PCA results presented here offer a convenient method to classify new spectra, to search for unusual spectra, to train various spectral classification methods, and to synthesize accurate colors in arbitrary optical bandpasses.

We present a sample of 723 optically selected BL Lac candidates from the SDSS DR7 spectroscopic database encompassing 8250 deg^2 of sky; our sample constitutes one of the largest uniform BL Lac samples yet derived. Each BL Lac candidate has a high-quality SDSS spectrum from which we determine spectroscopic redshifts for ~60% of the objects. Redshift lower limits are estimated for the remaining objects utilizing the lack of host galaxy flux contamination in their optical spectra; we find that objects lacking spectroscopic redshifts are likely at systematically higher redshifts. Approximately 80% of our BL Lac candidates match to a radio source in FIRST/NVSS, and ~40% match to a ROSAT X-ray source. The homogeneous multiwavelength coverage allows subdivision of the sample into 637 radio-loud BL Lac candidates and 86 weak-featured radio-quiet objects. The radio-loud objects broadly support the standard paradigm unifying BL Lac objects with beamed radio galaxies. We propose that the majority of the radio-quiet objects may be lower-redshift (z<2.2) analogs to high-redshift weak line quasars (i.e., AGN with unusually anemic broad emission line regions). These would constitute the largest sample of such objects, being of similar size and complementary in redshift to the samples of high-redshift weak line quasars previously discovered by the SDSS. However, some fraction of the weak-featured radio-quiet objects may instead populate a rare and extreme radio-weak tail of the much larger radio-loud BL Lac population. Serendipitous discoveries of unusual white dwarfs, high-redshift weak line quasars, and broad absorption line quasars with extreme continuum dropoffs blueward of rest-frame 2800 Angstroms are also briefly described.

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<i\lesssim19.1$ and with radio flux $S_{20}\geq1.25$mJy, from about 7000 deg$^2$ of sky. The selection criteria are positional coincidence within $1\arcsec$, radio and optical point source morphology, and an SDSS spectrum classified as stellar. The sample contamination is estimated by random matching to be $108\pm13$, suggesting that at most a small fraction of the selected candidates are genuine radio stars. Therefore, we rule out a very rare population of extremely radio-loud stars: no more than 1.2 of every million stars in the magnitude range $15<i<19.1$ stars has radio flux $S_{20}\geq1.25$ mJy. We investigate the optical and radio colors of the sample to find candidates that show the largest likelihood of being real radio stars. The significant outliers from the stellar locus, as well as the magnetically active stars, are the best candidates for follow-up radio observations. We conclude that, while the present wide-area radio surveys are not sensitive enough to provide homogeneous samples of the extremely rare radio stars, upcoming surveys which exploit the great sensitivity of current and planned telescopes do have sufficient sensitivity and will allow the properties of this class of object to be investigated in detail.

The four-image gravitationally lensed quasar QSO 2237+0305 is microlensed by stars in the lens galaxy. The amplitude of microlensing variability can be used to infer the relative size of the quasar as a function of wavelength; this provides a test of quasar models. Toward this end, we present Spitzer Space Telescope Infrared Spectrograph and Infrared Array Camera (IRAC) observations of QSO 2237+0305, finding the following. (1) The infrared (IR) spectral energy distribution (SED) is similar to that of other bright radio-quiet quasars, contrary to an earlier claim. (2) A dusty torus model with a small opening angle fits the overall shape of the IR SED well, but the quantitative agreement is poor due to an offset in wavelength of the silicate feature. (3) The flux ratios of the four lensed images can be derived from the IRAC data despite being unresolved. We find that the near-IR fluxes are increasingly affected by microlensing toward shorter wavelengths. (4) The wavelength dependence of the IRAC flux ratios is consistent with the standard quasar model in which an accretion disk and a dusty torus both contribute near 1 micron in the rest frame. This is also consistent with recent IR spectropolarimetry of nearby quasars.

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

Using effective temperature and metallicity derived from SDSS spectra for ~60,000 F and G type main sequence stars (0.2<g-r<0.6), we develop polynomial models for estimating these parameters from the SDSS u-g and g-r colors. We apply this method to SDSS photometric data for about 2 million F/G stars and measure the unbiased metallicity distribution for a complete volume-limited sample of stars at distances between 500 pc and 8 kpc. The metallicity distribution can be exquisitely modeled using two components with a spatially varying number ratio, that correspond to disk and halo. The two components also possess the kinematics expected for disk and halo stars. The metallicity of the halo component is spatially invariant, while the median disk metallicity smoothly decreases with distance from the Galactic plane from -0.6 at 500 pc to -0.8 beyond several kpc. The absence of a correlation between metallicity and kinematics for disk stars is in a conflict with the traditional decomposition in terms of thin and thick disks. We detect coherent substructures in the kinematics--metallicity space, such as the Monoceros stream, which rotates faster than the LSR, and has a median metallicity of [Fe/H]=-0.96, with an rms scatter of only ~0.15 dex. We extrapolate our results to the performance expected from the Large Synoptic Survey Telescope (LSST) and estimate that the LSST will obtain metallicity measurements accurate to 0.2 dex or better, with proper motion measurements accurate to ~0.2 mas/yr, for about 200 million F/G dwarf stars within a distance limit of ~100 kpc (g<23.5). [abridged]

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.

In addition to optical photometry of unprecedented quality, the Sloan Digital Sky Survey (SDSS) is also producing a massive spectroscopic database. We discuss determination of stellar parameters, such as effective temperature, gravity and metallicity from SDSS spectra, describe correlations between kinematics and metallicity, and study their variation as a function of the position in the Galaxy. We show that stellar parameter estimates by Beers et al. show a good correlation with the position of a star in the g-r vs. u-g color-color diagram, thereby demonstrating their robustness as well as a potential for photometric parameter estimation methods. Using Beers et al. parameters, we find that the metallicity distribution of the Milky Way stars at a few kpc from the galactic plane is bimodal with a local minimum at [Z/Zo]~ -1.3. The median metallicity for the low-metallicity [Z/Zo]<-1.3 subsample is nearly independent of Galactic cylindrical coordinates R and z, while it decreases with z for the high-metallicity [Z/Zo]> -1.3 sample. We also find that the low-metallicity sample has ~2.5 times larger velocity dispersion and that it does not rotate (at the ~10 km/s level), while the rotational velocity of the high-metallicity sample decreases smoothly with the height above the galactic plane.

Imaging data from the Sloan Digital Sky Survey is used to measure the empirical size-richness relation for a large sample of galaxy clusters. Using population subtraction methods, we determine the radius at which the cluster galaxy number density is 200/Omega_m times the mean galaxy density, without assuming a model for the radial distribution of galaxies in clusters. If these galaxies are unbiased on Mpc scales, this galaxy-density-based R_200 reflects the characteristic radii of clusters. We measure the scaling of this characteristic radius with richness over an order of magnitude in cluster richness, from rich clusters to poor groups. We use this information to examine the radial profiles of galaxies in clusters as a function of cluster richness, finding that the concentration of the galaxy distribution decreases with richness and is systematically lower than the concentrations measured for dark matter profiles in N-body simulations. Using these scaled radii, we investigate the behavior of the cluster luminosity function, and find that it is well matched by a Schechter function for galaxies brighter than M_r = -18 only after the central galaxy has been removed. We find that the luminosity function varies with richness and with distance from the cluster center, underscoring the importance of using an aperture that scales with cluster mass to compare physically equivalent regions of these different systems. We note that the lowest richness systems in our catalog have properties consistent with those expected of the earliest-forming halos; our cluster-finding algorithm, in addition to reliably finding clusters, may be efficient at finding fossil groups.

We discuss the UV, optical, and IR properties of the SDSS sources detected by GALEX as part of its All-sky Imaging Survey Early Release Observations. Virtually all of the GALEX sources in the overlap region are detected by SDSS. GALEX sources represent ~2.5% of all SDSS sources within these fields and about half are optically unresolved. Most unresolved GALEX/SDSS sources are bright blue turn-off thick disk stars and are typically detected only in the GALEX near-UV band. The remaining unresolved sources include low-redshift quasars, white dwarfs, and white dwarf/M dwarf pairs, and these dominate the optically unresolved sources detected in both GALEX bands. Almost all the resolved SDSS sources detected by GALEX are fainter than the SDSS 'main' spectroscopic limit. These sources have colors consistent with those of blue (spiral) galaxies (u-r<2.2), and most are detected in both GALEX bands. Measurements of their UV colors allow much more accurate and robust estimates of star-formation history than are possible using only SDSS data. Indeed, galaxies with the most recent (<20 Myr) star formation can be robustly selected from the GALEX data by requiring that they be brighter in the far-UV than in the near-UV band. However, older starburst galaxies have UV colors similar to AGN, and thus cannot be selected unambiguously on the basis of GALEX fluxes alone. With the aid of 2MASS data, we construct and discuss median 10 band UV-optical-IR spectral energy distributions for turn-off stars, hot white dwarfs, low-redshift quasars, and spiral and elliptical galaxies. We point out the high degree of correlation between the UV color and the contribution of the UV flux to the UV-optical-IR flux of galaxies detected by GALEX.