REVIEW 2 major objections 3 minor 43 references
Gravitational Lensing in the Schwarzschild Spacetime: Photon Rings in Vacuum and in the Presence of a Plasma
T0 review · 2 major / 3 minor · reviewed 2026-05-22 · grok-4.3
Pith's one-line read Inhomogeneous plasma around a Schwarzschild black hole changes the size and shape of photon rings in a frequency-dependent way that multifrequency observations can use to measure plasma properties.
desk verdict The paper supplies analytic expressions for photon rings under one specific inhomogeneous plasma profile around Schwarzschild, but the separability of the ray equations needs explicit confirmation. 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
Analytic integration of the null geodesic equations for the specific inhomogeneous plasma density profile, which yields explicit relations between impact parameters and observer angles.
What would settle it
Multifrequency images of photon rings around a known-mass black hole that show no measurable change in ring radius or shape between frequencies would contradict the predicted effect of the inhomogeneous plasma model.
Extended reading notes
Core claim
For the chosen analytic plasma model the geodesic equations integrate in closed form. An orthonormal tetrad at the observer converts the constants of motion into latitude-longitude coordinates on the celestial sphere. The direct image and the first- and second-order photon rings therefore occupy frequency-dependent locations and shapes. These locations are compared with the vacuum and homogeneous-plasma limits, and the redshift and travel time are computed explicitly. The structural differences that appear are shown to be usable for extracting plasma properties from multifrequency data.
Load-bearing premise
The plasma electron density is given by a functional form that allows the geodesic equations to be integrated completely in closed form.
Editorial extensions
If this is right
- The first-order photon ring contracts or expands with observing frequency in a manner distinct from the vacuum case.
- The second-order ring exhibits an even larger relative shift, amplifying the observable signature of the density gradient.
- Redshift and travel-time differences between frequencies provide independent observables that can be combined with ring geometry.
- Comparison with homogeneous-plasma results isolates the effect of the equatorial density increase.
Reading between the lines
- The same analytic approach could be tested against numerical ray-tracing codes that use more realistic density profiles.
- If the model is approximately correct, existing black-hole images could be re-analyzed at multiple frequencies to place limits on plasma gradients.
- Extension to Kerr spacetime would be a natural next step to assess whether spin changes the frequency dependence.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript examines gravitational lensing by a Schwarzschild black hole surrounded by an inhomogeneous, pressureless, non-magnetized plasma whose electron density increases toward the equatorial plane near the horizon. It derives the plasma-modified equations of motion, claims to solve them analytically, places an observer in the outer domain, introduces an orthonormal tetrad to map constants of motion to celestial-sphere coordinates, computes the direct image and first- and second-order photon rings, writes a lens equation, and evaluates redshift and travel time. These quantities are compared with the vacuum and homogeneous-plasma cases to argue that multifrequency observations of ring-structure changes can extract plasma properties.
Significance. If the analytic solutions are valid, the work supplies an explicit, parameter-free example in which plasma inhomogeneity produces traceable modifications to photon-ring geometry on the observer’s sky. This is directly relevant to interpreting multifrequency EHT or ngEHT data on black-hole shadows and rings. The absence of fitted parameters and the explicit comparison to vacuum and homogeneous cases are strengths that enhance the falsifiability of the proposed extraction method.
major comments (2)
- [§3] §3 (derivation of equations of motion): the central claim that the chosen n_e(r,θ) permits fully analytic integration rests on the separability of the Hamilton-Jacobi equation once the ω_p²(r,θ)/ω² term is included. The manuscript must demonstrate explicitly that no residual r–θ coupling survives after the two Killing constants are imposed; otherwise the reported closed-form r(λ) and θ(λ) are at best approximate and the subsequent lens equation and ring-coordinate expressions in §4 lose their analytic traceability.
- [§4] §4 (photon rings and lens equation): the argument that structural changes in the first- and second-order rings can be used to extract plasma parameters via multifrequency observations presupposes that the analytic solutions of §3 remain valid across the relevant impact-parameter range. Without explicit checks against limiting cases (equatorial plane, large-r asymptotics, or homogeneous-plasma reduction) or error estimates on the integration, the extraction claim is not yet load-bearing.
minor comments (3)
- [Abstract] The abstract states that the equations are “derived and solved analytically” but does not indicate the domain of validity or any approximations; this should be stated explicitly.
- [§2] Clarify the precise definition of the orthonormal tetrad and the mapping from conserved quantities to observer latitude-longitude coordinates; a short appendix or figure would help.
- [§5] Travel-time and redshift expressions should be compared numerically to the vacuum case for at least one concrete set of plasma parameters to illustrate the magnitude of the effect.
Simulated Author's Rebuttal
We thank the referee for the careful reading of the manuscript and for highlighting its potential relevance to multifrequency EHT observations. We address each major comment below and will incorporate the requested clarifications and validations into a revised version.
read point-by-point responses
-
Referee: [§3] §3 (derivation of equations of motion): the central claim that the chosen n_e(r,θ) permits fully analytic integration rests on the separability of the Hamilton-Jacobi equation once the ω_p²(r,θ)/ω² term is included. The manuscript must demonstrate explicitly that no residual r–θ coupling survives after the two Killing constants are imposed; otherwise the reported closed-form r(λ) and θ(λ) are at best approximate and the subsequent lens equation and ring-coordinate expressions in §4 lose their analytic traceability.
Authors: We agree that an explicit verification of separability strengthens the presentation. The plasma density profile n_e(r,θ) was specifically constructed so that the term ω_p²(r,θ)/ω², when combined with the two Killing constants, produces an additive separation of the Hamilton-Jacobi equation into purely r-dependent and θ-dependent parts with no cross terms remaining. In the revised manuscript we will insert a short subsection in §3 that substitutes the conserved quantities, cancels all mixed r–θ contributions term by term, and arrives at the decoupled ordinary differential equations whose solutions are the reported closed-form expressions for r(λ) and θ(λ). revision: yes
-
Referee: [§4] §4 (photon rings and lens equation): the argument that structural changes in the first- and second-order rings can be used to extract plasma parameters via multifrequency observations presupposes that the analytic solutions of §3 remain valid across the relevant impact-parameter range. Without explicit checks against limiting cases (equatorial plane, large-r asymptotics, or homogeneous-plasma reduction) or error estimates on the integration, the extraction claim is not yet load-bearing.
Authors: We accept that additional validation improves the robustness of the extraction argument. The revised §4 will contain direct comparisons of the analytic ring coordinates and lens equation against three limiting cases: (i) the equatorial-plane restriction, (ii) the large-r asymptotic expansion, and (iii) the exact reduction to the homogeneous-plasma solution. We will also report the maximum relative deviation obtained when the analytic expressions are cross-checked against numerical integration of the geodesic equations over the impact-parameter interval relevant to the first- and second-order rings, thereby quantifying the accuracy of the closed-form results. revision: yes
Circularity Check
No significant circularity; derivations are self-contained from metric and plasma ansatz
full rationale
The paper begins with the Schwarzschild metric and a postulated analytic electron-density profile for the pressureless plasma, derives the modified geodesic equations, and states that this specific profile permits explicit integration in quadratures. It then constructs the observer tetrad, lens equation, and redshift expressions directly from those solutions without any fitted parameters, self-citation chains, or uniqueness theorems that reduce the outputs to the inputs by construction. The reported ring coordinates and structural changes are computed quantities, not redefinitions or statistical fits of the same data. This matches the default expectation of an independent derivation chain.
Assumptions & free parameters
assumptions (2)
- domain assumption Light propagation obeys the geometric-optics limit in a cold, pressureless, non-magnetised plasma whose refractive index depends only on electron density.
- standard math The background spacetime is exactly the Schwarzschild metric.
Cite this review
Pith. "Pith review of Gravitational Lensing in the Schwarzschild Spacetime: Photon Rings in Vacuum and in the Presence of a Plasma." pith.science (2026). https://pith.science/paper/XJRFTEYR
@misc{pith2026250800624,
author = {Pith},
title = {Pith review of: Gravitational Lensing in the Schwarzschild Spacetime: Photon Rings in Vacuum and in the Presence of a Plasma},
year = {2026},
howpublished = {\url{https://pith.science/paper/XJRFTEYR}},
note = {Machine review of arXiv:2508.00624}
}
read the original abstract
Astrophysical black holes are usually surrounded by an accretion disk. At least parts of these accretion disks consist of a plasma in which light rays with different energies are dispersed. However, we usually do not know the exact configurations of these plasmas. In this paper we will now use the example of a Schwarzschild black hole embedded in an inhomogeneous pressureless and nonmagnetised plasma to investigate how the structural changes of the photon rings can help us to determine the properties of a plasma surrounding a black hole using multifrequency observations. For this purpose we will use a simple analytic model which describes a plasma whose electron density increases towards the equatorial plane when we approach the event horizon. For the chosen model we will first derive and then analytically solve the equations of motion. Then we will place an observer in the domain of outer communication and introduce an orthonormal tetrad to relate the constants of motion to latitude-longitude coordinates on the observer's celestial sphere. In the next step we will use the analytic solutions to investigate the geometric structures of the direct image as well as the photon rings of first and second order on the observer's celestial sphere. We will write down a lens equation and calculate the redshift and the travel time. We will compare the obtained results to results for photon rings in vacuum and in the presence of a homogeneous plasma. Finally, we will discuss which of these quantities can be used to extract information about the properties of the plasma.
Figures
Figures from the paper (10 more)
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
For the chosen model we will first derive and then analytically solve the equations of motion... using Jacobi’s elliptic functions and Legendre’s elliptic integrals
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Schwarzschild spacetime... line element gµνdxµdxν = −P(r)dt² + dr²/P(r) + r²(dϑ² + sin²ϑ dφ²)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Roots of the r Motion For calculating the roots of the r motion we start with writing down the conditional equation. In the presence of a plasma with the distribution specified in (6) it reads r4 + 2mE2 C E2 − E2 C r3 − K E2 − E2 C r2 + 2mK E2 − E2 C r = 0. (A1) We can easily see that in the vacuum case with EC = 0 the second term vanishes and otherwise t...
-
[2]
As first step we take (40) and write down the conditional equation for the roots
Roots of the ϑ Motion for Light Rays in the Inhomogeneous Plasma Calculating the roots of the ϑ motion is straight for- ward. As first step we take (40) and write down the conditional equation for the roots. For light rays propa- gating through the inhomogeneous plasma described by the distribution (6) it is a biquadratic polynomial and reads x4 + K − 2ω2...
-
[3]
Elementary Integrals In total we need to evaluate four different elementary integrals. After applying the coordinate transformation (26) we can rewrite all of them in one of the following forms I1 = Z y yi dy′ (y′ − ya) √y′ − y1 , (B1) I2 = Z y yi dy′ (y′ − ya)2 √y′ − y1 , (B2) where in the first integral the parameter ya can be −K/12, yH, or yph, while i...
-
[4]
II C we encountered several elliptic integrals
Elliptic Integrals When we solved the equations of motion in Sec. II C we encountered several elliptic integrals. In this paper we will rewrite them in terms of elementary functions and Legendre’s elliptic integrals of the first, second, and third kind. However, before we use them to solve the equations of motion we need to define them. This is the purpos...
-
[5]
, First M87 Event Horizon Telescope results
The Event Horizon Telescope Collaboration et al. , First M87 Event Horizon Telescope results. I. The shadow of the supermassive black hole, Astrophys. J. Lett. 875, L1 (2019)
work page 2019
-
[6]
, First Sagittarius A ∗ Event Horizon Telescope results
The Event Horizon Telescope Collaboration et al. , First Sagittarius A ∗ Event Horizon Telescope results. I. The shadow of the supermassive black hole in the center of the Milky Way, Astrophys. J. Lett. 930, L12 (2022)
work page 2022
-
[7]
R. P. Kerr, Gravitational field of a spinning mass as an example of algebraically special metrics, Phys. Rev. Lett 11, 237 (1963)
work page 1963
-
[8]
, First M87 Event Horizon Telescope results
The Event Horizon Telescope Collaboration et al. , First M87 Event Horizon Telescope results. V. Physical ori- gin of the asymmetric ring, Astrophys. J. Lett. 875, L5 (2019)
work page 2019
Show all 43 references
-
[9]
, First M87 Event Horizon Telescope results
The Event Horizon Telescope Collaboration et al. , First M87 Event Horizon Telescope results. VIII. Magnetic field structure near the event horizon, Astrophys. J. Lett. 910, L13 (2021)
2021
-
[10]
R. A. Breuer and J. Ehlers, Propagation of high- frequency electromagnetic waves through a magnetized plasma in curved space-time. I, Proc. R. Soc. A 370, 389–406 (1980)
1980
-
[11]
R. A. Breuer and J. Ehlers, Propagation of high- frequency electromagnetic waves through a magnetized plasma in curved spaces-time. II. Application of the asymptotic approximation, Proc. R. Soc. A 374, 65–86 (1981)
1981
-
[12]
Perlick, Ray Optics, Fermat’s Principle, and Appli- cations to General Relativity , Lecture Notes in Physics Monographs (Springer, Berlin, 2000)
V. Perlick, Ray Optics, Fermat’s Principle, and Appli- cations to General Relativity , Lecture Notes in Physics Monographs (Springer, Berlin, 2000)
2000
-
[13]
O. Y. Tsupko and G. S. Bisnovatyi-Kogan, Gravitational lensing in plasma: Relativistic images at homogeneous plasma, Phys. Rev. D 87, 124009 (2013)
2013
-
[14]
Er and S
X. Er and S. Mao, Effects of plasma on gravitational lensing, Mon. Not. R. Astron. Soc. 437, 2180 (2014)
2014
-
[15]
Feleppa, V
F. Feleppa, V. Bozza, and O. Y. Tsupko, Strong deflec- tion limit analysis of black hole lensing in inhomogeneous plasma, Phys. Rev. D 110, 064031 (2024)
2024
-
[16]
Perlick, O
V. Perlick, O. Y. Tsupko, and G. S. Bisnovatyi-Kogan, Influence of a plasma on the shadow of a spherically sym- metric black hole, Phys. Rev. D 92, 104031 (2015)
2015
-
[17]
Perlick and O
V. Perlick and O. Y. Tsupko, Light propagation in a plasma on Kerr spacetime: Separation of the Hamilton- Jacobi equation and calculation of the shadow, Phys. Rev. D 95, 104003 (2017)
2017
-
[18]
Bezdˇ ekov´ a, V
B. Bezdˇ ekov´ a, V. Perlick, and J. Biˇ c´ ak, Light propaga- tion in a plasma on an axially symmetric and stationary spacetime: Separability of the Hamilton–Jacobi equation and shadow, J. Math. Phys. 63, 092501 (2022)
2022
-
[19]
Perlick and O
V. Perlick and O. Y. Tsupko, Light propagation in a plasma on Kerr spacetime. II. Plasma imprint on pho- ton orbits, Phys. Rev. D 109, 064063 (2024)
2024
-
[20]
O. Y. Tsupko, Deflection of light rays by a spherically symmetric black hole in a dispersive medium, Phys. Rev. D 103, 104019 (2021)
2021
-
[21]
Luminet, Image of a spherical black hole with thin accretion disk, Astron
J.-P. Luminet, Image of a spherical black hole with thin accretion disk, Astron. Astrophys. 75, 228 (1979)
1979
-
[22]
Johannsen and D
T. Johannsen and D. Psaltis, Testing the no-hair theorem with observations in the electromagnetic spectrum. II. Black hole images, Astrophys. J. 718, 446 (2010)
2010
-
[23]
Aratore, O
F. Aratore, O. Y. Tsupko, and V. Perlick, Constraining spherically symmetric metrics by the gap between photon rings, Phys. Rev. D 109, 124057 (2024)
2024
-
[24]
A. E. Broderick et al. , The photon ring in M87 ∗, Astro- phys. J. 935, 61 (2022). 39
2022
-
[25]
Lockhart and S
W. Lockhart and S. E. Gralla, How narrow is the M87 ∗ ring – II. A new geometric model, Mon. Not. R. Astron. Soc. 517, 2462–2470 (2022)
2022
-
[26]
M. D. Johnson et al. , The Black Hole Explorer: moti- vation and vision, in Space Telescopes and Instrumenta- tion 2024: Optical, Infrared, and Millimeter Wave , Vol. 13092, edited by L. E. Coyle, S. Matsuura, and M. D. Perrin, International Society for Optics and Photonics (S...
2024
-
[27]
Falcke, F
H. Falcke, F. Melia, and E. Agol, Viewing the shadow of the black hole at the galactic center, Astrophys. J. 528, L13 (2000)
2000
-
[28]
Desire, A
T. Desire, A. C´ ardenas-Avenda˜ no, and A. Chael, Mul- tifrequency models of black hole photon rings from low- luminosity accretion disks, Astrophys. J.980, 262 (2025)
2025
-
[29]
Kobialko, D
K. Kobialko, D. Gal’tsov, and A. Molchanov, Gravita- tional shadow and emission spectrum of thin accretion disks in a plasma medium, arXiv:2505.07993 (2025)
2025
-
[30]
Perlick, Influence of a plasma on the shadow of a black hole, Astron
V. Perlick, Influence of a plasma on the shadow of a black hole, Astron. Rep. 67, S102 (2023)
2023
-
[31]
Mino, Perturbative approach to an orbital evolution around a supermassive black hole, Phys
Y. Mino, Perturbative approach to an orbital evolution around a supermassive black hole, Phys. Rev. D 67, 084027 (2003)
2003
-
[32]
S. E. Gralla and A. Lupsasca, Null geodesics of the Kerr exterior, Phys. Rev. D 101, 044032 (2020)
2020
-
[33]
T. C. Frost, Gravitational lensing of massive particles in the charged NUT spacetime, Phys. Rev. D 108, 124019 (2023)
2023
-
[34]
Hancock, Elliptic Integrals, 1st ed., edited by M
H. Hancock, Elliptic Integrals, 1st ed., edited by M. Mer- riman and R. S. Woodward, Mathematical Monographs (John Wiley & Sons, New York, 1917)
1917
-
[35]
Grenzebach, V
A. Grenzebach, V. Perlick, and C. L¨ ammerzahl, Photon regions and shadows of accelerated black holes, Int. J. Mod. Phys. D 24, 1542024 (2015)
2015
-
[36]
R. M. Wald, General Relativity (The University of Chicago Press, Chicago, 1984)
1984
-
[37]
A. Bohn, W. Throwe, F. H´ ebert, K. Henriksson, D. Bunandar, M. A. Scheel, and N. W. Taylor, What does a binary black hole merger look like?, Classical Quantum Gravity 32, 065002 (2015)
2015
-
[38]
Frittelli and E
S. Frittelli and E. T. Newman, Exact universal gravita- tional lensing equation, Phys. Rev. D 59, 124001 (1999)
1999
-
[39]
Perlick, Exact gravitational lens equation in spheri- cally symmetric and static spacetimes, Phys
V. Perlick, Exact gravitational lens equation in spheri- cally symmetric and static spacetimes, Phys. Rev. D 69, 064017 (2004)
2004
-
[40]
Bezanson, A
J. Bezanson, A. Edelman, S. Karpinski, and V. B. Shah, Julia: A fresh approach to numerical computing, SIAM Rev. 59, 65 (2017)
2017
-
[41]
M¨ uller and M
A. M¨ uller and M. Camenzind, Relativistic emission lines from accreting black holes - The effect of disk truncation on line profiles, Astron. Astrophys. 413, 861 (2004)
2004
-
[42]
L. M. Milne-Thomson, Elliptic Integrals, in Handbook of Mathematical Functions With Formulas, Graphs, and Mathematical Tables, Applied Mathematics Series, edited by M. Abramowitz and I. A. Stegun (U. S. Dept. of Com- merce, National Bureau of Standards, Washington D. C.,
-
[43]
10th ed., pp. 587–607
Reviewed May 22, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.