Gravitational Lensing in the Schwarzschild Spacetime: Photon Rings in Vacuum and in the Presence of a Plasma
Pith reviewed 2026-05-22 00:32 UTC · model grok-4.3
The pith
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.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
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.
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.
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.
Where Pith is 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.
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.
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.
Figures
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.
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.
Axiom & Free-Parameter Ledger
axioms (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.
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
-
[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)
work page 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)
work page 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)
work page 1981
-
[12]
V. Perlick, Ray Optics, Fermat’s Principle, and Appli- cations to General Relativity , Lecture Notes in Physics Monographs (Springer, Berlin, 2000)
work page 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)
work page 2013
- [14]
-
[15]
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)
work page 2024
-
[16]
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)
work page 2015
-
[17]
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)
work page 2017
-
[18]
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)
work page 2022
-
[19]
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)
work page 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)
work page 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)
work page 1979
-
[22]
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)
work page 2010
-
[23]
F. Aratore, O. Y. Tsupko, and V. Perlick, Constraining spherically symmetric metrics by the gap between photon rings, Phys. Rev. D 109, 124057 (2024)
work page 2024
-
[24]
A. E. Broderick et al. , The photon ring in M87 ∗, Astro- phys. J. 935, 61 (2022). 39
work page 2022
-
[25]
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)
work page 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 (SPIE, 2024) p. 130922D
work page 2024
- [27]
- [28]
-
[29]
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)
-
[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)
work page 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)
work page 2003
-
[32]
S. E. Gralla and A. Lupsasca, Null geodesics of the Kerr exterior, Phys. Rev. D 101, 044032 (2020)
work page 2020
-
[33]
T. C. Frost, Gravitational lensing of massive particles in the charged NUT spacetime, Phys. Rev. D 108, 124019 (2023)
work page 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)
work page 1917
-
[35]
A. Grenzebach, V. Perlick, and C. L¨ ammerzahl, Photon regions and shadows of accelerated black holes, Int. J. Mod. Phys. D 24, 1542024 (2015)
work page 2015
-
[36]
R. M. Wald, General Relativity (The University of Chicago Press, Chicago, 1984)
work page 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)
work page 2015
-
[38]
S. Frittelli and E. T. Newman, Exact universal gravita- tional lensing equation, Phys. Rev. D 59, 124001 (1999)
work page 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)
work page 2004
-
[40]
J. Bezanson, A. Edelman, S. Karpinski, and V. B. Shah, Julia: A fresh approach to numerical computing, SIAM Rev. 59, 65 (2017)
work page 2017
-
[41]
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)
work page 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
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.