We study the stellar population properties of the IRAC-detected $6 \lesssim z \lesssim 10$ galaxy candidates from the Spitzer UltRa Faint SUrvey Program (SURFS UP). Using the Lyman Break selection technique, we find a total of 16 new galaxy candidates at $6 \lesssim z \lesssim 10$ with $S/N \geq 3$ in at least one of the IRAC $3.6\mu$m and $4.5\mu$m bands. According to the best mass models available for the surveyed galaxy clusters, these IRAC-detected galaxy candidates are magnified by factors of $\sim 1.2$--$5.5$. We find that the IRAC-detected $6 \lesssim z \lesssim 10$ sample is likely not a homogeneous galaxy population: some are relatively massive (stellar mass as high as $4 \times 10^9\,M_{\odot}$) and evolved (age $\lesssim 500$ Myr) galaxies, while others are less massive ($M_{\text{stellar}}\sim 10^8\,M_{\odot}$) and very young ($\sim 10$ Myr) galaxies with strong nebular emission lines that boost their rest-frame optical fluxes. We identify two Ly$\alpha$ emitters in our sample from the Keck DEIMOS spectra, one at $z_{\text{Ly}\alpha}=6.76$ (in RXJ1347) and one at $z_{\text{Ly}\alpha}=6.32$ (in MACS0454). We show that IRAC $[3.6]-[4.5]$ color, when combined with photometric redshift, can be used to identify galaxies likely with strong nebular emission lines within certain redshift windows.

We present a reconstructions of galaxy-cluster-scale mass distributions from simulated gravitational lensing data sets including strong lensing, weak lensing shear, and measurements of quadratic image distortions -- flexion. The lensing data is constructed to make a direct comparison between mass reconstructions with and without flexion. We show that in the absence of flexion measurements, significant galaxy-group scale substructure can remain undetected in the reconstructed mass profiles, and that the resulting profiles underestimate the aperture mass in the substructure regions by $\sim25-40\%$. When flexion is included, subhaloes down to a mass of $\sim3\times10^{12}$ M$_\odot$ can be detected at an angular resolution smaller than 10\arcsec. Aperture masses from profiles reconstructed with flexion match the input distribution values to within an error of $\sim13\%$, including both statistical error and scatter. This demonstrates the important constraint that flexion measurements place on substructure in galaxy clusters and its utility for producing high-fidelity mass reconstructions.

We present a detailed gravitational lens model of the galaxy cluster RCS2 J232727.6-020437. Due to cosmological dimming of cluster members and ICL, its high redshift ($z=0.6986$) makes it ideal for studying background galaxies. Using new ACS and WFC3/IR HST data, we identify 16 multiple images. From MOSFIRE follow up, we identify a strong emission line in the spectrum of one multiple image, likely confirming the redshift of that system to $z=2.083$. With a highly magnified ($\mu\gtrsim2$) source plane area of $\sim0.7$ arcmin$^2$ at $z=7$, RCS2 J232727.6-020437 has a lensing efficiency comparable to the Hubble Frontier Fields clusters. We discover four highly magnified $z\sim7$ candidate Lyman-break galaxies behind the cluster, one of which may be multiply-imaged. Correcting for magnification, we find that all four galaxies are fainter than $0.5 L_{\star}$. One candidate is detected at ${>10\sigma}$ in both Spitzer/IRAC [3.6] and [4.5] channels. A spectroscopic follow-up with MOSFIRE does not result in the detection of the Lyman-alpha emission line from any of the four candidates. From the MOSFIRE spectra we place median upper limits on the Lyman-alpha flux of $5-14 \times 10^{-19}\, \mathrm{erg \,\, s^{-1} cm^{-2}}$ ($5\sigma$).

We present Spitzer/IRAC observations of nine $z'$-band dropouts highly magnified (2<mu<12) by the Bullet Cluster. We combine archival imaging with our Exploratory program (SURFS UP), which results in a total integration time of ~30 hr per IRAC band. We detect (>3sigma) in both IRAC bands the brightest of these high-redshift galaxies, with [3.6]=23.80+-0.28 mag, [4.5]=23.78+-0.25 mag, and (H-[3.6])=1.17+-0.32 mag. The remaining eight galaxies are undetected to [3.6]~26.4 mag and [4.5]~26.0 mag with stellar masses of ~5x10^7 M_sol. The detected galaxy has an estimated magnification of mu=12+-4, which implies this galaxy has an ultraviolet luminosity of L_1500~0.3 L*_{z=7} --- the lowest luminosity individual source detected in IRAC at z>7. By modeling the broadband photometry, we estimate the galaxy has an intrinsic star-formation rate of SFR~1.3 M_sol/yr and stellar mass of M~2x10^9 M_sol, which gives a specific star-formation rate of sSFR~0.7 Gyr^-1. If this galaxy had sustained this star-formation rate since z~20, it could have formed the observed stellar mass (to within a factor of ~2), we also discuss alternate star-formation histories and argue the exponentially-increasing model is unlikely. Finally, based on the intrinsic star-formation rate, we estimate this galaxy has a likely [C II] flux of <f_[C II]> = 10^{-17} erg/s/cm2.

SURFSUP is a joint Spitzer and HST Exploration Science program using 10 galaxy clusters as cosmic telescopes to study z >~ 7 galaxies at intrinsically lower luminosities, enabled by gravitational lensing, than blank field surveys of the same exposure time. Our main goal is to measure stellar masses and ages of these galaxies, which are the most likely sources of the ionizing photons that drive reionization. Accurate knowledge of the star formation density and star formation history at this epoch is necessary to determine whether these galaxies indeed reionized the universe. Determination of the stellar masses and ages requires measuring rest frame optical light, which only Spitzer can probe for sources at z >~ 7, for a large enough sample of typical galaxies. Our program consists of 550 hours of Spitzer/IRAC imaging covering 10 galaxy clusters with very well-known mass distributions, making them extremely precise cosmic telescopes. We combine our data with archival observations to obtain mosaics with ~30 hours exposure time in both 3.6$\mu$m and 4.5$\mu$m in the central 4 arcmin x 4 arcmin field and ~15 hours in the flanking fields. This results in 3-$\sigma$ sensitivity limits of ~26.6 and ~26.2AB magnitudes for the central field in the IRAC 3.6 and 4.5$\mu$m bands, respectively. To illustrate the survey strategy and characteristics we introduce the sample, present the details of the data reduction and demonstrate that these data are sufficient for in-depth studies of z >~ 7 sources (using a z=9.5 galaxy behind MACSJ1149.5+2223 as an example). For the first cluster of the survey (the Bullet Cluster) we have released all high-level data mosaics and IRAC empirical PSF models. In the future we plan to release these data products for the entire survey.

Measuring dark matter substructure within galaxy cluster haloes is a fundamental probe of the Lambda-CDM model of structure formation. Gravitational lensing is a technique for measuring the total mass distribution which is independent of the nature of the gravitating matter, making it a vital tool for studying these dark-matter dominated objects. We present a new method for measuring weak gravitational lensing flexions, the gradients of the lensing shear field, to measure mass distributions on small angular scales. While previously published methods for measuring flexions focus on measuring derived properties of the lensed images, such as shapelet coefficients or surface brightness moments, our method instead fits a mass-sheet-transformation-invariant Analytic Image Model (AIM) to the each galaxy image. This simple parametric model traces the distortion of lensed image isophotes and constrains the flexion fields. We test the AIM method using simulated data images with realistic noise and a variety of unlensed image properties, and show that it successfully reproduces the input flexion fields. We also apply the AIM method for flexion measurement to Hubble Space Telescope observations of Abell 1689, and detect mass structure in the cluster using flexions measured with the AIM method.

We present new Chandra observations of a high redshift (z~1) galaxy cluster discovered in the Red-Sequence Cluster Survey (RCS): RCS043938-2904.7. X-ray luminosity measurements and mass estimates are consistent with L_X-T_X and M_delta-T_X relationships obtained from low-redshift data. Assuming a single cluster, X-ray mass estimates are a factor of ~10-100 below the red-sequence optical richness mass estimate. Optical spectroscopy reveals that this cluster comprises two components which are close enough to perhaps be physically associated. We present simple modeling of this two-component system which then yields an X-ray mass and optical richness consistent with expectations from statistical samples of lower redshift clusters. An unexpectedly high gas mass fraction is measured assuming a single cluster, which independently supports this interpretation. Additional observations will be necessary to confirm the excess gas mass fraction and to constrain the mass distribution.

We present the results of Chandra observations of 13 optically-selected clusters with 0.6<z< 1.1, discovered via the Red-sequence Cluster Survey (RCS). All but one are detected at S/N>3; though 3 were not observed long enough to support detailed analysis. Surface brightness profiles are fit to beta-models. Integrated spectra are extracted within R(2500), and Tx and Lx information is obtained. We derive gas and total masses within R(2500) and R(500). Cosmologically corrected scaling relations are investigated, and we find the RCS clusters to be consistent with self-similar scaling expectations. However discrepancies exist between the RCS sample and lower-z X-ray selected samples for relationships involving Lx, with the higher-z RCS clusters having lower Lx for a given Tx. In addition, we find that gas mass fractions within R(2500) for the high-z RCS sample are lower than expected by a factor of ~2. This suggests that the central entropy of these high-z objects has been elevated by processes such as pre-heating, mergers, and/or AGN outbursts, that their gas is still infalling, or that they contain comparatively more baryonic matter in the form of stars. Finally, relationships between red-sequence optical richness (Bgc) and X-ray properties are fit to the data. For systems with measured Tx, we find that optical richness correlates with both Tx and mass, having a scatter of ~30% with mass for both X-ray and optically-selected clusters. However we also find that X-ray luminosity is not well correlated with richness, and that several of our sample appear to be significantly X-ray faint.