Showing 1–50 of 79 results for author: Cullen, F

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

Inferring the Distribution of the Ionising Photon Escape Fraction

Authors: Kimi C. Kreilgaard, Charlotte A. Mason, Fergus Cullen, Ryan Begley, Ross J. McLure

Abstract: The escape fraction of ionising photons from galaxies ($f_\mathrm{esc}$) is a key parameter for understanding how intergalactic hydrogen became reionised, but it remains mostly unconstrained. Measurements have been limited to the average value in galaxy ensembles and handfuls of individual detections. To help understand which mechanisms govern ionising photon escape, here we infer the distribution… ▽ More The escape fraction of ionising photons from galaxies ($f_\mathrm{esc}$) is a key parameter for understanding how intergalactic hydrogen became reionised, but it remains mostly unconstrained. Measurements have been limited to the average value in galaxy ensembles and handfuls of individual detections. To help understand which mechanisms govern ionising photon escape, here we infer the distribution of $f_\mathrm{esc}$. We develop a hierarchical Bayesian inference technique to estimate the population distribution of $f_\mathrm{esc}$ from the ratio of Lyman Continuum to non-ionising UV flux measured from broadband photometry. We apply it to a sample of 148 z ~ 3.5 star-forming galaxies from the VANDELS spectroscopic survey. We explore four physically motivated distributions: constant, log-normal, exponential and bimodal, recovering $\langle f_\mathrm{esc} \rangle \approx5\%$ for most models. We find the observations are best described by an exponential $f_\mathrm{esc}$ distribution with scale factor $μ=0.05^{+0.01}_{-0.02}$. This indicates most galaxies in our sample exhibit very low escape fractions while predicting substantial ionising photon leakage for only a few galaxies, implying a range of optical depths in the ISM and/or time variability in ionising photon escape. We rule out a bimodal distribution at high significance, indicating that a purely bimodal model of ionising photon escape (due to very strong sightline and/or time variability) is not favoured. We compare our recovered exponential distribution with the SPHINX simulations and find that, while the simulation also predicts an exponential-like distribution, it significantly underpredicts our inferred mean. The distribution of $f_\mathrm{esc}$ can be a vital test for simulations in understanding ionising photon leakage and is important to consider to gain a complete picture of reionisation. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 16 pages, 10 figures, main result in Figure 7, submitted to A&A

Comparing the VANDELS sample to a zoom-in Radiative Hydrodynamical Simulation: using the Si II and C II line spectra as tracers of galaxy evolution and Lyman Continuum leakage

Authors: Simon Gazagnes, Fergus Cullen, Valentin Mauerhofer, Ryan Begley, Danielle Berg, Jeremy Blaizot, John Chisholm, Thibault Garel, Floriane Leclercq, Ross J. McLure, Anne Verhamme

Abstract: We compare mock ultraviolet C II and Si II absorption and emission line features generated using a ~10$^9$ $M_\odot$ virtual galaxy with observations of 131 $z~3$ galaxies from the VANDELS survey. We find that the mock spectra reproduce reasonably well a large majority (83%) of the \vandels\ spectra ($χ^2<2$), but do not resemble the most massive objects ($>10^{10}M_\odot$) which exhibit broad abs… ▽ More We compare mock ultraviolet C II and Si II absorption and emission line features generated using a ~10$^9$ $M_\odot$ virtual galaxy with observations of 131 $z~3$ galaxies from the VANDELS survey. We find that the mock spectra reproduce reasonably well a large majority (83%) of the \vandels\ spectra ($χ^2<2$), but do not resemble the most massive objects ($>10^{10}M_\odot$) which exhibit broad absorption features. Interestingly, the best-matching mock spectra originate from periods of intense star formation in the virtual galaxy, where its luminosity is four times higher than in periods of relative quiescence. Furthermore, for each galaxy, we predict the Lyman Continuum (LyC) escape fractions using the environment of the virtual galaxy. We derive an average escape fraction of 0.01$\pm$0.02, consistent with other estimates from the literature. The predicted escape fractions are tightly correlated with the Lyman-$α$ escape fractions and highly consistent with observed empirical trends. Additionally, galaxies with larger predicted escape fractions exhibit bluer $β$ slopes, more Lyman-$α$ flux, and weaker low-ionization absorption lines. Building upon the good agreement between our predictions and observationally established LyC diagnostics, we examine the LyC leakage mechanisms in the simulation. We find that LyC photon leakage is enhanced in directions where the observed flux dominantly emerges from compact regions depleted of neutral gas and dust, mirroring the scenario inferred from observational data. In general, this study further highlights the potential of high-resolution radiation hydrodynamics simulations in analyzing UV absorption and emission line features and providing valuable insights into the LyC leakage of star-forming galaxies. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 15 pages, 10 figures, 0 birds. Accepted for publication in ApJ

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

The JWST-PRIMAL Legacy Survey. A JWST/NIRSpec reference sample for the physical properties and Lyman-$α$ absorption and emission of $\sim 500$ galaxies at $z=5.5-13.4$

Authors: K. E. Heintz, G. B. Brammer, D. Watson, P. A. Oesch, L. C. Keating, M. J. Hayes, Abdurro'uf, K. Z. Arellano-Córdova, A. C. Carnall, C. R. Christiansen, F. Cullen, R. Davé, P. Dayal, A. Ferrara, K. Finlator, J. P. U. Fynbo, S. R. Flury, V. Gelli, S. Gillman, R. Gottumukkala, K. Gould, T. R. Greve, S. E. Hardin, T. Y. -Y Hsiao, A. Hutter , et al. (23 additional authors not shown)

Abstract: One of the surprising early findings with JWST has been the discovery of a strong "roll-over" or a softening of the absorption edge of Ly$α$ in a large number of galaxies at ($z\gtrsim 6$), in addition to systematic offsets from photometric redshift estimates and fundamental galaxy scaling relations. This has been interpreted as damped Ly$α$ absorption (DLA) wings from high column densities of neu… ▽ More One of the surprising early findings with JWST has been the discovery of a strong "roll-over" or a softening of the absorption edge of Ly$α$ in a large number of galaxies at ($z\gtrsim 6$), in addition to systematic offsets from photometric redshift estimates and fundamental galaxy scaling relations. This has been interpreted as damped Ly$α$ absorption (DLA) wings from high column densities of neutral atomic hydrogen (HI), signifying major gas accretion events in the formation of these galaxies. To explore this new phenomenon systematically, we assemble the JWST/NIRSpec PRImordial gas Mass AssembLy (PRIMAL) legacy survey of 494 galaxies at $z=5.5-13.4$. We characterize this benchmark sample in full and spectroscopically derive the galaxy redshifts, metallicities, star-formation rates, and ultraviolet slopes. We define a new diagnostic, the Ly$α$ damping parameter $D_{\rm Lyα}$ to measure and quantify the Ly$α$ emission strength, HI fraction in the IGM, or local HI column density for each source. The JWST-PRIMAL survey is based on the spectroscopic DAWN JWST Archive (DJA-Spec). All the software, reduced spectra, and spectroscopically derived quantities and catalogs are made publicly available in dedicated repositories. The fraction of strong galaxy DLAs are found to be in the range $65-95\%$ at $z>5.5$. The fraction of strong Ly$α$ emitters (LAEs) is found to increase with decreasing redshift, in qualitative agreement with previous observational results, and are predominantly associated with low-metallicity and UV faint galaxies. By contrast, strong DLAs are observed in galaxies with a variety of intrinsic physical properties. Our results indicate that strong DLAs likely reflect a particular early assembly phase of reionization-era galaxies, at which point they are largely dominated by pristine HI gas accretion. [abridged] △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 20 pages, 18 figures, 3 tables. Submitted to A&A. Comments welcome! All data and catalogs are available through the DAWN JWST Archive (DJA): https://dawn-cph.github.io/dja/ and https://github.com/keheintz/jwst-primal

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

CLASSY IX: The Chemical Evolution of the Ne, S, Cl, and Ar Elements

Authors: Karla Z. Arellano-Córdova, Danielle A. Berg, Matilde Mingozzi, Bethan L. James, Noah S. J. Rogers, Evan D. Skillman, Fergus Cullen, Ryan Alexander, Ricardo O. Amorín, John Chisholm, Matthew Hayes, Timothy Heckman, Svean Hernandez, Nimisha Kumari, Claus Leitherer, Crystal L. Martin, Michael Maseda, Themiya Nanayakkara, Kaelee Parker, Swara Ravindranath, Alisson L. Strom, Fiorenzo Vincenzo, Aida Wofford

Abstract: To study the chemical evolution across cosmic epochs, we investigate Ne, S, Cl, and Ar abundance patterns in the COS Legacy Archive Spectroscopic SurveY (CLASSY). CLASSY comprises local star-forming galaxies (0.02 < z < 0.18) with enhanced star-formation rates, making them strong analogues to high-z star-forming galaxies. With direct measurements of electron temperature, we derive accurate ionic a… ▽ More To study the chemical evolution across cosmic epochs, we investigate Ne, S, Cl, and Ar abundance patterns in the COS Legacy Archive Spectroscopic SurveY (CLASSY). CLASSY comprises local star-forming galaxies (0.02 < z < 0.18) with enhanced star-formation rates, making them strong analogues to high-z star-forming galaxies. With direct measurements of electron temperature, we derive accurate ionic abundances for all elements and assess ionization correction factors (ICFs) to account for unseen ions and derive total abundances. We find Ne/O, S/O, Cl/O, and Ar/O exhibit constant trends with gas-phase metallicity for 12+log(O/H) < 8.5 but significant correlation for Ne/O and Ar/O with metallicity for 12+log(O/H) > 8.5, likely due to ICFs. Thus, applicability of the ICFs to integrated spectra of galaxies could bias results, underestimating true abundance ratios. Using CLASSY as a local reference, we assess the evolution of Ne/O, S/O, and Ar/O in galaxies at z>3, finding no cosmic evolution of Ne/O, while the lack of direct abundance determinations for S/O and Ar/O can bias the interpretation of the evolution of these elements. We determine the fundamental metallicity relationship (FMR) for CLASSY and compare to the high-redshift FMR, finding no evolution. Finally, we perform the first mass-neon relationship analysis across cosmic epochs, finding a slight evolution to high Ne at later epochs. The robust abundance patterns of CLASSY galaxies and their broad range of physical properties provide essential benchmarks for interpreting the chemical enrichment of the early galaxies observed with the JWST. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: Accepted for publication in The Astrophysical Journal. 20 pages (main body), 10 figures, 6 Tables

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

Two Distinct Classes of Quiescent Galaxies at Cosmic Noon Revealed by JWST PRIMER and UNCOVER

Authors: Sam E. Cutler, Katherine E. Whitaker, John R. Weaver, Bingjie Wang, Richard Pan, Rachel Bezanson, Lukas J. Furtak, Ivo Labbe, Joel Leja, Sedona H. Price, Yingjie Cheng, Maike Clausen, Fergus Cullen, Pratika Dayal, Anna de Graaff, Mark Dickinson, James S. Dunlop, Robert Feldmann, Marijn Franx, Mauro Giavalisco, Karl Glazebrook, Jenny E. Greene, Norman A. Grogin, Garth Illingworth, Anton M. Koekemoer , et al. (9 additional authors not shown)

Abstract: We present a measurement of the low-mass quiescent size-mass relation at Cosmic Noon (1<z<3) from the JWST PRIMER and UNCOVER treasury surveys, which highlights two distinct classes of quiescent galaxies. While the massive population is well studied at these redshifts, the low-mass end has been previously under-explored due to a lack of observing facilities with sufficient sensitivity and spatial… ▽ More We present a measurement of the low-mass quiescent size-mass relation at Cosmic Noon (1<z<3) from the JWST PRIMER and UNCOVER treasury surveys, which highlights two distinct classes of quiescent galaxies. While the massive population is well studied at these redshifts, the low-mass end has been previously under-explored due to a lack of observing facilities with sufficient sensitivity and spatial resolution. We select a conservative sample of low-mass quiescent galaxy candidates using rest-frame UVJ colors and specific star formation rate criteria and measure galaxy morphology in both rest-frame UV/optical wavelengths (F150W) and rest-frame near-infrared (F444W). We confirm an unambiguous flattening of the low-mass quiescent size-mass relation, which results from the separation of the quiescent galaxy sample into two distinct populations at $\log(M_\star/M_\odot)\sim10.3$: low-mass quiescent galaxies that are notably younger and have disky structures, and massive galaxies consistent with spheroidal morphologies and older median stellar ages. These separate populations imply mass quenching dominates at the massive end while other mechanisms, such as environmental or feedback-driven quenching, form the low-mass end. This stellar mass dependent slope of the quiescent size-mass relation could also indicate a shift from size growth due to star formation (low masses) to growth via mergers (massive galaxies). The transition mass between these two populations also corresponds with other dramatic changes and characteristic masses in several galaxy evolution scaling relations (e.g. star-formation efficiency, dust obscuration, and stellar-halo mass ratios), further highlighting the stark dichotomy between low-mass and massive galaxy formation. △ Less

Submitted 23 April, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

Comments: 16 pages, 6 figures, 1 table. Submitted to ApJL. Revised 2024 April 18

Quasar Sightline and Galaxy Evolution (QSAGE) -- III. The mass-metallicity and fundamental metallicity relation of $z \approx$ 2.2 galaxies

Authors: H. M. O. Stephenson, J. P. Stott, F. Cullen, R. M. Bielby, N. Amos, R. Dutta, M. Fumagalli, N. Tejos, J. N. Burchett, R. A. Crain, J. X. Prochaska

Abstract: We present analysis of the mass-metallicity relation (MZR) for a sample of 67 [OIII]-selected star-forming galaxies at a redshift range of $z=1.99 - 2.32$ ($z_{\text{med}} = 2.16$) using \emph{Hubble Space Telescope} Wide Field Camera 3 grism spectroscopy from the Quasar Sightline and Galaxy Evolution (QSAGE) survey. Metallicities were determined using empirical gas-phase metallicity calibrations… ▽ More We present analysis of the mass-metallicity relation (MZR) for a sample of 67 [OIII]-selected star-forming galaxies at a redshift range of $z=1.99 - 2.32$ ($z_{\text{med}} = 2.16$) using \emph{Hubble Space Telescope} Wide Field Camera 3 grism spectroscopy from the Quasar Sightline and Galaxy Evolution (QSAGE) survey. Metallicities were determined using empirical gas-phase metallicity calibrations based on the strong emission lines [OII]3727,3729, [OIII]4959,5007 and H$β$. Star-forming galaxies were identified, and distinguished from active-galactic nuclei, via Mass-Excitation diagrams. Using $z\sim0$ metallicity calibrations, we observe a negative offset in the $z=2.2$ MZR of $\approx -0.51$ dex in metallicity when compared to locally derived relationships, in agreement with previous literature analysis. A similar offset of $\approx -0.46$ dex in metallicity is found when using empirical metallicity calibrations that are suitable out to $z\sim5$, though our $z=2.2$ MZR, in this case, has a shallower slope. We find agreement between our MZR and those predicted from various galaxy evolution models and simulations. Additionally, we explore the extended fundamental metallicity relation (FMR) which includes an additional dependence on star formation rate (SFR). Our results consistently support the existence of the FMR, as well as revealing an offset of $0.28\pm0.04$ dex in metallicity compared to locally-derived relationships, consistent with previous studies at similar redshifts. We interpret the negative correlation with SFR at fixed mass, inferred from an FMR existing for our sample, as being caused by the efficient accretion of metal-poor gas fuelling SFR at cosmic noon. △ Less

Submitted 12 December, 2023; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: 15 pages, 10 figures. This paper has been accepted for publication in 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

Disappearing galaxies: the orientation dependence of JWST-bright, HST-dark, star-forming galaxy selection

Authors: R. K. Cochrane, D. Anglés-Alcázar, F. Cullen, C. C. Hayward

Abstract: Galaxies that are invisible in deep optical-NIR imaging but detected at longer wavelengths have been the focus of several recent observational studies, with speculation that they could constitute a substantial missing population and even dominate the cosmic star formation rate density at $z\gtrsim4$. The depths now achievable with JWST at the longest wavelengths probed by HST, coupled with the tra… ▽ More Galaxies that are invisible in deep optical-NIR imaging but detected at longer wavelengths have been the focus of several recent observational studies, with speculation that they could constitute a substantial missing population and even dominate the cosmic star formation rate density at $z\gtrsim4$. The depths now achievable with JWST at the longest wavelengths probed by HST, coupled with the transformative resolution at longer wavelengths, are already enabling detailed, spatially-resolved characterisation of sources that were invisible to HST, often known as `HST-dark' galaxies. However, until now, there has been little theoretical work to compare against. We present the first simulation-based study of this population, using highly-resolved galaxies from the Feedback in Realistic Environments (FIRE) project, with multi-wavelength images along several lines of sight forward-modelled using radiative transfer. We naturally recover a population of modelled sources that meet commonly-used selection criteria ($H_{\rm{AB}}>27\,\rm{mag}$ and $H_{\rm{AB}}-\rm{F444W}>2.3$). These simulated HST-dark galaxies lie at high redshifts ($z=4-7$), have high levels of dust attenuation ($A_{V}=2-4$), and display compact recent star formation ($R_{1/2,\,\rm{4.4\,μ\rm{m}}}\lesssim1\,\rm{kpc}$). Orientation is very important: for all but one of the 17 simulated galaxy snapshots with HST-dark sightlines, there exist other sightlines that do not meet the criteria. This result has important implications for comparisons between observations and models that do not resolve the detailed star-dust geometry, such as semi-analytic models or coarsely-resolved hydrodynamical simulations. Critically, we demonstrate that HST-dark sources are not an unexpected or exotic population, but a subset of high-redshift, highly-dust-attenuated sources viewed along certain lines of sight. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 12 pages, 8 figures. Accepted for publication in ApJ

The ALMA REBELS survey: obscured star formation in massive Lyman-break galaxies at z = 4-8 revealed by the IRX-$β$ and $M_{\star}$ relations

Authors: R. A. A. Bowler, H. Inami, L. Sommovigo, R. Smit, H. S. B. Algera, M. Aravena, L. Barrufet, R. Bouwens, E. da Cunha, F. Cullen, P. Dayal, I. de Looze, J. S. Dunlop, Y. Fudamoto, V. Mauerhofer, R. J. McLure, M. Stefanon, R. Schneider, A. Ferrara, L. Graziani, J. A. Hodge, T. Nanayakkara, M. Palla, S. Schouws, D. P. Stark , et al. (1 additional authors not shown)

Abstract: We investigate the degree of dust obscured star formation in 49 massive (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot})>9$) Lyman-break galaxies (LBGs) at $z = 6.5$-$8$ observed as part of the ALMA Reionization Era Bright Emission Line Survey (REBELS) large program. By creating deep stacks of the photometric data and the REBELS ALMA measurements we determine the average rest-frame UV, optical and far-… ▽ More We investigate the degree of dust obscured star formation in 49 massive (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot})>9$) Lyman-break galaxies (LBGs) at $z = 6.5$-$8$ observed as part of the ALMA Reionization Era Bright Emission Line Survey (REBELS) large program. By creating deep stacks of the photometric data and the REBELS ALMA measurements we determine the average rest-frame UV, optical and far-infrared (FIR) properties which reveal a significant fraction ($f_{\rm obs} = 0.4$-$0.7$) of obscured star formation, consistent with previous studies. From measurements of the rest-frame UV slope, we find that the brightest LBGs at these redshifts show bluer ($β\simeq -2.2$) colours than expected from an extrapolation of the colour-magnitude relation found at fainter magnitudes. Assuming a modified blackbody spectral-energy distribution (SED) in the FIR (with dust temperature of $T_{\rm d} = 46\,{\rm K}$ and $β_{\rm d} = 2.0$), we find that the REBELS sources are in agreement with the local ''Calzetti-like'' starburst Infrared-excess (IRX)-$β$ relation. By reanalysing the data available for 108 galaxies at $z \simeq 4$-$6$ from the ALPINE ALMA large program using a consistent methodology and assumed FIR SED, we show that from $z \simeq 4$-$8$, massive galaxies selected in the rest-frame UV have no appreciable evolution in their derived IRX-$β$ relation. When comparing the IRX-$M_{\star}$ relation derived from the combined ALPINE and REBELS sample to relations established at $z < 4$, we find a deficit in the IRX, indicating that at $z > 4$ the proportion of obscured star formation is lower by a factor of $\gtrsim 3$ at a given a $M_{\star}$. Our IRX-$β$ results are in good agreement with the high-redshift predictions of simulations and semi-analytic models for $z \simeq 7$ galaxies with similar stellar masses and SFRs. △ Less

Submitted 28 November, 2023; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 18 pages, 10 figures, 2 tables (plus 1 figure and 2 tables in the appendix). Updated to match MNRAS accepted version after minor corrections

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)

JWST observations of galaxy damping wings during reionization interpreted with cosmological simulations

Authors: Laura C. Keating, James S. Bolton, Fergus Cullen, Martin G. Haehnelt, Ewald Puchwein, Girish Kulkarni

Abstract: Spectra of the highest redshift galaxies taken with JWST are now allowing us to see into the heart of the reionization epoch. Many of these observed galaxies exhibit strong damping wing absorption redward of their Lyman-$α$ emission. These observations have been used to measure the redshift evolution of the neutral fraction of the intergalactic medium and sizes of ionized bubbles. However, these e… ▽ More Spectra of the highest redshift galaxies taken with JWST are now allowing us to see into the heart of the reionization epoch. Many of these observed galaxies exhibit strong damping wing absorption redward of their Lyman-$α$ emission. These observations have been used to measure the redshift evolution of the neutral fraction of the intergalactic medium and sizes of ionized bubbles. However, these estimates have been made using a simple analytic model for the intergalactic damping wing. We explore the recent observations with models of inhomogeneous reionization from the Sherwood-Relics simulation suite. We carry out a comparison between the damping wings calculated from the simulations and from the analytic model. We find that although the agreement is good on the red side of the Lyman-$α$ emission, there is a discrepancy on the blue side due to residual neutral hydrogen present in the simulations, which saturates the intergalactic absorption. For this reason, we find that it is difficult to reproduce the claimed observations of large bubble sizes at z ~ 7, which are driven by a detection of transmitted flux blueward of the Lyman-$α$ emission. We suggest instead that the observations can be explained by a model with smaller ionized bubbles and larger intrinsic Lyman-$α$ emission from the host galaxy. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: 13 pages, 7 figures. Submitted to MNRAS

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

Extreme broad He\2 emission at high and low redshifts: the dominant role of VMS in NGC 3125-A1 and CDFS131717

Authors: Aida Wofford, Andrés Sixtos, Stephane Charlot, Gustavo Bruzual, Fergus Cullen, Thomas M. Stanton, Svea Hernández, Linda J. Smith, Matthew Hayes

Abstract: Super star cluster (SSC) A1 (3.1E5 Msun) in NGC 3125 has one of the strongest (EW = 4.6 +/- 0.5 Ang) broad (FWHM = 1131 +\- 40 km/s) He II 1640 emission lines in the nearby Universe and constitutes an important template for interpreting observations of extreme He II emitters out to redshifts of z = 2-3. We use Cosmic Origins Spectrograph (COS) observations of A1 to show that there is no significan… ▽ More Super star cluster (SSC) A1 (3.1E5 Msun) in NGC 3125 has one of the strongest (EW = 4.6 +/- 0.5 Ang) broad (FWHM = 1131 +\- 40 km/s) He II 1640 emission lines in the nearby Universe and constitutes an important template for interpreting observations of extreme He II emitters out to redshifts of z = 2-3. We use Cosmic Origins Spectrograph (COS) observations of A1 to show that there is no significant contamination of the He II line with nebular emission and that the line is redshifted by 121 +/-17 km/s relative to ISM lines. We compare the COS G130M + G160M observations of A1 to recent binary BPASS and single-star Charlot & Bruzual (C&B) simple stellar population (SSP) models with Very Massive Stars (VMS) of up to 300 Msun. We suggest why BPASS models fail to reproduce A1's He II emission. On the other hand, a C&B model with Z = 0.008, age = 2.2 Myr, and VMS approaching the Eddington limit provides an excellent fit to the He II emission and fits reasonably well C III 1175, N V 1238,1241, and C IV 1548, 1551. We present O V 1371 line-profile predictions showing that this line constitutes an important tracer of youth and VMS in galaxies. Finally, we discuss the presence of VMS in CDFS131717, a highly star-forming low-metallicity galaxy located at z = 3.071, which has a tentative detection of O V absorption and strong broad He II emission. These features are rare and hint to the presence of short-lived VMS in the galaxy. Our results show the effect of the latest developments of stellar wind theory and the importance of accounting for VMS in models. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 17 pages, 17 figures, accepted in MNRAS

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

Insights into the reionization epoch from cosmic-noon-CIV emitters in the VANDELS survey

Authors: S. Mascia, L. Pentericci, A. Saxena, D. Belfiori, A. Calabrò, M. Castellano, A. Saldana-Lopez, M. Talia, R. Amorín, F. Cullen, B. Garilli, L. Guaita, M. Llerena, R. J. McLure, M. Moresco, P. Santini, D. Schaerer

Abstract: Recently, intense emission from nebular C III] and C IV emission lines have been observed in galaxies in the epoch of reionization ($z>6$) and have been proposed as the prime way of measuring their redshift and studying their stellar populations. These galaxies might represent the best examples of cosmic reionizers, as suggested by recent low-z observations of Lyman Continuum emitting galaxies, bu… ▽ More Recently, intense emission from nebular C III] and C IV emission lines have been observed in galaxies in the epoch of reionization ($z>6$) and have been proposed as the prime way of measuring their redshift and studying their stellar populations. These galaxies might represent the best examples of cosmic reionizers, as suggested by recent low-z observations of Lyman Continuum emitting galaxies, but it is hard to directly study the production and escape of ionizing photons at such high redshifts. The ESO spectroscopic public survey VANDELS offers the unique opportunity to find rare examples of such galaxies at cosmic noon ($z\sim 3$), thanks to the ultra deep observations available. We have selected a sample of 39 galaxies showing C IV emission, whose origin (after a careful comparison to photoionization models) can be ascribed to star formation and not to AGN. By using a multi-wavelength approach, we determine their physical properties including metallicity and ionization parameter and compare them to the properties of the parent population to understand what are the ingredients that could characterize the analogs of the cosmic reionizers. We find that C IV emitters are galaxies with high photons production efficiency and there are strong indications that they might have also large escape fraction: given the visibility of C IV in the epoch of reionization this could become the best tool to pinpoint the cosmic reioinzers. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 15 pages, 10 figures, resubmitted to A&A after addressing the referee comments

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

Being KLEVER at cosmic noon: ionised gas outflows are inconspicuous in low-mass star-forming galaxies but prominent in massive AGN hosts

Authors: Alice Concas, Roberto Maiolino, Mirko Curti, Connor Hayden-Pawson, Michele Cirasuolo, Gareth C. Jones, Amata Mercurio, Francesco Belfiore, Giovanni Cresci, Fergus Cullen, Filippo Mannucci, Alessandro Marconi, Michele Cappellari, Claudia Cicone, Yingjie Peng, Paulina Troncoso

Abstract: We investigate the presence of ionised gas outflows in a sample of 141 main-sequence star-forming galaxies at $1.2<z<2.6$ from the KLEVER (KMOS Lensed Emission Lines and VElocity Review) survey. Our sample covers an exceptionally wide range of stellar masses, $8.1<\log(M_\star/M_{\odot})<11.3$, pushing outflow studies into the dwarf regime thanks to gravitationally lensed objects. We stack optical… ▽ More We investigate the presence of ionised gas outflows in a sample of 141 main-sequence star-forming galaxies at $1.2<z<2.6$ from the KLEVER (KMOS Lensed Emission Lines and VElocity Review) survey. Our sample covers an exceptionally wide range of stellar masses, $8.1<\log(M_\star/M_{\odot})<11.3$, pushing outflow studies into the dwarf regime thanks to gravitationally lensed objects. We stack optical rest-frame emission lines (H$β$, [OIII], H$α$ and [NII]) in different mass bins and seek for tracers of gas outflows by using a novel, physically motivated method that improves over the widely used, simplistic double Gaussian fitting. We compare the observed emission lines with the expectations from a rotating disc (disc+bulge for the most massive galaxies) model, whereby significant deviations are interpreted as a signature of outflows. We find clear evidence for outflows in the most massive, $\log(M_\star/M_{\odot}) > 10.8$, AGN-dominated galaxies, suggesting that AGNs may be the primary drivers of these gas flows. Surprisingly, at $\log(M_\star/M_{\odot})\leq 9.6$, the observed line profiles are fully consistent with a rotating disc model, indicating that ionised gas outflows in dwarf galaxies might play a negligible role even during the peak of cosmic star-formation activity. Finally, we find that the observed mass loading factor scales with stellar mass as expected from the TNG50 cosmological simulation, but the ionised gas mass accounts for only 2$\%$ of the predicted value. This suggests that either the bulk of the outflowing mass is in other gaseous phases or the current feedback models implemented in cosmological simulations need to be revised. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS, 28 pages, 20 figures

The environmental dependence of the stellar and gas-phase mass-metallicity relation at 2 < z < 4

Authors: A. Calabro, L. Guaita, L. Pentericci, F. Fontanot, M. Castellano, G. De Lucia, T. Garofalo, P. Santini, F. Cullen, A. Carnall, B. Garilli, M. Talia, G. Cresci, M. Franco, J. P. U. Fynbo, N. P. Hathi, M. Hirschmann, A. Koekemoer, M. Llerena, L. Xie

Abstract: In the local universe, galaxies in clusters show different properties compared to more isolated systems. Understanding how this difference originates and whether it is already in place at high redshift is still a matter of debate. Thanks to uniquely deep optical spectra from the VANDELS survey, we investigate environmental effects on the stellar mass-metallicity relation (MZR) for a sample of ~100… ▽ More In the local universe, galaxies in clusters show different properties compared to more isolated systems. Understanding how this difference originates and whether it is already in place at high redshift is still a matter of debate. Thanks to uniquely deep optical spectra from the VANDELS survey, we investigate environmental effects on the stellar mass-metallicity relation (MZR) for a sample of ~1000 star-forming galaxies at 2<z<4. We complement our dataset with MOSFIRE follow-up of 21 galaxies to study the environmental dependence of the gas-phase MZR. Robust stellar and gas metallicities are derived, respectively, from well-calibrated photospheric absorptions features at 1501 and 1719 Åin the stacked spectra, and from optical emission lines ([OII]3726-3729, [OIII]5007, and Hbeta) in individual systems. We characterize the environment through multiple criteria by using the local galaxy density maps previously derived in VANDELS. We find that environmental effects are weak at these redshifts, and more important around the densest overdensity structures, where galaxies have a lower stellar Z (by 0.2 dex) and a lower gas-phase Z (by 0.1 dex) compared to the field, with a significance of 1 and 2 sigma, respectively. Crucially, this offset cannot be explained by a selection effect due to a higher SFR, a fainter UV continuum, or different dust attenuations and stellar ages. Despite the still low S/N of our results, we propose a combination of increased mergers and high-speed encounters, more efficient AGN feedback in dense cores, and cold gas inflows as viable mechanisms diluting the metal content of overdense galaxies or expelling their metals to the IGM. Finally, some tensions remain with semi-analytic models and hydrodynamical simulations, which predict no significant offset as a function of host halo mass, suggesting that an explicit implementation of environmental processes is needed. △ Less

Submitted 9 March, 2022; originally announced March 2022.

Comments: Accepted for publication in A&A ; 22 pages, 15 figures, 2 tables, and 1 Appendix

Journal ref: A&A 664, A75 (2022)

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 discovery of rest-frame UV colour gradients and a diversity of dust morphologies in bright z ~ 7 Lyman-break galaxies

Authors: R. A. A. Bowler, F. Cullen, R. J. McLure, J. S. Dunlop, A. Avison

Abstract: We present deep ALMA dust continuum observations for a sample of luminous ($M_{\rm UV} < -22$) star-forming galaxies at $z \simeq 7$. We detect five of the six sources in the far-infrared (FIR), providing key constraints on the obscured star-formation rate (SFR) and the infrared-excess-$β$ (IRX-$β$) relation without the need for stacking. Despite the galaxies showing blue rest-frame UV slopes (… ▽ More We present deep ALMA dust continuum observations for a sample of luminous ($M_{\rm UV} < -22$) star-forming galaxies at $z \simeq 7$. We detect five of the six sources in the far-infrared (FIR), providing key constraints on the obscured star-formation rate (SFR) and the infrared-excess-$β$ (IRX-$β$) relation without the need for stacking. Despite the galaxies showing blue rest-frame UV slopes ($β\simeq -2$) we find that 35-75 percent of the total SFR is obscured. We find the IRX-$β$ relation derived for these $z \simeq 7$ sources is consistent with that found for local star-burst galaxies. Using our relatively high-resolution (FWHM $\simeq 0.7\,{\rm arcsec}$) observations we identify a diversity of dust morphologies in the sample. We find both compact emission that appears offset relative to the unobscured components and extended dust emission that is co-spatial with the rest-frame UV light. In the majority of the sources we detect strong rest-frame UV colour gradients (with up to $Δβ\simeq 0.7$-$1.4$) as probed by the multi-band UltraVISTA ground-based data. The observed redder colours are spatially correlated with the location of the FIR detection. Our results show that even in bright Lyman-break galaxies at $z \simeq 7$ the peak of the star-formation is typically hosted by the fainter, redder, regions in the rest-frame UV, which have an obscured fraction of $f_{\rm obs} \ge 0.8$. As well as demonstrating the importance of dust obscured star-formation within the Epoch of Reionization, these observations provide an exciting taster of the rich spatially resolved datasets that will be obtained from JWST and high-resolution ALMA follow-up at these redshifts. △ Less

Submitted 5 February, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: 14 pages, 6 figures, 3 tables. Updated to MNRAS accepted version

The KLEVER survey: Nitrogen abundances at $z\sim$2 and probing the existence of a fundamental nitrogen relation

Authors: Connor Hayden-Pawson, Mirko Curti, Roberto Maiolino, Michele Cirasuolo, Francesco Belfiore, Michele Cappellari, Alice Concas, Giovanni Cresci, Fergus Cullen, Chiaki Kobayashi, Filippo Mannucci, Alessandro Marconi, Massimo Meneghetti, Amata Mercurio, Yingjie Peng, Mark Swinbank, Fiorenzo Vincenzo

Abstract: We present a comparison of the nitrogen-to-oxygen ratio (N/O) in 37 high-redshift galaxies at $z\sim$2 taken from the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey with a comparison sample of local galaxies, taken from the Sloan Digital Sky Survey (SDSS). The KLEVER sample shows only a mild enrichment in N/O of $+$0.1 dex when compared to local galaxies at a given gas-phase metall… ▽ More We present a comparison of the nitrogen-to-oxygen ratio (N/O) in 37 high-redshift galaxies at $z\sim$2 taken from the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey with a comparison sample of local galaxies, taken from the Sloan Digital Sky Survey (SDSS). The KLEVER sample shows only a mild enrichment in N/O of $+$0.1 dex when compared to local galaxies at a given gas-phase metallicity (O/H), but shows a depletion in N/O of $-$0.36 dex when compared at a fixed stellar mass (M$_*$). We find a strong anti-correlation in local galaxies between N/O and SFR in the M$_*$-N/O plane, similar to the anti-correlation between O/H and SFR found in the mass-metallicity relation (MZR). We use this anti-correlation to construct a fundamental nitrogen relation (FNR), analogous to the fundamental metallicity relation (FMR). We find that KLEVER galaxies are consistent with both the FMR and the FNR. This suggests that the depletion of N/O in high-$z$ galaxies when considered at a fixed M$_*$ is driven by the redshift-evolution of the mass-metallicity relation in combination with a near redshift-invariant N/O-O/H relation. Furthermore, the existence of an fundamental nitrogen relation suggests that the mechanisms governing the fundamental metallicity relation must be probed by not only O/H, but also N/O, suggesting pure-pristine gas inflows are not the primary driver of the FMR, and other properties such as variations in galaxy age and star formation efficiency must be important. △ Less

Submitted 9 March, 2022; v1 submitted 30 September, 2021; originally announced October 2021.

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

The MOSDEF Survey: Implications of the Lack of Evolution in the Dust Attenuation-Mass Relation to z~2

Authors: Alice E. Shapley, Ryan L. Sanders, Samir Salim, Naveen A. Reddy, Mariska Kriek, Bahram Mobasher, Alison Coil, Brian Siana, Sedona H. Price, Irene Shivaei, James S. Dunlop, Ross J. McLure, Fergus Cullen

Abstract: We investigate the relationship between dust attenuation and stellar mass ($M_*$) in star-forming galaxies over cosmic time. For this analysis, we compare measurements from the MOSFIRE Deep Evolution Field (MOSDEF) survey at $z\sim2.3$ and the Sloan Digital Sky Survey (SDSS) at $z\sim0$, augmenting the latter optical dataset with both UV Galaxy Evolution Explorer (GALEX) and mid-infrared Wide-fiel… ▽ More We investigate the relationship between dust attenuation and stellar mass ($M_*$) in star-forming galaxies over cosmic time. For this analysis, we compare measurements from the MOSFIRE Deep Evolution Field (MOSDEF) survey at $z\sim2.3$ and the Sloan Digital Sky Survey (SDSS) at $z\sim0$, augmenting the latter optical dataset with both UV Galaxy Evolution Explorer (GALEX) and mid-infrared Wide-field Infrared Survey Explorer (WISE) photometry from the GALEX-SDSS-WISE Catalog. We quantify dust attenuation using both spectroscopic measurements of H$α$ and H$β$ emission lines, and photometric measurements of the rest-UV stellar continuum. The H$α$/H$β$ ratio is used to determine the magnitude of attenuation at the wavelength of H$α$, $A_{{\rm H}α}$. Rest-UV colors and spectral-energy-distribution fitting are used to estimate $A_{1600}$, the magnitude of attenuation at a rest wavelength of 1600Å. As in previous work, we find a lack of significant evolution in the relation between dust attenuation and $M_*$ over the redshift range $z\sim0$ to $z\sim2.3$. Folding in the latest estimates of the evolution of $M_{\rm dust}$, $({M_{\rm dust}}/{M_{\rm gas}})$, and gas surface density at fixed $M_*$, we find that the expected $M_{\rm dust}$ and dust mass surface density are both significantly higher at $z\sim2.3$ than at $z\sim0$. These differences appear at odds with the lack of evolution in dust attenuation. To explain the striking constancy in attenuation vs. $M_*$, it is essential to determine the relationship between metallicity and $({M_{\rm dust}}/{M_{\rm gas}})$, the dust mass absorption coefficient, and dust geometry, and the evolution of these relations and quantities from $z\sim0$ to $z\sim2.3$. △ Less

Submitted 6 January, 2022; v1 submitted 29 September, 2021; originally announced September 2021.

Comments: 14 pages, 4 figures, accepted to ApJ

Metal-enriched halo gas across galaxy overdensities over the last 10 billion years

Authors: Rajeshwari Dutta, Michele Fumagalli, Matteo Fossati, Richard M. Bielby, John P. Stott, Emma K. Lofthouse, Sebastiano Cantalupo, Fergus Cullen, Robert A. Crain, Todd M. Tripp, J. Xavier Prochaska, Fabrizio Arrigoni Battaia, Joseph N. Burchett, Johan P. U. Fynbo, Michael T. Murphy, Joop Schaye, Nicolas Tejos, Tom Theuns

Abstract: We present a study of metal-enriched halo gas traced by MgII and CIV absorption at z<2 in the MUSE Analysis of Gas around Galaxies survey and the Quasar Sightline and Galaxy Evolution survey. Using these large and complete galaxy surveys in quasar fields, we study the dependence of the metal distribution on galaxy properties and overdensities, out to physical projected separations of 750 kpc. We f… ▽ More We present a study of metal-enriched halo gas traced by MgII and CIV absorption at z<2 in the MUSE Analysis of Gas around Galaxies survey and the Quasar Sightline and Galaxy Evolution survey. Using these large and complete galaxy surveys in quasar fields, we study the dependence of the metal distribution on galaxy properties and overdensities, out to physical projected separations of 750 kpc. We find that the cool, low-ionization gas is significantly affected by the environment across the full redshift range probed, with ~2-3 times more prevalent and stronger MgII absorption in higher overdensity group environments and in regions with greater overall stellar mass and star formation rates. Complementary to these results, we have further investigated the more highly ionized gas as traced by CIV absorption, and found that it is likely to be more extended than the MgII gas, with ~2 times higher covering fraction at a given distance. We find that the strength and covering fraction of CIV absorption show less significant dependence on galaxy properties and environment than the MgII absorption, but more massive and star-forming galaxies nevertheless also show ~2 times higher incidence of CIV absorption. The incidence of MgII and CIV absorption within the virial radius shows a tentative increase with redshift, being higher by a factor of ~1.5 and ~4, respectively, at z>1. It is clear from our results that environmental processes have a significant impact on the distribution of metals around galaxies and need to be fully accounted for when analyzing correlations between gaseous haloes and galaxy properties. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: 27 pages; accepted for publication in 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 evolution of the mass-metallicity relations from the VANDELS survey and the GAEA Semi-Analytic model

Authors: Fabio Fontanot, Antonello Calabrò, Margherita Talia, Filippo Mannucci, Marco Castellano, Giovanni Cresci, Gabriella De Lucia, Anna Gallazzi, Michaela Hirschmann, Laura Pentericci, Lizhi Xie, Ricardo Amorin, Micol Bolzonella, Angela Bongiorno, Olga Cucciati, Fergus Cullen, Johan P. U. Fynbo, Nimish Hathi, Pascale Hibon, Ross J. McLure, Lucia Pozzetti

Abstract: In this work, we study the evolution of the mass-metallicity relations (MZRs) as predicted by the GAlaxy Evolution and Assembly (GAEA) semi-analytic model. We contrast these predictions with recent results from the VANDELS survey, that allows us to expand the accessible redshift range for the stellar MZR up to $z\sim3.5$. We complement our study by considering the evolution of the gas-phase MZR in… ▽ More In this work, we study the evolution of the mass-metallicity relations (MZRs) as predicted by the GAlaxy Evolution and Assembly (GAEA) semi-analytic model. We contrast these predictions with recent results from the VANDELS survey, that allows us to expand the accessible redshift range for the stellar MZR up to $z\sim3.5$. We complement our study by considering the evolution of the gas-phase MZR in the same redshift range. We show that GAEA is able to reproduce the observed evolution of the $z<3.5$ gas-phase MZR and $z<0.7$ stellar MZR, while it overpredicts the stellar metallicity at $z\sim3.5$. Furthermore, GAEA also reproduces the so-called fundamental metallicity relation (FMR) between gas-phase metallicity, stellar mass and star formation rate (SFR). In particular, the gas-phase FMR in GAEA is already in place at $z\sim5$ and shows almost no evolution at lower redshift. GAEA predicts the existence of a stellar FMR, that is, however, characterized by a relevant redshift evolution, although its shape follows closely the gas-phase FMR. We also report additional unsolved tensions between model and data: the overall normalization of the predicted MZR agrees with observations only within $\sim$0.1 dex; the largest discrepancies are seen at $z\sim3.5$ where models tend to slightly overpredict observed metallicities; the slope of the predicted MZR at fixed SFR is too steep below a few ${\rm M}_\odot\ {\rm yr}^{-1}$. Finally, we provide model predictions for the evolution of the MZRs at higher redshifts, that would be useful in the context of future surveys, like those that will be performed with JWST. △ Less

Submitted 26 April, 2021; v1 submitted 16 April, 2021; originally announced April 2021.

Comments: 12 pages, 8 figures, MNRAS accepted, replaced to match the published version

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 VANDELS ESO public spectroscopic survey: final Data Release of 2087 spectra and spectroscopic measurements

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

Abstract: VANDELS is an ESO Public Spectroscopic Survey designed to build a sample of high signal to noise, medium resolution spectra of galaxies at redshift between 1 and 6.5. Here we present the final Public Data Release of the VANDELS Survey, comprising 2087 redshift measurements. We give a detailed description of sample selection, observations and data reduction procedures. The final catalogue reaches a… ▽ More VANDELS is an ESO Public Spectroscopic Survey designed to build a sample of high signal to noise, medium resolution spectra of galaxies at redshift between 1 and 6.5. Here we present the final Public Data Release of the VANDELS Survey, comprising 2087 redshift measurements. We give a detailed description of sample selection, observations and data reduction procedures. The final catalogue reaches a target selection completeness of 40% at iAB = 25. The high Signal to Noise ratio of the spectra (above 7 in 80% of the spectra) and the dispersion of 2.5Å allowed us to measure redshifts with high precision, the redshift measurement success rate reaching almost 100%. Together with the redshift catalogue and the reduced spectra, we also provide optical mid-IR photometry and physical parameters derived through SED fitting. The observed galaxy sample comprises both passive and star forming galaxies covering a stellar mass range 8.3< Log(M*/Msolar)<11.7. All catalogues and spectra are accessible through the survey database (http://vandels.inaf.it) where all information can be queried interactively, and via the ESO Archive (https://www.eso.org/qi/). △ Less

Submitted 19 January, 2021; originally announced January 2021.

Comments: 16 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 647, A150 (2021)

The VANDELS survey: the relation between UV continuum slope and stellar metallicity in star-forming galaxies at z~3

Authors: A. Calabrò, M. Castellano, L. Pentericci, F. Fontanot, N. Menci, F. Cullen, R. McLure, M. Bolzonella, A. Cimatti, F. Marchi, M. Talia, R. Amorín, G. Cresci, G. De Lucia, J. Fynbo, A. Fontana, M. Franco, N. P. Hathi, P. Hibon, M. Hirschmann, F. Mannucci, P. Santini, A. Saxena, D. Schaerer, L. Xie , et al. (1 additional authors not shown)

Abstract: The estimate of stellar metallicities (Z*) of high-z galaxies are of paramount importance in order to understand the complexity of dust effects and the reciprocal interrelations among stellar mass, dust attenuation, stellar age, and metallicity. Benefiting from uniquely deep FUV spectra of >500 star-forming galaxies at redshifts 2<z<5 extracted from the VANDELS survey and stacked in bins of stella… ▽ More The estimate of stellar metallicities (Z*) of high-z galaxies are of paramount importance in order to understand the complexity of dust effects and the reciprocal interrelations among stellar mass, dust attenuation, stellar age, and metallicity. Benefiting from uniquely deep FUV spectra of >500 star-forming galaxies at redshifts 2<z<5 extracted from the VANDELS survey and stacked in bins of stellar mass (M*) and UV continuum slope (beta), we estimate their stellar metallicities Z* from stellar photospheric absorption features at 1501 and 1719 Angstrom, which are calibrated with Starburst99 models and are largely unaffected by stellar age, dust, IMF, nebular continuum or interstellar absorption. Comparing them to photometric based spectral slopes in the range 1250-1750 Angstrom, we find that the stellar metallicity increases by ~0.5 dex from beta ~ -2 to beta ~ -1 (1 < A(1600) < 3.2), and a dependence with beta holds at fixed UV absolute luminosity M(UV) and stellar mass up to 10^(9.65) Msun. As a result, the metallicity is a fundamental ingredient for properly rescaling dust corrections based on M(UV) and M*. Using the same absorption features, we analyze the mass-metallicity relation (MZR) and find it is consistent with the previous VANDELS estimation based on a global fit of the FUV spectra. Similarly, we do not find a significant evolution between z=2 and z=3.5. Finally, the slopes of our MZR and Z*-beta relation are in agreement with the predictions of well-studied semi-analytic models of galaxy formation (SAM), while some tensions with observations remain as to the absolute metallicity normalization. The relation between UV slope and stellar metallicity is fundamental for the exploitation of large volume surveys with next-generation telescopes and for the physical characterization of galaxies in the first billion years of our Universe. △ Less

Submitted 12 November, 2020; originally announced November 2020.

Comments: 25 pages, 4 tables, 16 figures, and 5 Appendix subsections; accepted for publication in A&A

Journal ref: A&A 646, A39 (2021)

The First Robust Constraints on the Relationship Between Dust-to-Gas Ratio and Metallicity in Luminous Star-forming Galaxies at High Redshift

Authors: Alice E. Shapley, Fergus Cullen, James S. Dunlop, Ross J. McLure, Mariska Kriek, Naveen A. Reddy, Ryan L. Sanders

Abstract: We present rest-optical spectroscopic properties of a sample of four galaxies in the Atacama Large Millimeter/submillimeter Array Hubble Ultra Deep Field (ALMA HUDF). These galaxies span the redshift range $1.41 \leq z \leq 2.54$ and the stellar mass range $10.36\leq\log(M_*/{\rm M}_{\odot})\leq10.91$. They have existing far-infrared and radio measurements of dust-continuum and molecular gas emiss… ▽ More We present rest-optical spectroscopic properties of a sample of four galaxies in the Atacama Large Millimeter/submillimeter Array Hubble Ultra Deep Field (ALMA HUDF). These galaxies span the redshift range $1.41 \leq z \leq 2.54$ and the stellar mass range $10.36\leq\log(M_*/{\rm M}_{\odot})\leq10.91$. They have existing far-infrared and radio measurements of dust-continuum and molecular gas emission from which bolometric star-formation rates (SFRs), dust masses, and molecular gas masses have been estimated. We use new $H$- and $K$-band near-infrared spectra from the Keck/MOSFIRE spectrograph to estimate SFRs from dust-corrected H$α$ emission (SFR(H$α$)) and gas-phase oxygen abundances from the ratio [NII]$λ6584$/H$α$. We find that the dust-corrected SFR(H$α$) is systematically lower than the bolometric SFR by a factor of several, and measure gas-phase oxygen abundances in a narrow range, $12+\log(\mbox{O/H})=8.59-8.69$ ($0.8-1.0\: (\mbox{O/H})_{\odot}$). Relative to a large $z\sim 2$ comparison sample from the MOSDEF survey, the ALMA HUDF galaxies scatter roughly symmetrically around the best-fit linear mass-metallicity relation, providing tentative evidence for a flattening in the SFR dependence of metallicity at high stellar mass. Combining oxygen abundances with estimates of dust and molecular gas masses, we show that there is no significant evolution in the normalization of the dust-to-gas ratio DGR vs. metallicity relation from $z\sim0$ to $z\sim2$. This result is consistent with some semi-analytic models and cosmological simulations describing the evolution of dust in galaxies. Tracing the actual form of the DGR vs. metallicity relation at high redshift now requires combined measurements of dust, gas, and metallicity over a significantly wider range in metallicity. △ Less

Submitted 28 October, 2020; v1 submitted 21 September, 2020; originally announced September 2020.

Comments: 6 pages, 3 figures, accepted to ApJ Letters

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

Quasar Sightline and Galaxy Evolution (QSAGE) survey -- II. Galaxy overdensities around UV luminous quasars at z=1-2

Authors: J. P. Stott, R. M. Bielby, F. Cullen, J. N. Burchett, N. Tejos, M. Fumagalli, R. A. Crain, S. L. Morris, N. Amos, R. G. Bower, J. X. Prochaska

Abstract: We demonstrate that the UV brightest quasars at z=1-2 live in overdense environments. This is based on an analysis of deep Hubble Space Telescope WFC3 G141 grism spectroscopy of the galaxies along the lines-of-sight to UV luminous quasars in the redshift range z=1-2. This constitutes some of the deepest grism spectroscopy performed by WFC3, with 4 roll angles spread over a year of observations to… ▽ More We demonstrate that the UV brightest quasars at z=1-2 live in overdense environments. This is based on an analysis of deep Hubble Space Telescope WFC3 G141 grism spectroscopy of the galaxies along the lines-of-sight to UV luminous quasars in the redshift range z=1-2. This constitutes some of the deepest grism spectroscopy performed by WFC3, with 4 roll angles spread over a year of observations to mitigate the effect of overlapping spectra. Of the 12 quasar fields studied, 8 display evidence for a galaxy overdensity at the redshift of the quasar. One of the overdensities, PG0117+213 at z=1.50, has potentially 36 spectroscopically confirmed members, consisting of 19 with secure redshifts and 17 with single-line redshifts, within a cylinder of radius ~700 kpc. Its halo mass is estimated to be log (M/Msol)=14.7. This demonstrates that spectroscopic and narrow-band observations around distant UV bright quasars may be an excellent route for discovering protoclusters. Our findings agree with previous hints from statistical observations of the quasar population and theoretical works, as feedback regulated black hole growth predicts a correlation between quasar luminosity and halo mass. We also present the high signal-to-noise rest-frame optical spectral and photometric properties of the quasars themselves. △ Less

Submitted 19 August, 2020; v1 submitted 12 June, 2020; originally announced June 2020.

Comments: Published in MNRAS. This is the author accepted version. A small number of typos were corrected after journal proof stage

Journal ref: 2020, MNRAS, 497, 3083