Search

We present the first constraints on the prevalence of z>10 galaxies in the Hubble Ultra Deep Field (HUDF) leveraging new NIRCam observations from JEMS (JWST Extragalactic Medium-band Survey). These NIRCam observations probe redward of 1.6$μ$m, beyond the wavelength limit of HST, allowing us to search for galaxies to z>10. These observations indicate that the highest redshift candidate identified in the HUDF09 data with HST, UDFj-39546284, has a redshift of z>11.5, as had been suggested in analyses of the HUDF12/XDF data. This has now been confirmed with JWST NIRSpec. This source is thus the most distant galaxy discovered by HST in its >30 years of operation. Additionally, we identify nine other z~8-13 candidate galaxies over the HUDF, two of which are new discoveries that appear to lie at z~11-12. We use these results to characterize the evolution of the UV luminosity function (LF) from z~15 to z~8.7. While our LF results at z~8.7 and z~10.5 are consistent with previous findings over the HUDF, our new LF estimates at z~12.6 are higher than other results in the literature, potentially pointing to a milder evolution in the UV luminosity density from z~12.6. We emphasize that our LF results are uncertain given the small number of z~12.6 sources and limited volume probed. The new NIRCam data also indicate that the faint z~8-13 galaxies in the HUDF/XDF show blue UV-continuum slopes beta~-2.7, high specific star formation rates ~24.5 Gyr**-1, and high EW (~1300A) [OIII]+H$β$ emission, with two z~8.5 sources showing [OIII]+H$β$ EWs of ~2300 A., 21 pages, 14 figures, 7 tables, resubmitted to MNRAS, after replying to the referee report

We make use of ALMA continuum observations of $15$ luminous Lyman-break galaxies at $z$$\sim$$7$$-$$8$ to probe their dust-obscured star-formation. These observations are sensitive enough to probe to obscured SFRs of $20$ $M_{\odot}$$/$$yr$ ($3\sigma$). Six of the targeted galaxies show significant ($\geq$$3$$\sigma$) dust continuum detections, more than doubling the number of known dust-detected galaxies at $z$$>$$6.5$. Their IR luminosities range from $2.7$$\times$$10^{11}$ $L_{\odot}$ to $1.1$$\times$$10^{12}$ $L_{\odot}$, equivalent to obscured SFRs of $20$ to $105$ $M_{\odot}$$/$$yr$. We use our results to quantify the correlation of the infrared excess IRX on the UV-continuum slope $\beta_{UV}$ and stellar mass. Our results are most consistent with an SMC attenuation curve for intrinsic $UV$-slopes $\beta_{UV,intr}$ of $-2.63$ and most consistent with an attenuation curve in-between SMC and Calzetti for $\beta_{UV,intr}$ slopes of $-2.23$, assuming a dust temperature $T_d$ of $50$ K. Our fiducial IRX-stellar mass results at $z$$\sim$$7$$-$$8$ are consistent with marginal evolution from $z$$\sim$$0$. We then show how both results depend on $T_d$. For our six dust-detected sources, we estimate their dust masses and find that they are consistent with dust production from SNe if the dust destruction is low ($, ERC STFC

Deep rest-optical observations are required to accurately constrain the stellar populations of $z\sim8$ galaxies. Due to significant limitations in the availability of such data for statistically complete samples, observational results have been limited to modest numbers of bright or lensed sources. To revolutionize the present characterization of $z\sim8$ galaxies, we exploit the ultradeep ($\sim27$ mag, $3\sigma$) Spitzer/IRAC $3.6\mu$m and $4.5\mu$m data, probing the rest-frame optical at $z\sim8$, over $\sim200$ arcmin$^2$ of the GOODS fields from the recently completed GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) program, combined with observations in the CANDELS UDS and COSMOS fields. We stacked $\gtrsim100$ $z\sim8$ Lyman-Break galaxies in four bins of UV luminosity ($M_\mathrm{UV}\sim -20.7$ to $-18.4$) and study their $H_\mathrm{160}-[3.6]$ and $[3.6]-[4.5]$ colors. We find young ages ($\lesssim100$ Myr) for the three faintest stacks, inferred from their blue $H_\mathrm{160}-[3.6]\sim 0$ mag colors, consistent with a negative Balmer break. Meanwhile, the redder $H_\mathrm{160}-[3.6]$ color seen in the brightest stack is suggestive of slightly older ages. We explored the existence of a correlation between the UV luminosity and age, and find either no trend or fainter galaxies being younger. The stacked SEDs also exhibit very red $[3.6]-[4.5]\sim0.5$ mag colors, indicative of intense [OIII]+H$\beta$ nebular emission and SFR. The correspondingly high specific star-formation rates, sSFR$\gtrsim10$Gyr$^{-1}$, are consistent with recent determinations at similar redshifts and higher luminosities, and support the co-evolution between the sSFR and the specific halo mass accretion rate., Comment: Submitted to ApJ. Comments welcome. See Fig. 2 & 3 for the stacks, Fig. 4 & 6 for the colors/age and Fig. 8 for the sSFR evolution

The highly neutral inter-galactic medium (IGM) during the Epoch of Reionization (EoR) is expected to suppress Ly$\alpha$ emission with damping-wing absorption, causing nearly no Ly$\alpha$ detection from star-forming galaxies at $z{\sim}8$. However, spectroscopic observations of the 4 brightest galaxies (${\rm H}_{160}{\sim}25$ mag) at these redshifts do reveal prominent Ly$\alpha$ line, suggesting locally ionised IGM. In this paper, we explore the Ly$\alpha$ IGM transmission and environment of bright galaxies during the EoR using the Meraxes semi-analytic model. We find brighter galaxies to be less affected by damping-wing absorption as they are effective at ionizing surrounding neutral hydrogen. Specifically, the brightest sources (${\rm H}_{160}{\lesssim}25.5$ mag) lie in the largest ionized regions in our simulation, and have low attenuation of their Ly$\alpha$ from the IGM (optical depth ${, Comment: 10 pages, 3 figures, submitted to MNRAS, comments are welcome

Here we provide the most comprehensive determinations of the rest-frame $UV$ LF available to date with HST at z~2, 3, 4, 5, 6, 7, 8, and 9. Essentially all of the non-cluster extragalactic legacy fields are utilized, including the Hubble Ultra Deep Field (HUDF), the Hubble Frontier Field parallel fields, and all five CANDELS fields, for a total survey area of 1136 arcmin^2. Our determinations include galaxies at z~2-3 leveraging the deep HDUV, UVUDF, and ERS WFC3/UVIS observations available over a ~150 arcmin^2 area in the GOODS North and GOODS South regions. All together, our collective samples include >24,000 sources, >2.3x larger than previous selections with HST. 5766, 6332, 7240, 3449, 1066, 601, 246, and 33 sources are identified at z~2, 3, 4, 5, 6, 7, 8, and 9, respectively. Combining our results with an earlier z~10 LF determination by Oesch+2018a, we quantify the evolution of the $UV$ LF. Our results indicate that there is (1) a smooth flattening of the faint-end slope alpha from alpha~-2.4 at z~10 to -1.5 at z~2, (2) minimal evolution in the characteristic luminosity M* at z>~2.5, and (3) a monotonic increase in the normalization log_10 phi* from z~10 to z~2, which can be well described by a simple second-order polynomial, consistent with an "accelerated" evolution scenario. We find that each of these trends (from z~10 to z~2.5 at least) can be readily explained on the basis of the evolution of the halo mass function and a simple constant star formation efficiency model., 18 pages, 10 figures, 5 tables, in submission to ApJ, figures 9 and 10 show the main result

We present new stellar mass functions at z ∼ 6, z ∼ 7, z ∼ 8, z ∼ 9 and, for the first time, z ∼ 10, constructed from ∼800 Lyman-break galaxies previously identified over the eXtreme Deep Field and Hubble Ultra-Deep Field parallel fields and the five Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey fields. Our study is distinctive due to (1) the much deeper (∼200 hr) wide-area Spitzer/Infrared Array Camera (IRAC) imaging at 3.6 μm and 4.5 μm from the Great Observatories Origins Deep Survey Re-ionization Era Wide-area Treasury from Spitzer program (GREATS) and (2) consideration of z ∼ 6–10 sources over a 3× larger area than those of previous Hubble Space Telescope+Spitzer studies. The Spitzer/IRAC data enable ≥2σ rest-frame optical detections for an unprecedented 50% of galaxies down to a stellar mass limit of ∼ 10 8  ⊙ across all redshifts. Schechter fits to our volume densities suggest a combined evolution in the characteristic mass  * and normalization factor ϕ * between z ∼ 6 and z ∼ 8. The stellar mass density (SMD) increases by ∼1000× in the ∼500 Myr between z ∼ 10 and z ∼ 6, with indications of a steeper evolution between z ∼ 10 and z ∼ 8, similar to the previously reported trend of the star formation rate density. Strikingly, abundance matching to the Bolshoi–Planck simulation indicates halo mass densities evolving at approximately the same rate as the SMD between z ∼ 10 and z ∼ 4. Our results show that the stellar-to-halo mass ratios, a proxy for the star formation efficiency, do not change significantly over the huge stellar mass buildup occurred from z ∼ 10 to z ∼ 6, indicating that the assembly of stellar mass closely mirrors the buildup in halo mass in the first ∼1 Gyr of cosmic history. The James Webb Space Telescope is poised to extend these results into the “first galaxy” epoch at z ≳ 10.

We present the deepest Spitzer/IRAC $3.6$, $4.5$, $5.8$ and $8.0\mu$m wide-area mosaics yet over the GOODS-N and GOODS-S fields as part of the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) project. We reduced and mosaicked in a self-consistent way observations taken by the 11 different Spitzer/IRAC programs over the two GOODS fields from 12 years of Spitzer cryogenic and warm mission data. The cumulative depth in the $3.6\mu$m and $4.5\mu$m bands amounts to $\sim 4260$ hr, $\sim 1220$ hr of which are new very deep observations from the GREATS program itself. In the deepest area, the full-depth mosaics reach $\gtrsim200$ hr over an area of $\sim100$ arcmin$^2$, corresponding to a sensitivity of $\sim29$ AB magnitude at $3.6\mu$m ($1\sigma$ for point sources). Archival cryogenic $5.8\mu$m and $8.0\mu$m band data (a cumulative 976 hr) are also included in the release. The mosaics are projected onto the tangential plane of CANDELS/GOODS at a $0.3''$ pixel$^{-1}$ scale. This paper describes the methodology enabling, and the characteristics of, the public release of the mosaic science images, the corresponding coverage maps in the four IRAC bands, and the empirical Point-Spread Functions (PSFs). These PSFs enable mitigation of the source blending effects by taking into account the complex position-dependent variation in the IRAC images. The GREATS data products are in the Infrared Science Archive (IRSA). We also release the deblended $3.6$-to-$8.0\mu$m photometry for $9192$ Lyman-Break galaxies at $z\sim3.5-10$. GREATS will be the deepest mid-infrared imaging until JWST and, as such, constitutes a major resource for characterizing early galaxy assembly., Comment: Accepted for publication on ApJS. Mosaics, PSFs and catalogs are publicly available at https://home.strw.leidenuniv.nl/~stefanon/GREATS/

We present the Super Eight galaxies - a set of very luminous, high-redshift ($7.1, Comment: 21 pages, 17 figures, accepted for publication in ApJ

We report the results of an expanded search for z~9-10 candidates over the ~883 arcmin^2 CANDELS+ERS fields. This study adds 147 arcmin^2 to the search area we consider over the CANDELS COSMOS, UDS, and EGS fields, while expanding our selection to include sources with bluer J_{125}-H_{160} colors than our previous J_{125}-H_{160}>0.5 mag selection. In searching for new z~9-10 candidates, we make full use of all available HST, Spitzer/IRAC, and ground-based imaging data. As a result of our expanded search and use of broader color criteria, 3 new candidate z~9-10 galaxies are identified. We also find again the z=8.683 source previously confirmed by Zitrin+2015. This brings our sample of probable z~9-11 galaxy candidates over the CANDELS+ERS fields to 19 sources in total, equivalent to 1 candidate per 47 arcmin^2 (1 per 10 WFC3/IR fields). To be comprehensive, we also discuss 28 mostly lower likelihood z~9-10 candidates, including some sources that seem to be reliably at z>8 using the HST+IRAC data alone, but which the ground-based data show are much more likely at z8, with a factor of ~10 evolution seen in the luminosity density from z~10 to z~8., 22 pages, 12 figures, 6 tables, accepted for publication in the Astrophysical Journal

High-resolution multi-wavelength photometry is crucial to explore the spatial distribution of star formation in galaxies and understand how these evolve. To this aim, in this paper we exploit the deep, multi-wavelength Hubble Space Telescope (HST) data available in the central parts of the GOODS fields and study the distribution of star formation activity and mass in galaxies located at different positions with respect to the Main Sequence (MS) of star-forming galaxies. Our sample consists of galaxies with stellar mass $\geq 10^{9.5} M_{\odot}$ in the redshift range 0.2 $ \leq z \leq 1.2$. Exploiting 10-band photometry from the UV to the near-infrared at HST resolution, we derive spatially resolved maps of galaxies properties, such as stellar mass and star formation rate and specific star formation rate, with a resolution of $\sim 0.16$ arcsec. We find that the star formation activity is centrally enhanced in galaxies above the MS and centrally suppressed below the MS, with quiescent galaxies (1 dex below the MS) characterised by the highest suppression. The sSFR in the outer region does not show systematic trends of enhancement or suppression above or below the MS. The distribution of mass in MS galaxies indicates that bulges are growing when galaxies are still on the MS relation. Galaxies below the MS are more bulge-dominated with respect to MS counterparts at fixed stellar mass, while galaxies in the upper envelope are more extended and have S\'ersic indexes that are always smaller than or comparable to their MS counterparts. The suppression of star formation activity in the central region of galaxies below the MS hints at \textit{inside-out} quenching, as star formation is still ongoing in the outer regions., Comment: 17 pages, 12 figures, accepted for publication on A&A

We use a newly assembled large sample of 3,545 star-forming galaxies with secure spectroscopic, grism, and photometric redshifts at z=1.5-2.5 to constrain the relationship between UV slope (beta) and dust attenuation (L(IR)/L(UV)=IRX). Our sample benefits from the combination of deep Hubble WFC3/UVIS photometry from the Hubble Deep UV (HDUV) Legacy survey and existing photometric data compiled in the 3D-HST survey, and extends the range of UV luminosity and beta probed in previous UV-selected samples. IRX is measured using stacks of deep Herschel/PACS 100 and 160 micron data, and the results are compared with predictions of the IRX-beta relation for different assumptions of the stellar population model and obscuration curve. We find that z=1.5-2.5 galaxies have an IRX-beta relation that is consistent with the predictions for an SMC extinction curve if we invoke sub-solar metallicity models that are currently favored for high-redshift galaxies, while the commonly assumed starburst attenuation curve over-predicts the IRX at a given beta by a factor of ~3. The IRX of high-mass (M*>10^9.75 Msun) galaxies is a factor of >4 larger than that of low-mass galaxies, lending support for the use of stellar mass as a proxy for attenuation. The commonly observed trend of fainter galaxies having bluer beta may simply reflect bluer intrinsic UV slopes for such galaxies, rather than lower obscurations. The IRX-beta for young/low-mass galaxies implies a dust curve that is steeper than the SMC, suggesting a lower attenuation at a given beta relative to older/more massive galaxies. The lower attenuations and higher ionizing photon output implied by low metallicity stellar population models point to Lyman continuum production efficiencies, xi_ion, that may be elevated by a factor of ~2 relative to the canonical value for L* galaxies, aiding in their ability to keep the universe ionized at z~2. [Abridged], 22 pages, 3 tables, 11 figures, submitted to the Astrophysical Journal

The earliest galaxies are expected to emerge in the first billion years of the Universe during the Epoch of Reionization. However, both the spectroscopic confirmation of photometrically-selected galaxies at this epoch and the characterization of their early dynamical state has been hindered by the lack of bright, accessible lines to probe the velocity structure of their interstellar medium. We present the first ALMA spectroscopic confirmation of such sources at z > 6 using the far-infrared [C II]{\lambda}157.74{\mu}m emission line, and, for the first time, measurement of the velocity structure, for two galaxies at z = 6.8540+/-0.0003 and z = 6.8076+/-0.0002. Remarkably, the [C II] line luminosity from these galaxies is higher than previously found in `normal' star-forming galaxies at z > 6.5. This suggests that we are sampling a part of the galaxy population different from the galaxies found through detection of the Ly{\alpha} line. The luminous and extended [C II] detections reveal clear velocity gradients that, if interpreted as rotation, would suggest these galaxies have turbulent, yet rotation-dominated disks, with similar stellar-to-dynamical mass fractions as observed for H{\alpha} emitting galaxies 2 Gyr later at cosmic noon. Our novel approach for confirming galaxies during Reionization paves the way for larger studies of distant galaxies with spectroscopic redshifts from ALMA. Particularly important, this opens up opportunities for high angular-resolution [C II] dynamics in galaxies less than one billion years after the Big Bang., Comment: 14 pages, 5 figures, 1 table, published in Nature (submitted version)

We investigate the dependence of galaxy clustering at $z \sim 4 - 7$ on UV-luminosity and stellar mass. Our sample consists of $\sim$ 10,000 Lyman-break galaxies (LBGs) in the XDF and CANDELS fields. As part of our analysis, the $M_\star - M_{\rm UV}$ relation is estimated for the sample, which is found to have a nearly linear slope of $d\log_{10} M_\star / d M_{\rm UV} \sim 0.44$. We subsequently measure the angular correlation function and bias in different stellar mass and luminosity bins. We focus on comparing the clustering dependence on these two properties. While UV-luminosity is only related to recent starbursts of a galaxy, stellar mass reflects the integrated build-up of the whole star formation history, which should make it more tightly correlated with halo mass. Hence, the clustering segregation with stellar mass is expected to be larger than with luminosity. However, our measurements suggest that the segregation with luminosity is larger with $\simeq 90\%$ confidence (neglecting contributions from systematic errors). We compare this unexpected result with predictions from the \textsc{Meraxes} semi-analytic galaxy formation model. Interestingly, the model reproduces the observed angular correlation functions, and also suggests stronger clustering segregation with luminosity. The comparison between our observations and the model provides evidence of multiple halo occupation in the small scale clustering., Comment: 10 pages, 6 figures, 2 tables, accepted for publication in MNRAS

We present the Hubble Deep UV Legacy Survey (HDUV), a 132 orbit imaging program with the WFC3/UVIS camera onboard the Hubble Space Telescope (HST). The HDUV extends and builds on the few previous UV imaging surveys in the two GOODS/CANDELS-Deep fields to provide deep images over a total area of ~100 arcmin2 in the two filters F275W and F336W. Our release also includes all the F275W imaging data taken by the CANDELS survey, which were aligned using a novel approach and combined with the HDUV survey data. By reaching depths of 27.5-28.0 mag (5sigma, in 0.4" apertures), these are the deepest high-resolution UV data over such a large area taken to date. Such unique UV imaging enables a wide range of science by the community. Among the main goals of the HDUV survey are: (1) provide a complete sample of faint star-forming galaxies at z~1-3, (2) constrain the ionizing photon escape fraction from galaxies at z~2-3, and (3) track the build-up of bulges and the disappearance of clumpy disk galaxies through reliable internal stellar population properties at sub-kpc resolution out to z~3. The addition of the HDUV data further enhances the legacy value of the two GOODS/CANDELS-Deep fields, which now include deep 11-band HST imaging as well as very deep ancillary data from X-ray to radio, enabling unique multi-wavelength studies. Here, we provide an overview of the survey design, describe the data reduction, and highlight a few basic analyses on the images which are released to the community as high level science products via the Mikulski Archive for Space Telescopes (MAST)., Comment: 11 pages, 10 figures, accepted for publication in ApJS

We present six galaxies at z~2 that show evidence of Lyman continuum (LyC) emission based on the newly acquired UV imaging of the Hubble Deep UV legacy survey (HDUV) conducted with the WFC3/UVIS camera on the Hubble Space Telescope (HST). At the redshift of these sources, the HDUV F275W images partially probe the ionizing continuum. By exploiting the HST multi-wavelength data available in the HDUV/GOODS fields, models of the UV spectral energy distributions, and detailed Monte Carlo simulations of the intergalactic medium absorption, we estimate the absolute ionizing photon escape fractions of these galaxies to be very high -- typically >60% (>13% for all sources at 90% likelihood). Our findings are in broad agreement with previous studies that found only a small fraction of galaxies to show high escape fraction. These six galaxies comprise the largest sample yet of LyC leaking candidates at z~2 whose inferred LyC flux has been cleanly observed at HST resolution. While three of our six candidates show evidence of hosting an active galactic nucleus (AGN), two of these are heavily obscured and their LyC emission appears to originate from star-forming regions rather than the central nucleus. This suggests an AGN-aided pathway for LyC escape from these sources. Extensive multi-wavelength data in the GOODS fields, especially the near-IR grism spectra from the 3D-HST survey, enable us to study the candidates in detail and tentatively test some recently proposed indirect methods to probe LyC leakage -- namely, the [OIII]/[OII] line ratio and the H$��-$UV slope diagram. High-resolution spectroscopic followup of our candidates will help constrain such indirect methods which are our only hope of studying $f_{esc}$ at z~5-9 in the fast-approaching era of the James Webb Space Telescope., 11 pages, 6 figures, submitted to ApJ

We use the largest sample of z~6 galaxies to date from the first four Hubble Frontier Fields clusters to set constraints on the shape of the z~6 luminosity functions (LFs) to fainter than Muv=-14 mag. We quantify, for the first time, the impact of magnification uncertainties on LF results and thus provide more realistic constraints than other recent work. Our simulations reveal that for the highly-magnified sources the systematic uncertainties can become extremely large fainter than -14 mag, reaching several orders of magnitude at 95% confidence at ~-12 mag. Our new forward-modeling formalism incorporates the impact of magnification uncertainties into the LF results by exploiting the availability of many independent magnification models for the same cluster. One public magnification model is used to construct a mock high-redshift galaxy sample that is then analyzed using the other magnification models to construct a LF. Large systematic errors occur at high magnifications (mu>30) because of differences between the models. The volume densities we derive for faint (>-17 mag) sources are ~3-4x lower than one recent report and give a faint-end slope alpha=-1.92+/-0.04, which is 3.0-3.5sigma shallower (including or not including the size uncertainties, respectively). We introduce a new curvature parameter delta to model the faint end of the LF and demonstrate that the observations permit (at 68% confidence) a turn-over at z~6 in the range -15.3 to -14.2 mag, depending on the assumed lensing model. The present consideration of magnification errors and new size determinations raise doubts about previous reports regarding the form of the LF at >-14 mag. We discuss the implications of our turn-over constraints in the context of recent theoretical predictions., 33 pages, 19 figures, 5 tables, accepted for publication in the Astrophysical Journal, updated to match the version in press

We provide the first observational constraints on the sizes of the faintest galaxies lensed by the Hubble Frontier Fields (HFF) clusters. Ionizing radiation from faint galaxies likely drives cosmic reionization, and the HFF initiative provides a key opportunity to find such galaxies. Yet, we cannot really assess their ionizing emissivity without a robust measurement of their sizes, since this is key to quantifying both their prevalence and the faint-end slope to the UV luminosity function. Here we provide the first such size constraints with 2 new techniques. The first utilizes the fact that the detectability of highly-magnified galaxies as a function of shear is very dependent on a galaxy's size. Only the most compact galaxies will remain detectable in regions of high shear (vs. a larger detectable size range for low shear), a phenomenon we carefully quantify using simulations. Remarkably, however, no correlation is found between the surface density of faint galaxies and the predicted shear, using 87 faint high-magnification mu>10 z~2-8 galaxies seen behind the first 4 HFF clusters. This can only be the case if such faint (~-15 mag) galaxies have significantly smaller sizes than luminous galaxies. We constrain their half-light radii to be ~0.1-0.3)., Comment: 19 pages, 15 figures, 3 tables, accepted for publication in ApJ

The deep, wide-area (~800-900 arcmin**2) near-infrared/WFC3/IR + Spitzer/IRAC observations over the CANDELS fields have been a remarkable resource for constraining the bright end of high redshift UV luminosity functions (LFs). However, the lack of HST 1.05-micron observations over the CANDELS fields has made it difficult to identify z~9-10 sources robustly, since such data are needed to confirm the presence of an abrupt Lyman break at 1.2 microns. We report here on the successful identification of many such z~9-10 sources from a new HST program (z9-CANDELS) that targets the highest-probability z~9-10 galaxy candidates with observations at 1.05 microns, to search for a robust Lyman-break at 1.2 microns. The potential z~9-10 candidates are preselected from the full HST, Spitzer/IRAC S-CANDELS observations, and the deepest-available ground-based optical+near-infrared observations. We identified 15 credible z~9-10 galaxies over the CANDELS fields. Nine of these galaxies lie at z~9 and 5 are new identifications. Our targeted follow-up strategy has proven to be very efficient in making use of scarce HST time to secure a reliable sample of z~9-10 galaxies. Through extensive simulations, we replicate the selection process for our sample (both the preselection and follow-up) and use it to improve current estimates for the volume density of bright z~9 and z~10 galaxies. The volume densities we find are 5(-2)(+3)x and 8(-3)(+9)x lower, respectively, than found at z~8. When compared with the best-fit evolution (i.e., dlog_{10} rho(UV)/dz=-0.29+/-0.02) in the UV luminosities densities from z~8 to z~4 integrated to 0.3L*(z=3) (-20 mag), these luminosity densities are 2.6(-0.9)(+1.5)x and 2.2(-1.1)(+2.0)x lower, respectively, than the extrapolated trends. Our new results are broadly consistent with the "accelerated evolution" scenario at z>8, as seen in many theoretical models., 23 pages, 15 figures, 7 tables, updated to match the version in press, including some minor textual corrections identified at the proof stage

We present the details and early results from a deep near-infrared survey utilising the NICMOS instrument on the Hubble Space Telescope centred around massive M_* > 10^11 M_0 galaxies at 1.7 10^11 M_0, whereby we find an increase of a factor of eight between z = 3 and z = 1.5, demonstrating that this is an epoch when massive galaxies establish most of their mass.

[Abridged] We use the NEWFIRM Medium-Band Survey (NMBS) to characterize the properties of a mass-complete sample of 14 galaxies at 3.02.5x10^11 Msun, and to derive more accurate measurements of the high-mass end of the stellar mass function (SMF) of galaxies at z=3.5, with significantly reduced contributions from photometric redshift errors and cosmic variance to the total error budget of the SMF. The typical very massive galaxy at z=3.5 is red and faint in the observer's optical, with median r=26.1, and rest-frame U-V=1.6. About 60% of the sample have optical colors satisfying either the U- or the B-dropout color criteria, although ~50% of these galaxies have r>25.5. About 30% of the sample has SFRs from SED modeling consistent with zero. However, >80% of the sample is detected at 24 micron, with total infrared luminosities in the range (0.5-4.0)x10^13 Lsun. This implies the presence of either dust-enshrouded starburst activity (with SFRs of 600-4300 Msun/yr) and/or highly-obscured active galactic nuclei (AGN). The contribution of galaxies with M_star>2.5x10^11 Msun to the total stellar mass budget at z=3.5 is ~8%. We find an evolution by a factor of 2-7 and 3-22 from z~5 and z~6, respectively, to z=3.5. The previously found disagreement at the high-mass end between observed and model-predicted SMFs is now significant at the 3sigma level. However, systematic uncertainties dominate the total error budget, with errors up to a factor of ~8 in the densities, bringing the observed SMF in marginal agreement with the predicted SMF. Additional systematic uncertainties on the high-mass end could be introduced by either 1) the intense star-formation and/or the very common AGN activities as inferred from the MIPS 24 micron detections, and/or 2) contamination by a significant population of massive, old, and dusty galaxies at z~2.6., Comment: 20 pages, 11 figures. Accepted in ApJ. Minor changes to colors of figures to match accepted version

We investigate the Spitzer/IRAC properties of 36 z ~ 7 z 850-dropout galaxies and three z ~ 8 Y 098 galaxies derived from deep/wide-area WFC3/IR data of the Early Release Science, the ultradeep HUDF09, and wide-area NICMOS data. We fit stellar population synthesis models to the spectral energy distributions to derive mean redshifts, stellar masses, and ages. The z ~ 7 galaxies are best characterized by substantial ages (>100 Myr) and M/LV 0.2. The main trend with decreasing luminosity is that of bluing of the far-UV slope from ? ~ ?2.0 to ? ~ ?3.0. This can be explained by decreasing metallicity, except for the lowest luminosity galaxies (0.1L* z = 3), where low metallicity and smooth star formation histories (SFHs) fail to match the blue far-UV and moderately red H ? [3.6] color. Such colors may require episodic SFHs with short periods of activity and quiescence (on-off cycles) and/or a contribution from emission lines. The stellar mass of our sample of z ~ 7 star-forming galaxies correlates with star formation rate (SFR) according to log M* = 8.70(?0.09) + 1.06(?0.10)log SFR, implying that star formation may have commenced at z > 10. No galaxies are found with SFRs much higher or lower than the past averaged SFR suggesting that the typical star formation timescales are probably a substantial fraction of the Hubble time. We report the first IRAC detection of Y 098-dropout galaxies at z ~ 8. The average rest-frame U ? V 0.3 (AB) of the three galaxies are similar to faint z ~ 7 galaxies, implying similar M/L. The stellar mass density to M UV,AB < ?18 is ?*(z = 8) = 1.8+0.7 ?1.0 ? 106 M ? Mpc?3, following log ?*(z) = 10.6(?0.6) ? 4.4(?0.7) log(1 + z) [M ? Mpc?3] over 3 < z < 8.

We present an analysis of all prime HST legacy fields spanning >800 arcmin^2 for the search of z~10 galaxy candidates and the study of their UV luminosity function (LF). In particular, we present new z~10 candidates selected from the full Hubble Frontier Field (HFF) dataset. Despite the addition of these new fields, we find a low abundance of z~10 candidates with only 9 reliable sources identified in all prime HST datasets that include the HUDF09/12, the HUDF/XDF, all the CANDELS fields, and now the HFF survey. Based on this comprehensive search, we find that the UV luminosity function decreases by one order of magnitude from z~8 to z~10 at all luminosities over a four magnitude range. This also implies a decrease of the cosmic star-formation rate density by an order of magnitude within 170 Myr from z~8 to z~10. We show that this accelerated evolution compared to lower redshift can entirely be explained by the fast build-up of the dark matter halo mass function at z>8. Consequently, the predicted UV LFs from several models of galaxy formation are in good agreement with this observed trend, even though the measured UV LF lies at the low end of model predictions. In particular, the number of only 9 observed candidate galaxies is lower, by ~50%, than predicted by galaxy evolution models. The difference is generally still consistent within the Poisson and cosmic variance uncertainties. However, essentially all models predict larger numbers than observed. We discuss the implications of these results in light of the upcoming James Webb Space Telescope mission, which is poised to find much larger samples of z~10 galaxies as well as their progenitors at less than 400 Myr after the Big Bang., 13 pages, 6 figures, minor updates to match accepted version

We utilize the newly acquired, ultra-deep WFC3/IR observations over the Hubble Ultra Deep Field (HUDF) to search for star-forming galaxies at z ~ 8-8.5, only 600 million years from recombination, using a Y 105-dropout selection. The new 4.7 arcmin2 WFC3/IR observations reach to ~28.8 AB mag (5σ) in the Y 105 J 125 H 160 bands. These remarkable data reach ~1.5 AB mag deeper than the previous data over the HUDF, and now are an excellent match to the HUDF optical ACS data. For our search criteria, we use a two-color Lyman break selection technique to identify z ~ 8-8.5Y 105-dropouts. We find five likely z ~ 8-8.5 candidates. The sources have H 160-band magnitudes of ~28.3 AB mag and very blue UV-continuum slopes, with a median estimated β of –2.5 (where f λ λβ). This suggests that z ~ 8 galaxies are not only essentially dust free but also may have very young ages or low metallicities. The observed number of Y 105-dropout candidates is smaller than the 20 ± 6 sources expected assuming no evolution from z ~ 6, but is consistent with the five expected extrapolating the Bouwens et al. luminosity function (LF) results to z ~ 8. These results provide evidence that the evolution in the LF seen from z ~ 7 to z ~ 3 continues to z ~ 8. The remarkable improvement in the sensitivity of WFC3/IR has enabled Hubble Space Telescope to cross a threshold, revealing star-forming galaxies at z~ 8-9.

We present optical spectroscopic follow-up of a sample of Distant Red Galaxies (DRGs) with K < 22.5 (Vega), selected by J-K > 2.3, in the Hubble Deep Field South, the MS 1054-03 field, and the Chandra Deep Field South. Spectroscopic redshifts were obtained for 15 DRGs. Only 2 out of 15 DRGs are located at z < 2, suggesting a high efficiency to select high-redshift sources. From other spectroscopic surveys in the CDFS targeting intermediate to high redshift populations selected with different criteria, we find spectroscopic redshifts for a further 30 DRGs. We use the sample of spectroscopically confirmed DRGs to establish the high quality (scatter in ��z/(1+z) of ~ 0.05) of their photometric redshifts in the considered deep fields, as derived with EAZY (Brammer et al. 2008). Combining the spectroscopic and photometric redshifts, we find that 74% of DRGs with K < 22.5 lie at z > 2. The combined spectroscopic and photometric sample is used to analyze the distinct intrinsic and observed properties of DRGs at z < 2 and z > 2. In our photometric sample to K < 22.5, low-redshift DRGs are brighter in K than high-redshift DRGs by 0.7 mag, and more extincted by 1.2 mag in Av. Our analysis shows that the DRG criterion selects galaxies with different properties at different redshifts. Such biases can be largely avoided by selecting galaxies based on their rest-frame properties, which requires very good multi-band photometry and high quality photometric redshifts., Accepted for publication in the Astrophysical Journal, 13 pages, 8 figures, 5 tables

We select 25,000 galaxies from the NEWFIRM Medium Band Survey (NMBS) to study the rest-frame U-V color distribution of galaxies at 0 < z < 2.5. The five unique NIR filters of the NMBS enable the precise measurement of photometric redshifts and rest-frame colors for 9,900 galaxies at 1 < z < 2.5. The rest-frame U-V color distribution at all z, Comment: 5 pages + 4 figures, emulateapj format. Accepted for publication in ApJ Letters

We present new information on galaxies in the vicinity of luminous radio galaxies and quasars at z=4,5,6. These fields were previously found to contain overdensities of Lyman Break Galaxies (LBGs) or spectroscopic Lyman alpha emitters. We use HST and Spitzer data to infer stellar masses, and contrast our results with large samples of LBGs in more average environments as probed by the Great Observatories Origins Deep Survey (GOODS). The following results were obtained. First, LBGs in both overdense regions and in the field at z=4-5 lie on a very similar sequence in a z'-[3.6] versus [3.6] color-magnitude diagram. This is interpreted as a sequence in stellar mass (log[M*/Msun] = 9-11) in which galaxies become increasingly red due to dust and age as their star formation rate (SFR) increases. Second, the two radio galaxies are among the most massive objects (log[M*/Msun]~11) known to exist at z~4-5, and are extremely rare based on the low number density of such objects as estimated from the ~25x larger area GOODS survey. We suggest that the presence of these massive galaxies and supermassive black holes has been boosted through rapid accretion of gas or merging inside overdense regions. Third, the total stellar mass found in the z=4 ``proto-cluster'' TN1338 accounts for 4, based on a comparison with the massive X-ray cluster Cl1252 at z=1.2. Although future near-infrared observations should determine whether any massive galaxies are currently being missed, one possible explanation for this mass difference is that TN1338 evolves into a smaller cluster than Cl1252. This raises the interesting question of whether the most massive protocluster regions at z>4 remain yet to be discovered., The Astrophysical Journal, In Press (17 pages, 7 figures)

Several recent studies have shown that about half of the massive galaxies at z~2 are in a quiescent phase. Moreover, these galaxies are commonly found to be ultra-compact with half-light radii of ~1 kpc. We have obtained a ~29 hr spectrum of a typical quiescent, ultra-dense galaxy at z=2.1865 with the Gemini Near-Infrared Spectrograph. The spectrum exhibits a strong optical break and several absorption features, which have not previously been detected in z>2 quiescent galaxies. Comparison of the spectral energy distribution with stellar population synthesis models implies a low star formation rate (SFR) of 1-3 Msol/yr, an age of 1.3-2.2 Gyr, and a stellar mass of ~2x10^11 Msol. We detect several faint emission lines, with emission-line ratios of [NII]/Halpha, [SII]/Halpha and [OII]/[OIII] typical of low-ionization nuclear emission-line regions. Thus, neither the stellar continuum nor the nebular emission implies active star formation. The current SFR is, Accepted for publication in the Astrophysical Journal

We conduct a search for z>~7 dropout galaxies behind 11 massive lensing clusters using 21 arcmin**2 of deep HST NICMOS, ACS, and WFPC2 image data. In total, over this entire area, we find only one robust z~7 z-dropout candidate (previously reported around Abell 1689). Four less robust z-dropout and J-dropout candidates are also found. The nature of the four weaker candidates could not be precisely determined due to the limited depth of the available optical data, but detailed simulations suggest that all four could be low-redshift interlopers. We compare these numbers with what we might expect using the z~7 UV luminosity function (LF) determined from field searches. We predict 2.7 z~7 z-dropouts and 0.3 z~9 J-dropouts over our cluster search area, in reasonable agreement with our observational results, given the small numbers. The number of z>~7 candidates we find in the present search are much lower than has been reported in several previous studies of the prevalence of z>~7 galaxies behind lensing clusters. To understand these differences, we examined z>~7 candidates in other studies and conclude that only a small fraction are likely to be z>~7 galaxies. Our findings support models that show that gravitational lensing from clusters is of the most value for detecting galaxies at magnitudes brighter than L* (H, Comment: 9 pages, 3 figures, 5 tables, in press

We use all available deep optical ACS and near-IR data over both the HUDF and the two GOODS fields to search for star-forming galaxies at z>~7 and constrain the UV LF within the first 700 Myrs. Our data set includes ~23 arcmin^2 of deep NICMOS J+H data and ~248 arcmin^2 of ground-based (ISAAC+MOIRCS) data, coincident with ACS optical data of greater or equal depths. In total, we find 8 ~7.3 z-dropouts in our search fields, but no z~9 J-dropout candidates. A careful consideration of a wide variety of different contaminants suggest an overall contamination level of just ~12% for our z-dropout selection. After performing detailed simulations to accurately estimate the selection volumes, we derive constraints on the UV LFs at z~7 and z~9. For a faint-end slope alpha=-1.74, our most likely values for M*(UV) and phi* at z~7 are -19.8+/-0.4 mag and 1.1_{-0.7}^{+1.7} x 10^{-3} Mpc^{-3}, respectively. Our search results for z~9 J-dropouts set a 1 sigma lower limit on M*(UV) of -19.6 mag assuming that phi* and alpha are the same as their values at slightly later times. This lower limit on M*(UV) is 1.4 mag fainter than our best-fit value at z~4, suggesting that the UV LF has undergone substantial evolution over this time period. No evolution is ruled out at 99% confidence from z~7 to z~6 and at 80% confidence from z~9 to z~7. The inferred brightening in M*(UV) with redshift (i.e., M*(UV) = (-21.02+/-0.09) + (0.36+/-0.08)(z - 3.8)) matches the evolution expected in the halo mass function, if the mass-to-light ratio of halos evolves as ~(1+z)**{-1}. Finally, we consider the shape of the UV LF at z>~5 and discuss the implications of the Schechter-like form of the observed LFs, particularly the unexpected abrupt cut-off at the bright end.

The existence of massive galaxies with strongly suppressed star formation at z~2.3, identified in a previous paper, suggests that a red sequence may already be in place beyond z=2. In order to test this hypothesis, we study the rest-frame U-B color distribution of massive galaxies at 23 sigma) red sequence, which hosts ~60% of the stellar mass at the high-mass end. The red-sequence galaxies have little or no ongoing star formation, as inferred from both emission-line diagnostics and stellar continuum shapes. Their strong Balmer breaks and their location in the rest-frame U-B, B-V plane indicate that they are in a post-starburst phase, with typical ages of ~0.5-1.0 Gyr. In order to study the evolution of the red sequence, we compare our sample with spectroscopic massive galaxy samples at 0.0210^11 Msol) of the red sequence grow by factors of ~8 and ~6, respectively. We explore simple models to explain the observed evolution. Passive evolution models predict too strong d(U-B), and produce z~0 galaxies that are too red. More complicated models that include aging, galaxy transformations, and red mergers can explain both the number density and color evolution of the massive end of the red sequence between z~2.3 and the present., Comment: Accepted for publication in the Astrophysical Journal

We use the ACS BViz data from the HUDF and all other deep HST ACS fields (including the GOODS fields) to find large samples of star-forming galaxies at z~4 and z~5 and to extend our previous z~6 sample. These samples contain 4671, 1416, and 627 B, V, and i dropouts, respectively, and reach to extremely low luminosities (0.01-0.04 L* or M(UV)~-16 to -17), allowing us to determine the rest-frame UV luminosity function (LF) and faint-end slope alpha at z~4-6 to high accuracy. We find faint-end slopes alpha of -1.73+/-0.05 at z~4, -1.66+/-0.09 at z~5, and -1.74+/-0.16 at z~6 -- suggesting that the faint-end slope is very steep and shows little evolution with cosmic time. We find that M*(UV) brightens considerably in the 0.7 Gyr from z~6 to z~4 (by ~0.7 mag from M*=-20.24+/-0.19 to M*=-20.98+/-0.10). The observed increase in the characteristic luminosity over this range is almost identical to that expected for the halo mass function -- suggesting that the observed evolution is likely due to the hierarchical coalescence and merging of galaxies. The evolution in phi* is not significant. The UV luminosity density at z~6 is modestly lower (0.45+/-0.09 times) than that at z~4 (integrated to -17.5 AB mag) though a larger change is seen in the dust-corrected star-formation rate density. We thoroughly examine published LF results and assess the reasons for their wide dispersion. We argue that the results reported here are the most robust available. The extremely steep faint-end slopes alpha found here suggest that lower luminosity galaxies play a significant role in reionizing the universe. Finally, we consider recent search results for galaxies at z~7-8 and use them to extend our estimates of the evolution in M* from z~7-8 to z~4.

Using the Gemini Near-InfraRed Spectrograph (GNIRS), we have assembled a complete sample of 20 K-selected galaxies at 2.0, Comment: Accepted for publication in the Astrophysical Journal

Thomson optical depth tau measurements from Planck provide new insights into the reionization of the universe. In pursuit of model-independent constraints on the properties of the ionising sources, we determine the empirical evolution of the cosmic ionizing emissivity. We use a simple two-parameter model to map out the evolution in the emissivity at z>~6 from the new Planck optical depth tau measurements, from the constraints provided by quasar absorption spectra and from the prevalence of Ly-alpha emission in z~7-8 galaxies. We find the redshift evolution in the emissivity dot{N}_{ion}(z) required by the observations to be d(log Nion)/dz=-0.15(-0.11)(+0.08), largely independent of the assumed clumping factor C_{HII} and entirely independent of the nature of the ionising sources. The trend in dot{N}_{ion}(z) is well-matched by the evolution of the galaxy UV-luminosity density (dlog_{10} rho_UV/dz=-0.11+/-0.04) to a magnitude limit >~-13 mag, suggesting that galaxies are the sources that drive the reionization of the universe. The role of galaxies is further strengthened by the conversion from the UV luminosity density rho_UV to dot(N)_{ion}(z) being possible for physically-plausible values of the escape fraction f_{esc}, the Lyman-continuum photon production efficiency xi_{ion}, and faint-end cut-off $M_{lim}$ to the luminosity function. Quasars/AGN appear to match neither the redshift evolution nor normalization of the ionizing emissivity. Based on the inferred evolution in the ionizing emissivity, we estimate that the z~10 UV-luminosity density is 8(-4)(+15)x lower than at $z~6, consistent with the observations. The present approach of contrasting the inferred evolution of the ionizing emissivity with that of the galaxy UV luminosity density adds to the growing observational evidence that faint, star-forming galaxies drive the reionization of the universe., 20 pages, 12 figures, 5 tables, Astrophysical Journal, updated to match version in press, Figure 6 shows the main result of the paper

In this paper we study a key phase in the formation of massive galaxies: the transition of star forming galaxies into massive (M_stars~10^11 Msun), compact (r_e~1 kpc) quiescent galaxies, which takes place from z~3 to z~1.5. We use HST grism redshifts and extensive photometry in all five 3D-HST/CANDELS fields, more than doubling the area used previously for such studies, and combine these data with Keck MOSFIRE and NIRSPEC spectroscopy. We first confirm that a population of massive, compact, star forming galaxies exists at z~2, using K-band spectroscopy of 25 of these objects at 2.0, Comment: Accepted for publication in ApJ. Figure 18 summarizes the kinematics part of the paper. Figure 22 demonstrates how star forming galaxies grow in size and mass (in observations and simulations), and Figure 28 is a cartoon version of a model that describes the evolution of star forming and quiescent galaxies since z~3

Detection of Lyman-Break Galaxies (LBGs) at high-redshift can be affected by gravitational lensing induced by foreground deflectors not only in galaxy clusters, but also in blank fields. We quantify the impact of strong magnification in the samples of $B$, $V$, $i$, $z$ $\&$ $Y$ LBGs ($4\lesssim z \lesssim8$) observed in the XDF and GOODS/CANDELS fields, by investigating the proximity of dropouts to foreground objects. We find that $\sim6\%$ of bright LBGs ($m_{H_{160}}2$) by foreground objects. This fraction decreases from $\sim 3.5\%$ at $z\sim6$ to $\sim1.5\%$ at $z\sim4$. Since the observed fraction of strongly lensed galaxies is a function of the shape of the luminosity function (LF), it can be used to derive Schechter parameters, $\alpha$ and $M_{\star}$, independently from galaxy number counts. Our magnification bias analysis yields Schechter-function parameters in close agreement with those determined from galaxy counts albeit with larger uncertainties. Extrapolation of our analysis to $z\gtrsim 8$ suggests that future surveys with JSWT, WFIRST and EUCLID should find excess LBGs at the bright-end, even if there is an intrinsic exponential cutoff of number counts. Finally, we highlight how the magnification bias measurement near the detection limit can be used as probe of the population of galaxies too faint to be detected. Preliminary results using this novel idea suggest that the magnification bias at $M_{UV}\sim -18$ is not as strong as expected if $\alpha\lesssim -1.7$ extends well below the current detection limits in the XDF. At face value this implies a flattening of the LF at $M_{UV}\gtrsim-16.5$. However, selection effects and completeness estimates are difficult to quantify precisely. Thus, we do not rule out a steep LF extending to $M_{UV}\gtrsim -15$., Comment: Submitted to ApJ on 18/12/2014

We search the complete Hubble Frontier Field dataset of Abell 2744 and its parallel field for z~10 sources to further refine the evolution of the cosmic star-formation rate density (SFRD) at z>8. We independently confirm two images of the recently discovered triply-imaged z~9.8 source by Zitrin et al. (2014) and set an upper limit for similar z~10 galaxies with red colors of J_125-H_160>1.2 in the parallel field of Abell 2744. We utilize extensive simulations to derive the effective selection volume of Lyman-break galaxies at z~10, both in the lensed cluster field and in the adjacent parallel field. Particular care is taken to include position-dependent lensing shear to accurately account for the expected sizes and morphologies of highly-magnified sources. We show that both source blending and shear reduce the completeness at a given observed magnitude in the cluster, particularly near the critical curves. These effects have a significant, but largely overlooked, impact on the detectability of high-redshift sources behind clusters, and substantially reduce the expected number of highly-magnified sources. The detections and limits from both pointings result in a SFRD which is higher by 0.4+-0.4 dex than previous estimates at z~10 from blank fields. Nevertheless, the combination of these new results with all other estimates remain consistent with a rapidly declining SFRD in the 170 Myr from z~8 to z~10 as predicted by cosmological simulations and dark-matter halo evolution in LambdaCDM. Once biases introduced by magnification-dependent completeness are accounted for, the full six cluster and parallel Frontier Field program will be an extremely powerful new dataset to probe the evolution of the galaxy population at z>8 before the advent of the JWST., 10 pages, 7 figures, changed to match accepted version to appear in ApJ

We present a spectroscopic redshift measurement of a very bright Lyman break galaxy at z=7.7302+-0.0006 using Keck/MOSFIRE. The source was pre-selected photometrically in the EGS field as a robust z~8 candidate with H=25.0 mag based on optical non-detections and a very red Spitzer/IRAC [3.6]-[4.5] broad-band color driven by high equivalent width [OIII]+Hbeta line emission. The Lyalpha line is reliably detected at 6.1 sigma and shows an asymmetric profile as expected for a galaxy embedded in a relatively neutral inter-galactic medium near the Planck peak of cosmic reionization. The line has a rest-frame equivalent width of EW0=21+-4 A and is extended with V_FWHM=360+90-70 km/s. The source is perhaps the brightest and most massive z~8 Lyman break galaxy in the full CANDELS and BoRG/HIPPIES surveys, having assembled already 10^(9.9+-0.2) M_sol of stars at only 650 Myr after the Big Bang. The spectroscopic redshift measurement sets a new redshift record for galaxies. This enables reliable constraints on the stellar mass, star-formation rate, formation epoch, as well as combined [OIII]+Hbeta line equivalent widths. The redshift confirms that the IRAC [4.5] photometry is very likely dominated by line emission with EW0(OIII+Hbeta)= 720-150+180 A. This detection thus adds to the evidence that extreme rest-frame optical emission lines are a ubiquitous feature of early galaxies promising very efficient spectroscopic follow-up in the future with infrared spectroscopy using JWST and, later, ELTs., Comment: 6 pages, 4 figures, small updates to match ApJL accepted version

The IRAC ultradeep field (IUDF) and IRAC Legacy over GOODS (IGOODS) programs are two ultradeep imaging surveys at 3.6{\mu}m and 4.5{\mu}m with the Spitzer Infrared Array Camera (IRAC). The primary aim is to directly detect the infrared light of reionization epoch galaxies at z > 7 and to constrain their stellar populations. The observations cover the Hubble Ultra Deep Field (HUDF), including the two HUDF parallel fields, and the CANDELS/GOODS-South, and are combined with archival data from all previous deep programs into one ultradeep dataset. The resulting imaging reaches unprecedented coverage in IRAC 3.6{\mu}m and 4.5{\mu}m ranging from > 50 hour over 150 arcmin^2, > 100 hour over 60 sq arcmin2, to 200 hour over 5 - 10 arcmin$^2$. This paper presents the survey description, data reduction, and public release of reduced mosaics on the same astrometric system as the CANDELS/GOODS-South WFC3 data. To facilitate prior-based WFC3+IRAC photometry, we introduce a new method to create high signal-to-noise PSFs from the IRAC data and reconstruct the complex spatial variation due to survey geometry. The PSF maps are included in the release, as are registered maps of subsets of the data to enable reliability and variability studies. Simulations show that the noise in the ultradeep IRAC images decreases approximately as the square root of integration time over the range 20 - 200 hours, well below the classical confusion limit, reaching 1{\sigma} point source sensitivities as faint as of 15 nJy (28.5 AB) at 3.6{\mu}m and 18 nJy (28.3 AB) at 4.5{\mu}m. The value of such ultradeep IRAC data is illustrated by direct detections of z = 7 - 8 galaxies as faint as HAB = 28., Comment: Accepted for publication in ApJS

We have detected 506 i-dropouts (z~6 galaxies) in deep, wide-area HST ACS fields: HUDF, enhanced GOODS, and HUDF-Parallel ACS fields (HUDF-Ps). The contamination levels are =92% are at z~6). With these samples, we present the most comprehensive, quantitative analyses of z~6 galaxies yet and provide optimal measures of the UV luminosity function (LF) and luminosity density at z~6, and their evolution to z~3. We redetermine the size and color evolution from z~6 to z~3. Field-to-field variations (cosmic variance), completeness, flux, and contamination corrections are modelled systematically and quantitatively. After corrections, we derive a rest-frame continuum UV (~1350 A) LF at z~6 that extends to M_{1350,AB} ~ -17.5 (0.04L*(z=3)). There is strong evidence for evolution of the LF between z~6 and z~3, most likely through a brightening (0.6+/-0.2 mag) of M* (at 99.7% confidence) though the degree depends upon the faint-end slope. As expected from hierarchical models, the most luminous galaxies are deficient at z~6. Density evolution (phi*) is ruled out at >99.99% confidence. Despite large changes in the LF, the luminosity density at z~6 is similar (0.82+/-0.21x) to that at z~3. Changes in the mean UV color of galaxies from z~6 to z~3 suggest an evolution in dust content, indicating the true evolution is substantially larger: at z~6 the star formation rate density is just ~30% of the z~3 value. Our UV luminosity function is consistent with z~6 galaxies providing the necessary UV flux to reionize the universe.

We use the NICMOS Treasury and ACS HUDF images to measure the extinction corrected star formation history for 4681 galaxies in the region common to both images utilizing the star formation rate distribution function and other techniques similar to those employed with the NICMOS and WFPC2 images in the HDFN. Unlike the HDFN the NICMOS region of the HUDF appears to lack highly luminous and high star formation rate galaxies at redshifts beyond 3. The HUDF provides a region that is completely uncorrelated to the HDFN and therefore provides and independent measure of the star formation history of the universe. The combined HUDF and HDFN star formation rates show an average star formation rate of 0.2 solar masses per yer per cubic megaparsec. The average SFR of the combined fields at z = 1-3 is 0.29 solar masses per year per cubic megaparsec while the average at z = 4-6 is 1.2 solar masses per year per cubic megaparsec. The SFRs at all redshifts are within 3 sigma of the average over all redshifts.

We measured the angular clustering at z~6 from a large sample of i-dropout galaxies (293 with z 0.5L*) shows the presence of clustering at 94% significance. For this sample we derive a (co-moving) correlation length of r_0=4.5^{+2.1}_{-3.2} h_{72}^{-1} Mpc and bias b=4.1^{+1.5}_{-2.6}, using an accurate model for the redshift distribution. No clustering could be detected in the much deeper but significantly smaller UDF, yielding b 5 star formation might have occurred more efficiently compared to that at z=3-4. This may provide an explanation for the very mild evolution observed in the UV luminosity density between z=6 and z=3. Although our results are consistent with such a scenario, the errors are too large to find conclusive evidence for this.

HST ACS and NICMOS data are now of sufficient depth and areal coverage to place strong constraints on the formation and evolution of galaxies during the first 1–2 Gyrs of the universe. Of particular interest are galaxies at z ∼ 6 since they represent the earliest epoch accessible to current high-efficiency optical instrumentation. To this end, we have been involved in the systematic construction of a large sample of 346 z ∼ 6 objects from all the deepest wide-area HST data (UDF, UDF-Parallel, and GOODS fields). They have been used to construct an optimal determination of the rest-frame continuum UV LF at z ∼ 6. Our LF extends to over 3 magnitudes below L*, fainter than has been done at z ∼ 3. Over the interval z ∼ 6 to z ∼ 3, we find strong evidence for evolution in the UV LF. Though we can thus far make no strong claims on an evolution of the faint end slope, the characteristic luminosity appears to have approximately doubled over the interval z ∼ 6 to z ∼ 3, consistent with hierarchical expectations. Remarkably, this shift to lower luminosities extends to even higher redshifts. Using all deep J + H NICMOS observations (800 orbits in total), we have been able to demonstrate that the bright end of the LF ( > 0.3 L z = 3 ∗ ) is at least 5 times lower at z ∼ 10 than at z ∼ 4, with a similar deficit being established from our recent detections and first statistical sample of z ∼ 7–8 galaxies using our UDF NICMOS data. In these precedings, we discuss what is known about the UV LF and UV luminosity density at z ∼ 6–10 from current data and its evolution relative to z ∼ 3. We also describe several exciting prospects for advance in this area over the next year.

Hubble Space Telescope Advanced Camera for Surveys images of the young binary GG Tauri and its circumbinary disk in V and I bandpasses were obtained in 2002 and are the most detailed of this system to date. They confirm features previously seen in the disk including a "gap" apparently caused by shadowing from circumstellar material, an asymmetrical distribution of light about the line of sight on the near edge of the disk, enhanced brightness along the near edge of the disk due to forward scattering, and a compact reflection nebula near the secondary star. New features are seen in the ACS images: two short filaments along the disk, localized but strong variations in the disk intensity ("gaplets"), and a "spur" or filament extending from the reflection nebulosity near the secondary. The back side of the disk is detected in the V band for the first time. The disk appears redder than the combined light from the stars, which may be explained by a varied distribution of grain sizes. The brightness asymmetries along the disk suggest that it is asymmetrically illuminated by the stars due to extinction by nonuniform circumstellar material or that the illuminated surface of the disk is warped by tidal effects (or perhaps both). Localized, time-dependent brightness variations in the disk are also seen.

This paper describes the observations and data reduction techniques for the version 2.0 images and catalog of the NICMOS Ultra Deep Field Treasury program. All sources discussed in this paper are based on detections in the combined NICMOS F110W and F160W bands only. The NICMOS images are drizzled to 0.09 arc second pixels and aligned to the ACS UDF F850LP image which was rebinned to the same pixel scale. These form the NICMOS version 2.0 UDF images. The catalog sources are chosen with a conservative detection limit to avoid the inclusion of numerous spurious sources. The catalog contains 1293 objects in the 144 x 144 arc sececonds NICMOS subfield of the UDF. The 5 sigma signal to noise level is an average 0.6 arc second diameter aperture AB magnitude of ~27.7 at 1.1 and 1.6 microns. The catalog sources, listed in order of right ascension, satisfy a minimum signal to noise criterion of 1.4 sigma in at least 7 contiguous pixels of the combined F110W and F160W image