The Millimeter/X-ray Relation in Rapidly Accreting Supermassive Black Holes at z <0.16
Pith reviewed 2026-05-19 15:51 UTC · model grok-4.3
The pith
High-luminosity AGN lie above the millimeter-X-ray correlation defined by lower-luminosity sources.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
These high-luminosity AGN lie above the mm/X-ray correlation defined by lower-luminosity sources. A joint fit to both samples yields a second-degree polynomial with an intrinsic scatter of 0.32 dex. Furthermore, the mm emission correlates linearly with both the UV disk luminosity and L_bol, with intrinsic scatters of 0.45 and 0.35 dex, respectively. We propose that the deviation from the linear mm/X-ray relation arises from a two-component coronal electron population: thermal electrons that produce X-rays, but become less efficient at higher luminosities, and non-thermal electrons that produce mm emission and remain tied to L_bol.
What carries the argument
two-component coronal electron population: thermal electrons for X-ray production that become less efficient at high luminosities and non-thermal electrons for millimeter emission tied to bolometric luminosity
If this is right
- The X-ray corona becomes less efficient relative to millimeter output as bolometric luminosity rises.
- Millimeter emission remains a linear tracer of bolometric luminosity across a wide range of Eddington ratios.
- The combined sample is described by a quadratic rather than linear mm/X-ray relation.
- Additional millimeter contributions from outflows are possible but not required by the current SED and spectral data.
Where Pith is reading between the lines
- Future observations at multiple millimeter frequencies could confirm whether the emission remains unresolved and non-thermal across luminosities.
- The proposed shift in electron populations may alter predictions for X-ray variability amplitudes in luminous AGN.
- If the two-component picture is correct, millimeter surveys could provide independent estimates of bolometric luminosity for distant AGN where X-ray data are limited.
Load-bearing premise
The millimeter excess originates primarily from a non-thermal coronal electron population rather than from outflow-driven shocks or other extended structures.
What would settle it
High-resolution imaging or spectral index measurements showing resolved extended emission or a spectral index matching shock processes rather than coronal synchrotron would undermine the coronal two-component interpretation.
Figures
read the original abstract
A tight correlation between nuclear millimeter (mm) and X-ray emission has recently been found in nearby ($z < 0.01$) and low-Eddington ratio ($\rm \lambda_{Edd} < 0.1$) radio-quiet (RQ) Active Galactic Nuclei (AGN), suggesting a common origin in the hot X-ray corona. We test this relation in nine more distant RQ AGN ($z \sim 0.06-0.16$) with higher bolometric luminosities ($\log(L_{\rm bol}/\mathrm{erg\,s^{-1}})=45.3-46.3$), Eddington ratios ($\rm \lambda_{Edd} = 0.19-0.85$), and X-ray bolometric corrections ($\kappa_{2-10}=29-194$), selected from the Burst Alert Telescope (BAT) survey. We obtained quasi-simultaneous observations with Swift at 2-10 keV and the Atacama Large Millimeter/submillimeter Array (ALMA) at 100 GHz and with high angular resolution ($<0.14$"). We find that these high-luminosity AGN lie above the mm/X-ray correlation defined by lower-luminosity sources. A joint fit to both samples yields a second-degree polynomial with an intrinsic scatter of 0.32 dex. Furthermore, the mm emission correlates linearly with both the UV disk luminosity and $L_{\rm bol}$, with intrinsic scatters of 0.45 and 0.35 dex, respectively. We propose that the deviation from the linear mm/X-ray relation arises from a two-component coronal electron population: thermal electrons that produce X-rays, but become less efficient at higher luminosities, and non-thermal electrons that produce mm emission and remain tied to $L_{\rm bol}$. Additional mm emission from outflow-driven shocks may also contribute, though SED modeling and spectral index studies favor a coronal origin.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports quasi-simultaneous ALMA 100 GHz and Swift 2-10 keV observations of nine radio-quiet AGN at z = 0.06-0.16 with high bolometric luminosities (log L_bol = 45.3-46.3 erg s^{-1}) and Eddington ratios (0.19-0.85). These sources lie above the mm/X-ray correlation defined by lower-luminosity nearby AGN. A joint fit to both samples produces a second-degree polynomial with 0.32 dex intrinsic scatter. The mm emission correlates linearly with UV disk luminosity and L_bol (scatters 0.45 and 0.35 dex). The authors interpret the deviation via a two-component coronal electron population (thermal electrons for X-rays becoming less efficient at high L, non-thermal for mm tied to L_bol), while noting possible outflow contributions but favoring coronal origin from SED and spectral-index arguments.
Significance. If the observational deviation and polynomial fit hold, the result would indicate a change in coronal properties at high accretion rates, with implications for AGN corona models and multi-wavelength scaling relations. Strengths include the quasi-simultaneous high-resolution ALMA data isolating nuclear emission and the reported linear mm-L_bol relation as an empirical benchmark for luminous sources.
major comments (2)
- [Abstract] Abstract: the central claim that the nine high-luminosity sources lie above the existing mm/X-ray relation and that a joint second-degree polynomial fit yields 0.32 dex scatter is load-bearing, yet the manuscript provides no details on data reduction, error budgets, sample selection completeness, or the fitting procedure used to obtain the polynomial coefficients and intrinsic scatter values.
- [Abstract] Abstract: the proposed two-component coronal interpretation (thermal electrons less efficient at high luminosity, non-thermal electrons producing mm emission) is load-bearing for the physical model, but no quantitative bound is given on the allowable non-coronal contribution (e.g., upper limit on extended flux fraction or Bayes factor) from outflow-driven shocks or host structures; at z ~ 0.1 even a 10-20% extended component within the <0.14 arcsec beam would shift points enough to change the quadratic coefficient.
minor comments (1)
- The abstract states that SED modeling and spectral index studies favor a coronal origin; including the specific quantitative results or references from those analyses in the main text would improve clarity.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for identifying areas where additional clarity would strengthen the presentation. We address each major comment below and indicate the revisions we will make in the next version of the paper.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim that the nine high-luminosity sources lie above the existing mm/X-ray relation and that a joint second-degree polynomial fit yields 0.32 dex scatter is load-bearing, yet the manuscript provides no details on data reduction, error budgets, sample selection completeness, or the fitting procedure used to obtain the polynomial coefficients and intrinsic scatter values.
Authors: We note that the full manuscript contains dedicated sections describing the BAT-selected sample (Section 2), ALMA 100 GHz data reduction, calibration, and imaging with the <0.14 arcsec beam (Section 3.1), Swift 2-10 keV data processing and error estimation (Section 3.2), and the Bayesian MCMC fitting procedure used to derive the quadratic coefficients and 0.32 dex intrinsic scatter (Section 4.3). However, these details are not summarized in the abstract. We will revise the abstract to include a concise statement referencing the key methodological steps and will add a short methods overview paragraph early in the text to make the load-bearing claims more self-contained. revision: yes
-
Referee: [Abstract] Abstract: the proposed two-component coronal interpretation (thermal electrons less efficient at high luminosity, non-thermal electrons producing mm emission) is load-bearing for the physical model, but no quantitative bound is given on the allowable non-coronal contribution (e.g., upper limit on extended flux fraction or Bayes factor) from outflow-driven shocks or host structures; at z ~ 0.1 even a 10-20% extended component within the <0.14 arcsec beam would shift points enough to change the quadratic coefficient.
Authors: We agree that an explicit quantitative limit on possible non-coronal (outflow or host) contributions would strengthen the argument. The current manuscript relies on the sources remaining unresolved within the high-resolution ALMA beam, mm spectral indices, and SED modeling to favor a coronal origin, but does not provide a numerical upper bound on extended flux. We will add a new paragraph in the discussion section that uses the beam size and comparison with lower-resolution archival data to constrain any extended component to ≲15% for the sample; this level of contamination is insufficient to remove the quadratic deviation. A full Bayes-factor model comparison between pure-corona and corona-plus-outflow scenarios is not feasible with the existing dataset and would require additional multi-frequency modeling, which we will note as a limitation. revision: partial
Circularity Check
No significant circularity; empirical correlation from independent data
full rationale
The paper reports new quasi-simultaneous ALMA 100 GHz and Swift 2-10 keV observations of nine high-luminosity RQ AGN at z=0.06-0.16, measures their fluxes, places them relative to the existing lower-luminosity mm/X-ray relation, performs a joint second-degree polynomial fit to the combined samples, and reports linear mm-L_bol and mm-UV correlations. These steps are direct empirical measurements and statistical fits to fresh data; none reduce by the paper's own equations or self-citations to a quantity defined by the same inputs. The two-component corona proposal is presented as an interpretive suggestion supported by cited SED and spectral-index arguments, not as a mathematical derivation that loops back to the fitted parameters or sample selection by construction.
Axiom & Free-Parameter Ledger
free parameters (2)
- second-degree polynomial coefficients
- intrinsic scatter values (0.32, 0.45, 0.35 dex)
axioms (2)
- domain assumption Quasi-simultaneous ALMA and Swift observations capture the same coronal state without significant variability between epochs.
- domain assumption Bolometric luminosity and Eddington ratio estimates from BAT survey data are accurate enough to place sources in the high-luminosity regime.
invented entities (1)
-
two-component coronal electron population (thermal + non-thermal)
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
A joint fit to both samples yields a second-degree polynomial with an intrinsic scatter of 0.32 dex. ... two-component coronal electron population: thermal electrons that produce X-rays, but become less efficient at higher luminosities, and non-thermal electrons that produce mm emission and remain tied to L_bol.
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We find that these high-luminosity AGN lie above the mm/X-ray correlation defined by lower-luminosity sources.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
2012, The Astrophysical Journal, 747, 104, doi: 10.1088/0004-637X/747/2/104 Arnaud, K
Ackermann, M., Ajello, M., Allafort, A., et al. 2012, The Astrophysical Journal, 747, 104, doi: 10.1088/0004-637X/747/2/104 Arnaud, K. A. 1996, 101,
-
[2]
https://ui.adsabs.harvard.edu/abs/1996ASPC..101...17A Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481, doi: 10.1146/annurev.astro.46.060407.145222 Avni, Y ., & Tananbaum, H. 1982, The Astrophysical Journal, 262, L17, doi: 10.1086/183903 —. 1986, The Astrophysical Journal, 305, 83, doi: 10.1086/164230 Baldi, R. D., Behar, E., Lao...
-
[3]
https://ui.adsabs.harvard.edu/abs/1983QJRAS..24..267H Hwang, H.-C., Zakamska, N. L., Alexandroff, R. M., et al. 2018, Monthly Notices of the Royal Astronomical Society, 477, 830, doi: 10.1093/mnras/sty742 Ichikawa, K., Ricci, C., Ueda, Y ., et al. 2017, The Astrophysical Journal, 835, 74, doi: 10.3847/1538-4357/835/1/74 —. 2019, The Astrophysical Journal,...
-
[4]
2024, ApJ, 965, 68, doi: 10.3847/1538-4357/ad2fae
https://ui.adsabs.harvard.edu/abs/2007ASPC..376..127M Michiyama, T., Inoue, Y ., Doi, A., et al. 2024, The Astrophysical Journal, 965, 68, doi: 10.3847/1538-4357/ad2fae 17 Mullaney, J. R., Alexander, D. M., Goulding, A. D., & Hickox, R. C. 2011, Monthly Notices of the Royal Astronomical Society, 414, 1082, doi: 10.1111/j.1365-2966.2011.18448.x Nhat Ly, M....
-
[5]
Each image shows the beam in the lower left corner and the physical size of 0.′′1 in pc in the lower right corner. 19 A.2.Spectral index determination 1.96 1.98 2.00 2.02 log( / GHz) 1.20 1.15 1.10 1.05 log(S / mJy) Power-law fit = 0.16 ± 0.85 This work (a) Q 0119–286 1.96 1.98 2.00 2.02 log( / GHz) 0.09 0.10 0.11 0.12 0.13 0.14 0.15 log(S / mJy) Power-la...
work page 1949
-
[6]
The upper panel shows the spectra in counts
X-ray spectra of our nine sources at 2–10 keV . The upper panel shows the spectra in counts. The bottom panel shows the ratio of the data and the fitted model. C.ADDITIONAL TESTS ON THE MM/X-RAY RELATION C.1.The 0.3–2 keV energy range In addition to investigating the relation between the 2–10 keV and mm emission, we investigated the relation between 0.3– ...
work page 2017
-
[7]
2023 This work 42 43 44 45 log(L14 150 keV/erg s
Ricci et al. 2023 This work 42 43 44 45 log(L14 150 keV/erg s
work page 2023
- [8]
-
[9]
Left:The 100 GHz luminosityνL 100GHz in erg s−1 versus the intrinsic 0.3–2 keV luminosityL0.3−2keV in erg s−1 for sources by R23 and this work. We have fitted a second-degree polynomial relation to the mm and X-ray emission, presented in Equation C1, which has an intrinsic scatter of 0.46 dex.Right:The 100 GHz luminosityνL 100GHz in erg s−1 versus the 14–...
work page 2017
-
[10]
Finally, the SFR values for PG 0026+129 and PG 0052+251 were determined by Y
and were determined through IR SED decomposition. Finally, the SFR values for PG 0026+129 and PG 0052+251 were determined by Y . D´ıaz et al., (in preparation) from the strength of the PAH feature. 6.5 7.5 8.5 log(MBH/M ) 5.0 4.5 4.0 3.5 log( L100GHz/L2 10keV) This work Ricci et al. 2023 1 0 1 2 log(SFR/M yr
work page 2023
-
[11]
Left:The mm/X-ray luminosity ratio versus the black hole massM BH of the sources from this work and R23. We find no correlation between the two, obtaining a p-value of 0.03.Right:The mm/X-ray luminosity ratio versus the star-formation rate (SFR). Eight out of nine sources in this work only have upper limits on SFR. No correlation is found, obtaining a p-v...
work page 2016
-
[12]
power-law plus black body and ionized reflection usingXILLVER(Garc ´ıa & Kallman 2010), prominent atE>10 keV . Since a detailed fit is beyond the scope of this paper, we freeze the iron abundance and inclination of the disk in the model. The best-fit model has a statistic ofχ 2/dof=1043/913. Negative residuals are present in the soft X-ray band (E=0.5−2 k...
work page 2010
-
[13]
2 1 0 1 2 3 mm This work Ricci et al. 2023 1 0 1 2 log(SFR/M yr
work page 2023
-
[14]
Left:Spectral index (α mm) of the 100 GHz observations for this work, as obtained from fittingSmm ∝ν −αmm to the 100 GHz fluxes in mJy in the four separate SPWs, and R23 versus the Eddington ratioλEdd. We do not observe a correlation, obtaining a p-value of 0.17.Middle: αmm versus the bolometric luminosityL bol. Here, no correlation is observed either sin...
work page 2025
-
[15]
More details can be found in Section 5.3.3
indicates that the 100 GHz flux is dominated by a corona withr c =217±13 and logδ=−1.06±0.11. More details can be found in Section 5.3.3. PG 0052+251 has a less constrained SED, presented in Figure 4 as well, with VLA data in Bands C, X, and Q, VLBA data, and ALMA Band 5 data (187–202 GHz; Proposal ID 2023.1.01062.S; PI: F. Bauer). The 100 GHz flux could ...
work page 2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.