pith. sign in

arxiv: 2509.06933 · v2 · submitted 2025-09-08 · 🌌 astro-ph.HE · astro-ph.GA

Black-hole mass estimation through accretion disk spectral fitting for high-redshift blazars

G. Kyriopoulos , M. Petropoulou , G. Vasilopoulos , D. Hatzidimitriou This is my paper

Pith reviewed 2026-05-18 17:51 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.GA
keywords black hole masshigh-redshift blazarsaccretion diskIGM attenuationspectral energy distributionEddington ratioMCMC fittingShakura-Sunayev model
0
0 comments X

The pith

Neglecting intergalactic gas absorption overestimates black hole masses in high-redshift blazars.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper seeks to measure black hole masses and accretion rates for 23 high-redshift blazars by fitting their infrared-to-ultraviolet spectral energy distributions with a multi-temperature accretion disk model inside a single Markov Chain Monte Carlo framework. It explicitly includes the dimming effect of neutral hydrogen in the intergalactic medium that grows stronger at redshifts above 2. A sympathetic reader would care because these masses and Eddington ratios shape our picture of how the most massive black holes assembled and powered jets in the early universe. The fits show that skipping the IGM correction produces systematically higher masses and lower Eddington ratios, with the offset growing at larger distances. The work also finds that black-hole spin assumptions create a degeneracy in the results and that uniform methods reduce scatter when estimates from different studies are compared.

Core claim

By performing Bayesian MCMC fits of the Shakura-Sunayev multi-temperature accretion disk model to the compiled IR-UV photometry of 23 high-z blazars while accounting for IGM attenuation, the authors obtain median black-hole masses in the range 10^8 to 10^10 solar masses and a broad distribution of Eddington ratios from approximately 0.04 to 1. They demonstrate that omitting IGM attenuation leads to larger mass estimates and correspondingly lower Eddington ratios, with the bias becoming more pronounced at higher redshift. Individual mass values can differ by up to a factor of a few from literature values obtained with other techniques, and assumptions about black-hole spin introduce a clear,,

What carries the argument

The Shakura-Sunayev multi-temperature blackbody accretion disk model fitted via MCMC to the observed photometric SED after applying an IGM neutral-hydrogen attenuation correction.

If this is right

  • Black-hole masses for the sample span roughly 10^8 to 10^10 solar masses.
  • Eddington ratios range from sub-Eddington values near 0.04 up to near-Eddington values close to 1.
  • Omitting the IGM attenuation correction systematically inflates black-hole mass estimates and depresses Eddington ratios, with the size of the bias increasing at higher redshift.
  • Black-hole spin assumptions create a systematic degeneracy between mass and accretion rate in the fits.
  • Uniform statistical frameworks are required to avoid method-dependent offsets when black-hole estimates from different samples or techniques are compared.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same fitting approach could be applied to lower-redshift blazar samples to test whether accretion properties evolve with cosmic time.
  • Future ultraviolet spectroscopy that bypasses IGM absorption would tighten the mass and accretion-rate constraints for the highest-redshift objects.
  • If jet synchrotron contributes non-negligibly in the UV for some sources, the reported masses would shift downward once an explicit jet component is added to the model.

Load-bearing premise

The infrared-to-ultraviolet photometric data are dominated by thermal emission from the multi-temperature accretion disk rather than non-thermal emission from the relativistic jet.

What would settle it

An independent mass measurement for any one of the 23 blazars obtained through reverberation mapping or dynamical modeling that lies well outside the 1-sigma interval reported by the SED fit.

Figures

Figures reproduced from arXiv: 2509.06933 by D. Hatzidimitriou, G. Kyriopoulos, G. Vasilopoulos, M. Petropoulou.

Figure 1
Figure 1. Figure 1: Effective opacity (averaged over an ensemble of line of sights) due to neutral hydrogen gas in the intergalactic medium as a function of redshift for different wavelengths corresponding to the six Swift-UVOT filters. Results obtained from the analytical model of Inoue et al. (2014b) are plotted with solid lines. Results adopted from [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Infrared-to-ultraviolet SEDs of three high-z blazars (in the observer’s frame). Green points indicate publicly available pho￾tometric observations from the ASI/SSDC SED Builder (Stratta et al. 2011). Upper limits are plotted with black arrows. Model fits with and without IGM extinction are plotted with blue and orange colors, respectively. Dark and light shaded regions indicate the 68 and 90 per cent confi… view at source ↗
Figure 3
Figure 3. Figure 3: Decomposition of the total (jet+disk) model SEDs presented in [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Reduced corner plots for four blazars showing the posterior distributions for the pair of fitted parameters (MBH, M˙ ). The distribution of the Eddington ratio, which is a derived quantity, is also shown for comparison. Colors have the same meaning as in [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Comparison between our MCMC estimates (including extinction) and literature values for black hole mass MBH (left column) and Eddington ratio λEdd (right column). Top left panel: black hole mass scatter plot with points colored by literature (author+year) and error bars showing (asymmetric) MCMC 68% uncertainties. Top right panel: Eddington ratio scatter plot with points colored by literature (author+year) … view at source ↗
Figure 6
Figure 6. Figure 6: Top panel: SED of GB6 0805+614 with MCMC modeling re￾sults for a black hole with a = 0.998 and a = 0 overplotted in blue and orange colors, respectively. Bottom panel: Corner plot show￾ing the posterior distributions for M˙ , MBH and λEdd for a black hole with a = 0.998 and a = 0. Same color coding as in the top panel is used. Equation (17) thus returns the mathematically minimal seed mass consistent with … view at source ↗
Figure 7
Figure 7. Figure 7: Black hole growth history for the investigated objects. The shaded horizontal regions correspond to different black hole seeding channels: Pop III remnants (Mseed ≲ 102 M⊙, grey), stellar dynamical processes (Mseed ∼ 103–104 M⊙, cyan), and direct collapse (Mseed ∼ 104–106 M⊙, orange). These seed mass ranges are based on the classification described in Valiante et al. (2016). Parts of tracks with Mseed < 1 … view at source ↗
read the original abstract

High-redshift ($z>2$) blazars, with relativistic jets aligned toward us, probe the most powerful end of the active galactic nuclei (AGN) population. We aim at determining the black hole masses and mass accretion rates of high-$z$ blazars in a common framework that utilizes a Markov Chain Monte Carlo (MCMC) fitting method and the Shakura-Sunayev multi-temperature accretion disk model, accounting also for attenuation due to neutral hydrogen gas in the intergalactic medium (IGM). We compiled a sample of 23 high-redshift blazars from the literature with publicly available infrared-to-ultraviolet photometric data. We performed a Bayesian fit to the spectral energy distribution (SED) of the accretion disk, accounting for upper limits, and determined the black hole masses and mass accretion rates with their uncertainties. We also examined the impact of optical-ultraviolet attenuation due to gas in the IGM. We find that neglecting IGM attenuation in SED fits leads to systematically larger black-hole mass estimates and correspondingly lower Eddington ratios, with the bias becoming more severe at higher redshift. Our MCMC fits yield median black-hole masses in the range $\sim (10^{8}-10^{10})\,M_{\odot}$ and a broad distribution of median Eddington ratios ($\lambda_{\rm Edd}\sim0.04$ up to $\sim1$). Comparison with previous literature shows no clear method-dependent systematic offsets, although individual mass estimates can differ by up to a factor of a few. We also demonstrate that assumptions about black-hole spin introduce a systematic degeneracy. This work is to our knowledge the first systematic study to model the accretion-disk emission of a large sample of high-$z$ blazars within a single, consistent statistical framework. Our results emphasize the importance of accounting for IGM attenuation and of using uniform fitting methods when comparing disk-based black hole estimates across samples.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The paper presents a Bayesian MCMC analysis fitting the Shakura-Sunayev multi-temperature thin accretion disk model to the compiled IR-to-UV photometric SEDs of 23 high-redshift (z>2) blazars, explicitly including IGM attenuation as a multiplicative factor and handling upper limits. It reports that neglecting IGM attenuation produces systematically larger black-hole mass estimates and lower Eddington ratios, with the bias increasing at higher redshift; the fits yield median masses in the range ∼(10^8−10^10) M⊙ and Eddington ratios from ∼0.04 to ∼1. The work positions itself as the first systematic study applying this uniform statistical framework to a large high-z blazar sample and notes a spin-related degeneracy.

Significance. If the disk-dominance assumption holds, the results would usefully demonstrate the quantitative impact of IGM attenuation on high-z disk-based mass estimates and provide a consistent reference sample for comparing black-hole growth in the early universe. The MCMC approach with proper upper-limit treatment and the explicit comparison of fits with and without IGM correction are positive methodological features.

major comments (1)
  1. [Fitting procedure / abstract description] The central claim that omitting IGM attenuation leads to systematically inflated M_BH and reduced λ_Edd (with redshift-dependent severity) requires that every photometric point in the likelihood is attributable to the multi-temperature disk continuum. The fitting procedure (abstract and methods description) applies the Shakura-Sunayev model directly to the IR-UV photometry without an additional jet synchrotron component or quantitative test for non-thermal contamination. In blazars, jet emission commonly extends into the near-IR/optical; if present at even a subset of the data points, the MCMC will compensate by shifting M_BH and Ṁ, thereby altering both the absolute mass scale and the reported differential IGM bias. This assumption is load-bearing for the main result and needs explicit justification or modeling.
minor comments (1)
  1. [Abstract] The abstract states that assumptions about black-hole spin introduce a systematic degeneracy; a brief quantitative illustration of how the reported mass range changes under different spin priors would strengthen the presentation.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive feedback on our manuscript. We respond to the major comment below.

read point-by-point responses

  1. Referee: [Fitting procedure / abstract description] The central claim that omitting IGM attenuation leads to systematically inflated M_BH and reduced λ_Edd (with redshift-dependent severity) requires that every photometric point in the likelihood is attributable to the multi-temperature disk continuum. The fitting procedure (abstract and methods description) applies the Shakura-Sunayev model directly to the IR-UV photometry without an additional jet synchrotron component or quantitative test for non-thermal contamination. In blazars, jet emission commonly extends into the near-IR/optical; if present at even a subset of the data points, the MCMC will compensate by shifting M_BH and Ṁ, thereby altering both the absolute mass scale and the reported differential IGM bias. This assumption is load-bearing for the main result and needs explicit justification or modeling.

    Authors: We agree that the disk-dominance assumption is central and benefits from explicit discussion. Our analysis follows the standard approach in the literature for disk-based mass estimates in high-z AGN, where the UV portion of the SED is typically attributed to thermal disk emission (the big blue bump). We will revise the Methods and Discussion sections to include a dedicated justification of this assumption, supported by references to prior blazar SED studies that decompose thermal and non-thermal components and find disk dominance in the rest-frame UV for similar sources. We will also add a robustness test by repeating the MCMC fits for a subset of objects after excluding the near-IR photometric points most susceptible to synchrotron contamination, and we will report the resulting changes (or lack thereof) to the median M_BH values and the magnitude of the IGM-induced bias. These additions will be made without altering the primary conclusions on the importance of IGM attenuation. revision: yes

Circularity Check

0 steps flagged

No circularity: results obtained by direct MCMC fitting of external model to independent photometry

full rationale

The paper compiles independent IR-UV photometric data for 23 high-z blazars and performs Bayesian MCMC fits of the Shakura-Sunayev multi-temperature disk model, both with and without an IGM attenuation factor. The reported systematic offset in black-hole masses and Eddington ratios is the direct numerical outcome of these two comparative fits applied to the same external observations. No equation reduces a derived quantity to a fitted parameter by construction, no self-citation chain supplies a load-bearing premise, and no ansatz or uniqueness theorem is imported from prior author work. The central claim therefore rests on the application of a standard external model to independent data rather than on any self-referential loop.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central results rest on the applicability of the thin-disk model to high-redshift blazars and the assumption that the observed SED can be attributed primarily to disk emission after IGM correction.

free parameters (2)
  • Black hole mass M_BH
    Primary fitted parameter in the MCMC SED modeling that sets the overall normalization and temperature profile of the disk.
  • Mass accretion rate
    Fitted parameter that determines the disk luminosity and Eddington ratio for each source.
axioms (2)
  • domain assumption Shakura-Sunayev multi-temperature thin accretion disk model accurately describes the IR-UV SED
    Invoked when selecting the model for Bayesian fitting to the photometric data.
  • domain assumption Neutral hydrogen in the IGM produces the dominant attenuation that can be modeled as a simple absorption correction
    Used to modify the model SED before comparison with observations.

pith-pipeline@v0.9.0 · 5903 in / 1679 out tokens · 50993 ms · 2026-05-18T17:51:10.630614+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

45 extracted references · 45 canonical work pages · 1 internal anchor

  1. [1]

    M., Press, W

    Bardeen, J. M., Press, W. H., & Teukolsky, S. A. 1972, ApJ, 178, 347

    work page 1972

  2. [2]

    C., Volonteri, M., & Rees, M

    Begelman, M. C., Volonteri, M., & Rees, M. J. 2006, MNRAS, 370, 289

    work page 2006

  3. [3]

    2022, A&A, 660, A74

    Belladitta, S., Caccianiga, A., Diana, A., et al. 2022, A&A, 660, A74

    work page 2022

  4. [4]

    N., Hill, J

    Burrows, D. N., Hill, J. E., Nousek, J. A., et al. 2005, Space Sci. Rev., 120, 165

    work page 2005

  5. [5]

    2013, MNRAS, 431, 210

    Calderone, G., Ghisellini, G., Colpi, M., & Dotti, M. 2013, MNRAS, 431, 210

    work page 2013

  6. [6]

    A., Clayton, G

    Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245

    work page 1989

  7. [7]

    W., Lang, D., & Goodman, J

    Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP , 125, 306

    work page 2013

  8. [8]

    Frank, J., King, A., & Raine, D. J. 2002, Accretion Power in Astrophysics: Third Edition

    work page 2002

  9. [9]

    2016, Galaxies, 4, 36

    Ghisellini, G. 2016, Galaxies, 4, 36

    work page 2016

  10. [10]

    2010, MNRAS, 405, 387

    Ghisellini, G., Della Ceca, R., Volonteri, M., et al. 2010, MNRAS, 405, 387

    work page 2010

  11. [11]

    2013, MNRAS, 428, 1449

    Ghisellini, G., Nardini, M., Tagliaferri, J., G., et al. 2013, MNRAS, 428, 1449

    work page 2013

  12. [12]

    2011, MNRAS, 411, 901

    Ghisellini, G., Tagliaferri, G., Foschini, L., et al. 2011, MNRAS, 411, 901

    work page 2011

  13. [13]

    2015, MNRAS, 450, L34

    Ghisellini, G., Tagliaferri, G., Sbarrato, T., & Gehrels, N. 2015, MNRAS, 450, L34

    work page 2015

  14. [14]

    2014, Nature, 515, 376

    Ghisellini, G., Tavecchio, F., Maraschi, L., Celotti, A., & Sbarrato, T. 2014, Nature, 515, 376

    work page 2014

  15. [15]

    Gregory, P . C. 2005, Bayesian Logical Data Analysis for the Physical Sci- ences: A Comparative Approach with ‘Mathematica’ Support

    work page 2005

  16. [16]

    2008, PASP , 120, 405

    Greiner, J., Bornemann, W., Clemens, C., et al. 2008, PASP , 120, 405

    work page 2008

  17. [17]

    Gunn, J. E. & Peterson, B. A. 1965, ApJ, 142, 1633

    work page 1965

  18. [18]

    Haakonsen, C. B. & Rutledge, R. E. 2009, ApJS, 184, 138

    work page 2009

  19. [19]

    & Madau, P

    Haardt, F. & Madau, P . 2012, ApJ, 746, 125 HI4PI Collaboration, Ben Bekhti, N., Flöer, L., et al. 2016, A&A, 594, A116

    work page 2012

  20. [20]

    Hogg, D. W. 2000, Distance measures in cosmology

    work page 2000

  21. [21]

    Data analysis recipes: Fitting a model to data

    Hogg, D. W., Bovy, J., & Lang, D. 2010, arXiv e-prints, arXiv:1008.4686

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  22. [22]

    & Lindfors, E

    Hovatta, T. & Lindfors, E. 2019, New A Rev., 87, 101541

    work page 2019

  23. [23]

    2020, Annual Review of Astronomy and Astrophysics, 58, 27–97

    Inayoshi, K., Visbal, E., & Haiman, Z. 2020, Annual Review of Astronomy and Astrophysics, 58, 27–97

    work page 2020

  24. [24]

    H., Cohen, M., Masci, F., et al

    Jarrett, T. H., Cohen, M., Masci, F., et al. 2011, ApJ, 735, 112

    work page 2011

  25. [25]

    2016, MNRAS, 456, 2993

    Lupi, A., Haardt, F., Dotti, M., et al. 2016, MNRAS, 456, 2993

    work page 2016

  26. [26]

    2020, ApJ, 889, 164

    Marcotulli, L., Paliya, V ., Ajello, M., et al. 2020, ApJ, 889, 164

    work page 2020

  27. [27]

    Narayan, R. & Yi, I. 1995, ApJ, 452, 710

    work page 1995

  28. [28]

    Novikov, I. D. & Thorne, K. S. 1973, in Black Holes (Les Astres Occlus), ed. C. Dewitt & B. S. Dewitt, 343–450

    work page 1973

  29. [29]

    2016, MNRAS, 458, 3047 Planck Collaboration, Ade, P

    Pezzulli, E., Valiante, R., & Schneider, R. 2016, MNRAS, 458, 3047 Planck Collaboration, Ade, P . A. R., Aghanim, N., et al. 2016, A&A, 594, A13

    work page 2016

  30. [30]

    Reines, A. E. & Comastri, A. 2016, PASA, 33, e054

    work page 2016

  31. [31]

    Roming, P . W. A., Kennedy, T. E., Mason, K. O., et al. 2005, Space Sci. Rev., 120, 95

    work page 2005

  32. [32]

    Salpeter, E. E. 1964, ApJ, 140, 796 Sanchez Zaballa, J. M., Bottacini, E., & Tramacere, A. 2025, ApJ, 985, 179

    work page 1964

  33. [33]

    2012, PASP , 124, 1208

    Sawicki, M. 2012, PASP , 124, 1208

    work page 2012

  34. [34]

    2025, Space Sci

    Sbarrato, T., Ajello, M., Buson, S., et al. 2025, Space Sci. Rev., 221, 62

    work page 2025

  35. [35]

    2015, MNRAS, 446, 2483 Article number, page 10 of 15 G

    Sbarrato, T., Ghisellini, G., Tagliaferri, G., et al. 2015, MNRAS, 446, 2483 Article number, page 10 of 15 G. Kyriopoulos , M. Petropoulou , G. Vasilopoulos , and D. Hatzidimitriou : Black-hole masses of high-z blazars

    work page 2015

  36. [36]

    2013, ApJ, 777, 147

    Sbarrato, T., Tagliaferri, G., Ghisellini, G., et al. 2013, ApJ, 777, 147

    work page 2013

  37. [37]

    Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337

    work page 1973

  38. [38]

    T., Strauss, M

    Shen, Y., Richards, G. T., Strauss, M. A., et al. 2011, ApJS, 194, 45 S ˛ adowski, A. & Narayan, R. 2016, MNRAS, 456, 3929

    work page 2011

  39. [39]

    2011, arXiv e-prints, arXiv:1103.0749

    Stratta, G., Capalbi, M., Giommi, P ., et al. 2011, arXiv e-prints, arXiv:1103.0749

    work page arXiv 2011

  40. [40]

    2016, in Active Galactic Nuclei 12: A Multi-Messenger Perspective (AGN12), 4

    Valiante, R., Schneider, R., Volonteri, M., & Omukai, K. 2016, in Active Galactic Nuclei 12: A Multi-Messenger Perspective (AGN12), 4

    work page 2016

  41. [41]

    2010, A&A Rev., 18, 279

    Volonteri, M. 2010, A&A Rev., 18, 279

    work page 2010

  42. [42]

    2011, MNRAS, 416, 216

    Volonteri, M., Haardt, F., Ghisellini, G., & Della Ceca, R. 2011, MNRAS, 416, 216

    work page 2011

  43. [43]

    L., Eisenhardt, P

    Wright, E. L., Eisenhardt, P . R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

    work page 2010

  44. [44]

    2020, MNRAS, 491, 740

    Yang, G., Boquien, M., Buat, V ., et al. 2020, MNRAS, 491, 740

    work page 2020

  45. [45]

    & Miralda-Escudé, J

    Yoo, J. & Miralda-Escudé, J. 2004, ApJ, 614, L25 Article number, page 11 of 15 A&A proofs:manuscript no. paper Appendix A: Representative cases of problematic photometry To illustrate the types of photometric issues that motivated our visual inspection and filtering, we present two rep- resentative examples in Fig. A.1: (i) a catalog misassociation in the...

    work page 2004