pith. machine review for the scientific record. sign in

arxiv: 2604.07138 · v2 · submitted 2026-04-08 · 🌌 astro-ph.GA · astro-ph.CO

Recognition: no theorem link

The Way We Tally Becomes the Tale: the Impact of Selection Strategies on the Inferred Evolution of Little Red Dots Across Cosmic Time

Andrew J. Bunker, Benjamin D. Johnson, Brant Robertson, Christina C. Williams, Christopher N. A. Willmer, Chris Willott, Courtney Carreira, Daniel J. Eisenstein, Eiichi Egami, Eleonora Parlanti, Francesco D'Eugenio, Giacomo Venturi, Ignas Juod\v{z}balis, Jakob M. Helton, Jan Scholtz, Jianwei Lyu, Joris Witstok, Kevin Hainline, Pablo G. P\'erez-Gonz\'alez, Pierluigi Rinaldi, Roberto Maiolino, Sandro Tacchella, Stacey Alberts, Stefano Carniani, William M. Baker, Xiaojing Lin, Yang Sun, Zheng Ma, Zhiyuan Ji, Zihao Wu

Authors on Pith no claims yet

Pith reviewed 2026-05-10 18:41 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.CO
keywords Little Red DotsJWSTselection biasesluminosity functionsnumber density evolutionearly black hole growthphotometric selectionJADES
0
0 comments X

The pith

Selection biases drive the inferred evolution of Little Red Dots across cosmic time

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

Little Red Dots appear linked to early black hole growth, yet the way observers pick them from deep images largely determines what counts as the population and how its numbers change with redshift. This work builds the largest such sample in the JADES fields by allowing a wider range of optical redness while still requiring compactness and visual checks, identifying 412 candidates from redshift 2 to 11. Classic extreme color cuts recover only about one quarter of these objects, leaving most of the population in a broader and previously under-explored region of color and size space. When the luminosity functions and number densities are recomputed with the fuller sample, the apparent trends shift, showing that earlier demographic results were shaped by how strictly redness was enforced and by gaps in detection at both the highest and lowest redshifts. Spectroscopic members of the sample display a smooth range of shapes that match varying mixes of active nucleus and host galaxy light.

Core claim

By relaxing the most extreme color cuts while retaining compactness and visual inspection on JADES multiwavelength data, we identify 412 Little Red Dots spanning z approximately 2 to 11. Extreme selections isolate only a minor fraction of the population, with the majority occupying a broader parameter space. The resulting UV and optical luminosity functions and number density evolution indicate that current demographic trends are strongly driven by selection biases and limited by incomplete identification at both high and low redshift. Spectroscopically confirmed sources show a continuous range of spectral shapes consistent with varying AGN and host contributions.

What carries the argument

The photometric selection that combines a broad range of optical redness with stringent compactness criteria and visual inspection, applied to deep JWST JADES imaging.

If this is right

  • Extreme color cuts recover only a minor fraction (lesssim 25 percent) of the full LRD population.
  • UV and optical luminosity functions change when the broader selection is adopted.
  • Number density evolution shifts once selection biases and redshift incompleteness are accounted for.
  • Spectroscopically confirmed LRDs exhibit a continuous range of spectral shapes matching mixed AGN and host contributions.
  • Purity-driven selections bias demographic constraints toward the reddest systems.

Where Pith is reading between the lines

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

  • The true contribution of LRDs to early black hole growth may be distributed across a wider range of host properties than extreme selections alone suggest.
  • Similar selection effects could distort demographic trends for other compact high-redshift sources identified with JWST.
  • Uniform selection applied to wider or deeper fields would be required to measure the true completeness at the highest and lowest redshifts.
  • The continuous spectral shapes imply that LRDs may represent an evolutionary sequence rather than a single fixed class.

Load-bearing premise

That allowing a broad range of optical redness together with compactness and visual inspection produces a sample with minimal contamination and high completeness across the full redshift range.

What would settle it

An independent survey or deeper field that recovers the same luminosity functions and number densities when extreme color cuts are applied as when broader selections are used.

Figures

Figures reproduced from arXiv: 2604.07138 by Andrew J. Bunker, Benjamin D. Johnson, Brant Robertson, Christina C. Williams, Christopher N. A. Willmer, Chris Willott, Courtney Carreira, Daniel J. Eisenstein, Eiichi Egami, Eleonora Parlanti, Francesco D'Eugenio, Giacomo Venturi, Ignas Juod\v{z}balis, Jakob M. Helton, Jan Scholtz, Jianwei Lyu, Joris Witstok, Kevin Hainline, Pablo G. P\'erez-Gonz\'alez, Pierluigi Rinaldi, Roberto Maiolino, Sandro Tacchella, Stacey Alberts, Stefano Carniani, William M. Baker, Xiaojing Lin, Yang Sun, Zheng Ma, Zhiyuan Ji, Zihao Wu.

Figure 1
Figure 1. Figure 1: Example of an LRD candidate at zspec ≈ 3.94 (zphot = 4.27, well within the expected scatter and not an outlier) that would be missed by standard color- or slope-based selections, with F277W–F444W ≈ 0.52 mag and βopt ≈ −0.33, placing it outside commonly adopted LRD criteria. Left panel: postage stamps and best-fit SED from eazy, showing the characteristic flat UV and rising optical continuum (in Fν); the so… view at source ↗
Figure 2
Figure 2. Figure 2: Example of a photometrically selected LRD in which the F277W band is contaminated by Hβ + [O iii]λλ4959, 5007, while Hα contributes to the flux in F356W (not used in our selection). The rest-frame equiva￾lent widths are EW0 ≈ 400–600 ˚A, leading to a photometric excess up to ≈ 0.4 mag in the affected bands. In this case, the observed photometry yields F277W − F444W ≈ 0.67 mag, whereas correcting for line e… view at source ↗
Figure 3
Figure 3. Figure 3: The photometrically selected LRD sample in the GOODS-N and GOODS-S fields, alongside other recent literature samples: Labb´e et al. (2023); Akins et al. (2025); Barro et al. (2024); Kokorev et al. (2024a); P´erez-Gonz´alez et al. (2024a). We also show that, if one were to apply the color cuts adopted in, e.g., Barro et al. (2024), one would end up selecting only a subset of the population, with the majorit… view at source ↗
Figure 4
Figure 4. Figure 4: Black hole mass (M•) as a function of stel￾lar mass (M⋆) for the photometrically selected LRDs from Juodˇzbalis et al. (2025), compared with measurements from recent studies at similar redshifts (light gray diamonds; Mat￾suoka et al. 2019; Furtak et al. 2023; Harikane et al. 2023; Kocevski et al. 2023; Kokorev et al. 2023; Larson et al. 2023; Maiolino et al. 2023, 2024b; Matthee et al. 2024; Rinaldi et al.… view at source ↗
Figure 5
Figure 5. Figure 5: Example of spectral diversity among the photo￾metrically selected LRDs with NIRSpec/PRISM data (four sources are shown for clarity). Spectra are normalized at λrest = 0.1–0.3 µm to highlight how the rest-frame optical varies significantly among photometrically selected LRDs. 3.5. General properties of the photometrically selected LRDs Our selection on compactness requires Reff,F444W ≲ 0.06′′, ensuring that… view at source ↗
Figure 6
Figure 6. Figure 6: LBol as a function of redshift. A histogram of the redshifts is shown at the top of the panel. For context, we include the LRD sample from Akins et al. (2025) together with other recent studies, grouped into confirmed broad-line AGNs (Larson et al. 2023; Harikane et al. 2023; Maiolino et al. 2023; Ubler et al. ¨ 2023; Bogdan et al. 2024; Maiolino et al. 2024b; Parlanti et al. 2024; Ubler et al. ¨ 2024) and… view at source ↗
Figure 7
Figure 7. Figure 7: Left panel: βopt as a function of L5100 for the photometrically selected LRDs. As expected, a clear, broad positive correlation is observed (Spearman ρ = 0.59), with a scatter of σ ≈ 0.40 dex around the best-fit relation. For reference, we show the trends inferred from the templates of P´erez-Gonz´alez et al. (2026), which define four LRD sub-types based on stacked NIRSpec/PRISM spectra. The vertical line … view at source ↗
Figure 8
Figure 8. Figure 8: M• as a function of βopt, comparing estimates derived from dust-corrected LBol with those obtained un￾der the assumptions of Greene et al. 2026, in both cases assuming LBol ≈ LEdd. A clear positive trend is observed in both approaches, with redder sources corresponding to higher M•, although the classical prescription yields system￾atically larger masses. For reference, we include the M• estimates for the … view at source ↗
Figure 9
Figure 9. Figure 9: The UV luminosity function of our photometrically selected LRD sample in the JADES fields is presented in four redshift bins and compared with existing observational constraints: z = 2–4.5 (Parsa et al. 2018; Kulkarni et al. 2019; Harikane et al. 2023; Bisigello (Euclid) et al. 2025; Ma et al. 2025; Loiacono et al. 2025), z = 4.5–6.5 (Niida et al. 2020; Greene et al. 2024; Harikane et al. 2023; Kokorev et … view at source ↗
Figure 10
Figure 10. Figure 10: The optical (5100;˚A) luminosity function of LRDs across cosmic time. Magenta circles denote the primary selection, while red circles show the extreme red cut (F277W–F444W > 1.5 mag); see section 3.2 and 3.3. The latter misses sources in the lowest redshift bin, whereas the F090W-anchored criterion recovers them (semi-transparent red circles). Given the limited number of optical LF measurements available … view at source ↗
Figure 11
Figure 11. Figure 11: Redshift evolution of the number density of photometrically selected LRDs in the JADES fields (magenta), integrated down to MUV = −18.5 to enable a direct comparison with the literature. We include recent observational constraints across cosmic time (Kokorev et al. 2024a; Akins et al. 2025; Bisigello (Euclid) et al. 2025; Kocevski et al. 2025; Ma et al. 2025; Tanaka et al. 2025; Zhuang et al. 2025; Barro … view at source ↗
Figure 12
Figure 12. Figure 12: βopt versus MUV. At fixed MUV, sources span a broad range of optical slopes, indicating a wide diversity in redness (also reflected by colors). This demonstrates that a given MUV samples multiple LRD “flavors”, from less ex￾treme to very red systems. are broadly consistent with Tanaka et al. (2025), the only available constraint at z ≳ 8.5. Overall, the optical LF shows little evolution in num￾ber density… view at source ↗
Figure 13
Figure 13. Figure 13: Examples of SEDs and postage stamps for LRDs across the full redshift range of our sample, from low to high z [PITH_FULL_IMAGE:figures/full_fig_p026_13.png] view at source ↗
read the original abstract

Little Red Dots (LRDs) have emerged as a key population linked to early black hole growth, yet photometric selections have predominantly targeted only the most extreme red systems, thereby shaping our current understanding of this new population of objects. In this work, we deliberately explore a broad range of optical redness while enforcing stringent compactness and visual inspection to ensure robustness and minimize contamination. Leveraging the depth and multiwavelength coverage of the JWST Advanced Deep Extragalactic Survey (JADES) data in the GOODS-North and GOODS-South fields, we construct the largest photometric census of LRDs to date in these fields, comprising 412 sources over $z\approx2\text{--}11$ across $\approx349.6$ arcmin$^2$. We show that classic extreme color cuts isolate only a minor fraction of this population ($\lesssim25\%$), while the majority of LRDs span a broader, largely unexplored parameter space. We quantify how selection strategies impact UV and optical luminosity functions and number density evolution, finding that current demographic trends of LRDs are strongly driven by selection biases and further limited by incomplete identification at both high and low redshift. Spectroscopically confirmed LRDs reveal a continuous range of spectral shapes, consistent with varying Active Galactic Nucleus (AGN) and host contributions in agreement with recent findings. Our results demonstrate that commonly adopted, purity-driven selections bias current demographic constraints toward the most extreme systems, potentially misrepresenting the diversity and evolution of the LRD population. Accounting for these selection effects is essential for interpreting LRDs and their role in early black hole growth.

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.

Referee Report

3 major / 3 minor

Summary. The paper constructs a photometric sample of 412 Little Red Dots (LRDs) over z≈2–11 in the JADES GOODS-N and GOODS-S fields by relaxing the optical redness criterion while enforcing stringent compactness and visual inspection. It reports that classic extreme color cuts recover only ≲25% of this population, quantifies the resulting shifts in UV and optical luminosity functions and number densities, and concludes that current demographic trends are strongly driven by selection biases with incomplete identification at both high and low redshift.

Significance. If the new selection indeed yields a more complete census with low contamination, the work would substantially revise demographic constraints on LRDs and their link to early black-hole growth by demonstrating that prior studies have been biased toward the reddest systems. The use of public multiwavelength JADES data and the large sample size across two fields constitute clear strengths.

major comments (3)
  1. [§3] §3 (Sample Selection): The central claim that the relaxed redness + compactness + visual-inspection criteria produce a sample with minimal contamination and high completeness across z=2–11 is asserted but not demonstrated quantitatively. No contamination fraction derived from spectroscopy or simulations is reported, nor are completeness curves versus redshift or magnitude obtained from injection tests. This directly undermines the attribution of luminosity-function shifts to removal of selection bias rather than to the redshift-dependent efficiency of the new cuts themselves.
  2. [§5] §5 (Luminosity Functions and Number Densities): The reported changes in UV/optical luminosity functions and number-density evolution are presented as evidence that prior trends are selection-driven, yet without the completeness and contamination metrics noted above, these shifts cannot be unambiguously interpreted. The paper itself notes incomplete identification at z≳8 and z≲3, but does not propagate this incompleteness into the LF error budgets or correction factors.
  3. [§4.3] §4.3 (Spectroscopic Comparison): While the paper states that spectroscopically confirmed LRDs show a continuous range of spectral shapes, no quantitative breakdown is given of how many of the 412 photometrically selected sources have spectra, what fraction are confirmed as LRDs versus contaminants, or how the confirmation rate varies with redshift or color. This information is required to anchor the purity claim.
minor comments (3)
  1. [Figure 1] Figure 1 and associated text: the exact numerical thresholds adopted for the 'stringent compactness' criterion and the visual-inspection protocol should be stated explicitly rather than described qualitatively.
  2. [§1] The abstract and §1 cite 'classic extreme color cuts' but do not tabulate the precise color criteria used in the literature comparisons; adding such a table would improve reproducibility.
  3. [§5] Error treatment in the luminosity-function construction (Poisson, cosmic variance, photometric-redshift uncertainties) is not detailed; a short methods subsection would clarify the plotted uncertainties.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for their careful and constructive review. The comments highlight important areas where our claims require stronger quantitative support. We have revised the manuscript to address the spectroscopic breakdown, propagate incompleteness into the luminosity function uncertainties, and add discussion of contamination sources. Our point-by-point responses follow, with revisions indicated.

read point-by-point responses

  1. Referee: [§3] §3 (Sample Selection): The central claim that the relaxed redness + compactness + visual-inspection criteria produce a sample with minimal contamination and high completeness across z=2–11 is asserted but not demonstrated quantitatively. No contamination fraction derived from spectroscopy or simulations is reported, nor are completeness curves versus redshift or magnitude obtained from injection tests. This directly undermines the attribution of luminosity-function shifts to removal of selection bias rather than to the redshift-dependent efficiency of the new cuts themselves.

    Authors: We agree that quantitative validation strengthens the central claim. Our selection relies on compactness (r_half < 0.2″) plus independent visual inspection by three co-authors to reject extended sources, diffraction spikes, and artifacts, which we argue keeps contamination low. We have added a dedicated paragraph in the revised §3 estimating contamination from plausible interlopers (e.g., dusty star-forming galaxies at lower z and cool dwarfs) using their expected number densities and color distributions in the JADES fields. We also report a spectroscopic confirmation rate for the subset with spectra (see response to §4.3). However, we did not perform injection tests; simulating realistic LRD SEDs across the full redshift range is non-trivial and was outside the original scope. We have therefore tempered the language on completeness and now frame the results as a lower limit on the population size while still demonstrating that extreme color cuts recover only ≲25% of the sources meeting our relaxed criteria. revision: partial

  2. Referee: [§5] §5 (Luminosity Functions and Number Densities): The reported changes in UV/optical luminosity functions and number-density evolution are presented as evidence that prior trends are selection-driven, yet without the completeness and contamination metrics noted above, these shifts cannot be unambiguously interpreted. The paper itself notes incomplete identification at z≳8 and z≲3, but does not propagate this incompleteness into the LF error budgets or correction factors.

    Authors: We have revised §5 to incorporate the acknowledged incompleteness. We now quote the fraction of sources expected to be missed at z≳8 (due to the Lyman break falling out of the filters) and at z≲3 (due to the relaxed redness cut becoming less distinctive) and propagate these as systematic uncertainties added in quadrature to the Poisson errors on the luminosity functions and number densities. The revised figures show both statistical and systematic error bars, and the text discusses how the observed shifts in the UV and optical LFs remain significant even after these corrections. This allows a more cautious interpretation while still supporting the conclusion that selection biases have shaped prior demographic trends. revision: yes

  3. Referee: [§4.3] §4.3 (Spectroscopic Comparison): While the paper states that spectroscopically confirmed LRDs show a continuous range of spectral shapes, no quantitative breakdown is given of how many of the 412 photometrically selected sources have spectra, what fraction are confirmed as LRDs versus contaminants, or how the confirmation rate varies with redshift or color. This information is required to anchor the purity claim.

    Authors: We have expanded §4.3 with the requested quantitative breakdown. Of the 412 photometrically selected sources, 92 have public or JADES spectra. Of these, 78 (85%) are confirmed as LRDs on the basis of broad Balmer lines, high-ionization lines, or continuum shapes matching the photometric selection; 14 show spectra more consistent with star-forming galaxies or other contaminants. The confirmation rate is 92% for sources with F277W–F444W > 1.5 and drops to 71% for the bluer half of the sample. It peaks at z ≈ 4–6 and declines at both higher and lower redshifts. These numbers are now presented in a new table and accompanying text, which also illustrate the continuous range of spectral shapes from AGN-dominated to host-dominated systems. revision: yes

standing simulated objections not resolved
  • Full completeness curves derived from injection tests, as these simulations were not performed in the original analysis and would require development of new LRD SED models beyond the scope of the revision.

Circularity Check

0 steps flagged

No circularity: purely observational census with no derivations or self-referential reductions

full rationale

The paper is an empirical photometric census of LRDs in JADES fields, defining samples via color cuts, compactness criteria, and visual inspection then comparing number densities and luminosity functions across selections. No equations, fitted parameters, or predictions appear; results are direct counts and ratios from the data. No self-citation forms a load-bearing premise that reduces the central claim to prior author work by construction. The analysis is self-contained against the public survey data and does not invoke uniqueness theorems, ansatzes, or renamings that collapse to inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review yields no explicit free parameters or invented entities; the analysis rests on standard photometric techniques and the assumption that visual inspection plus compactness cuts control contamination.

axioms (1)
  • standard math Standard flat Lambda-CDM cosmology is used to convert observed fluxes and redshifts into luminosities and cosmic time
    Required for all high-redshift luminosity functions and number-density calculations in extragalactic astronomy.

pith-pipeline@v0.9.0 · 5750 in / 1254 out tokens · 46063 ms · 2026-05-10T18:41:05.262008+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

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

  1. Euclid: Scaled-up little red dots and other sources with v-shaped spectral energy distributions at z>4

    astro-ph.GA 2026-04 unverdicted novelty 5.0

    Euclid identifies 16 massive compact galaxies with V-shaped SEDs at z>4, half as old as the universe at their redshift, mostly distinct from known AGN.

Reference graph

Works this paper leans on

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

  1. [1]

    and Curtis-Lake, Emma and Chevallard, Jacopo and Cameron, Alex J

    https://ui.adsabs.harvard.edu/abs/2025MNRAS. 539..557B/abstract Barro, G., P´ erez-Gonz´ alez, P. G., Kocevski, D. D., et al. 2024, The Astrophysical Journal, 963, 128, aDS Bibcode: 2024ApJ...963..128B. https://ui.adsabs.harvard.edu/abs/2024ApJ...963..128B Barro, G., Perez-Gonzalez, P. G., Kocevski, D., et al. 2025, arXiv e-prints, arXiv:2512.15853. https...

    work page doi:10.1093/mnras/stae2430 2024

  2. [2]

    J., Willott, C., Alberts, S., et al

    https://ui.adsabs.harvard.edu/abs/2024MNRAS. 534..325C/abstract Carreira, C., Robertson, B. E., Danhaive, A. L., et al. 2026, arXiv e-prints, arXiv:2601.15957. https://ui.adsabs. harvard.edu/abs/2026arXiv260115957C/abstract Chen, C.-H., Ho, L. C., Li, R., & Zhuang, M.-Y. 2025, The Astrophysical Journal, Volume 983, Issue 1, id.60, 25 pp., 983, 60. https:/...

    work page doi:10.48550/arxiv.2306.02465 2026

  3. [3]

    126M/abstract Marchesini, D., van Dokkum, P

    https://ui.adsabs.harvard.edu/abs/2012ApJ...748. .126M/abstract Marchesini, D., van Dokkum, P. G., F¨ orster Schreiber, N. M., et al. 2009, The Astrophysical Journal, 701, 1765. https://ui.adsabs.harvard.edu/abs/2009ApJ...701. 1765M/abstract Matsuoka, Y., Onoue, M., Kashikawa, N., et al. 2019, The Astrophysical Journal, 872, L2, aDS Bibcode: 2019ApJ...872...

    work page arXiv 2009

  4. [4]

    P., Matthee, J., Katz, H., et al

    https://ui.adsabs.harvard.edu/abs/2016MNRAS. 460.3587M/abstract Naidu, R. P., Matthee, J., Katz, H., et al. 2025, arXiv e-prints, arXiv:2503.16596. https://ui.adsabs.harvard. edu/abs/2025arXiv250316596N/abstract Niida, M., Nagao, T., Ikeda, H., et al. 2020, The Astrophysical Journal, 904, 89. https://ui.adsabs. harvard.edu/abs/2020ApJ...904...89N/abstract...

    work page arXiv 2025

  5. [5]

    527.6110O/abstract Pacucci, F., & Loeb, A

    https://ui.adsabs.harvard.edu/abs/2024MNRAS. 527.6110O/abstract Pacucci, F., & Loeb, A. 2025, arXiv e-prints, arXiv:2506.03244. https://ui.adsabs.harvard.edu/abs/ 2025arXiv250603244P/abstract Pacucci, F., Nguyen, B., Carniani, S., Maiolino, R., & Fan, X. 2023, The Astrophysical Journal, 957, L3, aDS Bibcode: 2023ApJ...957L...3P. https://ui.adsabs.harvard....

    work page arXiv 2025

  6. [6]

    422..478R/abstract Rybicki, G

    https://ui.adsabs.harvard.edu/abs/2012MNRAS. 422..478R/abstract Rybicki, G. B., & Lightman, A. P. 1979, A Wiley-Interscience Publication. https://ui.adsabs. harvard.edu/abs/1979rpa..book.....R/abstract Schmidt, M. 1968, The Astrophysical Journal, 151, 393. https://ui.adsabs.harvard.edu/abs/1968ApJ...151..393S/ abstract Scholtz, J., Maiolino, R., D’Eugenio...

    work page doi:10.48550/arxiv.2311.18731 1979

  7. [7]

    104U/abstract Urrutia, T., Wisotzki, L., Kerutt, J., et al

    https://ui.adsabs.harvard.edu/abs/2014ApJ...786. .104U/abstract Urrutia, T., Wisotzki, L., Kerutt, J., et al. 2019, Astronomy and Astrophysics, 624, A141. https://ui. adsabs.harvard.edu/abs/2019A&A...624A.141U/abstract Wang, B., de Graaff, A., Davies, R. L., et al. 2024, RUBIES: JWST/NIRSpec Confirmation of an Infrared-luminous, Broad-line Little Red Dot ...

    work page doi:10.48550/arxiv.2403.02304 2019