Showing 1–50 of 50 results for author: McLeod, D J

The AURORA Survey: A New Era of Emission-line Diagrams with JWST/NIRSpec

Authors: Alice E. Shapley, Ryan L. Sanders, Michael W. Topping, Naveen A. Reddy, Danielle A. Berg, Rychard J. Bouwens, Gabriel Brammer, Adam C. Carnall, Fergus Cullen, Romeel Davé, James S. Dunlop, Richard S. Ellis, N. M. Förster Schreiber, Steven R . Furlanetto, Karl Glazebrook, Garth D. Illingworth, Tucker Jones, Mariska Kriek, Derek J. McLeod, Ross J. McLure, Desika Narayanan, Pascal Oesch, Anthony J. Pahl, Max Pettini, Daniel Schaerer , et al. (6 additional authors not shown)

Abstract: We present results on the emission-line properties of z=1.4-7.5 star-forming galaxies in the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) Cycle 1 JWST/NIRSpec program. Based on its depth, continuous wavelength coverage from 1--5 microns, and medium spectral resolution (R~1000), AURORA includes detections of a large suite of nebular emission lines spanning a broad… ▽ More We present results on the emission-line properties of z=1.4-7.5 star-forming galaxies in the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) Cycle 1 JWST/NIRSpec program. Based on its depth, continuous wavelength coverage from 1--5 microns, and medium spectral resolution (R~1000), AURORA includes detections of a large suite of nebular emission lines spanning a broad range in rest wavelength. We investigate the locations of AURORA galaxies in multiple different emission-line diagrams, including traditional "BPT" diagrams of [OIII]/Hbeta vs. [NII]/Halpha, [SII]/Halpha, and [OI]/Halpha, and the "ionization-metallicity" diagram of [OIII]/[OII] (O32) vs. ([OIII]+[OII])/Hbeta (R23). We also consider a bluer rest-frame "ionization-metallicity" diagram introduced recently to characterize z>10 galaxies: [NeIII]/[OII] vs. ([NeIII]+[OII])/Hdelta; as well as longer-wavelength diagnostic diagrams extending into the rest-frame near-IR: [OIII]/Hbeta vs. [SIII]/[SII] (S32); and HeI/Pagamma and [SIII]/Pagamma vs. [FeII]/Pabeta. With a significant boost in signal-to-noise and large, representative samples of individual galaxy detections, the AURORA emission-line diagrams presented here definitively confirm a physical picture in which chemically-young, alpha-enhanced, massive stars photoionize the ISM in distant galaxies with a harder ionizing spectrum at fixed nebular metallicity than in their z~0 counterparts. We also uncover previously unseen evolution prior to z~2 in the [OIII]/Hbeta vs. [NII]/Halpha diagram, which motivates deep NIRSpec observations at even higher redshift. Finally, we present the first statistical sample of rest-frame near-IR emission-line diagnostics in star-forming galaxies at high redshift. In order to truly interpret rest-frame near-IR line ratios including [FeII], we must obtain better constraints on dust depletion in the high-redshift ISM. △ Less

Submitted 2 July, 2024; v1 submitted 28 June, 2024; originally announced July 2024.

Comments: 17 pages, 8 figures, submitted to ApJ

JWST, ALMA, and Keck Spectroscopic Constraints on the UV Luminosity Functions at z~7-14: Clumpiness and Compactness of the Brightest Galaxies in the Early Universe

Authors: Yuichi Harikane, Akio K. Inoue, Richard S. Ellis, Masami Ouchi, Yurina Nakazato, Naoki Yoshida, Yoshiaki Ono, Fengwu Sun, Riku A. Sato, Seiji Fujimoto, Nobunari Kashikawa, Derek J. McLeod, Pablo G. Perez-Gonzalez, Marcin Sawicki, Yuma Sugahara, Yi Xu, Satoshi Yamanaka, Adam C. Carnall, Fergus Cullen, James S. Dunlop, Eiichi Egami, Norman Grogin, Yuki Isobe, Anton M. Koekemoer, Nicolas Laporte , et al. (10 additional authors not shown)

Abstract: We present the number densities and physical properties of the bright galaxies spectroscopically confirmed at $z\sim7-14$. Our sample is composed of 53 galaxies at $z_\mathrm{spec}\sim7-14$, including recently-confirmed galaxies at $z_\mathrm{spec}=12.34-14.32$ with JWST, as well as new confirmations at $z_\mathrm{spec}=6.583-7.643$ with $-24< M_\mathrm{UV}< -21$ mag using ALMA and Keck. Our JWST/… ▽ More We present the number densities and physical properties of the bright galaxies spectroscopically confirmed at $z\sim7-14$. Our sample is composed of 53 galaxies at $z_\mathrm{spec}\sim7-14$, including recently-confirmed galaxies at $z_\mathrm{spec}=12.34-14.32$ with JWST, as well as new confirmations at $z_\mathrm{spec}=6.583-7.643$ with $-24< M_\mathrm{UV}< -21$ mag using ALMA and Keck. Our JWST/NIRSpec observations have also revealed that very bright galaxy candidates at $z\sim10-13$ identified from ground-based telescope images before JWST are passive galaxies at $z\sim3-4$, emphasizing the necessity of strict screening and spectroscopy in the selection of the brightest galaxies at $z>10$. The UV luminosity functions derived from these spectroscopic results are consistent with a double power-law function, showing tensions with theoretical models at the bright end. To understand the origin of the overabundance of bright galaxies, we investigate their morphologies using JWST/NIRCam high-resolution images obtained in various surveys including PRIMER and COSMOS-Web. We find that $\sim70\%$ of the bright galaxies at $z\sim7$ exhibit clumpy morphologies with multiple sub-components, suggesting merger-induced starburst activity, which is consistent with SED fitting results showing bursty star formation histories. At $z\gtrsim10$, bright galaxies are classified into two types of galaxies; extended ones with weak high-ionization emission lines, and compact ones with strong high-ionization lines including NIV]$λ$1486, indicating that at least two different processes (e.g., merger-induced starburst and compact star formation/AGN) are shaping the physical properties of the brightest galaxies at $z\gtrsim10$ and are responsible for their overabundance. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 32 pages, 22 figures, submitted to ApJ

PRIMER: JWST/MIRI reveals the evolution of star-forming structures in galaxies at z<2.5

Authors: Yipeng Lyu, Benjamin Magnelli, David Elbaz, Pablo G. Pérez-González, Camila Correa, Emanuele Daddi, Carlos Gómez-Guijarro, James S. Dunlop, Norman A. Grogin, Anton M. Koekemoer, Derek J. McLeod, Shiying Lu

Abstract: The stellar structures of star-forming galaxies (SFGs) undergo significant size growth during their mass assembly and must pass through a compaction phase as they evolve into quiescent galaxies (QGs). To shed light on the mechanisms behind this structural evolution, we study the morphology of the star-forming components of 665 SFGs at 0<z<2.5 measured using JWST/MIRI observation and compare them w… ▽ More The stellar structures of star-forming galaxies (SFGs) undergo significant size growth during their mass assembly and must pass through a compaction phase as they evolve into quiescent galaxies (QGs). To shed light on the mechanisms behind this structural evolution, we study the morphology of the star-forming components of 665 SFGs at 0<z<2.5 measured using JWST/MIRI observation and compare them with the morphology of their stellar components taken from the literature. The stellar and star-forming components of most SFGs (66%) have extended disk-like structures that are aligned with each other and are of the same size. The star-forming components of these galaxies follow a mass-size relation, similar to that followed by their stellar components. At the highest mass, the optical Sérsic index of these SFGs increases to 2.5, suggesting the presence of a dominant stellar bulge. Because their star-forming components remain disk-like, these bulges cannot have formed by secular in-situ growth. We identify a second population of galaxies lying below the MIR mass-size relation, with compact star-forming components embedded in extended stellar components (EC galaxy). These galaxies are overall rare (15%) but become more dominant (30%) at high mass ($>10^{10.5}M_\odot$). The compact star-forming components of these galaxies are also concentrated and slightly spheroidal, suggesting that this compaction phase can build dense bulge in-situ. Finally, we identify a third population of SFGs (19%), with both compact stellar and star-forming components. The density of their stellar cores resemble those of QGs and are compatible with being the descendants of EC galaxy. Overall, the structural evolution of SFGs is mainly dominated by a secular inside-out growth, which can, however, be interrupted by violent compaction phase(s) that can build dominant stellar bulges like those in massive SFGs or QGs. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 24 pages, 17 figures, submitted to A&A, comments are welcome

The JWST EXCELS survey: Too much, too young, too fast? Ultra-massive quiescent galaxies at 3 < z < 5

Authors: A. C. Carnall, F. Cullen, R. J. McLure, D. J. McLeod, R. Begley, C. T. Donnan, J. S. Dunlop, A. E. Shapley, K. Rowlands, O. Almaini, K. Z. Arellano-Córdova, L. Barrufet, A. Cimatti, R. S. Ellis, N. A. Grogin, M. L. Hamadouche, G. D. Illingworth, A. M. Koekemoer, H. -H. Leung, C. C. Lovell, P. G. Pérez-González, P. Santini, T. M. Stanton, V. Wild

Abstract: We report ultra-deep, medium-resolution spectroscopic observations for 4 quiescent galaxies with log$_{10}(M_*/\mathrm{M_\odot})>11$ at $3 < z < 5$. These data were obtained with JWST NIRSpec as part of the Early eXtragalactic Continuum and Emission Line Science (EXCELS) survey, which we introduce in this work. The first pair of galaxies are newly selected from PRIMER UDS imaging, both at… ▽ More We report ultra-deep, medium-resolution spectroscopic observations for 4 quiescent galaxies with log$_{10}(M_*/\mathrm{M_\odot})>11$ at $3 < z < 5$. These data were obtained with JWST NIRSpec as part of the Early eXtragalactic Continuum and Emission Line Science (EXCELS) survey, which we introduce in this work. The first pair of galaxies are newly selected from PRIMER UDS imaging, both at $z=4.62$ and separated by $860$ pkpc on the sky, within a larger structure for which we confirm several other members. These galaxies formed at $z\simeq8-10$, and, despite their similar stellar masses, ages, and their proximity, they exhibit very different stellar metallicities, hinting at different formation pathways. These systems could plausibly merge by the present day to produce a local massive elliptical galaxy. The other 2 ultra-massive quiescent galaxies are previously known at $z=3.99$ and $3.19$, with the latter (ZF-UDS-7329) having been the subject of debate as potentially too old and too massive to be accommodated by the $Λ$-CDM halo-mass function. Both exhibit high stellar metallicities, and for ZF-UDS-7329 we are able to measure the $α-$enhancement, obtaining [Mg/Fe] = $0.42^{+0.19}_{-0.17}$. We finally evaluate whether these 4 galaxies are consistent with the $Λ$-CDM halo-mass function using an extreme value statistics approach. We find that the $z=4.62$ objects and the $z=3.19$ object are unlikely within our area under the assumption of standard stellar fractions ($f_*\simeq0.1-0.2$). However, these objects roughly align with the most massive galaxies expected under the assumption of 100 per cent conversion of baryons to stars ($f_*$=1). Our results suggest extreme galaxy formation physics during the first billion years, but no conflict with $Λ$-CDM cosmology. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 24 pages, 9 figures, submitted to MNRAS

The NIRVANDELS Survey: the stellar and gas-phase mass-metallicity relations of star-forming galaxies at z = 3.5

Authors: T. M. Stanton, F. Cullen, R. J. McLure, A. E. Shapley, K. Z. Arellano-Córdova, R. Begley, R. Amorín, L. Barrufet, A. Calabrò, A. C. Carnall, M. Cirasuolo, J. S. Dunlop, C. T. Donnan, M. L. Hamadouche, F. -Y. Liu, D. J. McLeod, L. Pentericci, L. Pozzetti, R. L. Sanders, D. Scholte, M. W. Topping

Abstract: We present determinations of the gas-phase and stellar metallicities of a sample of 65 star-forming galaxies at $z \simeq 3.5$ using rest-frame far-ultraviolet (FUV) spectroscopy from the VANDELS survey in combination with follow-up rest-frame optical spectroscopy from VLT/KMOS and Keck/MOSFIRE. We infer gas-phase oxygen abundances ($Z_{\mathrm{g}}$; tracing O/H) via strong optical nebular lines a… ▽ More We present determinations of the gas-phase and stellar metallicities of a sample of 65 star-forming galaxies at $z \simeq 3.5$ using rest-frame far-ultraviolet (FUV) spectroscopy from the VANDELS survey in combination with follow-up rest-frame optical spectroscopy from VLT/KMOS and Keck/MOSFIRE. We infer gas-phase oxygen abundances ($Z_{\mathrm{g}}$; tracing O/H) via strong optical nebular lines and stellar iron abundances ($Z_{\star}$; tracing Fe/H) from full spectral fitting to the FUV continuum. Our sample spans the stellar mass range $8.5 < \mathrm{log}(M_{\star}/\mathrm{M}_{\odot}) < 10.5$ and shows clear evidence for both a stellar and gas-phase mass-metallicity relation (MZR). We find that our O and Fe abundance estimates both exhibit a similar mass-dependence, such that $\mathrm{Fe/H}\propto M_{\star}^{0.30\pm0.11}$ and $\mathrm{O/H}\propto M_{\star}^{0.32\pm0.09}$. At fixed $M_{\star}$ we find that, relative to their solar values, O abundances are systematically larger than Fe abundances (i.e., $α$-enhancement).We estimate an average enhancement of $\mathrm{(O/Fe)} = 2.65 \pm 0.16 \times \mathrm{(O/Fe)_\odot}$ which appears to be independent of $M_{\star}$. We employ analytic chemical evolution models to place a constraint on the strength of galactic-level outflows via the mass-outflow factor ($η$). We show that outflow efficiencies that scale as $η\propto M_{\star}^{-0.32}$ can simultaneously explain the functional form of of the stellar and gas-phase MZR, as well as the degree of $α$-enhancement at fixed Fe/H. Our results add further evidence to support a picture in which $α$-enhanced abundance ratios are ubiquitous in high-redshift star-forming galaxies, as expected for young systems whose interstellar medium is primarily enriched by core-collapse supernovae. △ Less

Submitted 10 July, 2024; v1 submitted 1 May, 2024; originally announced May 2024.

Comments: 18 pages, 10 figures, 4 tables, accepted in MNRAS

Quiescent or dusty? Unveiling the nature of extremely red galaxies at $z>3$

Authors: L. Barrufet, P. Oesch, R. Marques-Chaves, K. Arellano-Cordova, J. F. W. Baggen, A. C. Carnall, F. Cullen, J. S. Dunlop, R. Gottumukkala, Y. Fudamoto, G. D. Illingworth, D. Magee, R. J. McLure, D. J. McLeod, M. J. Michałowski, M. Stefanon, P. G. van Dokkum, A. Weibel

Abstract: The advent of the JWST has revolutionised our understanding of high-redshift galaxies. In particular, the NIRCam instrument on-board JWST has revealed a population of Hubble Space Telescope (HST)-dark galaxies that had previously evaded optical detection, potentially due to significant dust obscuration, quiescence, or simply extreme redshift. Here, we present the first NIRSpec spectra of 23 HST-da… ▽ More The advent of the JWST has revolutionised our understanding of high-redshift galaxies. In particular, the NIRCam instrument on-board JWST has revealed a population of Hubble Space Telescope (HST)-dark galaxies that had previously evaded optical detection, potentially due to significant dust obscuration, quiescence, or simply extreme redshift. Here, we present the first NIRSpec spectra of 23 HST-dark galaxies ($\mathrm{H-F444W>1.75}$), unveiling their nature and physical properties. This sample includes both dusty and quiescent galaxies with spectroscopic data from NIRSpec/PRISM, providing accurate spectroscopic redshifts with $\mathrm{\overline{z}_{spec} = 4.1 \pm 0.7}$. The spectral features demonstrate that, while the majority of HST-dark galaxies are dusty, a substantial fraction, $\mathrm{13^{+9}_{-6} \%}$, are quiescent. For the dusty galaxies, we have quantified the dust attenuation using the Balmer decrement ($\mathrm{Hα/ Hβ}$), finding attenuations $\mathrm{A_{V} > 2\ mag}$. We find that HST-dark dusty galaxies are $\mathrm{Hα}$ emitters with equivalent widths spanning the range $\mathrm{ 68 A < EW_{Hα} < 550 A }$, indicative of a wide range of recent star-formation activity. Whether dusty or quiescent, we find that HST-dark galaxies are predominantly massive, with 85\% of the galaxies in the sample having masses $\mathrm{log(M_{*}/M_{\odot}) > 9.8}$. This pilot NIRSpec program reveals the diverse nature of HST-dark galaxies and highlights the effectiveness of NIRSpec/PRISM spectroscopic follow-up in distinguishing between dusty and quiescent galaxies and properly quantifying their physical properties. Upcoming research utilising higher-resolution NIRSpec data and combining JWST with ALMA observations will enhance our understanding of these enigmatic and challenging sources. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 13 pages, 11 figures

The Rise of Faint, Red AGN at $z>4$: A Sample of Little Red Dots in the JWST Extragalactic Legacy Fields

Authors: Dale D. Kocevski, Steven L. Finkelstein, Guillermo Barro, Anthony J. Taylor, Antonello Calabrò, Brivael Laloux, Johannes Buchner, Jonathan R. Trump, Gene C. K. Leung, Guang Yang, Mark Dickinson, Pablo G. Pérez-González, Fabio Pacucci, Kohei Inayoshi, Rachel S. Somerville, Elizabeth J. McGrath, Hollis B. Akins, Micaela B. Bagley, Laura Bisigello, Rebecca A. A. Bowler, Adam Carnall, Caitlin M. Casey, Yingjie Cheng, Nikko J. Cleri, Luca Costantin , et al. (32 additional authors not shown)

Abstract: We present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. These sources are likely heavily-reddened AGN that trace a previously-hidden phase of dust-obscured black hole growth in the early Universe. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifti… ▽ More We present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. These sources are likely heavily-reddened AGN that trace a previously-hidden phase of dust-obscured black hole growth in the early Universe. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifting bandpasses to sample the same rest-frame emission blueward and redward of the Balmer break. This approach allows us to identify LRDs over a wider redshift range and is less susceptible to contamination from galaxies with strong breaks that otherwise lack a rising red continuum. The redshift distribution of our sample increases at $z<8$ and then undergoes a rapid decline at $z\sim4.5$, which may tie the emergence, and obscuration, of these sources to the inside-out growth that galaxies experience during this epoch. We find that LRDs are 2-3 dex more numerous than bright quasars at $z\sim5-7$, but their number density is only 0.6-1 dex higher than X-ray and UV selected AGN at these redshifts. Within our sample, we have identified the first X-ray detected LRDs at $z=3.1$ and $z=4.66$. An X-ray spectral analysis confirms that these AGN are moderately obscured with $\log\,(N_{\rm H}/{\rm cm}^{2}$) of $23.3^{+0.4}_{-1.3}$ and $22.72^{+0.13}_{-0.16}$. Our analysis reveals that reddened AGN emission dominates their rest-optical light, while the rest-UV originates from their host galaxies. We also present NIRSpec follow-up spectroscopy of 17 LRDs that show broad emission lines consistent with AGN activity. The confirmed AGN fraction of our sample is $71\%$ for sources with F444W$<26.5$. In addition, we find three LRDs with narrow blue-shifted Balmer absorption features in their spectra, suggesting an outflow of high-density, low ionization gas from near the central engine of these faint, red AGN. △ Less

Submitted 19 April, 2024; v1 submitted 4 April, 2024; originally announced April 2024.

Comments: 23 pages, 17 figures, submitted to ApJ

Galaxy Build-up in the first 1.5 Gyr of Cosmic History: Insights from the Stellar Mass Function at $z\sim4-9$ from JWST NIRCam Observations

Authors: Andrea Weibel, Pascal A. Oesch, Laia Barrufet, Rashmi Gottumukkala, Richard S. Ellis, Paola Santini, John R. Weaver, Natalie Allen, Rychard Bouwens, Rebecca A. A. Bowler, Gabe Brammer, Adam C. Carnall, Fergus Cullen, Pratika Dayal, Callum T. Donnan, James S. Dunlop, Mauro Giavalisco, Norman A. Grogin, Garth D. Illingworth, Anton M. Koekemoer, Ivo Labbe, Danilo Marchesini, Derek J. McLeod, Ross J. McLure, Rohan P. Naidu , et al. (4 additional authors not shown)

Abstract: Combining the public JWST/NIRCam imaging programs CEERS, PRIMER and JADES, spanning a total area of $\sim$500 arcmin$^2$, we obtain a sample of $>$30,000 galaxies at $z\sim4-9$ that allows us to perform a complete, rest-optical selected census of the galaxy population at $z>3$. Comparing the stellar mass $M_*$ and the UV-slope $β$ distributions between JWST- and HST-selected samples, we generally… ▽ More Combining the public JWST/NIRCam imaging programs CEERS, PRIMER and JADES, spanning a total area of $\sim$500 arcmin$^2$, we obtain a sample of $>$30,000 galaxies at $z\sim4-9$ that allows us to perform a complete, rest-optical selected census of the galaxy population at $z>3$. Comparing the stellar mass $M_*$ and the UV-slope $β$ distributions between JWST- and HST-selected samples, we generally find very good agreement and no significant biases. Nevertheless, JWST enables us to probe a small population of UV-red galaxies that was missing from previous HST-based LBG samples. We measure galaxy stellar mass functions (SMFs) at $z\sim4-9$ and show that they are broadly consistent with existing literature results. However, UV-red galaxies dominate the high-mass end of the SMF at least out to $z\sim6$. In particular the most massive galaxies typically show very red colors between $λ_{obs}\sim1.5μ$m and $\sim4.5μ$m, and thus JWST's unprecedented resolution and sensitivity at these wavelengths yields more accurate constraints on their abundance and masses. The implied redshift evolution of the high-mass end of the SMF suggests a rapid build-up of massive dust-obscured as well as quiescent galaxies from $z\sim6$ to $z\sim4$ as well as an enhanced efficiency of star formation towards earlier times ($z\gtrsim6$). We find the SMFs to be steep over the entire redshift range, and slightly steepening with redshift from $z\sim 4-6$, reaching values of $\approx-2$ at $z\gtrsim6$. Finally, we show that the galaxy mass density grows by a factor $\sim20\times$ in the $\sim1$ Gyr of cosmic time from $z\sim9$ to $z\sim4$. Our results emphasize the importance of rest-frame optically-selected samples in inferring accurate distributions of physical properties and studying the mass build-up of galaxies in the first 1.5 Gyr of cosmic history. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 29 pages, 17 figures, submitted to MNRAS

Charting the main sequence of star-forming galaxies out to redshifts z~7

Authors: M. P. Koprowski, J. V. Wijesekera, J. S. Dunlop, D. J. McLeod, M. J. Michałowski, K. Lisiecki, R. J. McLure

Abstract: We present a new determination of the star-forming main sequence (MS), obtained through stacking 100k K-band-selected galaxies in the far-IR Herschel and James Clerk Maxwell Telescope (JCMT) imaging. By fitting the dust emission curve to the stacked far-IR photometry, we derive the IR luminosities (LIR) and, hence, star formation rates (SFR) out to z~7. The functional form of the MS is found, with… ▽ More We present a new determination of the star-forming main sequence (MS), obtained through stacking 100k K-band-selected galaxies in the far-IR Herschel and James Clerk Maxwell Telescope (JCMT) imaging. By fitting the dust emission curve to the stacked far-IR photometry, we derive the IR luminosities (LIR) and, hence, star formation rates (SFR) out to z~7. The functional form of the MS is found, with the linear SFR-M* relation that flattens at high stellar masses and the normalization that increases exponentially with redshift. We derive the corresponding redshift evolution of the specific star formation rate (sSFR) and compare our findings with the recent literature. We find our MS to be exhibiting slightly lower normalization at z<=2 and to flatten at larger stellar masses at high redshifts. By deriving the relationship between the peak dust temperature (Td) and redshift, where Td increases linearly from ~20K at z=0.5 to ~50 K at z=6, we conclude that the apparent inconsistencies in the shapes of the MS are most likely caused by the different dust temperatures assumed when deriving SFRs in the absence of far-IR data. Finally, we investigate the derived shape of the star-forming MS by simulating the time evolution of the observed galaxy stellar mass function (GSMF). While the simulated GSMF is in good agreement with the observed one, some inconsistencies persist. In particular, we find the simulated GSMF to be somewhat overpredicting the number density of low-mass galaxies at z>2. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 16 pages, 11 figures, 2 tables, submitted to A&A

JWST PRIMER: A new multi-field determination of the evolving galaxy UV luminosity function at redshifts $\mathbf{z \simeq 9-15}$

Authors: C. T. Donnan, R. J. McLure, J. S. Dunlop, D. J. McLeod, D. Magee, K. Z. Arellano-Córdova, L. Barrufet, R. Begley, R. A. A. Bowler, A. C. Carnall, F. Cullen, R. S. Ellis, A. Fontana, G. D. Illingworth, N. A. Grogin, M. L. Hamadouche, A. M. Koekemoer, F. -Y. Liu, C. Mason, P. Santini, T. M. Stanton

Abstract: We present a new determination of the evolving galaxy UV luminosity function (LF) over the redshift range $8.5<z<15.5$ using a combination of several major Cycle-1 JWST imaging programmes - PRIMER, JADES and NGDEEP. This multi-field approach yields a total of $\simeq370$ sq. arcmin of JWST/NIRCam imaging, reaching (5-$σ$) depths of $\simeq30$ AB mag in the deepest regions. We select a sample of 25… ▽ More We present a new determination of the evolving galaxy UV luminosity function (LF) over the redshift range $8.5<z<15.5$ using a combination of several major Cycle-1 JWST imaging programmes - PRIMER, JADES and NGDEEP. This multi-field approach yields a total of $\simeq370$ sq. arcmin of JWST/NIRCam imaging, reaching (5-$σ$) depths of $\simeq30$ AB mag in the deepest regions. We select a sample of 2548 galaxies with a significant probability of lying at high redshift ($p(z>8.5)>0.05$) to undertake a statistical calculation of the UV LF. Our new measurements span $\simeq4$ magnitudes in UV luminosity at $z=9-12.5$, placing new constraints on both the shape and evolution of the LF at early times. Our measurements yield a new estimate of the early evolution of cosmic star-formation rate density ($ρ_{\rm{SFR}}$) confirming the gradual decline deduced from early JWST studies, at least out to $z \simeq 12$. Finally we show that the observed early evolution of the galaxy UV LF (and $ρ_{\rm{SFR}}$) can be reproduced in a ${\rm Λ}$CDM Universe, with no change in dust properties or star-formation efficiency required out to $z \simeq 12$. Instead, a progressive trend towards younger stellar population ages can reproduce the observations, and the typical ages required at $z \simeq$ 8, 9, 10, and 11 all converge on $\simeq 380-330$ Myr after the Big Bang, indicative of a rapid emergence of early galaxies at $z \simeq 12 - 13$. This is consistent with the first indications of a steeper drop-off in $ρ_{\rm{SFR}}$ we find beyond $z \simeq 13$, possibly reflecting the rapid evolution of the halo mass function at earlier times. △ Less

Submitted 16 May, 2024; v1 submitted 5 March, 2024; originally announced March 2024.

Comments: 15 pages, 9 figures, submitted to MNRAS. v2: updated in response to reviewer comments

The sizes of bright Lyman-break galaxies at $z\simeq3-5$ with JWST PRIMER

Authors: R. G. Varadaraj, R. A. A. Bowler, M. J. Jarvis, N. J. Adams, N. Choustikov, A. M. Koekemoer, A. C. Carnall, D. J. McLeod, J. S. Dunlop, C. T. Donnan, N. A. Grogin

Abstract: We use data from the JWST Public Release IMaging for Extragalactic Research (PRIMER) survey to measure the size scaling relations of 1668 rest-frame UV-bright Lyman-break galaxies (LBGs) at $z=3-5$ with stellar masses $\mathrm{log}_{10}(M_{\star}/M_{\odot}) > 9$. The sample was selected from seeing-dominated ground-based data, presenting an unbiased sampling of the morphology and size distribution… ▽ More We use data from the JWST Public Release IMaging for Extragalactic Research (PRIMER) survey to measure the size scaling relations of 1668 rest-frame UV-bright Lyman-break galaxies (LBGs) at $z=3-5$ with stellar masses $\mathrm{log}_{10}(M_{\star}/M_{\odot}) > 9$. The sample was selected from seeing-dominated ground-based data, presenting an unbiased sampling of the morphology and size distributions of luminous sources. We fit Sérsic profiles to eight NIRCam bands and also measure a non-parametric half-light radius. We find that the size distributions with both measurements are well-fit by a log-normal distribution at all redshifts, consistent with disk formation models where size is governed by host dark-matter halo angular momentum. We find a size-redshift evolution of $R_{e} = 3.51(1+z)^{-0.60\pm0.22}$ kpc, in agreement with JWST studies. When considering the typical (modal) size over $z=3-5$, we find little evolution with bright LBGs remaining compact at $R_{e}\simeq0.7-0.9$ kpc. Simultaneously, we find evidence for a build-up of large ($R_{e} > 2$ kpc) galaxies by $z=3$. We find some evidence for a negatively sloped size-mass relation at $z=5$ when Sérsic profiles are used to fit the data in F200W. The intrinsic scatter in our size-mass relations increases at higher redshifts. Additionally, measurements probing the rest-UV (F200W) show larger intrinsic scatter than those probing the rest-optical (F356W). Finally, we leverage rest-UV and rest-optical photometry to show that disky galaxies are well established by $z=5$, but are beginning to undergo dissipative processes, such as mergers, by $z=3$. The agreement of our size-mass and size-luminosity relations with simulations provides tentative evidence for centrally concentrated star formation at high-redshift. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 15 pages, 9 figures, 3 tables. Submitted to MNRAS

The ultraviolet continuum slopes of high-redshift galaxies: evidence for the emergence of dust-free stellar populations at z > 10

Authors: F. Cullen, D. J. McLeod, R. J. McLure, J. S. Dunlop, C. T. Donnan, A. C. Carnall, L. C. Keating, D. Magee, K. Z. Arellano-Cordova, R. A. A. Bowler, R. Begley, S. R. Flury, M. L. Hamadouche, T. M. Stanton

Abstract: We present an analysis of the ultraviolet (UV) continuum slopes ($β$) for a sample of $172$ galaxy candidates at $8 < z_{\mathrm{phot}} < 16$ selected from a combination of JWST NIRCam imaging and COSMOS/UltraVISTA ground-based near-infrared imaging. Focusing primarily on a new sample of $121$ galaxies at $\langle z \rangle \simeq 11$ selected from $\simeq 320$ arcmin$^2$ of public JWST imaging da… ▽ More We present an analysis of the ultraviolet (UV) continuum slopes ($β$) for a sample of $172$ galaxy candidates at $8 < z_{\mathrm{phot}} < 16$ selected from a combination of JWST NIRCam imaging and COSMOS/UltraVISTA ground-based near-infrared imaging. Focusing primarily on a new sample of $121$ galaxies at $\langle z \rangle \simeq 11$ selected from $\simeq 320$ arcmin$^2$ of public JWST imaging data across $15$ independent data sets, we investigate the evolution of $β$ in the galaxy population at $z \geq 9$. We find a significant trend between $β$ and redshift, with the inverse-variance weighted mean UV slope evolving from $\langle β\rangle = -2.17 \pm 0.06$ at $z = 9.5$ to $\langle β\rangle = -2.59 \pm 0.06$ at $z = 11.5$. Based on a comparison with stellar population models including nebular continuum emission, we find that at $z>10.5$ the average UV continuum slope is consistent with the intrinsic blue limit of dust-free stellar populations $(β_{\mathrm{int}} \simeq -2.6)$. These results suggest that the moderately dust-reddened galaxy population at $z < 10$ was essentially unattenuated at $z \simeq 11$. The extremely blue galaxies being uncovered at $z>10$ place important constraints on dust attenuation in galaxies in the early Universe, and imply that the already observed galaxy population is likely supplying an ionising photon budget capable of maintaining ionised IGM fractions of $\gtrsim 5$ per cent at $z\simeq11$. △ Less

Submitted 7 May, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: 18 pages, 12 figures, MNRAS accepted

The VANDELS ESO public spectroscopic survey: The spectroscopic measurements catalogue

Authors: M. Talia, C. Schreiber, B. Garilli, L. Pentericci, L. Pozzetti, G. Zamorani, F. Cullen, M. Moresco, A. Calabrò, M. Castellano, J. P. U. Fynbo, L. Guaita, F. Marchi, S. Mascia, R. McLure, M. Mignoli, E. Pompei, E. Vanzella, A. Bongiorno, G. Vietri, R. O. Amorín, M. Bolzonella, A. C. Carnall, A. Cimatti, G. Cresci , et al. (14 additional authors not shown)

Abstract: VANDELS is a deep spectroscopic survey, performed with the VIMOS instrument at VLT, aimed at studying in detail the physical properties of high-redshift galaxies. VANDELS targeted about 2100 sources at 1<z<6.5 in the CANDELS Chandra Deep-Field South (CDFS) and Ultra-Deep Survey (UDS) fields. In this paper we present the public release of the spectroscopic measurement catalogues from this survey, f… ▽ More VANDELS is a deep spectroscopic survey, performed with the VIMOS instrument at VLT, aimed at studying in detail the physical properties of high-redshift galaxies. VANDELS targeted about 2100 sources at 1<z<6.5 in the CANDELS Chandra Deep-Field South (CDFS) and Ultra-Deep Survey (UDS) fields. In this paper we present the public release of the spectroscopic measurement catalogues from this survey, featuring emission and absorption line centroids, fluxes, and rest-frame equivalent widths obtained through a Gaussian fit, as well as a number of atomic and molecular indices (e.g. Lick) and continuum breaks (e.g. D4000), and including a correction to be applied to the error spectra. We describe the measurement methods and the validation of the codes that were used. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: 12 pages, 9 figures; accepted for publication in Astronomy & Astrophysics. The measurement catalogues are accessible through the survey database (http://vandels.inaf.it) where all information can be queried interactively, and at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/cgi-bin/qcat?J/A+A/

Journal ref: A&A 678, A25 (2023)

An optimal ALMA image of the Hubble Ultra Deep Field in the era of JWST: obscured star formation and the cosmic far-infrared background

Authors: Ryley Hill, Douglas Scott, Derek J. McLeod, Ross J. McLure, Scott C. Chapman, James S. Dunlop

Abstract: We combine archival ALMA data targeting the Hubble Ultra Deep Field (HUDF) to produce the deepest currently attainable 1-mm maps of this key region. Our deepest map covers 4.2arcmin^2, with a beamsize of 1.49''x1.07'' at an effective frequency of 243GHz (1.23mm). It reaches an rms of 4.6uJy/beam, with 1.5arcmin^2 below 9.0uJy/beam, an improvement of >5% (and up to 50% in some regions) over the bes… ▽ More We combine archival ALMA data targeting the Hubble Ultra Deep Field (HUDF) to produce the deepest currently attainable 1-mm maps of this key region. Our deepest map covers 4.2arcmin^2, with a beamsize of 1.49''x1.07'' at an effective frequency of 243GHz (1.23mm). It reaches an rms of 4.6uJy/beam, with 1.5arcmin^2 below 9.0uJy/beam, an improvement of >5% (and up to 50% in some regions) over the best previous map. We also make a wider, shallower map, covering 25.4arcmin^2. We detect 45 galaxies in the deep map down to 3.6sigma, 10 more than previously detected, and 39 of these galaxies have JWST counterparts. A stacking analysis on the positions of ALMA-undetected JWST galaxies with z<4 and stellar masses from 10^8.4 to 10^10.4 M_sun yields 10% more signal compared to previous stacking analyses, and we find that detected sources plus stacking contribute (10.0+/-0.5)Jy/deg^2 to the cosmic infrared background (CIB) at 1.23mm. Although this is short of the (uncertain) background level of about 20Jy/deg^2, we show that our measurement is consistent with the background if the HUDF is a mild (~2sigma) negative CIB fluctuation, and that the contribution from faint undetected objects is small and converging. In particular, we predict that the field contains about 60 additional 15uJy galaxies, and over 300 galaxies at the few uJy level. This suggests that JWST has detected essentially all of the galaxies that contribute to the CIB, as anticipated from the strong correlation between galaxy stellar mass and obscured star formation. △ Less

Submitted 1 February, 2024; v1 submitted 19 September, 2023; originally announced September 2023.

Comments: Published in MNRAS. All of the combined ALMA maps described in this paper are available at https://doi.org/10.5683/SP3/YWBVWH

Connecting the escape fraction of Lyman-alpha and Lyman-continuum photons in star-forming galaxies at $\mathbf{z\simeq 4-5}$

Authors: R. Begley, F. Cullen, R. J. McLure, A. E. Shapley, J. S. Dunlop, A. C. Carnall, D. J. McLeod, C. T. Donnan, M. L. Hamadouche, T. M. Stanton

Abstract: We present a study of the connection between the escape fraction of Lyman-alpha (Ly$α$) and Lyman-continuum (LyC) photons within a sample of N=152 star-forming galaxies selected from the VANDELS survey at $3.85<z_{spec}<4.95$. By combining measurements of H$α$ equivalent width $(W_λ(\rm{Hα}))$ derived from broad-band photometry with Ly$α$ equivalent width $(W_λ(Lyα))$ measurements from VANDELS spe… ▽ More We present a study of the connection between the escape fraction of Lyman-alpha (Ly$α$) and Lyman-continuum (LyC) photons within a sample of N=152 star-forming galaxies selected from the VANDELS survey at $3.85<z_{spec}<4.95$. By combining measurements of H$α$ equivalent width $(W_λ(\rm{Hα}))$ derived from broad-band photometry with Ly$α$ equivalent width $(W_λ(Lyα))$ measurements from VANDELS spectra, we individually estimate $f_{\rm{esc}}^{Lyα}$ for our full sample. In agreement with previous studies, we find a positive correlation between $W_λ(Lyα)$ and $f_{\rm{esc}}^{Lyα}$, increasing from $f_{\rm{esc}}^{Lyα}\simeq0.04$ at $W_λ(Lyα)=10$Å to $f_{\rm{esc}}^{Lyα}\simeq0.1$ at $W_λ(Lyα)=25$Å. For the first time at $z\sim4-5$, we investigate the relationship between $f_{\rm{esc}}^{Lyα}$ and $f_{\rm{esc}}^{\rm{LyC}}$ using $f_{\rm{esc}}^{\rm{LyC}}$ estimates derived using the equivalent widths of low-ionization, FUV absorption lines in composite VANDELS spectra. Our results indicate that $f_{\rm{esc}}^{\rm{LyC}}$ rises monotonically with $f_{\rm{esc}}^{Lyα}$, following the relation $f_{\rm{esc}}^{\rm{LyC}}\simeq 0.15^{+0.06}_{-0.04}f_{\rm{esc}}^{Lyα}$. Based on composite spectra of sub-samples with roughly constant $W_λ(Lyα)$, but very different $f_{\rm{esc}}^{Lyα}$, we show that the $f_{\rm{esc}}^{\rm{LyC}}-f_{\rm{esc}}^{Lyα}$ correlation is not driven by a secondary correlation between $f_{\rm{esc}}^{Lyα}$and $W_λ(Lyα)$. The $f_{\rm{esc}}^{\rm{LyC}}-f_{\rm{esc}}^{Lyα}$ correlation is in good qualitative agreement with theoretical predictions and provides further evidence that estimates of $f_{\rm{esc}}^{\rm{LyC}}$ within the Epoch of Reionization should be based on proxies sensitive to neutral gas density/geometry and dust attenuation. △ Less

Submitted 2 November, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Comments: 12 pages, 7 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

The galaxy UV luminosity function at $\mathbf{z \simeq 11}$ from a suite of public JWST ERS, ERO and Cycle-1 programs

Authors: D. J. McLeod, C. T. Donnan, R. J. McLure, J. S. Dunlop, D. Magee, R. Begley, A. C. Carnall, F. Cullen, R. S. Ellis, M. L. Hamadouche, T. M. Stanton

Abstract: We present a new determination of the evolving galaxy UV luminosity function (LF) over the redshift range $9.5<z<12.5$ based on a wide-area ($>250$ arcmin$^2$) data set of JWST NIRCam near-infrared imaging assembled from thirteen public JWST surveys. Our relatively large-area search allows us to uncover a sample of 61 robust $z>9.5$ candidates detected at $\geq 8σ$, and hence place new constraints… ▽ More We present a new determination of the evolving galaxy UV luminosity function (LF) over the redshift range $9.5<z<12.5$ based on a wide-area ($>250$ arcmin$^2$) data set of JWST NIRCam near-infrared imaging assembled from thirteen public JWST surveys. Our relatively large-area search allows us to uncover a sample of 61 robust $z>9.5$ candidates detected at $\geq 8σ$, and hence place new constraints on the intermediate-to-bright end of the UV LF. When combined with our previous JWST+UltraVISTA results, this allows us to measure the form of the LF over a luminosity range corresponding to four magnitudes ($M_{1500}$). At these early times we find that the galaxy UV LF is best described by a double power-law function, consistent with results obtained from recent ground-based and early JWST studies at similar redshifts. Our measurements provide further evidence for a relative lack of evolution at the bright-end of the UV LF at $z=9-11$, but do favour a steep faint-end slope ($α\leq-2$). The luminosity-weighted integral of our evolving UV LF provides further evidence for a gradual, smooth (exponential) decline in co-moving star-formation rate density ($ρ_{\mathrm{SFR}}$) at least out to $z\simeq12$, with our determination of $ρ_{\mathrm{SFR}}(z=11)$ lying significantly above the predictions of many theoretical models of galaxy evolution. △ Less

Submitted 8 November, 2023; v1 submitted 27 April, 2023; originally announced April 2023.

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

Confirmation and refutation of very luminous galaxies in the early universe

Authors: Pablo Arrabal Haro, Mark Dickinson, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Callum T. Donnan, Denis Burgarella, Adam Carnall, Fergus Cullen, James S. Dunlop, Vital Fernández, Seiji Fujimoto, Intae Jung, Melanie Krips, Rebecca L. Larson, Casey Papovich, Pablo G. Pérez-González, Ricardo O. Amorín, Micaela B. Bagley, Véronique Buat, Caitlin M. Casey, Katherine Chworowsky, Seth H. Cohen, Henry C. Ferguson, Mauro Giavalisco, Marc Huertas-Company , et al. (12 additional authors not shown)

Abstract: During the first 500 million years of cosmic history, the first stars and galaxies formed, seeding the Universe with heavy elements and eventually reionizing the intergalactic medium. Observations with JWST have uncovered a surprisingly high abundance of candidates for early star-forming galaxies, with distances (redshifts, $z$), estimated from multi-band photometry, as large as $z\approx 16$, far… ▽ More During the first 500 million years of cosmic history, the first stars and galaxies formed, seeding the Universe with heavy elements and eventually reionizing the intergalactic medium. Observations with JWST have uncovered a surprisingly high abundance of candidates for early star-forming galaxies, with distances (redshifts, $z$), estimated from multi-band photometry, as large as $z\approx 16$, far beyond pre-JWST limits. While generally robust, such photometric redshifts can suffer from degeneracies and occasionally catastrophic errors. Spectroscopic measurement is required to validate these sources and to reliably quantify physical properties that can constrain galaxy formation models and cosmology. Here we present JWST spectroscopy that confirms redshifts for two very luminous galaxies with $z > 11$, but also demonstrates that another candidate with suggested $z\approx 16$ instead has $z = 4.9$, with an unusual combination of nebular line emission and dust reddening that mimics the colors expected for much more distant objects. These results reinforce evidence for the early, rapid formation of remarkably luminous galaxies, while also highlighting the necessity of spectroscopic verification. The large abundance of bright, early galaxies may indicate shortcomings in current galaxy formation models, or deviation from physical properties (such as the stellar initial mass function) that are generally believed to hold at later times. △ Less

Submitted 15 August, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: Accelerated Article Preview (AAP) version published in Nature

The connection between stellar mass, age and quenching timescale in massive quiescent galaxies at $z \simeq 1$

Authors: M. L. Hamadouche, A. C. Carnall, R. J. McLure, J. S. Dunlop, R. Begley, F. Cullen, D. J. McLeod, C. T. Donnan, T. M. Stanton

Abstract: We present a spectro-photometric study of a mass-complete sample of quiescent galaxies at $1.0 < z < 1.3$ with $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}}) \geq 10.3$ drawn from the VANDELS survey, exploring the relationship between stellar mass, age and star-formation history. Within our sample of 114 galaxies, we derive a stellar-mass vs stellar-age relation with a slope of… ▽ More We present a spectro-photometric study of a mass-complete sample of quiescent galaxies at $1.0 < z < 1.3$ with $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}}) \geq 10.3$ drawn from the VANDELS survey, exploring the relationship between stellar mass, age and star-formation history. Within our sample of 114 galaxies, we derive a stellar-mass vs stellar-age relation with a slope of $1.20^{+0.28}_{-0.27}$ Gyr per decade in stellar mass. When combined with recent literature results, we find evidence that the slope of this relation remains consistent over the redshift interval $0<z<4$. The galaxies within the VANDELS quiescent display a wide range of star-formation histories, with a mean star-formation timescale of $1.5\pm{0.1}$ Gyr and a mean quenching timescale of $1.4\pm{0.1}$ Gyr. We also find a large scatter in the quenching timescales of the VANDELS quiescent galaxies, in agreement with previous evidence that galaxies at $z \sim 1$ cease star formation via multiple mechanisms. We then focus on the oldest galaxies in our sample, finding that the number density of galaxies that quenched before $z = 3$ with stellar masses $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}}) \geq 10.6$ is $ 1.12_{-0.72}^{+1.47} \times 10^{-5} \ \mathrm{Mpc}^{-3}$. Although uncertain, this estimate is in good agreement with the latest observational results at $3<z<4$, tentatively suggesting that neither rejuvenation nor merger events are playing a major role in the evolution of the oldest massive quiescent galaxies within the redshift interval $1<z<3$. △ Less

Submitted 10 March, 2023; originally announced March 2023.

Comments: Accepted for publication in MNRAS, 11 pages, 6 figures

A massive quiescent galaxy at redshift 4.658

Authors: A. C. Carnall, R. J. McLure, J. S. Dunlop, D. J. McLeod, V. Wild, F. Cullen, D. Magee, R. Begley, A. Cimatti, C. T. Donnan, M. L. Hamadouche, S. M. Jewell, S. Walker

Abstract: The extremely rapid assembly of the earliest galaxies during the first billion years of cosmic history is a major challenge for our understanding of galaxy formation physics. The advent of JWST has exacerbated this issue by confirming the existence of galaxies in significant numbers as early as the first few hundred million years. Perhaps even more surprisingly, in some galaxies, this initial high… ▽ More The extremely rapid assembly of the earliest galaxies during the first billion years of cosmic history is a major challenge for our understanding of galaxy formation physics. The advent of JWST has exacerbated this issue by confirming the existence of galaxies in significant numbers as early as the first few hundred million years. Perhaps even more surprisingly, in some galaxies, this initial highly efficient star formation rapidly shuts down, or quenches, giving rise to massive quiescent galaxies as little as 1.5 billion years after the Big Bang. However, due to their faintness and red colour, it has proven extremely challenging to learn about these extreme quiescent galaxies, or to confirm whether any exist at earlier times. Here we report the spectroscopic confirmation of a massive quiescent galaxy, GS-9209, at redshift $z=4.658$, just 1.25 billion years after the Big Bang, using JWST NIRSpec. From these data we infer a stellar mass of $M_* = 3.8\pm0.2\times10^{10}\ M_\odot$, which formed over a $\simeq200$ Myr period before this galaxy quenched its star formation activity at $z=6.5^{+0.2}_{-0.5}$, when the Universe was $\simeq800$ million years old. Based on the presence of broad H$α$ in the spectrum and a high narrow-line [NII]/H$α$ ratio, we infer the presence of an accreting supermassive black hole, with a mass of $M_\bullet = 5\pm1\times10^{8}\ M_\odot$. This large black hole mass relative to the stellar mass suggests that active galactic nucleus (AGN) feedback may have been responsible for quenching this galaxy. GS-9209 is also extremely compact, with an effective radius, $r_e=215\pm20$ parsecs. This galaxy is both a likely descendent of the highest-redshift submillimetre galaxies and quasars, and a likely progenitor for the dense, ancient cores of the most massive local galaxies. △ Less

Submitted 26 May, 2023; v1 submitted 26 January, 2023; originally announced January 2023.

Comments: 23 pages, 6 figures, 2 tables, accepted for publication in Nature

The abundance of $z \gtrsim 10$ galaxy candidates in the HUDF using deep JWST NIRCam medium-band imaging

Authors: C. T. Donnan, D. J. McLeod, R. J. McLure, J. S. Dunlop, A. C. Carnall, F. Cullen, D. Magee

Abstract: We utilise JWST NIRCam medium-band imaging to search for extreme redshift ($z \geq 9.5$) galaxy candidates in the Hubble Ultra Deep Field (HUDF) and the additional pointing within the GOODS-South field provided by the second NIRCam module. Our search reveals 6 robust candidates, 3 of which have recently been spectroscopically confirmed. One of these 3 is the previously controversial $z \simeq 12$… ▽ More We utilise JWST NIRCam medium-band imaging to search for extreme redshift ($z \geq 9.5$) galaxy candidates in the Hubble Ultra Deep Field (HUDF) and the additional pointing within the GOODS-South field provided by the second NIRCam module. Our search reveals 6 robust candidates, 3 of which have recently been spectroscopically confirmed. One of these 3 is the previously controversial $z \simeq 12$ galaxy candidate UDF-22980 which is now detected in five JWST NIRCam medium-band filters (F182M, F210M, F430M, F460M and F480M), efficiently excluding alternative low-redshift solutions and allowing us to now report a secure photometric redshift of $z = 11.6 \pm 0.2$. We also detect 2 galaxies at $z \geq 12.5$ including a newly-detected candidate in the imaging provided by the second NIRCam module (south-west of the HUDF) at $z = 12.6 \pm 0.6$. We determine the physical properties of the 6 galaxies by fitting the 14-band photometry with Bagpipes. We find stellar masses of $\log(M_{\star}/{\rm {M_{\odot}}}) \simeq 7.5 - 8.7$ and star-formation rates of $\log(\rm{SFR}/M_{\odot}^{-1} \rm{yr}^{-1}) \simeq 0.3 - 5.0$. Despite the relatively small cosmological volume covered by the HUDF itself and the second NIRCam module imaging, we find that the existence of these galaxies is fully consistent with the latest measurements of both the UV luminosity function and cosmic star-formation rate density at $z\simeq11$, supporting a gradual steady decline in the cosmic star-formation rate density out to at least $z\simeq15$. △ Less

Submitted 8 February, 2023; v1 submitted 20 December, 2022; originally announced December 2022.

Comments: 8 pages, 4 figures, accepted for publication in MNRAS

The VANDELS survey: the ionizing properties of star-forming galaxies at $3 \leq z \leq 5$ using deep rest-frame ultraviolet spectroscopy

Authors: A. Saldana-Lopez, D. Schaerer, J. Chisholm, A. Calabrò, L. Pentericci, F. Cullen, A. Saxena, R. Amorín, A. C. Carnall, F. Fontanot, J. P. U. Fynbo, L. Guaita, N. P. Hathi, P. Hibon Z. Ji D. J. McLeod, E. Pompei, G. Zamorani

Abstract: To better understand the ionizing properties of galaxies in the EoR, we investigate deep, rest-frame ultraviolet (UV) spectra of $\simeq 500$ star-forming galaxies at $3 \leq z \leq 5$ selected from the public ESO-VANDELS spectroscopic survey. The absolute ionizing photon escape fraction ($f_{\rm esc}^{\rm abs}$) is derived by combining absorption line measurements with estimates of the UV attenua… ▽ More To better understand the ionizing properties of galaxies in the EoR, we investigate deep, rest-frame ultraviolet (UV) spectra of $\simeq 500$ star-forming galaxies at $3 \leq z \leq 5$ selected from the public ESO-VANDELS spectroscopic survey. The absolute ionizing photon escape fraction ($f_{\rm esc}^{\rm abs}$) is derived by combining absorption line measurements with estimates of the UV attenuation. The ionizing production efficiency ($ξ_{ion}$) is calculated by fitting the far-UV (FUV) stellar continuum of the VANDELS galaxies. We find that the $f_{\rm esc}^{\rm abs}$ and $ξ_{ion}$ parameters increase towards low-mass, blue UV-continuum slopes and strong Ly$α$ emitting galaxies, and both are just slightly higher-than-average for the UV-faintest galaxies in the sample. Potential Lyman Continuum Emitters (LCEs) and selected Lyman Alpha Emitters (LAEs) show systematically higher $ξ_{ion}$ ($\log ξ_{ion}$ (Hz\erg) $\approx 25.38, 25.41$) than non-LCEs and non-LAEs ($\log ξ_{ion}$ (Hz\erg) $\approx 25.18, 25.14$) at similar UV magnitudes. This indicates very young underlying stellar populations ($\approx 10~{\rm Myr}$) at relatively low metallicities ($\approx 0.2~{\rm Z_{\odot}}$). The FUV non-ionizing spectra of potential LCEs is characterized by very blue UV slopes ($\leq -2$), enhanced Ly$α$ emission ($\leq -25$A), strong UV nebular lines (e.g., high CIV1550/CIII]1908 $\geq 0.75$ ratios), and weak absorption lines ($\leq 1$A). The latter suggests the existence of low gas-column-density channels in the interstellar medium which enables the escape of ionizing photons. By comparing our VANDELS results against other surveys in the literature, our findings imply that the ionizing budget in the EoR was likely dominated by UV-faint, low-mass and dustless galaxies. △ Less

Submitted 10 May, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 25 pages, 18 figures, 2 tables; accepted for publication in MNRAS on April 24, 2023

The ultraviolet continuum slopes ($\mathbfβ$) of galaxies at $\mathbf{z\simeq8-16}$ from JWST and ground-based near-infrared imaging

Authors: F. Cullen, R. J. McLure, D. J. McLeod, J. S. Dunlop, C. T. Donnan, A. C. Carnall, R. A. A. Bowler, R. Begley, M. L. Hamadouche, T. M. Stanton

Abstract: We study the rest-frame ultraviolet (UV) continuum slopes ($β$) of galaxies at redshifts $8 < z < 16$, using a combination of JWST ERO and ERS NIRcam imaging and ground-based near-infrared imaging of the COSMOS field. The combination of JWST and ground-based imaging provides a wide baseline in both redshift and absolute UV magnitude ($-22.6 < M_{\rm UV} < 17.9$), sufficient to allow a meaningful c… ▽ More We study the rest-frame ultraviolet (UV) continuum slopes ($β$) of galaxies at redshifts $8 < z < 16$, using a combination of JWST ERO and ERS NIRcam imaging and ground-based near-infrared imaging of the COSMOS field. The combination of JWST and ground-based imaging provides a wide baseline in both redshift and absolute UV magnitude ($-22.6 < M_{\rm UV} < 17.9$), sufficient to allow a meaningful comparison to previous results at lower redshift. Using a power-law fitting technique, we find that our full sample (median $M_{\rm UV}=-19.3\pm 1.3$) returns an inverse-variance weighted mean value of $\langle β\rangle = -2.10 \pm 0.05$, with a corresponding median value of $β=-2.29\pm 0.09$. These values imply that the UV colours of galaxies at $z>8$ are, on average, no bluer than the bluest galaxies in the local Universe. Moreover, we find evidence for a $β-M_{\rm UV}$ relation, such that brighter UV galaxies display redder UV slopes ($\rm{d}β/ \rm{d} M_{\rm UV} = -0.17 \pm 0.05$). Comparing to results at lower redshift, we find that the slope of our $β-M_{\rm UV}$ relation is consistent with the slope observed at $z\simeq 5$ and that, at a given $M_{\rm UV}$, our $8<z<16$ galaxies are somewhat bluer than their $z\simeq 5$ counterparts, with an inverse-variance weighted mean offset of $\langle Δβ\rangle = -0.38 \pm 0.09$. We do not find strong evidence that any objects in our sample display ultra-blue UV continuum slopes (i.e., $β\lesssim-3$) that would require their UV emission to be dominated by ultra-young, dust-free stellar populations with high Lyman-continuum escape fractions. Comparing our results to the predictions of theoretical galaxy formation models, we find that the galaxies in our sample are consistent with the young, metal-poor and moderately dust-reddened galaxies expected at $z>8$. △ Less

Submitted 6 January, 2023; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: 9 pages, 5 figures, MNRAS accepted version

A surprising abundance of massive quiescent galaxies at 3 < z < 5 in the first data from JWST CEERS

Authors: A. C. Carnall, D. J. McLeod, R. J. McLure, J. S. Dunlop, R. Begley, F. Cullen, C. T. Donnan, M. L. Hamadouche, S. M. Jewell, E. W. Jones, C. L. Pollock, V. Wild

Abstract: We report a robust sample of 10 massive quiescent galaxies at redshift, $z > 3$, selected using the first data from the JWST CEERS programme. Three of these galaxies are at $4 < z < 5$, constituting the best evidence to date for quiescent galaxies significantly before $z=4$. These extreme galaxies have stellar masses in the range log$_{10}(M_*/$M$_\odot) = 10.1-11.1$, and formed the bulk of their… ▽ More We report a robust sample of 10 massive quiescent galaxies at redshift, $z > 3$, selected using the first data from the JWST CEERS programme. Three of these galaxies are at $4 < z < 5$, constituting the best evidence to date for quiescent galaxies significantly before $z=4$. These extreme galaxies have stellar masses in the range log$_{10}(M_*/$M$_\odot) = 10.1-11.1$, and formed the bulk of their mass around $z \simeq 10$, with two objects having star-formation histories that suggest they had already reached log$_{10}(M_*/$M$_\odot) > 10$ by $z\gtrsim8$. We report number densities for our sample, demonstrating that, based on the small area of JWST imaging so far available, previous work appears to have underestimated the number of quiescent galaxies at $3 < z < 4$ by a factor of $3-5$, due to a lack of ultra-deep imaging data at $λ>2\,μ$m. This result deepens the existing tension between observations and theoretical models, which already struggle to reproduce previous estimates of $z>3$ quiescent galaxy number densities. Upcoming wider-area JWST imaging surveys will provide larger samples of such galaxies and more-robust number densities, as well as providing opportunities to search for quiescent galaxies at $z>5$. The galaxies we report are excellent potential targets for JWST NIRSpec spectroscopy, which will be required to understand in detail their physical properties, providing deeper insights into the processes responsible for forming massive galaxies and quenching star formation during the first billion years. △ Less

Submitted 6 February, 2023; v1 submitted 1 August, 2022; originally announced August 2022.

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

The evolution of the galaxy UV luminosity function at redshifts z ~ 8-15 from deep JWST and ground-based near-infrared imaging

Authors: C. T. Donnan, D. J. McLeod, J. S. Dunlop, R. J. McLure, A. C. Carnall, R. Begley, F. Cullen, M. L. Hamadouche, R. A. A. Bowler, D. Magee, H. J. McCracken, B. Milvang-Jensen, A. Moneti, T. Targett

Abstract: We reduce and analyse the available James Webb Space Telescope (JWST) ERO and ERS NIRCam imaging (SMACS0723, GLASS, CEERS) in combination with the latest deep ground-based near-infrared imaging in the COSMOS field (provided by UltraVISTA DR5) to produce a new measurement of the evolving galaxy UV luminosity function (LF) over the redshift range $z = 8 - 15$. This yields a new estimate of the evolu… ▽ More We reduce and analyse the available James Webb Space Telescope (JWST) ERO and ERS NIRCam imaging (SMACS0723, GLASS, CEERS) in combination with the latest deep ground-based near-infrared imaging in the COSMOS field (provided by UltraVISTA DR5) to produce a new measurement of the evolving galaxy UV luminosity function (LF) over the redshift range $z = 8 - 15$. This yields a new estimate of the evolution of UV luminosity density ($ρ_{\rm UV}$), and hence cosmic star-formation rate density ($ρ_{\rm SFR}$) out to within $< 300$\, Myr of the Big Bang. Our results confirm that the high-redshift LF is best described by a double power-law (rather than a Schechter) function up to $z\sim10$, and that the LF and the resulting derived $ρ_{\rm UV}$ (and thus $ρ_{\rm SFR}$), continues to decline gradually and steadily up to $z\sim15$ (as anticipated from previous studies which analysed the pre-existing data in a consistent manner to this study). We provide details of the 61 high-redshift galaxy candidates, 47 of which are new, that have enabled this new analysis. Our sample contains 6 galaxies at $z \ge 12$, one of which appears to set a new redshift record as an apparently robust galaxy candidate at $z \simeq 16.4$, the properties of which we therefore consider in detail. The advances presented here emphasize the importance of achieving high dynamic range in studies of early galaxy evolution, and re-affirm the enormous potential of forthcoming larger JWST programmes to transform our understanding of the young Universe. △ Less

Submitted 24 November, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 28 pages, 4 figures in main manuscript, accepted for publication in MNRAS. Updated zero-point corrections noted in Appendix C

A first look at the SMACS0723 JWST ERO: spectroscopic redshifts, stellar masses and star-formation histories

Authors: A. C. Carnall, R. Begley, D. J. McLeod, M. L. Hamadouche, C. T. Donnan, R. J. McLure, J. S. Dunlop, B. Milvang-Jensen, C. L. Bondestam, F. Cullen, S. M. Jewell, C. L. Pollock

Abstract: We present a first-look analysis of the JWST ERO data in the SMACS J0723.3-7327 cluster field. We begin by reporting 10 new spectroscopic redshifts from $λ_\mathrm{obs}=1.8-5.2μ$m NIRSpec medium-resolution ($R=λ/Δλ= 1000$) data. These are determined via multiple high-SNR emission line detections, with 5 objects at $1 < z < 3$ displaying multiple rest-frame near-infrared Hydrogen Paschen lines, and… ▽ More We present a first-look analysis of the JWST ERO data in the SMACS J0723.3-7327 cluster field. We begin by reporting 10 new spectroscopic redshifts from $λ_\mathrm{obs}=1.8-5.2μ$m NIRSpec medium-resolution ($R=λ/Δλ= 1000$) data. These are determined via multiple high-SNR emission line detections, with 5 objects at $1 < z < 3$ displaying multiple rest-frame near-infrared Hydrogen Paschen lines, and 5 objects at $5 < z < 9$ displaying rest-frame optical Oxygen and Hydrogen Balmer lines. For the 5 higher-redshift galaxies we extract fluxes in 6 NIRCam bands spanning $λ_\mathrm{obs}=0.8-5μ$m and perform spectral energy distribution fitting, in combination with existing HST photometry. The $7 < z < 9$ objects exhibit a U-shaped pattern across the F277W, F356W and F444W bands, indicating a Balmer break seen in emission (Balmer jump) and high-equivalent-width [O\,\textsc{iii}] emission. This indicates an extremely young stellar population, with the bulk of the current mass having formed within the past 10 Myr. We report robust stellar masses and mean stellar ages from our spectral fitting, with the four $z > 6$ galaxies exhibiting low stellar masses from log$_{10}(M_*/$M$_\odot)=7.1-8.2$ and correspondingly young mean stellar ages of only a few Myr. This work highlights the critical importance of combining large upcoming NIRCam surveys with NIRSpec follow-up to measure the spectroscopic redshifts necessary to robustly constrain physical parameters. △ Less

Submitted 24 October, 2022; v1 submitted 18 July, 2022; originally announced July 2022.

Comments: 5 pages, 2 figures, 2 tables, accepted for publication in MNRAS Letters

The VANDELS survey: a measurement of the average Lyman-continuum escape fraction of star-forming galaxies at z=3.5

Authors: R. Begley, F. Cullen, R. J. McLure, J. S. Dunlop, A. Hall, A. C. Carnall, M. L. Hamadouche, D. J. McLeod, R. Amorín, A. Calabrò, A. Fontana, J. P. U. Fynbo, L. Guaita, N. P. Hathi, P. Hibon, Z. Ji, M. Llerena, L. Pentericci, A. Saldana-Lopez, D. Schaerer, M. Talia, E. Vanzella, G. Zamorani

Abstract: We present a study designed to measure the average LyC escape fraction ($\langle f_{\rm esc}\rangle$) of star-forming galaxies at z=3.5. We assemble a sample of 148 galaxies from the VANDELS survey at $3.35\leq z_{\rm spec}\leq3.95$, selected to minimize line-of-sight contamination of their photometry. For this sample, we use ultra-deep, ground-based, $U-$band imaging and HST $V-$band imaging to r… ▽ More We present a study designed to measure the average LyC escape fraction ($\langle f_{\rm esc}\rangle$) of star-forming galaxies at z=3.5. We assemble a sample of 148 galaxies from the VANDELS survey at $3.35\leq z_{\rm spec}\leq3.95$, selected to minimize line-of-sight contamination of their photometry. For this sample, we use ultra-deep, ground-based, $U-$band imaging and HST $V-$band imaging to robustly measure the distribution of $\mathcal{R_{\rm obs}}$ $=(L_{\rm LyC}/L_{\rm UV})_{\rm obs}$. We then model the distribution as a function of $\langle f_{\rm esc}\rangle$, carefully accounting for attenuation by dust, and the IGM (and CGM). A maximum likelihood fit to the $\mathcal{R_{\rm obs}}$ distribution returns a best-fitting value of $\langle f_{\rm esc}\rangle =0.07\pm0.02$, a result confirmed using an alternative Bayesian inference technique (both exclude $\langle f_{\rm esc}\rangle=0.0$ at $> 3σ$). By splitting our sample in two, we find evidence that $\langle f_{\rm esc}\rangle$ is positively correlated with Ly$α$ equivalent width, with high and low sub-samples returning best fits of $\langle f_{\rm esc}\rangle=0.12^{+0.06}_{-0.04}$ and $\langle f_{\rm esc} \rangle=0.02^{+0.02}_{-0.01}$, respectively. In contrast, we find evidence that $\langle f_{\rm esc}\rangle$ is anti-correlated with intrinsic UV luminosity and UV dust attenuation; with low UV luminosity and dust attenuation sub-samples returning best fits in the range $0.10 \leq \langle f_{\rm esc}\rangle \leq 0.22$. We do not find evidence for a clear correlation between $f_{\rm esc}$ and galaxy stellar mass, suggesting it is not a primary indicator of leakage. Although larger samples are needed to further explore these trends, they suggest that it is entirely plausible that the low dust and metallicity galaxies found at z > 6 will display the $\langle f_{\rm esc}\rangle\geq0.1$ required to drive reionization. △ Less

Submitted 21 April, 2022; v1 submitted 8 February, 2022; originally announced February 2022.

Comments: 15 pages (+ appendix), 12 figures, 1 table, Accepted for publication in Monthly Notices of the Royal Astronomical Society

A combined VANDELS and LEGA-C study: the evolution of quiescent galaxy size, stellar mass and age from $\mathbf{\textit{z} = 0.6}$ to $\mathbf{\textit{z} = 1.3}$

Authors: M. L. Hamadouche, A. C. Carnall, R. J. McLure, J. S. Dunlop, D. J. McLeod, F. Cullen, R. Begley, M. Bolzonella, F. Buitrago, M. Castellano, O. Cucciati, A. Fontana, A. Gargiulo, M. Moresco, L. Pozzetti, G. Zamorani

Abstract: We study the relationships between stellar mass, size and age within the quiescent population, using two mass-complete spectroscopic samples with $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}})>10.3$, taken from VANDELS at $1.0<z<1.3$, and LEGA-C at $0.6<z<0.8$. Using robust D$_{n}$4000 values, we demonstrate that the well-known 'downsizing' signature is already in place by $z\simeq1.1$, with D… ▽ More We study the relationships between stellar mass, size and age within the quiescent population, using two mass-complete spectroscopic samples with $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}})>10.3$, taken from VANDELS at $1.0<z<1.3$, and LEGA-C at $0.6<z<0.8$. Using robust D$_{n}$4000 values, we demonstrate that the well-known 'downsizing' signature is already in place by $z\simeq1.1$, with D$_{n}$4000 increasing by $\simeq0.1$ across a $\simeq$ 1 dex mass interval for both VANDELS and LEGA-C. We then proceed to investigate the evolution of the quiescent galaxy stellar mass-size relation from $z\simeq1.1$ to $z\simeq0.7$. We find the median size increases by a factor of $1.9\pm{0.1}$ at $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}})=10.5$, and see tentative evidence for flattening of the relation, finding slopes of $α=0.72\pm0.06$ and $α=$ $0.56\pm0.04$ for VANDELS and LEGA-C respectively. We finally split our sample into galaxies above and below our fitted mass-size relations, to investigate how size and D$_{n}$4000 correlate. For LEGA-C, we see a clear difference, with larger galaxies found to have smaller D$_{n}$4000 at fixed stellar mass. Due to the faintness and smaller numbers of the VANDELS sample, we cannot confirm whether a similar relation exists at $z\simeq1.1$. We consider whether differences in stellar age or metallicity are most likely to drive this size-D$_{n}$4000 relation, finding that any metallicity differences are unlikely to fully explain the observed offset, meaning smaller galaxies must be older than their larger counterparts. We find the observed evolution in size, mass and D$_{n}$4000 across the $\simeq2$ Gyr from $z\sim1.1$ to $z\sim0.7$ can be explained by a simple toy model in which VANDELS galaxies evolve passively, whilst experiencing a series of minor mergers. △ Less

Submitted 24 February, 2022; v1 submitted 25 January, 2022; originally announced January 2022.

Comments: 14 pages, 9 figures, submitted to MNRAS

The stellar metallicities of massive quiescent galaxies at 1.0 < z < 1.3 from KMOS+VANDELS

Authors: A. C. Carnall, R. J. McLure, J. S. Dunlop, M. Hamadouche, F. Cullen, D. J. McLeod, R. Begley, R. Amorin, M. Bolzonella, M. Castellano, A. Cimatti, F. Fontanot, A. Gargiulo, B. Garilli, F. Mannucci, L. Pentericci, M. Talia, G. Zamorani, A. Calabro, G. Cresci, N. P. Hathi

Abstract: We present a rest-frame UV-optical stacked spectrum representative of massive quiescent galaxies at $1.0<z<1.3$ with log$(M_*/\rm{M_\odot})>10.8$. The stack is constructed using VANDELS survey data, combined with new KMOS observations. We apply two independent full-spectral-fitting approaches, measuring a total metallicity, [Z/H]=$-0.13\pm0.08$ with Bagpipes, and [Z/H]=$0.04\pm0.14$ with Alf, a fa… ▽ More We present a rest-frame UV-optical stacked spectrum representative of massive quiescent galaxies at $1.0<z<1.3$ with log$(M_*/\rm{M_\odot})>10.8$. The stack is constructed using VANDELS survey data, combined with new KMOS observations. We apply two independent full-spectral-fitting approaches, measuring a total metallicity, [Z/H]=$-0.13\pm0.08$ with Bagpipes, and [Z/H]=$0.04\pm0.14$ with Alf, a fall of $\sim0.2-0.3$ dex compared with the local Universe. We also measure an iron abundance, [Fe/H] =$-0.18\pm0.08$, a fall of $\sim0.15$ dex compared with the the local Universe. We measure the alpha enhancement via the magnesium abundance, obtaining [Mg/Fe]=$0.23\pm$0.12, consistent with similar-mass galaxies in the local Universe, indicating no evolution in the average alpha enhancement of log$(M_*/\rm{M_\odot})=11$ quiescent galaxies over the last $\sim8$ Gyr. This suggests the very high alpha enhancements recently reported for several bright $z\sim1-2$ quiescent galaxies are due to their extreme masses, log$(M_*/\rm{M_\odot})\gtrsim11.5$, rather than being typical of the $z\gtrsim1$ population. The metallicity evolution we observe with redshift (falling [Z/H], [Fe/H], constant [Mg/Fe]) is consistent with recent studies. We recover a mean stellar age of $2.5^{+0.6}_{-0.4}$ Gyr, corresponding to a formation redshift, $z_\rm{form}=2.4^{+0.6}_{-0.3}$. Recent studies have obtained varying average formation redshifts for $z\gtrsim1$ massive quiescent galaxies, and, as these studies report consistent metallicities, we identify different star-formation-history models as the most likely cause. Larger spectroscopic samples from upcoming ground-based instruments will provide precise constraints on ages and metallicities at $z\gtrsim1$. Combining these with precise JWST $z>2$ quiescent-galaxy stellar-mass functions will provide an independent test of formation redshifts derived from spectral fitting. △ Less

Submitted 4 March, 2022; v1 submitted 30 August, 2021; originally announced August 2021.

Comments: 16 pages, 3 figures, accepted for publication in ApJ

The VANDELS survey: global properties of CIII]$λ$1908Å emitting star-forming galaxies at z$\sim$3

Authors: M. Llerena, R. Amorín, F. Cullen, L. Pentericci, A. Calabrò, R. McLure, A. Carnall, E. Pérez-Montero, F. Marchi, A. Bongiorno, M. Castellano, A. Fontana, D. J. McLeod, M. Talia, N. P. Hathi, P. Hibon, F. Mannucci, A. Saxena, D. Schaerer, G. Zamorani

Abstract: We study the mean properties of a large representative sample of 217 galaxies showing CIII] emission at $2<z<4$, selected from a parent sample of $\sim$750 main-sequence star-forming galaxies in the VANDELS survey. These CIII] emitters have a broad range of UV luminosities, thus allowing a detailed stacking analysis to characterize their stellar mass, star formation rate (SFR) and stellar metallic… ▽ More We study the mean properties of a large representative sample of 217 galaxies showing CIII] emission at $2<z<4$, selected from a parent sample of $\sim$750 main-sequence star-forming galaxies in the VANDELS survey. These CIII] emitters have a broad range of UV luminosities, thus allowing a detailed stacking analysis to characterize their stellar mass, star formation rate (SFR) and stellar metallicity, as a function of the UV emission line ratios, EWs, and the carbon-to-oxygen (C/O) abundance ratio. Reliable CIII] detections represent $\sim$30% of the parent sample. Extreme CIII] emitters (EW(CIII])$\gtrsim$8Å) are exceedingly rare ($\sim$3%) in VANDELS. The UV line ratios of the sample suggest no ionization source other than massive stars. Stacks with larger EW(CIII]) show larger EW(Ly$α$) and lower metallicity, but not all CIII] emitters are Ly$α$ emitters. The stellar metallicities of CIII] emitters are not significantly different from that of the parent sample, increasing from $\sim$10% to $\sim$40% solar for stellar masses $\log$(M$_{\star}$/M$_{\odot})\sim$9-10.5. The stellar mass-metallicity relation of the CIII] emitters is consistent with previous works showing strong evolution from $z=0$ to $z\sim3$. The C/O abundances of the sample range 35%-150% solar, with a noticeable increase with FUV luminosity and a smooth decrease with the CIII] EW. We discuss the CIII] emitters in the C/O-Fe/H and the C/O-O/H planes and find they follow stellar and nebular abundance trends consistent with those of Milky Way halo and thick disc stars and local HII galaxies, respectively. A qualitative agreement is also found with chemical evolution models, which suggests that CIII] emitters at $z\sim$3 are experiencing an active phase of chemical enrichment. △ Less

Submitted 10 November, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 659, A16 (2022)

The VANDELS Survey: New constraints on the high-mass X-ray binary populations in normal star-forming galaxies at 3 < z < 5.5

Authors: A. Saxena, R. S. Ellis, P. U. Forster, A. Calabro, L. Pentericci, A. C. Carnall, M. Castellano, F. Cullen, A. Fontana, M. Franco, J. P. U. Fynbo, A. Gargiulo, B. Garilli, N. P. Hathi, D. J. McLeod, R. Amorin, G. Zamorani

Abstract: We use VANDELS spectroscopic data overlapping with the $\simeq$7 Ms Chandra Deep Field South survey to extend studies of high-mass X-ray binary systems (XRBs) in 301 normal star-forming galaxies in the redshift range $3 < z < 5.5$. Our analysis evaluates correlations between X-ray luminosities ($L_X$), star formation rates (SFR) and stellar metallicities ($Z_\star$) to higher redshifts and over a… ▽ More We use VANDELS spectroscopic data overlapping with the $\simeq$7 Ms Chandra Deep Field South survey to extend studies of high-mass X-ray binary systems (XRBs) in 301 normal star-forming galaxies in the redshift range $3 < z < 5.5$. Our analysis evaluates correlations between X-ray luminosities ($L_X$), star formation rates (SFR) and stellar metallicities ($Z_\star$) to higher redshifts and over a wider range in galaxy properties than hitherto. Using a stacking analysis performed in bins of both redshift and SFR for sources with robust spectroscopic redshifts without AGN signatures, we find convincing evolutionary trends in the ratio $L_X$/SFR to the highest redshifts probed, with a stronger trend for galaxies with lower SFRs. Combining our data with published samples at lower redshift, the evolution of $L_X$/SFR to $z\simeq5$ proceeds as $(1 + z)^{1.03 \pm 0.02}$. Using stellar metallicities derived from photospheric absorption features in our spectroscopic data, we confirm indications at lower redshifts that $L_X$/SFR is stronger for metal-poor galaxies. We use semi-analytic models to show that metallicity dependence of $L_X$/SFR alone may not be sufficient to fully explain the observed redshift evolution of X-ray emission from high-mass XRBs, particularly for galaxies with SFR $<30$ $M_\odot$ yr$^{-1}$. We speculate that the discrepancy may arise due to reduced overall stellar ages in the early Universe leading to higher $L_X$/SFR for the same metallicity. We use our data to define the redshift-dependent contribution of XRBs to the integrated X-ray luminosity density and, in comparison with models, find that the contribution of high-mass XRBs to the cosmic X-ray background at $z>6$ may be $\gtrsim 0.25$ dex higher than previously estimated. △ Less

Submitted 5 July, 2021; v1 submitted 6 April, 2021; originally announced April 2021.

Comments: 16 pages, 10 figures, Published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 505, Issue 4, August 2021, Pages 4798-4812

The NIRVANDELS Survey: a robust detection of $α$-enhancement in star-forming galaxies at $z\simeq3.4$

Authors: F. Cullen, A. E. Shapley, R. J. McLure, J. S. Dunlop, R. L. Sanders, M. W. Topping, N. A. Reddy, R. Amorin, R. Begley, M. Bolzonella, A. Calabro, A. C. Carnall, M. Castellano, A. Cimatti, M/ Cirasuolo, G. Cresci, A. Fontana, F. Fontanot, B. Garilli, L. Guaita, M. Hamadouche, N. P. Hathi, F. Mannucci, D. J. McLeod, L. Pentericci , et al. (3 additional authors not shown)

Abstract: We present results from the NIRVANDELS survey investigating the gas-phase metallicity ($\mathrm{Z}_{\mathrm{gas}}$, tracing O/H) and stellar metallicity ($Z_{\star}$, tracing Fe/H) of 33 star-forming galaxies at redshifts $2.95 < z < 3.80$. Based on a combined analysis of deep optical and near-IR spectra, tracing the rest-frame far ultraviolet and rest-frame optical respectively, we present the fi… ▽ More We present results from the NIRVANDELS survey investigating the gas-phase metallicity ($\mathrm{Z}_{\mathrm{gas}}$, tracing O/H) and stellar metallicity ($Z_{\star}$, tracing Fe/H) of 33 star-forming galaxies at redshifts $2.95 < z < 3.80$. Based on a combined analysis of deep optical and near-IR spectra, tracing the rest-frame far ultraviolet and rest-frame optical respectively, we present the first simultaneous determination of the stellar and gas-phase mass-metallicity relationships (MZRs) at $z\simeq3.4$. In both cases, we find that metallicity increases with increasing stellar mass ($M_{\star}$), and that the power-law slope at $M_{\star} \lesssim 10^{10} \mathrm{M}_{\odot}$ of both MZRs scales as $Z \propto M_{\star}^{0.3}$. Comparing the stellar and gas-phase MZRs, we present direct evidence for super-solar O/Fe ratios (i.e., $α$-enhancement) at $z>3$, finding $\mathrm{(O/Fe)}\simeq (2.54 \pm 0.38) \times \mathrm{(O/Fe)}_{\odot}$, with no clear dependence on $M_{\star}$. △ Less

Submitted 7 May, 2021; v1 submitted 10 March, 2021; originally announced March 2021.

Comments: 18 pages, 11 figures, MNRAS accepted

The evolution of the galaxy stellar mass function over the last twelve billion years from a combination of ground-based and HST surveys

Authors: D. J. McLeod, R. J. McLure, J. S. Dunlop, F. Cullen, A. C. Carnall, K. Duncan

Abstract: We present a new determination of the galaxy stellar mass function (GSMF) over the redshift interval $0.25 \leq z \leq 3.75$, derived from a combination of ground-based and Hubble Space Telescope (HST) imaging surveys. Based on a near-IR selected galaxy sample selected over a raw survey area of 3 deg$^{2}$ and spanning $\geq 4$ dex in stellar mass, we fit the GSMF with both single and double Schec… ▽ More We present a new determination of the galaxy stellar mass function (GSMF) over the redshift interval $0.25 \leq z \leq 3.75$, derived from a combination of ground-based and Hubble Space Telescope (HST) imaging surveys. Based on a near-IR selected galaxy sample selected over a raw survey area of 3 deg$^{2}$ and spanning $\geq 4$ dex in stellar mass, we fit the GSMF with both single and double Schechter functions, carefully accounting for Eddington bias to derive both observed and intrinsic parameter values. We find that a double Schechter function is a better fit to the GSMF at all redshifts, although the single and double Schechter function fits are statistically indistinguishable by $z=3.25$. We find no evidence for significant evolution in $M^{\star}$, with the intrinsic value consistent with $\log_{10}(M^{\star} / M_{\odot})=10.55\pm{0.1}$ over the full redshift range. Overall, our determination of the GSMF is in good agreement with recent simulation results, although differences persist at the highest stellar masses. Splitting our sample according to location on the UVJ plane, we find that the star-forming GSMF can be adequately described by a single Schechter function over the full redshift range, and has not evolved significantly since $z\simeq 2.5$. In contrast, both the normalization and functional form of the passive GSMF evolves dramatically with redshift, switching from a single to a double Schechter function at $z \leq 1.5$. As a result, we find that while passive galaxies dominate the integrated stellar-mass density at $z \leq 0.75$, they only contribute $\lesssim 10$ per cent by $z\simeq 3$. Finally, we provide a simple parameterization that provides an accurate estimate of the GSMF, both observed and intrinsic, at any redshift within the range $0 \leq z \leq 4$. △ Less

Submitted 5 March, 2021; v1 submitted 7 September, 2020; originally announced September 2020.

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

X-ray properties of He II 1640 emitting galaxies in VANDELS

Authors: A. Saxena, L. Pentericci, D. Schaerer, R. Schneider, R. Amorin, A. Bongiorno, A. Calabrò, M. Castellano, A. Cimatti, F. Cullen, A. Fontana, J. P. U. Fynbo, N. Hathi, D. J. McLeod, M. Talia, G. Zamorani

Abstract: We explore X-ray emission from a sample of 18 He II 1640 emitting star-forming galaxies at z ~ 2.3-3.6 from the VANDELS survey in the Chandra Deep Field South, to set constraints on the role of X-ray sources in powering the He II emission. We find that 4 He II emitters have tentative detections with S/N ~ 2 and have X-ray luminosities, L_X = 1.5-4.9 x 10^41 erg/s. The stacked luminosity of all 18… ▽ More We explore X-ray emission from a sample of 18 He II 1640 emitting star-forming galaxies at z ~ 2.3-3.6 from the VANDELS survey in the Chandra Deep Field South, to set constraints on the role of X-ray sources in powering the He II emission. We find that 4 He II emitters have tentative detections with S/N ~ 2 and have X-ray luminosities, L_X = 1.5-4.9 x 10^41 erg/s. The stacked luminosity of all 18 He II emitters is 2.6 x 10^41 erg/s, and that of a subset of 13 narrow He II emitters (FHWM(He II) < 1000 km/s) is 3.1 x 10^41 erg/s. We also measure stacked L_X for non-He II emitters through bootstrapping of matched samples, and find L_X = 2.5 x 10^41 erg/s, which is not significantly different from L_X measured for He II emitters. The L_X per star-formation rate for He II emitters (log (L_X/SFR) ~ 40.0) and non-emitters (log (L_X/SFR) ~ 39.9) are also comparable and in line with the redshift evolution and metallicity dependence predicted by models. Due to the non-significant difference between the X-ray emission from galaxies with and without He II, we conclude that X-ray binaries or weak or obscured AGNs are unlikely to be the dominant producers of He II ionising photons in VANDELS star-forming galaxies at z ~ 3. Given the comparable physical properties of both He II emitters and non-emitters reported previously, alternative He II ionising mechanisms such as localised low-metallicity stellar populations, Pop-III stars, etc. may need to be explored. △ Less

Submitted 8 July, 2020; v1 submitted 30 March, 2020; originally announced March 2020.

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

The role of galaxy mass on AGN emission: a view from the VANDELS survey

Authors: M. Magliocchetti, L. Pentericci, M. Cirasuolo, G. Zamorani, R. Amorin, A. Bongiorno, A. Cimatti, A. Fontana, B. Garilli, A. Gargiulo, N. P. Hathi, D. J. McLeod, R. J. McLure, M. Brusa, A. Saxena, M. Talia

Abstract: We present a comparative analysis of the properties of AGN emitting at radio and X-ray wavelengths. The study is performed on 907 X-ray AGN and 100 radio AGN selected on the CDFS and UDS fields and makes use of new and ancillary data available to the VANDELS collaboration. Our results indicate that the mass of the host galaxy is a fundamental quantity which determines the level of AGN activity at… ▽ More We present a comparative analysis of the properties of AGN emitting at radio and X-ray wavelengths. The study is performed on 907 X-ray AGN and 100 radio AGN selected on the CDFS and UDS fields and makes use of new and ancillary data available to the VANDELS collaboration. Our results indicate that the mass of the host galaxy is a fundamental quantity which determines the level of AGN activity at the various wavelengths. Indeed large stellar masses are found to be connected with AGN radio emission, as virtually all radio-active AGN reside within galaxies of M*>10^{10} Msun. Large stellar masses also seem to favour AGN activity in the X-ray, even though X-ray AGN present a mass distribution which is more spread out and with a non-negligible tail at M*<10^{9} Msun. Stellar mass alone is also observed to play a fundamental role in simultaneous radio and X-ray emission: the percentage of AGN active at both wavelengths increases from around 1% of all X-ray AGN residing within hosts of M*<10^{11} Msun to about 13% in more massive galaxies. In the case of radio-selected AGN, such a percentage moves from about 15% to about 45% (but up to 80% in the deepest fields). Neither cosmic epoch, nor radio luminosity, X-ray luminosity, Eddington ratio or star-formation rate of the hosts are found to be connected to an enhanced probability for joint radio+X-ray emission of AGN origin. Furthermore, only a loose relation is observed between X-ray and radio luminosity in those AGN which are simultaneously active at both frequencies. △ Less

Submitted 7 February, 2020; originally announced February 2020.

Comments: 18 pages, 18 figures, to appear on MNRAS

Timing the earliest quenching events with a robust sample of massive quiescent galaxies at 2 < z < 5

Authors: A. C. Carnall, S. Walker, R. J. McLure, J. S. Dunlop, D. J. McLeod, F. Cullen, V. Wild, R. Amorin, M. Bolzonella, M. Castellano, A. Cimatti, O. Cucciati, A. Fontana, A. Gargiulo, B. Garilli, M. J. Jarvis, L. Pentericci, L. Pozzetti, G. Zamorani, A. Calabro, N. P. Hathi, A. M. Koekemoer

Abstract: We present a sample of 151 massive ($M_* > 10^{10}\mathrm{M_\odot}$) quiescent galaxies at $2 < z < 5$, based on a sophisticated Bayesian spectral energy distribution fitting analysis of the CANDELS UDS and GOODS-South fields. Our sample includes a robust sub-sample of 61 objects for which we confidently exclude low-redshift and star-forming solutions. We identify 10 robust objects at $z>3$, of wh… ▽ More We present a sample of 151 massive ($M_* > 10^{10}\mathrm{M_\odot}$) quiescent galaxies at $2 < z < 5$, based on a sophisticated Bayesian spectral energy distribution fitting analysis of the CANDELS UDS and GOODS-South fields. Our sample includes a robust sub-sample of 61 objects for which we confidently exclude low-redshift and star-forming solutions. We identify 10 robust objects at $z>3$, of which 2 are at $z>4$. We report formation redshifts, demonstrating that the oldest objects formed at $z > 6$, however individual ages from our photometric data have significant uncertainties, typically $\sim0.5$ Gyr. We demonstrate that the UVJ colours of the quiescent population evolve with redshift at $z>3$, becoming bluer and more similar to post-starburst galaxies at lower redshift. Based upon this we construct a model for the time-evolution of quiescent galaxy UVJ colours, concluding that the oldest objects are consistent with forming the bulk of their stellar mass at $z\sim6-7$ and quenching at $z\sim5$. We report spectroscopic redshifts for two of our objects at $z=3.440$ and $3.396$, which exhibit extremely weak Ly$α$ emission in ultra-deep VANDELS spectra. We calculate star-formation rates based on these line fluxes, finding that these galaxies are consistent with our quiescent selection criteria, provided their Ly$α$ escape fractions are $>3$ and $>10$ per cent respectively. We finally report that our highest-redshift robust object exhibits a continuum break at $λ\sim7000$A in a spectrum from VUDS, consistent with our photometric redshift of $z_\mathrm{phot}=4.72^{+0.06}_{-0.04}$. If confirmed as quiescent this object would be the highest-redshift known quiescent galaxy. To obtain stronger constraints on the times of the earliest quenching events, high-SNR spectroscopy must be extended to $z\gtrsim3$ quiescent objects. △ Less

Submitted 29 May, 2020; v1 submitted 31 January, 2020; originally announced January 2020.

Comments: 14 pages, 8 figures, MNRAS accepted

The VANDELS survey: A strong correlation between Ly$α$ equivalent width and stellar metallicity at $\mathbf{3\leq z \leq 5}$

Authors: F. Cullen, R. J. McLure, J. S. Dunlop, A. C. Carnall, D. J. McLeod, A. E. Shapley, R. Amorín, M. Bolzonella, M. Castellano, A. Cimatti, M. Cirasuolo, O. Cucciati, A. Fontana, F. Fontanot, B. Garilli, L. Guaita, M. J. Jarvis, L. Pentericci, L. Pozzetti, M. Talia, G. Zamorani, A. Calabrò, G. Cresci, J. P. U. Fynbo, N. P. Hathi , et al. (4 additional authors not shown)

Abstract: We present the results of a new study investigating the relationship between observed Ly$α$ equivalent width ($W_λ$(Ly$α$)) and the metallicity of the ionizing stellar population ($Z_{\star}$) for a sample of $768$ star-forming galaxies at $3 \leq z \leq 5$ drawn from the VANDELS survey. Dividing our sample into quartiles of rest-frame $W_λ$(Ly$α$) across the range $-58 \unicode{xC5} \lesssim$… ▽ More We present the results of a new study investigating the relationship between observed Ly$α$ equivalent width ($W_λ$(Ly$α$)) and the metallicity of the ionizing stellar population ($Z_{\star}$) for a sample of $768$ star-forming galaxies at $3 \leq z \leq 5$ drawn from the VANDELS survey. Dividing our sample into quartiles of rest-frame $W_λ$(Ly$α$) across the range $-58 \unicode{xC5} \lesssim$ $W_λ$(Ly$α$) $\lesssim 110 \unicode{xC5}$ we determine $Z_{\star}$ from full spectral fitting of composite far-ultraviolet (FUV) spectra and find a clear anti-correlation between $W_λ$(Ly$α$) and $Z_{\star}$. Our results indicate that $Z_{\star}$ decreases by a factor $\gtrsim 3$ between the lowest $W_λ$(Ly$α$) quartile ($\langle$$W_λ$(Ly$α$)$\rangle=-18\unicode{xC5}$) and the highest $W_λ$(Ly$α$) quartile ($\langle$$W_λ$(Ly$α$)$\rangle=24\unicode{xC5}$). Similarly, galaxies typically defined as Lyman Alpha Emitters (LAEs; $W_λ$(Ly$α$) $>20\unicode{xC5}$) are, on average, metal poor with respect to the non-LAE galaxy population ($W_λ$(Ly$α$) $\leq20\unicode{xC5}$) with $Z_{\star}$$_{\rm{non-LAE}}\gtrsim 2 \times$ $Z_{\star}$$_{\rm{LAE}}$. Finally, based on the best-fitting stellar models, we estimate that the increasing strength of the stellar ionizing spectrum towards lower $Z_{\star}$ is responsible for $\simeq 15-25\%$ of the observed variation in $W_λ$(Ly$α$) across our sample, with the remaining contribution ($\simeq 75-85\%$) being due to a decrease in the HI/dust covering fractions in low $Z_{\star}$ galaxies. △ Less

Submitted 16 April, 2020; v1 submitted 29 January, 2020; originally announced January 2020.

Comments: 10 pages, 6 figures, MNRAS accepted

A lack of evolution in the very bright-end of the galaxy luminosity function from z = 8-10

Authors: R. A. A. Bowler, M. J. Jarvis, J. S. Dunlop, R. J. McLure, D. J. McLeod, N. J. Adams, B. Milvang-Jensen, H. J. McCracken

Abstract: We utilize deep near-infrared survey data from the UltraVISTA fourth data release (DR4) and the VIDEO survey, in combination with overlapping optical and Spitzer data, to search for bright star-forming galaxies at $z \gtrsim 7.5$. Using a full photometric redshift fitting analysis applied to the $\sim 6\,{\rm deg}^2$ of imaging searched, we find 27 Lyman-break galaxies (LBGs), including 20 new sou… ▽ More We utilize deep near-infrared survey data from the UltraVISTA fourth data release (DR4) and the VIDEO survey, in combination with overlapping optical and Spitzer data, to search for bright star-forming galaxies at $z \gtrsim 7.5$. Using a full photometric redshift fitting analysis applied to the $\sim 6\,{\rm deg}^2$ of imaging searched, we find 27 Lyman-break galaxies (LBGs), including 20 new sources, with best-fitting photometric redshifts in the range $7.4 < z < 9.1$. From this sample we derive the rest-frame UV luminosity function (LF) at $z = 8$ and $z = 9$ out to extremely bright UV magnitudes ($M_{\rm UV} \simeq -23$) for the first time. We find an excess in the number density of bright galaxies in comparison to the typically assumed Schechter functional form derived from fainter samples. Combined with previous studies at lower redshift, our results show that there is little evolution in the number density of very bright ($M_{\rm UV} \sim -23$) LBGs between $z \simeq 5$ and $z\simeq 9$. The tentative detection of an LBG with best-fit photometric redshift of $z = 10.9 \pm 1.0$ in our data is consistent with the derived evolution. We show that a double power-law fit with a brightening characteristic magnitude ($ΔM^*/Δz \simeq -0.5$) and a steadily steepening bright-end slope ($Δβ/Δz \simeq -0.5$) provides a good description of the $z > 5$ data over a wide range in absolute UV magnitude ($-23 < M_{\rm UV} < -17$). We postulate that the observed evolution can be explained by a lack of mass quenching at very high redshifts in combination with increasing dust obscuration within the first $\sim 1 \,{\rm Gyr}$ of galaxy evolution. △ Less

Submitted 30 January, 2020; v1 submitted 28 November, 2019; originally announced November 2019.

Comments: Updated to match MNRAS accepted version (accepted 27/01/2020) after minor revisions. 18 pages, 10 figures and 7 tables in the main text, plus additional figures in the appendix

The properties of He II 1640 emitters at z ~ 2.5-5 from the VANDELS survey

Authors: A. Saxena, L. Pentericci, M. Mirabelli, D. Schaerer, R. Schneider, F. Cullen, R. Amorin, M. Bolzonella, A. Bongiorno, A. C. Carnall, M. Castellano, O. Cucciati, A. Fontana, J. P. U. Fynbo, B. Garilli, A. Gargiulo, L. Guaita, N. P. Hathi, T. A. Hutchison, A. M. Koekemoer, F. Marchi, D. J. McLeod, R. J. McLure, C. Papovich, L. Pozzetti , et al. (2 additional authors not shown)

Abstract: Strong He II emission is produced by low-metallicity stellar populations. Here, we aim to identify and study a sample of He II $λ1640$-emitting galaxies at redshifts of $z \sim 2.5-5$ in the deep VANDELS spectroscopic survey.. We identified a total of 33 Bright He II emitters (S/N > 2.5) and 17 Faint emitters (S/N < 2.5) in the VANDELS survey and used the available deep multi-wavelength data to st… ▽ More Strong He II emission is produced by low-metallicity stellar populations. Here, we aim to identify and study a sample of He II $λ1640$-emitting galaxies at redshifts of $z \sim 2.5-5$ in the deep VANDELS spectroscopic survey.. We identified a total of 33 Bright He II emitters (S/N > 2.5) and 17 Faint emitters (S/N < 2.5) in the VANDELS survey and used the available deep multi-wavelength data to study their physical properties. After identifying seven potential AGNs in our sample and discarding them from further analysis, we divided the sample of \emph{Bright} emitters into 20 \emph{Narrow} (FWHM < 1000 km s$^{-1}$) and 6 \emph{Broad} (FWHM > 1000 km s$^{-1}$) He II emitters. We created stacks of Faint, Narrow, and Broad emitters and measured other rest-frame UV lines such as O III] and C III] in both individual galaxies and stacks. We then compared the UV line ratios with the output of stellar population-synthesis models to study the ionising properties of He II emitters. We do not see a significant difference between the stellar masses, star-formation rates, and rest-frame UV magnitudes of galaxies with He II and no He II emission. The stellar population models reproduce the observed UV line ratios from metals in a consistent manner, however they under-predict the total number of \heii ionising photons, confirming earlier studies and suggesting that additional mechanisms capable of producing He II are needed, such as X-ray binaries or stripped stars. The models favour subsolar metallicities ($\sim0.1Z_\odot$) and young stellar ages ($10^6 - 10^7$ years) for the He II emitters. However, the metallicity measured for He II emitters is comparable to that of non-He II emitters at similar redshifts. We argue that galaxies with He II emission may have undergone a recent star-formation event, or may be powered by additional sources of He II ionisation. △ Less

Submitted 23 March, 2020; v1 submitted 22 November, 2019; originally announced November 2019.

Comments: 21 pages, 11 figures (including appendix), Accepted for publication in A&A

Journal ref: A&A 636, A47 (2020)

High-velocity outflows in massive post-starburst galaxies at z > 1

Authors: David T. Maltby, Omar Almaini, Ross J. McLure, Vivienne Wild, James Dunlop, Kate Rowlands, William G. Hartley, Nina A. Hatch, Miguel Socolovsky, Aaron Wilkinson, Ricardo Amorin, Emma J. Bradshaw, Adam C. Carnall, Marco Castellano, Andrea Cimatti, Giovanni Cresci, Fergus Cullen, Stephane De Barros, Fabio Fontanot, Bianca Garilli, Anton M. Koekemoer, Derek J. McLeod, Laura Pentericci, Margherita Talia

Abstract: We investigate the prevalence of galactic-scale outflows in post-starburst (PSB) galaxies at high redshift ($1 < z < 1.4$), using the deep optical spectra available in the UKIDSS Ultra Deep Survey (UDS). We use a sample of $\sim40$ spectroscopically confirmed PSBs, recently identified in the UDS field, and perform a stacking analysis in order to analyse the structure of strong interstellar absorpt… ▽ More We investigate the prevalence of galactic-scale outflows in post-starburst (PSB) galaxies at high redshift ($1 < z < 1.4$), using the deep optical spectra available in the UKIDSS Ultra Deep Survey (UDS). We use a sample of $\sim40$ spectroscopically confirmed PSBs, recently identified in the UDS field, and perform a stacking analysis in order to analyse the structure of strong interstellar absorption features such as Mg ii ($\lambda2800$ Ang.). We find that for massive ($M_* > 10^{10}\rm\,M_{\odot}$) PSBs at $z > 1$, there is clear evidence for a strong blue-shifted component to the Mg ii absorption feature, indicative of high-velocity outflows ($v_{\rm out}\sim1150\pm160\rm\,km\,s^{-1}$) in the interstellar medium. We conclude that such outflows are typical in massive PSBs at this epoch, and potentially represent the residual signature of a feedback process that quenched these galaxies. Using full spectral fitting, we also obtain a typical stellar velocity dispersion $σ_*$ for these PSBs of $\sim200\rm\,km\,s^{-1}$, which confirms they are intrinsically massive in nature (dynamical mass $M_{\rm d}\sim10^{11}\rm\,M_{\odot}$). Given that these high-$z$ PSBs are also exceptionally compact ($r_{\rm e}\sim1$--$2\rm\,kpc$) and spheroidal (Sersic index $n\sim3$), we propose that the outflowing winds may have been launched during a recent compaction event (e.g. major merger or disc collapse) that triggered either a centralised starburst or active galactic nuclei (AGN) activity. Finally, we find no evidence for AGN signatures in the optical spectra of these PSBs, suggesting they were either quenched by stellar feedback from the starburst itself, or that if AGN feedback is responsible, the AGN episode that triggered quenching does not linger into the post-starburst phase. △ Less

Submitted 7 August, 2019; originally announced August 2019.

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

The VANDELS survey: the star-formation histories of massive quiescent galaxies at 1.0 < z < 1.3

Authors: A. C. Carnall, R. J. McLure, J. S. Dunlop, F. Cullen, D. J. McLeod, V. Wild, B. D. Johnson, S. Appleby, R. Davé, R. Amorin, M. Bolzonella, M. Castellano, A. Cimatti, O. Cucciati, A. Gargiulo, B. Garilli, F. Marchi, L. Pentericci, L. Pozzetti, C. Schreiber, M. Talia, G. Zamorani

Abstract: We present a Bayesian full-spectral-fitting analysis of 75 massive ($M_* > 10^{10.3} M_\odot$) UVJ-selected galaxies at redshifts of $1.0 < z < 1.3$, combining extremely deep rest-frame ultraviolet spectroscopy from VANDELS with multi-wavelength photometry. By the use of a sophisticated physical plus systematic uncertainties model, constructed within the Bagpipes code, we place strong constraints… ▽ More We present a Bayesian full-spectral-fitting analysis of 75 massive ($M_* > 10^{10.3} M_\odot$) UVJ-selected galaxies at redshifts of $1.0 < z < 1.3$, combining extremely deep rest-frame ultraviolet spectroscopy from VANDELS with multi-wavelength photometry. By the use of a sophisticated physical plus systematic uncertainties model, constructed within the Bagpipes code, we place strong constraints on the star-formation histories (SFHs) of individual objects. We firstly constrain the stellar mass vs stellar age relationship, finding a steep trend towards earlier average formation with increasing stellar mass of $1.48^{+0.34}_{-0.39}$ Gyr per decade in mass, although this shows signs of flattening at $M_* > 10^{11} M_\odot$. We show that this is consistent with other spectroscopic studies from $0 < z < 2$. This relationship places strong constraints on the AGN-feedback models used in cosmological simulations. We demonstrate that, although the relationships predicted by Simba and IllustrisTNG agree well with observations at $z=0.1$, they are too shallow at $z=1$, predicting an evolution of $<0.5$ Gyr per decade in mass. Secondly, we consider the connections between green-valley, post-starburst and quiescent galaxies, using our inferred SFH shapes and the distributions of galaxy physical properties on the UVJ diagram. The majority of our lowest-mass galaxies ($M_* \sim 10^{10.5} M_\odot$) are consistent with formation in recent ($z<2$), intense starburst events, with timescales of $\lesssim500$ Myr. A second class of objects experience extended star-formation epochs before rapidly quenching, passing through both green-valley and post-starburst phases. The most massive galaxies in our sample are extreme systems: already old by $z=1$, they formed at $z\sim5$ and quenched by $z=3$. However, we find evidence for their continued evolution through both AGN and rejuvenated star-formation activity. △ Less

Submitted 24 October, 2019; v1 submitted 26 March, 2019; originally announced March 2019.

Comments: 25 pages, 13 figures, MNRAS accepted

The VANDELS survey: the stellar metallicities of star-forming galaxies at 2.5 < z < 5.0

Authors: F. Cullen, R. J. McLure, J. S. Dunlop, S. Khochfar, R. Davé, R. Amorin, M. Bolzonella, A. C. Carnall, M. Castellano, A. Cimatti, M. Cirasuolo, G. Cresci, J. P. U. Fynbo, F. Fontanot, A. Gargiulo, B. Garilli, L. Guaita, N. Hathi, P. Hibon, F. Mannucci, F. Marchi, D. J. McLeod, L. Pentericci, L. Pozzetti, A. E. Shapley , et al. (2 additional authors not shown)

Abstract: We present the results of a study utilising ultra-deep, rest-frame UV, spectroscopy to quantify the relationship between stellar mass and stellar metallicity for 681 star-forming galaxies at $2.5<z<5.0$ ($\langle z \rangle = 3.5 \pm 0.6$) drawn from the VANDELS survey. Via a comparison with high-resolution stellar population models, we determine stellar metallicities for a set of composite spectra… ▽ More We present the results of a study utilising ultra-deep, rest-frame UV, spectroscopy to quantify the relationship between stellar mass and stellar metallicity for 681 star-forming galaxies at $2.5<z<5.0$ ($\langle z \rangle = 3.5 \pm 0.6$) drawn from the VANDELS survey. Via a comparison with high-resolution stellar population models, we determine stellar metallicities for a set of composite spectra formed from subsamples selected by mass and redshift. Across the stellar mass range $8.5 < \mathrm{log}(\langle M_{\ast} \rangle/\rm{M}_{\odot}) < 10.2$ we find a strong correlation between stellar metallicity and stellar mass, with stellar metallicity monotonically increasing from $Z_{\ast}/\mathrm{Z}_{\odot} < 0.09$ at $\langle M_{\ast} \rangle = 3.2 \times 10^{8} \rm{M}_{\odot}$ to $Z_{\ast}/Z_{\odot} = 0.27$ at $\langle M_{\ast} \rangle = 1.7 \times 10^{10} \rm{M}_{\odot}$. In contrast, at a given stellar mass, we find no evidence for significant metallicity evolution across the redshift range of our sample. However, comparing our results to the $z=0$ stellar mass-metallicity relation, we find that the $\langle z \rangle = 3.5$ relation is consistent with being shifted to lower metallicities by $\simeq 0.6$ dex. Contrasting our derived stellar metallicities with estimates of gas-phase metallicities at similar redshifts, we find evidence for enhanced $\rm{O}/\rm{Fe}$ ratios of the order (O/Fe) $\gtrsim 1.8$ $\times$ (O/Fe)$_{\odot}$. Finally, by comparing our results to simulation predictions, we find that the $\langle z \rangle = 3.5$ stellar mass-metallicity relation is consistent with current predictions for how outflow strength scales with galaxy mass. This conclusion is supported by an analysis of analytic models, and suggests that the mass loading parameter ($η=\dot{M}_{\mathrm{outflow}}/M_{\ast}$) scales as $η\propto M_{\ast}^β$ with $β\simeq -0.4$. △ Less

Submitted 22 May, 2019; v1 submitted 26 March, 2019; originally announced March 2019.

Comments: 21 pages (+ appendix), 13 figures, MNRAS accepted

The relationship between dust and [CI] at z=1 and beyond

Authors: N. Bourne, J. S. Dunlop, J. M. Simpson, K. E. Rowlands, J. E. Geach, D. J. McLeod

Abstract: Measuring molecular gas mass is vital for understanding the evolution of galaxies at high redshifts (z$\geq$1). Most measurements rely on CO as a tracer, but dependences on metallicity, dynamics and surface density lead to systematic uncertainties in high-z galaxies, where these physical properties are difficult to observe, and where the physical environments can differ systematically from those a… ▽ More Measuring molecular gas mass is vital for understanding the evolution of galaxies at high redshifts (z$\geq$1). Most measurements rely on CO as a tracer, but dependences on metallicity, dynamics and surface density lead to systematic uncertainties in high-z galaxies, where these physical properties are difficult to observe, and where the physical environments can differ systematically from those at z=0. Dust continuum emission provides a potential alternative assuming a known dust/gas ratio, but this must be calibrated on a direct gas tracer at z$\geq$1. In this paper we consider the [CI] 492-GHz emission line, which has been shown to trace molecular gas closely throughout Galactic clouds and has the advantages of being optically thin in typical conditions (unlike CO), and being observable at accessible frequencies at high redshifts (in contrast to the low-excitation lines of CO). We use the Atacama Large Millimetre/submillimetre Array (ALMA) to measure [CI], CO(4-3) and dust emission in a representative sample of star-forming galaxies at z=1, and combine these data with multi-wavelength spectral energy distributions to study relationships between dust and gas components of galaxies. We uncover a strong [CI]-dust correlation, suggesting that both trace similar phases of the gas. By incorporating other samples from the literature, we show that this correlation persists over a wide range of luminosities and redshifts up to z$\sim$4. Finally we explore the implications of our results as an independent test of literature calibrations for dust as a tracer of gas mass, and for predicting the CI abundance. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: MNRAS accepted. 24 pages, 12 figures

The VANDELS ESO public spectroscopic survey

Authors: R. J. McLure, L. Pentericci, A. Cimatti, J. S. Dunlop, D. Elbaz, A. Fontana, K. Nandra, R. Amorin, M. Bolzonella, A. Bongiorno, A. C. Carnall, M. Castellano, M. Cirasuolo, O. Cucciati, F. Cullen, S. De Barros, S. L. Finkelstein, F. Fontanot, P. Franzetti, M. Fumana, A. Gargiulo, B. Garilli, L. Guaita, W. G. Hartley, A. Iovino , et al. (70 additional authors not shown)

Abstract: VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric re… ▽ More VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric redshift pre-selection, 85% of the galaxies targeted by VANDELS were selected to be at z>=3. Exploiting the red sensitivity of the refurbished VIMOS spectrograph, the fundamental aim of the survey is to provide the high signal-to-noise ratio spectra necessary to measure key physical properties such as stellar population ages, masses, metallicities and outflow velocities from detailed absorption-line studies. Using integration times calculated to produce an approximately constant signal-to-noise ratio (20 < t_int < 80 hours), the VANDELS survey targeted: a) bright star-forming galaxies at 2.4 < z < 5.5, b) massive quiescent galaxies at 1.0 < z < 2.5, c) fainter star-forming galaxies at 3.0 < z < 7.0 and d) X-ray/Spitzer-selected active galactic nuclei and Herschel-detected galaxies. By targeting two extragalactic survey fields with superb multi-wavelength imaging data, VANDELS will produce a unique legacy data set for exploring the physics underpinning high-redshift galaxy evolution. In this paper we provide an overview of the VANDELS survey designed to support the science exploitation of the first ESO public data release, focusing on the scientific motivation, survey design and target selection. △ Less

Submitted 14 May, 2018; v1 submitted 20 March, 2018; originally announced March 2018.

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

Obscured star-formation in bright z ~ 7 Lyman-break galaxies

Authors: R. A. A. Bowler, N. Bourne, J. S. Dunlop, R. J. McLure, D. J. McLeod

Abstract: We present Atacama Large Millimeter/Submillimeter Array observations of the rest-frame far-infrared (FIR) dust continuum emission of six bright Lyman-break galaxies (LBGs) at $z \simeq 7$. One LBG is detected ($5.2σ$ at peak emission), while the others remain individually undetected at the $3σ$ level. The average FIR luminosity of the sample is found to be… ▽ More We present Atacama Large Millimeter/Submillimeter Array observations of the rest-frame far-infrared (FIR) dust continuum emission of six bright Lyman-break galaxies (LBGs) at $z \simeq 7$. One LBG is detected ($5.2σ$ at peak emission), while the others remain individually undetected at the $3σ$ level. The average FIR luminosity of the sample is found to be $L_{\rm FIR} \simeq 2 \times 10^{11}\,{\rm L}_{\odot}$, corresponding to an obscured star-formation rate (SFR) that is comparable to that inferred from the unobscured UV emission. In comparison to the infrared excess (IRX$\,=L_{\rm FIR}/L_{\rm UV}$)-$β$ relation, our results are consistent with a Calzetti-like attenuation law (assuming a dust temperature of T = 40-50 K). We find a physical offset of 3 kpc between the dust continuum emission and the rest-frame UV light probed by Hubble Space Telescope imaging for galaxy ID65666 at $z = 7.17^{+0.09}_{-0.06}$. The offset is suggestive of an inhomogeneous dust distribution, where 75% of the total star formation activity (SFR$ \,\simeq 70\,{\rm M}_{\odot}/{\rm yr}$) of the galaxy is completely obscured. Our results provide direct evidence that dust obscuration plays a key role in shaping the bright-end of the observed rest-frame UV luminosity function at $z \simeq 7$, in agreement with cosmological galaxy formation simulations. The existence of a heavily-obscured component of galaxy ID65666 indicates that dusty star-forming regions, or even entire galaxies, that are "UV-dark" are significant even in the $z \simeq 7$ galaxy population. △ Less

Submitted 29 August, 2018; v1 submitted 15 February, 2018; originally announced February 2018.

Comments: 15 pages, 8 figures, 3 tables. Updated to match MNRAS accepted version

Characterising the evolving $K$-band luminosity function using the UltraVISTA, CANDELS and HUDF surveys

Authors: Alice Mortlock, Ross J. McLure, Rebecca A. A. Bowler, Derek J. McLeod, Esther Marmol-Queralto, Shaghayegh Parsa, James S. Dunlop, Victoria A. Bruce

Abstract: We present the results of a new study of the K-band galaxy luminosity function (KLF) at redshifts z<3.75, based on a nested combination of the UltraVISTA, CANDELS and HUDF surveys. The large dynamic range in luminosity spanned by this new dataset (3-4 dex over the full redshift range) is sufficient to clearly demonstrate for the first time that the faint-end slope of the KLF at z>0.25 is relativel… ▽ More We present the results of a new study of the K-band galaxy luminosity function (KLF) at redshifts z<3.75, based on a nested combination of the UltraVISTA, CANDELS and HUDF surveys. The large dynamic range in luminosity spanned by this new dataset (3-4 dex over the full redshift range) is sufficient to clearly demonstrate for the first time that the faint-end slope of the KLF at z>0.25 is relatively steep (-1.3<alpha<-1.5 for a single Schechter function), in good agreement with recent theoretical and phenomenological models. Moreover, based on our new dataset we find that a double Schechter function provides a significantly improved description of the KLF at z<2. At redshifts z>0.25 the evolution of the KLF is remarkably smooth, with little or no evolution evident at faint (M_K>-20.5) or bright magnitudes (M_K<-24.5). Instead, the KLF is seen to evolve rapidly at intermediate magnitudes, with the number density of galaxies at M_K~-23 dropping by a factor of ~5 over the redshift interval 0.25<z<3.75. Motivated by this, we explore a simple description of the evolving KLF based on a double Schechter function with fixed faint-end slopes (alpha_1=-0.5, alpha_2=-1.5) and a shared characteristic magnitude (M_K*). According to this parameterisation, the normalisation of the component which dominates the faint-end of the KLF remains approximately constant, with phi*_2 decreasing by only a factor of ~2 between z~0 and z~3.25. In contrast, the component which dominates the bright end of the KLF at low redshifts evolves dramatically, becoming essentially negligible by z~3. Finally, we note that within this parameterisation, the observed evolution of M_K* between z~0 and z~3.25 is entirely consistent with M_K* corresponding to a constant stellar mass of M*~5x10^10 Msun at all redshifts. △ Less

Submitted 20 October, 2016; originally announced October 2016.

Comments: 18 pages, 10 figures, accepted for publication in MNRAS

No evidence for Population III stars or a Direct Collapse Black Hole in the z = 6.6 Lyman-$α$ emitter 'CR7'

Authors: R. A. A. Bowler, R. J. McLure, J. S. Dunlop, D. J. McLeod, E. R. Stanway, J. J. Eldridge, M. J. Jarvis

Abstract: The z = 6.6 Lyman-$α$ emitter 'CR7' has been claimed to have a Population III-like stellar population, or alternatively, be a candidate Direct Collapse Black Hole (DCBH). In this paper we investigate the evidence for these exotic scenarios using recently available, deeper, optical, near-infrared and mid-infrared imaging. We find strong Spitzer/IRAC detections for the main component of CR7 at 3.6 a… ▽ More The z = 6.6 Lyman-$α$ emitter 'CR7' has been claimed to have a Population III-like stellar population, or alternatively, be a candidate Direct Collapse Black Hole (DCBH). In this paper we investigate the evidence for these exotic scenarios using recently available, deeper, optical, near-infrared and mid-infrared imaging. We find strong Spitzer/IRAC detections for the main component of CR7 at 3.6 and 4.5 microns, and show that it has a blue colour ([3.6] - [4.5] $= -1.2\pm 0.3$). This colour cannot be reproduced by current Pop. III or pristine DCBH models. Instead, the results suggest that the [3.6] band is contaminated by the [OIII]4959,5007 emission line with an implied rest-frame equivalent width of EW_0 (H$β$ + [OIII]) $\gtrsim 2000$Å. Furthermore, we find that new near-infrared data from the UltraVISTA survey supports a weaker He II 1640 emission line than previously measured, with EW_0 $= 40 \pm 30$Å. For the fainter components of CR7 visible in Hubble Space Telescope imaging, we find no evidence that they are particularly red as previously claimed, and show that the derived masses and ages are considerably uncertain. In light of the likely detection of strong [OIII] emission in CR7 we discuss other more standard interpretations of the system that are consistent with the data. We find that a low-mass, narrow-line AGN can reproduce the observed features of CR7, including the lack of radio and X-ray detections. Alternatively, a young, low-metallicity (~1/200 solar) starburst, modelled including binary stellar pathways, can reproduce the inferred strength of the He II line and simultaneously the strength of the observed [OIII] emission, but only if the gas shows super-solar $α$-element abundances (O/Fe ~ 5 O/Fe solar). △ Less

Submitted 3 April, 2017; v1 submitted 2 September, 2016; originally announced September 2016.

Comments: 13 pages, 8 figures and 1 table. Updated to MNRAS accepted version

Unveiling the nature of bright z ~ 7 galaxies with the Hubble Space Telescope

Authors: R. A. A. Bowler, J. S. Dunlop, R. J. McLure, D. J. McLeod

Abstract: We present new Hubble Space Telescope/Wide Field Camera 3 imaging of 25 extremely luminous (-23.2 < M_ UV < -21.2) Lyman-break galaxies (LBGs) at z ~ 7. The sample was initially selected from 1.65 deg^2 of ground-based imaging in the UltraVISTA/COSMOS and UDS/SXDS fields, and includes the extreme Lyman-alpha emitters, `Himiko' and `CR7'. A deconfusion analysis of the deep Spitzer photometry availa… ▽ More We present new Hubble Space Telescope/Wide Field Camera 3 imaging of 25 extremely luminous (-23.2 < M_ UV < -21.2) Lyman-break galaxies (LBGs) at z ~ 7. The sample was initially selected from 1.65 deg^2 of ground-based imaging in the UltraVISTA/COSMOS and UDS/SXDS fields, and includes the extreme Lyman-alpha emitters, `Himiko' and `CR7'. A deconfusion analysis of the deep Spitzer photometry available suggests that these galaxies exhibit strong rest-frame optical nebular emission lines (EW_0(H_beta + [OIII]) > 600A). We find that irregular, multiple-component morphologies suggestive of clumpy or merging systems are common (f_multi > 0.4) in bright z ~ 7 galaxies, and ubiquitous at the very bright end (M_UV < -22.5). The galaxies have half-light radii in the range r_1/2 ~ 0.5-3 kpc. The size measurements provide the first determination of the size-luminosity relation at z ~ 7 that extends to M_UV ~ -23. We find the relation to be steep with r_1/2 ~ L^1/2. Excluding clumpy, multi-component galaxies however, we find a shallower relation that implies an increased star-formation rate surface density in bright LBGs. Using the new, independent, HST/WFC3 data we confirm that the rest-frame UV luminosity function at z ~ 7 favours a power-law decline at the bright-end, compared to an exponential Schechter function drop-off. Finally, these results have important implications for the Euclid mission, which we predict will detect > 1000 similarly bright galaxies at z ~ 7. Our new HST imaging suggests that the vast majority of these galaxies will be spatially resolved by Euclid, mitigating concerns over dwarf star contamination. △ Less

Submitted 19 December, 2016; v1 submitted 17 May, 2016; originally announced May 2016.

Comments: 26 pages, 11 figures and 5 tables. Updated to match MNRAS accepted version

The z = 9-10 galaxy population in the Hubble Frontier Fields and CLASH surveys: The z=9 LF and further evidence for a smooth decline in UV luminosity at z >= 8

Authors: D. J. McLeod, R. J. McLure, J. S. Dunlop

Abstract: We present the results of a search for z=9-10 galaxies within the first 8 pointings of the Hubble Frontier Fields (HFF) (4 clusters plus 4 parallel fields) and 20 cluster fields from the CLASH survey. Combined with our previous analysis of the Hubble Ultra-Deep field (HUDF), we have now completed a search for z=9-10 galaxies over ~130 sq. arcmin, across 29 HST WFC3/IR pointings. As in our recent s… ▽ More We present the results of a search for z=9-10 galaxies within the first 8 pointings of the Hubble Frontier Fields (HFF) (4 clusters plus 4 parallel fields) and 20 cluster fields from the CLASH survey. Combined with our previous analysis of the Hubble Ultra-Deep field (HUDF), we have now completed a search for z=9-10 galaxies over ~130 sq. arcmin, across 29 HST WFC3/IR pointings. As in our recent study of the first two HFF fields, we confine our primary search for high-redshift candidates in the HFF imaging to the uniformly deep (i.e. sigma_160>30 AB mag in 0.5-arcsec diameter apertures), relatively low magnification regions. In the CLASH fields our search was confined to uniformly deep regions where sigma_160>28.8 AB mag. Our SED fitting analysis unveils a sample of 33 galaxy candidates at z_phot>=8.4, five of which have primary photometric redshift solutions in the range 9.6<z_phot<11.2. By calculating a de-lensed effective volume for each candidate, the improved statistics and reduced cosmic variance provided by our new sample allows a more accurate determination of the UV-selected galaxy luminosity function (LF) at z~9. Our new results strengthen our previous conclusion that the LF appears to evolve smoothly from z=8 to z=9, an evolution which can be equally well modelled by a factor of ~2 drop in density, or a dimming of ~0.5 mag in M*. Moreover, based on our new sample, we are able to place initial constraints on the z=10 LF, finding that the number density at M_1500 ~ -19.7 is log(phi) = -4.1 (+0.2,-0.3), a factor of ~2 lower than at z=9. Finally, we use our new results to re-visit the issue of the decline in UV luminosity density at z>=8. We conclude that the data continue to support a smooth decline in rho_UV over the redshift interval 6<z<10, in agreement with simple models of early galaxy evolution driven by the growth in the underlying dark matter halo mass function. △ Less

Submitted 16 February, 2016; originally announced February 2016.

Comments: Submitted to MNRAS, 13 pages, 11 figures

The evolution of the equivalent width of the Ha emission line and specific star-formation rate in star-forming galaxies at 1<z<5

Authors: E. Marmol-Queralto, R. J. McLure, F. Cullen, J. S. Dunlop, A. Fontana, D. J. McLeod

Abstract: We present the results of a study which uses spectral energy distribution (SED) fitting to investigate the evolution of the equivalent width (EW) of the Halpha emission line in star-forming galaxies over the redshift interval 1<z<5. After first demonstrating the ability of our SED-fitting technique to recover EW(Ha) using a sample of galaxies at z~1.3 with EW(Ha) measurements from 3D-HST grism spe… ▽ More We present the results of a study which uses spectral energy distribution (SED) fitting to investigate the evolution of the equivalent width (EW) of the Halpha emission line in star-forming galaxies over the redshift interval 1<z<5. After first demonstrating the ability of our SED-fitting technique to recover EW(Ha) using a sample of galaxies at z~1.3 with EW(Ha) measurements from 3D-HST grism spectroscopy, we proceed to apply our technique to samples of spectroscopically confirmed and photometric-redshift selected star-forming galaxies at z>=1 in the CANDELS UDS and GOODS-S fields. Confining our analysis to a constant stellar mass range (9.5<log(M/Msun)<10.5), we find that the median EW(Ha) evolves only modestly with redshift, reaching a rest-frame value of EW(Ha)=301+/-30 Angs by redshift z~4.5. Furthermore, using estimates of star-formation rate (SFR) based on both UV luminosity and Ha line flux, we use our galaxy samples to compare the evolution of EW(Ha) and specific star-formation rate (sSFR). Our results indicate that over the redshift range 1<z<5, the evolution displayed by EW(Ha) and sSFR is consistent, and can be adequately parameterized as: propto (1+z)^(1.0+/-0.2). As a consequence, over this redshift range we find that the sSFR and rest-frame EW(Ha) of star-forming galaxies with stellar masses M~10^(10) Msun are related by: EW(Ha)/Ang=(63+/-7)sSFR/Gyr^(-1). Given the current uncertainties in measuring the SFRs of high-redshift galaxies, we conclude that EW(Ha) provides a useful independent tracer of sSFR for star-forming galaxies out to redshifts of z=5. △ Less

Submitted 18 May, 2016; v1 submitted 5 November, 2015; originally announced November 2015.

Comments: Version accepted to publication in MNRAS, 11 pages, 5 figures, 1 table

New redshift z ~ 9 galaxies in the Hubble Frontier Fields: Implications for early evolution of the UV luminosity density

Authors: Derek J. McLeod, Ross J. McLure, James S. Dunlop, Brant E. Robertson, Richard S. Ellis, Thomas T. Targett

Abstract: We present the results of a new search for galaxies at redshift z ~ 9 in the first two Hubble Frontier Fields with completed HST WFC3/IR and ACS imaging. To ensure robust photometric redshift solutions, and to minimize incompleteness, we confine our search to objects with H_{160} < 28.6 (AB mag), consider only image regions with an rms noise sigma_{160} > 30 mag (within a 0.5-arcsec diameter apert… ▽ More We present the results of a new search for galaxies at redshift z ~ 9 in the first two Hubble Frontier Fields with completed HST WFC3/IR and ACS imaging. To ensure robust photometric redshift solutions, and to minimize incompleteness, we confine our search to objects with H_{160} < 28.6 (AB mag), consider only image regions with an rms noise sigma_{160} > 30 mag (within a 0.5-arcsec diameter aperture), and insist on detections in both H_{160} and J_{140}. The result is a survey covering an effective area (after accounting for magnification) of 10.9 sq. arcmin, which yields 12 galaxies at 8.4 < z < 9.5. Within the Abell-2744 cluster and parallel fields we confirm the three brightest objects reported by Ishigaki et al. (2014), but recover only one of the four z > 8.4 sources reported by Zheng et al. (2014). In the MACSJ0416.1-240 cluster field we report five objects, and explain why each of these eluded detection or classification as z ~ 9 galaxies in the published searches of the shallower CLASH data. Finally, we uncover four z ~ 9 galaxies from the previously unsearched MACSJ0416.1-240 parallel field. Based on the published magnification maps we find that only one of these 12 galaxies is likely boosted by more than a factor of two by gravitational lensing. Consequently we are able to perform a fairly straightforward reanalysis of the normalization of the z ~ 9 UV galaxy luminosity function as explored previously in the HUDF12 programme. We conclude that the new data strengthen the evidence for a continued smooth decline in UV luminosity density (and hence star-formation rate density) from z ~ 8 to z ~ 9, contrary to recent reports of a marked drop-off at these redshifts. This provides further support for the scenario in which early galaxy evolution is sufficiently extended to explain cosmic reionization. △ Less

Submitted 15 May, 2015; v1 submitted 3 December, 2014; originally announced December 2014.

Comments: 14 pages, 7 figures, updated to match accepted version in MNRAS