{"id":"cb6c628b-686b-42ef-837b-40bd7a32eb22","arxiv_id":"2412.08359","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Explicit analytical formulas convert synchrotron Stokes parameters from the local disk frame to the observer frame around a Schwarzschild black hole, including disk thickness and radial accretion motion.","lead":"This paper derives fast analytical formulas for how synchrotron light from gas near a black hole changes its brightness, polarization degree, and polarization angle on the way to an observer. The new expressions cover matter that sits above the disk plane and moves inward as well as around, two effects earlier analytical models left out.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Main-text Eq. (39) for chi_SR does not match the general vector derivation in Appendix A; it is only valid for equatorial circular motion, contradicting the paper's claimed generalization.","rationale":"The reader's weakest assumption concerned the accuracy of the Poutanen (2020) light-bending approximation for elevated emitters. That is a legitimate external-validation concern, but it is secondary. The most load-bearing issue is internal: the main-text formula for chi_SR, which is part of the paper's central stated result, does not follow from the vector derivation in Appendix A except in the previously known equatorial circular limit. This is not a matter of numerical precision; it is an algebraic inconsistency in the claimed generalization. Because the correct derivation is present in the appendix, the paper may be repairable by replacing Eq. (39) with the coordinate form involving tilde{b}, and the figures may be unaffected if computed from the appendix expression. Nevertheless, as submitted, the central claim is misstated. I therefore retain the reader's CONDITIONAL verdict, but recommend that acceptance be conditioned on resolving this inconsistency and regenerating any affected figures, rather than only on a numerical error budget for the bending approximation. My disagreement with the reader is about which concern carries the most weight, not about whether the paper needs revision.","tokens_in":26190,"tokens_out":10479,"duration_ms":100754,"concrete_test":"Evaluate both expressions for a specific new-regime case, e.g. eta = 20 deg, sigma = 90 deg, r = 3, phi = 180 deg, i = 60 deg. Compute tan(chi_SR) from Eq. (39) (with cos(alpha)) and from Eq. (A.8) (with tilde{b}, Eq. (26)), using the same angles from Eqs. (5), (7), (8), (13). If the two values differ, the text is internally inconsistent; then check which value matches the contours in Fig. 2 or Fig. 5. A simpler check: for sigma = 180 deg and eta = 0, verify numerically that tilde{b} = -sin(alpha) sin(i) sin(phi) / sin(psi) is not equal to cos(alpha) = [sin(alpha) cos(psi) + sin(psi-alpha)] / sin(psi).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Eq. (39) gives the special-relativistic rotation of the polarization angle for arbitrary elevation eta and velocity direction sigma. However, Eq. (39) reads tan(chi_SR) = -beta cos(alpha) cos(zeta) / (sin^2(zeta) - beta cos(xi)), while the general vector expression derived in Appendix A, Eq. (A.8), is tan(chi_SR) = -beta (k0 dot (u x n)) (k0 dot n) / (1 - (k0 dot n)^2 - beta (k0 dot u)). Using the paper's definitions, k0 dot (u x n) is the quantity tilde{b} defined in Eq. (26), not cos(alpha). The equality tilde{b} = cos(alpha) holds only in the previously treated case eta = 0 and sigma = 90 deg, because then u x n = r. For the new regimes the paper claims to cover, such as sigma = 180 deg (radial infall) or eta != 0, tilde{b} and cos(alpha) differ. Thus Eq. (39) is not the specialization of Eq. (A.8) to the general geometry; the headline expression is internally inconsistent with the paper's own derivation. If the figures were computed with the correct appendix expression, the text requires a correction; if they were computed using Eq. (39), the results for elevated and non-circular cases are wrong. This is a more direct threat to the central claim than the external accuracy of the Poutanen bending approximation, which may still be adequate.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the analytical ray-tracing formalism artpol to synchrotron emission from matter at arbitrary elevation η above the equatorial plane and with non-circular velocity direction σ, in Schwarzschild spacetime. The central result is a decomposition of the total polarization-angle rotation into three additive contributions, χ_tot = χ_B + χ_GR + χ_SR, with explicit closed-form expressions given in main-text Eqs. (39), (41), and (45). The authors demonstrate limiting-case reductions to previous equatorial/circular results, present contour maps of the rotation angles, and apply the formalism to hot-spot polarimetric loops, disk-ring polarization, and images of a geometrically thick disk.","tokens_in":26429,"tokens_out":19861,"duration_ms":179389,"significance":"If correct, the formalism provides a very fast alternative to numerical ray tracing for computing Stokes images and time-dependent polarimetric signatures of Schwarzschild black hole accretion flows, filling a real gap: previous analytical models were limited to equatorial, circular motion. The vector-algebra framework in Appendix A is coherent, and the checks that the η=0, σ=90° limit reproduces Loktev et al. (2022) and that a vertical field yields constant PA are valuable sanity checks. The clear separation of GR, SR, and magnetic-field contributions is pedagogically useful. However, the significance is conditional on the correctness of the explicit main-text coordinate expressions, and I find that two of them (Eqs. 39 and 45) are inconsistent with the paper's own vector derivations in Appendix A.","major_comments":[{"comment":"Equation (39) does not follow from the general vector expression for χ_SR. Equation (A.8) contains in the numerator the factor k̂0·(û×n̂), which the authors themselves define as b̃ in Eq. (26). The main text instead uses cosα = k̂0·r̂. These two quantities are equal only for η = 0 and σ = 90°, because then û×n̂ = r̂. For the new regimes the paper claims to cover, e.g. η ≠ 0 or σ = 180°, b̃ and cosα differ; for σ = 180° at η = 0, b̃ = − sinα sin i sinφ / sinψ while cosα is positive. The correct specialization of Eq. (A.8) is tanχ_SR = −β b̃ cosζ / (sin²ζ − β cosξ). Because Eq. (39) is presented as the explicit expression for χ_SR and is used in the statements about the new effects, this is a load-bearing error that must be corrected and the affected calculations re-examined.","section":"§2.3, Eq. (39) and Appendix A, Eq. (A.8)"},{"comment":"Equation (45) appears to be the equatorial circular formula of Loktev et al. (2022) with cosα replaced by b̃, but that substitution is not sufficient for the general case and the expression is inconsistent with Eq. (A.6) even in a simple flat-space limit. Evaluating Eq. (A.6) for i = 30°, φ = 0°, η = 0, σ = 90°, β = 0.1, radial field (ϕ′ = 90°), gives tanχ_B ≈ +0.233, while Eq. (45) gives tanχ_B ≈ −0.173; the magnitudes disagree by about 30%. The direct geometric calculation of the angle between the bases B(k̂′₀, B̂′) and B(k̂′₀, n̂) confirms the vector-form result. The denominator and numerator of Eq. (45) are missing factors of δ and cosζ that appear in the derivation from Eq. (A.6). The correct coordinate expression must be re-derived from Eq. (A.6), and its reduction to the known η = 0, σ = 90° case must be verified.","section":"§2.3, Eq. (45)"},{"comment":"The paper lacks a numerical cross-check of the new terms (η ≠ 0, σ ≠ 90°) against exact geodesic integration. The only validations are limiting-case reductions and the reproduction of previous results for η = 0, σ = 90°. Given that the authors report PA differences of up to 30–90° for elevated matter and that the Poutanen (2020) bending formula was calibrated mainly for equatorial trajectories, an error budget for elevated and far-side photon paths is essential. The reader cannot assess whether the quoted accuracy (few degrees) survives for emission from above the plane or from the far side of the black hole. A comparison with a public numerical ray-tracing code for a set of representative (η, σ, i, r) points would resolve this concern.","section":"§3, Figs. 2, 5–7"},{"comment":"The same pattern of concern applies to Eq. (41) for χ_GR: the text says the expression is derived in Appendix A, but the reduction from the vector expression Eq. (A.10) is only demonstrated for η = 0. In light of the defects found in Eqs. (39) and (45), the derivation of Eq. (41) should be re-checked term-by-term against Eq. (A.10), and a numerical spot-check should be provided for η ≠ 0.","section":"§2.3, Eq. (41)"}],"minor_comments":[{"comment":"There is a stray text 'Loktev2024' on its own line that appears to be a leftover citation fragment; it should be removed or properly formatted.","section":"§2.2, after Eq. (35)"},{"comment":"The typesetting of Eq. (45) is ambiguous: the denominator contains an unmatched parenthesis and the fraction is hard to parse. Please rewrite with clear bracketing or display it in a more standard form.","section":"§2.3, Eq. (45)"},{"comment":"Table 1 lists σ = 90°, 120°, 150°, 180°, but the Fig. 7 caption says 'σ = 90° (black), 110° (red dashed), 150° (green dot-dashed), 180° (blue dotted)'. The value 110° vs 120° is inconsistent; please correct.","section":"Table 1 and Fig. 7 caption"},{"comment":"The sentence introducing the Luminet velocity law states that matter at height h moves with Keplerian velocity 'corresponding to its equatorial radius r_e'. This assumption is reasonable but should be justified, since it is an ad-hoc choice for elevated matter; the text later notes this is an assumption, but a brief comment on its physical motivation (e.g., vertical column rotating rigidly at the local Keplerian rate) would help.","section":"§2.1, Eq. (11)"}],"recommendation":"major_revision","confidential_remarks":"The main-text coordinate formulas (39) and (45) are internally inconsistent with the vector derivations in Appendix A, and the numerical spot-check I performed confirms that Eq. (45) gives incorrect values even in the flat-space equatorial case. This is more than a typo: it affects the central deliverable of the paper. The vector framework itself appears sound, so the paper is correctable, but it will require re-deriving the explicit formulas, re-running the calculations, and adding numerical validation. I recommend major revision with the requirement that the authors reconcile the main-text expressions with Appendix A and provide at least a spot-check against exact ray tracing for the new parameter regimes."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi — quick take on 2412.08359.\n\nThe paper is a useful extension of the artpol analytical ray-tracing formalism to synchrotron emission, adding non-zero disk thickness (η≠0) and non-circular velocity (σ≠90°). The vector-algebra setup in Appendix A is clean, the limiting-case checks back to Loktev et al. are sensible, and the applications to EHT/GRAVITY polarimetry are timely. For that, it deserves a real referee.\n\nBut there is a serious internal inconsistency in the central SR rotation formula. Eq. (39) claims tan χ_SR = −β cos α cos ζ / (sin^2 ζ − β cos ξ). The appendix derivation, Eq. (A.8), gives tan χ_SR = −β (k0·(u×n))(k0·n) / (1 − (k0·n)^2 − β k0·u). With your definitions, k0·(u×n) is exactly tilde{b} from Eq. (26), not cos α. These two agree only for the previously treated case η=0, σ=90°, because then u×n = r. For radial infall, or any elevated emitter, they differ. So the main-text equation does not follow from the paper’s own derivation. That is a load-bearing problem: either the figures were computed with Eq. (39) and the new results are wrong, or they were computed with Eq. (A.8) and the text needs a correction.\n\nThere is also the secondary issue the reader flagged: the Poutanen bending approximation is imported without a numerical error budget for the new elevated geometries. That is worth asking for, but the internal inconsistency is the sharper problem.\n\nThe paper is for people modeling polarized images and flares of Sgr A* and M87*. With the chi_SR issue fixed, the formalism would be a solid convenience tool. As it stands, I would not use Eq. (39) in my own work until the authors clarify which expression they actually used. I would still send it to peer review, with a request to reconcile the text with Appendix A and to add a numerical cross-check against exact geodesic integration for a few elevated, non-circular cases.\n\nNet: worth refereeing, but only after the authors resolve the Eq. (39)/(A.8) discrepancy.","headline":"Useful extension of the artpol ray-tracing formalism to synchrotron emission, but the main-text SR rotation formula conflicts with the appendix derivation and must be fixed before the new elevated/non-circular results can be trusted.","tokens_in":27033,"tokens_out":10520,"would_cite":false,"duration_ms":95165,"reading_group":"yes","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper derives closed-form formulas for the observed synchrotron polarization angle around a non-spinning black hole, as the sum of three independent rotations, valid for matter at any height and with any in-plane velocity direction.","keywords":["synchrotron polarization","analytical ray tracing","Schwarzschild black hole","polarization angle rotation","accretion disk thickness","Sgr A*","M87*","Stokes parameters"],"falsifier":"A direct numerical test would solve the Schwarzschild null geodesic and parallel-transport the electric-field vector for a point source at radius $r=3$ with elevation $\\eta=20^\\circ$, velocity direction $\\sigma=120^\\circ$, and a vertical magnetic field, then compare the two Stokes components over one full orbit with the closed-form expressions. A disagreement larger than the stated accuracy of the light-bending approximation would mark the radius and elevation range where the formulas stop being quantitative.","tokens_in":25884,"feed_emoji":"🌀","tokens_out":12499,"duration_ms":122480,"temperature":0.7,"pith_summary":"This paper aims to make polarimetric ray-tracing around a non-spinning black hole a matter of arithmetic rather than geodesic integration. It derives explicit formulas for the total rotation of the synchrotron polarization plane along the photon path, written as the sum of three separately computed rotations: one from the magnetic-field orientation, one from gravitational light bending, and one from the relativistic aberration caused by the matter's motion. The formulas cover matter elevated above or below the disk plane and matter with a non-circular velocity component, both left out of earlier analytical treatments. If the derivation is correct, the two observed linear-polarization components (Stokes Q and U) of optically thin synchrotron sources can be obtained from local emission by multiplying by a rotation matrix, which makes time-resolved flare modeling and image generation fast enough for systematic parameter searches.","feed_headline":"Polarization angle near black holes collapses to a three-term sum","feed_subtitle":"Closed-form formulas compute the full Stokes vector for elevated, non-circular flows near a Schwarzschild black hole.","key_machinery":"The working object is a chain of changes of reference axes for the photon's electric-field vector: from the emission frame tied to the magnetic field, to the disk normal, to the local fluid frame, and finally to the observer, with each step contributing $\\chi_B$, $\\chi_{\\rm SR}$, and $\\chi_{\\rm GR}$. The chain is closed by the paper's approximate light-bending relation (Eq. 18), which connects the emission angle $\\alpha$ to the observer's viewing angle $\\psi$ and supplies the lensing factor $L$ (a magnification factor from light bending), so no null geodesics are integrated. The three tangent identities (Eqs. 39, 41, 45) are the load-bearing formulas: they express aberration, light bending, and magnetic-field orientation in terms of the same set of angles ($\\psi$, $\\alpha$, $i$, $\\eta$, $\\sigma$, $\\varphi$), and the rotation matrix built from their sum converts local intensity, polarization degree, and polarization angle into observed Stokes parameters.","core_discovery":"The central claim is that around a Schwarzschild black hole the observed polarization angle of synchrotron emission is $\\chi = \\chi_0 + \\chi_B + \\chi_{\\rm GR} + \\chi_{\\rm SR}$, with $\\chi_0 = 90^\\circ$ set by the synchrotron mechanism, and with the three rotation terms given by closed trigonometric formulas (Eqs. 45, 41, and 39). The paper adapts the artpol technique (analytical ray tracing for spectro-polarimetry) to synchrotron radiation, parameterizing the matter by its elevation angle $\\eta$ above the disk plane, the direction $\\sigma$ of its in-plane velocity, the observer's inclination $i$, and the local magnetic-field components. Each term reduces to a function of these parameters, so the whole transformation from the local Stokes vector to the observer is a rotation matrix with no numerical geodesics. Because linear-polarization degree is unchanged by the transformation, only the orientation of the polarization plane needs to be traced. The formulas reduce to the known equatorial, circular-orbit results as special cases, and they quantify how disk thickness and accretion velocity direction change the rotation, with differences up to about 30 degrees for matter at radii $r \\lesssim 10$ and elevations of $\\pm 20^\\circ$.","pith_inferences":["Beyond the paper, the closed-form rotation angles make it feasible to fit observed flare loops and images with gradient-based optimization, since derivatives of the Stokes components with respect to radius, inclination, elevation, and velocity direction can be evaluated from the same formulas.","The paper states that the general expressions carry over to the Kerr metric, but the light-bending approximation it uses does not include the frame-dragging twist of photon trajectories; a practical spinning-black-hole version would need an additional spin-dependent term or a restriction to regions where the approximation error stays below observational precision.","Treating each layer of a thick disk as a Faraday screen along the line of sight, which the paper mentions as a possibility, would convert these intensity maps into wavelength-dependent polarization-rotation predictions that multi-frequency observations can test.","Because the net polarization of a thin ring depends strongly on inclination, crossing zero at face-on viewing and growing toward edge-on, the same formulas could track Lense-Thirring precession of an inner disk by mapping each precession phase to a point in the observed polarization plane."],"forward_implications":["Any optically thin synchrotron source around a non-spinning black hole can have its observed linear-polarization components computed by closed-form multiplication and rotation, so full images and light curves no longer require per-pixel numerical ray tracing.","Disk thickness is not a small correction: for radii within about 10 Schwarzschild radii, raising matter by 20 degrees changes the polarization-angle rotation by amounts comparable to the total rotation (up to about 30 degrees, with local excursions near 90 degrees), so razor-thin disk models mispredict the polarimetry of thick flows.","Changing the direction of matter motion shifts the zero-polarization critical point and can switch the polarization-angle behavior from one loop to two loops per orbit, which makes radial inflow degenerate with inclination in certain parameter ranges.","The zero-polarization critical point exists only for vertical magnetic fields and a limited range of observer inclinations (roughly 140 to 165 degrees for the Sgr A* geometry), so detecting a one-loop polarimetric flare simultaneously constrains the spot radius, viewing angle, and field topology.","Geometrically thick disk models produce ring-like images with a squashed shadow and spatially structured depolarization, giving analytic templates for horizon-scale interferometric images of M87* and Sgr A*."],"supporting_citations":[{"why":"Supplies the higher-order approximate light-bending relation and lensing factor that close the analytic system without numerical geodesics.","marker":"Poutanen (2020)"},{"why":"Establishes the basis-transformation decomposition of polarization rotation into the special-relativity and general-relativity terms that this paper extends to elevated and non-circular motion.","marker":"Loktev et al. (2022)"},{"why":"Validates the light-bending approximation against Kerr ray tracing and supplies the accuracy budget that the present paper inherits for its new geometries.","marker":"Loktev et al. (2024)"},{"why":"Provides the early simple light-bending and lensing-factor relations used in the flux and solid-angle expressions.","marker":"Beloborodov (2002)"},{"why":"Supplies the velocity and redshift law used to set the matter speed at each radius in the Schwarzschild metric.","marker":"Luminet (1979)"},{"why":"Gives the polarization degree for power-law synchrotron emission, which sets the local Stokes vector that the rotation matrix acts on.","marker":"Ginzburg & Syrovatskii (1965)"}],"fun_headline_variants":["Black hole polarization angle collapses to three-term sum","Closed-form sum explains black hole polarization angle","Analytical ray-tracing: polarization angle becomes three-term sum","Polarization angle around black holes: a three-term sum","Black hole synchrotron polarization: closed-form rotation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the paper's approximate formula connecting the photon's emission angle to the observer's viewing angle stays accurate to within a few degrees for trajectories that start above or below the disk plane and on the far side of the black hole; the paper quotes accuracy only for the equatorial, spinning-hole cases and does not supply a numerical error budget for the new geometry.","fun_headline_variants_meta":{"raw":{"variants":["Black hole polarization angle collapses to three-term sum","Closed-form sum explains black hole polarization angle","Analytical ray-tracing: polarization angle becomes three-term sum","Polarization angle around black holes: a three-term sum","Black hole synchrotron polarization: closed-form rotation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000643,"raw_usage":{"total_tokens":2974,"prompt_tokens":980,"completion_tokens":1994,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":596,"completion_tokens_details":{"reasoning_tokens":1916}},"tokens_in":596,"tokens_out":1994,"duration_ms":16129,"temperature":1.0,"reasoning_tokens":1916,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:53:30.618559+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct numerical test would solve the Schwarzschild null geodesic and parallel-transport the electric-field vector for a point source at radius $r=3$ with elevation $\\eta=20^\\circ$, velocity direction $\\sigma=120^\\circ$, and a vertical magnetic field, then compare the two Stokes components over one full orbit with the closed-form expressions. A disagreement larger than the stated accuracy of the light-bending approximation would mark the radius and elevation range where the formulas stop being quantitative.","supporting_citations":[{"cited_title":"2022, A&A, 660, A25","cited_arxiv_id":null,"evidence_quote":"Establishes the basis-transformation decomposition of polarization rotation into the special-relativity and general-relativity terms that this paper extends to elevated and non-circular motion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Validates the light-bending approximation against Kerr ray tracing and supplies the accuracy budget that the present paper inherits for its new geometries."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the early simple light-bending and lensing-factor relations used in the flux and solid-angle expressions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the polarization degree for power-law synchrotron emission, which sets the local Stokes vector that the rotation matrix acts on."}],"review_version":1}