# Photometric detection at 7.7 µm of a galaxy beyond redshift 14 with JWST/MIRI

Jakob M. Helton, George H. Rieke, Stacey Alberts, Zihao Wu, Daniel J.

Eisenstein, Kevin N. Hainline, Stefano Carniani, Zhiyuan Ji, William M.

Baker, Rachana Bhatawdekar, et al.

► To cite this version:

Jakob M. Helton, George H. Rieke, Stacey Alberts, Zihao Wu, Daniel J. Eisenstein, et al. Photometric detection at 7.7 µm of a galaxy beyond redshift 14 with JWST/MIRI. Nature Astronomy, 2025, 9, pp.729-740. 10.1038/s41550-025-02503-z. insu-05101651

HAL Id: insu-05101651

https://insu.hal.science/insu-05101651v1

Submitted on 6 Jun 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Photometric detection at 7.7 µm of a galaxy beyond redshift 14 with JWST/MIRI

Received: 17 May 2024

Accepted: 6 February 2025

Published online: 7 March 2025

<figure><img src="image_2.png" /></figure>

Check for updates

Jakob M. Helton Ⓓ1, George H. Rieke¹, Stacey Alberts¹, Zihao Wu Ⓓ2,
Daniel J. Eisenstein Ⓓ2, Kevin N. Hainline¹, Stefano Carniani Ⓓ3, Zhiyuan Ji¹,
William M. Baker Ⓓ4,5, Rachana Bhatawdekar Ⓓ6, Andrew J. Bunker Ⓓ7,
Phillip A. Cargile², Stéphane Charlot Ⓓ8, Jacopo Chevallard Ⓓ7,
Francesco D'Eugenio Ⓓ4,5, Eiichi Egami Ⓓ1, Benjamin D. Johnson Ⓓ2,
Gareth C. Jones Ⓓ7, Jianwei Lyu Ⓓ1, Roberto Maiolino Ⓓ4,5,9,
Pablo G. Pérez-González Ⓓ10, Marcia J. Rieke Ⓓ1, Brant Robertson Ⓓ11,
Aayush Saxena Ⓓ7,9, Jan Scholtz⁴,⁵, Irene Shivaei¹⁰, Fengwu Sun¹,
Sandro Tacchella Ⓓ4,5, Lily Whitler Ⓓ1, Christina C. Williams Ⓓ12,
Christopher N. A. Willmer Ⓓ1, Chris Willott Ⓓ13, Joris Witstok⁴,⁵ & Yongda Zhu Ⓓ1

The James Webb Space Telescope (JWST) has spectroscopically confirmed numerous galaxies at $z > 10$. While weak rest-frame ultraviolet emission lines have only been seen in a handful of sources, the stronger rest-frame optical emission lines are highly diagnostic and accessible at mid-infrared wavelengths with the Mid-Infrared Instrument (MIRI) of JWST. We report the photometric detection of the distant spectroscopically confirmed galaxy JADES-GS-z14-0 at $z = 14.32_{-0.20}^{+0.08}$ with MIRI at 7.7 µm. The most plausible solution for the stellar-population properties is that this galaxy contains half a billion solar masses in stars with a strong burst of star formation in the most recent few million years. For this model, at least one-third of the flux at 7.7 µm originates from the rest-frame optical emission lines $H\beta$ and/or [O III]λλ4959, 5007. The inferred properties of JADES-GS-z14-0 suggest rapid mass assembly and metal enrichment during the earliest phases of galaxy formation. This work demonstrates the unique power of mid-infrared observations in understanding galaxies at the redshift frontier.

With the launch of the James Webb Space Telescope (JWST), extragalactic astronomy fundamentally changed. The Near Infrared Camera (NIRCam) moved the photometric redshift frontier from $z \approx 10$ to $z \approx 14-16$ (see, for example, refs. 1–6), while the Near Infrared Spectrograph (NIRSpec) pushed the spectroscopic redshift frontier from $z \approx 8$ to $z \approx 12–14$ (see, for example, refs. 7–10). Crucially, JWST discovered an early period of galaxy formation that was more vigorous than expected, with a sizable population of luminous galaxies and supermassive black holes less than a billion years after the Big Bang.

A companion paper reports the spectroscopic confirmation of JADES-GS-z14-0 at redshift $z = 14.32_{-0.20}^{+0.08}$, which makes it the most distant galaxy with a spectroscopically confirmed redshift¹¹. This galaxy

is remarkably luminous, with a rest-frame UV absolute magnitude $M_{UV} \approx -20.81 \pm 0.16$, which may require a reassessment of ideas about early galaxy formation, suggesting a slow decline in the number density of galaxies at $z > 12$, with increasing efficiency of galaxy formation in halos at higher redshifts³. The rest-frame UV continuum slope $\beta_{UV} \approx -2.20 \pm 0.07$ is relatively red for a very young stellar population, suggesting that the UV emission is affected by a small amount of dust attenuation. The full-width at half-maximum (FWHM) of the intrinsic rest-frame UV light profile is $0.16 \pm 0.01$ arcsec (corresponding to a deconvolved half-light radius of 260 $\pm$ 20 pc). Given its spatial extent, the rest-frame UV emission appears not to be dominated by emission from an active galactic nucleus. The properties of JADES-GS-z14-0 add
<figure><img src="image_3.png" /><figcaption>**Fig. 1 | A distant galaxy spectroscopically confirmed by JADES.** This galaxy was initially selected from ultradeep NIRCam and MIRI imaging with JWST (F770W–F277W–F115W shown as an RGB false-colour mosaic in the lower right). It was targeted for NIRSpec micro shutter array follow-up observations and has been spectroscopically confirmed at redshift $z = 14.32^{+0.08}_{-0.20}$. JADES-GS-z14-0 is to the right and the foreground galaxy NIRCam ID 183349 to the left. The apparent</figcaption></figure>

colour of JADES-GS-z14-0 is caused by the absorption of the NIRCam/F115W flux by the intervening IGM and the rest-frame optical nebular emission-line excess in MIRI/F770W relative to NIRCam/F277W. A scale bar of 0.6" is provided which corresponds to roughly 1.9 physical kpc (pkpc) at the observed redshift ($z=14.32$).

to the evidence that a population of luminous and massive galaxies was already in place less than 300 Myr after the Big Bang, with number densities more than ten times higher than extrapolations based on pre-JWST observations, as demonstrated in ref. 5.

The rest-frame optical nebular emission lines are one of the primary means to characterize the physical conditions in galaxies. However, the redshift of JADES-GS-z14-0 has moved these lines into the wavelength coverage of the Mid-Infrared Instrument (MIRI), beyond the wavelength coverage of NIRCam and NIRSpec. The superb performance of JWST¹², alongside the remarkable brightness of the most extreme high-redshift galaxies, will allow MIRI to explore the rest-frame optical regime and provide important insights into the nature of the earliest galaxies. MIRI has been used at the redshift frontier with the recent spectroscopic identification of the rest-frame optical emission lines [O III]λλ4959, 5007 and Hα in the galaxy GHZ2/GLASS-z12 at $z = 12.33 \pm 0.02$ (ref. 13). These results highlight the power of combining observations from NIRCam, NIRSpec and MIRI to understand the properties of the very first galaxies.

In this work, we present the robust photometric detection of JADES-GS-z14-0 with MIRI at an observed wavelength of $\lambda_{obs} \approx 7.7$ µm, corresponding to rest-frame wavelengths of $\lambda_{rest} \approx 4,400–5,700$ Å. These ultradeep observations with MIRI/F770W provide more information about the nature of this remarkable galaxy, and are among the deepest mid-infrared integrations to date, with an on-source integration time of $t_{obs} \approx 23.8$ h for JADES-GS-z14-0.

Figure 1 presents in the lower right the JWST F770W–F277W–F115W false-colour image for JADES-GS-z14-0. The apparent colour of this galaxy is caused by (1) the complete absorption of emission by the intergalactic medium (IGM) in the F115W filter and (2) the excess rest-frame

<table><caption>**Table 1 | Photometry for JADES-GS-z14-0**</caption><thead><tr><td>Instrument/Filter</td><td>Model fitting (nJy)</td></tr></thead><tbody><tr><td>NIRCam/F090W</td><td>-2.1 ± 0.6</td></tr><tr><td>NIRCam/F115W</td><td>-0.8 ± 0.4</td></tr><tr><td>NIRCam/F150W</td><td>1.2 ± 0.5</td></tr><tr><td>NIRCam/F162M</td><td>-1.5 ± 0.9</td></tr><tr><td>NIRCam/F182M</td><td>13.9 ± 0.4</td></tr><tr><td>NIRCam/F200W</td><td>34.8 ± 0.5</td></tr><tr><td>NIRCam/F210M</td><td>46.5 ± 0.6</td></tr><tr><td>NIRCam/F250M</td><td>47.2 ± 0.5</td></tr><tr><td>NIRCam/F277W</td><td>55.1 ± 0.6</td></tr><tr><td>NIRCam/F300M</td><td>49.8 ± 0.5</td></tr><tr><td>NIRCam/F335M</td><td>43.4 ± 0.5</td></tr><tr><td>NIRCam/F356W</td><td>47.3 ± 0.5</td></tr><tr><td>NIRCam/F410M</td><td>46.1 ± 0.8</td></tr><tr><td>NIRCam/F444W</td><td>46.9 ± 0.6</td></tr><tr><td>MIRI/F770W</td><td>74.4 ± 5.6</td></tr></tbody></table>

Columns: (1) instrument and filter combinations and (2) fiducial model-fitting photometry assuming an extended morphology.

optical emission in the F770W filter relative to the rest-frame UV emission in the F277W filter. JADES-GS-z14-0 is close in projection to a foreground galaxy at a separation of roughly 0.4 arcsec to the east, which we refer to as NIRCam ID 183349 (ref. 11). The lensing magnification caused
<figure><img src="image_4.png" /></figure>

**Fig. 2 | SED modelling for JADES-GS-z14-0.** Upper panel: the measured spectral flux density and corresponding uncertainties are used to constrain the various SED models with BAGPIPES. As provided in Table 1, the mean and $1\sigma$ s.d. of the fiducial model-fitting photometry are the blue circles and error bars. The median of the SED models is the grey line and unfilled squares, while the $1\sigma$ confidence interval is the shaded region. Middle panel: the median model photometry is

compared with the measured fluxes and uncertainties ($\chi^2$), while the total $\chi^2$ value is reported on the right. Lower panels: the transmission curves for the various filters. These results suggest that the excess flux in F770W relative to F444W is from nebular emission-line contributions, while the underlying continuum is relatively flat at rest-frame optical wavelengths.

by 183349 and another low-redshift neighbouring object (roughly 2.2 arcsec to the south) is estimated to be small, with a lensing magnification factor of $\mu = 1.2$ (ref. 11). All of the analyses and results presented here have been corrected for this magnification factor.

Interpreting the ultradeep MIRI observations in the context of the ultradeep NIRCam observations requires measurement of the flux density in MIRI relative to the flux density in one of the long-wavelength NIRCam filters. The FWHMs of the F770W and F444W point spread functions (PSFs) are 0.269 arcsec and 0.145 arcsec, respectively, which makes separating JADES-GS-z14-0 from 183349 challenging, but essential. To meet this challenge, we measure photometry by performing model fitting on the individual exposures before mosaicing, illustrated by Extended Data Fig. 1. JADES-GS-z14-0 and 183349 are fitted simultaneously in all of the NIRCam and MIRI exposures with ForcePho (B.D.J. et al., in preparation) and GALFIT$^{14,15}$. We model the sources with Sérsic light profiles. An image of the difference between the observations and our model shows only small residuals, demonstrating the validity of the modelling. We measure the F444W flux density to be $f_{F444W} = 46.9 \pm 0.6$ nJy and the F770W flux density to be $f_{F770W} = 74.4 \pm 5.6$ nJy for JADES-GS-z14-0. The measurements correspond to an excess flux of $\Delta f = 27.5 \pm 5.6$ nJy in F770W with respect to F444W. We report these measurements in Table 1. Simple PSF photometry confirms these results (see Extended Data Fig. 2 for detailed results). The quoted photometric uncertainties indicate the signal-to-noise ratios (S/N) but are probably underestimates of the true errors, since we do not account for

systematic uncertainties related to, for example, photometric calibration, background subtraction and/or parametric assumptions for the intrinsic light profiles of JADES-GS-z14-0 and 183349.

To interpret the source of the excess flux in F770W relative to F444W, we model the spectral energy distribution (SED) of JADES-GS-z14-0 using two Bayesian fitting codes: BAGPIPES (Bayesian Analysis of Galaxies for Physical Inference and Parameter Estimation)$^{16}$ and Prospector$^{17}$. Generally, we find two types of solution for the stellar-population properties. To explain the red rest-frame UV continuum slope, BAGPIPES prefers solutions with relatively young stellar populations and non-zero diffuse dust attenuation, while Prospector prefers solutions with older stellar populations but roughly zero diffuse dust attenuation. The Prospector models contain a strong Balmer break (Extended Data Fig. 3), predicting stellar masses that are nearly an order of magnitude larger than the BAGPIPES models (Extended Data Table 1), rivalling the maximum expected halo mass for galaxies at these redshifts. Additionally, the Prospector models have the bulk of their stars forming at $z \approx 18-20$ (corresponding to median mass-weighted stellar ages of $t_{\star} \approx 80-100$ Myr), with no recent star formation, which exacerbates concerns about the predicted stellar mass. The lack of recent star formation would be unexpected given the predicted burstiness of star formation in early galaxies (see, for example, ref. 18). Reference 19 shows that galaxies at the redshift frontier should have stellar masses, stellar ages and diffuse dust attenuations that are consistent with the predictions from the BAGPIPES models (see discussion in Methods).
<figure><img src="image_5.png" /><figcaption>**Fig. 3 | Stellar-population synthesis (SPS) modelling using the NIRCam and MIRI photometry.** The measured fluxes and uncertainties are used to constrain the various SED models with BAGPIPES. Shown are the posterior distributions of $M_*, Z_*, t_*, SFR_{10}, A_V$ and $EW_{[OIII]+H\beta}$. We report the 16th, 50th and 84th percentiles</figcaption></figure>

after combining the posterior distributions from the various SED models. Three different SFHs are assumed: parametric constant SFH, parametric delayed-tau SFH and non-parametric continuity SFH.

Furthermore, ref. 20 shows that the strength of the Balmer break is reduced substantially with a top-heavy initial mass function (IMF), as is likely for JADES-GS-z14-0 (see discussion below). For these reasons, we consider the solutions from BAGPIPES to be more likely than the solutions from Prospector. Although we are unable to reject them on a formal basis, the Prospector fits would have radical implications for models of early galaxy evolution.

Figure 2 shows the measured photometry in all of the available NIRCam and MIRI filters, alongside the inferred SED determined with BAGPIPES. This fitting of the SED is self-consistently able to constrain the properties of the (1) stellar populations, (2) dust attenuation and (3) nebular gas. For simplicity, SEDs are modelled assuming a Kroupa stellar IMF$^{21}$. Absorption from the IGM and attenuation from diffuse interstellar dust are included, along with nebular contributions from continuum and emission lines. To understand the effects of differing star-formation histories (SFHs), we ran models for (1) a parametric constant SFH, (2) a parametric delayed-tau SFH and (3) a non-parametric continuity SFH. Comparing the measured photometry and corresponding uncertainties (filled circles) with the median of the inferred photometry from the SED modelling (unfilled squares) demonstrates the success in the modelling. The $1\sigma$ confidence interval (shaded regions) suggests that the majority (that is, at least 50%, and up to 100%) of the excess flux in the F770W filter is from the nebular emission lines Hβ and [O III]λλ4959, 5007, while the underlying continuum in the rest-frame optical regime is flat and consistent with the measured flux density in F444W.

Figure 3 and Extended Data Table 2 present the marginalized distributions for the BAGPIPES constraints on stellar mass ($M_*/M_\odot$, upper left), stellar metallicity ($Z_*/Z_\odot$, upper middle), $t_*$ (Myr, upper right), star-formation rate averaged over the previous 10 million years ($SFR_{10}$ ($M_\odot$ yr$^{-1}$), lower left), diffuse dust attenuation as measured in the V band ($A_V$ (mag), lower middle) and rest-frame equivalent width of [O III] + Hβ ($EW_{[OIII]+H\beta}$, Å), lower right). The equivalent width is measured for the combined rest-frame optical nebular emission lines Hβ and [O III]λλ4959, 5007. Similarly, Extended Data Fig. 4 presents

constraints on the joint posterior distributions and SFHs with BAGPIPES. The median derived $M_* \approx 10^{8.7} M_\odot$ (with a $1\sigma$ confidence interval of $M_* \approx 10^{8.3}-10^{9.2} M_\odot$) is nearly one-tenth of the current value for the Milky Way. The median star-formation rate $SFR \approx 25 M_\odot$ yr$^{-1}$ (with a $1\sigma$ confidence interval of $SFR \approx 20-31 M_\odot$ yr$^{-1}$) is consistent with expectations based on the empirical star-forming main sequence derived at lower redshifts ($z \approx 8$; for example, ref. 22). Taken together with the measured half-light radius, the star-formation rates imply star-formation rate surface densities $\Sigma_{SFR} \approx 64 M_\odot$ yr$^{-1}$ pc$^{-2}$ (with a $1\sigma$ confidence interval of $\Sigma_{SFR} \approx 50-78 M_\odot$ yr$^{-1}$ pc$^{-2}$), comparable to the most vigorous starbursts observed in the local Universe (see, for example, ref. 23).

The inferred stellar masses and star-formation rates are deduced from the emission produced by massive stars and an assumed 'local' IMF. However, it is likely that the formation of low-mass stars is strongly suppressed due to both the high temperature of the cosmic microwave background ($T \approx 60$ K at $z=20$; ref. 24) and the low metallicity of galaxies at $z>10$ (refs. 25,26). Recent work has found that stellar masses can be reduced by factors of three (or more) by changing assumptions about the IMF, without affecting the resulting SED$^{27}$, which our own analysis confirms.

Returning to the BAGPIPES constraints on physical properties, $t_*$ ≈ 15 Myr (with a $1\sigma$ confidence interval of $t_* \approx 3-43$ Myr) is consistent with values measured for some of the youngest galaxies observed at lower redshifts ($z \approx 8$; see, for example, refs. 28,29). Similarly, median $EW_{[OIII]+H\beta} \approx 370$ Å (with a $1\sigma$ confidence interval of $EW_{[OIII]+H\beta} \approx 240-730$ Å), consistent with values measured for those same galaxies observed at lower redshifts ($z \approx 8$; see, for example, refs. 28,29). The median $A_V \approx 0.56$ AB mag (as measured in the V band, with a $1\sigma$ confidence interval of $A_V \approx 0.43-0.64$ AB mag) suggests a small amount of attenuation at rest-frame optical wavelengths. A detailed discussion and interpretation of the dust content for JADES-GS-z14-0 is presented in ref. 11. Finally, the median $Z_* \approx 0.014 Z_\odot$ (with a $1\sigma$ confidence interval of $Z_* \approx 0.002-0.066 Z_\odot$) is largely unconstrained but consistent with metallicities that are less than 10% of the solar value. These results agree with those determined in ref. 11, which derives a stellar metallicity from
SED modelling with BEAGLE³⁰. Assuming that the stellar and gas-phase metallicities are the same, they find $log_{10}(O/H) + 12 = 7.2^{+0.7}_{-0.4}$.

Our own SED fitting finds an equivalent gas-phase metallicity of $log_{10}(O/H) + 12 = 6.9^{+0.7}_{-0.8}$.

The gas-phase metallicity encodes valuable information about the baryonic processes shaping the formation and evolution of galaxies. It has been found to be correlated with the emission-line ratio $[O \text{ III}]/H\beta$, which compares the strengths of the collisionally excited $[O \text{ III}]\lambda\lambda4959, 5007$ lines with the Balmer recombination $H\beta$ line. We can place loose constraints on this emission-line ratio by combining the stellar-population properties (and the predicted rest-frame optical continuum) with the measured excess flux in F770W relative to F444W. We proceed analogously to the fundamental relations derived in ref. 31 and compare inferred star-formation rates with observed $H\beta$ line luminosities. The relationship between these two quantities is derived for a sample of galaxies at $z \approx 8$ with $H\beta$ line flux measurements from NIRSpec/PRISM observations, which are from JWST Advanced Deep Extragalactic Survey (JADES) Data Release 3 (DR3)³². Combining our derived calibration with the inferred star-formation rates from BAGPIPES yields an $H\beta$ line flux $F_{H\beta} = 7.9^{+1.7}_{-1.8} \times 10^{-19}$ ergs s⁻¹ cm⁻². Comparing this derived quantity with two assumptions about the underlying continuum at rest-frame optical wavelengths yields $[O \text{ III}]/H\beta = 2.5^{+0.9}_{-0.6}$ and $1.9^{+2.6}_{-0.7}$ (Extended Data Fig. 5). For comparison, the typical value is $[O \text{ III}]/H\beta \approx 6$ for the aforementioned sample of galaxies at $z \approx 8$ from JADES DR3³². These results imply that the measured excess flux at 7.7 µm includes a substantial contribution from $[O \text{ III}]\lambda\lambda4959, 5007$, but with a gas-phase metallicity that is smaller than typical values at $z \approx 8$ (see discussion in Methods).

The detection of JADES-GS-z14-0 at $z > 14$ by MIRI demonstrates its power in understanding the properties of the earliest galaxies. The most plausible solution for the flux observed at 7.7 µm with MIRI suggests substantial contributions from the nebular emission lines $[O \text{ III}]\lambda\lambda4959, 5007$, which indicates metal enrichment for this galaxy. An alternative model is possible, but less likely, since it suggests an extreme stellar mass and a strong Balmer break from evolved stars. Deep spectroscopic follow-up observations with MIRI's low-resolution spectrometer are required to disentangle these interpretations by directly measuring the contributions to this flux from the nebular emission lines $H\beta$ and $[O \text{ III}]\lambda\lambda4959, 5007$. Such observations would also include $H\alpha$, producing a direct measurement of the star-formation rate. Given the size, luminosity and redshift of JADES-GS-z14-0, these measurements would build on a truly unique opportunity to study galaxy formation when the Universe was less than 300 Myr old.

## Methods

Throughout this work, we report wavelengths in air and adopt the standard flat $\Lambda$CDM cosmology from Planck18 with $H_0 = 67.4$ km s⁻¹ Mpc⁻¹ and $\Omega_m = 0.315$ (ref. 33). All magnitudes are in the AB system³⁴. Uncertainties are quoted as 68% (1σ) confidence intervals, unless otherwise stated.

### Observations

The observations used in this work consist of infrared imaging with NIRCam and MIRI in the Great Observatories Origins Deep Survey South (GOODS-S)³⁵ field, near the Hubble Ultra Deep Field³⁶ and the JADES Origins Field³⁷. The NIRCam data were primarily observed as part of JADES³⁸, but also as part of the First Reionization Epoch Spectroscopic Complete Survey³⁹. The MIRI data were observed as part of JADES.

JADES-GS-z14-0 was observed with NIRCam in four separate programmes, which we separate into three categories on the basis of observing time. Ultradeep observations were conducted with programme IDs 1210 (N. Lützgendorf) and 3215 (D.J.E.) across an area of roughly 10 arcmin² using 14 photometric filters, including seven wide bands (F090W, F115W, F150W, F200W, F277W, F356W and F444W) and seven medium bands (F162M, F182M, F210M, F250M, F300M, F335M and F410M). These were observed on 20–24 October 2022 (for

programme ID 1210) and 16–24 October 2023 (for programme ID 3215), with a range in integration times of 39–73 h for each of the wide bands and 55–165 h for each of the medium bands, reaching 5σ depths of 1.9–3.0 nJy and 1.6–2.6 nJy, respectively. Medium-depth observations were conducted with programme ID 1180 (D.J.E.) across an area of roughly 40 arcmin² using eight filters, including the same seven wide bands as programme ID 1210 plus one medium band (F410M). These were observed on 29 September–5 October 2022 and 28 September–3 October 2023, with a range in integration times of 6–8 h for each of the filters, reaching 5σ depths of 4.1–7.3 nJy. Shallow observations were conducted with programme ID 1895 (P. Oesch) across an area of roughly 60 arcmin² using three filters (F182M, F210M and F444W). These were observed on 13–18 November 2022, with a range in integration times of 0.25–1.25 h, reaching 5σ depths of 10–12 nJy. Depths are estimated using 0.2–arcsec–radius circular apertures assuming point–source morphologies. The aforementioned NIRCam observations have been presented and discussed extensively in the literature (for example, ⁴,⁵,³⁷–⁴⁰).

JADES-GS-z14-0 was also observed with MIRI as a coordinated parallel to NIRCam with programme ID 1180 (D.J.E.). The MIRI imaging includes four pointings near the JADES Origins Field at a position angle of 300°, which together produce an ultradeep contiguous mosaic of roughly 8.8 arcmin². In all four of the pointings, we conducted two separate nine-point dithers of 1,361 second individual exposures with MIRI for five different NIRCam filter pairs (2 × 9 × 5 total exposures), for a total exposure time of roughly 61.3 ks with MIRI for each of the two dithers. We conducted two additional four-point dithers of 1,361 second individual exposures with MIRI for three different NIRCam filter pairs (2 × 4 × 3 total exposures), for a total exposure time of roughly 16.3 ks for each of the two dithers. The readout mode SLOWRI was utilized with 57 groups to minimize data volume. These data were obtained on 29 September–5 October 2022 and 28 September–3 October 2023, with a typical on-source integration time of roughly 43.1 h (155.2 ks), reaching a 5σ depth of 21 nJy (28.1 AB mag) using 0.4–arcsec–radius circular apertures and assuming point–source morphologies as before. At the location of JADES-GS-z14-0, which falls near the edge of the MIRI imaging, the typical on-source integration time is roughly 23.8 h (85.7 ks) which corresponds to a 5σ depth of 28 nJy (27.8 AB mag). These MIRI observations have been presented and discussed previously at half depth in the literature (for example, ref. 41). This work is a presentation and discussion of the full-depth MIRI observations, which are some of the deepest images ever taken in the mid-infrared, alongside the ultradeep MIRI Deep Imaging Survey results in the Hubble Ultra Deep Field at 5.6 µm (refs. 42,43). The observations presented here will be described in more detail in a forthcoming paper from the JADES collaboration (S.A. et al., manuscript in preparation).

### Image reduction

A detailed description of the reduction, mosaicing, source detection and photometric measurements for the NIRCam data is provided in the first JADES data release in GOODS-S⁴⁰ and in ref. 11. Similarly, a detailed description of the reduction, mosaicing, source detection and photometric measurements for the MIRI data will be provided in an upcoming JADES data release in GOODS-S (S.A. et al., in preparation), but has already been partially introduced⁴¹,⁴⁴,⁴⁵. We briefly summarize the main steps of the reduction and mosaicing process for the MIRI data. Raw images are processed with the JWST Calibration Pipeline (v.1.12.5)⁴⁶ using the Calibration Reference Data System pipeline mapping 1188, similar to the procedure with NIRCam data. We run stage 1 of the JWST Calibration Pipeline using all of the default parameters, plus a correction for cosmic ray showers. Stage 2 is run using all of the default parameters, plus an additional custom subtraction of the background using our own super sky backgrounds. Following this step, we perform an additional astrometric correction. Finally, stage 3 is run using all of the default parameters, but without any further alignment or matching. The final image mosaic is registered to the Gaia DR3 frame⁴⁷ and
resampled onto the same world coordinate system as the NIRCam
image mosaics, but with a 0.060 arcsec per pixel grid.

Detection and photometry

To interpret the ultradeep MIRI observations in the context of the
ultradeep NIRCam observations, we measure the flux density in MIRI/
F770W relative to the flux density in one of the long-wavelength NIRCam
filters (that is, F444W). The MIRI detection is of modest S/N ≈ 13 with
respect to all of the NIRCam detections (S/N ≈ 25–100). Furthermore, the
MIRI/F770W diffraction-limited PSF (FWHM of 0.269 arcsec) is larger
than that of NIRCam/F444W (FWHM of 0.145 arcsec) by nearly a factor
of two. JADES-GS-z14-0 and the neighbouring foreground galaxy to the
east, NIRCam ID 183349, have similar measured flux densities in F444W
and are both morphologically extended (FWHMs of the deconvolved
light profiles are both roughly 0.16 arcsec)¹¹. These qualities make the
separation (or deblending) of JADES-GS-z14-0 from 183349 challeng-
ing, but essential for physical interpretation of the MIRI observations.

**Detailed model-fitting photometry.** Our primary approach for measuring photometry uses detailed model fitting in the vicinity of JADES-GS-z14-0 with the Bayesian fitting code ForcePho (B.D.J. et al., manuscript in preparation). The motivation for performing these more complicated photometric measurements is the complexity of the region surrounding JADES-GS-z14-0, with multiple bright foreground galaxies within a radius of a few arcseconds. It is important to properly disentangle the relative flux contributions from these galaxies. Furthermore, this region of the sky is located near the edge of the MIRI exposures, which may affect measurements of the photometric uncertainties. We briefly summarize here the main steps in measuring the detailed model-fitting photometry. Since the galaxies are better resolved with NIRCam, we determine their morphological properties (that is, half-light radius, Sérsic index, axis ratio, position angle and total flux for each NIRCam filter) by fitting the NIRCam data, and then use the results to constrain the MIRI photometry.

The first step for the NIRCam data is to construct 2.5 arc-
sec × 2.5 arcsec cutouts from all of the available exposures, with each
of these cutouts centred on the source location of JADES-GS-z14-0 from
the original photometric catalogue. Cutouts are performed on expo-
sures after reduction but before mosaicing. All of the detected sources
within these cutouts are simultaneously modelled to determine the rel-
ative flux contributions. These models assume an intrinsic Sérsic profile
that is identical for all the different NIRCam filters⁴⁸. We convolve these
intrinsic light profiles with Gaussian mixture approximations to the PSF
in each of the NIRCam exposures to produce observed light profiles. We
utilize the results of this detailed model-fitting approach by taking the
mean and s.d. of the resulting total flux posterior distributions, which
are sampled with Markov chain Monte Carlo techniques and account
for the covariance between parameters.

JADES-GS-z14-0 and 183349 are significantly detected (S/N > 5) in the vast majority of the individual NIRCam exposures, but they are only marginally detected (S/N = 2–3) in most of the individual MIRI exposures. This necessitates different detailed model-fitting photometric procedures for NIRCam and MIRI. Extended Data Figure 1 illustrates the photometric modelling process used to measure the flux density in F770W. The first step is once again to construct 2.5 arcsec × 2.5 arcsec cutouts, now with all of the available MIRI exposures. All of the detected sources within these cutouts are simultaneously modelled by adopting the inferred morphological properties from the ForcePho fitting to the NIRCam data (for example, half-light radius, Sérsic index, axis ratio and position angle). For JADES-GS-z14-0, the inferred morphological parameters are half-light radius r<sub>1/2</sub> = 0.0788 ± 0.0006 arcsec, Sérsic index *n* = 0.877 ± 0.027, axis ratio *b/a* = 0.425 ± 0.008 and position angle PA = 1.427 ± 0.008 rad. For 183349, the inferred parameters are r<sub>1/2</sub> = 0.0872 ± 0.0006 arcsec, *n* = 0.815 ± 0.020, *b/a* = 0.378 ± 0.006 and PA = 4.171 ± 0.003 rad. The adopted Sérsic

profile is used to forward-model the total flux in the F770W filter with GALFIT<sup>14,15</sup>. We adopt the mean and s.d. of the resulting total flux posterior distributions as fiducial for MIRI. For JADES-GS-z14-0, we measure $f_{\text{F444W}} = 46.9 \pm 0.6 \text{ nJy}$ and $f_{\text{F770W}} = 74.4 \pm 5.6 \text{ nJy}$, corresponding to $\Delta f = 27.5 \pm 5.6 \text{ nJy}$ in F770W with respect to F444W. We report these measurements in Table 1. The uncertainties are additionally validated using two different methodologies, producing consistent results within 10% of one another: (1) bootstrapping individual exposures and (2) measuring photometry in regions of the sky that are observed to be empty in the much deeper NIRCam images. For 183349, we measure $f_{\text{F444W}} = 48.9 \pm 0.6 \text{ nJy}$ and $f_{\text{F770W}} = 46.3 \pm 4.6 \text{ nJy}$.

The model-fitting approach assumes that JADES-GS-z14-0 and 183349 are well described by single Sérsic light profiles. To explore the impact of this assumption, we additionally adopt a non-parametric approach. We construct galaxy templates by deconvolving the NIR-Cam images in the F150W and F444W filters, where the F150W image corresponds to 183349 alone, while the F444W image corresponds to both JADES-GS-z14-0 and 183349. We conduct the deconvolution using the Wiener-Hunt method and Gaussian mixture approxima-tions for the PSFs to mitigate noise amplification at high frequencies. We then convolve the templates into the MIRI/F770W band and fit the amplitudes of the two templates simultaneously, where we measure on individual exposures and analytically derive the best estimated values and covariances as before. We assign zero weight to regions below a significance of 3σ to prevent overfitting. Finally, we convert the template amplitude to galaxy fluxes according to the ForcePho pho-tometry results in the NIRCam F150W and F444W bands. This method yields $f_{\text{F770W}} = 69.3 \pm 6.1 \text{ nJy}$ for JADES-GS-z14-0, which is consistent with the parametric model-fitting photometry. Moreover, we apply this non-parametric approach to measure fluxes across all of the NIRCam filters, obtaining results consistent with the reported photometry in Table 1. For 183349, this method yields $f_{\text{F770W}} = 47.8 \pm 4.0 \text{ nJy}$.

**PSF photometry.** Our primary measurements are based on the model-fitting approach described above. Alternatively, we could treat both galaxies as point sources, since their measured half-light angular sizes are notably smaller than the F770W PSF. Therefore, as a check on the model-fitting result, we also perform PSF photometry. To measure the observed flux density in F770W relative to F444W, we first convolve the F444W mosaic to the PSF of the F770W filter, which involves convolving the F444W mosaic with a kernel of the difference between the F444W and F770W PSFs. The adopted MIRI PSF is empirically measured from the final image mosaic and accounts for the ‘cruciform’ detector artefact⁴⁹. We simultaneously fit for the flux densities of JADES-GS-z14-0 and 183349 in both the convolved F444W and unconvolved F770W image mosaics. Extended Data Figure 2 illustrates the photometric modelling process used to measure the flux density in F770W. The residual images (shown in the right-hand panel) illustrate the validity of the point-source assumption and the success in the photometric modelling. For JADES-GS-z14-0, we measure $f_{\text{F444W}} = 46.4 \pm 1.2 \text{ nJy}$ and $f_{\text{F770W}} = 64.7 \pm 6.1 \text{ nJy}$, corresponding to $\Delta f = 18.3 \pm 6.2 \text{ nJy}$ in F770W with respect to F444W. The uncertainty on the F770W flux density is found by injecting fake point sources into the final image mosaic, extracting them with procedures identical to those used for JADES-GS-z14-0, and then calculating the sigma-clipped s.d. of the difference between the recovered and injected flux densities. These results are consistent with those obtained using both the parametric and non-parametric model-fitting approaches. The measured F770W flux density is slightly smaller for the PSF photometry, consistent with a small amount (roughly 10%) of galaxy light falling far enough outside the image cores that it is not included in this photometric approach.

**Conclusion.** To summarize, we obtain consistent results for the NIR-Cam and MIRI photometry of JADES-GS-z14-0 using three different methods. The PSF fitting approach provides a simple baseline that is
free from complex modelling assumptions. The model-fitting approach adopts both parametric and non-parametric galaxy models while measuring photometry from the individual exposures. This method more accurately accounts for the extended morphologies of JADES-GS-z14-0 and 183349, which improves measurements of the diffuse galaxy flux. This extended flux is the source of the larger F770W flux densities that are measured with the model-fitting approach when compared with the simpler PSF fitting approach.

## SPS modelling

The combination of ultradeep NIRCam, NIRSpec and MIRI observations provides an unprecedented opportunity to study the physical properties of JADES-GS-z14-0. We utilize the Bayesian SED fitting code BAGPIPES¹⁶ to self-consistently model the properties of the stellar populations, dust attenuation and nebular gas. We choose to sample the posterior distributions of these derived properties with the importance nested sampling code nautilus⁵⁰, assuming an effective sample size of 10⁴. Fits are performed on the fiducial model-fitting photometry after imposing an error floor of 5%. Such an error floor is imposed since the quoted photometric uncertainties are probably underestimates of the true errors. This is because we do not account for systematic uncertainties related to, for example, imperfect photometric calibration, background subtraction and/or parametric assumptions about the intrinsic light profiles (where we assume a single Sérsic profile for each source). We briefly summarize here the various components of the assumed physical model.

Stellar populations are derived using predefined SPS models³⁰, which are the 2016 updated version of previous models⁵¹, and are determined for a grid of simple stellar-population models with various ages and metallicities. We adopt the stellar library from the Medium-Resolution Isaac Newton Telescope Library of Empirical Spectra⁵², in addition to the stellar evolutionary tracks and isochrones from the Padova and Trieste Stellar Evolution Code⁵³. We assume a Kroupa IMF²¹ with a lower bound of 0.08 $M_\odot$ and an upper bound of 120 $M_\odot$. Absorption from the IGM is modelled after ref. 54, which is an updated version of the model from ref. 55. Dust attenuation is modelled after ref. 56 with one free parameter: $A_\nu$, assuming a uniform prior with min = 0.0 and max = 2.0. Nebular emission (from both emission lines and continuum) is modelled after ref. 57 using the photoionization code Cloudy⁵⁸ with one free parameter: the ionization parameter ($\log_{10}(U)$, assuming a uniform prior with min = -4.0 and max = -2.0). The gas-phase metallicity is fixed to the value of the stellar metallicity.

To explore the impact of assuming different SFHs on the results, we model the SEDs with two parametric SFHs (constant model and delayed-tau model) and one non-parametric SFH (continuity model). Each of the assumed SFHs has at least two free parameters: the total stellar mass formed ($\log_{10}(M_*/M_\odot)$, assuming a uniform prior with min = 6.0 and max = 12.0) and the stellar metallicity ($\log_{10}(Z_*/Z_\odot)$, assuming a uniform prior with min = -3.0 and max = 0.0). The constant-SFH model has one additional free parameter: the galaxy age ($t$ (Myr), assuming a uniform prior with min = 1.0 and max = $t_{\text{univ}}$ (Myr), where $t_{\text{univ}}$ is the age of the Universe measured with respect to the formation redshift $z_{\text{form}} = 20$). The delayed-tau SFH model has two additional free parameters: $t$ (Myr), assuming a uniform prior with min = 1.0 and max = $t_{\text{univ}}$ (Myr), and the e-folding time for the delayed-tau component ($\tau$ (Gyr), assuming a log-uniform prior with min = 0.001 and max = 30.0). Finally, the continuity SFH model has four additional free parameters, corresponding to the logarithm of the ratio of the star-formation rates in the five adjacent time-bins ($R_{\text{SFR}}$, assuming Student's t-distribution prior with $\mu = 0.0, \sigma = 0.3$). These time-bins are spaced at lookback times of 0–3, 3–10, 10–30, 30–100 and 100–$t_{\text{univ}}$ Myr, which assumes that the SFH starts at $z_{\text{form}} = 20$. These physical models have between five and eight free parameters, which should be compared with the 11 photometric detections that we have for JADES-GS-z14-0 (F182M, F200W, F210M, F250M, F277W, F300M, F335M, F356W, F410M, F444W and F770W).

Figure 3 presents the marginalized distributions for the BAGPIPES constraints on $M_*, Z_*, t_*$, $SFR_{10}$, $A_\nu$ and $EW_{[OIII]+H\beta}$. Results from each of the assumed SFHs are shown: the parametric constant-SFH model is represented by the purple solid lines, the parametric delayed-tau SFH model by the green dashed lines and the non-parametric continuity SFH model by the orange dotted lines. Throughout, we quote median derived properties and $1\sigma$ confidence intervals by concatenating the results from each of these SFHs. These measurements are reported in Extended Data Table 1, alongside measurements of the star-formation rate surface density ($\Sigma_{\text{SFR}_{10}}$), specific star-formation rate ($\text{sSFR}_{10}$) and $\log_{10}(U)$.

Similar to Fig. 3, Extended Data Fig. 4 shows the joint posterior distributions for some of the inferred physical parameters in the lower left, alongside the derived SFHs in the upper right. Results from each of the assumed SFHs are provided: the parametric constant-SFH model is shown in purple, the parametric delayed-tau SFH model in green and the non-parametric continuity SFH model in orange. The inferred stellar masses have two peaks in their posterior distributions, one at low $M_* \approx 10^{8.5} M_\odot$ and one at high $M_* \approx 10^{9.0} M_\odot$. The low-mass solutions suggest less stellar continuum at rest-frame optical wavelengths, and therefore larger equivalent widths of the rest-frame optical nebular emission lines $H\beta$ and [O III]$\lambda$4959, 5007. The opposite is true for the high-mass solutions, where the equivalent widths are smaller due to more stellar continuum in the rest-frame optical wavelengths. The delayed-tau and continuity SFHs prefer the low-mass solution, while the constant SFH finds equal weight for the two solutions. For all of the models, stellar mass is degenerate with stellar age, where the low-mass solution corresponds to younger stellar populations (with mass-weighted ages of a few million years), while the high-mass solution corresponds to older stellar populations (with ages of a few tens of millions of years). Some of the models suggest more extended periods of star formation out to lookback times of up to 100 Myr, which would be enough time to enrich this galaxy via type II supernovae. Despite these differences, the rest of the inferred physical parameters are similar for the two stellar-mass and stellar-age solutions. All of the inferred physical parameters are fairly well constrained by the existing observations, except for the stellar metallicity (and therefore also the gas-phase metallicity). The results of our SED modelling agree within the quoted uncertainties with the SED modelling of ref. 11, which includes JWST/NIRSpec data and adopts various physical models with BEAGLE³⁰.

To explore the impact of adopting different fitting codes, we additionally utilize the Bayesian SED fitting code Prospector (v.1.2.0)¹⁷ to self-consistently model the properties of the stellar populations, dust attenuation and nebular gas. We choose to sample the posterior distributions of these derived properties with the dynamic nested sampling code dynasty(v.1.2.3)⁵⁹, assuming an effective sample size of 10⁴. Fits are performed on the model-fitting photometry after imposing an error floor of 5%. The assumed physical model closely follows that of the BAGPIPES modelling, and we briefly summarize here the various components of this model.

Stellar populations are derived with the Flexible Stellar Population Synthesis code (FSPS)⁶⁰,⁶¹, which is accessed through the python-fsps bindings⁶². Stellar evolution is computed by the Modules for Experiments in Stellar Astrophysics package (MESA)⁶³⁻⁶⁶, while using the synthetic models from MESA Isochrones and Stellar Tracks⁶⁷,⁶⁸. The stellar libraries, in addition to the stellar evolutionary tracks and isochrones, are the primary differences between the physical models that we assumed for BAGPIPES and Prospector. Thus, we attribute any differences in the derived physical properties to these differences in the stellar libraries, stellar evolutionary tracks and isochrones. We assume a Kroupa IMF²¹ with a lower bound of 0.08 $M_\odot$ and an upper bound of 120 $M_\odot$. Absorption from the IGM is modelled after ref. 55. Dust attenuation is modelled after ref. 56 with one free parameter: $A_\nu$, assuming a uniform prior with min = 0.0 and max = 2.0. Nebular emission (from both emission lines and continuum) is modelled after ref. 57
using the photoionization code Cloudy⁵⁸ with one free parameter: log₁₀(U), assuming a uniform prior with min = -4.0 and max = -1.0. The gas-phase metallicity is fixed to the value of the stellar metallicity.

Finally, we assume the continuity model for the SFH. This is the only SFH that we assume for Prospector, since the two parametric SFHs (constant model and delayed-tau model) are unable to reproduce the observed photometry without invoking additional free parameters (for example, the escape fraction, or the shape of the diffuse dust attenuation curve). The inclusion of additional free parameters, such as the escape fraction, reduces the inferred stellar masses by roughly 0.5 dex without affecting any of the other inferred physical parameters. This SFH has two free parameters, corresponding to the total stellar mass formed (log₁₀(M<sub>⋆</sub>/M<sub>☉</sub>)), assuming a uniform prior with min = 6.0 and max = 12.0), and the stellar metallicity (log₁₀(Z<sub>⋆</sub>/Z<sub>☉</sub>), assuming a uniform prior with min = -3.0 and max = 0.0). It also has four additional free parameters, corresponding to the logarithm of R<sub>SFR</sub>, assuming Student's t-distribution prior with μ = 0.0, σ = 0.3. These time-bins are spaced at lookback times of 0-3, 3-10, 10-30, 30-100 and 100-t<sub>unv</sub> Myr, which assumes that the SFH starts at z<sub>form</sub> = 20. This physical model has eight free parameters, which should be compared with the 11 photometric detections that we have for JADES-GS-z14-0.

Similar to Fig. 2, Extended Data Fig. 3 shows the SED modelling for JADES-GS-z14-0 with Prospector. In the upper panel, the fiducial model-fitting photometry is used to constrain the various SED models. The median of these models is the grey line and unfilled squares, while the 1σ confidence interval is the shaded region. Extended Data Table 2 reports the 16th, 50th and 84th percentiles for the inferred physical properties. In the lower panel, the median model photometry is compared with the measured fluxes and uncertainties (χ), while the total χ² value is reported on the right. Prospector suggests that the excess flux in F770W relative to F444W is from stellar continuum rather than nebular emission-line contributions, producing strong Balmer breaks in the models, which we find unlikely for the following reasons.

1. Prospector predicts stellar masses of $M_{\star} \approx 10^{9.4} M_{\odot}$ that are nearly an order of magnitude larger than for the BAGPIPES models ($M_{\star} \approx 10^{8.7} M_{\odot}$). Reference 5 demonstrated that halos with virial masses of $M_{\star} \approx 10^{9.8-9.9} M_{\odot}$ have abundances comparable to those of nine galaxy candidates at $z=12-15$, including JADES-GS-z14-0. This implies that the Prospector-predicted stellar mass is only three times smaller than the halo mass predicted through abundance matching, which is smaller than the predicted stellar-to-halo-mass relation at $z \approx 14$ (see, for example, ref. 69).

2. Theoretical predictions from the First Light and Reionization Epoch Simulations (FLARES)¹⁹ suggest that galaxies with $M_{\star} \ge 10^{9.0} M_{\odot}$ only exist at $z < 14$, while galaxies with $M_{\star} \ge 10^{9.5} M_{\odot}$ only exist at $z < 13$. Additionally, theoretical predictions from the IllustrisTNG and THESAN projects⁷⁰ suggest that galaxies with $M_{\star} \ge 10^{9.0} M_{\odot}$ only exist at $z < 12$. Thus, cosmological galaxy formation simulations do not predict any galaxies at $z \approx 14$ with stellar masses that are comparable to those predicted by the Prospector models.

3. The issues related to the extreme stellar masses predicted by Prospector are exacerbated by the extreme stellar ages ($t_{\star} \approx 80-100$ Myr, corresponding to formation redshifts of $z \approx 18-20$). Given that cosmological simulations are unable to produce galaxies at $z \approx 14$ with stellar masses that are comparable to those predicted by the Prospector models, they certainly cannot produce those same galaxies at $z \approx 18-20$. On the other hand, the stellar masses, stellar ages, specific star-formation rates and diffuse dust attenuations predicted by the BAGPIPES models are consistent with theoretical predictions from FLARES¹⁹.

If they were correct, the physical properties predicted by the Prospector models would have radical implications for models of galaxy evolution

in the early Universe, since cosmological simulations do not predict any galaxies at the redshift frontier with these extreme inferred stellar masses and stellar ages. For the reasons outlined above, we consider the solutions from BAGPIPES to be more plausible than the solutions from Prospector, but we are unable to reject either of these solutions on a formal basis. As a reminder, we attribute any differences in the predicted physical properties between BAGPIPES and Prospector to differences in the assumed stellar libraries, in addition to the stellar evolutionary tracks and isochrones.

Bayesian SED fitting requires numerous assumptions to interpret observations and infer physical properties for galaxies. We have explored the impact of adopting different SFHs and fitting codes, where the former had very little impact on the derived stellar-population properties while the latter had substantial impact. Another assumption that can have substantial impact on the derived stellar-population properties is the IMF. In particular, it is likely that the formation of low-mass stars (that is, those with $M_{\star} < 1-3 M_{\odot}$) is strongly suppressed at high redshifts. This is caused by the high temperature of the cosmic microwave background²⁴ and the low metallicity of galaxies at $z > 10$ (refs. 25,26). To explore the impact of these physical processes, we utilize the Bayesian SED fitting code Prospector, since BAGPIPES does not allow changes to the assumed IMF. With regard to top-heavy IMFs, Prospector provides the more demanding test because it fits the stellar continuum with older, low-mass stars whose formation is likely to be suppressed. For simplicity, we will assume the same Kroupa IMF²¹ as before, but vary the lower and upper bounds of the mass range. Increasing the lower bound from 0.08 $M_{\odot}$ to 1 $M_{\odot}$ (3 $M_{\odot}$) decreases the inferred stellar masses by roughly 0.3 dex (0.4 dex), while increasing the upper bound from 120 $M_{\odot}$ to 300 $M_{\odot}$ has no effect on the inferred stellar masses. These results suggest that the reported stellar masses based on a local IMF can be overestimated by up to a factor of three without affecting the resulting SED. Recent work has found similar conclusions by changing assumptions about the IMF²⁷.

## Nebular emission-line flux predictions

The aforementioned stellar-population and dust properties provide predictions for the strength of the Balmer hydrogen recombination lines Hβ and Hα. With an estimate of the strength of Hβ, it is possible to determine the remaining excess flux in F770W relative to F444W, which can be attributed to the metallic collisionally excited lines [O III]λλ44959, 5007. These predictions require a series of steps and assumptions to convert the inferred star-formation rates into emission-line flux predictions then broadband flux predictions.

The first step in this process involves converting the BAGPIPES-derived star-formation rate into an Hβ line luminosity. To accomplish this, we proceed analogously to the fundamental relations derived in ref. 31. We compare BAGPIPES-derived star-formation rates with observed Hβ line luminosities for a sample of N = 27 galaxies at $z \approx 8$ with Hβ line flux measurements from NIRSpec/PRISM. These measurements, along with Kron photometry convolved to the F444W PSF, are from JADES DR3³². We select all sources in GOODS-S with NIRCam coverage and S/N > 3 for the Hβ line flux measurements. The measured photometry is used to constrain the various SED models with BAGPIPES, with procedures identical to those used for JADES-GS-z14-0, assuming the same three SFHs and concatenating the results from each. This produces a distribution of calibrations for each of the galaxies, which we can combine to obtain the full distribution of calibrations for the entire sample. They produce the following calibration:

$$C_{\mathrm{H\beta}} = \frac{L_{\mathrm{H\beta}} (\mathrm{ergs s}^{-1})}{\mathrm{SFR}_{10} (M_{\odot} \mathrm{yr}^{-1})} = 9.1 \pm 5.6 \times 10^{40}. \quad (1)$$

The reported calibration and corresponding uncertainty arise from taking the sigma-clipped mean and s.d. of the full distribution of calibrations. The quoted uncertainty reflects the scatter around the
average relation for individual galaxies. The dominant term in this uncertainty originates from the derived star-formation rates. Since the derived star-formation rate uncertainties are propagated throughout, we do not propagate the calibration uncertainty for the remaining calculations. The reported calibration results in a predicted Hβ line luminosity $L_{Hβ} \approx 230_{-150}^{+150} \times 10^{40}$ ergs s⁻¹.

The second step in this process involves converting the Hβ line luminosity into an Hβ line flux. This requires a luminosity distance, $d_L = 4.803 \times 10^{29}$ cm at z = 14.32, which is dependent on the assumed cosmology. This results in a predicted Hβ line flux $F_{Hβ} \approx 7.9_{-1.8}^{+1.7} \times 10^{-19}$ ergs s⁻¹ cm⁻².

The third step in this process involves converting the Hβ line flux into the equivalent signal in the F770W band. This requires the effective bandwidth of the F770W filter, $W_{v,F770W} = 1.00 \times 10^{13}$ Hz, as measured in frequency space and reported in the JWST User Documentation for the MIRI Filters and Dispersers. This results in a predicted Hβ F770W flux density $f_{Hβ,F770W} \approx 7.9_{-1.8}^{+1.7}$ nJy.

The final step in this process involves comparing the Hβ F770W flux density with the total nebular emission-line contribution to F770W. These nebular contributions are equal to the F770W flux density minus the total continuum contribution to F770W. For the continuum contribution to F770W, we first assume that the continuum is flat at rest-frame optical wavelengths and consistent with the measured flux density in F444W. We measure $\Delta f = 27.5 \pm 5.6$ nJy in F770W relative to F444W, which results in a predicted line ratio $[O III]/Hβ \approx 2.5_{-0.6}^{+0.9}$, as illustrated in the left-hand panel of Extended Data Fig. 5. The quoted uncertainties on the line ratio include uncertainties on the star-formation rates, but not on the rest-frame optical continuum. To account for uncertainties in the rest-frame optical continuum, we further assume the posterior distribution of continuum levels from BAGPIPES. This results in a predicted line ratio $[O III]/Hβ \approx 1.9_{-0.7}^{+2.6}$, as illustrated in the right-hand panel of Extended Data Fig. 5. Regardless of the assumption about the rest-frame optical continuum, we find that the excess flux at 7.7 µm includes a substantial contribution from the $[O III]\lambda4959, 5007$ lines. This line ratio is smaller than typical values observed for galaxies at $z \approx 8$, as shown in Extended Data Fig. 5 with individual galaxies that have Hβ line flux measurements from JADES DR3³².

The line ratio $[O III]/Hβ$ is found to be correlated but degenerate with the gas-phase metallicity, with an additional, secondary dependence on the ionization and excitation states of the nebular gas. The turnover in the relation between $[O III]/Hβ$ and gas-phase metallicity occurs for $Z_* \approx 25\% Z_⊙$ (see, for example, ref. 71), where lower values correspond to lower stellar masses, and higher values correspond to higher stellar masses. Given the inferred stellar masses and metallicities from the SED fitting for JADES-GS-z14-0 ($M_* \le 10^9 M_⊙$ and $Z_* \le 10\% Z_⊙$), we assume the lower value of the gas-phase metallicity to break the double-valued degeneracy between $[O III]/Hβ$. The inferred stellar masses and metallicities for the individual galaxies that have Hβ line flux measurements from JADES DR3³² are similar to those derived for JADES-GS-z14-0, so we can make the same assumption for those galaxies. Thus, given that the line ratio $[O III]/Hβ$ for JADES-GS-z14-0 is smaller than typical values observed for galaxies at $z \approx 8$, these results suggest a smaller gas-phase metallicity for JADES-GS-z14-0 when compared with galaxies at $z \approx 8$.

## Data availability

The NIRCam data that support the findings of this study are publicly available at https://archive.stsci.edu/hlsp/jades. The MIRI data that support the findings of this study will be made available in a future release; advanced access may be granted on reasonable request to the corresponding author.

## Code availability

The AstroPy⁷²,⁷³ software suite is publicly available, as is BAGPIPES¹⁶, Cloudy⁵⁷, dynasty³⁹, ForcePho (B.D.J. et al., manuscript in preparation), FSPS⁶⁰,⁶¹, GALFIT¹⁴,¹⁵, Matplotlib⁷⁴, nautilus⁵⁰, NumPy⁷⁵,⁷⁶, Pandas⁷⁷, photutils⁷⁸, Prospector¹⁷, python-fsps⁶², SciPy⁷⁹, seaborn⁸⁰, Tiny Tim⁸¹ and WebbPSF⁸².

## References

1. Castellano, M. et al. Early results from GLASS-JWST. III. Galaxy candidates at z 9–15. Astrophys. J. Lett. 938, L15 (2022).

2. Finkelstein, S. L. et al. CEERS Key Paper. I. An early look into the first 500 Myr of galaxy formation with JWST. Astrophys. J. Lett. 946, L13 (2023).

3. Harikane, Y. et al. A comprehensive study of galaxies at z 9–16 found in the early JWST data: ultraviolet luminosity functions and cosmic star formation history at the pre-reionization epoch. Astrophys. J. Suppl. Ser. 265, 5 (2023).

4. Robertson, B. E. et al. Identification and properties of intense star-forming galaxies at redshifts $z>10$. Nat. Astron. 7, 611–621 (2023).

5. Robertson, B. et al. Earliest galaxies in the JADES Origins Field: luminosity function and cosmic star formation rate density 300 Myr after the Big Bang. Astrophys. J. 970, 31 (2024).

6. Hainline, K. N. et al. The cosmos in its infancy: JADES galaxy candidates at $z>8$ in GOODS-S and GOODS-N. Astrophys. J. 964, 71 (2024).

7. Curtis-Lake, E. et al. Spectroscopic confirmation of four metal-poor galaxies at $z=10.3–13.2$. Nat. Astron. 7, 622–632 (2023).

8. Fujimoto, S. et al. UNCOVER: a NIRSpec census of lensed galaxies at $z=8.50–13.08$ probing a high-AGN fraction and ionized bubbles in the shadow. Astrophys. J. 977, 250 (2024).

9. Wang, B. et al. UNCOVER: illuminating the early Universe—JWST/NIRSpec confirmation of $z>12$ galaxies. Astrophys. J. Lett. 957, L34 (2023).

10. Castellano, M. et al. JWST NIRSpec spectroscopy of the remarkable bright galaxy GHZ2/GLASS-z12 at redshift 12.34. Astrophys. J. 972, 143 (2024).

11. Carniani, S. et al. Spectroscopic confirmation of two luminous galaxies at a redshift of 14. Nature 633, 318–322 (2024).

12. Rigby, J. et al. The science performance of JWST as characterized in commissioning. Publ. Astron. Soc. Pac. 135, 048001 (2023).

13. Zavala, J. A. et al. A luminous and young galaxy at $z=12.33$ revealed by a JWST/MIRI detection of Hα and [O III]. Nat. Astron. 9, 155–164 (2025).

14. Peng, C. Y., Ho, L. C., Impey, C. D. & Rix, H.-W. Detailed structural decomposition of galaxy images. Astron. J. 124, 266–293 (2002).

15. Peng, C. Y., Ho, L. C., Impey, C. D. & Rix, H.-W. Detailed decomposition of galaxy images. II. Beyond axisymmetric models. Astron. J. 139, 2097–2129 (2010).

16. Carnall, A. C., McLure, R. J., Dunlop, J. S. & Davé, R. Inferring the star formation histories of massive quiescent galaxies with BAGPIPES: evidence for multiple quenching mechanisms. Mon. Not. R. Astron. Soc. 480, 4379–4401 (2018).

17. Johnson, B. D., Leja, J., Conroy, C. & Speagle, J. S. Stellar population inference with Prospector. Astrophys. J. Suppl. Ser. 254, 22 (2021).

18. Sun, G. et al. Bursty star formation naturally explains the abundance of bright galaxies at cosmic dawn. Astrophys. J. Lett. 955, L35 (2023).

19. Wilkins, S. M. et al. First Light and Reionization Epoch Simulations (FLARES) V: the redshift frontier. Mon. Not. R. Astron. Soc. 519, 3118–3128 (2023).

20. Wilkins, S. M. et al. First Light and Reionization Epoch Simulations (FLARES)—XIV. The Balmer/4000 Å breaks of distant galaxies. Mon. Not. R. Astron. Soc. 527, 7965–7973 (2024).

21. Kroupa, P. & Boily, C. M. On the mass function of star clusters. Mon. Not. R. Astron. Soc. 336, 1188–1194 (2002).
22. Popesso, P. et al. The main sequence of star-forming galaxies across cosmic times. *Mon. Not. R. Astron. Soc.* **519**, 1526–1544 (2023).

23. Genzel, R. et al. A study of the gas-star formation relation over cosmic time. *Mon. Not. R. Astron. Soc.* **407**, 2091–2108 (2010).

24. Chon, S., Ono, H., Omukai, K. & Schneider, R. Impact of the cosmic background radiation on the initial mass function of metal-poor stars. *Mon. Not. R. Astron. Soc.* **514**, 4639–4654 (2022).

25. Steinhardt, C. L., Kokorev, V., Rusakov, V., Garcia, E. & Sneppen, A. Templates for fitting photometry of ultra-high-redshift galaxies. *Astrophys. J. Lett.* **951**, L40 (2023).

26. Chon, S., Hosokawa, T., Omukai, K. & Schneider, R. Impact of radiative feedback on the initial mass function of metal-poor stars. *Mon. Not. R. Astron. Soc.* **530**, 2453–2474 (2024).

27. Woodrum, C. et al. Using JADES NIRCam photometry to investigate the dependence of stellar mass inferences on the IMF in the early Universe. *Proc. Natl Acad. Sci. USA* **121**, e2317375121 (2024).

28. Endsley, R. et al. 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.* **524**, 2312–2330 (2023).

29. 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. *Mon. Not. R. Astron. Soc.* **533**, 1111–1142 (2024).

30. Chevallard, J. & Charlot, S. Modelling and interpreting spectral energy distributions of galaxies with BEAGLE. *Mon. Not. R. Astron. Soc.* **462**, 1415–1443 (2016).

31. Kennicutt, J. & Robert, C. Star formation in galaxies along the Hubble sequence. *Annu. Rev. Astron. Astrophys.* **36**, 189–232 (1998).

32. D'Eugenio, F. et al. JADES Data Release 3—NIRSpec/MSA spectroscopy for 4,000 galaxies in the GOODS fields. Preprint at https://arxiv.org/abs/2404.06531 (2024).

33. Planck Collaboration Planck 2018 results. VI. Cosmological parameters. *Astron. Astrophys.* **641**, A6 (2020).

34. Oke, J. B. & Gunn, J. E. Secondary standard stars for absolute spectrophotometry. *Astrophys. J.* **266**, 713–717 (1983).

35. Giavalisco, M. et al. The Great Observatories Origins Deep Survey: initial results from optical and near-infrared imaging. *Astrophys. J. Lett.* **600**, L93–L98 (2004).

36. Beckwith, S. V. W. et al. The Hubble Ultra Deep Field. *Astron. J.* **132**, 1729–1755 (2006).

37. Eisenstein, D. J. et al. The JADES Origins Field: a new JWST deep field in the JADES second NIRCam data release. Preprint at https://arxiv.org/abs/2310.12340 (2023).

38. Eisenstein, D. J. et al. Overview of the JWST Advanced Deep Extragalactic Survey (JADES). Preprint at https://arxiv.org/abs/2306.02465 (2023).

39. Oesch, P. A. et al. The JWST FRESCO survey: legacy NIRCam/ grism spectroscopy and imaging in the two GOODS fields. *Mon. Not. R. Astron. Soc.* **525**, 2864–2874 (2023).

40. Rieke, M. J. et al. JADES initial data release for the Hubble Ultra Deep Field: revealing the faint infrared sky with deep JWST NIRCam imaging. *Astrophys. J. Suppl. Ser.* **269**, 16 (2023).

41. Alberts, S. et al. To high redshift and low mass: exploring the emergence of quenched galaxies and their environments at $3<z<6$ in the ultradeep JADES MIRI F770W parallel. *Astrophys. J.* **975**, 85 (2024).

42. Rinaldi, P., et al. MIDIS: strong (Hβ+[O iii]) and Hα emitters at redshift $z \simeq 7-8$ unveiled with JWST NIRCam and MIRI imaging in the Hubble eXtreme Deep Field. *Astrophys. J.* **952**, 143 (2023).

43. Boogaard, L. A. et al. MIDIS: JWST/MIRI reveals the stellar structure of ALMA-selected galaxies in the Hubble Ultra Deep Field at cosmic noon. *Astrophys. J.* **969**, 27 (2024).

44. Lyu, J. et al. Active galactic nuclei selection and demographics: a new age with JWST/MIRI. *Astrophys. J.* **966**, 229 (2024).

45. Pérez-González, P. G. et al. What is the nature of little red dots and what is not, MIRI SMILES edition. *Astrophys. J.* **968**, 4 (2024).

46. Bushouse, H. et al. JWST Calibration Pipeline. Zenodo https://doi.org/10.5281/zenodo.10022973 (2023).

47. Gaia Collaboration Gaia Data Release 3. Summary of the content and survey properties. *Astron. Astrophys.* **674**, A1 (2023).

48. Sersic, J. L. *Atlas de Galaxias Australes* (Observatorio Astronómico de Córdoba, 1968).

49. Gáspár, A. et al. The quantum efficiency and diffractive image artifacts of Si:As IBC mid-IR detector arrays at 5–10 µm: implications for the JWST/MIRI detectors. *Publ. Astron. Soc. Pac.* **133**, 014504 (2021).

50. Lange, J. U. nautilus: boosting Bayesian importance nested sampling with deep learning. *Mon. Not. R. Astron. Soc.* **525**, 3181–3194 (2023).

51. Bruzual, G. & Charlot, S. Stellar population synthesis at the resolution of 2003. *Mon. Not. R. Astron. Soc.* **344**, 1000–1028 (2003).

52. Falcón-Barroso, J. et al. An updated MILES stellar library and stellar population models. *Astron. Astrophys.* **532**, A95 (2011).

53. Bressan, A. et al. PARSEC: stellar tracks and isochrones with the PAdova and TRIeste Stellar Evolution Code. *Mon. Not. R. Astron. Soc.* **427**, 127–145 (2012).

54. Inoue, A. K., Shimizu, I., Iwata, I. & Tanaka, M. An updated analytic model for attenuation by the intergalactic medium. *Mon. Not. R. Astron. Soc.* **442**, 1805–1820 (2014).

55. Madau, P. Radiative transfer in a clumpy Universe: the colors of high-redshift galaxies. *Astrophys. J.* **441**, 18–27 (1995).

56. Calzetti, D. et al. The dust content and opacity of actively star-forming galaxies. *Astrophys. J.* **533**, 682–695 (2000).

57. Byler, N., Dalcanton, J. J., Conroy, C. & Johnson, B. D. Nebular continuum and line emission in stellar population synthesis models. *Astrophys. J.* **840**, 44 (2017).

58. Ferland, G. J. et al. The 2013 release of Cloudy. *Rev. Mex. Astron. Astrofis.* **49**, 137–163 (2013).

59. Speagle, J. S. dynasty: a dynamic nested sampling package for estimating Bayesian posteriors and evidences. *Mon. Not. R. Astron. Soc.* **493**, 3132–3158 (2020).

60. Conroy, C., Gunn, J. E. & White, M. The propagation of uncertainties in stellar population synthesis modeling. I. The relevance of uncertain apects of stellar evolution and the initial mass function to the derived physical properties of galaxies. *Astrophys. J.* **699**, 486–506 (2009).

61. Conroy, C. & Gunn, J. E. The propagation of uncertainties in stellar population synthesis modeling. III. Model calibration, comparison, and evaluation. *Astrophys. J.* **712**, 833–857 (2010).

62. Foreman-Mackey, D., Sick, J. & Johnson, B. python-fsps: Python bindings to FSPS (v0.1.1). Zenodo https://zenodo.org/ records/12157 (2014).

63. Paxton, B. et al. Modules for Experiments in Stellar Astrophysics (MESA): Astrophys. J. Suppl. Ser. **192**, 3 (2011).

64. Paxton, B. et al. Modules for Experiments in Stellar Astrophysics (MESA): planets, oscillations, rotation, and massive stars. *Astrophys. J. Suppl. Ser.* **208**, 4 (2013).

65. Paxton, B. et al. Modules for Experiments in Stellar Astrophysics (MESA): binaries, pulsations, and explosions. *Astrophys. J. Suppl. Ser.* **220**, 15 (2015).

66. Paxton, B. et al. Modules for Experiments in Stellar Astrophysics (MESA): convective boundaries, element diffusion, and massive star explosions. *Astrophys. J. Suppl. Ser.* **234**, 34 (2018).

67. Dotter, A. MESA Isochrones and Stellar Tracks (MIST) O: methods for the construction of stellar isochrones. *Astrophys. J. Suppl. Ser.* **222**, 8 (2016).
68. Choi, J. et al. MESA Isochrones and Stellar Tracks (MIST). I. Solar-scaled models. *Astrophys. J.* **823**, 102 (2016).

69. Behroozi, P. et al. The Universe at $z>10$: predictions for JWST from the UNIVERSEMACHINE DR1. *Mon. Not. R. Astron. Soc.* **499**, 5702–5718 (2020).

70. Kannan, R. et al. The MillenniumTNG project: the galaxy population at $z\ge8$. *Mon. Not. R. Astron. Soc.* **524**, 2594–2605 (2023).

71. Wilkins, S. M. et al. First Light and Reionization Epoch Simulations (FLARES) XI: [O III] emitting galaxies at $5<z<10$. *Mon. Not. R. Astron. Soc.* **522**, 4014–4027 (2023).

72. Astropy Collaboration Astropy: a community Python package for astronomy. *Astron. Astrophys.* **558**, A33 (2013).

73. Astropy Collaboration et al. The Astropy Project: building an open-science project and status of the v2.0 core package. *Astron. J.* **156**, 123 (2018).

74. Hunter, J. D. Matplotlib: a 2D graphics environment. *Comput. Sci. Eng.* **9**, 90–95 (2007).

75. van der Walt, S., Colbert, S. C. & Varoquaux, G. The NumPy array: a structure for efficient numerical computation. *Comput. Sci. Eng.* **13**, 22–30 (2011).

76. Harris, C. R. et al. Array programming with NumPy. *Nature* **585**, 357–362 (2020).

77. The Pandas Development Team pandas-dev/pandas: Pandas. *Zenodo* https://zenodo.org/records/7093122 (2022).

78. Bradley, L. et al. astropy/photutils: 1.5.0. *Zenodo* https://zenodo.org/records/6825092 (2022).

79. Virtanen, P. et al. SciPy 1.0: fundamental algorithms for scientific computing in Python. *Nat. Methods* **17**, 261–272 (2020).

80. Waskom, M. seaborn: statistical data visualization. *J. Open Source Softw.* **6**, 3021 (2021).

81. Krist, J. E., Hook, R. N. & Stoehr, F. 20 years of Hubble Space Telescope optical modeling using Tiny Tim. *Proc. SPIE* **8127**, 81270J (2011).

82. Perrin, M. D. et al. Updated point spread function simulations for JWST with WebbPSF. *Proc. SPIE* **9143**, 91433X (2014).

## Acknowledgements

J.M.H., G.H.R., S.A., Z.H., D.J.E., P.A.C., E.E., B.D.J., M.J.R., B.R., F.S. and C.N.A.W. are supported by NASA contracts NAS5-02105 and NNX13AD82G to the University of Arizona. D.J.E. is supported as a Simons Investigator. S. Carniani acknowledges support by the European Union's HE ERC Starting Grant 101040227-WINGS. W.M.B. and J.S. acknowledge support by the Science and Technology Facilities Council (STFC), European Research Council (ERC) Advanced Grant 695671 'QUENCH'. A.J.B., J.C., G.C.J. and A.S. acknowledge funding from the 'FirstGalaxies' Advanced Grant from the ERC under the European Union's Horizon 2020 research and innovation programme (grant agreement 789056). F.D'E., R.M. and J.W. acknowledges support by the Science and Technology Facilities Council (STFC), by the ERC through Advanced Grant 695671 'QUENCH' and by the UKRI Frontier Research grant RISEandFALL. R.M. also acknowledges funding from a research professorship from the Royal Society. P.G.P.-G. acknowledges support from grant PID2022-139567NB-I00 funded by the Spanish Ministerio de Ciencia e Innovación MCIN/AEI/10.13039/501100011033, FEDER, UE. S.T. acknowledges support by the Royal Society Research Grant G125142. L.W. acknowledges support from the National Science Foundation Graduate Research Fellowship under grant DGE-2137419. 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 work made use of the lux supercomputer at UC Santa Cruz, which is funded by NSF MRI grant AST 1828315, as well as the High Performance Computing (HPC) resources at the University of Arizona, which are funded by the Office of Research Discovery and Innovation

(ORDI), Chief Information Officer (CIO) and University Information Technology Services (UITS). We respectfully acknowledge that the University of Arizona is on the land and territories of Indigenous peoples. At the time of writing, Arizona is home to 22 federally recognized tribes, with Tucson being home to the O'odham and the Yaqui. Committed to diversity and inclusion, the university strives to build sustainable relationships with sovereign Native Nations and Indigenous communities through education offerings, partnerships and community service.

## Author contributions

J.M.H. and G.H.R. led the writing of this paper. J.M.H., G.H.R., S.A., Z.W., D.J.E., Z.J., J.L., P.G.P.-G. and I.S. contributed to the MIRI imaging reduction and the modelling of the MIRI photometry. D.J.E., P.A.C., B.D.J., M.J.R., B.R., S.T. and C.N.A.W. contributed to the NIRCam imaging reduction and the modelling of the NIRCam photometry. S. Carniani, A.J.B., S. Charlot, J.C., F.D'E., R.M. and J.W. contributed to the NIRSpec data reduction and the modelling of the NIRSpec spectroscopy. J.M.H., G.H.R., S.A., K.N.H., S. Carniani, Z.J., S. Charlot, J.C., F.D'E., B.D.J., G.C.J., J.L., B.R., J.S., A.S., S.T., L.W. and J.W. contributed to the modelling and interpretation of the SED. K.N.H. contributed to the photometric redshift determination. F.S. contributed to the lensing magnification estimation. B.R. contributed to the imaging data visualization. K.N.H., M.J.R., B.R., C.C.W. and C.N.A.W contributed to the preflight NIRCam imaging data challenges. G.H.R., S.A., D.J.E., A.J.B., R.M., M.J.R., B.R., C.N.A.W. and C.W. contributed to the design of the MIRI, NIRCam and NIRSpec observations. G.H.R., S.A., K.N.H., R.B., E.E., M.J.R., F.S., C.C.W. and C.N.A.W. contributed to the design, construction and commissioning of MIRI, NIRCam and NIRSpec. W.M.B. and Y.Z. contributed comments to the paper.

## Competing interests

The authors declare no competing interests.

## Additional information

Extended data is available for this paper at
https://doi.org/10.1038/s41550-025-02503-z.

**Correspondence and requests for materials** should be addressed to Jakob M. Helton.

**Peer review information** *Nature Astronomy* thanks the anonymous reviewers for their contribution to the peer review of this work.

**Reprints and permissions information** is available at
www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**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/.

© The Author(s) 2025
¹Steward Observatory, University of Arizona, Tucson, AZ, USA. ²Center for Astrophysics, Harvard & Smithsonian, Cambridge, MA, USA.

³Scuola Normale Superiore, Pisa, Italy. ⁴Kavli Institute for Cosmology, University of Cambridge, Cambridge, UK. ⁵Cavendish Laboratory, University of Cambridge, Cambridge, UK. ⁶European Space Astronomy Centre (ESAC), European Space Agency (ESA), Madrid, Spain. ⁷Department of Physics, University of Oxford, Oxford, UK. ⁸Institut d'Astrophysique de Paris, Sorbonne Université, CNRS, UMR 7095, Paris, France. ⁹Department of Physics and Astronomy, University College London, London, UK. ¹⁰Centro de Astrobiologia (CAB), CSIC-INTA, Madrid, Spain. ¹¹Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA, USA. ¹²NSF's National Optical-Infrared Astronomy Research Laboratory, Tucson, AZ, USA.

¹³NRC Herzberg, Victoria, British Columbia, Canada. ☒e-mail: jakobhelton@arizona.edu
<figure><img src="image_6.png" /></figure>

**Extended Data Fig. 1 | An example of the model fitting photometry for JADES-GS-z14-0.** The individual exposures from the MIRI/F770W imaging are modeled assuming an extended morphology. These images are dominated by the diffraction width of the F770W PSF (FWHM = 0.269 arcseconds), which is nearly twice the measured half-light angular diameter with NIRCam (FWHM = 0.145 arcseconds). We measure photometry by fitting surface brightness profiles to the flux image of the object and its neighbors (left panel), allowing us

to construct accurate models (middle panel) which leave only slight residuals (right panel) relative to the data. The thumbnails are 2" × 2" which corresponds to roughly 6.5 pkpc × 6.5 pkpc at the observed redshift (z = 14.32). We note that the mean stacked flux image of the individual exposures is provided in the left panel for illustrative purposes, rather than the final image mosaic from Fig. 1 and Extended Data Figure 2.
<figure><img src="image_7.png" /><figcaption>**Extended Data Fig. 2** | An example of the PSF photometry for JADES-GS-z14-0.</figcaption></figure>

The final image mosaic for the MIRI/F770W imaging is modeled assuming a point source morphology for both galaxies. We measure photometry by fitting surface brightness profiles to the flux image of the object and its neighbors (left panel),

allowing us to construct an accurate model of the image (middle panel) that leaves only slight residuals (right panel) relative to the data. The thumbnails are 2" × 2" which corresponds to roughly 6.5 pkpc × 6.5 pkpc at the observed redshift.
<figure><img src="image_8.png" /></figure>

**Extended Data Fig. 3** | Constraints on joint posterior distributions and star formation histories (SFHs) with BAGPIPES. Various SED models with different SFH assumptions (constant in purple, delayed-tau in green, and continuity in orange) are implemented using BAGPIPES and the measured model fitting photometry. In the lower left, the joint posterior distributions are shown for stellar mass (M), stellar metallicity (Z), mass-weighted stellar age (t*), star-formation rate averaged over the previous 10 million years (SFR₁₀), diffuse

dust attenuation as measured in the V-band (A_V), and rest-frame equivalent width for the nebular emission lines Hβ and [OIII]λλ4959, 5007 (EW([OIII]+Hβ)). In the upper right, the derived SFHs are shown. The medians of the derived SFHs are the colored lines, while the 1σ confidence intervals are the shaded regions. This figure highlights degeneracies between some of the inferred physical parameters.
<figure><img src="image_9.png" /></figure>

**Extended Data Fig. 4** | SED modeling for JADES-GS-z14-0. Similar to Fig. 2, the measured spectral flux density and corresponding uncertainties are used to constrain the various SED models with Prospector. As provided in Table 1, the mean and 1σ standard deviation of the fiducial model fitting photometry are the red circles and error bars. The median of the SED models is the gray line and

unfilled squares, while the 1σ confidence interval is the shaded region. These results suggest that the excess flux in F770W relative to F444W is from stellar continuum rather than nebular emission line contributions, producing strong Balmer breaks in the models, which we find unlikely given the extreme inferred stellar masses and stellar ages.
<figure><img src="image_10.png" /></figure>

**Extended Data Fig. 5** | [OIII]/Hβ emission line ratio predictions for JADES-GS-z14-0. The physical properties inferred from SED modeling with BAGPIPES are used to predict the emission line ratio [OIII]/Hβ. This requires an assumption about the underlying continuum at rest-frame optical wavelengths. The left panel assumes the continuum is flat and consistent with the measured flux density in F444W, while the right panel assumes the posterior distribution of continuum levels from BAGPIPES. Various SED models with different SFH

assumptions (constant in purple, delayed-tau in green, and continuity in orange) are implemented. We report the 5th, 16th, 50th, 84th, and 95th percentiles after combining the posterior distributions from the various SED models. For comparison, we include measurements for a sample of galaxies at z ≈ 8 with Hβ line flux measurements. The median of this sample is the gray lines, while the 1σ confidence interval is the shaded regions.
<table><caption>Extended Data Table 1 | Inferred physical properties of JADES-GS-z14-0 from SED modeling with BAGPIPES</caption><thead><tr><th rowspan="2">Inferred Property</th><th colspan="3">Assumed Physical Model with BAGPIPES</th></tr><tr><th>Constant SFH</th><th>Delayed-Tau SFH</th><th>Continuity SFH</th></tr></thead><tbody><tr><td>log<sub>10</sub> (M<sub>*</sub>/M<sub>☉</sub>)</td><td>9.0<sup>+0.2</sup><sub>-0.7</sub></td><td>8.5<sup>+0.5</sup><sub>-0.2</sub></td><td>8.7<sup>+0.5</sup><sub>-0.2</sub></td></tr><tr><td>log<sub>10</sub> (Z<sub>*</sub>/Z<sub>☉</sub>)</td><td>-2.0<sup>+0.8</sup><sub>-0.7</sub></td><td>-1.5<sup>+0.4</sup><sub>-1.0</sub></td><td>-2.0<sup>+0.8</sup><sub>-0.7</sub></td></tr><tr><td>t<sub>*</sub>/Myr</td><td>22.0<sup>+23.5</sup><sub>-19.2</sub></td><td>6.1<sup>+22.6</sup><sub>-4.6</sub></td><td>22.1<sup>+29.0</sup><sub>-15.8</sub></td></tr><tr><td>SFR<sub>10</sub>/ [M<sub>☉</sub>/yr]</td><td>23.4<sup>+5.1</sup><sub>-3.8</sub></td><td>23.5<sup>+4.7</sup><sub>-4.2</sub></td><td>28.9<sup>+4.0</sup><sub>-7.7</sub></td></tr><tr><td>Σ<sub>SFR<sub>10</sub></sub>/ [M<sub>☉</sub>/yr/pc<sup>2</sup>]</td><td>59.6<sup>+12.9</sup><sub>-9.8</sub></td><td>59.8<sup>+12.0</sup><sub>-10.6</sub></td><td>73.5<sup>+10.3</sup><sub>-19.6</sub></td></tr><tr><td>sSFR<sub>10</sub>/Gyr<sup>-1</sup></td><td>26.9<sup>+74.8</sup><sub>-13.2</sub></td><td>86.4<sup>+15.2</sup><sub>-60.9</sub></td><td>58.5<sup>+44.2</sup><sub>-45.1</sub></td></tr><tr><td>A<sub>V</sub> / [AB mag]</td><td>0.60<sup>+0.06</sup><sub>-0.15</sub></td><td>0.52<sup>+0.11</sup><sub>-0.12</sub></td><td>0.55<sup>+0.06</sup><sub>-0.06</sub></td></tr><tr><td>log<sub>10</sub> (U)</td><td>-2.6<sup>+0.5</sup><sub>-0.8</sub></td><td>-2.5<sup>+0.3</sup><sub>-0.6</sub></td><td>-2.5<sup>+0.4</sup><sub>-0.8</sub></td></tr><tr><td>EW<sub>[OIII]+Hβ</sub>/Å</td><td>284<sup>+260</sup><sub>-83</sub></td><td>566<sup>+399</sup><sub>-226</sub></td><td>338<sup>+239</sup><sub>-90</sub></td></tr></tbody></table>

Columns: (1) name and units of the inferred property, (2) measurements using the constant SFH model, (3) measurements using the delayed-tau SFH model, and (4) measurements using the continuity SFH model. Reported measurements are median values and 1σ confidence intervals.
<table><caption>Extended Data Table 2 | Inferred physical properties of JADES-GS-z14-0 from SED modeling with Prospector</caption><thead><tr><th colspan="2">Assumed Physical Model with Prospector</th></tr><tr><th>Inferred Property</th><th>Continuity SFH</th></tr></thead><tbody><tr><td>log<sub>10</sub> (M*/M<sub>☉</sub>)</td><td>9.38<sup>+0.02</sup><sub>-0.02</sub></td></tr><tr><td>log<sub>10</sub> (Z*/Z<sub>☉</sub>)</td><td>0.08<sup>+0.02</sup><sub>-0.02</sub></td></tr><tr><td>t*/Myr</td><td>83.1<sup>+4.0</sup><sub>-3.0</sub></td></tr><tr><td>SFR<sub>10</sub> / [M<sub>☉</sub>/yr]</td><td>0.02<sup>+0.08</sup><sub>-0.02</sub></td></tr><tr><td>Σ<sub>SFR<sub>10</sub></sub> / [M<sub>☉</sub>/yr/pc<sup>2</sup>]</td><td>0.05<sup>+0.20</sup><sub>-0.05</sub></td></tr><tr><td>sSFR<sub>10</sub>/Gyr<sup>-1</sup></td><td>0.01<sup>+0.03</sup><sub>-0.01</sub></td></tr><tr><td>A<sub>V</sub> / [AB mag]</td><td>0.014<sup>+0.008</sup><sub>-0.008</sub></td></tr><tr><td>log<sub>10</sub> (U)</td><td>-2.1<sup>+0.4</sup><sub>-0.5</sub></td></tr><tr><td>EW<sub>[OIII]+Hβ</sub>/Å</td><td>2<sup>+16</sup><sub>-17</sub></td></tr></tbody></table>

Columns: (1) name and units of the inferred property and (2) measurements using the continuity SFH model. Reported measurements are median values and 1σ confidence intervals.