{"id":"28058ae7-f732-44de-83f0-c184530e3c34","arxiv_id":"2504.13304","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The Faraday-thick accretion flow around M87 should imprint a measurable, linear-in-wavelength-squared rotation and depolarization on the counter-jet's 43 GHz polarization, enabling accretion-flow constraints.","lead":"This paper predicts that the hot accretion flow around M87 should Faraday-rotate and depolarize the counter-jet's radio polarization at 43 GHz, even though the flow is optically thick to Faraday rotation. A smart generalist might read it because it offers a way to detect and measure the accretion flow onto a supermassive black hole, and it warns that previous rotation-measure-based mass estimates may need rethinking.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The RM fitted from the EVPA-λ² slope does not match the Faraday depth in the koral3D model, so the claim that a measured RM yields lower limits on density and accretion rate is not established.","rationale":"The reader's weakest assumption—compact blob size and density normalization—is a genuine threat to M87-specific applicability, and the paper's own §4.3 admits that a larger emission region can disrupt the linear λ² pattern. I agree that this warrants a conditional accept with softened claims and wider parameter exploration. However, I identify a more load-bearing concern that the reader did not emphasize: even granting the compact blob and the chosen normalization, the paper's most realistic koral3D model produces a fitted RM that is decoupled from the actual Faraday depth by a factor of 40–50. The abstract's promise that the linear relationship persists 'enabling us to constrain the flow's physical properties' rests on the standard thin-screen mapping RM ∝ ∫ n_e B∥ dl, but the paper's own §4.1 shows that mapping fails in the turbulent, Faraday-thick regime it claims is relevant to M87. The reader's proposed fixes—larger blobs and multiple snapshots—would not resolve this discrepancy, because it already occurs for a 2 rg blob within a single snapshot. The paper should either demonstrate a monotonic, calibratable mapping between the fitted RM and the model's Faraday depth (e.g., by density rescaling) or restrict its final claim to qualitative detection via depolarization rather than quantitative lower limits. Since the qualitative finding that linear EVPA can persist in a Faraday-thick screen remains supported, the appropriate verdict stays CONDITIONAL.","tokens_in":13606,"tokens_out":16294,"duration_ms":164501,"concrete_test":"Using the same koral3D snapshot and blob setup, rescale only the accretion-flow electron density by factors of 10 and 100 (keeping B, T, and blob fixed), recompute the EVPA-λ² relation, and compare the fitted RM-derived τ to the τρV map. In an external screen obeying Eq. (2), both quantities should scale as n_e^{3/2} and the fitted τ should equal the map value. If the fitted τ remains anomalously low or does not track the map, the lower-limit inversion in §4.2 fails in exactly the regime claimed to apply to M87. A complementary check is to compute the intensity-weighted mean and RMS of τρV over the blob for all four φblob positions and test whether 2RMλ² equals either; if it matches neither, the observed RM cannot be used with Eq. (2).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative promise—that the linear EVPA-λ² relation measured through a Faraday-thick, turbulent screen can be used to place lower limits on electron density, magnetic field, and mass accretion rate—is undercut by its own numerical results. In §4.1 the authors state that 'the τρV values from the RM fit do not match the code calculations shown in the maps.' Concretely, for the koral3D blob at φblob=π/2 (Fig. 4b), the fitted RM=1.69×10^5 rad m−2 implies τ=2RMλ²≈17 at 43 GHz, while the Faraday-depth map near the blob (Fig. 5b) shows τρV≈640–780, a factor of 40–50 discrepancy. This is not a small calibration offset. The cause is Faraday-depth fluctuations across even the compact 2 rg blob: the observed polarization is the vector sum of rays with many different rotation angles, so the slope of the net EVPA is not the line-of-sight-integrated Faraday depth. Therefore the inference in §4.2 that 'a RM measurement implies a lower limit on τρV, density, and in turn, a lower limit on Mdot' is not justified by the presented models. This is an internal inconsistency between the reported RM fits and the code's Faraday-depth output, not a disagreement with external consensus.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that resolved linear polarization observations of the M87 counter-jet at 43 GHz can probe the accretion flow via Faraday rotation. Using an analytic cylinder model and two numerical setups (a semi-analytic RIAF and a koral3D GRMHD snapshot) with the grtrans ray-tracing code, the authors find in all cases a Faraday-thick accretion flow (τ_ρV ≫ 1) with RM ~ 10^6 rad m^-2, yet the EVPA retains a linear wavelength-squared dependence. The koral3D model also yields strong depolarization due to Faraday-depth fluctuations. The authors propose that comparing counter-jet and forward-jet polarization can detect the accretion flow and place lower limits on electron density, magnetic field, and mass accretion rate.","tokens_in":13850,"tokens_out":4384,"duration_ms":40472,"significance":"The central result—that a linear EVPA–λ² relation can persist even when the Faraday screen is thick and depolarizing—is an important caution against the widespread assumption that linearity implies Faraday thinness. The paper is genuinely forward-modeling in spirit: the predicted RM is not fitted to M87's observed RM, the radiative transfer is done with an independent public code, and the density-scaling and blob-density tests cleanly isolate the accretion flow as the screen. If the claims hold, the work gives a concrete, falsifiable prediction for upcoming polarimetric VLBI observations. However, the quantitative inference from a measured RM to lower limits on density and Ṁ is weakened by an internal inconsistency between the fitted RM and the code's Faraday depths in the turbulent model, as detailed below.","major_comments":[{"comment":"For the koral3D blob at φ=π/2, the fitted RM of 1.69×10^5 rad m^-2 implies τ_ρV ≈ 17 at 43 GHz, whereas the Faraday-depth map near the blob shows τ_ρV ≈ 640–780, a factor of 40–50 discrepancy. The paper states in §4.1 that 'the τρV values from the RM fit do not match the code calculations shown in the maps,' but it then proceeds in §4.2 to claim that 'a RM measurement implies a lower limit on τρV, density, and in turn, a lower limit on Mdot.' This is not justified: in a turbulent screen with large Faraday-depth fluctuations, the slope of the net EVPA is not the path-integrated Faraday depth. The authors should either provide a quantitative relation between the fitted RM and the screen properties, or soften the lower-limit claim to be only qualitative.","section":"§4.1, Figs. 4b and 5b"},{"comment":"The numerical setup states 'The simulation is conducted for a BH mass of 6.5×10^5 M⊙.' This contradicts the M87 mass of (6.5±0.7)×10^9 M⊙ used elsewhere in the paper. If the code genuinely used 10^5 M⊙, all physical scales (lengths, densities, magnetic fields from the model parametrization) would be wrong by orders of magnitude, and the predicted RM values would not apply to M87. Please correct the typo or, if the simulations actually used a different mass, restate the value explicitly and check the resulting RM normalization.","section":"§3.1"},{"comment":"The text says the Faraday-thin limit is imposed 'by decreasing the accretion flow density by ∼ 100,' then reports that polarization is recovered 'while τρV ≫ 1.' The figure caption instead says the purple points correspond to 'imposing στρV = 1.' These are different statements: reducing n by 100 should reduce τ_ρV by roughly a factor of 100, while setting στρV = 1 acts on the fluctuations. Please clarify what was actually done, report the post-decrease τρV and στρV values, and reconcile the text and caption. This matters for the paper's conclusion that στρV rather than τρV controls depolarization.","section":"§3.3.2, Fig. 6"}],"minor_comments":[{"comment":"The running header contains 'F araday Rotation'; the spacing appears to be a typesetting error.","section":"Title page"},{"comment":"The mass '6 .5× 105M⊙' has irregular spacing and should be formatted as '6.5×10^5 M⊙' or, more likely, '6.5×10^9 M⊙' per the previous comment.","section":"§3.1"},{"comment":"The phrase 'when imposing στρV = 1, marked by the vertical dotted line' is helpful, but the main text should use the same notation and explain how στρV = 1 is achieved in practice.","section":"Fig. 6 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper's core idea is interesting and the forward-modeling approach is a strength, but the RM-versus-Faraday-depth discrepancy in the koral3D model is the key technical issue. I would not reject on this basis, because the authors openly acknowledge the mismatch and the qualitative finding (linearity plus depolarization in a thick screen) is novel and testable. However, the quantitative claim about lower limits on density and Ṁ needs to be either re-derived or removed. The BH mass typo must also be corrected."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it, and I mostly side with the reader's conditional verdict, with one emphasis added: the stress-test note lands, and it lands on the load-bearing part of the quantitative claims.\n\nWhat is new and good. Using the M87 counter-jet at 43 GHz as a background polarized source viewed through the accretion flow is genuinely new—prior RM work used the forward jet or the unresolved core. The analytic disk model is clean and gives a concrete, falsifiable prediction (RM ~ 10^6 rad/m^2, tau_rhoV ~ 10^3). The numerical tests are well executed: dropping the RIAF density by 100 changes the RM by 1000, changing the blob density does nothing, and the Faraday-off runs show no rotation. That confirms the screen interpretation about as cleanly as a forward model can. There is no circularity—no target RM is fitted, and grtrans is public. The authors also flag their own biggest assumptions in Section 4.3 (blob size and morphology, single snapshot, density normalization), which is more honest than most. And the central interpretive point—that a linear EVPA-lambda^2 relation does not imply a Faraday-thin screen—is correct, important, and underappreciated in the LLAGN RM literature. Burn 1966 made the argument; this paper gives it a concrete M87 application.\n\nThe soft spots, in proportion. First, 'in all cases' overstates it. Every case shares the same density normalization, and their own test shows that cutting n by ~100 kills the thick regime and restores high polarization. The RM also varies by ~50x with blob azimuth, so 'RM ~ 10^6' is not a single prediction. That is a wording problem more than a substance problem, because the physical conclusion survives.\n\nSecond—and this is the real issue—Section 4.1 admits the tau_rhoV values from the RM fit 'do not match the code calculations shown in the maps.' For the pi/2 blob, the fitted slope implies tau ~ 17 at 43 GHz while the Faraday-depth map shows ~640-780. Yet Section 4.2 still asserts that a measured RM implies a lower limit on tau_rhoV, density, and in turn Mdot. In a turbulent thick screen, the slope of the net EVPA is the slope of a vector sum, not the line-of-sight Faraday depth; the paper's own numbers show the two can differ by a factor of 40-50. The lower limit might survive in direction (the fitted tau is below the map tau), but nothing in the paper shows the bias is one-sided or bounded. As written, the quantitative inference chain from RM to density to Mdot is not established by the presented models.\n\nThird, the linearity result depends on a compact 2-rg blob. Their own 10-rg koral3D test breaks the pattern. The authors disclose this, but it means the robust observable is depolarization; the linear EVPA is conditional on counter-jet structure.\n\nBottom line. This is a solid forward-modeling paper with a real interpretive payoff. The qualitative prediction is testable with existing and future VLBI polarization observations of M87 or NGC 1052. It deserves a serious referee. My recommendation: send it to review, and require that the authors either calibrate the RM-to-tau relation against their own maps or soften the lower-limit claims in Section 4.2 and the abstract.","headline":"Solid new qualitative prediction—Faraday-thick M87 screen can still yield linear EVPA—but the paper's own koral3D fits miss the code's Faraday depth by up to ~50x, so the quantitative RM-to-density/Mdot lower-limit pipeline is not established.","tokens_in":14449,"tokens_out":6295,"would_cite":true,"duration_ms":56528,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"M87's accretion flow is Faraday-thick at 43 GHz yet still shows a linear polarization-angle signal.","keywords":["accretion flows","Faraday rotation","M87","counter-jet","radiatively inefficient accretion flow","polarization","supermassive black holes","mass accretion rate"],"falsifier":"Resolved 43 GHz linear polarization observations of M87's counter-jet over multiple epochs would settle it: the models predict $|\\mathrm{RM}|$ in the range roughly $10^5$–$10^7$ rad m$^{-2}$ with $\\tau_{\\rho V}\\gg1$ and low linear polarization degree, so measuring $|\\mathrm{RM}|$ well below $10^5$ rad m$^{-2}$ together with linear polarization above about 60 percent would rule out the assumed Faraday-thick screen.","tokens_in":13329,"feed_emoji":"🌀","tokens_out":11774,"duration_ms":95544,"temperature":0.7,"pith_summary":"The paper argues that the accretion flow onto M87's supermassive black hole is an external Faraday screen that is Faraday-thick at 43 GHz — the electric vector position angle (EVPA) rotates many times — yet it still produces the linear EVPA-versus-wavelength-squared relation usually taken as the hallmark of a Faraday-thin screen. Using an analytic cylindrical flow model and polarized ray-tracing radiative transfer through two accretion flow models, the paper finds rotation measures around $10^5$–$10^7$ rad m$^{-2}$ and Faraday depths far above unity in every case. The central consequence is that linearity alone cannot diagnose Faraday thinness in M87, and that resolved polarization measurements of the counter-jet, visible at 43 GHz, should reveal the accretion flow and place lower limits on electron density, magnetic field strength, and mass accretion rate. The same reasoning would apply to other low-luminosity active galaxies with resolved counter-jets, so a wrong assumption of Faraday thinness could bias RM-based accretion rate estimates elsewhere.","feed_headline":"Faraday-thick M87 flow still leaves a linear polarization signal","feed_subtitle":"If true, 43 GHz maps of the counter-jet could put lower limits on density, field, and accretion rate.","key_machinery":"The load-bearing element is the external Faraday screen: polarized synchrotron emission from a compact counter-jet blob at about 25 Schwarzschild radii passes through a geometrically thick, magnetized accretion flow whose Faraday depth $\\tau_{\\rho V}=2\\,\\mathrm{RM}\\,\\lambda^2$ is the line-of-sight integral of $n_e B_\\parallel$. The paper combines an analytic cylindrical disk model with polarized radiative transfer through a semi-analytic RIAF and a turbulent MHD snapshot. The result that carries the argument is that a large mean Faraday depth does not randomize the position angle; the EVPA keeps rotating with $\\lambda^2$, while the degree of depolarization is set by the fluctuation amplitude $\\sigma_{\\tau_{\\rho V}}$ rather than by $\\tau_{\\rho V}$ itself, as shown by recovery of high polarization when the flow is made Faraday thin.","core_discovery":"The paper's claim is that M87's accretion flow is Faraday thick at 7 mm but leaves a clean linear EVPA–$\\lambda^2$ signature. In the analytic cylindrical model, the screen gives $|\\mathrm{RM}| \\approx 9.93\\times 10^6\\,\\mathrm{rad\\,m^{-2}}$, i.e. Faraday depth $\\tau_{\\rho V}\\approx 973$ at 43 GHz. In a semi-analytic radiatively inefficient accretion flow with a compact polarized blob as the counter-jet, the fit gives $|\\mathrm{RM}| \\approx 1.13\\times 10^6\\,\\mathrm{rad\\,m^{-2}}$, $\\tau_{\\rho V}\\approx 111$. In a turbulent GRMHD snapshot, the screen remains Faraday thick with $\\tau_{\\rho V}\\gg 1$ for all tested blob positions, the EVPA still tracks $\\lambda^2$ linearly, and the emission is depolarized to a few percent by Faraday depth fluctuations. The paper reads this as showing that comparing the counter-jet and forward-jet polarization states can detect the accretion flow, and that linearity of EVPA with $\\lambda^2$ is not by itself evidence of a Faraday-thin screen.","pith_inferences":["A testable extension would be to run the same radiative-transfer experiment with a time sequence of turbulent MHD snapshots; the paper's single-snapshot result leaves open whether the linear $\\lambda^2$ relation persists across many realizations.","If the compact-blob assumption is relaxed, the linear pattern should break down for the turbulent screen, so multi-epoch, multi-frequency imaging could use the onset of nonlinearity to map the size of the counter-jet emission region.","The same Faraday-thick-but-linear behavior could occur along other lines of sight through radiatively inefficient flows, implying that unresolved RM measurements in low-luminosity active galactic nuclei may be more sensitive to screen fluctuations than to the mean magnetic field."],"forward_implications":["A linear EVPA versus $\\lambda^2$ relation measured in M87's counter-jet cannot by itself be taken as evidence for a Faraday-thin screen.","Comparing the counter-jet and forward-jet linear polarization at 43 GHz should reveal the accretion flow through rotation or depolarization and yield lower limits on electron density, magnetic field strength, and mass accretion rate.","Existing RM-based mass accretion rate limits for M87 that assumed a Faraday-thin spherical flow may instead be tracing the $\\tau_{\\rho V}=1$ surface, so the inferred rates would need reinterpretation.","The Faraday screen is time-variable in the turbulent model, so the measured RM and polarization degree should change with observing epoch and with position along the counter-jet.","Resolved polarized observations of other low-luminosity active galaxies with visible counter-jets, such as NGC 1052, could apply the same method."],"supporting_citations":[{"why":"Provides the theory that Faraday depth fluctuations depolarize the screen while leaving the mean EVPA–lambda^2 relation intact.","marker":"Burn 1966"},{"why":"Establishes the RM-to-mass-accretion-rate conversion for Sgr A* under a Faraday-thin spherical flow, the framework the paper contrasts with its Faraday-thick M87 prediction.","marker":"Marrone et al. (2006, 2007)"},{"why":"Gives the M87 230 GHz RM upper limit and accretion rate limit at 21 Schwarzschild radii that the paper reinterprets as possibly Faraday-thick.","marker":"Kuo et al. (2014)"},{"why":"Shows the counter-jet is detectable and resolved at 43 GHz, making it the observational target.","marker":"Walker et al. (2018)"},{"why":"Resolves counter-jet structure and linear polarization at 43 GHz, providing the comparison basis for forward-jet/counter-jet polarization.","marker":"Park et al. (2021)"},{"why":"Supplies the polarized radiative transfer machinery and the relation tau = 2 RM lambda^2 used to translate RM into Faraday depth.","marker":"Dexter 2016"},{"why":"Defines the semi-analytic RIAF density, temperature, and magnetic field profiles used as one accretion flow model.","marker":"Broderick et al. (2009)"},{"why":"Provides the turbulent GRMHD snapshot used as the more realistic accretion flow model.","marker":"Chael et al. (2019)"},{"why":"Supports neglecting the forward-jet contribution to RM because relativistic electrons suppress Faraday rotation.","marker":"Quataert & Gruzinov (2000)"}],"fun_headline_variants":["M87's Faraday-thick flow still shows linear polarization signature","Even when Faraday-thick, M87 flow leaves a clean EVPA","Counter-jet polarization can pin down M87 accretion flow","Linear Faraday signal despite thick flow in M87","M87 accretion flow: thick screen, but linear EVPA persists"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predictions assume the counter-jet is a small, compact patch of emission at a specific distance from the black hole, and that the surrounding gas has a particular density, temperature, and magnetic field strength; if the true emission region is larger or the gas thinner, the predicted Faraday-thick screen and linear signal weaken or disappear.","fun_headline_variants_meta":{"raw":{"variants":["M87's Faraday-thick flow still shows linear polarization signature","Even when Faraday-thick, M87 flow leaves a clean EVPA","Counter-jet polarization can pin down M87 accretion flow","Linear Faraday signal despite thick flow in M87","M87 accretion flow: thick screen, but linear EVPA persists"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000478,"raw_usage":{"total_tokens":2441,"prompt_tokens":1089,"completion_tokens":1352,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":705,"completion_tokens_details":{"reasoning_tokens":1267}},"tokens_in":705,"tokens_out":1352,"duration_ms":9637,"temperature":1.0,"reasoning_tokens":1267,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T12:10:46.143468+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Resolved 43 GHz linear polarization observations of M87's counter-jet over multiple epochs would settle it: the models predict $|\\mathrm{RM}|$ in the range roughly $10^5$–$10^7$ rad m$^{-2}$ with $\\tau_{\\rho V}\\gg1$ and low linear polarization degree, so measuring $|\\mathrm{RM}|$ well below $10^5$ rad m$^{-2}$ together with linear polarization above about 60 percent would rule out the assumed Faraday-thick screen.","supporting_citations":[{"cited_title":"2021, ApJ, 922, 180, doi: 10.3847/1538-4357/ac26bf","cited_arxiv_id":null,"evidence_quote":"Resolves counter-jet structure and linear polarization at 43 GHz, providing the comparison basis for forward-jet/counter-jet polarization."}],"review_version":1}