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.

This Astro2020 State of the Profession Consideration White Paper highlights the growth of astrostatistics and astroinformatics in astronomy, identifies key issues hampering the maturation of these new subfields, and makes recommendations for structural improvements at different levels that, if acted upon, will make significant positive impacts across astronomy.

Long gamma-ray bursts (GRBs) are the most luminous known electromagnetic radiation sources in the Universe for the 3 to 300 sec of their prompt flashes (isotropic X/ gamma-ray luminosities up to 10^53 ergs/sec). Their afterglows have first day rest-frame UV/optical absolute magnitudes AB = -30 to -23. This luminous continuum nUV-nIR back-light provides the ultimate probe of the SFR(z) back to the first Pop III to II.5 stars, expected to be massive and GRB progenitors. GRB afterglow spectra in the first 1 to 3 hours will directly measure their host galaxy ionization fraction x_i vs. z in the Epoch of Reionization (EOR), tracing the growth of structure. Only 28% of Swift GRBs have measured redshifts due to limited followup at R, J >21. Some ~25% of GRBs are optically dark due to dust absorption in their host galaxies, but those with low NH in their X-ray spectra are likely at z >7. Current 8-10m telescopes and coming ELTs cannot pursue optically dark GRBs promptly, nor can JWST or WFIRST slew within 0.5 to 1 days of a GRB. The Time-domain Spectroscopic Observatory (TSO) is a proposed Probe-class 1.3m telescope at L2, with imaging and spectroscopy (R = 200, 1800) in 4 bands (0.3 to 5um) and rapid slew capability to 90% of sky. TSO would finally utilize z > 6 to 12 GRBs as the most direct probe of the SFR(z), EOR(z), and possibly the first direct detection of the core collapse of the very first (Pop III) stars.

Time-domain Astrophysics (TDA), a foundation of Astronomy, has become a major part of current and projected (2020s) astrophysics. While much has been derived from temporal measures of flux and color, the real physics comes from spectroscopy. With LSST coming on line in 2022, with TDA as one of its original drivers, the deluge of Transients and new types of variables will be truly astronomical. With multi-wavelength targeted EM surveys and multi-messenger (e.g. LIGO-international and advanced ICECUBE), and the possibility of full-sky and full-time X-ray imagers, the discovery of new Transients and Variables will flood telescopes on the ground and in space, and this just for multi-band imaging without spectroscopy. In this White Paper we briefly summarize several long-standing major science objectives that can be realized with TDA imaging and spectroscopy (near UV to mid IR) from space. We provide a brief description of how these can be achieved with the Time-domain Spectroscopic Observatory (TSO), a Probe-Class mission concept that ELTs on the ground and Flagship missions in space can not achieve on their own.

With the epochal first detection of gravitational waves from a binary neutron star (NS) merger with the GW170817 event, and its direct confirmation that NS-NS mergers are significant sources of the of the r-process nucleosynthesis of heavy elements, an immense new arena for prompt EM (X-rays through IR and radio) studies of fundamental physics has been opened. Over the next decade, GW observatories will expand in scale and sensitivity so the need for facilities that can provide prompt, high sensitivity, broad-band EM followup becomes more urgent. NS-NS or NS-black hole (BH) mergers will be instantly recognized (and announced) by the LIGO-international collaboration. LSST will be a prime resource for rapid tiling of what will usually be large (~10-100 degree squared) error boxes. X-ray through IR Telescopes in space with (nearly) full-sky access that can rapidly image and tile are crucial for providing the earliest imaging and spectroscopic studies of the kilonova emission immediately following NS-NS mergers. The Time-domain Spectroscopic Observatory (TSO) is a proposed Probe-class 1.3 m telescope at L2, with imaging and spectroscopy (R = 200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky. TSO nUV-mid-IR spectra will enable new constraints on NS structure and nucleosynthesis.

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 sq. deg field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hbeta to z=8.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

We report on $g$, $r$ and $i$ band observations of the Interstellar Object 'Oumuamua (1I) taken on 2017 October 29 from 04:28 to 08:40 UTC by the Apache Point Observatory (APO) 3.5m telescope's ARCTIC camera. We find that 1I's colors are $g-r=0.41\pm0.24$ and $r-i=0.23\pm0.25$, consistent with the visible spectra of Masiero (2017), Ye et al. (2017) and Fitzsimmons et al. (2017), and most comparable to the population of Solar System C/D asteroids, Trojans, or comets. We find no evidence of any cometary activity at a heliocentric distance of 1.46 au, approximately 1.5 months after 1I's closest approach distance to the Sun. Significant brightness variability was seen in the $r$ observations, with the object becoming notably brighter towards the end of the run. By combining our APO photometric time series data with the Discovery Channel Telescope (DCT) data of Knight et al. (2017), taken 20 h later on 2017 October 30, we construct an almost complete light curve with a most probable lightcurve period of $P \simeq 4~{\rm h}$. Our results imply a double peaked rotation period of 8.1 $\pm$ 0.02 h, with a peak-to-peak amplitude of 1.5 - 2.1 mags. Assuming that 1I's shape can be approximated by an ellipsoid, the amplitude constraint implies that 1I has an axial ratio of 3.5 to 10.3, which is strikingly elongated. Assuming that 1I is rotating above its critical break up limit, our results are compatible with 1I having having modest cohesive strength and may have obtained its elongated shape during a tidal disruption event before being ejected from its home system. Astrometry useful for constraining 1I's orbit was also obtained and published in Weaver et al. (2017).

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

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.

Aims. The Sloan Digital Sky Survey (SDSS) offers a unique possibility of not only detecting Colour Induced Displacement (CID) double stars but also confirming these detections. Methods. Successive cuts are applied to the SDSS DR12 database in order to reduce the size of the sample to be considered. The resulting dataset is then screened with a criterion based on the distance and orientation of the photocentres in different photometric bands. Results. About 3 200 distinct objects are classified as CID double stars, 40 of which are confirmed with at least a second detection. A consistency check further validates these detections.

In area and depth, the Pan-STARRS1 (PS1) 3$\pi$ survey is unique among many-epoch, multi-band surveys and has enormous potential for all-sky identification of variable sources. PS1 has observed the sky typically seven times in each of its five bands ($grizy$) over 3.5 years, but unlike SDSS not simultaneously across the bands. Here we develop a new approach for quantifying statistical properties of non-simultaneous, sparse, multi-color lightcurves through light-curve structure functions, effectively turning PS1 into a $\sim 35$-epoch survey. We use this approach to estimate variability amplitudes and timescales $(\omega_r, \tau)$ for all point-sources brighter than $r_{\mathrm{P1}}=21.5$ mag in the survey. With PS1 data on SDSS Stripe 82 as ``ground truth", we use a Random Forest Classifier to identify QSOs and RR Lyrae based on their variability and their mean PS1 and WISE colors. We find that, aside from the Galactic plane, QSO and RR Lyrae samples of purity $\sim$75\% and completeness $\sim$92\% can be selected. On this basis we have identified a sample of $\sim 1,000,000$ QSO candidates, as well as an unprecedentedly large and deep sample of $\sim$150,000 RR Lyrae candidates with distances from $\sim$10 kpc to $\sim$120 kpc. Within the Draco dwarf spheroidal, we demonstrate a distance precision of 6\% for RR Lyrae candidates. We provide a catalog of all likely variable point sources and likely QSOs in PS1, a total of $25.8\times 10^6$ sources.

The Large Synoptic Survey Telescope (LSST) will be a ground-based, optical, all-sky, rapid cadence survey project with tremendous potential for discovering and characterizing asteroids. With LSST's large 6.5m diameter primary mirror, a wide 9.6 square degree field of view 3.2 Gigapixel camera, and rapid observational cadence, LSST will discover more than 5 million asteroids over its ten year survey lifetime. With a single visit limiting magnitude of 24.5 in r-band, LSST will be able to detect asteroids in the Main Belt down to sub-kilometer sizes. The current strawman for the LSST survey strategy is to obtain two visits (each visit being a pair of back-to-back 15s exposures) per field, separated by about 30 minutes, covering the entire visible sky every 3-4 days throughout the observing season, for ten years. The catalogs generated by LSST will increase the known number of small bodies in the Solar System by a factor of 10-100 times, among all populations. The median number of observations for Main Belt asteroids will be on the order of 200-300, with Near Earth Objects receiving a median of 90 observations. These observations will be spread among ugrizy bandpasses, providing photometric colors and allowing sparse lightcurve inversion to determine rotation periods, spin axes, and shape information. These catalogs will be created using automated detection software, the LSST Moving Object Processing System (MOPS), that will take advantage of the carefully characterized LSST optical system, cosmetically clean camera, and recent improvements in difference imaging. Tests with the prototype MOPS software indicate that linking detections (and thus discovery) will be possible at LSST depths with our working model for the survey strategy, but evaluation of MOPS and improvements in the survey strategy will continue. All data products and software created by LSST will be publicly available.

This paper introduces the multiband periodogram, a general extension of the well-known Lomb-Scargle approach for detecting periodic signals in time-domain data. In addition to advantages of the Lomb-Scargle method such as treatment of non-uniform sampling and heteroscedastic errors, the multiband periodogram significantly improves period finding for randomly sampled multiband light curves (e.g., Pan-STARRS, DES and LSST). The light curves in each band are modeled as arbitrary truncated Fourier series, with the period and phase shared across all bands. The key aspect is the use of Tikhonov regularization which drives most of the variability into the so-called base model common to all bands, while fits for individual bands describe residuals relative to the base model and typically require lower-order Fourier series. This decrease in the effective model complexity is the main reason for improved performance. We use simulated light curves and randomly subsampled SDSS Stripe 82 data to demonstrate the superiority of this method compared to other methods from the literature, and find that this method will be able to efficiently determine the correct period in the majority of LSST's bright RR Lyrae stars with as little as six months of LSST data. A Python implementation of this method, along with code to fully reproduce the results reported here, is available on GitHub.

We provide an overview of the science benefits of combining information from the Square Kilometre Array (SKA) and the Large Synoptic Survey Telescope (LSST). We first summarise the capabilities and timeline of the LSST and overview its science goals. We then discuss the science questions in common between the two projects, and how they can be best addressed by combining the data from both telescopes. We describe how weak gravitational lensing and galaxy clustering studies with LSST and SKA can provide improved constraints on the causes of the cosmological acceleration. We summarise the benefits to galaxy evolution studies of combining deep optical multi-band imaging with radio observations. Finally, we discuss the excellent match between one of the most unique features of the LSST, its temporal cadence in the optical waveband, and the time resolution of the SKA.

We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding 10 kpc and galactocentric distances exceeding 20 kpc. Our analysis utilizes Jeans equations to generate two-dimensional acceleration maps throughout the volume; this approach is thoroughly tested on a cosmologically derived N-body+SPH simulation of a MW-like galaxy. We show that the known accelerations (gradients of the gravitational potential) can be successfully recovered in such a realistic system. Leveraging the baryonic gravitational potential derived by Bovy & Rix (2013), we show that the gravitational potential implied by the SDSS observations cannot be explained, assuming Newtonian gravity, by visible matter alone: the gravitational force experienced by stars at galactocentric distances of 20 kpc is as much as three times stronger than what can be attributed to purely visible matter. We also show that the SDSS data provide a strong constraint on the shape of the dark matter halo potential. Within galactocentric distances of 20 kpc, the dark matter halo potential is well described as an oblate halo with axis ratio qDM=0.7+/-0.1; this corresponds to an axis ratio qDM=0.4+/-0.1 for the dark matter density distribution. Because of our precise two-dimensional measurements of the acceleration of the halo stars, we can reject several MOND models as an explanation of the observed behavior.

Astronomy and astrophysics are witnessing dramatic increases in data volume as detectors, telescopes and computers become ever more powerful. During the last decade, sky surveys across the electromagnetic spectrum have collected hundreds of terabytes of astronomical data for hundreds of millions of sources. Over the next decade, the data volume will enter the petabyte domain, and provide accurate measurements for billions of sources. Astronomy and physics students are not traditionally trained to handle such voluminous and complex data sets. In this paper we describe astroML; an initiative, based on Python and scikit-learn, to develop a compendium of machine learning tools designed to address the statistical needs of the next generation of students and astronomical surveys. We introduce astroML and present a number of example applications that are enabled by this package.

Through matches with the Sloan Digital Sky Survey (SDSS) catalogue we identify the location of various families of astronomical objects in WISE colour space. We identify reliable indicators that separate Galactic/local from extragalactic sources and concentrate here on the objects in our Galaxy and its closest satellites. We develop colour and magnitude criteria that are based only on WISE data to select asymptotic giant branch (AGB) stars with circumstellar dust shells, and separate them into O-rich and C-rich classes. With these criteria we produce an all-sky map for the count ratio of the two populations. The map reveals differences between the Galactic disc, the Magellanic Clouds and the Sgr Dwarf Spheroidal galaxy, as well as a radial gradient in the Large Magellanic Cloud (LMC) disc. We find that the C:O number ratio for dusty AGB stars increases with distance from the LMC centre about twice as fast as measured for near-IR selected samples of early AGB stars. Detailed radiative transfer models show that WISE colours are well explained by the emission of centrally heated dusty shells where the dust has standard properties of interstellar medium (ISM) grains. The segregation of different classes of objects in WISE colour space arises from differences in properties of the dust shells: those around young stellar objects have uniform density distributions while in evolved stars they have steep radial profiles.

We present the fiducial main sequence stellar locus traced by 10 photometric colors observed by SDSS, 2MASS, and WISE. Median colors are determined using 1,052,793 stars with r-band extinction less than 0.125. We use this locus to measure the dust extinction curve relative to the r-band, which is consistent with previous measurements in the SDSS and 2MASS bands. The WISE band extinction coefficients are larger than predicted by standard extinction models. Using 13 lines of sight, we find variations in the extinction curve in H, Ks, and WISE bandpasses. Relative extinction decreases towards Galactic anti-center, in agreement with prior studies. Relative extinction increases with Galactic latitude, in contrast to previous observations. This indicates a universal mid-IR extinction law does not exist due to variations in dust grain size and chemistry with Galactocentric position. A preliminary search for outliers due to warm circumstellar dust is also presented, using stars with high signal-to-noise in the W3-band. We find 199 such outliers, identified by excess emission in Ks-W3. Inspection of SDSS images for these outliers reveals a large number of contaminants due to nearby galaxies. Six sources appear to be genuine dust candidates, yielding a fraction of systems with infrared excess of 0.12$\pm$0.05%.

Context. The variations of solar activity over long time intervals using a solar activity reconstruction based on the cosmogenic radionuclide 10Be measured in polar ice cores are studied. Methods. By applying methods of nonlinear dynamics, the solar activity cycle is studied using solar activity proxies that have been reaching into the past for over 9300 years. The complexity of the system is expressed by several parameters of nonlinear dynamics, such as embedding dimension or false nearest neighbors, and the method of delay coordinates is applied to the time series. We also fit a damped random walk model, which accurately describes the variability of quasars, to the solar 10Be data and investigate the corresponding power spectral distribution. The periods in the data series were searched by the Fourier and wavelet analyses. The solar activity on the long-term scale is found to be on the edge of chaotic behavior. This can explain the observed intermittent period of longer lasting solar activity minima. Filtering the data by eliminating variations below a certain period (the periods of 380 yr and 57 yr were used) yields a far more regular behavior of solar activity. A comparison between the results for the 10Be data with the 14C data shows many similarities. Both cosmogenic isotopes are strongly correlated mutually and with solar activity. Finally, we find that a series of damped random walk models provides a good fit to the 10Be data with a fixed characteristic time scale of 1000 years, which is roughly consistent with the quasi-periods found by the Fourier and wavelet analyses.

The Wide-field Infrared Survey Explorer (WISE) has scanned the entire sky with unprecedented sensitivity in four infrared bands, at 3.4, 4.6, 12, and 22 micron. The WISE Point Source Catalog contains more than 560 million objects, among them hundreds of thousands of galaxies with Active Nuclei (AGN). While type 1 AGN, owing to their bright and unobscured nature, are easy to detect and constitute a rather complete and unbiased sample, their type 2 counterparts, postulated by AGN unification, are not as straightforward to identify. Matching the WISE catalog with known QSOs in the Sloan Digital Sky Survey we confirm previous identification of the type 1 locus in the WISE color space. Using a very large database of the popular CLUMPY torus models, we find the colors of the putative type 2 counterparts, and also, for the first time, predict their number vs. flux relation that can be expected to be observed in any given WISE color range. This will allow us to put statistically very significant constraints on the torus parameters. Our results are a successful test of the AGN unification scheme.

Over the last decade, quasar sample sizes have increased from several thousand to several hundred thousand, thanks mostly to SDSS imaging and spectroscopic surveys. LSST, the next-generation optical imaging survey, will provide hundreds of detections per object for a sample of more than ten million quasars with redshifts of up to about seven. We briefly review optical quasar selection techniques, with emphasis on methods based on colors, variability properties and astrometric behavior.

A damped random walk is a stochastic process, defined by an exponential covariance matrix that behaves as a random walk for short time scales and asymptotically achieves a finite variability amplitude at long time scales. Over the last few years, it has been demonstrated, mostly but not exclusively using SDSS data, that a damped random walk model provides a satisfactory statistical description of observed quasar variability in the optical wavelength range, for rest-frame timescales from 5 days to 2000 days. The best-fit characteristic timescale and asymptotic variability amplitude scale with the luminosity, black hole mass, and rest wavelength, and appear independent of redshift. In addition to providing insights into the physics of quasar variability, the best-fit model parameters can be used to efficiently separate quasars from stars in imaging surveys with adequate long-term multi-epoch data, such as expected from LSST.

The science goals for ground-based large-area surveys, such as the Dark Energy Survey, Pan-STARRS, and the Large Synoptic Survey Telescope, require calibration of broadband photometry that is stable in time and uniform over the sky to precisions of a per cent or better. This performance will need to be achieved with data taken over the course of many years, and often in less than ideal conditions. This paper describes a strategy to achieve precise internal calibration of imaging survey data taken in less than photometric conditions, and reports results of an observational study of the techniques needed to implement this strategy. We find that images of celestial fields used in this case study with stellar densities of order one per arcmin-squared and taken through cloudless skies can be calibrated with relative precision of 0.5 per cent (reproducibility). We report measurements of spatial structure functions of cloud absorption observed over a range of atmospheric conditions, and find it possible to achieve photometric measurements that are reproducible to 1 per cent in images that were taken through cloud layers that transmit as little as 25 per cent of the incident optical flux (1.5 magnitudes of extinction). We find, however, that photometric precision below 1 per cent is impeded by the thinnest detectable cloud layers. We comment on implications of these results for the observing strategies of future surveys.

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

Studies of stellar populations, understood to mean collections of stars with common spatial, kinematic, chemical, and/or age distributions, have been reinvigorated during the last decade by the advent of large-area sky surveys such as SDSS, 2MASS, RAVE, and others. We review recent analyses of these data that, together with theoretical and modeling advances, are revolutionizing our understanding of the nature of the Milky Way, and galaxy formation and evolution in general. The formation of galaxies like the Milky Way was long thought to be a steady process leading to a smooth distribution of stars. However, the abundance of substructure in the multi-dimensional space of various observables, such as position, kinematics, and metallicity, is by now proven beyond doubt, and demonstrates the importance of mergers in the growth of galaxies. Unlike smooth models that involve simple components, the new data reviewed here clearly show many irregular structures, such as the Sagittarius dwarf tidal stream and the Virgo and Pisces overdensities in the halo, and the Monoceros stream closer to the Galactic plane. These recent developments have made it clear that the Milky Way is a complex and dynamical structure, one that is still being shaped by the merging of neighboring smaller galaxies. We also briefly discuss the next generation of wide-field sky surveys, such as SkyMapper, Pan-STARRS, Gaia and LSST, which will improve measurement precision manyfold, and comprise billions of individual stars. The ultimate goal, development of a coherent and detailed story of the assembly and evolutionary history of the Milky Way and other large spirals like it, now appears well within reach.

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 investigate the use of optical photometric variability to select and identify blazars in large-scale time-domain surveys, in part to aid in the identification of blazar counterparts to the ~30% of gamma-ray sources in the Fermi 2FGL catalog still lacking reliable associations. Using data from the optical LINEAR asteroid survey, we characterize the optical variability of blazars by fitting a damped random walk model to individual light curves with two main model parameters, the characteristic timescales of variability (tau), and driving amplitudes on short timescales (sigma). Imposing cuts on minimum tau and sigma allows for blazar selection with high efficiency E and completeness C. To test the efficacy of this approach, we apply this method to optically variable LINEAR objects that fall within the several-arcminute error ellipses of gamma-ray sources in the Fermi 2FGL catalog. Despite the extreme stellar contamination at the shallow depth of the LINEAR survey, we are able to recover previously-associated optical counterparts to Fermi AGN with E > 88% and C = 88% in Fermi 95% confidence error ellipses having semimajor axis r < 8'. We find that the suggested radio counterpart to Fermi source 2FGL J1649.6+5238 has optical variability consistent with other gamma-ray blazars, and is likely to be the gamma-ray source. Our results suggest that the variability of the non-thermal jet emission in blazars is stochastic in nature, with unique variability properties due to the effects of relativistic beaming. After correcting for beaming, we estimate that the characteristic timescale of blazar variability is ~3 years in the rest-frame of the jet, in contrast with the ~320 day disk flux timescale observed in quasars. The variability-based selection method presented will be useful for blazar identification in time-domain optical surveys, and is also a probe of jet physics.

We search for evidence of dark matter in the Milky Way by utilizing the stellar number density distribution and kinematics measured by the Sloan Digital Sky Survey (SDSS) to heliocentric distances exceeding ~10 kpc. We employ the cylindrically symmetric form of Jeans equations and focus on the morphology of the resulting acceleration maps, rather than the normalization of the total mass as done in previous, mostly local, studies. Jeans equations are first applied to a mock catalog based on a cosmologically derived N-body + SPH simulation, and the known acceleration (gradient of gravitational potential) is successfully recovered. The same simulation is also used to quantify the impact of dark matter on the total acceleration. We use Galfast, a code designed to quantitatively reproduce SDSS measurements and selection effects, to generate a synthetic stellar catalog. We apply Jeans equations to this catalog and produce two-dimensional maps of stellar acceleration. These maps reveal that in a Newtonian framework, the implied gravitational potential cannot be explained by visible matter alone. The acceleration experienced by stars at galactocentric distances of ~20 kpc is three times larger than what can be explained by purely visible matter. The application of an analytic method for estimating the dark matter halo axis ratio to SDSS data implies an oblate halo with q_DM = 0.47 +/- 0.14 within the same distance range. These techniques can be used to map the dark matter halo to much larger distances from the Galactic center using upcoming deep optical surveys, such as LSST.

We investigate the dichotomy in the radio loudness distribution of quasars by modelling their radio emission and various selection effects using a Monte Carlo approach. The existence of two physically distinct quasar populations, the radio-loud and radio-quiet quasars, is controversial and over the last decade a bimodal distribution of radio loudness of quasars has been both affirmed and disputed. We model the quasar radio luminosity distribution with simple unimodal and bimodal distribution functions. The resulting simulated samples are compared to a fiducial sample of 8,300 quasars drawn from the SDSS DR7 Quasar Catalog and combined with radio observations from the FIRST survey. Our results indicate that the SDSS-FIRST sample is best described by a radio loudness distribution which consists of two components, with 12+/-1 % of sources in the radio-loud component. On the other hand, the evidence for a local minimum in the loudness distribution (bimodality) is not strong and we find that previous claims for its existence were probably affected by the incompleteness of the FIRST survey close to its faint limit. We also investigate the redshift and luminosity dependence of the radio loudness distribution and find tentative evidence that at high redshift radio-loud quasars were rarer, on average "louder", and exhibited a smaller range in radio loudness. In agreement with other recent work, we conclude that the SDSS-FIRST sample strongly suggests that the radio loudness distribution of quasars is not a universal function, and that more complex models than presented here are needed to fully explain available observations.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the original charge was to produce an interim design reference mission by mid-2011. That document was delivered to NASA and widely circulated within the astronomical community. In late 2011 the Astrophysics Division augmented its original charge, asking for two design reference missions. The first of these, DRM1, was to be a finalized version of the interim DRM, reducing overall mission costs where possible. The second of these, DRM2, was to identify and eliminate capabilities that overlapped with those of NASA's James Webb Space Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's ground-based Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall mission cost, while staying faithful to NWNH. This report presents both DRM1 and DRM2.

We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along Stripe 82 (S82), whereas the remaining ~25,000 have 2-3 observations due to scan overlaps. The observed time lags span the range from a day to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 4 years, and at rest-frame wavelengths from 1000A to 6000A. We publicly release a user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Survey, and we use it to analyze the ensemble properties of quasar variability. Based on a damped random walk (DRW) model defined by a characteristic time scale and an asymptotic variability amplitude that scale with the luminosity, black hole mass, and rest wavelength for individual quasars calibrated in S82, we can fully explain the ensemble variability statistics of the non-S82 quasars such as the exponential distribution of large magnitude changes. All available data are consistent with the DRW model as a viable description of the optical continuum variability of quasars on time scales of ~5-2000 days in the rest frame. We use these models to predict the incidence of quasar contamination in transient surveys such as those from PTF and LSST.

Broad Absorption Line (BAL) trough variability is predominantly due to cloud motion transverse to our line of sight. The rate at which the variability occurs indicates the velocity of the cloud, which can provide constraints on the cloud's distance from the central source. This requires detailed spectroscopy during a variability event. Such spectra have proven elusive, suggesting either the timescale of variability is too short to be caught, or too long to notice until a sufficient amount of time has passed. Photometric monitoring of BAL quasar colours may potentially be used as an early warning system to trigger time resolved spectroscopic monitoring of BAL variability. Towards this end, we are analyzing both BAL and non-BAL colour variability using time series photometry from Stripe 82 in the Sloan Digital Sky Survey.

We use the Eighth Data Release of Sloan Digital Sky Survey (SDSS DR8) catalog with its additional sky coverage of the southern Galactic hemisphere, to measure the extent and study the nature of the Virgo Overdensity (VOD; Juric et al. 2008). The data show that the VOD extends over no less than 2000 deg^2, with its true extent likely closer to 3000 deg^2. We test whether the VOD can be attributed to a tilt in the stellar halo ellipsoid with respect to the plane of the Galactic disk and find that the observed symmetry of the north-south Galactic hemisphere star counts excludes this possibility. We argue that the Virgo Overdensity, in spite of its wide area and cloud-like appearance, is still best explained by a minor merger. Its appearance and position is qualitatively similar to a near perigalacticon merger event and, assuming that the VOD and the Virgo Stellar Stream share the same progenitor, consistent with the VSS orbit determined by Casetti-Dinescu et al. (2009).

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u,g,r,i,z) photometry is used to determine the comets' colors, sizes, surface brightness profiles, and rates of dust production in terms of the Af{\rho} formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(< H) \propto 10(0.49\pm0.05)H for H < 18, with evidence of a much shallower fit N(< H) \propto 10(0.19\pm0.03)H for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 \pm 0.05 in g - r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets.

Two upcoming large scale surveys, the ESA Gaia and LSST projects, will bring a new era in astronomy. The number of binary systems that will be observed and detected by these projects is enormous, estimations range from millions for Gaia to several tens of millions for LSST. We review some tools that should be developed and also what can be gained from these missions on the subject of binaries and exoplanets from the astrometry, photometry, radial velocity and their alert systems.

We present details of the construction and characterization of the coaddition of the Sloan Digital Sky Survey Stripe 82 \ugriz\ imaging data. This survey consists of 275 deg$^2$ of repeated scanning by the SDSS camera of $2.5\arcdeg$ of $\delta$ over $-50\arcdeg \le \alpha \le 60\arcdeg$ centered on the Celestial Equator. Each piece of sky has $\sim 20$ runs contributing and thus reaches $\sim2$ magnitudes fainter than the SDSS single pass data, i.e. to $r\sim 23.5$ for galaxies. We discuss the image processing of the coaddition, the modeling of the PSF, the calibration, and the production of standard SDSS catalogs. The data have $r$-band median seeing of 1.1\arcsec, and are calibrated to $\le 1%$. Star color-color, number counts, and psf size vs modelled size plots show the modelling of the PSF is good enough for precision 5-band photometry. Structure in the psf-model vs magnitude plot show minor psf mis-modelling that leads to a region where stars are being mis-classified as galaxies, and this is verified using VVDS spectroscopy. As this is a wide area deep survey there are a variety of uses for the data, including galactic structure, photometric redshift computation, cluster finding and cross wavelength measurements, weak lensing cluster mass calibrations, and cosmic shear measurements.

Carollo et al. have recently resolved the stellar population of the Milky Way halo into at least two distinct components, an inner halo and an outer halo. This result has been criticized by Schoenrich et al., who claim that the retrograde signature associated with the outer halo is due to the adoption of faulty distances. We refute this claim, and demonstrate that the Schoenrich et al. photometric distances are themselves flawed because they adopted an incorrect main-sequence absolute magnitude relationship from the work of Ivezi\'c et al. When compared to the recommended relation from Ivezi\'c et al., which is tied to a Milky Way globular cluster distance scale and accounts for age and metallicity effects, the relation adopted by Schoenrich et al. yields up to 18% shorter distances for stars near the main-sequence turnoff (TO). Use of the correct relationship yields agreement between the distances assigned by Carollo et al. and Ivezi\'{c} et al. for low-metallicity dwarfs to within 6-10%. Schoenrich et al. also point out that intermediate-gravity stars (3.5 <= log g <= 4.0) with colors redder than the TO region are likely misclassified, with which we concur. We implement a new procedure to reassign luminosity classifications for the TO stars that require it. New derivations of the rotational behavior demonstrate that the retrograde signature and high velocity dispersion of the outer-halo population remains. We summarize additional lines of evidence for a dual halo, including a test of the retrograde signature based on proper motions alone, and conclude that the preponderance of evidence strongly rejects the single-halo interpretation.

The VAO (Virtual Astronomical Observatory) Science Council (VAO-SC) met on July 27-28, 2011 at the Harvard-Smithsonian Center for Astrophysics in Cambridge MA, to review the VAO performance during its first year of operations. In this meeting the VAO demonstrated the new tools for astronomers that are being released in September 2011 and presented plans for the second year of activities, resulting from studies conducted during the first year. This document contains the recommendations of the VAO-SC for the second year of activity of the VAO.

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.

We employ measurements of the [alpha/Fe] ratio derived from low-resolution (R~2000) spectra of 17,277 G-type dwarfs from the SEGUE survey to separate them into likely thin- and thick-disk subsamples. Both subsamples exhibit strong gradients of orbital rotational velocity with metallicity, of opposite signs, -20 to -30 km/s/dex for the thin-disk and +40 to +50 km/s/dex for the thick-disk population. The rotational velocity is uncorrelated with Galactocentric distance for the thin-disk subsample, and exhibits a small trend for the thick-disk subsample. The rotational velocity decreases with distance from the plane for both disk components, with similar slopes (-9.0 {\pm} 1.0 km/s/kpc). Thick-disk stars exhibit a strong trend of orbital eccentricity with metallicity (about -0.2/dex), while the eccentricity does not change with metallicity for the thin-disk subsample. The eccentricity is almost independent of Galactocentric radius for the thin-disk population, while a marginal gradient of the eccentricity with radius exists for the thick-disk population. Both subsamples possess similar positive gradients of eccentricity with distance from the Galactic plane. The shapes of the eccentricity distributions for the thin- and thick-disk populations are independent of distance from the plane, and include no significant numbers of stars with eccentricity above 0.6. Among several contemporary models of disk evolution we consider, radial migration appears to have played an important role in the evolution of the thin-disk population, but possibly less so for the thick disk, relative to the gas-rich merger or disk heating scenarios. We emphasize that more physically realistic models and simulations need to be constructed in order to carry out the detailed quantitative comparisons that our new data enable.

The luminosity changes of most types of variable stars are correlated in the different wavelengths, and these correlations may be exploited for several purposes: for variability detection, for distinction of microvariability from noise, for period search or for classification. Principal component analysis is a simple and well-developed statistical tool to analyze correlated data. We will discuss its use on variable objects of Stripe 82 of the Sloan Digital Sky Survey, with the aim of identifying new RR Lyrae and SX Phoenicis-type candidates. The application is not straightforward because of different noise levels in the different bands, the presence of outliers that can be confused with real extreme observations, under- or overestimated errors and the dependence of errors on the magnitudes. These particularities require robust methods to be applied together with the principal component analysis. The results show that PCA is a valuable aid in variability analysis with multi-band data.

We use Canada-France-Hawaii Telescope Legacy Survey data for 170 deg^2, recalibrated and transformed to the Sloan Digital Sky Survey ugri photometric system, to study the distribution of near-turnoff main-sequence stars in the Galactic halo along four lines of sight to heliocentric distances of ~35 kpc. We find that the halo stellar number density profile becomes steeper at Galactocentric distances greater than R_{gal}~28 kpc, with the power law index changing from n_{inner}=-2.62+-0.04 to n_{outer}=-3.8+-0.1. In particular, we test a series of single power law models and find them to be strongly disfavored by the data. The parameters for the best-fit Einasto profile are n=2.2+-0.2 and R_e=22.2+-0.4 kpc. We measure the oblateness of the halo to be q=c/a=0.70+-0.01 and detect no evidence of it changing across the range of probed distances. The Sagittarius stream is detected in the l=173 and b=-62 direction as an overdensity of [Fe/H] ~ -1.5 dex stars at R_{gal}~32 kpc, providing a new constraint for the Sagittarius stream and dark matter halo models. We also detect the Monoceros stream as an overdensity of [Fe/H] > -1.5 dex stars in the l=232 and b=26 direction at R_{gal}<25 kpc. In the two sightlines where we do not detect significant substructure, the median metallicity is found to be independent of distance within systematic uncertainties ([Fe/H] ~ -1.5+-0.1 dex).

The advent of deep, wide, accurate, digital photometric surveys exemplified by the Sloan Digital Sky Survey (SDSS) has had a profound impact on studies of the Milky Way. In the past decade, we have transitioned from a scarcity to an (over)abundance of precise, well calibrated, observations of stars over a large fraction of the Galaxy. The avalanche of data will continue throughout this decade, culminating with Gaia and LSST. This new reality will necessitate changes in methodology, habits, and expectations both on the side of the large survey projects as well as the astrophysics community at large. We argue, based on the experience with SDSS, that surveys should release data as early and often as possible incorporating incremental improvements in each subsequent release, as opposed to holding off for a single, big, final release. The scientific community will need to reciprocate by performing analyses and (re-analyses) appropriate to the current fidelity of the released data, understanding that these are continually evolving and improving products.

We present Spitzer infrared spectra and ultra-violet to mid-infrared spectral energy distributions (SEDs) of 25 luminous type 1 quasars at z \sim 2. In general, the spectra show a bump peaking around 3 {\mu}m, and the 10 {\mu}m silicate emission feature. The 3 {\mu}m emission is identified with hot dust emission at its sublimation temperature. We explore two approaches to modeling the SED: (i) using the Clumpy model SED from Nenkova et al. (2008a), and (ii) the Clumpy model SED, and an additional blackbody component to represent the 3 {\mu}m emission. In the first case, a parameter search of \sim 1.25 million Clumpy models shows: (i) if we ignore the UV-to-near-IR SED, models fit the 2-8 {\mu}m region well, but not the 10 {\mu}m feature; (ii) if we include the UV-to-near-IR SED in the fit, models do not fit the 2-8 {\mu}m region. The observed 10 {\mu}m features are broader and shallower than those in the best-fit models in the first approach. In the second case, the shape of the 10 {\mu}m feature is better reproduced by the Clumpy models. The additional blackbody contribution in the 2-8 {\mu}m range allows Clumpy models dominated by cooler temperatures (T < 800K) to better fit the 8-12{\mu}m SED. A centrally concentrated distribution of a small number of torus clouds is required in the first case, while in the second case the clouds are more spread out radially. The temperature of the blackbody component is ~ 1200 K as expected for graphite grains.

Periodic light curve behavior predicted for some binary black hole systems might be detected in large samples, such as the multi-million quasar sample expected from the Large Synoptic Survey Telescope (LSST). We investigate the false-alarm probability for the discovery of a periodic signal in light curves simulated using damped random walk (DRW) model. This model provides a good description of observed light curves, and does not include periodic behavior. We used the Lomb-Scargle periodogram to search for a periodic signal in a million simulated light curves that properly sample the DRW parameter space, and the LSST cadence space. We find that even a very conservative threshold for the false-alarm probability still yields thousands of "good" binary black hole candidates. We conclude that the future claims for binary black holes based on Lomb-Scargle analysis of LSST light curves will have to be interpreted with caution.

We develop a method for separating quasars from other variable point sources using SDSS Stripe 82 light curve data for ~10,000 variable objects. To statistically describe quasar variability, we use a damped random walk model parametrized by a damping time scale, tau, and an asymptotic amplitude (structure function), SF_inf. With the aid of an SDSS spectroscopically confirmed quasar sample, we demonstrate that variability selection in typical extragalactic fields with low stellar density can deliver complete samples with reasonable purity (or efficiency, E). Compared to a selection method based solely on the slope of the structure function, the inclusion of the tau information boosts E from 60% to 75% while maintaining a highly complete sample (98%) even in the absence of color information. For a completeness of C=90%, E is boosted from 80% to 85%. Conversely, C improves from 90% to 97% while maintaining E=80% when imposing a lower limit on tau. With the aid of color selection, the purity can be further boosted to 96%, with C= 93%. Hence, selection methods based on variability will play an important role in the selection of quasars with data provided by upcoming large sky surveys, such as Pan-STARRS and the Large Synoptic Survey Telescope (LSST). For a typical (simulated) LSST cadence over 10 years and a photometric accuracy of 0.03 mag (achieved at i~22), C is expected to be 88% for a simple sample selection criterion of tau>100 days. In summary, given an adequate survey cadence, photometric variability provides an even better method than color selection for separating quasars from stars.

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.

The separation of the Milky Way disk into a thin and thick component is supported by differences in the spatial, kinematic and metallicity distributions of their stars. These differences have led to the view that the thick disk formed early via a cataclysmic event and constitutes fossil evidence of the hierarchical growth of the Milky Way. We show here, using N-body simulations, how a double-exponential vertical structure, with stellar populations displaying similar dichotomies can arise purely through internal evolution. In this picture, stars migrate radially, while retaining nearly circular orbits, as described by Sellwood & Binney (2002). As stars move outwards they populate a thickened component. Such stars found at the present time in the solar neighborhood formed early in the disk's history at smaller radii where stars are more metal-poor and alpha-enhanced, leading to the properties observed for thick disk stars. Classifying stars as members of the thin or thick disk by either velocity or metallicity leads to an apparent separation in the other property as observed. This scenario is supported by the SDSS observation that stars in the transition region do not show any correlation between rotational velocity and metallicity. The good qualitative agreement between our simulation and observations in the Milky Way hint that the thick disk may be a ubiquitous galaxy feature generated by stellar migration. Nonetheless, we cannot exclude that some fraction of the thick disk is a fossil of a past more violent history, nor can this scenario explain thick disks in all galaxies.

Recommendations of the VAO-Science Council following the meeting of March 26-27, 2010. Refer to the webpage http://www.aui.edu/vao.php

We report spectroscopic observations with the Gemini South Telescope of 5 faint V~20 RR Lyrae stars associated with the Pisces overdensity. At a heliocentric and galactocentric distance of ~80 kpc, this is the most distant substructure in the Galactic halo known to date. We combined our observations with literature data and confirmed that the substructure is composed of two different kinematic groups. The main group contains 8 stars and has <V_{gsr}> = 50 km/s, while the second group contains four stars at a velocity of <V_{gsr}> = -52 km/s, where V_{gsr} is the radial velocity in the galactocentric standard of rest. The metallicity distribution of RR Lyrae stars in the Pisces overdensity is centered on [Fe/H]=-1.5 dex and has a width of 0.3 dex. The new data allowed us to establish that both groups are spatially extended making it very unlikely that they are bound systems, and are more likely to be debris of a tidally disrupted galaxy or galaxies. Due to small sky coverage, it is still unclear whether these groups have the same or different progenitors.

We report on new measurements of the luminosity function (LF) and mass function (MF) of field low-mass dwarfs derived from Sloan Digital Sky Survey (SDSS) Data Release 6 (DR6) photometry. The analysis incorporates ~15 million low-mass stars (0.1 Msun < M < 0.8 Msun), spread over 8,400 square degrees. Stellar distances are estimated using new photometric parallax relations, constructed from ugriz photometry of nearby low-mass stars with trigonometric parallaxes. We use a technique that simultaneously measures Galactic structure and the stellar LF from 7 < M_r < 16. We compare the LF to previous studies and convert to a MF using the mass-luminosity relations of Delfosse et al., 2000. The system MF, measured over -1.0 < log M/Msun < -0.1, is well-described by a log-normal distribution with Mo = 0.25 Msun. We stress that our results should not be extrapolated to other mass regimes. Our work generally agrees with prior low-mass stellar MFs and places strong constraints on future star-formation studies of the Milky Way.

We present the fifth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog, which is based upon the SDSS Seventh Data Release. The catalog, which contains 105,783 spectroscopically confirmed quasars, represents the conclusion of the SDSS-I and SDSS-II quasar survey. The catalog consists of the SDSS objects that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s or have interesting/complex absorption features, are fainter than i > 15.0 and have highly reliable redshifts. The catalog covers an area of 9380 deg^2. The quasar redshifts range from 0.065 to 5.46, with a median value of 1.49; the catalog includes 1248 quasars at redshifts greater than four, of which 56 are at redshifts greater than five. The catalog contains 9210 quasars with i < 18; slightly over half of the entries have i< 19. For each object the catalog presents positions accurate to better than 0.1" rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800-9200 Ang. at a spectral resolution R = 2000 the spectra can be retrieved from the SDSS public database using the information provided in the catalog. Over 96% of the objects in the catalog were discovered by the SDSS. We also include a supplemental list of an additional 207 quasars with SDSS spectra whose archive photometric information is incomplete.

The monochromatic illumination system is constructed to carry out in situ measurements of the response function of the mosaicked CCD imager used in the Sloan Digital Sky Survey (SDSS). The system is outlined and the results of the measurements, mostly during the first 6 years of the SDSS, are described. We present the reference response functions for the five colour passbands derived from these measurements, and discuss column to column variations and variations in time, and also their effects on photometry. We also discuss the effect arising from various, slightly different response functions of the associated detector systems that were used to give SDSS photometry. We show that the calibration procedures of SDSS remove these variations reasonably well with the resulting final errors from variant response functions being unlikely to be larger than 0.01 mag for g, r, i, and z bands over the entire duration of the survey. The considerable aging effect is uncovered in the u band, the response function showing a 30% decrease in the throughput in the short wavelength side during the survey years, which potentially causes a systematic error in photometry. The aging effect is consistent with variation of the instrumental sensitivity in u-band, which is calibrated out. The expected colour variation is consistent with measured colour variation in the catalog of repeated photometry. The colour variation is delta (u-g) ~ 0.01 for most stars, and at most delta (u-g) ~ 0.02 mag for those with extreme colours. We verified in the final catalogue that no systematic variations in excess of 0.01 mag are detected in the photometry which can be ascribed to aging and/or seasonal effects except for the secular u-g colour variation for stars with extreme colours.

The structure and kinematics of the recognized stellar components of the Milky Way are explored, based on well-determined atmospheric parameters and kinematic quantities for 32360 calibration stars from the Sloan Digital Sky Survey (SDSS) and its first extension, (SDSS-II), which included the sub-survey SEGUE: Sloan Extension for Galactic Understanding and Exploration. Full space motions for a sub-sample of 16920 stars, exploring a local volume within 4 kpc of the Sun, are used to derive velocity ellipsoids for the inner- and outer-halo components of the Galaxy, as well as for the canonical thick-disk and proposed metal-weak thick-disk populations. We first examine the question of whether the data require the presence of at least a two-component halo in order to account for the rotational behavior of likely halo stars in the local volume, and whether more than two components are needed. We also address the question of whether the proposed metal-weak thick disk is kinematically and chemically distinct from the canonical thick disk. In addition, we consider the fractions of each component required to understand the nature of the observed kinematic behavior of the stellar populations of the Galaxy as a function of distance from the plane. Scale lengths and scale heights for the thick-disk and metal-weak thick-disk components are determined. Spatial density profiles for the inner- and outer-halo populations are inferred from a Jeans Theorem analysis. The full set of calibration stars (including those outside the local volume) is used to test for the expected changes in the observed stellar metallicity distribution function with distance above the Galactic plane in-situ, due to the changing contributions from the underlying stellar populations. [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 compare observations of AGB stars and predictions of the Elitzur & Ivezic (2001) steady-state radiatively driven dusty wind model. The model results are described by a set of similarity functions of a single independent variable, and imply general scaling relations among the system parameters. We find that the model properly reproduces various correlations among the observed quantities and demonstrate that dust drift through the gas has a major impact on the structure of most winds. From data for nearby oxygen-rich and carbon-rich mass-losing stars we find that (1) the dispersion in grain properties within each group is rather small; (2) both the dust cross-section per gas particle and the dust-to-gas mass ratio are similar for the two samples even though the stellar atmospheres and grain properties are very different; (3) the dust abundance in both outflows is significantly below the Galactic average, indicating that most of the Galactic dust is not stardust - contrary to popular belief, but in support of Draine (2009). Our model results can be easily applied to recent massive data sets, such as the Spitzer SAGE survey of the Large Magellanic Cloud, and incorporated in galaxy evolution models.

We study Milky Way kinematics using a sample of 18.8 million main-sequence stars with r<20 and proper-motion measurements derived from SDSS and POSS astrometry, including ~170,000 stars with radial-velocity measurements from the SDSS spectroscopic survey. Distances to stars are determined using a photometric parallax relation, covering a distance range from ~100 pc to 10 kpc over a quarter of the sky at high Galactic latitudes (|b|>20 degrees). We find that in the region defined by 1 kpc <Z< 5 kpc and 3 kpc <R< 13 kpc, the rotational velocity for disk stars smoothly decreases, and all three components of the velocity dispersion increase, with distance from the Galactic plane. In contrast, the velocity ellipsoid for halo stars is aligned with a spherical coordinate system and appears to be spatially invariant within the probed volume. The velocity distribution of nearby ($Z<1$ kpc) K/M stars is complex, and cannot be described by a standard Schwarzschild ellipsoid. For stars in a distance-limited subsample of stars (<100 pc), we detect a multimodal velocity distribution consistent with that seen by HIPPARCOS. This strong non-Gaussianity significantly affects the measurements of the velocity ellipsoid tilt and vertex deviation when using the Schwarzschild approximation. We develop and test a simple descriptive model for the overall kinematic behavior that captures these features over most of the probed volume, and can be used to search for substructure in kinematic and metallicity space. We use this model to predict further improvements in kinematic mapping of the Galaxy expected from Gaia and LSST.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

The spectroscopic Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample represents the final set of galaxies observed using the original SDSS target selection criteria. We analyse the clustering of galaxies within this sample, including both the Luminous Red Galaxy (LRG) and Main samples, and also include the 2-degree Field Galaxy Redshift Survey (2dFGRS) data. Baryon Acoustic Oscillations are observed in power spectra measured for different slices in redshift; this allows us to constrain the distance--redshift relation at multiple epochs. We achieve a distance measure at redshift z=0.275, of r_s(z_d)/D_V(0.275)=0.1390+/-0.0037 (2.7% accuracy), where r_s(z_d) is the comoving sound horizon at the baryon drag epoch, D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3), D_A(z) is the angular diameter distance and H(z) is the Hubble parameter. We find an almost independent constraint on the ratio of distances D_V(0.35)/D_V(0.2)=1.736+/-0.065, which is consistent at the 1.1sigma level with the best fit Lambda-CDM model obtained when combining our z=0.275 distance constraint with the WMAP 5-year data. The offset is similar to that found in previous analyses of the SDSS DR5 sample, but the discrepancy is now of lower significance, a change caused by a revised error analysis and a change in the methodology adopted, as well as the addition of more data. Using WMAP5 constraints on Omega_bh^2 and Omega_ch^2, and combining our BAO distance measurements with those from the Union Supernova sample, places a tight constraint on Omega_m=0.286+/-0.018 and H_0 = 68.2+/-2.2km/s/Mpc that is robust to allowing curvature and non-Lambda dark energy. This result is independent of the behaviour of dark energy at redshifts greater than those probed by the BAO and supernova measurements. (abridged)

We utilize color information for an HI-selected sample of 195 galaxies to explore the star formation histories and physical conditions that produce the observed colors. We show that the HI selection creates a significant offset towards bluer colors that can be explained by enhanced recent bursts of star formation. There is also no obvious color bimodality, because the HI selection restricts the sample to bluer, actively star forming systems, diminishing the importance of the red sequence. Rising star formation rates are still required to explain the colors of galaxies bluer than g-r < 0.3. We also demonstrate that the colors of the bluest galaxies in our sample are dominated by emission lines and that stellar population synthesis models alone (without emission lines) are not adequate for reproducing many of the galaxy colors. These emission lines produce large changes in the r-i colors but leave the g-r color largely unchanged. In addition, we find an increase in the dispersion of galaxy colors at low masses that may be the result of a change in the star formation process in low-mass galaxies.

We present the optical data for 195 HI-selected galaxies that fall within both the Sloan Digital Sky Survey (SDSS) and the Parkes Equatorial Survey (ES). The photometric quantities have been independently recomputed for our sample using a new photometric pipeline optimized for large galaxies, thus correcting for SDSS's limited reliability for automatic photometry of angularly large or low surface brightness (LSB) galaxies. We outline the magnitude of the uncertainty in the SDSS catalog-level photometry and derive a quantitative method for correcting the over-sky subtraction in the SDSS photometric pipeline. The main thrust of this paper is to present the ES/SDSS sample and discuss the methods behind the improved photometry, which will be used in future scientific analysis. We present the overall optical properties of the sample and briefly compare to a volume-limited, optically-selected sample. Compared to the optically-selected SDSS sample (in the similar volume), HI-selected galaxies are bluer and more luminous (fewer dwarf ellipticals and more star formation). However, compared to typical SDSS galaxy studies, which have their own selection effects, our sample is bluer, fainter and less massive.

We present an improved analysis of halo substructure traced by RR Lyrae stars in the SDSS stripe 82 region. With the addition of SDSS-II data, a revised selection method based on new ugriz light curve templates results in a sample of 483 RR Lyrae stars that is essentially free of contamination. The main result from our first study persists: the spatial distribution of halo stars at galactocentric distances 5--100 kpc is highly inhomogeneous. At least 20% of halo stars within 30 kpc from the Galactic center can be statistically associated with substructure. We present strong direct evidence, based on both RR Lyrae stars and main sequence stars, that the halo stellar number density profile significantly steepens beyond a Galactocentric distance of ~30 kpc, and a larger fraction of the stars are associated with substructure. By using a novel method that simultaneously combines data for RR Lyrae and main sequence stars, and using photometric metallicity estimates for main sequence stars derived from deep co-added u-band data, we measure the metallicity of the Sagittarius dSph tidal stream (trailing arm) towards R.A.2h-3h and Dec~0 deg to be 0.3 dex higher ([Fe/H]=-1.2) than that of surrounding halo field stars. Together with a similar result for another major halo substructure, the Monoceros stream, these results support theoretical predictions that an early forming, smooth inner halo, is metal poor compared to high surface brightness material that have been accreted onto a later-forming outer halo. The mean metallicity of stars in the outer halo that are not associated with detectable clumps may still be more metal-poor than the bulk of inner-halo stars, as has been argued from other data sets.

It is possible to reliably identify white dwarfs (WDs) without recourse to spectra, instead using photometric and astrometric measurements to distinguish them from Main Sequence stars and quasars. WDs' colours can also be used to infer their intrinsic properties (effective temperature, surface gravity, etc.), but the results obtained must be interpreted with care. The difficulties stem from the existence of a solid angle degeneracy, as revealed by a full exploration of the likelihood, although this can be masked if a simple best-fit approach is used. Conversely, this degeneracy can be broken if a Bayesian approach is adopted, as it is then possible to utilise the prior information on the surface gravities of WDs implied by spectroscopic fitting. The benefits of such an approach are particularly strong when applied to outliers, such as the candidate halo and ultra-cool WDs identified by Vidrih et al. (2007). A reanalysis of these samples confirms their results for the latter sample but suggests that that most of the halo candidates are thick disk WDs in the tails of the photometric noise distribution.

We present the power spectrum of the reconstructed halo density field derived from a sample of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Seventh Data Release (DR7). The halo power spectrum has a direct connection to the underlying dark matter power for k <= 0.2 h/Mpc, well into the quasi-linear regime. This enables us to use a factor of ~8 more modes in the cosmological analysis than an analysis with kmax = 0.1 h/Mpc, as was adopted in the SDSS team analysis of the DR4 LRG sample (Tegmark et al. 2006). The observed halo power spectrum for 0.02 < k < 0.2 h/Mpc is well-fit by our model: chi^2 = 39.6 for 40 degrees of freedom for the best fit LCDM model. We find \Omega_m h^2 * (n_s/0.96)^0.13 = 0.141^{+0.009}_{-0.012} for a power law primordial power spectrum with spectral index n_s and \Omega_b h^2 = 0.02265 fixed, consistent with CMB measurements. The halo power spectrum also constrains the ratio of the comoving sound horizon at the baryon-drag epoch to an effective distance to z=0.35: r_s/D_V(0.35) = 0.1097^{+0.0039}_{-0.0042}. Combining the halo power spectrum measurement with the WMAP 5 year results, for the flat LCDM model we find \Omega_m = 0.289 +/- 0.019 and H_0 = 69.4 +/- 1.6 km/s/Mpc. Allowing for massive neutrinos in LCDM, we find \sum m_{\nu} < 0.62 eV at the 95% confidence level. If we instead consider the effective number of relativistic species Neff as a free parameter, we find Neff = 4.8^{+1.8}_{-1.7}. Combining also with the Kowalski et al. (2008) supernova sample, we find \Omega_{tot} = 1.011 +/- 0.009 and w = -0.99 +/- 0.11 for an open cosmology with constant dark energy equation of state w.

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.

We perform an extensive test of theoretical stellar models for main-sequence stars in ugriz, using cluster fiducial sequences obtained in the previous paper of this series. We generate a set of isochrones using the Yale Rotating Evolutionary Code (YREC) with updated input physics, and derive magnitudes and colors in ugriz from MARCS model atmospheres. These models match cluster main sequences over a wide range of metallicity within the errors of the adopted cluster parameters. However, we find a large discrepancy of model colors at the lower main sequence (Teff < ~4500 K) for clusters at and above solar metallicity. We also reach similar conclusions using the theoretical isochrones of Girardi et al. and Dotter et al., but our new models are generally in better agreement with the data. Using our theoretical isochrones, we also derive main-sequence fitting distances and turn-off ages for five key globular clusters, and demonstrate the ability to derive these quantities from photometric data in the Sloan Digital Sky Survey. In particular, we exploit multiple color indices (g - r, g - i, and g - z) in the parameter estimation, which allows us to evaluate internal systematic errors. Our distance estimates, with an error of sigma(m - M) = 0.03-0.11 mag for individual clusters, are consistent with Hipparcos-based subdwarf fitting distances derived in the Johnson-Cousins or Stromgren photometric systems.

Under the unified model for active galactic nuclei (AGNs), narrow-line (Type 2) AGNs are, in fact, broad-line (Type 1) AGNs but each with a heavily obscured accretion disk. We would therefore expect the optical continuum emission from Type 2 AGN to be composed mainly of stellar light and non-variable on the time-scales of months to years. In this work we probe the spectroscopic variability of galaxies and narrow-line AGNs using the multi-epoch data in the Sloan Digital Sky Survey (SDSS) Data Release 6. The sample contains 18,435 sources for which there exist pairs of spectroscopic observations (with a maximum separation in time of ~700 days) covering a wavelength range of 3900-8900 angstrom. To obtain a reliable repeatability measurement between each spectral pair, we consider a number of techniques for spectrophotometric calibration resulting in an improved spectrophotometric calibration of a factor of two. From these data we find no obvious continuum and emission-line variability in the narrow-line AGNs on average -- the spectroscopic variability of the continuum is 0.07+/-0.26 mag in the g band and, for the emission-line ratios log10([NII]/Halpha) and log10([OIII]/Hbeta), the variability is 0.02+/-0.03 dex and 0.06+/-0.08 dex, respectively. From the continuum variability measurement we set an upper limit on the ratio between the flux of varying spectral component, presumably related to AGN activities, and that of host galaxy to be ~30%. We provide the corresponding upper limits for other spectral classes, including those from the BPT diagram, eClass galaxy classification, stars and quasars.

The core of scientific research is turning new ideas into reality. From the school science fair to the search for the secrets of dark energy, high-quality research consists of scientific investigation constrained within the scope of a well-defined project. Large or small, generously funded or just scraping by,scientific projects use time, money, and information to turn ideas into plans, plans into action, and action into results. While we, as a community, do much to educate students in the techniques of research, we do not systematically train students in the nature and organization of scientific projects or in the techniques of project management. We propose a two-pronged attack to address this issue in the next decade. First, to generate a broad base of future scientists who have a basic familiarity with the ideas of projects, we propose that the community develop standards for the content of a project design and management course in astronomy and astrophysics. Second, to train future scientists to assume leadership roles in new investigations in astronomy and astrophysics, we propose that the community develop standards for graduate programs in the area of research project leadership.

Wide-angle surveys have been an engine for new discoveries throughout the modern history of astronomy, and have been among the most highly cited and scientifically productive observing facilities in recent years. This trend is likely to continue over the next decade, as many of the most important questions in astrophysics are best tackled with massive surveys, often in synergy with each other and in tandem with the more traditional observatories. We argue that these surveys are most productive and have the greatest impact when the data from the surveys are made public in a timely manner. The rise of the "survey astronomer" is a substantial change in the demographics of our field; one of the most important challenges of the next decade is to find ways to recognize the intellectual contributions of those who work on the infrastructure of surveys (hardware, software, survey planning and operations, and databases/data distribution), and to make career paths to allow them to thrive.

Tidal disruption events provide a unique probe of quiescent black holes in the nuclei of distant galaxies. The next generation of synoptic surveys will yield a large sample of flares from the tidal disruption of stars by massive black holes that will give insights to four key science questions: 1) What is the assembly history of massive black holes in the universe? 2) Is there a population of intermediate mass black holes that are the primordial seeds of supermassive black holes? 3) How can we increase our understanding of the physics of accretion onto black holes? 4) Can we localize sources of gravitational waves from the detection of tidal disruption events around massive black holes and recoiling binary black hole mergers?

The data obtained by the recent modern sky surveys enable detailed studies of the stellar distribution in the multi-dimensional space spanned by spatial coordinates, velocity and metallicity, from the solar neighborhood all the way out to the outer Milky Way halo. While these results represent exciting observational breakthroughs, their interpretation is not simple. For example, traditional decomposition of the thin and thick disks predicts a strong correlation in metallicity and kinematics at $\sim$1 kpc from the Galactic plane; however, recent SDSS--based work has demonstrated an absence of this correlation for disk stars. Instead, the variation of the metallicity and rotational velocity distributions can be modeled using non--Gaussian functions that retain their shapes and only shift as the distance from the mid--plane increases. To fully contextualize these recent observational results, a detailed comparison with sophisticated numerical models is necessary. Modern simulations have sufficient resolution and physical detail to study the formation of stellar disks and spheroids over a large baseline of masses and cosmic ages. We discuss preliminary comparisons of various observed maps and N--body model predictions and find them encouraging. In particular, the N--body disk models of Ro\v{s}kar et al. \cite{Roskar 2008} reproduce a change of disk scale height reminiscent of thin/thick disk decomposition, as well as metallicity and rotational velocity gradients, while not inducing a correlation of the latter two quantities, in qualitative agreement with SDSS observations.

The Large Synoptic Survey Telescope (LSST) will be a large, wide-field ground-based system designed to obtain, starting in 2015, multiple images of the sky that is visible from Cerro Pachon in Northern Chile. About 90% of the observing time will be devoted to a deep-wide-fast survey mode which will observe a 20,000 deg$^2$ region about 1000 times during the anticipated 10 years of operations (distributed over six bands, $ugrizy$). Each 30-second long visit will deliver 5$\sigma$ depth for point sources of $r\sim24.5$ on average. The co-added map will be about 3 magnitudes deeper, and will include 10 billion galaxies and a similar number of stars. We discuss various measurements that will be automatically performed for these 20 billion sources, and how they can be used for classification and determination of source physical and other properties. We provide a few classification examples based on SDSS data, such as color classification of stars, color-spatial proximity search for wide-angle binary stars, orbital-color classification of asteroid families, and the recognition of main Galaxy components based on the distribution of stars in the position-metallicity-kinematics space. Guided by these examples, we anticipate that two grand classification challenges for LSST will be 1) rapid and robust classification of sources detected in difference images, and 2) {\it simultaneous} treatment of diverse astrometric and photometric time series measurements for an unprecedentedly large number of objects.

Significant progress in the description of quasar variability has been recently made by employing SDSS and POSS data. Common to most studies is a fundamental assumption that photometric observations at two epochs for a large number of quasars will reveal the same statistical properties as well-sampled light curves for individual objects. We critically test this assumption using light curves for a sample of $\sim$2,600 spectroscopically confirmed quasars observed about 50 times on average over 8 years by the SDSS stripe 82 survey. We find that the dependence of the mean structure function computed for individual quasars on luminosity, rest-frame wavelength and time is qualitatively and quantitatively similar to the behavior of the structure function derived from two-epoch observations of a much larger sample. We also reproduce the result that the variability properties of radio and X-ray selected subsamples are different. However, the scatter of the variability structure function for fixed values of luminosity, rest-frame wavelength and time is similar to the scatter induced by the variance of these quantities in the analyzed sample. Hence, our results suggest that, although the statistical properties of quasar variability inferred using two-epoch data capture some underlying physics, there is significant additional information that can be extracted from well-sampled light curves for individual objects.

We report on a new measurement of the luminosity function (LF) and mass function (MF) of field low-mass dwarfs using Sloan Digital Sky Survey (SDSS) photometry. The final catalog is composed of ~15 million low-mass stars (0.1 Msun < M < 0.8 Msun), spread over 8,400 square degrees. Distances to the stars are estimated using new photometric parallax relations, constructed from ugriz photometry of nearby low-mass stars with trigonometric parallaxes. The LF is measured with a novel technique, which simultaneously measures Galactic structure and the stellar LF. The resulting LF is compared to previous studies and converted to a MF. The MF is well-described by a log-normal distribution, with Mo = 0.27 Msun.

Using SDSS Data Release 6, we construct two independent samples of candidate stellar wide binaries selected as i) pairs of unresolved sources with angular separation in the range 3'' - 16'', ii) common proper motion pairs with 5'' - 30'' angular separation, and make them publicly available. These samples are dominated by disk stars, and we use them to constrain the shape of the main-sequence photometric parallax relation M_r(r-i) and to study the properties of wide binary systems. We estimate M_r(r-i) by searching for a relation that minimizes the difference between distance moduli of primary and secondary components of wide binary candidates. We model M_r(r-i) by a fourth degree polynomial and determine the coefficients using Markov Chain Monte Carlo fitting, independently for each sample. Aided by the derived photometric parallax relation, we construct a series of high-quality catalogs of candidate main-sequence binary stars. Using these catalogs, we study the distribution of semi-major axes of wide binaries, a, in the 2,000 < a < 47,000 AU range. We find the observations to be well described by the Opik distribution, f(a)~1/a, for a<a_{break}, where a_{break} increases roughly linearly with the height Z above the Galactic plane (a_{break}~12,300 Z[kpc]^0.7 AU). The number of wide binary systems with 100 AU < a < a_{break}, as a fraction of the total number of stars, decreases from 0.9% at Z=0.5 kpc to 0.5% at Z=3 kpc. The probability for a star to be in a wide binary system is independent of its color. Given this color, the companions of red components seem to be drawn randomly from the stellar luminosity function, while blue components have a larger blue-to-red companion ratio than expected from luminosity function.

We present photometry for globular and open cluster stars observed with the Sloan Digital Sky Survey (SDSS). In order to exploit over 100 million stellar objects with r < 22.5 mag observed by SDSS, we need to understand the characteristics of stars in the SDSS ugriz filters. While star clusters provide important calibration samples for stellar colors, the regions close to globular clusters, where the fraction of field stars is smallest, are too crowded for the standard SDSS photometric pipeline to process. To complement the SDSS imaging survey, we reduce the SDSS imaging data for crowded cluster fields using the DAOPHOT/ALLFRAME suite of programs and present photometry for 17 globular clusters and 3 open clusters in a SDSS value-added catalog. Our photometry and cluster fiducial sequences are on the native SDSS 2.5-meter ugriz photometric system, and the fiducial sequences can be directly applied to the SDSS photometry without relying upon any transformations. Model photometry for red giant branch and main-sequence stars obtained by Girardi et al. cannot be matched simultaneously to fiducial sequences; their colors differ by ~0.02-0.05 mag. Good agreement (< ~0.02 mag in colors) is found with Clem et al. empirical fiducial sequences in u'g'r'i'z' when using the transformation equations in Tucker et al.

We present the observational results of a survey designed to target and detect asteroids whose colors are similar to those of Vesta family members and thus may be considered as candidates for having a basaltic composition. Fifty basaltic candidates were selected with orbital elements that lie outside of the Vesta dynamical family. Optical and near-infrared spectra were used to assign a taxonomic type to 11 of the 50 candidates. Ten of these were spectroscopically confirmed as V-type asteroids, suggesting that most of the candidates are basaltic and can be used to constrain the distribution of basaltic material in the Main Belt. Using our catalog of V-type candidates and the success rate of the survey, we calculate unbiased size-frequency and semi-major axis distributions of V-type asteroids. These distributions, in addition to an estimate for the total mass of basaltic material, suggest that Vesta was the predominant contributor to the basaltic asteroid inventory of the Main Belt, however scattered planetesimals from the inner Solar System (a < 2.0 AU) and other partially/fully differentiated bodies likely contributed to this inventory. In particular, we infer the presence of basaltic fragments in the vicinity of asteroid 15 Eunomia, which may be derived from a differentiated parent body in the middle Main Belt (2.5 < a < 2.8). We find no asteroidal evidence for a large number of previously undiscovered basaltic asteroids, which agrees with previous theories suggesting that basaltic fragments from the ~100 differentiated parent bodies represented in meteorite collections have been "battered to bits" [Burbine, T.H., Meibom, A., Binzel, R.P., 1996. Mantle material in the Main Belt: Battered to bits? Met. & Planet. Sci. 31, 607].

Asteroid families, traditionally defined as clusters of objects in orbital parameter space, often have distinctive optical colors. We show that the separation of family members from background interlopers can be improved with the aid of SDSS colors as a qualifier for family membership. Based on an ~88,000 object subset of the Sloan Digital Sky Survey Moving Object Catalog 4 with available proper orbital elements, we define 37 statistically robust asteroid families with at least 100 members using a simple Gaussian distribution model in both orbital and color space. The interloper rejection rate based on colors is typically ~10% for a given orbital family definition, with four families that can be reliably isolated only with the aid of colors. About 50% of all objects in this data set belong to families, and this fraction varies from about 35% for objects brighter than an H magnitude of 13 and rises to 60% for objects fainter than this. The fraction of C-type objects in families decreases with increasing H magnitude for H > 13, while the fraction of S-type objects above this limit remains effectively constant. This suggests that S-type objects require a shorter timescale for equilibrating the background and family size distributions via collisional processing. The size distributions for 15 families display a well-defined change of slope and can be modeled as a "broken" double power-law. Such "broken" size distributions are twice as likely for S-type familes than for C-type families, and are dominated by dynamically old families. The remaining families with size distributions that can be modeled as a single power law are dominated by young families. When size distribution requires a double power-law model, the two slopes are correlated and are steeper for S-type families.

According to unified schemes of Active Galactic Nuclei (AGN), the central engine is surrounded by dusty, optically thick clouds in a toroidal structure. We have recently developed a formalism that for the first time takes proper account of the clumpy nature of the AGN torus. We now provide a detailed report of our findings in a two-paper series. Here we present our general formalism for radiative transfer in clumpy media and construct its building blocks for the AGN problem -- the source functions of individual dusty clouds heated by the AGN radiation field. We show that a fundamental difference from smooth density distributions is that in a clumpy medium, a large range of dust temperatures coexist at the same distance from the radiation central source. This distinct property explains the low dust temperatures found close to the nucleus of NGC1068 in 10 \mic interferometric observations. We find that irrespective of the overall geometry, a clumpy dust distribution shows only moderate variation in its spectral energy distribution, and the 10\mic\ absorption feature is never deep. Furthermore, the X-ray attenuating column density is widely scattered around the column density that characterizes the IR emission. All of these properties are characteristic of AGN observations. The assembly of clouds into AGN tori and comparison with observations is presented in the companion paper.

From extensive radiative transfer calculations we find that clumpy torus models with \No \about 5--15 dusty clouds along radial equatorial rays successfully explain AGN infrared observations. The dust has standard Galactic composition, with individual cloud optical depth \tV \about 30--100 at visual. The models naturally explain the observed behavior of the 10\mic silicate feature, in particular the lack of deep absorption features in AGN of any type. The weak 10\mic emission feature tentatively detected in type 2 QSO can be reproduced if in these sources \No drops to \about 2 or \tV exceeds \about 100. The clouds angular distribution must have a soft-edge, e.g., Gaussian profile, the radial distribution should decrease as $1/r$ or $1/r^2$. Compact tori can explain all observations, in agreement with the recent interferometric evidence that the ratio of the torus outer to inner radius is perhaps as small as \about 5--10. Clumpy torus models can produce nearly isotropic IR emission together with highly anisotropic obscuration, as required by observations. In contrast with strict variants of unification schemes where the viewing-angle uniquely determines the classification of an AGN into type 1 or 2, clumpiness implies that it is only a probabilistic effect; a source can display type 1 properties even from directions close to the equatorial plane. The fraction of obscured sources depends not only on the torus angular thickness but also on the cloud number \No. The observed decrease of this fraction at increasing luminosity can be explained with a decrease of either torus angular thickness or cloud number, but only the latter option explains also the possible emergence of a 10\mic emission feature in QSO2.

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]

Abridged: We estimate the distances to ~48 million stars detected by the Sloan Digital Sky Survey and map their 3D number density distribution in 100 < D < 20 kpc range over 6,500 deg^2 of sky. The data show strong evidence for a Galaxy consisting of an oblate halo, a disk component, and a number of localized overdensities with exponential disk parameters (bias-corrected for an assumed 35% binary fraction) H_1 = 300 pc, L_1 = 2600 pc, H_2 = 900 pc, L_2 = 3600 pc, and local density normalization of 12%. We find the halo to be oblate, with best-fit axis ratio c/a = 0.64, r^{-2.8} profile, and the local halo-to-thin disk normalization of 0.5%. We estimate the errors of derived model parameters to be no larger than ~20% (disk scales) and ~10% (thick disk normalization). While generally consistent with the above model, the density distribution shows a number of statistically significant localized deviations. We detect two overdensities in the thick disk region at (R, Z) ~ (6.5, 1.5)kpc and (R, Z) ~ (9.5, 0.8) kpc, and a remarkable density enhancement in the halo covering >1000deg^2 of sky towards the constellation of Virgo, at distances of ~6-20 kpc. Compared to a region symmetric with respect to the l=0 line, the Virgo overdensity is responsible for a factor of 2 number density excess and may be a nearby tidal stream or a low-surface brightness dwarf galaxy merging with the Milky Way. After removal of the resolved overdensities, the remaining data are consistent with a smooth density distribution; we detect no evidence of further unresolved clumpy substructure at scales ranging from ~50pc in the disk, to ~1 - 2 kpc in the halo.

We have used data from the Sloan Digital Sky Survey (SDSS) Data Release 5 to explore the overall structure and substructure of the stellar halo of the Milky Way using about 4 million color-selected main sequence turn-off stars. We fit oblate and triaxial broken power-law models to the data, and found a `best-fit' oblateness of the stellar halo 0.5<c/a<0.8, and halo stellar masses between Galactocentric radii of 1 and 40kpc of (3.7+/-1.2)x10^8 M_sun. The density profile of the stellar halo is approximately r^{-3}; it is possible that the power law slope is shallower inside 20kpc and steeper outside that radius. Yet, we found that all smooth and symmetric models were very poor fits to the distribution of stellar halo stars because the data exhibit a great deal of spatial substructure. We quantified deviations from a smooth oblate/triaxial model using the RMS of the data around the model profile on scales >~100pc, after accounting for the (known) contribution of Poisson uncertainties. The fractional RMS deviation of the actual stellar distribution from any smooth, parameterized halo model is >~40%: hence, the stellar halo is highly structured. We compared the observations with simulations of galactic stellar halos formed entirely from the accretion of satellites in a cosmological context by analysing the simulations in the same way as the data. While the masses, overall profiles, and degree of substructure in the simulated stellar halos show considerable scatter, the properties and degree of substructure in the Milky Way's halo match well the properties of a `typical' stellar halo built exclusively out of the debris from disrupted satellite galaxies. Our results therefore point towards a picture in which an important fraction of the Milky Way's stellar halo has been accreted from satellite galaxies.

We explore the properties of the submillijansky radio population at 20 cm by applying a newly developed optical color-based method to separate star forming (SF) from AGN galaxies at intermediate redshifts (z<1.3). Although optical rest-frame colors are used, our separation method is shown to be efficient, and not biased against dusty starburst galaxies. This classification method has been calibrated and tested on a local radio selected optical sample. Given accurate multi-band photometry and redshifts, it carries the potential to be generally applicable to any galaxy sample where SF and AGN galaxies are the two dominant populations. In order to quantify the properties of the submillijansky radio population, we have analyzed ~2,400 radio sources, detected at 20 cm in the VLA-COSMOS survey. 90% of these have submillijansky flux densities. We classify the objects into 1) star candidates, 2) quasi stellar objects, 3) AGN, 4) SF, and 5) high redshift (z>1.3) galaxies. We find, for the composition of the submillijansky radio population, that SF galaxies are not the dominant population at submillijansky flux levels, as previously often assumed, but that they make up an approximately constant fraction of 30-40% in the flux density range of ~50 microJy to 0.7 mJy. In summary, based on the entire VLA-COSMOS radio population at 20 cm, we find that the radio population at these flux densities is a mixture of roughly 30-40% of SF and 50-60% of AGN galaxies, with a minor contribution (~10%) of QSOs.

AM CVn systems are a select group of ultracompact binaries with the shortest orbital periods of any known binary subclass; mass-transfer is likely from a low-mass (partially-)degenerate secondary onto a white dwarf primary, driven by gravitational radiation. In the past few years, the Sloan Digital Sky Survey (SDSS) has provided five new AM CVns. Here we report on two further candidates selected from more recent SDSS data. SDSS J1208+3550 is similar to the earlier SDSS discoveries, recognized as an AM CVn via its distinctive spectrum which is dominated by helium emission. From the expanded SDSS Data Release 6 (DR6) spectroscopic area, we provide an updated surface density estimate for such AM CVns of order 10^{-3.1} to 10^{-2.5} per deg^2 for 15<g<20.5. In addition, we present another new candidate AM CVn, SDSS J2047+0008, that was discovered in the course of followup of SDSS-II supernova candidates. It shows nova-like outbursts in multi-epoch imaging data; in contrast to the other SDSS AM CVn discoveries, its (outburst) spectrum is dominated by helium absorption lines, reminiscent of KL Dra and 2003aw. The variability selection of SDSS J2047+0008 from the 300 deg^2 of SDSS Stripe 82 presages further AM CVn discoveries in future deep, multicolor, and time-domain surveys such as LSST. The new additions bring the total SDSS yield to seven AM CVns thus far, a substantial contribution to this rare subclass, versus the dozen previously known.

We present a new public archive of light-motion curves in Sloan Digital Sky Survey (SDSS) Stripe 82, covering 99 deg in right ascension from RA = 20.7 h to 3.3 h and spanning 2.52 deg in declination from Dec = -1.26 to 1.26 deg, for a total sky area of ~249 sq deg. Stripe 82 has been repeatedly monitored in the u, g, r, i and z bands over a seven-year baseline. Objects are cross-matched between runs, taking into account the effects of any proper motion. The resulting catalogue contains almost 4 million light-motion curves of stellar objects and galaxies. The photometry are recalibrated to correct for varying photometric zeropoints, achieving ~20 mmag and ~30 mmag root-mean-square (RMS) accuracy down to 18 mag in the g, r, i and z bands for point sources and extended sources, respectively. The astrometry are recalibrated to correct for inherent systematic errors in the SDSS astrometric solutions, achieving ~32 mas and ~35 mas RMS accuracy down to 18 mag for point sources and extended sources, respectively. For each light-motion curve, 229 photometric and astrometric quantities are derived and stored in a higher-level catalogue. On the photometric side, these include mean exponential and PSF magnitudes along with uncertainties, RMS scatter, chi^2 per degree of freedom, various magnitude distribution percentiles, object type (stellar or galaxy), and eclipse, Stetson and Vidrih variability indices. On the astrometric side, these quantities include mean positions, proper motions as well as their uncertainties and chi^2 per degree of freedom. The here presented light-motion curve catalogue is complete down to r~21.5 and is at present the deepest large-area photometric and astrometric variability catalogue available.

We compare environmental effects in two analogous samples of galaxies, one from the Sloan Digital Sky Survey (SDSS) and the other from a semi-analytic model (SAM) based on the Millennium Simulation (MS), to test to what extent current SAMs of galaxy formation are reproducing environmental effects. We estimate the large-scale environment of each galaxy using a Bayesian density estimator based on distances to all ten nearest neighbors and compare broad-band photometric properties of the two samples as a function of environment. The feedbacks implemented in the semi-analytic model produce a qualitatively correct galaxy population with similar environmental dependence as that seen in SDSS galaxies. In detail, however, the colors of MS galaxies exhibit an exaggerated dependence on environment: the field contains too many blue galaxies while clusters contain too many red galaxies, compared to the SDSS sample. We also find that the MS contains a population of highly clustered, relatively faint red galaxies with velocity dispersions comparable to their Hubble flow. Such high-density galaxies, if they exist, would be overlooked in any low-redshift survey since their membership to a cluster cannot be determined due to the "Fingers of God" effect.

We have developed an end-to-end photometric data processing pipeline to compare current photometric algorithms commonly used on ground-based imaging data. This testbed is exceedingly adaptable, and enables us to perform many research and development tasks, including image subtraction and co-addition, object detection and measurements, the production of photometric catalogs, and the creation and stocking of database tables with time-series information. This testing has been undertaken to evaluate existing photometry algorithms for consideration by a next-generation image processing pipeline for the Large Synoptic Survey Telescope (LSST). We outline the results of our tests for four packages: The Sloan Digital Sky Survey's (SDSS) Photo package, Daophot and Allframe, DoPhot, and two versions of Source Extractor (SExtractor). The ability of these algorithms to perform point-source photometry, astrometry, shape measurements, star-galaxy separation, and to measure objects at low signal-to-noise is quantified. We also perform a detailed crowded field comparison of Daophot and Allframe, and profile the speed and memory requirements in detail for SExtractor. We find that both Daophot and Photo are able to perform aperture photometry to high enough precision to meet LSST's science requirements, and less adequately at PSF-fitting photometry. Photo performs the best at simultaneous point and extended-source shape and brightness measurements. SExtractor is the fastest algorithm, and recent upgrades in the software yield high-quality centroid and shape measurements with little bias towards faint magnitudes. Allframe yields the best photometric results in crowded fields.

We present an algorithm to photometrically calibrate wide field optical imaging surveys, that simultaneously solves for the calibration parameters and relative stellar fluxes using overlapping observations. The algorithm decouples the problem of "relative" calibrations, from that of "absolute" calibrations; the absolute calibration is reduced to determining a few numbers for the entire survey. We pay special attention to the spatial structure of the calibration errors, allowing one to isolate particular error modes in downstream analyses. Applying this to the Sloan Digital Sky Survey imaging data, we achieve ~1% relative calibration errors across 8500 sq.deg. in griz; the errors are ~2% for the u band. These errors are dominated by unmodelled atmospheric variations at Apache Point Observatory. These calibrations, dubbed "ubercalibration", are now public with SDSS Data Release 6, and will be a part of subsequent SDSS data releases.

We report the first determination of a distance bracket for the high-velocity cloud (HVC) complex C. Combined with previous measurements showing that this cloud has a metallicity of 0.15 times solar, these results provide ample evidence that complex C traces the continuing accretion of intergalactic gas falling onto the Milky Way. Accounting for both neutral and ionized hydrogen as well as He, the distance bracket implies a mass of 3-14x10^6 M_sun, and the complex represents a mass inflow of 0.1-0.25 M_sun/yr. We base our distance bracket on the detection of CaII absorption in the spectrum of the blue horizontal branch star SDSS J120404.78+623345.6, in combination with a significant non-detection toward the BHB star BS 16034-0114. These results set a strong distance bracket of 3.7-11.2 kpc on the distance to complex C. A more weakly supported lower limit of 6.7 kpc may be derived from the spectrum of the BHB star BS 16079-0017.

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.

The high- and intermediate-velocity interstellar clouds (HVCs/IVCs) are tracers of energetic processes in and around the Milky Way. Clouds with near-solar metallicity about one kpc above the disk trace the circulation of material between disk and halo (the Galactic Fountain). The Magellanic Stream consists of gas tidally extracted from the SMC, tracing the dark matter potential of the Milky Way. Several other HVCs have low-metallicity and appear to trace the continuing accretion of infalling intergalactic gas. These assertions are supported by the metallicities (0.1 to 1 solar) measured for about ten clouds in the past decade. Direct measurements of distances to HVCs have remained elusive, however. In this paper we present four new distance brackets, using VLT observations of interstellar \CaII H and K absorption toward distant Galactic halo stars. We derive distance brackets of 5.0 to 11.7 kpc for the Cohen Stream (likely to be an infalling low-metallicity cloud), 9.8 to 15.1 kpc for complex GCP (also known as the Smith Cloud or HVC40-15+100 and with still unknown origin), 1.0 to 2.7 kpc for an IVC that appears associated with the return flow of the Fountain in the Perseus Arm, and 1.8 to 3.8 kpc for cloud g1, which appears to be in the outflow phase of the Fountain. Our measurements further demonstrate that the Milky Way is accreting substantial amounts of gaseous material, which influences the Galaxy's current and future dynamical and chemical evolution.

We explore how well crowded field point-source photometry can be accomplished with SDSS data: We present a photometric pipeline based on DoPhot, and tuned for analyzing crowded-field images from the SDSS. Using Monte Carlo simulations we show that the completeness of source extraction is above 80% to i < 21 (AB) and a stellar surface density of about 200 sq.amin. Hence, a specialized data pipeline can efficiently be used for e.g. nearby resolved galaxies in SDSS images, where the standard SDSS photometric package Photo, when applied in normal survey mode, gives poor results. We apply our pipeline to an area of about 3.55sq.deg. around the dwarf spheroidal galaxy (dSph) Leo I, and construct a high S/N star-count map of Leo I via an optimized filter in color-magnitude space (g,r,i). Although the radial surface-density profile of the dwarf deviates from the best fit empirical King model towards outer radii, we find no evidence for tidal debris out to a stellar surface-density of 4*10^(-3) of the central value. We determine the total luminosity of Leo I, and model its mass using the spherical and isotropic Jeans equation. Assuming that 'mass follows light' we constrain a lower limit of the total mass of the dSph to be (1.7+/-0.2)*10^7 Msol. Contrary, if the mass in Leo I is dominated by a constant density dark-matter (DM) halo, then the mass within the central 12' is (2+/-0.6)*10^8 Msol. This leads to a mass-to-light ratio of >>6 (Ic_sol), and possibly >75 if the DM halo dominates the mass and extends further out than 12'. In summary, our results show that Leo I is a symmetric, relaxed and bound system; this supports the idea that Leo I is a dark-matter dominated system.

We present the fourth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog. The catalog contains 77,429 objects; this is an increase of over 30,000 entries since the previous edition. The catalog consists of the objects in the SDSS Fifth Data Release that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s, or have interesting/complex absorption features, are fainter than i=15.0, and have highly reliable redshifts. The area covered by the catalog is 5740 sq. deg. The quasar redshifts range from 0.08 to 5.41, with a median value of 1.48; the catalog includes 891 quasars at redshifts greater than four, of which 36 are at redshifts greater than five. Approximately half of the catalog quasars have i < 19; nearly all have i < 21. For each object the catalog presents positions accurate to better than 0.2 arcsec. rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains basic radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800--9200A at a spectral resolution of ~2000. The spectra can be retrieved from the public database using the information provided in the catalog. The average SDSS colors of quasars as a function of redshift, derived from the catalog entries, are presented in tabular form. Approximately 96% of the objects in the catalog were discovered by the SDSS.

We quantify the variability of faint unresolved optical sources using a catalog based on multiple SDSS imaging observations. The catalog covers SDSS Stripe 82, and contains 58 million photometric observations in the SDSS ugriz system for 1.4 million unresolved sources. In each photometric bandpass we compute various low-order lightcurve statistics and use them to select and study variable sources. We find that 2% of unresolved optical sources brighter than g=20.5 appear variable at the 0.05 mag level (rms) simultaneously in the g and r bands. The majority (2/3) of these variable sources are low-redshift (<2) quasars, although they represent only 2% of all sources in the adopted flux-limited sample. We find that at least 90% of quasars are variable at the 0.03 mag level (rms) and confirm that variability is as good a method for finding low-redshift quasars as is the UV excess color selection (at high Galactic latitudes). We analyze the distribution of lightcurve skewness for quasars and find that is centered on zero. We find that about 1/4 of the variable stars are RR Lyrae stars, and that only 0.5% of stars from the main stellar locus are variable at the 0.05 mag level. The distribution of lightcurve skewness in the g-r vs. u-g color-color diagram on the main stellar locus is found to be bimodal (with one mode consistent with Algol-like behavior). Using over six hundred RR Lyrae stars, we demonstrate rich halo substructure out to distances of 100 kpc. We extrapolate these results to expected performance by the Large Synoptic Survey Telescope and estimate that it will obtain well-sampled 2% accurate, multi-color lightcurves for ~2 million low-redshift quasars, and will discover at least 50 million variable stars.

We describe a standard star catalog constructed using multiple SDSS photometric observations (at least four per band, with a median of ten) in the $ugriz$ system. The catalog includes 1.01 million non-variable unresolved objects from the equatorial stripe 82 ($|\delta_{J2000}|<$ 1.266$^\circ$) in the RA range 20h 34m to 4h 00m, and with the corresponding $r$ band (approximately Johnson V band) magnitudes in the range 14--22. The distributions of measurements for individual sources demonstrate that the photometric pipeline correctly estimates random photometric errors, which are below 0.01 mag for stars brighter than (19.5, 20.5, 20.5, 20, 18.5) in $ugriz$, respectively (about twice as good as for individual SDSS runs). Several independent tests of the internal consistency suggest that the spatial variation of photometric zeropoints is not larger than $\sim$0.01 mag (rms). In addition to being the largest available dataset with optical photometry internally consistent at the $\sim$1% level, this catalog provides practical definition of the SDSS photometric system. Using this catalog, we show that photometric zeropoints for SDSS observing runs can be calibrated within nominal uncertainty of 2% even for data obtained through 1 mag thick clouds, and demonstrate the existence of He and H white dwarf sequences using photometric data alone. Based on the properties of this catalog, we conclude that upcoming large-scale optical surveys such as the Large Synoptic Survey Telescope will be capable of delivering robust 1% photometry for billions of sources.

We analyze 1187 observations of about 860 unique candidate Jovian Trojan asteroids listed in the 3rd release of Sloan Digital Sky Survey (SDSS) Moving Object Catalog. The sample is complete at the faint end to r=21.2 mag (apparent brightness) and H=13.8 (absolute brightness, approximately corresponding to 10 km diameter). A subset of 297 detections of previously known Trojans were used to design and optimize a selection method based on observed angular velocity that resulted in the remaining objects. Using a sample of objects with known orbits, we estimate that the candidate sample contamination is about 3%. The well-controlled selection effects, the sample size, depth and accurate five-band UV-IR photometry enabled several new findings and the placement of older results on a firmer statistical footing. We find that there are significantly more asteroids in the leading swarm (L4) than in the trailing swarm (L5): N(L4)/N(L5)=1.6$\pm$0.1, independently of limiting object's size. The overall counts normalization suggests that there are about as many Jovians Trojans as there are main-belt asteroids down to the same size limit, in agreement with earlier estimates. We find that Trojan asteroids have a remarkably narrow color distribution (root-mean-scatter of only $\sim$0.05 mag) that is significantly different from the color distribution of the main-belt asteroids. The color of Trojan asteroids is correlated with their orbital inclination, in a similar way for both swarms, but appears uncorrelated with the object's size. We extrapolate the results presented here and estimate that Large Synoptic Survey Telescope will determine orbits, accurate colors and measure light curves in six photometric bandpasses for about 100,000 Jovian Trojan asteroids.

The recently discovered Virgo stellar over-density, which expands over \~1000deg^2 perpendicularly to the Galactic disk plane (7< Z <15 kpc, R~7 kpc), is the largest clump of tidal debris ever detected in the outer halo and is likely related with the accretion of a nearby dwarf galaxy by the Milky Way. We carry out N-body simulations of the Sagittarius stream to show that this giant stellar over-density is a confirmation of theoretical model predictions for the leading tail of the Sagittarius stream to cross the Milky Way plane in the Solar neighborhood. Radial velocity measurements are needed to confirm this association and to further constrain the shape of the Milky Way dark matter halo through a new generation of theoretical models. If the identification of Virgo over-density and the Sagittarius leading arm is correct, we predict highly negative radial velocities for the stars of Virgo over-density. The detection of this new portion of the Sagittarius tidal stream would represent an excellent target for the on-going and future kinematic surveys and for dark matter direct detection experiments in the proximity of the Sun.

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.

(Abridged) The Large Synoptic Survey Telescope (LSST) is currently by far the most ambitious proposed ground-based optical survey. The main science themes that drive the LSST system design are Dark Energy and Matter, the Solar System Inventory, Transient Optical Sky and the Milky Way Mapping. The LSST system, with its 8.4m telescope and 3,200 Megapixel camera, will be sited at Cerro Pachon in northern Chile, with the first light scheduled for 2013. In a continuous observing campaign, LSST will cover the entire available sky every three nights in two photometric bands to a depth of V=25 per visit (two 15 second exposures), with exquisitely accurate astrometry and photometry. Over the proposed survey lifetime of 10 years, each sky location would be observed about 1000 times, with the total exposure time of 8 hours distributed over six broad photometric bandpasses (ugrizY). This campaign will open a movie-like window on objects that change brightness, or move, on timescales ranging from 10 seconds to 10 years. The survey will have a data rate of about 30 TB/night, and will collect over 60 PB of data over its lifetime, resulting in an incredibly rich and extensive public archive that will be a treasure trove for breakthroughs in many areas of astronomy. I describe how this archive will impact the AGB star research and speculate how the system could be further optimized by utilizing narrow-band TiO and CN filters.

We discuss several results made possible by accurate SDSS astrometric measurements in a large sky area, with emphasis on asteroids and stellar proper motions obtained by comparing POSS and SDSS. SDSS has observed over 200,000 moving objects in five photometric bands, corresponding to about two orders of magnitude increase over previous multi--color surveys. These data were used to extend the measurement of asteroid size distribution to a smaller size limit, to demonstrate that asteroid dynamical families, defined as clusters in orbital parameter space, also strongly segregate in color space, and to discover a correlation between asteroid age and colors. A preliminary analysis of SDSS-POSS proper motions for about 1 million M dwarf stars demonstrates that, in the 0.1-1 kpc distance range, the rotational velocity and its dispersion for disk stars increase with the distance from the Galactic plane.

(Abridged) The Large Synoptic Survey Telescope (LSST) is currently by far the most ambitious proposed ground-based optical survey. Solar System mapping is one of the four key scientific design drivers, with emphasis on efficient Near-Earth Object (NEO) and Potentially Hazardous Asteroid (PHA) detection, orbit determination, and characterization. In a continuous observing campaign of pairs of 15 second exposures of its 3,200 megapixel camera, LSST will cover the entire available sky every three nights in two photometric bands to a depth of V=25 per visit (two exposures), with exquisitely accurate astrometry and photometry. Over the proposed survey lifetime of 10 years, each sky location would be visited about 1000 times. The baseline design satisfies strong constraints on the cadence of observations mandated by PHAs such as closely spaced pairs of observations to link different detections and short exposures to avoid trailing losses. Equally important, due to frequent repeat visits LSST will effectively provide its own follow-up to derive orbits for detected moving objects. Detailed modeling of LSST operations, incorporating real historical weather and seeing data from LSST site at Cerro Pachon, shows that LSST using its baseline design cadence could find 90% of the PHAs with diameters larger than 250 m, and 75% of those greater than 140 m within ten years. However, by optimizing sky coverage, the ongoing simulations suggest that the LSST system, with its first light in 2013, can reach the Congressional mandate of cataloging 90% of PHAs larger than 140m by 2020.

We compare Stetson's photometric standards with measurements listed in a standard star catalog constructed using repeated SDSS imaging observations. The SDSS catalog includes over 700,000 candidate standard stars from the equatorial stripe 82 (|Dec|<1.266 deg) in the RA range 20h 34' to 4h 00', and with the $r$ band magnitudes in the range 14--21. The distributions of measurements for individual sources demonstrate that the SDSS photometric pipeline correctly estimates random photometric errors, which are below 0.01 mag for stars brighter than (19.5, 20.5, 20.5, 20, 18.5) in ugriz, respectively (about twice as good as for individual SDSS runs). We derive mean photometric transformations between the SDSS gri and the BVRI system using 1165 Stetson stars found in the equatorial stripe 82, and then study the spatial variation of the difference in zeropoints between the two catalogs. Using third order polynomials to describe the color terms, we find that photometric measurements for main-sequence stars can be transformed between the two systems with systematic errors smaller than a few millimagnitudes. The spatial variation of photometric zeropoints in the two catalogs typically does not exceed 0.01 magnitude. Consequently, the SDSS Standard Star Catalog for Stripe 82 can be used to calibrate new data in both the SDSS ugriz and the BVRI systems with a similar accuracy.

We present the Chandra Multiwavelength Project (ChaMP) X-ray point source catalog with ~6,800 X-ray sources detected in 149 Chandra observations covering \~10 deg^2. The full ChaMP catalog sample is seven times larger than the initial published ChaMP catalog. The exposure time of the fields in our sample ranges from 0.9 to 124 ksec, corresponding to a deepest X-ray flux limit of f_{0.5-8.0} = 9 x 10^{-16} erg/cm2/sec. The ChaMP X-ray data have been uniformly reduced and analyzed with ChaMP-specific pipelines, and then carefully validated by visual inspection. The ChaMP catalog includes X-ray photometric data in 8 different energy bands as well as X-ray spectral hardness ratios and colors. To best utilize the ChaMP catalog, we also present the source reliability, detection probability and positional uncertainty. To quantitatively assess those parameters, we performed extensive simulations. In particular, we present a set of empirical equations: the flux limit as a function of effective exposure time, and the positional uncertainty as a function of source counts and off axis angle. The false source detection rate is ~1% of all detected ChaMP sources, while the detection probability is better than ~95% for sources with counts >30 and off axis angle <5 arcmin. The typical positional offset between ChaMP X-ray source and their SDSS optical counterparts is 0.7+-0.4 arcsec, derived from ~900 matched sources.

Using a sample of optically-selected quasars from the Sloan Digital Sky Survey, we have determined the radio-loud fraction (RLF) of quasars as a function of redshift and optical luminosity. The sample contains more than 30,000 objects and spans a redshift range of 0<z<5 and a luminosity range of -30<M_i<-22. We use both the radio-to-optical flux ratio (R parameter) and radio luminosity to define radio-loud quasars. After breaking the correlation between redshift and luminosity due to the flux-limited nature of the sample, we find that the RLF of quasars decreases with increasing redshift and decreasing luminosity. The relation can be described in the form of log(RLF/(1-RLF)) = b_0 + b_z log(1+z) + b_M (M_{2500}+26), where M_{2500} is the absolute magnitude at rest-frame 2500A, and b_z, b_M<0. When using R>10 to define radio-loud quasars, we find that b_0=-0.132+/-0.116, b_z=-2.052+/-0.261, and b_M=-0.183+/-0.025. The RLF at z=0.5 declines from 24.3% to 5.6% as luminosity decreases from M_{2500}=-26 to M_{2500}=-22, and the RLF at M_{2500}=-26 declines from 24.3% to 4.1% as redshift increases from 0.5 to 3, suggesting that the RLF is a strong function of both redshift and luminosity. We also examine the impact of flux-related selection effects on the RLF determination using a series of tests, and find that the dependence of the RLF on redshift and luminosity is highly likely to be physical, and the selection effects we considered are not responsible for the dependence.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

We describe further results of a program aimed to yield ~10^4 fully characterized optical identifications of ROSAT X-ray sources. Our program employs X-ray data from the ROSAT All-Sky Survey (RASS), and both optical imaging and spectroscopic data from the Sloan Digital Sky Survey (SDSS). RASS/SDSS data from 5740 deg^2 of sky spectroscopically covered in SDSS Data Release 5 (DR5) provide an expanded catalog of 7000 confirmed quasars and other AGN that are probable RASS identifications. Again in our expanded catalog, the identifications as X-ray sources are statistically secure, with only a few percent of the SDSS AGN likely to be randomly superposed on unrelated RASS X-ray sources. Most identifications continue to be quasars and Seyfert 1s with 15<m<21 and 0.01<z<4; but the total sample size has grown to include very substantial numbers of even quite rare AGN, e.g., now including several hundreds of candidate X-ray emitting BL Lacs and narrow-line Seyfert 1 galaxies. In addition to exploring rare subpopulations, such a large total sample may be useful when considering correlations between the X-ray and the optical, and may also serve as a resource list from which to select the "best" object (e.g., X-ray brightest AGN of a certain subclass, at a preferred redshift or luminosity) for follow-on X-ray spectral or alternate detailed studies.

The photometric calibration of the Sloan Digital Sky Survey (SDSS) is a multi-step process which involves data from three different telescopes: the 1.0-m telescope at the US Naval Observatory (USNO), Flagstaff Station, Arizona (which was used to establish the SDSS standard star network); the SDSS 0.5-m Photometric Telescope (PT) at the Apache Point Observatory (APO), New Mexico (which calculates nightly extinctions and calibrates secondary patch transfer fields); and the SDSS 2.5-m telescope at APO (which obtains the imaging data for the SDSS proper). In this paper, we describe the Monitor Telescope Pipeline, MTPIPE, the software pipeline used in processing the data from the single-CCD telescopes used in the photometric calibration of the SDSS (i.e., the USNO 1.0-m and the PT). We also describe transformation equations that convert photometry on the USNO-1.0m u'g'r'i'z' system to photometry the SDSS 2.5m ugriz system and the results of various validation tests of the MTPIPE software. Further, we discuss the semi-automated PT factory, which runs MTPIPE in the day-to-day standard SDSS operations at Fermilab. Finally, we discuss the use of MTPIPE in current SDSS-related projects, including the Southern u'g'r'i'z' Standard Star project, the u'g'r'i'z' Open Star Clusters project, and the SDSS extension (SDSS-II).

We discuss the panchromatic properties of 99,088 galaxies selected from the SDSS Data Release 1 spectroscopic sample (a flux-limited sample for 1360 deg^2). These galaxies are positionally matched to sources detected by ROSAT, GALEX, 2MASS, IRAS, GB6, FIRST, NVSS and WENSS. We find strong correlations between the detection fraction at other wavelengths and optical properties such as flux, colors, and emission-line strengths. Using GALEX, SDSS, and 2MASS data, we construct the UV-IR broad-band spectral energy distributions for various types of galaxies, and find that they form a nearly one-parameter family. For example, based on SDSS u- and r-band data, supplemented with redshift, the K-band 2MASS magnitudes can be "predicted" with an rms scatter of only 0.2 mag. When a dust content estimate determined from SDSS data by Kauffmann et al. (2003) is also utilized, this scatter decreases to 0.1 mag. We demonstrate that this dust content is indeed higher for galaxies detected by IRAS and that it can be used to "predict" measured IRAS 60 micron flux density within a factor of two using only SDSS data. We also show that the position of a galaxy in the emission-line-based Baldwin-Phillips-Terlevich diagram is correlated with the optical light concentration index and u-r color determined from the SDSS broad-band imaging data, and discuss changes in the morphology of this diagram induced by requiring detections at other wavelengths. We study the IR-radio correlation and find evidence that its slope may be different for AGN and star-forming galaxies and related to the H_alpha/H_beta line strength ratio.

We synthesize the rest-frame Stroemgren colors using SDSS spectra for 99,088 galaxies selected from Data Release 1. This narrow-band ~200 AA photometric system (uz, vz, bz, yz), first designed for the determination of effective temperature, metallicity and gravity of stars, measures the continuum spectral slope of galaxies in the rest-frame 3200-5800 AA wavelength range. Galaxies form a remarkably narrow locus (~0.03 mag) in the resulting color-color diagram. The Bruzual & Charlot population synthesis models suggest that the position of a galaxy along this locus is controlled by a degenerate combination of metallicity and age of the dominant stellar population. Galaxy distribution along the locus is bimodal, with the local minimum corresponding to an ~1 Gyr old single stellar population. The position perpendicular to the locus is independent of metallicity and age, and reflects the galaxy's dust content, as implied by both the models and the statistics of IRAS detections. A comparison of this locus with the galaxy locus in the H_delta-D_n(4000) diagram, utilized by Kauffmann et al. (2003) to estimate stellar masses, reveals a tight correlation, although the two analyzed spectral ranges barely overlap. Overall, the galaxy spectral energy distribution in the entire UV to near-IR range can be described as a single-parameter family with an accuracy of 0.1 mag, or better. This nearly one-dimensional distribution of galaxies in the multi-dimensional space of measured parameters strongly supports the conclusion of Yip et al. (2004), based on a principal component analysis, that SDSS galaxy spectra can be described by a small number of eigenspectra. Apparently, the contributions of stellar populations that dominate the optical emission from galaxies are combined in a simple and well-defined way.

We present the 3D real space clustering power spectrum of a sample of \~600,000 luminous red galaxies (LRGs) measured by the Sloan Digital Sky Survey (SDSS), using photometric redshifts. This sample of galaxies ranges from redshift z=0.2 to 0.6 over 3,528 deg^2 of the sky, probing a volume of 1.5 (Gpc/h)^3, making it the largest volume ever used for galaxy clustering measurements. We measure the angular clustering power spectrum in eight redshift slices and combine these into a high precision 3D real space power spectrum from k=0.005 (h/Mpc) to k=1 (h/Mpc). We detect power on gigaparsec scales, beyond the turnover in the matter power spectrum, on scales significantly larger than those accessible to current spectroscopic redshift surveys. We also find evidence for baryonic oscillations, both in the power spectrum, as well as in fits to the baryon density, at a 2.5 sigma confidence level. The statistical power of these data to constrain cosmology is ~1.7 times better than previous clustering analyses. Varying the matter density and baryon fraction, we find \Omega_M = 0.30 \pm 0.03, and \Omega_b/\Omega_M = 0.18 \pm 0.04, The detection of baryonic oscillations also allows us to measure the comoving distance to z=0.5; we find a best fit distance of 1.73 \pm 0.12 Gpc, corresponding to a 6.5% error on the distance. These results demonstrate the ability to make precise clustering measurements with photometric surveys (abridged).

(abridged) We have investigated 136 Chandra extragalactic sources without broad optical emission lines, including 93 galaxies with narrow emission lines (NELG) and 43 with only absorption lines (ALG). Based on fx/fo, Lx, X-ray spectral hardness and optical emission line diagnostics, we have conservatively classified 36 normal galaxies (20 spirals and 16 ellipticals) and 71 AGNs. We found no statistically significant evolution in Lx/LB, within the limited z range. We have built log(N)-log(S), after correcting for completeness based on a series of simulations. The best-fit slope is -1.5 for both S and B energy bands, which is considerably steeper than that of the AGN-dominated cosmic background sources, but slightly flatter than the previous estimate, indicating normal galaxies will not exceed the AGN population until fx ~ 2 x 10-18 erg s-1 cm-2 (a factor of ~5 lower than the previous estimate). A group of NELGs appear to be heavily obscured in X-rays, i.e., a typical type 2 AGN. After correcting for intrinsic absorption, their X-ray luminosities could be Lx > 10^44 erg s-1, making them type 2 quasar candidates. While most X-ray luminous ALGs (XBONG - X-ray bright, optically normal galaxy candidates) do not appear to be significantly absorbed, we found two heavily obscured objects, which could be as luminous as an unobscured broad-line quasar. Among 43 ALGs, we found two E+A galaxy candidates with strong Balmer absorption lines, but no [OII] line. The X-ray spectra of both galaxies are soft and one of them has a nearby close companion galaxy, supporting the merger/interaction scenario rather than the dusty starburst hypothesis.

We determine the number counts and z=0-5 luminosity function for a well-defined, homogeneous sample of quasars from the Sloan Digital Sky Survey (SDSS). We conservatively define the most uniform statistical sample possible, consisting of 15,343 quasars within an effective area of 1622 deg^2 that was derived from a parent sample of 46,420 spectroscopically confirmed broad-line quasars in the 5282 deg^2 of imaging data from SDSS Data Release Three. The sample extends from i=15 to i=19.1 at z<3 and to i=20.2 for z>3. The number counts and luminosity function agree well with the results of the 2dF QSO Survey, but the SDSS data probe to much higher redshifts than does the 2dF sample. The number density of luminous quasars peaks between redshifts 2 and 3, although uncertainties in the selection function in this range do not allow us to determine the peak redshift more precisely. Our best fit model has a flatter bright end slope at high redshift than at low redshift. For z<2.4 the data are best fit by a redshift-independent slope of beta = -3.1 (Phi(L) propto L^beta). Above z=2.4 the slope flattens with redshift to beta=-2.37 at z=5. This slope change, which is significant at a >5-sigma level, must be accounted for in models of the evolution of accretion onto supermassive black holes.

We examine the ``puffed-up inner disk'' model (Dullemond, Dominik & Natta 2001), proposed for explaining the near-IR excess radiation from Herbig Ae/Be stars. Detailed model computations show that the observed near-IR excess requires more hot dust than is contained in the puffed-up disk rim. The rim can produce the observed near-IR excess only if its dust has perfectly gray opacity, but such dust is in conflict with the observed 10$\mu$m spectral feature. We find that a compact ($\sim$ 10 AU) tenuous ($\tau_V \la 0.4$) dusty halo around the disk inner regions contains enough dust to readily explain the observations. Furthermore, this model also resolves the puzzling relationship noted by Monnier and Millan-Gabet (2002) between luminosity and the interferometric inner radii of disks.

We study the luminosity and color dependence of the galaxy 2-point correlation function in the Sloan Digital Sky Survey, starting from a sample of 200,000 galaxies over 2500 deg^2. We concentrate on the projected correlation function w(r_p), which is directly related to the real space \xi(r). The amplitude of w(r_p) grows continuously with luminosity, rising more steeply above the characteristic luminosity L_*. Redder galaxies exhibit a higher amplitude and steeper correlation function at all luminosities. The correlation amplitude of blue galaxies increases continuously with luminosity, but the luminosity dependence for red galaxies is less regular, with bright red galaxies more strongly clustered at large scales and faint red galaxies more strongly clustered at small scales. We interpret these results using halo occupation distribution (HOD) models assuming concordance cosmological parameters. For most samples, an HOD model with two adjustable parameters fits the w(r_p) data better than a power-law, explaining inflections at r_p ~ 1-3 Mpc/h as the transition between the 1-halo and 2-halo regimes of \xi(r). The implied minimum mass for a halo hosting a central galaxy above a luminosity threshold L grows as M_min ~ L at low luminosities and more steeply above L_*. The mass at which an average halo has one satellite galaxy brighter than L is M_1 ~ 23 M_min(L). These results imply a conditional luminosity function (at fixed halo mass) in which central galaxies lie far above a Schechter function extrapolation of the satellite population. HOD models nicely explain the joint luminosity-color dependence of w(r_p) in terms of the color fractions of central and satellite populations as a function of halo mass. The inferred HOD properties are in good qualitative agreement with theoretical predictions.

We have examined the radial velocity data for stars spectroscopically observed by the Sloan Digital Sky Survey (SDSS) more than once to investigate the incidence of spectroscopic binaries, and to evaluate the accuracy of the SDSS stellar radial velocities. We find agreement between the fraction of stars with significant velocity variations and the expected fraction of binary stars in the halo and thick disk populations. The observations produce a list of 675 possible new spectroscopic binary stars and orbits for eight of them.

We analyze the SDSS Ly-alpha forest P_F(k,z) measurement to determine the linear theory power spectrum. Our analysis is based on fully hydrodynamic simulations, extended using hydro-PM simulations. We account for the effect of absorbers with damping wings, which leads to an increase in the slope of the linear power spectrum. We break the degeneracy between the mean level of absorption and the linear power spectrum without significant use of external constraints. We infer linear theory power spectrum amplitude Delta^2_L(k_p=0.009s/km,z_p=3.0)=0.452_{-0.057-0.116}^{+0.069+0.141} and slope n_eff=-2.321_{-0.047-0.102}^{+0.055+0.131} (possible systematic errors are included through nuisance parameters in the fit - a factor >~5 smaller errors would be obtained on both parameters if we ignored modeling uncertainties). The errors are correlated and not perfectly Gaussian, so we provide a chi^2 table to accurately describe the results. The result corresponds to sigma_8=0.85, n=0.94, for a LCDM model with Omega_m=0.3, Omega_b=0.04, and h=0.7, but is most useful in a combined fit with the CMB. The inferred curvature of the linear power spectrum and the evolution of its amplitude and slope with redshift are consistent with expectations for LCDM models, with the evolution of the slope, in particular, being tightly constrained. We use this information to constrain systematic contamination, e.g., fluctuations in the UV background. This paper should serve as a starting point for more work to refine the analysis, including technical improvements such as increasing the size and number of the hydrodynamic simulations, and improvements in the treatment of the various forms of feedback from galaxies and quasars.

The properties of the host galaxies of a well-defined sample of 2215 radio-loud AGN with redshifts 0.03 < z < 0.3, defined from the SDSS, are investigated. These are predominantly low radio luminosity sources, with 1.4GHz luminosities of 10^23 to 10^25 W/Hz. The fraction of galaxies that host radio-loud AGN with L(1.4GHz) > 10^23 W/Hz is a strong function of stellar mass, rising from nearly zero below a stellar mass of 10^10 Msun to more than 30% at 5x10^11 Msun. The integral radio luminosity function is derived in six ranges of stellar and black hole mass. Its shape is very similar in all of these ranges and can be well fitted by a broken power-law. Its normalisation varies strongly with mass, as M_*^2.5 or M_BH^1.6; this scaling only begins to break down when the predicted radio-loud fraction exceeds 20-30%. There is no correlation between radio and emission line luminosities for the radio-loud AGN in the sample and the probability that a galaxy of given mass is radio-loud is independent of whether it is optically classified as an AGN. The host galaxies of the radio-loud AGN have properties similar to those of ordinary galaxies of the same mass. All of these findings support the conclusion that the optical AGN and low radio luminosity AGN phenomena are independent and are triggered by different physical mechanisms. Intriguingly, the dependence on black hole mass of the radio-loud AGN fraction mirrors that of the rate at which gas cools from the hot atmospheres of elliptical galaxies. It is speculated that gas cooling provides a natural explanation for the origin of the radio-loud AGN activity, and it is argued that AGN heating could plausibly balance the cooling of the gas over time. [Abridged]

A sample of 2712 radio-luminous galaxies is defined from the second data release of the Sloan Digital Sky Survey (SDSS) by cross-comparing the main spectroscopic galaxy sample with two radio surveys: the NRAO VLA Sky Survey (NVSS) and the Faint Images of the Radio Sky at Twenty centimetres (FIRST) survey. The comparison is carried out in a multi-stage process and makes optimal use of both radio surveys by exploiting the sensitivity of the NVSS to extended and multi-component radio sources in addition to the high angular resolution of the FIRST images. A radio source sample with 95% completeness and 98.9% reliability is achieved, far better than would be possible for this sample if only one of the surveys was used. The radio source sample is then divided into two classes: radio-loud AGN and galaxies in which the radio emission is dominated by star formation. The division is based on the location of a galaxy in the plane of 4000Ang break strength versus radio luminosity per unit stellar mass and provides a sample of 2215 radio-loud AGN and 497 star forming galaxies brighter than 5mJy at 1.4GHz. A full catalogue of positions and radio properties is provided for these sources. The local radio luminosity function is then derived both for radio-loud AGN and for star-forming galaxies and is found to be in agreement with previous studies. By using the radio to far-IR correlation, the radio luminosity function of star forming galaxies is also compared to the luminosity function derived in the far-infrared. It is found to agree well at high luminosities but less so at lower luminosities, confirming that the linearity of the radio to far-IR correlation breaks down below about 10^22 W/Hz at 1.4GHz.

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.

We have used the 2dF instrument on the AAT to obtain redshifts of a sample of z<3, 18.0<g<21.85 quasars selected from SDSS imaging. These data are part of a larger joint programme: the 2dF-SDSS LRG and QSO Survey (2SLAQ). We describe the quasar selection algorithm and present the resulting luminosity function of 5645 quasars in 105.7 deg^2. The bright end number counts and luminosity function agree well with determinations from the 2dF QSO Redshift Survey (2QZ) data to g\sim20.2. However, at the faint end the 2SLAQ number counts and luminosity function are steeper than the final 2QZ results from Croom et al. (2004), but are consistent with the preliminary 2QZ results from Boyle et al. (2000). Using the functional form adopted for the 2QZ analysis, we find a faint end slope of beta=-1.78+/-0.03 if we allow all of the parameters to vary and beta=-1.45+/-0.03 if we allow only the faint end slope and normalization to vary. Our maximum likelihood fit to the data yields 32% more quasars than the final 2QZ parameterization, but is not inconsistent with other g>21 deep surveys. The 2SLAQ data exhibit no well defined ``break'' but do clearly flatten with increasing magnitude. The shape of the quasar luminosity function derived from 2SLAQ is in good agreement with that derived from type I quasars found in hard X-ray surveys. [Abridged]

We present a sample of 386 BL Lacertae (BL Lac) candidates identified from 2860 deg^2 of the Sloan Digital Sky Survey (SDSS) spectroscopic database. The candidates are primarily selected to have quasi-featureless optical spectra and low proper motions as measured from SDSS and USNO-B positions; however, our ability to separate Galactic from extragalactic quasi-featureless objects (QFOs) on the basis of proper motion alone is limited by the lack of reliable proper motion measurements for faint objects. Fortunately, high proper motion QFOs, mostly DC white dwarfs, populate a well defined region of color space, approximately corresponding to blackbodies with temperatures in the range 7000-12000 K. QFOs with measurable redshifts or X-ray or radio counterparts (i.e., evidence of an extragalactic/AGN nature) loosely follow a track in color space that corresponds to power-law continua plus host galaxy starlight, with typical power-law slopes in the range 1<alpha_{opt}<2 (f_{nu} propto nu^{-alpha}). Based largely on this remarkably clean color separation, we subdivide the sample into 240 probable candidates and 146 additional less probable (likely stellar) candidates. The probable BL Lac candidates have multi-wavelength properties consistent with the range of previously known BL Lacs, with an apparent preponderance of objects with synchrotron peaks at relatively high energies (HBL/ XBL-type). The majority of the 154 objects with measurable redshifts have z<1, with a median of 0.45; there are also a handful of high-redshift objects extending up to z=5.03. We identify a small number of potential radio-quiet BL Lac candidates, although more sensitive radio observations are needed to confirm their radio-quiet nature.

Here we present the New York University Value-Added Galaxy Catalog (NYU-VAGC), a catalog of local galaxies (mostly below a redshift of about 0.3) based on a set of publicly-released surveys (including the 2dFGRS, 2MASS, PSCz, FIRST, and RC3) matched to the Sloan Digital Sky Survey (SDSS) Data Release 2. Excluding areas masked by bright stars, the photometric sample covers 3514 square degrees and the spectroscopic sample covers 2627 square degrees (with about 85% completeness). Earlier, proprietary versions of this catalog have formed the basis of many SDSS investigations of the power spectrum, correlation function, and luminosity function of galaxies. We calculate and compile derived quantities (for example, K-corrections and structural parameters for galaxies). The SDSS catalog presented here is photometrically recalibrated, reducing systematic calibration errors across the sky from about 2% to about 1%. We include an explicit description of the geometry of the catalog, including all imaging and targeting information as a function of sky position. Finally, we have performed eyeball quality checks on a large number of objects in the catalog in order to flag deblending and other errors. This catalog is complementary to the SDSS Archive Servers, in that NYU-VAGC's calibration, geometrical description, and conveniently small size are specifically designed for studying galaxy properties and large-scale structure statistics using the SDSS spectroscopic catalog.

We present the large-scale correlation function measured from a spectroscopic sample of 46,748 luminous red galaxies from the Sloan Digital Sky Survey. The survey region covers 0.72 h^{-3} Gpc^3 over 3816 square degrees and 0.16<z<0.47, making it the best sample yet for the study of large-scale structure. We find a well-detected peak in the correlation function at 100h^{-1} Mpc separation that is an excellent match to the predicted shape and location of the imprint of the recombination-epoch acoustic oscillations on the low-redshift clustering of matter. This detection demonstrates the linear growth of structure by gravitational instability between z=1000 and the present and confirms a firm prediction of the standard cosmological theory. The acoustic peak provides a standard ruler by which we can measure the ratio of the distances to z=0.35 and z=1089 to 4% fractional accuracy and the absolute distance to z=0.35 to 5% accuracy. From the overall shape of the correlation function, we measure the matter density Omega_mh^2 to 8% and find agreement with the value from cosmic microwave background (CMB) anisotropies. Independent of the constraints provided by the CMB acoustic scale, we find Omega_m = 0.273 +- 0.025 + 0.123 (1+w_0) + 0.137 Omega_K. Including the CMB acoustic scale, we find that the spatial curvature is Omega_K=-0.010+-0.009 if the dark energy is a cosmological constant. More generally, our results provide a measurement of cosmological distance, and hence an argument for dark energy, based on a geometric method with the same simple physics as the microwave background anisotropies. The standard cosmological model convincingly passes these new and robust tests of its fundamental properties.

We summarize the Sloan Digital Sky Survey data acquisition and processing steps, and describe runQA, a pipeline designed for automated data quality assessment. In particular, we show how the position of the stellar locus in color-color diagrams can be used to estimate the accuracy of photometric zeropoint calibration to better than 0.01 mag in 0.03 deg2 patches. Using this method, we estimate that typical photometric zeropoint calibration errors for SDSS imaging data are not larger than ~0.01 mag in the g, r, and i bands, 0.02 mag in the z band, and 0.03 mag in the u band (root-mean-scatter for zeropoint offsets).

We present a catalog of 100,563 unresolved, UV-excess (UVX) quasar candidates to g=21 from 2099 deg^2 of the Sloan Digital Sky Survey (SDSS) Data Release One (DR1) imaging data. Existing spectra of 22,737 sources reveals that 22,191 (97.6%) are quasars; accounting for the magnitude dependence of this efficiency, we estimate that 95,502 (95.0%) of the objects in the catalog are quasars. Such a high efficiency is unprecedented in broad-band surveys of quasars. This ``proof-of-concept'' sample is designed to be maximally efficient, but still has 94.7% completeness to unresolved, g<~19.5, UVX quasars from the DR1 quasar catalog. This efficient and complete selection is the result of our application of a probability density type analysis to training sets that describe the 4-D color distribution of stars and spectroscopically confirmed quasars in the SDSS. Specifically, we use a non-parametric Bayesian classification, based on kernel density estimation, to parameterize the color distribution of astronomical sources -- allowing for fast and robust classification. We further supplement the catalog by providing photometric redshifts and matches to FIRST/VLA, ROSAT, and USNO-B sources. Future work needed to extend the this selection algorithm to larger redshifts, fainter magnitudes, and resolved sources is discussed. Finally, we examine some science applications of the catalog, particularly a tentative quasar number counts distribution covering the largest range in magnitude (14.2<g<21.0) ever made within the framework of a single quasar survey.

We present Sloan Digital Sky Survey and United States Naval Observatory observations of the V1647 Ori protostar and surrounding field near NGC 2068. V1647 Ori, the likely driving source for HH 23, brightened significantly in November 2003. Analysis of SDSS imaging acquired in November 1998 and February 2002 during the quiescent state, recent USNO photometry, and published 2MASS and Gemini data shows that the color changes associated with brightening suggest an EXor outburst rather than a simple dust clearing event.

We report the discovery of a new dwarf spheroidal satellite of M31, Andromeda IX, based on resolved stellar photometry from the Sloan Digital Sky Survey (SDSS). Using both SDSS and public archival data we have estimated its distance and other physical properties, and compared these to the properties of a previously known dwarf spheroidal companion, Andromeda V, also observed by SDSS. Andromeda IX is the lowest surface brightness galaxy found to date (mu_{V} ~ 26.8 mag arcsec^-2), and at the distance we estimate from the position of the tip of Andromeda IX's red giant branch, (m - M)_{0} ~ 24.5 (805 kpc), Andromeda IX would also be the faintest galaxy known (M_{V} ~ -8.3).

The Sloan Digital Sky Survey has revealed an overdensity of luminous red giant stars ~ 3 degrees (40 projected kpc) to the northeast of M31, which we have called Andromeda NE. The line-of-sight distance to Andromeda NE is within approximately 50 kpc of M31; Andromeda NE is not a physically unrelated projection. Andromeda NE has a g-band absolute magnitude of ~ -11.6 and central surface brightness of ~ 29 mag/sq.arcsec, making it nearly two orders of magnitude more diffuse than any known Local Group dwarf galaxy at that luminosity. Based on its distance and morphology, Andromeda NE is likely undergoing tidal disruption. Andromeda NE's red giant branch color is unlike that of M31's present-day outer disk or the stellar stream reported by Ibata et al. (2001), arguing against a direct link between Andromeda NE and these structures. However, Andromeda NE has a red giant branch color similar to that of the G1 clump; it is possible that these structures are both material torn off of M31's disk in the distant past, or that these are both part of one ancient stellar stream.

We combine the constraints from the recent Ly-alpha forest and bias analysis of the SDSS with previous constraints from SDSS galaxy clustering, the latest supernovae, and WMAP . Combining WMAP and the lya we find for the primordial slope n_s=0.98\pm 0.02. We see no evidence of running, dn/d\ln k=-0.003\pm 0.010, a factor of 3 improvement over previous constraints. We also find no evidence of tensors, r<0.36 (95% c.l.). A positive correlation between tensors and primordial slope disfavors chaotic inflation type models with steep slopes: V \propto \phi^4 is outside the 3-sigma contour. For the amplitude we find sigma_8=0.90\pm 0.03 from the lyaf and WMAP alone. We find no evidence of neutrino mass: for the case of 3 massive neutrino families with an inflationary prior, \sum m_{\nu}<0.42eV and the mass of lightest neutrino is m_1<0.13eV at 95% c.l. For the 3 massless + 1 massive neutrino case we find m_{\nu}<0.79eV for the massive neutrino, excluding at 95% c.l. all neutrino mass solutions compatible with the LSND results. We explore dark energy constraints in models with a fairly general time dependence of dark energy equation of state, finding Omega_lambda=0.72\pm 0.02, w(z=0.3)=-0.98^{+0.10}_{-0.12}, the latter changing to w(z=0.3)=-0.92^{+0.09}_{-0.10} if tensors are allowed. We find no evidence for variation of the equation of state with redshift, w(z=1)=-1.03^{+0.21}_{-0.28}. These results rely on the current understanding of the lyaf and other probes, which need to be explored further both observationally and theoretically, but extensive tests reveal no evidence of inconsistency among different data sets used here.

We perform an objective classification of 170,000 galaxy spectra in the Sloan Digital Sky Survey (SDSS) using the Karhunen-Lo\`eve (KL) transform. With about one-sixth of the total set of galaxy spectra which will be obtained by the survey, we are able to carry out the most extensive analysis of its kind to date. The formalism proposed by Connolly and Szalay is adopted to correct for gappy regions in the spectra, and to derive eigenspectra and eigencoefficients. From this analysis, we show that this gap-correction formalism leads to a converging set of eigenspectra and KL-repaired spectra. Furthermore, KL eigenspectra of galaxies are found to be convergent not only as a function of iteration, but also as a function of the number of randomly selected galaxy spectra used in the analysis. From these data a set of ten eigenspectra of galaxy spectra are constructed, with rest-wavelength coverage 3450-8350 angstrom. The eigencoefficients describing these galaxies naturally place the spectra into several classes defined by the plane formed by the first three eigencoefficients of each spectrum. Spectral types, corresponding to different Hubble-types and galaxies with extreme emission lines, are identified for the 170,000 spectra and are shown to be complementary to existing spectral classifications. From a non-parametric classification technique, we find that the population of galaxies can be divided into three classes which correspond to early late- through to intermediate late-types galaxies. (Abridged)

Type II quasars are luminous AGNs whose central engines and broad-line regions are obscured by intervening material; such objects only recently have been discovered in appreciable numbers. We study the multiwavelength properties of 291 type II AGN candidates (0.3 < z < 0.8) selected based on their optical emission line properties from the spectroscopic database of the Sloan Digital Sky Survey. This sample includes about 150 objects luminous enough to be classified as type II quasars. We matched the sample to the FIRST (20 cm), IRAS (12-100 micron), 2MASS (JHK_S) and RASS (0.1-2.4 keV) surveys. Roughly 10% of optically selected type II AGN candidates are radio-loud, comparable to the AGN population as a whole. About 40 objects are detected by IRAS at 60 micron and/or 100 micron, and the inferred mid/far-IR luminosities lie in the range nu L_nu=10^45-3x10^46 erg/sec. Average IR-to-[OIII]5007 ratios of objects in our sample are consistent with those of other AGNs. Objects from our sample are ten times less likely to have soft X-ray counterparts in RASS than type I AGNs with the same redshifts and [OIII]5007 luminosities. The few type II AGN candidates from our sample that are detected by RASS have harder X-ray spectra than those of type I AGNs. The multiwavelength properties of the type II AGN candidates from our sample are consistent with their interpretation as powerful obscured AGNs.

We summarize the detection rates at wavelengths other than optical for \~99,000 galaxies from the Sloan Digital Sky Survey (SDSS) Data Release 1 ``main'' spectroscopic sample. The analysis is based on positional cross-correlation with source catalogs from ROSAT, 2MASS, IRAS, GB6, FIRST, NVSS and WENSS surveys. We find that the rest-frame UV-IR broad-band galaxy SEDs form a remarkably uniform, nearly one parameter, family. As an example, the SDSS u and r band data, supplemented with redshift, can be used to predict K band magnitudes measured by 2MASS with an rms scatter of only 0.2 mag; when measurement uncertainties are taken into account, the astrophysical scatter appears not larger than ~0.1 mag.

We analyze the properties of quasar variability using repeated SDSS imaging data in five UV-to-far red photometric bands, accurate to 0.02 mag, for 13,000 spectroscopically confirmed quasars. The observed time lags span the range from 3 hours to over 3 years, and constrain the quasar variability for rest-frame time lags of up to two years, and at rest-frame wavelengths from 1000 Ang. to 6000 Ang. We demonstrate that 66,000 SDSS measurements of magnitude differences can be described within the measurement noise by a simple function of only three free parameters. The addition of POSS data constrains the long-term behavior of quasar variability and provides evidence for a turn-over in the structure function. This turn-over indicates that the characteristic time scale for optical variability of quasars is of the order 1 year.

We present a spectroscopic analysis of nearly 8000 late-type dwarfs in the Sloan Digital Sky Survey. Using the Halpha emission line as an activity indicator, we investigate the fraction of active stars as a function of spectral type and find a peak near type M8, confirming previous results. In contrast to past findings, we find that not all M7-M8 stars are active. We show that this may be a selection effect of the distance distributions of previous samples, as the active stars appear to be concentrated near the Galactic Plane. We also examine the activity strength (ratio of the luminosity emitted in Halpha to the bolometric luminosity) for each star, and find that the mean activity strength is constant over the range M0-M5 and declines at later types. The decline begins at a slightly earlier spectral type than previously found. We explore the effect that activity has on the broadband photometric colors and find no significant differences between active and inactive stars. We also carry out a search for subdwarfs using spectroscopic metallicity indicators, and find 60 subdwarf candidates. Several of these candidates are near the extreme subdwarf boundary. The spectroscopic subdwarf candidates are redder by \~0.2 magnitudes in g-r compared to disk dwarfs at the same r-i color.

We present a study of variable faint optical sources discovered by comparing the Sloan Digital Sky Survey (SDSS) and the Palomar Observatory Sky Survey (POSS) catalogs. We use SDSS measurements to photometrically recalibrate several publicly available POSS catalogs; a piecewise recalibration in 100 arcmin2 patches generally results in an improvement of photometric accuracy (rms) by nearly a factor of two, compared to the original data. The POSS I magnitudes can be improved to ~0.15 mag accuracy, and POSS II magnitudes to \~0.10 mag accuracy. We use the recalibrated catalogs for the ~2,000 deg2 of sky in the SDSS Data Release 1 to construct a catalog of ~60,000 sources variable on time scales 10-50 years. A series of statistical tests based on the morphology of SDSS color-color diagrams, as well as visual comparison of images and comparison with repeated SDSS observations, demonstrate the robustness of the selection methods. We quantify the distribution of variable sources in the SDSS color-color diagrams, and the variability characteristics of quasars. We detect a turn-over in quasar structure function which suggests that the characteristic time scale for quasar variability is of the order one year. The long-term (>1 year) quasar variability decreases with luminosity and rest-frame wavelength similarly to the short-term (<1 year) behavior. We also demonstrate that candidate RR Lyrae stars trace the same halo structures, such as the Sgr dwarf tidal stream, that were discovered using repeated SDSS observations. We utilize the POSS-SDSS selected candidates to constrain the halo structure in the parts of sky for which repeated SDSS observations do not exist. (abridged)

We positionally match the 6 cm GB6, 20 cm FIRST and NVSS, and 92 cm WENSS radio catalogs and find 16,500 matches in ~3,000 deg2 of sky. Using this unified radio database, we construct radio "color-magnitude-morphology" diagrams and find that they display a clear structure, rather than a random scatter. We propose a simple, yet powerful, method for morphological classification of radio sources based on FIRST and NVSS measurements. For a subset of matched sources, we find optical identifications using the SDSS Data Release 1 catalogs, and separate them into quasars and galaxies. Compact radio sources with flat radio spectra are dominated by quasars, while compact sources with steep spectra, and resolved radio sources, contain substantial numbers of both quasars and galaxies.

We report the discovery of a locus of stars in the SDSS g-r vs. u-g color-color diagram that connects the colors of white dwarfs and M dwarfs. While its contrast with respect to the main stellar locus is only ~1:2300, this previously unrecognized feature includes 863 stars from the SDSS Data Release 1. The position and shape of the feature are in good agreement with predictions of a simple binary star model that consists of a white dwarf and an M dwarf, with the components' luminosity ratio controlling the position along this binary system locus. SDSS DR1 spectra for 47 of these objects strongly support this model. The absolute magnitude--color distribution inferred for the white dwarf component is in good agreement with the models of Bergeron et al. (1995).

We report the first successful application of the astrometric color-induced displacement technique (CID, the displacement of the photocenter between different bandpasses due to a varying contribution of differently colored components to the total light), originally proposed by Wielen (1996) for discovering unresolved binary stars. Using the Sloan Digital Sky Survey (SDSS) Data Release 1 with 2.5 million stars brighter than 21m in the u and g bands, we select 419 candidate binary stars with CID greater than 0.5 arcsec. The SDSS colors of the majority of these candidates are consistent with binary systems including a white dwarf and any main sequence star with spectral type later than ~K7. The astrometric CID method discussed here is complementary to the photometric selection of binary stars in SDSS discussed by Smolcic et al. (2004), but there is considerable overlap (15%) between the two samples of selected candidates. This overlap testifies both to the physical soundness of both methods, as well as to the astrometric and photometric quality of SDSS data.

We present a detailed observational strategy for finding Near Earth Objects (NEOs) with the Sloan Digital Sky Survey (SDSS) telescope. We investigate strategies in normal, unbinned mode as well as binning the CCDs 2x2 or 3x3, which affects the sky coverage rate and the limiting apparent magnitude. We present results from 1 month, 3 year and 10 year simulations of such surveys. For each cadence and binning mode, we evaluate the possibility of achieving the Spaceguard goal of detecting 90% of 1 km NEOs (absolute magnitude H <= 18 for an albedo of 0.1). We find that an unbinned survey is most effective at detecting H <= 20 NEOs in our sample. However, a 3x3 binned survey reaches the Spaceguard Goal after only seven years of operation. As the proposed large survey telescopes (PanStarss; LSST) are at least 5-10 years from operation, an SDSS NEO survey could make a significant contribution to the detection and photometric characterization of the NEO population.

We present evidence for a ring of stars in the plane of the Milky Way, extending at least from l = 180 deg to l = 227 deg with turnoff magnitude $g \sim 19.5$; the ring could encircle the Galaxy. We infer that the low Galactic latitude structure is at a fairly constant distance of R = 18 +/- 2 kpc from the Galactic Center above the Galactic plane, and has R = 20 +/- 2 kpc in the region sampled below the Galactic plane. The evidence includes five hundred SDSS spectroscopic radial velocities of stars within 30 degrees of the plane. The velocity dispersion of the stars associated with this structure is found to be 27 km/s at (l,b) = (198, -27), 22 km/s at (l,b) = (225, 28), 30 km/s at (l,b) = (188, 24), and 30 km/s at (l,b) = (182, 27) degrees. The structure rotates in the same prograde direction as the Galactic disk stars, but with a circular velocity of 110+/-25 km/s. The narrow measured velocity dispersion is inconsistent with power law spheroid or thick disk populations. We compare the velocity dispersion in this structure with the velocity dispersion of stars in the Sagittarius dwarf galaxy tidal stream, for which we measure a velocity dispersion of 20 km/s at (l, b) = (165, -55) deg. We estimate a preliminary metallicity from the Ca II (K) line and color of the turnoff stars of [Fe/H] = -1.6 with a dispersion of 0.3 dex and note that the turnoff color is consistent with that of the spheroid population. We interpret our measurements as evidence for a tidally disrupted satellite of $2 \times 10^7$ to $5 \times 10^8 \rm M_\odot$ which rings the Galaxy.

We study the spatial variation of galaxy number counts using five band photometric images from the Sloan Digital Sky Survey. The spatial variation of this sample of 46 million galaxies collected from 2200 sq. degrees can be understood as the combination of Galactic extinction and large-scale clustering. With the use of the reddening map of Schlegel, Finkbeiner & Davis (1998), the standard extinction law is verified for the colour bands from u to z within 5% in the region of small extinction values, E(B-V)<0.15. The residual spatial variations of the number counts suggests that the error of global calibration for SDSS photometry is smaller than 0.02 mag.

Gravitational lensing is a powerful tool for the study of the distribution of dark matter in the Universe. The cold-dark-matter model of the formation of large-scale structures predicts the existence of quasars gravitationally lensed by concentrations of dark matter so massive that the quasar images would be split by over 7 arcsec. Numerous searches for large-separation lensed quasars have, however, been unsuccessful. All of the roughly 70 lensed quasars known, including the first lensed quasar discovered, have smaller separations that can be explained in terms of galaxy-scale concentrations of baryonic matter. Although gravitationally lensed galaxies with large separations are known, quasars are more useful cosmological probes because of the simplicity of the resulting lens systems. Here we report the discovery of a lensed quasar, SDSS J1004+4112, which has a maximum separation between the components of 14.62 arcsec. Such a large separation means that the lensing object must be dominated by dark matter. Our results are fully consistent with theoretical expectations based on the cold-dark-matter model.

We carry out a maximum-likelihood kinematic analysis of a sample of 1170 blue horizontal branch (BHB) stars from the Sloan Digital Sky Survey presented in Sirko et al. (2003) (Paper I). Monte Carlo simulations and resampling show that the results are robust to distance and velocity errors at least as large as the estimated errors from Paper I. The best-fit velocities of the Sun (circular) and halo (rotational) are 245.9 +/- 13.5 km/s and 23.8 +/- 20.1 km/s but are strongly covariant, so that v_0 - v_halo = 222.1 +/- 7.7 km/s. If one adopts standard values for the local standard of rest and solar motion, then the halo scarcely rotates. The velocity ellipsoid inferred for our sample is much more isotropic [(sigma_r,sigma_theta,sigma_phi) = (101.4 +/- 2.8, 97.7 +/- 16.4, 107.4 +/- 16.6) km/s] than that of halo stars in the solar neighborhood, in agreement with a recent study of the distant halo by Sommer-Larsen et al. (1997). The line-of-sight velocity distribution of the entire sample, corrected for the Sun's motion, is accurately gaussian with a dispersion of 101.6 +/- 3.0 km/s.

We isolate samples of 733 bright (g < 18) and 437 faint (g > 18) high-Galactic latitude blue horizontal branch stars with photometry and spectroscopy in the Sloan Digital Sky Survey (SDSS). Comparison of independent photometric and spectroscopic selection criteria indicates that contamination from F and blue-straggler stars is less than 10% for bright stars (g<18) and about 25% for faint stars (g>18), and this is qualitatively confirmed by proper motions based on the USNO-A catalog as first epoch. Analysis of repeated observations shows that the errors in radial velocity are approximately 26 km/s. A relation between absolute magnitude and color is established using the horizontal branches of halo globular clusters observed by SDSS. Bolometric corrections and colors are synthesized in the SDSS filters from model spectra. The redder stars agree well in absolute magitude with accepted values for RR Lyrae stars. The resulting photometric distances are accurate to about 0.2 magnitudes, with a median of about 25 kpc. Modest clumps in phase space exist and are consistent with the previously reported tidal stream of the Sagittarius dwarf galaxy. The sample is tabulated in electronic form in the online version of this article, or by request to the authors.

We utilize a complete sample of RR Lyrae stars discovered by the QUEST survey using light curves to design selection criteria based on SDSS colors. Thanks to the sensitivity of the u-g color to surface gravity and of g-r color to effective temperature, and to the small photometric errors (~0.02 mag) delivered by SDSS, RR Lyrae stars can be efficiently and robustly recognized even with single-epoch data. In a 100% complete color-selected sample, the selection efficiency (the fraction of RR Lyrae stars in the candidate sample) is 6%, and, by adjusting color cuts, it can be increased to 10% with a completeness of 80%, and to 60% with 28% completeness. Such color selection produces samples that are sufficiently clean for statistical studies of the Milky Way's halo substructure, and we utilize it to select 3,643 candidate RR Lyrae stars from SDSS Data Release 1. We demonstrate that this sample recovers known clumps of RR Lyrae stars associated with the Sgr dwarf tidal tail, and Pal 5 globular cluster, and use it to constrain the halo substructure away from the Sgr dwarf tidal tail. These results suggest that it will be possible to study the halo substructure out to ~70 kpc from the Galactic Center in the entire area imaged by the SDSS, and not only in the multiply observed regions.

We analyze the counts of low-redshift quasar candidates selected using nine-epoch SDSS imaging data. The co-added catalogs are more than 1 mag deeper than single-epoch SDSS data, and allow the selection of low-redshift quasar candidates using UV-excess and also variability techniques. The counts of selected candidates are robustly determined down to g=21.5. This is about 2 magnitudes deeper than the position of a change in the slope of the counts reported by Boyle et al. (1990, 2000) for a sample selected by UV-excess, and questioned by Hawkins & Veron (1995), who utilized a variability-selected sample. Using SDSS data, we confirm a change in the slope of the counts for both UV-excess and variability selected samples, providing strong support for the Boyle et al. results.

Recent claims by Ivezic et al. (2002) that the distribution of the radio-to-optical flux ratio, R, for quasars is bimodal (the so-called quasar radio dichotomy) were questioned on statistical grounds by Cirasuolo et al. (2003). We apply the approach suggested by Cirasuolo et al. to a sample of \~10,000 objects detected by SDSS and FIRST, and find support for the quasar radio dichotomy. The discrepancy between the claims by Cirasuolo et al. and the results presented here is most likely because 1) the ~100 times larger sample based on two homogeneous surveys that is used here allows a direct determination of the R distribution, rather than relying on indirect inferences based on Monte Carlo simulations of several heterogeneous surveys 2) the accurate SDSS colors and redshift information allow robust determination of the K-correction for R, which, if unaccounted for, introduces significant scatter that masks the intrinsic properties of the quasar R distribution.

Using a sample of over 25000 spectroscopically confirmed quasars from the Sloan Digital Sky Survey, we show how quasar variability in the rest frame optical/UV regime depends upon rest frame time lag, luminosity, rest wavelength, redshift, the presence of radio and X-ray emission, and the presence of broad absorption line systems. The time dependence of variability (the structure function) is well-fit by a single power law on timescales from days to years. There is an anti-correlation of variability amplitude with rest wavelength, and quasars are systematically bluer when brighter at all redshifts. There is a strong anti-correlation of variability with quasar luminosity. There is also a significant positive correlation of variability amplitude with redshift, indicating evolution of the quasar population or the variability mechanism. We parameterize all of these relationships. Quasars with RASS X-ray detections are significantly more variable (at optical/UV wavelengths) than those without, and radio loud quasars are marginally more variable than their radio weak counterparts. We find no significant difference in the variability of quasars with and without broad absorption line troughs. Models involving multiple discrete events or gravitational microlensing are unlikely by themselves to account for the data. So-called accretion disk instability models are promising, but more quantitative predictions are needed.

We analyse the bivariate distribution, in color versus absolute magnitude (u-r vs. M_r), of a low redshift sample of galaxies from the Sloan Digital Sky Survey (SDSS; 2400 deg^2, 0.004<z<0.08, -23.5<M_r<-15.5). We trace the bimodality of the distribution from luminous to faint galaxies by fitting double-Gaussians to the color functions separated in absolute magnitude bins. Color-magnitude (CM) relations are obtained for red and blue distributions (early- and late-type, predominantly field, galaxies) without using any cut in morphology. Instead, the analysis is based on the assumption of normal Gaussian distributions in color. We find that the CM relations are well fit by a straight line plus a tanh function. Both relations can be described by a shallow CM trend (slopes of about -0.04, -0.05) plus a steeper transition in the average galaxy properties over about two magnitudes. The midpoints of the transitions (M_r=-19.8 and -20.8 for the red and blue distributions, respectively) occur around 2x10^10 M_solar after converting luminosities to stellar mass. Separate luminosity functions are obtained for the two distributions. The red distribution has a more luminous characteristic magnitude and a shallower faint-end slope (M^*=-21.5, alpha=-0.8) compared to the blue distribution (alpha=-1.3 depending on the parameterization). These are approximately converted to galaxy stellar mass functions. The red distribution galaxies have a higher number density per magnitude for masses greater than about 3x10^10 M_solar. Using a simple merger model, we show that the differences between the two functions are consistent with the red distribution being formed from major galaxy mergers.

Type II quasars are the long-sought luminous analogs of type II (narrow emission line) Seyfert galaxies, suggested by unification models of active galactic nuclei (AGN) and postulated to account for an appreciable fraction of the cosmic hard X-ray background. We present a sample of 291 type II AGN at redshifts 0.3<Z<0.83 from the spectroscopic data of the Sloan Digital Sky Survey. These objects have narrow (FWHM<2000 km/s), high equivalent width emission lines with high-ionization line ratios. We describe the selection procedure and discuss the optical properties of the sample. About 50% of the objects have [OIII] \lambda_air 5007 line luminosities in the range 3\times 10^8-10^10 L_Sun, comparable to those of luminous (-27<M_B<-23) quasars; this, along with other evidence, suggests that the objects in the luminous subsample are type II quasars.

The Sloan Digital Sky Survey (SDSS) first data release provides a database of 106000 unique galaxies in the main galaxy sample with measured spectra. A sample of star-forming (SF) galaxies are identified from among the 3079 of these having 1.4 GHz luminosities from FIRST, by using optical spectral diagnostics. Using 1.4 GHz luminosities as a reference star formation rate (SFR) estimator insensitive to obscuration effects, the SFRs derived from the measured SDSS Halpha, [OII] and u-band luminosities, as well as far-infrared luminosities from IRAS, are compared. It is established that straightforward corrections for obscuration and aperture effects reliably bring the SDSS emission line and photometric SFR estimates into agreement with those at 1.4 GHz, although considerable scatter (~60%) remains in the relations. It thus appears feasible to perform detailed investigations of star formation for large and varied samples of SF galaxies through the available spectroscopic and photometric measurements from the SDSS. We provide herein exact prescriptions for determining the SFR for SDSS galaxies. The expected strong correlation between [OII] and Halpha line fluxes for SF galaxies is seen, but with a median line flux ratio F_[OII]/F_Halpha=0.23, about a factor of two smaller than that found in the sample of Kennicutt (1992). This correlation, used in deriving the [OII] SFRs, is consistent with the luminosity-dependent relation found by Jansen et al. (2001). The median obscuration for the SDSS SF systems is found to be A_Halpha=1.2 mag, while for the radio detected sample the median obscuration is notably higher, 1.6 mag, and with a broader distribution.

A new overdensity of A-colored stars in distant parts of the Milky Way's stellar halo, at a dereddened SDSS magnitude of g_0 = 20.3, is presented. Identification of associated variable RR Lyrae candidates supports the claim that these are blue horizontal branch stars. The inferred distance of these stars from the Galactic center is 90 kpc, assuming the absolute magnitude of these stars is M_g_0 = 0.7 and that the Sun is 8.5 kpc from the Galactic center. The new tidal debris is within 10 kpc of same plane as other confirmed tidal debris from the disruption of the Sagittarius dwarf galaxy, and could be associated with the trailing tidal arm. Distances to the Sagittarius stream estimated from M stars are about 13% smaller than our inferred distances. The tidal debris has a width of at least 10 degrees, and is traced for more than 20 degrees across the sky. The globular cluster NGC 2419 is located within the detected tidal debris, and may also have once been associated with the Sagittarius dwarf galaxy.

We report a detection of statistically significant color variations for a sample of 7,531 multiply observed asteroids that are listed in the Sloan Digital Sky Survey Moving Object Catalog. Using 5-band photometric observations accurate to ~0.02 mag, we detect color variations in the range 0.06-0.11 mag (rms). These variations appear uncorrelated with asteroids physical characteristics such as diameter (in the probed 1-10 km range), taxonomic class, and family membership. Despite such a lack of correlation, which implies a random nature for the cause of color variability, a suite of tests suggest that the detected variations are not instrumental effects. In particular, the observed color variations are incompatible with photometric errors, and, for objects observed at least four times, the color change in the first pair of observations is correlated with the color change in the second pair. These facts strongly suggest that the observed effect is real, and also indicate that for some asteroids color variations are larger than for others. The detected color variations can be explained as due to inhomogeneous albedo distribution over an asteroid surface. Although relatively small, these variations suggest that fairly large patches with different color than their surroundings exist on a significant fraction of asteroids. This conclusion is in agreement with spatially resolved color images of several large asteroids obtained by NEAR spacecraft and HST.

We discuss the density and radial velocity distributions of over 3000 candidate RR Lyrae stars selected by various methods using Sloan Digital Sky Survey data for about 1000 deg^2 of sky. This is more than 20 times larger sample than previously reported by SDSS (Ivezic et al. 2000), and includes candidate RR Lyrae stars out to the sample limit of 100 kpc. A cutoff in the radial distribution of halo RR Lyrae at ~50-60 kpc that was suggested by the early SDSS data appears to be a statistical anomaly confined to a small region (~100 deg^2). Despite the large increase in observed area, the most prominent features remain to be those associated with the Sgr dwarf tidal stream. We find multiple number density peaks along three lines of sight in the Sgr dwarf tidal stream plane, that may indicate several perigalactic passages of the Sgr dwarf galaxy.

We discuss a sample of over 3000 candidate RR Lyrae stars selected by various methods using Sloan Digital Sky Survey data for about 1000 deg^2 of sky. These stars probe the halo structure out to ~100 kpc from the Galactic center. Their spatial and radial velocity distributions are very inhomogeneous, with the most prominent features tracing the Sgr dwarf tidal stream. Outside the Sgr dwarf tidal stream, the spatial distribution in the 5-60 kpc range of Galactocentric radius R is well described by an R^{-3} power law.

Unified schemes of active galactic nuclei (AGN) require an obscuring dusty torus around the central source, giving rise to type 1 line spectrum for pole-on viewing and type 2 characteristics in edge-on sources. Infrared radiation at its different wavelengths is the best probe of the dust distribution, whether the torus orientation is edge-on or pole-on. The observed IR is in broad agreement with the unified scheme but serious problems remained in all early modeling efforts. In spite of a general awareness that the dust must be concentrated in clouds, clumpiness remained the one major ingredient missing in those radiative transfer studies because of the inherent difficulties it presents. We have recently developed the formalism to handle dust clumpiness and our results indicate that its inclusion may resolve the difficulties encountered by the previous theoretical efforts. We show that these problems find a natural explanation if the dust is contained in \about\ 5--10 clouds along radial rays through the torus. The spectral energy distributions (SED) of both type 1 and type 2 sources are properly reproduced from different viewpoints of the same object if the optical depth of each cloud is \ga 40 at visual wavelengths and the clouds' mean free path increases roughly in proportion to radial distance.

We study the IR emission from flared disks with and without additional optically thin halos. Flux calculations of a flared disk in vacuum can be considered a special case of the more general family of models in which the disk is imbedded in an optically thin halo. In the absence of such halo, flux measurements can never rule out its existence because the disk flaring surface defines a mathematically equivalent halo that produces the exact same flux at all IR wavelengths. When a flared disk with height H at its outer radius R is imbedded in a halo whose optical depth at visual wavelengths is \tau\sub{halo}, the system IR flux is dominated by the halo whenever \tau\sub{halo} > (1/4)H/R. Even when its optical depth is much smaller, the halo can still have a significant effect on the disk temperature profile. Imaging is the only way to rule out the existence of a potential halo, and we identify a decisive test that extracts a signature unique to flared disks from imaging observations.

The distribution of early-type galaxy velocity dispersions, phi(sigma), is measured using a sample drawn from the SDSS database. Its shape differs significantly from that which one obtains by simply using the mean correlation between luminosity, L, and velocity dispersion, sigma, to transform the luminosity function into a velocity function: ignoring the scatter around the mean sigma-L relation is a bad approximation. An estimate of the contribution from late-type galaxies is also made, which suggests that phi(sigma) is dominated by early-type galaxies at velocities larger than ~ 200 km/s.

We study the variation of the dark matter mass fraction of elliptical galaxies as a function of their luminosity, stellar mass, and size using a sample of 29,469 elliptical galaxies culled from the Sloan Digital Sky Survey. We model ellipticals as a stellar Hernquist profile embedded in an adiabatically compressed dark matter halo. This model allows us to estimate a dynamical mass ($M_{dynm}$) at the half-light radius from the velocity dispersion of the spectra, and to compare these to the stellar mass estimates ($M_{*}$) from Kauffmann et al (2003). We find that $M_{*}/L$ is independent of luminosity, while $M_{dynm}/L$ increases with luminosity, implying that the dark matter fraction increases with luminosity. We also observe that at a fixed luminosity or stellar mass, the dark matter fraction increases with increasing galaxy size or, equivalently, increases with decreasing surface brightness: high surface brightness galaxies show almost no evidence for dark matter, while in low surface brightness galaxies, the dark matter exceeds the stellar mass at the half light radius. We relate this to the fundamental plane of elliptical galaxies, suggesting that the tilt of this plane from simple virial predictions is due to the dark matter in galaxies. We find that a simple model where galaxies are embedded in dark matter halos and have a star formation efficiency independent of their surface brightness explains these trends. We estimate the virial mass of ellipticals as being approximately 7-30 times their stellar mass, with the lower limit suggesting almost all of the gas within the virial radius is converted into stars.

An initial assessment is made of white dwarf and hot subdwarf stars observed in the Sloan Digital Sky Survey. In a small area of sky (190 square degrees), observed much like the full survey will be, 269 white dwarfs and 56 hot subdwarfs are identified spectroscopically where only 44 white dwarfs and 5 hot subdwarfs were known previously. Most are ordinary DA (hydrogen atmosphere) and DB (helium) types. In addition, in the full survey to date, a number of WDs have been found with uncommon spectral types. Among these are blue DQ stars displaying lines of atomic carbon; red DQ stars showing molecular bands of C_2 with a wide variety of strengths; DZ stars where Ca and occasionally Mg, Na, and/or Fe lines are detected; and magnetic WDs with a wide range of magnetic field strengths in DA, DB, DQ, and (probably) DZ spectral types. Photometry alone allows identification of stars hotter than 12000 K, and the density of these stars for 15<g<20 is found to be ~2.2 deg^{-2} at Galactic latitudes 29-62 deg. Spectra are obtained for roughly half of these hot stars. The spectra show that, for 15<g<17, 40% of hot stars are WDs and the fraction of WDs rises to ~90% at g=20. The remainder are hot sdB and sdO stars.

We investigate the continuum and emission line properties of 4576 SDSS quasars as a function of their optical/UV SEDs. The optical/UV color distribution of our sample is roughly Gaussian, but with a red tail; we distinguish between 1) intrinsically blue (optically flat) quasars, 2) intrinsically red (optically steep) quasars, and 3) the 273 (6%) of our quasars whose continua are inconsistent with a single power-law and appear redder due to SMC-like dust reddening rather than synchrotron emission. The color distribution suggests that the population of moderately dust reddened broad-line quasars is smaller than that of unobscured quasars, but we estimate that a further 10% of the luminous quasar population is missing from the SDSS sample because of dust extinction with E(B-V)<0.5. We also investigate the emission and absorption line properties of these quasars as a function of color with regard to Boroson & Green type eigenvectors. Intrinsically red (optically steep) quasars tend to have narrower Balmer lines and weaker CIV, CIII], HeII and 3000A bump emission as compared with bluer (optically flatter) quasars. The change in strength of the 3000A bump appears to be dominated by the Balmer continuum and not by FeII emission. The dust reddened quasars have even narrower Balmer lines and weaker 3000A bumps, in addition to having considerably larger equivalent widths of [OII] and [OIII] emission. The fraction of broad absorption line quasars (BALQSOs) increases from ~3.4% for the bluest quasars to perhaps as large as 20% for the dust reddened quasars, but the intrinsic color distribution is affected by dust reddening. (abridged)

We examine the properties of the host galaxies of 22,623 narrow-line AGN with 0.02<z<0.3 selected from a complete sample of 122,808 galaxies from the Sloan Digital Sky Survey. We focus on the luminosity of the [OIII]$\lambda$5007 emission line as a tracer of the strength of activity in the nucleus. We study how AGN host properties compare to those of normal galaxies and how they depend on L[OIII]. We find that AGN of all luminosities reside almost exclusively in massive galaxies and have distributions of sizes, stellar surface mass densities and concentrations that are similar to those of ordinary early-type galaxies in our sample. The host galaxies of low-luminosity AGN have stellar populations similar to normal early-types. The hosts of high- luminosity AGN have much younger mean stellar ages. The young stars are not preferentially located near the nucleus of the galaxy, but are spread out over scales of at least several kiloparsecs. A significant fraction of high- luminosity AGN have strong H$\delta$ absorption-line equivalent widths, indicating that they experienced a burst of star formation in the recent past. We have also examined the stellar populations of the host galaxies of a sample of broad-line AGN. We conclude that there is no significant difference in stellar content between type 2 Seyfert hosts and QSOs with the same [OIII] luminosity and redshift. This establishes that a young stellar population is a general property of AGN with high [OIII] luminosities.

The potential of the Sloan Digital Sky Survey for wide-field variability studies is illustrated using multi-epoch observations for 3,000,000 point sources observed in 700 deg2 of sky, with time spans ranging from 3 hours to 3 years. These repeated observations of the same sources demonstrate that SDSS delivers ~0.02 mag photometry with well behaved and understood errors. We show that quasars dominate optically faint (r > 18) point sources that are variable on time scales longer than a few months, while for shorter time scales, and at bright magnitudes, most variable sources are stars.

[Abridged] We present here a new and homogeneous sample of 3340 galaxies selected from the Sloan Digital Sky Survey (SDSS) based solely on the observed strength of their Hdelta absorption line. These galaxies are commonly known as ``post-starburst'' or ``E+A'' galaxies, and the study of these galaxies has been severely hampered by the lack of a large, statistical sample of such galaxies. In this paper, we rectify this problem by selecting a sample of galaxies which possess an absorption Hdelta equivalent width of EW(Hdelta_max) - Delta EW(Hdelta_max) > 4A from 106682 galaxies in the SDSS. We have performed extensive tests on our catalog including comparing different methodologies of measuring the Hdelta absorption and studying the effects of stellar absorption, dust extinction, emission-filling and measurement error. The measured abundance of our Hdelta-selected (HDS) galaxies is 2.6 +/- 0.1% of all galaxies within a volume-limited sample of 0.05<z<0.1 and M(r*)<-20.5, which is consistent with previous studies of such galaxies in the literature. We find that only 25 of our HDS galaxies in this volume-limited sample (3.5+/-0.7%) show no evidence for OII and Halpha emission, thus indicating that true E+A (or k+a) galaxies are extremely rare objects at low redshift, i.e., only 0.09+/-0.02% of all galaxies in this volume-limited sample are true E+A galaxies. In contrast, 89+/-5% of our HDS galaxies in the volume-limited sample have significant detections of the OII and Halpha emission lines. We find 27 galaxies in our volume-limited HDS sample that possess no detectable OII emission, but do however possess detectable Halpha emission. These galaxies may be dusty star-forming galaxies. We provide the community with this new catalog of Hdelta-selected galaxies to aid in the understanding of these galaxies.

We present the discovery of three new quasars at z>6 in 1300 deg^2 of SDSS imaging data, J114816.64+525150.3 (z=6.43), J104845.05+463718.3 (z=6.23) and J163033.90+401209.6 (z=6.05). The first two objects have weak Ly alpha emission lines; their redshifts are determined from the positions of the Lyman break. They are only accurate to 0.05 and could be affected by the presence of broad absorption line systems. The last object has a Ly alpha strength more typical of lower redshift quasars. Based on a sample of six quasars at z>5.7 that cover 2870 deg^2 presented in this paper and in Paper I, we estimate the comoving density of luminous quasars at z 6 and M_{1450} < -26.8 to be (8 +/- 3)x10^{-10} Mpc^{-3} (for H_0 = 50 km/s/Mpc, Omega = 1). HST imaging of two z>5.7 quasars and high-resolution ground-based images (seeing 0.4'') of three additional z>5.7 quasars show that none of them is gravitationally lensed. The luminosity distribution of the high-redshfit quasar sample suggests the bright end slope of the quasar luminosity function at z 6 is shallower than Psi L^{-3.5} (2-sigma), consistent with the absence of strongly lensed objects.

We present a determination of the `Cosmic Optical Spectrum' of the Universe, i.e. the ensemble emission from galaxies, as determined from the red-selected Sloan Digital Sky Survey main galaxy sample and compare with previous results of the blue-selected 2dF Galaxy Redshift Survey. Broadly we find good agreement in both the spectrum and the derived star-formation histories. If we use a power-law star-formation history model where star-formation rate $\propto (1+z)^\beta$ out to z=1, then we find that $\beta$ of 2 to 3 is still the most likely model and there is no evidence for current surveys missing large amounts of star formation at high redshift. In particular `Fossil Cosmology' of the local universe gives measures of star-formation history which are consistent with direct observations at high redshift. Using the photometry of SDSS we are able to derive the cosmic spectrum in absolute units (i.e.$ W \AA$^{-1}$ Mpc$^{-3}$) at 2--5\AA resolution and find good agreement with published broad-band luminosity densities. For a Salpeter IMF the best fit stellar mass/light ratio is 3.7--7.5 $\Msun/\Lsun$ in the r-band (corresponding to $\omstars h = 0.0025$--0.0055) and from both the stellar emission history and the H$\alpha$ luminosity density independently we find a cosmological star-formation rate of 0.03--0.04 h $\Msun$ yr$^{-1}$ Mpc$^{-3}$ today.

We combine Sloan Digital Sky Survey spectra of 22,000 luminous, red, bulge-dominated galaxies to get high S/N average spectra in the rest-frame optical and ultraviolet (2600A to 7000A). The average spectra of these massive, quiescent galaxies are early-type with weak emission lines and with absorption lines indicating an apparent excess of alpha elements over solar abundance ratios. We make average spectra of subsamples selected by luminosity, environment and redshift. The average spectra are remarkable in their similarity. What variations do exist in the average spectra as a function of luminosity and environment are found to form a nearly one-parameter family in spectrum space. We present a high signal-to-noise ratio spectrum of the variation. We measure the properties of the variation with a modified version of the Lick index system and compare to model spectra from stellar population syntheses. The variation may be a combination of age and chemical abundance differences, but the conservative conclusion is that the quality of the data considerably exceeds the current state of the models.

We present a method for photometric selection of metal-poor halo giants from the imaging data of the Sloan Digital Sky Survey (SDSS). These stars are offset from the stellar locus in the (g-r) vs. (u-g) color-color diagram. Based on a sample of 29 candidates for which spectra were taken, we derive a selection efficiency of the order of 50%, for stars brighter than $r \sim 17^m$. The candidates selected in 400 deg$^2$ of sky from the SDSS Early Data Release trace the known halo structures (tidal streams from the Sagittarius dwarf galaxy, the Draco dwarf spheroidal galaxy), indicating that such a color-selected sample can be used to study the halo structure even without spectroscopic information. This method, and supplemental techniques for selecting halo stars, such as RR Lyrae stars and other blue horizontal branch stars, can produce an unprecedented three-dimensional map of the Galactic halo based on the SDSS imaging survey.

The astrometric calibration of the Sloan Digital Sky Survey is described. For point sources brighter than r ~ 20 the astrometric accuracy is 45 milliarcseconds (mas) rms per coordinate when reduced against the USNO CCD Astrograph Catalog, and 75 mas rms when reduced against Tycho-2, with an additional 20 - 30 mas systematic error in both cases. The rms errors are dominated by anomalous refraction and random errors in the primary reference catalogs. The relative astrometric accuracy between the r filter and each of the other filters (u g i z) is 25 - 35 mas rms. At the survey limit (r ~ 22), the astrometric accuracy is limited by photon statistics to approximately 100 mas rms for typical seeing. Anomalous refraction is shown to contain components correlated over two or more degrees on the sky.

We develop a new method to constrain the star formation histories, dust attenuation and stellar masses of galaxies. It is based on two stellar absorption line indices, the 4000 Angstrom break strength and the Balmer absorption line index Hdelta_A. Together, these indices allow us to constrain the mean stellar ages of galaxies and the fractional stellar mass formed in bursts over the past few Gyr. A comparison with broad band photometry then yields estimates of dust attenuation and of stellar mass. We generate a large library of Monte Carlo realizations of different star formation histories, including starbursts of varying strength and a range of metallicities. We use this library to generate median likelihood estimates of burst mass fractions, dust attenuation strengths, stellar masses and stellar mass-to-light ratios for a sample of 122,208 galaxies drawn from the Sloan Digital Sky Survey. The typical 95% confidence range in our estimated stellar masses is +-40%. We study how the stellar mass-to-light ratios of galaxies vary as a function of absolute magnitude, concentration index and photometric pass-band and how dust attenuation varies as a function of absolute magnitude and 4000 Angstrom break strength. We also calculate how the total stellar mass of the present Universe is distributed over galaxies as a function of their mass, size, concentration, colour, burst mass fraction and surface mass density. We find that most of the stellar mass in the local Universe resides in galaxies that have stellar masses \~5\times 10^10 M_sol, half light radii ~3 kpc, and half- light surface mass densities ~10^9 M_sol/kpc^2. The distribution of D(4000) is strongly bimodal, showing a clear division between galaxies dominated by old stellar populations and galaxies with more recent star formation.

Using photometry and spectroscopy of 144,609 galaxies from the Sloan Digital Sky Survey, we present bivariate distributions of pairs of seven galaxy properties: four optical colors, surface brightness, radial profile shape as measured by the Sersic index, and absolute magnitude. In addition, we present the dependence of local galaxy density (smoothed on 8 h^{-1} Mpc scales) on all of these properties. Several classic, well-known relations among galaxy properties are evident at extremely high signal-to-noise ratio: the color-color relations of galaxies, the color-magnitude relations, the magnitude-surface brightness relation, and the dependence of density on color and absolute magnitude. We show that most of the i-band luminosity density in the universe is in the absolute magnitude and surface brightness ranges used. Some of the relationships between parameters, in particular the color--magnitude relations, show stronger correlations for exponential galaxies and concentrated galaxies taken separately than for all galaxies taken together. We provide a simple set of fits of the dependence of galaxy properties on luminosity for these two sets of galaxies.

In this paper, we limit the stellar content of 13 high-velocity clouds (HVCs; including 1 compact HVC) using the detection limits of a new survey for resolved Milky Way satellite galaxies in the Sloan Digital Sky Survey Early Data Release (EDR). Our analysis is sensitive to stellar associations within the virial radius of the Milky Way that are up to 50 times fainter than the faintest known Milky Way satellites. Statistically, we find no stellar overdensity associated with any of the clouds. These non-detections suggest lower limits of M_{HI}/L_V that range between ~ 0.2 and 250, assuming cloud distances of 50 and 300 kpc and using fiducial optical scale lengths of 3' and 7'. We explore the implication of these non-detections on the origin of the HVCs.

We discuss optical colors of 10,592 asteroids with known orbits selected from a sample of 58,000 moving objects observed by the Sloan Digital Sky Survey (SDSS). This is more than ten times larger sample that includes both orbital parameters and multi-band photometric measurements than previously available. We confirm that asteroid dynamical families, defined as clusters in orbital parameter space, also strongly segregate in color space. In particular, we demonstrate that the three major asteroid families (Eos, Koronis, and Themis), together with the Vesta family, represent four main asteroid color types. Their distinctive optical colors indicate that the variations in chemical composition within a family are much smaller than the compositional differences between families, and strongly support earlier suggestions that asteroids belonging to a particular family have a common origin. We estimate that over 90% of asteroids belong to families.

We announce the first public release of the SDSS Moving Object Catalog, with SDSS observations for 58,117 asteroids. The catalog lists astrometric and photometric data for moving objects observed prior to Dec 15, 2001, and also includes orbital elements for 10,592 previously known objects. We analyze the correlation between the orbital parameters and optical colors for the known objects, and confirm that asteroid dynamical families, defined as clusters in orbital parameter space, also strongly segregate in color space. Their distinctive optical colors indicate that the variations in chemical composition within a family are much smaller than the compositional differences between families, and strongly support earlier suggestions that asteroids belonging to a particular family have a common origin.

Most dusty winds are described by a set of similarity functions of a single independent variable that can be chosen as the overall optical depth at visual. The self-similarity implies general scaling relations among the system parameters, in agreement with observations. Dust drift through the gas has a major impact on the structure of most winds.

We have investigated the mass loss processes in Long Period Variables, by combining radiative transfer modeling with mid-IR imaging and spectroscopic observations. We find a correlation between variability type and temporal variations in the mass loss rate. Mira variables are more likely to maintain constant mass loss rates over long period of time. Semiregular and Irregular variables are characterized by frequent interruptions in their mass loss rate, on timescales of ~100 yr. High resolution imaging of individual Long Period Variable sources show departures from spherical symmetry. We present the case of o Cet (Mira), in which the presence of a low mass companion plays an active role in shaping the circumstellar environment.

We positionally correlate asteroids from existing catalogs with a sample of $\about$18,000 asteroids detected by the Sloan Digital Sky Survey (SDSS, Ivezi\'{c} {\em et al.} 2001). We find 2641 unique matches, which represent the largest sample of asteroids with both accurate multi-color photometry and known orbital parameters. The matched objects are predominantly bright, and demonstrate that the SDSS photometric pipeline recovers \about90% of the known asteroids in the observed region. For the recovered asteroids we find a large offset (\about 0.4 mag) between Johnson V magnitudes derived from SDSS photometry and the predicted catalog-based visual magnitudes. This offset varies with the asteroid color from 0.34 mag for blue asteroids to 0.44 mag for red asteroids, and is probably caused by the use of unfiltered CCD observations in the majority of recent asteroid surveys. This systematic photometric error leads to an overestimate of the number of asteroids brighter than a given absolute magnitude limit by a factor of \about 1.7. The distribution of the matched asteroids in orbital parameter space indicates strong color segregation. We confirm that some families are dominated by a single asteroid type (e.g. the Koronis family by red asteroids and the Themis family by blue asteroids), while others appear to be a mixture of blue and red objects (e.g. the Nysa/Polana family). Asteroids with the bluest \iz colors, which can be associated with the Vesta family, show particularly striking localization in orbital parameter space.

We study the Galactic distribution of ~10,000 Asymptotic Giant Branch (AGB) stars selected by IRAS colors and variability index. The distance to each star is estimated by assuming a narrow luminosity function and a model-derived bolometric correction. The characteristic AGB star luminosity, L_AGB, is determined from the condition that the highest number density must coincide with the Galactic bulge. Assuming a bulge distance of 8 kpc, we determine L_AGB \~ 3,500 Lo, in close agreement with values obtained for nearby AGB stars using the HIPPARCOS data. We find that there are no statistically significant differences in the Galactic distribution of AGB stars with different IRAS colors, implying a universal density distribution. The direct determination of this distribution shows that it is separable in the radial, R, and vertical, z, directions. Perpendicular to the Galactic plane, the number density of AGB stars is well described by an exponential function with a vertical scale height of 300 pc. In the radial direction the number density of AGB stars is constant up to R ~ 5 kpc, and then it decreases exponentially with a scale length of \~1.6 kpc. This fall-off extends to at least 12 kpc, where the sample becomes too small. The overall normalization implies that there are about 200,000 AGB stars in the Galaxy. We estimate the [25]-[12] color distribution of AGB stars for an unbiased volume-limited sample. By using a model-dependent transformation between the color and mass-loss rate, Mdot, we constrain the time dependence of Mdot. The results suggest that for 10^-6 Mo/yr < Mdot < 10^-5 Mo/yr the mass-loss rate increases exponentially with time. We find only marginal evidence that the mass-loss rate increases with stellar mass.

Unified schemes of active galactic nuclei (AGN) require an obscuring dusty torus around the central source, giving rise to Seyfert 1 line spectrum for pole-on viewing and Seyfert 2 characteristics in edge-on sources. Although the observed IR is in broad agreement with this scheme, the behavior of the 10 micron silicate feature and the width of the far-IR emission peak remained serious problems in all previous modeling efforts. We show that these problems find a natural explanation if the dust is contained in about 5-10 clouds along radial rays through the torus. The spectral energy distributions (SED) of both type 1 and type 2 sources are properly reproduced from different viewpoints of the same object if the visual optical depth of each cloud is larger than about 60 and the clouds' mean free path increases roughly in proportion to radial distance.

We discuss the optical and radio properties of 30,000 FIRST sources positionally associated with an SDSS source in 1230 deg$^2$ of sky. The majority (83%) of the FIRST sources identified with an SDSS source brighter than r=21 are optically resolved. We estimate an upper limit of 5% for the fraction of quasars with broad-band optical colors indistinguishable from those of stars. The distribution of quasars in the radio flux -- optical flux plane supports the existence of the "quasar radio-dichotomy"; 8% of all quasars with i<18.5 are radio-loud and this fraction seems independent of redshift and optical luminosity. The radio-loud quasars have a redder median color by 0.08 mag, and a 3 times larger fraction of objects with red colors. FIRST galaxies represent 5% of all SDSS galaxies with r<17.5, and 1% for r<20, and are dominated by red galaxies. Magnitude and redshift limited samples show that radio galaxies have a different optical luminosity distribution than non-radio galaxies selected by the same criteria; when galaxies are further separated by their colors, this result remains valid for both blue and red galaxies. The distributions of radio-to-optical flux ratio are similar for blue and red galaxies in redshift-limited samples; this similarity implies that the difference in their luminosity functions, and resulting selection effects, are the dominant cause for the preponderance of red radio galaxies in flux-limited samples. We confirm that the AGN-to-starburst galaxy number ratio increases with radio flux, and find that radio emission from AGNs is more concentrated than radio emission from starburst galaxies (abridged).

We describe the algorithm for selecting quasar candidates for optical spectroscopy in the Sloan Digital Sky Survey. Quasar candidates are selected via their non-stellar colors in "ugriz" broad-band photometry, and by matching unresolved sources to the FIRST radio catalogs. The automated algorithm is sensitive to quasars at all redshifts lower than z=5.8. Extended sources are also targeted as low-redshift quasar candidates in order to investigate the evolution of Active Galactic Nuclei (AGN) at the faint end of the luminosity function. Nearly 95% of previously known quasars are recovered (based on 1540 quasars in 446 square degrees). The overall completeness, estimated from simulated quasars, is expected to be over 90%, whereas the overall efficiency (quasars:quasar candidates) is better than 65%. The selection algorithm targets ultraviolet excess quasars to i^*=19.1 and higher-redshift (z>3) quasars to i^*=20.2, yielding approximately 18 candidates per square degree. In addition to selecting ``normal'' quasars, the design of the algorithm makes it sensitive to atypical AGN such as Broad Absorption Line quasars and heavily reddened quasars.

We describe an automated method, the Cut & Enhance method (CE) for detecting clusters of galaxies in multi-color optical imaging surveys. This method uses simple color cuts, combined with a density enhancement algorithm, to up-weight pairs of galaxies that are close in both angular separation and color. The method is semi-parametric since it uses minimal assumptions about cluster properties in order to minimize possible biases. No assumptions are made about the shape of clusters, their radial profile or their luminosity function. The method is successful in finding systems ranging from poor to rich clusters of galaxies, of both regular and irregular shape. We determine the selection function of the CE method via extensive Monte Carlo simulations which use both the real, observed background of galaxies and a randomized background of galaxies. We use position shuffled and color shuffled data to perform the false positive test. We have also visually checked all the clusters detected by the CE method. We apply the CE method to the 350 deg^2 of the SDSS (Sloan Digital Sky Survey) commissioning data and construct a SDSS CE galaxy cluster catalog with an estimated redshift and richness for each cluster. The CE method is compared with other cluster selection methods used on SDSS data such as the Matched Filter (Postman et al. 1996, Kim et al. 2001), maxBCG technique (Annis et al. 2001) and Voronoi Tessellation (Kim et al. 2001). The CE method can be adopted for cluster selection in any multi-color imaging surveys.

We present multi-instrument optical observations of the High Energy Transient Explorer (HETE-2)/Interplanetary Network (IPN) error box of GRB 010921. This event was the first gamma ray burst (GRB) localized by HETE-2 which has resulted in the detection of an optical afterglow. In this paper we report the earliest known observations of the GRB010921 field, taken with the 0.11-m Livermore Optical Transient Imaging System (LOTIS) telescope, and the earliest known detection of the GRB010921 optical afterglow, using the 0.5-m Sloan Digital Sky Survey Photometric Telescope (SDSS PT). Observations with the LOTIS telescope began during a routine sky patrol 52 minutes after the burst. Observations were made with the SDSS PT, the 0.6-m Super-LOTIS telescope, and the 1.34-m Tautenburg Schmidt telescope at 21.3, 21.8, and 37.5 hours after the GRB, respectively. In addition, the host galaxy was observed with the USNOFS 1.0-m telescope 56 days after the burst. We find that at later times (t > 1 day after the burst), the optical afterglow exhibited a power-law decline with a slope of $\alpha = 1.75 \pm 0.28$. However, our earliest observations show that this power-law decline can not have extended to early times (t < 0.035 day).

We identify new structures in the halo of the Milky Way Galaxy from positions, colors and magnitudes of five million stars detected in the Sloan Digital Sky Survey. Most of these stars are within 1.26 degrees of the celestial equator. We present color-magnitude diagrams (CMDs) for stars in two previously discovered, tidally disrupted structures. The CMDs and turnoff colors are consistent with those of the Sagittarius dwarf galaxy, as had been predicted. In one direction, we are even able to detect a clump of red stars, similar to that of the Sagittarius dwarf, from stars spread across 110 square degrees of sky. Focusing on stars with the colors of F turnoff objects, we identify at least five additional overdensities of stars. Four of these may be pieces of the same halo structure, which would cover a region of the sky at least 40 degrees in diameter, at a distance of 11 kpc from the Sun (18 kpc from the center of the Galaxy). The turnoff is significantly bluer than that of thick disk stars, and closer to the Galactic plane than a power-law spheroid. We suggest two models to explain this new structure. One possibility is that this new structure could be a new dwarf satellite of the Milky Way, hidden in the Galactic plane, and in the process of being tidally disrupted. The other possibility is that it could be part of a disk-like distribution of stars which is metal-poor, with a scale height of approximately 2 kpc and a scale length of approximately 10 kpc. The fifth overdensity, which is 20 kpc away, is some distance from the Sagittarius dwarf streamer orbit and is not associated with any known structure in the Galactic plane. It is likely that there are many smaller streams of stars in the Galactic halo.

We present the detection limits of a new survey for resolved low surface brightness satellite galaxies to the Milky Way, based on the Sloan Digital Sky Survey (SDSS). Our survey exploits SDSS's major strengths (multi-color photometry, depth, large-scale, and uniformity) by combining filter smoothing with limits in both magnitude and color space to search for low surface brightness galaxies and stellar counterparts to the Compact High Velocity Clouds out to the Milky Way's virial radius (350 kpc). Our calculated detection limits for a purely old stellar population vary with galaxy size and distance between mu_{V,0} = 26.7 and 30.1 mag/sq ''. These limits will allow us to detect systems whose surface brightnesses are .5-3.9 mag/sq '' fainter than Sextans, the lowest surface brightness Local Group member known. Our survey not only is sensitive to lower surface brightness stellar populations than possible with previous Local Group surveys, but will also allow us to make an unbiased and well defined assessment of the completeness of the observed Local Group galaxy luminosity function, so that we may compare the results with the predictions of various structure formation scenarios.

We positionally match sources observed by the Sloan Digital Sky Survey (SDSS), the Two Micron All Sky Survey (2MASS), and the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. Practically all 2MASS sources are matched to an SDSS source within 2 arcsec; ~11% of them are optically resolved galaxies and the rest are dominated by stars. About 1/3 of FIRST sources are matched to an SDSS source within 2 arcsec; ~80% of these are galaxies and the rest are dominated by quasars. Based on these results, we project that by the completion of these surveys the matched samples will include about 10^7 stars and 10^6 galaxies observed by both SDSS and 2MASS, and about 250,000 galaxies and 50,000 quasars observed by both SDSS and FIRST. Here we present a preliminary analysis of the optical, infrared and radio properties for the extragalactic sources from the matched samples. In particular, we find that the fraction of quasars with stellar colors missed by the SDSS spectroscopic survey is probably not larger than ~10%, and that the optical colors of radio-loud quasars are ~0.05 mag. redder (with 4-sigma significance) than the colors of radio-quiet quasars.

We identify close pairs of galaxies from 278 deg^2 of Sloan Digital Sky Survey commissioning imaging data. The pairs are drawn from a sample of 330,041 galaxies with 18 < r^* < 20. We determine the angular correlation function of galaxy pairs, and find it to be stronger than the correlation function of single galaxies by a factor of 2.9 +/- 0.4. The two correlation functions have the same logarithmic slope of 0.77. We invert Limber's equation to estimate the three-dimensional correlation functions; we find clustering lengths of r_0= 4.2 +/- 0.4 h^{-1} Mpc for galaxies and 7.8 +/- 0.7 h^{-1} Mpc for galaxy pairs. These results agree well with the global richness dependence of the correlation functions of galaxy systems.

A sample of nearly 9000 early-type galaxies, in the redshift range 0.01<z<0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The sample was used to study how early-type galaxy observables, including luminosity L, effective radius R, surface brightness I, color, and velocity dispersion V, are correlated with one another. Measurement biases are understood with mock catalogs which reproduce all of the observed scaling relations. At any given redshift, the intrinsic distribution of luminosities, sizes and velocity dispersions in our sample are all approximately Gaussian. In the r* band L ~ V^3.91, L ~ R^1.58, R ~ I^(-0.75), and the Fundamental Plane relation is R ~ V^(1.49) I^(-0.75). These relations are approximately the same in the g*, i* and z* bands. At fixed luminosity, the mass-to-light ratio scales as M/L ~ L^0.14. The g*-r* color scales as V^0.25. Color also correlates with magnitude and size, but these correlations are entirely due to the L-V and R-V relations. Chemical evolution and star formation histories are investigated using co-added spectra of similar objects in our sample. Chemical abundances correlate primarily with velocity dispersion. At fixed V, the higher redshift population is bluer, is weaker in Mg2, and is stronger in Hbeta than the population nearby. In addition, the population at higher redshifts is slightly more luminous. These differences are consistent with that of a passively evolving population which formed the bulk of its stars about 9 Gyrs ago. The Fundamental Plane suggests that galaxies in dense regions are slightly different from those in less dense regions, but the co-added spectra and color--magnitude relations show no statistically significant dependence on environment.

Low dispersion optical spectra have been obtained with the Hobby-Eberly Telescope of 22 very red objects found in early imaging data from the Sloan Digital Sky Survey. The objects are assigned spectral types on the 2MASS system (Kirkpatrick et al. 1999) and are found to range from late M to late L. The red- and near-infrared colors from SDSS and 2MASS correlate closely with each other, and most of the colors are closely related to spectral type in this range; the exception is the (i^* - z^*) color, which appears to be independent of spectral type between about M7 and L4. The spectra suggest that this independence is due to the disappearance of the TiO and VO absorption in the i-band for later spectral types; to the presence of strong Na I and K I absorption in the i-band; and to the gradual disappearance of the 8400 Angstrom absorption of TiO and FeH in the z-band.

The commissioning year of the Sloan Digital Sky Survey has demonstrated that many cataclysmic variables have been missed in previous surveys with brighter limits. We report the identification of 22 cataclysmic variables, of which 19 are new discoveries and 3 are known systems (SW UMa, BH Lyn and Vir4). A compendium of positions, colors and characteristics of these systems obtained from the SDSS photometry and spectroscopy is presented along with data obtained during follow-up studies with the Apache Point Observatory (APO) and Manastash Ridge Observatory (MRO) telescopes. We have determined orbital periods for 3 of the new systems: two show dwarf nova outbursts, and the third is a likely magnetic system with eclipses of its region of line emission. Based on these results, we expect the completed survey to locate at least 400 new CVs. Most of these will be faint systems with low accretion rates that will provide new constraints on binary evolution models.

We present a comparison of three cluster finding algorithms from imaging data using Monte Carlo simulations of clusters embedded in a 25 deg^2 region of Sloan Digital Sky Survey (SDSS) imaging data: the Matched Filter (MF; Postman et al. 1996), the Adaptive Matched Filter (AMF; Kepner et al. 1999) and a color-magnitude filtered Voronoi Tessellation Technique (VTT). Among the two matched filters, we find that the MF is more efficient in detecting faint clusters, whereas the AMF evaluates the redshifts and richnesses more accurately, therefore suggesting a hybrid method (HMF) that combines the two. The HMF outperforms the VTT when using a background that is uniform, but it is more sensitive to the presence of a non-uniform galaxy background than is the VTT; this is due to the assumption of a uniform background in the HMF model. We thus find that for the detection thresholds we determine to be appropriate for the SDSS data, the performance of both algorithms are similar; we present the selection function for each method evaluated with these thresholds as a function of redshift and richness. For simulated clusters generated with a Schechter luminosity function (M_r^* = -21.5 and alpha = -1.1) both algorithms are complete for Abell richness >= 1 clusters up to z ~ 0.4 for a sample magnitude limited to r = 21. While the cluster parameter evaluation shows a mild correlation with the local background density, the detection efficiency is not significantly affected by the background fluctuations, unlike previous shallower surveys.

We present moderate resolution Keck spectroscopy of quasars at z=5.82, 5.99 and 6.28, discovered by the Sloan Digital Sky Survey (SDSS). We find that the Ly Alpha absorption in the spectra of these quasars evolves strongly with redshift. To z~5.7, the Ly Alpha absorption evolves as expected from an extrapolation from lower redshifts. However, in the highest redshift object, SDSSp J103027.10+052455.0 (z=6.28), the average transmitted flux is 0.0038+-0.0026 times that of the continuum level over 8450 A < lambda < 8710 A (5.95<z(abs)<6.16), consistent with zero flux. Thus the flux level drops by a factor of >150, and is consistent with zero flux in the Ly Alpha forest region immediately blueward of the Ly Alpha emission line, compared with a drop by a factor of ~10 at z(abs)~5.3. A similar break is seen at Ly Beta; because of the decreased oscillator strength of this transition, this allows us to put a considerably stronger limit, tau(eff) > 20, on the optical depth to Ly Alpha absorption at z=6. This is a clear detection of a complete Gunn-Peterson trough, caused by neutral hydrogen in the intergalactic medium. Even a small neutral hydrogen fraction in the intergalactic medium would result in an undetectable flux in the Ly Alpha forest region. Therefore, the existence of the Gunn-Peterson trough by itself does not indicate that the quasar is observed prior to the reionization epoch. However, the fast evolution of the mean absorption in these high-redshift quasars suggests that the mean ionizing background along the line of sight to this quasar has declined significantly from z~5 to 6, and the universe is approaching the reionization epoch at z~6.

We present the results from a survey of i-dropout objects selected from ~1550 deg^2 of multicolor imaging data from the Sloan Digital Sky Survey, to search for luminous quasars at z>5.8. Objects with i*-z*>2.2 and z*<20.2 are selected, and follow-up J band photometry is used to separate L and T type cool dwarfs from high-redshift quasars. We describe the discovery of three new quasars, at z=5.82, 5.99 and 6.28, respectively. Their spectra show strong and broad Ly alpha+NV emission lines, and very strong Ly alpha absorption, with a mean continuum decrement D_A > 0.90. The ARC 3.5m spectrum of the z=6.28 quasar shows that over a range of 300 A immediately blueward of the Ly alpha emission, the average transmitted flux is only 0.003 +/-0.020 times that of the continuum level, consistent with zero flux, and suggesting a tentative detection of the complete Gunn-Peterson trough. The existence of strong metal lines suggests early chemical enrichment in the quasar enviornment. The three new objects, together with the previously published z=5.8 quasar form a complete color-selected flux-limited sample at z>5.8. We estimate that at $z=6$, the comoving density of luminous quasars at M_1450 < -26.89 (h=0.5, Omega=1)is 1.1x10^-9 Mpc^-3. This is a factor of ~2 lower than that at z~5, and is consistent with an extrapolation of the observed quasar evolution at low-z. We discuss the contribution of quasars to the ionizing background at z~6. The luminous quasars discussed in the paper have central black hole masses of several times 10^9 M_sun by the Eddington argument. Their observed space density provides a sensitive test of models of quasar and galaxy formation at high redshift. (Abridged)

We address the dusty wind problem, from the point where dust formation has been completed and outward. Given grain properties, both radiative transfer and hydrodynamics components of the problem are fully defined by four additional input parameters. The wind radiative emission and the shape of its velocity profile are both independent of the actual magnitude of the velocity and are determined by just three dimensionless free parameters. Of the three, only one is always significant---for most of phase space the solution is described by a set of similarity functions of a single independent variable, which can be chosen as the overall optical depth at visual \tV. The self-similarity implies general scaling relations among mass loss rate (\Mdot), luminosity ($L$) and terminal velocity (\vf). Systems with different \Mdot, $L$ and \vf but the same combination $\Mdot/L^{3/4}$ necessarily have also the same $\Mdot\vf/L$. ... Eliminating \tV produces $\vf^3 = A \Mdot(1 + B \Mdot^{4/3}/L)^{-1.5}$, where $A$ and $B$ are coefficients that contain the only dependence of this universal correlation on chemical composition. At a given $L$, the maximal velocity of a dusty wind is $v\sub{max} \propto L^{1/4}$ attained at $\Mdot \propto L^{3/4}$, with proportionality coefficients derived from $A$ and $B$.

We describe the target selection and resulting properties of a spectroscopic sample of luminous, red galaxies (LRG) from the imaging data of the Sloan Digital Sky Survey (SDSS). These galaxies are selected on the basis of color and magnitude to yield a sample of luminous, intrinsically red galaxies that extends fainter and further than the main flux-limited portion of the SDSS galaxy spectroscopic sample. The sample is designed to impose a passively-evolving luminosity and rest-frame color cut to a redshift of 0.38. Additional, yet more luminous, red galaxies are included to a redshift of 0.5. Approximately 12 of these galaxies per square degree are targeted for spectroscopy, so the sample will number over 100,000 with the full survey. SDSS commissioning data indicate that the algorithm efficiently selects luminous (M_g=-21.4), red galaxies, that the spectroscopic success rate is very high, and that the resulting set of galaxies is approximately volume-limited out to z=0.38. When the SDSS is complete, the LRG spectroscopic sample will fill over 1h^-3 Gpc^3 with an approximately homogeneous population of galaxies and will therefore be well suited to studies of large-scale structure and clusters out to z=0.5.

We present new measurements of scaling laws relating the luminosity of galaxies to the amplitude and shape of their dark matter halos. Early imaging and spectroscopic data from the Sloan Digital Sky Survey are used to make weak lensing measurements of the surface mass density contrast Delta Sigma_+ around classes of lens objects. This surface mass density contrast as a function of radius is a measure of the galaxy-mass correlation function (GMCF). Because spectroscopic redshifts are available for all lens objects, the mass and distance scales are well constrained. The GMCF measured around ~31,000 lenses is well fit by a power law of the form Delta Sigma_+ = (2.5+0.7-0.6) (R/1 Mpc)^{-0.8+-0.2} h M_sun pc^-2. We compare this GMCF to galaxy luminosity, type, and environment, and find that it varies strongly with all three. We quantify these variations by comparing the normalization of a fit to the inner 260 h^-1 kpc, M_260, to the galaxy luminosity. While M_260 is not strongly related to luminosity in bluest band u', there is a simple, linear relation between M_260 and luminosity in redder bands (g', r', i', and z'). We test the universality of these mass-to-light scalings by independently measuring them for spiral and elliptical galaxies,and for galaxies in a variety of environments. We find remarkable consistency in these determinations in the red bands, especially i' and z'. This consistency across a wide range of systems suggests that the measured scaling represents an excellent cosmic average, and that the integrated star formation history of galaxies is strongly related to the dark matter environments in which they form.

We study the optical colors of 147,920 galaxies brighter than g* = 21, observed in five bands by the Sloan Digital Sky Survey (SDSS) over ~100 sq. deg. of high Galactic latitude sky along the Celestial Equator. The distribution of galaxies in the g*-r* vs. u*-g* color--color diagram is strongly bimodal, with an optimal color separator of u*-r* = 2.22. We use visual morphology and spectral classification of subsamples of 287 and 500 galaxies respectively, to show that the two peaks correspond roughly to early (E, S0, Sa) and late (Sb, Sc, Irr) type galaxies, as expected from their different stellar populations. We also find that the colors of galaxies are correlated with their radial profiles, as measured by the concentration index and by the likelihoods of exponential and de Vaucouleurs' profile fits. While it is well known that late type galaxies are bluer than early type galaxies, this is the first detection of a local minimum in their color distribution. In all SDSS bands, the counts vs. apparent magnitude relations for the two color types are significantly different, and indicate that the fraction of blue galaxies increases towards the faint end.

We present 13 Broad Absorption Line (BAL) quasars, including 12 new objects, which were identified in the Sloan Digital Sky Survey (SDSS) and matched within 2'' to sources in the FIRST radio survey catalog. The surface density of this sample of radio-detected BAL quasars is 4.5 +- 1.2 per 100 deg^2, i.e. approximately 4 times larger than previously found by the shallower FIRST Bright Quasar Survey (FBQS). A majority of these radio-detected BAL quasars are moderately radio-loud objects. The fraction of BAL quasars in the entire radio quasar sample, 4.8 +- 1.3 %, is comparable to the fraction of BAL quasars among the SDSS optical quasar sample (ignoring selection effects). We estimate that the true fraction of BAL quasars (mostly HiBALs) in the radio sample is 9.2 +- 2.6 % once selection effects are accounted for. We caution that the absorption troughs of 4 of the 13 radio-detected quasars considered do not strictly satisfy the standard BALnicity criterion. One or possibly two of the new radio-detected BAL quasars are of the rare ``iron LoBAL'' type. BAL quasars are generally redder than the median SDSS quasar at the same redshift.

The discovery of an optical counterpart to GRB010222 (detected by BeppoSAX; Piro 2001) was announced 4.4 hrs after the burst by Henden (2001a). The Sloan Digital Sky Survey's 0.5m photometric telescope (PT) and 2.5m survey telescope were used to observe the afterglow of GRB010222 starting 4.8 hours after the GRB. The 0.5m PT observed the afterglow in five, 300 sec g' band exposures over the course of half an hour, measuring a temporal decay rate in this short period of F_nu \propto t^{-1.0+/-0.5}. The 2.5m camera imaged the counterpart nearly simultaneously in five filters (u' g' r' i' z'), with r' = 18.74+/-0.02 at 12:10 UT. These multicolor observations, corrected for reddening and the afterglow's temporal decay, are well fit by the power-law F_nu \propto nu^{-0.90+/-0.03} with the exception of the u' band UV flux which is 20% below this slope. We examine possible interpretations of this spectral shape, including source extinction in a star forming region.

We present bright galaxy number counts in five broad bands ($u', g', r', i', z'$) from imaging data taken during the commissioning phase of the Sloan Digital Sky Survey (SDSS). The counts are derived from two independent stripes of imaging scans along the Celestial Equator, one each toward the North and the South Galactic cap, covering about 230 and 210 square degrees, respectively. A careful study is made to verify the reliability of the photometric catalog. For galaxies brighter than $r^* = 16$, the catalog produced by automated software is examined against eye inspection of all objects. Statistically meaningful results on the galaxy counts are obtained in the magnitude range $12 \le r^* \le 21$, using a sample of 900,000 galaxies. The counts from the two stripes differ by about 30% at magnitudes brighter than $r^*= 15.5$, consistent with a local $2\sigma$ fluctuation due to large scale structure in the galaxy distribution. The shape of the number counts-magnitude relation brighter than $r^* = 16$ is well characterized by $N \propto 10^{0.6m}$, the relation expected for a homogeneous galaxy distribution in a ``Euclidean'' universe. In the magnitude range $ 16 < r^* < 21$, the galaxy counts from both stripes agree very well, and follow the prediction of the no-evolution model, although the data do not exclude a small amount of evolution. We use empirically determined color transformations to derive the galaxy number counts in the $B$ and $I_{814}$ bands. We compute the luminosity density of the universe at zero redshift in the five SDSS bands and in the $B$ band. We find ${\cal L}_{B} = 2.4 \pm 0.4 \times 10^8L_\odot h $Mpc$^{-3}$, for a reasonably wide range of parameters of the Schechter luminosity function in the $B$ band.

We discuss measurements of the properties of about 10,000 asteroids detected in 500 deg2 of sky in the Sloan Digital Sky Survey (SDSS) commissioning data. The moving objects are detected in the magnitude range 14 < r < 21.5, with a baseline of 5 minutes. Extensive tests show that the sample is at least 98% complete, with the contamination rate of less than 3%. We find that the size distribution of asteroids resembles a broken power-law, independent of the heliocentric distance: D^{-2.3} for 0.4 km < D < 5 km, and D^{-4} for 5 km < D < 40 km. As a consequence of this break, the number of asteroids with r < 21.5 is ten times smaller than predicted by extrapolating the power-law relation observed for brighter asteroids (r < 18). The observed counts imply that there are about 530,000 objects with D>1 km in the asteroid belt, or about four times less than previous estimates. The distribution of main belt asteroids in the 4-dimensional SDSS color space is bimodal, and the two groups can be associated with S (rocky) and C (carbonaceous) asteroids. A strong bimodality is also seen in the heliocentric distribution of asteroids and suggests the existence of two distinct belts: the inner rocky belt, about 1 AU wide (FWHM) and centered at R~2.8 AU, and the outer carbonaceous belt, about 0.5 AU wide and centered at R~3.2 AU. The colors of Hungarias, Mars crossers, and near-Earth objects are more similar to the C-type than to S-type asteroids, suggesting that they originate in the outer belt. (abridged).

We summarise the properties of the Sloan Digital Sky Survey (SDSS) project, discuss our software infrastructure, and outline the architecture of the SDSS image processing pipelines. We then discuss two of the algorithms used in the SDSS image processing; the KL-transform based modelling of the spatial variation of the PSF, and the use of galaxy models in star/galaxy separation. We conclude with the first author's personal opinions on the challenges that the astronomical community faces with major software projects.

We present a lensing study of 42 galaxy clusters imaged in Sloan Digital Sky Survey (SDSS) commissioning data. Cluster candidates are selected optically from SDSS imaging data and confirmed for this study by matching to X-ray sources found independently in the ROSAT all sky survey (RASS). Five color SDSS photometry is used to make accurate photometric redshift estimates that are used to rescale and combine the lensing measurements. The mean shear from these clusters is detected to 2 h-1 Mpc at the 7-sigma level, corresponding to a mass within that radius of 4.2 +/- 0.6 x 10^14 h-1 M_sun. The shear profile is well fit by a power law with index -0.9 +/- 0.3, consistent with that of an isothermal density profile. This paper demonstrates our ability to measure ensemble cluster masses from SDSS imaging data.

We analyze the differences in infrared circumstellar dust emission between oxygen rich Mira and non-Mira stars, and find that they are statistically significant. In particular, we find that these stars segregate in the K-[12] vs. [12]-[25] color-color diagram, and have distinct properties of the IRAS LRS spectra, including the peak position of the silicate emission feature. We show that the infrared emission from the majority of non-Mira stars cannot be explained within the context of standard steady-state outflow models. The models can be altered to fit the data for non-Mira stars by postulating non-standard optical properties for silicate grains, or by assuming that the dust temperature at the inner envelope radius is significantly lower (300-400 K) than typical silicate grain condensation temperatures (800-1000 K). We argue that the latter is more probable and provide detailed model fits to the IRAS LRS spectra for 342 stars. These fits imply that 2/3 of non-Mira stars and 1/3 of Mira stars do not have hot dust ($>$ 500 K) in their envelopes. The absence of hot dust can be interpreted as a recent (order of 100 yr) decrease in the mass-loss rate. The distribution of best-fit model parameters agrees with this interpretation and strongly suggests that the mass loss resumes on similar time scales. Such a possibility appears to be supported by a number of spatially resolved observations (e.g. recent HST images of the multiple shells in the Egg Nebula) and is consistent with new dynamical models for mass loss on the Asymptotic Giant Branch.

For about 11,000 stars observed in the HIPPARCOS Survey and detected by IRAS we calculate bolometric luminosities by integrating their spectral energy distributions from the B band to far-IR wavelengths. We present an analysis of the dependence of dust emission on spectral type and correlations between the luminosity and dust emission for about 1000 sources with the best data (parallax error less than 30%, error in luminosity of about 50% or better). This subsample includes stars of all spectral types and is dominated by K and M giants. We use the IRAS [25]-[12] color to select stars with emission from circumstellar dust and show that they are found throughout the Hertzsprung-Russell diagram, including on the main sequence. Clear evidence is found that M giants with dust emission have luminosities about 3 times larger (about 3000 Lsun) than their counterparts without dust, and that mass loss on the asymptotic giant branch for both M and C stars requires a minimum luminosity of order 2000 Lsun. Above this threshold the mass-loss rate seems to be independent of, or only weakly dependent on, luminosity. We also show that the mass-loss rate for these stars is larger than the core mass growth rate, indicating that their evolution is dominated by mass loss.

Early data taken during commissioning of the SDSS have resulted in the discovery of a very cool white dwarf. It appears to have stronger collision induced absorption from molecular hydrogen than any other known white dwarf, suggesting it has a cooler temperature than any other. While its distance is presently unknown, it has a surprisingly small proper motion, making it unlikely to be a halo star. An analysis of white dwarf cooling times suggests that this object may be a low-mass star with a helium core. The SDSS imaging and spectroscopy also recovered LHS 3250, the coolest previously known white dwarf, indicating that the SDSS will be an effective tool for identifying these extreme objects.

We report the discovery of two well-defined tidal tails emerging from the sparse remote globular cluster Palomar 5. These tails stretch out symmetrically to both sides of the cluster in the direction of constant Galactic latitude and subtend an angle of 2.6 degrees on the sky. The tails have been detected in commissioning data of the Sloan Digital Sky Survey (SDSS), providing deep five-color photometry in a 2.5 degrees wide band along the equator. The stars in the tails make up a substantial part (~1/3) of the current total population of cluster stars in the magnitude interval 19.5 < i* < 22.0. This reveals that the cluster is subject to heavy mass loss. The orientation of the tails provides an important key for the determination of the cluster's Galactic orbit.

We discuss optical and infrared photometric properties of stars matched in the Two Micron All Sky Survey (2MASS) and the Sloan Digital Sky Survey (SDSS) commissioning data for 50 deg2 of sky. About 98% (63,000) of objects listed in the 2MASS Point Source Catalog in the analyzed area are matched within 2 arcsec to an SDSS source. The matched sources represent 8% of the about 800,000 SDSS sources in this area, and 15% of them are resolved in SDSS imaging data, although they are detected as point sources in 2MASS data. For about 14,000 stars with the smallest photometric errors (< 10%) in both surveys, we present optical and infrared color-magnitude and color-color diagrams. We use optical (SDSS) colors to identify the stellar spectral sequence and show that stars of different spectral types can have similar infrared colors, thus making the classification of stars based on only 2MASS data very difficult. However, a broad separation into ``early'' and ``late'' spectral types (relative to type K0) is possible with a reliability of about 95% even with 2MASS colors alone. The distributions of matched sources in color-magnitude and color-color diagrams are compared to the predictions of a stellar population synthesis code. We find that the models are in fair overall agreement with the data. The most significant discrepancies are found for the number ratio of ``early'' to ``late'' type stars (by about a factor of 2) and in the colors of M stars (up to 0.2 mag).

We report the discovery of a pair of z = 4.25 quasars with a separation of 33 arcseconds. The brighter of the two objects was identified as a high-redshift quasar candidate from Sloan Digital Sky Survey multicolor imaging data, and the redshift was measured from a spectrum obtained with the Hobby-Eberly Telescope. The slit orientation of this observation {\it by chance} included another quasar, approximately one magnitude fainter and having the same redshift as the target. This is the third serendipitous discovery of a z > 4 quasar. The differences in the relative strengths and profiles of the emission lines suggest that this is a quasar pair and not a gravitational lens. The two objects are likely to be physically associated; the projected physical separation is approximately 210 $h_{50}^{-1}$ kpc and the redshifts are identical to $\approx$ 0.01, implying a radial physical separation of 950 $h_{50}^{-1}$ kpc or less. The existence of this pair is strong circumstantial evidence that $z \sim 4$ quasars are clustered.

We report the discovery of five quasars with redshifts of 4.67 - 5.27 and z'-band magnitudes of 19.5-20.7 M_B ~ -27. All were originally selected as distant quasar candidates in optical/near-infrared photometry from the Sloan Digital Sky Survey (SDSS), and most were confirmed as probable high-redshift quasars by supplementing the SDSS data with J and K measurements. The quasars possess strong, broad Lyman-alpha emission lines, with the characteristic sharp cutoff on the blue side produced by Lyman-alpha forest absorption. Three quasars contain strong, broad absorption features, and one of them exhibits very strong N V emission. The amount of absorption produced by the Lyman-alpha forest increases toward higher redshift, and that in the z=5.27 object (D_A ~ 0.7) is consistent with a smooth extrapolation of the absorption seen in lower redshift quasars. The high luminosity of these objects relative to most other known objects at z >~ 5 makes them potentially valuable as probes of early quasar properties and of the intervening intergalactic medium.

We present a sample of 148 candidate RR Lyrae stars selected from Sloan Digital Sky Survey (SDSS) commissioning data for about 100 deg2 of sky surveyed twice with Delta(t) = 1.9946 days. Although the faint magnitude limit of the SDSS allows us to detect RR Lyrae stars to large galactocentric distances (about 100 kpc, or r = 21), we find no candidates fainter than r = 20, i.e. further than about 65 kpc from the Galactic center. On the assumption that all 148 candidates are indeed RR Lyrae stars (contamination by other species of variable star is probably less than 10%), we find that their volume density has roughly a power-law dependence on galactocentric radius, 1/R^(2.7+-0.2), between 10 and 50 kpc, and drops abruptly at R = 50-60 kpc, possibly indicating a sharp edge to the stellar halo as traced by RR Lyrae stars. The Galactic distribution of stars in this sample is very inhomogeneous and shows a clump of over 70 stars at about 45 kpc from the Galactic center. This clump is also detected in the distribution of nonvariable objects with RR Lyrae star colors. When sources in the clump are excluded, the best power-law fit becomes consistent with the 1/R^3 distribution found from surveys of bright RR Lyrae stars. These results imply that the halo contains clumpy overdensities inhomogeneously distributed within a smooth 1/R^3 background, with a possible cutoff at about 50 kpc.

A sample of 4208 objects with magnitude 15 < g* < 22 and colors of main sequence A stars has been selected from 370 square degrees of Sloan Digital Sky Survey (SDSS) commissioning observations. The data is from two long, narrow stripes, each with an opening angle of greater than 60 deg, at Galactic latitudes 36 < abs(b) < 63 on the celestial equator. An examination of the sample's distribution shows that these stars trace considerable substructure in the halo. Large overdensities of A-colored stars in the North at (l,b,R) = (350, 50, 46 kpc) and in the South at (157, -58, 33 kpc) and extending over tens of degrees are present in the halo of the Milky Way. Using photometry to separate the stars by surface gravity, both structures are shown to contain a sequence of low surface gravity stars consistent with identification as a blue horizontal branch (BHB). Both structures also contain a population of high surface gravity stars two magnitudes fainter than the BHB stars, consistent with their identification as blue stragglers (BSs). From the numbers of detected BHB stars, lower limits to the implied mass of the structures are 6x10^6 M_sun and 2x10^6 M_sun. The fact that two such large clumps have been detected in a survey of only 1% of the sky indicates that such structures are not uncommon in the halo. Simple spheroidal parameters are fit to a complete sample of the remaining unclumped BHB stars and yield (at r < 40 kpc) a fit to a halo distribution with flattening (c/a = 0.65+/-0.2) and a density falloff exponent of alpha = -3.2+/-0.3.

Unbiased analysis shows that IRAS data reliably differentiate between the early and late stages of stellar evolution because objects at these stages clearly segregate in infrared color-color diagrams. Structure in these diagrams is primarily controlled by the density distribution of circumstellar dust. The density profile around older objects is the steepest, declining as $r^{-2}$, while young objects have profiles that vary as $r^{-3/2}$ and flatter. The different density profiles reflect the different dynamics that govern the different environments. Our analysis also shows that high mass star formation is strongly concentrated within \about 5 kpc around the Galactic center, in support of other studies.

Atomic carbon has been detected in the envelopes of three carbon-rich evolved stars: HD 44179 (=AFGL 915, the `Red Rectangle'); HD 56126; and, tentatively, the carbon star V Hya. This brings to seven the number of evolved star envelopes in which CI has been detected. Upper limits were found for several other stars, including R CrB. CI was not detected in several oxygen rich post asymptotic giant branch (AGB) stars (OH231.8+4.2, for example), although it is detected in their carbon-rich analogues. Two trends are evident in the data. First, circumstellar envelopes with detectable CI are overwhelmingly carbon rich, suggesting that much of the CI is produced by the dissociation of molecules other than CO. Second, the more evolved the envelope away from the AGB, the higher is the CI/CO ratio. The oxygen-rich supergiant star Alpha Ori remains the only oxygen rich star with a wind containing detectable CI. These data suggest an evolutionary sequence for the CI/CO ratio in cool circumstellar envelopes. This ratio is small (a few %) while the star is on the AGB, and the CI is located in the outer envelope and produced by photodissociation. The ratio increases to about 0.5 as the star evolves away from the AGB because of the dissociation of CO and other carbon-bearing molecules by shocks caused by the fast winds which appear at the end of evolution on the AGB. Finally, the ratio becomes >> 1 as the central star becomes hot enough to photodissociate CO.

DUSTY solves the problem of radiation transport in a dusty environment. The code can handle both spherical and planar geometries. The user specifies the properties of the radiation source and dusty region, and the code calculates the dust temperature distribution and the radiation field in it. The solution method is based on a self-consistent equation for the radiative energy density, including dust scattering, absorption and emission, and does not introduce any approximations. The solution is exact to within the specified numerical accuracy. DUSTY has built in optical properties for the most common types of astronomical dust and comes with a library for many other grains. It supports various analytical forms for the density distribution, and can perform a full dynamical calculation for radiatively driven winds around AGB stars. The spectral energy distribution of the source can be specified analytically as either Planckian or broken power-law. In addition, arbitrary dust optical properties, density distributions and external radiation can be entered in user supplied files. Furthermore, the wavelength grid can be modified to accommodate spectral features. A single DUSTY run can process an unlimited number of models, with each input set producing a run of optical depths, as specified. The user controls the detail level of the output, which can include both spectral and imaging properties as well as other quantities of interest.

V Hya, an evolved carbon star with a complex circumstellar envelope, has two variability periods, 530d and 6000d (17 years). We analyze recent light curve data and show that both variations have benn present for at least 100 years and have been regular over this time. The 530d period and its 1.5 - 2 magnitude amplitude show that V Hya is a Mira variable. We suggest that the star is in a binary system (as also suspected from the structure of the circumstellar envelope) and that the 17-year variation is due to extinction by circumstellar dust orbiting with the companion. The properties of the envelope found from molecular line observations: the fast molecular wind, the relatively small size of the dense circumstellar envelope, and the high mass loss rate, all suggest that V Hya has entered its `superwind' phase. However, its spectral type, period, colors, and lack of ionizing radiation show that the star is still on the AGB. These properties add to the evidence that the complex structures of many planetary planetary nebulae, including fast stellar winds, originate during the final phases of mass loss on the AGB.

IR and mm-wave emission from Herbig Ae/Be stars has produced conflicting conclusions regarding the dust geometry in these objects. We show that the compact dimensions of the mm-wave emitting regions are a decisive indication for disks. But a disk cannot explain the spectral energy distribution (SED) unless it is embedded in an extended envelope that (1) dominates the IR emission and (2) provides additional disk heating on top of the direct stellar radiation. Detailed radiative transfer calculations based on the simplest model for envelope-embedded disks successfully fit the data from UV to mm wavelengths and show that the disks have central holes. This model also resolves naturally some puzzling results of IR imaging.

We find that AGB stars separate in the 25-12 vs. 12-K color-color diagram according to their chemistry (O, S vs. C) and variability type (Miras vs. SRb/Lb). While discrimination according to the chemical composition is not surprising, the separation of Miras from SRb/Lb variables is unexpected. We show that ``standard'' steady-state radiatively driven models provide excellent fits to the color distribution of Miras of all chemical types. However, these models are incapable of explaining the dust emission from O-rich SRb/Lb stars. The models can be altered to fit the data by postulating different optical properties for silicate grains, or by assuming that the dust temperature at the inner envelope radius is significantly lower (300-400 K) than typical condensation temperatures (800-1000 K), a possibility which is also supported by the detailed characteristics of LRS data. The absence of hot dust for SRb/Lb stars can be interpreted as a recent (order of 100 yr) decrease in the mass-loss rate. The distribution of O-rich SRb/Lb stars in the 25-12 vs. K-12 color-color diagram shows that the mass-loss rate probably resumes again, on similar time scales. It cannot be ruled out that the mass-loss rate is changing periodically on such time scales, implying that the stars might oscillate between the Mira and SRb/Lb phases during their AGB evolution as proposed by Kerschbaum et al. (1996). Such a possibility appears to be supported by numerous recent observations.

The Sloan Digital Sky Survey has obtained images in five broad-band colors for several hundred square degrees. We present color-color diagrams for stellar objects, and demonstrate that quasars are easily distinguished from stars by their distinctive colors. Follow-up spectroscopy in less than ten nights of telescope time has yielded 22 new quasars, 9 of them at $z> 3.65$, and one with $z = 4.75$, the second highest-redshift quasar yet known. Roughly 80% of the high-redshift quasar candidates selected by color indeed turn out to be high-redshift quasars.

We present a detailed, self-consistent model of radiatively driven stellar outflows which couples the radiative transfer and hydrodynamics equations. The circumstellar envelope, which consists of gas and dust, is described as a two-component fluid to account for relative drifts. Our results agree with both molecular line observations and infrared continuum spectra, and show that steady-state outflows driven by radiation pressure on dust grains adequately describe the surroundings of the majority of cool luminous evolved stars.

We have searched the IRAS Point Source Catalog and HIRES maps for young stellar objects (YSOs) in the direction of five \HI high-velocity clouds (HVCs). In agreement with optical searches in the halo, no evidence was found for extensive star-forming activity inside the high-latitude HVCs. Specifically, we have found no signs of star formation or YSOs in the direction of the A IV cloud or in the very-high-velocity clouds HVC~110-7-465 and HVC~114-10-440. We have identified only one young star in the direction of the M~I.1 cloud, which shows almost perfect alignment with a knot of \HI emission. Because of the small number of early-type stars observed in the halo, the probability for such a positional coincidence is low; thus, this young star appears to be physically associated with the M~I.1 cloud. We have also identified a good YSO candidate in the \HI shell-like structure observed in the core region of the low-latitude cloud complex H (HVC~131+1-200). This region could be a supernova remnant with several other YSO candidates formed along the shock front produced by the explosion. In agreement with recent theoretical estimates, these results point to a low but significant star-formation rate in intermediate and high Galactic latitude HVCs. For M~I.1 in particular, we estimate that the efficiency of the star-formation process is $M(YSO)/M(\HI)\ga 10^{-4}-10^{-3}$ by mass. Such efficiency is sufficient to account for (a) the existence of the few young blue stars whose ages imply that they were born in the Galactic halo, and (b) the nonprimordial metallicities inferred for some HVCs if their metal content proves to be low.

When verifying a sophisticated numerical code, it is a usual practice to compare the results with reliable solutions obtained by other means. This work provides such solutions for the wavelength dependent dust radiative transfer problem. We define a set of benchmark problems in spherical geometry and solve them by three radiative transfer codes which implement different numerical schemes. Results for the dust temperature and emerging spectra agree to better than 0.1%, and can be used as benchmark solutions for the verification of the dust radiative transfer codes.

Dust infrared emission possesses scaling properties. Overall luminosity is never an input parameter of the radiative transfer problem, spectral shape is the only relevant property of the heating radiation when the inner boundary of the dusty region is controlled by dust sublimation. Similarly, the absolute scales of densities and distances are irrelevant; the geometry enters only through angles, relative thicknesses and aspect ratios, and the actual magnitudes of densities and distances enter only through one independent parameter, the overall optical depth. Dust properties enter only through dimensionless, normalized distributions that describe the spatial variation of density and the wavelength dependence of scattering and absorption efficiencies. Scaling enables a systematic approach to modeling and classification of IR spectra. We develop a new, fully scale-free method for solving radiative transfer, present exact numerical results, and derive approximate analytical solutions for spherical geometry, covering the entire range of parameter space relevant to observations. Scaling implies tight correlations among the SEDs of various members of the same class of sources such as young stellar objects, late-type stars, etc. In particular, all members of the same class occupy common, well defined regions in color-color diagrams. The observational data corroborate the existence of these correlations.

The spectral shape of IR emission from Herbig Ae/Be stars has been invoked as evidence for accretion disks around high-mass protostars. Instead, we present here models based on spherical envelopes with $r^{-1.5}$ dust density profile that successfully explain the observed spectral shapes. The spectral energy distributions (SEDs) of eight primary candidates for protostellar disks are fitted in detail for all wavelengths available, from visual to far IR. The only envelope property adjusted in individual sources is the overall visual optical depth, and it ranges from 0.3 to 3. In each case, our models properly reproduce the data for both IR excess, visual extinction and reddening. The success of our models shows that accretion disks cannot make a significant contribution to the radiation observed in these pre-main sequence stars.

Infrared emission from radiatively heated dust possesses general scaling properties. The spectral shape is independent of overall luminosity when the inner boundary of the dusty region is controlled by dust sublimation; the only relevant property of the heating radiation is its frequency profile. For a given type of dust grains, the emission from red giants and supergiants is essentially controlled by a single parameter --- the overall dust optical depth. This leads to tight correlations among different spectral properties that explain many available observations and enable systematic studies of large data bases.

Infrared emission from the dust shell around IRC+10216 is analysed in detail, employing a self-consistent model for radiatively driven winds around late-type stars that couples the equations of motion and radiative transfer in the dust. The resulting model provides agreement with the wealth of available data, including the spectral energy distribution in the range 0.5--1000 \mic, and visibility and array observations. Previous conclusions about two dust shells, derived from modelling the data with a few single-temperature components of different radii, are not supported by our results. The extended, continuous temperature and density distributions derived from our model obviate the need for such discrete shells. The IR properties vary with the stellar phase, reflecting changes in both the dust condensation radius $r_1$ and the overall optical depth $\tau$ -- as the luminosity increases from minimum to maximum, $r_1$ increases while $\tau$ decreases. We find that the angular size of the dust condensation zone varies from 0.3 arcsec at minimum light to 0.5 arcsec at maximum. The shortage of flux at short wavelengths encountered in previous studies is resolved by employing a grain size distribution that includes grains larger than \about\ 0.1 \mic, required also for the visibility fits. This distribution is in agreement with the one recently proposed by Jura in a study that probed the outer regions of the envelope. Since our constraints on the size distribution mostly reflect the envelope's inner regions, the agreement of these independent studies is evidence against significant changes in grain sizes through effects like sputtering or grain growth after the initial formation at the dust condensation zone.

Infrared imaging properties of dusty winds around late-type stars are investigated in detail, employing a self-consistent model that couples the equations of motion and radiative transfer. Because of general scaling properties, the angular profiles of surface brightness are self-similar. In any given star, the profile shape is determined essentially by overall optical depth at each wavelength and it is self-similarly scaled by the size of the dust condensation zone. We find that the mid-IR is the best wavelength range to measure directly the angular size of this zone, and from {\it IRAS} data we identify the 15 best candidates for such future observations. We also show that the visibility function at short wavelengths ($\la$ 2 \mic) directly determines the scattering optical depth, and produce theoretical visibility curves for various characteristic wavelengths and the entire parameter range relevant to late-type stars. The infrared emission should display time variability because of cyclical changes in overall optical depth, reflecting luminosity-induced movement of the dust condensation point. Calculations of the wavelength dependence of photometric amplitudes and time variability of envelope sizes are in agreement with observations; envelopes are bigger and bluer at maximum light.

The dynamical structure and infrared emission of winds around late-type stars are studied in a self-consistent model that couples the equations of motion and radiative transfer. Both the dynamics and IR spectrum of the solution are fully characterized by the flux averaged optical depth. Five types of dust grains are considered: astronomical silicate, SiC, crystalline olivine, graphite and amorphous carbon as well as mixtures. Analysis of infrared signatures provides constraints on the grain chemical composition and indications for the simultaneous existence of silicate and carbon grains. The region of parameter space where radiation pressure can support a given mass-loss rate is identified and a new method to derive mass loss rates from IR data is presented. IRAS LRS classes are associated with optical depth for various grain materials and the regions of color-color diagrams expected to be populated by late-type stars are identified. We show that cirrus emission can severely affect point source measurements to the extent that their listed IRAS long wavelength fluxes are unreliable. After accounting for cirrus contamination, essentially all IRAS point sources located in the relevant regions of the color-color diagrams can be explained as late-type stars.