REVIEW 1 major objections 5 minor 98 references
Hamiltonian variational reconstruction of the 3D magnetic geometry of relativistic jets: accuracy across GRMHD models and the fundamental sign degeneracy
T0 review · 1 major / 5 minor · reviewed 2026-07-13 · grok-4.5
Pith's one-line read Jet width and linear polarization recover 3D magnetic geometry but never the field's sense of direction.
desk verdict Clean free-field reconstructor plus a formal non-identifiability proof: unsigned jet geometry is recoverable from W and p at high fidelity; poloidal sense is not, and three recovery routes fail for named physical reasons. 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
H-MOG, a Hamiltonian-regularized variational inversion that reconstructs amplitude and unit direction of the magnetic field on a Cartesian lattice from projected W(z) and p(z), optimized by automatic differentiation; its even prior and data terms make the sign degeneracy exact (Proposition 1).
What would settle it
Apply H-MOG to a GRMHD library that retains full non-axisymmetric structure and higher near-axis resolution, then check whether the unsigned alignment stays near 0.95–0.98 and whether any parity-even addition (or a parity-odd observable) can raise the signed alignment above chance.
Extended reading notes
Core claim
From jet width and linear polarization alone, H-MOG recovers the unsigned three-dimensional magnetic geometry of relativistic jets at ⟨|cos|⟩ ≃ 0.95–0.98 across ten GRMHD models; the poloidal sense cannot be recovered because those observables are invariant under B → −B, and three physically motivated attempts to break the degeneracy all fail for identifiable reasons.
Load-bearing premise
The truth fields used for scoring are axisymmetric mean fields after azimuthal averaging and resampling onto a uniform coarse Cartesian cube, so the high reported alignment measures recovery of large-scale mean geometry rather than the full native turbulent field.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces H-MOG, a variational method that reconstructs the 3D magnetic field of a relativistic jet on a Cartesian lattice from two projected observables, the jet width W(z) and linear polarization degree p(z), regularized by a Hamiltonian prior and optimized with automatic differentiation. Applied to ten GRMHD models spanning MAD and SANE states and five black-hole spins, it recovers the unsigned field orientation at ⟨|cos|⟩ ≃ 0.95–0.98. Proposition 1 proves that W and p (and the retained even prior) are invariant under B → −B, so the poloidal sense is formally non-identifiable from these data. Three routes to break the degeneracy—an RM term, full 3D sampling, and a Blandford–Znajek spin prior—are tested and each fails for a distinct physical reason; the spin–sense relation in the simulated jets is shown to be non-monotonic and state-dependent. The authors conclude that recovering sense requires a parity-odd observable such as Faraday tomography or circular polarization.
Significance. If the results hold, the paper cleanly separates what current VLBI/EHT linear-polarization and width data can determine (large-scale helical geometry) from what they cannot (poloidal sense). The formal invariance argument, the multi-model GRMHD campaign with reported snapshot scatter, the synthetic-jet validation across pitch and noise, and the explicit diagnosis of three failed recovery routes are all strengths. The finding that the spin–sense relation is not monotonic in these turbulent jets is an astrophysical result in its own right and cautions against silently imposing a BZ prior. The work is directly relevant to ongoing and next-generation EHT analyses of M87* and Sgr A*.
major comments (1)
- The high ⟨|cos|⟩ scores are measured against the same axisymmetrized, azimuthally averaged Cartesian cubes that the simplified forward model (ε ∝ |B|²) sees (Sect. 3.2). The paper states this limitation and does not claim recovery of native non-axisymmetric turbulence, so the central geometry claim remains sound for the large-scale mean field. A short, explicit statement in the abstract or conclusions that the reported accuracy is for the axisymmetric mean geometry would prevent over-reading by non-specialists.
minor comments (5)
- Table 1 and Table D.1 both report reconstruction statistics; a single consolidated table (or a clear pointer that D.1 is the complete version) would reduce redundancy.
- In Sect. 2.3 the plane-of-sky angle is written χ = arctan(Bz/Bx) + π/2; a brief note on the branch choice and how it is handled under automatic differentiation would aid reproducibility.
- Fig. 12 caption and text use “verso” for the mean poloidal sense; standard English “sense” or “sign” would be clearer for an international readership.
- Appendix C and Figs. D.1–D.2 usefully document convergence; a one-sentence statement of the final data residual level (already sub-percent) in the main text of Sect. 2.4 would help readers who skip the appendix.
- The frame-dragging term H_frame is correctly set to h = 0 for the main runs; stating the numerical value of h used in the failed Test 3 (Sect. 6.4) would make that experiment fully reproducible from the text alone.
Circularity Check
No significant circularity: geometry skill is scored against independent GRMHD truth never seen by the optimizer; the sign non-identifiability is a proved invariance of the observables, not a fitted tautology.
full rationale
The paper’s load-bearing chain is self-contained and externally scored. H-MOG minimizes a cost built only from projected W(z) and p(z) plus an even Hamiltonian prior; the reconstructed field is then compared cell-by-cell to GRMHD truth cubes that the optimizer never receives. That is ordinary inverse-problem validation, not a prediction forced by a fitted input. Proposition 1 is a short, definitional invariance proof (ε∝|B|² and 2χ → 2χ+2π leave I, W, Q, U, p unchanged; retained prior terms are even in S), establishing non-identifiability of sense rather than claiming to derive sense from even data. Prior-weight sweeps (Sect. 4.6) show ⟨|cos|⟩ stays high and ⟨cos⟩ stays near zero, so the prior does not smuggle the reported geometry or flip the sense. The three degeneracy-breaking tests (RM term, full-3D sampling, BZ spin prior) are empirical failures with named physical reasons, not self-citations or uniqueness theorems imported from the same author. There is no self-citation load-bearing chain (single author; library and theory citations are external), no fitted constant renamed as a prediction, and no ansatz smuggled via prior work that itself only assumes the target. The axisymmetric resampling of the truth cubes is a modeling limitation the paper states explicitly; it does not make the skill score circular, because the same forward model generates both the inverted observables and the scoring field. Score 0 is therefore the correct finding.
Assumptions & free parameters
free parameters (7)
- global prior weight α_prior =
0.5
- external-pressure weight α =
0.05
- shear-layer weight β =
0.3
- anti-kink weight λ =
1
- divergence penalty weight =
0.05
- Adam learning rate and iteration count =
0.05 / 400
- lattice size N and domain half-width L =
96 / 30 rg
assumptions (4)
- domain assumption Synchrotron emissivity is proportional to magnetic energy density, ε∝|B|².
- ad hoc to paper With frame-dragging coefficient h=0, every retained prior term is even under S→−S.
- domain assumption The axisymmetric mean field after azimuthal averaging faithfully represents the geometry that dominates the projected observables W(z) and p(z).
- standard math Jet width is the FWHM-equivalent second-moment width of the transverse brightness profile; polarization degree is formed from line-of-sight integrated Stokes Q and U.
Cite this review
Pith. "Pith review of Hamiltonian variational reconstruction of the 3D magnetic geometry of relativistic jets: accuracy across GRMHD models and the fundamental sign degeneracy." pith.science (2026). https://pith.science/paper/JLLR7WNZ
@misc{pith2026260709569,
author = {Pith},
title = {Pith review of: Hamiltonian variational reconstruction of the 3D magnetic geometry of relativistic jets: accuracy across GRMHD models and the fundamental sign degeneracy},
year = {2026},
howpublished = {\url{https://pith.science/paper/JLLR7WNZ}},
note = {Machine review of arXiv:2607.09569}
}
read the original abstract
Every resolved image of a relativistic jet encodes its magnetic field, yet no image records which way the field points along the axis. Synchrotron intensity and linear polarization are blind to this: both probe the field only through even combinations of its components. We introduce H-MOG, a variational method that reconstructs the 3D field on a lattice from two projected observables, the jet width W(z) and the linear polarization p(z), regularized by a Hamiltonian prior and optimized with automatic differentiation. We apply H-MOG to ten GRMHD simulations spanning MAD and SANE states and five black-hole spins (a* = -0.94, -0.5, 0, +0.5, +0.94), and test three routes to break the intrinsic sign degeneracy of the reconstruction: a Faraday rotation-measure term, full 3D sampling, and a Blandford-Znajek spin prior. The unsigned field orientation is recovered at <|cos|> ~ 0.95-0.98, far above the random expectation of 0.5. The sense of the poloidal field, however, is not recovered, and we prove it cannot be: W and p are invariant under B -> -B. All three routes to break this degeneracy fail for distinct physical reasons; the spin-sense relation in these turbulent jets is not monotonic, differing sharply between MAD and SANE. Recovering the sense requires a parity-odd observable: Faraday tomography or circular polarization.
Figures
Figures from the paper (8 more)
Reference graph
Works this paper leans on
-
[1]
and Olivares, H
Porth, O. and Olivares, H. and Mizuno, Y. and Younsi, Z. and Rezzolla, L. and others , title=. ApJS , year=
-
[2]
Etienne, Z. B. and Paschalidis, V. and Haas, R. and Moesta, P. and Shapiro, S. L. , title=. Phys. Rev. D , year=
-
[3]
McKinney, J. C. and Tchekhovskoy, A. and Blandford, R. D. , title=. MNRAS , year=
-
[4]
and Nakamura, M
Asada, K. and Nakamura, M. , title=. ApJL , year=
-
[5]
and Nakamura, M
Asada, K. and Nakamura, M. , title=. A&A , year=
-
[6]
and Asada, K
Lu, R.-S. and Asada, K. and Krichbaum, T. P. and others , title=. Nature , year=
-
[7]
McKinney, J. C. and Tchekhovskoy, A. and Blandford, R. D. , title=. MNRAS Lett. , year=
-
[8]
and Porth, O
Olivares, H. and Porth, O. and Mizuno, Y. and Rezzolla, L. , title=. Comput. Phys. Commun. , year=
Show all 98 references
-
[9]
Bellan, P. M. , title=. J. Plasma Phys. , year=
-
[10]
Blandford, R. D. and Znajek, R. L. , title=. MNRAS , year=
-
[11]
and Narayan, R
Tchekhovskoy, A. and Narayan, R. and McKinney, J. C. , title=. MNRAS Lett. , year=
-
[12]
Boroson, T. A. , title=. ApJL , year=
-
[13]
Nikonov, A. S. and Kovalev, Y. Y. and Kravchenko, E. V. and others , title=. MNRAS , year=
-
[14]
and Kino, M
Hada, K. and Kino, M. and Giovannini, G. and others , title=. arXiv e-prints , year=. 2403.18359 , archivePrefix=
-
[15]
and Zhao, G.-Y
Cui, Y. and Zhao, G.-Y. and Li, Z. and others , title=. Nat. Astron. , year=
-
[16]
and Philippov, A
Parfrey, K. and Philippov, A. A. and Cerutti, B. , title=. Phys. Rev. Lett. , year=
-
[17]
Blandford, R. D. and Payne, D. G. , title=. MNRAS , year=
-
[18]
Komissarov, S. S. , title=. MNRAS , year=
-
[19]
and Konigl, A
Vlahakis, N. and Konigl, A. , title=. ApJ , year=
-
[20]
, title=
Lyubarsky, Y. , title=. ApJ , year=
-
[21]
, title=
Toth, G. , title=. J. Comput. Phys. , year=
-
[22]
and Kemm, F
Dedner, A. and Kemm, F. and Kroner, D. and Munz, C.-D. and Schnitzer, T. and Wesenberg, M. , title=. J. Comput. Phys. , year=
-
[23]
and Bodo, G
Mignone, A. and Bodo, G. and Massaglia, S. and others , title=. ApJS , year=
-
[24]
Gammie, C. F. and McKinney, J. C. and Toth, G. , title=. ApJ , year=
-
[25]
and Narayan, R
Yuan, F. and Narayan, R. , title=. ARA&A , year=
-
[26]
Zavala, R. T. and Taylor, G. B. , title=. ApJ , year=
-
[27]
Lister, M. L. and. MOJAVE. XVII. Parsec-Scale AGN Jet Kinematics Analysis Based on 19 Years of VLBA Observations at 15 GHz , journal=. 2018 , volume=
2018
-
[28]
Homan, D. C. and Lister, M. L. and Aller, M. F. and others , title=. ApJ , year=
-
[29]
First M87 Event Horizon Telescope Results. VII. Polarization of the Ring , journal=. 2021 , volume=
2021
-
[30]
First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near the Event Horizon , journal=. 2021 , volume=
2021
-
[31]
Broderick, A. E. and McKinney, J. C. , title=. ApJ , year=
-
[32]
and Igumenshchev, I
Narayan, R. and Igumenshchev, I. V. and Abramowicz, M. A. , title=. PASJ , year=
-
[33]
and Asada, K
Nakamura, M. and Asada, K. and Hada, K. and others , title=. ApJ , year=
-
[34]
Ressler, S. M. and Tchekhovskoy, A. and Quataert, E. and Gammie, C. F. , title=. MNRAS , year=
-
[35]
and Hawley, J
Beckwith, K. and Hawley, J. F. and Krolik, J. H. , title=. ApJ , year=
-
[36]
and McKinney, J
Tchekhovskoy, A. and McKinney, J. C. and Narayan, R. , title=. ApJ , year=
-
[37]
and Shibata, K
Koide, S. and Shibata, K. and Kudoh, T. and Meier, D. L. , title=. Science , year=
-
[38]
and Goodman, J
Li, L.-X. and Goodman, J. and Narayan, R. , title=. ApJ , year=
-
[39]
Spruit, H. C. , title=. Lect. Notes Phys. , year=
-
[40]
Meier, D. L. and Koide, S. and Uchida, Y. , title=. Science , year=
-
[41]
White, C. J. and Stone, J. M. and Gammie, C. F. , title=. ApJS , year=
-
[42]
Fishbone, L. G. and Moncrief, V. , title=. ApJ , year=
-
[43]
and Sadowski, A
Narayan, R. and Sadowski, A. and Penna, R. F. and Kulkarni, A. K. , title=. MNRAS , year=
-
[44]
and Falcke, H
Moscibrodzka, M. and Falcke, H. and Shiokawa, H. , title=. A&A , year=
-
[45]
Walker, R. C. and Hardee, P. E. and Davies, F. B. and others , title=. ApJ , year=
-
[46]
and Krichbaum, T
Kim, J.-Y. and Krichbaum, T. P. and Lu, R.-S. and others , title=. A&A , year=
-
[47]
and Doi, A
Hada, K. and Doi, A. and Kino, M. and others , title=. Nature , year=
-
[48]
Doeleman, S. S. and Fish, V. L. and Schenck, D. E. and others , title=. Science , year=
-
[49]
Gralla, S. E. and Lupsasca, A. and Marrone, D. P. , title=. Phys. Rev. D , year=
-
[50]
Gralla, S. E. and Lupsasca, A. , title=. Class. Quantum Grav. , year=
-
[51]
First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring , journal=. 2019 , volume=
2019
-
[52]
and Gomez, J
Mizuno, Y. and Gomez, J. L. and Nishikawa, K.-I. and others , title=. ApJ , year=
-
[53]
and Lyubarsky, Y
Mizuno, Y. and Lyubarsky, Y. and Beckwith, K. and Hardee, P. E. , title=. ApJ , year=
-
[54]
and Lobanov, A
Mertens, F. and Lobanov, A. P. and Walker, R. C. and Hardee, P. E. , title=. A&A , year=
-
[55]
and Narayan, R
Sadowski, A. and Narayan, R. and McKinney, J. C. and Tchekhovskoy, A. , title=. MNRAS , year=
-
[56]
Balbus, S. A. and Hawley, J. F. , title=. ApJ , year=
-
[57]
Komissarov, S. S. and Barkov, M. V. and Vlahakis, N. and Konigl, A. , title=. MNRAS , year=
-
[58]
McKinney, J. C. , title=. MNRAS , year=
-
[59]
and Hada, K
Giroletti, M. and Hada, K. and Giovannini, G. and others , title=. A&A , year=
-
[60]
and Lu, R.-S
Akiyama, K. and Lu, R.-S. and Fish, V. L. and others , title=. ApJ , year=
-
[61]
First M87 Event Horizon Telescope Results. IX. Detection of Near-horizon Circular Polarization , journal=. 2023 , volume=
2023
-
[62]
First Sagittarius A* Event Horizon Telescope Results. VII. Polarization of the Ring , journal=. 2024 , volume=
2024
-
[63]
A&A , year=
Horizon-scale variability of. A&A , year=
-
[64]
and Camilloni, F
Meringolo, C. and Camilloni, F. and Rezzolla, L. , title=. ApJL , year=. 2507.08942 , archivePrefix=
-
[65]
and Rodriguez, M
Lupsasca, A. and Rodriguez, M. J. and Strominger, A. , title=. JHEP , year=
-
[66]
and others , title=
Lupsasca, A. and others , title=. ApJ , year=
-
[67]
2026 , volume=
Spatially resolved polarization swings in the binary black hole candidate OJ 287 with first EHT observations , journal=. 2026 , volume=
2026
-
[68]
and Prather, B
Dhruv, V. and Prather, B. and Wong, G. N. and Gammie, C. F. , title=. ApJS , year=
-
[69]
L., et al
Akiyama, K., Lu, R.-S., Fish, V. L., et al. 2015, ApJ, 807, 150
2015
-
[70]
& Nakamura, M
Asada, K. & Nakamura, M. 2011, ApJL, 745, L28
2011
-
[71]
Blandford, R. D. & Payne, D. G. 1982, MNRAS, 199, 883
1982
-
[72]
Blandford, R. D. & Znajek, R. L. 1977, MNRAS, 179, 433
1977
-
[73]
Broderick, A. E. & McKinney, J. C. 2010, ApJ, 725, 750
2010
-
[74]
N., & Gammie, C
Dhruv, V., Prather, B., Wong, G. N., & Gammie, C. F. 2025, ApJS, 277, 16
2025
-
[75]
S., Fish, V
Doeleman, S. S., Fish, V. L., Schenck, D. E., et al. 2012, Science, 338, 355
2012
-
[76]
2019, ApJL, 875, L1
Event Horizon Telescope Collaboration . 2019, ApJL, 875, L1
2019
-
[77]
2021 a , ApJL, 910, L12
Event Horizon Telescope Collaboration . 2021 a , ApJL, 910, L12
2021
-
[78]
2021 b , ApJL, 910, L13
Event Horizon Telescope Collaboration . 2021 b , ApJL, 910, L13
2021
-
[79]
2023, ApJL, 957, L20
Event Horizon Telescope Collaboration . 2023, ApJL, 957, L20
2023
-
[80]
2024, ApJL, 964, L25
Event Horizon Telescope Collaboration . 2024, ApJL, 964, L25
2024
-
[81]
F., McKinney, J
Gammie, C. F., McKinney, J. C., & Toth, G. 2003, ApJ, 589, 444
2003
-
[82]
2011, Nature, 477, 185
Hada, K., Doi, A., Kino, M., et al. 2011, Nature, 477, 185
2011
-
[83]
C., Lister, M
Homan, D. C., Lister, M. L., Aller, M. F., et al. 2009, ApJ, 706, 1253
2009
-
[84]
P., Lu, R.-S., et al
Kim, J.-Y., Krichbaum, T. P., Lu, R.-S., et al. 2018, A&A, 616, A188
2018
-
[85]
McKinney, J. C. 2006, MNRAS, 368, 1561
2006
-
[86]
C., Tchekhovskoy, A., & Blandford, R
McKinney, J. C., Tchekhovskoy, A., & Blandford, R. D. 2012, MNRAS, 423, 3083
2012
-
[87]
P., Walker, R
Mertens, F., Lobanov, A. P., Walker, R. C., & Hardee, P. E. 2016, A&A, 595, A54
2016
-
[88]
Mizuno, Y., Lyubarsky, Y., Beckwith, K., & Hardee, P. E. 2014, ApJ, 780, 87
2014
-
[89]
2018, ApJ, 868, 146
Nakamura, M., Asada, K., Hada, K., et al. 2018, ApJ, 868, 146
2018
-
[90]
V., & Abramowicz, M
Narayan, R., Igumenshchev, I. V., & Abramowicz, M. A. 2003, PASJ, 55, L69
2003
-
[91]
F., & Kulkarni, A
Narayan, R., Sadowski, A., Penna, R. F., & Kulkarni, A. K. 2012, MNRAS, 426, 3241
2012
-
[92]
2019, ApJS, 243, 26
Porth, O., Olivares, H., Mizuno, Y., et al. 2019, ApJS, 243, 26
2019
-
[93]
C., & Tchekhovskoy, A
Sadowski, A., Narayan, R., McKinney, J. C., & Tchekhovskoy, A. 2014, MNRAS, 439, 503
2014
-
[94]
Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2011, MNRAS Lett., 418, L79
2011
-
[95]
Toth, G. 2000, J. Comput. Phys., 161, 605
2000
-
[96]
& Konigl, A
Vlahakis, N. & Konigl, A. 2004, ApJ, 605, 656
2004
-
[97]
C., Hardee, P
Walker, R. C., Hardee, P. E., Davies, F. B., et al. 2018, ApJ, 855, 128
2018
-
[98]
Zavala, R. T. & Taylor, G. B. 2004, ApJ, 612, 749
2004
Reviewed July 13, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.