Source author record

Rychard Bouwens

Rychard Bouwens appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2026arXiv

The early Universe with JWST and ALMA

The Atacama Large Millimeter/submillimeter Array and the James Webb Space Telescope are transforming our understanding of galaxy formation and evolution in the early Universe. By combining their capabilities, these observatories provide unprecedented insights into the gas, dust, and stars of high-redshift galaxies at spatially resolved scales, unveiling the complexities of their interstellar medium, kinematics, morphology, active galactic nuclei, and star formation activity. This review summarizes recent breakthroughs in the study of galaxies during the first billion years of cosmic history, highlighting key discoveries, open questions, and current limitations. We discuss how observations, theoretical models, and simulations are shaping our understanding of early galaxy evolution and identify promising directions for future research. While significant progress can be achieved through optimized use of existing facilities and collaborative efforts, further advances will require enhanced angular resolution and sensitivity, motivating upgrades to current instruments and the development of next-generation observatories.

preprint2022arXiv

ALMA as a Redshift Machine: Using [CII] to Efficiently Confirm Galaxies in the Epoch of Reionization

The [CII]$_{158μm}$ line has long been proposed as a promising line to spectroscopically confirm galaxies in the epoch of reionization. In this paper we present the results of new ALMA observations spectral scanning for [CII] in six particularly luminous Lyman Break Galaxies at $z\sim7$. The six sources were drawn from a sample of bright $z\sim7$ galaxies identified using the wide-area optical, near-IR, and Spitzer/IRAC data over the COSMOS/UltraVISTA field and were targeted on the basis of tight constraints on their redshifts from their IRAC [3.6]-[4.5] colors. We detect significant ($>9σ$) [CII] lines in three of our six targets ($50\%$) co-spatial with the rest-$UV$ emission from the ground/space-based near-IR imaging. The luminosities of the [CII] lines lie in the range $5.6$ to $8.8\times10^{8}L_{\odot}$, consistent with the local [CII]-SFR relation. Meanwhile, their [CII]/$L_{IR}\sim1-3\times10^{-3}$ ratios are slightly elevated compared to local (U)LIRGS. This could be due to lower dust-to-gas or dust-to-metal ratios. We also find that our sources display a large kinematic diversity, with one source showing signs of rotation, one source a likely major merger and one dispersion dominated source that might contain a bright star-forming clump. Our results highlight the effectiveness of spectral scans with ALMA in spectroscopically confirming luminous galaxies in the epoch of reionization, something that is being be applied on a significantly larger sample in the on-going REBELS large program.

preprint2022arXiv

Dual constraints with ALMA: new [O III] 88 ${\rm μ}$m and dust-continuum observations reveal the ISM conditions of luminous LBGs at $z \sim 7$

We present new [O III] 88 ${\rm μ}$m observations of five bright $z \sim 7$ Lyman-break galaxies spectroscopically confirmed by ALMA through the [C II] 158 ${\rm μ}$m line, unlike recent [O III] detections where Lyman-${\rm α}$ was used. This nearly doubles the sample of Epoch of Reionisation galaxies with robust ($5 σ$) detections of [C II] and [O III]. We perform a multi-wavelength comparison with new deep HST images of the rest-frame UV, whose compact morphology aligns well with [O III] tracing ionised gas. By contrast, we find more spatially extended [C II] emission likely produced in neutral gas, as indicated by a [N II] 205 ${\rm μ}$m non-detection in one source. We find a positive correlation between the equivalent width of the optical [O III] and H${\rm β}$ lines and the [O III]/[C II] ratio, as seen in local metal-poor dwarf galaxies. Cloudy models of a nebula of typical density harbouring a young stellar population with a high ionisation parameter appear to adequately reproduce the far-infrared lines. Surprisingly, however, our models fail to reproduce the strength of [O III] 88 ${\rm μ}$m, unless we assume an ${\rm α}$/Fe enhancement and a near-solar nebular oxygen abundance. On spatially resolved scales, we find [O III]/[C II] shows a tentative anti-correlation with infrared excess, $L_{\rm IR}/L_{\rm UV}$, also seen on global scales in the local Universe. Finally, we introduce the far-infrared spectral energy distribution fitting code MERCURIUS to show that dust-continuum measurements of one source appear to favour a low dust temperature coupled with a high dust mass. This implies a high stellar metallicity yield and may point towards the need of dust production or grain-growth mechanisms beyond supernovae.

preprint2022arXiv

Radio and Far-IR Emission Associated with a Massive Star-forming Galaxy Candidate at z$\simeq$6.8: A Radio-Loud AGN in the Reionization Era?

We report the identification of radio (0.144-3 GHz), mid-IR, far-IR, and sub-mm (24-850$μ$m) emission at the position of one of 41 UV-bright (M$_\mathrm{UV}^{}\lesssim-21.25$) $z\simeq6.6-6.9$ Lyman-break galaxy candidates in the 1.5 deg$^2$ COSMOS field. This source, COS-87259, exhibits a sharp flux discontinuity (factor $>$3) between two narrow/intermediate bands at 9450 and 9700 Angstroms and is undetected in all nine bands blueward of 9600 Angstroms, as expected from a Lyman-alpha break at $z\simeq6.8$. The full multi-wavelength (X-ray through radio) data of COS-87529 can be self-consistently explained by a very massive (M$_{\ast}=10^{10.8}$ M$_{\odot}$) and extremely red (rest-UV slope $β=-0.59$) $z\simeq6.8$ galaxy with hyperluminous infrared emission (L$_{\mathrm{IR}}=10^{13.6}$ L$_{\odot}$) powered by both an intense burst of highly-obscured star formation (SFR$\approx$1800 M$_{\odot}$ yr$^{-1}$) and an obscured ($τ_{\mathrm{9.7μm}}=7.7\pm2.5$) radio-loud (L$_{\mathrm{1.4\ GHz}}\approx10^{25.4}$ W Hz$^{-1}$) AGN. The radio emission is compact (1.04$\pm$0.12 arcsec) and exhibits an ultra-steep spectrum between 1.32-3 GHz ($α=-1.57^{+0.22}_{-0.21}$) that flattens at lower frequencies ($α=-0.86^{+0.22}_{-0.16}$ between 0.144-1.32 GHz), consistent with known $z>4$ radio galaxies. We also demonstrate that COS-87259 may reside in a significant (11$\times$) galaxy overdensity at $z\simeq6.6-6.9$, as common for systems hosting radio-loud AGN. Nonetheless, a spectroscopic redshift will ultimately be required to establish the true nature of COS-87259 as we cannot yet completely rule out low-redshift solutions. If confirmed to lie at $z\simeq6.8$, the properties of COS-87259 would be consistent with a picture wherein AGN and highly-obscured star formation activity are fairly common among very massive (M$_{\ast}>10^{10}$ M$_{\odot}$) reionization-era galaxies.

preprint2020arXiv

Elevated ionizing photon production efficiency in faint high-equivalent-width Lyman-alpha emitters

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic Reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z~4-5 continuum-faint Lyman-alpha emitters from deep VLT/MUSE spectroscopy and directly measure their Halpha emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than -16 (~0.01 L_star), implying a median Lyman-alpha equivalent width of 249 Angstroms. By combining the Halpha measurement with the UV magnitude we determine the ionizing photon production efficiency, xi_ion, a first for such faint galaxies. The measurement of log (xi_ion [Hz/erg]) = 26.28 (+0.28; -0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower-luminosity, Lyman-alpha-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4e-4 and 0.008, with light-weighted ages less than 3 Myr.

preprint2020arXiv

Reconstructing the observed ionizing photon production efficiency at z~2 using stellar population models

The ionizing photon production efficiency, $ξ_{ion}$, is a critical parameter that provides a number of physical constraints to the nature of the early Universe, including the contribution of galaxies to the timely completion of the reionization of the Universe. Here we use KECK/MOSFIRE and ZFOURGE multi-band photometric data to explore the $ξ_{ion}$ of a population of galaxies at $z\sim2$ with $log_{10}(M_*/M_\odot)\sim9.0-11.5$. Our 130 \Halpha\ detections show a median $log_{10}(ξ_{ion}[Hz/erg])$ of $24.8\pm0.5$ when dust corrected using a Calzetti et al. (2000) dust prescription. Our values are typical of mass/magnitude selected $ξ_{ion}$ values observed in the $z\sim2$ Universe. Using BPASSv2.2.1 and Starburst99 stellar population models with simple parametric star-formation-histories (SFH), we find that even with models that account for effects of stellar evolution with binaries/stellar rotation, model galaxies at $log_{10}(ξ_{ion}[Hz/erg])\lesssim25.0$ have low H$α$ equivalent widths (EW) and redder colors compared to our $z\sim2$ observed sample. We find that introducing star-bursts to the SFHs resolve the tension with the models, however, due to the rapid time evolution of $ξ_{ion}$, H$α$ EWs, and rest-frame optical colors, our Monte Carlo simulations of star-bursts show that random distribution of star-bursts in evolutionary time of galaxies are unlikely to explain the observed distribution. Thus, either our observed sample is specially selected based on their past SFH or stellar models require additional mechanisms to reproduce the observed high UV luminosity of galaxies for a given production rate of hydrogen ionizing photons.

preprint2020arXiv

The ALMA Spectroscopic Survey in the HUDF: Deep 1.2 mm continuum number counts

We present the results from the 1.2 mm continuum image obtained as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin$^2$ within a primary beam response of 50% (10%) and a rms value of $9.3\thinspace{\rmμJy\thinspace beam^{-1}}$. We detect 35 sources at high significance (Fidelity $\geq0.5$), 32 of these have well characterized near-infrared HST counterparts. We estimate the 1.2 mm number counts to flux levels of $<30\thinspace{\rmμJy}$ in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below $S_ν\sim 0.1$ mJy. In a second approach, we constrain the number counts by using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jack-knifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting $\sim$93% ($\sim$100% if we include the lower fidelity detections) of the total continuum dust emission associated to galaxies in the HUDF. The ancillary data allows us to study the dependence of the 1.2 mm number counts on redshift ($z=0-4$), galaxy dust mass (${\rm M}_{\rm dust}=10^{7}-10^{9}{\rm M}_{\odot}$), stellar mass (${\rm M}_{*}=10^{9}-10^{12}{\rm M}_{\odot}$), and star-formation rate (${\rm SFR}=1-1000\thinspace{\rm M}_{\odot}\thinspace{\rm yr^{-1}}$). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift.

preprint2020arXiv

Towards Studying Hierarchical Assembly in Real Time: A Milky Way Progenitor Galaxy at z = 2.36 under the Microscope

We use Hubble Space Telescope (HST) imaging and near-infrared spectroscopy from Keck/MOSFIRE to study the sub-structure around the progenitor of a Milky Way-mass galaxy in the Hubble Frontier Fields (HFF). Specifically, we study an $r_e = 40^{+70}_{-30}$pc, $M_{\star} \sim 10^{8.2} M_{\odot}$ rest-frame ultra-violet luminous "clump" at a projected distance of $\sim$100~pc from a $M_{\star} \sim 10^{9.8}$M$_{\odot}$ galaxy at $z = 2.36$ with a magnification $μ= 5.21$. We measure the star formation history of the clump and galaxy by jointly modeling the broadband spectral energy distribution from HST photometry and H$α$ from MOSFIRE spectroscopy. Given our inferred properties (e.g., mass, metallicity, dust) of the clump and galaxy, we explore scenarios in which the clump formed \emph{in-situ} (e.g., a star forming complex) or \emph{ex-situ} (e.g., a dwarf galaxy being accreted). If it formed \emph{in-situ}, we conclude that the clump is likely a single entity as opposed to a aggregation of smaller star clusters, making it one of the most dense star clusters cataloged. If it formed \emph{ex-situ}, then we are witnessing an accretion event with a 1:40 stellar mass ratio. However, our data alone are not informative enough to distinguish between \emph{in-situ} and \emph{ex-situ} scenarios to a high level of significance. We posit that the addition of high-fidelity metallicity information, such as [OIII]4363Å, which can be detected at modest S/N with only a few hours of JWST/NIRSpec time, may be a powerful discriminant. We suggest that studying larger samples of moderately lensed sub-structures across cosmic time can provide unique insight into the hierarchical formation of galaxies like the Milky Way.

preprint2018arXiv

The Reionization Lensing Cluster Survey (RELICS) and the Brightest High-z Galaxies

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.

preprint2016arXiv

ALMA spectroscopic survey in the Hubble Ultra Deep Field: CO luminosity functions and the evolution of the cosmic density of molecular gas

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Molecular gas reservoirs in high-redshift galaxies

We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Survey Description

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $μ$Jy beam$^{-1}$ and 12.7 $μ$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: The Infrared Excess of UV-selected z=2-10 galaxies as a function of UV-continuum Slope and Stellar Mass

We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%.

preprint2016arXiv

Lyman-alpha and CIII] Emission in z=7-9 Galaxies: Accelerated Reionization Around Luminous Star Forming Systems?

We discuss new Keck/MOSFIRE spectroscopic observations of four luminous galaxies at z~7-9 selected to have intense optical line emission by Roberts-Borsani et al. (2016). Previous follow-up has revealed Lyman-alpha in two of the four galaxies. Our new MOSFIRE observations confirm that Lyman-alpha is present in the entire sample. We detect Lyman-alpha emission in COS-zs7-1, confirming its redshift as z=7.154, and we detect Lyman-alpha in EGS-zs8-2 at z=7.477, verifying a tentative detection presented in an earlier study. The ubiquity of Lyman-alpha in this sample is puzzling given that the IGM is likely significantly neutral over 7<z<9. To investigate this result in more detail, we have initiated a campaign to target UV metal emission in the four Lyman-alpha emitters as a probe of both the radiation field and the velocity offset of Lyman-alpha. Here we present the detection of intense CIII] emission in EGS-zs8-1, a galaxy from this sample previously shown to have Lyman-alpha at z=7.73. Photoionization models indicate that an intense radiation field and low metallicity are required to reproduce the intense CIII] and optical line emission. We argue that this extreme radiation field is likely to affect the local environment, increasing the transmission of Lyman-alpha through the galaxy. Moreover, the centroid of CIII] indicates that Lyman-alpha is redshifted from the systemic value by 340 km/s. This velocity offset is larger than that seen in less luminous systems, providing an additional explanation for the transmission of Lyman-alpha emission through the IGM. Since the transmission is further enhanced by the likelihood that such systems are also situated in the densest regions with the largest ionized bubbles, the visibility of Lyman-alpha at z>7 is expected to be strongly luminosity-dependent, with the most effective transmission occurring in systems with intense star formation.

preprint2016arXiv

The production and escape of Lyman-Continuum radiation from star-forming galaxies at z~2 and their redshift evolution

We study the production rate of ionizing photons of a sample of 588 H$α$ emitters (HAEs) and 160 Lyman-$α$ emitters (LAEs) at $z=2.2$ in the COSMOS field in order to assess the implied emissivity from galaxies, based on their UV luminosity. By exploring the rest-frame Lyman Continuum (LyC) with GALEX/$NUV$ data, we find f$_{\rm esc} < 2.8\, (6.4)$% through median (mean) stacking. By combining the H$α$ luminosity density with IGM emissivity measurements from absorption studies, we find a globally averaged $\langle$f$_{\rm esc}\rangle$ of $5.9^{+14.5}_{-4.2}$ % at $z=2.2$ if we assume HAEs are the only source of ionizing photons. We find similarly low values of the global $\langle$f$_{\rm esc}\rangle$ at $z\approx3-5$, also ruling out a high $\langle$f$_{\rm esc}\rangle$ at $z<5$. These low escape fractions allow us to measure $ξ_{ion}$, the number of produced ionizing photons per unit UV luminosity, and investigate how this depends on galaxy properties. We find a typical $ξ_{ion} \approx 10^{24.77\pm0.04}$ Hz erg$^{-1}$ for HAEs and $ξ_{ion} \approx 10^{25.14\pm0.09}$ Hz erg$^{-1}$ for LAEs. LAEs and low mass HAEs at $z=2.2$ show similar values of $ξ_{ion}$ as typically assumed in the reionization era, while the typical HAE is three times less ionizing. Due to an increasing $ξ_{ion}$ with increasing EW(H$α$), $ξ_{ion}$ likely increases with redshift. This evolution alone is fully in line with the observed evolution of $ξ_{ion}$ between $z\approx2-5$, indicating a typical value of $ξ_{ion} \approx 10^{25.4}$ Hz erg$^{-1}$ in the reionization era.

preprint2015arXiv

CLASH-X: A Comparison of Lensing and X-ray Techniques for Measuring the Mass Profiles of Galaxy Clusters

We present profiles of temperature (Tx), gas mass, and hydrostatic mass estimated from new and archival X-ray observations of CLASH clusters. We compare measurements derived from XMM and Chandra observations with one another and compare both to gravitational lensing mass profiles derived with CLASH HST and ground-based lensing data. Radial profiles of Chandra and XMM electron density and enclosed gas mass are nearly identical, indicating that differences in hydrostatic masses inferred from X-ray observations arise from differences in Tx measurements. Encouragingly, cluster Txs are consistent with one another at ~100-200 kpc radii but XMM Tx systematically decline relative to Chandra Tx at larger radii. The angular dependence of the discrepancy suggests additional investigation on systematics such as the XMM point spread function correction, vignetting and off-axis responses. We present the CLASH-X mass-profile comparisons in the form of cosmology-independent and redshift-independent circular-velocity profiles. Ratios of Chandra HSE mass profiles to CLASH lensing profiles show no obvious radial dependence in the 0.3-0.8 Mpc range. However, the mean mass biases inferred from the WL and SaWLens data are different. e.g., the weighted-mean value at 0.5 Mpc is <b> = 0.12 for the WL comparison and <b> = -0.11 for the SaWLens comparison. The ratios of XMM HSE mass profiles to CLASH lensing profiles show a pronounced radial dependence in the 0.3-1.0 Mpc range, with a weighted mean mass bias of value rising to <b>~0.3 at ~1 Mpc for the WL comparison and <b> of 0.25 for SaWLens comparison. The enclosed gas mass profiles from both Chandra and XMM rise to a value 1/8 times the total-mass profiles inferred from lensing at 0.5 Mpc and remain constant outside of that radius, suggesting that [8xMgas] profiles may be an excellent proxy for total-mass profiles at >0.5 Mpc in massive galaxy clusters.

preprint2015arXiv

Lyman-alpha Emission from a Luminous z=8.68 Galaxy: Implications for Galaxies as Tracers of Cosmic Reionization

We report the discovery of Lyman-alpha emission (Ly$α$) in the bright galaxy EGSY-2008532660 (hereafter EGSY8p7) using the MOSFIRE spectrograph at the Keck Observatory. First reported by Roberts-Borsani et al. (2015), it was selected for spectroscopic observations because of its photometric redshift ($z_{phot}=8.57^{+0.22}_{-0.43}$), apparent brightness (H$_{160}=25.26\pm0.09$) and red Spitzer/IRAC [3.6]-[4.5] color indicative of contamination by strong oxygen emission in the [4.5] band. With a total integration of $\sim$4.3 hours, our data reveal an emission line at $\simeq$11776 Å which we argue is likely Ly$α$ at a redshift $z_{spec}=8.683^{+0.001}_{-0.004}$, in good agreement with the photometric estimate. The line was detected independently on two nights using different slit orientations and its detection significance is $\sim7.5σ$. An overlapping skyline contributes significantly to the uncertainty on the total line flux although the significance of the detected line is robust to a variety of skyline-masking procedures. By direct addition and a Gaussian fit, we estimate a 95\% confidence range of 1.0--2.5$\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$, corresponding to a rest-frame equivalent width of 17--42 Å. EGSY8p7 is the most distant galaxy confirmed spectroscopically to date, and the third luminous source in the EGS field beyond $z_{phot}\gtrsim7.5$ with detectable Ly$α$ emission viewed at a time when the intergalactic medium is believed to be fairly neutral. Although the reionization process was probably patchy, we discuss whether luminous sources with prominent IRAC color excesses may harbor harder ionizing spectra than the dominant fainter population thereby creating earlier ionized bubbles. Further spectroscopic follow-up of such bright sources promises important insight into the early formation of galaxies.

preprint2015arXiv

Not In Our Backyard: Spectroscopic Support for the CLASH z=11 Candidate MACS0647-JD

We report on our first set of spectroscopic Hubble Space Telescope observations of the z~11 candidate galaxy strongly lensed by the MACSJ0647.7+7015 galaxy cluster. The three lensed images are faint and we show that these early slitless grism observations are of sufficient depth to investigate whether this high-redshift candidate, identified by its strong photometric break at ~1.5 micron, could possibly be an emission line galaxy at a much lower redshift. While such an interloper would imply the existence of a rather peculiar object, we show here that such strong emission lines would clearly have been detected. Comparing realistic, two-dimensional simulations to these new observations we would expect the necessary emission lines to be detected at >5 sigma while we see no evidence for such lines in the dispersed data of any of the three lensed images. We therefore exclude that this object could be a low redshift emission line interloper, which significantly increases the likelihood of this candidate being a bona fide z~11 galaxy.

preprint2015arXiv

The Sizes of Candidate $z\sim9-10$ Galaxies: confirmation of the bright CANDELS sample and relation with luminosity and mass

Recently, a small sample of six $z\sim9-10$ candidates was discovered in CANDELS that are $\sim10-20\times$ more luminous than any of the previous $z\sim9-10$ galaxies identified over the HUDF/XDF and CLASH fields. We measure the sizes of these candidates to map out the size evolution of galaxies from the earliest observable times. Their sizes are also used to provide a valuable constraint on whether these unusual galaxy candidates are at high redshift. Using galfit to derive sizes from the CANDELS F160W images of these candidates, we find a mean size of 0.13$\pm$0.02" (or 0.5$\pm$0.1 kpc at $z\sim9-10$). This handsomely matches the 0.6 kpc size expected extrapolating lower redshift measurements to $z\sim9-10$, while being much smaller than the 0.59" mean size for lower-redshift interlopers to $z\sim9-10$ photometric selections lacking the blue IRAC color criterion. This suggests that source size may be an effective constraint on contaminants from $z\sim9-10$ selections lacking IRAC data. Assuming on the basis of the strong photometric evidence that the Oesch et al. 2014 sample is entirely at $z\sim9-10$, we can use this sample to extend current constraints on the size-luminosity, size-mass relation, and size evolution of galaxies to $z\sim10$. We find that the $z\sim9-10$ candidate galaxies have broadly similar sizes and luminosities as $z\sim6$-8 counterparts with star-formation-rate surface densities in the range of $\rm Σ_{SFR}=1-20\, M_\odot~ yr^{-1}\, kpc^{-2}$. The stellar mass-size relation is uncertain, but shallower than those inferred for lower-redshift galaxies. In combination with previous size measurements at z=4-7, we find a size evolution of $(1+z)^{-m}$ with $m=1.0\pm0.1$ for $>0.3L^*_{z=3}$ galaxies, consistent with the evolution previously derived from $2 < z < 8$ galaxies.

preprint2014arXiv

Three Gravitationally Lensed Supernovae Behind CLASH Galaxy Clusters

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CLO12Car (z = 1.28), SN CLN12Did (z = 0.85), and SN CLA11Tib (z = 1.14), are located behind three different clusters, MACSJ1720.2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and Abell 383 (z = 0.187), respectively. Each SN was detected in Hubble Space Telescope (HST) optical and infrared images. Based on photometric classification, we find that SNe CLO12Car and CLN12Did are likely to be Type Ia supernovae (SNe Ia), while the classification of SN CLA11Tib is inconclusive. Using multi-color light-curve fits to determine a standardized SN Ia luminosity distance, we infer that SN CLO12Car was approximately 1.0 +/- 0.2 mag brighter than field SNe Ia at a similar redshift and ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is approximately 0.2 +/- 0.2 mag brighter than field SNe Ia. We derive independent estimates of the predicted magnification from CLASH strong+weak lensing maps of the clusters: 0.83 +/- 0.16 mag for SN CLO12Car, 0.28 +/- 0.08 mag for SN CLN12Did, and 0.43 +/- 0.11 mag for SN CLA11Tib. The two SNe Ia provide a new test of the cluster lens model predictions: we find that the magnifications based on the SN Ia brightness and those predicted by the lens maps are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe Ia (from the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct determination of absolute magnifications.

preprint2013arXiv

CLASH: Complete Lensing Analysis of the Largest Cosmic Lens MACS J0717.5+3745 and Surrounding Structures

The galaxy cluster MACS J0717.5+3745 (z=0.55) is the largest known cosmic lens, with complex internal structures seen in deep X-ray, Sunyaev-Zel'dovich effect and dynamical observations. We perform a combined weak and strong lensing analysis with wide-field BVRi'z' Subaru/Suprime-Cam observations and 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). We find consistent weak distortion and magnification measurements of background galaxies, and combine these signals to construct an optimally estimated radial mass profile of the cluster and its surrounding large-scale structure out to 5 Mpc/h. We find consistency between strong-lensing and weak-lensing in the region where these independent data overlap, <500 kpc/h. The two-dimensional weak-lensing map reveals a clear filamentary structure traced by distinct mass halos. We model the lensing shear field with 9 halos, including the main cluster, corresponding to mass peaks detected above 2.5σ_κ. The total mass of the cluster as determined by the different methods is M_{vir}=(2.8\pm0.4) \times 10^15 M_sun. Although this is the most massive cluster known at z>0.5, in terms of extreme value statistics we conclude that the mass of MACS J0717.5+3745 by itself is not in serious tension with LambdaCDM, representing only a ~2σ departure above the maximum simulated halo mass at this redshift.

preprint2013arXiv

Galaxy halo truncation and Giant Arc Surface Brightness Reconstruction in the Cluster MACSJ1206.2-0847

In this work we analyze the mass distribution of MACSJ1206.2-0847, especially focusing on the halo properties of its cluster members. The cluster appears relaxed in its X-ray emission, but has significant amounts of intracluster light which is not centrally concentrated, suggesting that galaxy-scale interactions are still ongoing despite the overall relaxed state. The cluster lenses 12 background galaxies into multiple images and one galaxy at $z=1.033$ into a giant arc and its counterimage. The multiple image positions and the surface brightness distribution (SFB) of the arc which is bent around several cluster members are sensitive to the cluster galaxy halo properties. We model the cluster mass distribution with a NFW profile and the galaxy halos with two parameters for the mass normalization and extent of a reference halo assuming scalings with their observed NIR--light. We match the multiple image positions at an r.m.s. level of $0.85\arcsec$ and can reconstruct the SFB distribution of the arc in several filters to a remarkable accuracy based on this cluster model. The length scale where the enclosed galaxy halo mass is best constrained is about 5 effective radii -- a scale in between those accessible to dynamical and field strong lensing mass estimates on one hand and galaxy--galaxy weak lensing results on the other hand. The velocity dispersion and halo size of a galaxy with $m_{\rm 160W,AB}=19.2$ or $M_{\rm B,Vega}=-20.7$ are $σ=150 \rm kms^{-1}$ and $r\approx 26\pm 6 \rm kpc$, indicating that the halos of the cluster galaxies are tidally stripped. We also reconstruct the unlensed source (which is smaller by a factor of $\sim5.8$ in area), demonstrating the increase of morphological information due to lensing and conclude that this galaxy has likely star--forming spiral arms with a red (older) central component.

preprint2013arXiv

Slow Evolution of the Specific Star Formation Rate at z>2: The Impact of Dust, Emission Lines, and A Rising Star Formation History

We measure the evolution of the specific star formation rate (sSFR = SFR / Mstellar) between redshift 4 and 6 to investigate the previous reports of "constant" sSFR at z>2. We obtain photometry on a large sample of galaxies at z~4-6 located in the GOODS-S field that have high quality imaging from HST and Spitzer. We have derived stellar masses and star formation rates (SFRs) through stellar population modeling of their spectral energy distributions (SEDs). We estimate the dust extinction from the observed UV colors. In the SED fitting process we have studied the effects of assuming a star formation history (SFH) both with constant SFR and one where the SFR rises exponentially with time. The latter SFH is chosen to match the observed evolution of the UV luminosity function. We find that neither the mean SFRs nor the mean stellar masses change significantly when the rising SFR (RSF) model is assumed instead of the constant SFR model. When focusing on galaxies with Mstar ~ 5x10^9 Msun, we find that the sSFR evolves weakly with redshift (sSFR(z) \propto (1+z)^(0.6+/-0.1) Gyr^-1), consistent with previous results and with recent estimates of the sSFR at z~2-3 using similar assumptions. We have also investigated the impact of optical emission lines on our results. We estimate that the contribution of emission lines to the rest-frame optical fluxes is only modest at z~4 and 5 but it could reach ~50% at z~6. When emission lines of this strength are taken into account, the sSFR shows somewhat higher values at high redshifts, according to the relation sSFR(z) \propto (1+z)^(1.0+/-0.1) Gyr^-1, i.e., ~2.3x higher at z~6 than at z~2. However, the observed evolution is substantially weaker than that found at z<2 or that expected from current models (which corresponds to sSFR(z) \propto (1+z)^(2.5) Gyr^-1). -abridged-

preprint2013arXiv

The Contribution of Halos with Different Mass Ratios to the Overall Growth of Cluster-Sized Halos

We provide a new observational test for a key prediction of the ΛCDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $σ$ agreement with ΛCDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated.

preprint2013arXiv

Tracing the Mass Growth and Star Formation Rate Evolution of Massive Galaxies from z~6 to z~1 in the Hubble Ultra-Deep Field

We present an analysis of $\sim$1500 H160-selected photometric galaxies detected to a limiting magnitude of 27.8 in the HUDF, using imaging from the HST WFC3/IR camera in combination with archival UV, optical, and NIR imaging. We fit photometric redshifts and stellar population estimates for all galaxies with well-determined Spitzer IRAC fluxes, allowing for the determination of the cumulative mass function within the range $1<z<6$. By selecting samples of galaxies at a constant cumulative number density, we explore the co-evolution of stellar masses and star formation rates from z$\sim$6. We find a steady increase in the SFRs of galaxies at constant number density from z$\sim$6 to z$\sim$3. The peak epoch of star formation is found to shift to later times for galaxies with increasing number densities, in agreement with the expectations from cosmic downsizing. The observed SFRs can fully account for the mass growth to z$\sim$2 amongst galaxies with cumulative number densities greater than 10$^{-3.5}$ Mpc$^{-3}$. For galaxies with a lower constant number density we find the observed stellar masses are $\sim$3 times greater than that which may be accounted for by the observed star formation alone at late times, implying that growth from mergers plays an important role at $z<2$. We additionally observe a decreasing sSFR, equivalent to approximately one order of magnitude, from z$\sim$6 to z$\sim$2 amongst galaxies with number densities less than 10$^{-3.5}$ Mpc$^{-3}$ along with significant evidence that at any redshift the sSFR is higher for galaxies at higher number density. The combination of these findings can qualitatively explain the previous findings of a sSFR plateau at high redshift. Tracing the evolution of the fraction of quiescent galaxies for samples matched in cumulative number density over this redshift range, we find no unambiguous examples of quiescent galaxies at $z>4$.

preprint2012arXiv

A highly magnified candidate for a young galaxy seen when the Universe was 500 Myrs old

The early Universe at redshift z\sim6-11 marks the reionization of the intergalactic medium, following the formation of the first generation of stars. However, those young galaxies at a cosmic age of \lesssim 500 million years (Myr, at z \gtrsim 10) remain largely unexplored as they are at or beyond the sensitivity limits of current large telescopes. Gravitational lensing by galaxy clusters enables the detection of high-redshift galaxies that are fainter than what otherwise could be found in the deepest images of the sky. We report the discovery of an object found in the multi-band observations of the cluster MACS1149+22 that has a high probability of being a gravitationally magnified object from the early universe. The object is firmly detected (12 sigma) in the two reddest bands of HST/WFC3, and not detected below 1.2 μm, matching the characteristics of z\sim9 objects. We derive a robust photometric redshift of z = 9.6 \pm 0.2, corresponding to a cosmic age of 490 \pm 15Myr (i.e., 3.6% of the age of the Universe). The large number of bands used to derive the redshift estimate make it one of the most accurate estimates ever obtained for such a distant object. The significant magnification by cluster lensing (a factor of \sim15) allows us to analyze the object's ultra-violet and optical luminosity in its rest-frame, thus enabling us to constrain on its stellar mass, star-formation rate and age. If the galaxy is indeed at such a large redshift, then its age is less than 200 Myr (at the 95% confidence level), implying a formation redshift of zf \lesssim 14. The object is the first z>9 candidate that is bright enough for detailed spectroscopic studies with JWST, demonstrating the unique potential of galaxy cluster fields for finding highly magnified, intrinsically faint galaxies at the highest redshifts.

preprint2012arXiv

CLASH: Mass Distribution in and around MACS J1206.2-0847 from a Full Cluster Lensing Analysis

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dlnΣ/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

preprint2012arXiv

CLASH: Precise New Constraints on the Mass Profile of Abell 2261

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with ΛCDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent ΛCDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

preprint2012arXiv

CLASH: Three Strongly Lensed Images of a Candidate z ~ 11 Galaxy

We present a candidate for the most distant galaxy known to date with a photometric redshift z = 10.7 +0.6 / -0.4 (95% confidence limits; with z < 9.5 galaxies of known types ruled out at 7.2-sigma). This J-dropout Lyman Break Galaxy, named MACS0647-JD, was discovered as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). We observe three magnified images of this galaxy due to strong gravitational lensing by the galaxy cluster MACSJ0647.7+7015 at z = 0.591. The images are magnified by factors of ~8, 7, and 2, with the brighter two observed at ~26th magnitude AB (~0.15 uJy) in the WFC3/IR F160W filter (~1.4 - 1.7 um) where they are detected at >~ 12-sigma. All three images are also confidently detected at >~ 6-sigma in F140W (~1.2 - 1.6 um), dropping out of detection from 15 lower wavelength HST filters (~0.2 - 1.4 um), and lacking bright detections in Spitzer/IRAC 3.6um and 4.5um imaging (~3.2 - 5.0 um). We rule out a broad range of possible lower redshift interlopers, including some previously published as high redshift candidates. Our high redshift conclusion is more conservative than if we had neglected a Bayesian photometric redshift prior. Given CLASH observations of 17 high mass clusters to date, our discoveries of MACS0647-JD at z ~ 10.8 and MACS1149-JD1 at z ~ 9.6 are consistent with a lensed luminosity function extrapolated from lower redshifts. This would suggest that low luminosity galaxies could have reionized the universe. However given the significant uncertainties based on only two galaxies, we cannot yet rule out the sharp drop off in number counts at z >~ 10 suggested by field searches.

preprint2012arXiv

Cluster Lensing And Supernova survey with Hubble (CLASH): An Overview

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated.

preprint2011arXiv

Evolution of Galaxy Stellar Mass Functions, Mass Densities, and Mass to Light Ratios from z 7 to z 4

We derive stellar masses from SED fitting to rest-frame optical and UV fluxes for 401 star-forming galaxies at z 4, 5, and 6 from Hubble-WFC3/IR observations of the ERS combined with the deep GOODS-S Spitzer/IRAC data (and include a previously-published z 7 sample). A mass-luminosity relation with strongly luminosity-dependent M/Luv ratios is found for the largest sample (299 galaxies) at z 4. The relation M \propto L_{UV,1500}^(1.7+/-0.2) has a well-determined intrinsic sample variance of 0.5 dex. This relation is also consistent with the more limited samples at z 5-7. This z 4 mass-luminosity relation, and the well-established faint UV luminosity functions at z 4-7, are used to derive galaxy mass functions (MF) to masses M~10^8 at z 4-7. A bootstap approach is used to derive the MFs to account for the large scatter in the M--Luv relation and the luminosity function uncertainties, along with an analytical crosscheck. The MFs are also corrected for the effects of incompleteness. The incompleteness-corrected MFs are steeper than previously found, with slopes α_M-1.4 to -1.6 at low masses. These slopes are, however, still substantially flatter than the MFs obtained from recent hydrodynamical simulations. We use these MFs to estimate the stellar mass density (SMD) of the universe to a fixed M_{UV,AB}<-18 as a function of redshift and find a SMD growth \propto(1+z)^{-3.4 +/-0.8} from z 7 to z 4. We also derive the SMD from the completeness-corrected MFs to a mass limit M~10^{8} Msun. Such completeness-corrected MFs and the derived SMDs will be particularly important for model comparisons as future MFs reach to lower masses.

preprint2011arXiv

The Rest Frame UV to Optical Colors and SEDs of z~4-7 Galaxies

We use the ultra-deep HUDF09 and the deep ERS data from the HST WFC3/IR camera, along with the wide area Spitzer/IRAC data from GOODS-S to derive SEDs of star-forming galaxies from the rest-frame UV to the optical over a wide luminosity range (M_1500 ~ -21 to M_1500 ~ -18) from z ~ 7 to z ~ 4. The sample contains ~ 400 z ~ 4, ~ 120 z ~ 5, ~ 60 z ~ 6, and 36 prior z ~ 7 galaxies. Median stacking enables the first comprehensive study of very faint high-z galaxies at multiple redshifts (e.g., [3.6] = 27.4 +/- 0.1 AB mag for the M_1500 ~ -18 sources at z ~ 4). At z ~ 4 our faint median-stacked SEDs reach to ~ 0.06 L*(z=4) and are combined with recently published results at high luminosity L > L* that extend to M_1500 ~ -23. We use the observed SEDs and template fits to derive rest frame UV-to-optical colors (U - V) at all redshifts and luminosities. We find that this color does not vary significantly with redshift at a fixed luminosity. The UV-to-optical color does show a weak trend with luminosity, becoming redder at higher luminosities. This is most likely due to dust. At z >~ 5 we find blue colors [3.6]-[4.5] ~ -0.3 mag that are most likely due to rest-frame optical emission lines contributing to the flux in the IRAC filter bandpasses. The scatter across our derived SEDs remains substantial, but the results are most consistent with a lack of any evolution in the SEDs with redshift at a given luminosity. The similarity of the SEDs suggests a self-similar mode of evolution over a timespan from 0.7 Gyr to 1.5 Gyr that encompasses very substantial growth in the stellar mass density in the universe (from ~ 4x10^6 to ~ 2x10^7 Msun Mpc^-3).

preprint2008arXiv

Confirmation of the remarkable compactness of massive quiescent galaxies at z~2.3: early-type galaxies did not form in a simple monolithic collapse

Using deep near-infrared spectroscopy Kriek et al. (2006) found that ~45% of massive galaxies at z~2.3 have evolved stellar populations and little or no ongoing star formation. Here we determine the sizes of these quiescent galaxies using deep, high-resolution images obtained with HST/NIC2 and laser guide star-assisted Keck/AO. Considering that their median stellar mass is 1.7x10^11 Solar masses the galaxies are remarkably small, with a median effective radius of 0.9 kpc. Galaxies of similar mass in the nearby Universe have sizes of ~5 kpc and average stellar densities which are two orders of magnitude lower than the z~2.3 galaxies. These results extend earlier work at z~1.5 and confirm previous studies at z>2 which lacked spectroscopic redshifts and imaging of sufficient resolution to resolve the galaxies. Our findings demonstrate that fully assembled early-type galaxies make up at most ~10% of the population of K-selected quiescent galaxies at z~2.3, effectively ruling out simple monolithic models for their formation. The galaxies must evolve significantly after z~2.3, through dry mergers or other processes, consistent with predictions from hierarchical models.

preprint2006arXiv

Spitzer IRAC confirmation of z_850-dropout galaxies in the Hubble Ultra Deep Field: stellar masses and ages at z~7

Using Spitzer IRAC mid-infrared imaging from the Great Observatories Origins Deep Survey, we study z_850-dropout sources in the Hubble Ultra Deep Field. After carefully removing contaminating flux from foreground sources, we clearly detect two z_850-dropouts at 3.6 micron and 4.5 micron, while two others are marginally detected. The mid-infrared fluxes strongly support their interpretation as galaxies at z~7, seen when the Universe was only 750 Myr old. The IRAC observations allow us for the first time to constrain the rest-frame optical colors, stellar masses, and ages of the highest redshift galaxies. Fitting stellar population models to the spectral energy distributions, we find photometric redshifts in the range 6.7-7.4, rest-frame colors U-V=0.2-0.4, V-band luminosities L_V=0.6-3 x 10^10 L_sun, stellar masses 1-10 x 10^9 M_sun, stellar ages 50-200 Myr, star formation rates up to ~25 M_sun/yr, and low reddening A_V<0.4. Overall, the z=7 galaxies appear substantially less massive and evolved than Lyman break galaxies or Distant Red Galaxies at z=2-3, but fairly similar to recently identified systems at z=5-6. The stellar mass density inferred from our z=7 sample is rho* = 1.6^{+1.6}_{-0.8} x 10^6 M_sun Mpc^-3 (to 0.3 L*(z=3)), in apparent agreement with recent cosmological hydrodynamic simulations, but we note that incompleteness and sample variance may introduce larger uncertainties. The ages of the two most massive galaxies suggest they formed at z>8, during the era of cosmic reionization, but the star formation rate density derived from their stellar masses and ages is not nearly sufficient to reionize the universe. The simplest explanation for this deficiency is that lower-mass galaxies beyond our detection limit reionized the universe.