Showing 1–50 of 197 results for author: McLure, R J

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

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

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

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

Comments: 17 pages, 8 figures, submitted to ApJ

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

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

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

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

Submitted 13 March, 2024; originally announced March 2024.

Comments: 29 pages, 17 figures, submitted to MNRAS

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

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

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

Submitted 11 March, 2024; originally announced March 2024.

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

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

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

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

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

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

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

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

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

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

Comments: 18 pages, 12 figures, MNRAS accepted

The 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

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

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

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

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

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

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

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

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

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

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

A new measurement of the expansion history of the Universe at z=1.26 with cosmic chronometers in VANDELS

Authors: E. Tomasetti, M. Moresco, N. Borghi, K. Jiao, A. Cimatti, L. Pozzetti, A. C. Carnall, R. J. McLure, L. Pentericci

Abstract: We derive a new constraint on the expansion history of the Universe by applying the cosmic chronometers method, studying the age evolution of high-redshift galaxies with a full-spectral-fitting approach. We select a sample of 39 massive ($log(M/M_\odot)>10.8$) and passive ($log(sSFR/yr^{-1})<-11$) galaxies from the data release 4 of the VANDELS survey at $1<z<1.5$, combining different selection cr… ▽ More We derive a new constraint on the expansion history of the Universe by applying the cosmic chronometers method, studying the age evolution of high-redshift galaxies with a full-spectral-fitting approach. We select a sample of 39 massive ($log(M/M_\odot)>10.8$) and passive ($log(sSFR/yr^{-1})<-11$) galaxies from the data release 4 of the VANDELS survey at $1<z<1.5$, combining different selection criteria to minimize the potential contamination by star-forming outliers. We perform full-spectral-fitting jointly on spectra and photometry of our sources with the code BAGPIPES, without any cosmological assumption on the age of the population. The derived physical properties of the selected galaxies are characteristic of a passive population, with short star formation timescales ($<τ>=0.28\pm0.02$ Gyr), low dust extinction ($<A_{V,dust}>=0.43\pm0.02$ mag), and sub-solar metallicities ($<Z/Z_{\odot}>=0.44\pm0.01$). The ages show a decreasing trend with redshift compatible with a standard cosmological model, even if no cosmological constraint is assumed in the fit, and a clear mass-downsizing pattern. Testing the impact of the star formation history on the results, we find only a maximum 2\% fluctuation in age and metallicity. By fitting the median age-redshift relation with a flat $Λ$CDM model and assuming a Gaussian prior on $Ω_{M,0}= 0.3\pm0.02$ from late-Universe probes, we obtain $H_0=67_{-15}^{+14}\:km\:s^{-1}\:Mpc^{-1}$. In the end, we derive a new estimate of the Hubble parameter with the cosmic chronometers method, $H(z=1.26)=135\pm65\:km\:s^{-1}\:Mpc^{-1}$ including statistical and systematic errors. While the error budget is currently dominated by the scarcity of the sample, this work proves the potential strength of the cosmic chronometers approach up to $z>1$, especially in view of incoming large spectroscopic surveys like Euclid. (abridged) △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: 18 pages, 12 figures, 4 tables. Submitted to A&A

Journal ref: A&A 679, A96 (2023)

The ionizing photon production efficiency of bright z$\sim$2-5 galaxies

Authors: M. Castellano, D. Belfiori, L. Pentericci, A. Calabrò, S. Mascia, L. Napolitano, F. Caro, S. Charlot, J. Chevallard, E. Curtis-Lake, M. Talia, A. Bongiorno, A. Fontana, J. P. U. Fynbo, B. Garilli, L. Guaita, R. J. McLure, E. Merlin, M. Mignoli, M. Moresco, E. Pompei, L. Pozzetti, A. Saldana Lopez, A. Saxena, P. Santini , et al. (5 additional authors not shown)

Abstract: We investigate the production efficiency of ionizing photons ($ξ_{ion}^*$) of 1174 galaxies with secure redshift at z=2-5 from the VANDELS survey to determine the relation between ionizing emission and physical properties of bright and massive sources. We constrain $ξ_{ion}^*$ and galaxy physical parameters by means of spectro-photometric fits performed with the BEAGLE code. The analysis exploits… ▽ More We investigate the production efficiency of ionizing photons ($ξ_{ion}^*$) of 1174 galaxies with secure redshift at z=2-5 from the VANDELS survey to determine the relation between ionizing emission and physical properties of bright and massive sources. We constrain $ξ_{ion}^*$ and galaxy physical parameters by means of spectro-photometric fits performed with the BEAGLE code. The analysis exploits the multi-band photometry in the VANDELS fields, and the measurement of UV rest-frame emission lines (CIII]$λ1909$, HeII$λ1640$, OIII]$λ1666$) from deep VIMOS spectra. We find no clear evolution of $ξ_{ion}^*$ with redshift within the probed range. The ionizing efficiency slightly increases at fainter $M_{UV}$, and bluer UV slopes, but these trends are less evident when restricting the analysis to a complete subsample at log(M$_{star}$/M$_{\odot}$)$>$9.5. We find a significant trend of increasing $ξ_{ion}^*$ with increasing EW(Ly$α$), with an average log($ξ_{ion}^*$/Hz erg$^{-1}$)$>$25 at EW$>$50Å, and a higher ionizing efficiency for high-EW CIII]$λ1909$ and OIII]$λ1666$ emitters. The most significant correlations are found with respect to stellar mass, specific star-formation rate (sSFR) and SFR surface density ($Σ_{SFR}$). The relation between $ξ_{ion}^*$ and sSFR shows a monotonic increase from log($ξ_{ion}^*$/Hz erg$^{-1}$) $\sim$24.5 at log(sSFR)$\sim$-9.5$yr^{-1}$ to $\sim$25.5 at log(sSFR)$\sim$-7.5$yr^{-1}$, a low scatter and little dependence on mass. The objects above the main-sequence of star-formation consistently have higher-than-average $ξ_{ion}^*$. A clear increase of $ξ_{ion}^*$ with $Σ_{SFR}$ is also found, with log($ξ_{ion}^*$/Hz erg$^{-1}$)$>$25 for objects at $Σ_{SFR}>$10 M$_{\odot}/yr/kpc^2$.(Abridged) △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: 11 pages, 9 figures, accepted for publication in A&A

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 ultraviolet continuum slopes ($\mathbfβ$) of galaxies at $\mathbf{z\simeq8-16}$ from JWST and ground-based near-infrared imaging

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

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

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

Comments: 9 pages, 5 figures, MNRAS accepted version

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Comments: 14 pages, 9 figures, submitted to MNRAS

The 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 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

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 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 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 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

MOONRISE: The Main MOONS GTO Extragalactic Survey

Authors: R. Maiolino, M. Cirasuolo, J. Afonso, F. E. Bauer, R. Bowler, O. Cucciati, E. Daddi, G. De Lucia, C. Evans, H. Flores, A. Gargiulo, B. Garilli, P. Jablonka, M. Jarvis, J. -P. Kneib, S. Lilly, T. Looser, M. Magliocchetti, Z. Man, F. Mannucci, S. Maurogordato, R. J. McLure, P. Norberg, P. Oesch, E. Oliva , et al. (12 additional authors not shown)

Abstract: The MOONS instrument possesses an exceptional combination of large multiplexing, high sensitivity, broad simultaneous spectral coverage (from optical to near-infrared bands), large patrol area and high fibre density. These properties provide the unprecedented potential of enabling, for the very first time, SDSS-like surveys around Cosmic Noon (z~1-2.5), when the star formation rate in the Universe… ▽ More The MOONS instrument possesses an exceptional combination of large multiplexing, high sensitivity, broad simultaneous spectral coverage (from optical to near-infrared bands), large patrol area and high fibre density. These properties provide the unprecedented potential of enabling, for the very first time, SDSS-like surveys around Cosmic Noon (z~1-2.5), when the star formation rate in the Universe peaked. The high-quality spectra delivered by MOONS will sample the same nebular and stellar diagnostics observed in extensive surveys of local galaxies, providing an accurate and consistent description of the evolution of various physical properties of galaxies, and hence a solid test of different scenarios of galaxy formation and transformation. Most importantly, by spectroscopically identifying hundreds of thousands of galaxies at high redshift, the MOONS surveys will be capable of determining the environments in which primeval galaxies lived and will reveal how such environments affected galaxy evolution. In this article, we specifically focus on the main Guaranteed Time Observation (GTO) MOONS extragalactic survey, MOONRISE, by providing an overview of its scientific goals and observing strategy. △ Less

Submitted 1 September, 2020; originally announced September 2020.

Comments: 6 pages, 6 figures, Published on the ESO Messenger No. 180

Journal ref: The Messenger, vol. 180, 2020, p. 24-29

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

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

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

Submitted 7 February, 2020; originally announced February 2020.

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

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

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

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

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

Comments: 14 pages, 8 figures, MNRAS accepted

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

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

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

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

Comments: 10 pages, 6 figures, MNRAS accepted

The rest-frame UV luminosity function at $z \simeq 4$: a significant contribution of AGN to the bright-end of the galaxy population

Authors: N. J. Adams, R. A. A. Bowler, M. J. Jarvis, B. Häußler, R. J. McLure, A. Bunker, J. S. Dunlop, A. Verma

Abstract: We measure the rest-frame UV luminosity function (LF) at $z \sim 4$ self-consistently over a wide range in absolute magnitude ($-27 \lesssim M_{\rm UV} \lesssim -20$). The LF is measured with 46,904 sources selected using a photometric redshift approach over $ \sim 6$ deg$^2$ of the combined COSMOS and XMM-LSS fields. We simultaneously fit for both AGN and galaxy LFs using a combination of Schecht… ▽ More We measure the rest-frame UV luminosity function (LF) at $z \sim 4$ self-consistently over a wide range in absolute magnitude ($-27 \lesssim M_{\rm UV} \lesssim -20$). The LF is measured with 46,904 sources selected using a photometric redshift approach over $ \sim 6$ deg$^2$ of the combined COSMOS and XMM-LSS fields. We simultaneously fit for both AGN and galaxy LFs using a combination of Schechter or Double Power Law (DPL) functions alongside a single power law for the faint-end slope of the AGN LF. We find a lack of evolution in the shape of the bright-end of the LBG component when compared to other studies at $z \simeq 5$ and evolutionary recipes for the UV LF. Regardless of whether the LBG LF is fit with a Schechter function or DPL, AGN are found to dominate at $M_{\rm UV} < -23.5$. We measure a steep faint-end slope of the AGN LF with $α_{AGN} = -2.09^{+0.35}_{-0.38}$ ($-1.66^{+0.29}_{-0.58}$) when fit alongside a Schechter function (DPL) for the galaxies. Our results suggest that if AGN are morphologically selected it results in a bias to lower number densities. Only by considering the full galaxy population over the transition region from AGN to LBG domination can an accurate measurement of the total LF be attained. △ Less

Submitted 3 March, 2020; v1 submitted 3 December, 2019; originally announced December 2019.

Comments: 14 pages, 9 figures, 3 tables, Updated to match MNRAS accepted version (accepted 3/3/2020) after minor revisions

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

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

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

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

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

The Intergalactic medium transmission towards z>4 galaxies with VANDELS and the impact of dust attenuation

Authors: R. Thomas, L. Pentericci, O. Le Fèvre, G. Zamorani, D. Schaerer, R. Amorin, M. Castellano, A. C. Carnall, S. Cristiani, F. Cullen, S. L. Finkelstein, F. Fontanot, L. Guaita, P. Hibon, N. Hathi, J. P. U. Fynbo, Y. Khusanova, A. M. Koekemoer, D. McLeod, R. J. McLure, F. Marchi, L. Pozzetti, A. Saxena, M. Talia, M. Bolzonella

Abstract: Aims. Our aim is to estimate the intergalactic medium transmission towards UV-selected star-forming galaxies at redshift 4 and above and study the effect of the dust attenuation on these measurements. Methods. The ultra-violet spectrum of high redshift galaxies is a combination of their intrinsic emission and the effect of the Inter-Galactic medium (IGM) absorption along their line of sight. Usi… ▽ More Aims. Our aim is to estimate the intergalactic medium transmission towards UV-selected star-forming galaxies at redshift 4 and above and study the effect of the dust attenuation on these measurements. Methods. The ultra-violet spectrum of high redshift galaxies is a combination of their intrinsic emission and the effect of the Inter-Galactic medium (IGM) absorption along their line of sight. Using data coming from the unprecedented deep spectroscopy from the VANDELS ESO public survey carried out with the VIMOS instrument we compute both the dust extinction and the mean transmission of the IGM as well as its scatter from a set of 281 galaxies at z>3.87. Because of a degeneracy between the dust content of the galaxy and the IGM, we first estimate the stellar dust extinction parameter E(B-V) and study the result as a function of the dust prescription. Using these measurements as constraint for the spectral fit we estimate the IGM transmission Tr(Lyalpha). Both photometric and spectroscopic SED fitting are done using the SPectroscopy And photometRy fiTting tool for Astronomical aNalysis (SPARTAN) that is able to fit the spectral continuum of the galaxies as well as photometric data. Results. Using the classical Calzetti's attenuation law we find that E(B-V) goes from 0.11 at z=3.99 to 0.08 at z=5.15. These results are in very good agreement with previous measurements from the literature. We estimate the IGM transmission and find that the transmission is decreasing with increasing redshift from Tr(Lyalpha)=0.53 at z=3.99 to 0.28 at z=5.15. We also find a large standard deviation around the average transmission that is more than 0.1 at every redshift. Our results are in very good agreement with both previous measurements from AGN studies and with theoretical models. △ Less

Submitted 28 November, 2019; originally announced November 2019.

Comments: Accepted for publication in A&A

Journal ref: A&A 634, A110 (2020)

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

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

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

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

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

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

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

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

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

Submitted 7 August, 2019; originally announced August 2019.

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

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

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

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

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

Comments: 25 pages, 13 figures, MNRAS accepted

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

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

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

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

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

The VANDELS survey: the role of ISM and galaxy physical properties on the escape of Lyα emission in z $\sim$ 3.5 star-forming galaxies

Authors: F. Marchi, L. Pentericci, L. Guaita, M. Talia, M. Castellano, N. Hathi, D. Schaerer, R. Amorin, A. C. Carnall, S. Charlot, J. Chevallard, F. Cullen, S. L. Finkelstein, A. Fontana, F. Fontanot, B. Garilli, P. Hibon, A. M. Koekemoer, D. Maccagni, R. J. McLure, C. Papovich, L. Pozzetti, A. Saxena

Abstract: We investigate the physical properties of a sample of 52 Lyα emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in the escape and spatial distribution of Lyα photons. For these galaxies, we derive different physical properties and characterize the Lyα emission in terms of kinematics, EW, FWHM and spatial extension and then e… ▽ More We investigate the physical properties of a sample of 52 Lyα emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in the escape and spatial distribution of Lyα photons. For these galaxies, we derive different physical properties and characterize the Lyα emission in terms of kinematics, EW, FWHM and spatial extension and then estimate the velocity of the neutral outflowing gas. We reproduce some of the well known correlations between Lyα EW and stellar mass, dust extinction and UV \b{eta} slope, in the sense that the emission line appears brighter in lower mass, less dusty and bluer galaxies. We do not find any correlation with the SED-derived star formation rate, while we find that galaxies with brighter Lyα tend to be more compact both in UV and in Lyα. Our data reveal a new interesting correlation between the Lyα velocity and the offset of the inter-stellar absorption lines with respect to the systemic redshift, in the sense that galaxies with larger inter-stellar medium (ISM) out-flow velocities show smaller Lyα velocity shifts. We interpret this relation in the context of the shell-model scenario, where the velocity of the ISM and the HI column density contribute together in determining the Lyα kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Lyα spatial profiles, a sign of increased scattering. However, we do not find any evidence that the HI column density is related to any other physical properties of the galaxies, although this might be due in part to the limited range of parameters that our sample spans. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: submitted to A&A

Journal ref: A&A 631, A19 (2019)

Maximising the power of deep extragalactic imaging surveys with the James Webb Space Telescope

Authors: T. W. Kemp, J. S. Dunlop, R. J. McLure, C. Schreiber, A. C. Carnall, F. Cullen

Abstract: We present a new analysis of the potential power of deep, near-infrared, imaging surveys with the James Webb Space Telescope (JWST) to improve our knowledge of galaxy evolution. In this work we properly simulate what can be achieved with realistic survey strategies, and utilise rigorous signal:noise calculations to calculate the resulting posterior constraints on the physical properties of galaxie… ▽ More We present a new analysis of the potential power of deep, near-infrared, imaging surveys with the James Webb Space Telescope (JWST) to improve our knowledge of galaxy evolution. In this work we properly simulate what can be achieved with realistic survey strategies, and utilise rigorous signal:noise calculations to calculate the resulting posterior constraints on the physical properties of galaxies. We explore a broad range of assumed input galaxy types (>20,000 models, including extremely dusty objects) across a wide redshift range (out to z~12), while at the same time considering a realistic mix of galaxy properties based on our current knowledge of the evolving population (as quantified through the Empirical Galaxy Generator: EGG). While our main focus is on imaging surveys with NIRCam, spanning lambda(obs) = 0.6-5.0 microns, an important goal of this work is to quantify the impact/added-value of: i) parallel imaging observations with MIRI at longer wavelengths, and ii) deeper supporting optical/UV imaging with HST (potentially prior to JWST launch) in maximising the power and robustness of a major extragalactic NIRCam survey. We show that MIRI parallel 7.7-micron imaging is of most value for better constraining the redshifts and stellar masses of the dustiest (A_V > 3) galaxies, while deep B-band imaging (reaching~28.5 AB mag) with ACS on HST is vital for determining the redshifts of the large numbers of faint/low-mass, z < 5 galaxies that will be detected in a deep JWST NIRCam survey. △ Less

Submitted 19 March, 2019; originally announced March 2019.

Comments: 19 Pages, 11 Figures, Submitted to MNRAS

An Ultra Deep Field survey with WFIRST

Authors: Anton M. Koekemoer, R. J. Foley, D. N. Spergel, M. Bagley, R. Bezanson, F. B. Bianco, R. Bouwens, L. Bradley, G. Brammer, P. Capak, I. Davidzon, G. De Rosa, M. E. Dickinson, O. Doré, J. S. Dunlop, R. S. Ellis, X. Fan, G. G. Fazio, H. C. Ferguson, A. V. Filippenko, S. Finkelstein, B. Frye, E. Gawiser, N. A. Grogin, N. P. Hathi , et al. (47 additional authors not shown)

Abstract: Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their siz… ▽ More Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.) △ Less

Submitted 19 March, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

VANDELS consortium public database

Authors: L. Pentericci, R. J. McLure, P. Franzetti, B. Garilli, the VANDELS team

Abstract: We provide access to the current VANDELS data release (DR2) through a dedicated database, accessible at the link http://vandels.inaf.it/db. The data release includes the spectra for all galaxies for which the scheduled integration time was completed during the first two observing seasons and the spectra for those galaxies for which the scheduled integration time was 50\% complete at the end of sea… ▽ More We provide access to the current VANDELS data release (DR2) through a dedicated database, accessible at the link http://vandels.inaf.it/db. The data release includes the spectra for all galaxies for which the scheduled integration time was completed during the first two observing seasons and the spectra for those galaxies for which the scheduled integration time was 50\% complete at the end of season two (i.e. 20/40 hours and 40/80 hours). The total number of spectra released is 1362 (586 in CDFS and 776 in UDS). Spectra are stored as multi-extension FITS files containing the 1D extracted spectrum, the 2D linearly re-sampled spectrum, the 1D sky spectrum, the 1D noise estimate and the image thumbnail of the object. △ Less

Submitted 13 November, 2018; originally announced November 2018.