Pith. sign in

REVIEW 3 major objections 6 minor 27 references

Detection of the cosmological evolution of the Doppler factor for blazars jets

T0 review · 3 major / 6 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Blazar jets are more strongly beamed at high redshift, according to gamma-ray spectra and optical variability of two blazar samples.

desk verdict A transparent measurement of redshift trends in Eb and tau_DRW that is being oversold as Doppler-factor evolution before the intrinsic-evolution degeneracy is tackled. read the letter →

arxiv 2501.15441 v1 pith:WWQQJT7H submitted 2025-01-26 astro-ph.HE

Dahai Yan , Haiyun Zhang , Jianeng Zhou , Lijuan Dong This is my paper · ORCID
classification astro-ph.HE
keywords blazarsDopplerfactorrelativisticjetscosmologicalevolutiongamma-rayspectradampedrandomwalkGaussianprocessregressionAGNcensus
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Blazar jets point nearly at Earth, so their observed brightness, photon energies, and variability timescales are all amplified by a Doppler factor $\delta$ set by the jet speed and viewing angle. This paper claims to detect, for the first time, that $\delta$ itself grows with cosmic epoch: $\delta\propto(1+z)^m$ with $m\approx 0.8$--$1$ on average. The claim comes from two independent datasets: the characteristic energies $E_b$ of the $\gamma$-ray spectra of 141 blazars give $m=0.81\pm0.12$, and the damping timescales of 89 optical light curves modeled as damped random walks give $m=1.09^{+0.25}_{-0.24}$. If the claim is right, high-redshift blazars are systematically more strongly beamed than nearby ones, which changes estimates of how common relativistic jets and obscured active galactic nuclei were in the early universe.

What carries the argument

The argument is carried by two redshift-dependent observables. The first is the characteristic energy $E_b$ of a log-parabolic gamma-ray spectrum, $N(E)\propto(E/E_b)^{-\alpha-b\log(E/E_b)}$, measured from 15-year Fermi-LAT spectra. The second is the damping timescale $\tau_{\rm DRW}$ of a damped random walk used in Gaussian-process fits to long-term optical light curves. Under relativistic beaming these transform as $E_b=E_b^{\rm int}\delta/(1+z)$ and $\tau_{\rm DRW}=\tau_{\rm int}(1+z)/\delta$. Writing $\delta=\delta_0(1+z)^m$ and taking logarithms turns both relations into linear regressions whose slopes are $(m-1)$ and $(1-m)$; a hierarchical Bayesian regression of $\log E_b$ and $\log\tau_{\rm DRW}$ against $\log(1+z)$ therefore yields $m$. The conversion from $m$ to the jet census uses $N_{\rm jetted}=2\Gamma^2 N_{\rm blazar}$ with $\Gamma\propto(1+z)^{0.8-1}$.

What would settle it

Take a complete, flux-limited blazar sample at $z\approx0.5$ and $z\approx3$ and measure the Doppler factor by an independent method such as the synchrotron self-absorption turnover or $\gamma\gamma$ opacity; if the median $\delta$ does not grow roughly as $(1+z)^{0.8}$, the claimed slope is an artifact of assuming constant intrinsic quantities or of selection effects.

Watch

Extended reading notes

Core claim

The paper reports a positive cosmological evolution of the Doppler factor in $\gamma$-ray bright blazars out to $z\approx3$. Using the observed relation $E_b=E_b^{\rm int}\delta/(1+z)$ for 141 blazars, the authors find $\log E_b$ grows with $\log(1+z)$ at a slope that implies $m=0.81\pm0.12$; using $\tau_{\rm DRW}=\tau_{\rm int}(1+z)/\delta$ for 89 blazars, they find $m=1.09^{+0.25}_{-0.24}$. They interpret the consistency as evidence that $\delta\propto(1+z)^{0.8-1}$ and that the jet Lorentz factor evolves as $\Gamma\propto(1+z)^{0.8-1}$. They also find that the index is larger for low-luminosity subsamples ($m\approx2$--$4$), indicating the evolution is not universal but stronger among fainter jets. On this basis they argue that at $z=3$ the co-moving density of jetted AGN is about 290 per Gpc$^3$, comparable to the total AGN density, so jets were much more common in the early universe than previously estimated.

Load-bearing premise

The argument assumes that the intrinsic characteristic energy and intrinsic damping timescale of a blazar's emission do not change with redshift; if those intrinsic quantities evolve, the measured slope would mix Doppler-factor evolution with that intrinsic evolution.

Editorial extensions

If this is right

  • At $z=3$, the co-moving space density of jetted AGN would be about 290 per Gpc$^3$, close to the total AGN density from the quasar luminosity function.
  • If local blazars have $\delta\sim10$, a $z=3$ blazar would have $\delta\sim30$--$40$, matching the hour-scale variability seen in the highest-redshift gamma-ray blazars.
  • The jet Lorentz factor would grow as $\Gamma\propto(1+z)^{0.8-1}$, implying that the most relativistic jets were more common at early cosmic times.
  • The current quasar luminosity function would underestimate the total AGN population at high redshift, with the missing sources likely in obscured phases.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The strong index found for low-luminosity subsamples ($m\approx2$--$4$) goes beyond the paper's average result; if real, it suggests the Doppler-factor evolution depends on jet power or accretion state, which the paper notes cannot be explored further with the current sample size.
  • A clean test would use Doppler factors measured by methods independent of $E_b$ and $\tau_{\rm DRW}$ (e.g., synchrotron self-absorption or gamma-ray opacity) on a redshift-matched sample; a null result would point to intrinsic evolution of $E_b^{\rm int}$ or $\tau_{\rm int}$ masquerading as beaming evolution.
  • The result implies that luminosity functions and black-hole growth models built from blazar samples should treat beaming as redshift-dependent; the paper leaves the full demographic revision implicit.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper claims to detect a cosmological evolution of the Doppler factor in Fermi-detected gamma-ray blazars, parameterized as delta ∝ (1+z)^m, by fitting the redshift dependence of two observed quantities: the characteristic energy Eb from gamma-ray spectra of 141 blazars and the optical damped-random-walk damping timescale tau_DRW from 89 blazars. The hierarchical Bayesian regressions of log10 Eb and log10 tau_DRW against log10(1+z) are interpreted through Eq. (2) as giving m = 0.81 ± 0.12 and m = 1.09 ± 0.25, respectively. The paper then argues that high-redshift blazars have larger Doppler factors, that the jet Lorentz factor evolves as Γ ∝ (1+z)^(0.8-1), and that the comoving density of jetted AGNs at z=3 could be about 290 per Gpc^3, comparable to the total AGN density.

Significance. If the detection is real, it has substantial implications for jet physics and for the census of high-redshift AGNs, and the use of two independent observables is a strength. The paper also reports its statistical errors transparently and uses a Bayesian regression approach. However, the central claim rests on an untested assumption that the intrinsic characteristic energy and damping timescale do not evolve with redshift, and the demographic extrapolation relies on additional assumptions about the origin of the Doppler-factor trend. The significance of the result is therefore conditional on breaking this degeneracy.

major comments (3)
  1. [Section II, Eq. (2)] The central inference is built on the assumption that the intrinsic characteristic energy Eb_int and the intrinsic damping timescale tau_int_DRW do not evolve with redshift. In Eq. (2), the observed slope of log10 Eb versus log10(1+z) is interpreted as m-1, but it is actually (m-1) + d log10 Eb_int / d log10(1+z), and similarly for log10 tau_DRW the slope is (1-m) + d log10 tau_int_DRW / d log10(1+z). The paper implicitly sets both intrinsic-evolution terms to zero and provides no test of this assumption. Because Eb_int plausibly depends on the external radiation field and the electron distribution, and tau_int_DRW plausibly scales with black hole mass or accretion rate, both of which evolve with cosmic time, the reported m values are degenerate with intrinsic evolution. Concretely, intrinsic-evolution slopes of p_E ≈ -0.19 and p_tau ≈ -0.09 would reconcile both measurements with a true m = 1. To support the detection claim, the authors should either constrain the intrinsic evolution with independent physical proxies (e.g., black hole mass, accretion rate, SED modeling) or compare with direct Doppler-factor estimates from VLBI kinematics or SED fitting.
  2. [Section IV] The demographic extrapolation to roughly 290 jetted AGN per Gpc^3 at z=3 depends on additional assumptions beyond the fitted m. First, the redshift trend in delta is attributed entirely to Γ ∝ (1+z)^m, whereas delta also depends on the viewing-angle distribution and on any redshift-dependent selection bias in the Fermi sample. Second, the relation Njetted = 2Γ^2 Nblazar is applied with a redshift-dependent Γ(z). The consistency with fast variability at z≥3 is suggestive but not a quantitative test of Doppler-factor evolution. The authors should demonstrate, or at least explicitly state, that the derived Γ(z) is consistent with independent constraints before using it to revise the high-redshift jetted AGN density.
  3. [Section III, Figure 1] The subsample test does not establish that the fitted m is free of selection effects. The three subsamples are small and defined post hoc in the Lγ-z plane, and the reported exponents (e.g., subsample 1: m ≈ 2.93 ± 0.97; subsample 3: m ≈ 1.80 ± 0.71) are mutually consistent within their large uncertainties and are subject to the same intrinsic-evolution degeneracy as the full-sample fit. The statement that 'the index itself appears to be evolving' is therefore not supported by the current data. At minimum, the authors should quantify the selection function of the Fermi sample or use a joint model that includes luminosity-dependent selection.
minor comments (6)
  1. [Section II, Eq. (1)] After the sentence 'Replacing δ with δ(1+z)^m', the equations still contain δ in the denominator and numerator; the authors should clarify that δ in the rewritten equations denotes the local normalization value at z=0, not the redshift-dependent quantity.
  2. [Abstract and Section IV] The text contains the typos 'straightly' where 'directly' is intended and 'Dopper' for 'Doppler'; please correct them.
  3. [Figures 2 and 4] The right panels label the intrinsic scatter as ' [dex]' without naming the variable; please add a parameter name such as σ_int on the axis.
  4. [Section III] The paper does not specify the prior distributions or convergence diagnostics used in the hierarchical Bayesian regression; please add this information in Section III or an appendix.
  5. [Section II] It is unclear whether the spectral sample of 141 blazars and the variability sample of 89 blazars overlap; please state the number of sources common to both samples and whether the two analyses are fully independent.
  6. [Section II] The description of the ZTF expansion of the optical light-curve sample omits details of the cadence, photometric bands, and quality cuts applied prior to the GP fit; please provide these details for reproducibility.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the Doppler-evolution index m is measured directly from the slopes in Eq. (2); the intrinsic-evolution degeneracy is an astrophysical assumption, not a definitional reduction.

full rationale

The paper's derivation chain is self-contained in the sense required here. The quantity m is introduced as a free parameter in the assumed redshift scaling δ ∝ δ0(1+z)^m, and it is then estimated by linear regression of log10 Eb and log10 τ_DRW against log10(1+z) (Eq. 2). The observed slopes m−1 and 1−m are the measurement, not a rename of the input. The density estimate at z=3 is a straightforward population extrapolation using the fitted Γ scaling, so it is not statistically forced by a fitted parameter masquerading as a prediction. The main caveat — that intrinsic Eb and τ_int_DRW are assumed redshift-independent — is a real astrophysical degeneracy that affects the interpretation, but it is an assumption about external physics, not a definitional identity or self-citation. The cited Zhang, Yan, and Zhang papers used for τ_DRW are the original data and method sources and are reproducible observational inputs; citing them is not load-bearing circularity. No uniqueness theorem is imported from the authors' prior work, and no ansatz is smuggled in solely by citation. The paper itself acknowledges residual selection effects. Under the hard-rule standard, no step reduces by construction to its inputs, so the circularity score is 0.

Assumptions & free parameters 3 free parameters · 3 assumptions · 0 invented entities

The central claim rests on the untested assumption that intrinsic Eb and τ are redshift-independent, plus standard beaming relations. No invented entities are introduced.

free parameters (3)
  • m (evolution index) = 0.81 ± 0.12 (spectral), 1.09+0.25−0.24 (variability)
    The slope of log Eb (or log τ) versus log(1+z), interpreted as the Doppler factor evolution index.
  • b (normalization intercept) = 2.93 ± 0.04 (spectral), 1.65 ± 0.07 (variability)
    Intercept of the linear regression, assumed to capture the constant intrinsic quantity times δ at z=0.
  • σ (intrinsic scatter) = 0.16 ± 0.01 (spectral), 0.28 ± 0.03 (variability)
    Intrinsic scatter about the regression line.
assumptions (3)
  • domain assumption Intrinsic Eb_int and τ_int are independent of redshift.
    This allows the redshift slope to be attributed to δ; it is not tested.
  • domain assumption Observed quantities scale as τ_obs = τ_int(1+z)/δ and E_obs = E_int δ/(1+z).
    Standard relativistic beaming relations used in Section II.
  • domain assumption Log-parabolic spectral fits and DRW/GP light-curve modeling are accurate.
    The values of Eb and τ depend on these functional forms, which are adopted from prior work without alternative testing.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Detection of the cosmological evolution of the Doppler factor for blazars jets." pith.science (2026). https://pith.science/paper/WWQQJT7H

@misc{pith2026250115441,
  author       = {Pith},
  title        = {Pith review of: Detection of the cosmological evolution of the Doppler factor for blazars jets},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WWQQJT7H}},
  note         = {Machine review of arXiv:2501.15441}
}
abstract

The Doppler factor ($\delta$) is a fundamental quality for the relativistic jets from active galactic nuclei (AGNs). It is also fundamental to assessing the number of the entire population of the jetted AGN at high redshift ($z$), and therefore is important for tracing the growth of supermassive black holes (SMBHs) across cosmic time. Here we present the identification of the positive cosmic evolution of the Doppler factor in {\it Fermi}-detected bright $\gamma-$ray blazars. The redshift dependence of the Doppler factor, $\delta\propto(1+z)^{0.8}$, is measured from the observed characteristic energies in the gamma-ray spectra of 141 blazars. Moreover, the analysis of the characteristic timescales derived from modeling the long-term optical light curves of 89 blazars with Gaussian process regression gives $\delta\propto(1+z)^{1.1}$, but with a larger scatter. \textbf{Note that each index is derived from the entire sample, representing an average evolution. Interestingly, the index itself also appears to evolve, with low-luminosity sources showing stronger evolution, as indicated by a larger index.} This detection straightly suggests that relativistic jets from AGNs are much more common at high redshifts than what is previously estimated.

Figures

Figures reproduced from arXiv: 2501.15441 by the authors.

Figure 1
Figure 1. FIG. 1. Blazar sample as a function of redshift and [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Fitting results for the entire spectral sample. Left: Linear regression result for the characteristic energies; Right: [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Blazar sample as a function of redshift and [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Fitting results for the entire variability sample. Left: Linear regression result for the characteristic variability timescales; [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

27 extracted references · 15 canonical work pages

  1. [1]

    2022, The As- trophysical Journal Supplement, 260, 53, doi: 10.3847/ 1538-4365/ac6751

    Abdollahi, S., Acero, F., Baldini, L., et al. 2022, The As- trophysical Journal Supplement, 260, 53, doi: 10.3847/ 1538-4365/ac6751

    work page 2022

  2. [2]

    S., Romani, R

    Ajello, M., Shaw, M. S., Romani, R. W., et al. 2012, The Astrophysical Journal, 751, 108, doi: 10.1088/ 0004-637X/751/2/108

    work page 2012

  3. [3]

    2020, The Astrophysical Journal, 892, 105, doi: 10.3847/ 1538-4357/ab791e

    Ajello, M., Angioni, R., Axelsson, M., et al. 2020, The Astrophysical Journal, 892, 105, doi: 10.3847/ 1538-4357/ab791e

    work page 2020

  4. [4]

    H., Lott, B., & The Fermi-LAT collaboration

    Ballet, J., Bruel, P., Burnett, T. H., Lott, B., & The Fermi-LAT collaboration. 2023, arXiv e-prints, arXiv:2307.12546, doi: 10.48550/arXiv.2307.12546

    arXiv doi:10.48550/arxiv.2307.12546 2023

  5. [5]

    2024, arXiv e-prints, arXiv:2407.07236, doi: 10.48550/arXiv.2407

    Banados, E., Momjian, E., Connor, T., et al. 2024, arXiv e-prints, arXiv:2407.07236, doi: 10.48550/arXiv.2407. 07236

    arXiv doi:10.48550/arxiv.2407.07236 2024

  6. [6]

    Belladitta, S., Caccianiga, A., Diana, A., et al. 2022, 5 0.0 0.1 0.2 0.3 0.4 0.5 0.6 log10(1 + z) 1.00 1.25 1.50 1.75 2.00 2.25 2.50log10( obs) m = 1.09+0.25 0.24 1.35 1.50 1.65 1.80 b b = 1.65+0.07 0.07 0.0 0.5 1.0 1.5 2.0 m 0.20 0.25 0.30 0.35 0.40 [dex] 1.35 1.50 1.65 1.80 b 0.20 0.25 0.30 0.35 0.40 [dex] [dex] = 0.28+0.03 0.03 FIG. 4. Fitting results...

    work page 2022

  7. [7]

    2019, Annu

    Blandford, R., Meier, D., & Readhead, A. 2019, Annu. Rev. Astron. Astrophys, 57, 467, doi: 10.1146/ annurev-astro-081817-051948

    work page 2019

  8. [8]

    P., Ros, E., & Zensus, J

    Boccardi, B., Krichbaum, T. P., Ros, E., & Zensus, J. A. 2017, aapr, 25, 4, doi: 10.1007/s00159-017-0105-6

    doi:10.1007/s00159-017-0105-6 2017

Show all 27 references
  1. [9]

    J., Shen, Y., Blaes, O., et al

    Burke, C. J., Shen, Y., Blaes, O., et al. 2021, Science, 373, 789, doi: 10.1126/science.abg9933

    2021 doi

  2. [10]

    2024, Astronomy and Astrophysics, 684, A98, doi: 10.1051/ 0004-6361/202348561

    Caccianiga, A., Ighina, L., Moretti, A., et al. 2024, Astronomy and Astrophysics, 684, A98, doi: 10.1051/ 0004-6361/202348561

    2024

  3. [11]

    2022, Mon

    Diana, A., Caccianiga, A., Ighina, L., et al. 2022, Mon. Not. R. Astron. Soc., 511, 5436, doi: 10.1093/mnras/ stac364

    2022 doi

  4. [12]

    2017, The Astronomical Journal, 154, 220, doi: 10.3847/1538-3881/aa9332

    Foreman-Mackey, D., Agol, E., Ambikasaran, S., & Angus, R. 2017, The Astronomical Journal, 154, 220, doi: 10.3847/1538-3881/aa9332

    2017 doi

  5. [13]

    2013, Mon

    Ghisellini, G., Haardt, F., Della Ceca, R., Volonteri, M., & Sbarrato, T. 2013, Mon. Not. R. Astron. Soc., 432, 2818, doi: 10.1093/mnras/stt637

    2013 doi

  6. [14]

    2014, Nature, 515, 376, doi: 10.1038/ nature13856

    Ghisellini, G., Tavecchio, F., Maraschi, L., Celotti, A., & Sbarrato, T. 2014, Nature, 515, 376, doi: 10.1038/ nature13856

    2014

  7. [15]

    C., Cohen, M

    Homan, D. C., Cohen, M. H., Hovatta, T., et al. 2021, The Astrophysical Journal, 923, 67, doi: 10.3847/ 1538-4357/ac27af

    2021

  8. [16]

    S., et al

    Lambrides, E., Chiaberge, M., Long, A. S., et al. 2024, The Astrophysical Journal Letter, 961, L25, doi: 10. 3847/2041-8213/ad11ee

    2024

  9. [17]

    2018, The Astrophysical Journal, 853, 159, doi: 10.3847/ 1538-4357/aaa3fb

    Li, S., Xia, Z.-Q., Liang, Y.-F., Liao, N.-H., & Fan, Y.-Z. 2018, The Astrophysical Journal, 853, 159, doi: 10.3847/ 1538-4357/aaa3fb

    2018

  10. [18]

    2019, The Astrophysical Journal Letter, 879, L9, doi: 10.3847/ 2041-8213/ab2893

    Liao, N.-H., Dou, L.-M., Jiang, N., et al. 2019, The Astrophysical Journal Letter, 879, L9, doi: 10.3847/ 2041-8213/ab2893

    2019

  11. [19]

    M., Assef, R

    Padovani, P., Alexander, D. M., Assef, R. J., et al. 2017, aapr, 25, 2, doi: 10.1007/s00159-017-0102-9

    2017 doi

  12. [20]

    2019, Mon

    Qu, Y., Zeng, H., & Yan, D. 2019, Mon. Not. R. Astron. Soc., 490, 758, doi: 10.1093/mnras/stz2651

    2019 doi

  13. [21]

    2015, Mon

    Sbarrato, T., Ghisellini, G., Tagliaferri, G., et al. 2015, Mon. Not. R. Astron. Soc., 446, 2483, doi: 10.1093/ mnras/stu2269

    2015

  14. [22]

    2022, Astronomy and Astrophysics, 663, A147, doi: 10.1051/0004-6361/202243569

    —. 2022, Astronomy and Astrophysics, 663, A147, doi: 10.1051/0004-6361/202243569

    2022 doi

  15. [23]

    F., Faucher-Gigu` ere, C.-A., et al

    Shen, X., Hopkins, P. F., Faucher-Gigu` ere, C.-A., et al. 2020, Mon. Not. R. Astron. Soc., 495, 3252, doi: 10. 1093/mnras/staa1381

    2020

  16. [24]

    2011, Mon

    Volonteri, M., Haardt, F., Ghisellini, G., & Della Ceca, R. 2011, Mon. Not. R. Astron. Soc., 416, 216, doi: 10. 1111/j.1365-2966.2011.19024.x

    2011

  17. [25]

    2017, The Astrophysical Journal, 846, 78, doi: 10.3847/ 1538-4357/aa8463

    Yuan, Z., Wang, J., Zhou, M., Qin, L., & Mao, J. 2017, The Astrophysical Journal, 846, 78, doi: 10.3847/ 1538-4357/aa8463

    2017

  18. [26]

    2022, The Astrophysical Journal, 930, 157, doi: 10.3847/1538-4357/ac679e

    Zhang, H., Yan, D., & Zhang, L. 2022, The Astrophysical Journal, 930, 157, doi: 10.3847/1538-4357/ac679e

    2022 doi

  19. [27]

    2023, The Astrophysical Journal, 944, 103, doi: 10

    —. 2023, The Astrophysical Journal, 944, 103, doi: 10. 3847/1538-4357/acafe5

    2023

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.