The JWST/AURORA Survey: Multiple Balmer and Paschen Emission Lines for Individual Star-forming Galaxies at z=1.5-4.4. I. A Diversity of Nebular Attenuation Curves and Evidence for Non-Unity Dust Covering Fractions

arxiv: 2506.17396 · v2 · pith:RMRKXZYAnew · submitted 2025-06-20 · 🌌 astro-ph.GA

The JWST/AURORA Survey: Multiple Balmer and Paschen Emission Lines for Individual Star-forming Galaxies at z=1.5-4.4. I. A Diversity of Nebular Attenuation Curves and Evidence for Non-Unity Dust Covering Fractions

Naveen A. Reddy , Alice E. Shapley , Ryan L. Sanders , Michael W. Topping , Richard S. Ellis , Max Pettini , Gabriel Brammer , Fergus Cullen

show 9 more authors

Natascha M. Forster Schreiber Ali A. Khostovan Derek J. McLeod Ross J. McLure Desika Narayanan Pascal A. Oesch Anthony J. Pahl Charles C. Steidel Danielle A. Berg

This is my paper

Pith reviewed 2026-05-19 07:53 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords nebular attenuation curvesdust covering fractionBalmer emission linesPaschen emission lineshigh-redshift star-forming galaxiesreddening measurementsISM geometry

0 comments p. Extension

pith:RMRKXZYA Add to your LaTeX paper

What is a Pith Number?

\usepackage{pith}
\pithnumber{RMRKXZYA}

Prints a linked pith:RMRKXZYA badge after your title and writes the identifier into PDF metadata. Compiles on arXiv with no extra files. Learn more

The pith

Star-forming galaxies at z=1.5-4.4 show diverse nebular attenuation curves explained by Milky Way-like dust grains with only partial covering of OB associations.

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

The paper measures multiple Balmer and Paschen emission lines in 24 individual star-forming galaxies at redshifts 1.5 to 4.4. It reports that Paschen-based reddening often exceeds Balmer-based values under a standard Galactic curve, indicating some heavily obscured star formation alongside less obscured regions. The derived attenuation curves display a wide range of normalizations, with Rv values from about 3.2 to 16.4. The authors show that this spread arises naturally when dust has Milky Way-like properties but covers only a fraction of the young stellar regions, with the covering fraction typically between 0.6 and 1.0. Such measurements of line ratios then offer a direct probe of how dust and gas obscure or allow escape of radiation from these distant galaxies.

Core claim

By analyzing multiple Balmer and Paschen lines, we find nebular reddening from Paschen lines exceeds that from Balmer lines for at least half the sample under the Galactic extinction curve, implying heavily reddened star formation. The nebular attenuation curves show a broad range of normalizations with Rv ~ 3.2-16.4. These offsets and high Rv values indicate variations in dust-stars geometry. We propose a model with subunity dust covering fraction fcov, and fitting to the HI recombination line ratios gives fcov ~ 0.6-1.0. The normalizations of the attenuation curves are driven primarily by fcov and the mix of reddened and unreddened OB associations. Thus the diversity of nebular attenuation

What carries the argument

Subunity dust covering fraction model fitted to HI recombination line ratios, which mixes contributions from reddened and unreddened OB associations to reproduce the observed range of attenuation curve normalizations.

If this is right

Integrated measurements of multiple Balmer and Paschen lines can place novel constraints on the dust covering fraction towards OB associations.
These covering fractions provide new avenues for exploring the impact of stellar feedback on ISM porosity.
The covering fraction influences the escape of Ly-alpha and Lyman continuum radiation from high-redshift galaxies.

Where Pith is reading between the lines

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

The geometry variations implied here may connect to how stellar feedback clears channels in the interstellar medium across different galaxy masses or environments.
Spatially resolved JWST observations could test whether the inferred covering fractions match the actual distribution of dust around young stars.
This line-ratio approach could be applied to lower-redshift samples to check whether the same partial-covering picture holds when more direct dust maps are available.

Load-bearing premise

The offsets between Balmer- and Paschen-derived reddening are due to variations in dust-stars geometry rather than differences in intrinsic dust properties or other unmodeled effects.

What would settle it

Finding that Paschen-derived reddening matches Balmer-derived reddening for the full sample when using a single fixed extinction curve, without any need for variable covering fractions, would falsify the proposed model.

read the original abstract

We present the nebular attenuation curves and dust covering fractions for 24 redshift z=1.5-4.4 star-forming galaxies using multiple Balmer and Paschen lines from the JWST/AURORA survey. Nebular reddening derived from Paschen lines exceeds that from Balmer lines for at least half the galaxies in the sample when assuming the commonly-adopted Galactic extinction curve, implying the presence of heavily reddened star formation. The nebular attenuation curves exhibit a broad range of normalizations (Rv ~ 3.2-16.4). Motivated by the offsets in reddening deduced from the Balmer and Paschen lines, and the high Rv values for the individual nebular attenuation curves, both of which suggest variations in the dust-stars geometry, we propose a model with a subunity dust covering fraction (fcov). Fitting such a model to the HI recombination line ratios indicates fcov ~ 0.6-1.0. The normalizations of the nebular attenuation curves, Rv, are driven primarily by fcov and the mix of reddened and unreddened OB associations. Thus, the diversity of nebular attenuation curves can be accommodated by assuming dust grain properties similar to that of Milky Way sightlines but with a subunity covering fraction of dust. Integrated measurements of multiple Balmer and Paschen lines can be used to place novel constraints on the dust covering fraction towards OB associations. These, in turn, provide new avenues for exploring the role of dust and gas covering fraction in a number of relevant aspects of high-redshift galaxies, including the impact of stellar feedback on ISM porosity and the escape of Ly-alpha and Lyman continuum radiation.

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 delivers new individual-galaxy constraints on high-redshift dust geometry from Balmer and Paschen lines, but leaves room for temperature effects to explain the key offsets.

read the letter

Colleague, This paper uses JWST observations from the AURORA survey to measure several Balmer and Paschen emission lines in 24 individual star-forming galaxies between z=1.5 and 4.4. The key finding is that for many of these objects, the reddening inferred from Paschen lines is higher than from Balmer lines under a standard Galactic curve. This leads to a wide range of Rv values from about 3 to 16, which they attribute to a mix of reddened and unreddened star-forming regions with a dust covering fraction below one. What the work does well is move beyond stacked spectra to look at individual galaxies. Having multiple lines per object allows them to map out the attenuation curve shape and normalization for each one. The data come from direct spectroscopy, so the line ratio measurements themselves look reliable. They then fit a two-component geometric model with fcov between 0.6 and 1.0, which reproduces the observed diversity while assuming Milky Way-like dust grains. The soft spot is in the interpretation of those Balmer-Paschen offsets. The stress-test note points out that electron temperature differences can change the intrinsic ratios by several percent per 1000 K, which would shift the derived E(B-V) without any change in dust geometry. The manuscript does not appear to test whether the offsets persist after allowing Te to vary or after using other line diagnostics to pin down temperature. If temperature variations contribute significantly, the need for the covering fraction model becomes less clear. The fitting of fcov directly to the same line ratios also introduces some dependence, though that's common in such models. This paper is aimed at researchers working on the interstellar medium and dust in high-redshift galaxies. Anyone thinking about how dust geometry affects star formation rate estimates, feedback, or the escape of ionizing radiation would find the constraints interesting. The new data on individual objects makes it worth a serious referee's time, even with the open questions around alternative explanations. I would send it out for peer review.

Referee Report

2 major / 2 minor

Summary. The paper reports JWST/AURORA spectroscopic measurements of multiple Balmer and Paschen recombination lines in 24 individual star-forming galaxies at z=1.5-4.4. It derives nebular attenuation curves that exhibit a wide range of normalizations (Rv ~ 3.2-16.4) and finds that Paschen-derived reddening often exceeds Balmer-derived reddening when assuming a Galactic extinction curve. Motivated by these offsets and high Rv values, the authors propose a two-component geometric model with Milky Way-like dust grains and a sub-unity dust covering fraction (fcov ~ 0.6-1.0) toward OB associations, attributing the diversity of attenuation curves to variations in dust-stars geometry rather than changes in intrinsic grain properties.

Significance. If the geometric interpretation holds, the work supplies direct observational constraints on dust covering fractions in high-redshift galaxies using multiple recombination lines, opening avenues to study ISM porosity, stellar feedback, and the escape of Lyα and Lyman-continuum radiation. The JWST-enabled access to Paschen lines alongside Balmer lines is a clear strength, providing a new empirical handle on geometry that is independent of continuum-based methods.

major comments (2)

[Abstract and proposed model section] Abstract and proposed model section: the central claim that Balmer-Paschen reddening offsets arise from sub-unity fcov with fixed case-B ratios at a single Te (~10^4 K) is load-bearing. Electron temperature variations of a few thousand K alter intrinsic Hα/Hβ, Hβ/Paβ, etc. ratios by several percent (and more at low density), which can shift apparent E(B-V) values without any change in geometry or grain properties. The manuscript does not report a systematic test that the offsets survive after allowing Te to vary (e.g., via [O III] or [S II] diagnostics) or after marginalizing over a Te prior in the fcov fit.
[Results on nebular attenuation curves] Results on nebular attenuation curves and fcov fitting: fcov is fitted directly to the same Balmer/Paschen line ratios used to measure the reddening offsets, creating moderate circularity. The manuscript should quantify the number of free parameters, show that the model has predictive power on held-out lines or independent diagnostics, and demonstrate that the derived Rv range (3.2-16.4) is not an artifact of the fitting procedure itself.

minor comments (2)

Figure captions and legends should explicitly label which curves correspond to individual galaxies versus the sample median or model predictions for clarity.
Add a short paragraph comparing the derived fcov values to independent constraints from UV continuum slopes or Lyα escape fractions in the same redshift range.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments, which have helped us identify areas where the manuscript can be strengthened. We address each major comment below, agreeing that additional tests are warranted to bolster the robustness of our conclusions regarding the geometric model.

read point-by-point responses

Referee: [Abstract and proposed model section] Abstract and proposed model section: the central claim that Balmer-Paschen reddening offsets arise from sub-unity fcov with fixed case-B ratios at a single Te (~10^4 K) is load-bearing. Electron temperature variations of a few thousand K alter intrinsic Hα/Hβ, Hβ/Paβ, etc. ratios by several percent (and more at low density), which can shift apparent E(B-V) values without any change in geometry or grain properties. The manuscript does not report a systematic test that the offsets survive after allowing Te to vary (e.g., via [O III] or [S II] diagnostics) or after marginalizing over a Te prior in the fcov fit.

Authors: We agree that electron temperature variations can affect the intrinsic recombination line ratios at the level of several percent, potentially contributing to apparent reddening offsets. Our analysis adopted the standard fixed Te ≈ 10^4 K assumption for case B recombination, but we did not perform a systematic marginalization or diagnostic-based test across the sample. In the revised manuscript, we will incorporate available [O III] λ4363 and [S II] temperature diagnostics where measurable, perform a sensitivity analysis allowing Te to vary by ±2000 K, and include a Te prior in the fcov fitting procedure to demonstrate that the observed Balmer-Paschen offsets and high Rv values persist under these variations. revision: yes
Referee: [Results on nebular attenuation curves] Results on nebular attenuation curves and fcov fitting: fcov is fitted directly to the same Balmer/Paschen line ratios used to measure the reddening offsets, creating moderate circularity. The manuscript should quantify the number of free parameters, show that the model has predictive power on held-out lines or independent diagnostics, and demonstrate that the derived Rv range (3.2-16.4) is not an artifact of the fitting procedure itself.

Authors: We acknowledge the risk of circularity since the same line ratios inform both the reddening offsets and the fcov fits. The model primarily varies fcov (one free parameter per galaxy) along with a single E(B-V) for the covered component, assuming Milky Way grain properties and case B ratios. To address this, the revised version will explicitly state the number of free parameters, implement a cross-validation test by fitting fcov to Balmer lines and predicting Paschen line ratios (and vice versa), and compare the resulting attenuation curves against independent continuum-based reddening estimates from the spectra. These additions will confirm that the reported Rv range reflects the data rather than fitting artifacts. revision: yes

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central model introduces fcov as a fitted parameter and assumes Milky Way-like dust properties; no new entities are postulated.

free parameters (1)

fcov = 0.6-1.0
Dust covering fraction fitted to match observed Balmer and Paschen line ratios for each galaxy.

axioms (1)

domain assumption Galactic extinction curve is a valid baseline for nebular attenuation
Used to compare Balmer and Paschen derived reddening and motivate the need for variable covering.

pith-pipeline@v0.9.0 · 5965 in / 1204 out tokens · 23710 ms · 2026-05-19T07:53:54.422942+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Fitting such a model to the HI recombination line ratios indicates fcov ~ 0.6-1.0. The normalizations of the nebular attenuation curves, Rv, are driven primarily by fcov and the mix of reddened and unreddened OB associations.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

The Internal Nebular Attenuation Curve of Three-Dimensional Turbulent HII regions
astro-ph.GA 2026-05 unverdicted novelty 6.0

Clumpy turbulent structures in HII regions leave the slope of the internal nebular attenuation curve unchanged because dust obscuration in dense clumps is offset by higher local emission-line production.
Calibrating Photometric Mid-Infrared Star Formation Rates for JWST
astro-ph.GA 2026-04 unverdicted novelty 4.0

Rest-frame 6-8um MIRI luminosity provides broken power-law SFR calibrations with 0.2-0.3 dex scatter and UV+IR composites at 0.15 dex, supporting robust use above log M* ~9 up to z~3.