We present J-H-K' photometry for a sample of 45 high redshift quasars found by the Sloan Digital Sky Survey. The sample was originally selected on the basis of optical colors and spans a redshift range from 3.6 to 5.03. Our photometry reflects the rest-frame SED longward of Ly alpha for all redshifts. The results show that the near-IR colors of high redshift quasars are quite uniform. We have modelled the continuum shape of the quasars (from just beyond Ly alpha to ~4000 A) with a power law of the form f_nu \propto nu^alpha, and find <alpha > =-0.57 with a scatter of 0.33. This value is similar to what is found for lower redshift quasars over the same restframe wavelength range, and we conclude that there is hardly any evolution in the continuum properties of optically selected quasars up to redshift 5. The spectral indices found by combining near-IR with optical photometry are in general consistent but slightly flatter than what is found for the same quasars using the optical spectra and photometry alone, showing that the continuum region used to determine the spectral indices can somewhat influence the results.

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.

The angular distribution of galaxies encodes a wealth of information about large scale structure. Ultimately, the Sloan Digital Sky Survey (SDSS) will record the angular positions of order 10^8 galaxies in five bands, adding significantly to the cosmological constraints. This is the first in a series of papers analyzing a rectangular stripe 2.5x90 degrees from early SDSS data. We present the angular correlation function for galaxies in four separate magnitude bins on angular scales ranging from 0.003 degrees to 15 degrees. Much of the focus of this paper is on potential systematic effects. We show that the final galaxy catalog -- with the mask accounting for regions of poor seeing, reddening, bright stars, etc. -- is free from external and internal systematic effects for galaxies brighter than r* = 22. Our estimator of the angular correlation function includes the effects of the integral constraint and the mask. The full covariance matrix of errors in these estimates is derived using mock catalogs with further estimates using a number of other methods.

We describe the algorithm that selects the main sample of galaxies for spectroscopy in the Sloan Digital Sky Survey from the photometric data obtained by the imaging survey. Galaxy photometric properties are measured using the Petrosian magnitude system, which measures flux in apertures determined by the shape of the surface brightness profile. The metric aperture used is essentially independent of cosmological surface brightness dimming, foreground extinction, sky brightness, and the galaxy central surface brightness. The main galaxy sample consists of galaxies with r-band Petrosian magnitude r < 17.77 and r-band Petrosian half-light surface brightness < 24.5 magnitudes per square arcsec. These cuts select about 90 galaxy targets per square degree, with a median redshift of 0.104. We carry out a number of tests to show that (a) our star-galaxy separation criterion is effective at eliminating nearly all stellar contamination while removing almost no genuine galaxies, (b) the fraction of galaxies eliminated by our surface brightness cut is very small (0.1%), (c) the completeness of the sample is high, exceeding 99%, and (d) the reproducibility of target selection based on repeated imaging scans is consistent with the expected random photometric errors. (abridged)

We present new VLT spectroscopic observations of the most distant quasar known, SDSS J1030+0524 at z=6.28 which was recently discovered by the Sloan Digital Sky Survey. We confirm the presence of a complete Gunn-Peterson trough caused by neutral hydrogen in the intergalactic medium. There is no detectable flux over the wavelength range from 8450 to 8710 A. We set an improved limit on the drop of the flux level blueward of the Ly alpha line: a factor of > 200. Below 8450 A the spectrum shows a rise in flux, with a large fraction (> 60 %) of the total emission produced by few narrow features of transmitted flux. We discuss the obvious proximity effect around this quasar, with the presence of transmitted flux with many absorption features in a region of about 23h^{-1} comoving Mpc. If assuming the surrounding medium were completely neutral, the size of this region would imply a quasar lifetime of \~1.3x10^7 years. We also present near-IR spectroscopy of both SDSS J1030+0524 and of SDSS J1306+05, the second most distant quasar known at redshift 6.0. We combine measurements of the CIV line and limits on the HeII emission with the NV line measurements from the optical spectra to derive line ratios, and by implication the abundances of these early quasar environments. The results are indistinguishable from those of lower redshift quasars and indicate little or no evolution in the abundances from z ~ 6 to z ~ 2. The line ratios suggest supersolar metallicities, implying that the first stars around the quasars must have formed at least a few hundreds of Myrs prior to the observation, i.e. at redshift higher than 8.

We compute the angular power spectrum C_l from 1.5 million galaxies in early SDSS data on large angular scales, l<600. The data set covers about 160 square degrees, with a characteristic depth of order 1 Gpc/h in the faintest (21<r<22) of our four magnitude bins. Cosmological interpretations of these results are presented in a companion paper by Dodelson et al (2001). The data in all four magnitude bins are consistent with a simple flat ``concordance'' model with nonlinear evolution and linear bias factors of order unity. Nonlinear evolution is particularly evident for the brightest galaxies. A series of tests suggest that systematic errors related to seeing, reddening, etc., are negligible, which bodes well for the sixtyfold larger sample that the SDSS is currently collecting. Uncorrelated error bars and well-behaved window functions make our measurements a convenient starting point for cosmological model fitting.

We study the process of cosmic reionization and estimate the ionizing background in the IGM using the Lyman series absorption in the spectra of the four quasars at 5.7<z<6.3 discovered by the SDSS. We derive the evolution of the ionizing background at high redshifts, using both semi-analytic techniques and cosmological simulations to model the density fluctuations in the IGM. The existence of the complete Ly alpha Gunn-Peterson trough in the spectrum of the z=6.28 quasar SDSS 1030+0524 indicates a photoionization rate Gamma_{-12} at z~6 lower than 0.08, at least a factor of 6 smaller than the value at z~3. The Ly beta and Ly gamma Gunn-Peterson troughs give an even stronger limit Gamma_{-12}<0.02 due to their smaller oscillator strengths, indicating that the ionizing background in the IGM at z~6 is more than 20 times lower than that at z~3. Meanwhile, the volume-averaged neutral hydrogen fraction increases from 10^{-5} at z~3 to >10^{-3} at z~6. At this redshift, the mass-averaged neutral hydrogen fraction is larger than 1%; the mildly overdense regions (delta > 3) are still mostly neutral and the comoving mean free path of ionizing photons is shorter than 8 Mpc. Comparison with simulations of cosmological reionization shows that the observed properties of the IGM at z~6 are typical of those in the era at the end of the overlap stage of reionization when the individual HII regions merge. Thus, z~6 marks the end of the reionization epoch. The redshift of reionization constrains the small scale power of the mass density fluctuations and the star forming efficiency of the first generation of objects.

We compare two methods for deriving distances to early-type galaxies: fundamental plane (FP) and surface brightness fluctuations (SBF) distances for 170 galaxies. A third set of distances is provided by predictions derived from the density field of the IRAS redshift survey. Overall there is good agreement. However, several nearby, low-luminosity, mainly S0 galaxies have systematically low FP distances; we conclude that this is at least partly due to recent star formation and consequently low mass-to-light ratios. The tie between the Cepheid-calibrated SBF distances (Mpc) and the far-field calibrated FP distances (km/s) yields a Hubble constant H_0 = 68 +/- 3 km/s/Mpc, while the comparison between SBF and the IRAS-reconstructed distances yields H_0 = 74 +/- 2 (internal errorbars). Possible explanations for the marginal inconsistency include systematic errors in the redshift survey completeness estimates or in the FP aperture corrections, but at this point, the best estimate of H_0 from early-type galaxies may be a simple average of the above two estimates. After revising the SBF distances to be in agreement with the final set of Key Project Cepheid distances, we conclude H_0 = 73 +/- 4 +/- 11, where the second errorbar represents the total systematic uncertainty in the distance zero point. We also discuss the recently introduced "fluctuation star count" parameter Nbar as a less demanding alternative to (V-I) for calibrating SBF distances. The Nbar-calibrated SBF method is akin to a hybrid SBF-FP distance indicator, and we find that the use of Nbar actually improves the SBF distances. Further study of Nbar would provide useful constraints for elliptical galaxy formation models. (abridged)

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 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 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 measurements of parameters of the 3-dimensional power spectrum of galaxy clustering from 222 square degrees of early imaging data in the Sloan Digital Sky Survey. The projected galaxy distribution on the sky is expanded over a set of Karhunen-Loeve eigenfunctions, which optimize the signal-to-noise ratio in our analysis. A maximum likelihood analysis is used to estimate parameters that set the shape and amplitude of the 3-dimensional power spectrum. Our best estimates are Gamma=0.188 +/- 0.04 and sigma_8L = 0.915 +/- 0.06 (statistical errors only), for a flat Universe with a cosmological constant. We demonstrate that our measurements contain signal from scales at or beyond the peak of the 3D power spectrum. We discuss how the results scale with systematic uncertainties, like the radial selection function. We find that the central values satisfy the analytically estimated scaling relation. We have also explored the effects of evolutionary corrections, various truncations of the KL basis, seeing, sample size and limiting magnitude. We find that the impact of most of these uncertainties stay within the 2-sigma uncertainties of our fiducial result.

The Sloan Digital Sky Survey is one of the first multicolor photometric and spectroscopic surveys designed to measure the statistical properties of galaxies within the local Universe. In this Letter we present some of the initial results on the angular 2-point correlation function measured from the early SDSS galaxy data. The form of the correlation function, over the magnitude interval 18<r*<22, is shown to be consistent with results from existing wide-field, photographic-based surveys and narrower CCD galaxy surveys. On scales between 1 arcminute and 1 degree the correlation function is well described by a power-law with an exponent of ~ -0.7. The amplitude of the correlation function, within this angular interval, decreases with fainter magnitudes in good agreement with analyses from existing galaxy surveys. There is a characteristic break in the correlation function on scales of approximately 1-2 degrees. On small scales, < 1', the SDSS correlation function does not appear to be consistent with the power-law form fitted to the 1'< theta <0.5 deg data. With a data set that is less than 2% of the full SDSS survey area, we have obtained high precision measurements of the power-law angular correlation function on angular scales 1' < theta < 1 deg, which are robust to systematic uncertainties. Because of the limited area and the highly correlated nature of the error covariance matrix, these initial results do not yet provide a definitive characterization of departures from the power-law form at smaller and larger angles. In the near future, however, the area of the SDSS imaging survey will be sufficient to allow detailed analysis of the small and large scale regimes, measurements of higher-order correlations, and studies of angular clustering as a function of redshift and galaxy type.

Early photometric data from the Sloan Digital Sky Survey (SDSS) contain angular positions for 1.5 million galaxies. In companion papers, the angular correlation function $w(\theta)$ and 2D power spectrum $C_l$ of these galaxies are presented. Here we invert Limber's equation to extract the 3D power spectrum from the angular results. We accomplish this using an estimate of $dn/dz$, the redshift distribution of galaxies in four different magnitude slices in the SDSS photometric catalog. The resulting 3D power spectrum estimates from $w(\theta)$ and $C_l$ agree with each other and with previous estimates over a range in wavenumbers $0.03 < k/{\rm h Mpc}^{-1} < 1$. The galaxies in the faintest magnitude bin ($21 < \rstar < 22$, which have median redshift $z_m=0.43$) are less clustered than the galaxies in the brightest magnitude bin ($18 < \rstar < 19$ with $z_m=0.17$), especially on scales where nonlinearities are important. The derived power spectrum agrees with that of Szalay et al. (2001) who go directly from the raw data to a parametric estimate of the power spectrum. The strongest constraints on the shape parameter $\Gamma$ come from the faintest galaxies (in the magnitude bin $21 < \rstar < 22$), from which we infer $\Gamma = 0.14^{+0.11}_{-0.06}$ (95% C.L.).

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.

The Sloan Digital Sky Survey has discovered a population of broad absorption line quasars with various extreme properties. Many show absorption from metastable states of FeII with varying excitations; several objects are almost completely absorbed bluewards of MgII; at least one shows stronger absorption from FeIII than FeII, indicating temperatures T>35000 K in the absorbing region; and one object even seems to have broad H-beta absorption. Many of these extreme BALs are also heavily reddened, though `normal' BALs (particularly LoBALs) from SDSS also show evidence for internal reddening.

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 compare the properties of recent samples of the lithium abundances in halo stars to one another and to the predictions of theoretical models including rotational mixing, and we examine the data for trends with metal abundance. We find from a KS test that in the absence of any correction for chemical evolution, the Ryan, Norris, & Beers (1999} sample is fully consistent with mild rotational mixing induced depletion and, therefore, with an initial lithium abundance higher than the observed value. Tests for outliers depend sensitively on the threshold for defining their presence, but we find a 10$--$45% probability that the RNB sample is drawn from the rotationally mixed models with a 0.2 dex median depletion (with lower probabilities corresponding to higher depletion factors). When chemical evolution trends (Li/H versus Fe/H) are treated in the linear plane we find that the dispersion in the RNB sample is not explained by chemical evolution; the inferred bounds on lithium depletion from rotational mixing are similar to those derived from models without chemical evolution. We find that differences in the equivalent width measurements are primarily responsible for different observational conclusions concerning the lithium dispersion in halo stars. The standard Big Bang Nucleosynthesis predicted lithium abundance which corresponds to the deuterium abundance inferred from observations of high-redshift, low-metallicity QSO absorbers requires halo star lithium depletion in an amount consistent with that from our models of rotational mixing, but inconsistent with no depletion.

The ubiquity of supermassive black holes (SMBHs) at the centers of nearby luminous galaxies can arise from the multiple mergers experienced by dark matter halos in hierarchical structure formation models, even if only a small fraction of these galaxies harbor SMBHs at high redshifts. We illustrate this possibility using cosmological Monte Carlo simulations of the merger history of dark matter halos and their associated SMBHs. In our most extreme models, in order to populate nearly every bright galaxy with a SMBH at z=0, only a few percent of the halos with virial temperatures above 10^4 K are required to harbor a SMBH at high redshift. This possibility must be included in studies of the luminosity function and the clustering properties of quasars. We predict the number of SMBH merger events that are detectable by the gravitational wave experiment LISA, as a function of redshift, out to z=5. Although the event rates can be significantly reduced in scenarios with rare SMBHs, a minimum of \~10 detectable merger events per year is predicted if SMBH binaries coalesce efficiently. The observed distribution of events with redshift could yield valuable information on the SMBH formation process. If SMBH binaries do not coalesce, we find that at least several SMBH slingshot ejections probably occurred from z=5 to the present in each galaxy more massive than ~10^11 Msun at z=0. Although our results are sensitive to the minimum cooling mass assumed for the formation of SMBHs, we expect the qualitative predictions of our models to be robust.

Measurements of the cosmological density parameter Omega using techniques that exploit the gravity-induced motions of galaxies constrain, in linear perturbation theory, the degenerate parameter combination beta = Omega^{0.6}/b, where the linear bias parameter b is the ratio of the fluctuation amplitudes of the galaxy and mass distributions. However, the relation between the galaxy and mass density fields depends on the complex physics of galaxy formation, and it can in general be non-linear, stochastic, and perhaps non-local. The one-parameter linear bias model is almost certainly oversimplified, which leads to the obvious question: What is the quantity beta that is actually measured by different techniques? To address this question, we estimate beta from galaxy distributions that are constructed by applying a variety of locally biased galaxy formation models to cosmological N-body simulations. We compare the values of beta estimated using three different techniques: a density-density comparison similar to the POTENT analysis, a velocity-velocity comparison similar to the VELMOD analysis, and an anisotropy analysis of the redshift-space power spectrum. In most cases, we find that beta estimated using all three methods is similar to the asymptotic value of Omega^{0.6}/b_{sigma}(R) at large R, where b_{sigma}(R) is the ratio of rms galaxy fluctuations to rms mass fluctuations on scale R. Thus, something close to the conventional interpretation of beta continues to hold even for complex bias models. Moreover, we find that beta estimates made using these three methods should, in principle, agree with each other. It is thus unlikely that non-linear or scale-dependent bias is responsible for the discrepancies that exist among current measurements of beta from different techniques.

High resolution N-body simulations of cold dark matter (CDM) models predict that galaxies and clusters have cuspy halos with excessive substructure. Observations reveal smooth halos with central density cores. One possible resolution of this conflict is that the dark matter is warm (WDM); this will suppress the power spectrum on small scales. The Lyman-alpha forest is a powerful probe of the linear power spectrum on these scales. We use collisionless N-body simulations to follow the evolution of structure in WDM models, and analyze artificial Lyman-alpha forest spectra extracted from them. By requiring that there is enough small-scale power in the linear power spectrum to reproduce the observed properties of the Lyman-alpha forest in quasar spectra, we derive a lower limit to the mass of the WDM particle of 750 eV. This limit is robust to reasonable uncertainties in our assumption about the temperature of the mean density gas (T0) at z=3. We argue that any model that suppresses the CDM linear theory power spectrum more severely than a 750 eV WDM particle cannot produce the Lyman-alpha forest.

We show that a regression of unsmoothed peculiar velocity measurements against peculiar velocities predicted from a smoothed galaxy density field leads to a biased estimate of the cosmological density parameter Omega, even when galaxies trace the underlying mass distribution and galaxy positions and velocities are known perfectly. The bias arises because the errors in the predicted velocities are correlated with the predicted velocities themselves. We investigate this bias using cosmological N-body simulations and analytic arguments. In linear perturbation theory, for cold dark matter power spectra and Gaussian or top hat smoothing filters, the bias in Omega is always positive, and its magnitude increases with increasing smoothing scale. This linear calculation reproduces the N-body results for Gaussian smoothing radii R_s > 10 Mpc/h, while non-linear effects lower the bias on smaller smoothing scales, and for R_s < 3 Mpc/h Omega is underestimated rather than overestimated. The net bias in Omega for a given smoothing filter depends on the underlying cosmological model. The effect on current estimates of Omega from velocity-velocity comparisons is probably small relative to other uncertainties, but taking full advantage of the statistical precision of future peculiar velocity data sets will require either equal smoothing of the predicted and measured velocity fields or careful accounting for the biases discussed here.

We present the results of reconstruction analysis of the galaxy distribution in a spherical region of radius 50 Mpc/h centered on the Local Group, as mapped by the IRAS Point Source Catalog Redshift Survey (PSCz). We reconstruct it using 15 different models for structure formation in the universe, each consisting of a set of assumptions regarding the value of the density parameter Omega_m, and the amplitude and nature of the biasing between IRAS galaxies and mass. For every model, we also reconstruct ten mock PSCz catalogs derived from the outputs of numerical simulations that have the appropriate assumptions about Omega_m and bias. We use a variety of statistics to compare the accuracy of each reconstruction of the PSCz catalog to the accuracy expected based on the corresponding mock catalog reconstructions. We find that gravitational instability of Gaussian mass density fluctuations can account for the galaxy distribution in the PSCz catalog, at least for some plausible assumptions about the value of Omega_m and the bias. Unbiased models in which IRAS galaxies trace mass fail to reconstruct the PSCz catalog accurately, both for Omega_m=0.4 and for Omega_m=1. Low Omega_m models in which IRAS galaxies are antibiased with respect to the mass are the most successful in reconstructing the PSCz catalog. In particular, a model with Omega_m=0.4 and IRAS galaxies related to the mass according to the predictions of a semi-analytic galaxy formation model is very successful in reconstructing the PSCz galaxy distribution.

Measurements of the cosmological density parameter (Omega) using techniques that exploit the gravity-induced motions of galaxies yield, in linear perturbation theory, the degenerate parameter combination beta=Omega^{0.6}/b, where the linear bias parameter b is the ratio of the fluctuation amplitudes of the galaxy and mass distributions. However, the relation between the mass and the galaxy density fields depends on the complex physics of galaxy formation, and it can in general be non-linear, non-local, and stochastic. There is a growing consensus that the one-parameter model for bias is oversimplified. This leaves us with the following question: What is the quantity that is actually being measured by the different techniques? In order to address this question, we present estimates of beta from galaxy distributions constructed by applying a variety of biased galaxy formation models to cosmological N-body simulations. We compare the values of beta estimated using two different techniques: the anisotropy of the redshift-space power spectrum and an idealized version of the POTENT method. In most cases, we find that the bias factor b=Omega^{0.6}/beta derived from redshift-space anisotropy or POTENT is similar to the large-scale value of b_{sigma}(R) defined by the ratio of rms fluctuation amplitudes of the galaxy and mass distributions. However, non-linearity of bias and residual effects of non-linear gravitational evolution both influence beta estimates at the 10-20% level.

We propose and apply a test that can detect any systematic errors in the Hipparcos parallaxes towards the Hyades cluster at the level of 0.3 mas. We show that the statistical parallax method subsumes the classical moving cluster methods and provides more accurate and robust estimates of the distance and the first two moments of the velocity distribution of the Hyades cluster namely, its bulk space velocity and the velocity dispersion tensor. We predict the parallaxes of Hyades cluster members using the common cluster space velocity derived from the statistical parallax method and their individual Hipparcos proper motions. We show that the parallaxes determined in this manner (pi_pm) are consistent at the 1-sigma level with the parallaxes (pi_orb) of three Hyades spectroscopic binary systems with orbital solutions. We find that <(pi_pm - pi_orb)> = 0.52 +- 0.47 mas, where the error is dominated by the errors in the orbital parallaxes. A reduction in this error would allow a test of the systematic errors in the Hipparcos parallaxes at the 0.3 mas level. If the Hyades distance scale is fixed by Hipparcos parallaxes, then its bulk velocity in equatorial coordinates is (Vx,Vy,Vz) = (-5.70 +- 0.20, 45.62 +- 0.11, 5.65 +- 0.08)km/s, its velocity dispersion is 320 +- 39 m/s, and the distance modulus to the centroid of our sample of 43 cluster members is 3.34 +- 0.02 mag.

We show that the errors in the Hipparcos parallaxes towards the Pleiades and the Hyades open clusters are spatially correlated over angular scales of 2 to 3 deg, with an amplitude of up to 2 mas. This correlation is stronger than expected based on the analysis of the Hipparcos catalog. We predict the parallaxes of individual cluster members, pi_pm, from their Hipparcos proper motions, assuming that all cluster members have the same space velocity. We compare pi_pm with their Hipparcos parallaxes, pi_Hip, and find that there are significant spatial correlations in pi_Hip. We derive a distance modulus to the Pleiades of 5.58 +- 0.18 mag using the radial-velocity gradient method. This value, agrees very well with the distance modulus of 5.60 +- 0.04 mag determined using the main-sequence fitting technique, compared with the value of 5.33 +- 0.06 inferred from the average of the Hipparcos parallaxes of the Pleiades members. We show that the difference between the main-sequence fitting distance and the Hipparcos parallax distance can arise from spatially correlated errors in the Hipparcos parallaxes of individual Pleiades members. Although the Hipparcos parallax errors towards the Hyades are spatially correlated in a manner similar to those of the Pleiades, the center of the Hyades is located on a node of this spatial structure. Therefore, the parallax errors cancel out when the average distance is estimated, leading to a mean Hyades distance modulus that agrees with the pre-Hipparcos value. We speculate that these spatial correlations are also responsible for the discrepant distances that are inferred using the mean Hipparcos parallaxes to some open clusters. Finally, we note that our conclusions are based on a purely geometric method and do not rely on any models of stellar isochrones.

We examine the influence of the morphology-density(MD) relation and a wide range of simple models for biased galaxy formation on statistical measures of large scale structure. We contrast the behavior of local biasing models, in which the efficiency of galaxy formation is determined by density, geometry, or velocity dispersion of the local mass distribution, with that of non-local biasing models, in which galaxy formation is modulated coherently over scales larger than the galaxy correlation length. If morphological segregation of galaxies is governed by a local MD relation, then the correlation function of E/S0 galaxies should be steeper and stronger than that of spiral galaxies on small scales, as observed, while on large scales the correlation functions of E/S0 and spiral galaxies should have the same shape but different amplitudes. Similarly, all of our local bias models produce scale-independent amplification of the correlation function and power spectrum in the linear and mildly non-linear regimes; only a non-local biasing mechanism can alter the shape of the power spectrum on large scales. Moments of the biased galaxy distribution retain the hierarchical pattern of the mass moments, but biasing alters the values and scale-dependence of the hierarchical amplitudes S3 and S4. Pair-weighted moments of the galaxy velocity distribution are sensitive to the details of the biasing prescription. The non-linearity of the relation between galaxy density and mass density depends on the biasing prescription and the smoothing scale, and the scatter in this relation is a useful diagnostic of the physical parameters that determine the bias. Although the sensitivity of galaxy clustering statistics to the details of biasing is an obstacle to testing cosmological models, it is an asset for testing galaxy formation theories.

We present a comparative study of six different methods for reversing the gravitational evolution of a cosmological density field to recover the primordial fluctuations: linear theory, the Gaussianization mapping scheme, two different quasi-linear dynamical schemes based on the Zel'dovich approximation, a Hybrid dynamical-Gaussianization method and the Path Interchange Zel'dovich Approximation (PIZA). The final evolved density field from an N-body simulation constitutes our test case. We use a variety of statistical measures to compare the initial density field recovered from it to the true initial density field, using each of the six different schemes. These include point-by-point comparisons of the density fields in real space, the individual modes in Fourier space, as well as global statistical properties such as the genus, the PDF of the density, and the distribution of peak heights and their shapes. We find linear theory to be the most inaccurate of all the schemes. The Gaussianization scheme is the least accurate after linear theory. The two quasi-linear dynamical schemes are more accurate than Gaussianization, although they break down quite drastically when used outside their range of validity - the quasi-linear regime. The complementary beneficial aspects of the dynamical and the Gaussianization schemes are combined in the Hybrid method. We find this Hybrid scheme to be more accurate and robust than either Gaussianization or the dynamical method alone. The PIZA scheme performs substantially better than the others in all point-by-point comparisons. However, it produces an oversmoothed initial density field, with a smaller number of peaks than expected, but recovers the PDF of the initial density with impressive accuracy on scales as small as 3Mpc/h.

In reconstruction analysis of galaxy redshift surveys, one works backwards from the observed galaxy distribution to the primordial density field in the same region, then evolves the primordial fluctuations forward in time with an N-body code. This incorporates assumptions about the cosmological parameters, the properties of primordial fluctuations, and the biasing relation between galaxies and mass. These can be tested by comparing the reconstruction to the observed galaxy distribution, and to peculiar velocity data. This paper presents a hybrid reconstruction method that combines the `Gaussianization'' technique of Weinberg(1992) with the dynamical schemes of Nusser & Dekel(1992) and Gramann(1993). We test the method on N-body simulations and on N-body mock catalogs that mimic the depth and geometry of the Point Source Catalog Redshift Survey and the Optical Redshift Survey. This method is more accurate than Gaussianization or dynamical reconstruction alone. Matching the observed morphology of clustering can limit the bias factor b, independent of Omega. Matching the cluster velocity dispersions and z-space distortions of the correlation function xi(s,mu) constrains the parameter beta=Omega^{0.6}/b. Relative to linear or quasi-linear approximations, a fully non-linear reconstruction makes more accurate predictions of xi(s,mu) for a given beta, thus reducing the systematic biases of beta measurements and offering further scope for breaking the degeneracy between Omega and b. It also circumvents the cosmic variance noise that limits conventional analyses of xi(s,mu). It can also improve the determination of Omega and b from joint analyses of redshift & peculiar velocity surveys as it predicts the fully non-linear peculiar velocity distribution at each point in z-space.

The depletion of lithium during the pre-main sequence and main sequence phases of stellar evolution plays a crucial role in the comparison of the predictions of big bang nucleosynthesis with the abundances observed in halo stars. In the past a wide range of possible depletion factors, ranging from minimal in standard (non-rotating) stellar models to as much as an order of magnitude in models including rotational mixing have been suggested. Using recent progress in the study of the angular momentum evolution of low mass stars, which now permits the construction of theoretical models that reproduce the angular momentum evolution of low mass open cluster stars, along with the observed distribution of initial angular momenta inferred from stellar rotation data in young open clusters, we study lithium depletion in main sequence halo stars. We predict a well-defined, nearly flat (in T_eff) halo lithium plateau with modest scatter and a small population of outliers. We also examine the relative depletions of 6Li and 7Li finding that the dispersion in the plateau and the 6Li/7Li depletion ratio scale with the absolute 7Li depletion. Using observational data to bound the 7Li depletion in main sequence halo stars, we find a maximum depletion of 0.4 dex is set by the observed dispersion and the 6Li/7Li depletion ratio, and a minimum depletion of 0.2 dex is required by both the presence of highly overdepleted halo stars and consistency with the solar and open cluster 7Li data. The cosmological implications of these bounds on the primordial abundance of 7Li are discussed. (Abridged)

There has been recent progress in the study of the angular momentum evolution of low mass stars (Krishnamurthi et al 1997a). Theoretical models can now be constructed which reproduce the angular momentum evolution of low mass open cluster stars and the distribution of initial conditions can be inferred from young clusters. In this poster we report on the application of these models to the problem of rotational mixing in halo stars. The distribution of initial conditions inferred from young clusters produces a well-defined halo lithium "plateau" with modest scatter and a small population of outliers. Different choices for the solar calibration produce a range of absolute depletion factors. We show that both the dispersion and the ratio of $^6$Li depletion to $^7$Li depletion increase as the absolute $^7$Li depletion increases. The measured $^6$Li in HD 84937 and the dispersion in the plateau set independent upper bounds on the $^7$Li depletion. Consistency with open clusters and the Sun, along with claims of an intrinsic dispersion in the plateau, set a lower bound. We derive a range of 0.2-0.4 dex $^7$Li depletion in halo field stars. Implications for cosmology are discussed.

We present H and K photometry of 42,000 stars in an area of 250 arcmin$^{2}$ centered on the Galactic center. We use the photometry to construct a dereddened K band luminosity function (LF) for this region, excluding the excessively crowded inner 2' of the Galaxy. This LF is intermediate between the LF of Baade's window and the LF of inner 2' of the Galactic center. We speculate that the bright stars in this region have an age which is intermediate between the starburst population in the Galactic center and the old bulge population. We present the coordinates and mags for 16 stars with K_{0} < 5 for spectroscopic follow up.

Lunar occultation can be used to measure the proper motions of some of the long time scale microlensing events, $t_{e} \gsim 70$ days, now being detected toward the Galactic bulge. The long events are difficult to explain within the context of standard models of the mass distribution and dynamics of the Galaxy. Han & Gould (1995b) have suggested that they may be due to a kinematically cold population near the Sun. To resolve the mass, distance, and velocity of individual events and so to determine their nature, one must measure parallaxes and proper motions. For long events, parallaxes can be often obtained from ground-based measurements, but proper motions can only rarely be determined using conventional methods. Lunar occultations are therefore key to the understanding of the long events. We carry out realistic simulations to estimate the uncertainty of these measurements and show that proper motions could be measured for about one long event per year.