Search

We present updated cosmological constraints from measurements of the gas mass fractions ($f_{gas}$) of massive, dynamically relaxed galaxy clusters. Our new data set has greater leverage on models of dark energy, thanks to the addition of the Perseus Cluster at low redshifts, two new clusters at redshifts $z>0.97$, and significantly longer observations of four clusters at $0.6, Comment: To be published in MNRAS. See https://github.com/abmantz/fgas-cosmo for code and https://github.com/abmantz/fgas-2021-paper for data presented in figures/tables

We develop a method to calibrate u-band photometry based on the observed color of blue galactic halo stars. The galactic halo stars belong to an old stellar population of the Milky Way and have relatively low metallicity. The "blue tip" of the halo population -- the main sequence turn-off (MSTO) stars -- is known to have a relatively uniform intrinsic edge u-g color with only slow spatial variation. In SDSS data, the observed variation is correlated with galactic latitude, which we attribute to contamination by higher-metallicity disk stars and fit with an empirical curve. This curve can then be used to calibrate u-band imaging if g-band imaging of matching depth is available. Our approach can be applied to single-field observations at $|b| > 30^\circ$, and removes the need for standard star observations or overlap with calibrated u-band imaging. We include in our method the calibration of g-band data with ATLAS-Refcat2. We test our approach on stars in KiDS DR 4, ATLAS DR 4, and DECam imaging from the NOIRLab Source Catalog (NSC DR2), and compare our calibration with SDSS. For this process, we use synthetic magnitudes to derive the color equations between these datasets, in order to improve zero-point accuracy. We find an improvement for all datasets, reaching a zero-point precision of 0.016 mag for KiDS (compared to the original 0.033 mag), 0.020 mag for ATLAS (originally 0.027 mag), and 0.016 mag for DECam (originally 0.041 mag). Thus, this method alone reaches the goal of 0.02 mag photometric precision in u-band for the Rubin Observatory's Legacy Survey of Space and Time (LSST)., Published in MNRAS. Code available at https://github.com/shuang92/blue_tip_calibration

We explore synergies between the space-based Wide-Field Infrared Survey Telescope (WFIRST) and the ground-based Rubin Observatory Legacy Survey of Space and Time (LSST). In particular, we consider a scenario where the currently envisioned survey strategy for WFIRST's High Latitude Survey (HLS), i.e., 2000 square degrees in four narrow photometric bands is altered in favor of a strategy that combines rapid coverage of the LSST area (to full LSST depth) in one band. We find that a 5-month WFIRST survey in the W-band can cover the full LSST survey area providing high-resolution imaging for >95% of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where WFIRST spends 1.5 years covering the LSST area. For this second scenario we quantify the constraining power on dark energy equation of state parameters from a joint weak lensing and galaxy clustering analysis, and compare it to an LSST-only survey and to the Reference WFIRST HLS survey. Our survey simulations are based on the WFIRST exposure time calculator and redshift distributions from the CANDELS catalog. Our statistical uncertainties account for higher-order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modeling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. Assuming the 5-month WFIRST wide scenario, we find a significant increase in constraining power for the joint LSST+WFIRST wide survey compared to LSST Y10 (FoM(Wwide)= 2.4 FoM(LSST)) and compared to LSST+WFIRST HLS (FoM(Wwide)= 5.5 FoM(HLS)).

Galaxy clusters have a triaxial matter distribution. The weak-lensing signal, an important part in cosmological studies, measures the projected mass of all matter along the line-of-sight, and therefore changes with the orientation of the cluster. Studies suggest that the shape of the brightest cluster galaxy (BCG) in the centre of the cluster traces the underlying halo shape, enabling a method to account for projection effects. We use 324 simulated clusters at four redshifts between 0.1 and 0.6 from `The Three Hundred Project' to quantify correlations between the orientation and shape of the BCG and the halo. We find that haloes and their embedded BCGs are aligned, with an average $\sim$20 degree angle between their major axes. The bias in weak lensing cluster mass estimates correlates with the orientation of both the halo and the BCG. Mimicking observations, we compute the projected shape of the BCG, as a measure of the BCG orientation, and find that it is most strongly correlated to the weak-lensing mass for relaxed clusters. We also test a 2-dimensional cluster relaxation proxy measured from BCG mass isocontours. The concentration of stellar mass in the projected BCG core compared to the total stellar mass provides an alternative proxy for the BCG orientation. We find that the concentration does not correlate to the weak-lensing mass bias, but does correlate with the true halo mass. These results indicate that the BCG shape and orientation for large samples of relaxed clusters can provide information to improve weak-lensing mass estimates., 15 pages, 9 figures, Figure 7 is key plot, Updated to match version accepted by MNRAS

Steinhardt et al., arXiv:2001.09999v2, The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope (HST) Treasury program taking data from 2018 to 2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in Wide Field Camera 3/IR F105W, F125W, and F160W and Advanced Camera for Surveys/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multiwavelength data sets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies that would otherwise lie below HST detection limits and model foreground clusters to study the properties of dark matter and galaxy assembly. The expanded area will provide the first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient follow-up. BUFFALO data including mosaics, value-added catalogs, and cluster mass distribution models will be released via MAST on a regular basis as the observations and analysis are completed for the six individual clusters., C.S. acknowledges support from the ERC Consolidator Grant funding scheme (project ConTExT, grant No. 648179). Y.M.B. acknowledges funding from the EU Horizon 2020 research and innovation program under Marie Skłodowska-Curie grant agreement 747645 (ClusterGal). J.M.D. acknowledges the support of project PGC2018-101814-B-100 (MCIU/AEI/MINECO/FEDER, UE), Ministerio de Ciencia, Investigación y Universidades.

Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05, Comment: 20 pages, 9 figures, 4 tables, submitted to MNRAS

Weak-lensing measurements of the masses of galaxy clusters are commonly based on the assumption of spherically symmetric density profiles. Yet, the cold dark matter model predicts the shapes of dark matter halos to be triaxial. Halo triaxiality, and the orientation of the major axis with respect to the line of sight, are expected to be the leading cause of intrinsic scatter in weak-lensing mass measurements. The shape of central cluster galaxies (Brightest Cluster Galaxies; BCGs) is expected to follow the shape of the dark matter halo. Here we investigate the use of BCG ellipticity as predictor of the weak-lensing mass bias in individual clusters compared to the mean. Using weak lensing masses $M^{\rm WL}_{500}$ from the Weighing the Giants project, and $M_{500}$ derived from gas masses as low-scatter mass proxy, we find that, on average, the lensing masses of clusters with the roundest / most elliptical 25% of BCGs are biased $\sim 20$% high / low compared to the average, as qualitatively predicted by the cold dark matter model. For cluster cosmology projects utilizing weak-lensing mass estimates, the shape of the BCG can thus contribute useful information on the effect of orientation bias in weak lensing mass estimates as well as on cluster selection bias., Comment: 9 pages, 3 figures, 3 tables

This paper presents multiwavelength photometric catalogues of the last two Hubble Frontier Fields (HFF), the massive galaxy clusters Abell 370 and RXC J2248.7-4431. The photometry ranges from imaging performed on the Hubble Space Telescope (HST) to ground based Very Large Telescope (VLT) and Spitzer/IRAC, in collaboration with the ASTRODEEP team, and using the ASTRODEEP pipeline. While the main purpose of this paper is to release the catalogues, we also perform, as a proof of concept, a brief analysis of z > 6 objects selected using drop-out method, as well as spectroscopically confirmed sources and multiple images in both clusters. While dropout methods yield a sample of high-z galaxies, the addition of longer wavelength data reveals that as expected the samples have substantial contamination at the ~30-45% level by dusty galaxies at lower redshifts. Furthermore, we show that spectroscopic redshifts are still required to unambiguously determine redshifts of multiply imaged systems. Finally, the now publicly available ASTRODEEP catalogues were combined for all HFFs and used to explore stellar properties of a large sample of 20,000 galaxies across a large photometric redshift range. The powerful magnification provided by the HFF clusters allows us an exploration of the properties of galaxies with intrinsic stellar masses as low as $M_* \gtrsim 10^7M_{\odot}$ and intrinsic star formation rates $\mbox{SFRs}\sim 0.1\mbox{-}1M_\odot/\mbox yr$ at z > 6., Replaced to match the accepted version

We describe the first spectroscopic tomographic (spectrotomographic) weak-lensing measurement for a galaxy cluster based only on background galaxies with spectroscopically determined redshifts. We use the massive cluster A2029 to demonstrate the power of combining spectroscopy and lensing to obtain accurate masses and to overcome biases from contamination and photometric redshift errors. We detect the shear signal from the cluster at > 3.9σ significance. The shear signal scales with source redshift in a way that is consistent with the angular diameter distance ratio variation in a ΛCDM universe. Furthermore, the amplitude of the measured signal is consistent with the X-ray mass. Upcoming spectroscopic instruments such as the Prime Focus Spectrograph on Subaru will permit spectrotomographic weak-lensing measurements with a signal-to-noise ratio comparable with current photometric-redshift-based weak-lensing measurements for hundreds of galaxy clusters. Thus, spectrotomography may enable sensitive cosmological constraints that complement and are independent of other measurement techniques.

We present here the weak gravitational lensing detection of four nearby galaxy clusters in the southern sky: Abell 2029, Abell 85, Abell 1606 and Abell 2457. The weak lensing detections of Abell 1606 and Abell 2457 are the first in the literature. This work capitalizes on the wide field of view of the Dark Energy Camera at the Cerro Tololo Inter-American Observatory, which we use to obtain deep, multi-wavelength imaging of all targets. We publish maps of the clusters' projected mass distributions, and obtain the $M_{200}$ of their clusters through NFW profile fits to the two-dimensional tangential ellipticity signal., 29 pages, 18 figures

Euclid and the Large Synoptic Survey Telescope (LSST) are poised to dramatically change the astronomy landscape early in the next decade. The combination of high cadence, deep, wide-field optical photometry from LSST with high resolution, wide-field optical photometry and near-infrared photometry and spectroscopy from Euclid will be powerful for addressing a wide range of astrophysical questions. We explore Euclid/LSST synergy, ignoring the political issues associated with data access to focus on the scientific, technical, and financial benefits of coordination. We focus primarily on dark energy cosmology, but also discuss galaxy evolution, transient objects, solar system science, and galaxy cluster studies. We concentrate on synergies that require coordination in cadence or survey overlap, or would benefit from pixel-level co-processing that is beyond the scope of what is currently planned, rather than scientific programs that could be accomplished only at the catalog level without coordination in data processing or survey strategies. We provide two quantitative examples of scientific synergies: the decrease in photo-z errors (benefitting many science cases) when high resolution Euclid data are used for LSST photo-z determination, and the resulting increase in weak lensing signal-to-noise ratio from smaller photo-z errors. We briefly discuss other areas of coordination, including high performance computing resources and calibration data. Finally, we address concerns about the loss of independence and potential cross-checks between the two missions and potential consequences of not collaborating., 33 pages, 7 figures, to appear in ApJS, revised with additional references and minor corrections

Exploiting the data of the Grism Lens-Amplified Survey from Space (GLASS), we characterize the spatial distribution of star formation in 76 high star forming galaxies in 10 clusters at 0.3< z, Comment: ApJ in press (16 pages, 8 figures)

We constrain cold dark energy of negligible sound speed using galaxy cluster abundance observations. In contrast to standard quasi-homogeneous dark energy, negligible sound speed implies clustering of the dark energy fluid at all scales, allowing us to measure the effects of dark energy perturbations at cluster scales. We compare those models and set the stage for using non-linear information from semi-analytical modelling in cluster growth data analyses. For this, we recalibrate the halo mass function with non-linear characteristic quantities, the spherical collapse threshold and virial overdensity, that account for model and redshift dependent behaviours, as well as an additional mass contribution for cold dark energy. We present the first constraints from this cold dark matter plus cold dark energy mass function using our cluster abundance likelihood, which self-consistently accounts for selection effects, covariances and systematic uncertainties. We combine cluster growth data with CMB, SNe Ia and BAO data, and find a shift between cold versus quasi-homogeneous dark energy of up to $1\sigma$. We make a Fisher matrix forecast of constraints attainable with cluster growth data from the on-going Dark Energy Survey (DES). For DES, we predict $\sim$50$\%$ tighter constraints on $\left(\Omega_\mathrm{m},w \right)$ for cold dark energy versus $w$CDM models, with the same free parameters. Overall, we show that cluster abundance analyses are sensitive to cold dark energy, an alternative, viable model that should be routinely investigated alongside the standard dark energy scenario., Comment: 14 pages, 6 figures, 3 tables. Accepted for publiction in MNRAS

This is the first in a series of papers in which we measure accurate weak-lensing masses for 51 of the most X-ray luminous galaxy clusters known at redshifts 0.15, Comment: 26 pages, 19 figures (Appendix C not included). Accepted after minor revision

We report weak-lensing masses for 51 of the most X-ray luminous galaxy clusters known. This cluster sample, introduced earlier in this series of papers, spans redshifts 0.15 < z_cl < 0.7, and is well suited to calibrate mass proxies for current cluster cosmology experiments. Cluster masses are measured with a standard `color-cut' lensing method from three-filter photometry of each field. Additionally, for 27 cluster fields with at least five-filter photometry, we measure high-accuracy masses using a new method that exploits all information available in the photometric redshift posterior probability distributions of individual galaxies. Using simulations based on the COSMOS-30 catalog, we demonstrate control of systematic biases in the mean mass of the sample with this method, from photometric redshift biases and associated uncertainties, to better than 3%. In contrast, we show that the use of single-point estimators in place of the full photometric redshift posterior distributions can lead to significant redshift-dependent biases on cluster masses. The performance of our new photometric redshift-based method allows us to calibrate `color-cut` masses for all 51 clusters in the present sample to a total systematic uncertainty of ~7% on the mean mass, a level sufficient to significantly improve current cosmology constraints from galaxy clusters. Our results bode well for future cosmological studies of clusters, potentially reducing the need for exhaustive spectroscopic calibration surveys as compared to other techniques, when deep, multi-filter optical and near-IR imaging surveys are coupled with robust photometric redshift methods., 25 pages, 14 figures. Accepted after minor revision. Lensing masses (table 4) and other data will be made available at http://www.slac.stanford.edu/~dapple/work/wtg.html

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., Comment: ApJ in press

We demonstrate that deep good-seeing VLT/HAWK-I $K_\mathrm{s}$ images complemented with $g$+$z$-band photometry can yield a sensitivity for weak lensing studies of massive galaxy clusters at redshifts \mbox{$0.7\lesssim z \lesssim 1.1$}, which is almost identical to the sensitivity of HST/ACS mosaics of single-orbit depth. Key reasons for this good performance are the excellent image quality frequently achievable for $K_\mathrm{s}$ imaging from the ground, a highly effective photometric selection of background galaxies, and a galaxy ellipticity dispersion that is noticeably lower than for optically observed high-redshift galaxy samples. Incorporating results from the 3D-HST and UltraVISTA surveys we also obtained a more accurate calibration of the source redshift distribution than previously achieved for similar optical weak lensing data sets. Here we studied the extremely massive galaxy cluster RCS2$J$232727.7$-$020437 (\mbox{$z=0.699$}), combining deep VLT/\mbox{HAWK-I} $K_\mathrm{s}$ images (point spread function with a 0\farcs35 full width at half maximum) with LBT/LBC photometry. The resulting weak lensing mass reconstruction suggests that the cluster consists of a single overdensity, which is detected with a peak significance of $10.1��$. We constrained the cluster mass to \mbox{$M_\mathrm{200c}/(10^{15} \mathrm{M}_\odot) =2.06^{+0.28}_{-0.26}(\mathrm{stat.})\pm 0.12 (\mathrm{sys.})$} assuming a spherical Navarro, Frenk \& White model and simulation-based priors on the concentration, making it one of the most massive galaxy clusters known in the \mbox{$z\gtrsim 0.7$} Universe. We also cross-checked the HAWK-I measurements through an analysis of overlapping HST/ACS images, yielding fully consistent estimates of the lensing signal., 17 pages, 14 figures, 1 table. Matches the version published in A&A. Includes minor changes to the text and corrections from the language editor

Post-starburst (E+A or k+a) spectra, characterized by their exceptionally strong Balmer lines in absorption and the lack of emission lines, belong to galaxies in which the star formation activity ended abruptly sometime during the past Gyr. We perform a spectral analysis of galaxies in clusters, groups, poor groups and the field at z=0.4-0.8 based on the ESO Distant Cluster Survey. The incidence of k+a's at these redshifts depends strongly on environment. K+a's reside preferentially in clusters and, unexpectedly, in a subset of the sigma = 200-400 km/s groups, those that have a low fraction of [OII] emitters. In these environments, 20-30% of the recently star-forming galaxies have had their star formation activity recently truncated. In contrast, there are proportionally fewer k+a's in the field, the poor groups and groups with a high [OII] fraction. The incidence of k+a galaxies correlates with the cluster velocity dispersion: more massive clusters have higher proportions of k+a's. Spectra of dusty starburst candidates, with strong Balmer absorption and emission lines, present a very different environmental dependence from k+a's. They are numerous in all environments at z=0.4-0.8, but they are especially numerous in all types of groups, favoring the hypothesis of triggering by a merger. Our observations are consistent with previous suggestions that cluster k+a galaxies are observed in a transition phase as massive S0 and Sa galaxies, evolving from star-forming later types to passively evolving early-type galaxies. The correlation between k+a fraction and cluster sigma supports the hypothesis that k+a galaxies in clusters originate from processes related to the intracluster medium, while several possibilities are discussed for the origin of the k+a frequency in low-[OII] groups.(abr.), Comment: accepted for publication in ApJ

We use the ROSAT All-Sky Survey to study the X-ray properties of a sample of 625 groups and clusters of galaxies selected from the Sloan Digital Sky Survey. We stack clusters with similar velocity dispersions and investigate whether their average X-ray luminosities and surface brightness profiles vary with the radio activity level of their central galaxies. We find that at a given value of $\sigma$, clusters with a central radio AGN have more concentrated X-ray surface brightness profiles, larger central galaxy masses, and higher X-ray luminosities than clusters with radio-quiet central galaxies. The enhancement in X-ray luminosity is more than a factor of two, is detected with better than 6$\sigma$ significance, and cannot be explained by X-ray emission from the radio AGN itself. This difference is largely due to a subpopulation of radio-quiet, high velocity dispersion clusters with low mass central galaxies. These clusters are underluminous at X-ray wavelengths when compared to otherwise similar clusters where the central galaxy is radio-loud, more massive, or both.

Using galaxy clusters from the ESO Distant Cluster Survey, we study how the distribution of galaxies along the colour-magnitude relation has evolved since z~0.8. While red-sequence galaxies in all these clusters are well described by an old, passively evolving population, we confirm our previous finding of a significant evolution in their luminosity distribution as a function of redshift. When compared to galaxy clusters in the local Universe, the high redshift EDisCS clusters exhibit a significant "deficit" of faint red galaxies. Combining clusters in three different redshift bins, and defining as `faint' all galaxies in the range 0.4 > L/L* > 0.1, we find a clear decrease in the luminous-to-faint ratio of red galaxies from z~0.8 to z~0.4. The amount of such a decrease appears to be in qualitative agreement with predictions of a model where the blue bright galaxies that populate the colour-magnitude diagram of high redshift clusters, have their star formation suppressed by the hostile cluster environment. Although model results need to be interpreted with caution, our findings clearly indicate that the red-sequence population of high-redshift clusters does not contain all progenitors of nearby red-sequence cluster galaxies. A significant fraction of these must have moved onto the red-sequence below z~0.8., 15 pages, 10 figures, accepted for publication in MNRAS

The abundance of massive galaxy clusters is a powerful probe of departures from General Relativity (GR) on cosmic scales. Despite current stringent constraints placed by stellar and galactic tests, on larger scales alternative theories of gravity such as $f(R)$ can still work as effective theories. Here we present constraints on two popular models of $f(R)$, Hu-Sawicki and "designer", derived from a fully self-consistent analysis of current samples of X-ray selected clusters and accounting for all the covariances between cosmological and astrophysical parameters. Using cluster number counts in combination with recent data from the cosmic microwave background (CMB) and the CMB lensing potential generated by large scale structures, as well as with other cosmological constraints on the background expansion history and its mean matter density, we obtain the upper bounds $\log_{10}|f_{R0}| < 4.79$ and $\log_{10}B_0 < 3.75$ at the 95.4 per cent confidence level, for the Hu-Sawicki (with $n=1$) and designer models, respectively. The robustness of our results derives from high quality cluster growth data for the most massive clusters known out to redshifts $z \sim 0.5$, a tight control of systematic uncertainties including an accurate and precise mass calibration from weak gravitational lensing data, and the use of the full shape of the halo mass function over the mass range of our data., 11 pages, 1 Table, 3 figures, matches published version on PRD

In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging, we have found four images of a single supernova forming an Einstein cross configuration around a redshift z=0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The cluster's gravitational potential also creates multiple images of the z=1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses., Published in the 6 March 2015 issue of Science; 17 pages, 7 figures, and 3 tables including Supplementary Materials

We analyse the spatially-resolved stellar populations of 9 local ($z 10^{11.3}$M$_{\odot}$) from the ATLAS$^{3D}$ survey (median [Z/H] $= 0.04\pm0.07$, $\Delta$[Z/H] $= -0.19\pm0.1$). However, massive early-type galaxies from ATLAS$^{3D}$ have consistently old ages (median Age $=12.0\pm3.8$Gyr). We also analyse the close massive companion galaxies of two of the BCGs. These galaxies have similar stellar populations to their respective BCGs., Comment: Accepted for publication, MNRAS, March 3, 2015

In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging, we have found four images of a single supernova forming an Einstein cross configuration around a redshift z = 0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The cluster's gravitational potential also creates multiple images of the z = 1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses.

This is the first in a series of papers studying the astrophysics and cosmology of massive, dynamically relaxed galaxy clusters. Here we present a new, automated method for identifying relaxed clusters based on their morphologies in X-ray imaging data. While broadly similar to others in the literature, the morphological quantities that we measure are specifically designed to provide a fair basis for comparison across a range of data quality and cluster redshifts, to be robust against missing data due to point-source masks and gaps between detectors, and to avoid strong assumptions about the cosmological background and cluster masses. Based on three morphological indicators - Symmetry, Peakiness and Alignment - we develop the SPA criterion for relaxation. This analysis was applied to a large sample of cluster observations from the Chandra and ROSAT archives. Of the 361 clusters which received the SPA treatment, 57 (16 per cent) were subsequently found to be relaxed according to our criterion. We compare our measurements to similar estimators in the literature, as well as projected ellipticity and other image measures, and comment on trends in the relaxed cluster fraction with redshift, temperature, and survey selection method. Code implementing our morphological analysis will be made available on the web., Comment: MNRAS, in press. 43 pages in total, of which 17 are tables (please think twice before printing). 18 figures, 4 tables. Machine-readable tables will be available from the journal and at the url below; code will be posted at http://www.slac.stanford.edu/~amantz/work/morph14/

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$��$m and 4.5$��$m in the central 4 arcmin x 4 arcmin field and ~15 hours in the flanking fields. This results in 3-$��$ sensitivity limits of ~26.6 and ~26.2AB magnitudes for the central field in the IRAC 3.6 and 4.5$��$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., 14 pages, 7 figures, accepted for publication in ApJ

We present improved methods for using stars found in astronomical exposures to calibrate both star and galaxy colors as well as to adjust the instrument flat field. By developing a spectroscopic model for the SDSS stellar locus in color-color space, synthesizing an expected stellar locus, and simultaneously solving for all unknown zeropoints when fitting to the instrumental locus, we increase the calibration accuracy of stellar locus matching. We also use a new combined technique to estimate improved flat-field models for the Subaru SuprimeCam camera, forming `star flats' based on the magnitudes of stars observed in multiple positions or through comparison with available SDSS magnitudes. These techniques yield galaxy magnitudes with reliable color calibration (< 0.01 - 0.02 mag accuracy) that enable us to estimate photometric redshift probability distributions without spectroscopic training samples. We test the accuracy of our photometric redshifts using spectroscopic redshifts z_s for ~5000 galaxies in 27 cluster fields with at least five bands of photometry, as well as galaxies in the COSMOS field, finding sigma((z_p - z_s)/(1 + z_s)) ~ 0.03 for the most probable redshift z_p. We show that the full posterior probability distributions for the redshifts of galaxies with five-band photometry exhibit good agreement with redshifts estimated from thirty-band photometry in the COSMOS field. The growth of shear with increasing distance behind each galaxy cluster shows the expected redshift-distance relation for a flat Lambda-CDM cosmology. Photometric redshifts and calibrated colors are used in subsequent papers to measure the masses of 51 galaxy clusters from their weak gravitational shear. We make our Python code for stellar locus matching available at http://big-macs-calibrate.googlecode.com; the code requires only a catalog and filter functions., Accepted by MNRAS with only minor revisions. Code available: http://big-macs-calibrate.googlecode.com (v2 latex symbols removed from abstract)

Using the VIMOS Integral Field Unit (IFU) spectrograph on the Very Large Telescope (VLT), we have spatially mapped the kinematic properties of 10 nearby Brightest Cluster Galaxies (BCGs) and 4 BCG companion galaxies located within a redshift of $z=0.1$. In the hierarchical formation model, these massive galaxies $(10^{10.5} M_{\odot} < M_{dyn} < 10^{11.9} M_{\odot})$ are expected to undergo more mergers than lower mass galaxies, and simulations show that dry minor mergers can remove angular momentum. We test whether BCGs have low angular momenta by using the $\lambda_{Re}$ parameter developed by the SAURON and ATLAS\textsuperscript{3D} teams and combine our kinematics with Sloan Digital Sky Survey (SDSS) photometry to analyze the BCGs' merger status. We find that 30% (3/10) of the BCGs and 100% of the companion galaxies (4/4) are fast rotators as defined by the ATLAS\textsuperscript{3D} criteria. Our fastest rotating BCG has a $\lambda_{Re}=0.35\pm0.05$. We increase the number of BCGs analyzed from 1 in the combined SAURON and ATLAS\textsuperscript{3D} surveys to 11 BCGs total and find that above $M_{dyn}\sim11.5 M_{\odot}$, virtually all galaxies regardless of environment are slow rotators. To search for signs of recent merging, we analyze the photometry of each system and use the $G-M_{20}$ selection criteria to identify mergers. We find that $40\pm20$% of our BCGs are currently undergoing or have recently undergone a merger (within 0.2 Gyrs). Surprisingly, we find no correlation between galaxies with high angular momentum and morphological signatures of merging., Comment: 14 pages, 9 figures, 3 tables. Accepted for publication in ApJ

We show that, observationally, the projected local density distribution in high-z clusters is shifted towards higher values compared to clusters at lower redshift. To search for the origin of this evolution, we analyze a sample of haloes selected from the Millennium Simulation and populated using semi-analytic models, investigating the relation between observed projected density and physical 3D density, using densities computed from the 10 and 3 closest neighbours. Both observationally and in the simulations, we study the relation between number of cluster members and cluster mass, and number of members per unit of cluster mass. We find that the observed evolution of projected densities reflects a shift to higher values of the physical 3D density distribution. In turn, this must be related with the globally higher number of galaxies per unit of cluster volume N/V in the past. We show that the evolution of N/V is due to a combination of two effects: a) distant clusters were denser in dark matter (DM) simply because the DM density within R_{200} (~the cluster virial radius) is defined to be a fixed multiple of the critical density of the Universe, and b) the number of galaxies per unit of cluster DM mass is remarkably constant both with redshift and cluster mass if counting galaxies brighter than a passively evolving magnitude limit. Our results highlight that distant clusters were much denser environments than today's clusters, both in galaxy number and mass, and that the density conditions felt by galaxies in virialized systems do not depend on the system mass., accepted for publication in MNRAS

We investigate the recent and current star formation activity of galaxies as function of distance from the cluster center in a sample of 521 SDSS clusters at z0.1 R_200., 16 pages, 12 figures, submitted to MNRAS

We present the rest-frame optical luminosity function (LF) of red sequence galaxies in 16 clusters at 0.4, Comment: Accepted for publication in the Astrophysical Journal. 36 pages, 16 figures, 10 tables

We investigate how the [OII] properties and the morphologies of galaxies in clusters and groups at z=0.4-0.8 depend on projected local galaxy density, and compare with the field at similar redshifts and clusters at low-z. In both nearby and distant clusters, higher-density regions contain proportionally fewer star-forming galaxies, and the average [OII] equivalent width of star-forming galaxies is independent of local density. However, in distant clusters the average current star formation rate (SFR) in star-forming galaxies seems to peak at densities ~15-40 galaxies Mpc^{-2}. At odds with low-z results, at high-z the relation between star-forming fraction and local density varies from high- to low-mass clusters. Overall, our results suggest that at high-z the current star formation (SF) activity in star-forming galaxies does not depend strongly on global or local environment, though the possible SFR peak seems at odds with this conclusion. We find that the cluster SFR normalized by cluster mass anticorrelates with mass and correlates with the star-forming fraction. These trends can be understood given a) that the average star-forming galaxy forms about 1 Msun/yr in all clusters; b) that the total number of galaxies scales with cluster mass and c) the dependence of star-forming fraction on cluster mass. We present the morphology-density (MD) relation for our z=0.4-0.8 clusters, and uncover that the decline of the spiral fraction with density is entirely driven by galaxies of types Sc or later. For galaxies of a given Hubble type, we see no evidence that SF properties depend on local environment. In contrast with recent findings at low-z, in our distant clusters the SF-density relation and the MD-relation are equivalent, suggesting that neither of the two is more fundamental than the other.(abr.), 21 pages, 14 figures, accepted for publication in ApJ

We use the Sloan Digital Sky Survey to construct a sample of 625 brightest group and cluster galaxies (BCGs) together with control samples of non-BCGs matched in stellar mass, redshift, and color. We investigate how the systematic properties of BCGs depend on stellar mass and on their privileged location near the cluster center. The groups and clusters that we study are drawn from the C4 catalogue of Miller et al. (2005) but we have developed improved algorithms for identifying the BCG and for measuring the cluster velocity dispersion. Since the SDSS photometric pipeline tends to underestimate the luminosities of large galaxies in dense environments, we have developed a correction for this effect which can be readily applied to the published catalog data. We find that BCGs are larger and have higher velocity dispersions than non-BCGs of the same stellar mass, which implies that BCGs contain a larger fraction of dark matter. In contrast to non-BCGs, the dynamical mass-to-light ratio of BCGs does not vary as a function of galaxy luminosity. Hence BCGs lie on a different fundamental plane than ordinary elliptical galaxies. BCGs also follow a steeper Faber-Jackson relation than non-BCGs, as suggested by models in which BCGs assemble via dissipationless mergers along preferentially radial orbits. We find tentative evidence that this steepening is stronger in more massive clusters. BCGs have similar mean stellar ages and metallicities to non-BCGs of the same mass, but they have somewhat higher alpha/Fe ratios, indicating that star formation may have occurred over a shorter timescale in the BCGs. Finally, we find that BCGs are more likely to host radio-loud active galactic nuclei than other galaxies of the same mass, but are less likely to host an optical AGN. The differences we find are more pronounced for the less massive BCGs., Comment: Replaced with slightly modified version accepted by MNRAS. 28 pages, 25 figures. Version with full resolution figures available at http://www.mpa-garching.mpg.de/~anja/bcgs_avdl.pdf

We present the first study of the spatial distribution of star formation in z~0.5 cluster galaxies. The analysis is based on data taken with the Wide Field Camera 3 as part of the Grism Lens-Amplified Survey from Space (GLASS). We illustrate the methodology by focusing on two clusters (MACS0717.5+3745 and MACS1423.8+2404) with different morphologies (one relaxed and one merging) and use foreground and background galaxies as field control sample. The cluster+field sample consists of 42 galaxies with stellar masses in the range 10^8-10^11 M_sun, and star formation rates in the range 1-20 M_sun/yr. Both in clusters and in the field, H{\alpha} is more extended than the rest-frame UV continuum in 60% of the cases, consistent with diffuse star formation and inside out growth. In ~20% of the cases, the H{\alpha} emission appears more extended in cluster galaxies than in the field, pointing perhaps to ionized gas being stripped and/or star formation being enhanced at large radii. The peak of the H{\alpha} emission and that of the continuum are offset by less than 1 kpc. We investigate trends with the hot gas density as traced by the X-ray emission, and with the surface mass density as inferred from gravitational lens models and find no conclusive results. The diversity of morphologies and sizes observed in H_alpha illustrates the complexity of the environmental process that regulate star formation. Upcoming analysis of the full GLASS dataset will increase our sample size by almost an order of magnitude, verifying and strengthening the inference from this initial dataset., Comment: 18 pages, 15 figures, accepted for publication in ApJ

White paper for the "Dark Energy and CMB" working group for the American Physical Society's Division of Particles and Fields long-term planning exercise ("Snowmass"); International audience; Ongoing and near-future imaging-based dark energy experiments are critically dependent upon photometric redshifts (a.k.a. photo-z's): i.e., estimates of the redshifts of objects based only on flux information obtained through broad filters. Higher-quality, lower-scatter photo-z's will result in smaller random errors on cosmological parameters; while systematic errors in photometric redshift estimates, if not constrained, may dominate all other uncertainties from these experiments. The desired optimization and calibration is dependent upon spectroscopic measurements for secure redshift information; this is the key application of galaxy spectroscopy for imaging-based dark energy experiments. Hence, to achieve their full potential, imaging-based experiments will require large sets of objects with spectroscopically-determined redshifts, for two purposes: * Training: Objects with known redshift are needed to map out the relationship between object color and z (or, equivalently, to determine empirically-calibrated templates describing the rest-frame spectra of the full range of galaxies, which may be used to predict the color-z relation). The ultimate goal of training is to minimize each moment of the distribution of differences between photometric redshift estimates and the true redshifts of objects, making the relationship between them as tight as possible. The larger and more complete our "training set" of spectroscopic redshifts is, the smaller the RMS photo-z errors should be, increasing the constraining power of imaging experiments. Requirements: Spectroscopic redshift measurements for ∼30,000 objects over >∼15 widely-separated regions, each at least ∼20 arcmin in diameter, and reaching the faintest objects used in a given experiment, will likely be necessary if photometric redshifts are to be trained and calibrated with conventional techniques. Larger, more complete samples (i.e., with longer exposure times) can improve photo-z algorithms and reduce scatter further, enhancing the science return from planned experiments greatly (increasing the Dark Energy Task Force figure of merit by up to ∼50%). Options: This spectroscopy will most efficiently be done by covering as much of the optical and near-infrared spectrum as possible at modestly high spectral resolution (λ/Δλ > ∼3000), while maximizing the telescope collecting area, field of view on the sky, and multiplexing of simultaneous spectra. The most efficient instrument for this would likely be either the proposed GMACS/MANIFEST spectrograph for the Giant Magellan Telescope or the OPTIMOS spectrograph for the European Extremely Large Telescope, depending on actual properties when built. The PFS spectrograph at Subaru would be next best and available considerably earlier, c. 2018; the proposed ngCFHT and SSST telescopes would have similar capabilities but start later. Other key options, in order of increasing total time required, are the WFOS spectrograph at TMT, MOONS at the VLT, and DESI at the Mayall 4 m telescope (or the similar 4MOST and WEAVE projects); of these, only DESI, MOONS, and PFS are expected to be available before 2020. Table 3 of this white paper summarizes the observation time required at each facility for strawman training samples. To attain secure redshift measurements for a high fraction of targeted objects and cover the full redshift span of future experiments, additional near-infrared spectroscopy will also be required; this is best done from space, particularly with WFIRST-2.4 and JWST. Calibration: The first several moments of redshift distributions (the mean, RMS redshift dispersion, etc.), must be known to high accuracy for cosmological constraints not to be systematics-dominated (equivalently, the moments of the distribution of differences between photometric and true redshifts could be determined instead). The ultimate goal of calibration is to characterize these moments for every subsample used in analyses - i.e., to minimize the uncertainty in their mean redshift, RMS dispersion, etc. - rather than to make the moments themselves small. Calibration may be done with the same spectroscopic dataset used for training if that dataset is extremely high in redshift completeness (i.e., no populations of galaxies to be used in analyses are systematically missed). Accurate photo-z calibration is necessary for all imaging experiments. Requirements: If extremely low levels of systematic incompleteness (