Homogeneous measurements of proximity zone sizes for 59 quasars in the Epoch of Reionization

Anna-Christina Eilers; Bram P. Venemans; Elia Pizzati; Joseph F. Hennawi; Silvia Onorato

arxiv: 2505.09676 · v2 · submitted 2025-05-14 · 🌌 astro-ph.GA · astro-ph.CO

Homogeneous measurements of proximity zone sizes for 59 quasars in the Epoch of Reionization

Silvia Onorato , Joseph F. Hennawi , Elia Pizzati , Bram P. Venemans , Anna-Christina Eilers This is my paper

Pith reviewed 2026-05-22 15:19 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.CO

keywords proximity zonesquasarsepoch of reionizationintergalactic mediumhigh-redshiftluminosity dependenceredshift evolutionquasar lifetimes

0 comments

The pith

Proximity zone sizes around high-redshift quasars shrink more steeply with redshift than earlier studies found.

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

This paper measures proximity zone sizes for 59 quasars at redshifts from 5.77 to 7.54, including 15 measured for the first time. It fits how these sizes scale with both quasar brightness and redshift across a combined sample of 100 objects, finding the luminosity trend matches theoretical expectations while the redshift trend is steeper. The work also flags 13 quasars with unusually small zones and attributes them to active lifetimes shorter than about 10,000 years once dense gas absorbers are ruled out. A reader would care because these zones show how early supermassive black holes ionize their surroundings and how long those black holes have been shining during the epoch of reionization.

Core claim

The central claim is that homogeneous measurements of proximity zone sizes R_p for 59 quasars spanning 5.77 ≤ z ≤ 7.54 yield a bivariate power-law fit in which R_p scales with luminosity as 10 to the power of -0.4 M_1450 / 2.87, agreeing with models, but scales with redshift as (1 + z) to the power of -2.44, steeper than previous results. Thirteen quasars show small R_p values best explained by quasar lifetimes t_Q less than or equal to 10^4 years after associated dense absorbers are excluded as the cause.

What carries the argument

Proximity zone size R_p, the distance from the quasar at which Lyman-alpha transmission falls below a fixed threshold, together with the bivariate power-law model fitted jointly to redshift and absolute magnitude M_1450.

If this is right

The steeper redshift evolution implies that the neutral fraction of the intergalactic medium or related environmental factors must change more rapidly near redshift 6 than earlier measurements suggested.
Short quasar lifetimes of order 10,000 years or less for some objects limit the total ionizing photons they can contribute during reionization.
Scatter in R_p at fixed luminosity and redshift is dominated by differences in quasar lifetimes, local density, and ionization state rather than measurement noise alone.
Homogeneous measurements across the sample reduce systematic offsets that previously complicated comparisons between different quasar studies.

Where Pith is reading between the lines

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

If short lifetimes turn out to be typical, models of supermassive black hole assembly at early times would need to incorporate even briefer active phases or higher duty-cycle bursts.
Higher-resolution spectra or larger samples could separate lifetime effects from local density variations by correlating small zones with other observables such as metal lines.
These results suggest that quasar-driven ionization bubbles grow more slowly or start later than assumed in many reionization simulations.

Load-bearing premise

The small proximity zones are produced by genuinely short quasar lifetimes rather than unmodeled intergalactic density fluctuations or systematic errors in the residuals from the bivariate fit.

What would settle it

Independent measurements of quasar activity timescales, such as through multi-epoch variability monitoring or direct mapping of local neutral hydrogen patches, that show the 13 small-zone objects have been active for longer than 10,000 years would falsify the short-lifetime explanation.

read the original abstract

The overionized regions surrounding high-redshift quasars, known as proximity zones, provide a window into the interaction between supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). We present new homogeneous measurements of proximity zone sizes ($R_{\mathrm{p}}$) for a sample of $59$ quasars spanning redshifts $5.77 \leq z \leq 7.54$ (median $z = 6.59$). For $15$ of these sources, we measure $R_{\mathrm{p}}$ for the first time. The quasars in our catalog have absolute magnitudes at rest-frame $1450$ \r{A} in the range $-29.13 \leq M_{1450} \leq -25.20$ (median $M_{1450} \simeq -26.49$), providing one of the most extensive data sets for exploring $R_{\mathrm{p}}$ at these epochs. The distribution of $R_{\mathrm{p}}$ values shows a large scatter at fixed redshift and luminosity, likely reflecting variations in quasar lifetimes ($t_{\mathrm{Q}}$), IGM density fluctuations, and IGM neutral fraction. We fit a bivariate power-law model to a large sample of $100$ objects to study the dependence of $R_{\mathrm{p}}$ with both $M_{1450}$ and $z$: we find that the evolution of $R_{\mathrm{p}}$ with luminosity is in agreement with the models ($R_{\mathrm{p}} \propto 10^{-0.4 M_{1450}/2.87}$), while the evolution of $R_{\mathrm{p}}$ with $z$ is steeper than previous works ($R_{\mathrm{p}} \propto (1+z)^{-2.44}$). We identify $13$ quasars with small proximity zone size, defined using the residuals of our fit. In all cases, except for J2211$-$6320, we rule out the presence of associated dense absorbers that prematurely truncate $R_{\mathrm{p}}$, and suggest a short $t_{\mathrm{Q}}$ ($\lesssim 10^4$ yr) as a possible explanation for their small proximity zone sizes.

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 paper adds a solid homogeneous catalog of 59 proximity zone sizes with 15 new measurements and a steeper redshift trend than prior work, but the short-lifetime interpretation for the 13 small-zone objects rests on an unquantified assumption about residuals.

read the letter

This paper's main contribution is a set of homogeneous proximity zone size measurements for 59 quasars at z from 5.77 to 7.54, including 15 new ones, plus a bivariate power-law fit on 100 objects total. The redshift slope comes out steeper than previous work at (1+z)^-2.44, while the luminosity dependence matches models. They do a good job keeping the measurements consistent across the sample. That reduces method-to-method differences that have plagued earlier compilations. The scatter at fixed luminosity and redshift is shown clearly, and they flag 13 small zones after ruling out absorbers in most cases. The weak point is the leap to short quasar lifetimes for those 13. They use residuals from the fit to pick them out, but the paper does not show whether typical IGM density fluctuations at these redshifts could produce similar small zones. Without that comparison, the lifetime claim is plausible but not yet required by the data. Readers studying the epoch of reionization and early supermassive black holes will find the catalog and the new slope useful. It is not a big theoretical advance, but the observational result is solid enough to merit referee time. I would send it out for review.

Referee Report

2 major / 2 minor

Summary. The manuscript presents homogeneous measurements of proximity zone sizes R_p for 59 quasars at 5.77 ≤ z ≤ 7.54, including 15 new measurements. Combining these with literature data for a total of 100 objects, the authors fit a bivariate power-law model R_p(M_1450, z) and report a luminosity dependence R_p ∝ 10^{-0.4 M_1450 / 2.87} that matches models, together with a steeper redshift evolution R_p ∝ (1+z)^{-2.44} than found in prior work. They identify 13 sources whose R_p values lie significantly below the fitted relation; after ruling out dense absorbers for 12 of them, the authors interpret the small zones as evidence for short quasar lifetimes t_Q ≲ 10^4 yr.

Significance. If the central interpretation holds, the work supplies one of the largest homogeneous R_p catalogs at z > 5.8 and strengthens the case for a steep redshift dependence that may constrain the neutral fraction or IGM thermal state during reionization. The explicit bivariate fit and the homogeneous measurement approach are clear strengths that facilitate future comparisons. The identification of short-lifetime candidates, if robust against IGM fluctuations, would directly inform SMBH growth timescales at early epochs.

major comments (2)

[abstract and small proximity zones section] Abstract and section on small proximity zones: the 13 small-R_p candidates are defined via negative residuals from the bivariate power-law fit. The manuscript notes large scatter at fixed luminosity and redshift but does not quantify the expected distribution of residuals arising from IGM density fluctuations at z ≈ 6–7.5 and the relevant neutral fractions; without this comparison it remains possible that the observed tail is produced by stochastic IGM structure rather than short t_Q, which is load-bearing for the lifetime interpretation.
[results section on the bivariate fit] Results section on the bivariate fit: the paper does not report the uncertainties or covariance on the fitted indices (-2.44 for redshift and the luminosity index), nor does it describe the error propagation or fitting method (e.g., whether measurement errors on individual R_p are included). This limits assessment of whether the steeper redshift slope is statistically significant relative to previous determinations.

minor comments (2)

A table or figure panel explicitly listing the 13 small-zone objects, their measured R_p, residuals, and the specific tests used to rule out dense absorbers would improve clarity and reproducibility.
The definition of R_p (e.g., the precise transmission threshold and wavelength range) should be restated in the main text even if it follows a prior reference, to make the homogeneous aspect self-contained.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their detailed and constructive comments on our manuscript. We address each of the major comments below and outline the revisions we plan to make.

read point-by-point responses

Referee: Abstract and section on small proximity zones: the 13 small-R_p candidates are defined via negative residuals from the bivariate power-law fit. The manuscript notes large scatter at fixed luminosity and redshift but does not quantify the expected distribution of residuals arising from IGM density fluctuations at z ≈ 6–7.5 and the relevant neutral fractions; without this comparison it remains possible that the observed tail is produced by stochastic IGM structure rather than short t_Q, which is load-bearing for the lifetime interpretation.

Authors: We agree that a quantitative comparison of the observed residuals to those expected from IGM density fluctuations would strengthen the case for interpreting the smallest zones as short-lifetime candidates. The manuscript already states that the large scatter at fixed luminosity and redshift arises from a combination of quasar lifetime variations, IGM density fluctuations, and neutral fraction. For the 13 sources we have ruled out dense absorbers in all but one case. While performing new Monte Carlo simulations of IGM structure at these redshifts and neutral fractions is beyond the present scope, we will revise the relevant section to cite existing cosmological simulations that characterize the expected scatter in proximity zone sizes due to IGM fluctuations. This will help readers evaluate whether the most extreme negative residuals are more readily explained by short t_Q. revision: partial
Referee: Results section on the bivariate fit: the paper does not report the uncertainties or covariance on the fitted indices (-2.44 for redshift and the luminosity index), nor does it describe the error propagation or fitting method (e.g., whether measurement errors on individual R_p are included). This limits assessment of whether the steeper redshift slope is statistically significant relative to previous determinations.

Authors: We thank the referee for highlighting this omission. Reporting the uncertainties, covariance, and details of the fitting procedure is necessary to allow readers to judge the statistical significance of the steeper redshift dependence relative to earlier studies. We will revise the results section to describe the bivariate power-law fitting method in full, specify whether measurement uncertainties on individual R_p values were incorporated, and provide the uncertainties and covariance matrix for the fitted indices. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results are direct measurements and empirical fits.

full rationale

The paper reports homogeneous R_p measurements for 59 quasars (15 new), combines with prior data to reach 100 objects, and fits a bivariate power-law in M_1450 and z. The reported indices (R_p ∝ 10^{-0.4 M_1450 / 2.87} and R_p ∝ (1+z)^{-2.44}) are direct outputs of that fit to observed values. The 13 small-R_p candidates are flagged by negative residuals from the same fit; this is a standard outlier definition, not a self-definitional loop or a prediction that reduces to the input by construction. No load-bearing self-citation, uniqueness theorem, or ansatz smuggling is present in the provided text. The derivation chain consists of measurement, fitting, and residual-based selection, all of which remain externally falsifiable against independent IGM simulations or larger samples.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The analysis rests on standard assumptions about quasar spectra and IGM transmission but introduces no new free parameters beyond the two power-law indices fitted to data; no invented entities are postulated.

free parameters (2)

luminosity index
The exponent -0.4 M_1450 / 2.87 is obtained by fitting the combined sample of 100 objects.
redshift index
The exponent -2.44 is obtained by fitting the combined sample of 100 objects.

axioms (1)

domain assumption Proximity zone size is set by the quasar lifetime, local IGM density, and global neutral fraction.
Invoked when interpreting scatter and small-zone candidates.

pith-pipeline@v0.9.0 · 5975 in / 1498 out tokens · 32048 ms · 2026-05-22T15:19:00.803634+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/Foundation/AbsoluteFloorClosure.lean reality_from_one_distinction unclear

?

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

We fit a bivariate power-law model to a large sample of 100 objects... R_p ∝ (1+z)^(-2.44) ... identify 13 quasars with small proximity zone size, defined using the residuals of our fit.
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

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

The distribution of R_p values shows a large scatter at fixed redshift and luminosity, likely reflecting variations in quasar lifetimes (t_Q), IGM density fluctuations...

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 1 Pith paper

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

QUEST (Quasar Unsupervised Encoder and Synthesis Tool): A machine learning framework to generate quasar spectra
astro-ph.GA 2025-10 unverdicted novelty 5.0

A variational autoencoder learns to generate and reconstruct quasar spectra from SDSS data, reproducing median and variance properties while enabling photometry synthesis and absorption-line interpolation without ad-h...