
==== Front
Nature
Nature
Nature
0028-0836
1476-4687
Nature Publishing Group UK London

38447669
7227
10.1038/s41586-024-07227-0
Article
A recently quenched galaxy 700 million years after the Big Bang
http://orcid.org/0000-0002-3642-2446
Looser Tobias J. tjl54@cam.ac.uk

12
D’Eugenio Francesco 12
http://orcid.org/0000-0002-4985-3819
Maiolino Roberto 123
http://orcid.org/0000-0002-7595-121X
Witstok Joris 12
http://orcid.org/0000-0001-9276-7062
Sandles Lester 12
http://orcid.org/0000-0002-9551-0534
Curtis-Lake Emma 4
Chevallard Jacopo 5
http://orcid.org/0000-0002-8224-4505
Tacchella Sandro 12
http://orcid.org/0000-0002-9280-7594
Johnson Benjamin D. 6
http://orcid.org/0000-0003-0215-1104
Baker William M. 12
http://orcid.org/0000-0002-1714-1905
Suess Katherine A. 789
http://orcid.org/0000-0002-6719-380X
Carniani Stefano 10
http://orcid.org/0000-0001-8895-0606
Ferruit Pierre 11
http://orcid.org/0000-0001-7997-1640
Arribas Santiago 12
Bonaventura Nina 1314
http://orcid.org/0000-0002-8651-9879
Bunker Andrew J. 5
http://orcid.org/0000-0002-0450-7306
Cameron Alex J. 5
http://orcid.org/0000-0003-3458-2275
Charlot Stephane 15
http://orcid.org/0000-0002-2678-2560
Curti Mirko 1216
http://orcid.org/0000-0002-2380-9801
de Graaff Anna 17
http://orcid.org/0000-0003-0695-4414
Maseda Michael V. 18
http://orcid.org/0000-0002-7028-5588
Rawle Tim 19
http://orcid.org/0000-0003-4996-9069
Rix Hans-Walter 17
Del Pino Bruno Rodríguez 12
http://orcid.org/0000-0001-8034-7802
Smit Renske 20
http://orcid.org/0000-0003-4891-0794
Übler Hannah 12
http://orcid.org/0000-0002-4201-7367
Willott Chris 21
Alberts Stacey 22
http://orcid.org/0000-0003-1344-9475
Egami Eiichi 22
http://orcid.org/0000-0002-2929-3121
Eisenstein Daniel J. 6
Endsley Ryan 23
http://orcid.org/0000-0002-8543-761X
Hausen Ryan 24
http://orcid.org/0000-0002-7893-6170
Rieke Marcia 22
http://orcid.org/0000-0002-4271-0364
Robertson Brant 7
Shivaei Irene 22
http://orcid.org/0000-0003-2919-7495
Williams Christina C. 25
http://orcid.org/0000-0003-4109-304X
Boyett Kristan 26
http://orcid.org/0000-0002-2178-5471
Chen Zuyi 22
Ji Zhiyuan 22
http://orcid.org/0000-0002-0267-9024
Jones Gareth C. 5
http://orcid.org/0000-0002-5320-2568
Kumari Nimisha 27
http://orcid.org/0000-0002-7524-374X
Nelson Erica 28
Perna Michele 12
http://orcid.org/0000-0001-5333-9970
Saxena Aayush 35
Scholtz Jan 12
1 https://ror.org/013meh722 grid.5335.0 0000 0001 2188 5934 Kavli Institute for Cosmology, University of Cambridge, Cambridge, UK
2 https://ror.org/013meh722 grid.5335.0 0000 0001 2188 5934 Cavendish Laboratory - Astrophysics Group, University of Cambridge, Cambridge, UK
3 https://ror.org/02jx3x895 grid.83440.3b 0000 0001 2190 1201 Department of Physics and Astronomy, University College London, London, UK
4 https://ror.org/0267vjk41 grid.5846.f 0000 0001 2161 9644 Centre for Astrophysics Research, Department of Physics, Astronomy and Mathematics, University of Hertfordshire, Hatfield, UK
5 https://ror.org/052gg0110 grid.4991.5 0000 0004 1936 8948 Department of Physics, University of Oxford, Oxford, UK
6 https://ror.org/03c3r2d17 grid.455754.2 Center for Astrophysics | Harvard & Smithsonian, Cambridge, MA USA
7 grid.205975.c 0000 0001 0740 6917 Department of Astronomy and Astrophysics, University of California, Santa Cruz, Santa Cruz, CA USA
8 grid.168010.e 0000000419368956 Kavli Institute for Particle Astrophysics and Cosmology (KIPAC), Stanford University, Stanford, CA USA
9 https://ror.org/00f54p054 grid.168010.e 0000 0004 1936 8956 Department of Physics, Stanford University, Stanford, CA USA
10 https://ror.org/03aydme10 grid.6093.c Scuola Normale Superiore, Pisa, Italy
11 https://ror.org/00kw1sm04 grid.450273.7 0000 0004 0623 7009 European Space Astronomy Centre (ESAC), European Space Agency (ESA), Madrid, Spain
12 grid.4711.3 0000 0001 2183 4846 Centro de Astrobiología (CAB), Spanish National Research Council (CSIC)–National Institute of Aerospace Technology (INTA), Madrid, Spain
13 grid.5254.6 0000 0001 0674 042X Cosmic Dawn Center (DAWN), Copenhagen, Denmark
14 https://ror.org/035b05819 grid.5254.6 0000 0001 0674 042X Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark
15 grid.435813.8 0000 0001 0540 8249 Sorbonne Université, CNRS, UMR 7095, Institut d’Astrophysique de Paris, Paris, France
16 https://ror.org/01qtasp15 grid.424907.c 0000 0004 0645 6631 European Southern Observatory, Garching bei Muenchen, Germany
17 https://ror.org/01vhnrs90 grid.429508.2 0000 0004 0491 677X Max-Planck-Institut für Astronomie, Heidelberg, Germany
18 https://ror.org/01y2jtd41 grid.14003.36 0000 0001 2167 3675 Department of Astronomy, University of Wisconsin-Madison, Madison, WI USA
19 https://ror.org/036f5mx38 grid.419446.a 0000 0004 0591 6464 European Space Agency (ESA) Office, Space Telescope Science Institute (STScI), Baltimore, MD USA
20 https://ror.org/04zfme737 grid.4425.7 0000 0004 0368 0654 Astrophysics Research Institute, Liverpool John Moores University, Liverpool, UK
21 grid.469915.6 0000 0001 1945 2224 NRC Herzberg, Victoria, British Columbia Canada
22 https://ror.org/03m2x1q45 grid.134563.6 0000 0001 2168 186X Steward Observatory, University of Arizona, Tucson, AZ USA
23 https://ror.org/00hj54h04 grid.89336.37 0000 0004 1936 9924 Department of Astronomy, University of Texas at Austin, Austin, TX USA
24 https://ror.org/00za53h95 grid.21107.35 0000 0001 2171 9311 Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD USA
25 https://ror.org/03zmsge54 grid.510764.1 NSF’s National Optical-Infrared Astronomy Research Laboratory (NOIRLab), Tucson, AZ USA
26 https://ror.org/01ej9dk98 grid.1008.9 0000 0001 2179 088X School of Physics, University of Melbourne, Parkville, Victoria Australia
27 https://ror.org/036f5mx38 grid.419446.a 0000 0004 0591 6464 AURA for European Space Agency, Space Telescope Science Institute, Baltimore, MD USA
28 https://ror.org/02ttsq026 grid.266190.a 0000 0000 9621 4564 Department of Astrophysical and Planetary Sciences, University of Colorado Boulder, Boulder, CO USA
6 3 2024
6 3 2024
2024
629 8010 5357
28 2 2023
22 2 2024
© The Author(s) 2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Local and low-redshift (z < 3) galaxies are known to broadly follow a bimodal distribution: actively star-forming galaxies with relatively stable star-formation rates and passive systems. These two populations are connected by galaxies in relatively slow transition. By contrast, theory predicts that star formation was stochastic at early cosmic times and in low-mass systems1–4. These galaxies transitioned rapidly between starburst episodes and phases of suppressed star formation, potentially even causing temporary quiescence—so-called mini-quenching events5,6. However, the regime of star-formation burstiness is observationally highly unconstrained. Directly observing mini-quenched galaxies in the primordial Universe is therefore of utmost importance to constrain models of galaxy formation and transformation7,8. Early quenched galaxies have been identified out to redshift z < 5 (refs. 9–12) and these are all found to be massive (M⋆ > 1010 M⊙) and relatively old. Here we report a (mini-)quenched galaxy at z = 7.3, when the Universe was only 700 Myr old. The JWST/NIRSpec spectrum is very blue (U–V = 0.16 ± 0.03 mag) but exhibits a Balmer break and no nebular emission lines. The galaxy experienced a short starburst followed by rapid quenching; its stellar mass (4–6 × 108 M⊙) falls in a range that is sensitive to various feedback mechanisms, which can result in perhaps only temporary quenching.

Analysis of the JWST/NIRSpec spectrum of the recently observed Lyman-break galaxy JADES-GS+53.15508-27.80178 revealed a redshift of z = 7.3, a Balmer break and a complete absence of nebular emission lines, indicating that quenching occurred only 700 million years after the Big Bang.

Subject terms

Early universe
Galaxies and clusters
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcMain

The galaxy was first described as a Lyman-break galaxy13 and was recently observed as part of our JWST Advanced Deep Extragalactic Survey (JADES; galaxy ID: JADES-GS+53.15508-27.80178; hereafter simply JADES-GS-z7-01-QU) through deep (28-h) NIRSpec-MSA observations with the prism. The galaxy was pre-selected with the photometric Lyman dropout technique and a blue rest-frame ultraviolet (UV) colour.

The spectrum of JADES-GS-z7-01-QU is shown in Fig. 1. The redshift z = 7.29 ± 0.01 is unambiguously determined (using the BEAGLE code; see Methods) from the combined observed wavelengths of the characteristic Lyα drop and Balmer break.Fig. 1 NIRSpec R100/prism spectrum of JADES-GS-z7-01-QU.

The absence of emission lines, together with the Balmer break, reveals that this is a—temporarily or permanently—(mini-)quenched, post-starburst galaxy. The clearly detected Lyα drop and the Balmer break unambiguously give a redshift of z = 7.3. The vertical dashed lines indicate the rest-frame wavelengths of the strongest nebular emission lines. The red line indicates the pPXF spectral fit. The upper panel shows the signal-to-noise ratio (S/N) in the 2D prism spectrum. The bottom panel shows the ratio between the residuals of the fit and the noise. For reference, the flux in the F200W NIRCam filter is 3.33 ± 0.08 × 10−17 erg cm−2 s−1, fully consistent with the spectrum.

The 3σ upper limit on the Hβ emission-line flux, F(Hβ) < 6.1 × 10−20 erg cm−2 s−1, implies an upper limit on the star-formation rate (SFR) of <0.65 M⊙ yr−1 over the past 3–10 Myr (even accounting for dust attenuation; see Methods). Even stronger constraints come from the [O III]λ5008 line: we find F([O III]λ5008) < 6.5 × 10−20 erg s−1 cm−2, which—combined with a conservative assumption on the [O III]λ5008/Hβ ratios in high-z galaxies14,15—implies a 3σ limit on the SFR five times lower than the Hβ-derived value. The absence of emission lines is independently confirmed by the medium-resolution spectrum (see Methods).

We measure a UV slope β = −2.09 ± 0.09, typical for galaxies at 6 < z < 10 (refs. 16,17), indicating strong star-formation activity during the past 100 Myr before observation. In the rest-frame visible, we detect a clear Balmer break and Hδ absorption with equivalent width EWHδA=4.8±1.0Å. This value, combined with the absence of emission lines, means that JADES-GS-z7-01-QU meets the most common spectroscopic definition of a post-starburst galaxy18,19, that is, a galaxy that has only recently stopped forming stars.

Previous high-redshift works have identified Balmer-break galaxies in the epoch of reionization20–22, indicating the existence of evolved stellar populations and even proposing quiescent phases in these objects20,21. However, without spectroscopy, one cannot rule out the presence of emission lines with low equivalent width or that strong emission lines masquerade as Balmer breaks. Furthermore, because of the lack of atmospheric transmission at wavelengths longer than 2.5 microns, it is impossible to investigate Balmer breaks at z > 5 from the ground. Therefore, before JWST, it was impossible to confirm the absence of continuing star formation.

Crucially, on the basis of colours alone, this (mini-)quenched galaxy would have been identified as ‘star forming’ by the colour selection criteria23, even if including the extension to fast-quenched galaxies24. Indeed, its rest-frame U–V colour of 0.16 ± 0.03 mag places it outside the quiescent region of the UVJ diagram, regardless of V–J colour25, as it is the case for other quiescent galaxies at high redshift7. However, thanks to JWST/NIRSpec, we can place stringent upper limits on the nebular emission-line fluxes.

Are there potential alternatives to the quenched interpretation? A very high escape fraction of ionizing Lyman-continuum (LyC) photons with fesc > 0.9 could strongly suppress nebular emission26. However, if fesc is high, this would be because nearly all of the interstellar medium (ISM) was ejected or consumed by star formation27; yet, if the ISM is absent, there is no fuel for star formation and the galaxy must be quenched. This makes the galaxy highly interesting in the context of reionization, as a remnant leaker28. The question is whether the object is (still) a remnant leaker at the epoch of observation. In other words, whether there are still very young stellar populations (a few Myr old) that would still be producing ionizing photons associated with O-type stars and which would largely escape the galaxy, as fesc ≈ 1. This scenario is disfavoured by the normal UV slope β (ref. 29), the Balmer break and by the strong Hδ absorption.

Statistically, a very recently (<10 Myr) star-forming solution with high fesc is also disfavoured by our further analysis. Indeed, by making use of the flexibility of the software BEAGLE to model the observed spectrum, we find that a high-fesc, recently star-forming solution—although possible—is strongly disfavoured compared with the quenched (>3–10 Myr) solution (see Methods). Furthermore, as we will discuss below, both the pPXF and Prospector codes, which can optionally decouple the continuum from the nebular lines (which are degenerate with fesc), do not favour a solution with very recent star formation. The second alternative that we cannot completely rule out is the presence of completely obscured star formation, as advocated for some post-starburst galaxies in the local Universe30. However, we note that high dust masses and high dust extinction in such low-mass systems, at such high redshift, have never been observed31.

To estimate the physical properties of the galaxy including stellar mass M⋆, SFR, star-formation history (SFH), dust attenuation and stellar metallicity, we apply joint spectrophotometric modelling of its spectral energy distribution (SED). To marginalize over model assumptions and implementation, we use four different SED-fitting codes (pPXF, BAGPIPES, Prospector and BEAGLE; see Methods). Figure 1 shows, as an example, the best-fit pPXF model in red, overlaid on the spectrum.

The methods agree on a low stellar mass of M⋆ = 4–6 × 108 M⊙ (Table 1); in other words, this is an object in the dwarf-galaxy regime—essentially the same mass as the nearby, actively star-forming Small Magellanic Cloud, but at z = 7.3 and quenched.Table 1 Key physical quantities inferred by the four full spectral fitting codes pPXF, BAGPIPES, BEAGLE and Prospector

Key inferred properties	pPXF	BAGPIPES	BEAGLE	Prospector	
log10(M⋆/M⊙)	–	8.5 ± 0.1	8.8−0.2+0.1	8.7−0.1+0.1	
log10[SFR (M⊙ yr−1)]	–	<−1.0	−2.5−1.0+1.0	−2.6−2.7+1.5	
log10(Z/Z⊙)	<−2.0	−0.7 ± 0.1	−1.9−0.2+0.4	−1.7−0.2+0.2	
tquench (Myr)	About 50	18−5+5	16−4+7	38−10+9	
tform (Myr)	About 150	37−5+8	93−47+69	116−45+85	
AV (mag)	0.4 ± 0.1	0.32−0.23+0.26	0.51−0.04+0.03	0.1−0.0+0.1	
M⋆, stellar mass; SFR, star-formation rate; Z, metallicity; tquench, quenching lookback time; tform, formation lookback time; AV, effective dust attenuation optical depth.

Figure 2 shows the SFH of the galaxy, as inferred by the four codes. All models agree that JADES-GS-z7-01-QU is quenched and give similar stellar population parameters. The oldest notable population of stars is 40–150 Myr old, corresponding to a formation redshift z = 7.6–8.8, whereas the youngest stars have ages 20–50 Myr, corresponding to a quenching redshift of z = 7.4–7.7. These numbers imply that JADES-GS-z7-01-QU formed in a burst of star formation lasting only 20–100 Myr, consistent with the formation timescales of star-forming galaxies at similar redshifts2.Fig. 2 The SFH of the galaxy as inferred by four different full spectral fitting codes, which use different (effective) priors on the SFH of the galaxy.

All four codes confirm that the galaxy is quenched at the epoch of observation and reconstruct comparable SFHs. a, The stellar age–metallicity grid resulting from the pPXF fit. The code reconstructs dominant metal-poor populations forming from approximately 100 Myr to approximately 20 Myr before observation. The colour bar represents the fractional mass distribution over the SSP grid. b, The SFH inferred by BAGPIPES. The solid green line shows the median posterior, the shaded region shows the 1σ range, indicating a single star-formation burst lasting approximately 20 Myr and quenching approximately 20 Myr before observation. c, The SFH inferred by BEAGLE, which suggests that the galaxy formed approximately 100 Myr before the epoch of observation and quenched approximately 10–20 Myr before observation. d, The SFH inferred by Prospector, which suggests that the galaxy quenched approximately 20–30 Myr before observation after a starburst lasting approximately 50 Myr.

The SFR at the time of observation inferred by BAGPIPES, BEAGLE and Prospector are extremely low, between 10−2.6 and 10−1.4 M⊙ yr−1, yielding specific SFRs ranging between 10−2.3 Gyr−1 and 0.1 Gyr−1. These values are between 2 and 3 orders of magnitude below the main sequence of star-forming galaxies at this redshift32–36 and below the widely used threshold sSFR10 < 0.2/tH = 0.29 Gyr−1, on 10-Myr timescales, hence confirming that the galaxy is quenched at the epoch of observation. Crucially, the four codes agree that the galaxy has been strongly star forming between 10 and 100 Myr before the epoch of observation.

Three of the four codes infer a tentative low average stellar metallicity of the galaxy of log10(Z/Z⊙) ≈ −2 (in which Z⊙ is the solar metallicity), whereas BAGPIPES infers log10(Z/Z⊙) ≈ −0.7. pPXF indicates the presence of a weak enriched population representing only 5% of the total stellar mass of the galaxy, which formed last before quenching. However, we note that stellar metallicity measurements are uncertain with the low-resolution prism spectroscopy.

Which physical mechanism(s) quenched the galaxy?

The inferred mass of this galaxy rules out that it has been quenched by the UV background37; indeed, numerical simulations predict that this quenching mechanism works only for very-low-mass galaxies with M⋆ ≈ 105–107 M⊙ (maximally < 108 M⊙)38.

In the local Universe, galaxies in the mass range of our target are quenched primarily by environment39,40. It has been postulated that some satellite galaxies may experience environment-driven quenching already during the epoch of reionization41. However, we do not find any massive galaxies nearby (see Methods), disfavouring environmental effects as the quenching mechanism for this target.

Given the short inferred duration of the SFH and the rapidity of the transition to quiescence, it seems more reasonable to speculate that JADES-GS-z7-01-QU may have experienced a powerful outflow, driven by either star-formation feedback (radiation-pressure, supernovae might act too slowly) or accretion on a primeval supermassive black hole, which rapidly ejected most of the star-forming gas42. This scenario is supported by the tentative low average stellar metallicity inferred by three of the codes. Indeed, ejective feedback mechanisms might have rapidly removed gas from the galaxy and quenched it, before the ISM could be substantially enriched with new metals. A slower quenching process (such as the starvation scenarios) would have probably resulted in a longer transition between star forming and quenched and into higher-metallicity stellar populations, formed out of recycled gas produced by stellar evolution and returned to the ISM by means of supernovae43,44.

These outflow events, driven by either star formation or active galactic nucleus, might have mini-quenched star formation only temporarily45, until new or re-accreted material replenishes the supply of gas available for star formation and rejuvenates the galaxy. The latter picture may be qualitatively in agreement with a wide range of cosmological simulations predicting that a population of galaxies in the early Universe goes through periodic bursts of star formation, interspersed with periods of suppressed star formation45–47. Although the expected SFHs are very ‘bursty’, these recent simulations struggle to achieve the complete quenching observed by us for galaxies with mass similar to our system.

More generally, interpreting these observations with existing simulations is complicated because, according to current theories47, this object occupies the transition region between bursty and stable SFHs. Moreover, it is important to note that these models do not include active galactic nucleus feedback, which recent observations have shown to be important in local galaxies of this mass range48. These difficulties mean that JADES-GS-z7-01-QU provides the community with the opportunity to shed light on this pivotal mass range.

We conclude by emphasizing that the discovery and spectroscopic analysis of a (mini-)quenched galaxy at redshift z = 7.3 by our JADES collaboration ushers the era in which we can constrain theoretical feedback models using direct observations of the primordial Universe. However, this is just the starting point for the JWST mission: upcoming and future observations will start the transition from the ‘discovery’ phase to the statistical characterization of the properties of the first (mini-)quenched galaxies.

Methods

JWST/NIRSpec spectra

The NIRSpec49 prism/R100 and gratings/R1000 spectra of JADES-GS-z7-01-QU presented in this work were obtained as part of our JADES GTO programme (PI: N. Lützgendorf, ID: 1210) observations in the Great Observatories Origins Deep Survey South (GOODS-S) field between 21 and 25 October 2022. The R100 observations were obtained using the disperser/filter configuration PRISM/CLEAR, which covers the wavelength range between 0.6 μm and 5.3 μm and provides spectra with a wavelength-dependent spectral resolution of R ≈ 30–330. The R100 spectrum of JADES-GS-z7-01-QU is presented in Fig. 1.

The medium-resolution R1000 observations, with a spectral resolution of R ≈ 500–1,340, used the disperser/filter configurations G140M/F070LP, G235M/F170LP and G395M/F290LP, which were exposed for 14 h, 7 h and 7 h. A zoom-in on the R1000 spectrum (into the region with spectral lines best tracing star-formation activity) is shown in Extended Data Fig. 1. Finally, high-resolution R2700 observations used G395H/F290LP and were exposed for 7 h (like the R1000 spectrum, the R2700 spectrum of JADES-GS-z7-01-QU contains no detections, hence is not shown).

The programme observed a total of 253 galaxies over three dither pointings, with JADES-GS-z7-01-QU being observed in each of the three pointings. Each dither pointing had a different microshutter array (MSA) configuration to place the spectra at different positions on the detector to decrease the impact of detector gaps, mitigate detector artefacts and improve the signal-to-noise ratio for high-priority targets, while increasing the density of observed targets. Within each individual dither pointing, the telescope executed a three-nod pattern (by slightly reorienting the telescope by the length of one microshutter, keeping the same MSA configuration). In each of the three nodding pointings, three microshutters were opened for each target, with the targets in the central shutter. Each three-point nodding was executed within 8,403 s. The nodding pattern has been repeated four times in the PRISM/CLEAR configuration, two times in the G140M/F070LP combination, once in the G235M/F170LP combination and once in the G395M/F290LP combination. This resulted in a total exposure time for JADES-GS-z7-01-QU of 28 h in R100, 14 h in G140M and 7 h in each of G235M, G395M and G395H.

The flux-calibrated spectra were extracted using a customized pipeline developed by the NIRSpec GTO team, which builds on the publicly available ESA NIRSpec Science Operations Team (SOT) pipeline50. A detailed description of the custom pipeline will be presented in a forthcoming technical paper (Carniani et al., in preparation) and more information can be found in ref. 51. We summarize here the main steps and the differences to the publicly available pipeline. For each exposure, we extract the count rate for each pixel, removing cosmic rays and flagging saturation. The 2D spectrum is background subtracted on the basis of the two other exposures in the three-nod pattern. The individual 2D spectra are flat-fielded and illumination-corrected, taking into account the wavelength-dependent throughput. The wavelength and flux calibration was then applied, with each pixel of the 2D spectrum having an associated wavelength and position along the shutter. We applied a wavelength-dependant path-loss correction to account for flux falling outside the microshutter, taking into account the considerable point spread function variation of NIRSpec, treating the target as a point source. For the prism, we used an irregular spectral wavelength grid, taking into account the resolution (R) as a function of wavelength. The 1D spectra for the three nod positions from each of the three pointings are combined by a weighted average into a single 1D spectrum. Outliers are rejected with a sigma-clipping algorithm. The presented 1D spectra come from a combination of the 1D individual spectra and are not an extraction from the presented combined 2D spectra.

JWST/NIRCam image and morphology

A JWST/NIRCam F444W-F200W-F090W rgb (red-green-blue) colour image of JADES-GS-z7-01-QU from our JADES programme (PI: Daniel J. Eisenstein, ID: 1180), created from cutouts of the mosaics in each filter, at wavelengths λ ≈ 0.8–5 μm, is shown in Extended Data Fig. 2.

For the spectrophotometric modelling of JADES-GS-z7-01-QU, we used the photometry from the JADES and JEMS52 NIRCam53,54 surveys. In particular, the modelling included deep infrared NIRCam observations with the following filters: F090W, F115W, F150W, F182M, F200W, F210M, F277W, F335M, F356W, F410M, F430M, F444W, F460M and F480M. The JADES photometry reduction pipeline made use of the JWST Calibration Pipeline (JWSTCP, v1.9.2) with the CRDS pmap context 1039. The raw images were transformed into count-rate images, making use of JWSTCP stage 1, for which detector-level corrections and ‘snowballs’ were accounted. The count-rate images were then flat-fielded and flux-calibrated with a customized methodology, using JWSTCP stage 2. Finally, the mosaics were created using stage 3 of the pipeline. For further details on the JADES photometry data reduction pipeline, see refs. 55,56.

To obtain the morphological parameters of JADES-GS-z7-01-QU, we fit the NIRCam photometry with Forcepho (Johnson et al., in preparation). Forcepho models galaxies and substructures (for example, clumps or blended companions) as several Sérsic profiles convolved with the instrument point spread functions as mixtures of Gaussians by forward-modelling the light distribution in all individual exposures and filters and sampling the joint posterior probability distribution of all parameters through Markov chain Monte Carlo. For more details on the multicomponent modelling procedure, see ref. 56. JADES-GS-z7-01-QU appears as a compact, discy galaxy (half-light radius Re = 36 ± 1 mas ≙ 0.2 kpc ≙ 0.04 arcsec, Sérsic index n = 0.95 ± 0.03; Extended Data Fig. 2). The images also show a distinct, fainter source 0.13 arcsec to the east. This secondary source could not be deblended in the spectroscopy but we obtained deblended photometry using Forcepho. The contribution of the secondary source to the total flux ranges from a maximum of 27% (in the F115W band) to 17% (in the F444W band), therefore its SED is much bluer than that of the main source. Its photometric redshift z = 7.50 ± 0.13 (1σ) is consistent with the spectroscopic redshift of the main source. At a redshift of z = 7.3, this secondary source would lie within 0.7 kpc (or 3Re) from the centre of JADES-GS-z7-01-QU; its interpretation as a clump or satellite is unclear. To attempt removing its contribution from the spectrum of the main source, we extracted a spectrum from the central three pixels (0.3 arcsec) from the NIRSpec 2-d spectrum; using this spectrum does not change the interpretation of our results, that is, JADES-GS-z7-01-QU is still quenched.

As discussed in the main text, quenching by environment is ruled out for JADES-GS-z7-01-QU, as no other galaxy resides nearby. This can be verified with JADES NIRCam imaging on our publicly available website and more specifically the interactive tool FitsMap: https://jades.idies.jhu.edu/public/?ra=53.1554497&dec=-27.8018917&zoom=9 at the coordinates RA = 53.1551 and dec. = −27.8018.

Full spectral fitting

pPXF

The red model fit of the stellar continuum in Fig. 1 was performed with the χ2-minimization Penalized PiXel-Fitting code pPXF (refs. 57,58), using a library of single stellar population (SSP) templates spectra obtained combining the synthetic C3K model atmospheres59 with MIST isochrones60 and solar abundances. The SSP spectra span a full 2D logarithmic grid of 62 ages and 10 metallicities from ageSSP = 106.0 years to 109.2 years (generously older than the age of the Universe at z = 7.3) and log10(Z/Z⊙)SSP = −2.5 to 0.5. Owing to the low resolution of the R100 spectrum, we fix the stellar velocity dispersion to its virial estimate σ*≈σvir≡GM*/(5Re)=50kms−1. To account for dust reddening, the fitted SSP are multiplicatively coupled to the dust attenuation curve in ref. 61. To infer the stellar population weight-grid shown in Fig. 2a, following ref. 62, we first convolve the SSP templates to match the wavelength-dependant spectral resolution of the prism spectrum. Then, to avoid numerical problems, both the spectrum and the templates are renormalized by the median flux per spectral pixel. Then we run an initial fit with pPXF and we σ-clip outliers in the spectrum. Finally, we perform a residual-based bootstrapping of the initial pPXF best fit, without regularization57,58, over 1,000 iterations. The inferred bootstrapped SSP grids are averaged to recover the non-parametric SFH, consistent with the intrinsic noise of the spectrum, presented in Fig. 2a.

We infer a dust attenuation of the stars in this galaxy of AV = 0.4 ± 0.1. It should be noted that the presence of dust in the pPXF fit is mainly driven by the UV slope. The complex physics of the Lyα drop is not included in the SSP templates. Masking this part of the spectrum returns a nearly dust-free fit with older and metal-richer stellar populations, which would make JADES-GS-z7-01-QU even more quenched. As stated in the main text, we infer an extremely low average stellar metallicity of log10(Z/Z⊙) ≈ −2 with pPXF. It should be noted that the dominant reconstructed stellar populations lie at log10(Z/Z⊙) ≈ −2.5, at the boundary of the available grid of synthetic spectra. This suggests that model SSP spectra of even lower metallicity might be needed in the future to accurately model the stellar populations in galaxies at high redshift. However, we note that the metallicity measurements are uncertain, owing to the low resolution of the prism. We infer that the oldest notable population of stars (that is, indicating the start of the star formation) in the galaxy is 150 Myr old, whereas the youngest is 50 Myr old, resulting in an extremely short duration of the star formation of just 100 Myr between the formation of the galaxy and its quenching.

BAGPIPES

We used the Bayesian Analysis of Galaxies for Physical Inference and Parameter EStimation (BAGPIPES) code63 to simultaneously fit the NIRSpec PRISM measurements and NIRCam photometry. Following ref. 64, we used the updated BC03 stellar population models65,66 combined with the stellar MILES library67 and the updated stellar evolution tracks68,69. For the presented BAGPIPES fit, we assumed two bins of constant SFH, one fixed bin over the past 10 Myr and one variable bin spanning a range beyond 10 Myr (minimum age ranging between 10 Myr and 0.5 Gyr, maximum age between 11 Myr and the age of the Universe). We varied the total stellar mass formed between 0 and 1015 M⊙ and the stellar metallicity of the variable SFH bin between 0.01 Z⊙ and 1.5 Z⊙ (the 10-Myr bin having a metallicity of 0.2 Z⊙ to match the inferred metallicity of the variable-SFH bin). Nebular emission is modelled self-consistently with a grid of CLOUDY70 models with the ionization parameter (−3 < log10U < −0.5) as a free parameter. We included a flexible dust attenuation prescription71 with visual extinction and power-law slope freely varying (0 < AV < 7, 0.4 < n < 1.5) while fixing the fraction of attenuation from stellar birth clouds to 60% (the remaining fraction arising in the diffuse ISM; ref. 72). A first-order correction polynomial73 is fitted to the spectroscopic data to account for aperture and flux calibration effects. The spectrophotometric fit and the corresponding corner plot are shown in Extended Data Fig. 3. We find that nearly no wavelength-dependant correction is necessary at the blue end of the spectrum, whereas at the red end, a correction of 15% is applied. Crucially, we find a very low SFR (consistent with 0) in the past 10 Myr for JADES-GS-z7-01-QU, noting that other tested SFH parametrizations, namely the double-power-law SFH described in ref. 74 and a single-bin constant SFH with flexible beginning and end of star formation, return consistent results and most crucially agree that the galaxy is quenched. We infer that the oldest stellar population is 40 Myr old, which is equivalent to a formation redshift of z = 7.6. The galaxy has been quenched for 10 Myr, resulting in a short duration of star formation of 20 Myr from the formation of the galaxy to its quenching.

BEAGLE

We use the Bayesian analysis tool BEAGLE (ref. 66) to fit to the R100/prism spectrum of JADES-GS-z7-01-QU. The BEAGLE code incorporates a consistent modelling of stellar radiation and its transfer through the interstellar and intergalactic media. We model the SFH as an initial delayed exponential with maximum stellar age, tform (years), and location of the peak of star formation as free parameters. To disentangle the current SFR from the integrated property of total stellar mass, we allow for the most recent episode of star formation to be modelled as a constant with free parameters SFR (M⊙ yr−1) and duration, tquench (years) (which can vary between 107 and 108 years). The nebular emission is characterized by the interstellar metallicity, the ionization parameter, the mass fraction of interstellar metals locked within dust grains and, crucially, fesc (which can vary between 0 and 1). Dust attenuation follows the two-component prescription of ref. 71, in which we fit for the total effective V-band attenuation optical depth (fixing the ratio of V-band ISM attenuation to the V-band ISM + birth cloud attenuation to 0.4). We also fit for stellar metallicity, stellar mass formed and redshift, totalling 12 free parameters. A list of the free parameters and the adopted priors is presented in Extended Data Table 1.

The corner plot in Extended Data Fig. 4 shows the BEAGLE posterior probability distributions of the BEAGLE fit. The 2D (off-diagonal) and 1D (along the main diagonal) subplots show the posterior distributions on stellar mass M⋆, metallicity Z, SFR, maximum age of stars tform, minimum age of stars tquench, redshift z, effective dust attenuation optical depth in the V-band AV and the escape fraction of ionizing photons fesc. The dark, medium and light blue contours show the extents of the 1σ, 2σ and 3σ credible regions.

BEAGLE gives a current SFR of less than 10−1.5 M⊙ yr−1, a formation time of less than 160 Myr before observation and a quenching time of roughly 15 Myr before observation.

We also note that BEAGLE, as for the other three codes, requires some degree of dust attenuation, which suggests that some cold gas is still present, which—in turn—is incompatible with fesc ≈ 1.

Prospector

We use the Bayesian SED fitting code Prospector75 to model the spectrophotometric data of JADES-GS-z7-01-QU. The posterior corner plot for several key parameters from Prospector is shown in Extended Data Fig. 5. The code uses a flexible spectroscopic calibration model, combined with forward modelling of spectra and photometry, to infer physical properties. Following the setup in ref. 76, we include a flexible SFH (ten bins with the bursty continuity prior), a flexible attenuation law (diffuse dust optical depth with a power-law modifier to shape the attenuation curve of the diffuse dust in ref. 61) and fit for the stellar metallicity. Notably, Prospector infers a low-dust attenuation with AV=0.1−0.0+0.1 with a rather steep attenuation law AUV/AV=2.6−0.8+1.4. This is consistent with the idea that the galaxy has a low gas content and the low SFR in the past 30 Myr before observation. Prospector infers that the oldest stellar population (as defined by the lookback time when the first 10% of the stellar mass formed) has an age of about 100 Myr, which means a nominal formation redshift of z = 8.8. The SFR increases markedly approximately 80 Myr before observation. After this final burst, lasting around 50 Myr, the galaxy quenched on a short timescale.

We have also experimented with the standard continuity prior77, which weights against sharp transition in the SFH. The overall shape of the SFH is the same, indicating that the data strongly prefer a decreasing SFH in the past roughly 50 Myr. Quantitatively, the recent SFR (averaged over the past 10 Myr) increases with this prior to log10SFRM⊙year−1=−0.4−0.9+0.4, which is still consistent with being quenched and within the uncertainties of the fiducial value obtained with the bursty continuity prior. The quenching time is slightly more recent 24−9+6Myr, but consistent within the uncertainties quoted in Table 1.

Star-forming, high-fesc interpretation

It should be noted that the complete absence of nebular lines always allows, by construction, a solution with fesc ≈ 1 (regardless of whether the galaxy has been recently star forming or quiescent)—the question is whether this solution is accompanied by the production of ionizing photons associated with continuing star formation.

The fiducial BEAGLE posterior distribution does not highlight a solution with high fesc and very recent star formation26,29,78. By contrast, although fesc is unconstrained, even a value approaching unity indicates a low SFR < 0.1 M⊙ yr−1 at the 3σ level (fifth subplot from the left at the bottom of Extended Data Fig. 4).

To assess the very recently star-forming and high-fesc scenario quantitatively, we use BEAGLE to compare two SED models. The model already described (see the ‘BEAGLE’ section) formally allows a star-forming solution with high fesc. The alternative model has a simplified SFH consisting of a constant SFR; in this way, low-SFR solutions are effectively removed by the constraint to form sufficient stellar mass of the appropriate age to reproduce the observed spectrum. This alternative model gives fesc=0.98−0.04+0.01 and SFR=0.63−0.05+0.05M⊙year−1, which is a much higher SFR than the alternative solution. To select the preferred model, we use the Bayes ratio, that is, the ratio between the evidence of the models. The log difference between the evidences, that is, the Bayes factor, is ln(K) = 4.1 ± 0.3; according to Jeffreys’ criterion79, this is strong evidence for the quenched solution and we adopt it as our fiducial model.

As an extra test, we assumed a model with the same setup as the fiducial run, but forcing the escape fraction to fesc > 0.9. We find that the result is equal to the fiducial run and the galaxy remains quenched.

Empirical measurements

To estimate the flux upper limits on Hβ and [O III]λ5008, we sum the formal variance over three pixels. For EWHδA we use the bands in the Lick definition80 but without any further correction owing to spectral resolution.

We derive an upper limit on the SFR from the 3σ upper limit on the Hβ emission-line flux, F(Hβ) < 6.1 × 10−20 erg cm−2 s−1. To correct this flux for dust attenuation, we assume the Milky Way attenuation law81, which seems appropriate for galaxies at least up until z = 2.5 (refs. 82,83). Given that Hβ is not detected, we cannot measure the Balmer decrement. We therefore derive the nebular AV from the continuum AV = 0.51 mag inferred from BEAGLE (the highest value between all models) and upscale this value by 0.64, the median continuum-to-nebular AV ratio inferred from local galaxies84 (of stellar mass comparable with JADES-GS-z7-01-QU). The flux is converted to a luminosity assuming the Planck18 cosmology85. To convert the Hβ attenuation-corrected luminosity to a SFR, we use the conversion factor 2.1 × 10−42 M⊙ yr−1 erg−1 s, appropriate for a Chabrier initial mass function with a high-mass cutoff of 100 M⊙ and metallicity Z = 0.27 Z⊙ (ref. 83) (note that this value of the metallicity is higher than that inferred from the data; this provides a conservative estimate). This gives a SFR of 0.57 M⊙ yr−1. Even stronger constraints come from the [O III]λ5008 line: we find F([O III]λ5008) < 6.5 × 10−20 erg s−1 cm−2, which, combined with a conservative assumption of [O III]λ5008/Hβ ratios in high-z galaxies14,15, implies a 3σ limit on the SFR roughly five times lower than the Hβ-derived value (SFR = 0.12 M⊙ yr−1).

Alternatively, assuming the median (and the extreme) observed Balmer decrement 3.5 (5.5) from ref. 83, we would obtain nebular AV values of 0.63 and 2.05 mag, respectively. These translate into [O III]λ5008-derived SFRs of 0.10 and 0.34 M⊙ yr−1, respectively. As a comparison, the SFR threshold for quiescence at z = 7.3 is 0.18 M⊙ yr−1 (obtained from the threshold in sSFR defined by 0.2/tH(z) times the BEAGLE stellar mass (ref. 86)). Thus, in all but the most extreme scenario, JADES-GS-z7-01-QU would meet the formal threshold for quiescence. The absence of emission lines is independently confirmed by the medium-resolution spectrum (see Extended Data Fig. 1).

Online content

Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at 10.1038/s41586-024-07227-0.

Extended data figures and tables

Extended Data Fig. 1 NIRSpec R1000/grating spectrum of the (mini-)quenched galaxy JADES-GS-z7-01-QU at z = 7.3.

The spectrum confirms the absence of emission lines. The blue-shaded region shows the 1D noise level. The upper panel shows the signal-to-noise ratio (S/N) in the 2D grating spectrum. The spectrum is median-smoothed, for visualization.

Extended Data Fig. 2 JWST/NIRCam image covering JADES-GS-z7-01-QU and its nearby projected environment.

The NIRCam F444W-F200W-F090W rgb colour image is created from cutouts of the mosaics at wavelengths λ ≈ 0.8–5 μm. The five NIRSpec microshutter positions used for this target are overlaid in white.

Extended Data Fig. 3 Summary of key outputs by BAGPIPES.

Bottom left, corner plot. Top right, spectrophotometric BAGPIPES fit of the JADES-GS-z7-01-QU R100/prism spectrum.

Extended Data Fig. 4 Summary of key outputs by BEAGLE.

Bottom left, corner plot with free fesc. Top right, BEAGLE maximum a posteriori model of the R100 spectrum.

Extended Data Fig. 5 Summary of key outputs by Prospector.

Bottom left, corner plot with stellar mass M⋆, SFR, tform, tquench, dust attenuation AV, AUV/AV and stellar metallicity Z. Top right, Prospector maximum a posteriori model of the R100 spectrum.

Extended Data Table 1 Parameters and associated priors set in BEAGLE

Parameters and associated priors set in BEAGLE

N(a,b) is the normal distribution with mean a and standard deviation b. †In practice, BEAGLE will not allow the age of the oldest stars to be greater than the time between z = 20 and the sampled redshift.

Acknowledgements

We thank the referees for their very helpful comments, based on which we have substantially improved the discussion and presentation of the results. T.J.L., F.D.E., R.M., J.W., W.M.B., L.S. and J.S. acknowledge support by the Science and Technology Facilities Council (STFC), by the European Research Council (ERC) through Advanced Grant 695671 ‘QUENCH’ and by the UKRI Frontier Research grant RISEandFALL. T.J.L. acknowledges support by the STFC Center for Doctoral Training in Data Intensive Science. R.M. also acknowledges funding from a research professorship from the Royal Society. J.W. further acknowledges support from the Fondation MERAC. This study made use of the Prospero high-performance computing facility at Liverpool John Moores University. B.D.J., E.E., M.R. and B.R. acknowledge support from the JWST/NIRCam Science Team contract to the University of Arizona, NAS5-02015. E.C.-L. acknowledges support of a STFC Webb Fellowship (ST/W001438/1). S. Carniani acknowledges support by the European Union’s HE ERC Starting grant no. 101040227 – WINGS. A.J.B., J.C., A.J.C., A.S. and G.C.J. acknowledge funding from the ‘FirstGalaxies’ Advanced Grant from the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 789056). R.S. acknowledges support from a STFC Ernest Rutherford Fellowship (ST/S004831/1). S. Arribas, B.R.D.P. and M.P. acknowledge support from the research project PID2021-127718NB-I00 of the Spanish Ministry of Science and Innovation/State Agency of Research (MICIN/AEI). The Cosmic Dawn Center (DAWN) is financed by the Danish National Research Foundation under grant no. 140. H.Ü. gratefully acknowledges support by the Isaac Newton Trust and by the Kavli Foundation through a Newton-Kavli Junior Fellowship. D.J.E. is supported as a Simons Investigator and by JWST/NIRCam contract to the University of Arizona, NAS5-02015. Funding for this research was provided by the Johns Hopkins University, Institute for Data Intensive Engineering and Science (IDIES). The research of C.C.W. is supported by NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. This research is supported in part by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. We thank C. Lovell, A. Pillepich, D. Sijacki and N. Laporte for helpful comments and discussions.

Author contributions

T.J.L., F.D.E. and R.M. led the writing of the paper. All authors have contributed to the interpretation of the results. T.J.L., J.W., L.S., S.T., F.D.E., E.C.-L., J.C., B.D.J., W.M.B. and K.A.S. led the spectro(photometric) modelling of JADES-GS-z7-01-QU and the data visualization. A.J.B., A.d.G., C.C.W., C.W., D.J.E., H.-W.R., M.R., P.F., R.M., S. Alberts and S. Arribas contributed to the design of the JADES survey. C.W. contributed to the design of the spectroscopic observations and MSA configurations. E.C.-L., R.E. and R.H. contributed to the photometric redshift determination and target selection. A.J.C., A.J.B., C.W., E.C.-L., H.Ü. and K.B. contributed to the selection, prioritization and visual inspection of the targets. S. Carniani, M.C., J.W., P.F., S. Arribas and B.R.D.P. contributed to the NIRSpec data reduction and to the development of the NIRSpec pipeline. S. Arribas, S. Charlot, J.C., M.C., F.D.E., A.d.G., E.C.-L., M.V.M., R.M., B.R.D.P., T.J.L., A.S., L.S., J.S., R.S. and J.W. contributed to the development of the JADES tools for the spectroscopic data analysis. B.R., S.T., B.D.J., C.C.W., D.J.E., I.S., M.R., R.E. and Z.C. contributed to the JADES imaging data reduction. C.C.W., S.T., M.V.M., B.R., B.D.J., C.W., D.J.E., Z.J., A.S., K.B., A.J.B., S. Carniani, S. Charlot, J.C., E.C.-L., A.d.G., E.E., N.K., R.M., E.N., M.R., L.S., I.S., R.S., K.A.S. and H.Ü. contributed to the JEMS survey. R.H. and B.R. contributed to the JADES imaging data visualization. B.D.J., S.T., A.d.G. and R.E. contributed to the modelling of galaxy photometry. N.B. and S. Arribas contributed to the design and optimization of the MSA configurations. P.F., T.R., N.K. and B.R.D.P. contributed to the design, construction and commissioning of NIRSpec. M.R., C.W., E.E. and C.C.W. contributed to the design, construction and commissioning of NIRCam. B.R., C.W., D.J.E., M.R. and R.M. serve as the JADES Steering Committee.

Peer review

Peer review information

Nature thanks the anonymous reviewers for their contribution to the peer review of this work.

Data availability

The reduced spectra that support the findings of this study are publicly available on GitHub: https://github.com/tobiaslooser/JWST-reveals-a-recently-mini-quenched-galaxy-at-z-7.3. See MAST at Space Telescope Science Institute for the original data: https://archive.stsci.edu/hlsp/jades.

Code availability

The pPXF, BAGPIPES and Prospector codes are publicly available. BEAGLE is available through a Docker image (distributed through Docker Hub) on request at https://www.iap.fr/beagle/.

Competing interests

The authors declare no competing interests.

Extended data

is available for this paper at 10.1038/s41586-024-07227-0.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Endsley R A JWST/NIRCam study of key contributors to reionization: the star-forming and ionizing properties of UV-faint z ~ 7–8 galaxies Mon. Not. R. Astron. Soc. 2023 524 2312 2330 10.1093/mnras/stad1919
2. Tacchella S JWST NIRCam + NIRSpec: interstellar medium and stellar populations of young galaxies with rising star formation and evolving gas reservoirs Mon. Not. R. Astron. Soc. 2023 522 6236 6249 10.1093/mnras/stad1408
3. Mason CA Trenti M Treu T The brightest galaxies at cosmic dawn Mon. Not. R. Astron. Soc. 2023 521 497 503 10.1093/mnras/stad035
4. Whitler L Star formation histories of UV-luminous galaxies at z ≃ 6.8: implications for stellar mass assembly at early cosmic times Mon. Not. R. Astron. Soc. 2023 519 5859 5881 10.1093/mnras/stad004
5. Endsley, R. et al. The star-forming and ionizing properties of dwarf z~6-9 galaxies in JADES: insights on bursty star formation and ionized bubble growth. Preprint at https://arxiv.org/abs/2306.05295 (2023).
6. Dome, T. et al. Mini-quenching of z = 4 – 8 galaxies by bursty star formation. Mon. Not. R. Astron. Soc. 527, 2139–2151 (2024).
7. Merlin E Red and dead CANDELS: massive passive galaxies at the dawn of the Universe Mon. Not. R. Astron. Soc. 2019 490 3309 3328 10.1093/mnras/stz2615
8. Lovell, C. C. et al. First Light And Reionisation Epoch Simulations (FLARES) VIII. The emergence of passive galaxies at z ≥ 5. Preprint at https://arxiv.org/abs/2211.07540 (2023).
9. Glazebrook K A massive, quiescent galaxy at a redshift of 3.717 Nature 2017 544 71 74 10.1038/nature21680 28382981
10. Carnall, A. C. et al. A surprising abundance of massive quiescent galaxies at 3 < z < 5 in the first data from JWST CEERS. Preprint at https://arxiv.org/abs/2208.00986 (2023).
11. Valentino F Quiescent galaxies 1.5 billion years after the Big Bang and their progenitors Astrophys. J. 2020 889 93 10.3847/1538-4357/ab64dc
12. Nanayakkara T A population of faint, old, and massive quiescent galaxies at 3 < z < 4 revealed by JWST NIRSpec Spectroscopy Sci. Rep. 2024 14 3724 10.1038/s41598-024-52585-4 38355772
13. Oesch PA z ~ 7 galaxies in the HUDF: first epoch WFC3/IR results Astrophys. J. Lett. 2010 709 L16 L20 10.1088/2041-8205/709/1/L16
14. Cameron AJ JADES: probing interstellar medium conditions at z ~ 5.5–9.5 with ultra-deep JWST/NIRSpec spectroscopy Astron. Astrophys. 2023 677 A115 10.1051/0004-6361/202346107
15. Sanders RL Shapley AE Topping MW Reddy NA Brammer GB Excitation and ionization properties of star-forming galaxies at z = 2.0–9.3 with JWST/NIRSpec Astrophys. J. 2023 955 54 10.3847/1538-4357/acedad
16. Bhatawdekar R Conselice CJ UV spectral slopes at z = 6–9 in the Hubble Frontier Fields: lack of evidence for unusual or population III stellar populations Astrophys. J. 2021 909 144 10.3847/1538-4357/abdd3f
17. Bunker AJ JADES NIRSpec spectroscopy of GN-z11: Lyman-α emission and possible enhanced nitrogen abundance in a z = 10.60 luminous galaxy Astron. Astrophys. 2023 677 A88 10.1051/0004-6361/202346159
18. Goto T A catalogue of local E+A (post-starburst) galaxies selected from the Sloan Digital Sky Survey Data Release 5 Mon. Not. R. Astron. Soc. 2007 381 187 193 10.1111/j.1365-2966.2007.12227.x
19. Wu P-F Fast and slow paths to quiescence: ages and sizes of 400 quiescent galaxies from the LEGA-C survey Astrophys. J. 2018 868 37 10.3847/1538-4357/aae822
20. Hashimoto T The onset of star formation 250 million years after the Big Bang Nature 2018 557 392 395 10.1038/s41586-018-0117-z 29769675
21. Roberts-Borsani GW Ellis RS Laporte N Interpreting the Spitzer/IRAC colours of 7 ≤ z ≤ 9 galaxies: distinguishing between line emission and starlight using ALMA Mon. Not. R. Astron. Soc. 2020 497 3440 3450 10.1093/mnras/staa2085
22. Laporte N Probing cosmic dawn: ages and star formation histories of candidate z ≥ 9 galaxies Mon. Not. R. Astron. Soc. 2021 505 3336 3346 10.1093/mnras/stab1239
23. Williams RJ Quadri RF Franx M van Dokkum P Labbé I Detection of quiescent galaxies in a bicolor sequence from Z = 0–2 Astrophys. J. 2009 691 1879 1895 10.1088/0004-637X/691/2/1879
24. Belli S Newman AB Ellis RS MOSFIRE spectroscopy of quiescent galaxies at 1.5 < z < 2.5. II. Star formation histories and galaxy quenching Astrophys. J. 2019 874 17 10.3847/1538-4357/ab07af
25. Park M Rapid quenching of galaxies at cosmic noon Astrophys. J 2023 953 119 10.3847/1538-4357/acd54a
26. Zackrisson E The spectral evolution of the first galaxies. III. Simulated James Webb Space Telescope spectra of reionization-epoch galaxies with Lyman-continuum Leakage Astrophys. J. 2017 836 78 10.3847/1538-4357/836/1/78
27. Rosdahl J LyC escape from SPHINX galaxies in the Epoch of Reionization Mon. Not. R. Astron. Soc. 2022 515 2386 2414 10.1093/mnras/stac1942
28. Katz H Two modes of LyC escape from bursty star formation: implications for [C II] deficits and the sources of reionization Mon. Not. R. Astron. Soc. 2023 518 270 285 10.1093/mnras/stac3019
29. Topping MW Searching for extremely blue UV continuum slopes at z = 7–11 in JWST/NIRCam imaging: implications for stellar metallicity and ionizing photon escape in early galaxies Astrophys. J. 2022 941 153 10.3847/1538-4357/aca522
30. Baron, D. et al. Star formation and molecular gas properties of post-starburst galaxies. Preprint at https://arxiv.org/abs/2204.11881 (2022).
31. Sandles, L. et al. JADES: Balmer decrement measurements at redshifts 4 < z < 7. Preprint at https://arxiv.org/abs/2306.03931 (2023).
32. Steinhardt CL Star formation at 4 < z < 6 from the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH) Astrophys. J. Lett. 2014 791 L25 10.1088/2041-8205/791/2/L25
33. Pearson WJ Main sequence of star forming galaxies beyond the Herschel confusion limit Astron. Astrophys. 2018 615 A146 10.1051/0004-6361/201832821
34. Sandles L Curtis-Lake E Charlot S Chevallard J Maiolino R Bayesian hierarchical modelling of the M*–SFR relation from 1 ≲ z ≲ 6 in ASTRODEEP Mon. Not. R. Astron. Soc. 2022 515 2951 2969 10.1093/mnras/stac1999
35. Santini P The star formation main sequence in the Hubble Space Telescope Frontier Fields Astrophys. J. 2017 847 76 10.3847/1538-4357/aa8874
36. Laporte N Ellis RS Witten CEC Roberts-Borsani G Resolving ambiguities in the inferred star formation histories of intense [O III] emitters in the reionization era Mon. Not. R. Astron. Soc. 2023 523 3018 3024 10.1093/mnras/stad1597
37. Efstathiou G Suppressing the formation of dwarf galaxies via photoionization Mon. Not. R. Astron. Soc. 1992 256 43P 47P 10.1093/mnras/256.1.43P
38. Katz H How to quench a dwarf galaxy: the impact of inhomogeneous reionization on dwarf galaxies and cosmic filaments Mon. Not. R. Astron. Soc. 2020 494 2200 2220 10.1093/mnras/staa639
39. Peng Y-j Mass and environment as drivers of galaxy evolution in SDSS and zCOSMOS and the origin of the Schechter function Astrophys. J. 2010 721 193 221 10.1088/0004-637X/721/1/193
40. Bluck AFL How do central and satellite galaxies quench? – Insights from spatially resolved spectroscopy in the MaNGA survey Mon. Not. R. Astron. Soc. 2020 499 230 268 10.1093/mnras/staa2806
41. Gelli V Salvadori S Ferrara A Pallottini A Carniani S Dwarf satellites of high-z Lyman break galaxies: a free lunch for JWST Astrophys. J. Lett. 2021 913 L25 10.3847/2041-8213/abfe6c
42. Koudmani S Sijacki D Bourne MA Smith MC Fast and energetic AGN-driven outflows in simulated dwarf galaxies Mon. Not. R. Astron. Soc. 2019 484 2047 2066 10.1093/mnras/stz097
43. Peng Y Maiolino R Cochrane R Strangulation as the primary mechanism for shutting down star formation in galaxies Nature 2015 521 192 195 10.1038/nature14439 25971510
44. Trussler J Both starvation and outflows drive galaxy quenching Mon. Not. R. Astron. Soc. 2020 491 5406 5434 10.1093/mnras/stz3286
45. Ceverino D Klessen RS Glover SCO FirstLight – II. Star formation rates of primeval galaxies from z = 5–15 Mon. Not. R. Astron. Soc. 2018 480 4842 4850
46. Kimm T Cen R Devriendt J Dubois Y Slyz A Towards simulating star formation in turbulent high-z galaxies with mechanical supernova feedback Mon. Not. R. Astron. Soc. 2015 451 2900 2921 10.1093/mnras/stv1211
47. Ma X Simulating galaxies in the reionization era with FIRE-2: galaxy scaling relations, stellar mass functions, and luminosity functions Mon. Not. R. Astron. Soc. 2018 478 1694 1715 10.1093/mnras/sty1024
48. Penny SJ SDSS-IV MaNGA: evidence of the importance of AGN feedback in low-mass galaxies Mon. Not. R. Astron. Soc. 2018 476 979 998 10.1093/mnras/sty202
49. Jakobsen P The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope. I. Overview of the instrument and its capabilities Astron. Astrophys. 2022 661 A80 10.1051/0004-6361/202142663
50. Ferruit P The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope. II. Multi-object spectroscopy (MOS) Astron. Astrophys. 2022 661 A81 10.1051/0004-6361/202142673
51. Bunker, A. J. et al. JADES NIRSpec Initial Data Release for the Hubble Ultra Deep Field: redshifts and line fluxes of distant galaxies from the deepest JWST Cycle 1 NIRSpec Multi-Object Spectroscopy. Preprint at https://arxiv.org/abs/2306.02467 (2023).
52. Williams CC JEMS: a deep medium-band imaging survey in the Hubble Ultra Deep Field with JWST NIRCam and NIRISS Astrophys. J. Suppl. Ser. 2023 268 64 10.3847/1538-4365/acf130
53. Rieke, M. J., Kelly, D. & Horner, S. in Cryogenic Optical Systems and Instruments XI Vol. 5904 (eds Heaney, J. B. & Burriesci, L. G.) 590401 (SPIE, 2005).
54. Rieke MJ Performance of NIRCam on JWST in flight Publ. Astron. Soc. Pac. 2023 135 028001 10.1088/1538-3873/acac53
55. Robertson BE Identification and properties of intense star-forming galaxies at redshifts z > 10 Nat. Astron. 2023 7 611 621 10.1038/s41550-023-01921-1
56. Tacchella S JADES imaging of GN-z11: revealing the morphology and environment of a luminous galaxy 430 Myr after the Big Bang Astrophys. J. 2023 952 74 10.3847/1538-4357/acdbc6
57. Cappellari M Improving the full spectrum fitting method: accurate convolution with Gauss–Hermite functions Mon. Not. R. Astron. Soc. 2017 466 798 811 10.1093/mnras/stw3020
58. Cappellari, M. Full spectrum fitting with photometry in pPXF: stellar population versus dynamical masses, non-parametric star formation history and metallicity for 3200 LEGA-C galaxies at redshift z~0.8. Preprint at https://arxiv.org/abs/2208.14974 (2023).
59. Conroy C Resolving the metallicity distribution of the stellar halo with the H3 survey Astrophys. J. 2019 887 237 10.3847/1538-4357/ab5710
60. Choi J Mesa Isochrones and Stellar Tracks (MIST). I. Solar-scaled models Astrophys. J. 2016 823 102 10.3847/0004-637X/823/2/102
61. Calzetti D The dust content and opacity of actively star-forming galaxies Astrophys. J. 2000 533 682 695 10.1086/308692
62. Looser, T. J. et al. JADES: differing assembly histories of galaxies – observational evidence for bursty SFHs and (mini-)quenching in the first billion years of the Universe. Preprint at https://arxiv.org/abs/2306.02470 (2023).
63. Carnall AC McLure RJ Dunlop JS Davé R Inferring the star formation histories of massive quiescent galaxies with BAGPIPES: evidence for multiple quenching mechanisms Mon. Not. R. Astron. Soc. 2018 480 4379 4401 10.1093/mnras/sty2169
64. Witstok J Carbonaceous dust grains seen in the first billion years of cosmic time Nature 2023 621 267 270 10.1038/s41586-023-06413-w 37467786
65. Bruzual G Charlot S Stellar population synthesis at the resolution of 2003 Mon. Not. R. Astron. Soc. 2003 344 1000 1028 10.1046/j.1365-8711.2003.06897.x
66. Chevallard J Charlot S Modelling and interpreting spectral energy distributions of galaxies with BEAGLE Mon. Not. R. Astron. Soc. 2016 462 1415 1443 10.1093/mnras/stw1756
67. Sánchez-Blázquez P Medium-resolution Isaac Newton Telescope library of empirical spectra Mon. Not. R. Astron. Soc. 2006 371 703 718 10.1111/j.1365-2966.2006.10699.x
68. Bressan A PARSEC: stellar tracks and isochrones with the PAdova and TRieste Stellar Evolution Code Mon. Not. R. Astron. Soc. 2012 427 127 145 10.1111/j.1365-2966.2012.21948.x
69. Marigo P Bressan A Nanni A Girardi L Pumo ML Evolution of thermally pulsing asymptotic giant branch stars – I. The COLIBRI code Mon. Not. R. Astron. Soc. 2013 434 488 526 10.1093/mnras/stt1034
70. Ferland GJ The 2017 release Cloudy Rev. Mex. Astron. Astrofis. 2017 53 385 438
71. Charlot S Fall SM A simple model for the absorption of starlight by dust in galaxies Astrophys. J. 2000 539 718 731 10.1086/309250
72. Chevallard J Simulating and interpreting deep observations in the Hubble Ultra Deep Field with the JWST/NIRSpec low-resolution ‘prism’ Mon. Not. R. Astron. Soc. 2019 483 2621 2640 10.1093/mnras/sty2426
73. Carnall AC The VANDELS survey: the star-formation histories of massive quiescent galaxies at 1.0 < z < 1.3 Mon. Not. R. Astron. Soc. 2019 490 417 439 10.1093/mnras/stz2544
74. Carnall AC A massive quiescent galaxy at redshift 4.658 Nature 2023 619 716 719 10.1038/s41586-023-06158-6 37216978
75. Johnson BD Leja J Conroy C Speagle JS Stellar population inference with Prospector Astrophys. J. Suppl. Ser. 2021 254 22 10.3847/1538-4365/abef67
76. Tacchella S Fast, slow, early, late: quenching massive galaxies at z ~ 0.8 Astrophys. J. 2022 926 134 10.3847/1538-4357/ac449b
77. Leja J Carnall AC Johnson BD Conroy C Speagle JS How to measure galaxy star formation histories. II. Nonparametric models Astrophys. J. 2019 876 3 10.3847/1538-4357/ab133c
78. Trebitsch M Blaizot J Rosdahl J Devriendt J Slyz A Fluctuating feedback-regulated escape fraction of ionizing radiation in low-mass, high-redshift galaxies Mon. Not. R. Astron. Soc. 2017 470 224 239 10.1093/mnras/stx1060
79. Jeffreys, H. Theory of Probability (Clarendon Press, 1939).
80. Worthey G Faber SM Gonzalez JJ Burstein D Old stellar populations. 5: absorption feature indices for the complete LICK/IDS sample of stars Astrophys. J. Suppl. Ser. 1994 94 687 722 10.1086/192087
81. Cardelli JA Clayton GC Mathis JS The relationship between infrared, optical, and ultraviolet extinction Astrophys. J. 1989 345 245 256 10.1086/167900
82. Reddy NA The HDUV Survey: a revised assessment of the relationship between UV slope and dust attenuation for high-redshift galaxies Astrophys. J. 2018 853 56 10.3847/1538-4357/aaa3e7
83. Shapley AE Sanders RL Redd NA Topping MW Bramme GB JWST/NIRSpec Balmer-line measurements of star formation and dust attenuation at z ~ 3–6 Astrophys. J. 2023 954 157 10.3847/1538-4357/acea5a
84. Zahid HJ Kudritzki R-P Conroy C Andrews B Ho IT Stellar absorption line analysis of local star-forming galaxies: the relation between stellar mass, metallicity, dust attenuation, and star formation rate Astrophys. J. 2017 847 18 10.3847/1538-4357/aa88ae
85. Planck Collaboration Planck 2018 results. VI. Cosmological parameters Astron. Astrophys. 2020 641 A6 10.1051/0004-6361/201833910
86. Pacifici C The evolution of star formation histories of quiescent galaxies Astrophys. J. 2016 832 79 10.3847/0004-637X/832/1/79
