{"id":"72eb26da-af14-478a-8121-ffc85593e24b","arxiv_id":"2505.16846","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Using ray-traced GR-MHD thin-disk models scaled to M87 and SgrA*, the paper shows that at 230 GHz high-mass systems are optically thin while low-mass systems are optically thick, making outflow emission more visible in lower-mass AGNs.","lead":"The paper simulates thin accretion disks around black holes with resistive GR-MHD and ray-traces the emission to predict how M87-like and SgrA*-like systems would look at horizon scales. It finds that at 230 GHz, high-mass systems are optically thin and dominated by the photon ring, while low-mass systems are optically thick and show stronger outflow emission, which may guide future VLBI observations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"High-mass optically thin at 230 GHz rests on an unconstrained SED peak; the paper's own Sect. 4 says observed M87/SgrA* SEDs peak above 230 GHz.","rationale":"I read the paper as an exploratory numerical experiment whose central claim is explicitly conditional on a toy model. The reader's weakest assumption—that a thin, resistive, axisymmetric Keplerian disk represents the accretion-ejection structure of real LLAGNs such as M87 and SgrA*—is reasonable and is acknowledged by the authors. However, there is a sharper, text-internal problem: the central high-mass optically-thin result depends on the model's SED peak position relative to 230 GHz, and that peak position is not constrained by the fitting procedure. In fact, Sect. 4 states that the literature SEDs for M87 and SgrA* peak above 230 GHz and 345 GHz, which would make both systems optically thick at these frequencies. The paper is transparent about this tension, which is why I do not move the verdict to reject: the work is a legitimate toy-model study, and the authors explicitly say they are not proposing alternative models for M87 and SgrA*. Still, the advertised conclusion that high-mass systems are optically thin and photon-ring-dominated at 230 GHz is not supported by the observational constraints as presented; it is a consequence of the unconstrained, model-dependent peak location. Conditional acceptance is therefore the right level, with the concrete test above as the condition that would settle whether the claimed dichotomy survives once the observed SED peak is enforced. The paper does have independent support in the form of the publicly available GRTRANS code and previously published simulation data, but the radiative postprocessing parameters (R_low = 1, R_high = 80) and the low-frequency-only fitting procedure are exactly where the load-bearing vulnerability sits.","tokens_in":19307,"tokens_out":7610,"duration_ms":63179,"concrete_test":"Refit the HIGH-mass SIM20 (and at least one other model) including the published M87 sub-mm/far-infrared SED points from Prieto et al. (2016), rather than only low-frequency radio points, or equivalently impose nu_peak > 230 GHz as a prior. If the best-fit peak stays above 230 GHz, compute tau_230 along the image-plane line of sight; tau > 1 would directly falsify the claim that M87-like high-mass systems are optically thin at 230 GHz and dominated by the lensed photon ring rather than outflow emission.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.2's central dichotomy—230 GHz (and 345 GHz) on the optically thin branch for the HIGH-mass system but on the self-absorbed branch for LOW-mass—requires the model's high-mass synchrotron peak to lie below 230 GHz. The SED fits in Sect. 3.1, however, deliberately constrain only the low-frequency data points ('an approximate fit to the data points at low frequencies is assumed to be correct'), so the position of the peak is not tested against the M87 data. Section 4 then concedes that the literature SEDs for M87 and SgrA* have peaks above 230 GHz and 345 GHz for both systems, which would put both sources on the self-absorbed slope at the frequencies of interest. The optically thin high-mass regime is therefore an output of the toy thin-disk model, not a consequence of the observational constraints the paper uses to label the systems. If the true or data-constrained peak for the HIGH-mass case lies above 230 GHz, the predicted LOW/HIGH contrast disappears and the 'photon-ring ideal' conclusion for high-mass systems is not supported. This is the most load-bearing point because everything downstream—outflow prominence, ring dominance, resolvability—follows from placing 230 GHz on opposite sides of the SED peak.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents six axisymmetric resistive GR-MHD simulations of thin Keplerian accretion disks, post-processed with the GRTRANS ray-tracing code. The authors scale the scale-free simulations to the masses and SEDs of M87 (HIGH-mass) and SgrA* (LOW-mass), adjusting the Eddington ratio to fit low-frequency radio data. They find that at 230 GHz and 345 GHz the HIGH-mass models lie on the optically thin synchrotron branch (yielding a lensed-photon-ring-dominated image), while the LOW-mass models lie on the self-absorbed branch (yielding relatively brighter outflow emission). They conclude that high-mass systems are ideal for probing photon rings and low-mass systems for probing outflows, and that current EHT resolution cannot distinguish the models.","tokens_in":19604,"tokens_out":7754,"duration_ms":58463,"significance":"If the central mass-dependent optical-thickness dichotomy were robustly established, it would be a useful organizing principle for interpreting EHT/ngEHT images of low-luminosity AGNs and for planning space-VLBI observations. The paper's strengths are a clearly described numerical pipeline using public codes (rHARM3D, GRTRANS), an explicit enumeration of six dynamical models with different floor densities, spins, and magnetic field strengths, and falsifiable predictions (e.g., photon-ring dominance for high-mass systems at 230 GHz, outflow brightness for low-mass systems). The models also reproduce Eddington-ratio magnitudes consistent with EHTC MAD/SANE fits. However, the central result depends on the position of the synthetic SED peak, which is not constrained by the fitted data and is in tension with the observed SEDs cited by the authors themselves.","major_comments":[{"comment":"The claim in Sect. 3.2 that 230 GHz and 345 GHz lie on the optically thin branch for the HIGH-mass system rests on the unconstrained location of the SED peak: Section 2.2 states that the fits target 'data points at low frequencies,' and Section 4 concedes that the observed M87 and SgrA* SEDs peak above 230/345 GHz, which would put both sources on the self-absorbed slope at these frequencies. The high-mass optically thin result is therefore an output of the toy model's assumed SED shape rather than a consequence of the observational constraints used to label the systems. Because the outflow-versus-photon-ring dichotomy follows entirely from this placement, this issue is load-bearing.","section":"Sections 2.2, 3.2, and 4"},{"comment":"The analysis uses models that the authors themselves rule out. The text states that SIM26 and SIM20 'can be completely ruled out' for the M87 SED and that SIM26 and SIM23 'can be completely ruled out' for the SgrA* SED, yet the subsequent image analysis and the general conclusions in Sections 3.3-3.5 and 5 draw on the full set of six models. Since the optically-thick/outflow-prominent behavior may be driven by the excluded models, the authors should either restrict the analysis to the viable subset (SIM21, SIM22, SIM24) or demonstrate that the conclusions are unchanged when the ruled-out models are removed.","section":"Section 3.1"},{"comment":"The generalization of the conclusions to real sources is not secure because the models are thin, axisymmetric, resistive Keplerian disks with density floors, whereas M87 and SgrA* are normally modeled as thick, optically thin ADAF-like MAD/SANE flows. The authors acknowledge this ('toy model'), but the conclusions are phrased for 'HIGH-mass systems' and 'LOW-mass systems' and specifically mention M87 as an example where photon-ring probing is favorable. A quantitative discussion of how a thick-disk geometry would shift the SED peak for the same mass and Eddington ratio is needed to support the applicability of the mass dichotomy to the EHT targets.","section":"Sections 2.1 and 4"},{"comment":"No uncertainties are provided for the Eddington ratios from the SED fits, and the fits fix Rlow=1, Rhigh=80, and i=17° without exploring the degeneracy of these parameters with the accretion rate. Because the position of 230 GHz relative to the peak is the key diagnostic, the absence of a parameter study or error estimate weakens the quantitative strength of the central claim. At minimum, the authors should identify how much the peak frequency and the 230-GHz optical depth vary over the plausible ranges of these postprocessing parameters.","section":"Section 3.1 and Table 1"}],"minor_comments":[{"comment":"There is a stray closing bracket after the citation: 'SED data of M87 and SgrA* (Narayan et al. 1998; Prieto et al. 2016)] as the constraint parameters.'","section":"Introduction"},{"comment":"In 'The thermal synchrotron spectra generated thus for HIGH-mass and LOW-mass systems is shown in Fig. 1,' the verb should agree with the plural subject 'spectra.'","section":"Section 3.1"},{"comment":"The sentence 'The reference data at different frequencies may have been obtained with a resolution and field of view (FOV) greater than or less than those assumed to obtain the total SED in this work, but for a broad comparison between the different models, an approximate fit to the data points at low frequencies is assumed to be correct' is ambiguous: it is not clear whether the comparison to data is normalized to the total SED or to the fitted frequency range, and this should be clarified.","section":"Section 2.2"},{"comment":"The statement that 'the peak frequency being the same for a given black hole mass irrespective of the simulated model implies that the peak frequency is determined by the emission from the shape of the accreting disk' is not fully explained; the relation to the disk shape rather than to the magnetic field strength should be clarified.","section":"Section 3.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is an exploratory modeling study with a clear, reproducible pipeline, but its central claim is more conditional than the abstract suggests. The authors should be asked to rework the analysis to separate model predictions from observational constraints and to be explicit about the toy-model dependence. The self-contradiction regarding ruled-out models is a correctable flaw."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the arXiv paper. The thing to know: this is a cleanly executed exploratory ray-tracing study of resistive GR-MHD thin-disk models, and the central mass-dependent optical-thickness claim is a model output, not something the SED fits to M87/SgrA* actually constrain. The paper is honest about that—Section 4 concedes the observed SEDs peak above 230/345 GHz for both sources, which would put both on the self-absorbed slope—but the abstract and conclusions state the LOW/HIGH contrast in more general terms than the models support.\n\nThe new bit is applying GRTRANS to these specific rHARM3D thin-disk simulations, separating disk from outflow with a polar-angle cut, and showing that in these models 230 GHz falls on opposite sides of the synchrotron peak for high- and low-mass scalings. That is a legitimate extension of the dynamical work from Bandyopadhyay et al. (2021), and it gives EHT/ngEHT target selectors a concrete hypothesis to check. The pipeline is described well: public ray-tracing code, six models, transparent parameter choices.\n\nThe load-bearing weakness is that the peak placement for the HIGH-mass case is not actually fit to the M87 data—the fits target low-frequency points only, and the peak position is essentially a consequence of the thin-disk model's density and temperature structure. The paper itself admits the real sources' SEDs peak above 230 GHz. So the 'photon-ring ideal' conclusion for high-mass systems is conditional on the toy model, not on the data. There are also no uncertainties on the Eddington ratios, the Rlow=1/Rhigh=80 and 10-degree decomposition are hand-picked, and the analysis deliberately retains models the SED fits rule out (SIM26 and SIM20 for M87; SIM26 and SIM23 for SgrA*). Those are real flaws, but the authors flag them; it is more overreach in framing than hidden error.\n\nWho this is for: someone working on EHT/ngEHT image interpretation and accretion-state contrasts. It won't settle whether M87 or SgrA* is thick or thin, but it lays out a crisp prediction that can be checked with multi-frequency VLBI.\n\nSend it to peer review. The methods are reproducible, the limitations are mostly acknowledged, and the mass-dependent SED-peak argument is worth airing even if the current data don't pin the peak. A good referee should push for the authors to either fit the peak frequency to the observed SEDs or soften the abstract.","headline":"Useful exploratory ray-tracing study with a mass-dependent optical-thickness claim that is more model-dependent than the abstract suggests; deserves peer review with pressure to fit or caveat the SED peak.","tokens_in":20190,"tokens_out":3018,"would_cite":false,"duration_ms":21919,"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":"The paper claims that a black hole's mass and accretion rate determine whether its horizon-scale emission shows outflow or the lensed photon ring, purely through where the observing frequency sits on the synchrotron spectrum.","keywords":["black hole accretion","thin accretion disks","GR-MHD simulations","ray tracing","synchrotron self-absorption","photon ring","AGN jets and outflows","Event Horizon Telescope"],"falsifier":"Take a low-luminosity AGN of roughly $10^{6}$ to $10^{7}$ solar masses with a measured SED peak above 230 GHz, and image it at 230 GHz with a space-VLBI baseline resolving about 0.16 microarcsecond: the paper predicts the emission should be optically thick and dominated by outflow/jet structure rather than the lensed photon ring, so observing a ring-dominated image would falsify the central claim.","tokens_in":19109,"feed_emoji":"🕳️","tokens_out":7516,"duration_ms":55379,"temperature":0.7,"pith_summary":"This paper argues that for accreting black holes, the observability of outflows versus the lensed photon ring at a given frequency is set by the position of that frequency on the thermal synchrotron spectrum, which in turn is fixed by the black hole mass and accretion rate. Using resistive GR-MHD simulations of thin Keplerian disks, ray-traced under the mass and spectral constraints of M87 and SgrA*, the authors find that at 230 GHz a low-mass system sits on the self-absorbed side of the spectrum and therefore appears optically thick, with disk-wind and jet emission dominating the image. A high-mass system with similar physics sits on the optically thin side, so the lensed photon ring dominates. The authors conclude that probing outflow emission on horizon scales favors low-mass AGNs, while high-mass AGNs are better targets for gravitational-lensing and photon-ring tests.","feed_headline":"Mass decides: outflow glow or photon ring at 230 GHz","feed_subtitle":"A mass-based rule tells which AGNs reveal outflows and which reveal the lensed ring.","key_machinery":"The carrying mechanism is synchrotron self-absorption relative to the SED peak: below the peak the source is optically thick and brightness traces density, while above it the source is optically thin and brightness traces the electron distribution. The peak frequency itself is shown to depend on black hole mass (shifting lower for higher mass) rather than on the details of the simulation, so the observer's frequency relative to the peak determines whether the disk, outflow, or photon ring dominates the image. This is implemented by postprocessing resistive GR-MHD snapshots of a thin Keplerian disk with the GRTRANS ray-tracing code, using an electron-temperature prescription with a plasma-beta-dependent ion-to-electron temperature ratio.","core_discovery":"The central claim is that, for a given dynamical model, the emission properties at an observing frequency are completely determined by where that frequency lies on the synchrotron SED, and the SED's peak frequency is set by the black hole mass and accretion rate. Because synchrotron self-absorption makes the spectrum optically thick below the peak, a low-mass system like SgrA* at 230 GHz is self-absorbed and its low-density outflow regions brighten, whereas a high-mass system like M87 at the same frequency is on the optically thin tail and exhibits mainly the lensed photon ring. The paper demonstrates this mass-dependent contrast through synthetic intensity maps and radial profiles, separates disk from outflow contributions, and notes that the Eddington ratios inferred from thin-disk SED fits match those derived by the EHT with thicker MAD/SANE models.","pith_inferences":["If the mass-position-on-SED logic holds beyond thin disks, it suggests a selection rule for future VLBI targets: low-mass LLAGNs are the natural laboratories for outflow-launching physics, and high-mass ones for strong-gravity tests, regardless of whether the underlying flow is thin or thick.","The same reasoning predicts that pushing to higher observing frequencies (e.g., 345 GHz) should make high-mass systems even more ring-dominated, while low-mass systems only become optically thin at still higher frequencies; this is testable with ngEHT-class arrays.","A direct observational discriminator would be a low-mass AGN with a known SED peak above 230 GHz: if sub-microarcsecond imaging still shows a ring-dominated image, the thin-disk assumption or the SED-position logic would need revision.","The axisymmetric thin-disk setup cannot capture turbulent, non-axisymmetric structures, so a natural next test is a 3D resistive GR-MHD run with slow-light ray tracing to see whether the mass-dependent outflow/ring contrast survives realistic variability."],"forward_implications":["At 230 GHz, low-mass AGNs (SgrA*-like) should show outflow-dominated, optically thick emission, while high-mass AGNs (M87-like) should show a dominant lensed photon ring.","The outflow component generally contributes the maximum to total emission at low inclination angles for most models, except the strongest-magnetic-field run.","Doppler beaming strongly modulates outflow brightness with viewing angle, especially in low-mass systems.","Current EHT resolution (about 20 microarcseconds) cannot distinguish the different thin-disk models; a space-VLBI baseline (Geo- or L2-class) is needed to separate them.","Because the inferred Eddington ratios match those from thicker MAD/SANE models, the thin-disk geometry requires higher density and optical depth to produce the same accretion power."],"supporting_citations":[{"why":"Supplies the six resistive GR-MHD thin-disk models and the earlier finding that black hole mass and accretion rate set the synchrotron SED peak.","marker":"Bandyopadhyay et al. (2021)"},{"why":"Provides the reference simulation setup and the disk/outflow separation by cutting 10 degrees above and below the equatorial plane.","marker":"Vourellis et al. (2019)"},{"why":"Provides the GRTRANS code used to ray-trace the synchrotron emission and solve radiative transfer along photon paths.","marker":"Dexter (2016)"},{"why":"Supplies the observed M87 SED used to fit the Eddington ratio for the HIGH-mass system.","marker":"Prieto et al. (2016)"},{"why":"Supplies the observed SgrA* SED used to fit the Eddington ratio for the LOW-mass system.","marker":"Narayan et al. (1998)"},{"why":"Establishes the resistive GR-MHD launching of disk winds that the thin-disk models rely on.","marker":"Qian et al. (2018)"},{"why":"Gives the plasma-beta-dependent electron-temperature prescription used to compute thermal synchrotron emission.","marker":"Mo´scibrodzka et al. (2016)"}],"fun_headline_variants":["Mass rule for AGNs at 230 GHz: outflow vs photon ring","Thin-disk AGNs: mass sets whether outflows or ring shine","Low-mass AGNs show outflows, high-mass show photon rings","Self-absorption makes AGN mass the key at 230 GHz","AGN mass flips emission: from self-absorbed outflows to ring"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central result assumes that a thin, axisymmetric, resistive Keplerian disk with numerical density floors is a faithful stand-in for the real accretion-ejection flow near the horizons of low-luminosity AGNs such as M87 and SgrA*, whose observed emission is usually modeled with thicker, optically thin flows.","fun_headline_variants_meta":{"raw":{"variants":["Mass rule for AGNs at 230 GHz: outflow vs photon ring","Thin-disk AGNs: mass sets whether outflows or ring shine","Low-mass AGNs show outflows, high-mass show photon rings","Self-absorption makes AGN mass the key at 230 GHz","AGN mass flips emission: from self-absorbed outflows to ring"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000496,"raw_usage":{"total_tokens":2495,"prompt_tokens":1073,"completion_tokens":1422,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":689,"completion_tokens_details":{"reasoning_tokens":1325}},"tokens_in":689,"tokens_out":1422,"duration_ms":8661,"temperature":1.0,"reasoning_tokens":1325,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:54:00.569662+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a low-luminosity AGN of roughly $10^{6}$ to $10^{7}$ solar masses with a measured SED peak above 230 GHz, and image it at 230 GHz with a space-VLBI baseline resolving about 0.16 microarcsecond: the paper predicts the emission should be optically thick and dominated by outflow/jet structure rather than the lensed photon ring, so observing a ring-dominated image would falsify the central claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the six resistive GR-MHD thin-disk models and the earlier finding that black hole mass and accretion rate set the synchrotron SED peak."},{"cited_title":"2016, MNRAS, 462, 115","cited_arxiv_id":null,"evidence_quote":"Provides the GRTRANS code used to ray-trace the synchrotron emission and solve radiative transfer along photon paths."},{"cited_title":"E., Popham, R","cited_arxiv_id":null,"evidence_quote":"Supplies the observed SgrA* SED used to fit the Eddington ratio for the LOW-mass system."},{"cited_title":"2018, ApJ, 859, 28","cited_arxiv_id":null,"evidence_quote":"Establishes the resistive GR-MHD launching of disk winds that the thin-disk models rely on."}],"review_version":1}