We present deep Chandra observations of A3411-12, a remarkable merging cluster that hosts the most compelling evidence for electron re-acceleration at cluster shocks to date. Using the $Y_X-M$ scaling relation, we find $r_{500} \sim 1.3$ Mpc, $M_{500} = (7.1 \pm 0.7) \times 10^{14} \ M_{\rm{\odot}}$, $kT=6.5\pm 0.1$ keV, and a gas mass of $M_{\rm g,500} = (9.7 \pm 0.1) \times 10^{13} M_\odot$. The gas mass fraction within $r_{500}$ is $f_{\rm g} = 0.14 \pm 0.01$. We compute the shock strength using density jumps to conclude that the Mach number of the merging subcluster is small ($M \leq 1.15_{-0.09}^{+0.14}$). We also present pseudo-density, projected temperature, pseudo-pressure, and pseudo-entropy maps. Based on the pseudo-entropy map we conclude that the cluster is undergoing a mild merger, consistent with the small Mach number. On the other hand, radio relics extend over Mpc scale in the A3411-12 system, which strongly suggests that a population of energetic electrons already existed over extended regions of the cluster.

Merging galaxy clusters may provide a unique window into the behavior of dark matter and the evolution of member galaxies. To interpret these natural collider experiments we must account for how much time has passed since pericenter passage (TSP), the maximum relative speed of the merging subclusters, merger phase (outbound after first pericenter or returning for second pericenter), and other dynamical parameters that are not directly observable. These quantities are often inferred from staged simulations or analytical timing arguments that include neither substructure, large-scale structure, nor a cosmologically motivated range of impact parameters. We include all these effects by extracting dynamical parameters from analog systems in a cosmological n-body simulation, and we present constraints for 11 observed systems. The TSP and viewing angles we derive are consistent with those of staged hydrodynamical simulations, but we find lower maximum speeds. Compared to the analytical MCMAC method we find lower TSP, and viewing angles that put the separation vector closer to the plane of the sky; we attribute this to the MCMAC assumption of zero pericenter distance. We discuss potential extensions to the basic analog method as well as complementarities between methods.

Weak gravitational lensing measurements based on photometry are limited by shape noise, the variance in the unknown pre-lensing orientations of the source galaxies. If the source is a disk galaxy with a well-ordered velocity field, however, velocity field data can support simultaneous inference of the shear, inclination, and position angle, virtually eliminating shape noise. We use the Fisher Information Matrix formalism to forecast the precision of this method in the idealized case of a perfectly ordered velocity field defined on an infinitesimally thin disk. For nearly face-on targets one shear component, $\gamma_\times$, can be constrained to $0.003\frac{90}{I_0}\frac{25}{n_{\rm pix}}$ where $I_0$ is the S/N of the central intensity pixel and $n_{\rm pix}$ is the number of pixels across a diameter enclosing enclosing 80% of the light. The other shear component, $\gamma_+$, is degenerate with the magnification $\mu$ but with a loose prior ($\sigma_\mu=1$) it can reach $0.009\frac{90}{I_0}\frac{25}{n_{\rm pix}}$. These constraints, however, degrade quickly at higher inclinations. We show that in general these constraints apply not to the shear components themselves but to two eigenvectors in the $(\gamma_\times,\gamma_+,\mu)$ space, with the third (usually $\mu$-like) eigenvector nearly unconstrained. We also forecast the potential of less expensive partial observations of the velocity field such as slit spectroscopy. We conclude by outlining some ways in which real galaxies depart from our idealized model, which may present barriers to practical implementation.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

A115 is a merging galaxy cluster at $z\sim0.2$ with a number of remarkable features including a giant ($\sim2.5$ Mpc) radio relic, two asymmetric X-ray peaks with trailing tails, and a peculiar line-of-sight velocity structure. We present a multi-wavelength study of A115 using optical imaging data from Subaru, X-ray data from $Chandra$, and spectroscopic data from the Keck/DEIMOS and MMT/Hectospec instruments. Our weak-lensing analysis shows that the cluster is comprised of two subclusters whose mass centroids are in excellent agreement with the two BCG positions ($\lesssim10$"). By modeling A115 with a superposition of two Navarro-Frenk-White halos, we determine the masses of the northern and southern subclusters to be $M_{200}=1.58_{-0.49}^{+0.56}\times 10^{14} \text{M}_{\odot}$ and $3.15_{-0.71}^{+0.79}\times 10^{14} \text{M}_{\odot}$, respectively. Combining the two halos, we estimate the total cluster mass to be $M_{200}=6.41_{-1.04}^{+1.08}\times10^{14} \text{M}_{\odot}$ at $R_{200}=1.67_{-0.09}^{+0.10}$ Mpc. These weak-lensing masses are significantly (a factor of 3-10) lower than what is implied by the X-ray and optical spectroscopic data. We attribute the difference to the gravitational and hydrodynamic disruption caused by the collision between the two subclusters.

The orientations of brightest cluster galaxies (BCGs) and their host clusters tend to be aligned, but the mechanism driving this is not clear. To probe the role of cluster mergers in this process, we quantify alignments of 38 BCGs in 22 clusters undergoing major mergers (up to $\sim 1$ Gyr after first pericenter). We find alignments entirely consistent with those of clusters in general. This suggests that alignments are robust against major cluster mergers. If, conversely, major cluster mergers actually help orient the BCG, such a process is acting quickly because the orientation is in place within $\sim 1$ Gyr after first pericenter.

In an influential recent paper, Harvey et al (2015) derive an upper limit to the self-interaction cross section of dark matter ($\sigma_{\rm DM} < 0.47$ cm$^2$/g at 95\% confidence) by averaging the dark matter-galaxy offsets in a sample of merging galaxy clusters. Using much more comprehensive data on the same clusters, we identify several substantial errors in their offset measurements. Correcting these errors relaxes the upper limit on $\sigma_{\rm DM}$ to $\lesssim 2$ cm$^2$/g, following the Harvey et al prescription for relating offsets to cross sections in a simple solid body scattering model. Furthermore, many clusters in the sample violate the assumptions behind this prescription, so even this revised upper limit should be used with caution. Although this particular sample does not tightly constrain self-interacting dark matter models when analyzed this way, we discuss how merger ensembles may be used more effectively in the future. We conclude that errors inherent in using single-band imaging to identify mass and light peaks do not necessarily average out in a sample of this size, particularly when a handful of substructures constitute a majority of the weight in the ensemble.

We present cosmological parameter measurements from the Deep Lens Survey (DLS) using galaxy-mass and galaxy-galaxy power spectra in the multipole range $\ell=250\sim2000$. We measure galaxy-galaxy power spectra from two lens bins centered at $z\sim0.27$ and $0.54$ and galaxy-mass power spectra by cross-correlating the positions of galaxies in these two lens bins with galaxy shapes in two source bins centered at $z\sim0.64$ and $1.1$. We marginalize over a baryonic feedback process using a single-parameter representation and a sum of neutrino masses, as well as photometric redshift and shear calibration systematic uncertainties. For a flat $\Lambda$CDM cosmology, we determine $S_8\equiv\sigma_8\sqrt{\Omega_m/0.3}=0.810^{+0.039}_{-0.031}$, in good agreement with our previous DLS cosmic shear and the Planck Cosmic Microwave Background (CMB) measurements. Without the baryonic feedback marginalization, $S_8$ decreases by $\sim0.05$ because the dark matter-only power spectrum lacks the suppression at the highest $\ell$'s due to Active Galactic Nuclei (AGN) feedback. Together with the Planck CMB measurement, we constrain the baryonic feedback parameter to $A_{baryon}=1.07^{+0.31}_{-0.39}$, which suggests an interesting possibility that the actual AGN feedback might be stronger than the recipe used in the OWLS simulations. The interpretation is limited by the validity of the baryonic feedback simulation and the one-parameter representation of the effect.

A key uncertainty in interpreting observations of bimodal merging galaxy clusters is the unknown angle between the subcluster separation vector and the plane of the sky. We present a new method for constraining this key parameter. We find analogs of observed systems in cosmological n-body simulations and quantify their likelihood of matching the observed projected separation and relative radial velocities between subclusters, as a function of viewing angle. We derive constraints on the viewing angle of many observed bimodal mergers including the Bullet Cluster (1E 0657-558) and El Gordo (ACT-CL J0102-4915). We also present more generic constraints as a function of projected separation and relative radial velocity, which can be used to assess additional clusters as information about them becomes available. The constraints from these two observables alone are weak (typically $\gtrsim 70-75^\circ$ at 68\% confidence and $\gtrsim 55-60^\circ$ at 95\% confidence) but incorporate much more cosmological context than the classical timing argument, marginalizing over many realizations of substructure, peculiar velocities, and so on. Compared to the MCMAC code, which implements the timing argument on NFW halos, our constraints generally predict subcluster separation vectors closer to the plane of the sky. This is because in realistic mergers the subcluster velocity vectors are not entirely parallel to the separation vector (i.e, the mergers are not perfectly head-on). As a result, observation of a nonzero relative radial velocity does not exclude a separation vector in the plane of the sky, as it does in the head-on timing argument employed by MCMAC.

Golovich et al. 2017b presents an optical imaging and spectroscopic survey of 29 radio relic merging galaxy clusters. In this paper, we study this survey to identify substructure and quantify the dynamics of the mergers. Using a combined photometric and spectroscopic approach, we identify the minimum number of substructures in each system to describe the galaxy populations and estimate the line of sight velocity difference between likely merging subclusters. We find that the line-of-sight velocity components of the mergers are typically small compared with the maximum three dimensional relative velocity (usually $<1000$ km s$^{-1}$ and often consistent with zero). This suggests that the merger axes of these systems are generally in or near the plane of the sky matching findings in magneto-hydrodynamical simulations. In 28 of the 29 systems we identify substructures in the galaxy population aligned with the radio relic(s) and presumed associated merger induced shock. From this ensemble, we identify eight systems to include in a `gold' sample that is prime for further observation, modeling, and simulation study. Additional papers will present weak lensing mass maps and dynamical modeling for each merging system, ultimately leading to new insight into a wide range of astrophysical phenomena at some of the largest scales in the universe.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

The second most significant detection of the Planck Sunyaev Zel'dovich survey, PLCK~G287.0+32.9 ($z=0.385$) boasts two similarly bright radio relics and a radio halo. One radio relic is located $\sim 400$ kpc northwest of the X-ray peak and the other $\sim 2.8$ Mpc to the southeast. This large difference suggests that a complex merging scenario is required. A key missing puzzle for the merging scenario reconstruction is the underlying dark matter distribution in high resolution. We present a joint Subaru Telescope and {\it Hubble Space Telescope} weak-lensing analysis of the cluster. Our analysis shows that the mass distribution features four significant substructures. Of the substructures, a primary cluster of mass $M_{200\text{c}}=1.59^{+0.25}_{-0.22}\times 10^{15} \ h^{-1}_{70} \ \text{M}_{\odot}$ dominates the weak-lensing signal. This cluster is likely to be undergoing a merger with one (or more) subcluster whose mass is approximately a factor of 10 lower. One candidate is the subcluster of mass $M_{200\text{c}}=1.16^{+0.15}_{-0.13}\times 10^{14} \ h^{-1}_{70} \ \text{M}_{\odot}$ located $\sim 400$ kpc to the southeast. The location of this subcluster suggests that its interaction with the primary cluster could be the source of the NW radio relic. Another subcluster is detected $\sim 2$ Mpc to the SE of the X-ray peak with mass $M_{200\text{c}}=1.68^{+0.22}_{-0.20}\times 10^{14} \ h^{-1}_{70} \ \text{M}_{\odot}$. This SE subcluster is in the vicinity of the SE radio relic and may have created the SE radio relic during a past merger with the primary cluster. The fourth subcluster, $M_{200\text{c}}=1.87^{+0.24}_{-0.22}\times 10^{14} \ h^{-1}_{70} \ \text{M}_{\odot}$, is northwest of the X-ray peak and beyond the NW radio relic.

Multi-band photometric and multi-object spectroscopic surveys of merging galaxy clusters allow for the characterization of the distributions of constituent dark matter and galaxy populations, constraints on the dynamics of the merging subclusters, and an understanding of galaxy evolution of member galaxies. We present deep photometric observations from Subaru/SuprimeCam and a catalog of $\sim$5400 spectroscopic cluster members from Keck/DEIMOS across 29 merging galaxy clusters ranging in redshift from $z=0.07$ to $0.55$. The ensemble is compiled based on the presence of radio relics, which highlight cluster scale collisionless shocks in the intra-cluster medium. Together with the spectroscopic and photometric information, the velocities, timescales, and geometries of the respective merging events may be tightly constrained. In this preliminary analysis, the velocity distributions of 28 of the 29 clusters are shown to be well fit by single Gaussians. This indicates that radio relic mergers largely occur transverse to the line of sight and/or near apocenter. In this paper, we present our optical and spectroscopic surveys, preliminary results, and a discussion of the value of radio relic mergers for developing accurate dynamical models of each system.

ZwCl 2341.1+0000, a merging galaxy cluster with disturbed X-ray morphology and widely separated ($\sim$3 Mpc) double radio relics, was thought to be an extremely massive ($10-30 \times 10^{14} M_\odot$) and complex system with little known about its merger history. We present JVLA 2-4 GHz observations of the cluster, along with new spectroscopy from our Keck/DEIMOS survey, and apply Gaussian Mixture Modeling to the three-dimensional distribution of 227 confirmed cluster galaxies. After adopting the Bayesian Information Criterion to avoid overfitting, which we discover can bias total dynamical mass estimates high, we find that a three-substructure model with a total dynamical mass estimate of $9.39 \pm 0.81 \times 10^{14} M_\odot$ is favored. We also present deep Subaru imaging and perform the first weak lensing analysis on this system, obtaining a weak lensing mass estimate of $5.57 \pm 2.47 \times 10^{14} M_\odot$. This is a more robust estimate because it does not depend on the dynamical state of the system, which is disturbed due to the merger. Our results indicate that ZwCl 2341.1+0000 is a multiple merger system comprised of at least three substructures, with the main merger that produced the radio relics occurring near to the plane of the sky, and a younger merger in the North occurring closer to the line of sight. Dynamical modeling of the main merger reproduces observed quantities (relic positions and polarizations, subcluster separation and radial velocity difference), if the merger axis angle of $\sim$10$^{+34}_{-6}$ degrees and the collision speed at pericenter is $\sim$1900$^{+300}_{-200}$ km/s.

Pairs of radio relics are believed to form during cluster mergers, and are best observed when the merger occurs in the plane of the sky. Mergers can also produce radio halos, through complex processes likely linked to turbulent re-acceleration of cosmic-ray electrons. However, only some clusters with double relics also show a radio halo. Here, we present a novel method to derive upper limits on the radio halo emission, and analyse archival X-ray Chandra data, as well as galaxy velocity dispersions and lensing data, in order to understand the key parameter that switches on radio halo emission. We place upper limits on the halo power below the $P_{\rm 1.4 \, GHz}\, M_{500}$ correlation for some clusters, confirming that clusters with double relics have different radio properties. Computing X-ray morphological indicators, we find that clusters with double relics are associated with the most disturbed clusters. We also investigate the role of different mass-ratios and time-since-merger. Data do not indicate that the merger mass ratio has an impact on the presence or absence of radio halos (the null hypothesis that the clusters belong to the same group cannot be rejected). However, the data suggests that the absence of radio halos could be associated with early and late mergers, but the sample is too small to perform a statistical test. Our study is limited by the small number of clusters with double relics. Future surveys with LOFAR, ASKAP, MeerKat and SKA will provide larger samples to better address this issue.

We analyze a rich dataset including Subaru/SuprimeCam, HST/ACS and WFC3, Keck/DEIMOS, Chandra/ACIS-I, and JVLA/C and D array for the merging galaxy cluster ZwCl 0008.8+5215. With a joint Subaru/HST weak gravitational lensing analysis, we identify two dominant subclusters and estimate the masses to be M$_{200}=\text{5.7}^{+\text{2.8}}_{-\text{1.8}}\times\text{10}^{\text{14}}\,\text{M}_{\odot}$ and 1.2$^{+\text{1.4}}_{-\text{0.6}}\times10^{14}$ M$_{\odot}$. We estimate the projected separation between the two subclusters to be 924$^{+\text{243}}_{-\text{206}}$ kpc. We perform a clustering analysis on confirmed cluster member galaxies and estimate the line of sight velocity difference between the two subclusters to be 92$\pm$164 km s$^{-\text{1}}$. We further motivate, discuss, and analyze the merger scenario through an analysis of the 42 ks of Chandra/ACIS-I and JVLA/C and D polarization data. The X-ray surface brightness profile reveals a remnant core reminiscent of the Bullet Cluster. The X-ray luminosity in the 0.5-7.0 keV band is 1.7$\pm$0.1$\times$10$^{\text{44}}$ erg s$^{-\text{1}}$ and the X-ray temperature is 4.90$\pm$0.13 keV. The radio relics are polarized up to 40$\%$. We implement a Monte Carlo dynamical analysis and estimate the merger velocity at pericenter to be 1800$^{+\text{400}}_{-\text{300}}$ km s$^{-\text{1}}$. ZwCl 0008.8+5215 is a low-mass version of the Bullet Cluster and therefore may prove useful in testing alternative models of dark matter. We do not find significant offsets between dark matter and galaxies, as the uncertainties are large with the current lensing data. Furthermore, in the east, the BCG is offset from other luminous cluster galaxies, which poses a puzzle for defining dark matter -- galaxy offsets.

We continue the study of the first sample of shear-selected clusters (Wittman et al. 2006) from the initial 8.6 square degrees of the Deep Lens Survey (DLS, Wittman et al. 2002); a sample with well-defined selection criteria corresponding to the highest ranked shear peaks in the survey area. We aim to characterize the weak lensing selection by examining the sample's X-ray properties. There are multiple X-ray clusters associated with nearly all the shear peaks: 14 X-ray clusters corresponding to seven DLS shear peaks. An additional three X-ray clusters cannot be definitively associated with shear peaks, mainly due to large positional offsets between the X-ray centroid and the shear peak. Here we report on the X-ray properties of the 17 X-ray clusters. The X-ray clusters display a wide range of luminosities and temperatures; the Lx-Tx relation we determine for the shear-associated X-ray clusters is consistent with X-ray cluster samples selected without regard to dynamical state, while it is inconsistent with self-similarity. For a subset of the sample, we measure X-ray masses using temperature as a proxy, and compare to weak lensing masses determined by the DLS team (Abate et al. 2009; Wittman et al. 2014). The resulting mass comparison is consistent with equality. The X-ray and weak lensing masses show considerable intrinsic scatter (~48%), which is consistent with X-ray selected samples when their X-ray and weak lensing masses are independently determined.

Dark matter with a non-zero self-interacting cross section ($\sigma_{\rm SIDM}$) has been posited as a solution to a number of outstanding astrophysical mysteries. Many studies of merging galaxy clusters have given constraints on $\sigma_{\rm SIDM}$ based on the spatial offset between the member galaxy population and the dark matter distribution. Assuming $\sigma_{\rm SIDM} = 0$, how likely is it for us to see the galaxy-DM offset values observed in merging clusters of galaxies? To answer this question, we formulate a hypothesis test using data from Illustris, a $\Lambda$CDM cosmological simulation. We select 43 Illustris clusters and their galaxy members at z~0 and examine the accuracy of commonly used galaxy summary statistics, including kernel-density-estimation (KDE) luminosity peak, KDE number density peak, shrinking aperture, centroid and the location of the brightest cluster galaxy (BCG). We use the dark-matter particles to reproduce commonly adopted methods to identify dark-matter peaks based on gravitational lensing cluster maps. By analysing each cluster in 768 projections, we determine the optimistic noise floor in the measurements of the galaxy-DM offsets. We find that the choice of the galaxy summary statistics affects the inferred offset values substantially, with the BCG and the luminosity peak giving the tightest 68-th percentile offset levels, $\lesssim$ 4 kpc and $\lesssim$ 32 kpc, respectively. Shrinking aperture, number density and centroid give a large offset scatter of about 50-100 kpc at the 68-th percentile level, even for clusters with only one dominant mass component. Out of the 15 reported offsets from observed merging clusters that we examined, 13 of them are consistent with Illustris unrelaxed cluster offsets at the 2-sigma (95-th percentile) level, i.e. consistent with the hypothesis that $\Lambda$CDM is the true underlying physical model.

On the largest scales, the Universe consists of voids and filaments making up the cosmic web. Galaxy clusters are located at the knots in this web, at the intersection of filaments. Clusters grow through accretion from these large-scale filaments and by mergers with other clusters and groups. In a growing number of galaxy clusters, elongated Mpc-size radio sources have been found, so-called radio relics. These relics are thought to trace relativistic electrons in the intracluster plasma accelerated by low-Mach number collisionless shocks generated by cluster-cluster merger events. A long-standing problem is how low-Mach number shocks can accelerate electrons so efficiently to explain the observed radio relics. Here we report on the discovery of a direct connection between a radio relic and a radio galaxy in the merging galaxy cluster Abell 3411-3412. This discovery indicates that fossil relativistic electrons from active galactic nuclei are re-accelerated at cluster shocks. It also implies that radio galaxies play an important role in governing the non-thermal component of the intracluster medium in merging clusters.

We investigate star formation in DLSCL J0916.2+2953, a dissociative merger of two clusters at z=0.53 that has progressed $1.1^{+1.3}_{-0.4}$ Gyr since first pass-through. We attempt to reveal the effects a collision may have had on the evolution of the cluster galaxies by tracing their star formation history. We probe current and recent activity to identify a possible star formation event at the time of the merger using EW(Hd), EW([OII]), and D$_{n}$4000 measured from the composite spectra of 64 cluster and 153 coeval field galaxies. We supplement $Keck$ DEIMOS spectra with DLS and $HST$ imaging to determine the color, stellar mass, and morphology of each galaxy and conduct a comprehensive study of the populations in this complex structure. Spectral results indicate the average cluster and cluster red sequence galaxies experienced no enhanced star formation relative to the surrounding field during the merger, ruling out a predominantly merger-quenched population. We find that the average blue galaxy in the North cluster is currently active and in the South cluster is currently post-starburst having undergone a recent star formation event. While the North activity could be latent or long-term merger effects, a young blue stellar population and irregular geometry suggest the cluster was still forming prior the collision. While the South activity coincides with the time of the merger, the blue early-type population could be a result of secular cluster processes. The evidence suggests that the dearth or surfeit of activity is indiscernible from normal cluster galaxy evolution.

Merging galaxy clusters have been touted as one of the best probes for constraining self-interacting dark matter, but few simulations exist to back up this claim. We simulate equal mass mergers of 10$^{15}$ M$_\odot$ halos, like the El Gordo and Sausage clusters, with cosmologically-motivated halo and merger parameters, and with velocity-independent dark-matter self-interactions. Although the standard lore for merging clusters is that self-interactions lead to large separations between the galaxy and dark-matter distributions, we find that maximal galaxy-dark matter offsets of $\lesssim~20$~kpc form for a self-interaction cross section of $\sigma_\text{SI}/m_\chi$ = 1 cm$^2$/g. This is an order of magnitude smaller than those measured in observed equal mass and near equal mass mergers, and is likely to be even smaller for lower-mass systems. While competitive cross-section constraints are thus unlikely to emerge from offsets, we find other signatures of self-interactions which are more promising. Intriguingly, we find that after dark matter halos coalesce, the collisionless galaxies (and especially the Brightest Cluster Galaxy [BGC]) oscillate around the center of the merger remnant on stable orbits of 100 kpc for $\sigma_\text{SI}/m_\chi = 1$~cm$^2$/g for at least several Gyr, well after the clusters have relaxed. If BCG miscentering in relaxed clusters remains a robust prediction of SIDM under the addition of gas physics, substructure, merger mass ratios (e.g., 10:1 like the Bullet Cluster), and complex cosmological merger histories, the observed BCG offsets may constrain $\sigma_\text{SI}/m_\chi \lesssim$ 0.1 cm$^2$/g---the tightest constraint yet.

We present an analysis of the merging cluster MACS J1149.5+2223 using archival imaging from Subaru/Suprime-Cam and multi-object spectroscopy from Keck/DEIMOS and Gemini/GMOS. We employ two and three dimensional substructure tests and determine that MACS J1149.5+2223 is composed of two separate mergers between three subclusters occurring $\sim$1 Gyr apart. The primary merger gives rise to elongated X-ray morphology and a radio relic in the southeast. The brightest cluster galaxy is a member of the northern subcluster of the primary merger. This subcluster is very massive (16.7$^{+\text{1.25}}_{-\text{1.60}}\times\text{10}^{\text{14}}$ M$_{\odot}$). The southern subcluster is also very massive (10.8$^{+\text{3.37}}_{-\text{3.54}}\times\text{10}^{\text{14}}$ M$_{\odot}$), yet it lacks an associated X-ray surface brightness peak, and it has been unidentified previously despite the detailed study of this \emph{Frontier Field} cluster. A secondary merger is occurring in the north along the line of sight with a third, less massive, subcluster (1.20$^{+\text{0.19}}_{-\text{0.34}}\times\text{10}^{\text{14}}$ M$_{\odot}$). We perform a Monte Carlo dynamical analysis on the main merger and estimate a collision speed at pericenter of 2770$^{+\text{610}}_{-\text{310}}$ km s$^{-\text{1}}$. We show the merger to be returning from apocenter with core passage occurring 1.16$^{+\text{0.50}}_{-\text{0.25}}$ Gyr before the observed state. We identify the line of sight merging subcluster in a strong lensing analysis in the literature and show that it is likely bound to MACS J1149 despite having reached an extreme collision velocity of $\sim$4000 km s$^{-\text{1}}$.

We present a tomographic cosmic shear study from the Deep Lens Survey (DLS), which, providing a limiting magnitude r_{lim}~27 (5 sigma), is designed as a pre-cursor Large Synoptic Survey Telescope (LSST) survey with an emphasis on depth. Using five tomographic redshift bins, we study their auto- and cross-correlations to constrain cosmological parameters. We use a luminosity-dependent nonlinear model to account for the astrophysical systematics originating from intrinsic alignments of galaxy shapes. We find that the cosmological leverage of the DLS is among the highest among existing >10 sq. deg cosmic shear surveys. Combining the DLS tomography with the 9-year results of the Wilkinson Microwave Anisotropy Probe (WMAP9) gives Omega_m=0.293_{-0.014}^{+0.012}, sigma_8=0.833_{-0.018}^{+0.011}, H_0=68.6_{-1.2}^{+1.4} km/s/Mpc, and Omega_b=0.0475+-0.0012 for LCDM, reducing the uncertainties of the WMAP9-only constraints by ~50%. When we do not assume flatness for LCDM, we obtain the curvature constraint Omega_k=-0.010_{-0.015}^{+0.013} from the DLS+WMAP9 combination, which however is not well constrained when WMAP9 is used alone. The dark energy equation of state parameter w is tightly constrained when Baryonic Acoustic Oscillation (BAO) data are added, yielding w=-1.02_{-0.09}^{+0.10} with the DLS+WMAP9+BAO joint probe. The addition of supernova constraints further tightens the parameter to w=-1.03+-0.03. Our joint constraints are fully consistent with the final Planck results and also the predictions of a LCDM universe.

The galaxy cluster RX J0603.3+4214 at z=0.225 is one of the rarest clusters boasting an extremely large (~2 Mpc) radio-relic. Because of the remarkable morphology of the relic, the cluster is nicknamed "Toothbrush Cluster". Although the cluster's underlying mass distribution is one of the critical pieces of information needed to reconstruct the merger scenario responsible for the puzzling radio-relic morphology, its proximity to the Galactic plane b~10 deg has imposed significant observational challenges. We present a high-resolution weak-lensing study of the cluster with Subaru/Suprime Cam and Hubble Space Telescope imaging data. Our mass reconstruction reveals that the cluster is comprised of complicated dark matter substructures closely tracing the galaxy distribution, however in contrast with the relatively simple binary X-ray morphology. Nevertheless, we find that the cluster mass is still dominated by the two most massive clumps aligned north-south with a ~3:1 mass ratio (M_{200}=6.29_{-1.62}^{+2.24} x 10^{14} Msun and 1.98_{-0.74}^{+1.24} x 10^{14} Msun for the northern and southern clumps, respectively). The southern mass peak is ~2' offset toward the south with respect to the corresponding X-ray peak, which has a "bullet"-like morphology pointing south. Comparison of the current weak-lensing result with the X-ray, galaxy, and radio-relic suggests that perhaps the dominant mechanism responsible for the observed relic may be a high-speed collision of the two most massive subclusters, although the peculiarity of the morphology necessitates involvement of additional sub-clusters. Careful numerical simulations should follow in order to obtain more complete understanding of the merger scenario utilizing all existing observations.

We describe a new test of photometric redshift performance given a spectroscopic redshift sample. This test complements the traditional comparison of redshift {\it differences} by testing whether the probability density functions $p(z)$ have the correct {\it width}. We test two photometric redshift codes, BPZ and EAZY, on each of two data sets and find that BPZ is consistently overconfident (the $p(z)$ are too narrow) while EAZY produces approximately the correct level of confidence. We show that this is because EAZY models the uncertainty in its spectral energy distribution templates, and that post-hoc smoothing of the BPZ $p(z)$ provides a reasonable substitute for detailed modeling of template uncertainties. Either remedy still leaves a small surplus of galaxies with spectroscopic redshift very far from the peaks. Thus, better modeling of low-probability tails will be needed for high-precision work such as dark energy constraints with the Large Synoptic Survey Telescope and other large surveys.

Diffuse radio emission in the form of radio halos and relics has been found in a number of merging galaxy clusters. These structures indicate that shock and turbulence associated with the merger accelerate electrons to relativistic energies. We report the discovery of a radio relic + radio halo system in PSZ1 G108.18-11.53 (z=0.335). This cluster hosts the second most powerful double radio relic system ever discovered. We observed PSZ1 G108.18-11.53 with the Giant Meterwave Radio Telescope (GMRT) and the Westerbork Synthesis Radio Telescope (WSRT). We obtained radio maps at 147, 323, 607 and 1380 MHz. We also observed the cluster with the Keck telescope, obtaining the spectroscopic redshift for 42 cluster members. From the injection index we obtained the Mach number of the shocks generating the two radio relics. For the southern shock we found M = 2.33^{+0.19}_{-0.26}, while the northern shock Mach number goes from M = 2.20^{+0.07}_{-0.14} in the north part down to M = 2.00^{+0.03}_{-0.08} in the southern region. If the relation between the injection index and the Mach number predicted by diffusive shock acceleration (DSA) theory holds, this is the first observational evidence for a gradient in the Mach number along a galaxy cluster merger shock.

Merging galaxy clusters with radio relics provide rare insights to the merger dynamics as the relics are created by the violent merger process. We demonstrate one of the first uses of the properties of the radio relic to reduce the uncertainties of the dynamical variables and determine the 3D configuration of a cluster merger, ACT-CL J0102-4915, nicknamed El Gordo. From the double radio relic observation and the X-ray observation of a comet-like gas morphology induced by motion of the cool core, it is widely believed that El Gordo is observed shortly after the first core-passage of the subclusters. We employ a Monte Carlo simulation to investigate the three-dimensional (3D) configuration and dynamics of El Gordo. Using the polarization fraction of the radio relic, we constrain the estimate of the angle between the plane of the sky and the merger axis to be $\alpha = 21~{\rm degree} \pm^9_{11}$. We find the relative 3D merger speed of El Gordo to be $2400\pm^{400}_{200}~{\rm km}~{\rm s}^{-1}$ at pericenter. The two possible estimates of the time-since-pericenter are $0.46\pm^{0.09}_{0.16}$ Gyr and $0.91\pm^{0.22}_{0.39}$ Gyr for the outgoing and returning scenario respectively. We put our estimates of the time-since-pericenter into context by showing that if the time-averaged shock velocity is approximately equal to or smaller than the pericenter velocity of the corresponding subcluster in the center of mass frame, the two subclusters are more likely to be moving towards, rather than away, from each other, post apocenter. We compare and contrast the merger scenario of El Gordo with that of the Bullet Cluster, and show that this late-stage merging scenario explains why the southeast dark matter lensing peak of El Gordo is closer to the merger center than the southeast cool core.

Cluster mergers may play a fundamental role in the formation and evolution of cluster galaxies. Stroe et al. (2014) revealed unexpected over-densities of candidate H$\alpha$ emitters near the ~1 Mpc-wide shock fronts of the massive (~2x10$^{15}$M$_{\odot}$) "Sausage" merging cluster, CIZA J2242.8+5301. We used Keck/DEIMOS and WHT/AF2 to confirm 83 H$\alpha$ emitters in and around the merging cluster. We find that cluster star-forming galaxies in the hottest X-ray gas and/or in the cluster sub-cores (away from the shock fronts) show high [SII]6716/[SII]6761 and high [SII]6716/H$\alpha$, implying very low electron densities (<30x lower than all other star-forming galaxies outside the cluster) and significant contribution from supernovae, respectively. All cluster star-forming galaxies near the cluster centre show evidence of significant outflows (blueshifted Na D~200-300km/s), likely driven by supernovae. Strong outflows are also found for the cluster H$\alpha$ AGN. H$\alpha$ star-forming galaxies in the merging cluster follow the z~0 mass-metallicity relation, showing systematically higher metallicity (~0.15-0.2 dex) than H$\alpha$ emitters outside the cluster (projected R>2.5 Mpc). This suggests that the shock front may have triggered remaining metal-rich gas which galaxies were able to retain into forming stars. Our observations show that the merger of impressively massive (~10$^{15}$M$_\odot$) clusters can provide the conditions for significant star-formation and AGN activity, but, as we witness strong feedback by star-forming galaxies and AGN (and given how massive the merging cluster is), such sources will likely quench in a few 100 Myrs.

CIZA J2242.8+5301 (`Sausage') and 1RXS J0603.3+4213 (`Toothbrush') are two low-redshift ($z\sim0.2$), massive ($\sim2\times10^{15}M_\odot$), post-core passage merging clusters, which host shock waves traced by diffuse radio emission. To study their star-formation properties, we uniformly survey the `Sausage' and `Toothbrush' clusters in broad and narrow band filters and select a sample of $201$ and $463$ line emitters, down to a rest-frame equivalent width ($13${\AA}). We robustly separate between H$\alpha$ and higher redshift emitters using a combination of optical multi-band (B, g, V, r, i, z) and spectroscopic data. We build H$\alpha$ luminosity functions for the entire cluster region, near the shock fronts, and away from the shock fronts and find striking differences between the two clusters. In the dynamically younger, $1$ Gyr old `Sausage' cluster we find numerous ($59$) H$\alpha$ emitters above a star-formation rate (SFR) of $0.17$ M$_{\sun}$ yr$^{-1}$ surprisingly located in close proximity to the shock fronts, embedded in very hot intra-cluster medium plasma. The SFR density for the cluster population is at least at the level of typical galaxies at $z\sim2$. Down to the same star-formation rate, the possibly dynamically more evolved `Toothbrush' cluster has only $9$ H$\alpha$ galaxies. The cluster H$\alpha$ galaxies fall on the SFR-stellar mass relation $z\sim0.2$ for the field. However, the `Sausage' cluster has an H$\alpha$ emitter density $>20$ times that of blank fields. If the shock passes through gas-rich cluster galaxies, the compressed gas could collapse into dense clouds and excite star-formation for a few $100$ Myr. This process ultimately leads to a rapid consumption of the molecular gas, accelerating the transformation of gas-rich field spirals into cluster S0s or ellipticals.

CIZA J2242.8+5301 is a merging system with a prominent (~2 Mpc long) radio relic, which together with the morphology of the X-ray emission provides strong evidence for a violent collision along the N-S axis. We present our constraints on the dark matter distribution of this unusual system using Subaru and CFHT imaging data. Measuring a high S/N lensing signal from this cluster is potentially a challenging task because of its proximity to the Milky Way plane (|b|~5 deg). We overcome this challenge with careful observation planning and systematics control, which enables us to successfully map the dark matter distribution of the cluster with high fidelity. The resulting mass map shows that the mass distribution is highly elongated along the N-S merger axis inferred from the orientation of the radio relics. Based on our mass reconstruction, we identify two sub-clusters, which coincide with the cluster galaxy distributions. We determine their masses using MCMC analysis by simultaneously fitting two NFW halos without fixing their centroids. The resulting masses of the northern and southern systems are $M_{200}=11.0_{-3.2}^{+3.7}\times10^{14} M_{\sun}$ and $9.8_{-2.5}^{+3.8}\times10^{14} M_{\sun}$, respectively, indicating that we are witnessing a post-collision of two giant systems of nearly equal mass. When the mass and galaxy centroids are compared in detail, we detect ~ 1' (~190 kpc) offsets in both northern and southern sub-clusters. We find that the galaxy luminosity-mass offset for the northern clump is statistically significant at the ~2 sigma level whereas the detection is only marginal for the southern sub-cluster in part because of a relatively large mass centroid error. We conclude that it is yet premature to uniquely attribute the galaxy-mass misalignment to SIDM and discuss caveats.

X-ray and radio observations of CIZA J2242.8+5301 suggest that it is a major cluster merger. Despite being well studied in the X-ray, and radio, little has been presented on the cluster structure and dynamics inferred from its galaxy population. We carried out a deep (i<25) broad band imaging survey of the system with Subaru SuprimeCam (g & i bands) and the Canada France Hawaii Telescope (r band) as well as a comprehensive spectroscopic survey of the cluster area (505 redshifts) using Keck DEIMOS. We use this data to perform a comprehensive galaxy/redshift analysis of the system, which is the first step to a proper understanding the geometry and dynamics of the merger, as well as using the merger to constrain self-interacting dark matter. We find that the system is dominated by two subclusters of comparable richness with a projected separation of 6.9'^{+0.7}_{-0.5} (1.3^{+0.13}_{-0.10} Mpc). We find that the north and south subclusters have similar redshifts of z=0.188 with a relative line-of-sight velocity difference of 69+/-190 km/s. We also find that north and south subclusters have velocity dispersions of 1160^{+100}_{-90} km/s and 1080^{+100}_{-70} km/s, respectively. These correspond to masses of 16.1^{+4.6}_{-3.3}x10^14 M_sun and 13.0^{+4.0}_{-2.5}x10^14 M_sun, respectively. While velocity dispersion measurements of merging clusters can be biased we believe the bias in this system to be minor due to the large projected separation and nearly plane-of-sky merger configuration. CIZA J2242.8+5301 is a relatively clean dissociative cluster merger with near 1:1 mass ratio, which makes it an ideal merger for studying merger associated physical phenomena.

We present the first sample of 882 optically selected galaxy clusters in the Deep Lens Survey (DLS), selected with the Bayesian Cluster Finder. We create mock DLS data to assess completeness and purity rates, and find that both are at least $70\%$ within 0.1$\le z \le$ 1.2 for clusters with $M_{200}\ge 1.2\times 10^{14}M_{\odot}$. We verified the integrity of the sample by performing several comparisons with other optical, weak lensing, X-ray and spectroscopic surveys which overlap the DLS footprint: the estimated redshifts are consistent with the spectroscopic redshifts of known clusters (for $z>0.25$ where saturation in the DLS is not an issue); our richness estimates in combination with a previously calibrated richness-mass relation yields individual cluster mass estimates consistent with available SHeLS dynamical mass estimates; synthetic mass maps made from the optical mass estimates are correlated ($>3\sigma$ significance) with the weak lensing mass maps; and the mass function thus derived is consistent with theoretical predictions for the CDM scenario. With the verified sample we investigated correlations between the brightest cluster galaxies (BCG) properties and the host cluster properties within a broader range in redshift (0.25 $\le z \le$ 0.8) and mass ($\ge2.4\times 10^{14}M_{\odot}$) than in previous work. We find that the slope of the BCG magnitude-redshift relation throughout this redshift range is consistent with that found at lower redshifts. This result supports an extrapolation to higher redshift of passive evolution of the BCG within the hierarchical scenario.

The galaxy cluster Abell 781 West has been viewed as a challenge to weak gravitational lensing mass calibration, as Cook and dell'Antonio (2012) found that the weak lensing signal-to-noise in three independent sets of observations was consistently lower than expected from mass models based on X-ray and dynamical measurements. We correct some errors in statistical inference in Cook and dell'Antonio (2012) and show that their own results agree well with the dynamical mass and exhibit at most 2.2--2.9$\sigma$ low compared to the X-ray mass, similar to the tension between the dynamical and X-ray masses. Replacing their simple magnitude cut with weights based on source photometric redshifts eliminates the tension between lensing and X-ray masses; in this case the weak lensing mass estimate is actually higher than, but still in agreement with, the dynamical estimate. A comparison of lensing analyses with and without photometric redshifts shows that a 1--2$\sigma$ chance alignment of low-redshift sources lowers the signal-to-noise observed by all previous studies which used magnitude cuts rather than photometric redshifts. The fluctuation is unexceptional, but appeared to be highly significant in Cook and dell'Antonio (2012) due to the errors in statistical interpretation.

We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz multi-band imaging survey of five 4 sq. degree fields with two National Optical Astronomy Observatory (NOAO) 4-meter telescopes at Kitt Peak and Cerro Tololo. For both telescopes, the change of the point-spread-function (PSF) shape across the focal plane is complicated, and the exposure-to-exposure variation of this position-dependent PSF change is significant. We overcome this challenge by modeling the PSF separately for individual exposures and CCDs with principal component analysis (PCA). We find that stacking these PSFs reproduces the final PSF pattern on the mosaic image with high fidelity, and the method successfully separates PSF-induced systematics from gravitational lensing effects. We calibrate our shears and estimate the errors, utilizing an image simulator, which generates sheared ground-based galaxy images from deep Hubble Space Telescope archival data with a realistic atmospheric turbulence model. For cosmological parameter constraints, we marginalize over shear calibration error, photometric redshift uncertainty, and the Hubble constant. We use cosmology-dependent covariances for the Markov Chain Monte Carlo analysis and find that the role of this varying covariance is critical in our parameter estimation. Our current non-tomographic analysis alone constrains the Omega_M-sigma_8 likelihood contour tightly, providing a joint constraint of Omega_M=0.262+-0.051 and sigma_8=0.868+-0.071. We expect that a future DLS weak-lensing tomographic study will further tighten these constraints because explicit treatment of the redshift dependence of cosmic shear more efficiently breaks the Omega_M-sigma_8 degeneracy. Combining the current results with the Wilkinson Microwave Anisotropy Probe 7-year (WMAP7) likelihood data, we obtain Omega_M=0.278+-0.018 and sigma_8=0.815+-0.020.

The NASA Exoplanet Science Institute (NExScI) hosts the annual Sagan Workshops, thematic meetings aimed at introducing researchers to the latest tools and methodologies in exoplanet research. The theme of the Summer 2012 workshop, held from July 23 to July 27 at Caltech, was to explore the use of exoplanet light curves to study planetary system architectures and atmospheres. A major part of the workshop was to use hands-on sessions to instruct attendees in the use of three open source tools for the analysis of light curves, especially from the Kepler mission. Each hands-on session involved the 160 attendees using their laptops to follow step-by-step tutorials given by experts. We describe how we used the Amazon Elastic Cloud 2 to run these applications.

Using about 450,000 galaxies in the Deep Lens Survey, we present a detection of the gravitational magnification of z > 4 Lyman Break Galaxies by massive foreground galaxies with 0.4 < z < 1.0, grouped by redshift. The magnification signal is detected at S/N greater than 20, and rigorous checks confirm that it is not contaminated by any galaxy sample overlap in redshift. The inferred galaxy mass profiles are consistent with earlier lensing analyses at lower redshift. We then explore the tomographic lens magnification signal by splitting our foreground galaxy sample into 7 redshift bins. Combining galaxy-magnification cross-correlations and galaxy angular auto-correlations, we develop a bias-independent estimator of the tomographic signal. As a diagnostic of magnification tomography, the measurement of this estimator rejects a flat dark matter dominated Universe at > 7.5{\sigma} with a fixed \sigma_8 and is found to be consistent with the expected redshift-dependence of the WMAP7 {\Lambda}CDM cosmology.

We examine the projected correlation of galaxies with mass from small scales (<few hundred kpc) where individual dark matter halos dominate, out to 15 Mpc where correlated large-scale structure dominates. We investigate these profiles as a function of galaxy luminosity and redshift. Selecting 0.8 million galaxies in the Deep Lens Survey, we use photometric redshifts and stacked weak gravitational lensing shear tomography out to radial scales of 1 degree from the centers of foreground galaxies. We detect correlated mass density from multiple halos and large-scale structure at radii larger than the virial radius, and find the first observational evidence for growth in the galaxy-mass correlation on 10 Mpc scales with decreasing redshift and fixed range of luminosity. For a fixed range of redshift, we find a scaling of projected halo mass with rest-frame luminosity similar to previous studies at lower redshift. We control systematic errors in shape measurement and photometric redshift, enforce volume completeness through absolute magnitude cuts, and explore residual sample selection effects via simulations.

We introduce a new method for constraining the redshift distribution of a set of galaxies, using weak gravitational lensing shear. Instead of using observed shears and redshifts to constrain cosmological parameters, we ask how well the shears around clusters can constrain the redshifts, assuming fixed cosmological parameters. This provides a check on photometric redshifts, independent of source spectral energy distribution properties and therefore free of confounding factors such as misidentification of spectral breaks. We find that ~40 massive ($\sigma_v=1200$ km/s) cluster lenses are sufficient to determine the fraction of sources in each of six coarse redshift bins to ~11%, given weak (20%) priors on the masses of the highest-redshift lenses, tight (5%) priors on the masses of the lowest-redshift lenses, and only modest (20-50%) priors on calibration and evolution effects. Additional massive lenses drive down uncertainties as $N_{lens}^0.5$, but the improvement slows as one is forced to use lenses further down the mass function. Future large surveys contain enough clusters to reach 1% precision in the bin fractions if the tight lens mass priors can be maintained for large samples of lenses. In practice this will be difficult to achieve, but the method may be valuable as a complement to other more precise methods because it is based on different physics and therefore has different systematic errors.

We present DLSCL J0916.2+2951 (z=0.53), a newly discovered major cluster merger in which the collisional cluster gas has become dissociated from the collisionless galaxies and dark matter. We identified the cluster using optical and weak lensing observations as part of the Deep Lens Survey. Our follow-up observations with Keck, Subaru, Hubble Space Telescope, and Chandra show that the cluster is a dissociative merger and constrain the dark matter self-interaction cross-section {\sigma}_{DM}m_{DM}^{-1}\leq7 cm^2g^{-1}. The system is observed at least 0.7\pm0.2 Gyr since first pass-through, thus providing a picture of cluster mergers 2-5 times further progressed than similar systems observed to date. This improved temporal leverage has implications for our understanding of merging clusters and their impact on galaxy evolution.

We describe the internal photometric calibration of the Deep Lens Survey, which consists of five widely separated fields observed by two different observatories. Adopting the global linear least-squares ("ubercal") approach developed for the Sloan Digital Sky Survey (SDSS), we derive flatfield corrections for all observing runs, which indicate that the original sky flats were nonuniform by up to 0.13 mag peak to valley in $\z$ band, and by up to half that amount in {\it BVR}. We show that application of these corrections reduces spatial nonuniformities in corrected exposures to the 0.01-0.02 mag level. We conclude with some lessons learned in applying ubercal to a survey structured very differently from SDSS, with isolated fields, multiple observatories, and shift-and-stare rather than drift-scan imaging. Although the size of the error caused by using sky or dome flats is instrument- and wavelength-dependent, users of wide-field cameras should not assume that it is small. Pipeline developers should facilitate routine application of this procedure, and surveys should include it in their plans from the outset.

GRavitational lEnsing Accuracy Testing 2010 (GREAT10) is a public image analysis challenge aimed at the development of algorithms to analyze astronomical images. Specifically, the challenge is to measure varying image distortions in the presence of a variable convolution kernel, pixelization and noise. This is the second in a series of challenges set to the astronomy, computer science and statistics communities, providing a structured environment in which methods can be improved and tested in preparation for planned astronomical surveys. GREAT10 extends upon previous work by introducing variable fields into the challenge. The "Galaxy Challenge" involves the precise measurement of galaxy shape distortions, quantified locally by two parameters called shear, in the presence of a known convolution kernel. Crucially, the convolution kernel and the simulated gravitational lensing shape distortion both now vary as a function of position within the images, as is the case for real data. In addition, we introduce the "Star Challenge" that concerns the reconstruction of a variable convolution kernel, similar to that in a typical astronomical observation. This document details the GREAT10 Challenge for potential participants. Continually updated information is also available from http://www.greatchallenges.info.

The detection of galaxy clusters in present and future surveys enables measuring mass-to-light ratios, clustering properties, galaxy cluster abundances and therefore, constraining cosmological parameters. We present a new technique for detecting galaxy clusters, which is based on the Matched Filter Algorithm from a Bayesian point of view. The method is able to determine the position, redshift and richness of the cluster through the maximization of a filter depending on galaxy luminosity, density and photometric redshift combined with a galaxy cluster prior that accounts for color-magnitude relations and BCG-redshift relation. We tested the algorithm through realistic mock galaxy catalogs, revealing that the detections are 100% complete and 80% pure for clusters up to z $<$1.2 and richer than $\Lambda_{CL}>$20 (Abell Richness $\sim$0, M$\sim4\times10^{14} M_{\odot}$). The completeness and purity remains approximately the same if we do not include the prior information, implying that this method is able to detect galaxy cluster with and without a well defined red sequence. We applied the algorithm to the CFHTLS Archive Research Survey (CARS) data, recovering similar detections as previously published using the same or deeper data plus additional clusters which appear to be real.

We present a sample of 17 newly discovered ultracool dwarf candidates later than ~M8, drawn from 231.90 arcmin2 of {\it Hubble Space Telescope} Wide Field Camera 3 infrared imaging. By comparing the observed number counts for 17.5<J_125<25.5 AB mag to an exponential disk model, we estimate a vertical scale height of z_scl=290 +- 25 (random) +- 30 (systematic) pc for a binarity fraction of f_b=0. While our estimate is roughly consistent with published results, we suggest that the differences can be attributed to sample properties, with the present sample containing far more substellar objects than previous work. We predict the object counts should peak at J_{125}~24 AB mag due to the exponentially-declining number density at the edge of the disc. We conclude by arguing that trend in scale height with spectral type may breakdown for brown dwarfs since they do not settle onto the main sequence.

The detection of galaxy clusters in present and future surveys enables measuring mass-to-light ratios, clustering properties or galaxy cluster abundances and therefore, constraining cosmological parameters. We present a new technique for detecting galaxy clusters, which is based on the Matched Filter Algorithm from a Bayesian point of view. The method is able to determine the position, redshift and richness of the cluster through the maximization of a filter depending on galaxy luminosity, density and photometric redshift combined with a galaxy cluster prior. We tested the algorithm through realistic mock galaxy catalogs, revealing that the detections are 100% complete and 80% pure for clusters up to z <1.2 and richer than \Lambda > 25 (Abell Richness > 0). We applied the algorithm to the CFHTLS Archive Research Survey (CARS) data, recovering similar detections as previously published using the same data plus additional clusters that are very probably real. We also applied this algorithm to the Deep Lens Survey (DLS), obtaining the first sample of optical-selected galaxy in this survey. The sample is complete up to redshift 0.7 and we detect more than 780 cluster candidates up to redshift 1.2. We conclude by discussing the differences between previous weak lensing detections in this survey and optical detections in both samples.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We present the results of the GREAT08 Challenge, a blind analysis challenge to infer weak gravitational lensing shear distortions from images. The primary goal was to stimulate new ideas by presenting the problem to researchers outside the shear measurement community. Six GREAT08 Team methods were presented at the launch of the Challenge and five additional groups submitted results during the 6 month competition. Participants analyzed 30 million simulated galaxies with a range in signal to noise ratio, point-spread function ellipticity, galaxy size, and galaxy type. The large quantity of simulations allowed shear measurement methods to be assessed at a level of accuracy suitable for currently planned future cosmic shear observations for the first time. Different methods perform well in different parts of simulation parameter space and come close to the target level of accuracy in several of these. A number of fresh ideas have emerged as a result of the Challenge including a re-examination of the process of combining information from different galaxies, which reduces the dependence on realistic galaxy modelling. The image simulations will become increasingly sophisticated in future GREAT challenges, meanwhile the GREAT08 simulations remain as a benchmark for additional developments in shear measurement algorithms.

We present weak lensing mass estimates of seven shear-selected galaxy cluster candidates from the Deep Lens Survey. The clusters were previously identified as mass peaks in convergence maps of 8.6 sq. deg of R band imaging, and followed up with X-ray and spectroscopic confirmation, spanning a redshift range 0.19 - 0.68. Most clusters contained multiple X-ray peaks, yielding 17 total mass concentrations. In this paper, we constrain the masses of these X-ray sources with weak lensing, using photometric redshifts from the full set of BVRz' imaging to properly weight background galaxies according to their lensing distance ratios. We fit both NFW and singular isothermal sphere profiles, and find that the results are insensitive to the assumed profile. We also show that the results do not depend significantly on the assumed prior on the position of the mass peak, but that this may become an issue in future larger samples. The inferred velocity dispersions for the extended X-ray sources range from 250-800 km/s, with the exception of one source for which no lensing signal was found. This work further establishes shear selection as a viable technique for finding clusters, but also highlights some unresolved issues such as determination of the mass profile center without biasing the mass estimate, and fully accounting for line-of-sight projections. A follow-up paper will examine the mass-X-ray scaling relations of these clusters.

We use simulations to demonstrate that photometric redshift "errors" can be greatly reduced by using the photometric redshift probability distribution p(z) rather than a one-point estimate such as the most likely redshift. In principle this involves tracking a large array of numbers rather than a single number for each galaxy. We introduce a very simple estimator that requires tracking only a single number for each galaxy, while retaining the systematic-error-reducing properties of using the full p(z) and requiring only very minor modifications to existing photometric redshift codes. We find that using this redshift estimator (or using the full p(z)) can substantially reduce systematics in dark energy parameter estimation from weak lensing, at no cost to the survey.

We present a regularized maximum likelihood weak lensing reconstruction of the Deep Lens Survey F2 field (4 deg^2). High signal-to-noise ratio peaks in our lensing significance map appear to be associated with possible projected filamentary structures. The largest apparent structure extends for over a degree in the field and has contributions from known optical clusters at three redshifts (z ~ 0.3, 0.43, 0.5). Noise in weak lensing reconstructions is known to potentially cause "false positives"; we use Monte Carlo techniques to estimate the contamination in our sample, and find that 10-25% of the peaks are expected to be false detections. For significant lensing peaks we estimate the total signal-to-noise ratio of detection using a method that accounts for pixel-to-pixel correlations in our reconstruction. We also report the detection of a candidate relative underdensity in the F2 field with a total signal-to-noise ratio of ~ 5.5.

The GRavitational lEnsing Accuracy Testing 2008 (GREAT08) Challenge focuses on a problem that is of crucial importance for future observations in cosmology. The shapes of distant galaxies can be used to determine the properties of dark energy and the nature of gravity, because light from those galaxies is bent by gravity from the intervening dark matter. The observed galaxy images appear distorted, although only slightly, and their shapes must be precisely disentangled from the effects of pixelisation, convolution and noise. The worldwide gravitational lensing community has made significant progress in techniques to measure these distortions via the Shear TEsting Program (STEP). Via STEP, we have run challenges within our own community, and come to recognise that this particular image analysis problem is ideally matched to experts in statistical inference, inverse problems and computational learning. Thus, in order to continue the progress seen in recent years, we are seeking an infusion of new ideas from these communities. This document details the GREAT08 Challenge for potential participants. Please visit http://www.great08challenge.info for the latest information.

The Cold Dark Matter theory of gravitationally-driven hierarchical structure formation has earned its status as a paradigm by explaining the distribution of matter over large spans of cosmic distance and time. However, its central tenet, that most of the matter in the universe is dark and exotic, is still unproven; the dark matter hypothesis is sufficiently audacious as to continue to warrant a diverse battery of tests. While local searches for dark matter particles or their annihilation signals could prove the existence of the substance itself, studies of cosmological dark matter in situ are vital to fully understand its role in structure formation and evolution. We argue that gravitational lensing provides the cleanest and farthest-reaching probe of dark matter in the universe, which can be combined with other observational techniques to answer the most challenging and exciting questions that will drive the subject in the next decade: What is the distribution of mass on sub-galactic scales? How do galaxy disks form and bulges grow in dark matter halos? How accurate are CDM predictions of halo structure? Can we distinguish between a need for a new substance (dark matter) and a need for new physics (departures from General Relativity)? What is the dark matter made of anyway? We propose that the central tool in this program should be a wide-field optical imaging survey, whose true value is realized with support in the form of high-resolution, cadenced optical/infra-red imaging, and massive-throughput optical spectroscopy.

The coming decade will be an exciting period for dark energy research, during which astronomers will address the question of what drives the accelerated cosmic expansion as first revealed by type Ia supernova (SN) distances, and confirmed by later observations. The mystery of dark energy poses a challenge of such magnitude that, as stated by the Dark Energy Task Force (DETF), "nothing short of a revolution in our understanding of fundamental physics will be required to achieve a full understanding of the cosmic acceleration." The lack of multiple complementary precision observations is a major obstacle in developing lines of attack for dark energy theory. This lack is precisely what next-generation surveys will address via the powerful techniques of weak lensing (WL) and baryon acoustic oscillations (BAO) -- galaxy correlations more generally -- in addition to SNe, cluster counts, and other probes of geometry and growth of structure. Because of their unprecedented statistical power, these surveys demand an accurate understanding of the observables and tight control of systematics. This white paper highlights the opportunities, approaches, prospects, and challenges relevant to dark energy studies with wide-deep multiwavelength photometric redshift surveys. Quantitative predictions are presented for a 20000 sq. deg. ground-based 6-band (ugrizy) survey with 5-sigma depth of r~27.5, i.e., a Stage 4 survey as defined by the DETF.

Next generation probes of dark matter and dark energy require high precision reconstruction of faint galaxy shapes from hundreds of dithered exposures. Current practice is to stack the images. While valuable for many applications, this stack is a highly compressed version of the data. Future weak lensing studies will require analysis of the full dataset using the stack and its associated catalog only as a starting point. We describe a "Multi-Fit" algorithm which simultaneously fits individual galaxy exposures to a common profile model convolved with each exposure's point spread function at that position in the image. This technique leads to an enhancement of the number of usable small galaxies at high redshift and, more significantly, a decrease in systematic shear error.

We compare X-ray and weak-lensing masses for four galaxy clusters that comprise the top-ranked shear-selected cluster system in the Deep Lens Survey. The weak-lensing observations of this system, which is associated with A781, are from the Kitt Peak Mayall 4-m telescope, and the X-ray observations are from both Chandra and XMM-Newton. For a faithful comparison of masses, we adopt the same matter density profile for each method, which we choose to be an NFW profile. Since neither the X-ray nor weak-lensing data are deep enough to well constrain both the NFW scale radius and central density, we estimate the scale radius using a fitting function for the concentration derived from cosmological hydrodynamic simulations and an X-ray estimate of the mass assuming isothermality. We keep this scale radius in common for both X-ray and weak-lensing profiles, and fit for the central density, which scales linearly with mass. We find that for three of these clusters, there is agreement between X-ray and weak-lensing NFW central densities, and thus masses. For the other cluster, the X-ray central density is higher than that from weak-lensing by 2 sigma. X-ray images suggest that this cluster may be undergoing a merger with a smaller cluster. This work serves as an additional step towards understanding the possible biases in X-ray and weak-lensing cluster mass estimation methods. Such understanding is vital to efforts to constrain cosmology using X-ray or weak-lensing cluster surveys to trace the growth of structure over cosmic time.

We examine the impact of non-Gaussian photometry errors on photometric redshift performance. We find that they greatly increase the scatter, but this can be mitigated to some extent by incorporating the correct noise model into the photometric redshift estimation process. However, the remaining scatter is still equivalent to that of a much shallower survey with Gaussian photometry errors. We also estimate the impact of non-Gaussian errors on the spectroscopic sample size required to verify the photometric redshift rms scatter to a given precision. Even with Gaussian {\it photometry} errors, photometric redshift errors are sufficiently non-Gaussian to require an order of magnitude larger sample than simple Gaussian statistics would indicate. The requirements increase from this baseline if non-Gaussian photometry errors are included. Again the impact can be mitigated by incorporating the correct noise model, but only to the equivalent of a survey with much larger Gaussian photometry errors. However, these requirements may well be overestimates because they are based on a need to know the rms, which is particularly sensitive to tails. Other parametrizations of the distribution may require smaller samples.

We investigate the impact of photometric signal-to-noise (S/N) on the precision of photometric redshifts in multi-band imaging surveys, using both simulations and real data. We simulate the optical 4-band (BVRz) Deep Lens Survey (DLS, Wittman etal 2002), and use the publicly available Bayesian Photometric Redshift code BPZ by Benitez (2000). The simulations include a realistic range of magnitudes and colors and vary from infinite S/N to S/N=5. The real data are from DLS photometry and two spectroscopic surveys, and explore a range of S/N by adding noise to initially very high S/N photometry. Precision degrades steadily as S/N drops, both because of direct S/N effects and because lower S/N is linked to fainter galaxies with a weaker magnitude prior. If a simple S/N cut were used, S/N>17 in R (corresponding, in the DLS, to lower S/N in other bands) would be required to keep the scatter in Deltaz = (zspec-zphot)/(1+zspec) to less than 0.1. However, cutting on ODDS (a measure of the peakiness of the probability density function provided by BPZ) greater than 0.4 provides roughly double the number of usable galaxies with the same scatter. Ellipticals form the tightest zspec-zphot relation, and cutting on type=elliptical provides better precision than the ODDS>0.9 cut, but this eliminates the vast majority of galaxies in a deep survey. In addition to being more efficient than a type cut, ODDS also has the advantages of working with all types of galaxies (although ellipticals are overrepresented) and of being a continuous parameter for which the severity of the cut can be adjusted as desired.

We present the first sample of galaxy clusters selected on the basis of their weak gravitational lensing shear. The shear induced by a cluster is a function of its mass profile and its redshift relative to the background galaxies being sheared; in contrast to more traditional methods of selecting clusters, shear selection does not depend on the cluster's star formation history, baryon content, or dynamical state. Because mass is the property of clusters which provides constraints on cosmological parameters, the dependence on these other parameters could induce potentially important biases in traditionally-selected samples. Comparison of a shear-selected sample with optically and X-ray selected samples is therefore of great importance. Here we present the first step toward a new shear-selected sample: the selection of cluster candidates from the first 8.6 deg$^2$ of the 20 deg$^2$ Deep Lens Survey (DLS), and tabulation of their basic properties such as redshifts and optical and X-ray counterparts.

We report the serendipitous discovery of two strong gravitational lens candidates (ACS J160919+6532 and ACS J160910+6532) in deep images obtained with the Advanced Camera for Surveys on the Hubble Space Telescope, each less than 40 arcsec from the previously known gravitational lens system CLASS B1608+656. The redshifts of both lens galaxies have been measured with Keck and Gemini: one is a member of a small galaxy group at z~0.63, which also includes the lensing galaxy in the B1608+656 system, and the second is a member of a foreground group at z~0.43. By measuring the effective radii and surface brightnesses of the two lens galaxies, we infer their velocity dispersions based on the passively evolving Fundamental Plane (FP) relation. Elliptical isothermal lens mass models are able to explain their image configurations within the lens hypothesis, with a velocity dispersion compatible with that estimated from the FP for a reasonable source-redshift range. Based on the large number of massive early-type galaxies in the field and the number-density of faint blue galaxies, the presence of two additional lens systems around CLASS B1608+656 is not unlikely in hindsight. Gravitational lens galaxies are predominantly early-type galaxies, which are clustered, and the lensed quasar host galaxies are also clustered. Therefore, obtaining deep high-resolution images of the fields around known strong lens systems is an excellent method of enhancing the probability of finding additional strong gravitational lens systems.

The Shear TEsting Programme, STEP, is a collaborative project to improve the accuracy and reliability of all weak lensing measurements in preparation for the next generation of wide-field surveys. In this first STEP paper we present the results of a blind analysis of simulated ground-based observations of relatively simple galaxy morphologies. The most successful methods are shown to achieve percent level accuracy. From the cosmic shear pipelines that have been used to constrain cosmology, we find weak lensing shear measured to an accuracy that is within the statistical errors of current weak lensing analyses, with shear measurements accurate to better than 7%. The dominant source of measurement error is shown to arise from calibration uncertainties where the measured shear is over or under-estimated by a constant multiplicative factor. This is of concern as calibration errors cannot be detected through standard diagnostic tests. The measured calibration errors appear to result from stellar contamination, false object detection, the shear measurement method itself, selection bias and/or the use of biased weights. Additive systematics (false detections of shear) resulting from residual point-spread function anisotropy are, in most cases, reduced to below an equivalent shear of 0.001, an order of magnitude below cosmic shear distortions on the scales probed by current surveys. Our results provide a snapshot view of the accuracy of current ground-based weak lensing methods and a benchmark upon which we can improve. To this end we provide descriptions of each method tested and include details of the eight different implementations of the commonly used Kaiser, Squires and Broadhurst (1995) method (KSB+) to aid the improvement of future KSB+ analyses.

The Smithsonian Hectospec Lensing Survey (SHELS) combines a large deep complete redshift survey with a weak lensing map from the Deep Lens Survey (Wittman et al. 2002; 2005). We use maps of the velocity dispersion based on systems identified in the redshift survey to compare the three-dimensional matter distribution with the two-dimensional projection mapped by weak lensing. We demonstrate directly that the lensing map images the three-dimensional matter distribution obtained from the kinematic data.

Weak lensing observations have the potential to be even more powerful than cosmic microwave background (CMB) observations in constraining cosmological parameters. However, the practical limits to weak lensing observations are not known. Most theoretical studies of weak lensing constraints on cosmology assume that the only limits are shot noise on small scales, and cosmic variance on large scales. For future large surveys, shot noise will be so low that other, systematic errors will likely dominate. Here we examine a potential source of additive systematic error for ground-based observations: spurious power induced by the atmosphere. We show that this limit will not be a significant factor even in future massive surveys such as LSST.

We present an analysis of weak lensing observations for RXJ1347-1145 over a 43' X 43' field taken in B and R filters on the Blanco 4m telescope at CTIO. RXJ1347-1145 is a massive cluster at redshift z=0.45. Using a population of galaxies with 20<R<26, we detect a weak lensing signal at the p<0.0005 level, finding best-fit parameters of \sigma_v=1400^{+130}_{-140} km s^{-1} for a singular isothermal sphere model and r_{200} = 3.5^{+0.8}_{-0.2} Mpc with c = 15^{+64}_{-10} for a NFW model in an \Omega_m = 0.3, \Omega_\Lambda = 0.7 cosmology. In addition, a mass to light ratio M/L_R =90 \pm 20 M_\odot / L_{R\odot} was determined. These values are consistent with the previous weak lensing study of RXJ1347--1145 by Fischer and Tyson, 1997, giving strong evidence that systemic bias was not introduced by the relatively small field of view in that study. Our best-fit parameter values are also consistent with recent X-ray studies by Allen et al, 2002 and Ettori et al, 2001, but are not consistent with recent optical velocity dispersion measurements by Cohen and Kneib, 2002.

We present a weak lensing analysis of the high-redshift cluster Cl 1604+4304. At z=0.90, this is the highest-redshift cluster yet detected with weak lensing. It is also one of a sample of high-redshift, optically-selected clusters whose X-ray temperatures are lower than expected based on their velocity dispersions. Both the gas temperature and galaxy velocity dispersion are proxies for its mass, which can be determined more directly by a lensing analysis. Modeling the cluster as a singular isothermal sphere, we find that the mass contained within projected radius R is 3.69+-1.47 * (R/500 kpc) 10^14 M_odot. This corresponds to an inferred velocity dispersion of 1004+-199 km/s, which agrees well with the measured velocity dispersion of 989+98-76 km/s (Gal & Lubin 2004). These numbers are higher than the 575+110-85 km/s inferred from Cl 1604+4304 X-ray temperature, however all three velocity dispersion estimates are consistent within ~ 1.9 sigma.

We report on the methodology and first results from the Deep Lens Survey transient search. We utilize image subtraction on survey data to yield all sources of optical variability down to 24th magnitude. Images are analyzed immediately after acquisition, at the telescope and in near-real time, to allow for followup in the case of time-critical events. All classes of transients are posted to the web upon detection. Our observing strategy allows sensitivity to variability over several decades in timescale. The DLS is the first survey to classify and report all types of photometric and astrometric variability detected, including solar system objects, variable stars, supernovae, and short timescale phenomena. Three unusual optical transient events were detected, flaring on thousand-second timescales. All three events were seen in the B passband, suggesting blue color indices for the phenomena. One event (OT 20020115) is determined to be from a flaring Galactic dwarf star of spectral type dM4. From the remaining two events, we find an overall rate of \eta = 1.4 events deg-2 day-1 on thousand-second timescales, with a 95% confidence limit of \eta < 4.3. One of these events (OT 20010326) originated from a compact precursor in the field of galaxy cluster Abell 1836, and its nature is uncertain. For the second (OT 20030305) we find strong evidence for an extended extragalactic host. A dearth of such events in the R passband yields an upper 95% confidence limit on short timescale astronomical variability between 19.5 < R < 23.4 of \eta_R < 5.2. We report also on our ensemble of astrometrically variable objects, as well as an example of photometric variability with an undetected precursor.

We report on preliminary results of our X-ray survey of the most massive clusters currently identified from the Deep Lens Survey (DLS). The DLS cluster sample is selected based on weak lensing shear, which makes it possible for the first time to study clusters in a baryon-independent way. In this article we present X-ray properties of a subset of the shear-selected cluster sample.

We report the weak lensing discovery, spectroscopic confirmation, and weak lensing tomography of a massive cluster of galaxies at $z=0.68$, demonstrating that shear selection of clusters works at redshifts high enough to be cosmologically interesting. The mass estimate from weak lensing, $11.1 +- 2.8 x 10^{14} (r/Mpc)$ solar masses within projected radius r, agrees with that derived from the spectroscopy ($\sigma_v = 980 km s^{-1}$), and with the position of an arc which is likely to be a strongly lensed background galaxy. The redshift estimate from weak lensing tomography is consistent with the spectroscopy, demonstrating the feasibility of baryon-unbiased mass surveys. This tomographic technique will be able to roughly identify the redshifts of any dark clusters which may appear in shear-selected samples, up to z ~ 1.

We measure seeing-corrected ellipticities for 2 x 10^6 galaxies with magnitude R<23 in 12 widely separated fields totalling 75 deg^2 of sky. At angular scales >30\arcmin, ellipticity correlations are detected at high significance and exhibit nearly the pure "E-mode" behavior expected of weak gravitational lensing. Even when smoothed to the full field size of 2.5 degrees, which is ~25h^-1 Mpc at the lens distances, an rms shear variance of <\gamma^2>^1/2 = 0.0012 +- 0.0003 is detected. At smaller angular scales there is significant "B-mode" power, an indication of residual uncorrected PSF distortions. The >30\arcmin data constrain the power spectrum of matter fluctuations on comoving scales of ~10h^-1 Mpc to have \sigma_8 (\Omega_m/0.3)^{0.57} = 0.71^{+0.12}_{-0.16} (95% CL, \LambdaCDM, \Gamma=0.21), where the systematic error includes statistical and calibration uncertainties, cosmic variance, and a conservative estimate of systematic contamination based upon the detected B-mode signal. This normalization of the power spectrum is lower than previous weak-lensing results but generally consistent them, is at the lower end of the \sigma_8 range from various analyses of galaxy cluster abundances, and agrees with recent determinations from CMB and galaxy clustering. The large and dispersed sky coverage of our survey reduces random errors and cosmic variance, while the relatively shallow depth allows us to use existing redshift-survey data to reduce systematic uncertainties in the N(z) distribution to insignificance. Reanalysis of the data with more sophisticated algorithms will hopefully reduce the systematic (B-mode) contamination, and allow more precise, multidimensional constraint of cosmological parameters.

The Deep Lens Survey (DLS) is a deep BVRz' imaging survey of seven 2x2 degree fields, with all data to be made public. The primary scientific driver is weak gravitational lensing, but the survey is also designed to enable a wide array of other astrophysical investigations. A unique feature of this survey is the search for transient phenomena. We subtract multiple exposures of a field, detect differences, classify, and release transients on the Web within about an hour of observation. Here we summarize the scientific goals of the DLS, field and filter selection, observing techniques and current status, data reduction, data products and release, and transient detections. Finally, we discuss some lessons which might apply to future large surveys such as LSST.

The number of mass clusters and their distribution in redshift are very sensitive to the density of matter Omega_m and the equation of state of dark energy w. Using weak lens gravitational tomography one can detect clusters of dark matter, weigh them, image their projected mass distribution, and determine their 3-D location. The degeneracy curve in the Omega_m - w plane is nearly orthogonal to that from CMB or SN measurements. Thus, a combination of CMB data with weak lens tomography of clusters can yield precision measurements of Omega_m and w, independently of the SN observations. The Large Synoptic Survey Telescope (LSST) will repeatedly survey 30,000 square degrees of the sky in multiple wavelengths. LSST will create a 3-D tomographic assay of mass overdensities back to half the age of the universe by measuring the shear and color-redshift of billions of high redshift galaxies. By simultaneously measuring several functions of cosmic shear and mass cluster abundance, LSST will provide a number of independent constraints on the dark energy density and the equation of state.

Weak lensing review at the first-year graduate student level.

We report the discovery of a cluster of galaxies via its weak gravitational lensing effect on background galaxies, the first spectroscopically confirmed cluster to be discovered through its gravitational effects rather than by its electromagnetic radiation. This fundamentally different selection mechanism promises to yield mass-selected, rather than baryon or photon-selected, samples of these important cosmological probes. We have confirmed this cluster with spectroscopic redshifts of fifteen members at z=0.276, with a velocity dispersion of 615 km/s. We use the tangential shear as a function of source photometric redshift to estimate the lens redshift independently and find z_l = 0.30 +- 0.08. The good agreement with the spectroscopy indicates that the redshift evolution of the mass function may be measurable from the imaging data alone in shear-selected surveys.

We present a survey of mass profiles and mass-to-light ratios of eight typical galaxy clusters at a common redshift (z ~ 0.2). We use weak gravitational lensing as a probe because it is unique in avoiding any assumptions about the dynamical state of the clusters. To avoid bias toward the rare and spectacular clusters that are easy targets for lensing work, we selected an ensemble of much more common clusters with moderate X-ray luminosity. Although the survey is still in progress, two conclusions are emerging: (1) within a cluster, mass follows light very closely on the angular scales that we can measure, $0.2-2h^{-1}$ Mpc, and (2) there is a significant cluster-to-cluster scatter in mass-to-light (M/L) ratios despite uniformity of observing, reduction, and analysis procedures. We also derive an estimate of $\Omega_{matter}$ based on extrapolation from the mass properties of these typical clusters. Finally, we discuss the discovery of other clusters in our fields through their lensing signal.

Weak gravitational lensing enables direct reconstruction of dark matter maps over cosmologically significant volumes. This research is currently telescope-limited. The Dark Matter Telescope (DMT) is a proposed 8.4 m telescope with a 3 degree field of view, with an etendue of 260 $(m. degree)^2$, ten times greater than any other current or planned telescope. With its large etendue and dedicated observational mode, the DMT fills a nearly unexplored region of parameter space and enables projects that would take decades on current facilities. The DMT will be able to reach 10-sigma limiting magnitudes of 27-28 magnitude in the wavelength range .3 - 1 um over a 7 square degree field in 3 nights of dark time. Here we review its unique weak lensing cosmology capabilities and the design that enables those capabilities.

Most of the matter in the universe is not luminous and can be observed directly only through its gravitational effect. An emerging technique called weak gravitational lensing uses background galaxies to reveal the foreground dark matter distribution on large scales. Light from very distant galaxies travels to us through many intervening overdensities which gravitationally distort their apparent shapes. The observed ellipticity pattern of these distant galaxies thus encodes information about the large-scale structure of the universe, but attempts to measure this effect have been inconclusive due to systematic errors. We report the first detection of this ``cosmic shear'' using 145,000 background galaxies to reveal the dark matter distribution on angular scales up to half a degree in three separate lines of sight. The observed angular dependence of this effect is consistent with that predicted by two leading cosmological models, providing new and independent support for these models.

We determine the luminosity function (LF) of galaxies in the core of the Coma cluster for M_R<=-11.4 (assuming H_0=75 km/s/Mpc), a magnitude regime previously explored only in the Local Group. Objects are counted in a deep CCD image of Coma having RMS noise of 27.7 R mag~arcsec$^{-2}$. A correction for objects in the foreground or background of the Coma cluster---and the uncertainty in this correction---are determined from images of five other high-latitude fields, carefully matched to the Coma image in both resolution and noise level. Accurate counts of Coma cluster members are obtained as faint as R=25.5, or M_R=-9.4. The LF for galaxies is well fit by a power law dN/dL\propto L^\alpha, with \alpha=-1.42\pm0.05, over the range -19.4<M_R<-11.4; faintward of this range, the galaxies are unresolved and indistinguishable from globular clusters, but the data are consistent with an extrapolation of the power law. Surface brightness biases are minimized since galaxies are not subjected to morphological selection, and the limiting detection isophote is 27.6 R mag~arcsec$^{-2}$. We find the typical $M_R\approx-12$ Coma cluster galaxy to have an exponential scale length $\approx200$~pc, similar to Local Group galaxies of comparable magnitude. These extreme dwarf galaxies show a surface density increasing towards the giant elliptical NGC~4874 as $r^{-1.3}$, similar to the diffuse light and globular cluster distributions. The luminosity in the detected dwarf galaxies is at most a few percent of the total diffuse light of the giant galaxies in the cluster, and the contribution of the dwarfs to the mass of the cluster is likely negligible as well.