Showing 51–100 of 175 results for author: Bouwens, R J

Unveiling the dynamical state of massive clusters through the ICL fraction

Authors: Y. Jiménez-Teja, R. Dupke, N. Benítez, A. M. Koekemoer, A. Zitrin, K. Umetsu, B. L. Ziegler, B. L. Frye, H. Ford, R. J. Bouwens, L. D. Bradley, T. Broadhurst, D. Coe, M. Donahue, G. J. Graves, C. Grillo, L. Infante, S. Jouvel, D. D. Kelson, O. Lahav, R. Lazkoz, D. Lemze, D. Maoz, E. Medezinski, P. Melchior , et al. (12 additional authors not shown)

Abstract: We have selected a sample of eleven massive clusters of galaxies observed by the Hubble Space Telescope in order to study the impact of the dynamical state on the IntraCluster Light (ICL) fraction, the ratio of total integrated ICL to the total galaxy member light. With the exception of the Bullet cluster, the sample is drawn from the Cluster Lensing and Supernova Survey and the Frontier Fields pr… ▽ More We have selected a sample of eleven massive clusters of galaxies observed by the Hubble Space Telescope in order to study the impact of the dynamical state on the IntraCluster Light (ICL) fraction, the ratio of total integrated ICL to the total galaxy member light. With the exception of the Bullet cluster, the sample is drawn from the Cluster Lensing and Supernova Survey and the Frontier Fields program, containing five relaxed and six merging clusters. The ICL fraction is calculated in three optical filters using the CHEFs IntraCluster Light Estimator, a robust and accurate algorithm free of a priori assumptions. We find that the ICL fraction in the three bands is, on average, higher for the merging clusters, ranging between $\sim7-23\%$, compared with the $\sim 2-11\%$ found for the relaxed systems. We observe a nearly constant value (within the error bars) in the ICL fraction of the regular clusters at the three wavelengths considered, which would indicate that the colors of the ICL and the cluster galaxies are, on average, coincident and, thus, their stellar populations. However, we find a higher ICL fraction in the F606W filter for the merging clusters, consistent with an excess of lower-metallicity/younger stars in the ICL, which could have migrated violently from the outskirts of the infalling galaxies during the merger event. △ Less

Submitted 13 March, 2018; originally announced March 2018.

Comments: 25 pages, 13 figures, 2 tables. Accepted for publication in ApJ

The ALMA Frontier Fields Survey - IV. Lensing-corrected 1.1 mm number counts in Abell 2744, MACSJ0416.1-2403 and MACSJ1149.5+2223

Authors: A. M. Muñoz Arancibia, J. González-López, E. Ibar, F. E. Bauer, M. Carrasco, N. Laporte, T. Anguita, M. Aravena, F. Barrientos, R. J. Bouwens, R. Demarco, L. Infante, R. Kneissl, N. Nagar, N. Padilla, C. Romero-Cañizales, P. Troncoso, A. Zitrin

Abstract: [abridged] Characterizing the number counts of faint, dusty star-forming galaxies is currently a challenge even for deep, high-resolution observations in the FIR-to-mm regime. They are predicted to account for approximately half of the total extragalactic background light at those wavelengths. Searching for dusty star-forming galaxies behind massive galaxy clusters benefits from strong lensing, en… ▽ More [abridged] Characterizing the number counts of faint, dusty star-forming galaxies is currently a challenge even for deep, high-resolution observations in the FIR-to-mm regime. They are predicted to account for approximately half of the total extragalactic background light at those wavelengths. Searching for dusty star-forming galaxies behind massive galaxy clusters benefits from strong lensing, enhancing their measured emission while increasing spatial resolution. Derived number counts depend, however, on mass reconstruction models that properly constrain these clusters. We estimate the 1.1 mm number counts along the line of sight of three galaxy clusters, i.e. Abell 2744, MACSJ0416.1-2403 and MACSJ1149.5+2223, which are part of the ALMA Frontier Fields Survey. We perform detailed simulations to correct these counts for lensing effects. We use several publicly available lensing models for the galaxy clusters to derive the intrinsic flux densities of our sources. We perform Monte Carlo simulations of the number counts for a detailed treatment of the uncertainties in the magnifications and adopted source redshifts. We find an overall agreement among the number counts derived for the different lens models, despite their systematic variations regarding source magnifications and effective areas. Our number counts span ~2.5 dex in demagnified flux density, from several mJy down to tens of uJy. Our number counts are consistent with recent estimates from deep ALMA observations at a 3$σ$ level. Below $\approx$ 0.1 mJy, however, our cumulative counts are lower by $\approx$ 1 dex, suggesting a flattening in the number counts. In our deepest ALMA mosaic, we estimate number counts for intrinsic flux densities $\approx$ 4 times fainter than the rms level. This highlights the potential of probing the sub-10 uJy population in larger samples of galaxy cluster fields with deeper ALMA observations. △ Less

Submitted 14 September, 2018; v1 submitted 11 December, 2017; originally announced December 2017.

Comments: 19 pages, 14 figures, 3 tables. Accepted for publication in A&A

Journal ref: A&A 620, A125 (2018)

Very low-luminosity galaxies in the early universe have observed sizes similar to single star cluster complexes

Authors: R. J. Bouwens, G. D. Illingworth, P. A. Oesch, M. Maseda, B. Ribeiro, M. Stefanon, D. Lam

Abstract: We compare the sizes and luminosities of 307 faint z=6-8 sources revealed by the Hubble Frontier Fields (HFF) program with sources in the nearby universe. Making use of the latest lensing models and data from the first four HFF clusters with an extensive suite of public lens models, we measure both the sizes and luminosities for 153 z~6, 101 z~7, and 53 z~8 galaxies. The sizes range over more than… ▽ More We compare the sizes and luminosities of 307 faint z=6-8 sources revealed by the Hubble Frontier Fields (HFF) program with sources in the nearby universe. Making use of the latest lensing models and data from the first four HFF clusters with an extensive suite of public lens models, we measure both the sizes and luminosities for 153 z~6, 101 z~7, and 53 z~8 galaxies. The sizes range over more than a decade from ~500 to <50 pc. Extremely small sizes are inferred for many of our lowest luminosity sources, reaching individual sizes as small as 10-30 pc (the smallest is 11(-6)(+28) pc). The uncertainty in these measures ranges from 80 pc for the largest sources to typically about 20 pc for the smallest. Such sizes are smaller than extrapolations of the size-luminosity relation, and expectations for the completeness of our faint samples, suggesting a likely break in the size-luminosity relation at ~-17 mag with size proportional to L**(0.50(-0.11)(+0.10)). The sizes and luminosities of the lowest-luminosity sources are similar to those of single star cluster complexes like 30 Doradus in the lower-redshift universe and -- in a few cases -- super star clusters. Remarkably, our identification of these compact, faint star-forming sources in the z~6-8 universe also allow us to set upper limits on the proto-globular cluster LF at z~6. Comparisons with recent models allow us to rule out (with some caveats) some scenarios for proto-globular cluster formation and set useful upper limits on other less extreme ones. Our results suggest we may be very close to discovering a bona-fide population of forming globular clusters at high redshift. △ Less

Submitted 28 August, 2018; v1 submitted 6 November, 2017; originally announced November 2017.

Comments: The authors have decided to withdraw this manuscript from the Astrophysical Journal given the unprofessional footnotes in the original astro-ph posting, offering our sincere apologies

The Dearth of z~10 Galaxies in all HST Legacy Fields -- The Rapid Evolution of the Galaxy Population in the First 500 Myr

Authors: P. A. Oesch, R. J. Bouwens, G. D. Illingworth, I. Labbe, M. Stefanon

Abstract: 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… ▽ More 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. △ Less

Submitted 27 February, 2018; v1 submitted 30 October, 2017; originally announced October 2017.

Comments: 13 pages, 6 figures, minor updates to match accepted version

The MUSE Hubble Ultra Deep Field Survey: IV. Global properties of C III] emitters

Authors: Michael V. Maseda, Jarle Brinchmann, Marijn Franx, Roland Bacon, Rychard J. Bouwens, Kasper B. Schmidt, Leindert A. Boogaard, Thierry Contini, Anna Feltre, Hanae Inami, Wolfram Kollatschny, Raffaella A. Marino, Johan Richard, Anne Verhamme, Lutz Wisotzki

Abstract: The C III] 1907,1909 emission doublet has been proposed as an alternative to Lyman-alpha in redshift confirmations of galaxies at z > 6 since it is not attenuated by the largely neutral intergalactic medium at these redshifts and is believed to be strong in the young, vigorously star-forming galaxies present at these early cosmic times. We present a statistical sample of 17 C III]-emitting galaxie… ▽ More The C III] 1907,1909 emission doublet has been proposed as an alternative to Lyman-alpha in redshift confirmations of galaxies at z > 6 since it is not attenuated by the largely neutral intergalactic medium at these redshifts and is believed to be strong in the young, vigorously star-forming galaxies present at these early cosmic times. We present a statistical sample of 17 C III]-emitting galaxies beyond z~1.5 using 30 hour deep VLT/MUSE integral field spectroscopy covering 2 square arcminutes in the Hubble Deep Field South (HDFS) and Ultra Deep Field (UDF), achieving C III] sensitivities of ~2e-17 erg/s/cm^2 in the HDFS and ~7e-18 erg/s/cm^2 in the UDF. The rest-frame equivalent widths range from 2 to 19 Angstroms. These 17 galaxies represent ~3% of the total sample of galaxies found between 1.5 < z < 4. They also show elevated star formation rates, lower dust attenuation, and younger mass-weighted ages than the general population of galaxies at the same redshifts. Combined with deep slitless grism spectroscopy from the HST/WFC3 in the UDF, we can tie the rest-frame ultraviolet C III] emission to rest-frame optical emission lines, namely [O III] 5007, finding a strong correlation between the two. Down to the flux limits that we observe (~1e-18 erg/s/cm^2 with the grism data in the UDF), all objects with a rest-frame [O III] 4959,5007 equivalent width in excess of 250 Angstroms, the so-called Extreme Emission Line Galaxies, have detections of C III] in our MUSE data. More detailed studies of the C III]-emitting population at these intermediate redshifts will be crucial to understand the physical conditions in galaxies at early cosmic times and to determine the utility of C III] as a redshift tracer. △ Less

Submitted 17 October, 2017; originally announced October 2017.

Comments: 16 pages and 13 figures, accepted for publication in A&A

Dust masses of $z>5$ galaxies from SED fitting and ALMA upper limits

Authors: Hiroyuki Hirashita, Denis Burgarella, Rychard J. Bouwens

Abstract: We aim at constraining the dust mass in high-redshift ($z\gtrsim 5$) galaxies using the upper limits obtained by ALMA in combination with the rest-frame UV--optical spectral energy distributions (SEDs). For SED fitting, because of degeneracy between dust extinction and stellar age, we focus on two extremes: continuous star formation (Model A) and instantaneous star formation (Model B). We apply th… ▽ More We aim at constraining the dust mass in high-redshift ($z\gtrsim 5$) galaxies using the upper limits obtained by ALMA in combination with the rest-frame UV--optical spectral energy distributions (SEDs). For SED fitting, because of degeneracy between dust extinction and stellar age, we focus on two extremes: continuous star formation (Model A) and instantaneous star formation (Model B). We apply these models to Himiko (as a representative UV-bright object) and a composite SED of $z>5$ Lyman break galaxies (LBGs). For Himiko, Model A requires a significant dust extinction, which leads to a high dust temperature $>70$ K for consistency with the ALMA upper limit. This high dust temperature puts a strong upper limit on the total dust mass $M_\mathrm{d}\lesssim 2\times 10^6$ M$_{\odot}$, and the dust mass produced per supernova (SN) $m_\mathrm{d,SN}\lesssim 0.1$ M$_{\odot}$. Such a low $m_\mathrm{d,SN}$ suggests significant loss of dust by reverse shock destruction or outflow, and implies that SNe are not the dominant source of dust at high $z$. Model B allows $M_\mathrm{d}\sim 2\times 10^7$ M$_{\odot}$ and $m_\mathrm{d,SN}\sim 0.3$ M$_{\odot}$.} We could distinguish between Models A and B if we observe Himiko at {wavelength $<$ 1.2 mm by ALMA. For the LBG sample, we obtain $M_\mathrm{d}\lesssim 2\times 10^6$ M$_{\odot}$ for a typical LBG at $z>5$, but this only puts an upper limit for $m_\mathrm{d,SN}$ as $\sim 2$ M$_{\odot}$. This clarifies the importance of observing UV-bright objects (like Himiko) to constrain the dust production by SNe. △ Less

Submitted 7 September, 2017; originally announced September 2017.

Comments: 12 pages, 6 figures, accepted for publication in MNRAS

Rotation in [CII]-emitting gas in two galaxies at a redshift of 6.8

Authors: Renske Smit, Rychard J. Bouwens, Stefano Carniani, Pascal A. Oesch, Ivo Labbé, Garth D. Illingworth, Paul van der Werf, Larry D. Bradley, Valentino Gonzalez, Jacqueline A. Hodge, Benne W. Holwerda, Roberto Maiolino, Wei Zheng

Abstract: 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 pre… ▽ More 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]λ157.74μ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α 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α 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. △ Less

Submitted 26 January, 2018; v1 submitted 14 June, 2017; originally announced June 2017.

Comments: 14 pages, 5 figures, 1 table, published in Nature (submitted version)

HST imaging of the brightest z~8-9 galaxies from UltraVISTA: the extreme bright end of the UV luminosity function

Authors: Mauro Stefanon, Ivo Labbé, Rychard J. Bouwens, Gabriel B. Brammer, Pascal Oesch, Marijn Franx, Johan P. U. Fynbo, Bo Milvang-Jensen, Adam Muzzin, Garth D. Illingworth, Olivier Le Fèvre, Karina I. Caputi, Benne W. Holwerda, Henry J. McCracken, Renske Smit, Dan Magee

Abstract: We report on the discovery of three especially bright candidate $z_{phot} \gtrsim 8$ galaxies. Five sources were targeted for follow-up with HST/WFC3, selected from a larger sample of 16 bright ($24.8 \lesssim H\lesssim25.5$~mag) candidate $z\gtrsim 8$ LBGs identified over the 1.6 degrees$^2$ of the COSMOS/UltraVISTA field. These were identified as Y and J dropouts by leveraging the deep (Y-to-… ▽ More We report on the discovery of three especially bright candidate $z_{phot} \gtrsim 8$ galaxies. Five sources were targeted for follow-up with HST/WFC3, selected from a larger sample of 16 bright ($24.8 \lesssim H\lesssim25.5$~mag) candidate $z\gtrsim 8$ LBGs identified over the 1.6 degrees$^2$ of the COSMOS/UltraVISTA field. These were identified as Y and J dropouts by leveraging the deep (Y-to-$K_{S} \sim 25.3-24.8$~mag, $5σ$) NIR data from the UltraVISTA DR3 release, deep ground based optical imaging from the CFHTLS and Subaru Suprime Cam programs and Spitzer/IRAC mosaics combining observations from the SMUVS and SPLASH programs. Through the refined spectral energy distributions, which now also include new HyperSuprime Cam g, r, i, z and Y band data, we confirm that 3/5 galaxies have robust $z_{phot}\sim8.0-8.7$, consistent with the initial selection. The remaining 2/5 galaxies have a nominal $z_{phot}\sim2$. However, if we use the HST data alone, these objects have increased probability of being at $z\sim9$. Furthermore, we measure mean UV continuum slopes $β=-1.91\pm0.26$ for the three $z\sim8-9$ galaxies, marginally bluer than similarly luminous $z\sim4-6$ in CANDELS but consistent with previous measurements of similarly luminous galaxies at $z\sim7$. The circularized effective radius for our brightest source is $0.9\pm0.2$ kpc, similar to previous measurements for a bright $z\sim11$ galaxy and bright $z\sim7$ galaxies. Finally, enlarging our sample to include the six brightest $z\sim8$ LBGs identified over UltraVISTA (i.e., including three other sources from Labbe et al. 2017, in prep.) we estimate for the first time the volume density of galaxies at the extreme bright ($M_{UV}\sim-22$~mag) end of the $z\sim8$ UV LF. Despite this exceptional result, the still large statistical uncertainties do not allow us to discriminate between a Schechter and a double power-law form. △ Less

Submitted 16 October, 2017; v1 submitted 14 June, 2017; originally announced June 2017.

Comments: Resubmitting after addressing the referee's comments. Appendix B includes the additional 3 sources from Labbe+2017, used for the LF estimate

The HDUV Survey: A Revised Assessment of the Relationship between UV Slope and Dust Attenuation for High-Redshift Galaxies

Authors: Naveen A. Reddy, Pascal A. Oesch, Rychard J. Bouwens, Mireia Montes, Garth D. Illingworth, Charles C. Steidel, Pieter G. van Dokkum, Hakim Atek, Marcella C. Carollo, Anna Cibinel, Brad Holden, Ivo Labbe, Dan Magee, Laura Morselli, Erica J. Nelson, Steve Wilkins

Abstract: 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… ▽ More 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] △ Less

Submitted 25 May, 2017; originally announced May 2017.

Comments: 22 pages, 3 tables, 11 figures, submitted to the Astrophysical Journal

The rest-frame optical (900nm) galaxy luminosity function at z~4-7: abundance matching points to limited evolution in the Mstar/Mhalo ratio at z>4

Authors: Mauro Stefanon, Rychard J. Bouwens, Ivo Labbé, Adam Muzzin, Danilo Marchesini, Pascal Oesch, Valentino Gonzalez

Abstract: We present the first determination of the galaxy luminosity function (LF) at z~4, 5, 6 and 7 in the rest-frame optical at lambda_rest~900 nm (z' band). The rest-frame optical light traces the content in low-mass evolved stars (~stellar mass - Mstar), minimizing potential measurement biases for Mstar. Moreover it is less affected by nebular line emission contamination and dust attenuation, is indep… ▽ More We present the first determination of the galaxy luminosity function (LF) at z~4, 5, 6 and 7 in the rest-frame optical at lambda_rest~900 nm (z' band). The rest-frame optical light traces the content in low-mass evolved stars (~stellar mass - Mstar), minimizing potential measurement biases for Mstar. Moreover it is less affected by nebular line emission contamination and dust attenuation, is independent of stellar population models, and can be probed up to z~8 through Spitzer/IRAC. Our analysis leverages the unique full depth Spitzer/IRAC 3.6um-to-8.0um data over the CANDELS/GOODS-N, CANDELS/GOODS-S and COSMOS/UltraVISTA fields. We find that at absolute magnitudes M_z' fainter than >-23 mag, M_z' linearly correlates with M_UV,1600. At brighter M_z', M_UV,1600 presents a turnover, suggesting that the stellar mass-to-light ratio Mstar/L_UV,1600 could be characterised by a very broad range of values at high stellar masses. Median-stacking analysis recovers a Mstar/L_z' roughly independent on M_z' for M_z'>-23 mag, but exponentially increasing at brighter magnitudes. We find that the evolution of the LF marginally prefers a pure evolution in luminosity over a pure evolution in density, with the characteristic luminosity decreasing by a factor ~5x between z~4 and z~7. Direct application of the recovered Mstar/L_z' generates stellar mass functions consistent with average measurements from the literature. Measurements of the stellar-to-halo mass ratio at fixed cumulative number density show that it is roughly constant with redshift for Mh>10^12Msun. This is also supported by the fact that the evolution of the LF at 4<z<7 can be accounted for by a rigid displacement in luminosity corresponding to the evolution of the halo mass from abundance matching. △ Less

Submitted 15 June, 2017; v1 submitted 28 November, 2016; originally announced November 2016.

Comments: New sect. 2.3 discussing selection biases; new panel in fig.9 with Mstar vs. M_z' . Updated to reflect the accepted version

The HDUV Survey: Six Lyman Continuum Emitter Candidates at z~2 Revealed by HST UV Imaging

Authors: R. P. Naidu, P. A. Oesch, N. Reddy, B. Holden, C. C. Steidel, M. Montes, H. Atek, R. J. Bouwens, C. M. Carollo, A. Cibinel, G. D. Illingworth, I. Labbe, D. Magee, L. Morselli, E. J. Nelson, P. G. van Dokkum, S. Wilkins

Abstract: 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… ▽ More 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. △ Less

Submitted 21 November, 2016; originally announced November 2016.

Comments: 11 pages, 6 figures, submitted to ApJ

The z~6 Luminosity Function Fainter than -15 mag from the Hubble Frontier Fields: The Impact of Magnification Uncertainties

Authors: R. J. Bouwens, P. A. Oesch, G. D. Illingworth, R. S. Ellis, M. Stefanon

Abstract: 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… ▽ More 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. △ Less

Submitted 6 June, 2017; v1 submitted 2 October, 2016; originally announced October 2016.

Comments: 33 pages, 19 figures, 5 tables, accepted for publication in the Astrophysical Journal, updated to match the version in press

Extremely Small Sizes for Faint z~2-8 Galaxies in the Hubble Frontier Fields: A Key Input For Establishing their Volume Density and UV Emissivity

Authors: R. J. Bouwens, G. D. Illingworth, P. A. Oesch, H. Atek, D. Lam, M. Stefanon

Abstract: 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 k… ▽ More 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 <~30 mas (<160-240 pc). As a 2nd size probe, we rotate and stack 26 faint high-magnification sources along the major shear axis. Less elongation is found than even for objects with an intrinsic half-light radius of 10 mas. Together these results indicate that extremely faint z~2-8 galaxies have near point-source profiles in the HFF dataset (half-light radii conservatively <30 mas and likely 5-10 mas). These results suggest smaller completeness corrections and hence much lower volume densities for faint z~2-8 galaxies and shallower faint-end slopes than have been derived in many recent studies (by factors of ~2-3 and by dalpha>~0.1-0.3). △ Less

Submitted 6 June, 2017; v1 submitted 2 August, 2016; originally announced August 2016.

Comments: 19 pages, 15 figures, 3 tables, accepted for publication in ApJ

The ALMA Frontier Fields survey I: 1.1 mm continuum detections in Abell 2744, MACSJ0416.1-2403 and MACSJ1149.5+2223

Authors: J. González-López, F. E. Bauer, C. Romero-Cañizales, R. Kneissl, E. Villard, R. Carvajal, S. Kim, N. Laporte, T. Anguita, M. Aravena, R. J. Bouwens, L. Bradley, M. Carrasco, R. Demarco, H. Ford, E. Ibar, L. Infante, H. Messias, A. M. Muñoz Arancibia, N. Nagar, N. Padilla, E. Treister, P. Troncoso, A. Zitrin

Abstract: Dusty star-forming galaxies are among the most prodigious systems at high redshift (z>1), characterized by high star formation rates and huge dust reservoirs. The bright end of this population has been well characterized in recent years, but considerable uncertainties remain for fainter dusty star-forming galaxies, which are responsible for the bulk of star formation at high redshift and thus play… ▽ More Dusty star-forming galaxies are among the most prodigious systems at high redshift (z>1), characterized by high star formation rates and huge dust reservoirs. The bright end of this population has been well characterized in recent years, but considerable uncertainties remain for fainter dusty star-forming galaxies, which are responsible for the bulk of star formation at high redshift and thus play a key role in galaxy growth and evolution. In this first paper of our series, we describe our methods for finding high redshift faint dusty galaxies using millimeter observations with ALMA. We obtained ALMA 1.1 mm mosaic images for three strong-lensing galaxy clusters from the Frontier Fields survey. The 2'x2' mosaics overlap with the deep HST WFC3/IR footprints and encompass the high magnification regions of each cluster. The combination of extremely high ALMA sensitivity and the magnification power of these clusters allows us to systematically probe the sub-mJy population of dusty star-forming galaxies over a large surveyed area. We present a description of the reduction and analysis of the ALMA continuum observations for the galaxy clusters Abell 2744 (z=0.308), MACSJ0416.1-2403 (z=0.396) and MACSJ1149.5+2223 (z=0.543), for which we reach observed rms sensitivities of 55, 59 and 71 $μ$Jy beam$^{-1}$ respectively. We detect 12 dusty star-forming galaxies at S/N>5.0 across the three clusters, all of them presenting coincidence with NIR detected counterparts in the HST images. None of the sources fall close to the lensing caustics, thus they are not strongly lensed. The observed 1.1 mm flux densities for the total sample of galaxies range from 0.41 to 2.82 mJy, with observed effective radii spanning <0.05" to 0.37"$\pm$0.21". The lensing-corrected sizes of the detected sources appear to be in the same range as those measured in brighter samples, albeit with possibly larger dispersion. △ Less

Submitted 24 October, 2016; v1 submitted 13 July, 2016; originally announced July 2016.

Comments: Accepted for publication in A&A. New version addressing last referee comments, 17 pages and 12 figures

Journal ref: A&A 597, A41 (2017)

All NIRspec needs is HST/WFC3 pre-imaging? The use of Milky Way Stars in WFC3 Imaging to Register NIRspec MSA Observations

Authors: B. W. Holwerda, R. J. Bouwens, M. Trenti, M. A. Kenworthy

Abstract: The James Webb Space Telescope (JWST) will be an exquisite new near-infrared observatory with imaging and multi-object spectroscopy through ESA's NIRspec instrument with its unique Micro-Shutter Array (MSA), allowing for slits to be positioned on astronomical targets by opening specific 0.002"-wide micro shutter doors. To ensure proper target acquisition, the on-sky position of the MSA needs to… ▽ More The James Webb Space Telescope (JWST) will be an exquisite new near-infrared observatory with imaging and multi-object spectroscopy through ESA's NIRspec instrument with its unique Micro-Shutter Array (MSA), allowing for slits to be positioned on astronomical targets by opening specific 0.002"-wide micro shutter doors. To ensure proper target acquisition, the on-sky position of the MSA needs to be verified before spectroscopic observations start. An onboard centroiding program registers the position of pre-identified guide stars in a Target Acquisition (TA) image, a short pre-spectroscopy exposure without dispersion (image mode) through the MSA with all shutters open. The outstanding issue is the availability of Galactic stars in the right luminosity range for TA relative to typical high redshift targets. We explore this here using the stars and $z\sim8$ candidate galaxies identified in the source extractor catalogs of Brightest of Reionizing Galaxies survey (BoRG[z8]), a pure-parallel program with Hubble Space Telescope Wide-Field Camera 3. We find that (a) a single WFC3 field contains enough Galactic stars to satisfy the NIRspec astrometry requirement (20 milli-arcseconds), provided its and the NIRspec TA's are $m_{lim}>24.5$ AB in WFC3 F125W, (b) a single WFC3 image can therefore serve as the pre-image if need be, (c) a WFC3 mosaic and accompanying TA image satisfy the astrometry requirement at $\sim23$ AB mag in WFC3 F125W, (d) no specific Galactic latitude requires deeper TA imaging due to a lack of Galactic stars, and (e) a depth of $\sim24$ AB mag in WFC3 F125W is needed if a guide star in the same MSA quadrant as a target is required. We take the example of a BoRG identified $z\sim8$ candidate galaxy and require a Galactic star within 20" of it. In this case, a depth of 25.5 AB in F125W is required (with $\sim$97% confidence). △ Less

Submitted 28 June, 2016; originally announced June 2016.

Comments: 17 pages, 15 figures, to appear in the Journal of Astronomical Instrumentation

Galaxy candidates at z ~ 10 in archival data from the Brightest of Reionizing Galaxies (BoRG[z8]) survey

Authors: S. R. Bernard, D. Carrasco, M. Trenti, P. A. Oesch, J. F. Wu, L. D. Bradley, K. B. Schmidt, R. J. Bouwens, V. Calvi, C. A. Mason, M. Stiavelli, T. Treu

Abstract: The Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) enabled the search for the first galaxies observed at z ~ 8 - 11 (500 - 700 Myr after the Big Bang). To continue quantifying the number density of the most luminous galaxies (M_AB ~ -22.0) at the earliest epoch observable with HST, we search for z ~ 10 galaxies (F125W-dropouts) in archival data from the Brightest of Reionizing Gala… ▽ More The Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) enabled the search for the first galaxies observed at z ~ 8 - 11 (500 - 700 Myr after the Big Bang). To continue quantifying the number density of the most luminous galaxies (M_AB ~ -22.0) at the earliest epoch observable with HST, we search for z ~ 10 galaxies (F125W-dropouts) in archival data from the Brightest of Reionizing Galaxies (BoRG[z8]) survey, originally designed for detection of z ~ 8 galaxies (F098M-dropouts). By focusing on the deepest 293 arcmin^2 of the data along 62 independent lines of sight, we identify six z ~ 10 candidates satisfying the color selection criteria, detected at S/N > 8 in F160W with M_AB = -22.8 to -21.1 if at z = 10. Three of the six sources, including the two brightest, are in a single WFC3 pointing (~ 4 arcmin^2), suggestive of significant clustering, which is expected from bright galaxies at z ~ 10. However, the two brightest galaxies are too extended to be likely at z ~ 10, and one additional source is unresolved and possibly a brown dwarf. The remaining three candidates have m_AB ~ 26, and given the area and completeness of our search, our best estimate is a number density of sources that is marginally higher but consistent at 2σ with searches in legacy fields. Our study highlights that z ~ 10 searches can yield a small number of candidates, making tailored follow-ups of HST pure-parallel observations viable and effective. △ Less

Submitted 19 June, 2016; originally announced June 2016.

Comments: Accepted by ApJ. 10 pages, 5 figures

The Size and Shape of the Milky Way Disk and Halo from M-type Brown Dwarfs in the BoRG Survey

Authors: Isabel van Vledder, Dieuwertje van der Vlugt, B. W. Holwerda, M. A. Kenworthy, R. J. Bouwens, M. Trenti

Abstract: We have identified 274 M-type Brown Dwarfs in the Hubble Space Telescope's Wide Field Camera 3 (WFC3) pure parallel fields from the Brightest of Reionizing Galaxies (BoRG) survey for high redshift galaxies. These are near-infrared observations with multiple lines-of-sight out of our Milky Way. Using these observed M-type Brown Dwarfs we fitted a Galactic disk and halo model with a Markov chain Mon… ▽ More We have identified 274 M-type Brown Dwarfs in the Hubble Space Telescope's Wide Field Camera 3 (WFC3) pure parallel fields from the Brightest of Reionizing Galaxies (BoRG) survey for high redshift galaxies. These are near-infrared observations with multiple lines-of-sight out of our Milky Way. Using these observed M-type Brown Dwarfs we fitted a Galactic disk and halo model with a Markov chain Monte Carlo (MCMC) analysis. This model worked best with the scale length of the disk fixed at $h$ = 2.6 kpc. For the scale height of the disk, we found $z_0 = 0.29^{+0.02}_{-0.019}$ kpc and for the central number density $ρ_0 = 0.29^{+0.20}_{-0.13}$ \#/pc$^3$. For the halo we derived a flattening parameter $κ$ = 0.45$\pm{0.04}$ and a power-law index $p$ = 2.4$\pm{0.07}$. We found the fraction of M-type brown dwarfs in the local density that belong to the halo to be $f_{h}$ = 0.0075$^{+0.0025}_{-0.0019}$. We found no correlation between subtype of M-dwarf and any model parameters. The total number of M-type Brown Dwarfs in the disk and halo was determined to be $58.2^{+9.81}_{-6.70} \times10^{9}$. We found an upper limit for the fraction of M-type Brown Dwarfs in the halo of 7$^{+5}_{-4}$\%. The upper limit for the total Galactic Disk mass in M-dwarfs is $4.34^{+0.73}_{-0.5}\times10^{9}$ $M_{\odot}$, assuming all M-type Brown Dwarfs have a mass of $80 M_J$. △ Less

Submitted 16 March, 2016; originally announced March 2016.

Comments: 9 Figures, 8 Tables, MNRAS

A Remarkably Luminous Galaxy at z=11.1 Measured with Hubble Space Telescope Grism Spectroscopy

Authors: P. A. Oesch, G. Brammer, P. G. van Dokkum, G. D. Illingworth, R. J. Bouwens, I. Labbe, M. Franx, I. Momcheva, M. L. N. Ashby, G. G. Fazio, V. Gonzalez, B. Holden, D. Magee, R. E. Skelton, R. Smit, L. R. Spitler, M. Trenti, S. P. Willner

Abstract: We present Hubble WFC3/IR slitless grism spectra of a remarkably bright $z\gtrsim10$ galaxy candidate, GN-z11, identified initially from CANDELS/GOODS-N imaging data. A significant spectroscopic continuum break is detected at $λ=1.47\pm0.01~μ$m. The new grism data, combined with the photometric data, rule out all plausible lower redshift solutions for this source. The only viable solution is that… ▽ More We present Hubble WFC3/IR slitless grism spectra of a remarkably bright $z\gtrsim10$ galaxy candidate, GN-z11, identified initially from CANDELS/GOODS-N imaging data. A significant spectroscopic continuum break is detected at $λ=1.47\pm0.01~μ$m. The new grism data, combined with the photometric data, rule out all plausible lower redshift solutions for this source. The only viable solution is that this continuum break is the Ly$α$ break redshifted to ${z_\mathrm{grism}=11.09^{+0.08}_{-0.12}}$, just $\sim$400 Myr after the Big Bang. This observation extends the current spectroscopic frontier by 150 Myr to well before the Planck (instantaneous) cosmic reionization peak at z~8.8, demonstrating that galaxy build-up was well underway early in the reionization epoch at z>10. GN-z11 is remarkably and unexpectedly luminous for a galaxy at such an early time: its UV luminosity is 3x larger than L* measured at z~6-8. The Spitzer IRAC detections up to 4.5 $μ$m of this galaxy are consistent with a stellar mass of ${\sim10^{9}~M_\odot}$. This spectroscopic redshift measurement suggests that the James Webb Space Telescope (JWST) will be able to similarly and easily confirm such sources at z>10 and characterize their physical properties through detailed spectroscopy. Furthermore, WFIRST, with its wide-field near-IR imaging, would find large numbers of similar galaxies and contribute greatly to JWST's spectroscopy, if it is launched early enough to overlap with JWST. △ Less

Submitted 1 March, 2016; originally announced March 2016.

Comments: 12 pages, 11 figures, accepted for publication in ApJ

Bright galaxies at Hubble's redshift detection frontier: Preliminary results and design from the redshift z~9-10 BoRG pure-parallel HST survey

Authors: V. Calvi, M. Trenti, M. Stiavelli, P. Oesch, L. D. Bradley, K. B. Schmidt, D. Coe, G. Brammer, S. Bernard, R. J. Bouwens, D. Carrasco, C. M. Carollo, B. W. Holwerda, J. W. MacKenty, C. A. Mason, J. M. Shull, T. Treu

Abstract: We present the first results and design from the redshift z~9-10 Brightest of the Reionizing Galaxies {\it Hubble Space Telescope} survey BoRG[z9-10], aimed at searching for intrinsically luminous unlensed galaxies during the first 700 Myr after the Big Bang. BoRG[z9-10] is the continuation of a multi-year pure-parallel near-IR and optical imaging campaign with the Wide Field Camera 3. The ongoing… ▽ More We present the first results and design from the redshift z~9-10 Brightest of the Reionizing Galaxies {\it Hubble Space Telescope} survey BoRG[z9-10], aimed at searching for intrinsically luminous unlensed galaxies during the first 700 Myr after the Big Bang. BoRG[z9-10] is the continuation of a multi-year pure-parallel near-IR and optical imaging campaign with the Wide Field Camera 3. The ongoing survey uses five filters, optimized for detecting the most distant objects and offering continuous wavelength coverage from λ=0.35μm to λ=1.7μm. We analyze the initial ~130 arcmin$^2$ of area over 28 independent lines of sight (~25% of the total planned) to search for z>7 galaxies using a combination of Lyman break and photometric redshift selections. From an effective comoving volume of (5-25) $times 10^5$ Mpc$^3$ for magnitudes brighter than $m_{AB}=26.5-24.0$ in the $H_{160}$-band respectively, we find five galaxy candidates at z~8.3-10 detected at high confidence (S/N>8), including a source at z~8.4 with mAB=24.5 (S/N~22), which, if confirmed, would be the brightest galaxy identified at such early times (z>8). In addition, BoRG[z9-10] data yield four galaxies with $7.3 \lesssim z \lesssim 8$. These new Lyman break galaxies with m$\lesssim26.5$ are ideal targets for follow-up observations from ground and space based observatories to help investigate the complex interplay between dark matter growth, galaxy assembly, and reionization. △ Less

Submitted 16 December, 2015; originally announced December 2015.

Comments: Accepted for publication on ApJ. 21 pages, 11 figures, 4 tables

The Lyman-continuum photon production efficiency in the high-redshift Universe

Authors: Stephen M. Wilkins, Yu Feng, Tiziana Di Matteo, Rupert Croft, Elizabeth R. Stanway, Rychard J. Bouwens, Peter Thomas

Abstract: The Lyman Continuum photon production efficiency ($ξ_{\rm ion}$) is a critical ingredient for inferring the number of photons available to reionise the intergalactic medium. To estimate the theoretical production efficiency in the high-redshift Universe we couple the BlueTides cosmological hydrodynamical simulation with a range of stellar population synthesis models. We find Lyman Continuum photon… ▽ More The Lyman Continuum photon production efficiency ($ξ_{\rm ion}$) is a critical ingredient for inferring the number of photons available to reionise the intergalactic medium. To estimate the theoretical production efficiency in the high-redshift Universe we couple the BlueTides cosmological hydrodynamical simulation with a range of stellar population synthesis models. We find Lyman Continuum photon production efficiencies of $\log_{10}(ξ_{\rm ion}/{\rm erg^{-1}\, Hz})\approx 25.1-25.5$ depending on the choice of stellar population synthesis model. These results are broadly consistent with recent observational constraints at high-redshift though favour a model incorporating the effects of binary evolution △ Less

Submitted 14 March, 2016; v1 submitted 10 December, 2015; originally announced December 2015.

Comments: 4 pages, 2 figures, accepted to MNRAS. Minor changes reflecting referee comments and community feedback

Inferred Hα Flux as a Star-Formation Rate Indicator at z ~ 4-5: Implications for Dust Properties, Burstiness, and the z = 4-8 Star-Formation-Rate Functions

Authors: Renske Smit, Rychard J. Bouwens, Ivo Labbé, Marijn Franx, Stephen M. Wilkins, Pascal A. Oesch

Abstract: We derive Hα fluxes for a large spectroscopic and photometric-redshift-selected sample of sources over GOODS-North and South in the redshift range z = 3.8-5.0 with deep HST, Spitzer/IRAC, and ground-based observations. The Hα flux is inferred based on the offset between the IRAC 3.6 μm flux and that predicted from the best-fit SED. We demonstrate that the Hα flux correlates well with dust- correct… ▽ More We derive Hα fluxes for a large spectroscopic and photometric-redshift-selected sample of sources over GOODS-North and South in the redshift range z = 3.8-5.0 with deep HST, Spitzer/IRAC, and ground-based observations. The Hα flux is inferred based on the offset between the IRAC 3.6 μm flux and that predicted from the best-fit SED. We demonstrate that the Hα flux correlates well with dust- corrected UV star-formation rate (SFR) and therefore can serve as an independent SFR indicator. However, we also find a systematic offset in the SFR_Hα/SFR_UV ratios for z ~ 4-5 galaxies relative to local relations (assuming the same dust corrections for nebular regions and stellar light). We show that we can resolve the modest tension in the inferred SFRs by assuming bluer intrinsic UV slopes (increasing the dust correction), a rising star-formation history or assuming a low metallicity stellar population with a hard ionizing spectrum (increasing the L_Hα/SFR ratio). Using Hα as a SFR indicator, we find a normalization of the star formation main sequence in good agreement with recent SED-based determinations and also derive the SFR functions at z ~ 4-8. In addition, we assess for the first time the burstiness of star formation in z ~ 4 galaxies on <100 Myr time scales by comparing UV and Hα-based sSFRs; their one-to-one relationship argues against significantly bursty star-formation histories. Further progress will be made on these results, by incorporating new results from ALMA to constrain the dust-obscured star formation in high-redshift UV-selected samples. △ Less

Submitted 1 November, 2016; v1 submitted 27 November, 2015; originally announced November 2015.

Comments: 22 pages, 13 figures, 9 tables, ApJ, in press

The Lyman-Continuum Photon Production Efficiency xi_{ion} of z~4-5 Galaxies from IRAC-based Halpha Measurements: Implications for the Escape Fraction and Cosmic Reionization

Authors: R. J. Bouwens, R. Smit, I. Labbe, M. Franx, J. Caruana, P. Oesch, M. Stefanon, N. Rasappu

Abstract: Galaxies represent one of the preferred candidate sources to drive the reionization of the universe. Even as gains are made in mapping the galaxy UV luminosity density to z>6, significant uncertainties remain regarding the conversion to the implied ionizing emissivity. The relevant unknowns are the Lyman-continuum (LyC) photon production efficiency xi_{ion} and the escape fraction f_{esc}. As we s… ▽ More Galaxies represent one of the preferred candidate sources to drive the reionization of the universe. Even as gains are made in mapping the galaxy UV luminosity density to z>6, significant uncertainties remain regarding the conversion to the implied ionizing emissivity. The relevant unknowns are the Lyman-continuum (LyC) photon production efficiency xi_{ion} and the escape fraction f_{esc}. As we show here, the first of these unknowns is directly measureable in z=4-5 galaxies, based on the impact the Halpha line has on the observed IRAC fluxes. By computing a LyC photon production rate from the implied Halpha luminosities for a broad selection of z=4-5 galaxies and comparing this against the dust-corrected UV-continuum luminosities, we provide the first-ever direct estimates of the LyC photon production efficiency xi_{ion} for the z>~4 galaxy population. We find log_{10} xi_{ion}/[Hz/ergs] to have a mean value of 25.27_{-0.03}^{+0.03} and 25.34_{-0.02}^{+0.02} for sub-L* z=4-5 galaxies adopting Calzetti and SMC dust laws, respectively. Reassuringly, both values are consistent with standardly assumed xi_{ion}'s in reionization models, with a slight preference for higher xi_{ion}'s (by ~0.1 dex) adopting the SMC dust law. A modest ~0.03-dex increase in these estimates would result if the escape fraction for ionizing photons is non-zero and galaxies dominate the ionizing emissivity at z~4.4. High values of xi_{ion} (~25.5-25.8 dex) are derived for the bluest galaxies (beta<-2.3) in our samples, independent of dust law and consistent with results for a z=7.045 galaxy. Such elevated values of xi_{ion} would have important consequences, indicating that f_{esc} cannot be in excess of 13% unless the galaxy UV luminosity function does not extend down to -13 mag or the clumping factor is greater than 3. A low escape fraction would fit well with the low rate of LyC leakage observed at z~3. △ Less

Submitted 1 June, 2016; v1 submitted 26 November, 2015; originally announced November 2015.

Comments: 12 pages, 6 figures, 3 tables, updated to match the version in press, both statistical and systematic errors given for xi_{ion}

The clustering and halo occupation distribution of Lyman-break galaxies at $z\sim4$

Authors: Jaehong Park, Han-Seek Kim, J. Stuart B. Wyithe, C. G. Lacey, C. M. Baugh, R. L. Barone-Nugent, M. Trenti, R. J. Bouwens

Abstract: We investigate the clustering of Lyman-break galaxies (LBGs) at $z\sim4$. Using the hierarchical galaxy formation model GALFORM, we predict, for the first time using a semi-analytical model with feedback from active galactic nuclei (AGN), the angular correlation function (ACF) of LBGs and find agreement within $3\,σ$ with new measurements of the ACF from surveys including the Hubble eXtreme Deep F… ▽ More We investigate the clustering of Lyman-break galaxies (LBGs) at $z\sim4$. Using the hierarchical galaxy formation model GALFORM, we predict, for the first time using a semi-analytical model with feedback from active galactic nuclei (AGN), the angular correlation function (ACF) of LBGs and find agreement within $3\,σ$ with new measurements of the ACF from surveys including the Hubble eXtreme Deep Field (XDF) and CANDELS field. Our simulations confirm the conclusion reached using independent models that although the predicted ACFs reproduce the trend of increased clustering with luminosity, the dependence is less strong than observed. We find that for the detection limits of the XDF field central LBGs at $z\sim 4$ predominantly reside in haloes of mass $\sim 10^{11}-10^{12}h^{-1}M_{\rm \odot}$ and that satellites reside in larger haloes of mass $\sim 10^{12}-10^{13}h^{-1}M_{\rm \odot}$. The model predicts fewer bright satellite LBGs at $z\sim4$ than is inferred from measurements of the ACF at small scales. By analysing the halo occupation distribution (HOD) predicted by the model, we find evidence that AGN feedback affects the HOD of central LBGs in massive haloes. This is a new high-redshift test of this important feedback mechanism. We investigate the effect of photometric errors in the observations on the ACF predictions. We find that the observational uncertainty in the galaxy luminosity reduces the clustering amplitude and that this effect increases towards faint galaxies, particularly on small scales. To compare properties of model with observed LBGs this uncertainty must be considered. △ Less

Submitted 8 June, 2016; v1 submitted 5 November, 2015; originally announced November 2015.

Comments: 16 pages, 13 figures, accepted for publication in MNRAS

High-Redshift Galaxy Surveys and the Reionization of the Universe

Authors: Rychard J. Bouwens

Abstract: Star-forming galaxies in the early universe provide us with perhaps the most natural way of explaining the reionization of the universe. Current observational results are sufficiently comprehensive, as to allow us to approximately calculate how the ionizing radiation from galaxies varies as a function of cosmic time. Important uncertainties in modeling reionization by galaxies revolve around the e… ▽ More Star-forming galaxies in the early universe provide us with perhaps the most natural way of explaining the reionization of the universe. Current observational results are sufficiently comprehensive, as to allow us to approximately calculate how the ionizing radiation from galaxies varies as a function of cosmic time. Important uncertainties in modeling reionization by galaxies revolve around the escape fraction and its luminosity and redshift dependence, a possible truncation of the galaxy luminosity function at the faint end, and an evolution in the production efficiency of Lyman-continuum photons with cosmic time. Despite these uncertainties, plausible choices for these parameters naturally predict a cosmic ionizing emissivity at z~6-10 whose evolution and overall normalization is in excellent agreement with that derived from current observational constraints. This strongly suggests that galaxies provide the necessary photons to reionize the universe. △ Less

Submitted 15 November, 2015; v1 submitted 3 November, 2015; originally announced November 2015.

Comments: 37 pages, 14 figures, Book chapter in "Understanding the Epoch of Cosmic Reionization: Challenges and Progress", Springer International Publishing, Ed. Andrei Mesinger, ISBN 978-3-319-21956-1, corrected a few minor typos

Quantifying the UV continuum slopes of galaxies to z~10 using deep Hubble and Spitzer/IRAC observations

Authors: Stephen M. Wilkins, Rychard J. Bouwens, Pascal A. Oesch, Ivo Labbe, Mark Sargent, Joseph Caruana, Julie Wardlow, Scott Clay

Abstract: Measurements of the UV-continuum slopes provide valuable information on the physical properties of galaxies forming in the early universe, probing the dust reddening, age, metal content, and even the escape fraction. While constraints on these slopes generally become more challenging at higher redshifts as the UV continuum shifts out of the Hubble Space Telescope bands (particularly at z>7), such… ▽ More Measurements of the UV-continuum slopes provide valuable information on the physical properties of galaxies forming in the early universe, probing the dust reddening, age, metal content, and even the escape fraction. While constraints on these slopes generally become more challenging at higher redshifts as the UV continuum shifts out of the Hubble Space Telescope bands (particularly at z>7), such a characterisation actually becomes abruptly easier for galaxies in the redshift window z=9.5-10.5 due to the Spitzer/IRAC 3.6um-band probing the rest-UV continuum and the long wavelength baseline between this Spitzer band and the Hubble H-band. Higher S/N constraints on the UV slope are possible at z~10 than at z=8. Here we take advantage of this opportunity and five recently discovered bright z=9.5-10.5 galaxies to present the first measurements of the mean slope for a multi-object sample of galaxy candidates at z~10. We find the measured observed slopes of these candidates are $-2.1\pm0.3\pm0.2$ (random and systematic), only slightly bluer than the measured slopes at 3.5<z<7.5 for galaxies of similar luminosities. Small increases in the stellar ages, metallicities, and dust content of the galaxy population from z~10 to z~7 could easily explain the apparent evolution in slopes. △ Less

Submitted 6 October, 2015; originally announced October 2015.

Comments: 10 pages, 10 figures, accepted in MNRAS

Ultradeep IRAC Imaging Over The HUDF And GOODS-South: Survey Design And Imaging Data Release

Authors: I. Labbe, P. A. Oesch, G. D. Illingworth, P. G. van Dokkum, R. J. Bouwens, M. Franx, C. M. Carollo, M. Trenti, B. Holden, R. Smit, V. Gonzalez, D. Magee, M. Stiavelli, M. Stefanon

Abstract: The IRAC ultradeep field (IUDF) and IRAC Legacy over GOODS (IGOODS) programs are two ultradeep imaging surveys at 3.6μm and 4.5μ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… ▽ More The IRAC ultradeep field (IUDF) and IRAC Legacy over GOODS (IGOODS) programs are two ultradeep imaging surveys at 3.6μm and 4.5μ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μm and 4.5μ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σ point source sensitivities as faint as of 15 nJy (28.5 AB) at 3.6μm and 18 nJy (28.3 AB) at 4.5μm. The value of such ultradeep IRAC data is illustrated by direct detections of z = 7 - 8 galaxies as faint as HAB = 28. △ Less

Submitted 29 July, 2015; originally announced July 2015.

Comments: Accepted for publication in ApJS

The Bright End of the z~9 and z~10 UV Luminosity Functions using all five CANDELS Fields

Authors: R. J. Bouwens, P. A. Oesch, I. Labbe, G. D. Illingworth, G. G. Fazio, D. Coe, B. Holwerda, R. Smit, M. Stefanon, P. G. van Dokkum, M. Trenti, M. L. N. Ashby, J. -S. Huang, L. Spitler, C. Straatman, L. Bradley, D. Magee

Abstract: 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 th… ▽ More 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. △ Less

Submitted 11 August, 2016; v1 submitted 2 June, 2015; originally announced June 2015.

Comments: 23 pages, 15 figures, 7 tables, updated to match the version in press, including some minor textual corrections identified at the proof stage

Z > 7 galaxies with red Spitzer/IRAC [3.6]-[4.5] colors in the full CANDELS data set: the brightest-known galaxies at Z ~ 7-9 and a probable spectroscopic confirmation at Z=7.48

Authors: G. W. Roberts-Borsani, R. J. Bouwens, P. A. Oesch, I. Labbe, R. Smit, G. D. Illingworth, P. van Dokkum, B. Holden, V. Gonzalez, M. Stefanon, B. Holwerda, S. Wilkins

Abstract: We identify 4 unusually bright (H < 25.5) galaxies from HST and Spitzer CANDELS data with probable redshifts z ~ 7-9. These identifications include the brightest-known galaxies to date at z > 7.5. As Y-band observations are not available over the full CANDELS program to perform a standard Lyman-break selection of z > 7 galaxies, we employ an alternate strategy using deep Spitzer/IRAC data. We iden… ▽ More We identify 4 unusually bright (H < 25.5) galaxies from HST and Spitzer CANDELS data with probable redshifts z ~ 7-9. These identifications include the brightest-known galaxies to date at z > 7.5. As Y-band observations are not available over the full CANDELS program to perform a standard Lyman-break selection of z > 7 galaxies, we employ an alternate strategy using deep Spitzer/IRAC data. We identify z ~ 7.1 - 9.1 galaxies by selecting z >~ 6 galaxies from the HST CANDELS data that show quite red IRAC [3.6]-[4.5] colors, indicating strong [OIII]+Hbeta lines in the 4.5 micron band. This selection strategy was validated using a modest sample for which we have deep Y-band coverage, and subsequently used to select the brightest z > 7 sources. Applying the IRAC criteria to all HST-selected optical-dropout galaxies over the full ~900 arcmin**2 of the CANDELS survey revealed four unusually bright z ~ 7.1, 7.6, 7.9 and 8.6 candidates. The median [3.6]-[4.5] color of our selected z ~ 7.1-9.1 sample is consistent with rest-frame [OIII]+Hbeta EWs of ~1500A, in the [4.5] band. Keck/MOSFIRE spectroscopy has been independently reported for two of our selected sources, showing Ly-alpha at redshifts of 7.7302+/-0.0006 and 8.683^+0.001_-0.004, respectively. We present similar Keck/MOSFIRE spectroscopy for a third selected galaxy with a probable 4.7sigma Ly-alpha line at z_spec=7.4770+/-0.0008. All three have H-band magnitudes of ~25 mag and are ~0.5 mag more luminous (M(UV) ~ -22.0) than any previously discovered z ~ 8 galaxy, with important implications for the UV LF. Our 3 brightest, highest redshift z > 7 galaxies all lie within the CANDELS EGS field, providing a dramatic illustration of the potential impact of field-to-field variance. △ Less

Submitted 1 April, 2016; v1 submitted 2 June, 2015; originally announced June 2015.

Comments: 19 pages, 13 figures, 5 tables, accepted for publication in the Astrophysical Journal. Updated to match the version in press and reflect spectroscopic confirmation of the highest redshift candidate in our sample by Zitrin+2015

Report number: ApJ99302

Reionization after Planck: The Derived Growth of the Cosmic Ionizing Emissivity now matches the Growth of the Galaxy UV Luminosity Density

Authors: R. J. Bouwens, G. D. Illingworth, P. A. Oesch, J. Caruana, B. Holwerda, R. Smit, S. Wilkins

Abstract: 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 measu… ▽ More 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. △ Less

Submitted 5 September, 2015; v1 submitted 27 March, 2015; originally announced March 2015.

Comments: 20 pages, 12 figures, 5 tables, Astrophysical Journal, updated to match version in press, Figure 6 shows the main result of the paper

A Spectroscopic Redshift Measurement for a Luminous Lyman Break Galaxy at z=7.730 using Keck/MOSFIRE

Authors: P. A. Oesch, P. G. van Dokkum, G. D. Illingworth, R. J. Bouwens, I. Momcheva, B. Holden, G. W. Roberts-Borsani, R. Smit, M. Franx, I. Labbe, V. Gonzalez, D. Magee

Abstract: 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 reli… ▽ More 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. △ Less

Submitted 3 May, 2015; v1 submitted 18 February, 2015; originally announced February 2015.

Comments: 6 pages, 4 figures, small updates to match ApJL accepted version

High-precision Photometric Redshifts from Spitzer/IRAC: Extreme [3.6]-[4.5] Colors Identify Galaxies in the Redshift Range z~6.6-6.9

Authors: Renske Smit, Rychard J. Bouwens, Marijn Franx, Pascal A. Oesch, Matthew L. N. Ashby, S. P. Willner, Ivo Labbe, Benne Holwerda, Giovanni G. Fazio, J. -S. Huang

Abstract: One of the most challenging aspects of studying galaxies in the z>~7 universe is the infrequent confirmation of their redshifts through spectroscopy, a phenomenon thought to occur from the increasing opacity of the intergalactic medium to Lya photons at z>6.5. The resulting redshift uncertainties inhibit the efficient search for [C II] in z~7 galaxies with sub-mm instruments such as ALMA, given th… ▽ More One of the most challenging aspects of studying galaxies in the z>~7 universe is the infrequent confirmation of their redshifts through spectroscopy, a phenomenon thought to occur from the increasing opacity of the intergalactic medium to Lya photons at z>6.5. The resulting redshift uncertainties inhibit the efficient search for [C II] in z~7 galaxies with sub-mm instruments such as ALMA, given their limited scan speed for faint lines. One means by which to improve the precision of the inferred redshifts is to exploit the potential impact of strong nebular emission lines on the colors of z~4-8 galaxies as observed by Spitzer/IRAC. At z~6.8, galaxies exhibit IRAC colors as blue as [3.6]-[4.5] ~-1, likely due to the contribution of [O III]+Hb to the 3.6 mum flux combined with the absence of line contamination in the 4.5 mum band. In this paper we explore the use of extremely blue [3.6]-[4.5] colors to identify galaxies in the narrow redshift window z~6.6-6.9. When combined with an I-dropout criterion, we demonstrate that we can plausibly select a relatively clean sample of z~6.8 galaxies. Through a systematic application of this selection technique to our catalogs from all five CANDELS fields, we identify 20 probable z~6.6-6.9 galaxies. We estimate that our criteria select the ~50% strongest line emitters at z~6.8 and from the IRAC colors we estimate a typical [O III]+Hb rest-frame equivalent width of 1085A for this sample. The small redshift uncertainties on our sample make it particularly well suited for follow-up studies with facilities such as ALMA. △ Less

Submitted 1 December, 2014; originally announced December 2014.

Comments: In submission to the Astrophysical Journal, updated in response to the referee report, 13 pages, 11 figures, 1 table

First Frontier Field Constraints on the Cosmic Star-Formation Rate Density at z~10 - The Impact of Lensing Shear on Completeness of High-Redshift Galaxy Samples

Authors: P. A. Oesch, R. J. Bouwens, G. D. Illingworth, M. Franx, S. M. Ammons, P. G. van Dokkum, M. Trenti, I. Labbe

Abstract: 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… ▽ More 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. △ Less

Submitted 21 July, 2015; v1 submitted 3 September, 2014; originally announced September 2014.

Comments: 10 pages, 7 figures, changed to match accepted version to appear in ApJ

Measurement of Galaxy Clustering at z~7.2 and the Evolution of Galaxy Bias from 3.8<z<8 in the XDF, GOODS-S AND GOODS-N

Authors: R. L. Barone-Nugent, M. Trenti, J. S. B. Wyithe, R. J. Bouwens, P. A. Oesch, G. D. Illingworth, C. M. Carollo, J. Su, M. Stiavelli, I. Labbe, P. G. van Dokkum

Abstract: Lyman-Break Galaxy (LBG) samples observed during reionization ($z\gtrsim6$) with Hubble Space Telescope's Wide Field Camera 3 are reaching sizes sufficient to characterize their clustering properties. Using a combined catalog from the Hubble eXtreme Deep Field and CANDELS surveys, containing $N=743$ LBG candidates at z>6.5 at a mean redshift of $z=7.2$, we detect a clear clustering signal in the a… ▽ More Lyman-Break Galaxy (LBG) samples observed during reionization ($z\gtrsim6$) with Hubble Space Telescope's Wide Field Camera 3 are reaching sizes sufficient to characterize their clustering properties. Using a combined catalog from the Hubble eXtreme Deep Field and CANDELS surveys, containing $N=743$ LBG candidates at z>6.5 at a mean redshift of $z=7.2$, we detect a clear clustering signal in the angular correlation function (ACF) at $\sim4σ$, corresponding to a real-space correlation length $r_{0}=6.7^{+0.9}_{-1.0}h^{-1}$cMpc. The derived galaxy bias $b=8.6^{+0.9}_{-1.0}$ is that of dark-matter halos of $M=10^{11.1^{+0.2}_{-0.3}}$M$_{\odot}$ at $z=7.2$, and highlights that galaxies below the current detection limit ($M_{AB}\sim-17.7$) are expected in lower-mass halos ($M\sim10^{8}-10^{10.5}$M$_{\odot}$). We compute the ACF of LBGs at $z\sim3.8-z\sim5.9$ in the same surveys. A trend of increasing bias is found from $z=3.8$ ($b\sim3.0$) to $z=7.2$ ($b\sim8.6$), broadly consistent with galaxies at fixed luminosity being hosted in dark-matter halos of similar mass at $4<z<6$, followed by a slight rise in halo masses at $z\sim7$ ($\sim2σ$ confidence). Separating the data at the median luminosity of the $z=7.2$ sample ($M_{UV}=-19.4$) shows higher clustering at $z=5.9$ for bright galaxies ($r_{0}=5.5^{+1.4}_{-1.5}h^{-1}$cMpc, $b=6.2^{+1.2}_{-1.5}$) compared to faint galaxies ($r_{0}=1.9^{+1.1}_{-1.0}h^{-1}$cMpc, $b=2.7\pm1.2$) implying a constant mass-to-light ratio $\frac{dlogM}{dlogL}\sim1.2^{+1.8}_{-0.8}$. A similar trend is present in the $z=7.2$ sample with larger uncertainty. Finally, our bias measurements allow us to investigate the fraction of dark-matter halos hosting UV-bright galaxies (the duty-cycle, $ε_{DC}$). At $z=7.2$ values near unity are preferred, which may be explained by the shortened halo assembly time at high-redshift. △ Less

Submitted 28 July, 2014; originally announced July 2014.

Comments: 7 pages, 2 tables, 5 figures, accepted for publication in The Astrophysical Journal

UV Luminosity Functions at redshifts z~4 to z~10: 10000 Galaxies from HST Legacy Fields

Authors: R. J. Bouwens, G. D. Illingworth, P. A. Oesch, M. Trenti, I. Labbe', L. Bradley, M. Carollo, P. G. van Dokkum, V. Gonzalez, B. Holwerda, M. Franx, L. Spitler, R. Smit, D. Magee

Abstract: The remarkable HST datasets from the CANDELS, HUDF09, HUDF12, ERS, and BoRG/HIPPIES programs have allowed us to map out the evolution of the UV LF from z~10 to z~4. We have identified 5859, 3001, 857, 481, 217, and 6 galaxy candidates at z~4, z~5, z~6, z~7, z~8, and z~10, respectively from the ~1000 arcmin**2 area probed. The selection of z~4-8 galaxies over the five CANDELS fields allows us to as… ▽ More The remarkable HST datasets from the CANDELS, HUDF09, HUDF12, ERS, and BoRG/HIPPIES programs have allowed us to map out the evolution of the UV LF from z~10 to z~4. We have identified 5859, 3001, 857, 481, 217, and 6 galaxy candidates at z~4, z~5, z~6, z~7, z~8, and z~10, respectively from the ~1000 arcmin**2 area probed. The selection of z~4-8 galaxies over the five CANDELS fields allows us to assess the cosmic variance; the largest variations are apparent at z>=7. Our new LF determinations at z~4 and z~5 span a 6-mag baseline (-22.5 to -16 AB mag). These determinations agree well with previous estimates, but the larger samples and volumes probed here result in a more reliable sampling of >L* galaxies and allow us to reassess the form of the UV LFs. Our new LF results strengthen our earlier findings to 3.4 sigma significance for a steeper faint-end slope to the UV LF at z>4, with alpha evolving from alpha=-1.64+/-0.04 at z~4 to alpha=-2.06+/-0.13 at z~7 (and alpha = -2.02+/-0.23 at z~8), consistent with that expected from the evolution of the halo mass function. With our improved constraints at the bright end, we find less evolution in the characteristic luminosity M* over the redshift range z~4 to z~7; the observed evolution in the LF is now largely represented by changes in phi*. No evidence for a non-Schechter-like form to the z~4-8 LFs is found. A simple conditional LF model based on halo growth and evolution in the M/L ratio of halos ((1+z)**-1.5) provides a good representation of the observed evolution. △ Less

Submitted 22 January, 2015; v1 submitted 17 March, 2014; originally announced March 2014.

Comments: 53 pages, 28 figures, 11 tables, accepted for publication in the Astrophysical Journal

The Most Luminous z~9-10 Galaxy Candidates yet Found: The Luminosity Function, Cosmic Star-Formation Rate, and the First Mass Density Estimate at 500 Myr

Authors: P. A. Oesch, R. J. Bouwens, G. D. Illingworth, I. Labbe, R. Smit, P. G. van Dokkum, I. Momcheva, M. L. N. Ashby, G. G. Fazio, J. Huang, S. P. Willner, V. Gonzalez, D. Magee, M. Trenti, G. B. Brammer, R. E. Skelton, L. R. Spitler

Abstract: [abridged] We present the discovery of four surprisingly bright (H_160 ~ 26 - 27 mag AB) galaxy candidates at z~9-10 in the complete HST CANDELS WFC3/IR GOODS-N imaging data, doubling the number of z~10 galaxy candidates that are known, just ~500 Myr after the Big Bang. Two similarly bright sources are also detected in a systematic re-analysis of the GOODS-S data set. Three of the four galaxies in… ▽ More [abridged] We present the discovery of four surprisingly bright (H_160 ~ 26 - 27 mag AB) galaxy candidates at z~9-10 in the complete HST CANDELS WFC3/IR GOODS-N imaging data, doubling the number of z~10 galaxy candidates that are known, just ~500 Myr after the Big Bang. Two similarly bright sources are also detected in a systematic re-analysis of the GOODS-S data set. Three of the four galaxies in GOODS-N are significantly detected at 4.5-6.2sigma in the very deep Spitzer/IRAC 4.5 micron data, as is one of the GOODS-S candidates. Furthermore, the brightest of our candidates (at z=10.2+-0.4) is robustly detected also at 3.6 micron (6.9sigma), revealing a flat UV spectral energy distribution with a slope beta=-2.0+-0.2, consistent with demonstrated trends with luminosity at high redshift. The abundance of such luminous candidates suggests that the luminosity function evolves more significantly in phi_* than in L_* at z>~8 with a higher number density of bright sources than previously expected. Despite the discovery of these luminous candidates, the cosmic star formation rate density for galaxies with SFR >0.7 M_sun/yr shows an order-of-magnitude increase in only 170 Myr from z ~ 10 to z ~ 8, consistent with previous results. Based on the IRAC detections, we derive galaxy stellar masses at z~10, finding that these luminous objects are typically 10^9 M_sun. The cosmic stellar mass density at z~10 is log10 rho_* = 4.7^+0.5_-0.8 M_sun Mpc^-3 for galaxies brighter than M_UV~-18. The remarkable brightness, and hence luminosity, of these z~9-10 candidates highlights the opportunity for deep spectroscopy to determine their redshift and nature, demonstrates the value of additional search fields covering a wider area to understand star-formation in the very early universe, and highlights the opportunities for JWST to map the buildup of galaxies at redshifts much earlier than z~10. △ Less

Submitted 6 January, 2014; v1 submitted 9 September, 2013; originally announced September 2013.

Comments: 20 pages, 12 figures, changed to match resubmitted version to ApJ

Evidence for Ubiquitous, High-EW Nebular Emission in z~7 Galaxies: Towards a Clean Measurement of the Specific Star Formation Rate using a Sample of Bright, Magnified Galaxies

Authors: R. Smit, R. J. Bouwens, I. Labbe, W. Zheng, L. Bradley, M. Donahue, D. Lemze, J. Moustakas, K. Umetsu, A. Zitrin, D. Coe, M. Postman, V. Gonzalez, M. Bartelmann, N. Benitez, T. Broadhurst, H. Ford, C. Grillo, L. Infante, Y. Jimenez-Teja, S. Jouvel, D. D. Kelson, O. Lahav, D. Maoz, E. Medezinski , et al. (8 additional authors not shown)

Abstract: Growing observational evidence now indicates that nebular line emission has a significant impact on the rest-frame optical fluxes of z~5-7 galaxies observed with Spitzer. This line emission makes z~5-7 galaxies appear more massive, with lower specific star formation rates. However, corrections for this line emission have been very difficult to perform reliably due to huge uncertainties on the over… ▽ More Growing observational evidence now indicates that nebular line emission has a significant impact on the rest-frame optical fluxes of z~5-7 galaxies observed with Spitzer. This line emission makes z~5-7 galaxies appear more massive, with lower specific star formation rates. However, corrections for this line emission have been very difficult to perform reliably due to huge uncertainties on the overall strength of such emission at z>~5.5. Here, we present the most direct observational evidence yet for ubiquitous high-EW [OIII]+Hbeta line emission in Lyman-break galaxies at z~7, while also presenting a strategy for an improved measurement of the sSFR at z~7. We accomplish this through the selection of bright galaxies in the narrow redshift window z~6.6-7.0 where the IRAC 4.5 micron flux provides a clean measurement of the stellar continuum light. Observed 4.5 micron fluxes in this window contrast with the 3.6 micron fluxes which are contaminated by the prominent [OIII]+Hbeta lines. To ensure a high S/N for our IRAC flux measurements, we consider only the brightest (H_{160}<26 mag) magnified galaxies we have identified in CLASH and other programs targeting galaxy clusters. Remarkably, the mean rest-frame optical color for our bright seven-source sample is very blue, [3.6]-[4.5]=-0.9+/-0.3. Such blue colors cannot be explained by the stellar continuum light and require that the rest-frame EW of [OIII]+Hbeta be greater than 637 Angstroms for the average source. The bluest four sources from our seven-source sample require an even more extreme EW of 1582 Angstroms. Our derived lower limit for the mean [OIII]+Hbeta EW could underestimate the true EW by ~2x based on a simple modeling of the redshift distribution of our sources. We can also set a robust lower limit of >~4 Gyr^-1 on the specific star formation rates based on the mean SED for our seven-source sample. (abridged) △ Less

Submitted 22 July, 2013; originally announced July 2013.

Comments: 9 pages, 6 figures, 1 table, submitted to the Astrophysical Journal

UV-Continuum Slopes of >4000 z~4-8 Galaxies from the HUDF/XDF, HUDF09, ERS, CANDELS-South, and CANDELS-North Fields

Authors: R. J. Bouwens, G. D. Illingworth, P. A. Oesch, I. Labbe, P. G. van Dokkum, M. Trenti, M. Franx, R. Smit, V. Gonzalez, D. Magee

Abstract: We measure the UV-continuum slope beta for over 4000 high-redshift galaxies over a wide range of redshifts z~4-8 and luminosities from the HST HUDF/XDF, HUDF09-1, HUDF09-2, ERS, CANDELS-N, and CANDELS-S data sets. Our new beta results reach very faint levels at z~4 (-15.5 mag: 0.006 L*(z=3)), z~5 (-16.5 mag: 0.014L*(z=3)), and z~6 and z~7 (-17 mag: 0.025 L*(z=3)). Inconsistencies between previous… ▽ More We measure the UV-continuum slope beta for over 4000 high-redshift galaxies over a wide range of redshifts z~4-8 and luminosities from the HST HUDF/XDF, HUDF09-1, HUDF09-2, ERS, CANDELS-N, and CANDELS-S data sets. Our new beta results reach very faint levels at z~4 (-15.5 mag: 0.006 L*(z=3)), z~5 (-16.5 mag: 0.014L*(z=3)), and z~6 and z~7 (-17 mag: 0.025 L*(z=3)). Inconsistencies between previous studies led us to conduct a comprehensive review of systematic errors and develop a new technique for measuring beta that is robust against biases that arise from the impact of noise. We demonstrate, by object-by-object comparisons, that all previous studies, including our own and those done on the latest HUDF12 dataset, suffer from small systematic errors in beta. We find that after correcting for the systematic errors (typically d(beta) ~0.1-0.2) all beta results at z~7 from different groups are in excellent agreement. The mean beta we measure for faint (-18 mag: 0.1L*(z=3)) z~4, z~5, z~6, and z~7 galaxies is -2.03+/-0.03+/-0.06 (random and systematic errors), -2.14+/-0.06+/-0.06, -2.24+/-0.11+/-0.08, and -2.30+/-0.18+/-0.13, respectively. Our new beta values are redder than we have reported in the past, but bluer than other recent results. Our previously reported trend of bluer beta's at lower luminosities is confirmed, as is the evolution to bluer beta's at high redshifts. beta appears to show only a mild luminosity dependence faintward of M(UV,AB) ~ -19 mag, suggesting that the mean beta asymptotes to ~ -2.2 to -2.4 for faint z>~4 galaxies. At z~7, the observed beta's suggest non-zero, but low dust extinction, and they agree well with values predicted in cosmological hydrodynamical simulations. △ Less

Submitted 3 August, 2014; v1 submitted 12 June, 2013; originally announced June 2013.

Comments: 41 pages, 28 figures, 9 tables, accepted for publication in the Astrophysical Journal, makes use of the full 10-epochs of CANDELS North data

The HST eXtreme Deep Field XDF: Combining all ACS and WFC3/IR Data on the HUDF Region into the Deepest Field Ever

Authors: G. D. Illingworth, D. Magee, P. A. Oesch, R. J. Bouwens, I. Labbe, M. Stiavelli, P. G. van Dokkum, M. Franx, M. Trenti, C. M. Carollo, V. Gonzalez

Abstract: The eXtreme Deep Field (XDF) combines data from ten years of observations with the HST Advanced Camera for Surveys (ACS) and the Wide-Field Camera 3 Infra-Red (WFC3/IR) into the deepest image of the sky ever in the optical/near-IR. Since the initial observations on the Hubble Ultra-Deep Field (HUDF) in 2003, numerous surveys and programs, including supernova followup, HUDF09, CANDELS, and HUDF12 h… ▽ More The eXtreme Deep Field (XDF) combines data from ten years of observations with the HST Advanced Camera for Surveys (ACS) and the Wide-Field Camera 3 Infra-Red (WFC3/IR) into the deepest image of the sky ever in the optical/near-IR. Since the initial observations on the Hubble Ultra-Deep Field (HUDF) in 2003, numerous surveys and programs, including supernova followup, HUDF09, CANDELS, and HUDF12 have contributed additional imaging data across the HUDF region. Yet these have never been combined and made available as one complete ultra-deep optical and near-infrared image dataset. We do so now for the eXtreme Deep Field (XDF) program. Our new and improved processing techniques provide higher quality reductions of the total dataset. All WFC3 near-IR and optical ACS data sets have been fully combined and accurately matched, resulting in the deepest imaging ever taken at these wavelengths ranging from 29.1 to 30.3 AB mag (5sigma in a 0.35" diameter aperture) in 9 filters. The gains in the optical for the four filters done in the original ACS HUDF correspond to a typical improvement of 0.15 mag, with gains of 0.25 mag in the deepest areas. Such gains are equivalent to adding ~130 to ~240 orbits of ACS data to the HUDF. Improved processing alone results in a typical gain of ~0.1 mag. Our 5sigma (optical+near-IR) SExtractor catalogs reveal about 14140 sources in the full field and about 7121 galaxies in the deepest part of the XDF (the HUDF09 region). The XDF is the deepest image of the universe ever taken, reaching, in the combined image for a flat f_nu source, to 31.2 AB mag 5sigma (32.9 at 1sigma) in a 0.35" diameter aperture. △ Less

Submitted 11 September, 2013; v1 submitted 8 May, 2013; originally announced May 2013.

Comments: 13 pages, 12 figures, ApJS in press; changed to match accepted version; XDF imaging data available on MAST Archive at http://archive.stsci.edu/prepds/xdf/, for more information visit also http://xdf.ucolick.org/

Probing the Dawn of Galaxies at z~9-12: New Constraints from HUDF12/XDF and CANDELS Data

Authors: P. A. Oesch, R. J. Bouwens, G. D. Illingworth, I. Labbe, M. Franx, P. G. van Dokkum, M. Trenti, M. Stiavelli, V. Gonzalez, D. Magee

Abstract: We present a comprehensive analysis of z>8 galaxies based on ultra-deep WFC3/IR data. We constrain the evolution of the UV luminosity function (LF) and luminosity densities from z~11 to z~8 by exploiting all the WFC3/IR data over the Hubble Ultra-Deep Field from the HUDF09 and the new HUDF12 program, in addition to the HUDF09 parallel field data, as well as wider area WFC3/IR imaging over GOODS-So… ▽ More We present a comprehensive analysis of z>8 galaxies based on ultra-deep WFC3/IR data. We constrain the evolution of the UV luminosity function (LF) and luminosity densities from z~11 to z~8 by exploiting all the WFC3/IR data over the Hubble Ultra-Deep Field from the HUDF09 and the new HUDF12 program, in addition to the HUDF09 parallel field data, as well as wider area WFC3/IR imaging over GOODS-South. Galaxies are selected based on the Lyman Break Technique in three samples centered around z~9, z~10 and z~11, with seven z~9 galaxy candidates, and one each at z~10 and z~11. We confirm a new z~10 candidate (with z=9.8+-0.6) that was not convincingly identified in our first z~10 sample. The deeper data over the HUDF confirms all our previous z>~7.5 candidates as genuine high-redshift candidates, and extends our samples to higher redshift and fainter limits (H_160~29.8 mag). We perform one of the first estimates of the z~9 UV LF and improve our previous constraints at z~10. Extrapolating the lower redshift UV LF evolution should have revealed 17 z~9 and 9 z~10 sources, i.e., a factor ~3x and 9x larger than observed. The inferred star-formation rate density (SFRD) in galaxies above 0.7 M_sun/yr decreases by 0.6+-0.2 dex from z~8 to z~9, in good agreement with previous estimates. The low number of sources found at z>8 is consistent with a very rapid build-up of galaxies across z~10 to z~8. From a combination of all current measurements, we find a best estimate of a factor 10x decrease in the SFRD from z~8 to z~10, following (1+z)^(-11.4+-3.1). Our measurements thus confirm our previous finding of an accelerated evolution beyond z~8, and signify a rapid build-up of galaxies with M_UV<-17.7 within only ~200 Myr from z~10 to z~8, in the heart of cosmic reionization. △ Less

Submitted 25 January, 2013; originally announced January 2013.

Comments: 21 pages, 13 figures, 6 tables; submitted to ApJ

A Tentative Detection of an Emission Line at 1.6 μm for the z~12 Candidate UDFj-39546284

Authors: Gabriel B. Brammer, Pieter G. van Dokkum, Garth D. Illingworth, Rychard J. Bouwens, Ivo Labbe, Marijn Franx, Ivelina Momcheva, Pascal A. Oesch

Abstract: We present deep WFC3 grism observations of the candidate z~12 galaxy UDFj-39546284 in the HST Ultra Deep Field (UDF), by combining spectroscopic data from the 3D-HST and CANDELS surveys. The total exposure time is 40.5 ks and the spectrum covers 1.10 < λ < 1.65μm. We search for faint emission lines by cross-correlating the 2D G141 spectrum with the observed H160 morphology, a technique that is uni… ▽ More We present deep WFC3 grism observations of the candidate z~12 galaxy UDFj-39546284 in the HST Ultra Deep Field (UDF), by combining spectroscopic data from the 3D-HST and CANDELS surveys. The total exposure time is 40.5 ks and the spectrum covers 1.10 < λ < 1.65μm. We search for faint emission lines by cross-correlating the 2D G141 spectrum with the observed H160 morphology, a technique that is unique to slitless spectroscopy at HST resolution. We find a 2.7σ detection of an emission line at 1.599 μm---just redward of the JH140 filter---with flux 3.5 +/- 1.3x10-18 erg/s/cm2. Assuming the line is real, it contributes 110 +/- 40% of the observed H160 flux and has an observed equivalent width > 7300 Å. If the line is confirmed, it could be Ly-α at z = 12.12. However, a more plausible interpretation, given current results, could be a lower redshift feature such as [O III]λ4959,5007 at z = 2.19. We find two other 3D-HST [O III] emitters within 1000 km/s of that redshift in the GOODS-South field. Additional support for this interpretation comes from the discovery of a bright "[O III] blob" with a secure G141 grism redshift of z = 1.605. This object has a strikingly large observed equivalent width of nearly 9000 Å that results in similar "dropout" colors as UDFj-39546284. △ Less

Submitted 23 January, 2013; v1 submitted 2 January, 2013; originally announced January 2013.

Comments: Accepted for publication in the Astrophysical Journal Letters

Stellar Kinematics of z~2 Galaxies and the Inside-Out Growth of Quiescent Galaxies

Authors: Jesse van de Sande, Mariska Kriek, Marijn Franx, Pieter G. van Dokkum, Rachel Bezanson, Rychard J. Bouwens, Ryan F. Quadri, Hans-Walter Rix, Rosalind E. Skelton

Abstract: Using stellar kinematics measurements, we investigate the growth of massive, quiescent galaxies from z~2 to today. We present X-Shooter spectra from the UV to NIR and dynamical mass measurements of 5 quiescent massive (>10^11 Msun) galaxies at z~2. This triples the sample of z>1.5 galaxies with well constrained (dsigma <100 km/s) velocity dispersion measurements. From spectral population synthesis… ▽ More Using stellar kinematics measurements, we investigate the growth of massive, quiescent galaxies from z~2 to today. We present X-Shooter spectra from the UV to NIR and dynamical mass measurements of 5 quiescent massive (>10^11 Msun) galaxies at z~2. This triples the sample of z>1.5 galaxies with well constrained (dsigma <100 km/s) velocity dispersion measurements. From spectral population synthesis modeling we find that these galaxies have stellar ages that range from 0.5-2 Gyr, with no signs of ongoing star formation. We measure velocity dispersions (290-450 km/s) from absorption lines and find that they are 1.6-2.1 times higher than those of galaxies in the SDSS at the same mass. Sizes are measured using GALFIT from HST-WFC3 H160 and UDS K-band images. The dynamical masses correspond well to the SED based stellar masses, with dynamical masses that are ~15% higher. We find that M_*/M_dyn may decrease slightly with time, which could reflect the increase of the dark matter fraction within an increasing effective radius. We combine different stellar kinematic studies from the literature, and examine the structural evolution from z~2 to z~0: we confirm that at fixed dynamical mass, the effective radius increases by a factor of ~2.8, and the velocity dispersion decreases by a factor of ~1.7. The mass density within one effective radius decreases by a factor of ~20, while within a fixed physical radius (1 kpc) it decreases only mildly (~2). When we allow for an evolving mass limit by selecting a population of galaxies at fixed number density, a stronger size growth with time is found (~4), velocity dispersion decreases by a factor of ~1.4, and interestingly, the mass density within 1 kpc is consistent with no evolution. This finding suggests that massive quiescent galaxies at ~2 grow inside-out, consistent with the expectations from minor mergers. △ Less

Submitted 6 June, 2013; v1 submitted 14 November, 2012; originally announced November 2012.

Comments: 25 pages, 20 figures. Accepted for publication in the Astrophysical Journal

Photometric Constraints on the Redshift of z~10 candidate UDFj-39546284 from deeper WFC3/IR+ACS+IRAC observations over the HUDF

Authors: R. J. Bouwens, P. A. Oesch, G. D. Illingworth, I. Labbe, P. G. van Dokkum, G. Brammer, D. Magee, L. Spitler, M. Franx, R. Smit, M. Trenti, V. Gonzalez, C. M. Carollo

Abstract: Ultra-deep WFC3/IR observations on the HUDF from the HUDF09 program revealed just one plausible z~10 candidate UDFj-39546284. UDFj-39546284 had all the properties expected of a galaxy at z~10 showing (1) no detection in the deep ACS+WFC3 imaging data blueward of the F160W band, exhibiting (2) a blue spectral slope redward of the break, and showing (3) no prominent detection in deep IRAC observatio… ▽ More Ultra-deep WFC3/IR observations on the HUDF from the HUDF09 program revealed just one plausible z~10 candidate UDFj-39546284. UDFj-39546284 had all the properties expected of a galaxy at z~10 showing (1) no detection in the deep ACS+WFC3 imaging data blueward of the F160W band, exhibiting (2) a blue spectral slope redward of the break, and showing (3) no prominent detection in deep IRAC observations. The new, similarly deep WFC3/IR HUDF12 F160W observations over the HUDF09/XDF allow us to further assess this candidate. These observations show that this candidate, previously only detected at ~5.9 sigma in a single band, clearly corresponds to a real source. It is detected at ~5.3 sigma in the new H-band data and at ~7.8 sigma in the full 85-orbit H-band stack. Interestingly, the non-detection of the source (<1 sigma) in the new F140W observations suggests a higher redshift. Formally, the best-fit redshift of the source utilizing all the WFC3+ACS (and IRAC+K-band) observations is 11.8+/-0.3. However, we consider the z~12 interpretation somewhat unlikely, since the source would either need to be ~20x more luminous than expected or show very high-EW Ly-alpha emission (which seems improbable given the extensive neutral gas prevalent early in the reionization epoch). Lower-redshift solutions fail if only continuum models are allowed. Plausible lower-redshift solutions require that the H-band flux be dominated by line emission such as Halpha or [OIII] with extreme EWs. The tentative detection of line emission at 1.6 microns in UDFj-39546284 in a companion paper suggests that such emission may have already been found. △ Less

Submitted 14 February, 2013; v1 submitted 13 November, 2012; originally announced November 2012.

Comments: 6 pages, 4 figures, 1 table, accepted for publication in ApJ Letters, updated to match the version in press

A Rest-frame Optical View on z~4 Galaxies I: Color and Age Distributions from Deep IRAC Photometry of the IUDF10 and GOODS Surveys

Authors: P. A. Oesch, I. Labbe, R. J. Bouwens, G. D. Illingworth, V. Gonzalez, M. Franx, M. Trenti, B. P. Holden, P. G. van Dokkum, D. Magee

Abstract: We present a study of rest-frame UV-to-optical color distributions for z~4 galaxies based on the combination of deep HST/ACS+WFC3/IR data with Spitzer/IRAC imaging. In particular, we use new, ultra-deep data from the IRAC Ultradeep Field program (IUDF10). Our sample contains a total of ~2600 galaxies selected as B-dropout Lyman Break Galaxies (LBGs) in the HUDF and one of its deep parallel fields,… ▽ More We present a study of rest-frame UV-to-optical color distributions for z~4 galaxies based on the combination of deep HST/ACS+WFC3/IR data with Spitzer/IRAC imaging. In particular, we use new, ultra-deep data from the IRAC Ultradeep Field program (IUDF10). Our sample contains a total of ~2600 galaxies selected as B-dropout Lyman Break Galaxies (LBGs) in the HUDF and one of its deep parallel fields, the HUDF09-2, as well as GOODS-North and South. This sample is used to investigate the UV continuum slopes beta and Balmer break colors (J_125-[4.5]) as a function of rest-frame optical luminosity. The [4.5] filter is chosen to avoid potential contamination by strong rest-frame optical emission lines. We find that galaxies at M_z<-21.5 (roughly corresponding to L*[z~4]) are significantly redder than their lower luminosity counterparts. The UV continuum slopes and the J_125-[4.5] colors are well correlated. The most simple explanation for this correlation is that the dust reddening at these redshifts is better described by an SMC-like extinction curve, rather than the typically assumed Calzetti reddening. After correcting for dust, we find that the galaxy population shows mean stellar population ages in the range 10^8.5 to 10^9 yr, with a dispersion of ~0.5 dex, and only weak trends as a function of luminosity. In contrast to some results from the literature, we find that only a small fraction of galaxies shows Balmer break colors which are consistent with extremely young ages, younger than 100 Myr. Under the assumption of smooth star-formation histories, this fraction is only 12-19% for galaxies at M_z<-19.75. Our results are consistent with a gradual build-up of stars and dust in galaxies at z>4, with only a small fraction of stars being formed in short, intense bursts of star-formation. △ Less

Submitted 5 November, 2012; originally announced November 2012.

Comments: 11 pages, 10 figures; submitted to ApJ

The X-ray Star Formation Story as Told by Lyman Break Galaxies in the 4 Ms CDF-S

Authors: Antara R. Basu-Zych, Bret D. Lehmer, Ann E. Hornschemeier, Rychard J. Bouwens, Tassos Fragos, Pascal A. Oesch, Krzysztof Belczynski, W. N. Brandt, Vassiliki Kalogera, Bin Luo, Neal Miller, James R. Mullaney, Panayiotis Tzanavaris, Yongquan Xue, Andreas Zezas

Abstract: We present results from deep X-ray stacking of >4000 high redshift galaxies from z~1 to 8 using the 4 Ms Chandra Deep Field South (CDF-S) data, the deepest X-ray survey of the extragalactic sky to date. The galaxy samples were selected using the Lyman break technique based primarily on recent HST ACS and WFC3 observations. Based on such high specific star formation rates (sSFRs): log SFR/M* > -8.7… ▽ More We present results from deep X-ray stacking of >4000 high redshift galaxies from z~1 to 8 using the 4 Ms Chandra Deep Field South (CDF-S) data, the deepest X-ray survey of the extragalactic sky to date. The galaxy samples were selected using the Lyman break technique based primarily on recent HST ACS and WFC3 observations. Based on such high specific star formation rates (sSFRs): log SFR/M* > -8.7, we expect that the observed properties of these LBGs are dominated by young stellar populations. The X-ray emission in LBGs, eliminating individually detected X-ray sources (potential AGN), is expected to be powered by X-ray binaries and hot gas. We find, for the first time, evidence of evolution in the X-ray/SFR relation. Based on X-ray stacking analyses for z<4 LBGs (covering ~90% of the Universe's history), we find that the 2-10 keV X-ray luminosity evolves weakly with redshift (z) and SFR as log LX = 0.93 log (1+z) + 0.65 log SFR + 39.80. By comparing our observations with sophisticated X-ray binary population synthesis models, we interpret that the redshift evolution of LX/SFR is driven by metallicity evolution in HMXBs, likely the dominant population in these high sSFR galaxies. We also compare these models with our observations of X-ray luminosity density (total 2-10 keV luminosity per Mpc^3) and find excellent agreement. While there are no significant stacked detections at z>5, we use our upper limits from 5<z<8 LBGs to constrain the SMBH accretion history of the Universe around the epoch of reionization. △ Less

Submitted 11 October, 2012; originally announced October 2012.

Comments: 16 pages, 9 figures. Accepted for publication in ApJ

The Spectral Energy Distributions Of z~8 Galaxies From The IRAC Ultra Deep Fields: Emission Lines, Stellar Masses, And Specific Star Formation Rates At 650 Myr

Authors: I. Labbe, P. A. Oesch, R. J. Bouwens, G. D. Illingworth, D. Magee, V. Gonzalez, C. M. Carollo, M. Franx, M. Trenti, P. G. van Dokkum, M. Stiavelli

Abstract: Using new ultradeep Spitzer/IRAC photometry from the IRAC Ultradeep Field program (IUDF), we investigate the stellar populations of a sample of 63 Y-dropout galaxy candidates at z~8, only 650Myr after the Big Bang. The sources are selected from HST/ACS+WFC3/IR data over the Hubble Ultra Deep Field (HUDF), two HUDF parallel fields, and wide area data over the CANDELS/GOODS-South. The new Spitzer/IR… ▽ More Using new ultradeep Spitzer/IRAC photometry from the IRAC Ultradeep Field program (IUDF), we investigate the stellar populations of a sample of 63 Y-dropout galaxy candidates at z~8, only 650Myr after the Big Bang. The sources are selected from HST/ACS+WFC3/IR data over the Hubble Ultra Deep Field (HUDF), two HUDF parallel fields, and wide area data over the CANDELS/GOODS-South. The new Spitzer/IRAC data increase the coverage at 3.6 micron and 4.5 micron to ~120h over the HUDF reaching depths of ~28 (AB,1 sigma). The improved depth and inclusion of brighter candidates result in direct >3 sigma IRAC detections of 20/63 sources, of which 11/63 are detected at > 5 sigma. The average [3.6]-[4.5] colors of IRAC detected galaxies at z~8 are markedly redder than those at z~7, observed only 130Myr later. The simplest explanation is that we witness strong rest-frame optical emission lines (in particular [OIII]4959,5007+Hbeta) moving through the IRAC bandpasses with redshift. Assuming that the average rest-frame spectrum is the same at both z~7 and z~8 we estimate a rest-frame equivalent width of W([OIII]4959,5007+Hbeta) = 670 (+260,-170) Angstrom contributing 0.56 (+0.16,-0.11) mag to the 4.5 micron filter at z~8. The corresponding W(Halpha) = 430 (+160,-110) Angstrom implies an average specific star formation rate of sSFR = 11 (+11,-5) Gyr^-1 and a stellar population age of 100 (+100,-50) Myr. Correcting the spectral energy distribution for the contribution of emission lines lowers the average best-fit stellar masses and mass-to-light ratios by x3, decreasing the integrated stellar mass density to rho*(z=8,MUV<-18)=0.6 (+0.4,-0.3) x 10^6 Msun Mpc^-3. △ Less

Submitted 20 September, 2013; v1 submitted 13 September, 2012; originally announced September 2012.

Comments: 6 pages, 5 figures, accepted for publication in ApJL

The Evolution of Mass-size Relation for Lyman Break Galaxies From z=1 to z=7

Authors: Moein Mosleh, Rik J. Williams, Marijn Franx, Valentino Gonzalez, Rychard J. Bouwens, Pascal Oesch, Ivo Labbe, Garth D. Illingworth, Michele Trenti

Abstract: For the first time, we study the evolution of the stellar mass-size relation for star-forming galaxies from z ~ 4 to z ~ 7 from Hubble-WFC3/IR camera observations of the HUDF and Early Release Science (ERS) field. The sizes are measured by determining the best fit model to galaxy images in the rest-frame 2100 Å\ with the stellar masses estimated from SED fitting to rest-frame optical (from Spitzer… ▽ More For the first time, we study the evolution of the stellar mass-size relation for star-forming galaxies from z ~ 4 to z ~ 7 from Hubble-WFC3/IR camera observations of the HUDF and Early Release Science (ERS) field. The sizes are measured by determining the best fit model to galaxy images in the rest-frame 2100 Å\ with the stellar masses estimated from SED fitting to rest-frame optical (from Spitzer/IRAC) and UV fluxes. We show that the stellar mass-size relation of Lyman-break galaxies (LBGs) persists, at least to z ~ 5, and the median size of LBGs at a given stellar mass increases towards lower redshifts. For galaxies with stellar masses of 9.5<Log(M*/Msun)<10.4 sizes evolve as $(1+z)^{-1.20\pm0.11}$. This evolution is very similar for galaxies with lower stellar masses of 8.6<Log(M*/Msun)<9.5 which is $r_{e} \propto (1+z)^{-1.18\pm0.10}$, in agreement with simple theoretical galaxy formation models at high z. Our results are consistent with previous measurements of the LBGs mass-size relation at lower redshifts (z ~ 1-3). △ Less

Submitted 27 July, 2012; originally announced July 2012.

Comments: Accepted for publication in The Astrophysical Journal Letters

The Brightest of Reionizing Galaxies Survey: Constraints on the Bright End of the z~8 Luminosity Function

Authors: L. D. Bradley, M. Trenti, P. A. Oesch, M. Stiavelli, T. Treu, R. J. Bouwens, J. M. Shull, B. W. Holwerda, N. Pirzkal

Abstract: We report the discovery of 33 Lyman break galaxy (LBG) candidates at z~8 detected in Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging as part of the Brightest of Reionizing Galaxies (BoRG) pure-parallel survey. The ongoing BoRG survey currently has the largest area (274 arcmin^2) with Y_098 (or Y_105), J_125, and H_160 band coverage needed to search for z~8 galaxies, about three tim… ▽ More We report the discovery of 33 Lyman break galaxy (LBG) candidates at z~8 detected in Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging as part of the Brightest of Reionizing Galaxies (BoRG) pure-parallel survey. The ongoing BoRG survey currently has the largest area (274 arcmin^2) with Y_098 (or Y_105), J_125, and H_160 band coverage needed to search for z~8 galaxies, about three times the current CANDELS area, and slightly larger than what will be the final CANDELS wide component with Y_105 data. Our sample of 33 relatively bright Y_098-dropout galaxies have J_125 band magnitudes between 25.5 and 27.4 mag. This is the largest sample of bright (J_125 <~ 27.4) z~8 galaxy candidates presented to date. Combining our dataset with the Hubble Ultra-Deep Field (HUDF09) dataset, we constrain the rest-frame ultraviolet galaxy luminosity function at z~8 over the widest dynamic range currently available. The combined datasets are well fitted by a Schechter function, i.e. φ(L) = φ_* (L/L_*)^α\ e^{-(L/L_*)}, without evidence for an excess of sources at the bright end. At 68% confidence, we derive φ_* = (4.3^{+3.5}_{-2.1}) \times 10^{-4} Mpc^{-3}, M_* = -20.26^{+0.29}_{-0.34}, and a very steep faint-end slope α= -1.98^{+0.23}_{-0.22}. While the best-fit parameters still have a strong degeneracy, especially between φ_* and M_*, our improved coverage at the bright end has reduced the uncertainty of the faint-end power-law slope at z~8 compared to the best previous determination at +/-0.4. With a future expansion of the BoRG survey, combined with planned ultradeep WFC3/IR observations, it will be possible to further reduce this uncertainty and clearly demonstrate the steepening of the faint-end slope compared to measurements at lower redshift, thereby confirming the key role played by small galaxies in the reionization of the universe. △ Less

Submitted 9 October, 2012; v1 submitted 16 April, 2012; originally announced April 2012.

Comments: 13 pages, 9 figures, 9 tables, accepted for publication in the Astrophysical Journal

The Star Formation Rate Function for Redshift z~4-7 Galaxies: Evidence for a Uniform Build-Up of Star-Forming Galaxies During the First 3 Gyr of Cosmic Time

Authors: Renske Smit, Rychard J. Bouwens, Marijn Franx, Garth D. Illingworth, Ivo Labbé, Pacal A. Oesch, Pieter G. van Dokkum

Abstract: We combine recent estimates of dust extinction at z~4-7 with UV luminosity function (LF) determinations to derive star formation rate (SFR) functions at z~4, 5, 6 and 7. SFR functions provide a more physical description of galaxy build-up at high redshift and allow for direct comparisons to SFRs at lower redshifts determined by a variety of techniques. Our SFR functions are derived from well-estab… ▽ More We combine recent estimates of dust extinction at z~4-7 with UV luminosity function (LF) determinations to derive star formation rate (SFR) functions at z~4, 5, 6 and 7. SFR functions provide a more physical description of galaxy build-up at high redshift and allow for direct comparisons to SFRs at lower redshifts determined by a variety of techniques. Our SFR functions are derived from well-established z~4-7 UV LFs, UV-continuum slope trends with redshift and luminosity, and IRX-beta relations. They are well-described by Schechter relations. We extend the comparison baseline for SFR functions to z~2 by considering recent determinations of the Hα and mid-IR luminosity functions. The low-end slopes of the SFR functions are flatter than for the UV LFs, Δα\sim+0.13, and show no clear evolution with cosmic time (z~0-7). In addition, we find that the characteristic value SFR* from the Schechter fit to SFR function exhibits consistent, and substantial, linear growth as a function of redshift from ~5 M_sun/yr at z~8, 650 Myr after the Big Bang, to ~100 M_sun/yr at z~2, ~2.5 Gyr later. Recent results at z~10, close to the onset of galaxy formation, are consistent with this trend. The uniformity of this evolution is even greater than seen in the UV LF over the redshift range z~2-8, providing validation for our dust corrections. These results provide strong evidence that galaxies build up uniformly over the first 3 Gyr of cosmic time. △ Less

Submitted 10 July, 2012; v1 submitted 16 April, 2012; originally announced April 2012.

Comments: Added an appendix, 1 figure and 3 tables: 9 pages, 5 figures, 4 tables, ApJ, in press

The Bright End of the UV Luminosity Function at z~8: New Constraints from CANDELS Data in GOODS-South

Authors: P. A. Oesch, R. J. Bouwens, G. D. Illingworth, V. Gonzalez, M. Trenti, P. G. van Dokkum, M. Franx, I. Labbe, C. M. Carollo, D. Magee

Abstract: We present new z~8 galaxy candidates from a search over ~95 arcmin^2 of WFC3/IR data, tripling the previous search area for bright z~8 galaxies. Our analysis uses newly acquired WFC3/IR imaging data from the CANDELS Multi-Cycle Treasury program over the GOODS South field. These new data are combined with existing deep optical ACS imaging to search for relatively bright (M_UV < -19.5 mag) z~8 galax… ▽ More We present new z~8 galaxy candidates from a search over ~95 arcmin^2 of WFC3/IR data, tripling the previous search area for bright z~8 galaxies. Our analysis uses newly acquired WFC3/IR imaging data from the CANDELS Multi-Cycle Treasury program over the GOODS South field. These new data are combined with existing deep optical ACS imaging to search for relatively bright (M_UV < -19.5 mag) z~8 galaxy candidates using the Lyman Break technique. These new candidates are used to determine the bright end of the UV luminosity function (LF) of star-forming galaxies at z~8. To minimize contamination from lower redshift galaxies, we make full use of all optical ACS data and impose strict non-detection criteria based on an optical chi^2_opt flux measurement. In the whole search area we identify 16 candidate z~8 galaxies, spanning a magnitude range H_160 = 25.7-27.9 mag. The new data show that the UV LF is a factor ~1.7x lower at M_UV < -19.5 mag than determined from the HUDF09 and ERS data alone. Combining this new sample with the previous candidates from the HUDF09 and ERS data allows us to perform the most accurate measurement of the z~8 UV LF yet. Schechter function fits to the combined data result in a best-fit characteristic magnitude of M*(z=8) = -20.04+-0.46 mag. The faint-end slope is very steep, though quite uncertain, with alpha = -2.06+-0.32. A combination of wide area data with additional ultra-deep imaging will be required to significantly reduce the uncertainties on these parameters in the future. △ Less

Submitted 2 November, 2012; v1 submitted 3 January, 2012; originally announced January 2012.

Comments: 16 pages, 12 figures; ApJ in press

Overdensities of Y-dropout Galaxies from the Brightest-of-Reionizing Galaxies Survey: A Candidate Protocluster at Redshift z~8

Authors: M. Trenti, L. D. Bradley, M. Stiavelli, J. M. Shull, P. Oesch, R. J. Bouwens, J. A. Munoz, E. Romano-Diaz, T. Treu, I. Shlosman, C. M. Carollo

Abstract: Theoretical and numerical modeling of dark-matter halo assembly predicts that the most luminous galaxies at high redshift are surrounded by overdensities of fainter companions. We test this prediction with HST observations acquired by our Brightest of Reionizing Galaxies (BoRG) survey, which identified four very bright z~8 candidates as Y-dropout sources in four of the 23 non-contiguous WFC3 field… ▽ More Theoretical and numerical modeling of dark-matter halo assembly predicts that the most luminous galaxies at high redshift are surrounded by overdensities of fainter companions. We test this prediction with HST observations acquired by our Brightest of Reionizing Galaxies (BoRG) survey, which identified four very bright z~8 candidates as Y-dropout sources in four of the 23 non-contiguous WFC3 fields observed. We extend here the search for Y-dropouts to fainter luminosities (M_* galaxies with M_AB\sim-20), with detections at >5sigma confidence (compared to >8sigma confidence adopted earlier) identifying 17 new candidates. We demonstrate that there is a correlation between number counts of faint and bright Y-dropouts at >99.84% confidence. Field BoRG58, which contains the best bright z\sim8 candidate (M_AB=-21.3), has the most significant overdensity of faint Y-dropouts. Four new sources are located within 70arcsec (corresponding to 3.1 comoving Mpc at z=8) from the previously known brighter z\sim8 candidate. The overdensity of Y-dropouts in this field has a physical origin to high confidence (p>99.975%), independent of completeness and contamination rate of the Y-dropout selection. We modeled the overdensity by means of cosmological simulations and estimate that the principal dark matter halo has mass M_h\sim(4-7)x10^11Msun (\sim5sigma density peak) and is surrounded by several M_h\sim10^11Msun halos which could host the fainter dropouts. In this scenario, we predict that all halos will eventually merge into a M_h>2x10^14Msun galaxy cluster by z=0. Follow-up observations with ground and space based telescopes are required to secure the z\sim8 nature of the overdensity, discover new members, and measure their precise redshift. △ Less

Submitted 21 November, 2011; v1 submitted 3 October, 2011; originally announced October 2011.

Comments: Minor revision: ApJ accepted [17 pages (emulateapj style), 7 figures, 2 tables]