We present strong gravitational lensing models for 37 galaxy clusters from the SDSS Giant Arcs Survey. We combine data from multi-band Hubble Space Telescope WFC3imaging, with ground-based imaging and spectroscopy from Magellan, Gemini, APO, and MMT, in order to detect and spectroscopically confirm new multiply-lensed background sources behind the clusters. We report spectroscopic or photometric redshifts of sources in these fields, including cluster galaxies and background sources. Based on all available lensing evidence, we construct and present strong lensing mass models for these galaxy clusters.

We describe redMaPPer, a new red-sequence cluster finder specifically designed to make optimal use of ongoing and near-future large photometric surveys. The algorithm has multiple attractive features: (1) It can iteratively self-train the red-sequence model based on minimal spectroscopic training sample, an important feature for high redshift surveys; (2) It can handle complex masks with varying depth; (3) It produces cluster-appropriate random points to enable large-scale structure studies; (4) All clusters are assigned a full redshift probability distribution P(z); (5) Similarly, clusters can have multiple candidate central galaxies, each with corresponding centering probabilities; (6) The algorithm is parallel and numerically efficient: it can run a Dark Energy Survey-like catalog in ~500 CPU hours; (7) The algorithm exhibits excellent photometric redshift performance, the richness estimates are tightly correlated with external mass proxies, and the completeness and purity of the corresponding catalogs is superb. We apply the redMaPPer algorithm to ~10,000 deg^2 of SDSS DR8 data, and present the resulting catalog of ~25,000 clusters over the redshift range 0.08<z<0.55. The redMaPPer photometric redshifts are nearly Gaussian, with a scatter \sigma_z ~ 0.006 at z~0.1, increasing to \sigma_z~0.02 at z~0.5 due to increased photometric noise near the survey limit. The median value for |\Delta z|/(1+z) for the full sample is 0.006. The incidence of projection effects is low (<=5%). Detailed performance comparisons of the redMaPPer DR8 cluster catalog to X-ray and SZ catalogs are presented in a companion paper (Rozo & Rykoff 2014).

We present the discovery and a detailed multi-wavelength study of a strongly-lensed luminous infrared galaxy at z=0.816. Unlike most known lensed galaxies discovered at optical or near-infrared wavelengths this lensed source is red, r-Ks = 3.9 [AB], which the data presented here demonstrate is due to ongoing dusty star formation. The overall lensing magnification (a factor of 17) facilitates observations from the blue optical through to 500micron, fully capturing both the stellar photospheric emission as well as the re-processed thermal dust emission. We also present optical and near-IR spectroscopy. These extensive data show that this lensed galaxy is in many ways typical of IR-detected sources at z~1, with both a total luminosity and size in accordance with other (albeit much less detailed) measurements in samples of galaxies observed in deep fields with the Spitzer telescope. Its far-infrared spectral energy distribution is well-fit by local templates that are an order of magnitude less luminous than the lensed galaxy; local templates of comparable luminosity are too hot to fit. Its size (D~7kpc) is much larger than local luminous infrared galaxies, but in line with sizes observed for such galaxies at z~1. The star formation appears uniform across this spatial scale. In this source, the luminosity of which is typical of sources that dominate the cosmic infrared background, we find that star formation is spatially extended and well organised, quite unlike the compact merger-driven starbursts which are typical for sources of this luminosity at z~0.

We report the discovery of a unique gravitational lens system, SDSSJ2222+2745, producing five spectroscopically confirmed images of a z_s=2.82 quasar lensed by a foreground galaxy cluster at z_l=0.49. We also present photometric and spectroscopic evidence for a sixth lensed image of the same quasar. The maximum separation between the quasar images is 15.1". Both the large image separations and the high image multiplicity of the lensed quasar are in themselves exceptionally rare, and observing the combination of these two factors is an exceptionally unlikely occurrence in present datasets. This is only the third known case of a quasar lensed by a cluster, and the only one with six images. The lens system was discovered in the course of the Sloan Giant Arcs Survey, in which we identify candidate lenses in the Sloan Digital Sky Survey and target these for follow up and verification with the 2.56m Nordic Optical Telescope. Multi-band photometry obtained over multiple epochs from September 2011 to September 2012 reveal significant variability at the ~10-30% level in some of the quasar images, indicating that measurements of the relative time delay between quasar images will be feasible. In this lens system we also identify a bright (g = 21.5) giant arc corresponding to a strongly lensed background galaxy at z_s=2.30. We fit parametric models of the lens system, constrained by the redshift and positions of the quasar images and the redshift and position of the giant arc. The predicted time delays between different pairs of quasar images range from ~100 days to ~6 years.

We present new HST/WFC3 imaging data of RCSGA 032727-132609, a bright lensed galaxy at z=1.7 that is magnified and stretched by the lensing cluster RCS2 032727-132623. Using this new high-resolution imaging, we modify our previous lens model (which was based on ground-based data) to fully understand the lensing geometry, and use it to reconstruct the lensed galaxy in the source plane. This giant arc represents a unique opportunity to peer into 100-pc scale structures in a high-redshift galaxy. This new source reconstruction will be crucial for a future analysis of the spatially-resolved rest-UV and rest-optical spectra of the brightest parts of the arc.

Using the sample of Type Ia supernovae (SNe Ia) discovered by the Hubble Space Telescope (HST) Cluster Supernova Survey and augmented with HST-observed SNe Ia in the GOODS fields, we search for correlations between the properties of SNe and their host galaxies at high redshift. We use galaxy color and quantitative morphology to determine the red sequence in 25 clusters and develop a model to distinguish passively evolving early-type galaxies from star-forming galaxies in both clusters and the field. With this approach, we identify six SN Ia hosts that are early-type cluster members and eleven SN Ia hosts that are early-type field galaxies. We confirm for the first time at z>0.9 that SNe Ia hosted by early-type galaxies brighten and fade more quickly than SNe Ia hosted by late-type galaxies. We also show that the two samples of hosts produce SNe Ia with similar color distributions. The relatively simple spectral energy distributions (SEDs) expected for passive galaxies enable us to measure stellar masses of early-type SN hosts. In combination with stellar mass estimates of late-type GOODS SN hosts from Thomson & Chary (2011), we investigate the correlation of host mass with Hubble residual observed at lower redshifts. Although the sample is small and the uncertainties are large, a hint of this relation is found at z>0.9. By simultaneously fitting the average cluster galaxy formation history and dust content to the red-sequence scatters, we show that the reddening of early-type cluster SN hosts is likely E(B-V) <~ 0.06. The similarity of the field and cluster early-type host samples suggests that field early-type galaxies that lie on the red sequence may also be minimally affected by dust. Hence, the early-type hosted SNe Ia studied here occupy a more favorable environment to use as well-characterized high-redshift standard candles than other SNe Ia.

We present the discovery of an extremely bright and extended lensed source from the second Red Sequence Cluster Survey (RCS2). RCSGA 032727-132609 is spectroscopically confirmed as a giant arc and counter-image of a background galaxy at $z=1.701$, strongly-lensed by the foreground galaxy cluster RCS2 032727-132623 at $z=0.564$. The giant arc extends over $\sim 38$\,\arcsec and has an integrated $g$-band magnitude of 19.15, making it $\sim 20$ times larger and $\sim 4$ times brighter than the prototypical lensed galaxy MS1512-cB58. This is the brightest distant lensed galaxy in the Universe known to date. Its location in the `redshift desert' provides unique opportunities to connect between the large samples of galaxies known at $z\sim3$ and $z\sim1$. We have collected photometry in 9 bands, ranging from $u$ to $K_s$, which densely sample the rest-frame UV and optical light, including the age-sensitive 4000\AA\ break. A lens model is constructed for the system, and results in a robust total magnification of $2.04 \pm 0.16$ for the counter-image; we estimate an average magnification of $17.2 \pm 1.4$ for the giant arc based on the relative physical scales of the arc and counter-image. Fits of single-component spectral energy distribution (SED) models to the photometry result in a moderately young age, $t = 115 \pm 65$\,Myr, small amounts of dust, $E(B-V) \le 0.035$, and an exponentially declining star formation history with \textit{e}-folding time $\tau = 10-100$\,Myr. After correcting for the lensing magnification, we find a stellar mass of $\log(\mathrm{M}/\mathrm{M}_\odot)=10.0 \pm 0.1$. Allowing for episodic star formation, an underlying old burst could contain up to twice the mass inferred from single-component modeling. This stellar mass estimate is consistent with the average stellar mass of a sample of `BM' galaxies ($1.4 < z < 2.0$) studied by Reddy et al. (2006).

We report the discovery of 12 new fossil groups of galaxies, systems dominated by a single giant elliptical galaxy and cluster-scale gravitational potential, but lacking the population of bright galaxies typically seen in galaxy clusters. These fossil groups (FGs), selected from the maxBCG optical cluster catalog, were detected in snapshot observations with the Chandra X-ray Observatory. We detail the highly successful selection method, with an 80% success rate in identifying 12 FGs from our target sample of 15 candidates. For 11 of the systems, we determine the X-ray luminosity, temperature, and hydrostatic mass, which do not deviate significantly from expectations for normal systems, spanning a range typical of rich groups and poor clusters of galaxies. A small number of detected FGs are morphologically irregular, possibly due to past mergers, interaction of the intra-group medium (IGM) with a central AGN, or superposition of multiple massive halos. Two-thirds of the X-ray-detected FGs exhibit X-ray emission associated with the central BCG, although we are unable to distinguish between AGN and extended thermal galaxy emission using the current data. This sample, a large increase in the number of known FGs, will be invaluable for future planned observations to determine FG temperature, gas density, metal abundance, and mass distributions, and to compare to normal (non-fossil) systems. Finally, the presence of a population of galaxy-poor systems may bias mass function determinations that measure richness from galaxy counts. When used to constrain power spectrum normalization and {\Omega}_m, these biased mass functions may in turn bias these results.

We study the mass distribution of a sample of 28 galaxy clusters using strong and weak lensing observations. The clusters are selected via their strong lensing properties as part of the Sloan Giant Arcs Survey (SGAS) from the Sloan Digital Sky Survey (SDSS). Mass modelling of the strong lensing information from the giant arcs is combined with weak lensing measurements from deep Subaru/Suprime-cam images to primarily obtain robust constraints on the concentration parameter and the shape of the mass distribution. We find that the concentration c_vir is a steep function of the mass, c_vir \propto M_vir^-0.59\pm0.12, with the value roughly consistent with the lensing-bias-corrected theoretical expectation for high mass (10^15 h^-1 M_sun) clusters. However, the observationally inferred concentration parameters appear to be much higher at lower masses (10^14 h^-1 M_sun), possibly a consequence of the modification to the inner density profiles provided by baryon cooling. The steep mass-concentration relation is also supported from direct stacking analysis of the tangential shear profiles. In addition, we explore the two-dimensional shape of the projected mass distribution by stacking weak lensing shear maps of individual clusters with prior information on the position angle from strong lens modelling, and find significant evidence for a large mean ellipticity with the best-fit value of e = 0.47 \pm 0.06 for the mass distribution of the stacked sample. We find that the luminous cluster member galaxy distribution traces the overall mass distribution very well, although the distribution of fainter cluster galaxies appears to be more extended than the total mass.

We report a measurement of the Type Ia supernova (SN Ia) rate in galaxy clusters at 0.9 < z < 1.45 from the Hubble Space Telescope (HST) Cluster Supernova Survey. This is the first cluster SN Ia rate measurement with detected z > 0.9 SNe. Finding 8 +/- 1 cluster SNe Ia, we determine a SN Ia rate of 0.50 +0.23-0.19 (stat) +0.10-0.09 (sys) SNuB (SNuB = 10^-12 SNe L_{sun,B}^-1 yr^-1). In units of stellar mass, this translates to 0.36 +0.16-0.13 (stat) +0.07-0.06 (sys) SNuM (SNuM = 10^-12 SNe M_sun^-1 yr^-1). This represents a factor of approximately 5 +/- 2 increase over measurements of the cluster rate at z < 0.2. We parameterize the late-time SN Ia delay time distribution with a power law (proportional to t^s). Under the assumption of a cluster formation redshift of z_f = 3, our rate measurement in combination with lower-redshift cluster SN Ia rates constrains s = -1.41 +0.47/-0.40, consistent with measurements of the delay time distribution in the field. This measurement is generally consistent with expectations for the "double degenerate" scenario and inconsistent with some models for the "single degenerate" scenario predicting a steeper delay time distribution at large delay times. We check for environmental dependence and the influence of younger stellar populations by calculating the rate specifically in cluster red-sequence galaxies and in morphologically early-type galaxies, finding results similar to the full cluster rate. Finally, the upper limit of one host-less cluster SN Ia detected in the survey implies that the fraction of stars in the intra-cluster medium is less than 0.47 (95% confidence), consistent with measurements at lower redshifts.

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate out to z ~ 1.6 from the Hubble Space Telescope Cluster Supernova Survey. In observations spanning 189 orbits with the Advanced Camera for Surveys we discovered 29 SNe, of which approximately 20 are SNe Ia. Twelve of these SNe Ia are located in the foregrounds and backgrounds of the clusters targeted in the survey. Using these new data, we derive the volumetric SN Ia rate in four broad redshift bins, finding results consistent with previous measurements at z > 1 and strengthening the case for a SN Ia rate that is equal to or greater than ~0.6 x 10^-4/yr/Mpc^3 at z ~ 1 and flattening out at higher redshift. We provide SN candidates and efficiency calculations in a form that makes it easy to rebin and combine these results with other measurements for increased statistics. Finally, we compare the assumptions about host-galaxy dust extinction used in different high-redshift rate measurements, finding that different assumptions may induce significant systematic differences between measurements.

We present weak gravitational lensing analysis of 22 high-redshift (z >~1) clusters based on Hubble Space Telescope images. Most clusters in our sample provide significant lensing signals and are well detected in their reconstructed two-dimensional mass maps. Combining the current results and our previous weak-lensing studies of five other high-z clusters, we compare gravitational lensing masses of these clusters with other observables. We revisit the question whether the presence of the most massive clusters in our sample is in tension with the current LambdaCDM structure formation paradigm. We find that the lensing masses are tightly correlated with the gas temperatures and establish, for the first time, the lensing mass-temperature relation at z >~ 1. For the power law slope of the M-TX relation (M propto T^{\alpha}), we obtain \alpha=1.54 +/- 0.23. This is consistent with the theoretical self-similar prediction \alpha=3/2 and with the results previously reported in the literature for much lower redshift samples. However, our normalization is lower than the previous results by 20-30%, indicating that the normalization in the M-TX relation might evolve. After correcting for Eddington bias and updating the discovery area with a more conservative choice, we find that the existence of the most massive clusters in our sample still provides a tension with the current Lambda CDM model. The combined probability of finding the four most massive clusters in this sample after marginalization over current cosmological parameters is less than 1%.

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.

Maximizing the utility of upcoming photometric cluster surveys requires a thorough understanding of the richness-mass relation of galaxy clusters. We use Monte Carlo simulations to study the impact of various sources of observational scatter on this relation. Cluster ellipticity, photometric errors, photometric redshift errors, and cluster-to-cluster variations in the properties of red-sequence galaxies contribute negligible noise. Miscentering, however, can be important, and likely contributes to the scatter in the richness-mass relation of galaxy maxBCG clusters at the low mass end, where centering is more difficult. We also investigate the impact of projection effects under several empirically motivated assumptions about cluster environments. Using SDSS data and the maxBCG cluster catalog, we demonstrate that variations in cluster environments can rarely (\approx 1% - 5% of the time) result in significant richness boosts. Due to the steepness of the mass/richness function, the corresponding fraction of optically selected clusters that suffer from these projection effects is \approx 5% - 15%. We expect these numbers to be generic in magnitude, but a precise determination requires detailed, survey-specific modeling.

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 place constraints on the average density (Omega_m) and clustering amplitude (sigma_8) of matter using a combination of two measurements from the Sloan Digital Sky Survey: the galaxy two-point correlation function, w_p, and the mass-to-galaxy-number ratio within galaxy clusters, M/N, analogous to cluster M/L ratios. Our w_p measurements are obtained from DR7 while the sample of clusters is the maxBCG sample, with cluster masses derived from weak gravitational lensing. We construct non-linear galaxy bias models using the Halo Occupation Distribution (HOD) to fit both w_p and M/N for different cosmological parameters. HOD models that match the same two-point clustering predict different numbers of galaxies in massive halos when Omega_m or sigma_8 is varied, thereby breaking the degeneracy between cosmology and bias. We demonstrate that this technique yields constraints that are consistent and competitive with current results from cluster abundance studies, even though this technique does not use abundance information. Using w_p and M/N alone, we find Omega_m^0.5*sigma_8=0.465+/-0.026, with individual constraints of Omega_m=0.29+/-0.03 and sigma_8=0.85+/-0.06. Combined with current CMB data, these constraints are Omega_m=0.290+/-0.016 and sigma_8=0.826+/-0.020. All errors are 1-sigma. The systematic uncertainties that the M/N technique are most sensitive to are the amplitude of the bias function of dark matter halos and the possibility of redshift evolution between the SDSS Main sample and the maxBCG sample. Our derived constraints are insensitive to the current level of uncertainties in the halo mass function and in the mass-richness relation of clusters and its scatter, making the M/N technique complementary to cluster abundances as a method for constraining cosmology with future galaxy surveys.

We present results from a spectroscopic program targeting 26 strong lensing cluster cores that were visually identified in the Sloan Digital Sky Survey (SDSS) and the Second Red-Sequence Cluster Survey (RCS-2). The 26 galaxy cluster lenses span a redshift range of 0.2 < z < 0.65, and our spectroscopy reveals 69 unique background sources with redshifts as high as z=5.200. We also identify redshifts for 262 cluster member galaxies and measure the velocity dispersions and dynamical masses for 18 clusters where we have redshifts for N \geq 10 cluster member galaxies. We include an accounting for the expected biases in dynamical masses of strong lensing selected clusters as predicted by results from numerical simulations and discuss possible sources of bias in our observations. The median dynamical mass of the 18 clusters with N \geq 10 spectroscopic cluster members is M_{Vir} = 7.84 x 10^14 M_sun h_{0.7}^{-1}, which is somewhat higher than predictions for strong lensing selected clusters in simulations. The disagreement is not significant considering the large uncertainty in our dynamical data, systematic uncertainties in the velocity dispersion calibration, and limitations of the theoretical modeling. Nevertheless our study represents an important first step toward characterizing large samples of clusters that are identified in a systematic way as systems exhibiting dramatic strong lensing features.

We measure the redshift distribution of a sample of 28 giant arcs discovered as a part of the Sloan Giant Arcs Survey (SGAS). Gemini/GMOS-North spectroscopy provides precise redshifts for 24 arcs, and "redshift desert" constraints for the remaining four. This is a direct measurement of the redshift distribution of a uniformly selected sample of bright giant arcs, which is an observable that can be used to inform efforts to predict giant arc statistics. Our primary giant arc sample has a median redshift z=1.821 and nearly two thirds of the arcs - 64% - are sources at z \gtrsim 1.4, indicating that the population of background sources that are strongly lensed into bright giant arcs resides primarily at high redshift. We also analyze the distribution of redshifts for 19 secondary strongly lensed background sources that are not visually apparent in SDSS imaging, but were identified in deeper follow-up imaging of the lensing cluster fields. Our redshift sample for the secondary sources is not spectroscopically complete, but combining it with our primary giant arc sample suggests that a large fraction of all background galaxies which are strongly lensed by foreground clusters reside at z \gtrsim 1.4. Kolmogorov-Smirnov (KS) tests indicate that our well-selected, spectroscopically complete primary giant arc redshift sample can be reproduced with a model distribution that is constructed from a combination of results from studies of strong lensing clusters in numerical simulations, and observational constraints on the galaxy luminosity function.

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.

We have measured the Sunyaev Zel'dovich (SZ) effect for a sample of ten strong lensing selected galaxy clusters using the Sunyaev Zel'dovich Array (SZA). The SZA is sensitive to structures on spatial scales of a few arcminutes, while the strong lensing mass modeling constrains the mass at small scales (typically < 30"). Combining the two provides information about the projected concentrations of the strong lensing clusters. The Einstein radii we measure are twice as large as expected given the masses inferred from SZ scaling relations. A Monte Carlo simulation indicates that a sample randomly drawn from the expected distribution would have a larger median Einstein radius than the observed clusters about 3% of the time. The implied overconcentration has been noted in previous studies with smaller samples of lensing clusters. It persists for this sample, with the caveat that this could result from a systematic effect such as if the gas fractions of the strong lensing clusters are substantially below what is expected.

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 observations of two strongly lensed $z\sim5$ Lyman-$\alpha$ Emitting (LAE) galaxies that were discovered in the Sloan Giant Arcs Survey (SGAS). We identify the two sources as SGAS J091541+382655, at $z=5.200$, and SGAS J134331+415455 at $z=4.994$. We measure their AB magnitudes at $(i,z)=(23.34\pm0.09,23.29\pm0.13$) mags and $(i,z)=(23.78\pm0.18,24.24^{+0.18}_{-0.16}$) mags, and the rest-frame equivalent widths of the Lyman-$\alpha$ emission at $25.3\pm4.1$\AA~and $135.6\pm20.3$\AA~for SGAS J091541+382655 and SGAS J134331+415455, respectively. Each source is strongly lensed by a massive galaxy cluster in the foreground, and the magnifications due to gravitational lensing are recovered from strong lens modeling of the foreground lensing potentials. We use the magnification to calculate the intrinsic, unlensed Lyman-$\alpha$ and UV continuum luminosities for both sources, as well as the implied star formation rates (SFR). We find SGAS J091541+382655 and SGAS J134341+415455 to be galaxies with (L$_{Ly-\alpha}$, L$_{UV})\leq(0.6$L$_{Ly-\alpha}^{*}, 2$L$_{UV}^{*}$) and (L$_{Ly-\alpha}$, L$_{UV})=(0.5$L$_{Ly-\alpha}^{*}, 0.9$L$_{UV}^{*}$), respectively. Comparison of the spectral energy distributions (SEDs) of both sources against stellar population models produces estimates of the mass in young stars in each galaxy: we report an upper limit of M$_{stars} \leq 7.9^{+3.7}_{-2.5} \times 10^{7}$ M$_{\sun} h_{0.7}^{-1}$ for SGAS J091531+382655, and a range of viable masses for SGAS J134331+415455 of $2\times10^{8}$ M$_{\sun} h_{0.7}^{-1} <$ M$_{stars} < 6\times10^{9}$ M$_{\sun} h_{0.7}^{-1}$.

We report the discovery of two strongly-lensed z ~ 3 Lyman Break Galaxies (LBGs) discovered as u-band dropouts as part of the SDSS Giant Arcs Survey (SGAS). The first, SGAS J122651.3+215220 at z = 2.9233 is lensed by one of several sub-clusters, SDSS J1226+2152, in a complex massive cluster at z = 0.43. Its (g, r, i) magnitudes are (21.14, 20.60, 20.51) which translate to surface brightnesses, mu_{g,r,i}, of (23.78, 23.11, 22.81). The second, SGAS J152745.1+065219, is an LBG at z = 2.7593 lensed by the foreground SDSS J1527+0652 at z = 0.39, with (g, r, z)=(20.90, 20.52, 20.58) and mu_{g,r,z}=(25.15, 24.52, 24.12). Moderate resolution spectroscopy confirms the redshifts suggested by photometric breaks, and shows both absorption and emission features typical of LBGs. Lens mass models derived from combined imaging and spectroscopy reveal that SGAS J122651.3+215220 is a highly magnified source (M ~40), while SGAS J152745.1+065219 is magnified by no more than M ~ 15. Compared to LBG survey results (Steidel et al. 2003), the luminosities and lensing-corrected magnitudes suggest that SGAS J122651.3+215220 is among the faintest 20% of LBGs in that sample. SGAS J152745.1+065219, on the other hand, appears to be more representative of the average LBG, similar to the "Cosmic Eye".

We examine the alignment between Brightest Cluster Galaxies (BCGs) and their host clusters in a sample of 7031 clusters with 0.08<z<0.44 found using a matched-filter algorithm and an independent sample of 5744 clusters with 0.1<z<0.3 selected with the maxBCG algorithm, both extracted from the Sloan Digital Sky Survey Data Release 6 imaging data. We confirm that BCGs are preferentially aligned with the cluster's major axis; clusters with dominant BCGs (>0.65 mag brighter than the mean of the second and third ranked galaxies) show stronger alignment than do clusters with less dominant BCGs at the 4.4 sigma level. Rich clusters show a stronger alignment than do poor clusters at the 2.3 sigma level. Low redshift clusters (z<0.26) show more alignment than do high redshift (z>0.26) clusters, with a difference significant at the 3.0 sigma level. Our results do not depend on the algorithm used to select the cluster sample, suggesting that they are not biased by systematics of either algorithm. The correlation between BCG dominance and cluster alignment may be a consequence of the hierarchical merging process which forms the cluster. The observed redshift evolution may follow from secondary infall at late redshifts.

We present Chandra snapshot observations of the first large X-ray sample of optically identified fossil groups. For 9 of 14 candidate groups, we are able to determine the X-ray luminosity and temperature, which span a range typical of large ellipticals to rich groups of galaxies. We discuss these initial results in the context of group IGM and central galaxy ISM evolution, and we also describe plans for a deep X-ray follow-up program.

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.

Using deep two-band imaging from the Hubble Space Telescope, we measure the color-magnitude relations (CMR) of E/S0 galaxies in a set of 9 optically-selected clusters principally from the Red-Sequence Cluster Survey (RCS) at 0.9 < z < 1.23. We find that the mean scatter in the CMR in the observed frame of this set of clusters is 0.049 +/- 0.008, as compared to 0.031 +/- 0.007 in a similarly imaged and identically analyzed X-ray sample at similar redshifts. Single-burst stellar population models of the CMR scatter suggest that the E/S0 population in these RCS clusters truncated their star-formation at z~1.6, some 0.9 Gyrs later than their X-ray E/S0 counterparts which were truncated at z~2.1. The notion that this is a manifestation of the differing evolutionary states of the two populations of cluster galaxies is supported by comparison of the fraction of bulge-dominated galaxies found in the two samples which shows that optically-selected clusters contain a smaller fraction of E/S0 galaxies at the their cores

We derive radial mass profiles of four strong lensing selected clusters which show prominent giant arcs (Abell 1703, SDSS J1446+3032, SDSS J1531+3414, and SDSS J2111-0115), by combining detailed strong lens modeling with weak lensing shear measured from deep Subaru Suprime-cam images. Weak lensing signals are detected at high significance for all four clusters, whose redshifts range from z=0.28 to 0.64. We demonstrate that adding strong lensing information with known arc redshifts significantly improves constraints on the mass density profile, compared to those obtained from weak lensing alone. While the mass profiles are well fitted by the universal form predicted in N-body simulations of the LCDM model, all four clusters appear to be slightly more centrally concentrated (the concentration parameters c_vir \sim 8) than theoretical predictions, even after accounting for the bias toward higher concentrations inherent in lensing selected samples. Our results are consistent with previous studies which similarly detected a concentration excess, and increases the total number of clusters studied with the combined strong and weak lensing technique to ten. Combining our sample with previous work, we find that clusters with larger Einstein radii are more anomalously concentrated. We also present a detailed model of the lensing cluster A1703 with constraints from multiple image families, and find the dark matter inner density profile to be cuspy with the slope consistent with -1, in agreement with expectations.

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.

SN 2008D was discovered while following up an unusually bright X-ray transient in the nearby spiral galaxy NGC 2770. We present early optical spectra (obtained 1.75 days after the X-ray transient) which allowed the first identification of the object as a supernova at redshift z = 0.007. These spectra were acquired during the initial declining phase of the light curve, likely produced in the stellar envelope cooling after shock breakout, and rarely observed. They exhibit a rather flat spectral energy distribution with broad undulations, and a strong, W-shaped feature with minima at 3980 and 4190 AA (rest frame). We also present extensive spectroscopy and photometry of the supernova during the subsequent photospheric phase. Unlike supernovae associated with gamma-ray bursts, SN 2008D displayed prominent He features and is therefore of Type Ib.

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 on the detection of reddening toward z ~ 0.2 galaxy clusters. This is measured by correlating the Sloan Digital Sky Survey cluster and quasar catalogs and by comparing the photometric and spectroscopic properties of quasars behind the clusters to those in the field. We find mean E(B-V) values of a few times 10^-3 mag for sight lines passing ~Mpc from the clusters' center. The reddening curve is typical of dust but cannot be used to distinguish between different dust types. The radial dependence of the extinction is shallow near the cluster center suggesting that most of the detected dust lies at the outskirts of the clusters. Gravitational magnification of background z ~ 1.7 sources seen on Mpc (projected) scales around the clusters is found to be of order a few per cent, in qualitative agreement with theoretical predictions. Contamination by different spectral properties of the lensed quasar population is unlikely but cannot be excluded.

Imaging data from the Sloan Digital Sky Survey are used to characterize the population of galaxies in groups and clusters detected with the MaxBCG algorithm. We investigate the dependence of Brightest Cluster Galaxy (BCG) luminosity, and the distributions of satellite galaxy luminosity and satellite color, on cluster properties over the redshift range 0.1 < z < 0.3. The size of the dataset allows us to make measurements in many bins of cluster richness, radius and redshift. We find that, within r_200 of clusters with mass above 3e13 h-1 M_sun, the luminosity function of both red and blue satellites is only weakly dependent on richness. We further find that the shape of the satellite luminosity function does not depend on cluster-centric distance for magnitudes brighter than ^{0.25}M_i - 5log(h) < -19. However, the mix of faint red and blue galaxies changes dramatically. The satellite red fraction is dependent on cluster-centric distance, galaxy luminosity and cluster mass, and also increases by ~5% between redshifts 0.28 and 0.2, independent of richness. We find that BCG luminosity is tightly correlated with cluster richness, scaling as L_{BCG} ~ M_{200}^{0.3}, and has a Gaussian distribution at fixed richness, with sigma_{log L} ~ 0.17 for massive clusters. The ratios of BCG luminosity to total cluster luminosity and characteristic satellite luminosity scale strongly with cluster richness: in richer systems, BCGs contribute a smaller fraction of the total light, but are brighter compared to typical satellites. This study demonstrates the power of cross-correlation techniques for measuring galaxy populations in purely photometric data.

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.

We introduce a framework for describing the halo selection function of optical cluster finders. We treat the problem as being separable into a term that describes the intrinsic galaxy content of a halo (the Halo Occupation Distribution, or HOD) and a term that captures the effects of projection and selection by the particular cluster finding algorithm. Using mock galaxy catalogs tuned to reproduce the luminosity dependent correlation function and the empirical color-density relation measured in the SDSS, we characterize the maxBCG algorithm applied by Koester et al. to the SDSS galaxy catalog. We define and calibrate measures of completeness and purity for this algorithm, and demonstrate successful recovery of the underlying cosmology and HOD when applied to the mock catalogs. We identify principal components -- combinations of cosmology and HOD parameters -- that are recovered by survey counts as a function of richness, and demonstrate that percent-level accuracies are possible in the first two components, if the selection function can be understood to ~15% accuracy.

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)

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

The localization of the short-duration, hard-spectrum GRB 050509b was a watershed event. Thanks to the nearly immediate relay of the GRB position by Swift, we began imaging the GRB field 8 minutes after the burst and continued for the following 8 days. No convincing optical/infrared candidate afterglow or supernova was found for the object. We present a re-analysis of the XRT afterglow and find an absolute position that is ~4" to the west of the XRT position reported previously. Close to this position is a bright elliptical galaxy with redshift z=0.2248, about 1' from the center of a rich cluster of galaxies. Based on positional coincidences, the GRB and the bright elliptical are likely to be physically related. We thus have discovered evidence that at least some short-duration, hard-spectra GRBs arise at cosmological distances. However, while GRB 050509b was underluminous compared to long-duration GRBs, we demonstrate that the ratio of the blast-wave energy to the gamma-ray energy is consistent with that of long-duration GRBs. Based on this analysis, on the location of the GRB (40 +- 13 kpc from a bright galaxy), on the galaxy type (elliptical), and the lack of a coincident supernova, we suggest that there is now observational consistency with the hypothesis that short-hard bursts arise during the merger of a compact binary. We limit the properties of a Li-Paczynski ''mini-supernova.'' Other progenitor models are still viable, and additional rapidly localized bursts from the Swift mission will undoubtedly help to further clarify the progenitor picture. (abridged)