arxiv: 2601.06968 · v1 · submitted 2026-01-11 · 🌌 astro-ph.GA

Recognition: no theorem link

Nitrogen abundances in star-forming galaxies 2.2 Gyr after the Big Bang are not elevated

D. Schaerer , Y.I. Izotov , R. Marques-Chaves , C. C. Steidel , N. Reddy , A. E. Shapley , S. Mascia , J. Chisholm

show 17 more authors

S. R. Flury N. Guseva T. Heckman A. Henry A.K. Inoue I. Jung H. Kusakabe K. Mawatari P. Oesch G. Oestlin L. Pentericci N. Roy A. Saldana-Lopez R. Sato E. Vanzella A. Verhamme B. Wang

Authors on Pith no claims yet

Pith reviewed 2026-05-16 15:21 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords nitrogen abundanceN/O ratiohigh-redshift galaxiesJWST spectroscopychemical evolutionstar-forming galaxiesmetallicity

0 comments

The pith

Nitrogen-to-oxygen ratios in typical z~3 galaxies match those in local galaxies across the same metallicity range.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper measures nitrogen and oxygen abundances in 25 typical star-forming galaxies at redshift approximately 3 using deep JWST spectra and the direct method based on auroral line detections. It finds that the average N/O ratio and its scatter are the same as in low-redshift galaxies and HII regions over metallicities from 7.5 to 8.44 in 12+log(O/H). This result indicates that nitrogen enrichment has not evolved with redshift for typical galaxies over a large fraction of cosmic time. The finding matters because it shows that standard chemical enrichment processes were already operating similarly early in the universe's history, rather than requiring special high-redshift conditions for nitrogen production. The study also links this to BPT diagram positions, confirming the offset from local galaxies is not caused by higher nitrogen.

Core claim

In a sample of typical star-forming galaxies at z ~ 3 observed with JWST, direct-method abundances yield a mean log(N/O) = -1.29 with uncertainties, over a metallicity range 12+log(O/H) = 7.5 to 8.44. This value and its scatter match those seen in local galaxies and HII regions, demonstrating no redshift evolution of the N/O ratio for these objects. The galaxies also exhibit a BPT diagram offset similar to low-z Lyman continuum leakers, but this is unrelated to nitrogen abundance.

What carries the argument

Direct T_e method for nebular abundance determination using the [O III] λ4363 auroral line in rest-optical spectra.

If this is right

Establishes a local-like baseline for N/O at z~3 for future high-z comparisons.
Indicates that chemical evolution models need to produce similar N/O ratios by 2.2 Gyr after the Big Bang.
Shows that BPT offsets in high-redshift galaxies arise from factors other than enhanced nitrogen.
Supports using local scaling relations for interpreting early galaxy spectra in many cases.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Nitrogen production mechanisms, such as from AGB stars, must have been efficient early on to match local ratios.
The lack of evolution may point to similar star formation efficiencies or initial mass functions in typical galaxies across time.
This baseline allows testing for deviations in more extreme high-z populations like those with different star formation rates.
Could be extended by checking if oxygen and nitrogen yields are consistent with standard supernova and stellar models at these epochs.

Load-bearing premise

The subsample of galaxies with auroral line detections is representative of the broader population of typical star-forming galaxies at z~3, and the assumptions underlying the direct abundance method apply equally at high redshift.

What would settle it

A larger sample of z~3 galaxies without auroral lines or using alternative abundance methods like strong-line calibrations or photoionization models showing systematically higher N/O ratios would contradict the no-evolution claim.

Figures

Figures reproduced from arXiv: 2601.06968 by A. E. Shapley, A. Henry, A.K. Inoue, A. Saldana-Lopez, A. Verhamme, B. Wang, C. C. Steidel, D. Schaerer, E. Vanzella, G. Oestlin, H. Kusakabe, I. Jung, J. Chisholm, K. Mawatari, L. Pentericci, N. Guseva, N. Reddy, N. Roy, P. Oesch, R. Marques-Chaves, R. Sato, S. Mascia, S. R. Flury, T. Heckman, Y.I. Izotov.

**Figure 1.** Figure 1: JWST spectra from LyC22 showing object 1005 (top panel, LACES z = 3.676 LAE in SSA22) and 300009 (bottom, z = 2.767 LBG in the Westphal field). The blue line shows the adopted fit to the continuum. The main emission and absorption lines are indicated. 1.8 1.6 1.4 1.2 1.0 0.8 log(N/O) - this work 1.8 1.6 1.4 1.2 1.0 0.8 log(N/O) - other methods Strom+2017 Cataldi+2025 Cataldi - SDSS [PITH_FULL_IMAGE:figure… view at source ↗

**Figure 2.** Figure 2: Comparison of the N/O abundances derived following the direct method and the assumptions described in this work (x-axis) with those adopting the strong line calibrations of Cataldi et al. (2025b) based on JWST samples (in black; applicable for z > 1) and an SDSS sample (red, applicable at low-z), and the calibration from Strom et al. (2017) (blue). All strong line calibrations use the extinction-correcte… view at source ↗

**Figure 3.** Figure 3: Left: Classical emission line diagnostic diagram showing the galaxies from the LyC22 sample (blue and red symbols). LyC22 objects with significant [O iii] λ4363 detections are surrounded by a black square. LyC22 galaxies with log(N/O) > −1.2 are shown in red. We only show objects where the involved emission lines are detected at ≥ 3σ. Typical uncertainties are comparable to the size of the symbols, hence n… view at source ↗

**Figure 4.** Figure 4: Derived abundance ratio N/O from rest-optical lines of the LyC22 galaxies (blue circles) and low-z samples as a function of O/H. The low-z star-forming galaxies and Hii regions from the compilation of Izotov et al. (2023) and the LzLCS are shown by small black and grey symbols, respectively. Dash-dotted and dotted lines show the average trend observed in low-z star-forming galaxies, as parametrized by Vil… view at source ↗

read the original abstract

Using deep medium-resolution JWST rest-optical spectra of a sample of typical star-forming galaxies (Lyman break galaxies and Lyman-$\alpha$ emitters) from the LyC22 survey at $z \sim 3$, we determined the nebular abundances of N, O, and Ne relative to H for a subsample of 25 objects with the direct method, based on auroral [OIII]4363 line detections. Our measurements increases the number of accurate N/O determinations at $z \sim 2-4$ using a homogeneous approach. We found a mean value of $\log({\rm N/O})=-1.29^{+0.25}_{-0.21} $ over a metallicity range 12+log(O/H)=7.5 to 8.44. The observed N/O ratio and scatter are indistinguishable from that observed in low-z galaxies and HII regions over the same metallicity range, showing thus no redshift evolution of N/O for typical galaxies over a significant fraction of cosmic time. We also show that typical $z \sim 3$ galaxies show a similar offset in the BPT diagram as galaxies from the low-z Lyman Continuum Survey (LzLCS), when compared to the average of SDSS galaxies, and show that this offset is not due to enhanced nitrogen abundances. Our results establish a basis for future studies of the evolution of N and O at higher redshifts.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper adds 25 new direct-method N/O points at z~3 that line up with local values over the same metallicity range.

read the letter

The main result here is straightforward: typical star-forming galaxies at z~3 show a mean log(N/O) of -1.29 with scatter that matches local galaxies and HII regions across 12+log(O/H) from 7.5 to 8.44. The authors used JWST rest-optical spectra from the LyC22 survey and pulled direct abundances from the 25 objects where the auroral [OIII]4363 line was detected. This nearly doubles the number of homogeneous high-redshift direct N/O measurements in that range, which is the concrete addition worth noting. They also check that the BPT offset in these galaxies is not driven by nitrogen enhancement, consistent with the abundance numbers themselves. The analysis sticks to standard direct-method assumptions and tests them against local analogs, with no obvious internal contradictions in the ionization corrections or temperature handling once the full text is read. A soft spot is whether the auroral-line subsample fully represents the broader z~3 population. Galaxies bright enough in [OIII]4363 could favor higher-excitation conditions, and while the paper discusses selection, the generalization to all typical galaxies rests on that assumption holding. This is useful for people working on chemical evolution models or building baselines for even higher-redshift work. The data are observational, the comparison is direct, and the claim follows without circularity. I would send it to peer review. The sample increase is real and the central comparison is testable enough to benefit from referee input.

Referee Report

2 major / 3 minor

Summary. The manuscript reports direct-method nebular abundance measurements of N, O, and Ne for a subsample of 25 typical star-forming galaxies (Lyman break galaxies and Lyman-α emitters) at z ≈ 3 from the LyC22 survey, using deep JWST rest-optical spectra with auroral [O III] λ4363 detections. The authors derive a mean log(N/O) = −1.29^{+0.25}_{-0.21} over 12 + log(O/H) = 7.5–8.44 and find the N/O ratio and scatter indistinguishable from local galaxies and H II regions over the same metallicity range, concluding no redshift evolution in N/O for typical galaxies. They further show that the BPT diagram offset relative to SDSS is unrelated to nitrogen enhancement and comparable to the low-z Lyman Continuum Survey sample.

Significance. If the measurements and comparison hold, the result is significant because it supplies a homogeneous, direct-method anchor for N/O at z ∼ 3, substantially increasing the number of accurate high-redshift determinations and demonstrating that nitrogen enrichment in typical star-forming galaxies has not evolved detectably over 2.2 Gyr. This provides a firm empirical baseline for chemical-evolution models and future JWST studies at higher redshift.

major comments (2)

[§3 and §5] The central claim that the 25 auroral-line detections are representative of the full population of typical z ∼ 3 star-forming galaxies (and thus that N/O shows no evolution) rests on the assumption that selection effects tied to [O III] λ4363 detectability do not bias the N/O distribution. Section 3 (sample and observations) and the comparison in §5 should include a quantitative test (e.g., Kolmogorov–Smirnov or Anderson–Darling statistic on the parent-sample metallicity and excitation distributions) to demonstrate that the subsample is unbiased within the quoted uncertainties.
[§4 (abundance derivation)] The direct-method N/O values rely on standard local-calibrated ionization correction factors and the assumption of a single-zone temperature structure. The manuscript should report explicit sensitivity tests (varying T_e gradients or ICF prescriptions by ±20 %) in the abundance analysis section to show that the reported mean log(N/O) and its overlap with the local relation remain unchanged; without these, the no-evolution conclusion is vulnerable to systematic offsets that are not yet quantified at z ∼ 3.

minor comments (3)

[Abstract] Abstract: the sentence 'Our measurements increases the number...' contains a subject-verb agreement error and should read 'increase'.
[Figures 4–5] Figure 4 or 5 (N/O vs. O/H comparison): the local comparison sample should be explicitly labeled with its reference (e.g., 'Izotov et al. 2006' or equivalent) in the legend and caption for immediate clarity.
[§6 (BPT analysis)] The BPT offset discussion would benefit from a short statement of the exact line-ratio definitions used (e.g., [O III]/Hβ and [N II]/Hα) to avoid any ambiguity with alternative calibrations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their positive assessment of our manuscript and for the constructive suggestions that will improve the statistical robustness of our claims. We address each major comment below and will incorporate the requested quantitative tests and sensitivity analyses into the revised manuscript.

read point-by-point responses

Referee: [§3 and §5] The central claim that the 25 auroral-line detections are representative of the full population of typical z ∼ 3 star-forming galaxies (and thus that N/O shows no evolution) rests on the assumption that selection effects tied to [O III] λ4363 detectability do not bias the N/O distribution. Section 3 (sample and observations) and the comparison in §5 should include a quantitative test (e.g., Kolmogorov–Smirnov or Anderson–Darling statistic on the parent-sample metallicity and excitation distributions) to demonstrate that the subsample is unbiased within the quoted uncertainties.

Authors: We agree that a quantitative statistical comparison is needed to strengthen the representativeness argument. In the revised manuscript we will add Kolmogorov–Smirnov tests (and, if appropriate, Anderson–Darling tests) comparing the distributions of gas-phase metallicity (12 + log(O/H)) and excitation (e.g., [O III] λ5007 / Hβ) between the 25 auroral-line galaxies and the full parent LyC22 sample at z ≈ 3. The p-values and any implications for possible selection bias in the N/O distribution will be reported in Sections 3 and 5. revision: yes
Referee: [§4 (abundance derivation)] The direct-method N/O values rely on standard local-calibrated ionization correction factors and the assumption of a single-zone temperature structure. The manuscript should report explicit sensitivity tests (varying T_e gradients or ICF prescriptions by ±20 %) in the abundance analysis section to show that the reported mean log(N/O) and its overlap with the local relation remain unchanged; without these, the no-evolution conclusion is vulnerable to systematic offsets that are not yet quantified at z ∼ 3.

Authors: We acknowledge that explicit quantification of systematic uncertainties arising from the adopted ICFs and temperature structure is important, especially at high redshift. In the revised §4 we will include a dedicated sensitivity analysis in which (i) the nitrogen ICF is varied by ±20 % around the standard prescription and (ii) a two-zone temperature model is adopted with T_e differences of ±20 % between the high- and low-ionization zones. We will recompute the mean log(N/O) and its 16–84 percentile range under each variation and demonstrate that the overlap with the local N/O–metallicity relation is preserved within the quoted uncertainties. revision: yes

Circularity Check

0 steps flagged

No significant circularity; pure observational measurement

full rationale

The paper reports direct-method N/O and O/H measurements from auroral-line detections in JWST spectra of 25 z~3 galaxies. The central claim (no redshift evolution of N/O) follows from a straightforward statistical comparison of these observed values to local HII-region and galaxy samples over the same metallicity range. No equations, fitted parameters, or self-citations reduce the reported mean log(N/O) or its scatter to the input data by construction; the abundance derivations rely on standard atomic physics and ionization corrections that are independent of the redshift-evolution conclusion. The selection and temperature-structure assumptions are explicitly discussed but do not create a self-referential loop.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review; full details of nebular-analysis assumptions unavailable. Standard assumptions in direct-method abundance work are presumed but unverified here.

axioms (1)

domain assumption Direct-method abundance calculation using auroral [OIII]4363 is valid for z~3 galaxies
Invoked to derive N/O from observed line ratios.

pith-pipeline@v0.9.0 · 5686 in / 1139 out tokens · 26261 ms · 2026-05-16T15:21:35.731153+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

46 extracted references · 46 canonical work pages · 2 internal anchors

[1]

Z., Berg, D

Arellano-Córdova, K. Z., Berg, D. A., Mingozzi, M., et al. 2025, arXiv e-prints, arXiv:2507.11658 Arellano-Córdova, K. Z., Cullen, F., Carnall, A. C., et al. 2024, arXiv e-prints, arXiv:2412.10557

work page arXiv 2025
[2]

J., & Scott, P

Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481

work page 2009
[3]

A., Phillips, M

Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5

work page 1981
[4]

A., Erb, D

Berg, D. A., Erb, D. K., Henry, R. B. C., Skillman, E. D., & McQuinn, K. B. W. 2019, ApJ, 874, 93

work page 2019
[5]

A Fleeting GLIMPSE of N/O Enrichment at Cosmic Dawn: Evidence for Wolf Rayet N Stars in a z = 6.1 Galaxy

Berg, D. A., Naidu, R. P., Chisholm, J., et al. 2025, arXiv e-prints, arXiv:2511.13591

work page internal anchor Pith review Pith/arXiv arXiv 2025
[6]

2023, msaexp: NIRSpec analyis tools

Brammer, G. 2023, msaexp: NIRSpec analyis tools

work page 2023
[7]

2017, MNRAS, 466, 798

Cappellari, M. 2017, MNRAS, 466, 798

work page 2017
[8]

A., Clayton, G

Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245

work page 1989
[9]

Tracing nitrogen enrichment across cosmic time with JWST

Cataldi, E., Belfiore, F., Curti, M., et al. 2025b, arXiv e-prints, arXiv:2512.07955

work page internal anchor Pith review Pith/arXiv arXiv
[10]

2012, MNRAS, 427, 1953

Christensen, L., Richard, J., Hjorth, J., et al. 2012, MNRAS, 427, 1953

work page 2012
[11]

& Gunn, J

Conroy, C. & Gunn, J. E. 2010, FSPS: Flexible Stellar Population Synthesis

work page 2010
[12]

Edmunds, M. G. & Pagel, B. E. J. 1978, MNRAS, 185, 77P Fernández, V., Amorín, R., Firpo, V., & Morisset, C. 2024, A&A, 688, A69

work page 1978
[13]

J., Tang, M., Robertson, B

Fletcher, T. J., Tang, M., Robertson, B. E., et al. 2019, ApJ, 878, 87

work page 2019
[14]

G., Izotov, Y

Guseva, N. G., Izotov, Y. I., Stasińska, G., et al. 2011, A&A, 529, A149

work page 2011
[15]

J., Saldana-Lopez, A., Citro, A., et al

Hayes, M. J., Saldana-Lopez, A., Citro, A., et al. 2025, ApJ, 982, 14

work page 2025
[16]

Henry, R. B. C., Edmunds, M. G., & Köppen, J. 2000, ApJ, 541, 660

work page 2000
[17]

I., Orlitová, I., Schaerer, D., et al

Izotov, Y. I., Orlitová, I., Schaerer, D., et al. 2016, Nature, 529, 178

work page 2016
[18]

I., Schaerer, D., Worseck, G., et al

Izotov, Y. I., Schaerer, D., Worseck, G., et al. 2023, MNRAS, 522, 1228

work page 2023
[19]

I., Stasińska, G., Meynet, G., Guseva, N

Izotov, Y. I., Stasińska, G., Meynet, G., Guseva, N. G., & Thuan, T. X. 2006, A&A, 448, 955

work page 2006
[20]

I., Worseck, G., Schaerer, D., et al

Izotov, Y. I., Worseck, G., Schaerer, D., et al. 2018, MNRAS, 478, 4851

work page 2018
[21]

E., Silveyra, A

Jaskot, A. E., Silveyra, A. C., Plantinga, A., et al. 2024, ApJ, 972, 92

work page 2024
[22]

I., Stasińska, G., Meynet, G., Guseva, N

Ji, X., Belokurov, V., Maiolino, R., et al. 2025, arXiv e-prints, arXiv:2505.12505 Article number, page 6 of 7 Schaerer et al.: Normal N/O abundances at redshiftz∼3

work page arXiv 2025
[23]

M., White, S

Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003, MN- RAS, 341, 33

work page 2003
[24]

J., Heisler, C

Kewley, L. J., Heisler, C. A., Dopita, M. A., & Lumsden, S. 2001, ApJS, 132, 37

work page 2001
[25]

J., Maier, C., Yabe, K., et al

Kewley, L. J., Maier, C., Yabe, K., et al. 2013b, ApJL, 774, L10 Le Fèvre, O., Lemaux, B. C., Nakajima, K., et al. 2019, A&A, 625, A51

work page 2019
[26]

2025, arXiv e-prints, arXiv:2510.12411

Marques-Chaves, R., Schaerer, D., Dessauges-Zavadsky, M., et al. 2025, arXiv e-prints, arXiv:2510.12411

work page arXiv 2025
[27]

2016, The Astrophysical Journal, 828, 18

Masters, D., Faisst, A., & Capak, P. 2016, The Astrophysical Journal, 828, 18

work page 2016
[28]

2025, arXiv e- prints, arXiv:2511.20484

Morel, I., Schaerer, D., Marques-Chaves, R., et al. 2025, arXiv e- prints, arXiv:2511.20484

work page arXiv 2025
[29]

Groves, B. A. 2017, MNRAS, 466, 4403

work page 2017
[30]

J., Shapley, A., Steidel, C

Pahl, A. J., Shapley, A., Steidel, C. C., Chen, Y., & Reddy, N. A. 2021, MNRAS, 505, 2447

work page 2021
[31]

F., & Tsang, B

Pascale, M., Dai, L., McKee, C. F., & Tsang, B. T. H. 2023, ApJ, 957, 77 Pérez-Montero, E. & Contini, T. 2009, MNRAS, 398, 949

work page 2023
[32]

2023, ApJ, 948, 83

Brammer, G. 2023, ApJ, 948, 83

work page 2023
[33]

L., Shapley, A

Sanders, R. L., Shapley, A. E., Clarke, L., et al. 2023, ApJ, 943, 75

work page 2023
[34]

L., Shapley, A

Sanders, R. L., Shapley, A. E., Topping, M. W., et al. 2025, arXiv e-prints, arXiv:2508.10099 Schaerer,D.,Marques-Chaves,R.,Atek,H.,etal.2025,arXive-prints, arXiv:2512.16549

work page arXiv 2025
[35]

C., et al

Scholte, D., Cullen, F., Carnall, A. C., et al. 2025, arXiv e-prints, arXiv:2502.10499

work page arXiv 2025
[36]

Shapley, A. E. 2011, ARA&A, 49, 525

work page 2011
[37]

E., Reddy, N

Shapley, A. E., Reddy, N. A., Kriek, M., et al. 2015, ApJ, 801, 88

work page 2015
[38]

C., Bogosavljević, M., Shapley, A

Steidel, C. C., Bogosavljević, M., Shapley, A. E., et al. 2018, ApJ, 869, 123

work page 2018
[39]

C., Rudie, G

Steidel, C. C., Rudie, G. C., Strom, A. L., et al. 2014, ApJ, 795, 165

work page 2014
[40]

C., Strom, A

Steidel, C. C., Strom, A. L., Pettini, M., et al. 2016, ApJ, 826, 159

work page 2016
[41]

2025, ApJ, 981, 136

Stiavelli, M., Morishita, T., Chiaberge, M., et al. 2025, ApJ, 981, 136

work page 2025
[42]

L., Steidel, C

Strom, A. L., Steidel, C. C., Rudie, G. C., et al. 2017, ApJ, 836, 164

work page 2017
[43]

W., Stark, D

Topping, M. W., Stark, D. P., Senchyna, P., et al. 2024, MNRAS, 529, 3301

work page 2024
[44]

Vila-Costas, M. B. & Edmunds, M. G. 1993, MNRAS, 265, 199

work page 1993
[45]

& Kobayashi, C

Vincenzo, F. & Kobayashi, C. 2018, MNRAS, 478, 155

work page 2018
[46]

E., Rigby, J

Welch, B., Rivera-Thorsen, T. E., Rigby, J. R., et al. 2025, ApJ, 980, 33 1 DepartmentofAstronomy,UniversityofGeneva,CheminPe- gasi 51, 1290 Versoix, Switzerland 2 CNRS, IRAP, 14 Avenue E. Belin, 31400 Toulouse, France 3 Bogolyubov Institute for Theoretical Physics, National Academy of Sciences of Ukraine, 14-b Metrolohichna str., Kyiv, 03143, Ukraine 4 C...

work page 2025