REVIEW 2 major objections 3 minor 25 references
Variable jet Lorentz factors can explain soft self-absorbed radio spectra of accreting black-holes
T0 review · 2 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Soft radio core spectra in accreting black holes can arise from jets whose Doppler factor grows outward, removing the need for energy injection at large distances.
desk verdict A plausible new mechanism for soft core radio spectra — variable Doppler factor — but the acceleration branch currently leans on an acknowledged non-self-consistent simplification, so the claim needs tempering or a self-consistent follow-up. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the partially self-absorbed synchrotron spectrum of a conical jet, computed by integrating the radiative-transfer equation along the projected jet (equations 3–4). The controlling identity is the spectral-index formula $\alpha = [5a+3b+2(b-1)p-13]/[2a-2+b(p+2)]$ from Königl (1981): with the standard Blandford–Königl scalings $N\propto\xi^{-2}$, $B\propto\xi^{-1}$, $r\propto\xi$ it gives $\alpha=0$ for constant $\Gamma$. Replacing constant $\Gamma$ by $\Gamma=\Gamma_{\rm fin}(\xi/\xi_{\rm max})^q$ lets the Doppler factor $\delta(\xi)$ vary, and because the flux in the self-absorbed regime scales roughly as $\delta^2$, a Doppler factor that rises outward enhances the large-distance contribution and softens the spectrum. The sign of the $\delta(\Gamma)$ derivative with respect to $\xi$ is set by whether the viewing angle is below or above $\arcsin(1/\Gamma_0)$.
What would settle it
A self-consistent calculation of a conical jet with variable Lorentz factor, in which the electron-density and magnetic-field profiles are transformed in the comoving frame rather than held at their constant-$\Gamma$ forms, would settle the point: if every such solution gives $\alpha\ge0$, the mechanism fails. Observationally, one could measure the Doppler-factor gradient along a jet with multi-epoch very-long-baseline interferometry in a source showing $\alpha<0$ and check that $\delta$ indeed rises with distance.
Extended reading notes
Core claim
The central claim is that the spectral index of partially synchrotron self-absorbed jet emission, $\alpha$, is controlled not only by the radial profiles of electron density and magnetic field but also by the gradient of the jet Doppler factor. For the standard conical-jet scalings, which give $\alpha=0$ at constant Lorentz factor, an outward increase of $\delta(\xi)$ boosts the emission from large distances and drives $\alpha$ below zero without any increase of the electron or magnetic energy flux. The sign of the effect depends on the viewing angle: near the characteristic angle $i=\arcsin(1/\Gamma_0)$, both acceleration and deceleration lower $\delta$, whereas smaller angles make $\delta$ rise with acceleration and larger angles make it rise with deceleration. Applied to Mrk 421 and Mrk 501, the model reproduces their soft radio-to-IR spectra with slowly varying Lorentz factors, $q=0.1$ and $q=-0.15$ respectively.
Load-bearing premise
The argument assumes that when the jet accelerates or decelerates, the density, magnetic field, and radius still follow the same simple power laws along the jet that they would in a constant-speed jet; if changing the speed reshapes those profiles, the predicted spectral index could change.
Editorial extensions
If this is right
- Observed soft core spectra with $\alpha<0$ no longer force the conclusion that relativistic electrons or magnetic flux are being injected at large distances; a Doppler-factor gradient can do the same work.
- A single source can have the same sign of spectral softening for opposite kinematic behaviours: acceleration in jets seen at small angles and deceleration in jets seen near or above $1/\Gamma_0$, so the spectral index alone does not distinguish acceleration from deceleration without knowing the viewing angle.
- The model reproduces the quiescent radio-to-X-ray continua of Mrk 421 with an accelerating jet and Mrk 501 with a decelerating jet, showing the mechanism is viable for BL Lac objects.
- Because the mechanism operates through partially self-absorbed emission, it preserves the core-shift phenomenon that argues for synchrotron self-absorption, unlike optically thin alternatives.
- The same variable-$\Gamma$ effect should apply to hard-state black-hole binaries, extending the explanation beyond AGN cores.
Reading between the lines
- A natural next step is to treat $N$, $B$, and $r$ self-consistently under a variable $\Gamma$; the paper retains the constant-$\Gamma$ power laws as an illustrative simplification, so the quantitative value of $\alpha$ in a fully relativistic jet model could differ significantly.
- If the mechanism holds, the distribution of core spectral indices ($\langle\alpha\rangle\approx 0$, $\sigma_\alpha\approx 0.4$) could be mapped onto a distribution of Doppler-factor gradients, offering an indirect kinematic census of jets that complements VLBI proper-motion measurements.
- The same Doppler-gradient logic applies to the spectral hardening ($\alpha>0$) seen in many cores: some of those sources may be jets in which $\delta$ decreases outward, rather than systems with genuinely dissipative electron or magnetic energy profiles.
- A focused test would compare fitted $\delta(\xi)$ profiles with direct very-long-baseline interferometry measurements of jet acceleration in a small sample of nearby BL Lacs; the mechanism predicts agreement in sign between the spectral softening and the measured Doppler-factor gradient.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes that soft partially synchrotron self-absorbed radio spectra (α < 0) from the cores of radio-loud AGN can be produced by jets whose Doppler factor increases with distance, without the large energy deposition at large radii that constant-Γ models would require. It adopts the Blandford–König framework (Eqs. 3–4), retains the standard power-law scalings for the electron density, magnetic field, and conical radius (Eq. 7), and allows a power-law varying bulk Lorentz factor Γ ∝ ξ^q (Eq. 8). The paper argues that δ(ξ) increases with distance either for accelerating jets seen at small angles or for decelerating jets seen at larger angles, and it presents fits to the quiescent radio-to-X-ray spectra of Mrk 421 (q = 0.1, acceleration) and Mrk 501 (q = −0.15, deceleration). The author explicitly states that the fits are illustrative, non-unique, and based on simplified scalings that neglect the effect of a variable Γ on the comoving-frame quantities.
Significance. If the mechanism is correct, it offers an observationally motivated alternative to the energetic requirement for soft core spectra and directly connects the observed α < 0 population to measured jet acceleration/deceleration. The paper is transparent in its assumptions, uses a full radiative-transfer integral rather than only the δ^2 approximation, and its two example fits demonstrate numerical viability under the stated assumptions. The main weakness is that the adopted scalings are not self-consistent with a varying Γ; in particular, the accelerating branch is subject to a competing hardening effect, so the general 'either acceleration or deceleration' claim is not yet established. This is a fixable but load-bearing gap, and the deceleration branch is substantially more robust than the acceleration branch.
major comments (2)
- [Section 2, Eq. (7); Section 3, Fig. 4] The paper allows Γ ∝ ξ^q (Eq. 8) while keeping the constant-velocity scalings N ∝ ξ^−2 and B ∝ ξ^−1 (Eq. 7). For a steady conical jet with conserved particle number and magnetic flux, a varying Γ changes these scalings; substituting the self-consistent scalings N ∝ (Γξ^2)^−1 and B ∝ (Γξ)^−1 into the paper's own spectral-index formula (Eq. 1) gives α = 2q/(1 + q), which is positive for an accelerating jet (q > 0). The Mrk 421 fit (Fig. 4) uses q = 0.1 and does not include this correction, so it does not demonstrate that the δ(ξ) gradient can overcome the baseline hardening. The caveat in Section 2 about Eq. (7) is candid, but the abstract and conclusions still claim the acceleration branch as part of the explanation. The author should either provide a self-consistent calculation (analytical or numerical) showing that the Doppler effect dominates, or restrict the main claim to the deceleration branch.
- [Section 2, Eqs. (3)–(4); Section 4, Conclusions] The paper's general statement that soft spectra 'can be obtained if the Doppler factor ... increases with distance' is inferred from the δ^2 scaling of the optically thin case, but in the partially self-absorbed regime the δ dependence enters inside the ξ-integral of Eq. (3) and is convolved with the τ_sa structure. No asymptotic or parameter-space analysis is given for the resulting spectral index as a function of q and of the viewing-angle offset from 1/Γ. The two fits are point examples rather than a demonstration of the claimed dichotomy between acceleration and deceleration. A short analytic treatment or a parameter scan for α(q, i) would considerably strengthen the central claim.
minor comments (3)
- [Section 3, Figs. 4–5] The fits have many free parameters (Γ_fin, q, i, Θ, p, N0, B0, z0, zmax, η_acc, M) and no quoted uncertainties or degeneracy analysis; the author acknowledges this, but the abstract wording 'we find our model can explain' is stronger than the illustrative fits support. Consider adding a sentence in the abstract or conclusions that the fits are existence proofs rather than unique models.
- [Section 3, Fig. 5] The caption states that the feature around 1 eV is the host galaxy and is not subtracted from the average spectrum; please clarify whether this feature is included in the fit or is excluded from the model comparison, since an unmodeled component in the data can affect how the fit is assessed.
- [Section 2, Fig. 1] The validation of Eqs. (3)–(4) against the exact angular dependence in Fig. 1 is given for Γ = 20 and p = 2.5 only; a sentence stating the expected range of validity in Γ and p would help readers judge the approximation for the larger and smaller Γ values used in the fits.
Circularity Check
No significant circularity: the spectral softening is a computed output of the radiative-transfer model, not an assumed or fitted result.
full rationale
The central claim is not circular. The paper computes the partially self-absorbed spectrum from the radiative-transfer integral Eqs. (3)-(4), with the Doppler factor δ(ξ) entering through Eq. (2); the constant-Γ case (q=0) yields α=0 and fails to reproduce the radio-to-IR data (Fig. 3), while the variable-Γ cases with q=0.1 and q=-0.15 are integrated explicitly and give α≈-0.2 (Figs. 4, 5). Thus the softening is a computed output, not an input or a renamed parameter. The fits to Mrk 421 and Mrk 501 use free parameters, but the paper explicitly disclaims uniqueness and says the fits 'are intended to show the viability of the general model,' so the conclusion is an existence/demonstration claim, not a fitted quantity presented as a prediction. The self-citations (Zdziarski, Lubiński & Sikora 2012; Zdziarski, Stawarz & Sikora 2019) supply the standard synchrotron source-function and absorption-coefficient formalism; that formalism is not the target result and is independent of the variable-Γ mechanism, so it is not load-bearing circular self-citation. No uniqueness theorem from the authors is invoked to forbid alternatives; in fact alternatives (Potter & Cotter 2013a,b; optically thin emission) are discussed. The one substantive caveat is acknowledged by the paper itself in Section 2, Eq. (7): the BK scalings N∝ξ^-2, B∝ξ^-1 are retained 'while these dependencies neglect the effect of variable Γ on the quantities in the comoving frame,' with a self-consistent treatment deferred to work in preparation. This is a physical-consistency limitation that may affect the acceleration branch, but it is not circular: the predicted softening follows from the stated equations rather than being identical to the assumed δ(ξ) profile.
Assumptions & free parameters
free parameters (11)
- Γ_fin (terminal Lorentz factor) =
30 (Mrk 421), 10 (Mrk 501)
- q (Lorentz factor power-law index) =
0.1 (Mrk 421), -0.15 (Mrk 501)
- i (viewing angle) =
0.5/Γ_fin (Mrk 421), 3/Γ_fin (Mrk 501)
- Θ (jet half-opening angle) =
0.2/Γ_fin (Mrk 421), 0.1/Γ_fin (Mrk 501)
- p (electron power-law index) =
2.4 (Mrk 421), 2.54 (Mrk 501)
- N0 (electron density normalization) =
1.7e10 cm^-3 (Mrk 421), 3.5e8 cm^-3 (Mrk 501)
- B0 (magnetic field at z0) =
40 G (Mrk 421), 400 G (Mrk 501)
- z0 (onset distance) =
20 rg (Mrk 421), 60 rg (Mrk 501)
- zmax (terminal distance) =
100 pc (both)
- η_acc (acceleration efficiency) =
7e-7 (Mrk 421), 5e-5 (Mrk 501)
- M (black hole mass) =
2e8 M_sun (Mrk 421), 1e9 M_sun (Mrk 501)
assumptions (5)
- standard math The Blandford-Konigl model of partially self-absorbed synchrotron emission from a conical jet is valid.
- domain assumption The electron distribution in the jet frame is a power law in γ with index p.
- domain assumption The magnetic field, electron density, and jet radius follow the standard BK79 dependencies (N∝ξ^-2, B∝ξ^-1, r∝ξ) even when Γ varies.
- domain assumption The observed radio core spectra are dominated by partially self-absorbed synchrotron emission, not optically thin emission.
- domain assumption The quiescent spectra of Mrk 421 and Mrk 501 are representative and the data are reliable.
Cite this review
Pith. "Pith review of Variable jet Lorentz factors can explain soft self-absorbed radio spectra of accreting black-holes." pith.science (2026). https://pith.science/paper/WNPBYK6J
@misc{pith2026190804216,
author = {Pith},
title = {Pith review of: Variable jet Lorentz factors can explain soft self-absorbed radio spectra of accreting black-holes},
year = {2026},
howpublished = {\url{https://pith.science/paper/WNPBYK6J}},
note = {Machine review of arXiv:1908.04216}
}
abstract
We study the effect of variable jet bulk Lorentz factors, i.e., either jet acceleration or deceleration, on partially synchrotron self-absorbed radio spectra from cores of radio-loud active galactic nuclei and black-hole binaries in the hard state. In about a half of quasars and radio galaxies, their core radio spectra are observed to be soft, i.e., have the spectral index of $\alpha<0$. If they are emitted by jets with constant Lorentz factors, that softness implies deposition of large amounts of energy at large distances from the centre. We show here that such soft spectra can be explained without that energetic requirement by emission of jets with the Doppler factor increasing with the distance. This can happen for either jet acceleration or deceleration, depending on the jet viewing angle. We find our model can explain the quiescent radio to X-ray spectra of the BL Lac objects Mrk 421 and Mrk 501.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
- [1]
- [2]
- [3]
- [4]
-
[5]
D., K \"o nigl A., 1979, ApJ, 232, 34 (BK79)
Blandford R. D., K \"o nigl A., 1979, ApJ, 232, 34 (BK79)
work page 1979
-
[6]
Cawthorne T. V., 1991, in P. A. Hughes, ed., Beams and Jets in Astrophysics. Cambridge Univ.\ Press, Cambridge, p.\ 187
work page 1991
-
[7]
Clausen-Brown E., Savolainen T., Pushkarev A. B., Kovalev Y. Y., Zensus J. A., 2013, A&A, 558, A144
work page 2013
-
[8]
Fender, R. P., Pooley, G. G., Durouchoux, P., Tilanus, R. P. J., Brocksopp, C., 2000, MNRAS, 312, 853
work page 2000
Show all 25 references
-
[9]
W., Fabian A
Guilbert P. W., Fabian A. C., Rees M. J., 1983, MNRAS, 205, 593
1983
-
[10]
S., 2011, Mem.\ Soc.\ Astron.\ Ital., 82, 95
Komissarov S. S., 2011, Mem.\ Soc.\ Astron.\ Ital., 82, 95
2011
-
[11]
K \"o nigl A., 1981, ApJ, 243, 700
1981
-
[12]
A., Bridle A
Laing R. A., Bridle A. H., 2004, MNRAS, 348, 1459
2004
-
[13]
R., Blandford R
Lind K. R., Blandford R. D., 1985, ApJ, 295, 358
1985
-
[14]
L., et al., 2019, ApJ, 874, 43
Lister M. L., et al., 2019, ApJ, 874, 43
2019
-
[15]
J., Cotter G., 2013a, MNRAS, 429, 1189
Potter W. J., Cotter G., 2013a, MNRAS, 429, 1189
-
[16]
J., Cotter G., 2013b, MNRAS, 436, 304
Potter W. J., Cotter G., 2013b, MNRAS, 436, 304
-
[17]
B., Kovalev Y
Pushkarev A. B., Kovalev Y. Y., Lister M. L., Savolainen T., 2009, A&A, 507, L33
2009
-
[18]
B., Hovatta T., Kovalev Y
Pushkarev A. B., Hovatta T., Kovalev Y. Y., Lister M. L., Lobanov A. P., Savolainen T., Zensus J. A., 2012, A&A, 545, A113
2012
-
[19]
M., Mannheim K., 2003, A&A, 397, 121
Rieger F. M., Mannheim K., 2003, A&A, 397, 121
2003
-
[20]
Sikora M., Madejski G., Moderski R., Poutanen J., 1997, ApJ, 484, 108
1997
-
[21]
M., Zhang B.-B., Mao J., 2018, ApJS, 239, 33
Yuan Z., Wang J., Worrall D. M., Zhang B.-B., Mao J., 2018, ApJS, 239, 33
2018
-
[22]
A., Lubi \'n ski P., Sikora M., 2012, MNRAS, 423, 663
Zdziarski A. A., Lubi \'n ski P., Sikora M., 2012, MNRAS, 423, 663
2012
-
[23]
A., Paul D., Osborne R., Rao A
Zdziarski A. A., Paul D., Osborne R., Rao A. R., 2016, MNRAS, 463, 1153
2016
-
[24]
A., Stawarz ., Sikora M., 2019, MNRAS, 485, 1210
Zdziarski A. A., Stawarz ., Sikora M., 2019, MNRAS, 485, 1210
2019
-
[25]
\.Z ywucka N., Goyal A., Jamrozy M., Ostrowski M., Stawarz ., 2014, A&A, 563, 135
2014
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.