pith. sign in

arxiv: 2606.12076 · v1 · pith:D6IID2NWnew · submitted 2026-06-10 · 🌌 astro-ph.GA

X-ray luminosity function of Compton-thick AGN in the early Universe (z > 3). Robustness and biases of the CTK population

Pith reviewed 2026-06-27 09:29 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords Compton-thick AGNX-ray luminosity functionhigh-redshift AGNobscured AGNactive galactic nuclei24-micron photometryblack hole growthcosmic X-ray background
0
0 comments X

The pith

Compton-thick AGN make up a steady 17 percent of the AGN population from redshift 3 to 6.

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

This paper examines the X-ray luminosity function of Compton-thick AGN at redshifts above 3 using large X-ray selected samples from XMM-Newton and Chandra surveys. The authors first identify high-probability CTK candidates via X-ray spectral fitting and then use multiwavelength SED modeling to check their properties. For most sources the X-ray luminosities appear overestimated relative to infrared luminosities, so the team applies 24-micron photometry to derive infrared luminosities and update the hydrogen column density posteriors. This step lowers the number of sources classified as CTK and produces a more conservative luminosity function. The final result is that CTK AGN comprise 17 percent of the total AGN population between redshifts 3 and 6, with no statistically significant change in that fraction up to redshift 6.

Core claim

After updating hydrogen column density posteriors with infrared-informed priors for the full sample, the analysis finds that Compton-thick AGN constitute 17 per cent of the total AGN population at 3 < z < 6. This fraction shows no statistically significant evolution up to z approximately 6. The overall obscured AGN fraction increases at higher redshifts, but the stable CTK fraction suggests that typical host galaxy interstellar media at these epochs cannot produce Compton-thick obscuration levels.

What carries the argument

Infrared luminosity priors from 24-micron photometry used to update hydrogen column density posteriors from X-ray spectral fitting.

If this is right

  • The CTK phase remains a persistent component of black hole growth throughout cosmic history.
  • Incorporating infrared priors yields lower and more physically consistent CTK number densities than X-ray-only methods.
  • The increase in the overall obscured AGN fraction toward higher redshifts does not extend to the CTK regime.
  • Host-galaxy interstellar media at z greater than 3 lack the density needed to produce Compton-thick obscuration.

Where Pith is reading between the lines

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

  • Black-hole growth models must incorporate a constant contribution from the CTK phase across the first few billion years.
  • Deeper multiwavelength follow-up could test whether the 17 percent fraction persists at redshifts beyond 6.
  • X-ray-only CTK identifications at high redshift may carry systematic overestimation biases that infrared priors help correct.

Load-bearing premise

Infrared luminosities derived from 24-micron photometry provide unbiased priors that correctly update the hydrogen column density posteriors without introducing new systematic errors in CTK classification.

What would settle it

An independent high-redshift sample with direct hard X-ray spectroscopy or submillimeter data that measures a CTK fraction significantly different from 17 percent would falsify the central claim.

Figures

Figures reproduced from arXiv: 2606.12076 by Angel Ruiz, Ektoras Pouliasis, Ioannis Georgantopoulos.

Figure 1
Figure 1. Figure 1: Intrinsic X-ray luminosity in the 2–10 keV band versus [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Spectral energy distribution of source lid_1278 ( [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Intrinsic X-ray luminosity in the 2–10 keV band versus the monochromatic 6 [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Joint posterior distribution in the log NH −log L6 µm plane (blue shaded region) for source lid_1278, one of our CTK can￾didates. The infrared luminosity L6 µm is inferred from the X￾ray luminosity posterior via the Stern (2015) relation. The black dashed line and grey shaded region show the SED-inferred L6 µm value and its uncertainty, while the red dotted contours indicate the posterior distribution afte… view at source ↗
Figure 5
Figure 5. Figure 5: Observed distributions of hydrogen column density and X-ray luminosity for the full [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Probability of the selected CTK candidates of having [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: X-ray luminosity and absorption function for two redshift intervals ( [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: X-ray luminosity function for the CTK population (24 [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
read the original abstract

The population of Compton-thick (CTK) AGN represents a critical yet elusive phase in the growth of supermassive black holes. Constraining their abundance and evolution at high z is essential for understanding both SMBH growth and the origin of the cosmic X-ray background. We investigate the X-ray luminosity function (XLF) of CTK AGN at z > 3 using one of the largest available samples of X-ray-selected AGN at high z, containing 811 sources from XMM-Newton XXL-N and Chandra CCLS and CDF-S/N surveys. We first selected a subsample of ten high-probability CTK candidates, identified through x-ray spectral fitting. Their multiwavelength properties are examined through SED modelling to assess the robustness of their CTK classification. For most sources, the inferred X-ray luminosities appear overestimated when compared with their IR luminosities. After updating the NH posteriors with IR-informed priors, only three sources remain consistent with the CTK regime. To compute the XLF for the entire CTK AGN population, we used 24 microns photometry to estimate IR luminosities and update the X-ray posteriors for all the remaining sources. Incorporating IR priors systematically reduces the inferred CTK number densities, yielding a more conservative and physically consistent estimate of the XLF. We find that CTK AGN constitute 17 per cent of the total AGN population at 3 < z < 6, consistent with results at lower z. Our analysis reveals no statistically significant evolution in the CTK fraction up to z about 6, suggesting that the most heavily obscured accretion phase remains a persistent component of black hole growth throughout cosmic history. While the overall obscured AGN fraction (NH > 1e23 cm-2) increases toward higher redshifts, the stable CTK fraction supports the interpretation that, at these epochs, the interstellar medium in typical host galaxies cannot produce CTK levels of obscuration.

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

2 major / 2 minor

Summary. The manuscript analyzes the X-ray luminosity function of Compton-thick (CTK) AGN at z > 3 using a sample of 811 X-ray-selected sources from XMM-Newton XXL-N and Chandra surveys. Ten high-probability CTK candidates are first identified via X-ray spectral fitting; multiwavelength SED modeling and IR priors from 24μm photometry are then used to refine NH posteriors, reducing the CTK count to three. The same IR-prior update is applied to the full sample to derive the XLF, yielding the result that CTK AGN constitute 17% of the AGN population at 3 < z < 6 with no statistically significant evolution up to z ≈ 6.

Significance. If the result holds after validation of the IR priors, the work provides a valuable constraint on the persistent fraction of heavily obscured accretion at high redshift, with implications for SMBH growth models and the cosmic X-ray background. The large sample size and the conservative approach of updating posteriors with independent IR data are positive features that strengthen the analysis relative to purely X-ray-based studies.

major comments (2)
  1. [Results / posterior-update procedure] The derivation of the headline 17% CTK fraction (abstract and results section) rests on updating NH posteriors for all 811 sources using IR luminosities inferred from 24μm photometry. At 3 < z < 6 the observed 24μm band samples rest-frame ≲ 8 μm; the paper does not quantify how deviations in the AGN-torus versus star-formation contribution or in the IR–X-ray scaling relation from the low-z calibration affect the CTK classification probability, yet this step reduces the initial candidate list from ten to three and is applied uniformly to the remaining sources.
  2. [Discussion / evolution test] The claim of “no statistically significant evolution” in the CTK fraction up to z ≈ 6 (abstract and discussion) is not accompanied by the explicit statistical test, posterior comparison, or completeness-corrected number-density uncertainties that would allow the reader to assess whether the 17% value is distinguishable from lower-z measurements within the reported errors.
minor comments (2)
  1. [Abstract] The abstract states that “X-ray luminosities appear overestimated when compared with their IR luminosities” for most of the ten candidates, but does not report the quantitative offset or the exact IR–X-ray relation adopted.
  2. Table or figure captions should explicitly state the redshift bins used for the 3 < z < 6 CTK fraction and the total AGN comparison sample.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments, which highlight important aspects of our methodology and statistical presentation. We address each major comment point-by-point below. Revisions have been made to incorporate additional quantification and explicit tests as detailed in the responses.

read point-by-point responses
  1. Referee: [Results / posterior-update procedure] The derivation of the headline 17% CTK fraction (abstract and results section) rests on updating NH posteriors for all 811 sources using IR luminosities inferred from 24μm photometry. At 3 < z < 6 the observed 24μm band samples rest-frame ≲ 8 μm; the paper does not quantify how deviations in the AGN-torus versus star-formation contribution or in the IR–X-ray scaling relation from the low-z calibration affect the CTK classification probability, yet this step reduces the initial candidate list from ten to three and is applied uniformly to the remaining sources.

    Authors: We acknowledge that the 24 μm photometry at z > 3 probes rest-frame wavelengths where both torus and star-formation emission can contribute, and that the IR–X-ray scaling is calibrated at lower redshifts. Our IR luminosities were derived via full multiwavelength SED modeling that explicitly separates AGN and host contributions before applying the IR–X-ray prior to update the NH posteriors. While the original manuscript did not include a dedicated sensitivity analysis, we have added a new subsection (Section 3.4) that quantifies the impact of plausible variations: (i) ±0.4 dex scatter in the IR–X-ray relation, (ii) different torus templates, and (iii) increased star-formation contamination. The resulting CTK fraction remains within 14–20 % across these tests, confirming that the reduction from ten to three candidates and the final 17 % value are robust within the reported uncertainties. We have also clarified in the text that the prior is applied uniformly precisely because it provides a conservative, independent constraint. revision: yes

  2. Referee: [Discussion / evolution test] The claim of “no statistically significant evolution” in the CTK fraction up to z ≈ 6 (abstract and discussion) is not accompanied by the explicit statistical test, posterior comparison, or completeness-corrected number-density uncertainties that would allow the reader to assess whether the 17% value is distinguishable from lower-z measurements within the reported errors.

    Authors: We agree that an explicit statistical comparison strengthens the no-evolution statement. In the revised manuscript we have added: (i) a new panel in Figure 8 showing the CTK fraction versus redshift with completeness-corrected number-density uncertainties derived from the XLF posterior; (ii) a direct Bayesian comparison of the CTK-fraction posterior in the 3 < z < 6 bin against the lower-redshift literature values, yielding P(consistent with no evolution) = 0.68; and (iii) a Kolmogorov–Smirnov test on the binned fractions (p = 0.42). These additions are now reported in Section 5.2 and allow the reader to evaluate the claim quantitatively. revision: yes

Circularity Check

0 steps flagged

No circularity; observational update with independent IR priors

full rationale

The derivation computes the CTK fraction (17% at 3<z<6) after updating NH posteriors for the full sample using 24μm-derived IR luminosities as priors. This step is an external data incorporation, not a self-definitional loop, fitted-input prediction, or self-citation chain. The initial candidate selection (10 sources) and reduction to 3 are data-driven; applying the same prior to the remaining 801 sources does not equate the output fraction to the input X-ray fit by construction. No equations or citations reduce the central claim to its own inputs. The analysis remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the assumption that 24-micron photometry supplies reliable, unbiased priors for NH and that the X-ray selected sample remains representative after the prior update; no free parameters are explicitly fitted beyond the standard spectral model parameters.

axioms (2)
  • domain assumption X-ray spectra of AGN can be adequately described by an absorbed power-law continuum plus reflection component for the purpose of NH estimation.
    Invoked when performing the initial spectral fitting to identify CTK candidates.
  • domain assumption Infrared luminosity at 24 microns traces the intrinsic AGN power without significant host-galaxy contamination at these redshifts.
    Required to convert 24-micron flux into an IR prior for the X-ray NH posterior.

pith-pipeline@v0.9.1-grok · 5901 in / 1627 out tokens · 32410 ms · 2026-06-27T09:29:14.039501+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

117 extracted references · 20 canonical work pages · 2 internal anchors

  1. [1]

    Astropy: A Community Python Package for Astronomy

    , title =. 2013 , month = oct, pages =. doi:10.1051/0004-6361/201322068 , eid =. arXiv , bdsk-url-1 =:1307.6212 , groups =

  2. [2]

    2007 , month = feb, pages =

    , title =. 2007 , month = feb, pages =. doi:10.1051/0004-6361:20066334 , eprint =

  3. [3]

    2005 , month = oct, pages =

    , title =. 2005 , month = oct, pages =. doi:10.1051/0004-6361:20042134 , eprint =

  4. [4]

    2009 , month = feb, pages =

    , title =. 2009 , month = feb, pages =. arXiv , doi =:0803.0547 , groups =

  5. [5]

    2008 , month = sep, pages =

    , title =. 2008 , month = sep, pages =. arXiv , doi =:0806.0511 , groups =

  6. [6]

    2008 , month = sep, pages =

    , title =. 2008 , month = sep, pages =. arXiv , doi =:0806.0512 , groups =

  7. [7]

    2012 , month = jun, pages =

    , title =. 2012 , month = jun, pages =. arXiv , doi =:1204.4359 , groups =

  8. [8]

    2003 , month = dec, pages =

    , title =. 2003 , month = dec, pages =. doi:10.1086/378940 , eprint =

  9. [9]

    2015 , month = aug, number =

    , title =. 2015 , month = aug, number =. arXiv , doi =:1503.01120 , groups =

  10. [10]

    2015 , month = dec, pages =

    , title =. 2015 , month = dec, pages =. arXiv , doi =:1511.04184 , groups =

  11. [11]

    2021 , month = nov, pages =

    , title =. 2021 , month = nov, pages =. arXiv , doi =:2108.11337 , groups =

  12. [12]

    2012 , month =

    , title =. 2012 , month =. arXiv , doi =:1209.5398 , groups =

  13. [13]

    2012 , month = jun, number =

    , title =. 2012 , month = jun, number =. arXiv , doi =:1112.1949 , groups =

  14. [14]

    2008 , month =

    , title =. 2008 , month =. arXiv , doi =:0803.0636 , groups =

  15. [15]

    2018 , month = sep, number =

    , title =. 2018 , month = sep, number =. arXiv , doi =:1801.02634 , groups =

  16. [16]

    2005 , month = feb, number =

    , title =. 2005 , month = feb, number =. arXiv , doi =:astro-ph/0410527 , groups =

  17. [17]

    2015 , month = jan, pages =

    , title =. 2015 , month = jan, pages =. arXiv , doi =:1501.01982 , groups =

  18. [18]

    2014 , month = sep, pages =

    , title =. 2014 , month = sep, pages =. arXiv , doi =:1406.4502 , groups =

  19. [19]

    2003 , month =

    , title =. 2003 , month =. arXiv , doi =:astro-ph/0309134 , groups =

  20. [20]

    2021 , month = apr, number =

    The Journal of Open Source Software , title =. 2021 , month = apr, number =. arXiv , doi =:2101.09604 , groups =

  21. [21]

    2016 , month = jan, number =

    Statistics and Computing , title =. 2016 , month = jan, number =. arXiv , doi =:1407.5459 , groups =

  22. [22]

    2017 , month = mar, number =

    , title =. 2017 , month = mar, number =. arXiv , doi =:1610.09380 , groups =

  23. [23]

    2014 , month = apr, pages =

    , title =. 2014 , month = apr, pages =. arXiv , doi =:1402.0004 , groups =

  24. [24]

    2019 , month = sep, pages =

    , title =. 2019 , month = sep, pages =. arXiv , doi =:1907.13137 , groups =

  25. [25]

    2015 , month = apr, number =

    , title =. 2015 , month = apr, number =. arXiv , doi =:1501.02805 , groups =

  26. [26]

    2011 , month = feb, pages =

    , title =. 2011 , month = feb, pages =. arXiv , doi =:1012.0302 , groups =

  27. [27]

    2000 , month = aug, number =

    , title =. 2000 , month = aug, number =. arXiv , doi =:astro-ph/0003128 , groups =

  28. [28]

    2016 , month = mar, number =

    , title =. 2016 , month = mar, number =. arXiv , doi =:1601.00941 , groups =

  29. [29]

    1995 , month = apr, pages =

    , title =. 1995 , month = apr, pages =. arXiv:astro-ph/9409067 , groups =

  30. [30]

    2014 , month =

    , title =. 2014 , month =. arXiv , doi =:1402.1495 , groups =

  31. [31]

    2013 , month = dec, pages =

    The Messenger , title =. 2013 , month = dec, pages =

  32. [32]

    2016 , month = mar, pages =

    , title =. 2016 , month = mar, pages =. arXiv , doi =:1601.06002 , groups =

  33. [33]

    2015 , month = oct, number =

    , title =. 2015 , month = oct, number =. arXiv , doi =:1507.07558 , groups =

  34. [34]

    2017 , month = aug, number =

    , title =. 2017 , month = aug, number =. arXiv , doi =:1704.08296 , groups =

  35. [35]

    2019 , month = jan, pages =

    , title =. 2019 , month = jan, pages =. arXiv , doi =:1809.03747 , groups =

  36. [36]

    2008 , month = jun, number =

    , title =. 2008 , month = jun, number =. arXiv , doi =:0804.4007 , groups =

  37. [37]

    2011 , month = oct, pages =

    , title =. 2011 , month = oct, pages =. arXiv , doi =:1107.3777 , groups =

  38. [38]

    2013 , month = jul, pages =

    , title =. 2013 , month = jul, pages =. doi:10.1051/0004-6361/201220828 , eid =

  39. [39]

    2004 , month = jan, pages =

    , title =. 2004 , month = jan, pages =. doi:10.1086/379232 , eprint =

  40. [40]

    2022 , month = oct, pages =

    , title =. 2022 , month = oct, pages =. arXiv , doi =:2206.03508 , groups =

  41. [41]

    2014 , month = feb, pages =

    , title =. 2014 , month = feb, pages =. arXiv , doi =:1312.1248 , groups =

  42. [42]

    2013 , month = jun, pages =

    , title =. 2013 , month = jun, pages =. arXiv , doi =:1301.7307 , groups =

  43. [43]

    2007 , month = jan, number =

    Computing in Science and Engineering , title =. 2007 , month = jan, number =. doi:10.1109/MCSE.2007.55 , groups =

  44. [44]

    2013 , month = jan, number =

    , title =. 2013 , month = jan, number =. arXiv , doi =:1206.4263 , groups =

  45. [45]

    2023 , month = jan, number =

    , title =. 2023 , month = jan, number =. arXiv , doi =:2211.03800 , groups =

  46. [46]

    2015 , month = jan, pages =

    , title =. 2015 , month = jan, pages =. arXiv , doi =:1409.1867 , groups =

  47. [47]

    2007 , month = aug, number =

    , title =. 2007 , month = aug, number =. arXiv , doi =:astro-ph/0604426 , groups =

  48. [48]

    2004 , editor =

    Bayesian Inference and Maximum Entropy Methods in Science and Engineering , title =. 2004 , editor =. doi:10.1063/1.1835214 , eprint =

  49. [49]

    2017 , month = jan, pages =

    , title =. 2017 , month = jan, pages =. arXiv , doi =:1611.03501 , groups =

  50. [50]

    1969 , month = aug, number =

    , title =. 1969 , month = aug, number =. doi:10.1038/223690a0 , groups =

  51. [51]

    1983 , month = mar, pages =

    , title =. 1983 , month = mar, pages =. doi:10.1086/183981 , groups =

  52. [52]

    1984 , month = aug, pages =

    , title =. 1984 , month = aug, pages =. doi:10.1086/162331 , groups =

  53. [53]

    2016 , month = jan, number =

    , title =. 2016 , month = jan, number =. arXiv , doi =:1512.01105 , groups =

  54. [54]

    2018 , month = mar, number =

    , title =. 2018 , month = mar, number =. arXiv , doi =:1801.01881 , groups =

  55. [55]

    2015 , month = may, number =

    , title =. 2015 , month = may, number =. arXiv , doi =:1501.04335 , groups =

  56. [56]

    2013 , month =

    The Messenger , title =. 2013 , month =

  57. [57]

    2000 , month = jan, pages =

    , title =. 2000 , month = jan, pages =

  58. [58]

    2015 , month = may, pages =

    , title =. 2015 , month = may, pages =. arXiv , doi =:1503.00056 , groups =

  59. [59]

    2018 , month = jan, pages =

    , title =. 2018 , month = jan, pages =. doi:10.1093/pasj/psx063 , eid =

  60. [60]

    2015 , month = aug, pages =

    , title =. 2015 , month = aug, pages =. arXiv , doi =:1505.00811 , groups =

  61. [61]

    2023 , month = feb, number =

    , title =. 2023 , month = feb, number =. arXiv , doi =:2210.08030 , groups =

  62. [62]

    2016 , month = jun, note =

    , title =. 2016 , month = jun, note =. arXiv , doi =:1512.04317 , groups =

  63. [63]

    2022 , month = nov, pages =

    , title =. 2022 , month = nov, pages =. arXiv , doi =:2206.13935 , groups =

  64. [64]

    title Compton-thick Accretion in the Local Universe

    , title =. 2015 , month = dec, pages =. doi:10.1088/2041-8205/815/1/L13 , eid =

  65. [65]

    2017 , month =

    , title =. 2017 , month =. arXiv , doi =:1709.03989 , groups =

  66. [66]

    1968 , month = feb, pages =

    , title =. 1968 , month = feb, pages =. doi:10.1086/149446 , groups =

  67. [67]

    2017 , month = mar, number =

    , title =. 2017 , month = mar, number =. arXiv , doi =:1702.04729 , groups =

  68. [68]

    2021 , month = oct, number =

    , title =. 2021 , month = oct, number =. arXiv , doi =:2105.05659 , groups =

  69. [69]

    2023 , month = aug, pages =

    , title =. 2023 , month = aug, pages =. arXiv , doi =:2305.13368 , groups =

  70. [70]

    2012 , month = mar, number =

    , title =. 2012 , month = mar, number =. arXiv , doi =:1109.1286 , groups =

  71. [71]

    2016 , month = may, number =

    , title =. 2016 , month = may, number =. arXiv , doi =:1602.06954 , groups =

  72. [72]

    2015 , month = jul, pages =

    , title =. 2015 , month = jul, pages =. arXiv , doi =:1506.00162 , groups =

  73. [73]

    2016 , month = oct, number =

    , title =. 2016 , month = oct, number =. arXiv , doi =:1608.07579 , groups =

  74. [74]

    2014 , month = may, number =

    , title =. 2014 , month = may, number =. arXiv , doi =:1402.1836 , groups =

  75. [75]

    2022 , month = dec, number =

    , title =. 2022 , month = dec, number =. arXiv , doi =:2209.07797 , groups =

  76. [76]

    2018 , month = jan, number =

    , title =. 2018 , month = jan, number =. arXiv , doi =:1709.07892 , groups =

  77. [77]

    2016 , month = jun, pages =

    , title =. 2016 , month = jun, pages =. arXiv , doi =:1602.06299 , groups =

  78. [78]

    2019 , month =

    , title =. 2019 , month =. arXiv , doi =:1811.03094 , groups =

  79. [79]

    1991 , month = mar, pages =

    , title =. 1991 , month = mar, pages =. doi:10.1093/mnras/249.2.352 , groups =

  80. [80]

    title A two-phase model for the X-ray emission from Seyfert galaxies

    , title =. 1991 , month = oct, pages =. doi:10.1086/186171 , groups =

Showing first 80 references.