Showing 1–50 of 71 results for author: González, V

ANDES, the high-resolution spectrograph for the ELT: RIZ Spectrograph preliminary design

Authors: Bruno Chazelas, Yevgeniy Ivanisenko, Audrey Lanotte, Pablo Santos Diaz, Ludovic Genolet, Michael Sordet, Ian Hughes, Christophe Lovis, Tobias M. Schmidt, Manuel Amate, José Peñate Castro, Afrodisio Vega Moreno, Fabio Tenegi, Roberto Simoes, Jonay I. González Hernández, María Rosa Zapatero Osorio, Javier Piqueras, Tomás Belenguer Dávila, Rocío Calvo Ortega, Roberto Varas González, Luis Miguel González Fernández, Pedro J. Amado, Jonathan Kern, Frank Dionies, Svend-Marian Bauer , et al. (22 additional authors not shown)

Abstract: We present here the preliminary design of the RIZ module, one of the visible spectrographs of the ANDES instrument 1. It is a fiber-fed high-resolution, high-stability spectrograph. Its design follows the guidelines of successful predecessors such as HARPS and ESPRESSO. In this paper we present the status of the spectrograph at the preliminary design stage. The spectrograph will be a warm, vacuum-… ▽ More We present here the preliminary design of the RIZ module, one of the visible spectrographs of the ANDES instrument 1. It is a fiber-fed high-resolution, high-stability spectrograph. Its design follows the guidelines of successful predecessors such as HARPS and ESPRESSO. In this paper we present the status of the spectrograph at the preliminary design stage. The spectrograph will be a warm, vacuum-operated, thermally controlled and fiber-fed echelle spectrograph. Following the phase A design, the huge etendue of the telescope will be reformed in the instrument with a long slit made of smaller fibers. We discuss the system design of the spectrographs system. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: Paper submitted to the SPIE astronomical telescope and instrumentation 2024, conference title : Ground-based and Airborne Instrumentation for Astronomy X, paper reference number : 13096-171

Understanding the Universal Dust Attenuation Scaling Relation of Star-Forming Galaxies

Authors: J. Qin, X. Z. Zheng, S. Wuyts, Z. Lv, M. Qiao, J. -S. Huang, F. S. Liu, A. Katsianis, V. Gonzalez, F. Bian, H. Xu, Z. Pan, W. Liu, Q. -H. Tan, F. X. An, D. D. Shi, Y. Zhang, R. Wen, S. Liu, C. Yang

Abstract: Star-forming galaxies (SFGs) adhere to a surprisingly tight scaling relation of dust attenuation parameterized by the infrared excess (IRX=$L_{\rm IR}/L_{\rm UV}$), being jointly determined by the star formation rate (SFR), galaxy size ($R_{\rm e}$), metallicity ($Z$/Z$_\odot$) and axial ratio ($b/a$). We examine how these galaxy parameters determine the effective dust attenuation and give rise to… ▽ More Star-forming galaxies (SFGs) adhere to a surprisingly tight scaling relation of dust attenuation parameterized by the infrared excess (IRX=$L_{\rm IR}/L_{\rm UV}$), being jointly determined by the star formation rate (SFR), galaxy size ($R_{\rm e}$), metallicity ($Z$/Z$_\odot$) and axial ratio ($b/a$). We examine how these galaxy parameters determine the effective dust attenuation and give rise to the universal IRX relation, utilizing a simple two-component star-dust geometry model in which dust in the dense and diffuse interstellar medium (ISM) follows exponential mass density profiles, connected with but not necessarily identical to the stellar mass profiles. Meanwhile, empirical relations are adopted to link galaxy properties, including the gas--star formation relation, the dust-to-stellar size relation, as well as the dust-to-gas ratio versus metallicity relation. By fitting a large sample of local SFGs with the model, we obtain the best-fitting model parameters as a function of metallicity, showing that the two-component geometry model is able to successfully reproduce the dependence of IRX on SFR, $R_{\rm e}$, $b/a$ at given $Z$/Z$_\odot$, as well as the dependence of power-law indices on metallicity. Moreover, we also retrieve constraints on the model geometry parameters, including the optical depth of birth clouds (BCs), BC-to-total dust mass fraction, BC covering factor of UV-emitting stars, and star-to-total dust disc radius ratio, which all evolve with galaxy metallicity. Finally, a consistent picture of how the star-dust geometry in SFGs evolves with galaxy metallicity is discussed. △ Less

Submitted 30 January, 2024; v1 submitted 27 December, 2023; originally announced December 2023.

Comments: 20 pages, 10 figures, published in MNRAS (2024, Volume 528, Issue 1, pp.658-675); A PHTHON package IRX_TAU_TOT is available at https://github.com/LvZF/irx_tau_tot/ to calculate the total dust optical depth of a galaxy with given metallicity and best-fitting geometry parameters

Journal ref: MNRAS, 528, 658 (2024)

Unveiling the hidden universe with JWST: The contribution of dust-obscured galaxies to the stellar mass function at $z\sim3-8$

Authors: R. Gottumukkala, L. Barrufet, P. A. Oesch, A. Weibel, N. Allen, B. Alcalde Pampliega, E. J. Nelson, C. C. Williams, G. Brammer, Y. Fudamoto, V. González, K. E. Heintz, G. Illingworth, D. Magee, R. P. Naidu, M. Shuntov, M. Stefanon, S. Toft, F. Valentino, M. Xiao

Abstract: With the advent of JWST, we can probe the rest-frame optical emission of galaxies at $z>3$ with high sensitivity and spatial resolution, making it possible to accurately characterise red, optically-faint galaxies and thus move towards a more complete census of the galaxy population at high redshifts. To this end, we present a sample of 148 massive, dusty galaxies from the JWST/CEERS survey, colour… ▽ More With the advent of JWST, we can probe the rest-frame optical emission of galaxies at $z>3$ with high sensitivity and spatial resolution, making it possible to accurately characterise red, optically-faint galaxies and thus move towards a more complete census of the galaxy population at high redshifts. To this end, we present a sample of 148 massive, dusty galaxies from the JWST/CEERS survey, colour-selected using solely JWST bands. With deep JWST/NIRCam data from 1.15$μ$m to 4.44$μ$m and ancillary HST/ACS and WFC3 data, we determine the physical properties of our sample using spectral energy distribution fitting with BAGPIPES. We demonstrate that our selection method efficiently identifies massive ($\mathrm{\langle \log M_\star/M_\odot \rangle \sim 10}$) and dusty ($\mathrm{\langle A_V\rangle \sim 2.7\ mag}$) sources, with a majority at $z>3$ and predominantly lying on the galaxy main-sequence. The main results of this work are the stellar mass functions (SMF) of red, optically-faint galaxies from redshifts between $3<z<8$: these galaxies make up a significant relative fraction of the pre-JWST total SMF at $3<z<4$ and $4<z<6$, and dominate the high-mass end of the pre-JWST SMF at $6<z<8$, suggesting that our census of the galaxy population needs amendment at these epochs. While larger areas need to be surveyed in the future, our results suggest already that the integrated stellar mass density at $\mathrm{\log M_\star/M_\odot\geq9.25}$ may have been underestimated in pre-JWST studies by up to $\sim$15-20\% at $z\sim3-6$, and up to $\sim$45\% at $z\sim6-8$, indicating the rapid onset of obscured stellar mass assembly in the early universe. △ Less

Submitted 13 June, 2024; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: 19 pages, 10 figures, published in MNRAS; view published article at https://academic.oup.com/mnras/article/530/1/966/7630203

The ALMA Reionization Era Bright Emission Line Survey (REBELS): The molecular gas content of galaxies at z~7

Authors: M. Aravena, K. E. Heintz, M. Dessauges-Zavadsky, P. A. Oesch, H. S. B. Algera, R. J. Bouwens, E. Da Cunha, P. Dayal, I. De Looze, A. Ferrara, Y. Fudamoto, V. Gonzalez, L. Graziani, H. Inami, A. Pallottini, R. Schneider, S. Schouws, L. Sommovigo, M. Topping, P. van der Werf, M. Palla

Abstract: A key to understanding the formation of the first galaxies is to quantify the content of the molecular gas as the fuel for star formation activity through the epoch of reionization. In this paper, we use the 158$μ$m [CII] fine-structure emission line as a tracer of the molecular gas in the interstellar medium (ISM) in a sample of $z=6.5-7.5$ galaxies recently unveiled by the Reionization Era Brigh… ▽ More A key to understanding the formation of the first galaxies is to quantify the content of the molecular gas as the fuel for star formation activity through the epoch of reionization. In this paper, we use the 158$μ$m [CII] fine-structure emission line as a tracer of the molecular gas in the interstellar medium (ISM) in a sample of $z=6.5-7.5$ galaxies recently unveiled by the Reionization Era Bright Line Emission Survey, REBELS, with the Atacama Large Millimeter/submillimeter Array. We find substantial amounts of molecular gas ($\sim10^{10.5}\ M_\odot$) comparable to those found in lower redshift galaxies for similar stellar masses ($\sim10^{10}\ M_\odot$). The REBELS galaxies appear to follow the standard scaling relations of molecular gas to stellar mass ratio ($μ_{\rm mol}$) and gas depletion timescale ($t_{\rm dep}$) with distance to the star-forming main-sequence expected from extrapolations of $z\sim1-4$ observations. We find median values at $z\sim7$ of $μ_{\rm mol}=2.6_{-1.4}^{4.1}$ and $t_{\rm dep}=0.5_{-0.14}^{+0.26}$ Gyr, indicating that the baryonic content of these galaxies is gas-phase dominated and little evolution from $z\sim7$ to 4. Our measurements of the cosmic density of molecular gas, log$(ρ_{\rm mol}/(M_\odot {\rm Mpc}^{-3}))=6.34^{+0.34}_{-0.31}$, indicate a steady increase by an order of magnitude from $z\sim7$ to 4. △ Less

Submitted 29 September, 2023; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 11 pages, 7 figures. Typo fixed in an author's name

The ALMA REBELS Survey: Discovery of a massive, highly star-forming and morphologically complex ULIRG at $z =7.31$

Authors: A. P. S. Hygate, J. A. Hodge, E. da Cunha, M. Rybak, S. Schouws, H. Inami, M. Stefanon, L. Graziani, R. Schneider, P. Dayal, R. J. Bouwens, R. Smit, R. A. A. Bowler, R. Endsley, V. Gonzalez, P. A. Oesch, D. P. Stark, H. S. B. Algera, M. Aravena, L. Barrufet, A. Ferrara, Y. Fudamoto, J. H. A, I. De Looze, T. Nanayakkara , et al. (5 additional authors not shown)

Abstract: We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rmμm$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a sig… ▽ More We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rmμm$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a significant stellar mass of $M_{*}=8^{+4}_{-2}\times10^{9}$ M$_\odot$. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR $=199^{+101}_{-63}$ M$_\odot$ yr$^{-1}$. This is about four times the SFR estimated from an extrapolated main-sequence. We also infer a [CII]-based molecular gas mass of $M_{\rm H_2}=5.1^{+5.1}_{-2.6}\times10^{10}$ $M_\odot$, implying a molecular gas depletion time of $ t_{\rm depl, H_2}=0.3^{+0.3}_{-0.2}$ Gyr. We observe a [CII] velocity gradient consistent with disc rotation, but given the current resolution we cannot rule out a more complex velocity structure such as a merger. The spectrum exhibits excess [CII] emission at large positive velocities ($\sim500$ km s$^{-1}$), which we interpret as either a merging companion or an outflow. In the outflow scenario, we derive a lower limit of the mass outflow rate of 200 M$_\odot$ yr$^{-1}$, which is consistent with expectations for a star formation-driven outflow. Given its large stellar mass, SFR and molecular gas reservoir $\sim700$ Myr after the Big Bang, we explore the future evolution of REBELS-25. Considering a simple, conservative model assuming an exponentially declining star formation history, constant star formation efficiency, and no additional gas inflow, we find that REBELS-25 has the potential to evolve into a galaxy consistent with the properties of high-mass quiescent galaxies recently observed at $z\sim4$. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: Accepted for publication in MNRAS. 21 pages, 8 figures

The Hα Luminosity Function of Galaxies at z {\sim} 4.5

Authors: Victoria Bollo, Valentino González, Mauro Stefanon, Pascal A. Oesch, Rychard J. Bouwens, Renske Smit, Garth D. Illingworth, Ivo Labbé

Abstract: We present the Hα luminosity function (LF) derived from a large sample of Lyman break galaxies at z {\sim} 4.5 over the GOODS-South and North fields. This study makes use of the new, full-depth Spitzer/IRAC [3.6] and [4.5] imaging from the GOODS Re-ionization Era wide-Area Treasury from the Spitzer program. The Hα flux is derived from the offset between the continuum flux estimated from the best-f… ▽ More We present the Hα luminosity function (LF) derived from a large sample of Lyman break galaxies at z {\sim} 4.5 over the GOODS-South and North fields. This study makes use of the new, full-depth Spitzer/IRAC [3.6] and [4.5] imaging from the GOODS Re-ionization Era wide-Area Treasury from the Spitzer program. The Hα flux is derived from the offset between the continuum flux estimated from the best-fit spectral energy distribution, and the observed photometry in IRAC [3.6]. From these measurements, we build the Hα LF and study its evolution providing the best constraints of this property at high redshift, where spectroscopy of Hα is not yet available. Schechter parameterizations of the Hα LF show a decreasing evolution of {Φ^\star} with redshift, increasing evolution in L{^\star}, and no significant evolution in the faint-end slope at high z. We find that star formation rates (SFRs) derived from Hα are higher than those derived from the rest-frame UV for low SFR galaxies but the opposite happens for the highest SFRs. This can be explained by lower mass galaxies (also lower SFR) having, on average, rising star formation histories (SFHs), while at the highest masses the SFHs may be declining. The SFR function is steeper, and because of the excess SFR(Hα) compared to SFR(UV) at low SFRs, the SFR density estimated from Hα is higher than the previous estimates based on UV luminosities. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: Astrophysical Journal, Accepted, 17 pages, 16 figures

Journal ref: The Astrophysical Journal, 946:117 (17pp), 2023 April 1

Performance of the MALTA Telescope

Authors: Milou van Rijnbach, Giuliano Gustavino, Phil Allport, Igancio Asensi, Dumitru Vlad Berlea, Daniela Bortoletto, Craig Buttar, Edoardo Charbon, Florian Dachs, Valerio Dao, Dominik Dobrijevic, Leyre Flores Sanz de Acedo, Andrea Gabrielli, Martin Gazi, Laura Gonella, Vicente Gonzalez, Stefan Guindon, Matt LeBlanc, Heinz Pernegger, Francesco Piro, Petra Riedler, Heidi Sandaker, Abhishek Sharma, Carlos Solans Sanchez, Walter Snoeys , et al. (5 additional authors not shown)

Abstract: MALTA is part of the Depleted Monolithic Active Pixel sensors designed in Tower 180nm CMOS imaging technology. A custom telescope with six MALTA planes has been developed for test beam campaigns at SPS, CERN, with the ability to host several devices under test. The telescope system has a dedicated custom readout, online monitoring integrated into DAQ with realtime hit map, time distribution and ev… ▽ More MALTA is part of the Depleted Monolithic Active Pixel sensors designed in Tower 180nm CMOS imaging technology. A custom telescope with six MALTA planes has been developed for test beam campaigns at SPS, CERN, with the ability to host several devices under test. The telescope system has a dedicated custom readout, online monitoring integrated into DAQ with realtime hit map, time distribution and event hit multiplicity. It hosts a dedicated fully configurable trigger system enabling to trigger on coincidence between telescope planes and timing reference from a scintillator. The excellent time resolution performance allows for fast track reconstruction, due to the possibility to retain a low hit multiplicity per event which reduces the combinatorics. This paper reviews the architecture of the system and its performance during the 2021 and 2022 test beam campaign at the SPS North Area. △ Less

Submitted 21 April, 2023; v1 submitted 3 April, 2023; originally announced April 2023.

The ALMA REBELS Survey: The First Infrared Luminosity Function Measurement at $\mathbf{z \sim 7}

Authors: L. Barrufet, P. A. Oesch, R. Bouwens, H. Inami, L. Sommovigo, H. Algera, E. da Cunha, M. Aravena, P. Dayal, A. Ferrara, Y. Fudamoto, V. Gonzalez, L. Graziani, A. Hygate, I. de Looze, T. Nanayakkara, A. Pallottini, R. Schneider, M. Stefanon, M. Topping, P. van Der Werf

Abstract: We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a UV-selected galaxy sample with ALMA spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at $\mathrm{z=6.4-7.7}$ with [CII] 158$\micron$ line scans. 16 sources exhibit… ▽ More We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a UV-selected galaxy sample with ALMA spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at $\mathrm{z=6.4-7.7}$ with [CII] 158$\micron$ line scans. 16 sources exhibit a dust detection, 15 of which are also spectroscopically confirmed through the [CII] line. The IR luminosities of the sample range from $\log L_{IR}/L_\odot=11.4$ to 12.2. Using the UVLF as a proxy to derive the effective volume for each of our target sources, we derive IRLF estimates, both for detections and for the full sample including IR luminosity upper limits. The resulting IRLFs are well reproduced by a Schechter function with the characteristic luminosity of $\log L_{*}/L_\odot=11.6^{+0.2}_{-0.1}$. Our observational results are in broad agreement with the average of predicted IRLFs from simulations at $z\sim7$. Conversely, our IRLFs lie significantly below lower redshift estimates, suggesting a rapid evolution from $z\sim4$ to $z\sim7$, into the reionization epoch. The inferred obscured contribution to the cosmic star-formation rate density at $z\sim7$ amounts to $\mathrm{log(SFRD/M_{\odot}/yr/Mpc^{3}) = -2.66^{+0.17}_{-0.14} }$ which is at least $\sim$10\% of UV-based estimates. We conclude that the presence of dust is already abundant in the EoR and discuss the possibility of unveiling larger samples of dusty galaxies with future ALMA and JWST observations. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 9 pages, 5 figures

Search for $^{22}$Na in novae supported by a novel method for measuring femtosecond nuclear lifetimes

Authors: C. Fougères, F. de Oliveira Santos, J. José, C. Michelagnoli, E. Clément, Y. H. Kim, A. Lemasson, V. Guimaraes, D. Barrientos, D. Bemmerer, G. Benzoni, A. J. Boston, R. Bottger, F. Boulay, A. Bracco, I. Celikovic, B. Cederwall, M. Ciemala, C. Delafosse, C. Domingo-Pardo, J. Dudouet, J. Eberth, Z. Fulop, V. Gonzalez, J. Goupil , et al. (36 additional authors not shown)

Abstract: Classical novae are thermonuclear explosions in stellar binary systems, and important sources of $^{26}$Al and $^{22}$Na. While gamma rays from the decay of the former radioisotope have been observed throughout the Galaxy, $^{22}$Na remains untraceable. The half-life of $^{22}$Na (2.6 yr) would allow the observation of its 1.275 MeV gamma-ray line from a cosmic source. However, the prediction of s… ▽ More Classical novae are thermonuclear explosions in stellar binary systems, and important sources of $^{26}$Al and $^{22}$Na. While gamma rays from the decay of the former radioisotope have been observed throughout the Galaxy, $^{22}$Na remains untraceable. The half-life of $^{22}$Na (2.6 yr) would allow the observation of its 1.275 MeV gamma-ray line from a cosmic source. However, the prediction of such an observation requires good knowledge of the nuclear reactions involved in the production and destruction of this nucleus. The $^{22}$Na($p,γ$)$^{23}$Mg reaction remains the only source of large uncertainty about the amount of $^{22}$Na ejected. Its rate is dominated by a single resonance on the short-lived state at 7785.0(7) keV in $^{23}$Mg. In the present work, a combined analysis of particle-particle correlations and velocity-difference profiles is proposed to measure femtosecond nuclear lifetimes. The application of this novel method to the study of the $^{23}$Mg states, combining magnetic and highly-segmented tracking gamma-ray spectrometers, places strong limits on the amount of $^{22}$Na produced in novae, explains its non-observation to date in gamma rays (flux < 2.5x$10^{-4}$ ph/(cm$^2$s)), and constrains its detectability with future space-borne observatories. △ Less

Submitted 12 December, 2022; originally announced December 2022.

Comments: 18 pages, 3 figures, 1 table

Journal ref: Nature Communications 14, 2023

The ALMA REBELS Survey: The Dust-obscured Cosmic Star Formation Rate Density at Redshift 7

Authors: Hiddo Algera, Hanae Inami, Pascal Oesch, Laura Sommovigo, Rychard Bouwens, Michael Topping, Sander Schouws, Mauro Stefanon, Daniel Stark, Manuel Aravena, Laia Barrufet, Elisabete da Cunha, Pratika Dayal, Ryan Endsley, Andrea Ferrara, Yoshinobu Fudamoto, Valentino Gonzalez, Luca Graziani, Jacqueline Hodge, Alexander Hygate, Ilse de Looze, Themiya Nanayakkara, Rafaella Schneider, Paul van der Werf

Abstract: Cosmic dust is an essential component shaping both the evolution of galaxies and their observational signatures. How quickly dust builds up in the early Universe remains an open question that requires deep observations at (sub-)millimeter wavelengths to resolve. Here we use Atacama Large Millimeter Array observations of 45 galaxies from the Reionization Era Bright Emission Line Survey (REBELS) and… ▽ More Cosmic dust is an essential component shaping both the evolution of galaxies and their observational signatures. How quickly dust builds up in the early Universe remains an open question that requires deep observations at (sub-)millimeter wavelengths to resolve. Here we use Atacama Large Millimeter Array observations of 45 galaxies from the Reionization Era Bright Emission Line Survey (REBELS) and its pilot programs, designed to target [CII] and dust emission in UV-selected galaxies at $z\sim7$, to investigate the dust content of high-redshift galaxies through a stacking analysis. We find that the typical fraction of obscured star formation $f_\mathrm{obs} = \mathrm{SFR}_\mathrm{IR} / \mathrm{SFR}_\mathrm{UV + IR}$ depends on stellar mass, similar to what is observed at lower redshift, and ranges from $f_\mathrm{obs} \approx 0.3 - 0.6$ for galaxies with $\log_{10}\left(M_\star / M_\odot\right) = 9.4 - 10.4$. We further adopt the $z\sim7$ stellar mass function from the literature to extract the obscured cosmic star formation rate density (SFRD) from the REBELS survey. Our results suggest only a modest decrease in the SFRD between $3\lesssim z \lesssim 7$, with dust-obscured star formation still contributing $\sim30\%$ at $z\sim7$. While we extensively discuss potential caveats, our analysis highlights the continued importance of dust-obscured star formation even well into the epoch of reionization. △ Less

Submitted 9 November, 2022; v1 submitted 17 August, 2022; originally announced August 2022.

Comments: MNRAS Accepted (01 Nov 2022)

Unveiling the Nature of Infrared Bright, Optically Dark Galaxies with Early JWST Data

Authors: L. Barrufet, P. A. Oesch, A. Weibel, G. Brammer, R. Bezanson, R. Bouwens, Y. Fudamoto, V. Gonzalez, R. Gottumukkala, G. Illingworth, K. E. Heintz, B. Holden, I. Labbe, D. Magee, R. P. Naidu, E. Nelson, M. Stefanon, R. Smit, P. van Dokkum, J. Weaver, C. Williams

Abstract: Over the last few years, both ALMA and Spitzer/IRAC observations have revealed a population of likely massive galaxies at $z>3$ that was too faint to be detected in HST rest-frame ultraviolet imaging. However, due to the very limited photometry for individual galaxies, the true nature of these so-called HST-dark galaxies has remained elusive. Here, we present the first sample of such galaxies obse… ▽ More Over the last few years, both ALMA and Spitzer/IRAC observations have revealed a population of likely massive galaxies at $z>3$ that was too faint to be detected in HST rest-frame ultraviolet imaging. However, due to the very limited photometry for individual galaxies, the true nature of these so-called HST-dark galaxies has remained elusive. Here, we present the first sample of such galaxies observed with very deep, high-resolution NIRCam imaging from the Early Release Science Program CEERS. 30 HST-dark sources are selected based on their red colours across 1.6 $μ$m to 4.4 $μ$m. Their physical properties are derived from 12-band multi-wavelength photometry, including ancillary HST imaging. We find that these galaxies are generally heavily dust-obscured ($A_{V}\sim2$ mag), massive ($\log (M/M_{\odot}) \sim10$), star-forming sources at $z\sim2-8$ with an observed surface density of $\sim0.8$ arcmin$^{-2}$. This suggests that an important fraction of massive galaxies may have been missing from our cosmic census at $z>3$ all the way into the Reionization epoch. The HST-dark sources lie on the main sequence of galaxies and add an obscured star formation rate density (SFRD) of $\mathrm{3.2^{+1.8}_{-1.3} \times 10^{-3} M_{\odot}/yr/Mpc^{3}}$ at $z\sim7$ showing likely presence of dust in the Epoch of Reionization. Our analysis shows the unique power of JWST to reveal this previously missing galaxy population and to provide a more complete census of galaxies at $z=2-8$ based on rest-frame optical imaging. △ Less

Submitted 15 March, 2023; v1 submitted 29 July, 2022; originally announced July 2022.

Comments: 9 pages, 7 figures, accepted to MNRAS

Deep Spitzer/IRAC Data for z~10 galaxies Reveal Blue Balmer Break Colors: Young Stellar Populations at ~500 Myr of Cosmic Time

Authors: Mauro Stefanon, Rychard J. Bouwens, Ivo Labbé, Garth D. Illingworth, Valentino Gonzalez, Pascal A. Oesch

Abstract: We present the deepest constraints yet on the median rest-UV+optical SED of $z\sim10$ galaxies, prior to JWST science operations. We constructed stacks based on four robust $J_{125}$-dropouts, previously identified across the GOODS fields. We used archival HST/WFC3 data and the full depth Spitzer/IRAC mosaics from the GREATS program, the deepest coverage at $\sim3-5μ$m to date. The most remarkable… ▽ More We present the deepest constraints yet on the median rest-UV+optical SED of $z\sim10$ galaxies, prior to JWST science operations. We constructed stacks based on four robust $J_{125}$-dropouts, previously identified across the GOODS fields. We used archival HST/WFC3 data and the full depth Spitzer/IRAC mosaics from the GREATS program, the deepest coverage at $\sim3-5μ$m to date. The most remarkable feature of the SED is a blue IRAC $[3.6]-[4.5]=-0.18\pm0.25$ mag color. We also find a nearly flat $H_{160}-[3.6]=0.07\pm0.22$ mag color, corresponding to a UV slope $β= -1.92\pm0.25$. This is consistent with previous studies, and indicative of minimal dust absorption. The observed blue IRAC color and SED fitting suggest that $z\sim10$ galaxies have very young (few $\times10$ Myr) stellar populations, with $80\%$ of stars being formed in the last $\lesssim 160$ Myr ($2σ$). While an exciting result, the uncertainties on the SED are too large to allow us to place strong constraints on the presence of a nebular continuum in $z\sim10$ galaxies (as might be suggested by the blue $[3.6]-[4.5] < 0$ mag color). The resulting sSFR is consistent with the specific accretion rate of dark matter halos, indicative of a star-formation efficiency showing quite limited evolution at such early epochs. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: Submitted to ApJ. 9 pages, 3 figures. Comments welcome

The ALMA REBELS Survey: Average [CII] $158\,\rm{μm}$ sizes of Star-Forming Galaxies from $z\sim 7$ to $z\sim 4$

Authors: Y. Fudamoto, R. Smit, R. A. A. Bowler, P. A. Oesch, R. Bouwens, M. Stefanon, H. Inami, R. Endsley, V. Gonzalez, S. Schouws, D. Stark, H. S. B. Algera, M. Aravena, L. Barrufet, E. da Cunha, P. Dayal, A. Ferrara, L. Graziani, J. A. Hodge, A. P. S. Hygate, A. K. Inoue, T. Nanayakkara, A. Pallottini, E. Pizzati, R. Schneider , et al. (14 additional authors not shown)

Abstract: We present the average [CII] $158\,\rm{μm}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{μm}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective ra… ▽ More We present the average [CII] $158\,\rm{μm}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{μm}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4<z<6$ galaxies observed by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). By analysing [CII] sizes of $4<z<6$ galaxies in two redshift bins, we find an average [CII] size of $r_{\rm e}=2.2\pm0.2\,\rm{kpc}$ and $r_{\rm e}=2.5\pm0.2\,\rm{kpc}$ for $z\sim5.5$ and $z\sim4.5$ galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the [CII] $158\,{\rm μm}$ emitting regions at redshift between $z\sim7$ and $z\sim4$. This finding suggest that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range. △ Less

Submitted 15 June, 2022; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 9pages, 4 figures, submitted to ApJ; minor corrections were made for typos in the author list and the title

The Physical Properties of Star-Forming Galaxies with Strong [O III] Lines at z=3.25

Authors: Run Wen, FangXia An, Xian Zhong Zheng, Dong Dong Shi, Jianbo Qin, Valentino Gonzalez, Fuyan Bian, Haiguang Xu, Zhizheng Pan, Qing-Hua Tan, Wenhao Liu, Min Fang, Jian Ren, Yu Heng Zhang, Man Qiao, Shuang Liu

Abstract: We present an analysis of physical properties of 34 [O III] emission-line galaxies (ELGs) at z=3.254$\pm$0.029 in the Extended Chandra Deep Field South (ECDFS). These ELGs are selected from deep narrow H2S(1) and broad Ks imaging of 383 arcmin$^{2}$ obtained with CFHT/WIRCam. We construct spectral energy distributions (SEDs) from U to Ks to derive the physical properties of ELGs. These [O III] ELG… ▽ More We present an analysis of physical properties of 34 [O III] emission-line galaxies (ELGs) at z=3.254$\pm$0.029 in the Extended Chandra Deep Field South (ECDFS). These ELGs are selected from deep narrow H2S(1) and broad Ks imaging of 383 arcmin$^{2}$ obtained with CFHT/WIRCam. We construct spectral energy distributions (SEDs) from U to Ks to derive the physical properties of ELGs. These [O III] ELGs are identified as starburst galaxies with strong [O III] lines of L([O III]) ~ 10$^{42.6}$ - 10$^{44.2}$ erg s$^{-1}$, and have stellar masses of M* ~ 10$^{9.0}$-10$^{10.6}$ M$_\odot$ and star formation rates of ~ 10-210 M$_\odot$ yr$^{-1}$. Our results show that 24% of our sample galaxies are dusty with Av > 1 mag and EW(OIII)$_{rest}$ ~ 70-500 $Å$, which are often missed in optically selected [O III] ELG samples. Their rest-frame UV and optical morphologies from HST/ACS and HST/WFC3 deep imaging reveal that these [O III] ELGs are mostly multiple-component systems (likely mergers) or compact. And 20% of them are nearly invisible in the rest-frame UV owing to heavy dust attenuation. Interestingly, we find that our samples reside in an overdensity consisting of two components: one southeast (SE) with an overdensity factor of $δ_{gal}$ ~ 41 over a volume of 13$^{3}$ cMpc$^{3}$ and the other northwest (NW) with $δ_{gal}$ ~ 38 over a volume of 10$^{3}$ cMpc$^{3}$. The two overdense substructures are expected to be virialized at z=0 with a total mass of ~ 1.1 x 10$^{15}$ M$_\odot$ and ~ 4.8 x 10$^{14}$ M$_\odot$, and probably merge into a Coma-like galaxy cluster. △ Less

Submitted 22 July, 2022; v1 submitted 26 May, 2022; originally announced May 2022.

Comments: 22 pages, 11 figures, 3 tables. Accepted for publication in ApJ

Journal ref: Astrophysics Journal, year:2022, month:july, volume:933, pages:50

High Equivalent Width of Hα+[N II] Emission in z~8 Lyman-break Galaxies from IRAC 5.8μm Observations: Evidence for Efficient Lyman-continuum Photon production in the Epoch of Re-ionization

Authors: Mauro Stefanon, Rychard J. Bouwens, Garth D. Illingworth, Ivo Labbé, Pascal A. Oesch, Valentino Gonzalez

Abstract: We measure, for the first time, the median equivalent width (EW) of H$α$+[N II] in star-forming galaxies at $z\sim8$. Our estimate leverages the unique photometric depth of the Spitzer/IRAC $5.8μ$m-band mosaics (probing $\approx 5500 - 7100$ A at $z\sim8$) of the GOODS Reionization Era Wide Area Treasury from Spitzer (GREATS) program. We median stacked the stamps of $102$ Lyman-break galaxies in t… ▽ More We measure, for the first time, the median equivalent width (EW) of H$α$+[N II] in star-forming galaxies at $z\sim8$. Our estimate leverages the unique photometric depth of the Spitzer/IRAC $5.8μ$m-band mosaics (probing $\approx 5500 - 7100$ A at $z\sim8$) of the GOODS Reionization Era Wide Area Treasury from Spitzer (GREATS) program. We median stacked the stamps of $102$ Lyman-break galaxies in the $3.6, 4.5, 5.8$ and $8.0μ$m bands, after carefully removing potential contamination from neighbouring sources. We infer an extreme rest-frame EW$_0$(H$α$+[N II])$=2328^{+1326}_{-1127}$ A from the measured red $[3.6]-[5.8]=0.82\pm0.27$ mag, consistent with young ($\lesssim10^7$ yr) average stellar population ages at $z\sim8$. This implies an ionizing photon production efficiency of $\log(ξ_{\mathrm{ion},0}/\mathrm{erg\ Hz}^{-1})=25.97^{+0.18}_{-0.28}$. Such a high value for photo production, similar to the highest values found at $z\lesssim4$, indicates that only modest escape fractions $f_\mathrm{esc}\lesssim0.3$ (at $2σ$) are sufficient for galaxies brighter than $M_\mathrm{UV}<-18$ mag to re-ionize the neutral Hydrogen at $z\sim8$. This requirement is relaxed even more to $f_\mathrm{esc}\le 0.1$ when considering galaxies brighter than $M_\mathrm{UV}\approx -13$ mag, consistent with recent luminosity functions and as typically assumed in studies addressing re-ionization. These exceptional results clearly indicate that galaxies can be the dominant source of reionizing photons, and provide us with an exciting glimpse into what we might soon learn about the early universe, and particularly about the Reionization Epoch, from forthcoming JWST/MIRI and NIRCam programs. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: Submitted to ApJ. 12 pages, 5 figures (main) + appendix. Comments welcome

The ALMA REBELS Survey: Dust Continuum Detections at z > 6.5

Authors: Hanae Inami, Hiddo S. B. Algera, Sander Schouws, Laura Sommovigo, Rychard Bouwens, Renske Smit, Mauro Stefanon, Rebecca A. A. Bowler, Ryan Endsley, Andrea Ferrara, Pascal Oesch, Daniel Stark, Manuel Aravena, Laia Barrufet, Elisabete da Cunha, Pratika Dayal, Ilse De Looze, Yoshinobu Fudamoto, Valentino Gonzalez, Luca Graziani, Jacqueline A. Hodge, Alexander P. S. Hygate, Themiya Nanayakkara, Andrea Pallottini, Dominik A. Riechers , et al. (3 additional authors not shown)

Abstract: We report 18 dust continuum detections ($\geq 3.3σ$) at $\sim88{\rm μm}$ and $158{\rm μm}$ out of 49 ultraviolet(UV)-bright galaxies ($M_{\rm UV} < -21.3$ mag) at $z>6.5$, observed by the Cycle-7 ALMA Large Program, REBELS and its pilot programs. This has more than tripled the number of dust continuum detections known at $z>6.5$. Out of these 18 detections, 12 are reported for the first time as pa… ▽ More We report 18 dust continuum detections ($\geq 3.3σ$) at $\sim88{\rm μm}$ and $158{\rm μm}$ out of 49 ultraviolet(UV)-bright galaxies ($M_{\rm UV} < -21.3$ mag) at $z>6.5$, observed by the Cycle-7 ALMA Large Program, REBELS and its pilot programs. This has more than tripled the number of dust continuum detections known at $z>6.5$. Out of these 18 detections, 12 are reported for the first time as part of REBELS. In addition, 15 of the dust continuum detected galaxies also show a [CII]$_{\rm 158{\rm μm}}$ emission line, providing us with accurate redshifts. We anticipate more line emission detections from six targets (including three continuum detected targets) where observations are still ongoing. The dust continuum detected sources in our sample tend to have a redder UV spectral slope than the ones without a dust continuum detection. We estimate that all of the sources have an infrared (IR) luminosity ($L_{\rm IR}$) in a range of $3-8 \times 10^{11} L_\odot$, except for one with $L_{\rm IR} = 1.5^{+0.8}_{-0.5} \times 10^{12}\,L_{\odot}$. Their fraction of obscured star formation is significant at $\gtrsim 50\%$. Some of the dust continuum detected galaxies show spatial offsets ($\sim 0.5-1.5''$) between the rest-UV and far-IR emission peaks. These separations appear to have an increasing trend against an indicator that suggests spatially decoupled phases of obscured and unobscured star formation. REBELS offers the best available statistical constraints on obscured star formation in UV-bright, massive galaxies at $z > 6.5$. △ Less

Submitted 26 September, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 17 pages, 9 figures, 4 tables, Accepted for publication in MNRAS

Journal ref: MNRAS, 515, 3126 (2022)

The ALMA REBELS Survey: Specific Star-Formation Rates in the Reionization Era

Authors: Michael W. Topping, Daniel P. Stark, Ryan Endsley, Rychard J. Bouwens, Sander Schouws, Renske Smit, Mauro Stefanon, Hanae Inami, Rebecca A. A. Bowler, Pascal Oesch, Valentino Gonzalez, Pratika Dayal, Elisabete da Cunha, Hiddo Algera, Paul van der Werf, Andrea Pallottini, Laia Barrufet De Soto, Raffaella Schneider, Ilse De Looze, Laura Sommovigo, Lily Whitler, Luca Graziani, Yoshinobu Fudamoto, Andrea Ferrara

Abstract: We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we de… ▽ More We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we derive stellar masses from SED fitting and total SFRs from calibrations of the UV and FIR emission. The median sSFR is $18_{-5}^{+7}$ Gyr$^{-1}$, significantly larger than literature measurements lacking constraints in the FIR. The increase in sSFR reflects the larger obscured SFRs we derive from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star formation history (SFH) in reionization-era galaxies. When large sSFR galaxies are modeled with non-parametric SFHs, the derived stellar masses can increase by an order of magnitude relative to constant star formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses that are easily able to accommodate the dominant old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties in the derived parameters, we find that the sSFR increases rapidly toward higher redshifts for massive galaxies ($9.6<\log(\rm M_*/M_{\odot})<9.8$), with a power law that goes as $(1+z)^{1.7\pm0.3}$, broadly consistent with expectations from the evolving baryon accretion rates. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 18 pages, 10 figures, submitted to MNRAS

The ALMA REBELS Survey: the dust content of $z \sim 7$ Lyman Break Galaxies

Authors: P. Dayal, A. Ferrara, L. Sommovigo, R. Bouwens, P. A. Oesch, R. Smit, V. Gonzalez, S. Schouws, M. Stefanon, C. Kobayashi, J. Bremer, H. S. B. Algera, M. Aravena, R. A. A. Bowler, E. da Cunha, Y. Fudamoto, L. Graziani, J. Hodge, H. Inami, I. De Looze, A. Pallottini, D. Riechers, R. Schneider, D. Stark, R. Endsley

Abstract: We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destructi… ▽ More We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard timescale $τ_0= 30$ Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07-0.1\%$ and a UV-to-total star formation rate relation such that $log (ψ_{\rm UV}) = -0.05 ~[log (ψ)]^{2} + 0.86 ~log(ψ) -0.05$ (implying that 55-80\% of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9-10} M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at $z=7$. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94-1.1\%$) that are up to $18\times$ larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing $τ_0$ to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor $ \sim 2$. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: Accepted to MNRAS

The ALMA REBELS Survey. Epoch of Reionization giants: properties of dusty galaxies at $z \approx 7$

Authors: A. Ferrara, L. Sommovigo, P. Dayal, A. Pallottini, R. J. Bouwens, V. Gonzalez, H. Inami, R. Smit, R. A. A. Bowler, R. Endsley, P. Oesch, S. Schouws, D. Stark, M. Stefanon, M. Aravena, E. da Cunha, I. De Looze, Y. Fudamoto, L. Graziani, J. Hodge, D. Riechers, R. Schneider, H. S. B. Algera, L. Barrufet, A. P. S. Hygate , et al. (5 additional authors not shown)

Abstract: We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $β$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158μ$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different d… ▽ More We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $β$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158μ$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have (28-90.5)% of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 < {\rm SFR}/ (M_\odot {\rm yr}^{-1}) < 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 M_\odot$ and $52 \pm 11$ K, respectively. In some galaxies dust is abundant (REBELS-14, $M'_d \approx 3.4 \times 10^7 M_\odot$), or hot (REBELS-18, $T'_d \approx 67$ K). The dust distribution is compact ($<0.3$ kpc for 70% of the galaxies). The dust yield per supernova is $0.1 \le y_d/M_\odot \le 3.3$, with 70% of the galaxies requiring $y_d < 0.25 M_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_d > 1 M_\odot$. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [CII] nicely follow the local $L_{\rm CII}-$SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of $0.11\, y_P^{-2}$ Gyr, where $y_P$ is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed ($β, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(β-β_{\rm int})$, where $β_{\rm int}$ is the intrinsic UV slope, exceeds $I_m^*\approx 1120$ for a MW curve. For these objects we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-$β$ relation. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 22 pages, 7 figures, accepted by MNRAS. Comments welcome

The REBELS ALMA Survey: cosmic dust temperature evolution out to z $\sim$ 7

Authors: L. Sommovigo, A. Ferrara, A. Pallottini, P. Dayal, R. J. Bouwens, R. Smit, E. da Cunha, I. De Looze, R. A. A. Bowler, J. Hodge, H. Inami, P. Oesch, R. Endsley, V. Gonzalez, S. Schouws, D. Stark, M. Stefanon, M. Aravena, L. Graziani, D. Riechers, R. Schneider, P. van der Werf, H. Algera, L. Barrufet, Y. Fudamoto , et al. (5 additional authors not shown)

Abstract: ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star format… ▽ More ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [CII] 158$μ$m line and underlying dust continuum measurements, we derive $T_ d$ in the continuum and [CII] detected $z\approx 7$ galaxies in the ALMA Large Project REBELS sample. We find $39\ \mathrm{K} < T_d < 58\ \mathrm{K}$, and dust masses in the narrow range $M_d = (0.9-3.6)\times 10^7 M_{\odot}$. These results allow us to extend for the first time the reported $T_d(z)$ relation into the Epoch of Reionization. We produce a new physical model that explains the increasing $T_ d(z)$ trend with the decrease of gas depletion time, $t_{dep}=M_g/\mathrm{SFR}$, induced by the higher cosmological accretion rate at early times; this hypothesis yields $T_d \propto (1+z)^{0.4}$. The model also explains the observed $T_d$ scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, $T_d$ only depends on the gas column density (metallicity), $T_d \propto N_H^{1/6}$ ($T_d \propto Z^{-1/6}$). REBELS galaxies are on average relatively transparent, with effective gas column densities around $N_H \simeq (0.03-1)\times 10^{21} \mathrm{cm}^{-2}$. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated $T_d \gg 60$ K, are significantly obscured, low-metallicity systems. In fact $T_d$ is higher in metal-poor systems due to their smaller dust content, which for fixed $L_{ IR}$ results in warmer temperatures. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS

The ALMA REBELS Survey: Efficient Ly$α$ Transmission of UV-Bright z$\simeq$7 Galaxies from Large Velocity Offsets and Broad Line Widths

Authors: Ryan Endsley, Daniel P. Stark, Rychard J. Bouwens, Sander Schouws, Renske Smit, Mauro Stefanon, Hanae Inami, Rebecca A. A. Bowler, Pascal Oesch, Valentino Gonzalez, Manuel Aravena, Elisabete da Cunha, Pratika Dayal, Andrea Ferrara, Luca Graziani, Themiya Nanayakkara, Andrea Pallottini, Raffaella Schneider, Laura Sommovigo, Michael Topping, Paul van der Werf, Anne Hutter

Abstract: Recent work has shown that UV-luminous reionization-era galaxies often exhibit strong Lyman-alpha emission despite being situated at redshifts where the IGM is thought to be substantially neutral. It has been argued that this enhanced Ly$α$ transmission reflects the presence of massive galaxies in overdense regions which power large ionized bubbles. An alternative explanation is that massive galax… ▽ More Recent work has shown that UV-luminous reionization-era galaxies often exhibit strong Lyman-alpha emission despite being situated at redshifts where the IGM is thought to be substantially neutral. It has been argued that this enhanced Ly$α$ transmission reflects the presence of massive galaxies in overdense regions which power large ionized bubbles. An alternative explanation is that massive galaxies shift more of their Ly$α$ profile to large velocities (relative to the systemic redshift) where the IGM damping wing absorption is reduced. Such a mass-dependent trend is seen at lower redshifts, but whether one exists at $z\sim7$ remains unclear owing to the small number of existing systemic redshift measurements in the reionization era. This is now changing with the emergence of [CII]-based redshifts from ALMA. Here we report MMT/Binospec Ly$α$ spectroscopy of eight UV-bright ($\mathrm{M_{UV}}^{}\sim-22$) galaxies at $z\simeq7$ selected from the ALMA REBELS survey. We detect Ly$α$ in 4 of 8 galaxies and use the [CII] systemic redshifts to investigate the Ly$α$ velocity profiles. The Ly$α$ lines are significantly redshifted from systemic (average velocity offset=223 km/s) and broad (FWHM$\approx$300$-$650 km/s), with two sources showing emission extending to $\approx$750 km/s. We find that the broadest Ly$α$ profiles are associated with the largest [CII] line widths, suggesting a potential link between the Ly$α$ FWHM and the dynamical mass. Since Ly$α$ photons at high velocities transmit efficiently through the $z=7$ IGM, our data suggest that velocity profiles play a significant role in boosting the Ly$α$ visibility of the most UV-luminous reionization-era galaxies. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: 17 pages, 7 figures, submitted to MNRAS

Systematic biases in determining dust attenuation curves through galaxy SED fitting

Authors: Jianbo Qin, Xian Zhong Zheng, Min Fang, Zhizheng Pan, Stijn Wuyts, Yong Shi, Yingjie Peng, Valentino Gonzalez, Fuyan Bian, Jia-Sheng Huang, Qiu-Sheng Gu, Wenhao Liu, Qinghua Tan, Dong Dong Shi, Jian Ren, Yuheng Zhang, Man Qiao, Run Wen, Shuang Liu

Abstract: While the slope of the dust attenuation curve ($δ$) is found to correlate with effective dust attenuation ($A_V$) as obtained through spectral energy distribution (SED) fitting, it remains unknown how the fitting degeneracies shape this relation. We examine the degeneracy effects by fitting SEDs of a sample of local star-forming galaxies (SFGs) selected from the Galaxy And Mass Assembly survey, in… ▽ More While the slope of the dust attenuation curve ($δ$) is found to correlate with effective dust attenuation ($A_V$) as obtained through spectral energy distribution (SED) fitting, it remains unknown how the fitting degeneracies shape this relation. We examine the degeneracy effects by fitting SEDs of a sample of local star-forming galaxies (SFGs) selected from the Galaxy And Mass Assembly survey, in conjunction with mock galaxy SEDs of known attenuation parameters. A well-designed declining starburst star formation history is adopted to generate model SED templates with intrinsic UV slope ($β_0$) spanning over a reasonably wide range. The best-fitting $β_0$ for our sample SFGs shows a wide coverage, dramatically differing from the limited range of $β_0<-2.2$ for a starburst of constant star formation. Our results show that strong degeneracies between $β_0$, $δ$, and $A_V$ in the SED fitting induce systematic biases leading to a false $A_V$--$δ$ correlation. Our simulation tests reveal that this relationship can be well reproduced even when a flat $A_V$--$δ$ relation is taken to build the input model galaxy SEDs. The variations in best-fitting $δ$ are dominated by the fitting errors. We show that assuming a starburst with constant star formation in SED fitting will result in a steeper attenuation curve, smaller degeneracy errors, and a stronger $A_V$--$δ$ relation. Our findings confirm that the $A_V$--$δ$ relation obtained through SED fitting is likely driven by the systematic biases induced by the fitting degeneracies between $β_0$, $δ$, and $A_V$. △ Less

Submitted 14 January, 2022; originally announced January 2022.

Comments: 21 pages, 13 figures, accepted for publication in the MNRAS, Comments welcome!

The cosmic environment overtakes the local density in shaping galaxy star formation

Authors: Jian Ren, Zhizheng Pan, XianZhong Zheng, Jianbo Qin, DongDong Shi, Valentino Gonzalez, Fuyan Bian, Jia-Sheng Huang, Min Fang, Wenhao Liu, Run Wen, Yuheng Zhang, Man Qiao, Shuang Liu

Abstract: The gas supply from the cosmic web is the key to sustain star formation in galaxies. It remains to be explored how the cosmic large-scale structure (LSS) effects on galaxy evolution at given local environments. We examine galaxy specific star formation rate as a function of local density in a LSS at $z=0.735$ in the Extended Chandra Deep Field South. The LSS is mapped by 732 galaxies with $R<24$\,… ▽ More The gas supply from the cosmic web is the key to sustain star formation in galaxies. It remains to be explored how the cosmic large-scale structure (LSS) effects on galaxy evolution at given local environments. We examine galaxy specific star formation rate as a function of local density in a LSS at $z=0.735$ in the Extended Chandra Deep Field South. The LSS is mapped by 732 galaxies with $R<24$\,mag and redshift at $0.72\le z \le 0.75$ collected from the literature and our spectroscopic observations with Magellan/IMACS, consisting of five galaxy clusters/groups and surrounding filaments over an area of $23.9 \times22.7$\,co-moving\,Mpc$^2$. The spread of spectroscopic redshifts corresponds a velocity dispersion of 494\,km\,s$^{-1}$, indicating the LSS likely to be a thin sheet with a galaxy density $\gtrsim 3.9$ times that of the general field. These clusters/groups in this LSS mostly exhibit elongated morphologies and multiple components connected with surrounding filaments. Strikingly, we find that star-forming galaxies in the LSS keep star formation at the same level as field, and show no dependence on local density but stellar mass. Meanwhile, an increasing fraction of quiescent galaxies is detected at increasing local density in both the LSS and the field, consistent with the expectation that galaxy mass and local dense environment hold the key to quench star formation. Combined together, we conclude that the cosmic environment of the LSS overtakes the local environment in remaining galaxy star formation to the level of the field. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Comments: 14 pages, 13 figures, Accepted for publication in MNRAS

The Spitzer/IRAC Legacy over the GOODS Fields: Full-Depth 3.6, 4.5, 5.8 and 8.0um Mosaics and Photometry for > 9000 Galaxies at z~3.5-10 from the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS)

Authors: Mauro Stefanon, Ivo Labbé, Pascal A. Oesch, Stephane de Barros, Valentino Gonzalez, Rychard J. Bouwens, Marijn Franx, Garth D. Illingworth, Brad Holden, Dan Magee, Renske Smit, Pieter van Dokkum

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

Submitted 12 October, 2021; originally announced October 2021.

Comments: Accepted for publication on ApJS. Mosaics, PSFs and catalogs are publicly available at https://home.strw.leidenuniv.nl/~stefanon/GREATS/

Normal, Dust-Obscured Galaxies in the Epoch of Reionization

Authors: Y. Fudamoto, P. A. Oesch, S. Schouws, M. Stefanon, R. Smit, R. J. Bouwens, R. A. A. Bowler, R. Endsley, V. Gonzalez, H. Inami, I. Labbe, D. Stark, M. Aravena, L. Barrufet, E. da Cunha, P. Dayal, A. Ferrara, L. Graziani, J. Hodge, A. Hutter, Y. Li, I. De Looze, T. Nanayakkara, A. Pallottini, D. Riechers , et al. (4 additional authors not shown)

Abstract: Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extrem… ▽ More Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extreme starburst galaxies, and companions of rare quasars. These studies conclude that the contribution of dust-obscured galaxies to the cosmic star formation rate density at $z>6$ is sub-dominant. Recent ALMA and Spitzer observations have identified a more abundant, less extreme population of obscured galaxies at $z=3-6$. However, this population has not been confirmed in the reionization epoch so far. Here, we report the discovery of two dust-obscured star forming galaxies at $z=6.6813\pm0.0005$ and $z=7.3521\pm0.0005$. These objects are not detected in existing rest-frame UV data, and were only discovered through their far-infrared [CII] lines and dust continuum emission as companions to typical UV-luminous galaxies at the same redshift. The two galaxies exhibit lower infrared luminosities and star-formation rates than extreme starbursts, in line with typical star-forming galaxies at $z\sim7$. This population of heavily dust-obscured galaxies appears to contribute 10-25 per cent to the $z>6$ cosmic star formation rate density. △ Less

Submitted 21 September, 2021; originally announced September 2021.

Comments: 16 pages, 7 figures, 4 tables, published in Nature

Reionization Era Bright Emission Line Survey: Selection and Characterization of Luminous Interstellar Medium Reservoirs in the z>6.5 Universe

Authors: R. J. Bouwens, R. Smit, S. Schouws, M. Stefanon, R. Bowler, R. Endsley, V. Gonzalez, H. Inami, D. Stark, P. Oesch, J. Hodge, M. Aravena, E. da Cunha, P. Dayal, I. de Looze, A. Ferrara, Y. Fudamoto, L. Graziani, C. Li, T. Nanayakkara, A. Pallotini, R. Schneider, L. Sommovigo, M. Topping, P. van der Werf , et al. (6 additional authors not shown)

Abstract: The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z>6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7-deg**2 area for bright 158-micron [CII] and 88… ▽ More The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z>6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7-deg**2 area for bright 158-micron [CII] and 88-micron [OIII] lines and dust-continuum emission. Selection of the 40 REBELS targets was done by combining our own and other photometric selections, each of which is subject to extensive vetting using three completely independent sets of photometry and template-fitting codes. Building on the observational strategy deployed in two pilot programs, we are increasing the number of massive interstellar medium (ISM) reservoirs known at z>6.5 by ~4-5x to >30. In this manuscript, we motivate the observational strategy deployed in the REBELS program and present initial results. Based on the 60.6 hours of ALMA observations taken in the first year of the program (November 2019 to January 2020), 18 highly significant >~7sigma [CII] lines have already been discovered, the bulk of which (13/18) also show >~3.3 sigma dust-continuum emission. These newly discovered lines more than triple the number of bright ISM-cooling lines known in the z>6.5 universe, such that the number of ALMA-derived redshifts at z>6.5 already rival Lya redshift discoveries. An analysis of the completeness of our search results vs. star formation rate (SFR) suggests an ~79% efficiency in scanning for [CII] when the SFR(UV+IR) is in excess of 28 M_sol/yr. These new LP results further demonstrate ALMA's efficiency as a "redshift machine", particularly in the epoch of reionization. △ Less

Submitted 2 March, 2022; v1 submitted 25 June, 2021; originally announced June 2021.

Comments: 30 pages, 4 tables, 18 figures, accepted for publication in ApJ

Galaxy Stellar Mass Functions from z~10 to z~6 using the Deepest Spitzer/IRAC Data: No Significant Evolution in the Stellar-to-Halo Mass Ratio of Galaxies in the First Gyr of Cosmic Time

Authors: Mauro Stefanon, Rychard J. Bouwens, Ivo Labbé, Garth D. Illingworth, Valentino Gonzalez, Pascal A. Oesch

Abstract: We present new stellar mass functions at $z\sim6$, $z\sim7$, $z\sim8$, $z\sim9$ and, for the first time, $z\sim10$, constructed from $\sim800$ Lyman-Break galaxies previously identified over the XDF/UDF, parallels and the five CANDELS fields. Our study is distinctive due to (1) the much deeper ($\sim200$ hour) wide-area Spitzer/IRAC imaging at $3.6μ$m and $4.5μ$m from the GOODS Re-ionization Era w… ▽ More We present new stellar mass functions at $z\sim6$, $z\sim7$, $z\sim8$, $z\sim9$ and, for the first time, $z\sim10$, constructed from $\sim800$ Lyman-Break galaxies previously identified over the XDF/UDF, parallels and the five CANDELS fields. Our study is distinctive due to (1) the much deeper ($\sim200$ hour) wide-area Spitzer/IRAC imaging at $3.6μ$m and $4.5μ$m from the GOODS Re-ionization Era wide Area Treasury from Spitzer (GREATS) program and (2) consideration of $z\sim6-10$ sources over a $3\times$ larger area than previous HST+Spitzer studies. The Spitzer/IRAC data enable $\ge2σ$ rest-frame optical detections for an unprecedented $50\%$ of galaxies down to a stellar mass limit of $\sim10^{8}\mathcal{M}_\odot$ across all redshifts. Schechter fits to our volume densities suggest a combined evolution in characteristic mass $\mathcal{M}^*$ and normalization factor $φ^*$ between $z\sim6$ and $z\sim8$. The stellar mass density (SMD) increases by $\sim1000\times$ in the $\sim500$ Myr between $z\sim10$ and $z\sim6$, with indications of a steeper evolution between $z\sim10$ and $z\sim8$, similar to the previously-reported trend of the star-formation rate density. Strikingly, abundance matching to the Bolshoi-Planck simulation indicates halo mass densities evolving at approximately the same rate as the SMD between $z\sim10$ and $z\sim4$. Our results show that the stellar-to-halo mass ratios, a proxy for the star-formation efficiency, do not change significantly over the huge stellar mass build-up occurred from $z\sim10$ to $z\sim6$, indicating that the assembly of stellar mass closely mirrors the build-up in halo mass in the first $\sim1$ Gyr of cosmic history. JWST is poised to extend these results into the "first galaxy" epoch at $z\gtrsim10$. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: Submitted to ApJ, comments welcome

Blue Rest-Frame UV-Optical Colors in z~8 Galaxies from GREATS: Very Young Stellar Populations at ~650 Myr of Cosmic Time

Authors: Mauro Stefanon, Rychard J. Bouwens, Ivo Labbé, Garth D. Illingworth, Pascal A. Oesch, Pieter van Dokkum, Valentino Gonzalez

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

Submitted 10 March, 2021; originally announced March 2021.

Comments: Submitted to ApJ. Comments welcome. See Fig. 2 & 3 for the stacks, Fig. 4 & 6 for the colors/age and Fig. 8 for the sSFR evolution

LATIS: The Ly$α$ Tomography IMACS Survey

Authors: Andrew B. Newman, Gwen C. Rudie, Guillermo A. Blanc, Daniel D. Kelson, Sunny Rhoades, Tyson Hare, Victoria Pérez, Andrew J. Benson, Alan Dressler, Valentino Gonzalez, Juna A. Kollmeier, Nicholas P. Konidaris, John S. Mulchaey, Michael Rauch, Olivier Le Fèvre, Brian C. Lemaux, Olga Cucciati, Simon J. Lilly

Abstract: We introduce LATIS, the Ly$α$ Tomography IMACS Survey, a spectroscopic survey at Magellan designed to map the z=2.2-2.8 intergalactic medium (IGM) in three dimensions by observing the Ly$α$ forest in the spectra of galaxies and QSOs. Within an area of 1.7 deg${}^2$, we will observe approximately half of $\gtrsim L^*$ galaxies at z=2.2-3.2 for typically 12 hours, providing a dense network of sightl… ▽ More We introduce LATIS, the Ly$α$ Tomography IMACS Survey, a spectroscopic survey at Magellan designed to map the z=2.2-2.8 intergalactic medium (IGM) in three dimensions by observing the Ly$α$ forest in the spectra of galaxies and QSOs. Within an area of 1.7 deg${}^2$, we will observe approximately half of $\gtrsim L^*$ galaxies at z=2.2-3.2 for typically 12 hours, providing a dense network of sightlines piercing the IGM with an average transverse separation of 2.5 $h^{-1}$ comoving Mpc (1 physical Mpc). At these scales, the opacity of the IGM is expected to be closely related to the dark matter density, and LATIS will therefore map the density field in the $z \sim 2.5$ universe at $\sim$Mpc resolution over the largest volume to date. Ultimately LATIS will produce approximately 3800 spectra of z=2.2-3.2 galaxies that probe the IGM within a volume of $4 \times 10^6 h^{-3}$ Mpc${}^3$, large enough to contain a representative sample of structures from protoclusters to large voids. Observations are already complete over one-third of the survey area. In this paper, we describe the survey design and execution. We present the largest IGM tomographic maps at comparable resolution yet made. We show that the recovered matter overdensities are broadly consistent with cosmological expectations based on realistic mock surveys, that they correspond to galaxy overdensities, and that we can recover structures identified using other tracers. LATIS is conducted in Canada-France-Hawaii Telescope Legacy Survey fields, including COSMOS. Coupling the LATIS tomographic maps with the rich data sets collected in these fields will enable novel studies of environment-dependent galaxy evolution and the galaxy-IGM connection at cosmic noon. △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: accepted to ApJ; Fig. 26 contains animated rendering of the IGM maps

The high redshift SFR-M* relation is sensitive to the employed star formation rate and stellar mass indicators: Towards addressing the tension between observations and simulations

Authors: A. Katsianis, V Gonzalez, D. Barrientos, X. Yang, C. D. P. Lagos, J. Schaye, P. Camps, A. Trčka, M. Baes, M. Stalevski, G. A. Blanc, T. Theuns

Abstract: There is a severe tension between the observed star formation rate (SFR) - stellar mass (${\rm M}_{\star}$) relations reported by different authors at $z = 1-4$. In addition, the observations have not been successfully reproduced by state-of-the-art cosmological simulations which tend to predict a factor of 2-4 smaller SFRs at a fixed ${\rm M}_{\star}$. We examine the evolution of the SFR… ▽ More There is a severe tension between the observed star formation rate (SFR) - stellar mass (${\rm M}_{\star}$) relations reported by different authors at $z = 1-4$. In addition, the observations have not been successfully reproduced by state-of-the-art cosmological simulations which tend to predict a factor of 2-4 smaller SFRs at a fixed ${\rm M}_{\star}$. We examine the evolution of the SFR$-{\rm M}_{\star}$ relation of $z = 1-4 $ galaxies using the SKIRT simulated spectral energy distributions of galaxies sampled from the EAGLE simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR$-{\rm M}_{\star}$ relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 ${\rm μm}$, MIPS - 24, 70 and 160 ${\rm μm}$ or SPIRE - 250, 350, 500 ${\rm μm}$) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by -0.5 dex at $z \sim 4$ but overpredict the measurements by 0.3 dex at $z \sim 1$. Relations relying on dust-corrected rest-frame UV luminosities, are flatter since they overpredict/underpredict SFRs for low/high star forming objects and yield deviations from the intrinsic relation from 0.10 dex to -0.13 dex at $z \sim 4$. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs. △ Less

Submitted 16 January, 2020; originally announced January 2020.

Comments: 15 pages, 6 Figures, 4 tables, accepted for publication in MNRAS, Figures 1, 2, 3, 4 summarize results

The Hubble Legacy Field GOODS-S Photometric Catalog

Authors: Katherine E. Whitaker, Mohammad Ashas, Garth Illingworth, Daniel Magee, Joel Leja, Pascal Oesch, Pieter van Dokkum, Lamiya Mowla, Rychard Bouwens, Marijn Franx, Bradford Holden, Ivo Labbé, Marc Rafelski, Harry Teplitz, Valentino Gonzalez

Abstract: This manuscript describes the public release of the Hubble Legacy Fields (HLF) project photometric catalog for the extended GOODS-South region from the Hubble Space Telescope (HST) archival program AR-13252. The analysis is based on the version 2.0 HLF data release that now includes all ultraviolet (UV) imaging, combining three major UV surveys. The HLF data combines over a decade worth of 7475 ex… ▽ More This manuscript describes the public release of the Hubble Legacy Fields (HLF) project photometric catalog for the extended GOODS-South region from the Hubble Space Telescope (HST) archival program AR-13252. The analysis is based on the version 2.0 HLF data release that now includes all ultraviolet (UV) imaging, combining three major UV surveys. The HLF data combines over a decade worth of 7475 exposures taken in 2635 orbits totaling 6.3 Msec with the HST Advanced Camera for Surveys Wide Field Channel (ACS/WFC) and the Wide Field Camera 3 (WFC3) UVIS/IR Channels in the greater GOODS-S extragalactic field, covering all major observational efforts (e.g., GOODS, GEMS, CANDELS, ERS, UVUDF and many other programs; see Illingworth et al 2019, in prep). The HLF GOODS-S catalogs include photometry in 13 bandpasses from the UV (WFC3/UVIS F225W, F275W and F336W filters), optical (ACS/WFC F435W, F606W, F775W, F814W and F850LP filters), to near-infrared (WFC3/IR F098M, F105W, F125W, F140W and F160W filters). Such a data set makes it possible to construct the spectral energy distributions (SEDs) of objects over a wide wavelength range from high resolution mosaics that are largely contiguous. Here, we describe a photometric analysis of 186,474 objects in the HST imaging at wavelengths 0.2--1.6$μ$m. We detect objects from an ultra-deep image combining the PSF-homogenized and noise-equalized F850LP, F125W, F140W and F160W images, including Gaia astrometric corrections. SEDs were determined by carefully taking the effects of the point-spread function in each observation into account. All of the data presented herein are available through the HLF website (https://archive.stsci.edu/prepds/hlf/). △ Less

Submitted 15 August, 2019; originally announced August 2019.

Comments: Hubble Legacy Fields GOODS-S public data release available at https://archive.stsci.edu/prepds/hlf/. Accepted for publication in ApJS (20 pages, 22 figures, 2 tables)

An evolving and mass dependent ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation for galaxies

Authors: Antonios Katsianis, Xianzhong Zheng, Valentino Gonzalez, Guillermo Blanc, Claudia del P. Lagos, Luke J. M. Davies, Peter Camps, Ana Trčka, Maarten Baes, Joop Schaye, James W. Trayford, Tom Theuns, Marko Stalevski

Abstract: The scatter (${\rmσ_{\text{sSFR}}}$) of the specific star formation rates (sSFRs) of galaxies is a measure of the diversity in their star formation histories (SFHs) at a given mass. In this paper we employ the EAGLE simulations to study the dependence of the ${\rm σ_{\text{sSFR}}}$ of galaxies on stellar mass (${\rm M_{\star}}$) through the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation in… ▽ More The scatter (${\rmσ_{\text{sSFR}}}$) of the specific star formation rates (sSFRs) of galaxies is a measure of the diversity in their star formation histories (SFHs) at a given mass. In this paper we employ the EAGLE simulations to study the dependence of the ${\rm σ_{\text{sSFR}}}$ of galaxies on stellar mass (${\rm M_{\star}}$) through the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation in $ {\rm z \sim 0-4}$. We find that the relation evolves with time, with the dispersion depending on both stellar mass and redshift. The models point to an evolving U-shape form for the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation with the scatter being minimal at a characteristic mass $M^{\star}$ of ${\rm 10^{9.5}}$ ${\rm M_{\odot}}$ and increasing both at lower and higher masses. This implication is that the diversity of SFHs increases towards both at the low- and high-mass ends. We find that active galactic nuclei feedback is important for increasing the ${\rm σ_{\text{sSFR}}}$ for high mass objects. On the other hand, we suggest that SNe feedback increases the ${\rm σ_{\text{sSFR}}}$ of galaxies at the low-mass end. We also find that excluding galaxies that have experienced recent mergers does not significantly affect the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation. Furthermore, we employ the combination of the EAGLE simulations with the radiative transfer code SKIRT to evaluate the effect of SFR/stellar mass diagnostics in the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation and find that the ${\rm SFR/M_{\star}}$ methodologies (e.g. SED fitting, UV+IR, UV+IRX-$β$) widely used in the literature to obtain intrinsic properties of galaxies have a large effect on the derived shape and normalization of the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation. △ Less

Submitted 6 May, 2019; originally announced May 2019.

Comments: 23 pages, 10 figures, Accepted to ApJ

The GREATS H$β$+[OIII] Luminosity Function and Galaxy Properties at $\mathbf{z\sim8}$: Walking the Way of JWST

Authors: S. De Barros, P. A. Oesch, I. Labbé, M. Stefanon, V. González, R. Smit, R. J. Bouwens, G. D. Illingworth

Abstract: The James Webb Space Telescope will allow to spectroscopically study an unprecedented number of galaxies deep into the reionization era, notably by detecting [OIII] and H$β$ nebular emission lines. To efficiently prepare such observations, we photometrically select a large sample of galaxies at $z\sim8$ and study their rest-frame optical emission lines. Combining data from the GOODS Re-ionization… ▽ More The James Webb Space Telescope will allow to spectroscopically study an unprecedented number of galaxies deep into the reionization era, notably by detecting [OIII] and H$β$ nebular emission lines. To efficiently prepare such observations, we photometrically select a large sample of galaxies at $z\sim8$ and study their rest-frame optical emission lines. Combining data from the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) survey and from HST, we perform spectral energy distribution (SED) fitting, using synthetic SEDs from a large grid of photoionization models. The deep Spitzer/IRAC data combined with our models exploring a large parameter space enables to constrain the [OIII]+H$β$ fluxes and equivalent widths for our sample, as well as the average physical properties of $z\sim8$ galaxies, such as the ionizing photon production efficiency with $\log(ξ_\mathrm{ion}/\mathrm{erg}^{-1}\hspace{1mm}\mathrm{Hz})\geq25.77$. We find a relatively tight correlation between the [OIII]+H$β$ and UV luminosity, which we use to derive for the first time the [OIII]+H$β$ luminosity function (LF) at $z\sim8$. The $z\sim8$ [OIII]+H$β$ LF is higher at all luminosities compared to lower redshift, as opposed to the UV LF, due to an increase of the [OIII]+H$β$ luminosity at a given UV luminosity from $z\sim3$ to $z\sim8$. Finally, using the [OIII]+H$β$ LF, we make predictions for JWST/NIRSpec number counts of $z\sim8$ galaxies. We find that the current wide-area extragalactic legacy fields are too shallow to use JWST at maximal efficiency for $z\sim8$ spectroscopy even at 1hr depth and JWST pre-imaging to $\gtrsim30$ mag will be required. △ Less

Submitted 22 March, 2019; originally announced March 2019.

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

MUSCAT: The Mexico-UK Sub-Millimetre Camera for AsTronomy

Authors: Thomas L. R. Brien, Peter A. R. Ade, Peter S. Barry, Edgar Castillo-Domìnguez, Daniel Ferrusca, Thomas Gascard, Victor Gómez, Peter C. Hargrave, Amber L. Hornsby, David Hughes, Enzo Pascale, Josie D. A. Parrianen, Abel Perez, Sam Rowe, Carole Tucker, Salvador Ventura González, Simon M. Doyle

Abstract: The Mexico-UK Sub-millimetre Camera for AsTronomy (MUSCAT) is a large-format, millimetre-wave camera consisting of 1,500 background-limited lumped-element kinetic inductance detectors (LEKIDs) scheduled for deployment on the Large Millimeter Telescope (Volcán Sierra Negra, Mexico) in 2018. MUSCAT is designed for observing at 1.1 mm and will utilise the full 40' field of view of the LMTs upgraded 5… ▽ More The Mexico-UK Sub-millimetre Camera for AsTronomy (MUSCAT) is a large-format, millimetre-wave camera consisting of 1,500 background-limited lumped-element kinetic inductance detectors (LEKIDs) scheduled for deployment on the Large Millimeter Telescope (Volcán Sierra Negra, Mexico) in 2018. MUSCAT is designed for observing at 1.1 mm and will utilise the full 40' field of view of the LMTs upgraded 50-m primary mirror. In its primary role, MUSCAT is designed for high-resolution follow-up surveys of both galactic and extra-galactic sub-mm sources identified by Herschel. MUSCAT is also designed to be a technology demonstrator that will provide the first on-sky demonstrations of novel design concepts such as horn-coupled LEKID arrays and closed continuous cycle miniature dilution refrigeration. Here we describe some of the key design elements of the MUSCAT instrument such as the novel use of continuous sorption refrigerators and a miniature dilutor for continuous 100-mK cooling of the focal plane, broadband optical coupling to Aluminium LEKID arrays using waveguide chokes and anti-reflection coating materials as well as with the general mechanical and optical design of MUSCAT. We explain how MUSCAT is designed to be simple to upgrade and the possibilities for changing the focal plane unit that allows MUSCAT to act as a demonstrator for other novel technologies such as multi-chroic polarisation sensitive pixels and on-chip spectrometry in the future. Finally, we will report on the current status of MUSCAT's commissioning. △ Less

Submitted 23 July, 2018; originally announced July 2018.

Comments: Presented at SPIE Astronomical Telescopes + Instrumentation, 2018, Austin, Texas, United States

The evolution of the star formation rate function in the EAGLE simulations: A comparison with UV, IR and H$α$ observations from z ~ 8 to z ~ 0

Authors: A. Katsianis, G. Blanc, C. P. Lagos, N. Tejos, R. G. Bower, A. Alavi, V. Gonzalez, T. Theuns, M. Schaller, S. Lopez

Abstract: We investigate the evolution of the galaxy Star Formation Rate Function (SFRF) and Cosmic Star Formation Rate Density (CSFRD) of $z\sim 0-8 $ galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. In addition, we present a compilation of UV, IR and H$α$ SFRFs and compare these with the predictions from the EAGLE suite of cosmological hydrodynamic simulations… ▽ More We investigate the evolution of the galaxy Star Formation Rate Function (SFRF) and Cosmic Star Formation Rate Density (CSFRD) of $z\sim 0-8 $ galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. In addition, we present a compilation of UV, IR and H$α$ SFRFs and compare these with the predictions from the EAGLE suite of cosmological hydrodynamic simulations. We find that the constraints implied by different indicators are inconsistent with each other for the highest star-forming objects at z < 2, a problem that is possibly related to selection biases and the uncertainties of dust attenuation effects. EAGLE's feedback parameters were calibrated to reproduce realistic galaxy sizes and stellar masses at z = 0.1. In this work we test if and why those choices yield realistic Star Formation Rates (SFRs) for $z \sim 0-8$ as well. We demonstrate that SNe feedback plays a major role at setting the abundance of galaxies at all star-forming regimes, especially at high redshifts. On the contrary, Active Galactic Nuclei (AGN) feedback becomes more prominent at lower redshifts and is a major mechanism that affects only the highest star-forming systems. Furthermore, we find that galaxies with SFR $\sim 1-10 \, {\rm M_{\odot} \, yr^{-1}}$ dominate the CSFRD at redshifts z < 5, while rare high star-forming galaxies (SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$) contribute significantly only briefly around the peak era ($z \sim 2$) and then are quenched by AGN feedback. In the absence of this prescription objects with SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$ would dominate the CSFRD, while the cosmic budget of star formation would be extremely high. Finally, we demonstrate that the majority of the cosmic star formation occurs in relatively rare high mass halos ($ {\rm M_{Halo}} \sim 10^{11-13} \, {\rm M_{\odot}}$) even at the earliest epochs. △ Less

Submitted 6 August, 2017; originally announced August 2017.

Comments: 21 pages, 12 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)

A Comprehensive Study of Ly$α$ Emission in the High-redshift Galaxy Population

Authors: Grecco A. Oyarzún, Guillermo A. Blanc, Valentino González, Mario Mateo, John I. Bailey III

Abstract: We present an exhaustive census of Lyman alpha (Ly$α$) emission in the general galaxy population at $3<z<4.6$. We use the Michigan/Magellan Fiber System (M2FS) spectrograph to study a stellar mass (M$_*$) selected sample of 625 galaxies homogeneously distributed in the range $7.6<\log{\mbox{M$_*$/M$_{\odot}$}}<10.6$. Our sample is selected from the 3D-HST/CANDELS survey, which provides the complem… ▽ More We present an exhaustive census of Lyman alpha (Ly$α$) emission in the general galaxy population at $3<z<4.6$. We use the Michigan/Magellan Fiber System (M2FS) spectrograph to study a stellar mass (M$_*$) selected sample of 625 galaxies homogeneously distributed in the range $7.6<\log{\mbox{M$_*$/M$_{\odot}$}}<10.6$. Our sample is selected from the 3D-HST/CANDELS survey, which provides the complementary data to estimate Ly$α$ equivalent widths ($W_{Lyα}$) and escape fractions ($f_{esc}$) for our galaxies. We find both quantities to anti-correlate with M$_*$, star-formation rate (SFR), UV luminosity, and UV slope ($β$). We then model the $W_{Lyα}$ distribution as a function of M$_{UV}$ and $β$ using a Bayesian approach. Based on our model and matching the properties of typical Lyman break galaxy (LBG) selections, we conclude that the $W_{Lyα}$ distribution in such samples is heavily dependent on the limiting M$_{UV}$ of the survey. Regarding narrowband surveys, we find their $W_{Lyα}$ selections to bias samples toward low M$_*$, while their line-flux limitations preferentially leave out low-SFR galaxies. We can also use our model to predict the fraction of Ly$α$-emitting LBGs at $4\leqslant z\leqslant 7$. We show that reported drops in the Ly$α$ fraction at $z\geqslant6$, usually attributed to the rapidly increasing neutral gas fraction of the universe, can also be explained by survey M$_{UV}$ incompleteness. This result does not dismiss reionization occurring at $z\sim7$, but highlights that current data is not inconsistent with this process taking place at $z>7$. △ Less

Submitted 8 June, 2017; v1 submitted 6 June, 2017; originally announced June 2017.

Comments: Accepted for publication in ApJ; Figure 3 updated

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

Galaxy-galaxy lensing estimators and their covariance properties

Authors: Sukhdeep Singh, Rachel Mandelbaum, Uroš Seljak, Anže Slosar, Jose Vazquez Gonzalez

Abstract: We study the covariance properties of real space correlation function estimators -- primarily galaxy-shear correlations, or galaxy-galaxy lensing -- using SDSS data for both shear catalogs and lenses (specifically the BOSS LOWZ sample). Using mock catalogs of lenses and sources, we disentangle the various contributions to the covariance matrix and compare them with a simple analytical model. We sh… ▽ More We study the covariance properties of real space correlation function estimators -- primarily galaxy-shear correlations, or galaxy-galaxy lensing -- using SDSS data for both shear catalogs and lenses (specifically the BOSS LOWZ sample). Using mock catalogs of lenses and sources, we disentangle the various contributions to the covariance matrix and compare them with a simple analytical model. We show that not subtracting the lensing measurement around random points from the measurement around the lens sample is equivalent to performing the measurement using the lens density field instead of the lens over-density field. While the measurement using the lens density field is unbiased (in the absence of systematics), its error is significantly larger due to an additional term in the covariance. Therefore, this subtraction should be performed regardless of its beneficial effects on systematics. Comparing the error estimates from data and mocks for estimators that involve the over-density, we find that the errors are dominated by the shape noise and lens clustering, that empirically estimated covariances (jackknife and standard deviation across mocks) are consistent with theoretical estimates, and that both the connected parts of the 4-point function and the super-sample covariance can be neglected for the current levels of noise. While the trade-off between different terms in the covariance depends on the survey configuration (area, source number density), the diagnostics that we use in this work should be useful for future works to test their empirically-determined covariances. △ Less

Submitted 23 August, 2017; v1 submitted 2 November, 2016; originally announced November 2016.

Comments: 18 pages. Summary in figure 7 and table 1. Matches accepted version in MNRAS

The Hubble Legacy Fields (HLF-GOODS-S) v1.5 Data Products: Combining 2442 Orbits of GOODS-S/CDF-S Region ACS and WFC3/IR Images

Authors: Garth Illingworth, Daniel Magee, Rychard Bouwens, Pascal Oesch, Ivo Labbe, Pieter van Dokkum, Katherine Whitaker, Bradford Holden, Marijn Franx, Valentino Gonzalez

Abstract: We have submitted to MAST the 1.5 version data release of the Hubble Legacy Fields (HLF) project covering a 25 x 25 arcmin area over the GOODS-S (ECDF-S) region from the HST archival program AR-13252. The release combines exposures from Hubble's two main cameras, the Advanced Camera for Surveys (ACS/WFC) and the Wide Field Camera 3 (WFC3/IR), taken over more than a decade between mid-2002 to the e… ▽ More We have submitted to MAST the 1.5 version data release of the Hubble Legacy Fields (HLF) project covering a 25 x 25 arcmin area over the GOODS-S (ECDF-S) region from the HST archival program AR-13252. The release combines exposures from Hubble's two main cameras, the Advanced Camera for Surveys (ACS/WFC) and the Wide Field Camera 3 (WFC3/IR), taken over more than a decade between mid-2002 to the end of 2016. The HLF includes essentially all optical (ACS/WFC F435W, F606W, F775W, F814W and F850LP filters) and infrared (WFC3/ IR F098M, F105W, F125W, F140W and F160W filters) data taken by Hubble over the original CDF-S region including the GOODS-S, ERS, CANDELS and many other programs (31 in total). The data has been released at https://archive.stsci.edu/prepds/hlf/ as images with a common astrometric reference frame, with corresponding inverse variance weight maps. We provide one image per filter of WFC3/IR images at 60 mas per pixel resolution and two ACS/WFC images per filter, at both 30 and 60 mas per pixel. Since this comprehensive dataset combines data from 31 programs on the GOODS-S/CDF-S, the AR proposal identified the MAST products by the global name "Hubble Legacy Field", with this region being identified by "HLF-GOODS-S". This dataset complements that of the Frontier Fields program. The total incorporated in the HLF-GOODS-S is 5.8 Msec in 7211 exposures from 2442 orbits. This is ~70% of a HST full cycle! △ Less

Submitted 10 May, 2017; v1 submitted 2 June, 2016; originally announced June 2016.

Comments: Version 1.5 release data can be found at: https://archive.stsci.edu/prepds/hlf/ 9 pages, 3 figures - minor update in number of orbits/exposures

How Lyman Alpha Emission Depends On Galaxy Stellar Mass

Authors: Grecco A. Oyarzún, Guillermo A. Blanc, Valentino González, Mario Mateo, John I. Bailey III, Steven L. Finkelstein, Paulina Lira, Jeffrey D. Crane, Edward W. Olszewski

Abstract: In this work, we show how the stellar mass (M) of galaxies affects the 3<z<4.6 Ly-alpha equivalent width (EW) distribution. To this end, we design a sample of 629 galaxies in the M range 7.6 < logM/Msun < 10.6 from the 3D-HST/CANDELS survey. We perform spectroscopic observations of this sample using the Michigan/Magellan Fiber System, allowing us to measure Ly-alpha fluxes and use 3D-HST/CANDELS a… ▽ More In this work, we show how the stellar mass (M) of galaxies affects the 3<z<4.6 Ly-alpha equivalent width (EW) distribution. To this end, we design a sample of 629 galaxies in the M range 7.6 < logM/Msun < 10.6 from the 3D-HST/CANDELS survey. We perform spectroscopic observations of this sample using the Michigan/Magellan Fiber System, allowing us to measure Ly-alpha fluxes and use 3D-HST/CANDELS ancillary data. In order to study the Ly-alpha EW distribution dependence on M, we split the whole sample in three stellar mass bins. We find that, in all bins, the distribution is best represented by an exponential profile of the form dN(M)/dEW= A(M)exp(-EW/W0(M))/W0(M). Through a Bayesian analysis, we confirm that lower M galaxies have higher Ly-alpha EWs. We also find that the fraction A of galaxies featuring emission and the e-folding scale W0 of the distribution anti- correlate with M, recovering expressions of the forms A(M)= -0.26(.13) logM/Msun+3.01(1.2) and W0(M)= -15.6(3.5) logM/Msun +166(34). These results are crucial for proper interpretation of Ly-alpha emission trends reported in the literature that may be affected by strong M selection biases. △ Less

Submitted 11 April, 2016; originally announced April 2016.

Comments: 4 pages, 5 figures

Journal ref: ApJL, 821, L14 (2016)

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

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

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

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

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

Rest-Frame Optical Emission Lines in z~3.5 Lyman Break selected Galaxies: The Ubiquity of Unusually High [OIII]/Hbeta Ratios at 2 Gyr

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

Abstract: We present K-band spectra of rest-frame optical emission lines for 24 star-forming galaxies at z~3.2-3.7 using MOSFIRE on the Keck 1 telescope. Strong rest-frame optical [O III] and Hbeta emission lines were detected in 18 LBGs. The median flux ratio of [O III]5007 to Hbeta is 5.1+/-0.5, a factor of 5-10x higher than in local galaxies with similar stellar masses. The observed Hbeta luminosities ar… ▽ More We present K-band spectra of rest-frame optical emission lines for 24 star-forming galaxies at z~3.2-3.7 using MOSFIRE on the Keck 1 telescope. Strong rest-frame optical [O III] and Hbeta emission lines were detected in 18 LBGs. The median flux ratio of [O III]5007 to Hbeta is 5.1+/-0.5, a factor of 5-10x higher than in local galaxies with similar stellar masses. The observed Hbeta luminosities are in good agreement with expectations from the estimated star-formation rates, and none of our sources are detected in deep X-ray stacks, ruling out significant contamination by active galactic nuclei. Combining our sample with a variety of LBGs from the literature, including 49 galaxies selected in a very similar manner, we find a high median ratio of [OIII]/Hbeta = 4.8+0.8-1.7. This high ratio seems to be an ubiquitous feature of z~3-4 LBGs, very different from typical local star-forming galaxies at similar stellar masses. The only comparable systems at z~0 are those with similarly high specific star-formation rates, though ~5x lower stellar masses. High specific star-formation rates either result in a much higher ionization parameter or other unusual conditions for the interstellar medium, which result in a much higher [OIII]/Hbeta line ratio. This implies a strong relation between a global property of a galaxy, the specific star-formation rate, and the local conditions of ISM in star-forming regions. △ Less

Submitted 23 March, 2016; v1 submitted 21 January, 2014; originally announced January 2014.

Comments: 14 pages, 8 figures, 5 color, published in ApJ, updated to reflect published version

Journal ref: ApJ, 820, 73 (2016)

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

Authors: Britt F. Lundgren, Pieter van Dokkum, Marijn Franx, Ivo Labbe, Michele Trenti, Rychard Bouwens, Valentino Gonzalez, Garth Illingworth, Daniel Magee, Pascal Oesch, Massimo Stiavelli

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

Submitted 30 October, 2013; v1 submitted 28 October, 2013; originally announced October 2013.

Comments: 21 pages, 16 figures; Accepted for publication in ApJ

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