{"id":"51e93e76-f074-4f1e-8088-a274f0bdbb51","arxiv_id":"2505.13242","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"In the changing-state AGN NGC 1566, the 230 GHz millimeter flux tracks the 14-150 keV X-ray flux during a state transition, indicating the millimeter emission likely originates in the X-ray corona.","lead":"Astronomers used archival ALMA and Swift observations to watch the compact millimeter emission of the changing-state active galaxy NGC 1566 brighten and fade along with its X-ray emission. The result suggests the millimeter light comes from the hot magnetized corona around the black hole, rather than from a jet or star formation, supporting the use of millimeter emission as a tracer of black hole activity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Figure 3 gives p=0.168 for the mm–X-ray Spearman correlation, but the text quotes p=0.0168; recomputation from Table A.2 confirms rho≈0.54, p≈0.17, so the 8-epoch correlation is not significant and cannot support the abstract's 'positive correlation' or the corona-origin conclusion.","rationale":"The reader's weakest assumption identifies the non-simultaneity and the p-value inconsistency; the p-value issue is the most load-bearing because it is an internal contradiction that can be settled immediately from the paper's own table. Recomputing the Spearman statistic from Table A.2 confirms rho ≈ 0.54 and p ≈ 0.17, matching Figure 3 and not the text's 0.0168. Thus the central empirical claim — a positive mm–X-ray correlation — is not significant at the 5% level with eight epochs. The corona-origin conclusion in Section 3.3.4 explicitly relies on 'the observed tight correlation between the mm and X-ray flux,' so it collapses with the correlation. The paper still contains useful ALMA flux measurements and a plausible discussion of source size and magnetic field, but those do not establish a common origin for the mm and X-ray emission. A revised paper could report the non-detection of a significant correlation, soften the abstract and conclusions, and reframe the work as a variability study; as written, the headline claim overstates the evidence and should be rejected.","tokens_in":14985,"tokens_out":5625,"duration_ms":50005,"concrete_test":"Recompute the Spearman rank correlation and an exact permutation p-value from the eight (log F230GHz, log F14-150keV) pairs in Table A.2, preferably from the unrounded values before the quoted 0.13 dex uncertainties are applied. If the two-sided p is approximately 0.17, the text's p = 0.0168 is a typo and the abstract and conclusions must be revised to state that no significant correlation is found. As a second check, rerun the correlation using only the pre-outburst and post-outburst epochs and also omitting the two 2018 outburst epochs (epochs 4 and 5), to test whether any apparent trend is driven solely by the state change rather than by a steady mm–X-ray coupling.","verdict_should_be":"REJECT","load_bearing_attack":"The central conclusion, that the mm emission originates in the magnetized X-ray corona because it tracks the X-ray flux, rests on the claimed mm–X-ray correlation in Section 3.2. That correlation is not statistically established. The manuscript contains two inconsistent p-values: the text reports p = 0.0168, while Figure 3 reports p-value: 1.68e-01. Recomputing the Spearman correlation from the eight log F230GHz and log F14-150keV pairs in Table A.2 gives rho ≈ 0.54 and a two-sided p ≈ 0.17, matching the figure rather than the text. With n = 8, rho = 0.54 is comfortably non-significant, so the null hypothesis of no correlation cannot be rejected. The abstract's 'positive correlation' and the conclusion's 'strong correlation' are therefore unsupported by the paper's own data. The quoted intrinsic scatter of 0.05 dex, used as evidence of tightness, is also hard to reconcile with the observed scatter and with the non-significant Spearman statistic; it depends on the fitting method and does not rescue significance. Because the corona-origin argument in Section 3.3.4 is built on this correlation, the headline claim inherits the statistical failure. Non-simultaneous X-ray sampling and the Gamma = 1.8 flux conversion compound the problem, but even taking those at face value, the claimed correlation does not reach significance.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports the first ALMA Band 6 (230 GHz) continuum variability study of a changing-state AGN, NGC 1566, using eight archival epochs from 2014–2023. The authors measure the nuclear mm flux and in-band spectral index, compare the mm flux to 14–150 keV X-ray fluxes, and argue that the mm emission tracks the X-ray emission during the 2018 changing-state event. They interpret the flat spectral index, the mm emitting size, and the exclusion of star formation, compact jets, and outflow shocks as evidence that the mm emission arises from synchrotron radiation in the magnetized X-ray corona. The central quantitative claim is a positive mm–X-ray correlation with a Spearman rho of 0.54, p = 0.0168, and an intrinsic scatter of 0.05 dex.","tokens_in":15244,"tokens_out":7432,"duration_ms":64273,"significance":"If the claimed correlation and coronal origin were robustly established, this would be a valuable single-object contribution to the debate on the origin of compact mm emission in radio-quiet AGNs, and it would support the use of mm emission as a tracer of AGN activity. The paper has genuine strengths: it analyzes archival ALMA data with attention to beam-size effects (Appendix C), provides variability information across an accretion-state transition, and discusses alternative emission mechanisms in a structured way. However, the headline correlation is not statistically significant on the data presented, and the reported intrinsic scatter is difficult to reconcile with the measurement uncertainties. The scientific value of the study is real, but the claims as currently written outrun the statistical evidence.","major_comments":[{"comment":"The text reports a Spearman p-value of 0.0168 for the mm–X-ray correlation, while Figure 3 shows p = 1.68e-01. Recomputing the Spearman correlation from the eight (log F14-150 keV, log F230 GHz) pairs in Table A.2 gives rho ≈ 0.54 and a two-sided p ≈ 0.17, consistent with the figure and not with the text. With n = 8, rho = 0.54 does not reach the 0.05 significance level, so the null hypothesis of no correlation cannot be rejected. Therefore the abstract's 'positive correlation', Section 3.3.4's 'tight correlation', and Summary item 1's 'strong correlation' are not supported by the paper's own data. Because the corona-origin conclusion in Section 3.3.4 is built on this correlation, the main claim needs substantial reanalysis or rephrasing.","section":"§3.2 and Figure 3"},{"comment":"The reported intrinsic scatter of 0.05 dex is inconsistent with the quoted flux uncertainties and with the data. The 230 GHz fluxes carry 0.13 dex uncertainties (10% errors per Table A.1 note), so the observed scatter around the best-fit line should be at least comparable to the measurement error unless the model absorbs it; an intrinsic scatter smaller than the errors is not plausible without a detailed fitting procedure. Computing residuals from the stated best-fit relation log F230 = 0.39 log F14-150 - 10.62 gives an RMS residual of about 0.10 dex, not 0.05 dex. The authors should specify the regression method, how the intrinsic scatter was derived, and whether measurement errors were included; as presented, the tightness claim is unsupported.","section":"§3.2, abstract, and Table A.2"},{"comment":"The X-ray fluxes used in the correlation are not simultaneous with the ALMA observations and are converted to the 14–150 keV band from different original bands (14–195 keV, 15–50 keV, 2–10 keV) assuming a fixed photon index Gamma = 1.8. During a rapidly changing state, a fixed spectral slope and month-long binning of BAT data can introduce systematic scatter and possibly bias the correlation. The authors should either repeat the correlation using only strictly contemporaneous measurements, or explicitly quantify how the non-simultaneity and the Gamma = 1.8 assumption affect the result. This is a secondary concern given the non-significant Spearman test, but it should be addressed before any revised correlation claim is made.","section":"§2 and §3.2"}],"minor_comments":[{"comment":"The bottom axis of Figure 3 appears to show MJD values (57000–60000) rather than log F14-150 keV, despite the caption describing the plot as a relation between mm and X-ray flux. The axis labels should be corrected so that the figure matches its caption.","section":"Figure 3"},{"comment":"The text refers to 'free-free emission from dust' and 'free-free emission from a dusty torus'; free-free emission originates from ionized gas, not from dust. The wording should be changed to, e.g., 'free-free emission from ionized gas in the dusty torus' or similar.","section":"§3.3.1 and §3.4"},{"comment":"In the outflow rate calculation, the substitution 'NH = nR' should read 'NH = ne R' or 'NH = ne R' with ne the electron density; the symbol n is ambiguous.","section":"§3.3.3"},{"comment":"The observatory name is misspelled as 'ASSAS-Sn'; it should be 'ASAS-SN'.","section":"Appendix B"},{"comment":"The entries in the alpha_mm column contain spaces between the integer and the decimal point (e.g., '0 .11± 0.26'). These should be formatted as decimal numbers.","section":"Table A.2"}],"recommendation":"major_revision","confidential_remarks":"The discrepancy between the p-value quoted in the text (0.0168) and the one shown in Figure 3 (0.168) is serious enough that the editor should ask the authors to verify all reported statistics. The regression-derived intrinsic scatter also needs scrutiny. If the correlation is genuinely not significant, the manuscript can still be publishable as a single-object variability case study, but the abstract and conclusions must be rewritten to report a tentative trend rather than a 'strong correlation'. The paper's acknowledgment of a specific reviewer by name is unusual but not scientifically problematic."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the first ALMA mm variability study of a changing-state AGN is genuinely new, but the central correlation between mm and X-ray flux does not survive contact with the paper's own statistics. Text says p=0.0168; Figure 3 says p=1.68e-01. Recomputing from Table A.2 gives rho≈0.54, p≈0.17, so with eight epochs the correlation is not significant. The abstract's 'positive correlation' and conclusion's 'strong correlation' go beyond the data.\n\nWhat earns credit: the ALMA reduction looks careful, including the beam-size consistency check in Appendix C. The eight-epoch light curve shows a clean factor-of-two rise around the 2018 outburst and decay afterward. The flat in-band spectral index is consistent with a compact optically thick source. The authors do useful work rejecting star formation, dust, free-free, jet, and outflow contributions; those arguments are not circular and the numbers are plausible.\n\nThe soft spots are real and load-bearing. The p-value mismatch is not a typo; it is the evidentiary foundation for the corona-origin claim. The 0.05 dex intrinsic scatter also sits oddly with the observed scatter around the fit, so the fitting method needs reporting and justification. The X-ray fluxes are binned over 30 days, converted from different bands assuming Gamma=1.8, and not simultaneous with ALMA; that compounds the significance problem. The linear fit slope of 0.39±0.13 may be parametrically significant, but the non-parametric Spearman test is not, and the paper never addresses that tension.\n\nWhom is it for? AGN observers working on mm emission and changing-look events. The light curve and spectral indices are useful reference data even if the interpretation is too strong. It deserves serious peer review, not desk rejection, but a good referee should demand corrected statistics, explicit fitting details, and a more restrained conclusion.","headline":"The new ALMA data on NGC 1566 are worth a look, but the paper's central mm-X-ray correlation does not survive contact with its own statistics, so the corona-origin conclusion is overstated.","tokens_in":637,"tokens_out":1260,"would_cite":false,"duration_ms":60527,"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 230 GHz glow of NGC 1566 comes from the magnetized X-ray corona.","keywords":["NGC 1566","changing-state AGN","radio-quiet AGN","millimeter continuum","X-ray corona","synchrotron emission","ALMA Band 6","accretion state"],"falsifier":"Recompute the Spearman test on the eight pairs in Table A.2; if the p-value is 0.168 rather than 0.0168, the claimed correlation is not significant. Alternatively, obtain simultaneous ALMA and X-ray observations across a changing-state transition and test whether the 0.05 dex scatter and flat spectral index persist on daily timescales.","tokens_in":14737,"feed_emoji":"📡","tokens_out":9043,"duration_ms":80161,"temperature":0.7,"pith_summary":"This paper uses eight epochs of ALMA 230 GHz observations of NGC 1566, a nearby radio-quiet galaxy that changed spectral state, to argue that the compact millimeter continuum is synchrotron radiation from the magnetized X-ray corona. The central evidence is a positive correlation between the 230 GHz flux and the 14-150 keV X-ray flux (Spearman $\\rho=0.54$, intrinsic scatter 0.05 dex) together with a flat in-band spectral index $\\alpha_{\\rm mm}$ that marks the source as compact and optically thick. The authors exclude star formation, outflow-driven shocks, and a compact jet as dominant contributors, and conclude that the corona is the most probable origin, making millimeter emission a practical tracer of accretion activity in radio-quiet active galactic nuclei.","feed_headline":"230 GHz light in NGC 1566 tracks its X-ray corona","feed_subtitle":"Eight ALMA epochs show mm flux tracking X-rays, tying radio-quiet AGN millimeter light to the magnetized corona.","key_machinery":"The mechanism at work is the magnetized X-ray corona: a compact plasma of hot electrons near the supermassive black hole in which magnetic reconnection heats the plasma, Compton scattering up-scatters UV seed photons to X-ray energies, and non-thermal electrons radiate synchrotron emission that becomes self-absorbed and peaks near 100-300 GHz. The paper identifies the 230 GHz emission as this coronal synchrotron component using three linked diagnostics: the flat in-band spectral index $\\alpha_{\\rm mm}$, the absence of a correlation between $\\alpha_{\\rm mm}$ and flux, and a self-absorbed size estimate from the Laor & Behar (2008) formula, which gives $R_{\\rm mm}\\sim100{-}200\\,R_g$ and an equipartition magnetic field $B_{\\rm eq}\\sim100{-}300$ G.","core_discovery":"The paper's central claim is that NGC 1566's unresolved 230 GHz emission is produced by self-absorbed synchrotron radiation in the same magnetized corona that generates its 14-150 keV X-rays, and not by a jet, outflow shock, or star formation. Across the 2014-2023 changing-state cycle the 230 GHz flux rose about a factor of two during the 2018 outburst and later faded, tracking the X-ray flux; the in-band spectral index stayed flat in the range $-0.3$ to $0.1$, as expected for an optically thick synchrotron source; and the inferred self-absorbed source size ($R_{\\rm mm}\\sim100{-}200\\,R_g$, or $3{-}7\\times10^{-5}$ pc) is comparable to the X-ray corona size obtained from earlier spectral fitting. The paper also finds that the millimeter-to-X-ray flux ratio is higher at low Eddington ratios, which it interprets as easier particle acceleration in lower-density accretion flows.","pith_inferences":["With only eight non-simultaneous epochs, the tight 0.05 dex scatter should be tested by simultaneous, daily-cadence ALMA and X-ray monitoring through a future state transition; time lags between bands could change the apparent correlation.","The paper text gives the correlation's p-value as 0.0168, while Figure 3 prints 1.68e-01; until this discrepancy is resolved, the statistical significance of the millimeter-X-ray correlation should be regarded as unsettled.","If the coronal origin holds, 230 GHz variability could become a practical way to measure coronal magnetic fields and particle acceleration efficiency in radio-quiet AGNs without deep X-ray spectroscopy."],"forward_implications":["Millimeter monitoring can serve as a tracer of accretion state in radio-quiet AGNs, even when X-ray monitoring is unavailable or the nucleus is obscured.","The corona's 230 GHz emission is a stable, optically thick component across the state change, so brightness variability reflects changes in source size, particle density, or magnetic field rather than a steepening electron spectrum.","The higher millimeter-to-X-ray ratio at low Eddington ratio implies the corona's non-thermal electron content is not simply proportional to accretion power, consistent with more efficient particle acceleration in low-density flows.","The opposite behavior from black hole X-ray binaries (which quench jet radio emission in high-accretion states) argues against a compact jet origin for this radio-quiet source's millimeter light."],"supporting_citations":[{"why":"Supplies the 98-source X-ray versus 230 GHz correlation, with slope ~1.19 and scatter 0.30 dex, that this paper compares against.","marker":"Kawamuro et al. (2022)"},{"why":"Establishes the 100 GHz X-ray correlation with 0.22 dex scatter used as the population-level baseline.","marker":"Ricci et al. (2023)"},{"why":"Provides the theoretical coronal synchrotron model predicting millimeter emission peaking near 100 GHz with a flat self-absorbed spectrum.","marker":"Inoue & Doi (2014)"},{"why":"Supplies the self-absorbed synchrotron source-size formula and equipartition estimate used to size the millimeter-emitting region.","marker":"Laor & Behar (2008)"},{"why":"Gives the X-ray spectral modeling of NGC 1566 that yields the corona size and documents the 2018 outburst.","marker":"Jana et al. (2021)"},{"why":"Provides the radio-to-far-infrared coronal SED model and magnetic field constraints for comparison with the equipartition estimate.","marker":"del Palacio et al. (2025)"}],"fun_headline_variants":["Millimeter glow from NGC 1566's corona, not its jet","230 GHz tracks X-rays: corona, not jet, in NGC 1566","Corona powers NGC 1566's millimeter light, ALMA shows","Compact synchrotron from corona: NGC 1566's mm origin","NGC 1566's mm radiation reveals its X-ray corona"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the X-ray fluxes, which are not simultaneous with the ALMA observations and are converted from Swift/BAT and Swift/XRT bands assuming a fixed photon index $\\Gamma=1.8$ with 30-day binning, faithfully represent the X-ray state at each millimeter epoch; if the conversion or the time lag is biased, the reported correlation could weaken or disappear.","fun_headline_variants_meta":{"raw":{"variants":["Millimeter glow from NGC 1566's corona, not its jet","230 GHz tracks X-rays: corona, not jet, in NGC 1566","Corona powers NGC 1566's millimeter light, ALMA shows","Compact synchrotron from corona: NGC 1566's mm origin","NGC 1566's mm radiation reveals its X-ray corona"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000155,"raw_usage":{"total_tokens":1257,"prompt_tokens":1030,"completion_tokens":227,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":646,"completion_tokens_details":{"reasoning_tokens":130}},"tokens_in":646,"tokens_out":227,"duration_ms":2885,"temperature":1.0,"reasoning_tokens":130,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:16:19.299709+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the Spearman test on the eight pairs in Table A.2; if the p-value is 0.168 rather than 0.0168, the claimed correlation is not significant. Alternatively, obtain simultaneous ALMA and X-ray observations across a changing-state transition and test whether the 0.05 dex scatter and flat spectral index persist on daily timescales.","supporting_citations":[{"cited_title":"2022, , 938, 87","cited_arxiv_id":null,"evidence_quote":"Supplies the 98-source X-ray versus 230 GHz correlation, with slope ~1.19 and scatter 0.30 dex, that this paper compares against."},{"cited_title":"& Doi , A","cited_arxiv_id":null,"evidence_quote":"Provides the theoretical coronal synchrotron model predicting millimeter emission peaking near 100 GHz with a flat self-absorbed spectrum."},{"cited_title":"2021, , 507, 687","cited_arxiv_id":null,"evidence_quote":"Gives the X-ray spectral modeling of NGC 1566 that yields the corona size and documents the 2018 outburst."},{"cited_title":"Millimeter emission from supermassive black hole coronae","cited_arxiv_id":"2504.07762","evidence_quote":"Provides the radio-to-far-infrared coronal SED model and magnetic field constraints for comparison with the equipartition estimate."}],"review_version":1}