arxiv: 2604.11879 · v1 · submitted 2026-04-13 · 🌌 astro-ph.GA

Recognition: unknown

A Deep ALMA Survey of the Redshift Distribution of Dusty Star-forming Galaxies

S. J. McKay , A. J. Barger , L. L. Cowie , F. E. Bauer

Authors on Pith no claims yet

Pith reviewed 2026-05-10 14:55 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords dusty star-forming galaxiesredshift distributionALMA spectroscopysubmillimeter galaxiesGOODS-Sphotometric redshiftshigh-redshift galaxiesearly universe

0 comments

The pith

ALMA spectroscopy shows fewer than 10 percent of bright 870-micron dusty galaxies lie at redshifts above 4.

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

The paper conducts an ALMA spectroscopic survey of dusty star-forming galaxies selected by their 870-micron emission in the GOODS-S field. By securing redshifts for most of the brighter sources, the authors model the redshift distribution and conclude that only a small fraction, less than or equal to 10 percent, are at redshifts above 4 and even fewer above 5. This points to a steep decline in the population of massive dusty galaxies in the early universe within the first 1.5 billion years. They also note that common photometric redshift techniques tend to overestimate the redshifts for these galaxies. The high spectroscopic completeness allows for a reliable view of how these galaxies are distributed across cosmic time.

Core claim

The authors use ALMA to obtain linescans that identify or confirm spectroscopic redshifts for 52 out of 75 DSFGs in their sample, achieving 97% completeness for sources brighter than 2.5 mJy. From this, they determine that the redshift distribution has at most 10% of the f870 ≳2 mJy DSFGs at z>4 and 2% at z>5. This reflects a steep decline in the abundance of massive dusty galaxies in the first 1.5 Gyr of the universe.

What carries the argument

The 870 μm flux-selected sample of DSFGs with ALMA linescan spectroscopy used to secure redshifts and model the distribution down to the confusion limit.

If this is right

The abundance of massive dusty galaxies drops sharply before the first 1.5 Gyr after the Big Bang.
Photometric redshift methods for DSFGs show outlier fractions above 20 percent and systematically overestimate high redshifts.
JWST spectroscopic coverage identifies nearly all targeted DSFGs with little dependence on flux or redshift.
The redshift distribution can be reliably modeled for the full sample down to near the 850-micron confusion limit.

Where Pith is reading between the lines

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

Galaxy formation models may need to delay the peak of dust-obscured star formation to match the low fraction of early massive dusty systems.
Cross-comparisons with selections at other submillimeter wavelengths could show whether 870-micron sources represent a distinct population.
Deeper ALMA observations of fainter sources would test if the same steep decline holds below the current flux limit.

Load-bearing premise

The 870 micron selection produces an unbiased sample of DSFGs whose redshifts can be securely identified without large incompleteness or errors at high redshift.

What would settle it

An independent complete redshift survey of the same field that finds substantially more than 10 percent of bright 870-micron DSFGs at z greater than 4.

Figures

Figures reproduced from arXiv: 2604.11879 by A. J. Barger, F. E. Bauer, L. L. Cowie, S. J. McKay.

**Figure 1.** Figure 1: RGB thumbnails of the ALMA sample. Unless otherwise noted, the filters are from JWST/NIRCam with R = F444W, G = F277W, and B = F150W. Where sources are not covered in these bands, we use other available HST and/or JWST bands and give the filter names on the panels. Cutouts are 4′′ × 4 ′′; north is up and east is to the left. We label the source name, 870 𝜇m flux density, and redshift on each panel [PITH_F… view at source ↗

**Figure 2.** Figure 2: All confirmed > 4.5𝜎 line detections (source names given in upper left corner of each panel). We show the spectra (black line) and error per channel (gray dotted line and shaded region) binned to a resolution of 80 km s−1 . The Gaussian fit to the line is shown (dark purple curve) with channels within ±2𝜎 shaded blue for clarity. We mark the adopted frequency of the line (cyan line) and label the transitio… view at source ↗

**Figure 3.** Figure 3: Cumulative number of ALMA sources above a given 𝑓870 𝜇m versus 𝑓870 𝜇m. We show the full sample (gray solid line), sources with secure or tentative spec-zs (light blue dotted line), and sources with secure spec-zs (dark blue dashed line). Above 𝑓870 𝜇m = 2.5 mJy (green dashed line), the spectroscopic completeness reaches 72% (97%) for secure (tentative) spec-zs. Our ALMA survey is flux-complete above 𝑓870… view at source ↗

**Figure 4.** Figure 4: Difference between photo-zs/grism-zs and spec-zs (Δ𝑧 ≡ 𝑧phot − 𝑧spec) vs. 𝑧spec for sources with both to quantify the effectiveness of various photo-z estimation methods. Each panel indicates the method and the photometric data used for that photo-z estimate, as well as the normalized median absolute deviation, 𝜎NMAD; the outlier fraction, 𝜂; and the bias for that method. Tentative spec-zs that we inferred… view at source ↗

**Figure 5.** Figure 5: Absolute difference (|Δ𝑧phot|) between photo-z/grism-z estimates vs. 𝐾𝑠 magnitude (top) and vs. spec-z, where available (bottom), for DSFGs with multiple photo-z/grism-z estimates. We use the S16 photo-zs as a reference, plotting the difference between these and the E26 photo-zs (purple stars), the MAGPHYS photo-zs (cyan diamonds), and the HST grism-zs (M16; red squares). Objects with |Δ𝑧phot| > 3, which … view at source ↗

**Figure 8.** Figure 8: Redshift distribution (scaled to unit probability) of the flux-complete ( 𝑓870 𝜇m ≥ 2.25 mJy) subset (gray histogram) and best-fit log-normal distribution with 5%–95% credible interval (blue curve and shaded region). For comparison, we show the distribution of redshifts (mainly MAGPHYS photo-zs) for the 707 AS2UDS DSFGs (dashed histogram; U. Dudzevici¯ut ˇ e et al. ˙ 2020). We mark several known overdensi… view at source ↗

**Figure 9.** Figure 9: Center: Best redshift vs. 𝑓870 𝜇m for our DSFG sample. We plot photo-zs (gray), tentative spec-zs identified from a single line (light blue), and secure spec-zs (dark blue). We show the best fit line to our data (black solid line), which agrees well with the best-fit relation from the 708 AS2UDS DSFGs (S. M. Stach et al. 2019; orange dashed-dotted line) and the spectroscopic sample of C.-C. Chen et al. (20… view at source ↗

**Figure 10.** Figure 10: 𝑓870 𝜇m vs. 𝐾𝑠 magnitude for our DSFG sample, showing the different methods used to obtain redshifts (see legend). Sources may have redshifts from multiple instruments, as indicated by the stacked symbols. At 𝑓870 𝜇m > 2.5 mJy (green dotted line and shaded region), the secure (tentative) spectroscopic completeness is 72% (97%); below, it is 40% (49%). We also mark 𝐾𝑠 = 23 mag (gray line) for illustrative … view at source ↗

**Figure 11.** Figure 11: Redshift distribution from [PITH_FULL_IMAGE:figures/full_fig_p016_11.png] view at source ↗

read the original abstract

We present an Atacama Large Millimeter/submillimeter Array (ALMA) spectroscopic follow-up survey of an 870 $\mu$m-selected sample of dusty star-forming galaxies (DSFGs) in the GOODS-S field. We use these linescans to identify or confirm spectroscopic redshifts (spec-zs) for 20 sources. Including spec-zs from the literature, there are now secure or tentative spec-zs for 52 out of 75 DSFGs (69%). At $f_{870}>2.5$ mJy, the sample is 97% spectroscopically complete, allowing us to model the full DSFG redshift distribution down to nearly the confusion limit for a 15-m telescope at 850 $\mu$m. This is the highest completeness for an unbiased sample at this flux limit to date. We find that nearly all of the DSFGs in our sample that were targeted with JWST/NIRSpec were spectroscopically identified, without much dependence on near-infrared or submillimeter flux or redshift. However, only 29% of our sample have JWST spectroscopic coverage. We use the spec-zs to evaluate various photometric redshift (photo-z) estimates, finding that all methods exhibit an outlier fraction of at least $>20$%. Nearly all of the photo-z methods tend to overshoot the redshifts, leading to overestimates of the number of DSFGs at high redshift ($z>4$). Our results suggest that $\lesssim10$% of $f_{870} \gtrsim 2$ mJy DSFGs lie at $z>4$ and $\lesssim2$% lie at $z>5$, reflecting a steep decline in the abundance of massive dusty galaxies in the first 1.5 Gyr.

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

This ALMA survey reaches 97% spectroscopic completeness for f870>2.5 mJy DSFGs and uses new linescans plus literature redshifts to cap the z>4 fraction at ≲10%.

read the letter

This paper's main contribution is a high-completeness spectroscopic redshift survey of 870 μm-selected dusty star-forming galaxies in GOODS-S. With ALMA linescans they secure redshifts for 20 new sources, bringing the total to 52 out of 75 with secure or tentative spec-zs. At fluxes above 2.5 mJy the sample reaches 97% completeness, which is the highest reported for an unbiased selection at this depth. They use this to argue that at most 10% of such DSFGs are at z>4 and 2% at z>5. The work stands out for directly comparing multiple photo-z methods on the identical objects and showing that all of them have outlier rates above 20%, with a tendency to push redshifts too high. This explains discrepancies with earlier estimates and supports a steeper drop-off in massive dusty galaxies at early times. The JWST cross-check, where almost all targeted sources were identified, adds confidence even though coverage is limited to 29% of the sample. The potential weak points are in the details of line identification and the exact construction of the redshift distribution model, which the abstract does not spell out. However, the full manuscript apparently provides those without internal contradictions or selection biases that would invalidate the upper limits. The numbers are not derived from fitted parameters on the same data, so there's little circularity. This is a solid data paper for anyone studying the redshift distribution of DSFGs or the early universe's star-forming population. It gives updated constraints that models of galaxy assembly and dust production will need to accommodate. The evidence from the spec-zs looks reliable enough to deserve referee attention. I would recommend sending it for peer review. The completeness and photo-z comparison are worth the time to verify.

Referee Report

2 major / 2 minor

Summary. The paper presents ALMA 870 μm linescan follow-up of a sample of 75 DSFGs in GOODS-S, securing spectroscopic redshifts for 52 sources (69% overall) including 20 new identifications from the linescans. At f870 > 2.5 mJy the sample reaches 97% spectroscopic completeness, enabling modeling of the redshift distribution down to near the confusion limit. The work evaluates multiple photo-z methods (all showing >20% outlier fractions that tend to overestimate z) and concludes that ≲10% of f870 ≳2 mJy DSFGs lie at z>4 and ≲2% at z>5.

Significance. If the identification security and completeness hold, the result supplies one of the highest-completeness spectroscopic constraints on the high-redshift tail of an unbiased 870 μm-selected DSFG sample. The direct use of ALMA linescans plus literature cross-matches to derive the low z>4 and z>5 fractions, rather than relying on fitted parameters from the same data, strengthens the claim of a steep decline in massive dusty galaxies within the first 1.5 Gyr.

major comments (2)

[ALMA linescan results and redshift assignments] The criteria used to identify lines in the ALMA linescans, assign secure versus tentative redshifts, and quantify the associated error budget are not described in sufficient detail. These choices directly underpin the reported 97% completeness at f870 > 2.5 mJy and the resulting upper limits on the high-z fractions.
[Redshift distribution modeling] The construction of the redshift distribution model, including any selection-function corrections, handling of sources below the confusion limit, and propagation of identification uncertainties, is not explicitly laid out. This modeling step is load-bearing for the claimed ≲10% (z>4) and ≲2% (z>5) fractions.

minor comments (2)

[JWST comparison] The statement that 'nearly all' JWST/NIRSpec-targeted sources were identified should be quantified with the exact fraction and any dependence on flux or redshift.
[Data presentation] Figure or table presenting the full list of redshift assignments (secure/tentative, method, and quality flags) would improve traceability of the completeness calculation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of our work and for the constructive comments, which help improve the clarity of the manuscript. We address each major comment below and will revise the paper accordingly to provide the requested details.

read point-by-point responses

Referee: [ALMA linescan results and redshift assignments] The criteria used to identify lines in the ALMA linescans, assign secure versus tentative redshifts, and quantify the associated error budget are not described in sufficient detail. These choices directly underpin the reported 97% completeness at f870 > 2.5 mJy and the resulting upper limits on the high-z fractions.

Authors: We agree that additional detail is needed on the line identification process. In the revised manuscript we will insert a new subsection (Section 3.2) that explicitly states the S/N thresholds applied to candidate lines, the requirement of either a second line or a secure cross-match to literature data for a 'secure' redshift versus a single-line 'tentative' assignment, and the Monte Carlo procedure used to propagate the linescan noise properties into the final redshift uncertainties. These additions will directly support the quoted 97% completeness figure at f870 > 2.5 mJy. revision: yes
Referee: [Redshift distribution modeling] The construction of the redshift distribution model, including any selection-function corrections, handling of sources below the confusion limit, and propagation of identification uncertainties, is not explicitly laid out. This modeling step is load-bearing for the claimed ≲10% (z>4) and ≲2% (z>5) fractions.

Authors: We accept that the modeling procedure requires more explicit description. The revised text will expand Section 4 to detail (i) the flux-dependent selection function derived from the ALMA primary beam and noise maps, (ii) the treatment of the 23 sources below the nominal confusion limit by restricting the high-completeness analysis to the f870 > 2.5 mJy subsample while retaining the full sample for the overall distribution, and (iii) the bootstrap/Monte Carlo resampling that folds in both secure/tentative redshift assignments and photo-z outlier rates to produce the final upper limits. These clarifications will make the derivation of the ≲10% (z>4) and ≲2% (z>5) fractions fully reproducible. revision: yes

Circularity Check

0 steps flagged

No significant circularity; empirical fractions from direct spec-z identifications

full rationale

The paper's central results on the redshift distribution and high-z fractions (≲10% at z>4, ≲2% at z>5 for f870 ≳2 mJy DSFGs) are obtained by direct counting of secure spectroscopic redshifts from ALMA linescans (20 new) plus literature (total 52/75, 97% complete above 2.5 mJy). No equations, fitted parameters, or self-citations reduce these observed fractions to quantities defined from the same data; the completeness assessment and photo-z evaluation are separate and do not loop back into the reported upper limits. The derivation chain is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that the 870 μm flux-limited sample is representative and that ALMA linescans yield reliable redshifts; no free parameters or invented entities are described in the abstract.

axioms (1)

domain assumption The 870 μm selection provides an unbiased sample of DSFGs whose redshift distribution can be modeled to the confusion limit.
Invoked to claim the sample allows modeling of the full DSFG redshift distribution.

pith-pipeline@v0.9.0 · 5651 in / 1469 out tokens · 47694 ms · 2026-05-10T14:55:29.083762+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

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

Extending the ALMA survey of the SCUBA-2 CLS UDS field: Tracing the obscured formation of spheroids across z~1-4
astro-ph.GA 2026-04 unverdicted novelty 5.0

Faint 870-um galaxies at z<~2.5 show extended star-forming discs with lower gas fractions while brighter higher-z systems are compact unstable starbursts, marking a shift in obscured spheroid formation.

Reference graph

Works this paper leans on

3 extracted references · 1 canonical work pages · cited by 1 Pith paper

[1]

, " * write output.state after.block = add.period write newline

ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...
[2]

write newline

" write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...
[3]

The connection between dusty star-forming galaxies and the first massive quenched galaxies

thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

work page doi:10.17909/8tdj-8n28 2017