Microlensing surveys have discovered thousands of events with almost all events discovered within the Galactic bulge or toward the Magellanic clouds. The Zwicky Transient Facility (ZTF), while not designed to be a microlensing campaign, is an optical time-domain survey that observes the entire northern sky every few nights including the Galactic plane. ZTF observes $\sim10^9$ stars in g-band and r-band and can significantly contribute to the observed microlensing population. We predict that ZTF will observe $\sim$1100 microlensing events in three years of observing within $10^\circ$ degrees latitude of the Galactic plane, with $\sim$500 events in the outer Galaxy ($\ell \geq 10^\circ$). This yield increases to $\sim$1400 ($\sim$800) events by combining every three ZTF exposures, $\sim$1800 ($\sim$900) events if ZTF observes for a total of five years, and $\sim$2400 ($\sim$1300) events for a five year survey with post-processing image stacking. Using the microlensing modeling software PopSyCLE, we compare the microlensing populations in the Galactic bulge and the outer Galaxy. We also present an analysis of the microlensing event ZTF18abhxjmj to demonstrate how to leverage these population statistics in event modeling. ZTF will constrain Galactic structure, stellar populations, and primordial black holes through photometric microlensing.

Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $\theta_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $\theta_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $\theta_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties.

We present a new Milky Way microlensing simulation code, dubbed PopSyCLE (Population Synthesis for Compact object Lensing Events). PopSyCLE is the first resolved microlensing simulation to include a compact object distribution derived from numerical supernovae explosion models and both astrometric and photometric microlensing effects. We demonstrate the capabilities of PopSyCLE by investigating the optimal way to find black holes (BHs) with microlensing. Candidate BHs have typically been selected from wide-field photometric microlensing surveys, such as OGLE, by selecting events with long Einstein crossing times ($t_E>120$ days). These events can be selected at closest approach and monitored astrometrically in order to constrain the mass of each lens; PopSyCLE predicts a BH detection rate of ~40% for such a program. We find that the detection rate can be enhanced to ~85% by selecting events with both $t_E>120$ days and a microlensing parallax of $\pi_E<0.08$. Unfortunately, such a selection criterion cannot be applied during the event as $\pi_E$ requires both pre- and post-peak photometry. However, historical microlensing events from photometric surveys can be revisited using this new selection criteria in order to statistically constrain the abundance of BHs in the Milky Way. The future WFIRST microlensing survey provides both precise photometry and astrometry and will yield individual masses of $\mathcal{O}(100-1000) black holes, which is at least an order of magnitude more than is possible with individual candidate follow-up with current facilities. The resulting sample of BH masses from WFIRST will begin to constrain the shape of the black hole present-day mass function, BH multiplicity, and BH kick velocity distributions.

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.

Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper highlights the scientific motivation, computational and algorithmic needs to build joint pixel level processing capabilities, which the individual projects by themselves will not be able to support. Through this white paper, we request that the Astro2020 decadal committee recognize the JSP effort as a multi-agency project with the natural outcome being a collaborative effort among groups which are normally supported by a single agency. JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also result in the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone.

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/ .

One of the key phases of stellar evolution that remains poorly understood is stellar death. We lack a predictive model for how a star of a given mass explodes and what kind of remnant it leaves behind (i.e. the initial-final mass relation, IFMR). Progress has been limited due to the difficulty in finding and weighing black holes and neutron stars in large numbers. Technological advances that allow for sub-milliarcsecond astrometry in crowded fields have opened a new window for finding black holes and neutron stars: astrometric gravitational lensing. Finding and weighing a sample of compact objects with astrometric microlensing will allow us to place some of the first constraints on the present-day mass function of isolated black holes and neutron stars, their multiplicity, and their kick velocities. All of these are fundamental inputs into understanding the death phase of stellar evolution, improving supernovae models, and interpreting LIGO detections in an astrophysical context. To achieve these goals, we require large area surveys, such as the WFIRST exoplanet microlensing survey, to photometrically identify long-duration (>120 day), un-blended microlensing events as candidate compact objects. We also require high-precision astrometric follow-up monitoring using extremely large telescopes, equipped with adaptive optics, such as TMT and GMT.

Gravitational microlensing is one of the few means of finding primordial black holes (PBHs), if they exist. Recent LIGO detections of 30 Msun black holes have re-invigorated the search for PBHs in the 10-100 Msun mass regime. Unfortunately, individual PBH microlensing events cannot easily be distinguished from stellar lensing events from photometry alone. However, the distribution of microlensing timescales (tE, the Einstein radius crossing time) can be analyzed in a statistical sense using models of the Milky Way with and without PBHs. While previous works have presented both theoretical models and observational constrains for PBHs (e.g. Calcino et al. 2018; Niikura et al. 2019), surprisingly, they rarely show the observed quantity -- the tE distribution -- for different abundances of PBHs relative to the total dark matter mass (fPBH). We present a simple calculation of how the tE distribution changes between models with and without PBHs.

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.

Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST).

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

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.

Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by RXJ0152.7-1357, for which we identify an additional multiple image, guided by a light-traces-mass approach we identify seven new sets of multiply imaged background sources lensed by this cluster, spanning the redshift range [1.79-3.93]. A total of 25 multiple images are seen over a small area of ~0.4 $arcmin^2$, allowing us to put relatively high-resolution constraints on the inner matter distribution. Although modestly massive, the high degree of substructure together with its very elongated shape make RXJ0152.7-1357 a very efficient lens for its size. This cluster also comprises the third-largest sample of z~6-7 candidates in the RELICS survey. Finally, we present a comparison of our resulting mass distribution and magnification estimates with those from a Lenstool model. These models are made publicly available through the MAST archive.

We present the results from new deep {\it Chandra} ($\sim410$~ks) and {\it Suzaku} ($\sim180$ ks) observations of the merging galaxy cluster ZwCl\,0008.8+5215 ($z=0.104$). Previous radio observations revealed the presence of a double radio relic located diametrically west and east of the cluster center. Using our new {\it Chandra} data, we find evidence for the presence of a shock at the location of the western relic, RW, with a Mach number $\mathcal{M}_{S_X}=1.48^{+0.50}_{-0.32}$ from the density jump. We also measure $\mathcal{M}_{T_X}=2.35^{+0.74}_{-0.55}$ and $\mathcal{M}_{T_X}=2.02^{+0.74}_{-0.47}$ from the temperature jump, with {\it Chandra} and {\it Suzaku} respectively. These values are consistent with the Mach number estimate from a previous study of the radio spectral index, under the assumption of diffusive shock acceleration ($\mathcal{M}_{\rm RW}=2.4^{+0.4}_{-0.2}$). Interestingly, the western radio relic does not entirely trace the X-ray shock. A possible explanation is that the relic traces fossil plasma from nearby radio galaxies which is re-accelerated at the shock. For the eastern relic we do not detect an X-ray surface brightness discontinuity, despite the fact that radio observations suggest a shock with $\mathcal{M}_{\rm RE}=2.2^{+0.2}_{-0.1}$. The low surface brightness and reduced integration time for this region might have prevented the detection. {\it Chandra} surface brightness profile suggests $\mathcal{M}\lesssim1.5$, while {\it Suzaku} temperature measurements found $\mathcal{M}_{T_X}=1.54^{+0.65}_{-0.47}$. Finally, we also detect a merger induced cold front on the western side of the cluster, behind the shock that traces the western relic.

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.

The overwhelming majority of objects visible to LSST lie within the Galactic Plane. Though many previous surveys have avoided this region for fear of stellar crowding, LSST's spatial resolution combined with its state-of-the-art Difference Image Analysis mean that it can conduct a high cadence survey of most of the Galaxy for the first time. Here we outline the many areas of science that would greatly benefit from an LSST survey that included the Galactic Plane, Magellanic Clouds and Bulge at a cadence of 2-3 d. Particular highlights include measuring the mass spectrum of black holes, and mapping the population of exoplanets in the Galaxy in relation to variations in star forming environments. But the same survey data will provide a goldmine for a wide range of science, and we explore possible survey strategies which maximize the scientific return for a number of fields including young stellar objects, cataclysmic variables and Neptune Trojans.

NASA's WFIRST mission will perform a wide-field, NIR survey of the Galactic Bulge to search for exoplanets via the microlensing techniques. As the mission is due to launch in the mid-2020s, around half-way through the LSST Main Survey, we have a unique opportunity to explore synergistic science from two landmark programs. LSST can survey the entire footprint of the WFIRST microlensing survey in a single Deep Drilling Field. Here we explore the great scientific potential of this proposal and recommend the most effective observing strategies.

Despite progress in understanding radio relics, there are still open questions regarding the underlying particle acceleration mechanisms. In this paper we present deep 1--4 GHz VLA observations of CIZA\,J2242.8+5301 ($z=0.1921$), a double radio relic cluster characterized by small projection on the plane of the sky. Our VLA observations reveal, for the first time, the complex morphology of the diffuse sources and the filamentary structure of the northern relic. We discover new faint diffuse radio emission extending north of the main northern relic. Our Mach number estimates for the northern and southern relics, based on the radio spectral index map obtained using the VLA observations and existing LOFAR and GMRT data, are consistent with previous radio and X-ray studies ($\mathcal{M}_{\rm RN}=2.58\pm0.17$ and $\mathcal{M}_{\rm RS}=2.10\pm0.08$). However, color-color diagrams and modelings suggest a flatter injection spectral index than the one obtained from the spectral index map, indicating that projection effects might be not entirely negligible. The southern relic consists of five "arms". Embedded in it, we find a tailed radio galaxy which seems to be connected to the relic. A spectral index flattening, where the radio tail connects to the relic, is also measured. We propose that the southern relic may trace AGN fossil electrons that are re-accelerated at a shock, with an estimated strength of $\mathcal{M}=2.4$. High-resolution mapping of other tailed radio galaxies also supports a scenario where AGN fossil electrons are revived by the merger event and could be related to the formation of some diffuse cluster radio emission.

We present a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the Reionization Lensing Cluster Survey. We use a Light-Traces-Mass technique to uncover sets of multiply images and constrain the mass distribution of the clusters. These mass models are the first published for Abell S295 and MACS J0159.8-0849, and are improvements over previous models for Abell 697 and MACS J0025.4-1222. Our analysis for MACS J0025.4-1222 and Abell S295 shows a bimodal mass distribution supporting the merger scenarios proposed for these clusters. The updated model for MACS J0025.4-1222 suggests a substantially smaller critical area than previously estimated. For MACS J0159.8-0849 and Abell 697 we find a single peak and relatively regular morphology, revealing fairly relaxed clusters. Despite being less prominent lenses, three of these clusters seem to have lensing strengths, i.e. cumulative area above certain magnification, similar to the Hubble Frontier Fields clusters (e.g., A($\mu>5$) $\sim 1-3$ arcmin$^2$, A($\mu>10$) $\sim 0.5-1.5$ arcmin$^2$), which in part can be attributed to their merging configurations. We make our lens models publicly available through the Mikulski Archive for Space Telescopes. Finally, using Gemini-N/GMOS spectroscopic observations we detect a single emission line from a high-redshift $J_{125}\simeq25.7$ galaxy candidate lensed by Abell 697. While we cannot rule out a lower-redshift solution, we interpret the line as Ly$\alpha$ at $z=5.800\pm 0.001$, in agreement with its photometric redshift and dropout nature. Within this scenario we measure a Ly$\alpha$ rest-frame equivalent width of $52\pm22$ \AA, and an observed Gaussian width of $117\pm 15$ km/s.

We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use as constraints for the model. We find a foreground structure at $z\sim0.4$, which we include as a second cluster-sized halo in one of our models; however two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within $r=26.7"$ (the region where the strong lensing constraints exist) to be $M=2.51^{+0.15}_{-0.09}\times 10^{14}$~M$_{\odot}$. If we extrapolate out to $r_{500}$, our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev-Zel'dovich effect ($M\sim10^{15}$~M$_{\odot}$). This cluster is lensing a previously reported $z\sim10$ galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known.

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.

Strong gravitational lensing by galaxy clusters has become a powerful tool for probing the high-redshift Universe, magnifying distant and faint background galaxies. Reliable strong lensing (SL) models are crucial for determining the intrinsic properties of distant, magnified sources and for constructing their luminosity function. We present here the first SL analysis of MACS J0308.9+2645 and PLCK G171.9-40.7, two massive galaxy clusters imaged with the Hubble Space Telescope in the framework of the Reionization Lensing Cluster Survey (RELICS). We use the Light-Traces-Mass modeling technique to uncover sets of multiply imaged galaxies and constrain the mass distribution of the clusters. Our SL analysis reveals that both clusters have particularly large Einstein radii ($\theta_E>30"$ for a source redshift of $z_s=2$), providing fairly large areas with high magnifications, useful for high-redshift galaxy searches ($\sim2$ arcmin$^{2}$ with $\mu>5$ to $\sim1$ arcmin$^{2}$ with $\mu>10$, similar to a typical \textit{Hubble Frontier Fields} cluster). We also find that MACS J0308.9+2645 hosts a promising, apparently bright (J$\sim23.2-24.6$ AB), multiply imaged high-redshift candidate at $z\sim6.4$. These images are amongst the brightest high-redshift candidates found in RELICS. Our mass models, including magnification maps, are made publicly available for the community through the Mikulski Archive for Space Telescopes.

We present LOFAR $120-168$ MHz images of the merging galaxy cluster Abell 1240 that hosts double radio relics. In combination with the GMRT $595-629$ MHz and VLA $2-4$ GHz data, we characterised the spectral and polarimetric properties of the radio emission. The spectral indices for the relics steepen from their outer edges towards the cluster centre and the electric field vectors are approximately perpendicular to the major axes of the relics. The results are consistent with the picture that these relics trace large-scale shocks propagating outwards during the merger. Assuming diffusive shock acceleration (DSA), we obtain shock Mach numbers of $\mathcal{M}=2.4$ and $2.3$ for the northern and southern shocks, respectively. For $\mathcal{M}\lesssim3$ shocks, a pre-existing population of mildly relativistic electrons is required to explain the brightness of the relics due to the high ($>10$ per cent) particle acceleration efficiency required. However, for $\mathcal{M}\gtrsim4$ shocks the required efficiency is $\gtrsim1\%$ and $\gtrsim0.5\%$, respectively, which is low enough for shock acceleration directly from the thermal pool. We used the fractional polarization to constrain the viewing angle to $\geqslant(53\pm3)^\circ$ and $\geqslant(39\pm5)^\circ$ for the northern and southern shocks, respectively. We found no evidence for diffuse emission in the cluster central region. If the halo spans the entire region between the relics ($\sim1.8\,\text{Mpc}$) our upper limit on the power is $P_\text{1.4 GHz}=(1.4\pm0.6)\times10^{23}\,\text{W}\text{Hz}^{-1}$ which is approximately equal to the anticipated flux from a cluster of this mass. However, if the halo is smaller than this, our constraints on the power imply that the halo is underluminous.

Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the z ~ 6 - 8 candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ~200 arcmin^2. These clusters were selected to be excellent lenses and we find similar high-redshift sample sizes and magnitude distributions as CLASH. We discover 321 candidate galaxies with photometric redshifts between z ~ 6 to z ~ 8, including extremely bright objects with H-band magnitudes of m_AB ~ 23 mag. As a sample, the observed (lensed) magnitudes of these galaxies are among the brightest known at z> 6, comparable to much wider, blank-field surveys. RELICS demonstrates the efficiency of using strong gravitational lenses to produce high-redshift samples in the epoch of reionization. These brightly observed galaxies are excellent targets for follow-up study with current and future observatories, including the James Webb Space Telescope.

Strong gravitational lensing by galaxy clusters magnifies background galaxies, enhancing our ability to discover statistically significant samples of galaxies at z>6, in order to constrain the high-redshift galaxy luminosity functions. Here, we present the first five lens models out of the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury Program, based on new HST WFC3/IR and ACS imaging of the clusters RXC J0142.9+4438, Abell 2537, Abell 2163, RXC J2211.7-0349, and ACT-CLJ0102-49151. The derived lensing magnification is essential for estimating the intrinsic properties of high-redshift galaxy candidates, and properly accounting for the survey volume. We report on new spectroscopic redshifts of multiply imaged lensed galaxies behind these clusters, which are used as constraints, and detail our strategy to reduce systematic uncertainties due to lack of spectroscopic information. In addition, we quantify the uncertainty on the lensing magnification due to statistical and systematic errors related to the lens modeling process, and find that in all but one cluster, the magnification is constrained to better than 20% in at least 80% of the field of view, including statistical and systematic uncertainties. The five clusters presented in this paper span the range of masses and redshifts of the clusters in the RELICS program. We find that they exhibit similar strong lensing efficiencies to the clusters targeted by the Hubble Frontier Fields within the WFC3/IR field of view. Outputs of the lens models are made available to the community through the Mikulski Archive for Space Telescopes

The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lensing. Discovered in the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury program and companion S-RELICS Spitzer program, this candidate has a lensed H-band magnitude of 25.7+/-0.1 AB mag. With a magnification of \mu~4-7 estimated from our lens models, the de-lensed intrinsic magnitude is 27.6+/-0.3 AB mag, and the half-light radius is r_e<0.8 kpc, both consistent with other z>9 candidates. The inferred stellar mass (log [M* /M_Sun]=9.7^{+0.7}_{-0.5}) and star formation rate (\log [SFR/M_Sun yr^{-1}]=1.3^{+0.2}_{-0.3}) indicate that this candidate is a typical star-forming galaxy on the z>6 SFR-M* relation. We note that three independent lens models predict two counterimages, at least one of which should be of a similar magnitude to the arc, but these counterimages are not yet detected. Counterimages would not be expected if the arc were at lower redshift. However, the only spectral energy distributions capable of fitting the Hubble and Spitzer photometry well at lower redshifts require unphysical combinations of z~2 galaxy properties. The unprecedented lensed size of this z~10 candidate offers the potential for the James Webb Space Telescope to study the geometric and kinematic properties of a galaxy observed 500 Myr after the Big Bang.

We report L-band VLA observations of 1RXS J0603.3+4214, a cluster that hosts a bright radio relic, known as the Toothbrush, and an elongated giant radio halo. These new observations allow us to study the surface brightness distribution down to one arcsec resolution with very high sensitivity. Our images provide an unprecedented detailed view of the Toothbrush, revealing enigmatic filamentary structures. To study the spectral index distribution, we complement our analysis with published LOFAR and GMRT observations. The bright `brush' of the Toothbrush shows a prominent narrow ridge to its north with a sharp outer edge. The spectral index at the ridge is in the range $-0.70\leq\alpha\leq-0.80$. We suggest that the ridge is caused by projection along the line of sight. With a simple toy model for the smallest region of the ridge, we conclude that the magnetic field is below $5\,\rm\mu G$ and varies significantly across the shock front. Our model indicates that the actual Mach number is higher than that obtained from the injection index and agrees well with the one derived from the overall spectrum, namely ${\cal M}=3.78^{+0.3}_{-0.2}$. The radio halo shows an average spectral index of $\alpha=-1.16\pm0.05$ and a slight gradient from north to south. The southernmost part of the halo is steeper and possibly related to a shock front. Excluding the southernmost part, the halo morphology agrees very well with the X-ray morphology. A power-law correlation is found between the radio and X-ray surface brightness

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.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

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.

We examine the effects of an impending cluster merger on galaxies in the large scale structure (LSS) RX J0910 at $z =1.105$. Using multi-wavelength data, including 102 spectral members drawn from the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey and precise photometric redshifts, we calculate star formation rates and map the specific star formation rate density of the LSS galaxies. These analyses along with an investigation of the color-magnitude properties of LSS galaxies indicate lower levels of star formation activity in the region between the merging clusters relative to the outskirts of the system. We suggest that gravitational tidal forces due to the potential of the merging halos may be the physical mechanism responsible for the observed suppression of star formation in galaxies caught between the merging clusters.

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.

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.

To investigate the relationship between thermal and non-thermal components in merger galaxy clusters, we present deep JVLA and Chandra observations of the HST Frontier Fields cluster MACS J0717.5+3745. The Chandra image shows a complex merger event, with at least four components belonging to different merging subclusters. NW of the cluster, $\sim 0.7$ Mpc from the center, there is a ram-pressure-stripped core that appears to have traversed the densest parts of the cluster after entering the ICM from the direction of a galaxy filament to the SE. We detect a density discontinuity NNE of this core which we speculate is associated with a cold front. Our radio images reveal new details for the complex radio relic and radio halo in this cluster. In addition, we discover several new filamentary radio sources with sizes of 100-300 kpc. A few of these seem to be connected to the main radio relic, while others are either embedded within the radio halo or projected onto it. A narrow-angled-tailed (NAT) radio galaxy, a cluster member, is located at the center of the radio relic. The steep spectrum tails of this AGN leads into the large radio relic where the radio spectrum flattens again. This morphological connection between the NAT radio galaxy and relic provides evidence for re-acceleration (revival) of fossil electrons. The presence of hot $\gtrsim 20$ keV ICM gas detected by Chandra near the relic location provides additional support for this re-acceleration scenario.

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.

We infer gravitational lensing shear and convergence fields from galaxy ellipticity catalogs under a spatial process prior for the lensing potential. We demonstrate the performance of our algorithm with simulated Gaussian-distributed cosmological lensing shear maps and a reconstruction of the mass distribution of the merging galaxy cluster Abell 781 using galaxy ellipticities measured with the Deep Lens Survey. Given interim posterior samples of lensing shear or convergence fields on the sky, we describe an algorithm to infer cosmological parameters via lens field marginalization. In the most general formulation of our algorithm we make no assumptions about weak shear or Gaussian distributed shape noise or shears. Because we require solutions and matrix determinants of a linear system of dimension that scales with the number of galaxies, we expect our algorithm to require parallel high-performance computing resources for application to ongoing wide field lensing surveys.

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

We have developed a maximum likelihood source detection method capable of detecting ultra-faint streaks with surface brightnesses approximately an order of magnitude fainter than the pixel level noise. Our maximum likelihood detection method is a model based approach that requires no a priori knowledge about the streak location, orientation, length, or surface brightness. This method enables discovery of typically undiscovered objects, and enables the utilization of low-cost sensors (i.e., higher-noise data). The method also easily facilitates multi-epoch co-addition. We will present the results from the application of this method to simulations, as well as real low earth orbit observations.

We present a Bayesian algorithm to combine optical imaging of unresolved objects from distinct epochs and observation platforms for orbit determination and tracking. By propagating the non-Gaussian uncertainties we are able to optimally combine imaging of arbitrary signal-to-noise ratios, allowing the integration of data from low-cost sensors. Our Bayesian approach to image characterization also allows large compression of imaging data without loss of statistical information. With a computationally efficient algorithm to combine multiple observation epochs and multiple telescopes, we show statistically optimal orbit inferences.

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 report on high-resolution JVLA and Chandra observations of the HST Frontier Cluster MACS J0717.5+3745. MACS J0717.5+3745 offers the largest contiguous magnified area of any known cluster, making it a promising target to search for lensed radio and X-ray sources. With the high-resolution 1.0-6.5 GHz JVLA imaging in A and B configuration, we detect a total of 51 compact radio sources within the area covered by the HST imaging. Within this sample we find 7 lensed sources with amplification factors larger than $2$. None of these sources are identified as multiply-lensed. Based on the radio luminosities, the majority of these sources are likely star forming galaxies with star formation rates of 10-50 M$_\odot$ yr$^{-1}$ located at $1 \lesssim z \lesssim 2$. Two of the lensed radio sources are also detected in the Chandra image of the cluster. These two sources are likely AGN, given their $2-10$ keV X-ray luminosities of $\sim 10^{43-44}$ erg s$^{-1}$. From the derived radio luminosity function, we find evidence for an increase in the number density of radio sources at $0.6<z<2.0$, compared to a $z < 0.3$ sample. Our observations indicate that deep radio imaging of lensing clusters can be used to study star forming galaxies, with star formation rates as low as $\sim10$ M$_{\odot}$ yr$^{-1}$, at the peak of cosmic star formation history.

We present results based on X-ray, optical, and radio observations of the massive galaxy cluster CIZA J0107.7+5408. We find that this system is a post core passage, dissociative, binary merger, with the optical galaxy density peaks of each subcluster leading their associated X-ray emission peaks. This separation occurs because the diffuse gas experiences ram pressure forces while the effectively collisionless galaxies (and presumably their associated dark matter halos) do not. This system contains double peaked diffuse radio emission, possibly a double radio relic with the relics lying along the merger axis and also leading the X-ray cores. We find evidence for a temperature peak associated with the SW relic, likely created by the same merger shock that is powering the relic radio emission in this region. Thus, this system is a relatively rare clean example of a dissociative binary merger, which can in principle be used to place constraints on the self-interaction cross-section of dark matter. Low frequency radio observations reveal ultra-steep spectrum diffuse radio emission that is not correlated with the X-ray, optical, or high frequency radio emission. We suggest that these sources are radio phoenixes, which are preexisting non-thermal particle populations that have been re-energized through adiabatic compression by the same merger shocks that power the radio relics. Finally, we place upper limits on inverse Compton emission from the SW radio relic.

The HST Frontier Fields cluster MACS J1149.6+2223 is one of the most complex merging clusters, believed to consist of four dark matter halos. We present results from deep (365 ks) Chandra observations of the cluster, which reveal the most distant cold front (z=0.544) discovered to date. In the cluster outskirts, we also detect hints of a surface brightness edge that could be the bow shock preceding the cold front. The substructure analysis of the cluster identified several components with large relative radial velocities, thus indicating that at least some collisions occur almost along the line of sight. The inclination of the mergers with respect to the plane of the sky poses significant observational challenges at X-ray wavelengths. MACS J1149.6+2223 possibly hosts a steep-spectrum radio halo. If the steepness of the radio halo is confirmed, then the radio spectrum, combined with the relatively regular ICM morphology, could indicate that MACS J1149.6+2223 is an old merging cluster.

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 present deep LOFAR observations between 120-181 MHz of the "Toothbrush" (RX J0603.3+4214), a cluster that contains one of the brightest radio relic sources known. Our LOFAR observations exploit a new and novel calibration scheme to probe 10 times deeper than any previous study in this relatively unexplored part of the spectrum. The LOFAR observations, when combined with VLA, GMRT, and Chandra X-ray data, provide new information about the nature of cluster merger shocks and their role in re-accelerating relativistic particles. We derive a spectral index of $\alpha = -0.8 \pm 0.1$ at the northern edge of the main radio relic, steepening towards the south to $\alpha \approx - 2$. The spectral index of the radio halo is remarkably uniform ($\alpha = -1.16$, with an intrinsic scatter of $\leq 0.04$). The observed radio relic spectral index gives a Mach number of $\mathcal{M} = 2.8^{+0.5}_{-0.3}$, assuming diffusive shock acceleration (DSA). However, the gas density jump at the northern edge of the large radio relic implies a much weaker shock ($\mathcal{M} \approx 1.2$, with an upper limit of $\mathcal{M} \approx 1.5$). The discrepancy between the Mach numbers calculated from the radio and X-rays can be explained if either (i) the relic traces a complex shock surface along the line of sight, or (ii) if the radio relic emission is produced by a re-accelerated population of fossil particles from a radio galaxy. Our results highlight the need for additional theoretical work and numerical simulations of particle acceleration and re-acceleration at cluster merger shocks.

LOFAR, the Low-Frequency Array, is a powerful new radio telescope operating between 10 and 240 MHz. LOFAR allows detailed sensitive high-resolution studies of the low-frequency radio sky. At the same time LOFAR also provides excellent short baseline coverage to map diffuse extended emission. However, producing high-quality deep images is challenging due to the presence of direction dependent calibration errors, caused by imperfect knowledge of the station beam shapes and the ionosphere. Furthermore, the large data volume and presence of station clock errors present additional difficulties. In this paper we present a new calibration scheme, which we name facet calibration, to obtain deep high-resolution LOFAR High Band Antenna images using the Dutch part of the array. This scheme solves and corrects the direction dependent errors in a number of facets that cover the observed field of view. Facet calibration provides close to thermal noise limited images for a typical 8 hr observing run at $\sim$ 5arcsec resolution, meeting the specifications of the LOFAR Tier-1 northern survey.

Using overlapping fields with space-based Hubble Space Telescope (HST) and ground-based Subaru Telescope imaging we identify a population of blended galaxies that are blended to such a large degree that they are detected as single objects in the ground-based monochromatic imaging, which we label as 'ambiguous blends'. For deep imaging data, such as the depth targeted with the Large Synoptic Survey Telescope (LSST), the ambiguous blend population is both large ($\sim 14$%) and has a distribution of ellipticities that is different from that of unblended objects in a way that will likely be important for the weak lensing measurements. Most notably, for a limiting magnitude of $i \sim 27$ we find that ambiguous blending results in a ~14% increase in shear noise (or ~12% decrease in the effective projected number density of lensed galaxies; neff) due to 1) larger intrinsic ellipticity dispersion, 2) a scaling with the galaxy number density $N_{gal}$ that is shallower than 1/$\sqrt{N_{gal}}$. For the LSST Gold Sample ($i < 25.3$) there is a ~7% increase in shear noise (or ~7% decrease in $n_{eff}$). More importantly than these increases in the shear noise, we find that the ellipticity distribution of ambiguous blends has an RMS 13% larger than that of non-blended galaxies. Given the need of future weak lensing surveys to constrain the ellipticity distribution of galaxies to better than a percent in order to mitigate cosmic shear multiplicative biases, the different ellipticity distribution of ambiguous blends could be a dominant systematic if unaccounted for.

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.

CIZA J2242.8+5301 ($z = 0.188$, nicknamed 'Sausage') is an extremely massive ($M_{200}\sim 2.0 \times 10^{15}M_\odot$ ), merging cluster with shock waves towards its outskirts, which was found to host numerous emission-line galaxies. We performed extremely deep Westerbork Synthesis Radio Telescope HI observations of the 'Sausage' cluster to investigate the effect of the merger and the shocks on the gas reservoirs fuelling present and future star formation (SF) in cluster members. By using spectral stacking, we find that the emission-line galaxies in the 'Sausage' cluster have, on average, as much HI gas as field galaxies (when accounting for the fact cluster galaxies are more massive than the field galaxies), contrary to previous studies. Since the cluster galaxies are more massive than the field spirals, they may have been able to retain their gas during the cluster merger. The large HI reservoirs are expected to be consumed within $\sim0.75-1.0$ Gyr by the vigorous SF and AGN activity and/or driven out by the out-flows we observe. We find that the star-formation rate in a large fraction of H$\alpha$ emission-line cluster galaxies correlates well with the radio broad band emission, tracing supernova remnant emission. This suggests that the cluster galaxies, all located in post-shock regions, may have been undergoing sustained SFR for at least 100 Myr. This fully supports the interpretation proposed by Stroe et al. (2015) and Sobral et al. (2015) that gas-rich cluster galaxies have been triggered to form stars by the passage of the shock.

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the S\'{e}rsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

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.

We present results from Chandra, XMM-Newton, and ROSAT observations of the Planck SZ-detected cluster A3716 (PLCKG345.40-39.34 - G345). We show that G345 is, in fact, two subclusters separated on the sky by 400 kpc. We measure the subclusters' gas temperatures (~ 2-3 keV), total (~ 1-2 x 10^14 solar masses) and gas (~ 1-2 x 10^13 solar masses) masses, gas mass fraction within r500, entropy profiles, and X-ray luminosities (~ 10^43 erg/s). Using the gas density and temperature profiles for both subclusters, we show that there is good (0.8 sigma) agreement between the expected Sunyaev-Zel'dovich signal predicted from the X-ray data and that measured from the Planck mission, and better agreement within 0.6 sigma when we re-computed the Planck value assuming a two component cluster model, with relative amplitudes fixed based on the X-ray data. Dynamical analysis shows that the two galaxy subclusters are very likely (> 97% probability) gravitationally bound, and in the most likely scenario, the subclusters will undergo core passage in 500 +- 200 Myr. The northern subcluster is centrally peaked and has a low entropy core, while the southern subcluster has a high central entropy. The high central entropy in the southern subcluster can be explained either by the mergers of several groups, as suggested by the presence of five giant ellipticals or by AGN energy injection, as suggested by the presence of a strong radio source in one of its massive elliptical galaxies, or by a combination of both processes.

Extended steep-spectrum radio emission in a galaxy cluster is usually associated with a recent merger. However, given the complex scenario of galaxy cluster mergers, many of the discovered sources hardly fit into the strict boundaries of a precise taxonomy. This is especially true for radio phoenixes that do not have very well defined observational criteria. Radio phoenixes are aged radio galaxy lobes whose emission is reactivated by compression or other mechanisms. Here, we present the detection of a radio phoenix close to the moment of its formation. The source is located in Abell 1033, a peculiar galaxy cluster which underwent a recent merger. To support our claim, we present unpublished Westerbork Synthesis Radio Telescope and Chandra observations together with archival data from the Very Large Array and the Sloan Digital Sky Survey. We discover the presence of two sub-clusters displaced along the N-S direction. The two sub-clusters probably underwent a recent merger which is the cause of a moderately perturbed X-ray brightness distribution. A steep-spectrum extended radio source very close to an AGN is proposed to be a newly born radio phoenix: the AGN lobes have been displaced/compressed by shocks formed during the merger event. This scenario explains the source location, morphology, spectral index, and brightness. Finally, we show evidence of a density discontinuity close to the radio phoenix and discuss the consequences of its presence.

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.

Point estimators for the shearing of galaxy images induced by gravitational lensing involve a complex inverse problem in the presence of noise, pixelization, and model uncertainties. We present a probabilistic forward modeling approach to gravitational lensing inference that has the potential to mitigate the biased inferences in most common point estimators and is practical for upcoming lensing surveys. The first part of our statistical framework requires specification of a likelihood function for the pixel data in an imaging survey given parameterized models for the galaxies in the images. We derive the lensing shear posterior by marginalizing over all intrinsic galaxy properties that contribute to the pixel data (i.e., not limited to galaxy ellipticities) and learn the distributions for the intrinsic galaxy properties via hierarchical inference with a suitably flexible conditional probabilitiy distribution specification. We use importance sampling to separate the modeling of small imaging areas from the global shear inference, thereby rendering our algorithm computationally tractable for large surveys. With simple numerical examples we demonstrate the improvements in accuracy from our importance sampling approach, as well as the significance of the conditional distribution specification for the intrinsic galaxy properties when the data are generated from an unknown number of distinct galaxy populations with different morphological characteristics.

The statistics of shear peaks have been shown to provide valuable cosmological information beyond the power spectrum, and will be an important constraint of models of cosmology with the large survey areas provided by forthcoming astronomical surveys. Surveys include masked areas due to bright stars, bad pixels etc, which must be accounted for in producing constraints on cosmology from shear maps. We advocate a forward-modeling approach, where the impact of masking (and other survey artifacts) are accounted for in the theoretical prediction of cosmological parameters, rather than removed from survey data. We use masks based on the Deep Lens Survey, and explore the impact of up to 37% of the survey area being masked on LSST and DES-scale surveys. By reconstructing maps of aperture mass, the masking effect is smoothed out, resulting in up to 14% smaller statistical uncertainties compared to simply reducing the survey area by the masked area. We show that, even in the presence of large survey masks, the bias in cosmological parameter estimation produced in the forward-modeling process is ~1%, dominated by bias caused by limited simulation volume. We also explore how this potential bias scales with survey area and find that small survey areas are more significantly impacted by the differences in cosmological structure in the data and simulated volumes, due to cosmic variance.

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.

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.

Context. While searching the NRAO VLA Sky Survey (NVSS) for diffuse radio emission, we have serendipitously discovered extended radio emission close to the Galactic plane. The radio morphology suggests the presence of a previously unknown Galactic supernova remnant. An unclassified {\gamma}-ray source detected by EGRET (3EG J1744-3934) is present in the same location and may stem from the interaction between high-speed particles escaping the remnant and the surrounding interstellar medium. Aims. Our aim is to confirm the presence of a previously unknown supernova remnant and to determine a possible association with the {\gamma}-ray emission 3EG J1744-3934. Methods. We have conducted optical and radio follow-ups of the target using the Dark Energy Camera (DECam) on the Blanco telescope at Cerro Tololo Inter-American Observatory (CTIO) and the Giant Meterwave Radio Telescope (GMRT). We then combined these data with archival radio and {\gamma}-ray observations. Results. While we detected the extended emission in four different radio bands (325, 1400, 2417, and 4850 MHz), no optical counterpart has been identified. Given its morphology and brightness, it is likely that the radio emission is caused by an old supernova remnant no longer visible in the optical band. Although an unclassified EGRET source is co-located with the supernova remnant, Fermi-LAT data do not show a significant {\gamma}-ray excess that is correlated with the radio emission. However, in the radial distribution of the {\gamma}-ray events, a spatially extended feature is related with SNR at a confidence level $\sim 1.5$ {\sigma}. Conclusions. We classify the newly discovered extended emission in the radio band as the old remnant of a previously unknown Galactic supernova: SNR G351.0-5.4.

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.

Merging galaxy clusters have become one of the most important probes of dark matter, providing evidence for dark matter over modified gravity and even constraints on the dark matter self-interaction cross-section. To properly constrain the dark matter cross-section it is necessary to understand the dynamics of the merger, as the inferred cross-section is a function of both the velocity of the collision and the observed time since collision. While the best understanding of merging system dynamics comes from N-body simulations, these are computationally intensive and often explore only a limited volume of the merger phase space allowed by observed parameter uncertainty. Simple analytic models exist but the assumptions of these methods invalidate their results near the collision time, plus error propagation of the highly correlated merger parameters is unfeasible. To address these weaknesses I develop a Monte Carlo method to discern the properties of dissociative mergers and propagate the uncertainty of the measured cluster parameters in an accurate and Bayesian manner. I introduce this method, verify it against an existing hydrodynamic N-body simulation, and apply it to two known dissociative mergers: 1ES 0657-558 (Bullet Cluster) and DLSCL J0916.2+2951 (Musket Ball Cluster). I find that this method surpasses existing analytic models - providing accurate (10% level) dynamic parameter and uncertainty estimates throughout the merger history. This coupled with minimal required a priori information (subcluster mass, redshift, and projected separation) and relatively fast computation (~6 CPU hours) makes this method ideal for large samples of dissociative merging clusters.

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.