Quantifying biases in stellar masses of JWST high-z quasar host galaxies caused by quasar subtraction
Pith reviewed 2026-05-19 09:13 UTC · model grok-4.3
The pith
Correcting for quasar subtraction biases leaves the mismatch between high-redshift black hole and host galaxy masses intact.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors construct corrected host magnitude posteriors from simulated JWST images and find that quasar subtraction either recovers the true host magnitude or underestimates it by at most 0.2 mag. After performing SED fitting on the corrected photometry, stellar mass estimates prove robust or change by less than 0.3 dex. The mismatch between black hole and host galaxy growth therefore persists, with some quasars showing overmassive black holes and others residing in relatively massive galaxies.
What carries the argument
Statistical construction of corrected host magnitude posteriors from BlueTides mock images to isolate biases introduced by quasar point source removal.
If this is right
- Stellar mass estimates from JWST photometry remain reliable within 0.3 dex after accounting for quasar subtraction.
- The apparent mismatch between black hole masses and host stellar masses at high redshift is not produced by these observational biases.
- A subset of hosts, including J0844-0132, J0911+0152, and J1146-0005, stay unconstrained because key photometric bands provide only upper limits.
- Accounting for the bias does not eliminate the observed diversity in host galaxy masses at early times.
Where Pith is reading between the lines
- The persistence of the mismatch implies that explanations must be sought in black hole seeding or growth physics rather than in data processing artifacts.
- The statistical correction method can be applied to additional JWST targets to produce more uniform host mass catalogs.
- If the simulation light profiles match reality, then future observations at longer wavelengths may further tighten the mass constraints.
Load-bearing premise
The correction assumes the point spread function is modeled perfectly and that the BlueTides simulation accurately represents the light profiles of real high-redshift quasar host galaxies.
What would settle it
Independent stellar mass measurements of the same host galaxies obtained without point source subtraction, for example from higher-resolution imaging or different wavelengths, that differ by more than 0.3 dex from the corrected values.
read the original abstract
JWST has enabled dozens of high-$z$ quasar host galaxy detections. Many of these observations imply galaxies with black holes that are overmassive compared to their low-$z$ counterparts. However, the bright quasar point source removal can cause significant biases in recovered host magnitudes and stellar mass measurements due to the degeneracy in host galaxy and quasar light. We develop a statistical method to disentangle the quasar host galaxy stellar mass measurements from observational biases during the point source removal assuming the PSF is modelled perfectly. We use the BlueTides simulation to generate mock images and perform point source removal on thousands of simulated high-$z$ quasar host galaxies, constructing corrected host magnitude posteriors. We find that removing a bright quasar in JWST photometry tends to either correctly recover or modestly misestimate host magnitudes, with a maximum magnitude underestimate of 0.2 mag. With our corrected magnitude posteriors, we perform SED fitting on each quasar host galaxy and compare the stellar mass measurement before and after the correction. We find that stellar mass estimates are generally robust, or misestimated by <0.3 dex. We also find that the stellar masses of a subset of hosts (J0844-0132, J0911+0152, and J1146-0005) remain unconstrained, as key photometric bands provide only flux upper limits. Accounting for observational biases does not resolve the apparent mismatch between black hole and host galaxy growth at high-$z$, where some quasars appear to host overmassive black holes while others reside in relatively massive galaxies.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims to quantify biases in stellar masses of JWST high-z quasar host galaxies due to quasar subtraction using a statistical method based on BlueTides simulation mocks. They generate thousands of mock images, perform point source removal, and derive corrected magnitude posteriors showing at most 0.2 mag underestimate. SED fitting then shows stellar mass shifts <0.3 dex, and they conclude that this does not resolve the black hole-host galaxy growth mismatch at high redshift.
Significance. The result, if valid, is significant for the interpretation of JWST high-redshift quasar observations. The simulation-driven forward modeling provides specific, testable bias values (0.2 mag, <0.3 dex) and avoids circularity by using independent mocks rather than fitting to the real data. This supports the robustness of the observed mismatch between overmassive black holes and host masses, advancing understanding of early galaxy and BH co-evolution.
major comments (2)
- [§3] §3 (PSF modeling assumption): The bias quantification assumes the PSF is modelled perfectly with no residuals or wavelength-dependent errors. No sensitivity analysis to realistic PSF imperfections is presented, which directly affects the reliability of the reported 0.2 mag maximum underestimate and <0.3 dex mass shift bounds when applied to actual JWST photometry.
- [§4.2] §4.2 (simulation validation): The central claim that accounting for biases does not resolve the mismatch depends on BlueTides faithfully reproducing the distribution of real high-z host sizes, Sérsic indices, and quasar-host flux ratios. No direct comparison of simulated structural parameters to the observed sample (e.g., J0844-0132) or other high-z simulations is provided, which is load-bearing for generalizing the correction.
minor comments (2)
- The abstract is dense; splitting the description of the method and results into shorter sentences would improve clarity for readers.
- [Figure captions] Figure captions should explicitly define the range of quasar-host flux ratios and Sérsic indices used in the mocks to aid interpretation of the posterior distributions.
Simulated Author's Rebuttal
We thank the referee for their constructive comments on our manuscript. We address the major comments point by point below, providing clarifications and indicating revisions where we agree that additional discussion strengthens the paper.
read point-by-point responses
-
Referee: [§3] §3 (PSF modeling assumption): The bias quantification assumes the PSF is modelled perfectly with no residuals or wavelength-dependent errors. No sensitivity analysis to realistic PSF imperfections is presented, which directly affects the reliability of the reported 0.2 mag maximum underestimate and <0.3 dex mass shift bounds when applied to actual JWST photometry.
Authors: We appreciate the referee pointing out this limitation. As explicitly stated in the abstract and throughout Section 3, our analysis assumes a perfectly modeled PSF to isolate the bias arising from the quasar subtraction procedure under ideal conditions. This provides a baseline quantification of the degeneracy between quasar and host light. We agree that realistic PSF imperfections could add further uncertainty in real data applications. In the revised manuscript, we will expand the discussion in Section 3 to include a qualitative assessment of how PSF residuals might affect the results and emphasize that our reported maximum bias of 0.2 mag is a lower bound under the perfect PSF assumption. No changes to the quantitative results are needed, as they are conditional on this assumption. revision: partial
-
Referee: [§4.2] §4.2 (simulation validation): The central claim that accounting for biases does not resolve the mismatch depends on BlueTides faithfully reproducing the distribution of real high-z host sizes, Sérsic indices, and quasar-host flux ratios. No direct comparison of simulated structural parameters to the observed sample (e.g., J0844-0132) or other high-z simulations is provided, which is load-bearing for generalizing the correction.
Authors: We thank the referee for this observation. The BlueTides simulation has been validated in multiple prior works for reproducing high-redshift galaxy properties, including sizes and luminosity functions. Our approach uses a large statistical sample from the simulation to perform forward modeling of the observational biases, rather than fitting to specific observed hosts. This avoids circularity with the real data. We will revise the manuscript to add references to BlueTides validation papers and clarify in Section 4.2 that while direct comparisons to individual targets like J0844-0132 are not performed (as the simulation provides a population-level distribution), the structural parameters are consistent with expectations from other high-z simulations. This bolsters the applicability of our bias corrections. revision: partial
Circularity Check
No significant circularity in forward-modeling bias correction
full rationale
The paper develops a statistical correction by generating mock JWST images from the BlueTides simulation, performing quasar subtraction on those mocks, and deriving magnitude posteriors that are then applied to real observations. This forward-modeling approach uses independent simulation inputs to quantify biases rather than fitting parameters to the target JWST data or re-deriving results from the observations themselves. No equations or steps reduce the final stellar-mass corrections or mismatch conclusion to the input data by construction, and the explicit assumption of perfect PSF modeling is stated as a boundary condition rather than a self-referential definition. The derivation remains self-contained against external simulation benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The PSF is modelled perfectly
- domain assumption BlueTides simulation accurately captures the light profiles and properties of real high-z quasar hosts
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.