We present preliminary trigonometric parallaxes of 184 late-T and Y dwarfs using observations from Spitzer (143), USNO (18), NTT (14), and UKIRT (9). To complete the 20-pc census of $\ge$T6 dwarfs, we combine these measurements with previously published trigonometric parallaxes for an additional 44 objects and spectrophotometric distance estimates for another 7. For these 235 objects, we estimate temperatures, sift into five 150K-wide $T_{\rm eff}$ bins covering the range 300-1050K, determine the completeness limit for each, and compute space densities. To anchor the high-mass end of the brown dwarf mass spectrum, we compile a list of early- to mid-L dwarfs within 20 pc. We run simulations using various functional forms of the mass function passed through two different sets of evolutionary code to compute predicted distributions in $T_{\rm eff}$. The best fit of these predictions to our L, T, and Y observations is a simple power-law model with $\alpha \approx 0.6$ (where $dN/dM \propto M^{-\alpha}$), meaning that the slope of the field substellar mass function is in rough agreement with that found for brown dwarfs in nearby star forming regions and young clusters. Furthermore, we find that published versions of the log-normal form do not predict the steady rise seen in the space densities from 1050K to 350K. We also find that the low-mass cutoff to formation, if one exists, is lower than $\sim$5 $M_{Jup}$, which corroborates findings in young, nearby moving groups and implies that extremely low-mass objects have been forming over the lifetime of the Milky Way.

This paper presents new trigonometric parallaxes and proper motions for 214 stars. The measurements were made at the US Naval Observatory Flagstaff Station (NOFS) between 1989 and 2017, and the average uncertainty in the parallax values is 0.6 mas. We find good agreement with Gaia Data Release 2 measurements for the stars in common, although there may be a small systematic offset similar to what has been found by other investigators. The sample is matched to catalogs and the literature to create a photometric dataset which spans the ultraviolet to the mid-infrared. New mid-infrared photometry is obtained for nineteen stars from archived Spitzer mosaics. New optical spectroscopy is presented for seven systems and additional spectra were obtained from the literature. We identify a sub-sample of 179 white dwarfs (WDs) at distances of 25 - 200 pc. Their spectral energy distributions (SEDs) are analyzed using model atmospheres. The models reproduce the entire flux-calibrated SED very well and provide the atmospheric chemical composition, temperature, surface gravity, mass and cooling age of each WD. Twenty-six WDs are newly classified and twelve systems are presented as candidate unresolved binaries. We confirm one WD+red dwarf system and identify two WDs as candidate dust disk systems. Twelve old and high-velocity systems are identified as candidate thick disk or halo objects. The WDs in the sample generally have Galactic disk-like ages of < 8 Gyr and masses close to the canonical 0.6 M_Sun.

Parallaxes are presented for a sample of 20 nearby dwarf carbon stars. The inferred luminosities cover almost two orders of magnitude. Their absolute magnitudes and tangential velocities confirm prior expectations that some originate in the Galactic disk, although more than half of this sample are halo stars. Three stars are found to be astrometric binaries, and orbital elements are determined; their semimajor axes are 1 -- 3 AU, consistent with the size of an AGB mass-transfer donor star.

APOGEE-2 is a high-resolution, near-infrared spectroscopic survey observing roughly 300,000 stars across the entire sky. It is the successor to APOGEE and is part of the Sloan Digital Sky Survey IV (SDSS-IV). APOGEE-2 is expanding upon APOGEE's goals of addressing critical questions of stellar astrophysics, stellar populations, and Galactic chemodynamical evolution using (1) an enhanced set of target types and (2) a second spectrograph at Las Campanas Observatory in Chile. APOGEE-2 is targeting red giant branch (RGB) and red clump (RC) stars, RR Lyrae, low-mass dwarf stars, young stellar objects, and numerous other Milky Way and Local Group sources across the entire sky from both hemispheres. In this paper, we describe the APOGEE-2 observational design, target selection catalogs and algorithms, and the targeting-related documentation included in the SDSS data releases.

New, updated, and/or revised CCD parallaxes determined with the Strand Astrometric Reflector at the Naval Observatory Flagstaff Station (NOFS) are presented. Included are results for 309 late-type dwarf and subdwarf stars observed over the 30+ years that the program operated. For 124 of the stars, parallax determinations from other investigators have already appeared in the literature and we compare the different results. Also included here is new or updated $VI$ photometry on the Johnson-Kron-Cousins system for all but a few of the faintest targets. Together with 2MASS $JHK_s$ near-infrared photometry, a sample of absolute magnitude versus color and color versus color diagrams are constructed. Since large proper motion was a prime criterion for targeting the stars, the majority turn out to be either M-type subdwarfs or late M-type dwarfs. The sample also includes 50 dwarf or subdwarf L-type stars, and four T dwarfs. Possible halo subdwarfs are identified in the sample based on tangential velocity, subluminosity, and spectral type. Residuals from the solutions for parallax and proper motion for several stars show evidence of astrometric perturbations.

We present a detailed analysis of the white dwarf luminosity functions derived from the local 40 pc sample and the deep proper motion catalog of Munn et al (2014, 2017). Many of the previous studies ignored the contribution of thick disk white dwarfs to the Galactic disk luminosity function, which results in an erronous age measurement. We demonstrate that the ratio of thick/thin disk white dwarfs is roughly 20\% in the local sample. Simultaneously fitting for both disk components, we derive ages of 6.8-7.0 Gyr for the thin disk and 8.7 $\pm$ 0.1 Gyr for the thick disk from the local 40 pc sample. Similarly, we derive ages of 7.4-8.2 Gyr for the thin disk and 9.5-9.9 Gyr for the thick disk from the deep proper motion catalog, which shows no evidence of a deviation from a constant star formation rate in the past 2.5 Gyr. We constrain the time difference between the onset of star formation in the thin disk and the thick disk to be $1.6^{+0.3}_{-0.4}$ Gyr. The faint end of the luminosity function for the halo white dwarfs is less constrained, resulting in an age estimate of $12.5^{+1.4}_{-3.4}$ Gyr for the Galactic inner halo. This is the first time ages for all three major components of the Galaxy are obtained from a sample of field white dwarfs that is large enough to contain significant numbers of disk and halo objects. The resultant ages agree reasonably well with the age estimates for the oldest open and globular clusters.

The dwarf carbon (dC) star SDSS J112801.67+004034.6 has an unusually high radial velocity, 531$\pm 4$ km s$^{-1}$. We present proper motion and new spectroscopic observations which imply a large Galactic rest frame velocity, 425$\pm 9$ km s$^{-1}$. Several other SDSS dC stars are also inferred to have very high galactocentric velocities, again each based on both high heliocentric radial velocity and also confidently detected proper motions. Extreme velocities and the presence of $C_2$ bands in the spectra of dwarf stars are both rare. Passage near the Galactic center can accelerate stars to such extreme velocities, but the large orbital angular momentum of SDSS J1128 precludes this explanation. Ejection from a supernova in a binary system or disruption of a binary by other stars are possibilities, particularly as dC stars are thought to obtain their photospheric $C_2$ via mass transfer from an evolved companion.

A catalog of 8472 white dwarf (WD) candidates is presented, selected using reduced proper motions from the deep proper motion catalog of Munn et al. 2014. Candidates are selected in the magnitude range 16 < r < 21.5 over 980 square degrees, and 16 < r < 21.3 over an additional 1276 square degrees, within the Sloan Digital Sky Survey (SDSS) imaging footprint. Distances, bolometric luminosities, and atmospheric compositions are derived by fitting SDSS ugriz photometry to pure hydrogen and helium model atmospheres (assuming surface gravities log g = 8). The disk white dwarf luminosity function (WDLF) is constructed using a sample of 2839 stars with 5.5 < M_bol < 17, with statistically significant numbers of stars cooler than the turnover in the luminosity function. The WDLF for the halo is also constructed, using a sample of 135 halo WDs with 5 < M_bol < 16. We find space densities of disk and halo WDs in the solar neighborhood of 5.5 +- 0.1 x 10^-3 pc^-3 and 3.5 +- 0.7 x 10^-5 pc^-3, respectively. We resolve the bump in the disk WDLF due to the onset of fully convective envelopes in WDs, and see indications of it in the halo WDLF as well.

We use extensive gravity-sensitive DDO 51 photometry over 5100 square degrees, combined with SDSS broadband photometry, to select a catalog of $\sim 4,000$ giant stars covering a large fraction of the high Galactic latitude sky and reaching out to $\sim 80$ kpc in the Galactic halo. This sample of bright and unbiased tracers enables us to measure the radial profile and 3D structure of the stellar halo to large distance which had previously only been measured with sparse tracers or small samples. Using population synthesis models to reproduce the observed giant star luminosity function, we find that the halo maintains a $r^{-3.5}$ profile from $30$ to $80$ kpc with no signs of a truncation or sharp break over this range. The radial profile measurement is largely insensitive to individual halo substructure components, but we find that attempting to measure the shape of the halo is overwhelmed by the Sagittarius stream such that no ellipsoidal shape is a satisfactory description in this region. These measurements allow us to begin placing the Milky Way in context with the growing sample of external galaxies where similar halo profile measurements are available, with the goal of further linking the properties of stellar halos to the accretion histories that formed them.

We present optical spectroscopy and near-infrared photometry of 57 faint ($g= 19-22$) high proper motion white dwarfs identified through repeat imaging of $\approx3100$ square degrees of the Sloan Digital Sky Survey footprint by Munn et al. (2014). We use $ugriz$ and $JH$ photometry to perform a model atmosphere analysis, and identify ten ultracool white dwarfs with $T_{\rm eff}<4000$ K, including the coolest pure H atmosphere white dwarf currently known, J1657+2638, with $T_{\rm eff}= 3550 \pm 100$K. The majority of the objects with cooling ages larger than 9 Gyr display thick disc kinematics and constrain the age of the thick disc to $\geq11$ Gyr. There are four white dwarfs in our sample with large tangential velocities ($v_{tan} > 120$ km $s^{-1}$) and UVW velocities that are more consistent with the halo than the Galactic disc. For typical $0.6 M_{\odot}$ white dwarfs, the cooling ages for these halo candidates range from 2.3 to 8.5 Gyr. However, the total main-sequence + white dwarf cooling ages of these stars would be consistent with the Galactic halo if they are slightly undermassive. Given the magnitude limits of the current large scale surveys, many of the coolest and oldest white dwarfs remain undiscovered in the solar neighborhood, but upcoming surveys such as GAIA and the Large Synoptic Survey Telescope (LSST) should find many of these elusive thick disc and halo white dwarfs.

We present the first results from a minute cadence survey of a three square degree field obtained with the Dark Energy Camera. We imaged part of the Canada-France-Hawaii Telescope Legacy Survey area over eight half-nights. We use the stacked images to identify 111 high proper motion white dwarf candidates with g $\leq24.5$ mag and search for eclipse-like events and other sources of variability. We find a new g = 20.64 mag pulsating ZZ Ceti star with pulsation periods of 11-13 min. However, we do not find any transiting planetary companions in the habitable zone of our target white dwarfs. Given the probability of eclipses of 1% and our observing window from the ground, the non-detection of such companions in this first field is not surprising. Minute cadence DECam observations of additional fields will provide stringent constraints on the frequency of planets in the white dwarf habitable zone.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE), one of the programs in the Sloan Digital Sky Survey III (SDSS-III), has now completed its systematic, homogeneous spectroscopic survey sampling all major populations of the Milky Way. After a three year observing campaign on the Sloan 2.5-m Telescope, APOGEE has collected a half million high resolution (R~22,500), high S/N (>100), infrared (1.51-1.70 microns) spectra for 146,000 stars, with time series information via repeat visits to most of these stars. This paper describes the motivations for the survey and its overall design---hardware, field placement, target selection, operations---and gives an overview of these aspects as well as the data reduction, analysis and products. An index is also given to the complement of technical papers that describe various critical survey components in detail. Finally, we discuss the achieved survey performance and illustrate the variety of potential uses of the data products by way of a number of science demonstrations, which span from time series analysis of stellar spectral variations and radial velocity variations from stellar companions, to spatial maps of kinematics, metallicity and abundance patterns across the Galaxy and as a function of age, to new views of the interstellar medium, the chemistry of star clusters, and the discovery of rare stellar species. As part of SDSS-III Data Release 12, all of the APOGEE data products are now publicly available.

URAT1 is an observational, astrometric catalog covering most of the Dec >= -15 deg area and a magnitude range of about R = 3 to 18.5. Accurate positions (typically 10 to 30 mas standard error) are given for over 228 million objects at a mean epoch around 2013.5. For the over 188 million objects matched with the 2MASS point source catalog proper motions (typically 5 to 7 mas/yr std. errors) are provided. These data are supplemented by 2MASS and APASS photometry. Observations, reductions and catalog construction are described together with results from external data verifications. The catalog data are served by CDS, Starsbourg (I/329). There is no DVD release.

A new proper motion catalog is presented, combining the Sloan Digital Sky Survey (SDSS) with second epoch observations in the r band within a portion of the SDSS imaging footprint. The new observations were obtained with the 90prime camera on the Steward Observatory Bok 90 inch telescope, and the Array Camera on the U.S. Naval Observatory, Flagstaff Station, 1.3 meter telescope. The catalog covers 1098 square degrees to r = 22.0, an additional 1521 square degrees to r = 20.9, plus a further 488 square degrees of lesser quality data. Statistical errors in the proper motions range from 5 mas/year at the bright end to 15 mas/year at the faint end, for a typical epoch difference of 6 years. Systematic errors are estimated to be roughly 1 mas/year for the Array Camera data, and as much as 2 - 4 mas/year for the 90prime data (though typically less). The catalog also includes a second epoch of r band photometry.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

A large fraction of white dwarfs (WDs) may host planets in their habitable zones. These planets may provide our best chance to detect bio-markers on a transiting exoplanet, thanks to the diminished contrast ratio between the Earth-sized WD and its Earth-sized planets. The JWST is capable of obtaining the first spectroscopic measurements of such planets, yet there are no known planets around WDs. Here we propose to take advantage of the unique capability of the Kepler spacecraft in the 2-Wheels mode to perform a transit survey that is capable of identifying the first planets in the habitable zone of a WD. We propose to obtain Kepler time-series photometry of 10,000 WDs in the SDSS imaging area to search for planets in the habitable zone. Thanks to the large field of view of Kepler, for the first time in history, a large number of WDs can be observed at the same time, which is essential for discovering transits. Our proposed survey requires a total of 200 days of observing time, and will find up to 100 planets in the WD habitable zone. This survey will maintain Kepler's spirit of searching for habitable Earths, but near new hosts. With few-day observations and minute-cadences per field, it will also open up a completely unexplored discovery space. In addition to planets, this survey is sensitive to pulsating WDs, as well as eclipsing short period stellar and substellar companions. These have important implications for constraining the double WD merger rate and their contribution to Type Ia supernovae and the gravitational wave foreground. Given the relatively low number density of our targets, this program can be combined with other projects that would benefit from high cadence and many-fields observations with Kepler, e.g. a transit survey of a magnitude-limited, complete sample of nearby M dwarfs or asteroseismology of variable stars (e.g. RR Lyrae) in the same fields.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a high-resolution infrared spectroscopic survey spanning all Galactic environments (i.e., bulge, disk, and halo), with the principal goal of constraining dynamical and chemical evolution models of the Milky Way. APOGEE takes advantage of the reduced effects of extinction at infrared wavelengths to observe the inner Galaxy and bulge at an unprecedented level of detail. The survey's broad spatial and wavelength coverage enables users of APOGEE data to address numerous Galactic structure and stellar populations issues. In this paper we describe the APOGEE targeting scheme and document its various target classes to provide the necessary background and reference information to analyze samples of APOGEE data with awareness of the imposed selection criteria and resulting sample properties. APOGEE's primary sample consists of ~100,000 red giant stars, selected to minimize observational biases in age and metallicity. We present the methodology and considerations that drive the selection of this sample and evaluate the accuracy, efficiency, and caveats of the selection and sampling algorithms. We also describe additional target classes that contribute to the APOGEE sample, including numerous ancillary science programs, and we outline the targeting data that will be included in the public data releases.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

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 Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

We present optical spectroscopy and near-infrared photometry of 126 cool white dwarfs in the Sloan Digital Sky Survey (SDSS). Our sample includes high proper motion targets selected using the SDSS and USNO-B astrometry and a dozen previously known ultracool white dwarf candidates. Our optical spectroscopic observations demonstrate that a clean selection of large samples of cool white dwarfs in the SDSS (and the SkyMapper, Pan-STARRS, and the Large Synoptic Survey Telescope datasets) is possible using a reduced proper motion diagram and a tangential velocity cut-off (depending on the proper motion accuracy) of 30 km/s. Our near-infrared observations reveal eight new stars with significant absorption. We use the optical and near-infrared photometry to perform a detailed model atmosphere analysis. More than 80% of the stars in our sample are consistent with either pure hydrogen or pure helium atmospheres. However, the eight stars with significant infrared absorption and the majority of the previously known ultracool white dwarf candidates are best explained with mixed hydrogen and helium atmosphere models. The age distribution of our sample is consistent with a Galactic disk age of 8 Gyr. A few ultracool white dwarfs may be as old as 12-13 Gyr, but our models have problems matching the spectral energy distributions of these objects. There are only two halo white dwarf candidates in our sample. However, trigonometric parallax observations are required for accurate mass and age determinations and to confirm their membership in the halo.

We report the discovery of three nearby old halo white dwarf candidates in the Sloan Digital Sky Survey (SDSS), including two stars in a common proper motion binary system. These candidates are selected from our 2800 square degree proper motion survey on the Bok and U.S. Naval Observatory Flagstaff Station 1.3m telescopes, and they display proper motions of 0.4-0.5 arcsec/yr. Follow-up MMT spectroscopy and near-infrared photometry demonstrate that all three objects are hydrogen-dominated atmosphere white dwarfs with Teff = 3700 - 4100 K. For average mass white dwarfs, these temperature estimates correspond to cooling ages of 9-10 Gyr, distances of 70-80 pc, and tangential velocities of 140-200 km/s. Based on the UVW space velocities, we conclude that they most likely belong to the halo. Furthermore, the combined main-sequence and white dwarf cooling ages are 10-11 Gyr. Along with SDSS J1102+4113, they are the oldest field white dwarfs currently known. These three stars represent only a small fraction of the halo white dwarf candidates in our proper motion survey, and they demonstrate that deep imaging surveys like the Pan-STARRS and Large Synoptic Survey Telescope should find many old thick disk and halo white dwarfs that can be used to constrain the age of the Galactic thick disk and halo.

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.

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]

We identify seven new ultracool white dwarfs discovered in the Sloan Digital Sky Survey (SDSS). The SDSS photometry, spectra, and proper motions are presented, and additional BVRI data are given for these and other previously discovered ultracool white dwarfs. The observed colors span a remarkably wide range, qualitatively similar to colors predicted by models for very cool white dwarfs. One of the new stars (SDSS J1251+44) exhibits strong collision-induced absorption (CIA) in its spectra, while the spectra and colors of the other six are consistent with mild CIA. Another of the new discoveries (SDSS J2239+00A) is part of a binary system -- its companion is also a cool white dwarf, and other data indicate that the companion exhibits an infrared flux deficiency, making this the first binary system composed of two CIA white dwarfs. A third discovery (SDSS J0310-00) has weak Balmer emission lines. The proper motions of all seven stars are consistent with membership in the disk or thick disk.

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.

The halo of the Milky Way provides unique elemental abundance and kinematic information on the first objects to form in the Universe, which can be used to tightly constrain models of galaxy formation and evolution. Although the halo was once considered a single component, evidence for its dichotomy has slowly emerged in recent years from inspection of small samples of halo objects. Here we show that the halo is indeed clearly divisible into two broadly overlapping structural components -- an inner and an outer halo -- that exhibit different spatial density profiles, stellar orbits and stellar metallicities (abundances of elements heavier than helium). The inner halo has a modest net prograde rotation, whereas the outer halo exhibits a net retrograde rotation and a peak metallicity one-third that of the inner halo. These properties indicate that the individual halo components probably formed in fundamentally different ways, through successive dissipational (inner) and dissipationless (outer) mergers and tidal disruption of proto-Galactic clumps.

A 2.5 x 100 degree region along the celestial equator (Stripe 82) has been imaged repeatedly from 1998 to 2005 by the Sloan Digital Sky Survey. A new catalogue of ~4 million light-motion curves, together with over 200 derived statistical quantities, for objects in Stripe 82 brighter than r~21.5 has been constructed by combining these data by Bramich et al. (2007). This catalogue is at present the deepest catalogue of its kind. Extracting the ~130000 objects with highest signal-to-noise ratio proper motions, we build a reduced proper motion diagram to illustrate the scientific promise of the catalogue. In this diagram disk and halo subdwarfs are well-separated from the cool white dwarf sequence. Our sample of 1049 cool white dwarf candidates includes at least 8 and possibly 21 new ultracool H-rich white dwarfs (T_eff < 4000K) and one new ultracool He-rich white dwarf candidate identified from their SDSS optical and UKIDSS infrared photometry. At least 10 new halo white dwarfs are also identified from their kinematics.

We present the first parallax and luminosity measurements for an L subdwarf, the sdL7 2MASS J05325346+8246465. Observations conducted over three years by the USNO infrared astrometry program yield an astrometric distance of 26.7+/-1.2 pc and a proper motion of 2.6241+/-0.0018"/yr. Combined with broadband spectral and photometric measurements, we determine a luminosity of log(Lbol/Lsun) = -4.24+/-0.06 and Teff = 1730+/-90 K (the latter assuming an age of 5-10 Gyr), comparable to mid-type L field dwarfs. Comparison of the luminosity of 2MASS J05325346+8246465 to theoretical evolutionary models indicates that its mass is just below the sustained hydrogen burning limit, and is therefore a brown dwarf. Its kinematics indicate a ~110 Myr, retrograde Galactic orbit which is both eccentric (3 <~ R <~ 8.5 kpc) and extends well away from the plane (Delta_Z = +/-2 kpc), consistent with membership in the inner halo population. The relatively bright J-band magnitude of 2MASS J05325346+8246465 implies significantly reduced opacity in the 1.2 micron region, consistent with inhibited condensate formation as previously proposed. Its as yet unknown subsolar metallicity remains the primary limitation in constraining its mass; determination of both parameters would provide a powerful test of interior and evolutionary models for low-mass stars and brown dwarfs.

We have assembled a spectroscopic sample of low-mass dwarfs observed as part of the Sloan Digital Sky Survey along one Galactic sightline, designed to investigate the observable properties of the thin and thick disks. This sample of ~7400 K and M stars also has measured ugriz photometry, proper motions, and radial velocities. We have computed UVW space motion distributions, and investigate their structure with respect to vertical distance from the Galactic Plane. We place constraints on the velocity dispersions of the thin and thick disks, using two-component Gaussian fits. We also compare these kinematic distributions to a leading Galactic model. Finally, we investigate other possible observable differences between the thin and thick disks, such as color, active fraction and metallicity.

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.

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.

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.

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

A sample of white dwarfs is selected from SDSS DR3 imaging data using their reduced proper motions, based on improved proper motions from SDSS plus USNO-B combined data. Numerous SDSS and followup spectra (Kilic et al. 2005) are used to quantify completeness and contamination of the sample; kinematic models are used to understand and correct for velocity-dependent selection biases. A luminosity function is constructed covering the range 7 < M_bol < 16, and its sensitivity to various assumptions and selection limits is discussed. The white dwarf luminosity function based on 6000 stars is remarkably smooth, and rises nearly monotonically to M_bol = 15.3. It then drops abruptly, although the small number of low-luminosity stars in the sample and their unknown atmospheric composition prevent quantitative conclusions about this decline. Stars are identified that may have high tangential velocities, and a preliminary luminosity function is constructed for them.

We present the third edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog. The catalog consists of the 46,420 objects in the SDSS Third Data Release that have luminosities larger than M_i = -22 (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 are unambiguously broad absorption line quasars, are fainter than i = 15.0, and have highly reliable redshifts. The area covered by the catalog is 4188 sq. deg. The quasar redshifts range from 0.08 to 5.41, with a median value of 1.47; the high-redshift sample includes 520 quasars at redshifts greater than four, of which 17 are at redshifts greater than five. 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 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 about 2000; the spectra can be retrieved from the public database using the information provided in the catalog. A total of 44,221 objects in the catalog were discovered by the SDSS; 28,400 of the SDSS discoveries are reported here for the first time.

A reduced proper motion diagram utilizing Sloan Digital Sky Survey (SDSS) photometry and astrometry and USNO-B plate astrometry is used to separate cool white dwarf candidates from metal-weak, high-velocity main sequence Population II stars (subdwarfs) in the SDSS Data Release 2 imaging area. Follow-up spectroscopy using the Hobby-Eberly Telescope, the MMT, and the McDonald 2.7m Telescope is used to demonstrate that the white dwarf and subdwarf loci separate cleanly in the reduced proper motion diagram, and that the contamination by subdwarfs is small near the cool white dwarf locus. This enables large statistically complete samples of white dwarfs, particularly the poorly understood cool white dwarfs, to be created from the SDSS imaging survey, with important implications for white dwarf luminosity function studies. SDSS photometry for our sample of cool white dwarfs is compared to current white dwarf models.

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 describe a new survey for cool white dwarfs that supplements Sloan Digital Sky Survey photometry with USNO proper motions and follow-up spectroscopy. To date we have discovered and spectroscopically confirmed 80 new moderate temperature and cool white dwarfs. We have also found a handful of high-velocity white dwarfs and we expect a sizable fraction of these to be thick disk or possibly halo objects. Our survey is designed to find ~10,000 new white dwarfs, although only ~60 will be among the faintest white dwarfs (M(V) >= 16), where most of the age-sensitivity resides. We discuss an extension of our survey to V ~ 22.

We have compiled an extensive dataset on potential parts of the Monoceros tidal stream and performed an exhaustive survey of dwarf galaxy semi-analytic orbits in order to constrain its orbital properties. The best-fit orbits are subsequently realized as self-consistent N-body simulations in order to reproduce the spatial and velocity distribution of satellite debris. We find that all kinematic and geometric constraints can be fit by a single stream allowing for multiple wraps. The orbital eccentricity and inclination of the progenitor are strongly constrained to be e = 0.10 \pm 0.05 and i = 25 \pm 5 deg. Ten new estimates of proper motions from the Sloan Digital Sky Survey (SDSS) clearly exclude all retrograde orbits. Particles lost by the satellite populate two, nearly concentric rings naturally explaining the detection of stream stars at both 6--8 kpc (Ibata et al. 2003, Newberg et al. 2003) and 12--18 kpc (the Tri/And stream; Rocha-Pinto et al. 2003) from the Sun. We have attempted to predict the present location of the Monoceros stream progenitor using different information: (i) the kinematical and spatial distribution of detections and (ii) the different mean metallicity in the innner and the outer rings. A current location at l ~ 245, b ~ -18., with a distance to the Sun r_s~15 kpc, reproduces the available constraints. The mass of the progenitor has been poorly determined due to the slow orbital decay. Similar fits have been obtained for masses (3-9)10^8 M_sun. Additionally, we have analyzed to possible common origin of the Canis Major dwarf and the Monoceros stream, unfortunately, without a conclusive result.

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 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 have compiled L' and M' photometry of 63 single and binary M, L,and T dwarfs obtained at UKIRT using the MKO filter set. This compilation includes new L' data for 8 L dwarfs and 13 T dwarfs and new M' data for 7 L dwarfs, 5 T dwarfs, and Gl 229A. We compute L_bol, BC_K, and T_eff for 42 dwarfs whose flux-calibrated JHK spectra, L' photometry, and trigonometric parallaxes are available, and we estimate these quantities for 9 other dwarfs whose parallaxes and flux-calibrated spectra have been obtained. BC_K is a well-behaved function of near-IR spectral type with a dispersion of ~0.1 mag for types M6-T5; it is significantly more scattered for types T5-T9. T_eff declines steeply and monotonically for types M6-L7 and T4-T9, but is nearly constant at ~1450 K for types L7-T4 with assumed ages of ~3 Gyr. This constant T_eff is evidenced by nearly unchanging values of L'-M' between types L6 and T3. We compare the observed K, L', and M' luminosities of L and T dwarfs in our sample with those predicted by model atmospheres with varying surface gravities and sedimentation efficiencies. The models indicate that the L3-T4.5 dwarfs generally have higher gravities (log g = 5.0-5.5) than the T6-T9 dwarfs (log g = 4.5-5.0). The predicted M' luminosities of T dwarfs are 1.5-2.5 times larger than those derived empirically for the T dwarfs in our sample. This discrepancy is attributed to absorption at 4.5-4.9 um by CO, which is not expected under the condition of chemical equilibrium assumed in the models. We determine that the L3 dwarf Kelu-1 and the T0 dwarf SDSS J0423-0414 are probable binary systems. We compute log(L_bol/L_sun) = -5.73 +/- 0.05 and T_eff = 600-750 K for the T9 dwarf 2MASS J0415-0935, which supplants Gl 570D as the least luminous and coolest brown dwarf presently known.

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 present the full spectroscopic white dwarf and hot subdwarf sample from the SDSS first data release, DR1. We find 2551 white dwarf stars of various types, 240 hot subdwarf stars, and an additional 144 objects we have identified as uncertain white dwarf stars. Of the white dwarf stars, 1888 are non-magnetic DA types and 171, non-magnetic DBs. The remaining (492) objects consist of all different types of white dwarf stars: DO, DQ, DC, DH, DZ, hybrid stars like DAB, etc., and those with non-degenerate companions. We fit the DA and DB spectra with a grid of models to determine the Teff and log(g) for each object. For all objects, we provide coordinates, proper motions, SDSS photometric magnitudes, and enough information to retrieve the spectrum/image from the SDSS public database. This catalog nearly doubles the known sample of spectroscopically-identified white dwarf stars. In the DR1 imaged area of the sky, we increase the known sample of white dwarf stars by a factor of 8.5. We also comment on several particularly interesting objects in this sample.

We present preliminary trigonometric parallaxes and proper motions for 22 L dwarfs and 18 T dwarfs measured using the ASTROCAM infrared imager. Relative to absolute parallax corrections are made by employing 2MASS and/or SDSS photometry for reference frame stars. We combine USNO infrared and optical parallaxes with the best available CIT system photometry to determine M_J, M_H, and M_K values for 37 L dwarfs between spectral types L0 to L8 and 19 T dwarfs between spectral types T0.5 and T8 and present selected absolute magnitude versus spectral type and color diagrams, based on these results. Luminosities and temperatures are estimated for these objects. Of special interest are the distances of several objects which are at or near the L-T dwarf boundary so that this important transition can be better understood. The previously reported early-mid T dwarf luminosity excess is clearly confirmed and found to be present at J, H, and K. The large number of objects that populate this luminosity excess region indicates that it cannot be due entirely to selection effects. The T dwarf sequence is extended to M_J~16.9 by 2MASS J041519-0935 which, at d = 5.74 pc, is found to be the least luminous [log(L/L_sun)=-5.58] and coldest (T_eff~760 K) brown dwarf known. Combining results from this paper with earlier USNO CCD results we find that, in contrast to the L dwarfs, there are no examples of low velocity (V_tan < 20 km/s) T dwarfs. We briefly discuss future directions for the USNO infrared astrometry program.

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.

Many open questions in X-ray astronomy are limited by the relatively small number of objects in uniform optically-identified samples, especially when rare subclasses are considered, or subsets isolated to search for evolution or correlations between wavebands. We describe initial results of a program aimed to ultimately yield 10^4 X-ray source identifications--a sample about an order of magnitude larger than earlier efforts. The technique employs X-ray data from the ROSAT All-Sky Survey (RASS), and optical imaging and spectroscopic followup from the Sloan Digital Sky Survey (SDSS). Optical objects in the SDSS catalogs are automatically cross-correlated with RASS X-ray source positions; then priorities for follow-on SDSS optical spectra of candidate counterparts are automatically assigned using an algorithm based on the known fx/fopt ratios for various classes of X-ray emitters. SDSS parameters for optical morphology, magnitude, colors, plus FIRST radio data, serve as proxies for object class. Initial application of this approach to 1400 deg^2 of sky provides a catalog of 1200 spectroscopically confirmed quasars/AGN that are probable RASS identifications. Most of the IDs are new, and only a few percent of the AGN are likely to be random superpositions. The magnitude and redshift ranges of the counterparts extend over 15<m<21 and 0.03<z<3.6. Although most IDs are quasars and Sy 1s, a variety of other AGN subclasses are also sampled. Substantial numbers of rare AGN are found, including more than 130 narrow-line Seyfert 1s and 45 BL Lac candidates. These results already provide a sizeable set of new IDs, show utility of the sample in multi-waveband studies, and demonstrate the capability of the RASS/SDSS approach to efficiently proceed towards the largest homogeneously selected/observed sample of X-ray emitting AGN. Abridged Abstract

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

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.

USNO-B is an all-sky catalog that presents positions, proper motions, magnitudes in various optical passbands, and star/galaxy estimators for 1,042,618,261 objects derived from 3,643,201,733 separate observations. The data were obtained from scans of 7,435 Schmidt plates taken for the various sky surveys during the last 50 years. USNO-B1.0 is believed to provide all-sky coverage, completeness down to V = 21, 0.2 arcsecond astrometric accuracy at J2000, 0.3 magnitude photometric accuracy in up to five colors, and 85% accuracy for distinguishing stars from non-stellar objects. A brief discussion of various issues is given here, but the actual data are available from http://www.nofs.navy.mil and other sites.

Using a catalog of 147,986 galaxy redshifts and fluxes from the Sloan Digital Sky Survey (SDSS) we measure the galaxy luminosity density at z=0.1 in five optical bandpasses corresponding to the SDSS bandpasses shifted to match their restframe shape at z=0.1. We denote the bands {0.1}{u}, {0.1}{g}, {0.1}{r}, {0.1}{i}, {0.1}{z}, with \lambda_{eff} = [3216, 4240, 5595, 6792, 8111] Angstroms respectively. We use a maximum likelihood method which allows for a general form for the shape of the luminosity function, simple luminosity and number evolution, incorporates flux uncertainties, and accounts for the flux limits of the survey. We find luminosity densities at z=0.1 in absolute AB magnitudes in a Mpc^3 of [-14.10 \pm 0.15, -15.18 \pm 0.03, -15.90 \pm 0.03, -16.24 \pm 0.03, -16.56 \pm 0.02] in [{0.1}{u}, {0.1}{g}, {0.1}{r}, {0.1}{i}, {0.1}{z}], respectively, using \Omega_0 =0.3, \Omega_\Lambda=0.7, and h=1, and using Petrosian magnitudes. Similar results are obtained using Sersic model magnitudes, suggesting that flux from outside the Petrosian apertures is not a major correction. In the {0.1}{r} band, the best fit Schechter function to our results has \phi_\ast = (1.49 \pm 0.04) \times 10^{-2} h^3 Mpc^{-3}, M_\ast - 5\log_{10} h = -20.44 \pm 0.01, and \alpha = -1.05\pm 0.01. In solar luminosities, the luminosity density in {0.1}{r} is (1.84 \pm 0.04) h 10^8 L_{{0.1}{r},\odot} Mpc^{-3}. Our results are consistent with other estimates of the luminosity density, from the Two-degree Field Galaxy Redshift Survey and the Millenium Galaxy Catalog. They represent a substantial change (\sim 0.5 mag) from earlier SDSS luminosity density results based on commissioning data, almost entirely because of the inclusion of evolution in the luminosity function model.

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.

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.

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

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.

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.

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.

Through use of multiband (U, B, R, I) photometry we have isolated high redshift (3.0<z<3.5) galaxy candidates in a survey of 1.27 deg^2 to R = 21.25 and a survey of 0.02 deg^2 to R = 23.5. Our pool of candidates constrains the nature of the 3.0 < z < 3.5 luminosity function over the range L* < L < 100 L*, if we grant a similar level of completeness to these data as for very faint samples (to R = 25.5) selected in a similar fashion. Our constraints agree with the high redshift sky density at R = 20.5 estimated from Yee et al.'s (1996) serendipitous discovery of a bright, z = 2.7 galaxy, as well as the density at R ~ 23 by Steidel et al. (1996b). We strongly rule out -- by more than two orders of magnitude at M(R) = -25 -- the L > L* luminosity function for z = 3-5 galaxies obtained by a photometric redshift analysis of the Hubble Deep Field (HDF) by Gwyn & Hartwick (1996). Our results at R ~ 23 are more consistent with the photometric redshift analysis of the faint HDF galaxies by Sawicki & Yee (1996), but our present upper limits at the brightest magnitudes (R < 21.5, M(R) < -24) allow more generous volume densities of these super-L* galaxies.

As a means of better understanding the evolution of optically selected galaxies we consider the distribution of galaxies within the multicolor space $U$, $B_J$, $R_F$ and $I_N$. We find that they form an almost planar distribution out to $B_J =22.5$ and $z<0.3$. The position of a galaxy within this plane is dependent on its redshift, luminosity and spectral type. While in the original $U$, $B_J$, $R_F$ and $I_N$ space these properties are highly correlated we can define an optimal rotation of the photometric axes that makes much of this information orthogonal. Fitting the observed spectroscopic redshifts with a quadratic function of the four magnitudes we show that redshifts for galaxies can be estimated to an accuracy better than $\Delta z =0.05$. This dispersion is due to the photometric uncertainties within the photographic data. Assuming no galaxy evolution we derive a set of simulated galaxy fluxes in the U, J, F and N passbands. Using these data we investigate how the redshift is encoded within the broadband magnitudes and the intrinsic dispersion of the photometric-redshift relation. We find that the signal that defines a galaxy's photometric redshift is not related to specific absorption or emission lines but comes from the break in the overall shape of the galaxy continuum at around 4000 \AA. Using high signal-to-noise photometric data we estimate that it is possible to achieve an intrinsic dispersion of less than $\Delta z =0.02$.