Large surveys of galaxy clusters with the Hubble Space Telescope (HST) and Spitzer, including the Cluster Lensing And Supernova survey with Hubble 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, described here. Our 188-orbit Hubble Treasury Program observed 41 clusters at 0.182 ≤ z ≤ 0.972 with Advanced Camera for Surveys (ACS) and WFC3/IR imaging spanning 0.4-1.7 m. We selected 21 of the most massive clusters known based on Planck PSZ2 estimates and 20 additional clusters based on observed or inferred lensing strength. RELICS observed 46 WFC3/IR pointings (∼200 arcmin2) each 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 over 300 z ∼ 6-10 candidates, including the brightest z ∼ 6 candidates known, and the most distant spatially resolved lensed arc known at z ∼ 10. Spitzer IRAC imaging (945 hr awarded, plus 100 archival, spanning 3.0-5.0 m) 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. Reduced HST images, catalogs, and lens models are available on MAST, and reduced Spitzer images are available on IRSA.
ABSTRACT We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity functions at z = 4-8. We use deep Hubble Space Telescope imaging over the Cosmic Assembly ...Near-infrared Deep Extragalactic Legacy Survey/GOODS fields, the Hubble Ultra Deep Field, and the Hubble Frontier Field deep parallel observations near the Abell 2744 and MACS J0416.1-2403 clusters. The combination of these surveys provides an effective volume of 0.6-1.2 × 106 Mpc3 over this epoch, allowing us to perform a robust search for faint 18) and bright (M 21) high-redshift galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 candidate galaxies at 3.5 8.5, with >1000 galaxies at 6-8. We measure both a stepwise luminosity function for candidate galaxies in our redshift samples, and a Schechter function, using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end, our UV luminosity functions agree with previous studies, yet we find a higher abundance of UV-bright candidate galaxies at 6. Our best-fit value of the characteristic magnitude is consistent with −21 at 5, which is different than that inferred based on previous trends at lower redshift, and brighter at ∼2 significance than previous measures at z = 6 and 7. At z = 8, a single power law provides an equally good fit to the UV luminosity function, while at z = 6 and 7 an exponential cutoff at the bright end is moderately preferred. We compare our luminosity functions to semi-analytical models, and find that the lack of evolution in is consistent with models where the impact of dust attenuation on the bright end of the luminosity function decreases at higher redshift, although a decreasing impact of feedback may also be possible. We measure the evolution of the cosmic star-formation rate (SFR) density by integrating our observed luminosity functions to , correcting for dust attenuation, and find that the SFR density declines proportionally to (1 ) at 4, which is consistent with observations at 9. Our observed luminosity functions are consistent with a reionization history that starts at 10, completes at 6, and reaches a midpoint (x 0.5) at 6.7 9.4. Finally, using a constant cumulative number density selection and an empirically derived rising star-formation history, our observations predict that the abundance of bright z = 9 galaxies is likely higher than previous constraints, although consistent with recent estimates of bright 10 galaxies.
We explore star formation histories (SFHs) of galaxies based on the evolution of the star formation rate stellar mass relation (SFR-M). Using data from the FourStar Galaxy Evolution Survey (ZFOURGE) ...in combination with far-IR imaging from the Spitzer and Herschel observatories we measure the SFR-M relation at 0.5 < z <. Similar to recent works we find that the average infrared spectral energy distributions of galaxies are roughly consistent with a single infrared template across a broad range of redshifts and stellar masses, with evidence for only weak deviations. We find that these two estimates are in broad qualitative agreement, but that there is room for improvement at a more detailed level. At early times the SFHs suggest mass growth rates that are as much as 10 x higher than inferred from the SMF. However, at later times the SFHs under-predict the inferred evolution, as is expected in the case of additional growth due to mergers.
Abstract
The study of galaxy evolution hinges on our ability to interpret multiwavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy ...distribution (SED) models, which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to wide and deep multiwave band galaxy surveys, the volume of high-quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED-fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply 14 of the most commonly used SED-fitting codes on samples from the CANDELS photometric catalogs at
z
∼ 1 and
z
∼ 3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust-attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust and active galactic nucleus models) on the derived stellar masses, SFRs, and
A
V
values. We then assess the difference among the codes on the SFR–stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (∼0.1 dex), SFR (∼0.3 dex), and dust attenuation (∼0.3 mag). Finally, we present some resources summarizing best practices in SED fitting.
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 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 arcmin2. These clusters were selected to be excellent lenses, and we find similar high-redshift sample sizes and magnitude distributions as the Cluster Lensing And Supernova survey with Hubble (CLASH). We discover 257, 57, and eight candidate galaxies at z ∼ 6, 7, and 8 respectively, (322 in total). The observed (lensed) magnitudes of the z ∼ 6 candidates are as bright as AB mag ∼23, making them among the brightest known at these redshifts, comparable with discoveries from 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.
This is the first in a series of papers examining the demographics of star-forming (SF) galaxies at 0.2 < z < 2.5 in CANDELS. We study 9100 galaxies from GOODS-S and UDS, having published values of ...redshifts, masses, star formation rates (SFRs), and dust attenuation (AV) derived from UV-optical spectral energy distribution fitting. In agreement with previous works, we find that the UVJ colors of a galaxy are closely correlated with its specific star formation rate (SSFR) and AV. We define rotated UVJ coordinate axes, termed SSED and CSED, that are parallel and perpendicular to the SF sequence and derive a quantitative calibration that predicts SSFR from CSED with an accuracy of ∼0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 m agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the SF main sequence is identified. These objects are located between SF and quiescent galaxies in UVJ space, and have lower AV and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and SF galaxies are given as a function of mass and redshift.
ABSTRACT Although there has been much progress in understanding how galaxies evolve, we still do not understand how and when they stop forming stars and become quiescent. We address this by applying ...our galaxy spectral energy distribution models, which incorporate physically motivated star formation histories (SFHs) from cosmological simulations, to a sample of quiescent galaxies at . A total of 845 quiescent galaxies with multi-band photometry spanning rest-frame ultraviolet through near-infrared wavelengths are selected from the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) data set. We compute median SFHs of these galaxies in bins of stellar mass and redshift. At all redshifts and stellar masses, the median SFHs rise, reach a peak, and then decline to reach quiescence. At high redshift, we find that the rise and decline are fast, as expected, because the universe is young. At low redshift, the duration of these phases depends strongly on stellar mass. Low-mass galaxies ( ) grow on average slowly, take a long time to reach their peak of star formation ( Gyr), and then the declining phase is fast ( Gyr). Conversely, high-mass galaxies ( ) grow on average fast ( Gyr), and, after reaching their peak, decrease the star formation slowly ( ). These findings are consistent with galaxy stellar mass being a driving factor in determining how evolved galaxies are, with high-mass galaxies being the most evolved at any time (i.e., downsizing). The different durations we observe in the declining phases also suggest that low- and high-mass galaxies experience different quenching mechanisms, which operate on different timescales.
Measurements of stellar properties of galaxies when the universe was less than one billion years old yield some of the only observational constraints on the onset of star formation. We present here ...the inclusion of Spitzer/IRAC imaging in the fitting of the spectral energy distribution of the seven highest-redshift galaxy candidates selected from the Hubble Space Telescope (HST) imaging of the Reionization Lensing Cluster Survey. We find that for six out of eight HST-selected z ∼ 8 sources, the z ∼ 8 solutions are still strongly preferred over z ∼ 1-2 solutions after the inclusion of Spitzer fluxes, and two prefer a z ∼ 7 solution, which we defer to a later analysis. We find a wide range of intrinsic stellar masses (5 × 106-4 × 109 M ), star formation rates (0.2-14 M yr−1), and ages (30-600 Myr) among our sample. Of particular interest is A1763-1434, which shows evidence of an evolved stellar population (∼500 Myr) at z ∼ 8, implying that its first generation of star formation occurred <100 Myr after the Big Bang. SPT0615-JD, a spatially resolved z ∼ 10 candidate, remains at its high redshift, supported by deep Spitzer/IRAC data, and also shows some evidence for an evolved stellar population. Even with the lensed, bright apparent magnitudes of these z 8 candidates (H = 26.1-27.8 AB mag), only the James Webb Space Telescope will be able to exclude the possibility of abnormally strong nebular emission, large dust content, or some combination thereof, and confirm the presence of evolved stellar populations early in the universe.
Strong gravitational lensing by galaxy clusters magnifies background galaxies, enhancing our ability to discover statistically significant samples of galaxies at , 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.
Abstract
We investigate the properties of galaxies as they shut off star formation over the 4 billion years surrounding peak cosmic star formation. To do this, we categorize ∼7000 galaxies from 1 <
...z
< 4 into 90 groups based on the shape of their spectral energy distributions (SEDs) and build composite SEDs with
R
∼ 50 resolution. These composite SEDs show a variety of spectral shapes and also show trends in parameters such as color, mass, star formation rate, and emission-line equivalent width. Using emission-line equivalent widths and strength of the 4000 Å break,
, we categorize the composite SEDs into five classes: extreme emission line, star-forming, transitioning, post-starburst, and quiescent galaxies. The transitioning population of galaxies shows modest H
α
emission (EW
REST
∼ 40 Å) compared to more typical star-forming composite SEDs at log
10
(
M
/
M
⊙
) ∼ 10.5 (EW
REST
∼ 80 Å). Together with their smaller sizes (3 kpc vs. 4 kpc) and higher Sérsic indices (2.7 vs. 1.5), this indicates that morphological changes initiate before the cessation of star formation. The transitional group shows a strong increase of over 1 dex in number density from
z
∼ 3 to
z
∼ 1, similar to the growth in the quiescent population, while post-starburst galaxies become rarer at
z
≲ 1.5. We calculate average quenching timescales of 1.6 Gyr at
z
∼ 1.5 and 0.9 Gyr at
z
∼ 2.5 and conclude that a fast-quenching mechanism producing post-starbursts dominated the quenching of galaxies at early times, while a slower process has become more common since
z
∼ 2.