{"id":"9a9f1e3b-ebd1-44b7-9770-c5c13ab00d50","arxiv_id":"2501.15441","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"Blazar jets appear to have Doppler factors that increase with redshift, δ ∝ (1+z)^{0.8-1.1}, implying more numerous high-redshift jetted AGN.","lead":"This paper claims that the Doppler factor of blazar jets grows with redshift, roughly as (1+z)^0.8 to (1+z)^1.1, based on gamma-ray spectral energies of 141 blazars and optical variability timescales of 89 blazars. If true, high-redshift jets would be more common and more powerful than previously thought, changing estimates of supermassive black hole growth.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The inferred Doppler evolution index m is degenerate with the assumed constancy of the intrinsic characteristic energy Eb_int and damping timescale tau_int; the paper provides no independent calibration of that assumption, so the detection is not yet isolated.","rationale":"The paper's statistical analysis is transparent, the hierarchical Bayesian regression is appropriate, and using two independent observables (gamma-ray Eb and optical tau_DRW) is a genuine strength. The subsample analysis in the L_gamma-z plane is a reasonable attempt to address luminosity-related selection, and the final paragraph honestly concedes that selection cannot be fully removed. However, the single most load-bearing condition for the central claim is that Eb_int and tau_int_DRW do not evolve with redshift. Without an independent calibration of that assumption, the observed slopes in Eq. (2) cannot uniquely identify Doppler-factor evolution. This is not a formal inconsistency in the fitting procedure; it is a physical degeneracy that the data alone cannot break. The reader's weakest_assumption identifies exactly this issue, and the proposed test using direct Doppler measurements would settle whether the concern lands. Given that the paper is otherwise publishable with a clear caveat or reframing, the conditional verdict remains appropriate and no adjustment is needed.","tokens_in":7060,"tokens_out":8534,"duration_ms":97575,"concrete_test":"For the overlapping FSRQs, obtain independent Doppler factors from MOJAVE 15 GHz apparent jet speeds (or, where unavailable, from one-zone SED/gamma-gamma opacity modeling; cf. Homan et al. 2021). Fit log10(delta_ind) = a + m_direct log10(1+z) over the same redshift range and compare m_direct with m=0.8-1. Then compute the residual intrinsic quantities Eb_obs(1+z)/delta_ind and tau_obs delta_ind/(1+z) and test whether either retains a significant log10(1+z) trend. If m_direct is consistent with 0.8-1 and the residuals are flat, the assumption of constant intrinsic parameters is supported and the claim stands. If m_direct is consistent with zero or the residual intrinsic trends remain significant, the reported m is not isolated from intrinsic redshift evolution and the central claim would need to be reframed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (2) is the load-bearing step. The regression attributes the entire redshift slope of log10(Eb) and log10(tau_DRW) to the Doppler factor: the fitted slopes are interpreted as (m-1) and (1-m), respectively. But the actual observed slopes are (m-1)+p_E and (1-m)+p_tau, where p_E = d log10(Eb_int)/d log10(1+z) and p_tau = d log10(tau_int_DRW)/d log10(1+z). The paper implicitly sets p_E = p_tau = 0 and never tests this. For FSRQs, the intrinsic Compton peak energy depends on the external radiation field, the electron distribution, and the accretion/black-hole state, while the optical damping timescale plausibly scales with black hole mass or accretion rate; both are known to evolve with cosmic time. The two reported values, m_E = 0.81 +/- 0.12 and m_tau = 1.09 +/- 0.25, agree within about 1 sigma, but small nonzero intrinsic-evolution slopes (e.g., p_E ~ -0.19 and p_tau ~ -0.09) could reconcile both with a true m = 1. The paper also does not compare with direct Doppler-factor estimates from VLBI kinematics or SED modeling. The final paragraph acknowledges that selection effects cannot be fully eliminated, but the intrinsic-evolution degeneracy is not flagged. The headline conclusion, including the estimate of ~290 jetted AGN per Gpc^3 at z=3, therefore rests on an unverified assumption that Eb_int and tau_int are redshift-independent.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":7387,"tokens_out":4850,"duration_ms":43708,"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":[{"comment":"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.","section":"Section II, Eq. (2)"},{"comment":"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.","section":"Section IV"},{"comment":"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.","section":"Section III, Figure 1"}],"minor_comments":[{"comment":"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.","section":"Section II, Eq. (1)"},{"comment":"The text contains the typos 'straightly' where 'directly' is intended and 'Dopper' for 'Doppler'; please correct them.","section":"Abstract and Section IV"},{"comment":"The right panels label the intrinsic scatter as ' [dex]' without naming the variable; please add a parameter name such as σ_int on the axis.","section":"Figures 2 and 4"},{"comment":"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.","section":"Section III"},{"comment":"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.","section":"Section II"},{"comment":"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.","section":"Section II"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a timely question and contains useful statistical work, but the central claim of a detection is not yet isolated from the assumed redshift-independence of the intrinsic quantities. I recommend major revision rather than rejection because the issue is identifiable and potentially addressable with additional analysis or a substantially more cautious interpretation. The self-citation to the authors' previous work on tau_DRW measurements is appropriate given the continuity of the methods."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Things you should know: the paper is a clean, transparent measurement of how the gamma-ray characteristic energy and optical DRW damping timescale scale with (1+z), and it interprets the resulting slopes as evidence that the blazar Doppler factor evolves as (1+z)^{0.8-1}. What is actually new is the measurement, not the formalism. Two independent observables, a clear hierarchical Bayesian regression, and an explicit average index that two probes roughly agree on—that is worth taking seriously.\n\nThe soft spot is the step from observed slopes to Doppler factor. Equation (2) writes log10 Eb = const + (m-1) log10(1+z), and analogously for tau_DRW, with the intercept treated as a redshift-independent constant. But the intercept contains Eb_int and tau_int. If those evolve with z, the fitted slope measures (m-1)+p_E and (1-m)+p_tau, not (m-1) and (1-m). The paper never justifies p_E=p_tau=0. In FSRQs, the Compton peak energy and the accretion-related damping timescale are both expected to have some cosmic evolution. Small nonzero intrinsic slopes (p_E ~ -0.19, p_tau ~ -0.09) could reconcile both probes with a true m=1. So the agreement between the two fitted indices, within about 1 sigma, does not isolate Doppler evolution.\n\nA second concern is the luminosity-binned subsamples, which give m ~ 1.3-2.9 for the spectral sample and m ~ 1.7 for variability. The paper reads this as the index itself evolving. The natural alternative is that binning in the L_gamma-z plane selects on beaming and leaves residual selection effects. The final paragraph acknowledges that selection effects cannot be entirely eliminated, but it does not flag the more basic intrinsic-evolution degeneracy. There is also no comparison with direct Doppler-factor estimates from VLBI kinematics or SED modeling.\n\nI agree with the reader's conditional verdict. If the claim survives, it matters for the high-z jetted AGN census—the ~290 Gpc^-3 estimate at z=3 is a direct but model-dependent extrapolation. As written, the detection is underdetermined. This is an addressable flaw, not a fatal one. I would send the paper to a serious referee who can force the authors to test the intrinsic-evolution assumption, compare with independent Doppler estimates, and reframe the claim as a measurement of apparent redshift trends if necessary. The descriptive part—how Eb and tau_DRW depend on redshift—is worth preserving regardless.","headline":"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.","tokens_in":7891,"tokens_out":2253,"would_cite":false,"duration_ms":21458,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Blazar jets are more strongly beamed at high redshift, according to gamma-ray spectra and optical variability of two blazar samples.","keywords":["blazars","Doppler factor","relativistic jets","cosmological evolution","gamma-ray spectra","damped random walk","Gaussian process regression","AGN census"],"falsifier":"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.","tokens_in":6869,"feed_emoji":"🔭","tokens_out":9288,"duration_ms":75393,"temperature":0.7,"pith_summary":"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.","feed_headline":"Blazar jets grow more beamed with redshift, data show","feed_subtitle":"Gamma-ray spectra and optical variability both point to stronger beaming at high redshift, changing the jet census.","key_machinery":"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}$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"supplies the blazar sample with reliable redshifts spanning 0.065-3.033.","marker":"[14]"},{"why":"provides the 15-year Fermi-LAT spectra and the log-parabolic fits from which Eb is measured.","marker":"[3]"},{"why":"provides the optical light-curve modeling and damping timescales for the variability sample, extended here with ZTF data.","marker":"[27]"},{"why":"supplies the Gaussian-process regression method used to fit damped random walks to the optical light curves.","marker":"[12]"},{"why":"gives the z=3 blazar space density and jetted-fraction estimates used to quantify the census implication.","marker":"[11]"},{"why":"provides the bolometric quasar luminosity function used to estimate the total AGN density at z=3.","marker":"[23]"},{"why":"documents the positive correlation between gamma-ray luminosity and Doppler factor used to assess selection effects.","marker":"[15]"}],"fun_headline_variants":["Doppler factor rises with redshift in blazar jets","Blazar jet beaming intensifies across cosmic time","Cosmic evolution of blazar Doppler factor detected","High-redshift blazars show stronger jet beaming","Blazar jets increasingly beamed at high redshift"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Doppler factor rises with redshift in blazar jets","Blazar jet beaming intensifies across cosmic time","Cosmic evolution of blazar Doppler factor detected","High-redshift blazars show stronger jet beaming","Blazar jets increasingly beamed at high redshift"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000243,"raw_usage":{"total_tokens":1575,"prompt_tokens":1041,"completion_tokens":534,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":657,"completion_tokens_details":{"reasoning_tokens":459}},"tokens_in":657,"tokens_out":534,"duration_ms":5075,"temperature":1.0,"reasoning_tokens":459,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:17:06.763070+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"2014, Nature, 515, 376, doi: 10.1038/ nature13856","cited_arxiv_id":null,"evidence_quote":"supplies the blazar sample with reliable redshifts spanning 0.065-3.033."},{"cited_title":"2020, The Astrophysical Journal, 892, 105, doi: 10.3847/ 1538-4357/ab791e","cited_arxiv_id":null,"evidence_quote":"provides the 15-year Fermi-LAT spectra and the log-parabolic fits from which Eb is measured."},{"cited_title":"2023, The Astrophysical Journal, 944, 103, doi: 10","cited_arxiv_id":null,"evidence_quote":"provides the optical light-curve modeling and damping timescales for the variability sample, extended here with ZTF data."},{"cited_title":"F., Faucher-Gigu` ere, C.-A., et al","cited_arxiv_id":null,"evidence_quote":"provides the bolometric quasar luminosity function used to estimate the total AGN density at z=3."},{"cited_title":"C., Cohen, M","cited_arxiv_id":null,"evidence_quote":"documents the positive correlation between gamma-ray luminosity and Doppler factor used to assess selection effects."}],"review_version":1}