The 2019 Nobel Prize in Physics honors three pioneering scientists for their fundamental contributions to basic cosmic questions - Professor James Peebles (Princeton University), Michel Mayor (University of Geneva), and Didier Queloz (University of Geneva and the University of Cambridge) - "for contributions to our understanding of the evolution of the universe and Earth's place in the cosmos," with one half to James Peebles "for theoretical discoveries in physical cosmology," and the other half jointly to Michel Mayor and Didier Queloz "for the discovery of an exoplanet orbiting a solar-type star." We summarize the historical and scientific backdrop to this year's Physics Nobel.

We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$\sigma$ significance. PSF also slightly flattens the measured colour profiles.

Aims. We wish to determine the distribution of dark matter halo masses as a function of the stellar mass and the stellar mass profile, for massive galaxies in the BOSS CMASS sample. Methods. We use grizy photometry from HSC to obtain S\'ersic fits and stellar masses of CMASS galaxies for which HSC weak lensing data is available, visually selected to have spheroidal morphology. We apply a cut in stellar mass, $\log{M_*/M_\odot} > 11.0$,selecting $\sim$10, 000 objects. Using a Bayesian hierarchical inference method, we first investigate the distribution of S\'ersic index and size as a function of stellar mass. Then, making use of shear measurements from HSC, we measure the distribution of halo mass as a function of stellar mass, size and S\'ersic index. Results. Our data reveals a steep stellar mass-size relation $R_e \propto M_*^{\beta_R}$, with $\beta_R$ larger than unity, and a positive correlation between S\'ersic index and stellar mass: $n \propto M_*^{0.46}$. Halo mass scales approximately with the 1.7 power of the stellar mass. We do not find evidence for an additional dependence of halo mass on size or S\'ersic index at fixed stellar mass. Conclusions. Our results disfavour galaxy evolution models that predict significant differences in the size growth efficiency of galaxies living in low and high mass halos.

Using deep images from the Hyper Suprime-Cam (HSC) survey and taking advantage of its unprecedented weak lensing capabilities, we reveal a remarkably tight connection between the stellar mass distribution of massive central galaxies and their host dark matter halo mass. Massive galaxies with more extended stellar mass distributions tend to live in more massive dark matter haloes. We explain this connection with a phenomenological model that assumes, (1) a tight relation between the halo mass and the total stellar content in the halo, (2) that the fraction of in-situ and ex-situ mass at $r<10$ kpc depends on halo mass. This model provides an excellent description of the stellar mass functions (SMF) of total stellar mass ($M_{\star}^{\rm Max}$) and stellar mass within inner 10 kpc ($M_{\star}^{10}$) and also reproduces the HSC weak lensing signals of massive galaxies with different stellar mass distributions. The best-fit model shows that halo mass varies significantly at fixed total stellar mass (as much as 0.4 dex) with a clear dependence on $M_{\star}^{10}$. Our two-parameter $M_{\star}^{\rm Max}$-$M_{\star}^{10}$ description provides a more accurate picture of the galaxy-halo connection at the high-mass end than the simple stellar-halo mass relation (SHMR) and opens a new window to connect the assembly history of halos with those of central galaxies. The model also predicts that the ex-situ component dominates the mass profiles of galaxies at $r< 10$ kpc for $\log M_{\star} \ge 11.7$). The code used for this paper is available online: https://github.com/dr-guangtou/asap

We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

Using $\sim$140 deg$^2$ Subaru Hyper Suprime-Cam (HSC) survey data, we stack the weak lensing (WL) signal around five Planck clusters found within the footprint. This yields a 15$\sigma$ detection of the mean Planck cluster mass density profile. The five Planck clusters span a relatively wide mass range, $M_{\rm WL,500c} = (2-30)\times10^{14}\,M_\odot/h$ with a mean mass of $M_{\rm WL,500c} = (4.15\pm0.61)\times10^{14}\,M_\odot/h$. The ratio of the stacked Planck Sunyaev-Zel'dovich (SZ) mass to the stacked WL mass is $ \langle M_{\rm SZ}\rangle/\langle M_{\rm WL}\rangle = 1-b = 0.80\pm0.14$. This mass bias is consistent with previous WL mass calibrations of Planck clusters within the errors. We discuss the implications of our findings for the calibration of SZ cluster counts and the much discussed tension between Planck SZ cluster counts and Planck $\Lambda$CDM cosmology.

We present a statistical study of the redshift evolution of the cluster galaxy population over a wide redshift range from 0.1 to 1.1, using $\sim 1900$ optically-selected CAMIRA clusters from $\sim 232$~deg$^2$ of the Hyper Suprime-Cam (HSC) Wide S16A data. Our stacking technique with a statistical background subtraction reveals color-magnitude diagrams of red-sequence and blue cluster galaxies down to faint magnitudes of $m_z\sim 24$. We find that the linear relation of red-sequence galaxies in the color-magnitude diagram extends down to the faintest magnitudes we explore with a small intrinsic scatter $\sigma_{\rm int}(g-r)<0.1$. The scatter does not evolve significantly with redshift. The stacked color-magnitude diagrams are used to define red and blue galaxies in clusters for studying their radial number density profiles without resorting to photometric redshifts of individual galaxies. We find that red galaxies are significantly more concentrated toward cluster centers and blue galaxies dominate the outskirt of clusters. We explore the fraction of red galaxies in clusters as a function of redshift, and find that the red fraction decreases with increasing distances from cluster centers. The red fraction exhibits a moderate decrease with increasing redshift. The radial number density profiles of cluster member galaxies are also used to infer the location of the steepest slope in the three dimensional galaxy density profiles. For a fixed threshold in richness, we find little redshift evolution in this location.

We present a measurement of the linear growth rate of structure, \textit{f} from the Sloan Digital Sky Survey III (SDSS III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) using Convolution Lagrangian Perturbation Theory (CLPT) with Gaussian Streaming Redshift-Space Distortions (GSRSD) to model the two point statistics of BOSS galaxies in DR12. The BOSS-DR12 dataset includes 1,198,006 massive galaxies spread over the redshift range $0.2 < z < 0.75$. These galaxy samples are categorized in three redshift bins. Using CLPT-GSRSD in our analysis of the combined sample of the three redshift bins, we report measurements of $f \sigma_8$ for the three redshift bins. We find $f \sigma_8 = 0.430 \pm 0.054$ at $z_{\rm eff} = 0.38$, $f \sigma_8 = 0.452 \pm 0.057$ at $z_{\rm eff} = 0.51$ and $f \sigma_8 = 0.457 \pm 0.052$ at $z_{\rm eff} = 0.61$. Our results are consistent with the predictions of Planck $\Lambda$CDM-GR. Our constraints on the growth rates of structure in the Universe at different redshifts serve as a useful probe, which can help distinguish between a model of the Universe based on dark energy and models based on modified theories of gravity. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. 2016 to produce the final cosmological constraints from BOSS.

We study the effects of filaments on galaxy properties in the Sloan Digital Sky Survey (SDSS) Data Release 12 using filaments from the `Cosmic Web Reconstruction' catalogue (Chen et al. 2016), a publicly available filament catalogue for SDSS. Since filaments are tracers of medium-to-high density regions, we expect that galaxy properties associated with the environment are dependent on the distance to the nearest filament. Our analysis demonstrates that a red galaxy or a high-mass galaxy tend to reside closer to filaments than a blue or low-mass galaxy. After adjusting the effect from stellar mass, on average, early-forming galaxies or large galaxies have a shorter distance to filaments than late-forming galaxies or small galaxies. For the Main galaxy sample (MGS), all signals are very significant ($>6\sigma$). For the LOWZ and CMASS sample, the stellar mass and size are significant ($>2 \sigma$). The filament effects we observe persist until $z = 0.7$ (the edge of the CMASS sample). Comparing our results to those using the galaxy distances from redMaPPer galaxy clusters as a reference, we find a similar result between filaments and clusters. Moreover, we find that the effect of clusters on the stellar mass of nearby galaxies depends on the galaxy's filamentary environment. Our findings illustrate the strong correlation of galaxy properties with proximity to density ridges, strongly supporting the claim that density ridges are good tracers of filaments.

The Sloan Digital Sky Survey has surveyed 14,555 square degrees of the sky, and delivered over a trillion pixels of imaging data. We present the large-scale clustering of 1.6 million quasars between z = 0.5 and z = 2.5 that have been classified from this imaging, representing the highest density of quasars ever studied for clustering measurements. This data set spans ~11,000 square degrees and probes a volume of 80(Gpc/h)^3. In principle, such a large volume and medium density of tracers should facilitate high-precision cosmological constraints. We measure the angular clustering of photometrically classified quasars using an optimal quadratic estimator in four redshift slices with an accuracy of ~25% over a bin width of l ~10 - 15 on scales corresponding to matter-radiation equality and larger (l ~ 2 - 30). Observational systematics can strongly bias clustering measurements on large scales, which can mimic cosmologically relevant signals such as deviations from Gaussianity in the spectrum of primordial perturbations. We account for systematics by employing a new method recently proposed by Agarwal et al. (2014) to the clustering of photometrically classified quasars. We carefully apply our methodology to mitigate known observational systematics and further remove angular bins that are contaminated by unknown systematics. Combining quasar data with the photometric luminous red galaxy (LRG) sample of Ross et al. (2011) and Ho et al. (2012), and marginalizing over all bias and shot noise-like parameters, we obtain a constraint on local primordial non-Gaussianity of fNL = -113+/-154 (1\sigma error). [Abridged]

Intrinsic alignments of galaxies have been shown to contaminate weak gravitational lensing observables on linear scales, $r>$ 10 $h^{-1}$Mpc, but studies of alignments in the non-linear regime have thus far been inconclusive. We present an estimator for extracting the intrinsic alignment signal of galaxies around stacked clusters of galaxies from multiband imaging data. Our estimator removes the contamination caused by galaxies that are gravitationally lensed by the clusters and scattered in redshift space due to photometric redshift uncertainties. It uses posterior probability distributions for the redshifts of the galaxies in the sample and it is easily extended to obtain the weak gravitational lensing signal while removing the intrinsic alignment contamination. We apply this algorithm to groups and clusters of galaxies identified in the Sloan Digital Sky Survey `Stripe 82' coadded imaging data over $\sim 150$ deg$^2$. We find that the intrinsic alignment signal around stacked clusters in the redshift range $0.1<z<0.4$ is consistent with zero. In terms of the tidal alignment model of Catelan et al. (2001), we set joint constraints on the strength of the alignment and the bias of the lensing groups and clusters on scales between 0.1 and $10\,h^{-1}$ Mpc, $b_LC_1\rho_{\rm crit} = -2_{-14}^{+14} \times 10^{-4}$. This constrains the contamination fraction of alignment to lensing signal to the range between $[-18,23]$ per cent below scales of 1 $h^{-1}$ Mpc at 95 per cent confidence level, and this result depends on our photometric redshift quality and selection criteria used to identify background galaxies. Our results are robust to the choice of photometric band in which the shapes are measured ($i$ and $r$) and to centring on the Brightest Cluster Galaxy or on the geometrical centre of the clusters.

We investigate the luminosity and colour dependence of clustering of CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The halo occupation distribution framework is adopted to model the projected two-point correlation function measurements on small and intermediate scales (from $0.02$ to $60\,h^{-1}{\rm {Mpc}}$) and to interpret the observed trends and infer the connection of galaxies to dark matter halos. We find that luminous red galaxies reside in massive halos of mass $M{\sim}10^{13}$--$10^{14}\,h^{-1}{\rm M_\odot}$ and more luminous galaxies are more clustered and hosted by more massive halos. The strong small-scale clustering requires a fraction of these galaxies to be satellites in massive halos, with the fraction at the level of 5--8 per cent and decreasing with luminosity. The characteristic mass of a halo hosting on average one satellite galaxy above a luminosity threshold is about a factor $8.7$ larger than that of a halo hosting a central galaxy above the same threshold. At a fixed luminosity, progressively redder galaxies are more strongly clustered on small scales, which can be explained by having a larger fraction of these galaxies in the form of satellites in massive halos. Our clustering measurements on scales below $0.4\,h^{-1}{\rm {Mpc}}$ allow us to study the small-scale spatial distribution of satellites inside halos. While the clustering of luminosity-threshold samples can be well described by a Navarro-Frenk-White (NFW) profile, that of the reddest galaxies prefers a steeper or more concentrated profile. Finally, we also use galaxy samples of constant number density at different redshifts to study the evolution of luminous galaxies, and find the clustering to be consistent with passive evolution in the redshift range of $0.5 \lesssim z \lesssim 0.6$.

We use high-resolution cosmological hydrodynamic galaxy formation simulations to gain insights on how galaxies lose their cold gas at low redshift as they migrate from the field to the high density regions of clusters of galaxies. We find that beyond three cluster virial radii, the fraction of gas-rich galaxies is constant, representing the field. Within three cluster-centric radii, the fraction of gas-rich galaxies declines steadily with decreasing radius, reaching $\mathbf{<10%}$ near the cluster center. Our results suggest that galaxies start to feel the impact of the cluster environment on their gas content well beyond the cluster virial radius. We show that almost all gas-rich galaxies at the cluster virial radius are falling in for the first time at nearly radial orbits. Furthermore, we find that almost no galaxy moving outward at the cluster virial radius is gas-rich (with gas to baryon ratio greater than $\mathbf{1%}$). These results suggest that galaxies that fall into clusters lose their cold gas within a single radial round-trip.

How is mass distributed in the Universe? How does it compare with the distribution of light and stars? We address these questions by examining the distribution of mass, determined from weak lensing observations, and starlight, around $>10^5$ SDSS MaxBCG groups and clusters as a function of environment and scale, from deep inside clusters to large cosmic scales of $22 h^{-1}$ Mpc. The observed cumulative mass-to-light profile, $M/L (< r)$, rises on small scales, reflecting the increasing $M/L$ of the central bright galaxy of the cluster, then flattens to a nearly constant ratio on scales above $\sim 300 h^{-1}$ kpc, where light follows mass on all scales and in all environments. A trend of slightly decreasing $M/L (r)$ with scale is shown to be consistent with the varying stellar population following the morphology-density relation. This suggests that stars trace mass remarkably well even though they represent only a few percent of the total mass. We determine the stellar mass fraction and find it to be nearly constant on all scales above $\sim 300 h^{-1}$ kpc, with $M_{*}/M_{tot} \simeq 1.0\pm0.4\%$. We further suggest that most of the dark matter in the Universe is located in the large halos of individual galaxies ($\sim 300$ kpc for $L^{*}$ galaxies); we show that the entire $M/L (r)$ profile -- from groups and clusters to large-scale structure -- can be accounted for by the aggregate masses of the individual galaxies (whose halos may be stripped off but still remain in the clusters), plus gas. We use the observed mass-to-light ratio on large scales to determine the mass density of the Universe: $\Omega_{m} = 0.24 \pm 0.02 \times b_{M/L}^{2} = 0.26 \pm 0.02.$

Photometric large scale structure (LSS) surveys probe the largest volumes in the Universe, but are inevitably limited by systematic uncertainties. Imperfect photometric calibration leads to biases in our measurements of the density fields of LSS tracers such as galaxies and quasars, and as a result in cosmological parameter estimation. Earlier studies have proposed using cross-correlations between different redshift slices or cross-correlations between different surveys to reduce the effects of such systematics. In this paper we develop a method to characterize unknown systematics. We demonstrate that while we do not have sufficient information to correct for unknown systematics in the data, we can obtain an estimate of their magnitude. We define a parameter to estimate contamination from unknown systematics using cross-correlations between different redshift slices and propose discarding bins in the angular power spectrum that lie outside a certain contamination tolerance level. We show that this method improves estimates of the bias using simulated data and further apply it to photometric luminous red galaxies in the Sloan Digital Sky Survey as a case study.

We measure the luminosity and color dependence and the redshift evolution of galaxy clustering in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Ninth Data Release. We focus on the projected two-point correlation function (2PCF) of subsets of its CMASS sample, which includes about 260,000 galaxies over ~3,300 sq. deg in the redshift range 0.43<z<0.7. To minimize the selection effect on galaxy clustering, we construct well-defined luminosity and color subsamples by carefully accounting for the CMASS galaxy selection cuts. The 2PCF of the whole CMASS sample, if approximated by a power-law, has a correlation length of r_0=7.93\pm0.06Mpc/h and an index of \gamma=1.85\pm0.01. Clear dependences on galaxy luminosity and color are found for the projected 2PCF in all redshift bins, with more luminous and redder galaxies generally exhibiting stronger clustering and steeper 2PCF. The color dependence is also clearly seen for galaxies within the red sequence, consistent with the behavior of SDSS-II main sample galaxies at lower redshifts. At a given luminosity (k+e corrected), no significant evolution of the projected 2PCFs with redshift is detected for red sequence galaxies. We also construct galaxy samples of fixed number density at different redshifts, using redshift-dependent magnitude thresholds. The clustering of these galaxies in the CMASS redshift range is found to be consistent with that predicted by passive evolution. Our measurements of the luminosity and color dependence and redshift evolution of galaxy clustering will allow for detailed modeling of the relation between galaxies and dark matter halos and new constraints on galaxy formation and evolution.

We obtain constraints on cosmological parameters from the spherically averaged redshift-space correlation function of the CMASS Data Release 9 (DR9) sample of the Baryonic Oscillation Spectroscopic Survey (BOSS). We combine this information with additional data from recent CMB, SN and BAO measurements. Our results show no significant evidence of deviations from the standard flat-Lambda CDM model, whose basic parameters can be specified by Omega_m = 0.285 +- 0.009, 100 Omega_b = 4.59 +- 0.09, n_s = 0.96 +- 0.009, H_0 = 69.4 +- 0.8 km/s/Mpc and sigma_8 = 0.80 +- 0.02. The CMB+CMASS combination sets tight constraints on the curvature of the Universe, with Omega_k = -0.0043 +- 0.0049, and the tensor-to-scalar amplitude ratio, for which we find r < 0.16 at the 95 per cent confidence level (CL). These data show a clear signature of a deviation from scale-invariance also in the presence of tensor modes, with n_s <1 at the 99.7 per cent CL. We derive constraints on the fraction of massive neutrinos of f_nu < 0.049 (95 per cent CL), implying a limit of sum m_nu < 0.51 eV. We find no signature of a deviation from a cosmological constant from the combination of all datasets, with a constraint of w_DE = -1.033 +- 0.073 when this parameter is assumed time-independent, and no evidence of a departure from this value when it is allowed to evolve as w_DE(a) = w_0 + w_a (1 - a). The achieved accuracy on our cosmological constraints is a clear demonstration of the constraining power of current cosmological observations.

We present the final statistical sample of lensed quasars from the Sloan Digital Sky Survey (SDSS) Quasar Lens Search (SQLS). The well-defined statistical lens sample consists of 26 lensed quasars brighter than i=19.1 and in the redshift range of 0.6<z<2.2 selected from 50,836 spectroscopically confirmed quasars in the SDSS Data Release 7 (DR7), where we restrict the image separation range to 1"<\theta<20" and the i-band magnitude differences in two image lenses to be smaller than 1.25 mag. The SDSS DR7 quasar catalog also contains 36 additional lenses identified with various techniques. In addition to these lensed quasars, we have identified 81 pairs of quasars from follow-up spectroscopy, 26 of which are physically associated binary quasars. The statistical lens sample covers a wide range of image separations, redshifts, and magnitudes, and therefore is suitable for systematic studies of cosmological parameters and surveys of the structure and evolution of galaxies and quasars.

We present a statistical analysis of the final lens sample from the Sloan Digital Sky Survey Quasar Lens Search (SQLS). The number distribution of a complete subsample of 19 lensed quasars selected from 50,836 source quasars is compared with theoretical expectations, with particular attention to the selection function. Assuming that the velocity function of galaxies does not evolve with redshift, the SQLS sample constrains the cosmological constant to \Omega_\Lambda=0.79^{+0.06}_{-0.07}(stat.)^{+0.06}_{-0.06}(syst.) for a flat universe. The dark energy equation of state is found to be consistent with w=-1 when the SQLS is combined with constraints from baryon acoustic oscillation (BAO) measurements or results from the Wilkinson Microwave Anisotropy Probe (WMAP). We also obtain simultaneous constraints on cosmological parameters and redshift evolution of the galaxy velocity function, finding no evidence for redshift evolution at z<1 in any combinations of constraints. For instance, number density evolution quantified as \nu_n=d\ln\phi_*/d\ln(1+z) and the velocity dispersion evolution \nu_\sigma=d\ln\sigma_*/d\ln(1+z) are constrained to \nu_n=1.06^{+1.36}_{-1.39}(stat.)^{+0.33}_{-0.64}(syst.) and \nu_\sigma=-0.05^{+0.19}_{-0.16}(stat.)^{+0.03}_{-0.03}(syst.) respectively when the SQLS result is combined with BAO and WMAP for flat models with a cosmological constant. We find that a significant amount of dark energy is preferred even after fully marginalizing over the galaxy evolution parameters. Thus the statistics of lensed quasars robustly confirm the accelerated cosmic expansion.

We measure the large-scale intrinsic alignments of galaxy clusters in the Sloan Digital Sky Survey (SDSS) using subsets of two cluster catalogues: 6625 clusters with 0.1<z<0.3 from the maxBCG cluster catalogue (Koester et al. 2007, 7500 sq. deg.), and 8081 clusters with 0.08<z<0.44 from the Adaptive Matched Filter catalogue (Dong et al. 2008, 6500 sq. deg.). We search for two types of cluster alignments using pairs of clusters: the alignment between the projected major axes of the clusters (`correlation' alignment), and the alignment between one cluster major axis and the line connecting it to the other cluster in the pair (`pointing' alignment). In each case, we use the cluster member galaxy distribution as a tracer of the cluster shape. All measurements are carried out with each catalogue separately, to check for dependence on cluster selection procedure. We find a strong detection of the pointing alignment on scales up to 100 Mpc/h, at the 6 or 10-sigma level depending on the cluster selection algorithm used. The correlation alignment is only marginally detected up to ~20 Mpc/h, at the 2 or 2.5-sigma level. These results support our current theoretical understanding of galaxy cluster intrinsic alignments in the LCDM paradigm, although further work will be needed to understand the impact of cluster selection effects and observational measurement errors on the amplitude of the detection.

The Sloan Digital Sky Survey (SDSS) surveyed 14,555 square degrees, and delivered over a trillion pixels of imaging data. We present a study of galaxy clustering using 900,000 luminous galaxies with photometric redshifts, spanning between $z=0.45$ and $z=0.65$, constructed from the SDSS using methods described in Ross et al. (2011). This data-set spans 11,000 square degrees and probes a volume of $3h^{-3} \rm{Gpc}^3$, making it the largest volume ever used for galaxy clustering measurements. We present a novel treatment of the observational systematics and its applications to the clustering signals from the data set. In this paper, we measure the angular clustering using an optimal quadratic estimator at 4 redshift slices with an accuracy of ~15% with bin size of delta_l = 10 on scales of the Baryon Acoustic Oscillations (BAO) (at l~40-400). We derive cosmological constraints using the full-shape of the power-spectra. For a flat Lambda CDM model, when combined with Cosmic Microwave Background Wilkinson Microwave Anisotropy Probe 7 (WMAP7) and H_0 constraints from 600 Cepheids observed by HST, we find \Omega_\Lambda = 0.73 +/- 0.019 and H_0 to be 70.5 +/- 1.6 km/s/Mpc. For an open Lambda CDM model, when combined with WMAP7 + HST, we find $\Omega_K = 0.0035 +/- 0.0054, improved over WMAP7+HST alone by 40%. For a wCDM model, when combined with WMAP7+HST+SN, we find w = -1.071 +/- 0.078, and H_0 to be 71.3 +/- 1.7 km/s/Mpc, which is competitive with the latest large scale structure constraints from large spectroscopic surveys such as SDSS Data Release 7 (DR7) (Reid et al. 2010, Percival et al. 2010, Montesano et al. 2011) and WiggleZ (Blake et al. 2011). The SDSS-III Data Release 8 (SDSS-III DR8) Angular Clustering Data allows a wide range of investigations into the cosmological model, cosmic expansion (via BAO), Gaussianity of initial conditions and neutrino masses. (abridged)

We measure the luminosity and color dependence of galaxy clustering in the SDSS DR7 main galaxy sample, focusing on the projected correlation function w_p(r_p) of volume-limited samples. We interpret our measurements using halo occupation distribution (HOD) modeling assuming a Lambda-CDM cosmology. The amplitude of w_p(r_p) grows slowly with luminosity for L < L_* and increases sharply at higher luminosities, with bias factor b(>L)=1.06+0.23(L/L_*)^{1.12}. At fixed luminosity, redder galaxies have a stronger and steeper w_p(r_p), a trend that runs steadily from the bluest galaxies to the reddest galaxies. The individual luminosity trends for the red and blue galaxy populations are strikingly different. Blue galaxies show a slow but steady increase of w_p(r_p) with luminosity, at all scales. The large-scale clustering of red galaxies shows little luminosity dependence until a sharp increase at L > 4L_*, but the lowest luminosity red galaxies (0.04-0.25 L_*) show very strong clustering on scales r_p < 2 Mpc/h. Most of the observed trends can be naturally understood within the LCDM+HOD framework. The growth of w_p(r_p) with luminosity reflects an overall shift in the halo mass scale, in particular an increase in the minimum host halo mass Mmin. The mass at which a halo has, on average, one satellite galaxy brighter than L is M_1 ~ 17 Mmin(L) over most of the luminosity range. The growth and steepening of w_p(r_p) for redder galaxies reflects the increasing fraction of galaxies that are satellite systems in high mass halos instead of central systems in low mass halos, a trend that is especially marked at low luminosities. Our extensive measurements, provided in tabular form, will allow detailed tests of theoretical models of galaxy formation, a firm grounding of semi-empirical models of the galaxy population, and new cosmological tests.

Observations of clusters of galaxies suggest that they contain significantly fewer baryons (gas plus stars) than the cosmic baryon fraction. This `missing baryon' puzzle is especially surprising for the most massive clusters which are expected to be representative of the cosmic matter content of the universe (baryons and dark matter). Here we show that the baryons may not actually be missing from clusters, but rather are extended to larger radii than typically observed. The baryon deficiency is typically observed in the central regions of clusters (~0.5 the virial radius). However, the observed gas-density profile is significantly shallower than the mass-density profile, implying that the gas is more extended than the mass and that the gas fraction increases with radius. We use the observed density profiles of gas and mass in clusters to extrapolate the measured baryon fraction as a function of radius and as a function of cluster mass. We find that the baryon fraction reaches the cosmic value near the virial radius for all groups and clusters above 5e13 solar masses. This suggests that the baryons are not missing, they are simply located in cluster outskirts. Heating processes (shock-heating of the intracluster gas, plus supernovae and AGN feedback) that cause the gas to expand are likely explanations for these results. Upcoming observations should be able to detect these baryons.

We present the second report of our systematic search for strongly lensed quasars from the data of the Sloan Digital Sky Survey (SDSS). From extensive follow-up observations of 136 candidate objects, we find 36 lenses in the full sample of 77,429 spectroscopically confirmed quasars in the SDSS Data Release 5. We then define a complete sample of 19 lenses, including 11 from our previous search in the SDSS Data Release 3, from the sample of 36,287 quasars with i<19.1 in the redshift range 0.6<z<2.2, where we require the lenses to have image separations of 1"<\theta<20" and i-band magnitude differences between the two images smaller than 1.25 mag. Among the 19 lensed quasars, 3 have quadruple-image configurations, while the remaining 16 show double images. This lens sample constrains the cosmological constant to be \Omega_\Lambda=0.84^{+0.06}_{-0.08}(stat.)^{+0.09}_{-0.07}(syst.) assuming a flat universe, which is in good agreement with other cosmological observations. We also report the discoveries of 7 binary quasars with separations ranging from 1.1" to 16.6", which are identified in the course of our lens survey. This study concludes the construction of our statistical lens sample in the full SDSS-I data set.

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

We develop an improved mass tracer for clusters of galaxies from optically observed parameters, and calibrate the mass relation using weak gravitational lensing measurements. We employ a sample of ~ 13,000 optically-selected clusters from the Sloan Digital Sky Survey (SDSS) maxBCG catalog, with photometric redshifts in the range 0.1-0.3. The optical tracers we consider are cluster richness, cluster luminosity, luminosity of the brightest cluster galaxy (BCG), and combinations of these parameters. We measure the weak lensing signal around stacked clusters as a function of the various tracers, and use it to determine the tracer with the least amount of scatter. We further use the weak lensing data to calibrate the mass normalization. We find that the best mass estimator for massive clusters is a combination of cluster richness, N_{200}, and the luminosity of the brightest cluster galaxy, L_{BCG}: M_{200\bar{\rho}} = (1.27 +/- 0.08) (N_{200}/20)^{1.20 +/- 0.09} (L_{BCG}/\bar{L}_{BCG}(N_{200}))^{0.71 +/- 0.14} \times 10^{14} h^{-1} M_sun, where $\bar{L}_{BCG}(N_{200})$ is the observed mean BCG luminosity at a given richness. This improved mass tracer will enable the use of galaxy clusters as a more powerful tool for constraining cosmological parameters.

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

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

We examine the correlation function \xi of the Sloan Digital Sky Survey (SDSS) Luminous Red Galaxy sample (LRG) at large scales (60<s<400 Mpc/h) using the final data release (DR7; 105,831 LRGs between 0.16<z<0.47). Using mock catalogs, we demonstrate that the observed baryonic acoustic peak and larger scale signal are consistent with LCDM at the 1.5\sigma level. The signal at 155<s<200 Mpc/h tends to be high relative to theoretical expectations; this slight deviation can be attributed to a bright subsample of the LRGs. Fitting data to a non-linear, redshift-space, template based-model, we constrain the peak position at s_p=103.6+3.6-2.4 Mpc/h when fitting the range 60<s<150 Mpc/h (1\sigma uncertainties measured from the mocks. This redshift-space distance s_p is related to the comoving sound horizon scale r_s after taking into account matter clustering non-linearities, redshift distortions and galaxy clustering bias. Mock catalogs show that the probability that a DR7-sized sample would not have an identifiable peak is at least 10%. As a consistency check of a fiducial cosmology, we use the observed s_p to obtain the distance D_V=[(1+z)^2D_A^2cz/H(z)]^(1/3) relative to the acoustic scale. We find r_s/D_V(z=0.278)=0.1394+-0.0049. This result is in excellent agreement with Percival et. al (2009), who examine roughly the same data set, but using the power spectrum. Comparison with other determinations in the literature are also in very good agreement. We have tested our results against a battery of possible systematic effects, finding all effects are smaller than our estimated sample variance.

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

We analyze the cross-correlation of 2,705 unambiguously intervening Mg II (2796,2803A) quasar absorption line systems with 1,495,604 luminous red galaxies (LRGs) from the Fifth Data Release of the Sloan Digital Sky Survey within the redshift range 0.36<=z<=0.8. We confirm with high precision a previously reported weak anti-correlation of equivalent width and dark matter halo mass, measuring the average masses to be log M_h(M_[solar]h^-1)=11.29 [+0.36,-0.62] and log M_h(M_[solar]h^-1)=12.70 [+0.53,-1.16] for systems with W[2796A]>=1.4A and 0.8A<=W[2796A]<1.4A, respectively. Additionally, we investigate the significance of a number of potential sources of bias inherent in absorber-LRG cross-correlation measurements, including absorber velocity distributions and the weak lensing of background quasars, which we determine is capable of producing a 20-30% bias in angular cross-correlation measurements on scales less than 2'. We measure the Mg II - LRG cross-correlation for 719 absorption systems with v<60,000 km s^-1 in the quasar rest frame and find that these associated absorbers typically reside in dark matter haloes that are ~10-100 times more massive than those hosting unambiguously intervening Mg II absorbers. Furthermore, we find evidence for evolution of the redshift number density, dN/dz, with 2-sigma significance for the strongest (W>2.0A) absorbers in the DR5 sample. This width-dependent dN/dz evolution does not significantly affect the recovered equivalent width-halo mass anti-correlation and adds to existing evidence that the strongest Mg II absorption systems are correlated with an evolving population of field galaxies at z<0.8, while the non-evolving dN/dz of the weakest absorbers more closely resembles that of the LRG population.

We present measurements of the quasar two-point correlation function, \xi_{Q}, over the redshift range z=0.3-2.2 based upon data from the SDSS. Using a homogeneous sample of 30,239 quasars with spectroscopic redshifts from the DR5 Quasar Catalogue, our study represents the largest sample used for this type of investigation to date. With this redshift range and an areal coverage of approx 4,000 deg^2, we sample over 25 h^-3 Gpc^3 (comoving) assuming the current LCDM cosmology. Over this redshift range, we find that the redshift-space correlation function, xi(s), is adequately fit by a single power-law, with s_{0}=5.95+/-0.45 h^-1 Mpc and \gamma_{s}=1.16+0.11-0.16 when fit over s=1-25 h^-1 Mpc. Using the projected correlation function we calculate the real-space correlation length, r_{0}=5.45+0.35-0.45 h^-1 Mpc and \gamma=1.90+0.04-0.03, over scales of rp=1-130 h^-1 Mpc. Dividing the sample into redshift slices, we find very little, if any, evidence for the evolution of quasar clustering, with the redshift-space correlation length staying roughly constant at s_{0} ~ 6-7 h^-1 Mpc at z<2.2 (and only increasing at redshifts greater than this). Comparing our clustering measurements to those reported for X-ray selected AGN at z=0.5-1, we find reasonable agreement in some cases but significantly lower correlation lengths in others. We find that the linear bias evolves from b~1.4 at z=0.5 to b~3 at z=2.2, with b(z=1.27)=2.06+/-0.03 for the full sample. We compare our data to analytical models and infer that quasars inhabit dark matter haloes of constant mass M ~2 x 10^12 h^-1 M_Sol from redshifts z~2.5 (the peak of quasar activity) to z~0. [ABRIDGED]

We report the discovery of five gravitationally lensed quasars from the Sloan Digital Sky Survey (SDSS). All five systems are selected as two-image lensed quasar candidates from a sample of high-redshift (z>2.2) SDSS quasars. We confirmed their lensing nature with additional imaging and spectroscopic observations. The new systems are SDSS J0819+5356 (source redshift z_s=2.237, lens redshift z_l=0.294, and image separation \theta=4.04"), SDSS J1254+2235 (z_s=3.626, \theta=1.56"), SDSS J1258+1657 (z_s=2.702, \theta=1.28"), SDSS J1339+1310 (z_s=2.243, \theta=1.69"), and SDSS J1400+3134 (z_s=3.317, \theta=1.74"). We estimate the lens redshifts of the latter four systems to be z_l=0.2-0.8 from the colors and magnitudes of the lensing galaxies. We find that the image configurations of all systems are well reproduced by standard mass models. Although these lenses will not be included in our statistical sample of z_s<2.2 lenses, they expand the number of lensed quasars which can be used for high-redshift galaxy and quasar studies.

Using a homogenous sample of 38,208 quasars with a sky coverage of $4000 {\rm deg^2}$ drawn from the SDSS Data Release Five quasar catalog, we study the dependence of quasar clustering on luminosity, virial black hole mass, quasar color, and radio loudness. At $z<2.5$, quasar clustering depends weakly on luminosity and virial black hole mass, with typical uncertainty levels $\sim 10%$ for the measured correlation lengths. These weak dependences are consistent with models in which substantial scatter between quasar luminosity, virial black hole mass and the host dark matter halo mass has diluted any clustering difference, where halo mass is assumed to be the relevant quantity that best correlates with clustering strength. However, the most luminous and most massive quasars are more strongly clustered (at the $\sim 2\sigma$ level) than the remainder of the sample, which we attribute to the rapid increase of the bias factor at the high-mass end of host halos. We do not observe a strong dependence of clustering strength on quasar colors within our sample. On the other hand, radio-loud quasars are more strongly clustered than are radio-quiet quasars matched in redshift and optical luminosity (or virial black hole mass), consistent with local observations of radio galaxies and radio-loud type 2 AGN. Thus radio-loud quasars reside in more massive and denser environments in the biased halo clustering picture. Using the Sheth et al.(2001) formula for the linear halo bias, the estimated host halo mass for radio-loud quasars is $\sim 10^{13} h^{-1}M_\odot$, compared to $\sim 2\times 10^{12} h^{-1}M_\odot$ for radio-quiet quasar hosts at $z\sim 1.5$.

We cross-correlate large scale structure (LSS) observations from a number of surveys with CMB anisotropies from WMAP to investigate the Integrated Sachs-Wolfe (ISW) effect as a function of redshift, covering z~0.1-2.5. Our main goal is to go beyond reporting detections towards developing a reliable likelihood analysis that allows one to determine cosmological constraints from ISW observations. With this in mind we spend a considerable amount of effort in determining the redshift-dependent bias and redshift distribution b(z)*dN/dz of these samples by matching with spectroscopic observations where available, and analyzing auto-power spectra and cross-power spectra between the samples. The data sets we use are 2-Micron All Sky Survey (2MASS) samples, Sloan Digital Sky Survey (SDSS) photometric Luminous Red Galaxies, SDSS photometric quasars and NRAO VLA Sky Survey (NVSS) radio sources. We make a joint analysis of all samples constructing a full covariance matrix, which we subsequently use for cosmological parameter fitting. We report a 3.7 sigma detection of ISW combining all the datasets. We combine the ISW likelihood function with weak lensing of CMB and CMB power spectrum to constrain the equation of state of dark energy and the curvature of the Universe. While ISW does not significantly improve the constraints in the simplest 6-parameter flat Lambda CDM model, it improves constraints on 7-parameter models with curvature by a factor of 3.2 (relative to WMAP alone) to Omega_K=-0.004^{+0.014}_{-0.020}, and with dark energy equation of state by 15% to w=-1.01^{+0.30}_{-0.40}. (Abridged.)

We investigate the correlation of gravitational lensing of the cosmic microwave background (CMB) with several tracers of large-scale structure, including luminous red galaxies (LRGs), quasars, and radio sources. The lensing field is reconstructed based on the CMB maps from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite; the LRGs and quasars are observed by the Sloan Digital Sky Survey (SDSS); and the radio sources are observed in the NRAO VLA Sky Survey (NVSS). Combining all three large-scale structure samples, we find evidence for a positive cross-correlation at the $2.5\sigma$ level ($1.8\sigma$ for the SDSS samples and $2.1\sigma$ for NVSS); the cross-correlation amplitude is $1.06\pm 0.42$ times that expected for the WMAP cosmological parameters. Our analysis extends other recent analyses in that we carefully determine bias weighted redshift distribution of the sources, which is needed for a meaningful cosmological interpretation of the detected signal. We investigate contamination of the signal by Galactic emission, extragalactic radio and infrared sources, thermal and kinetic Sunyaev-Zel'dovich effects, and the Rees-Sciama effect, and find all of them to be negligible.

This is the first in a series of papers on the weak lensing effect caused by clusters of galaxies in Sloan Digital Sky Survey. The photometrically selected cluster sample, known as MaxBCG, includes ~130,000 objects between redshift 0.1 and 0.3, ranging in size from small groups to massive clusters. We split the clusters into bins of richness and luminosity and stack the surface density contrast to produce mean radial profiles. The mean profiles are detected over a range of scales, from the inner halo (25 kpc/h) well into the surrounding large scale structure (30 Mpc/h), with a significance of 15 to 20 in each bin. The signal over this large range of scales is best interpreted in terms of the cluster-mass cross-correlation function. We pay careful attention to sources of systematic error, correcting for them where possible. The resulting signals are calibrated to the ~10% level, with the dominant remaining uncertainty being the redshift distribution of the background sources. We find that the profiles scale strongly with richness and luminosity. We find the signal within a given richness bin depends upon luminosity, suggesting that luminosity is more closely correlated with mass than galaxy counts. We split the samples by redshift but detect no significant evolution. The profiles are not well described by power laws. In a subsequent series of papers we invert the profiles to three-dimensional mass profiles, show that they are well fit by a halo model description, measure mass-to-light ratios and provide a cosmological interpretation.

We select a sample of $\sim 4200$ traditionally defined broad absorption line quasars (BALQs) from the Fifth Data Release quasar catalog of the Sloan Digital Sky Survey. For a statistically homogeneous quasar sample with $1.7\le z\le 4.2$, the BAL quasar fraction is $\sim 14%$ and is almost constant with redshift. We measure the auto-correlation of non-BAL quasars (nonBALQs) and the cross-correlation of BALQs with nonBALQs using this statistically homogeneous sample, both in redshift space and using the projected correlation function. We find no significant difference between the clustering strengths of BALQs and nonBALQs. Assuming a power-law model for the real space correlation function $\xi(r)=(r/r_0)^{-1.8}$, the correlation length for nonBALQs is $r_0=7.6\pm 0.8 h^{-1}{\rm Mpc}$; for BALQs, the cross-correlation length is $r_0=7.4\pm 1.1 h^{-1}{\rm Mpc}$. Our clustering results suggest that BALQs live in similar large-scale environments as do nonBALQs.

We report the first results of our systematic search for strongly lensed quasars using the spectroscopically confirmed quasars in the Sloan Digital Sky Survey (SDSS). Among 46,420 quasars from the SDSS Data Release 3 (~4188 deg^2), we select a subsample of 22,683 quasars that are located at redshifts between 0.6 and 2.2 and are brighter than the Galactic extinction corrected i-band magnitude of 19.1. We identify 220 lens candidates from the quasar subsample, for which we conduct extensive and systematic follow-up observations in optical and near-infrared wavebands, in order to construct a complete lensed quasar sample at image separations between 1'' and 20'' and flux ratios of faint to bright lensed images larger than 10^{-0.5}. We construct a statistical sample of 11 lensed quasars. Ten of these are galaxy-scale lenses with small image separations (~1''-2'') and one is a large separation (15'') system which is produced by a massive cluster of galaxies, representing the first statistical sample of lensed quasars including both galaxy- and cluster-scale lenses. The Data Release 3 spectroscopic quasars contain an additional 11 lensed quasars outside the statistical sample.

We present cosmological results from the statistics of lensed quasars in the Sloan Digital Sky Survey (SDSS) Quasar Lens Search. By taking proper account of the selection function, we compute the expected number of quasars lensed by early-type galaxies and their image separation distribution assuming a flat universe, which is then compared with 7 lenses found in the SDSS Data Release 3 to derive constraints on dark energy under strictly controlled criteria. For a cosmological constant model (w=-1) we obtain \Omega_\Lambda=0.74^{+0.11}_{-0.15}(stat.)^{+0.13}_{-0.06}(syst.). Allowing w to be a free parameter we find \Omega_M=0.26^{+0.07}_{-0.06}(stat.)^{+0.03}_{-0.05}(syst.) and w=-1.1\pm0.6(stat.)^{+0.3}_{-0.5}(syst.) when combined with the constraint from the measurement of baryon acoustic oscillations in the SDSS luminous red galaxy sample. Our results are in good agreement with earlier lensing constraints obtained using radio lenses, and provide additional confirmation of the presence of dark energy consistent with a cosmological constant, derived independently of type Ia supernovae.

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

We present a catalog of galaxy clusters selected using the maxBCG redsequence method from Sloan Digital Sky Survey photometric data. This catalog includes 13,823 clusters with velocity dispersions greater than 400 km/s, and is the largest galaxy cluster catalog assembled to date. They are selected in an approximately volume-limited way from a 0.5 Gpc^3 region covering 7500 square degrees of sky between redshifts 0.1 and 0.3. (ABRIGDED)

(Abridged) We study the two-point correlation function of a uniformly selected sample of 4,426 luminous optical quasars with redshift $2.9 \le z\le 5.4$ selected over 4041 deg$^2$ from the Fifth Data Release of the Sloan Digital Sky Survey. For a real-space correlation function of the form $\xi(r)=(r/r_0)^{-\gamma}$, the fitted parameters in comoving coordinates are $r_0 = 15.2 \pm 2.7 h^{-1}$ Mpc and $\gamma = 2.0 \pm 0.3$, over a scale range $4\le r_p\le 150 h^{-1}$ Mpc. Thus high-redshift quasars are appreciably more strongly clustered than their $z \approx 1.5$ counterparts, which have a comoving clustering length $r_0 \approx 6.5 h^{-1}$ Mpc. Dividing our sample into two redshift bins: $2.9\le z\le 3.5$ and $z\ge 3.5$, and assuming a power-law index $\gamma=2.0$, we find a correlation length of $r_0 = 16.9 \pm 1.7 h^{-1}$ Mpc for the former, and $r_0 = 24.3 \pm 2.4 h^{-1}$ Mpc for the latter. Following Martini & Weinberg, we relate the clustering strength and quasar number density to the quasar lifetimes and duty cycle. Using the Sheth & Tormen halo mass function, the quasar lifetime is estimated to lie in the range $4\sim 50$ Myr for quasars with $2.9\le z\le 3.5$; and $30\sim 600$ Myr for quasars with $z\ge 3.5$. The corresponding duty cycles are $0.004\sim 0.05$ for the lower redshift bin and $0.03\sim 0.6$ for the higher redshift bin. The minimum mass of halos in which these quasars reside is $2-3\times 10^{12}\ h^{-1}M_\odot$ for quasars with $2.9\le z\le 3.5$ and $4-6\times 10^{12}\ h^{-1}M_\odot$ for quasars with $z\ge 3.5$.

We determine the shape, multiplicity, size, and radial structure of superclusters in the LambdaCDM concordance cosmology from z = 0 to z = 2. Superclusters are defined as clusters of clusters in our large-scale cosmological simulation. We find that superclusters are triaxial in shape; many have flattened since early times to become nearly two-dimensional structures at present, with a small fraction of filamentary systems. The size and multiplicity functions are presented at different redshifts. Supercluster sizes extend to scales of ~ 100 - 200 Mpc/h. The supercluster multiplicity (richness) increases linearly with supercluster size. The density profile in superclusters is approximately isothermal (~ R^{-2}) and steepens on larger scales. These results can be used as a new test of the current cosmology when compared with upcoming observations of large-scale surveys.

Using ~300,000 photometrically classified quasars, by far the largest quasar sample ever used for such analyses, we study the redshift and luminosity evolution of quasar clustering on scales of ~50 kpc/h to ~20 Mpc/h from redshifts of z~0.75 to z~2.28. We parameterize our clustering amplitudes using realistic dark matter models, and find that a LCDM power spectrum provides a superb fit to our data with a redshift-averaged quasar bias of b_Q = 2.41+/-0.08 ($P_{<\chi^2}=0.847$) for $\sigma_8=0.9$. This represents a better fit than the best-fit power-law model ($\omega = 0.0493\pm0.0064\theta^ {-0.928\pm0.055}$; $P_{<\chi^2}=0.482$). We find b_Q increases with redshift. This evolution is significant at >99.6% using our data set alone, increasing to >99.9999% if stellar contamination is not explicitly parameterized. We measure the quasar classification efficiency across our full sample as a = 95.6 +/- ^{4.4}_{1.9}%, a star-quasar separation comparable with the star-galaxy separation in many photometric studies of galaxy clustering. We derive the mean mass of the dark matter halos hosting quasars as MDMH=(5.2+/-0.6)x10^{12} M_solar/h. At z~1.9 we find a $1.5\sigma$ deviation from luminosity-independent quasar clustering; this suggests that increasing our sample size by a factor of 1.8 could begin to constrain any luminosity dependence in quasar bias at z~2. Our results agree with recent studies of quasar environments at z < 0.4, which detected little luminosity dependence to quasar clustering on proper scales >50 kpc/h. At z < 1.6, our analysis suggests that b_Q is constant with luminosity to within ~0.6, and that, for g < 21, angular quasar autocorrelation measurements are unlikely to have sufficient statistical power at z < 1.6 to detect any luminosity dependence in quasars' clustering.

We study quasar clustering on small scales, modeling clustering amplitudes using halo-driven dark matter descriptions. From 91 pairs on scales <35 kpc/h, we detect only a slight excess in quasar clustering over our best-fit large-scale model. Integrated across all redshifts, the implied quasar bias is b_Q = 4.21+/-0.98 (b_Q = 3.93+/-0.71) at ~18 kpc/h (~28 kpc/h). Our best-fit (real-space) power index is ~-2 (i.e., $\xi(r) \propto r^{-2}$), implying steeper halo profiles than currently found in simulations. Alternatively, quasar binaries with separation <35 kpc/h may trace merging galaxies, with typical dynamical merger times t_d~(610+/-260)m^{-1/2} Myr/h, for quasars of host halo mass m x 10^{12} Msolar/h. We find UVX quasars at ~28 kpc/h cluster >5 times higher at z > 2, than at z < 2, at the $2.0\sigma$ level. However, as the space density of quasars declines as z increases, an excess of quasar binaries (over expectation) at z > 2 could be consistent with reduced merger rates at z > 2 for the galaxies forming UVX quasars. Comparing our clustering at ~28 kpc/h to a $\xi(r)=(r/4.8\Mpch)^{-1.53}$ power-law, we find an upper limit on any excess of a factor of 4.3+/-1.3, which, noting some caveats, differs from large excesses recently measured for binary quasars, at $2.2\sigma$. We speculate that binary quasar surveys that are biased to z > 2 may find inflated clustering excesses when compared to models fit at z < 2. We provide details of 111 photometrically classified quasar pairs with separations <0.1'. Spectroscopy of these pairs could significantly constrain quasar dynamics in merging galaxies.

We measure the cosmological matter density by observing the positions of baryon acoustic oscillations in the clustering of galaxies in the Sloan Digital Sky Survey (SDSS). We jointly analyse the main galaxies and LRGs in the SDSS DR5 sample, using over half a million galaxies in total. The oscillations are detected with 99.74% confidence (3.0sigma assuming Gaussianity) compared to a smooth power spectrum. When combined with the observed scale of the peaks within the CMB, we find a best-fit value of Omega_m=0.256+0.029-0.024 (68% confidence interval), for a flat Lambda cosmology when marginalising over the Hubble parameter and the baryon density. This value of the matter density is derived from the locations of the baryon oscillations in the galaxy power spectrum and in the CMB, and does not include any information from the overall shape of the power spectra. This is an extremely clean cosmological measurement as the physics of the baryon acoustic oscillation production is well understood, and the positions of the oscillations are expected to be independent of systematics such as galaxy bias.

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

We report on the first results of an imaging survey to detect strong gravitational lensing targeting the richest clusters selected from the photometric data of the Sloan Digital Sky Survey (SDSS) with follow-up deep imaging observations from the Wisconsin Indiana Yale NOAO (WIYN) 3.5m telescope and the University of Hawaii 88-inch telescope (UH88). The clusters are selected from an area of 8000 deg^2 using the Red Cluster Sequence technique and span the redshift range 0.1 < z < 0.6, corresponding to a comoving cosmological volume of ~ 2 Gpc^3. Our imaging survey thus targets a volume more than an order of magnitude larger than any previous search. A total of 240 clusters were imaged of which 141 had sub-arcsecond image quality. Our survey has uncovered16 new lensing clusters with definite giant arcs, an additional 12 systems for which the lensing interpretation is very likely, and 9 possible lenses which contain shorter arclets or candidate arcs which are less certain and will require further observations to confirm their lensing origin. The number of new cluster lenses detected in this survey is likely > 30. Among these new systems are several of the most dramatic examples of strong gravitational lensing ever discovered with multiple bright arcs at large angular separation. These will likely become 'poster-child' gravitational lenses similar to Abell 1689 and CL0024+1654. The new lenses discovered in this survey will enable future sysetmatic studies of the statistics of strong lensing and its implications for cosmology and our structure formation paradigm.

We present a Fourier analysis of the clustering of galaxies in the combined Main galaxy and Luminous Red Galaxy (LRG) Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5) sample. The aim of our analysis is to consider how well we can measure the cosmological matter density using the signature of the horizon at matter-radiation equality embedded in the large-scale power spectrum. The new data constrains the power spectrum on scales 100--600h^-1Mpc with significantly higher precision than previous analyses of just the SDSS Main galaxies, due to our larger sample and the inclusion of the LRGs. This improvement means that we can now reveal a discrepancy between the shape of the measured power and linear CDM models on scales 0.01<k<0.15hMpc^-1, with linear model fits favouring a lower matter density (Omega_m=0.22+/-0.04) on scales 0.01<k<0.06hMpc^-1 and a higher matter density (Omega_m=0.32+/-0.01) when smaller scales are included, assuming a flat LCDM model with h=0.73 and n_s=0.96. This discrepancy could be explained by scale-dependent bias and, by analysing subsamples of galaxies, we find that the ratio of small-scale to large-scale power increases with galaxy luminosity, so all of the SDSS galaxies cannot trace the same power spectrum shape over 0.01<k<0.2hMpc^-1. However, the data are insufficient to clearly show a luminosity-dependent change in the largest scale at which a significant increase in clustering is observed, although they do not rule out such an effect. Significant scale-dependent galaxy bias on large-scales, which changes with the r-band luminosity of the galaxies, could potentially explain differences in our Omega_m estimates and differences previously observed between 2dFGRS and SDSS power spectra and the resulting parameter constraints.

We identify galaxy groups and clusters in volume-limited samples of the SDSS redshift survey, using a redshift-space friends-of-friends algorithm. We optimize the friends-of-friends linking lengths to recover galaxy systems that occupy the same dark matter halos, using a set of mock catalogs created by populating halos of N-body simulations with galaxies. Extensive tests with these mock catalogs show that no combination of perpendicular and line-of-sight linking lengths is able to yield groups and clusters that simultaneously recover the true halo multiplicity function, projected size distribution, and velocity dispersion. We adopt a linking length combination that yields, for galaxy groups with ten or more members: a group multiplicity function that is unbiased with respect to the true halo multiplicity function; an unbiased median relation between the multiplicities of groups and their associated halos; a spurious group fraction of less than ~1%; a halo completeness of more than ~97%; the correct projected size distribution as a function of multiplicity; and a velocity dispersion distribution that is ~20% too low at all multiplicities. These results hold over a range of mock catalogs that use different input recipes of populating halos with galaxies. We apply our group-finding algorithm to the SDSS data and obtain three group and cluster catalogs for three volume-limited samples that cover 3495.1 square degrees on the sky. We correct for incompleteness caused by fiber collisions and survey edges, and obtain measurements of the group multiplicity function, with errors calculated from realistic mock catalogs. These multiplicity function measurements provide a key constraint on the relation between galaxy populations and dark matter halos.

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

We present measurements of the normalised redshift-space three-point correlation function (Q_z) of galaxies from the Sloan Digital Sky Survey (SDSS) main galaxy sample. We have applied our "npt" algorithm to both a volume-limited (36738 galaxies) and magnitude-limited sample (134741 galaxies) of SDSS galaxies, and find consistent results between the two samples, thus confirming the weak luminosity dependence of Q_z recently seen by other authors. We compare our results to other Q_z measurements in the literature and find it to be consistent within the full jack-knife error estimates. However, we find these errors are significantly increased by the presence of the ``Sloan Great Wall'' (at z ~ 0.08) within these two SDSS datasets, which changes the 3-point correlation function (3PCF) by 70% on large scales (s>=10h^-1 Mpc). If we exclude this supercluster, our observed Q_z is in better agreement with that obtained from the 2dFGRS by other authors, thus demonstrating the sensitivity of these higher-order correlation functions to large-scale structures in the Universe. This analysis highlights that the SDSS datasets used here are not ``fair samples'' of the Universe for the estimation of higher-order clustering statistics and larger volumes are required. We study the shape-dependence of Q_z(s,q,theta) as one expects this measurement to depend on scale if the large scale structure in the Universe has grown via gravitational instability from Gaussian initial conditions. On small scales (s <= 6h^-1 Mpc), we see some evidence for shape-dependence in Q_z, but at present our measurements are consistent with a constant within the errors (Q_z ~ 0.75 +/- 0.05). On scales >10h^-1 Mpc, we see considerable shape-dependence in Q_z.

Where do quasars reside? Are quasars located in environments similar to those of typical L* galaxies, and, if not, how do they differ? An answer to this question will help shed light on the triggering process of quasar activity. We use the Sloan Digital Sky Survey to study the environment of quasars and compare it directly with the environment of galaxies. We find that quasars (M_i < -22, z < 0.4) are located in higher local overdensity regions than are typical L* galaxies. The enhanced environment around quasars is a local phenomenon; the overdensity relative to that around L* galaxies is strongest within 100 kpc of the quasars. In this region, the overdensity is a factor of 1.4 larger than around L* galaxies. The overdensity declines monotonically with scale to nearly unity at ~1 Mpc, where quasars inhabit environments comparable to those of L* galaxies. The small-scale density enhancement depends on quasar luminosity, but only at the brightest end: the most luminous quasars reside in higher local overdensity regions than do fainter quasars. The mean overdensity around the brightest quasars (M_i < -23.3) is nearly three times larger than around L* galaxies while the density around dimmer quasars (M_i = -22.0 to -23.3) is ~1.4 times that of L* galaxies. By ~0.5 Mpc, the dependence on quasar luminosity is no longer significant. The overdensity on all scales is independent of redshift to z = 0.4. The results suggest a picture in which quasars typically reside in L* galaxies, but have a local excess of neighbors within ~0.1 - 0.5 Mpc; this local density excess likely contributes to the triggering of quasar activity through mergers and other interactions.

We present the discovery of seven quasars at z>5.7, selected from ~2000 deg^2 of multicolor imaging data of the Sloan Digital Sky Survey (SDSS). The new quasars have redshifts z from 5.79 to 6.13. Five are selected as part of a complete flux-limited sample in the SDSS Northern Galactic Cap; two have larger photometric errors and are not part of the complete sample. One of the new quasars, SDSS J1335+3533 (z=5.93), exhibits no emission lines; the 3-sigma limit on the rest-frame equivalent width of Ly alpha+NV line is 5 A. It is the highest redshift lineless quasar known, and could be a gravitational lensed galaxy, a BL Lac object or a new type of quasar. Two new z>6 quasars, SDSS 1250+3130 (z=6.13) and SDSS J1137+3549 (z=6.01), show deep Gunn-Peterson absorption gaps in Ly alpha. These gaps are narrower the complete Gunn-Peterson absorption troughs observed among quasars at z>6.2 and do not have complete Ly beta absorption.

We investigate the dependence of ellipticities of clusters of galaxies on cosmological parameters using large-scale cosmological simulations. We determine cluster ellipticities out to redshift unity for LCDM models with different mean densities $\Omega_m$ and amplitudes of mass fluctuation $\sigma_{8,0}$. The mean ellipticity increases monotonically with redshift for all models. Larger values of $\sigma_{8,0}$, i.e., earlier cluster formation time, produce lower ellipticities. The dependence of ellipticity on $\Omega_m$ is relatively weak in the range $0.2 \leq \Omega_m \leq 0.5$ for high mass clusters. The mean ellipticity $\bar{e}(z)$ decreases linearly with the amplitude of fluctuations at the cluster redshift $z$, nearly independent of $\Omega_m$; on average, older clusters are more relaxed and are thus less elliptical. The distribution of ellipticities about the mean is approximated by a Gaussian, allowing a simple characterization of the evolution of ellipticity with redshift as a function of cosmological parameters. At $z=0$, the mean ellipticity of high mass clusters is approximated by $\bar{e}(z=0) = 0.248-0.069 \sigma_{8,0} + 0.013 \Omega_{m,0}$. This relation opens up the possibility that, when compared with future observations of large cluster samples, the mean cluster ellipticity and its evolution could be used as a new, independent tool to constrain cosmological parameters, especially the amplitude of mass fluctuations, $\sigma_{8,0}$.

We have produced a new conformal map of the universe illustrating recent discoveries, ranging from Kuiper belt objects in the Solar system, to the galaxies and quasars from the Sloan Digital Sky Survey. This map projection, based on the logarithm map of the complex plane, preserves shapes locally, and yet is able to display the entire range of astronomical scales from the Earth's neighborhood to the cosmic microwave background. The conformal nature of the projection, preserving shapes locally, may be of particular use for analyzing large scale structure. Prominent in the map is a Sloan Great Wall of galaxies 1.37 billion light years long, 80% longer than the Great Wall discovered by Geller and Huchra and therefore the largest observed structure in the universe.

We report on i-band snapshot observations of 157 Sloan Digital Sky Survey (SDSS) quasars at 4<z<5.4 using the Advanced Camera for Surveys on the Hubble Space Telescope (HST) to search for evidence of gravitational lensing of these sources. None of the quasars appear to be strongly lensed and multiply imaged at the angular resolution (~0.1") and sensitivity of HST. The non-detection of strong lensing in these systems constrains the z=4-5 luminosity function to an intrinsic slope of beta>-3.8 (3 sigma), assuming a break in the quasar luminosity function at M_1450^star=-24.5. This constraint is considerably stronger than the limit of beta>-4.63 obtained from the absence of lensing in four z>5.7 quasars. Such constraints are important for our understanding of the true space density of high-redshift quasars and the ionization state of the early universe.

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

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

The amount and nature of dark energy (DE) can be tightly constrained by measuring the spatial correlation features and evolution of a sample of ~ 100,000 galaxy clusters over the redshift range 0<z < 1.5. Such an X-ray survey will discover all collapsed structures with mass above 3.5e14 solar masss at redshifts z<2 (i.e. the full range where such objects are expected) in the high Galactic latitude sky. Above this mass threshold the tight correlations between X-ray observables and mass allow direct interpretation of the data. We describe the constraints on Dark Energy that can be inferred from such a survey, using powerful self-calibration techniques to relate the X-ray observables (luminosity and temperature) to the underlying mass.

We report the first result of the clustering analysis of Sloan Digital Sky Survey (SDSS) quasars. We compute the two-point correlation function (2PCF) of SDSS quasars in redshift space at $8h^{-1}{\rm Mpc} < s < 500h^{-1}{\rm Mpc}$, with particular attention to its baryonic signature. Our sample consists of 19986 quasars extracted from the SDSS Data Release 4 (DR4). The redshift range of the sample is $0.72 \le z \le 2.24$ (the mean redshift is $\bar z = 1.46$) and the reddening-corrected $i$-band apparent magnitude range is $15.0 \le m_{i,{\rm rc}} \le 19.1$. Due to the relatively low number density of the quasar sample, the bump in the power spectrum due to the baryon density, $\Omega_{\rm b}$, is not clearly visible. The effect of the baryon density is, however, to distort the overall shape of the 2PCF.The degree of distortion makes it an interesting alternate measure of the baryonic signature. Assuming a scale-independent linear bias and the spatially flat universe, i.e., $\Omega_{\rm b} + \Omega_{\rm d} + \Omega_\Lambda =1$, where $\Omega_{\rm d}$ and $\Omega_\Lambda$ denote the density parameters of dark matter and the cosmological constant, we combine the observed quasar 2PCF and the predicted matter 2PCF to put constraints on $\Omega_{\rm b}$ and $\Omega_\Lambda$. Our result is fitted as $0.80- 2.8\Omega_{\rm b} < \Omega_\Lambda < 0.90 - 1.4\Omega_{\rm b}$ at the 2$\sigma$ confidence level, which is consistent with results from other cosmological observations such as WMAP. (abridged)

We present a sample of 218 new quasar pairs with proper transverse separations R_prop < 1 Mpc/h over the redshift range 0.5 < z < 3.0, discovered from an extensive follow up campaign to find companions around the Sloan Digital Sky Survey and 2dF Quasar Redshift Survey quasars. This sample includes 26 new binary quasars with separations R_prop < 50 kpc/h (theta < 10 arcseconds), more than doubling the number of such systems known. We define a statistical sample of binaries selected with homogeneous criteria and compute its selection function, taking into account sources of incompleteness. The first measurement of the quasar correlation function on scales 10 kpc/h < R_prop < 400 kpc/h is presented. For R_prop < 40 kpc/h, we detect an order of magnitude excess clustering over the expectation from the large scale R_prop > 3 Mpc/h quasar correlation function, extrapolated down as a power law to the separations probed by our binaries. The excess grows to ~ 30 at R_prop ~ 10 kpc/h, and provides compelling evidence that the quasar autocorrelation function gets progressively steeper on sub-Mpc scales. This small scale excess can likely be attributed to dissipative interaction events which trigger quasar activity in rich environments. Recent small scale measurements of galaxy clustering and quasar-galaxy clustering are reviewed and discussed in relation to our measurement of small scale quasar clustering.

We present an 8 sigma detection of cosmic magnification measured by the variation of quasar density due to gravitational lensing by foreground large scale structure. To make this measurement we used 3800 square degrees of photometric observations from the Sloan Digital Sky Survey (SDSS) containing \~200,000 quasars and 13 million galaxies. Our measurement of the galaxy-quasar cross-correlation function exhibits the amplitude, angular dependence and change in sign as a function of the slope of the observed quasar number counts that is expected from magnification bias due to weak gravitational lensing. We show that observational uncertainties (stellar contamination, Galactic dust extinction, seeing variations and errors in the photometric redshifts) are well controlled and do not significantly affect the lensing signal. By weighting the quasars with the number count slope, we combine the cross-correlation of quasars for our full magnitude range and detect the lensing signal at >4 sigma in all five SDSS filters. Our measurements of cosmic magnification probe scales ranging from 60 kpc/h to 10 Mpc/h and are in good agreement with theoretical predictions based on the WMAP concordance cosmology. As with galaxy-galaxy lensing, future measurements of cosmic magnification will provide useful constraints on the galaxy-mass power spectrum.

We present the "C4 Cluster Catalog", a new sample of 748 clusters of galaxies identified in the spectroscopic sample of the Second Data Release (DR2) of the Sloan Digital Sky Survey (SDSS). The C4 cluster--finding algorithm identifies clusters as overdensities in a seven-dimensional position and color space, thus minimizing projection effects which plagued previous optical clusters selection. The present C4 catalog covers ~2600 square degrees of sky with groups containing 10 members to massive clusters having over 200 cluster members with redshifts. We provide cluster properties like sky location, mean redshift, galaxy membership, summed r--band optical luminosity (L_r), velocity dispersion, and measures of substructure. We use new mock galaxy catalogs to investigate the sensitivity to the various algorithm parameters, as well as to quantify purity and completeness. These mock catalogs indicate that the C4 catalog is ~90% complete and 95% pure above M_200 = 1x10^14 solar masses and within 0.03 <=z <= 0.12. The C4 algorithm finds 98% of X-ray identified clusters and 90% of Abell clusters within 0.03 <= z <= 0.12. We show that the L_r of a cluster is a more robust estimator of the halo mass (M_200) than the line-of-sight velocity dispersion or the richness of the cluster. L_r. The final SDSS data will provide ~2500 C4 clusters and will represent one of the largest and most homogeneous samples of local clusters.

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.

We combine the measurements of luminosity dependence of bias with the luminosity dependent weak lensing analysis of dark matter around galaxies to derive the galaxy bias and constrain nonlinear mass and cosmological parameters. We take advantage of theoretical and simulation predictions that predict that while halo bias is rapidly increasing with mass for high masses, it is nearly constant in low mass halos. We use a new weak lensing analysis around the same SDSS galaxies to determine their halo mass probability distribution. These halo mass probability distributions are used to predict the bias for each luminosity subsample and we find an excellent agreement with observed values. The required nonlinear mass suggests slightly lower matter density than usually assumed, Omegam=0.25+/- 0.03 for the simplest models. We combine the bias constraints with those from the WMAP and the SDSS power spectrum analysis to derive new constraints on bias and sigma_8. For the most general parameter space we find sigma_8=0.88+/- 0.06 and b_*=0.99+/- 0.07. In the context of spatially flat models we improve the limit on the neutrino mass for the case of 3 degenerate families from m_nu<0.6eV without bias to m_nu<0.18eV with bias (95% c.l.), which is weakened to m_nu<0.24eV if running is allowed. The corresponding limit for 3 massless + 1 massive neutrino is 1.37eV.

We use the projected correlation function w_p(r_p) of a volume-limited subsample of the Sloan Digital Sky Survey (SDSS) main galaxy redshift catalogue to measure the halo occupation distribution (HOD) of the galaxies of the sample. Simultaneously, we allow the cosmology to vary within cosmological constraints imposed by cosmic microwave background experiments in a Lambda-CDM model. We find that combining w_p(r_p) for this sample alone with the observations by WMAP, ACBAR, CBI and VSA can provide one of the most precise techniques available to measure cosmological parameters. For a minimal flat six-parameter Lambda-CDM model with an HOD with three free parameters, we find Omega_m=0.278^{+0.027}_{-0.027}, sigma_8=0.812^{+0.028}_{-0.027}, and H_0=69.8^{+2.6}_{-2.6}km s^{-1} Mpc^{-1}; these errors are significantly smaller than from CMB alone and similar to those obtained by combining CMB with the large-scale galaxy power spectrum assuming scale-independent bias. The corresponding HOD parameters describing the minimum halo mass and the normalization and cut-off of the satellite mean occupation are M_min=(3.03^{+0.36}_{-0.36})x 10^{12} h^{-1} M_sun, M_1 = (4.58^{+0.60}_{-0.60})x 10^{13} h^{-1} M_sun, and kappa=4.44^{+0.51}_{-0.69}. When more parameters are added to the HOD model, the error bars on the HOD parameters increase because of degeneracies, but the error bars on the cosmological parameters do not increase greatly. Similar modeling for other galaxy samples could reduce the statistical errors on these results, while more thorough investigations of the cosmology dependence of nonlinear halo bias and halo mass functions are needed to eliminate remaining systematic uncertainties, which may be comparable to statistical uncertainties.

The evolution of the cluster mass function and the cluster correlation function from z = 0 to z = 3 are determined using 10^6 clusters obtained from high-resolution simulations of the current best-fit LCDM cosmology (\Omega_m = 0.27, \sigma_8 = 0.84, h = 0.7). The results provide predictions for comparisons with future observations of high redshift clusters. A comparison of the predicted mass function of low redshift clusters with observations from early Sloan Digital Sky Survey data, and the predicted abundance of massive distant clusters with observational results, favor a slightly larger amplitude of mass fluctuations (\sigma_8 = 0.9) and lower density parameter (\Omega_m = 0.2); these values are consistent within 1-\sigma with the current observational and model uncertainties. The cluster correlation function strength increases with redshift for a given mass limit; the clusters were more strongly correlated in the past, due to their increasing bias with redshift - the bias reaches b = 100 at z = 2 for M > 5 x 10^13 h^-1 M_sun. The richness-dependent cluster correlation function, represented by the correlation scale versus cluster mean separation relation, R0-d, is generally consistent with observations. This relation can be approximated as R_0 = 1.7 d^0.6 h^-1 Mpc for d = 20 - 60 h^-1 Mpc. The R0-d relation exhibits surprisingly little evolution with redshift for z < 2; this can provide a new test of the current LCDM model when compared with future observations of high redshift clusters.

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 construct and analyze a u-band selected galaxy sample from the SDSS Southern Survey, which covers 275 sq. deg. The sample includes 43223 galaxies with spectroscopic redshifts in the range 0.005<z<0.3 and with 14.5<u<20.5. The S/N in the u-band Petrosian aperture is improved by coadding multiple epochs of imaging data and by including sky-subtraction corrections. Luminosity functions for the near-UV ^{0.1}u band (322+-26 nm) are determined in redshift slices of width 0.02, which show a highly significant evolution in M^{star} of -0.8+-0.1 mag between z=0 and 0.3; with M^{star} = -18.84+-0.05 (AB mag), log phi^{star} = -2.06+-0.03 (Mpc^{-3}) and log rho_L = 19.11+-0.02 (W Hz^{-1} Mpc^{-3}) at z=0.1. The faint-end slope determined for z<0.06 is given by alpha=-1.05+-0.08. This is in agreement with recent determinations from GALEX at shorter wavelengths. Comparing our z<0.3 luminosity density measurements with 0.2<z<1.2 from COMBO-17, we find that the 280-nm density evolves as rho_L proportional to (1+z)^{beta} with beta=2.1+-0.2; and find no evidence for any change in slope over this redshift range. By comparing with other measurements of cosmic star formation history, we estimate that the effective dust attenuation at 280 nm has increased by 0.8+-0.3 mag between z=0 and 1.

We report the discovery of two doubly-imaged quasars, SDSS J100128.61+502756.9 and SDSS J120629.65+433217.6, at redshifts of 1.838 and 1.789 and with image separations of 2.86'' and 2.90'', respectively. The objects were selected as lens candidates from the Sloan Digital Sky Survey (SDSS). Based on the identical nature of the spectra of the two quasars in each pair and the identification of the lens galaxies, we conclude that the objects are gravitational lenses. The lenses are complicated; in both systems there are several galaxies in the fields very close to the quasars, in addition to the lens galaxies themselves. The lens modeling implies that these nearby galaxies contribute significantly to the lens potentials. On larger scales, we have detected an enhancement in the galaxy density near SDSS J100128.61+502756.9. The number of lenses with image separation of ~3'' in the SDSS already exceeds the prediction of simple theoretical models based on the standard Lambda-dominated cosmology and observed velocity function of galaxies.

We present the small-scale (0.2h^{-1} to 7h^{-1} Mpc) cross-correlations between 32,000 luminous early-type galaxies and a reference sample of 16 million normal galaxies from the Sloan Digital Sky Survey. Our method allows us to construct the spherically averaged, real-space cross-correlation function between the spectroscopic LRG sample and galaxies from the SDSS imaging. We report the cross-correlation as a function of scale, luminosity, and redshift. We find very strong luminosity dependences in the clustering amplitudes, up to a factor of 4 over a factor of 4 in luminosity, and measure this dependence with high signal-to-noise ratio. The luminosity dependence of bias is found to depend on scale, with more variation on smaller scales. The clustering as a function of scale is not a power law, but instead has a dip at 1h^{-1} Mpc and an excess on small scales. The fraction of red galaxies within the L* sample surrounding LRGs is a strong function of scale, as expected. However, the fraction of red galaxies evolves in redshift similarly on small and large scales, suggesting that cluster and field populations are changing in the same manner. The results highlight the advantage on small scales of using cross-correlation methods as a means of avoiding shot noise in samples of rare galaxies.

We test the homogeneity of the Universe at $z\sim 0.3$ with the Luminous Red Galaxy (LRG) spectroscopic sample of the Sloan Digital Sky Survey. First, the mean number $N(R)$ of LRGs within completely surveyed LRG-centered spheres of comoving radius $R$ is shown to be proportional to $R^3$ at radii greater than $R\sim 70 h^{-1} \mathrm{Mpc}$. The test has the virtue that it does not rely on the assumption that the LRG sample has a finite mean density; its results show, however, that there \emph{is} such a mean density. Secondly, the survey sky area is divided into 10 disjoint solid angular regions and the fractional rms density variations of the LRG sample in the redshift range $0.2<z<0.35$ among these ($\sim 2\times10^7 h^{-3} \mathrm{Mpc^3}$) regions is found to be 7 percent of the mean density. This variance is consistent with typical biased \lcdm models and puts very strong constraints on the quality of SDSS photometric calibration.

We present H-alpha rotation curves of two low luminosity field galaxies with r-band absolute magnitudes M_r=-13.9^{+0.8}_{-0.5} and M_r=-14.7^{+0.3}_{-0.2} (for h=0.7; the large error bars reflect distance uncertainties). Most previously studied galaxies in this luminosity range are members of groups defined by brighter galaxies, but these two systems, selected from Blanton et al.'s (2004) sample of low luminosity galaxies in the Sloan Digital Sky Survey (SDSS), appear to have no bright companions. The measured rotation speeds at the outer extent of the H-alpha rotation curves, 34.8+/-3.8 km/s and 30.9+/-7.2 km/s, are larger than the values of 16.0^{+6.1}_{-5.4} km/s and 20.9^{+6.2}_{-5.2} km/s predicted by extrapolating the inverse Tully-Fisher relation of luminous SDSS galaxies to these faint luminosities. However, a previous HI measurement of the first galaxy shows that it has a gas mass similar to its stellar mass, and the total baryonic mass is consistent with that predicted by McGaugh et al.'s (2000) "baryonic Tully-Fisher relation." We find r-band dynamical mass-to-light ratios within the radii of the last H-alpha data points (about 1.8 disk scale lengths in each case) of 12.6^{+4.7}_{-4.5} Msun/Lsun and 4.8^{+2.5}_{-2.1} Msun/Lsun, much higher than the values ~ 1 Msun/Lsun expected for the stellar populations. The dynamical properties of these galaxies, including the rotation speeds and evidence for high gas fractions and dark matter domination within the luminous extent of the galaxy, are consistent with those of previously studied faint galaxies in nearby groups. Further studies of the SDSS sample will allow characterization of low luminosity galaxies over the full range of environments in which they reside.

The ellipticities and alignments of clusters of galaxies, and their evolution with redshift, are examined in the context of a Lambda-dominated cold dark matter cosmology. We use a large-scale, high-resolution N-body simulation to model the matter distribution in a light cone containing ~10^6 clusters out to redshifts of z=3. Cluster ellipticities are determined as a function of mass, radius, and redshift, both in 3D and in projection. We find strong cluster ellipticities: the mean ellipticity increases with redshift from 0.3 at z=0 to 0.5 at z=3, for both 3D and 2D ellipticities; the evolution is well-fit by e=0.33+0.05z. The ellipticities increase with cluster mass and with cluster radius; the main cluster body is more elliptical than the cluster cores, but the increase of ellipticities with redshift is preserved. Using the fitted cluster ellipsoids, we determine the alignment of clusters as a function of their separation. We find strong alignment of clusters for separations <100 Mpc/h; the alignment increases with decreasing separation and with increasing redshift. The evolution of clusters from highly aligned and elongated systems at early times to lower alignment and elongation at present reflects the hierarchical and filamentary nature of structure formation. These measures of cluster ellipticity and alignment will provide a new test of the current cosmological model when compared with upcoming cluster surveys.

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 present estimates of cosmological parameters from the application of the Karhunen-Loeve transform to the analysis of the 3D power spectrum of density fluctuations using Sloan Digital Sky Survey galaxy redshifts. We use Omega_m*h and f_b = Omega_b/Omega_m to describe the shape of the power spectrum, sigma8 for the (linearly extrapolated) normalization, and beta to parametrize linear theory redshift space distortions. On scales k < 0.16 h/Mpc, our maximum likelihood values are Omega_m*h = 0.264 +/-0.043, f_b = 0.286 +/- 0.065, sigma8 = 0.966 +/- 0.048, and beta = 0.45 +/- 0.12. When we take a prior on Omega_b from WMAP, we find Omega_m*h = 0.207 +/- 0.030, which is in excellent agreement with WMAP and 2dF. This indicates that we have reasonably measured the gross shape of the power spectrum but we have difficulty breaking the degeneracy between Omega_m*h and f_b because the baryon oscillations are not resolved in the current spectroscopic survey window function.

We present galaxy-galaxy lensing measurements over scales 0.025 to 10 Mpc/h in the Sloan Digital Sky Survey. Using a flux-limited sample of 127,001 lens galaxies with spectroscopic redshifts and mean luminosity <L> = L_* and 9,020,388 source galaxies with photometric redshifts, we invert the lensing signal to obtain the galaxy-mass correlation function xi_{gm}. We find xi_{gm} is consistent with a power-law, xi_{gm} = (r/r_0)^{-gamma}, with best-fit parameters gamma = 1.79 +/- 0.06 and r_0 = (5.4+/-0.7)(0.27/Omega_m)^{1/gamma} Mpc/h. At fixed separation, the ratio xi_{gg}/xi_{gm} = b/r where b is the bias and r is the correlation coefficient. Comparing to the galaxy auto-correlation function for a similarly selected sample of SDSS galaxies, we find that b/r is approximately scale independent over scales 0.2-6.7 Mpc/h, with mean <b/r> = (1.3+/-0.2)(Omega_m/0.27). We also find no scale dependence in b/r for a volume limited sample of luminous galaxies (-23.0 < M_r < -21.5). The mean b/r for this sample is <b/r>_{Vlim} = (2.0+/-0.7)(Omega_m/0.27). We split the lens galaxy sample into subsets based on luminosity, color, spectral type, and velocity dispersion, and see clear trends of the lensing signal with each of these parameters. The amplitude and logarithmic slope of xi_{gm} increases with galaxy luminosity. For high luminosities (L ~5 L_*), xi_{gm} deviates significantly from a power law. These trends with luminosity also appear in the subsample of red galaxies, which are more strongly clustered than blue galaxies.

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.

We study the recently discovered gravitational lens SDSS J1004+4112, the first quasar lensed by a cluster of galaxies. It consists of four images with a maximum separation of 14.62''. The system has been confirmed as a lensed quasar at z=1.734 on the basis of deep imaging and spectroscopic follow-up observations. We present color-magnitude relations for galaxies near the lens plus spectroscopy of three central cluster members, which unambiguously confirm that a cluster at z=0.68 is responsible for the large image separation. We find a wide range of lens models consistent with the data, but they suggest four general conclusions: (1) the brightest cluster galaxy and the center of the cluster potential well appear to be offset by several kpc; (2) the cluster mass distribution must be elongated in the North--South direction, which is consistent with the observed distribution of cluster galaxies; (3) the inference of a large tidal shear (~0.2) suggests significant substructure in the cluster; and (4) enormous uncertainty in the predicted time delays between the images means that measuring the delays would greatly improve constraints on the models. We also compute the probability of such large separation lensing in the SDSS quasar sample, on the basis of the CDM model. The lack of large separation lenses in previous surveys and the discovery of one in SDSS together imply a mass fluctuation normalization \sigma_8=1.0^{+0.4}_{-0.2} (95% CL), if cluster dark matter halos have an inner slope -1.5. Shallower profiles would require higher values of \sigma_8. Although the statistical conclusion might be somewhat dependent on the degree of the complexity of the lens potential, the discovery is consistent with the predictions of the abundance of cluster-scale halos in the CDM scenario. (Abridged)

Compact groups (CGs) of galaxies -- relatively poor groups of galaxies in which the typical separations between members is of the order of a galaxy diameter -- offer an exceptional laboratory for the study of dense galaxian environments with short (<1Gyr) dynamical time-scales. In this paper, we present an objectively defined catalog of CGs in 153 sq deg of the Sloan Digital Sky Survey Early Data Release (SDSS EDR). To identify CGs, we applied a modified version of Hickson's (1982) criteria aimed at finding the highest density CGs and thus reducing the number of chance alignments. Our catalog contains 175 CGs down to a limiting galaxy magnitude of r* = 21. The resulting catalog has a median depth of approximately z = 0.13, substantially deeper than previous CG catalogs. Since the SDSS will eventually image up to one quarter of the celestial sphere, we expect our final catalog, based upon the completed SDSS, will contain on the order of 5,000 - 10,000 CGs. This catalog will be useful for conducting studies of the general characteristics of CGs, their environments, and their component galaxies.

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

We study the evolution of the cluster correlation function and its richness-dependence from z = 0 to z = 3 using large-scale cosmological simulations. A standard flat LCDM model with \Omega_m = 0.3 and, for comparison, a tilted \Omega_m = 1 model, TSCDM, are used. The evolutionary predictions are presented in a format suitable for direct comparisons with observations. We find that the cluster correlation strength increases with redshift: high redshift clusters are clustered more strongly (in comoving scale) than low redshift clusters of the same mass. The increased correlations with redshift, in spite of the decreasing mass correlation strength, is caused by the strong increase in cluster bias with redshift: clusters represent higher density peaks of the mass distribution as the redshift increases. The richness-dependent cluster correlation function, presented as the correlation-scale versus cluster mean separation relation, R_0 - d, is found to be, remarkably, independent of redshift to z <~ 2 for LCDM and z <~ 1 for TCDM (for a fixed correlation function slope and cluster mass within a fixed comoving radius). The non-evolving R_0 - d relation implies that both the comoving clustering scale and the cluster mean separation increase with redshift for the same mass clusters so that the R_0 - d relation remains essentially unchanged. The evolution of the R_0 - d relation from z ~ 0 to z ~ 3 provides an important new tool in cosmology; it can be used to break degeneracies that exist at z ~ 0 and provide precise determination of cosmological parameters.

The cluster correlation function and its richness dependence are determined from 1108 clusters of galaxies -- the largest sample of clusters studied so far -- found in 379 deg^2 of Sloan Digital Sky Survey early data. The results are compared with previous samples of optically and X-ray selected clusters. The richness-dependent correlation function increases monotonically from an average correlation scale of ~ 12 h^{-1} Mpc for poor clusters to ~ 25 h^{-1} Mpc for the richer, more massive clusters with a mean separation of ~ 90 h^{-1} Mpc. X-ray selected clusters suggest slightly stronger correlations than optically selected clusters (~ 2-\sigma). The results are compared with large-scale cosmological simulations. The observed richness-dependent cluster correlation function is well represented by the standard flat LCDM model (\Omega_m ~= 0.3, h ~= 0.7), and is inconsistent with the considerably weaker correlations predicted by \Omega_m = 1 models. An analytic relation for the correlation scale versus cluster mean separation, r_0 - d, that best describes the observations and the LCDM prediction is r_0 ~= 2.6 \sqrt{d} (for d ~= 20 - 90 h^{-1} Mpc). Data from the complete Sloan Digital Sky Survey, when available, will greatly enhance the accuracy of the results and allow a more precise determination of cosmological parameters.

Over the last few years, the Sloan Digital Sky Survey (SDSS) has discovered several hundred quasars with redshift between 4.0 and 6.4. Including the effects of magnification bias, one expects a priori that an appreciable fraction of these objects are gravitationally lensed. We have used the Advanced Camera for Surveys on the Hubble Space Telescope to carry out a snapshot imaging survey of high-redshift SDSS quasars to search for gravitationally split lenses. This paper, the first in a series reporting the results of the survey, describes snapshot observations of four quasars at z = 5.74, 5.82, 5.99 and 6.30, respectively. We find that none of these objects has a lensed companion within 5 magnitudes with a separation larger than 0.3 arcseconds; within 2.5 magnitudes, we can rule out companions within 0.1 arcseconds. Based on the non-detection of strong lensing in these four systems, we constrain the z~6 luminosity function to a slope of beta>-4.63 (3 sigma), assuming a break in the quasar luminosity function at M_{1450}^*=-24.0. We discuss the implications of this constraint on the ionizing background due to quasars in the early universe. Given that these quasars are not highly magnified, estimates of the masses of their central engines by the Eddington argument must be taken seriously, possibly challenging models of black hole formation.

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

We report the discovery of a new gravitationally lensed quasar from the Sloan Digital Sky Survey, SDSS J090334.92+502819.2. This object was targeted for SDSS spectroscopy as a Luminous Red Galaxy (LRG), but manual examination of the spectrum showed the presence of a quasar at z= 3.6 in addition to a red galaxy at z=0.388, and the SDSS image showed a second possible quasar image nearby. Follow-up imaging and spectroscopy confirmed the lensing hypothesis. In images taken at the ARC 3.5-meter telescope, two quasars are separated by 2.8 arc-seconds; the lensing galaxy is clearly seen and is blended with one of the quasar images. Spectroscopy taken at the Keck II telescope shows that the quasars have identical redshifts of z=3.6 and both show the presence of the same broad absorption line-like troughs. We present simple lens models which account for the geometry and magnifications. The lens galaxy lies near two groups of galaxies and may be a part of them. The models suggest that the groups may contribute considerable shear and may have a strong effect on the lens configuration.

The distribution in color and absolute magnitude is presented for 55,158 galaxies taken from the Sloan Digital Sky Survey in the redshift range $0.08<z<0.12$, as a function of galaxy number overdensity in a line-of-sight cylinder of transverse radius $1 h^{-1} \mathrm{Mpc}$. In all environments, bulge-dominated galaxies (defined to be those with radial profiles best fit with large S\'ersic indices) form a narrow, well defined color--magnitude relation. Although the most luminous galaxies reside preferentially in the highest density regions, there is only a barely detectable variation in the color (zero-point) or slope of the color--magnitude relation ($<0.02 \mathrm{mag}$ in $^{0.1}[g-r]$ or $[B-V]$). These results constrain variations with environmental density in the ages or metallicities of typical bulge-dominated galaxies to be under 20 percent.

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

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

We investigate the evolution of the fractions of late type cluster galaxies as a function of redshift, using one of the largest, most uniform cluster samples available. The sample consists of 514 clusters of galaxies in the range 0.02<z<0.3 from the Sloan Digital Sky Survey Cut & Enhance galaxy cluster catalog. This catalog was created using a single automated cluster finding algorithm on uniform data from a single telescope, with accurate CCD photometry, thus, minimizing selection biases. We use four independent methods to analyze the evolution of the late type galaxy fraction. Specifically, we select late type galaxies based on: restframe g-r color, u-r color, galaxy profile fitting and concentration index. The first criterion corresponds to the one used in the classical Butcher-Oemler analyses. The last three criteria are more sensitive to the morphological type of the galaxies. In all four cases, we find an increase in the fraction of late type galaxies with increasing redshift, significant at the 99.9% level. The results confirm that cluster galaxies do change colors with redshift (the Butcher-Oemler effect) and, in addition, they change their morphology to later-type toward higher redshift -- indicating a morphological equivalent of the Butcher-Oemler effect. We also find a tendency of richer clusters to have lower fractions of late type galaxies. The trend is consistent with a ram pressure stripping model, where richer clusters have more effective ram pressure due to their higher temperature.

We present here results on the composite luminosity functions of galaxies in the clusters of galaxies selected from the SDSS Cut and Enhance cluster catalog (CE; Goto et al. 2001). We construct the composite luminosity function in the five SDSS bands, u,g,r,i and z, using 204 CE clusters ranging from z=0.02 to z=0.25. We use photometric redshifts to construct composite luminosity functions. The robustness of the weighting scheme was tested using Monte Carlo simulation. We find the slope of composite LFs become flatter toward redder color band. Comparing with field LFs of SDSS (Blanton et al. 2001), cluster LFs have brighter characteristic magnitude and flatter slopes in g, r, i and z bands. These results are consistent with the hypothesis that the cluster LF has two distinct underlying populations i.e. the bright end of the LF is dominated by bright early types that follow a gaussian--like luminosity distribution, while the faint--end of the cluster LF is a steep power--law like function dominated by star--forming (bluer) galaxies. We also study the composite luminosity functions for early--type and late--type galaxies using profile fits, a concentration parameter and u-r color to classify galaxy morphology. The strong dependence of LF on galaxy morphology is found.The faint end slope of the LF is always flatter for early--type galaxies than late--type regardless of passband and methodology. These results are consistent with the hypothesis that the cluster region are dominated by bright elliptical galaxies.

Understanding the clustering of galaxies has long been a goal of modern observational cosmology. Utilizing our photometric redshift technique a volume limited sample containing more than 2 million galaxies is constructed from the SDSS galaxy catalog. In the largest such analysis to date, we study the angular clustering as a function of luminosity and spectral type. Using Limber's equation we calculate the clustering length for the full data set as r0=5.77+/-0.10 Mpc/h. We find that r0 increases with luminosity by a factor of 1.6 over the sampled luminosity range, in agreement with previous redshift surveys. We also find that both the clustering length and the slope of the correlation function depend on the galaxy type. In particular, by splitting the galaxies in four groups by their rest-frame type we find a bimodal behavior in their clustering properties. Galaxies with spectral types similar to elliptical galaxies have a correlation length of 6.59 +/- 0.17 Mpc/h and a slope of the angular correlation function of 0.96 +/- 0.05 while blue galaxies have a clustering length of 4.51 +/- 0.19 Mpc/h and a slope of 0.68 +/- 0.09. The two intermediate color groups behave like their more extreme 'siblings', rather than showing a gradual transition in slope. We discuss these correlations in the context of current cosmological models for structure formation.

We present a catalog of 799 clusters of galaxies in the redshift range z_est = 0.05 - 0.3 selected from ~400 deg^2 of early SDSS commissioning data along the celestial equator. The catalog is based on merging two independent selection methods -- a color-magnitude red-sequence maxBCG technique (B), and a Hybrid Matched-Filter method (H). The BH catalog includes clusters with richness \Lambda >= 40 (Matched-Filter) and N_gal >= 13 (maxBCG), corresponding to typical velocity dispersion of \sigma_v >~ 400 km s^{-1} and mass (within 0.6 h^{-1) Mpc radius) >~ 5*10^{13} h^{-1} M_sun. This threshold is below Abell richness class 0 clusters. The average space density of these clusters is 2*10^{-5} h^3 Mpc^{-3}. All NORAS X-ray clusters and 53 of the 58 Abell clusters in the survey region are detected in the catalog; the 5 additional Abell clusters are detected below the BH catalog cuts. The cluster richness function is determined and found to exhibit a steeply decreasing cluster abundance with increasing richness. We derive observational scaling relations between cluster richness and observed cluster luminosity and cluster velocity dispersion; these scaling relations provide important physical calibrations for the clusters. The catalog can be used for studies of individual clusters, for comparisons with other sources such as X-ray clusters and AGNs, and, with proper correction for the relevant selection functions, also for statistical analyses of clusters.

We determine the linear amplitude of mass fluctuations in the universe, sigma_8, from the abundance of massive clusters at redshifts z=0.5 to 0.8. The evolution of massive clusters depends exponentially on the amplitude of mass fluctuations and thus provides a powerful measure of this important cosmological parameter. The relatively high abundance of massive clusters observed at z>0.5, and the relatively slow evolution of their abundance with time, suggest a high amplitude of mass fluctuations: sigma_8=0.9 +-10% for Omega_m=0.4, increasing slightly to sigma_8=0.95 for Omega_m=0.25 and sigma_8=1.0 for Omega_m=0.1 (flat CDM models). We use the cluster abundance observed at z=0.5 to 0.8 to derive a normalization relation from the high-redshift clusters, which is only weakly dependent on Omega_m: sigma_8*Omega_m^0.14 = 0.78 +-0.08. When combined with recent constraints from the present-day cluster mass function (sigma_8*Omega_m^0.6=0.33 +-0.03) we find sigma_8=0.98 +-0.1 and Omega_m=0.17 +-0.05. Low sigma_8 values (<0.7) are unlikely; they produce an order of magnitude fewer massive clusters than observed.

We present constraints on several cosmological parameters from a combined analysis of the most recent Cosmic Microwave Background anisotropy data and the Sloan Digital Sky Survey cluster mass function. We find that the combination of the two data sets breaks several degeneracies among the parameters and provides the following constraints: $\sigma_8=0.76\pm0.09$, $\Omega_m=0.26^{+0.06}_{-0.07}$, $h=0.66^{+0.05}_{-0.06}$, $n=0.96 \pm 0.05$, $\tau_c=0.07^{+0.07}_{-0.05}$.

We study the mean environments of galaxies in the Sloan Digital Sky Survey as a function of rest-frame luminosity and color. Overdensities in galaxy number are estimated in $8 h^{-1} \mathrm{Mpc}$ and $1 h^{-1} \mathrm{Mpc}$ spheres centered on $125,000$ galaxies taken from the SDSS spectroscopic sample. We find that, at constant color, overdensity is independent of luminosity for galaxies with the blue colors of spirals. This suggests that, at fixed star-formation history, spiral-galaxy mass is a very weak function of environment. Overdensity does depend on luminosity for galaxies with the red colors of early types; both low-luminosity and high-luminosity red galaxies are found to be in highly overdense regions.

The mass function of clusters of galaxies is determined from 400 deg^2 of early commissioning imaging data of the Sloan Digital Sky Survey; ~300 clusters in the redshift range z = 0.1 - 0.2 are used. Clusters are selected using two independent selection methods: a Matched Filter and a red-sequence color magnitude technique. The two methods yield consistent results. The cluster mass function is compared with large-scale cosmological simulations. We find a best-fit cluster normalization relation of sigma_8*omega_m^0.6 = 0.33 +- 0.03 (for 0.1 ~< omega_m ~< 0.4), or equivalently sigma_8 = (0.16/omega_m)^0.6. The amplitude of this relation is significantly lower than the previous canonical value, implying that either omega_m is lower than previously expected (omega_m = 0.16 if sigma_8 = 1) or sigma_8 is lower than expected (sigma_8 = 0.7 if omega_m = 0.3). The best-fit mass function parameters are omega_m = 0.19 (+0.08,-0.07) and sigma_8 = 0.9 (+0.3,-0.2). High values of omega_m (>= 0.4) and low sigma_8 (=< 0.6) are excluded at >~ 2 sigma.

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.

We present the topology of a volume-limited sample of 11,884 galaxies, selected from an apparent-magnitude limited sample of over 100,000 galaxies observed as part of the Sloan Digital Sky Survey (SDSS). The data currently cover three main regions on the sky: one in the Galactic north and one in the south, both at zero degrees declination, and one area in the north at higher declination. Each of these areas covers a wide range of survey longitude but a narrow range of survey latitude, allowing the two dimensional genus to be measured. The genus curves of the SDSS sub-samples are similar, after appropriately normalizing these measurements for the different areas. We sum the genus curves from the three areas to obtain the total genus curve of the SDSS. The total curve has a shape similar to the genus curve derived from mock catalogs drawn from the Hubble Volume Lambda CDM simulation and is similar to that of a Gaussian random field. Likewise, comparison with the genus of the 2dFGRS, after normalization for the difference in area, reveals remarkable similarity in the topology of these samples. We test for the effects of galaxy type segregation by splitting the SDSS data into thirds, based on the u^*- r^* colors of the galaxies, and measure the genus of the reddest and bluest sub-samples. This red/blue split in u^*- r^* is essentially a split by morphology (Strateva et al. 2001). We find that the genus curve for the reddest galaxies exhibits a ``meatball'' shift of the topology -- reflecting the concentration of red galaxies in high density regions -- compared to the bluest galaxies and the full sample, in agreement with predictions from simulations.

Galaxy masses can be estimated by a variety of methods; each applicable in different circumstances, and each suffering from different systematic uncertainties. Confirmation of results obtained by one technique with analysis by another is particularly important. Recent SDSS weak lensing measurements of the projected-mass correlation function reveal a linear relation between galaxy luminosities and the depth of their dark matter halos (measured on 260 \hinv kpc scales). In this work we use an entirely independent dynamical method to confirm these results. We begin by assembling a sample of 618 relatively isolated host galaxies, surrounded by a total of 1225 substantially fainter satellites. We observe the mean dynamical effect of these hosts on the motions of their satellites by assembling velocity difference histograms. Dividing the sample by host properties, we find significant variations in satellite velocity dispersion with host luminosity. We quantify these variations using a simple dynamical model, measuring \mtsd a dynamical mass within 260 \hinv kpc. The appropriateness of this mass reconstruction is checked by conducting a similar analysis within an N-body simulation. Comparison between the dynamical and lensing mass-to-light scalings shows reasonable agreement, providing some quantitative confirmation for the lensing results.

We report on exploratory Chandra observations of the three highest redshift quasars known (z = 5.82, 5.99, and 6.28), all found in the Sloan Digital Sky Survey. These data, combined with a previous XMM-Newton observation of a z = 5.74 quasar, form a complete set of color-selected, z > 5.7 quasars. X-ray emission is detected from all of the quasars at levels that indicate that the X-ray to optical flux ratios of z ~ 6 optically selected quasars are similar to those of lower redshift quasars. The observations demonstrate that it will be feasible to obtain quality X-ray spectra of z ~ 6 quasars with current and future X-ray missions.

We show that the observed mass-to-light ratio of galaxy clusters increases with cluster temperature as expected from cosmological simulations. Contrary to previous observational suggestions, we find a mild but robust increase of M/L from poor (T = 1 - 2 keV) to rich (T = 12 keV) clusters; over this range, the mean M/L_V increases by a factor of about 2. The best fit relation satisfies M/L_V = (170 +- 30) T(keV)^(0.3 +- 0.1) h at z=0, with a large scatter. This trend confirms predictions from cosmological simulations which show that the richest clusters are antibiased, with a higher ratio of mass per unit light than average. The antibias increases with cluster temperature. The effect is caused by the relatively older age of the high-density clusters, where light has declined more significantly than average since their earlier formation time. Combining the current observations with simulations, we find a global value of M/L_V = 240 +- 50 h, and a corresponding mass density of the universe of Omega_m = 0.17 +- 0.05.

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

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

One of the most vital observational clues for unraveling the origin of Brightest Cluster Galaxies (BCG) is the observed alignment of the BCGs with their host cluster and its surroundings. We have examined the BCG-cluster alignment effect, using clusters of galaxies detected from the Sloan Digital Sky Survey (SDSS). We find that the BCGs are preferentially aligned with the principal axis of their hosts, to a much higher redshift (z >~ 0.3) than probed by previous studies (z <~ 0.1). The alignment effect strongly depends on the magnitude difference of the BCG and the second and third brightest cluster members: we find a strong alignment effect for the dominant BCGs, while less dominant BCGs do not show any departure from random alignment with respect to the cluster. We therefore claim that the alignment process originates from the same process that makes the BCG grow dominant, be it direct mergers in the early stage of cluster formation, or a later process that resembles the galactic cannibalism scenario. We do not find strong evidence for (or against) redshift evolution between 0<z<0.45, largely due to the insufficient sample size (< 200 clusters). However, we have developed a framework by which we can examine many more clusters in an automated fashion for the upcoming SDSS cluster catalogs, which will provide us with better statistics for systematic investigations of the alignment with redshift, richness and morphology of both the cluster and the BCG.

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

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

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

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

High-resolution N-body simulations of four popular Cold Dark Matter cosmologies (LCDM, OCDM, QCDM, and tilted SCDM), each containing 10^5 clusters of galaxies in a cubic gigaparsec volume, are used to determine the evolution of the cluster mass function from z=3 to z=0. The large volume and high resolution of these simulations allow an accurate measure of the evolution of cosmologically important (but rare) massive clusters at high redshift. The simulated mass function is presented for cluster masses within several radii typically used observationally (R=0.5, 1.0, and 1.5 Mpc/h, both comoving and physical) in order to enable direct comparison with current and future observations. The simulated evolution is compared with current observations of massive clusters at redshifts 0.3<z<0.8. The Omega_m=1 tilted SCDM model, which exhibits very rapid evolution of the cluster abundance, produces too few clusters at z>0.3 and no massive clusters at z>0.5, in stark contradiction with observations. The Omega_m=0.3 models- LCDM, OCDM, and QCDM- all exhibit considerably weaker evolution and are consistent with current data. Among these low density models, OCDM evolves the least. These trends are enhanced at high redshift and can be used to discriminate between flat and open low density models. The simulated mass functions are compared with the Press-Schechter approximation. Standard Press-Schechter predicts too many low mass clusters at z=0, and too few clusters at higher redshift. We modify the approximation by a simple parameterization of the density contrast threshold for collapse, which has a redshift dependence. This modified Press-Schechter approximation provides a good fit to the simulated mass functions.

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

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

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

We use large-scale cosmological simulations to estimate the mass-to-light ratio of galaxy systems as a function of scale, and compare the results with observations of galaxies, groups, clusters, and superclusters of galaxies. We find remarkably good agreement between observations and simulations. Specifically, we find that the simulated mass-to-light ratio increases with scale on small scales and flattens to a constant value on large scales, as suggested by observations. We find that while mass typically follows light on large scales, high overdensity regions --- such as rich clusters and superclusters of galaxies --- exhibit higher M/L_B values than average, while low density regions exhibit lower M/L_B values; high density regions are thus antibiased in M/L_B, with mass more strongly concentrated than blue light. This is true despite the fact that the galaxy mass density is unbiased or positively biased relative to the total mass density in these regions. The M/L_B antibias is likely due to the relatively old age of the high density regions, where light has declined significantly since their early formation time, especially in the blue band which traces recent star formation. Comparing the simulated results with observations, we place a powerful constraint on the mass density of the universe; using, for the first time, the entire observed mass-to-light function, from galaxies to superclusters, we find Omega =0.16+/-0.05.

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

The Sloan Digital Sky Survey has obtained 450 square degrees of photometric scan data, in five bands (u',g',r',i',z'), which we use to identify clusters of galaxies. We illustrate how we do star-galaxy separation, and present a simple and elegant method of detecting overdensities in the galaxy distribution, using the Voronoi Tessellation.

The "cosmic triangle" is introduced as a way of representing the past, present, and future status of the universe. Our current location within the cosmic triangle is determined by the answers to three questions: How much matter is in the universe? Is the expansion rate slowing down or speeding up? And, is the universe flat? A review of recent observations suggests a universe that is lightweight (matter density about one-third the critical value), is accelerating, and is flat. The acceleration implies the existence of cosmic dark energy that overcomes the gravitational self-attraction of matter and causes the expansion to speed up.

We derive a simple empirical photometric redshift estimator using a training set of galaxies with multiband photometry and measured redshifts in the Hubble Deep Field (HDF). This estimator is model-independent; it does not use spectral templates. The dispersion between the estimated redshifts and the spectroscopically measured ones is small; the dispersions range from $\sigma_z\simeq 0.03$ to 0.1 for $z\la 2$ galaxies, and from $\sigma_z\simeq 0.14$ to 0.25 for $z\ga 2$ galaxies. The predictions provided by our empirical redshift estimator agree well with recently measured galaxy redshifts. We illustrate how our empirical redshift estimator can be modified to include flat spectrum galaxies with $1.4\la z \la 2$.

We describe an automated method for detecting clusters of galaxies in imaging and redshift galaxy surveys. The Adaptive Matched Filter (AMF) method utilizes galaxy positions, magnitudes, and---when available---photometric or spectroscopic redshifts to find clusters and determine their redshift and richness. The AMF can be applied to most types of galaxy surveys: from two-dimensional (2D) imaging surveys, to multi-band imaging surveys with photometric redshifts of any accuracy (2.5D), to three-dimensional (3D) redshift surveys. The AMF can also be utilized in the selection of clusters in cosmological N-body simulations. The AMF identifies clusters by finding the peaks in a cluster likelihood map generated by convolving a galaxy survey with a filter based on a model of the cluster and field galaxy distributions. In tests on simulated 2D and 2.5D data with a magnitude limit of r' ~ 23.5, clusters are detected with an accuracy of Delta z ~ 0.02 in redshift and ~10% in richness to z < 0.5. Detecting clusters at higher redshifts is possible with deeper surveys. In this paper we present the theory behind the AMF and describe test results on synthetic galaxy catalogs.

I show that three independent methods utilizing clusters of galaxies - cluster dynamics and mass-to-light ratio, baryon fraction in clusters, and cluster evolution - all indicate the same robust result: the mass-density of the universe is low, Omega ~ 0.2, and the mass approximately traces light on large scales.

We derive simple empirical color-redshift relations for $z\la 4$ galaxies in the Hubble Deep Field (HDF) using a linear function of three photometric colors ($U-B$, $B-V$, $V-I$). The dispersion between the estimated redshifts and the spectroscopically observed ones is small for relations derived in several separate color regimes; the dispersions range from $\sigma_z\simeq 0.03$ to 0.1 for $z\la 2$ galaxies, and from $\sigma_z\simeq 0.14$ to 0.25 for $z\ga 2 $ galaxies. We apply the color-redshift relations to the HDF photometric catalog and obtain estimated redshifts that are consistent with those derived from spectral template fitting methods. The advantage of these color-redshift relations is that they are simple and easy to use and do not depend on the assumption of any particular spectral templates; they provide model independent redshift estimates for $z\la 4$ galaxies using only multi-band photometry, and they apply to about 90% of all galaxies. We provide a color-based estimated redshift catalog of HDF galaxies to $z\la 4$. We use the estimated redshifts to investigate the redshift distribution of galaxies in the HDF; we find peaks in the redshift distribution that suggest large-scale clustering of galaxies to at least $z\sim 1$ and that are consistent with those identified in spectroscopic probes of the HDF.

How much matter is there in the universe? Does the universe have the critical density needed to stop its expansion, or is the universe underweight and destined to expand forever? We show that several independent measures, especially those utilizing the largest bound systems known - clusters of galaxies - all indicate that the mass-density of the universe is insufficient to halt the expansion. A promising new method, the evolution of the number density of clusters with time, provides the most powerful indication so far that the universe has a sub-critical density. We show that different techniques reveal a consistent picture of a lightweight universe with only ~20-30% of the critical density. Thus, the universe may expand forever.

The existence of the three most massive clusters of galaxies observed so far at z>0.5 is used to constrain the mass density parameter of the universe, Omega, and the amplitude of mass fluctuations, sigma_8. We find Omega=0.2 (+0.3,-0.1), and sigma_8=1.2 (+0.5,-0.4) (95 %). We show that the existence of even the single most distant cluster at z=0.83, MS1054-03, with its large gravitational lensing mass, high temperature, and large velocity dispersion, is sufficient to establish powerful constraints. High-density, Omega=1 (sigma_8 ~ 0.5-0.6) Gaussian models are ruled out by these data (< 10^{-6} probability); the Omega=1 models predict only ~10^{-5} massive clusters at z > 0.65 (~10^{-3} at z > 0.5) instead of the 1 (3) clusters observed.

Clusters of galaxies, the most massive virialized systems known, provide a powerful tool for studying the structure, the mass density, and the cosmology of our universe. Clusters furnish one of the best estimates of the dynamical mass density parameter on 1 Mpc scale, $\Omega_{dyn}$; the best measure of the baryon density fraction in the universe; an excellent tracer of the large-scale structure of the universe; and, most recently, a powerful tracer of the evolution of structure in the universe and its unique cosmological implications. I review the above measures and show that they portray a consistent picture of the universe and place stringent constraints on cosmology. Each of these independent measures suggests a low-density universe, $\Omega_{m}$ $\simeq $ 0.3 $\pm $ 0.1, with mass approximately tracing light on large scales.

We present a method for determining the rms mass fluctuations on 8 h^-1 Mpc scale, sigma8. The method utilizes the rate of evolution of the abundance of rich clusters of galaxies. Using the Press-Schechter approximation, we show that the cluster abundance evolution is a strong function of sigma8: d log n/dz ~ -1/sigma8^2; low sigma8 models evolve exponentially faster than high sigma8 models, for a given mass cluster. For example, the number density of Coma-like clusters decreases by a factor of ~10^3$ from z = 0 to z ~ 0.5 for sigma8=0.5 models, while the decrease is only a factor of ~5 for sigma8 ~ 1. The strong exponential dependence on sigma8 arises because clusters represent rarer density peaks in low sigma8 models. We show that the evolution rate at z < 1 is insensitive to the density parameter $\Omega$ or to the exact shape of the power spectrum. Cluster evolution therefore provides a powerful constraint on sigma8. Using available cluster data to z ~ 0.8, we find sigma8 = 0.83 +/- 0.15. This amplitude implies a bias parameter b ~ 1/sigma8 = 1.2 +/- 0.2, i.e., a nearly unbiased universe with mass approximately tracing light on large scales.

We show that the evolution of the number density of rich clusters of galaxies breaks the degeneracy between Omega (the mass density ratio of the universe) and sigma_{8} (the normalization of the power spectrum), sigma_{8}Omega^{0.5} \simeq 0.5, that follows from the observed present-day abundance of rich clusters. The evolution of high-mass (Coma-like) clusters is strong in Omega=1, low-sigma_{8} models (such as the standard biased CDM model with sigma_{8} \simeq 0.5), where the number density of clusters decreases by a factor of \sim 10^{3} from z = 0 to z \simeq 0.5; the same clusters show only mild evolution in low-Omega, high-sigma_{8} models, where the decrease is a factor of \sim 10. This diagnostic provides a most powerful constraint on Omega. Using observations of clusters to z \simeq 0.5-1, we find only mild evolution in the observed cluster abundance. We find Omega = 0.3 \pm 0.1 and sigma_{8} = 0.85 \pm 0.15 (for Lambda = 0 models; for Omega + Lambda = 1 models, Omega = 0.34 \pm 0.13). These results imply, if confirmed by future surveys, that we live in a low-den sity, low-bias universe.

How is the universe organized on large scales? How did this structure evolve from the unknown initial conditions to the present time? The answers to these questions will shed light on the cosmology we live in, the amount, composition and distribution of matter in the universe, the initial spectrum of density fluctuations, and the formation and evolution of galaxies, clusters of galaxies, and larger scale structures. I review observational studies of large-scale structure, describe the peculiar velocity field on large scales, and present the observational constraints placed on cosmological models and the mass density of the universe. While major progress has been made over the last decade, most observational problems in the field of large-scale structure remain open. I will highlight some of the unsolved problems that seem likely to be solved in the next decade.

Rich clusters of galaxies are the most massive virialized systems known. Even though they contain only a small fraction of all galaxies, rich clusters provide a powerful tool for the study of galaxy formation, dark matter, large-scale structure, and cosmology. Superclusters, the largest known systems of galaxies, extend to (approx) 100h^{-1} Mpc in size and highlight the large-scale structure of the universe. This large-scale structure reflects initial conditions in the early universe and places strong contraints on models of galaxy formation and on cosmology. Some of the questions that can be addressed with clusters and superclusters of galaxies include: How did galaxies and larger structures form and evolve? What is the amount, composition, and distribution of matter in clusters and larger structures? How does the cluster mass density relate to the matter density in the universe? What constraints can the cluster and supercluster data place on cosmology? I will discuss in these lectures some of the properties of clusters and superclusters of galaxies that can be used to investigate these topics.

How is the universe organized on large scales? How did this structure evolve from the unknown initial conditions of a rather smooth early universe to the present time? The answers to these questions will shed light on the cosmology we live in, the amount, composition and distribution of matter in the universe, the initial spectrum of density fluctuations that gave rise to this structure, and the formation and evolution of galaxies, clusters of galaxies, and larger scale structures. To address these fundamental questions, large and accurate sky surveys are needed---in various wavelengths and to various depths. In this presentation I review current observational studies of large scale structure, present the constraints these observations place on cosmological models and on the amount of dark matter in the universe, and highlight some of the main unsolved problems in the field of large-scale structure that could be solved over the next decade with the aid of current and future surveys. I briefly discuss some of

Observations in the optical, in X-rays, and gravitational lensing of galaxies, clusters of galaxies, and large-scale structure are beginning to provide clues to the dark matter problem. I review the impact of these observations on some of the main questions relating to dark matter: How much dark matter is there? Where is it located? What is the nature of the dark matter? and what is the amount of baryonic dark matter.

The peculiar velocity function of clusters of galaxies is determined using an accurate sample of cluster velocities based on Tully-Fisher distances of Sc galaxies (Giovanelli et al 1995b). In contrast with previous results based on samples with considerably larger velocity uncertainties, the observed velocity function does not exhibit a tail of high velocity clusters. The results indicate a low probability of $\lesssim$\,5\% of finding clusters with one-dimensional velocities greater than $\sim$ 600 {\kms}. The root-mean-square one-dimensional cluster velocity is 293$\pm$28 {\kms}. The observed cluster velocity function is compared with expectations from different cosmological models. The absence of a high velocity tail in the observed function is most consistent with a low mass-density ($\Omega \sim$0.3) CDM model, and is inconsistent at $\gtrsim 3 \sigma$ level with $\Omega$= 1.0 CDM and HDM models. The root-mean-square one-dimensional cluster velocities in these models correspond, respectively, to 314, 516, and 632 {\kms} (when convolved with the observational uncertainties). Comparison with the observed RMS cluster velocity of 293$\pm$28 {\kms} further supports the low-density CDM model.

The observed baryon fraction and velocity--temperature relation in clusters of galaxies are compared with hydrodynamic simulations in two cosmological models : standard (Omega = 1) and a low-density flat (Omega=0.45 and \lambda=0.55) CDM models, normalized to the COBE background fluctuations. The observed properties of clusters include the velocity dispersion versus temperature relation, the gas mass versus total mass relation, and the gas mass fraction versus velocity dispersion relation. We find that, while both cosmological models reproduce well the shape of these observed functions, only low-density CDM can reproduce the observed amplitudes. The cluster gas mass fraction reflects approximately the baryon fraction in the models, with a slight anti-bias. Therefore, due to the low baryon density given by nucleosynthesis, Omega = 1 models produce too few baryons in clusters compared with observations. Scaling our results as a function of Omega, we find that a low-density CDM model, with Omega approximately 0.3 - 0.4, best reproduces the observed mean baryon fraction in clusters. The observed beta parameter of clusters, beta = sigma^2/(k T/mu m_p) = 0.94 \pm 0.08 discriminates less well between the models; it is consistent with that produced by low-density CDM (1.10 \pm 0.22), while it is slightly larger than expected but still consistent with Omega = 1 (0.70 \pm 0.14).

How much dark matter is there in the universe and where is it located? These are two of the most fundamental questions in cosmology. We use in this paper optical and x-ray mass determinations of galaxies, groups, and clusters of galaxies to suggest that most of the dark matter may reside in very large halos around galaxies, typically extending to ~200 kpc for bright galaxies. We show that the mass-to-light ratio of galaxy systems does not increase significantly with linear scale beyond the very large halos suggested for individual galaxies. Rather, the total mass of large scale systems such as groups and rich clusters of galaxies, even superclusters, can on average be accounted for by the total mass of their member galaxies, including their large halos (which may be stripped-off in the dense cluster environment but still remain in the clusters) plus the mass of the hot intracluster gas. This conclusion also suggests that we may live in a low-density universe with Omega ~ 0.2 - 0.3.

Rich clusters of galaxies, the largest virialized systems known, provide a powerful tool for the study of cosmology. Some of the fundamental questions that can be addressed with clusters of galaxies include: how did galaxies and large-scale structure form and evolve? What is the amount, composition and distribution of matter in the universe? I review some of the studies utilizing clusters of galaxies to investigare, among others: - The dark matter on clusters scale and the mean mass-density of the universe; - The large-scale structure of the universe; - The peculiar velocity field on large scales; - The mass-function of groups and clusters of galaxies; - The constraints placed on specific cosmological models using the cluster data.

The Sloan Digital Sky Survey (SDSS) will provide a complete imaging and spectroscopic survey of the high-latitude northern sky. The 2D survey will image the sky in five colors and will contain nearly 5 x 107 galaxies to g ~ 23m. The spectroscopic survey will obtain spectra of the brightest 106 galaxies, 105 quasars, and 103.5 rich clusters of galaxies (to g~18.3-19.3m, respectively). I summarize some of the science opportunities that will be made possible by this survey for studying the clustering and large-scale structure of the universe. The survey will identify a complete sample of several thousand rich clusters of galaxies, both in 2D and 3D - the largest automated sample yet available. The extensive cluster sample can be used to determine critical clustering properties such as the luminosity-function, velocity-function, and mass-function of clusters of galaxies (a critical test for cosmological models), detailed cluster dynamics and W(dyn), the cluster correlation function and its dependence on richness, cluster evolution, superclustering and voids to the largest scales yet observed, the motions of clusters and their large-scale peculiar velocity field, as well as detailed correlations between x-ray and optical properties of clusters, the density-morphology relation, and cluster-quasar associations. The large redshift survey, reaching to a depth of 600h-1 Mpc, will accurately map the largest scales yet observed, determine the power-spectrum and correlation function on these large scales for different type galaxies, and study the clustering of quasars to high redshifts (z 4). The implications of the survey for cosmological models, the dark matter, and W are also discussed.

I summarize the successes and failures of the most popular current cosmological models in accounting for the available observations. This evaluation, presented as the summary of the 11th Potsdam Workshop on Large-Scale Structure, was discussed by all the participants. I present the ratings given to each model based on a confrontation with the data, and the final vote taken.

The distributions of peculiar velocities of rich clusters and of groups of galaxies are investigated for different cosmological models and are compared with observations. Four cosmological models are studied: standard ($\Omega=1$) CDM, low-density ($\Omega=0.3$) CDM, HDM ($\Omega=1$), and PBI ($\Omega=0.3$), all COBE normalized. We find that rich clusters of galaxies exhibit a Maxwellian distribution of peculiar velocities in all models,as expected from a Gaussian initial density field. The cluster 3-D velocity distribution typically peaks at $v \sim 600$ km s$^{-1}$, and extends to high cluster velocities of $v \sim 2000$ km s$^{-1}$. The low-density CDM model exhibits somewhat lower velocities: it peaks at $\sim 400$ km s$^{-1}$, and extends to $\sim 1200$ km s$^{-1}$. Approximately 10\% ($\sim$1\% for low- density CDM) of all model rich clusters move with high peculiar velocities of $v \ge 10^3$ km s$^{-1}$. The highest velocity clusters frequently originate in dense superclusters. The model velocity distributions of groups and clusters of galaxies are compared with observations. The data exhibit a larger high-velocity tail, to $v_r \geq 2000$ km s$^{-1}$, than seen in the model simulations (except HDM). Due to the large observational uncertainties, however, the data are consistent at a $\sim 1$ to $3\sigma$ level with the model predictions and with a Gaussian initial density field. Accurate observations of cluster peculiar velocities, especially at the high-velocity tail, should provide powerful constraints on the cosmological models.

We determine the velocity correlation function, pairwise peculiar velocity difference, and root-mean-square pairwise peculiar velocity dispersion of rich clusters of galaxies, as a function of pair separation, for three cosmological models: Omega=1 and Omega=0.3 CDM, and Omega=0.3 PBI models (all flat and COBE-normalized). We find that close cluster pairs, with separation r<10Mpc/h, exhibit strong attractive peculiar velocities in all models; the cluster pairwise velocities depend sensitively on the model. The mean pairwise attractive velocity of clusters on 5Mpc/h scale ranges from 1700 km/s for Omega=1 CDM, to 1000 km/s for PBI, to 700 km/s for Omega=0.3 CDM. The small-scale pairwise velocities depend also on cluster mass: richer, more massive clusters exhibit stronger attractive velocities than less massive clusters. On large scales, from 20 to 200Mpc/h, the cluster peculiar velocities are increasingly dominated by bulk and random motions; they are independent of cluster mass. The cluster velocity correlation function, which reflects the bulk motion minus the relative motion of pairs, is negative on small scales for Omega=1 and Omega=0.3 CDM, and positive for PBI; this indicates stronger pairwise motion than bulk motion on small scales for CDM, and relatively larger bulk motions for PBI. The cluster velocity correlation function is positive on very large scales, from 10 to 200Mpc/h, for all models. These positive correlations, which decrease monotonically with scale, indicate significant bulk motions of clusters up to 200Mpc/h. The strong dependence of the cluster velocity functions on models, especially at small separations, makes them useful tools in constraining cosmological models when compared with observations.

What is the role of clusters of galaxies in probing the large-scale velocity field of the universe? We investigate the distribution of peculiar velocities of clusters of galaxies in the popular low-density ($\Omega=0.3$) flat Cold-Dark-Matter (CDM) cosmological model, which best fits many large-scale structure observations. An $\Omega=1$ CDM model is also studied for comparison. We find that clusters of galaxies are efficient tracers of the large-scale velocity field. The clusters exhibit a Maxwellian distribution of peculiar velocities, as expected from Gaussian initial density fluctuations. The cluster 3-D velocity distribution for the $\Omega=0.3$ model peaks at $v \sim 400$ km s$^{-1}$, and extends to high velocities of $v \sim 1200$ km s$^{-1}$. The rms peculiar velocity of the clusters is $440$ km s$^{-1}$. Approximately 10\% of all model clusters move with high peculiar velocities of $v \ge 700$ km s$^{-1}$. The observed velocity distribution of clusters of galaxies is compared with the predictions from cosmological models. The observed data exhibit a larger velocity tail than seen in the model simulations; however, due to the large observational uncertainties, the data are consistent at a $\sim 3\sigma$ level with the model predictions, and with a Gaussian initial density field. The large peculiar velocities reported for some clusters of galaxies ($v \geq 3000$ km s$^{-1}$) are likely to be overestimated, if the current model is viable.

The distributions of peculiar velocities of rich clusters and of groups of galaxies are investigated for different cosmological models and are compared with observations. Four cosmological models are studied: standard ($\Omega=1$) CDM, low-density CDM, HDM ($\Omega=1$), and PBI. We find that rich clusters of galaxies exhibit a Maxwellian distribution of peculiar velocities in all models, as expected from a Gaussian initial density fluctuation field. The cluster 3-D velocity distribution is generally similar in the models: it peaks at $v \sim 500$ km s$^{-1}$, and extends to high cluster velocities of $v \sim 1500$ km s$^{-1}$. Approximately 10\% of all model rich clusters move with high peculiar velocities of $v \ge 10^3$ km s$^{-1}$. The highest velocity clusters frequently originate in dense superclusters. The group velocity distribution is, in general, similar to the velocity distribution of the rich clusters. In all but the low-density CDM model, the mass exhibits a longer tail of high velocities than do the clusters. This high-velocity tail originates mostly from the high velocities that exist within rich clusters. The model velocity distributions of groups and clusters of galaxies are compared with observations. The data are generally consistent with the models, but exhibit a somewhat larger high-velocity tail, to $v_r \sim 3000$ km s$^{-1}$. While this high-velocity tail is similar to the HDM model predictions, the data are consistent with the other models studied, including the low-density CDM model, which best fits most other large-scale structure observations. The observed velocity

We show that the spatial correlation function of a flux-limited sample of X-ray selected clusters of galaxies will exhibit a correlation scale that is smaller than the correlation scale of a volume-limited, richness-limited sample of comparable apparent spatial density. The flux-limited sample contains clusters of different richnesses at different distances: poor groups are found nearby and rich clusters at greater distances. Since the cluster correlation strength is known to increase with richness, the flux-limited sample averages over the correlations of poor and rich clusters. On the other hand, a volume-limited, richness-limited sample has a minimum richness threshold, and a constant mixture of richnesses with redshift. Using the observed correlation scale for rich ($R\ge 1$) clusters, $r_o (R\ge 1) = 21\pm 2 h^{-1}$Mpc that was determined from previous volume-limited studies, we derive for the ROSAT flux-limited X-ray cluster sample $r_o({\rm flux}{\rm -}{\rm limited}) \approx 14h^{-1}$Mpc, in agreement with the recently observed value of $13.7\pm 2.3 h^{-1}$Mpc.

Previous comparisons of optical and X-ray observations of clusters of galaxies have lead to the so-called ``$\beta$ - discrepancy'' that has persisted for the last decade. The standard hydrostatic-isothermal model for clusters predicts that the parameter $\beta_{spec} = \sigma_{r}^{2}/(kT/\mu m_{p})$, which describes the ratio of energy per unit mass in galaxies to that in the gas, should equal the parameter $\beta_{fit}$ (where $\rho_{gas}(r) \propto \rho_{gal}(r)^{\beta_{fit}}$) determined from the X-ray surface brightness distribution. The observations suggest an apparent discrepancy : $\beta_{spec} \sim 1.2$ (i.e., the galaxies are ``hotter'' than the gas) while $\beta_{fit} \sim 0.65$ (i.e., the gas is ``hotter'' and more extended than the galaxies). Here we show that the discrepancy is resolved when the actual observed galaxy distribution in clusters is used, $\rho_{gal}(r) \propto r^{-2.4 \pm 0.2}$, instead of the previously assumed steeper King approximation, $\rho_{gal}(r) \propto r^{-3}$. Using a large sample of clusters, we find best-fit mean values of $\beta_{spec} = 0.94 \pm 0.08$ and $\beta_{fit}^{corrected} = 1.25 \times \beta_{fit} = 0.84 \pm 0.10$. These results resolve the $\beta$ - discrepancy and provide support for the hydrostatic cluster model.