JWST/NIRSpec EXCELS Survey Overview
- The JWST/NIRSpec EXCELS Survey is a pioneering deep spectroscopic program using medium-resolution data to study galaxies across z ≃1.6–8 with unprecedented metallicity diagnostics.
- It employs a multi-shutter assembly with three grating/filter combinations, providing simultaneous rest-frame UV and optical coverage with 4–5.5 hr integrations per pointing.
- Observations reveal detailed trends in ionization, chemical enrichment, and feedback processes, setting new benchmarks for understanding galaxy evolution from reionization to cosmic noon.
The JWST/NIRSpec EXCELS Survey (“EXtraordinarily deep spectroscopic survey of the Cosmic EarLy SystEms”) is a comprehensive program exploiting the sensitivity and spectral multiplexing capabilities of JWST/NIRSpec to investigate the physical conditions, chemical enrichment, ionization, feedback, and assembly pathways of galaxies from the epoch of reionization through cosmic noon. Using medium-resolution (R≃1000) spectroscopy across 1–5.3 μm, the survey targets a broad redshift (1.6<z<8.0) and stellar-mass range, enabling homogeneous studies of both star-forming and quiescent systems. EXCELS delivers rest-frame optical diagnostics for the highest-redshift galaxies accessed to date, with a suite of new direct-method metallicity measurements, strong-line calibrations, and insights into the co-evolution of stellar populations, gas flows, and supermassive black holes.
1. Survey Architecture and Observational Strategy
EXCELS is structured around the use of JWST/NIRSpec in its Multi-Shutter Assembly (MSA) mode, configured with three medium-resolution grating/filter combinations: G140M/F100LP (1.0–1.8 μm), G235M/F170LP (1.7–3.1 μm), and G395M/F290LP (2.9–5.3 μm). Each grating provides R ∼ 1000 spectroscopy, allowing simultaneous coverage of rest-frame UV and optical features for galaxies over a wide redshift span. Typical integration times are ≈4 hr (G140M, G395M) and ≈5.5 hr (G235M) per pointing.
Target selection leverages deep HST and JWST/NIRCam photometry in the PRIMER UDS field, with robust photometric and/or spectroscopic pre-selection spanning massive quiescent galaxies (log₁₀(Mₛₜₐᵣ/Mₒ)>10 at 3<z<5), post-starburst (PSB) galaxies, and star-forming systems from the main sequence to the low-mass regime (8.1<log₁₀(Mₛₜₐᵣ/Mₒ)<10.3). The MSA is configured to maximize overlap and multiplex for high-priority targets, resulting in four MSA pointings covering ≃160 arcmin².
Raw NIRSpec data are reduced using the JWST pipeline with numerous customizations: 1/f noise removal, snowball/cosmic-ray rejection, precise dither combination, and optimal extraction techniques that mitigate spectrally undersampled PSF “wiggles.” Absolute flux calibration is achieved through synthetic integration of the spectra over observed photometric filters, ensuring global consistency with multi-band photometry (Scholte et al., 14 Feb 2025, Leung et al., 5 Feb 2026, Carnall et al., 2024).
2. Direct Metallicity Measurements and Strong-Line Diagnostics
A cornerstone of EXCELS is the high-reliability detection of auroral [O III] λ4363 in 22 galaxies (1.65<z\<7.92, median z=4.05), doubling the census of Tₑ-based abundance measurements at z\>3. The Tₑ-method employs the [O III] λ4363/λ5007 flux ratio, modeled with PyNeb in a two-zone approximation (O²⁺—high ionization; O⁺, N⁺—low ionization). Ionic abundances are computed directly and converted into total O/H and N/O using standard empirical temperature relations.
Single-ratio strong-line metallicity calibrations (e.g. R₂, R₃, O₃₂, N2, O3N2) suffer from systematic evolution with ionization parameter (U) and relative abundance (notably N/O). EXCELS implements and recalibrates the redshift-invariant oxygen-only diagnostic [Laseter et al. 2024], a linear combination of R₂ = log([O II]/Hβ) and R₃ = log([O III]/Hβ) with:
with a fitted cubic/quartic relation to O/H, achieving 0.11–0.14 dex scatter, independent of U and redshift. Additionally, the diagnostic extends robust metallicity estimation to z > 9.5 by exploiting accessible [O II] and [Ne III] lines:
with excellent performance and minimal dust or U sensitivity, valid over (Scholte et al., 14 Feb 2025).
3. Chemical Enrichment: Trends, Abundance Ratios, and Enrichment Pathways
EXCELS establishes precise mass-metallicity relations (MZR) and fundamental metallicity relation (FMR) evolution at 2<z\<8, anchored by both direct and strong-line measurements across 65 star-forming galaxies. The MZR exhibits a modest decrease in normalization (ΔZ₉≈0.10±0.04 dex from z=3.2 to z=5.5) and unchanged slope (γ≈0.3), with galaxies at fixed mass attaining ~50% of their present-day O/H by z=3 (first 15% of cosmic time). The FMR, as parameterized by Andrews & Martini (2013) and Curti et al. (2020), displays systematic offsets (Δ=–0.17 to –0.37 dex at z=3–7), reflecting non-equilibrium conditions in high-sSFR galaxies at z\>3 (Stanton et al., 1 Nov 2025).
Abundance patterns reveal elevated N/O (log(N/O)≈–1.1 to –1.3 at 12+log(O/H)≈8.0), exceeding local medians by 0.2–0.3 dex, and high N/O outliers at z≈5 relative to typical star-forming galaxies. Direct Tₑ-based C/O and C/N in two z≃5 galaxies indicate low C/O (–1.02≤log(C/O)≤–0.82) and depressed C/N, aligning with core-collapse SN yields and chemical evolution models favoring bottom-heavy or standard IMFs with rapid N enrichment from 4–7 Mₒ AGB stars (Arellano-Córdova et al., 2024). In extreme cases, such as EXCELS-63107 at z=8.271, metallicities of 12+log(O/H)=6.89 are reached, requiring a top-heavy IMF or Population III-like sources to achieve the hard ionizing spectrum and nebular excitation ([O III] λ4363/λ5007, β=–3.3) (Cullen et al., 19 Jan 2025).
4. Ionization Conditions, Ionizing Efficiency, and Star Formation Feedback
Multi-grating NIRSpec data enable robust determinations of ionizing photon production efficiencies (ξ₍ᵢₒₙ₎) across a wide dynamic range in Hα and [O III] equivalent widths (EWs), UV slopes (β_UV), dust attenuation, and SFRs. A sample of 159 SFGs at 1<z<8 reveals that ξ₍ᵢₒₙ,₀₎ is dominantly regulated by W(Hα), W([O III]), β_UV, and E(B–V)ₙₑb, with multivariate regression yielding tight predictive relations (scatter σ_int~0.04 dex, MAE~0.15 dex). SFGs with W(Hα) ≳ 250 Å or β < –1.6 consistently surpass the canonical reionization threshold (⟨log ξ₍ᵢₒₙ,₀₎⟩ ≳ 25.2). High E(B–V)ₙₑb sources form a population of outliers, likely due to heavily obscured star formation or differential dust geometry, emphasizing the need for caution in using single indicators in highly reddened galaxies (Begley et al., 30 Sep 2025).
5. Gas Flows, Outflows, and Feedback Across Galaxy Types
High-S/N, medium-resolution spectra enable detailed studies of both neutral and low-ionization gas flows via Na D, Fe II, and Mg II absorption. In star-forming galaxies at z∼3, outflows traced by Fe II and Mg II show mean velocities of –50 to –80 km s⁻¹, with detection likelihood strongly increasing for higher M*, SFR, and Σ_SFR, though no statistically significant scaling of outflow velocity with galaxy property is found individually. Mg II emission (“emitters”) is preferentially seen in low-mass, low-A_V, and high-sSFR galaxies, consistent with a porous neutral ISM promoting resonant-line and LyC escape. Stacking by inclination demonstrates that outflows manifest most strongly along galaxy minor axes (face-on view), while inflows predominantly occur along the disk plane, matching predictions from simulations and local outflow observations (Kehoe et al., 10 Jun 2026).
In quiescent and PSB galaxies at 1.8<z<4.6, ~23% display strong, blueshifted NaD absorption (v_flow=300–1300 km s⁻¹), with outflow mass-loading (Ṁ_out=20–200 Mₒ yr⁻¹) and energy/momentum rates inconsistent with ongoing SFR or detected AGN, but in line with “fossil” outflows from short-lived (<5 Myr), episodic AGN activity. Simulations (EAGLE) and incidence statistics are consistent with duty cycles of t_cycle≃40 Myr, with AGN “on” for ≃5 Myr, outflow visible ≃10 Myr, and a ∼20 Myr inflow phase, supporting a “flicker feedback” maintenance mechanism for quenching (Taylor et al., 5 Jan 2026).
6. Galaxy Assembly Histories, Downsizing, and Implications for Early Structure Formation
Full spectral fitting of quiescent galaxies at 3<z\<5 (N=14; S/N\>10 Å⁻¹) confirms that stellar mass assembly is highly accelerated in the early universe: the t_form–M* relation has slope Δt_form/Δlog₁₀M ≃ 1.5 Gyr dex⁻¹, with the most massive galaxies forming by z≳7, in concordance with archaeological downsizing established at z<3. Stellar metallicities reach solar or even super-solar values at log₁₀(M*/Mₒ) > 10.5. Attempts to measure α/Fe find large model-dependent uncertainties and require substantially higher S/N (≥100 per R∼1000 element). Robust assembly durations (τ₁₀₋₉₀ ≈ 200 Myr) imply formation SFRs in excess of 300 Mₒ yr⁻¹ (Leung et al., 5 Feb 2026).
Ultra-massive (log₁₀(M*/Mₒ)>11) quiescent galaxies at 3<z<5, among the earliest such systems known, possess exceptionally high stellar metallicities (Z_\approx1–2 Z_⊙), compact sizes (r_e<1 kpc), and in at least one case, strong α-enhancement ([Mg/Fe]=0.42{+0.19}_{-0.17}). Extreme value statistics show their abundance is only consistent with ΛCDM cosmology if stellar baryon conversion fractions f_=1 are permitted, implying strikingly efficient baryonic assembly and a fundamentally different regime of early galaxy formation (Carnall et al., 2024).
7. Broader Impact and Future Directions
The EXCELS program delivers a multi-faceted dataset that connects ionized gas-phase diagnostics, stellar metallicities, abundance ratios, star formation rates, and feedback-induced gas flows at the highest redshifts yet accessible with rest-frame optical spectroscopy. The redshift-invariant and calibrations provide methodological advances now central in interpreting high-z spectra where traditional diagnostics fail. Chemical enrichment results highlight significant departures in CNO patterns at fixed O/H versus the local universe, tied to IMF variations and enrichment timescales. The identification of episodic AGN-driven outflows as a general feature of post-starburst galaxies, and the direct mapping of baryon cycling geometry via inclination-resolved absorption, set new benchmarks for models of galaxy feedback and quenching. Finally, the census of quiescent massive galaxies constrains the tail of the halo mass function and informs on the earliest assembly of the red sequence.
Future work will require both greater sample sizes and higher S/N to constrain abundance patterns at the ≲0.1 dex level, leveraging advances in α-enhanced stellar models and deeper spectral integration. Integration of direct-method abundances, multi-phase gas diagnostics, and resolved kinematic IFU follow-up will further refine constraints on both galaxy formation physics and the cosmological context for the emergence of structure in the first billion years (Scholte et al., 14 Feb 2025, Skarbinski et al., 22 Sep 2025, Leung et al., 5 Feb 2026, Carnall et al., 2024, Cullen et al., 19 Jan 2025, Taylor et al., 5 Jan 2026, Stanton et al., 1 Nov 2025, Arellano-Córdova et al., 2024, Begley et al., 30 Sep 2025, Kehoe et al., 10 Jun 2026).