{"id":"2bba67ab-486d-446b-ba63-332235c8dd82","arxiv_id":"2505.07993","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A blackbody accretion disk viewed through transparent, rotating plasma produces frequency-dependent shadows and emission maps, with the frequency of the brightest total flux controlled mainly by viewing angle.","lead":"This paper computes what a glowing disk around a black hole would look like when the black hole is wrapped in thin plasma that bends and blocks light differently at different frequencies. It shows that the disk's image and brightness peaks change strongly with observing frequency, which could help future multi-frequency black hole telescopes measure the viewing angle.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Transparent-plasma assumption (eta=0, J=0) is the load-bearing weakness: the claimed ngEHT diagnostics in Figs. 15-17 are pure-transmission results, yet the paper targets frequencies where sources become optically thick.","rationale":"The central mathematical claim - that Eq. (17) gives the observed specific intensity of a blackbody disk seen through a non-absorbing, non-emitting, co-rotating cold plasma - is coherent; I checked the derivation from T conservation, the v_ph v_g=1 identity, and the redshift factor, and found no algebraic error that would break it. The typo in Eq. (77) (r^sigma should be sigma r) is real but non-load-bearing: the positivity intervals and Fig. 3 use the correct form, and the same correction is implicit in Eqs. (78)-(81). What is genuinely load-bearing is the physical use of the formula: the paper presents secondary images and flux-maximum positions as ngEHT diagnostics, but these predictions are computed entirely in pure transmission. At the frequencies the paper names in its introduction, the surrounding plasma is known to become optically thick, so eta=0 and J=0 are not safe defaults. This does not make the paper internally inconsistent; it makes the central observable claim conditional on an unquantified assumption. The reader already assigned CONDITIONAL for essentially this reason, and I see no reason to move the verdict, so I keep UNCHANGED. The proposed concrete test, adding a simple kappa_nu to the same pipeline, would settle whether the absorption concern actually shifts the claimed diagnostics.","tokens_in":1006,"tokens_out":1381,"duration_ms":155133,"concrete_test":"Re-run the authors' ray tracer with a simple two-fluid absorption: set the invariant extinction coefficient kappa_nu(r)=kappa_0 [omega_p(r)^2 / omega_s^2] with omega_s=(1+z)omega_o the local frequency and leave J=0 initially, tuning kappa_0 so the line-of-sight optical depth is about 1 at 230 GHz for M87*-like parameters; repeat for sigma=2,4,6, bar_theta=84 and 45 degrees, dot_M=0.1, and compare F(omega_o)/F_Sch(omega_o) with Figs. 15-17. If the peak frequency shifts by more than the peak width or the F>1 enhancement disappears, the inclination-angle diagnostic is absorption-sensitive; if the changes are within line width, the transparency idealization is benign. A second run with J chosen so the plasma re-emits locally (tau much less than 1 but nonzero) would test shadow filling.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing weakness is the transparency assumption used to derive Eq. (17) and to produce every image and flux curve in Section V. Equation (17) follows from dT/dlambda = -eta*T + J (Eq. 12) only when eta=0 and J=0, and the paper sets both to zero without quantitative justification. The paper itself (Introduction) states that at the ngEHT frequencies 86/230/345 GHz real sources 'transit from optically thin to optically thick', so the band invoked is exactly where tau can be order unity. The claimed observables - secondary-image visibility (Figs. 7-10), shadow emergence as a function of omega_o, and the F(omega_o)/F_Sch(omega_o) maxima whose peak position is said to measure inclination (Figs. 15-17) - are all pure-transmission ray-tracing results. Adding absorption would add an exp(-tau) suppression to secondary images, which are the origin of the F/F_Sch>1 peaks, and would add intrinsic plasma emission J that can fill the shadow, so the central diagnostic is not established for the target observations. This is a correctness risk for the astrophysical claim, not an internal inconsistency: the mathematical derivation of Eq. (17) under eta=J=0 is coherent.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper develops a ray-tracing framework for thermal radiation from a Novikov-Thorne thin accretion disk in Schwarzschild spacetime, surrounded by a cold, non-magnetized, pressureless, transparent plasma in stationary rotation. The central formal result is Eq. (17), an observed specific intensity that contains a plasma factor (1 - \\bar{\\omega}_p^2/\\omega_o^2) and a redshifted Planck factor. The authors derive analytical shadow radii for power-law plasma profiles, compute frequency-resolved images, total integrated intensities, and total-flux spectra for different inclinations, accretion rates, and disk sizes, and propose that the frequency of the total-flux maximum can serve as an inclination diagnostic for ngEHT observations.","tokens_in":17804,"tokens_out":17021,"duration_ms":149987,"significance":"If the results hold, the paper provides a useful semi-analytical framework that connects plasma-modified ray propagation with thermal disk emission, and the proposed observables (secondary-image visibility near the plasma equilibrium frequency and flux-maximum frequency shifts) are interesting. The derivation of Eq. (17) from the invariant transport equation is internally coherent under the stated assumptions, and the sigma=2 shadow formula (84) is an elegant closed-form result. The numerical study covers a wide parameter space and the paper ships clearly presented figures. The main caveats are the assumed plasma transparency and the errors in two analytical displays; these affect the astrophysical interpretation and must be corrected before the ngEHT claims are accepted.","major_comments":[{"comment":"The shadow radius is printed as R^2 = r^3(r^sigma - 2M(sigma+1)) / ((sigma-2)(r-2M)^2), but the correct algebra from Eq. (62) gives R^2 = r^3(sigma r - 2M(sigma+1)) / ((sigma-2)(r-2M)^2). The typo is not cosmetic: for sigma=0 the printed numerator is dimensionally inconsistent (r^sigma=1), and the allowed intervals (78)-(81) as well as Figure 3 are wrong for sigma values other than 1 and 2. Please correct the formula and recompute the affected analytical panels.","section":"Section V, Eq. (77)"},{"comment":"The displayed chain of inequalities is algebraically false. For alpha=gamma=1, rho=-5, and Omega_p=0.2, the second inequality in (55) states -5.76 <= -6.51, which is false. The intended conclusion (omega_p^2/omega^2 <= 1) is nevertheless correct and follows directly from (1-alpha gamma^{-1} rho^2)(1-gamma Omega_p^2/alpha) <= (1-Omega_p rho)^2, so the proof should be rewritten; as printed, the claim that the condition omega >= omega_p is automatically satisfied is not established.","section":"Section III, Eq. (55)"},{"comment":"The model sets the absorption and emission coefficients to zero (eta=0, J=0), yet the Introduction targets the 86-345 GHz band where Sgr A* and M87* are said to transit from optically thin to optically thick. Absorption would introduce a factor exp(-tau) into the secondary-image contributions that produce the F/F_Sch>1 peaks, and intrinsic plasma emission J would tend to fill the shadow. No estimate of the optical depth is given anywhere in the paper. The formal derivation is internally consistent, but the ngEHT diagnostic claim is not yet established for the targeted sources; please either include optical-depth estimates that justify transparency in the relevant regions or explicitly restrict the conclusions to optically thin configurations.","section":"Section II (after Eq. (16)) and Section V (Figs. 15-17)"}],"minor_comments":[{"comment":"The text reads 'at frequency omega'_o = h omega'_o'; this is a typo and should be 'omega'_o = h omega_o'.","section":"Section V, after Eq. (75)"},{"comment":"The sentence 'T'_D = T_D' conflicts with the rescaling F'_D = h^4 F_D; from Eq. (67) one obtains T'_D = h T_D. The ratios in Eq. (76) are consistent with the latter statement, so the sentence should be corrected.","section":"Section V, after Eq. (75)"},{"comment":"T_D is an intrinsic disk temperature, so writing T_D(omega_o) is misleading; the observed-frequency dependence enters through the redshift and through ray tracing, not through T_D itself. Please clarify the notation.","section":"Section IV, Eq. (69)"},{"comment":"No convergence study is reported for the RK4 step size or for the angular resolution of the stereographic grid. Since the positions of flux maxima are quantitative claims, a resolution check should be added.","section":"Section V"},{"comment":"Since Eq. (17) inherits the validity of the blackbody-in-plasma distribution, a short derivation of Eq. (10), or at least a statement of its regime of validity, would make the paper more self-contained.","section":"Section II, Eq. (10)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the central transport derivation is sound under its stated assumptions. However, the typo in Eq. (77), the invalid inequality in Eq. (55), and the unquantified transparency assumption are substantive enough that the paper should not be accepted in its current form. The corrections are local and the framework is worth publishing after revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a coherent, genuinely new calculation of frequency-resolved disk images and spectra in a dispersive plasma, and the central equation (17) is a clean transport result. The combination of a Novikov-Thorne disk with full radiative transfer in a cold, rotating plasma is not in the cited literature, and the paper delivers it cleanly. The coordinate-rotation trick that reduces all trajectories to the equatorial plane is elegant and makes the numerics fast, and the proof that the condition ω ≥ ω_p is automatic on the equations of motion is a useful add. The analytical shadow formulas, including the σ = 2 case and the equilibrium frequency, are solid under the stated plasma model. The paper is also honest about its main simplification: absorption and intrinsic emission are set to zero.\n\nThe load-bearing soft spot is exactly that transparency assumption. Equation (17) follows only for η = J = 0, and the introduction itself notes that at the ngEHT band (86, 230, 345 GHz) sources like Sgr A* and M87* transit from optically thin to optically thick. The secondary images that produce the F/F_Sch > 1 peaks in Figs. 15–17 are pure-transmission features; they will be exponentially suppressed if absorption is included, and intrinsic plasma emission can fill the shadow. So the claim that the flux-maximum frequency provides a direct inclination diagnostic is not established for the target observations. This is not a flaw in the mathematics, but it is a real flaw in the leap to an observational prediction. The paper needs at least a quantitative estimate of optical depth at those frequencies, or a clear statement that the results are a proof of concept only.\n\nMechanical issues: Eq. (77) as printed looks dimensionally off for general σ – probably r^σ should be σ r or something similar. There are no convergence tests and no code or data shipped, so the numerical results are not independently checkable as they stand. These are minor to moderate and fixable.\n\nWho this is for: people computing black hole images in dispersive media, and theorists planning multi-frequency ngEHT analyses. It deserves a serious referee. I would send it to review with a request to address the optical-depth question and fix the typo. My own verdict would be a conditional accept, not a reject; the math under the stated assumptions holds up, but the astrophysical interpretation needs more work.","headline":"Useful new calculation of plasma-transported disk spectra, but the ngEHT diagnostic rests on an unexamined transparency assumption.","tokens_in":18429,"tokens_out":3402,"would_cite":true,"duration_ms":34591,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","83C10","85A25"],"pacs":["04.70.Bw","95.30.Jx","98.62.Mw"],"model":"deepseek-v4-flash","headline":"A thin thermal disk around a Schwarzschild black hole, seen through transparent plasma, produces frequency-dependent shadow images whose total-flux peak is set mainly by the observer's inclination angle.","keywords":["black hole shadow","thin accretion disk","cold plasma","gravitational lensing in plasma","radiative transfer","frequency-dependent images","multi-frequency black hole imaging"],"falsifier":"Point one of the currently imaged supermassive black holes at the three planned radio bands (86, 230, and 345 GHz). The model predicts a specific frequency scaling of the shadow radius (for the $\\sigma=2$ plasma profile, $R^2(\\omega_o)=M^2(27-\\omega_o^{-2})$), no shadow at all below a calculable equilibrium frequency, and a total-flux peak whose frequency tracks the viewing angle rather than the accretion rate; observing a shadow below the equilibrium frequency, or a flux peak that moves with the accretion rate, would break the claim. The cleanest numerical check is to rerun the same ray tracing with absorption and intrinsic emission switched on and see whether the peak-frequency diagnostic survives.","tokens_in":17314,"feed_emoji":"🕳️","tokens_out":19522,"duration_ms":154588,"temperature":0.7,"pith_summary":"The paper builds a unified model of black-hole shadow and disk emission in a plasma: a thin thermal accretion disk around a Schwarzschild black hole, immersed in cold, pressureless, transparent plasma that co-rotates with the disk. It derives a modified blackbody intensity formula for light that has travelled through the plasma, and shows that the observed image depends strongly on frequency, so the same disk looks different at different wavelengths. If the model holds, multi-frequency observations of black holes can read the inclination angle of the disk off the frequency at which the total received flux peaks, and can spot secondary and relativistic images that a single-frequency observation would miss. This is the kind of signal the next-generation horizon-resolving observations at 86, 230, and 345 GHz are being built to see.","feed_headline":"One spectral peak reveals a black hole's tilt","feed_subtitle":"Photons bend differently in disk plasma, so the peak of the received flux pins down the viewing angle.","key_machinery":"The carrying mechanism is the invariant distribution function of the radiation, $T=\\omega^{-3}v_g^{-2}I(\\omega)$, which the transport equation $dT/d\\lambda=-\\eta T+J$ keeps constant along a ray whenever the plasma neither absorbs nor emits ($\\eta=0$, $J=0$). Inserting the equilibrium distribution of blackbody radiation in a plasma and using the identity $v_{ph}v_g=1$ that follows from the cold-plasma dispersion relation $k^2=\\omega^2-\\omega_p^2$, the plasma-frequency factors cancel and the transport law becomes the observed-intensity formula (17). The frequency-dependent ray paths themselves are governed by the effective potential $V=\\beta^{-1}(\\alpha^{-1}-\\gamma^{-1}\\rho'^2-\\omega_p^2/\\omega_\\infty^2)$: it is the $\\omega_o$-dependence entering through this potential that makes the shadow, the secondary images, and the flux spectra vary with wavelength. Around that core, the paper adds two analytic tools: coordinate rotations that reduce every photon orbit to an equatorial trajectory, and the equilibrium radius $r_e=2M(\\sigma+1)/\\sigma$ at which gravitational attraction and plasma repulsion balance, erasing the shadow.","core_discovery":"The paper's central claim is that the specific intensity received from a blackbody thin accretion disk through transparent cold plasma is $$I_o(\\omega_o)=\\frac{\\$omega_o^{3}$}{4\\$pi^{3}$}\\,\\frac{1-\\bar\\$omega_p^{2}$/\\$omega_o^{2}$}{$e^{{(1+z)\\omega_o/T_D}}$-1},$$ with $\\bar\\omega_p$ the plasma frequency at the observer, $1+z$ the redshift accumulated along the ray, and $T_D$ the local disk temperature. The plasma density at the emission point cancels out of this formula, because the product of phase and group velocity in a cold plasma equals one ($v_{ph}v_g=1$), yet the observed spectrum is not a rescaled blackbody: rays of different frequencies move in different effective potentials, so the redshift, the source temperature, and even whether a ray hits the disk at all depend on $\\omega_o$. Ray tracing in Schwarzschild spacetime then shows the shadow radius shrinking as frequency rises, the shadow vanishing completely below an equilibrium frequency at which the plasma's repulsive force balances gravity, and secondary images of the near disk edge appearing near that frequency; the computed flux at each frequency, normalized to the vacuum value, develops relative maxima that can exceed one, and the position of the dominant maximum is controlled mainly by the observer's inclination angle.","pith_inferences":["Editorial inference: the same transport machinery should transfer to a rotating black hole, where the flux-maximum frequency would plausibly encode the spin as well as the inclination angle, giving multi-frequency campaigns a second clean observable.","Editorial inference: the factor $(1-\\bar\\omega_p^2/\\omega_o^2)$ inside the intensity formula makes the low-frequency turnover of the observed spectrum a direct probe of the plasma frequency at the observer, so it could be used to estimate the line-of-sight electron density toward the source.","Editorial inference: the absorption-free predictions should be most trustworthy on the optically thin side of the band; near the frequencies where the source turns optically thick, rays linger in dense plasma, so the secondary-image and peak-frequency diagnostics would need to be recomputed with absorption included."],"forward_implications":["Below an equilibrium frequency $\\omega_e$, the gravitational shadow disappears entirely: plasma repulsion reflects every geodesic before it reaches the horizon, so a low-frequency observation shows disk light but no shadow.","Secondary images of the near edge of the disk become observable close to $\\omega_e$, especially for steep plasma profiles ($\\sigma=4,6$), whereas at higher frequencies the images approach their vacuum shapes.","The frequency position of the total-flux maximum is set mainly by the observer's inclination angle, giving a geometric diagnostic that does not depend on the poorly constrained mass accretion rate, which instead controls the peak's amplitude.","Normalized total intensity is independent of the accretion rate in vacuum but not in plasma, because dispersion makes the redshift and the source temperature frequency-dependent; the plasma spreads the same total energy over a larger image area, suppressing the peak intensity.","Because the shadow is frequency-dependent in plasma but frequency-independent in vacuum, multi-frequency imaging in the 86\\u2013345 GHz band can separate plasma effects from spacetime-geometry effects on the shadow boundary."],"supporting_citations":[{"why":"Supplies the equilibrium distribution function of blackbody radiation in a plasma (Eq. 10), the external input from which the observed-intensity formula (17) is derived.","marker":"[43, 44]"},{"why":"Supplies the relativistic transport equation in dispersive media, including the invariant distribution function that the paper propagates along rays.","marker":"[46]"},{"why":"Supplies the Hamiltonian for photon motion in cold plasma whose equations of motion the ray tracing integrates.","marker":"[11, 46]"},{"why":"Defines the plasma frequency in terms of electron number density, fixing the frequency-dependent refractive properties of the medium.","marker":"[9]"},{"why":"Provides the standard thin accretion disk model whose radial flux profile sets the local blackbody temperature in the intensity formula.","marker":"[41, 42]"},{"why":"The earlier discussion of plasma reflection of light that the paper's equilibrium-radius and shadow-erasure analysis builds on.","marker":"[21]"},{"why":"Shows the cold-plasma Hamiltonian does not depend on the plasma velocity, which licenses the co-rotating plasma model used here.","marker":"[18]"},{"why":"Defines the relativistic and secondary images of the disk that the paper tracks through the plasma at different frequencies.","marker":"[50, 56]"},{"why":"The vacuum shadow construction the plasma shadow formula is compared with and reduces to in the no-plasma limit.","marker":"[54]"}],"fun_headline_variants":["Frequency-dependent shadows reveal black hole tilt","Plasma shifts black hole spectrum, exposing inclination","Peak in disk spectrum traces observer's angle","Plasma controls black hole shadow size and flux peak"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on assuming the plasma around the disk is perfectly transparent, absorbing no light and emitting none of its own, so that its only effect on the radiation is to bend rays and shift frequencies; once realistic absorption in the 86\\u2013345 GHz band is switched on, the predicted images, secondary-image visibility, and flux-peak frequencies all change.","fun_headline_variants_meta":{"raw":{"variants":["Frequency-dependent shadows reveal black hole tilt","Plasma shifts black hole spectrum, exposing inclination","Peak in disk spectrum traces observer's angle","Plasma controls black hole shadow size and flux peak"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00093,"raw_usage":{"total_tokens":3977,"prompt_tokens":937,"completion_tokens":3040,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":553,"completion_tokens_details":{"reasoning_tokens":2982}},"tokens_in":553,"tokens_out":3040,"duration_ms":20514,"temperature":1.0,"reasoning_tokens":2982,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:08:18.379331+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Point one of the currently imaged supermassive black holes at the three planned radio bands (86, 230, and 345 GHz). The model predicts a specific frequency scaling of the shadow radius (for the $\\sigma=2$ plasma profile, $R^2(\\omega_o)=M^2(27-\\omega_o^{-2})$), no shadow at all below a calculable equilibrium frequency, and a total-flux peak whose frequency tracks the viewing angle rather than the accretion rate; observing a shadow below the equilibrium frequency, or a flux peak that moves with the accretion rate, would break the claim. The cleanest numerical check is to rerun the same ray tracing with absorption and intrinsic emission switched on and see whether the peak-frequency diagnostic survives.","supporting_citations":[{"cited_title":"Quasi-local studies of the particle surfaces and their stability in general spacetimes","cited_arxiv_id":"2208.03661","evidence_quote":"The earlier discussion of plasma reflection of light that the paper's equilibrium-radius and shadow-erasure analysis builds on."}],"review_version":1}