REVIEW 3 major objections 3 minor 30 references
Relativistic polytrope from the collimation and acceleration profiles of the M87 jet at subparsec scales and thermodynamic evidence for the Blandford-Znajek mechanism
T0 review · 3 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The M87 jet's measured acceleration curve pins its plasma to the hot 4/3 polytropic index, no matter what balances the jet.
desk verdict A clean analytic result — Gamma=4/3 from beta=1 across three force-balance regimes — but the inference hangs on identifying the observed VLBI gamma with the bulk-flow Lorentz factor, and the Blandford-Znajek and pair-plasma conclusions outrun the evidence. 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 set of conserved quantities along magnetic flux tubes in stationary axisymmetric ideal MHD, together with the radial force-balance equation. The conserved integrals are magnetic flux, the Ferraro isorotation frequency $\Omega_F$, the mass-flux integral $\eta = \gamma \rho v_p / B_p$, energy $\mathcal{E} = \gamma h \eta - \Omega_F I / 2\pi$, angular momentum $L = \gamma h \eta r v_\varphi - I/2\pi$, and entropy $S = p/\rho^\Gamma$. The author feeds the observed power laws $\gamma \propto r^\beta$ with $\beta = 1$ and $r \propto z^\alpha$ into three limiting equilibria. In the transverse electromagnetic case the key identity is the near cancellation $E_r \approx B_\varphi$ from ideal conductivity, which leaves a subdominant pressure $p \sim E_r^2 / 8\pi \gamma^2 \propto r^{-4}$; in the hot centrifugal case entropy conservation plus $\gamma r v_\varphi = \text{const}$ forces the same $r^{-4}$. With $\rho \propto r^{-2-\beta}$, entropy conservation then requires $\Gamma = 4/3$ independently of the chosen equilibrium.
What would settle it
A precise re-measurement of the exponent $\beta$ in $\gamma \propto r^\beta$ on subparsec scales would settle the claim: a value clearly different from 1 changes the predicted index (e.g., $\Gamma_{\rm LM} = 4/(2+\beta)$), while near 8 marcsec the electron-proton cooling scenario predicts a distinct $\gamma \propto r^{0.46}$ segment whose presence or absence would directly test the pair-plasma interpretation.
Extended reading notes
Core claim
The paper's central claim is that the measured relation $\gamma \propto r$ (with $r \propto z^{0.58}$) at subparsec scales fixes the equation of state of the M87 jet plasma. Starting from stationary, axisymmetric ideal MHD, the author treats the jet as a stack of magnetic flux tubes on which magnetic flux, Ferraro isorotation, mass, energy, angular momentum, and entropy $S = p/\rho^{\Gamma}$ are conserved. Combining the transverse force-balance equation with the density scaling $\rho \propto \gamma^{-1} r^{-2}$ gives three pressure scalings – $p \propto r^{-4}$ for longitudinal magnetic pressure, for transverse electromagnetic pressure after the near cancellation $E_r \approx B_\varphi$, and for hot centrifugal pressure – and every route yields $\Gamma = 4/(2+\beta) = 4/3$ at the observed $\beta = 1$. The same machinery rules out a cold centrifugal flow, whose index would be sub-4/3. Consequently the flow is hot, and energy conservation forces the longitudinal electric current to be an exact integral of motion, $I = \text{const}$. The author then extends the argument beyond 8 marcsec, where the acceleration exponent drops to $\approx 0.28$: the inferred indices are consistent with a cold $\Gamma = 5/3$ flow, and the sharp jump from $4/3$ to $5/3$ without lingering at $1.44$ is taken as thermodynamic evidence for an electron-positron plasma and for partial Blandford-Znajek launching.
Load-bearing premise
The load-bearing premise is that the observed Lorentz factor is the true bulk-flow speed of a steady, axisymmetric, infinitely conducting plasma with entropy conserved along each magnetic flux tube; if the radio images track pattern speeds or if the jet dissipates energy significantly, the inferred $4/3$ index does not follow.
Editorial extensions
If this is right
- At subparsec scales the M87 jet is hot in the thermodynamic sense: random particle motion is relativistic, not just the bulk flow.
- The longitudinal electric current is exactly conserved along each magnetic tube, implying a stable internal electromagnetic-current structure that does not dissipate or redistribute over the observed region.
- The acceleration break at 8 marcsec can be interpreted as plasma cooling from $\Gamma = 4/3$ to $\Gamma = 5/3$, with no need for a change in jet collimation or force balance.
- A direct jump from $4/3$ to $5/3$ during cooling favours an electron-positron plasma, supporting models in which the Blandford-Znajek mechanism contributes to launching M87's jet.
- The same profile-inversion method could be applied to other well-resolved relativistic jets to infer their polytropic index, composition, and launch mechanism from kinematic measurements alone.
Reading between the lines
- Beyond the paper, the invariance of $\Gamma = 4/3$ across all three force-balance regimes suggests the conclusion is essentially kinematic: any stationary axisymmetric ideal-MHD jet with $\gamma \propto r$ and parabolic collimation would have to be hot, so the result may survive even if M87's internal magnetic structure is more complex than the limiting equilibria considered.
- A testable extension: near 8 marcsec, an electron-proton plasma cooling through the $\Gamma \approx 1.44$ stage should produce a distinct acceleration segment $\gamma \propto r^{0.46} \propto z^{0.27}$; searching for it with high longitudinal resolution would directly discriminate pair plasma from normal plasma.
- Beyond 8 marcsec the method gives force-balance-dependent indices ($\Gamma_{\rm LM} \approx 1.76$, $\Gamma_{\rm TEM} \approx 1.12$, $\Gamma_{\rm CC} \approx 1.64$), so a natural next step is to use independent measurements — polarization maps or internal velocity structure — to identify which transverse force actually dominates there.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This Letter derives a polytropic equation of state for the M87 jet at subparsec scales from the observed collimation profile r ∝ z^0.58 and acceleration profile γ ∝ r (Eq. 2). Under the assumptions of stationary, axisymmetric, ideal MHD flow with conserved entropy along each magnetic flux tube, the author shows that the observed index β = 1 in γ ∝ r^β implies a polytropic index Γ = 4/3, independent of whether the transverse pressure is balanced by longitudinal magnetic pressure (Eq. 24), transverse electromagnetic pressure (Eq. 35), or centrifugal pressure in a hot flow (Eq. 52). The paper further argues that Γ = 4/3 and β = 1 imply exact conservation of the longitudinal electric current (Eq. 54) and that the change of acceleration profile at 8 marcsec, with β ≈ 0.28, corresponds to plasma cooling from Γ = 4/3 to Γ = 5/3 without an intermediate value 1.44, which the author interprets as evidence for an electron-positron plasma and at least partial operation of the Blandford-Znajek mechanism.
Significance. If the central inference holds, the paper makes a substantive contribution: it converts two observational power laws into a thermodynamic constraint on the jet plasma, a conserved current, and a falsifiable prediction about the absence of a 1.44 polytropic plateau. The algebra in Section 3 is clean and internally consistent: the three force-balance regimes all yield Γ = 4/3 for β = 1, and the later derivation of current conservation follows analytically once the density and enthalpy scalings are accepted. The paper is appropriately cautious in stating that the result does not prove a pure Blandford-Znajek jet. The main risk is not the mathematics but the physical identification of the observed VLBI Lorentz factor with the bulk-flow Lorentz factor on a fixed magnetic tube, together with the assumption of strictly adiabatic, shock-free flow. These assumptions are load-bearing and need explicit defense or quantitative qualification before the central claim can be accepted.
major comments (3)
- [§1, §3.2 (Eqs. 2, 21)] The central inference treats the measured VLBI Lorentz factor as the bulk-flow Lorentz factor on a fixed magnetic flux tube, but this identification is not justified. Mertens et al. (2016) measure apparent motions of radio features; these can be pattern speeds (waves or recollimation shocks) rather than the local thermodynamic bulk speed, particularly given the expansion/recollimation cycles cited in the paper's own §1. Equation (21), ρ ∝ γ^{-1} r^{-2}, follows from η = γρv_p/B_p only if γ is the bulk Lorentz factor of the same plasma crossing the flux tube at radius r. If γ_obs is a pattern speed or samples a different streamline, then ρ(r), p(r), and the entropy matching in Eqs. (24)-(26), (35)-(36), and (52) do not follow, and β = 1 in the observed frame does not force Γ = 4/3. The manuscript should either argue that the radio features trace the bulk flow in this region, or quantify the systematic error introduced by this assumption, for example by comparing with independent bulk-speed estimates or by modeling the pattern-to-bulk ratio.
- [§2, Eqs. (14)-(15)] Entropy conservation S = p/ρ^Γ along each magnetic tube is assumed throughout, but the same region is described in §1 as containing 'a sequence of subsequent expansions and contractions' and 'three consecutive expansion/recollimation cycles' (Walker et al. 2018). Recollimation shocks and dissipation violate entropy conservation; under such conditions, the polytropic relation used to infer Γ from the scaling of p and ρ is not guaranteed. The derivation needs an explicit argument that dissipation is negligible at subparsec scales in M87, or the conclusion should be restricted to shock-free segments of the jet.
- [§3.2, Eq. (21)] The approximation v_p ≈ 1 used to obtain ρ ∝ γ^{-1} r^{-2} may not hold over the full range of z < 8 marcsec, where the measured Lorentz factors are moderate (apparent speeds of the order of a few c). For γ ≈ 2, v_p ≈ 0.87, and the density scaling acquires a γ-dependent correction that is not a pure power law; this correction propagates into the entropy-matching step. The authors should show either that this correction is negligible for the inferred Γ, or quantify the range of β and γ over which the conclusion Γ ≈ 4/3 is robust.
minor comments (3)
- [Eq. (17)] The fraction is typeset ambiguously; it should be written as (p + (B_z^2 + B_φ^2 - E_r^2)/(8π))′ to make the grouping clear.
- [§3.3] The statement that theoretical estimates give B_φ ≫ B_z and v_φ ≪ 1 relies on earlier work (Sob'yanin 2017); an independent observational or numerical check would strengthen the case.
- [Throughout] The unit 'marcsec' is used; consider using 'mas' for consistency with the standard astronomical literature.
Circularity Check
No circularity: Gamma=4/3 is inferred from the observed beta=1 scaling via conservation laws, not fitted or assumed.
full rationale
The derivation takes the observed profile gamma proportional to r (beta=1, Eq. 2) as an input and combines it with mass conservation and entropy conservation to infer the polytropic index. Specifically, Eq. (21) gives rho proportional to gamma^{-1} r^{-2} from the conserved integral eta, then with p proportional to r^{-4} in the longitudinal-magnetic branch, entropy conservation forces Gamma = 4/(2+beta), which becomes 4/3 at beta=1 (Eqs. 24-26). The transverse-electromagnetic and hot-centrifugal branches use different pressure scalings but give the same Gamma=4/3 at beta=1 (Eqs. 35-36 and 51-52). This is a deductive inference from observations, not a relabeling of a fitted parameter. The later conclusion I=const (Eq. 54) follows analytically once gamma h=const and Gamma=4/3 are established, and is not used earlier as an input. Citations to Sob'yanin 2017 and 2018 supply the radial momentum form and some interpretive estimates (B_phi >> B_z, v_phi << 1, jet-in-jet, precession), but the same polytropic result is obtained independently through other force-balance branches, so those self-citations are not load-bearing for the central claim. The main vulnerability, namely whether the observed VLBI Lorentz factor is the bulk-flow Lorentz factor on a fixed magnetic tube, is a physical identification assumption and a correctness risk, not a circularity.
Assumptions & free parameters
free parameters (2)
- Lorentz-factor-radius index beta (z < 8 marcsec) =
1
- Lorentz-factor-radius index beta (z > 8 marcsec) =
approx 0.28
assumptions (5)
- domain assumption Ideal MHD with infinite conductivity (frozen-in condition, Eq. 4).
- domain assumption Stationary, axisymmetric jet with conserved integrals of motion along magnetic tubes (particle flux eta, energy E, angular momentum L, entropy S).
- domain assumption The plasma obeys a polytropic equation of state p=(Gamma-1) rho epsilon (Eq. 13) with a single index Gamma.
- domain assumption Poloidal velocity is highly relativistic, vp approximately 1, so eta conservation gives rho proportional to gamma^-1 r^-2 (Eq. 21).
- domain assumption Entropy is conserved along each magnetic tube (Eq. 15), so no dissipation or shocks are present.
Cite this review
Pith. "Pith review of Relativistic polytrope from the collimation and acceleration profiles of the M87 jet at subparsec scales and thermodynamic evidence for the Blandford-Znajek mechanism." pith.science (2026). https://pith.science/paper/KNYYDP7T
@misc{pith2026190805485,
author = {Pith},
title = {Pith review of: Relativistic polytrope from the collimation and acceleration profiles of the M87 jet at subparsec scales and thermodynamic evidence for the Blandford-Znajek mechanism},
year = {2026},
howpublished = {\url{https://pith.science/paper/KNYYDP7T}},
note = {Machine review of arXiv:1908.05485}
}
read the original abstract
Recent Very Long Baseline Interferometry observations of the relativistic jet in the M87 radio galaxy at 43 GHz show gradual relativistic acceleration of the plasma and suggest a linear dependence of Lorentz factor on jet radius at scales up to 8 marcsec (0.65 pc) from the core (2.5 marcsec in projection). General analysis of integrals of motion being unaltered along the jet and reflecting fundamental conservation laws shows that the above dependence implies a polytropic equation of state with index 4/3. The inferred value of the polytropic index appears independent of the exact nature of forces sustaining the transverse balance of the jet and indicates exact conservation of the longitudinal electric current and hence the existence of a stable internal electromagnetic structure at the scales under consideration. At this index the flow is hot and corresponds to relativistic thermodynamic motion of particles. Considerable weakening of the acceleration efficiency after 8 marcsec with the jet form being unchanged can be related to the plasma cooling, when the polytropic index becomes 5/3. Such a sharp change in the index without intermediate delay at 1.44 during cooling favours the existence of an electron-positron plasma and requires at least partial participation of the Blandford-Znajek mechanism in the launching of the M87 jet.
Reference graph
Works this paper leans on
-
[1]
Asada K., Nakamura M., 2012, ApJ, 745, L28
work page 2012
-
[2]
Asada K., Nakamura M., Doi A., Nagai H., Inoue M., 2014, ApJ, 781, L2
work page 2014
- [3]
- [4]
- [5]
- [6]
-
[7]
Chattopadhyay I., Ryu D., 2009, ApJ, 694, 492
work page 2009
- [8]
Show all 30 references
-
[9]
J., 2010, ApJS, 191, 113
Choi E., Wiita P. J., 2010, ApJS, 191, 113
2010
-
[10]
D., 1918, Publ
Curtis H. D., 1918, Publ. Lick Obser., 13, 9
1918
-
[11]
Falle S. A. E. G., Komissarov S. S., 1996, MNRAS, 278, 586s
1996
-
[12]
Fendt C., 1997, A&A, 323, 999
1997
-
[13]
Hada K., 2017, Galaxies, 5, 2
2017
-
[14]
Hada K., Doi A., Kino M., Nagai H., Hagiwara Y., Kawaguchi N., 2011, Nature, 477, 185
2011
-
[15]
Hada K., et al., 2013, ApJ, 775, 70
2013
-
[16]
Hada K., et al., 2017, PASJ, 69, 71
2017
-
[17]
P., Bremer M., Zensus J
Kim J.-Y., Lu R.-S., Krichbaum T. P., Bremer M., Zensus J. A., Walker R. C., 2016, Galaxies, 4, 39
2016
-
[18]
Y., Lister M
Kovalev Y. Y., Lister M. L., Homan D. C., Kellermann K. I., 2007, ApJ, 668, L27
2007
-
[19]
B., Van Moer M., 2019, preprint (arXiv:1904.08428)
Liska M., Hesp C., Tchekhovskoy A., Ingram A., van der Klis M., Markoff S. B., Van Moer M., 2019, preprint (arXiv:1904.08428)
2019 arXiv
-
[20]
E., 2009, ApJ, 698, 1570
Lyubarsky Y. E., 2009, ApJ, 698, 1570
2009
-
[21]
G., 1971, ApJ, 165, 147
Mathews W. G., 1971, ApJ, 165, 147
1971
-
[22]
Meliani Z., Sauty C., Tsinganos K., Vlahakis N., 2004, A&A, 425, 773
2004
-
[23]
P., Walker R
Mertens F., Lobanov A. P., Walker R. C., Hardee P. E., 2016, A&A, 595, A54
2016
-
[24]
Nakamura M., et al., 2018, ApJ, 868, 146
2018
-
[25]
Ogihara T., Takahashi K., Toma K., 2019, ApJ, 877, 19
2019
-
[26]
N., 2017, MNRAS, 471, 4121
Sob'yanin D. N., 2017, MNRAS, 471, 4121
2017
-
[27]
N., 2018, MNRAS, 479, L65
Sob'yanin D. N., 2018, MNRAS, 479, L65
2018
-
[28]
L., 1957, The Relativistic Gas
Synge J. L., 1957, The Relativistic Gas. North-Holland, Amsterdam
1957
-
[29]
C., Hardee P
Walker R. C., Hardee P. E., Davies F. B., Ly C., Junor W., 2018, ApJ, 855, 128
2018
-
[30]
EHT Collaboration 2019, ApJ, 875, L1
2019
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.