Constraints on the mass distribution in high-redshift clusters of galaxies are not currently very strong. We aim to constrain the mass profile, M(r), and dynamical status of the $z \sim 0.8$ LCDCS 0504 cluster of galaxies characterized by prominent giant gravitational arcs near its center. Our analysis is based on deep X-ray, optical, and infrared imaging, as well as optical spectroscopy. We model the mass distribution of the cluster with three different mass density profiles, whose parameters are constrained by the strong lensing features of the inner cluster region, by the X-ray emission from the intra-cluster medium, and by the kinematics of 71 cluster members. We obtain consistent M(r) determinations from three methods (dispersion-kurtosis, caustics and MAMPOSSt), out to the cluster virial radius and beyond. The mass profile inferred by the strong lensing analysis in the central cluster region is slightly above, but still consistent with, the kinematics estimate. On the other hand, the X-ray based M(r) is significantly below both the kinematics and strong lensing estimates. Theoretical predictions from $\Lambda$CDM cosmology for the concentration--mass relation are in agreement with our observational results, when taking into account the uncertainties in both the observational and theoretical estimates. There appears to be a central deficit in the intra-cluster gas mass fraction compared to nearby clusters. Despite the relaxed appearance of this cluster, the determinations of its mass profile by different probes show substantial discrepancies, the origin of which remains to be determined. The extension of a similar dynamical analysis to other clusters of the DAFT/FADA survey will allow to shed light on the possible systematics that affect the determination of mass profiles of high-z clusters, possibly related to our incomplete understanding of intracluster baryon physics.

Stripe 82 in the Sloan Digital Sky Survey was observed multiple times, allowing deeper images to be constructed by coadding the data. Here we analyze the ellipticities of background galaxies in this 275 square degree region, searching for evidence of distortions due to cosmic shear. The E-mode is detected in both real and Fourier space with $>5$-$\sigma$ significance on degree scales, while the B-mode is consistent with zero as expected. The amplitude of the signal constrains the combination of the matter density $\Omega_m$ and fluctuation amplitude $\sigma_8$ to be $\Omega_m^{0.7}\sigma_8 = 0.252^{+0.032}_{-0.052}$.

The shapes of distant galaxies are sheared by intervening galaxy clusters. We examine this effect in Stripe 82, a 275 square degree region observed multiple times in the Sloan Digital Sky Survey and coadded to achieve greater depth. We obtain a mass-richness calibration that is similar to other SDSS analyses, demonstrating that the coaddition process did not adversely affect the lensing signal. We also propose a new parameterization of the effect of tomography on the cluster lensing signal which does not require binning in redshift, and we show that using this parameterization we can detect tomography for stacked clusters at varying redshifts. Finally, due to the sensitivity of the tomographic detection to accurately marginalizing over the effect of the cluster mass, we show that tomography at low redshift (where dependence on exact cosmological models is weak) can be used to constrain mass profiles in clusters.

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

We present and describe a catalog of galaxy photometric redshifts (photo-z's) for the Sloan Digital Sky Survey (SDSS) Coadd Data. We use the Artificial Neural Network (ANN) technique to calculate photo-z's and the Nearest Neighbor Error (NNE) method to estimate photo-z errors for $\sim$ 13 million objects classified as galaxies in the coadd with $r < 24.5$. The photo-z and photo-z error estimators are trained and validated on a sample of $\sim 83,000$ galaxies that have SDSS photometry and spectroscopic redshifts measured by the SDSS Data Release 7 (DR7), the Canadian Network for Observational Cosmology Field Galaxy Survey (CNOC2), the Deep Extragalactic Evolutionary Probe Data Release 3(DEEP2 DR3), the VIsible imaging Multi-Object Spectrograph - Very Large Telescope Deep Survey (VVDS) and the WiggleZ Dark Energy Survey. For the best ANN methods we have tried, we find that 68% of the galaxies in the validation set have a photo-z error smaller than $\sigma_{68} =0.031$. After presenting our results and quality tests, we provide a short guide for users accessing the public data.

The study of intracluster light can help us to understand the mechanisms taking place in galaxy clusters, and to place constraints on the cluster formation history and physical properties. However, owing to the intrinsic faintness of ICL emission, most searches and detailed studies of ICL have been limited to redshifts z<0.4.We search for ICL in a subsample of ten clusters detected by the ESO Distant Cluster Survey (EDisCS), at redshifts 0.4<z<0.8, that are also part of our DAFT/FADA Survey. We analyze the ICL by applying the OV WAV package, a wavelet-based technique, to deep HST ACS images in the F814W filter and to V-band VLT/FORS2 images of three clusters. Detection levels are assessed as a function of the diffuse light source surface brightness using simulations. In the F814W filter images, we detect diffuse light sources in all the clusters, with typical sizes of a few tens of kpc (assuming that they are at the cluster redshifts). The ICL detected by stacking the ten F814W images shows an 8sigma detection in the source center extending over a ~50x50kpc2 area, with a total absolute magnitude of -21.6 in the F814W filter, equivalent to about two L* galaxies per cluster. We find a weak correlation between the total F814W absolute magnitude of the ICL and the cluster velocity dispersion and mass. There is no apparent correlation between the cluster mass-to-light ratio (M/L) and the amount of ICL, and no evidence for any preferential orientation in the ICL source distribution. We find no strong variation in the amount of ICL between z=0 and z=0.8. In addition, we find wavelet-detected compact objects (WDCOs) in the three clusters for which data in two bands are available; these objects are probably very faint compact galaxies that in some cases are members of the respective clusters. We have shown that ICL is important in clusters at least up to z=0.8.

We present measurements of two types of cluster galaxy alignments based on a volume limited and highly pure ($\ge$ 90%) sample of clusters from the GMBCG catalog derived from SDSS DR7. We detect a clear BCG alignment (the alignment of major axis of the BCG toward the distribution of cluster satellite galaxies). We find that the BCG alignment signal becomes stronger as the redshift and BCG absolute magnitude decrease, and becomes weaker as BCG stellar mass decreases. No dependence of the BCG alignment on cluster richness is found. We can detect a statistically significant ($\ge$ 3 sigma) satellite alignment (the alignment of the major axes of the cluster satellite galaxies toward the BCG) only when we use the isophotal fit position angles (PAs, hereafter), and the satellite alignment depends on the apparent magnitudes rather than the absolute magnitudes of the BCGs. This suggests the detected satellite alignment based on isophotoal PAs from the SDSS pipeline is possibly due to the contamination from the diffuse light of nearby BCGs. We caution that this should not be simply interpreted as non-existence of the satellite alignment, but rather that we cannot detect them with our current photometric SDSS data. We perform our measurements on both SDSS $r$ band and $i$ band data, but did not observe a passband dependence of the alignments.

Future orbiting observatories will survey large areas of sky in order to constrain the physics of dark matter and dark energy using weak gravitational lensing and other methods. Lossy compression of the resultant data will improve the cost and feasibility of transmitting the images through the space communication network. We evaluate the consequences of the lossy compression algorithm of Bernstein et al. (2010) for the high-precision measurement of weak-lensing galaxy ellipticities. This square-root algorithm compresses each pixel independently, and the information discarded is by construction less than the Poisson error from photon shot noise. For simulated space-based images (without cosmic rays) digitized to the typical 16 bits per pixel, application of the lossy compression followed by image-wise lossless compression yields images with only 2.4 bits per pixel, a factor of 6.7 compression. We demonstrate that this compression introduces no bias in the sky background. The compression introduces a small amount of additional digitization noise to the images, and we demonstrate a corresponding small increase in ellipticity measurement noise. The ellipticity measurement method is biased by the addition of noise, so the additional digitization noise is expected to induce a multiplicative bias on the galaxies' measured ellipticities. After correcting for this known noise-induced bias, we find a residual multiplicative ellipticity bias of m ~ -4x10^{-4}. This bias is small when compared to the many other issues that precision weak lensing surveys must confront, and furthermore we expect it to be reduced further with better calibration of ellipticity measurement methods.

We present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density \rho_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find \Omega_\Lambda = 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

Reducing the scatter between cluster mass and optical richness is a key goal for cluster cosmology from photometric catalogs. We consider various modifications to the red-sequence matched filter richness estimator of Rozo et al. (2009), and evaluate their impact on the scatter in X-ray luminosity at fixed richness. Most significantly, we find that deeper luminosity cuts can reduce the recovered scatter, finding that sigma_lnLX|lambda=0.63+/-0.02 for clusters with M_500c >~ 1.6e14 h_70^-1 M_sun. The corresponding scatter in mass at fixed richness is sigma_lnM|lambda ~ 0.2-0.3 depending on the richness, comparable to that for total X-ray luminosity. We find that including blue galaxies in the richness estimate increases the scatter, as does weighting galaxies by their optical luminosity. We further demonstrate that our richness estimator is very robust. Specifically, the filter employed when estimating richness can be calibrated directly from the data, without requiring a-priori calibrations of the red-sequence. We also demonstrate that the recovered richness is robust to up to 50% uncertainties in the galaxy background, as well as to the choice of photometric filter employed, so long as the filters span the 4000 A break of red-sequence galaxies. Consequently, our richness estimator can be used to compare richness estimates of different clusters, even if they do not share the same photometric data. Appendix 1 includes "easy-bake" instructions for implementing our optimal richness estimator, and we are releasing an implementation of the code that works with SDSS data, as well as an augmented maxBCG catalog with the lambda richness measured for each cluster.

We present a large catalog of optically selected galaxy clusters from the application of a new Gaussian Mixture Brightest Cluster Galaxy (GMBCG) algorithm to SDSS Data Release 7 data. The algorithm detects clusters by identifying the red sequence plus Brightest Cluster Galaxy (BCG) feature, which is unique for galaxy clusters and does not exist among field galaxies. Red sequence clustering in color space is detected using an Error Corrected Gaussian Mixture Model. We run GMBCG on 8240 square degrees of photometric data from SDSS DR7 to assemble the largest ever optical galaxy cluster catalog, consisting of over 55,000 rich clusters across the redshift range from 0.1 < z < 0.55. We present Monte Carlo tests of completeness and purity and perform cross-matching with X-ray clusters and with the maxBCG sample at low redshift. These tests indicate high completeness and purity across the full redshift range for clusters with 15 or more members.

The red sequence is an important feature of galaxy clusters and plays a crucial role in optical cluster detection. Measurement of the slope and scatter of the red sequence are affected both by selection of red sequence galaxies and measurement errors. In this paper, we describe a new error corrected Gaussian Mixture Model for red sequence galaxy identification. Using this technique, we can remove the effects of measurement error and extract unbiased information about the intrinsic properties of the red sequence. We use this method to select red sequence galaxies in each of the 13,823 clusters in the maxBCG catalog, and measure the red sequence ridgeline location and scatter of each. These measurements provide precise constraints on the variation of the average red galaxy populations in the observed frame with redshift. We find that the scatter of the red sequence ridgeline increases mildly with redshift, and that the slope decreases with redshift. We also observe that the slope does not strongly depend on cluster richness. Using similar methods, we show that this behavior is mirrored in a spectroscopic sample of field galaxies, further emphasizing that ridgeline properties are independent of environment.

We present the 'simage' software suite for the simulation of artificial extragalactic images, based empirically around real observations of the Hubble Ultra Deep Field (UDF). The simulations reproduce galaxies with realistic and complex morphologies via the modeling of UDF galaxies as shapelets. Images can be created in the B, V, i and z bands for both space- and ground-based telescopes and instruments. The simulated images can be produced for any required field size, exposure time, Point Spread Function (PSF), telescope mirror size, pixel resolution, field star density, and a variety of detector noise sources. It has the capability to create images with both a pre-determined number of galaxies or one calibrated to the number counts of pre-existing data sets such as the HST COSMOS survey. In addition, simple options are included to add a known weak gravitational lensing (both shear and flexion) to the simulated images. The software is available in Interactive Data Language (IDL) and can be freely downloaded for scientific, developmental and teaching purposes.

As a contribution to the understanding of the dark energy concept, the Dark energy American French Team (DAFT, in French FADA) has started a large project to characterize statistically high redshift galaxy clusters, infer cosmological constraints from Weak Lensing Tomography, and understand biases relevant for constraining dark energy and cluster physics in future cluster and cosmological experiments. The purpose of this paper is to establish the basis of reference for the photo-z determination used in all our subsequent papers, including weak lensing tomography studies. This project is based on a sample of 91 high redshift (z>0.4), massive clusters with existing HST imaging, for which we are presently performing complementary multi-wavelength imaging. This allows us in particular to estimate spectral types and determine accurate photometric redshifts for galaxies along the lines of sight to the first ten clusters for which all the required data are available down to a limit of I_AB=24/24.5 with the LePhare software. The accuracy in redshift is of the order of 0.05 for the range 0.2<z<1.5. We verified that the technique applied to obtain photometric redshifts works well by comparing our results to with previous works. In clusters, photoz accuracy is degraded for bright absolute magnitudes and for the latest and earliest type galaxies. The photoz accuracy also only slightly varies as a function of the spectral type for field galaxies. As a consequence, we find evidence for an environmental dependence of the photoz accuracy, interpreted as the standard used Spectral Energy Distributions being not very well suited to cluster galaxies. Finally, we modeled the LCDCS 0504 mass with the strong arcs detected along this line of sight.

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.

Measurements of X-ray scaling laws are critical for improving cosmological constraints derived with the halo mass function and for understanding the physical processes that govern the heating and cooling of the intracluster medium. In this paper, we use a sample of 206 X-ray selected galaxy groups to investigate the scaling relation between X-ray luminosity (Lx) and halo mass (M00) where M200 is derived via stacked weak gravitational lensing. This work draws upon a broad array of multi-wavelength COSMOS observations including 1.64 square degrees of contiguous imaging with the Advanced Camera for Surveys (ACS) and deep XMM-Newton/Chandra imaging. The combined depth of these two data-sets allows us to probe the lensing signals of X-ray detected structures at both higher redshifts and lower masses than previously explored. Weak lensing profiles and halo masses are derived for nine sub-samples, narrowly binned in luminosity and redshift. The COSMOS data alone are well fit by a power law, M200 ~ Lx^a, with a slope of a=0.66+-0.14. These results significantly extend the dynamic range for which the halo masses of X-ray selected structures have been measured with weak gravitational lensing. As a result, tight constraints are obtained for the slope of the M-Lx relation. The combination of our group data with previously published cluster data demonstrates that the M-Lx relation is well described by a single power law, a=0.64+-0.03, over two decades in mass, 10^13.5-10^15.5 h72^-1 Msun. These results are inconsistent at the 3.7 level with the self-similar prediction of a=0.75. We examine the redshift dependence of the M-Lx relation and find little evidence for evolution beyond the rate predicted by self-similarity from z ~ 0.25 to z ~ 0.8.

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

We present a new survey strategy to discover and study high redshift Type Ia supernovae (SNe Ia) using the Hubble Space Telescope (HST). By targeting massive galaxy clusters at 0.9<z<1.5, we obtain a twofold improvement in the efficiency of finding SNe compared to an HST field survey and a factor of three improvement in the total yield of SN detections in relatively dust-free red-sequence galaxies. In total, sixteen SNe were discovered at z>0.95, nine of which were in galaxy clusters. This strategy provides a SN sample that can be used to decouple the effects of host galaxy extinction and intrinsic color in high redshift SNe, thereby reducing one of the largest systematic uncertainties in SN cosmology.

We use the abundance and weak lensing mass measurements of the SDSS maxBCG cluster catalog to simultaneously constrain cosmology and the richness--mass relation of the clusters. Assuming a flat \LambdaCDM cosmology, we find \sigma_8(\Omega_m/0.25)^{0.41} = 0.832\pm 0.033 after marginalization over all systematics. In common with previous studies, our error budget is dominated by systematic uncertainties, the primary two being the absolute mass scale of the weak lensing masses of the maxBCG clusters, and uncertainty in the scatter of the richness--mass relation. Our constraints are fully consistent with the WMAP five-year data, and in a joint analysis we find \sigma_8=0.807\pm 0.020 and \Omega_m=0.265\pm 0.016, an improvement of nearly a factor of two relative to WMAP5 alone. Our results are also in excellent agreement with and comparable in precision to the latest cosmological constraints from X-ray cluster abundances. The remarkable consistency among these results demonstrates that cluster abundance constraints are not only tight but also robust, and highlight the power of optically-selected cluster samples to produce precision constraints on cosmological parameters.

We measure the logarithmic scatter in mass at fixed richness for clusters in the maxBCG cluster catalog, an optically selected cluster sample drawn from SDSS imaging data. Our measurement is achieved by demanding consistency between available weak lensing and X-ray measurements of the maxBCG clusters, and the X-ray luminosity--mass relation inferred from the 400d X-ray cluster survey, a flux limited X-ray cluster survey. We find \sigma_{\ln M|N_{200}}=0.45^{+0.20}_{-0.18} (95% CL) at N_{200} ~ 40, where N_{200} is the number of red sequence galaxies in a cluster. As a byproduct of our analysis, we also obtain a constraint on the correlation coefficient between \ln Lx and \ln M at fixed richness, which is best expressed as a lower limit, r_{L,M|N} >= 0.85 (95% CL). This is the first observational constraint placed on a correlation coefficient involving two different cluster mass tracers. We use our results to produce a state of the art estimate of the halo mass function at z=0.23 -- the median redshift of the maxBCG cluster sample -- and find that it is consistent with the WMAP5 cosmology. Both the mass function data and its covariance matrix are presented.

Minimizing the scatter between cluster mass and accessible observables is an important goal for cluster cosmology. In this work, we introduce a new matched filter richness estimator, and test its performance using the maxBCG cluster catalog. Our new estimator significantly reduces the variance in the L_X-richness relation, from \sigma_{\ln L_X}^2=(0.86\pm0.02)^2 to \sigma_{\ln L_X}^2=(0.69\pm0.02)^2. Relative to the maxBCG richness estimate, it also removes the strong redshift dependence of the richness scaling relations, and is significantly more robust to photometric and redshift errors. These improvements are largely due to our more sophisticated treatment of galaxy color data. We also demonstrate the scatter in the L_X-richness relation depends on the aperture used to estimate cluster richness, and introduce a novel approach for optimizing said aperture which can be easily generalized to other mass tracers.

We present a new measurement of the scaling relation between X-ray luminosity and total mass for 17,000 galaxy clusters in the maxBCG cluster sample. Stacking sub-samples within fixed ranges of optical richness, N_200, we measure the mean 0.1-2.4 keV X-ray luminosity, <L_X>, from the ROSAT All-Sky Survey. The mean mass, <M_200>, is measured from weak gravitational lensing of SDSS background galaxies (Johnston et al. 2007). For 9 <= N_200 < 200, the data are well fit by a power-law, <L_X>/10^42 h^-2 erg/s = (12.6+1.4-1.3 (stat) +/- 1.6 (sys)) (<M_200>/10^14 h^-1 M_sun)^1.65+/-0.13. The slope agrees to within 10% with previous estimates based on X-ray selected catalogs, implying that the covariance in L_X and N_200 at fixed halo mass is not large. The luminosity intercent is 30%, or 2\sigma, lower than determined from the X-ray flux-limited sample of Reiprich & Bohringer (2002), assuming hydrostatic equilibrium. This difference could arise from a combination of Malmquist bias and/or systematic error in hydrostatic mass estimates, both of which are expected. The intercept agrees with that derived by Stanek et al. (2006) using a model for the statistical correspondence between clusters and halos in a WMAP3 cosmology with power spectrum normalization sigma_8 = 0.85. Similar exercises applied to future data sets will allow constraints on the covariance among optical and hot gas properties of clusters at fixed mass.

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 present measurements of the excess mass-to-light ratio measured aroundMaxBCG galaxy clusters observed in the SDSS. This red sequence cluster sample includes objects from small groups with masses ranging from ~5x10^{12} to ~10^{15} M_{sun}/h. Using cross-correlation weak lensing, we measure the excess mass density profile above the universal mean \Delta \rho(r) = \rho(r) - \bar{\rho} for clusters in bins of richness and optical luminosity. We also measure the excess luminosity density \Delta l(r) = l(r) - \bar{l} measured in the z=0.25 i-band. For both mass and light, we de-project the profiles to produce 3D mass and light profiles over scales from 25 kpc/ to 22 Mpc/h. From these profiles we calculate the cumulative excess mass M(r) and excess light L(r) as a function of separation from the BCG. On small scales, where \rho(r) >> \bar{\rho}, the integrated mass-to-light profile may be interpreted as the cluster mass-to-light ratio. We find the M/L_{200}, the mass-to-light ratio within r_{200}, scales with cluster mass as a power law with index 0.33+/-0.02. On large scales, where \rho(r) ~ \bar{\rho}, the M/L approaches an asymptotic value independent of cluster richness. For small groups, the mean M/L_{200} is much smaller than the asymptotic value, while for large clusters it is consistent with the asymptotic value. This asymptotic value should be proportional to the mean mass-to-light ratio of the universe <M/L>. We find <M/L>/b^2_{ml} = 362+/-54 h (statistical). There is additional uncertainty in the overall calibration at the ~10% level. The parameter b_{ml} is primarily a function of the bias of the L <~ L_* galaxies used as light tracers, and should be of order unity. Multiplying by the luminosity density in the same bandpass we find \Omega_m/b^2_{ml} = 0.02+/-0.03, independent of the Hubble parameter.

Determining the scaling relations between galaxy cluster observables requires large samples of uniformly observed clusters. We measure the mean X-ray luminosity--optical richness (L_X--N_200) relation for an approximately volume-limited sample of more than 17,000 optically-selected clusters from the maxBCG catalog spanning the redshift range 0.1<z<0.3. By stacking the X-ray emission from many clusters using ROSAT All-Sky Survey data, we are able to measure mean X-ray luminosities to ~10% (including systematic errors) for clusters in nine independent optical richness bins. In addition, we are able to crudely measure individual X-ray emission from ~800 of the richest clusters. Assuming a log-normal form for the scatter in the L_X--N_200 relation, we measure \sigma_\ln{L}=0.86+/-0.03 at fixed N_200. This scatter is large enough to significantly bias the mean stacked relation. The corrected median relation can be parameterized by L_X = (e^\alpha)(N_200/40)^\beta 10^42 h^-2 ergs/s, where \alpha = 3.57+/-0.08 and \beta = 1.82+/-0.05. We find that X-ray selected clusters are significantly brighter than optically-selected clusters at a given optical richness. This selection bias explains the apparently X-ray underluminous nature of optically-selected cluster catalogs.

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

The distribution of galaxies in position and velocity around the centers of galaxy clusters encodes important information about cluster mass and structure. Using the maxBCG galaxy cluster catalog identified from imaging data obtained in the Sloan Digital Sky Survey, we study the BCG-galaxy velocity correlation function. By modeling its non-Gaussianity, we measure the mean and scatter in velocity dispersion at fixed richness. The mean velocity dispersion increases from 202+/-10 km/s for small groups to more than 854+/-102 km/s for large clusters. We show the scatter to be at most 40.5+/-3.5%, declining to 14.9+/-9.4% in the richest bins. We test our methods in the C4 cluster catalog, a spectroscopic cluster catalog produced from the Sloan Digital Sky Survey DR2 spectroscopic sample, and in mock galaxy catalogs constructed from N-body simulations. Our methods are robust, measuring the scatter to well within one-sigma of the true value, and the mean to within 10%, in the mock catalogs. By convolving the scatter in velocity dispersion at fixed richness with the observed richness space density function, we measure the velocity dispersion function of the maxBCG galaxy clusters. Although velocity dispersion and richness do not form a true mass-observable relation, the relationship between velocity dispersion and mass is theoretically well characterized and has low scatter. Thus our results provide a key link between theory and observations up to the velocity bias between dark matter and galaxies.

We interpret and model the statistical weak lensing measurements around 130,000 groups and clusters of galaxies in the Sloan Digital Sky Survey presented by Sheldon et al. 2007 (Paper I). We present non-parametric inversions of the 2D shear profiles to the mean 3D cluster density and mass profiles in bins of both optical richness and cluster i-band luminosity. We correct the inferred 3D profiles for systematic effects, including non-linear shear and the fact that cluster halos are not all precisely centered on their brightest galaxies. We also model the measured cluster shear profile as a sum of contributions from the brightest central galaxy, the cluster dark matter halo, and neighboring halos. We infer the relations between mean cluster virial mass and optical richness and luminosity over two orders of magnitude in cluster mass; the virial mass at fixed richness or luminosity is determined with a precision of 13% including both statistical and systematic errors. We also constrain the halo concentration parameter and halo bias as a function of cluster mass; both are in good agreement with predictions of LCDM models. The methods employed here will be applicable to deeper, wide-area optical surveys that aim to constrain the nature of the dark energy, such as the Dark Energy Survey, the Large Synoptic Survey Telescope and space-based surveys.

We perform a maximum likelihood analysis of the cluster abundance measured in the SDSS using the maxBCG cluster finding algorithm. Our analysis is aimed at constraining the power spectrum normalization $\sigma_8$, and assumes flat cosmologies with a scale invariant spectrum, massless neutrinos, and CMB and supernova priors Omega_m*h^2=0.128+/-0.01 and h=0.72+/-0.05 respectively. Following the method described in the companion paper Rozo et al. 2007, we derive \sigma_8=0.92+/-0.10$ (1-sigma) after marginalizing over all major systematic uncertainties. We place strong lower limits on the normalization, sigma_8>0.76 (95% CL) (>0.68 at 99% CL). We also find that our analysis favors relatively low values for the slope of the Halo Occupation Distribution (HOD), alpha=0.83+/-0.06. The uncertainties of these determinations will substantially improve upon completion of an ongoing campaign to estimate dynamical, weak lensing, and X-ray cluster masses in the SDSS maxBCG cluster sample.

With a primary goal of conducting precision weak lensing measurements from space, the COSMOS survey has imaged the largest contiguous area observed by the Hubble Space Telescope (HST) to date using the Advanced Camera for Surveys (ACS). This is the first paper in a series where we describe our strategy for addressing the various technical challenges in the production of weak lensing measurements from the COSMOS data. The COSMOS ACS catalog is constructed from 575 ACS/WFC tiles (1.64 deg^2) and contains a total 1.2x10^6 objects to a limiting magnitude of F814W=26.5. This catalog is made publicly available. The shapes of galaxies have been measured and corrected for the distortion induced by the time varying ACS Point Spread Function and for Charge Transfer Efficiency effects. Next, simulated images are used to derive the shear susceptibility factors that are necessary in order to transform shape measurements into unbiased shear estimators. Finally, for each galaxy, we derive a shape measurement error and utilize this quantity to extract the intrinsic shape noise of the galaxy sample. Interestingly, our results indicate that the intrinsic shape noise varies little with either size, magnitude or redshift. Representing a number density of 66 galaxies per arcmin^2, the final COSMOS weak lensing catalog contains 3.9x10^5 galaxies with accurate shape measurements. The properties of the COSMOS weak lensing catalog described throughout this paper will provide key input numbers for the preparation and design of next-generation wide field space missions.

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)

We present a three dimensional cosmic shear analysis of the Hubble Space Telescope COSMOS survey, the largest ever optical imaging program performed in space. We have measured the shapes of galaxies for the tell-tale distortions caused by weak gravitational lensing, and traced the growth of that signal as a function of redshift. Using both 2D and 3D analyses, we measure cosmological parameters Omega_m, the density of matter in the universe, and sigma_8, the normalization of the matter power spectrum. The introduction of redshift information tightens the constraints by a factor of three, and also reduces the relative sampling (or "cosmic") variance compared to recent surveys that may be larger but are only two dimensional. From the 3D analysis, we find sigma_8*(Omega_m/0.3)^0.44=0.866+^0.085_-0.068 at 68% confidence limits, including both statistical and potential systematic sources of error in the total budget. Indeed, the absolute calibration of shear measurement methods is now the dominant source of uncertainty. Assuming instead a baseline cosmology to fix the geometry of the universe, we have measured the growth of structure on both linear and non-linear physical scales. Our results thus demonstrate a proof of concept for tomographic analysis techniques that have been proposed for future weak lensing surveys by a dedicated wide-field telescope in space.

Measurements of galaxy cluster abundances, clustering properties, and mass to- light ratios in current and future surveys can provide important cosmological constraints. Digital wide-field imaging surveys, the recently-demonstrated fidelity of red-sequence cluster detection techniques, and a new generation of realistic mock galaxy surveys provide the means for construction of large, cosmologicallyinteresting cluster samples, whose selection and properties can be understood in unprecedented depth. We present the details of the "maxBCG" algorithm, a cluster-detection technique tailored to multi-band CCD-imaging data. MaxBCG primarily relies on an observational cornerstone of massive galaxy clusters: they are marked by an overdensity of bright, uniformly red galaxies. This detection scheme also exploits classical brightest cluster galaxies (BCGs), which are often found at the center of these same massive clusters. (ABRIDGED)

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

The Shear TEsting Programme (STEP) is a collaborative project to improve the accuracy and reliability of weak lensing measurement, in preparation for the next generation of wide-field surveys. We review sixteen current and emerging shear measurement methods in a common language, and assess their performance by running them (blindly) on simulated images that contain a known shear signal. We determine the common features of algorithms that most successfully recover the input parameters. We achieve previously unattained discriminatory precision in our analysis, via a combination of more extensive simulations, and pairs of galaxy images that have been rotated with respect to each other, thus removing noise from their intrinsic ellipticities. The robustness of our simulation approach is also confirmed by testing the relative calibration of methods on real data. Weak lensing measurement has improved since the first STEP paper. Several methods now consistently achieve better than 2% precision, and are still being developed. However, the simulations can now distinguish all methods from perfect performance. Our main concern continues to be the potential for a multiplicative shear calibration bias: not least because this can not be internally calibrated with real data. We determine which galaxy populations are responsible and, by adjusting the simulated observing conditions, we also investigate the effects of instrumental and atmospheric parameters. We have isolated several previously unrecognised aspects of galaxy shape measurement, in which focussed development could provide further progress towards the sub-percent level of precision desired for future surveys. [ABRIDGED]

In this paper we present the first results of a deep spectroscopic survey of faint quasars in the Sloan Digital Sky Survey (SDSS) Southern Survey, a deep survey carried out by repeatedly imaging a 270 deg^2 area. Quasar candidates were selected from the deep data with good completeness over 0<z<5, and 2 to 3 magnitudes fainter than the SDSS main survey. Spectroscopic follow-up was carried out on the 6.5m MMT with Hectospec. The preliminary sample of this SDSS faint quasar survey (hereafter SFQS) covers ~ 3.9 deg^2, contains 414 quasars, and reaches g=22.5. The overall selection efficiency is ~ 66% (~ 80% at g<21.5); the efficiency in the most difficult redshift range (2<z<3) is better than 40%. We use the 1/V_{a} method to derive a binned estimate of the quasar luminosity function (QLF) and model the QLF using maximum likelihood analysis. The best model fits confirm previous results showing that the QLF has steep slopes at the bright end and much flatter slopes (-1.25 at z<2.0 and -1.55 at z>2.0) at the faint end, indicating a break in the QLF slope. Using a luminosity-dependent density evolution model, we find that the quasar density at M_{g}<-22.5 peaks at z~2, which is later in cosmic time than the peak of z~2.5 found from surveys of more luminous objects. The SFQS QLF is consistent with the results of the 2dF QSO Redshift Survey, the SDSS, and the 2dF-SDSS LRG and QSO Survey, but probes fainter quasars. We plan to obtain more quasars from future observations and establish a complete faint quasar sample with more than 1000 objects over 10 deg^2.

Using a sample of nearly 20,000 massive early-type galaxies selected from the Sloan Digital Sky Survey, we study the color-magnitude relation for the most luminous (L > 2.2 L^{*}) field galaxies in the redshift range 0.1<z<0.4 in several colors. The intrinsic dispersion in galaxy colors is quite small in all colors studied, but the 40 milli-mag scatter in the bluest colors is a factor of two larger than the 20 milli-mag measured in the reddest bands. While each of three simple models constructed for the star formation history in these systems can satisfy the constraints placed by our measurements, none of them produce color distributions matching those observed. Subdividing by environment, we find the dispersion for galaxies in clusters to be about 11% smaller than that of more isolated systems. Finally, having resolved the red sequence, we study the color dependence of the composite spectra. Bluer galaxies on the red sequence are found to have more young stars than red galaxies; the extent of this spectral difference is marginally better described by passive evolution of an old stellar population than by a model consisting of a recent trace injection of young stars.

In this paper we describe a detailed analysis of the photometric uncertainties present within the Sloan Digital Sky Survey (SDSS) imaging survey based on repeat observations of approximately 200 square degrees of the sky. We show that, for the standard SDSS aperture systems (petrocounts, counts_model, psfcounts and cmodel_counts), the errors generated by the SDSS photometric pipeline under-estimate the observed scatter in the individual bands. The degree of disagreement is a strong function of aperture and magnitude (ranging from 20% to more than a factor of 2). We also find that the photometry in the five optical bands can be highly correlated for both point sources and galaxies, although the correlation for point sources is almost entirely due to variable objects. Without correcting for this covariance the SDSS color errors could be in over-estimated by a factor of two to three. Combining these opposing effects, the SDSS errors on the colors differ from the observed color variation by approximately 10-20% for most apertures and magnitudes. We provide a prescription to correct the errors derived from the SDSS photometric pipeline as a function of magnitude and a semi-analytic method for generating the appropriate covariance between the different photometric passbands. Given the intrinsic nature of these correlations, we expect that all current and future multi-band surveys will also observe strongly covariant magnitudes. The ability of these surveys to complete their science goals is largely dependent on color-based target selection and photometric redshifts; these results show the importance of spending a significant fraction of early survey operations on re-imaging to empirically determine the photometric covariance of any observing/reduction pipeline.

We report follow-up observations of two gravitational lens candidates identified in the Sloan Digital Sky Survey (SDSS) dataset. We have confirmed that SDSS J102111.02+491330.4 is a previously unknown gravitationally lensed quasar. This lens system exhibits two images of a $z = 1.72$ quasar, with an image separation of $1{\farcs}14 \pm 0.04$. Optical and near-IR imaging of the system reveals the presence of the lensing galaxy between the two quasar images. Observations of SDSS J112012.12+671116.0 indicate that it is more likely a binary quasar than a gravitational lens. This system has two quasars at a redshift of $z = 1.49$, with an angular separation of $1{\farcs}49 \pm 0.02$. However, the two quasars have markedly different SEDs and no lens galaxy is apparent in optical and near-IR images of this system. We also present a list of 31 SDSS lens candidates which follow-up observations have confirmed are \textit{not} gravitational lenses.

We present a new non-parametric method for determining mean 3D density and mass profiles from weak lensing measurements around stacked samples of galaxies or clusters, that is, from measurement of the galaxy-shear or cluster-shear correlation functions. Since the correlation function is statistically isotropic, this method evades problems, such as projection of large-scale structure along the line of sight or halo asphericity, that complicate attempts to infer masses from weak lensing measurements of individual objects. We demonstrate the utility of this method in measuring halo profiles, galaxy-mass and cluster-mass cross-correlation functions, and cluster virial masses. We test this method on an N-body simulation and show that it correctly and accurately recovers the 3D density and mass profiles of halos. We find no evidence of problems due to a mass sheet degeneracy in the simulation. Cross-correlation lensing provides a powerful method for calibrating the mass-observable relation for use in measurement of the cluster mass function in large surveys. It can also be used on large scales to measure and remove the halo bias and thereby provide a direct measurement of Omega_m * sigma_8.

We report the discovery of a new two-image gravitationally lensed quasar, SDSS J024634.11-082536.2 (SDSS J0246-0825). This object was selected as a lensed quasar candidate from the Sloan Digital Sky Survey (SDSS) by the same algorithm that was used to discover other SDSS lensed quasars (e.g., SDSS J0924+0219). Multicolor imaging with the Magellan Consortium's Walter Baade 6.5-m telescope and the spectroscopic observations using the W. M. Keck Observatory's Keck II telescope confirm that SDSS J0246-0825 consists of two lensed images ($\Delta{\theta}=$1\farcs04) of a source quasar at z=1.68. Imaging observations with the Keck telescope and the Hubble Space Telescope reveal an extended object between the two quasar components, which is likely to be a lensing galaxy of this system. From the absorption lines in the spectra of quasar components and the apparent magnitude of the galaxy, combined with the expected absolute magnitude from the Faber-Jackson relation, we estimate the redshift of the lensing galaxy to be z=0.724. A highly distorted ring is visible in the Hubble Space Telescope images, which is likely to be the lensed host galaxy of the source quasar. Simple mass modeling predicts the possibility that there is a small (faint) lensing object near the primary lensing galaxy.

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.

In preparation for the Supernova Survey of the Sloan Digital Sky Survey (SDSS) II, a proposed 3-year extension to the SDSS, we have conducted an early engineering and science run during the fall of 2004, which consisted of approximately 20 scheduled nights of repeated imaging of half of the southern equatorial stripe. Transient supernova-like events were detected in near real-time and photometric measurements were made in the five SDSS filter bandpasses with a cadence of ~2 days. Candidate type Ia supernovae (SNe) were pre-selected based on their colors, light curve shape, and the properties of the host galaxy. Follow-up spectroscopic observations were performed with the Astrophysical Research Consortium 3.5m telescope and the 9.2m Hobby-Eberly Telescope to confirm their types and measure the redshifts. The 2004 campaign resulted in 22 spectroscopically confirmed SNe, which includes 16 type Ia, 5 type II, and 1 type Ib/c. These SN Ia will help fill in the sparsely sampled redshift interval of z = 0.05 - 0.35, the so-called 'redshift desert', in the Hubble diagram. Detailed investigation of the spectral properties of these moderate-redshift SNe Ia will also provide a bridge between local SNe and high-redshift objects, and will help us understand the systematics for future cosmological applications that require high photometric precision. Finally, the large survey volume also provides the opportunity to select unusual supernovae for spectroscopic study that are poorly sampled in other surveys. We report on some of the early results from this program and discuss potential future applications.

We report the discovery of a fifth image in the large separation lensed quasar system SDSS J1004+4112. A faint point source located 0.2'' from the center of the brightest galaxy in the lensing cluster is detected in images taken with the Advanced Camera for Surveys (ACS) and the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) on the Hubble Space Telescope. The flux ratio between the point source and the brightest lensed component in the ACS image is similar to that in the NICMOS image. The location and brightness of the point source are consistent with lens model predictions for a lensed image. We therefore conclude that the point source is likely to be a fifth image of the source quasar. In addition, the NICMOS image reveals the lensed host galaxy of the source quasar, which can strongly constrain the structure of the lensing critical curves and thereby the mass distribution of the lensing cluster.

We discuss the galaxy-galaxy-mass three-point correlation function and show how to measure it with weak gravitational lensing. The method entails choosing a large of pairs of foreground lens galaxies and constructing a mean shear map with respect to their axis, by averaging the ellipticities of background source galaxies. An average mass map can be reconstructed from this shear map and this will represent the average mass distribution around pairs of galaxies. We show how this mass map is related to the projected galaxy-galaxy-mass three-point correlation function. Using a large N-body dark matter simulation populated with galaxies using the Halo Occupation Distribution (HOD) bias prescription, we compute these correlation functions, mass maps, and shear maps. The resultant mass maps are distinctly bimodal, tracing the galaxy centers and remaining anisotropic up to scales much larger than the galaxy separation. At larger scales, the shear is approximately tangential about the center of the pair but with small azimuthal variation in amplitude. We estimate the signal-to-noise ratio of the reconstructed mass maps for a survey of similar depth to the Sloan Digital Sky Survey (SDSS) and conclude that the galaxy-galaxy-mass three-point function should be measurable with the current SDSS weak lensing data. Measurements of this three-point function, along with galaxy-galaxy lensing and galaxy auto-correlation functions, will provide new constraints on galaxy bias models. The anisotropic shear profile around close pairs of galaxies is a prediction of cold dark matter models and may be difficult to reconcile with alternative theories of gravity without dark matter.

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

We present 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 seven epochs of spectroscopy on the quadruply imaged quasar SDSS J1004+4112, spanning observed-frame time delays from 1 to 322 days. The spectra reveal differences in the emission lines between the lensed images. Specifically, component A showed a strong enhancement in the blue wings of several high-ionization lines relative to component B, which lasted at least 28 days (observed frame) then faded. Since the predicted time delay between A and B is <30 days, our time coverage suggests that the event was not intrinsic to the quasar. We attribute these variations to microlensing of part of the broad emission line region of the quasar, apparently resolving structure in the source plane on a scale of ~10^{16} cm at z=1.734. In addition, we observed smaller differences in the emission line profiles between components A and B that persisted throughout the time span, which may also be due to microlensing or millilensing. Further spectroscopic monitoring of this system holds considerable promise for resolving the structure of the broad emission line region in quasars.

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 report the discovery of the two-image gravitationally lensed quasar SDSS J1335+0118. The object was selected as a lens candidate from the Sloan Digital Sky Survey. The imaging and spectroscopic follow-up observations confirm that the system exhibits two gravitationally lensed images of a quasar at z=1.57. The image separation is 1.56''. We also detect an extended component between the two quasar images, likely the lensing galaxy. Preliminary mass modeling predicts the differential time delay \Delta t ~ 30 h^{-1} day assuming the redshift of the lens galaxy is 0.5.

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 report the discovery of SDSSJ115517.35+634622.0, a previously unknown gravitationally lensed quasar. The lens system exhibits two images of a $z = 2.89$ quasar, with an image separation of $1{\farcs}832 \pm 0.007$ . Near-IR imaging of the system reveals the presence of the lensing galaxy between the two quasar images. Based on absorption features seen in the Sloan Digital Sky Survey (SDSS) spectrum, we determine a lens galaxy redshift of $z = 0.1756$. The lens is rather unusual in that one of the quasar images is only $0{\farcs}22\pm0{\farcs}07$ ($\sim 0.1 R_{\rm eff}$) from the center of the lens galaxy and photometric modeling indicates that this image is significantly brighter than predicted by a SIS model. This system was discovered in the course of an ongoing search for strongly lensed quasars in the dataset from the SDSS.

We measure the large-scale real-space power spectrum P(k) using a sample of 205,443 galaxies from the Sloan Digital Sky Survey, covering 2417 square degrees with mean redshift z~0.1. We employ a matrix-based method using pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 22 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.02 h/Mpc < k < 0.3h/Mpc. We pay particular attention to modeling, quantifying and correcting for potential systematic errors, nonlinear redshift distortions and the artificial red-tilt caused by luminosity-dependent bias. Our final result is a measurement of the real-space matter power spectrum P(k) up to an unknown overall multiplicative bias factor. Our calculations suggest that this bias factor is independent of scale to better than a few percent for k<0.1h/Mpc, thereby making our results useful for precision measurements of cosmological parameters in conjunction with data from other experiments such as the WMAP satellite. As a simple characterization of the data, our measurements are well fit by a flat scale-invariant adiabatic cosmological model with h Omega_m =0.201+/- 0.017 and L* galaxy sigma_8=0.89 +/- 0.02 when fixing the baryon fraction Omega_b/Omega_m=0.17 and the Hubble parameter h=0.72; cosmological interpretation is given in a companion paper.

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

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

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

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.

We present measurements of the angular cross-correlation between luminous red galaxies from the Sloan Digital Sky Survey and the cosmic microwave background temperature maps from the Wilkinson Microwave Anisotropy Probe. We find a statistically significant achromatic positive correlation between these two data sets, which is consistent with the expected signal from the late Integrated Sachs-Wolfe (ISW) effect. We do not detect any anti-correlation on small angular scales as would be produced from a large Sunyaev-Zel'dovich (SZ) effect, although we do see evidence for some SZ effect for our highest redshift samples. Assuming a flat universe, our preliminary detection of the ISW effect provides independent physical evidence for the existence of dark energy.

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 report the discovery of a new gravitationally lensed quasar from the Sloan Digital Sky Survey, SDSS J092455.87+021924.9 (SDSS J0924+0219). This object was selected from among known SDSS quasars by an algorithm that was designed to select another known SDSS lensed quasar (SDSS 1226-0006A,B). Five separate components, three of which are unresolved, are identified in photometric follow-up observations obtained with the Magellan Consortium's 6.5m Walter Baade telescope at Las Campanas Observatory. Two of the unresolved components (designated A and B) are confirmed to be quasars with z=1.524; the velocity difference is less than 100 km sec^{-1} according to spectra taken with the W. M. Keck Observatory's Keck II telescope on Mauna Kea. A third stellar component, designated C, has the colors of a quasar with redshift similar to components A and B. The maximum separation of the point sources is 1.78". The other two sources, designated G and D, are resolved. Component G appears to be the best candidate for the lensing galaxy. Although component D is near the expected position of the fourth lensed component in a four image lens system, its properties are not consistent with being the image of a quasar at z~1.5. Nevertheless, the identical redshifts of components A and B and the presence of component C strongly suggest that this object is a gravitational lens. Our observations support the idea that a foreground object reddens the fourth lensed component and that another unmodeled effect (such as micro- or milli-lensing) demagnificates it, but we cannot rule out the possibility that SDSS0924+0219 is an example of the relatively rare class of ``three component'' lens systems.

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

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

We present initial results for counts in cells statistics of the angular distribution of galaxies in early data from the Sloan Digital Sky Survey (SDSS). We analyze a rectangular stripe $2.5^\circ$ wide, covering approximately 160 sq. degrees, containing over $10^6$ galaxies in the apparent magnitude range $18 < r^\prime < 22$, with areas of bad seeing, contamination from bright stars, ghosts, and high galactic extinction masked out. This survey region, which forms part of the SDSS Early Data Release, is the same as that for which two-point angular clustering statistics have recently been computed. The third and fourth moments of the cell counts, $s_3$ (skewness) and $s_4$ (kurtosis), constitute the most accurate measurements to date of these quantities (for $r^\prime < 21$) over angular scales $0.015^\circ-0.3^\circ$. They display the approximate hierarchical scaling expected from non-linear structure formation models and are in reasonable agreement with the predictions of $\Lambda$-dominated cold dark matter models with galaxy biasing that suppresses higher order correlations at small scales. The results are in general consistent with previous measurements in the APM, EDSGC, and Deeprange surveys. These results suggest that the SDSS imaging data are free of systematics to a high degree and will therefore enable determination of the skewness and kurtosis to 1% and less then 10%, as predicted by Colombi, Szapudi, & Szalay (1998).

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

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

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

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

We 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 present two weak lensing reconstructions of the nearby ($z_{cl}=0.055$) merging cluster Abell 3667, based on observations taken $\sim 1$ year apart under different seeing conditions. This is the lowest redshift cluster with a weak lensing mass reconstruction to date. The reproducibility of features in the two mass maps demonstrate that weak lensing studies of low-redshift clusters are feasible. These data constitute the first results from an X-ray luminosity-selected weak lensing survey of 19 low-redshift ($z< 0.1$) southern clusters.

We describe first results of a project to create weak lensing mass maps for a complete, X-ray luminosity-limited sample of 19 nearby (z < 0.1) southern galaxy clusters scheduled for Sunyaev-Zel'dovich observations by the Viper Telescope at the South Pole. We have collected data on 1/3 of the sample and present motivation for the project as well as projected mass maps of two clusters.