Black Hole Mass Estimates Research Articles

The SDSS-V Black Hole Mapper Reverberation Mapping Project: A Kinematically Variable Broad-line Region and Consequences for the Masses of Luminous Quasars

We present a velocity-resolved reverberation mapping analysis of the hypervariable quasar RM160 (SDSS J141041.25+531849.0) at z = 0.359 with 153 spectroscopic epochs of data representing a 10 yr baseline (2013–2023). We split the baseline into two regimes based on the 3× flux increase in the light curve: a “low state” phase during the years 2013–2019 and a “high state” phase during the years 2022–2023. The velocity-resolved lag profiles (VRLPs) indicate that gas with different kinematics dominates the line emission in different states. The Hβ VRLP begins with a signature of inflow onto the broad-line region (BLR) in the low state, while in the high state it is flatter with less signature of inflow. The Hα VRLP begins consistent with a virialized BLR in the low state, while in the high state shows a signature of inflow. The differences in the kinematics between the Balmer lines and between the low state and the high state suggests complex BLR dynamics. We find that the BLR radius and velocity (both FWHM and σ) do not obey a constant virial product throughout the monitoring period. We find that the BLR lags and continuum luminosity are correlated, consistent with rapid response of the BLR gas to the illuminating continuum. The BLR kinematic profile changes in unpredictable ways that are not related to continuum changes and reverberation lag. Our observations indicate that nonvirial kinematics can significantly contribute to observed line profiles, suggesting caution for black hole mass estimation in luminous and highly varying quasars like RM160.

Stripe 82X Data Release 3: Multiwavelength Catalog with New Spectroscopic Redshifts and Black Hole Masses

We present the third catalog release of the wide-area (31.3 deg2) Stripe 82 X-ray survey. This catalog combines previously published X-ray source properties with multiwavelength counterparts and photometric redshifts, presents 343 new spectroscopic redshifts, and provides black hole masses for 1297 Type 1 active galactic nuclei (AGN). With spectroscopic redshifts for 3457 out of 6181 Stripe 82X sources, the survey has a spectroscopic completeness of 56%. This completeness rises to 90% when considering the contiguous portions of the Stripe 82X survey with homogeneous X-ray coverage at an optical magnitude limit of r < 22. Within that portion of the survey, 23% of AGN can be considered obscured by being either a Type 2 AGN, reddened (R − K > 4, Vega), or X-ray obscured with a column density of N H > 1022 cm−2. Unlike other surveys, there is only an 18% overlap between Type 2 and X-ray obscured AGN. We calculated black hole masses for Type 1 AGN that have Sloan Digital Sky Survey spectra using virial mass estimators calibrated on the Hβ, Mg ii, Hα, and C iv emission lines. We find wide scatter in these black hole mass estimates, indicating that statistical analyses should use black hole masses calculated from the same formula to minimize bias. We find that the AGN with the highest X-ray luminosities are accreting at the highest Eddington ratios, consistent with the picture that most black hole mass accretion happens in the phase when the AGN is luminous (L 2−10keV > 1045 erg s−1).

Red Type-1 Quasars after Cosmic Noon and Impact on L UV-related Quasar Statistics

Over the past decades, nearly a million quasars have been explored to shed light on the evolution of supermassive black holes and galaxies. The ultraviolet-to-optical spectra of type-1 quasars particularly offer insights into their black hole activities. Recent findings, however, raise questions about the prevalence of red type-1 quasars of which colors might be due to dust obscuration and their potential influence on luminosity-related properties of quasars. We examine the fraction of red type-1 quasars within the redshift range of 0.68 ≤ z < 2.20, applying a spectral energy distribution (SED) fitting using optical-to-mid-infrared (MIR) photometric data of Sloan Digital Sky Survey Data Release 14 quasars. Approximately 10% of the type-1 quasars exhibit red colors suggestive of dust obscuration. There is an association between the brightness of the MIR luminosity and a higher fraction of red type-1 quasars, albeit with negligible redshift evolution. By employing E(B − V) values from the SED fitting, we obtained dereddened luminosity of the red type-1 quasars and reassess the quasar luminosity function (QLF) and black hole mass (M BH) estimates. Result shows a modest increase in the number density of bright quasars, linking to more flatten bright-end slope of QLFs, while M BH adjustments are minimal. Current SDSS selections with optical colors could miss a significant population of heavily dust-obscured quasars. As future MIR surveys like Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer expand, they may reveal enough obscured quasars to prompt a more profound revision of fundamental quasar properties.

DAVOS: Dwarf Active Galactic Nuclei from Variability for the Origins of Seeds: Properties of Variability-selected Active Galactic Nuclei in the COSMOS Field and Expectations for the Rubin Observatory

We study the black hole mass–host galaxy stellar mass relation, M BH–M *, of a sample of z < 4 optically variable active galactic nuclei (AGNs) in the COSMOS field. The parent sample of 491 COSMOS AGNs were identified by optical variability from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) program. Using publicly available catalogs and spectra, we consolidate their spectroscopic redshifts and estimate virial black hole masses using broad-line widths and luminosities. We show that variability searches with deep, high-precision photometry like the HSC-SSP can identity AGNs in low-mass galaxies up to z ∼ 1. However, their black holes are more massive given their host galaxy stellar masses than predicted by the local relation for active galaxies. We report that z ∼ 0.5–4 variability-selected AGNs are meanwhile more consistent with the M BH–M * relation for local inactive early-type galaxies. This result is in agreement with most previous studies of the M BH–M * relation at similar redshifts and indicates that AGNs selected from variability are not intrinsically different from the broad-line Type 1 AGN population at similar luminosities. Our results demonstrate the need for robust black hole and stellar mass estimates for intermediate-mass black hole candidates in low-mass galaxies at similar redshifts to anchor this scaling relation. Assuming that these results do not reflect a selection bias, they appear to be consistent with self-regulated feedback models wherein the central black hole and stars in galaxies grow in tandem.

Quasi-stars as a Means of Rapid Black Hole Growth in the Early Universe

JWST observations demonstrate that supermassive black holes (SMBHs) exist by redshifts z ≳ 10, providing further evidence for “direct collapse” black hole (BH) formation, whereby massive (∼103–5 M ⊙) SMBH seeds are generated within a few million years as a byproduct of the rapid inflow of gas into the centers of protogalaxies. Here we analyze the intermediate “quasi-star” phase that accompanies some direct-collapse models, during which a natal BH accretes mass from and energetically sustains (through accretion) an overlying gaseous envelope. We argue that previous estimates of the maximum BH mass that can be reached during this stage, ∼1% of the total quasi-star mass, are unphysical, and arise from underestimating the efficiency with which energy can be transported outward from regions close to the BH. We construct new quasi-star models that consist of an inner, “saturated convection” region (which conforms to a convection-dominated accretion flow near the BH) matched to an outer, adiabatic envelope. These solutions exist up to a BH mass of ∼60% of the total quasi-star mass, at which point the adiabatic envelope contains only 2% of the mass (with the remaining ∼38% in the saturated-convection region), and this upper limit is reached within a time of 20–40 Myr. We conclude that quasi-stars remain a viable route for producing SMBHs at large redshifts, which is consistent with recent JWST observations.

The black hole masses of high-redshift QSOs

ABSTRACT Observations of high-redshift quasars frequently promote suggestions of large black hole masses, whose presence so early in cosmic time is not easily explicable. I consider the parallel with ultraluminous X-ray sources (ULXs) – now known to be stellar-mass black hole (and neutron star) binaries apparently radiating far above their Eddington luminosities LEdd. The true luminosity in ULXs is actually only of order LEdd, for stellar-mass accretors, but has a very anisotropic (‘beamed’) component, plus a near-isotropic component of similar luminosity but much lower specific intensity. Observers viewing ULXs from within the beam but assuming spherical symmetry deduce a luminosity ≫LEdd. These features appear because the accretors are fed mass at highly super-Eddington rates, most of it expelled in high-speed (v &amp;gt; 0.2c) outflows from the accretion disc. I show that in similarly beamed AGN, emission-line properties would be essentially the same as in unbeamed sources, but standard virial mass indicators unusable because velocity widths are dominated by the outflows, not bound motions about the black holes. In an ensemble of this kind the apparently most luminous systems are always the most distant, but have the lowest black hole masses. Interpreting observations of this ensemble without knowing that they are beamed leads instead to very high black hole mass estimates. The analogy with ULXs therefore suggests that high-redshift quasars might actually have central black hole masses which could have grown from stellar values within the lookback time. I consider how one might test these ideas observationally.

Exploring Low-mass Black Holes through Tidal Disruption Events in the Early Universe: Perspectives in the Era of the JWST, Roman Space Telescope, and LSST Surveys

The James Webb Space Telescope (JWST) has recently uncovered the presence of low-luminosity active galactic nuclei (AGNs) at z = 4–11. Spectroscopic observations have provided estimates of the nuclear black hole (BH) masses for these sources, extending the low-mass boundary down to M • ∼ 106–7 M ⊙. Despite this breakthrough, the observed lowest mass of BHs is still ≳1–2 orders of magnitude heavier than the predicted mass range of their seed population, thereby leaving the initial mass distribution of massive BHs poorly constrained. In this paper, we focus on UV-to-optical (in the rest frame) flares of stellar tidal disruption events (TDEs) embedded in low-luminosity AGNs as a tool for exploring low-mass BH populations with ≲104–6 M ⊙. We provide an estimate of the TDE rate over z = 4–11, associated with the properties of JWST-detected AGN host galaxies, and we find that deep and wide survey programs with JWST and the Roman Space Telescope (RST) can detect and identify TDEs up to z ≃ 4–7. The predicted detection numbers of TDEs at z > 4 in 1 yr are NTDE∼2–10(0.2–2) for the JADES-Medium (and COSMOS-Web) survey with JWST and NTDE∼2–10(8–50) for the deep (and wide) tiers of the High Latitude Time Domain Survey with RST. We further discuss survey strategies for hunting for transient high-redshift TDEs in wide-field surveys with RST, as well as a joint observation campaign with the Vera C. Rubin Observatory for enhancing the detection number. The high-redshift TDE search will give us a unique opportunity to probe the mass distribution of early BH populations.

The Seoul National University AGN Monitoring Project. III. Hβ Lag Measurements of 32 Luminous Active Galactic Nuclei and the High-luminosity End of the Size–Luminosity Relation

We present the main results from a long-term reverberation mapping campaign carried out for the Seoul National University AGN Monitoring Project (SAMP). High-quality data were obtained during 2015–2021 for 32 luminous active galactic nuclei (AGNs; i.e., continuum luminosity in the range of 1044–46 erg s−1) at a regular cadence, of 20–30 days for spectroscopy and 3–5 days for photometry. We obtain time lag measurements between the variability in the Hβ emission and the continuum for 32 AGNs; 25 of those have the best lag measurements based on our quality assessment, examining correlation strength and the posterior lag distribution. Our study significantly increases the current sample of reverberation-mapped AGNs, particularly at the moderate-to-high-luminosity end. Combining our results with literature measurements, we derive an Hβ broadline region size–luminosity relation with a shallower slope than reported in the literature. For a given luminosity, most of our measured lags are shorter than the expectations, implying that single-epoch black hole mass estimators based on previous calibrations could suffer large systematic uncertainties.

Subaru High-z Exploration of Low-luminosity Quasars (SHELLQs). XVII. Black Hole Mass Distribution at z ∼ 6 Estimated via Spectral Comparison with Low-z Quasars

We report the distribution of black hole (BH) masses and Eddingont ratios estimated for a sample of 131 low luminosity quasars in the early cosmic epoch (5.6 < z < 7.0). Our work is based on the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which has constructed a low luminosity quasar sample down to M 1450 ∼ − 21 mag, exploiting the survey data of Hyper Suprime-Cam installed on Subaru Telescope. The discovery spectra of these quasars are limited to the rest-frame wavelengths of ∼1200–1400 Å, which contain no emission lines that can be used as BH mass estimators. In order to overcome this problem, we made use of low-z counterpart spectra from the Sloan Digital Sky Survey, which are spectrally matched to the high-z spectra in overlapping wavelengths. We then combined the C iv emission line widths of the counterparts with the continuum luminosity from the SHELLQs data to estimate BH masses. The resulting BH mass distribution has a range of ∼107–10 M ⊙, with most of the quasars having BH masses ∼108 M ⊙ with sub-Eddington accretion. The present study provides not only a new insight into normal quasars in the reionization epoch, but also a new promising way to estimate BH masses of high-z quasars without near-infrared spectroscopy.

Spectroscopic and Photometric Study of the Low-luminous AGN NGC 4395

The study of intermediate mass black holes (IMBHs) and estimating the black hole mass could potentially elucidate the origin of supermassive black hole (SMBH) seeds present at high redshifts and how they co-evolve with their host galaxy. The difficulty in detecting these IMBHs and the requirement of high spatial resolution to measure BHs mass makes this task challenging. In this work, we have performed photometric and spectroscopic analyses of an extremely low luminosity Seyfert 1 galaxy NGC4395. Our observations revealed strong emission lines in the spectra, and we observed a fractional variability in the V-band of 12%. Utilizing the broad-line region (BLR) size-luminosity relation for active galactic nuclei (AGNs), we estimated the size of the BLR to be 74±14 light minutes, based on a continuum luminosity of 1038 erg s−1 at 5100 Å. Additionally, the gas clouds are found to be revolving with a velocity dispersion of 384 km s−1 around the central engine, leading to a black hole mass estimation of ∼3.87−1.1+1.1 ×105 M⊙.

BASS. XXXV. The M BH–σ* Relation of 105 Month Swift-BAT Type 1 AGNs

We present two independent measurements of stellar velocity dispersions (σ ⋆) from the Ca ii H+K λ3969, 3934 and Mg i b λ 5183, 5172, 5167 region (3880–5550 Å) and the calcium triplet region (8350–8750 Å) for 173 hard X-ray-selected Type 1 active galactic nuclei (AGNs; z ≤ 0.08) from the 105 month Swift-BAT catalog. We construct one of the largest samples of local Type 1 AGNs that have both single-epoch virial black hole mass (M BH) estimates and σ ⋆ measurements obtained from high spectral resolution data, allowing us to test the usage of such methods for supermassive black hole studies. We find that the two independent σ ⋆ measurements are highly consistent with each other, with an average offset of only 0.002 ± 0.001 dex. Comparing M BH estimates based on broad emission lines and stellar velocity dispersion measurements, we find that the former is systematically lower by ≈0.12 dex. Consequently, Eddington ratios estimated through broad-line M BH determinations are similarly biased (but in the opposite way). We argue that the discrepancy is driven by extinction in the broad-line region. We also find an anticorrelation between the offset from the M BH–σ ⋆ relation and the Eddington ratio. Our sample of Type 1 AGNs shows a shallower M BH–σ ⋆ relation (with a power-law exponent of ≈3.5) compared with that of inactive galaxies (with a power-law exponent of ≈4.5), confirming earlier results obtained from smaller samples.

Estimators of Bolometric Luminosity and Black Hole Mass with Mid-infrared Continuum Luminosities for Dust-obscured Quasars: Prevalence of Dust-obscured SDSS Quasars

We present bolometric luminosity (L bol) and black hole (BH) mass (M BH) estimators based on mid-infrared (MIR) continuum luminosity (hereafter, L MIR) that are measured from infrared (IR) photometric data. The L MIR-based estimators are relatively immune from dust extinction effects, hence they can be used for dust-obscured quasars. To derive the L bol and M BH estimators, we use unobscured quasars selected from the Sloan Digital Sky Survey (SDSS) quasar catalog, which have wide ranges of L bol (1044.62–1046.16 erg s−1) and M BH (107.14–109.69 M ⊙). We find empirical relations between (i) continuum luminosity at 5100 Å (hereafter, L5100) and L MIR; (ii) L bol and L MIR. Using these relations, we derive the L MIR-based L bol and M BH estimators. We find that our estimators allow the determination of L bol and M BH at an accuracy of ∼0.2 dex against the fiducial estimates based on the optical properties of the unobscured quasars. We apply the L MIR-based estimators to SDSS quasars at z ≲ 0.5 including obscured ones. The ratios of L bol from the L MIR-based estimators to those from the optical luminosity-based estimators become larger with the amount of the dust extinction, and a non-negligible fraction (∼15%) of the SDSS quasars exhibits ratios greater than 1.5. This result suggests that dust extinction can significantly affect physical parameter derivations even for SDSS quasars, and that dust extinction needs to be carefully taken into account when deriving quasar properties.

Black Hole Accretion with Saturated Magnetic Pressure and Disk Wind

We construct an analytical black hole accretion disk model that incorporates both magnetic pressure and disk wind, which are found to be important from numerical simulations. A saturated magnetic pressure that relates the Alfvén velocity with local Keplerian velocity and gas sound speed is assumed in addition to radiation and gas pressures. The mass accretion rate is assumed to have a power-law form in response to mass loss in the wind. We find three sets of self-consistent solutions that are thermally stable and satisfy the model assumptions. At high accretion rates, the disk is geometrically and optically thick, resembling the slim disk solution. At relatively low accretion rates, our model predicts an accretion flow consisting of a geometrically thin and optically thick outer disk (similar to the standard disk), and a geometrically thick and optically thin inner disk (similar to the advection-dominated accretion flow, or ADAF). Thus, this is a natural solution for a truncated disk connected with an inner ADAF, which has been proposed to explain some observations. The magnetic pressure plays a more important role than the outflow in shaping the disk structure. The observed disk luminosity tends to saturate around 8 times the Eddington limit, suggesting that supercritical accretion onto black holes can be used for a black hole mass estimate, or a standard candle with known black hole masses.

Uncertainty quantification of the virial black hole mass with conformal prediction

ABSTRACT Precise measurements of the black hole mass are essential to gain insight on the black hole and host galaxy co-evolution. A direct measure of the black hole mass is often restricted to nearest galaxies and instead, an indirect method using the single-epoch virial black hole mass estimation is used for objects at high redshifts. However, this method is subjected to biases and uncertainties as it is reliant on the scaling relation from a small sample of local active galactic nuclei. In this study, we propose the application of conformalized quantile regression (CQR) to quantify the uncertainties of the black hole predictions in a machine learning setting. We compare CQR with various prediction interval techniques and demonstrated that CQR can provide a more useful prediction interval indicator. In contrast to baseline approaches for prediction interval estimation, we show that the CQR method provides prediction intervals that adjust to the black hole mass and its related properties. That is it yields a tighter constraint on the prediction interval (hence more certain) for a larger black hole mass, and accordingly, bright and broad spectral line width source. Using a combination of neural network model and CQR framework, the recovered virial black hole mass predictions and uncertainties are comparable to those measured from the Sloan Digital Sky Survey. The code is publicly available.

Gemini Near Infrared Spectrograph–Distant Quasar Survey: Prescriptions for Calibrating UV-based Estimates of Supermassive Black Hole Masses in High-redshift Quasars

The most reliable single-epoch supermassive black hole mass (M BH) estimates in quasars are obtained by using the velocity widths of low-ionization emission lines, typically the Hβ λ4861 line. Unfortunately, this line is redshifted out of the optical band at z ≈ 1, leaving M BH estimates to rely on proxy rest-frame ultraviolet (UV) emission lines, such as C iv λ1549 or Mg ii λ2800, which contain intrinsic challenges when measuring, resulting in uncertain M BH estimates. In this work, we aim at correcting M BH estimates derived from the C iv and Mg ii emission lines based on estimates derived from the Hβ emission line. We find that employing the equivalent width of C iv in deriving M BH estimates based on Mg ii and C iv provides values that are closest to those obtained from Hβ. We also provide prescriptions to estimate M BH values when only C iv, only Mg ii, and both C iv and Mg ii are measurable. We find that utilizing both emission lines, where available, reduces the scatter of UV-based M BH estimates by ∼15% when compared to previous studies. Lastly, we discuss the potential of our prescriptions to provide more accurate and precise estimates of M BH given a much larger sample of quasars at 3.20 ≲ z ≲ 3.50, where both Mg ii and Hβ can be measured in the same near-infrared spectrum.

Shedding New Light on Weak Emission-line Quasars in the C iv–Hβ Parameter Space

Weak emission-line quasars (WLQs) are a subset of type 1 quasars that exhibit extremely weak Lyα + N v λ1240 and/or C iv λ1549 emission lines. We investigate the relationship between emission-line properties and accretion rate for a sample of 230 “ordinary” type 1 quasars and 18 WLQs at z < 0.5 and 1.5 < z < 3.5 that have rest-frame ultraviolet and optical spectral measurements. We apply a correction to the Hβ-based black hole mass (M BH) estimates of these quasars using the strength of the optical Fe ii emission. We confirm previous findings that WLQs’ M BH values are overestimated by up to an order of magnitude using the traditional broad-emission-line region size–luminosity relation. With this M BH correction, we find a significant correlation between Hβ-based Eddington luminosity ratios and a combination of the rest-frame C iv equivalent width and C iv blueshift with respect to the systemic redshift. This correlation holds for both ordinary quasars and WLQs, which suggests that the two-dimensional C iv parameter space can serve as an indicator of accretion rate in all type 1 quasars across a wide range of spectral properties.

What Does the Geometry of the Hβ BLR Depend On?

We combine our dynamical modeling black-hole mass measurements from the Lick AGN Monitoring Project 2016 sample with measured cross-correlation time lags and line widths to recover individual scale factors, f, used in traditional reverberation-mapping analyses. We extend our sample by including prior results from Code for AGN Reverberation and Modeling of Emission Lines (caramel) studies that have utilized our methods. Aiming to improve the precision of black-hole mass estimates, as well as uncover any regularities in the behavior of the broad-line region (BLR), we search for correlations between f and other AGN/BLR parameters. We find (i) evidence for a correlation between the virial coefficient and black-hole mass, (ii) marginal evidence for a similar correlation between and black-hole mass, (iii) marginal evidence for an anticorrelation of BLR disk thickness with and , and (iv) marginal evidence for an anticorrelation of inclination angle with , , and . Last, we find marginal evidence for a correlation between line-profile shape, when using the root-mean-square spectrum, , and the virial coefficient, , and investigate how BLR properties might be related to line-profile shape using caramel models.

Using rational filters to uncover the first ringdown overtone in GW150914

There have been debates in the literature about the existence of the first overtone in the ringdown of GW150914. We develop a novel Bayesian framework to reanalyze the data of this event, by incorporating a new technique, the ``rational filter'' that can clean particular modes from the ringdown signal. We examine the existence of the first overtone in GW150914 from multiple novel perspectives. First, we confirm that the estimates of the remnant black hole mass and spin are more consistent with those obtained from the full inspiral-merger-ringdown signal when including the first overtone at an early stage of the ringdown (right after the inferred signal peak); such improvement fades away at later times. Second, we formulate a new way to compare the ringdown models with and without the first overtone by calculating the Bayes factor at different times during the ringdown. We obtain a Bayes factor of 600 at the time when the signal amplitude reaches its peak. The Bayes factor decreases sharply when moving away from the peak time and eventually oscillates around a small value when the overtone signal is expected to have decayed. Third, we clean the fundamental mode from the ringdown of GW150914 and estimate the amplitudes of the modes using the filtered data with Markov chain Monte Carlo (MCMC). The inferred amplitude of the fundamental mode is $\ensuremath{\sim}0$ whereas the amplitude of the first overtone remains almost unchanged, implying that the filtered data is consistent with a first-overtone-only template. Similarly, if we remove the first overtone from the GW150914 data, the filtered data are consistent with a fundamental-mode-only template. Finally, after removing the fundamental mode, we use MCMC to infer the remnant black hole mass and spin from the first overtone alone. We find the posteriors are still informative and consistent with those inferred from the fundamental mode. The conclusions are also verified through simulations in Gaussian noise using a GW150914-like numerical relativity waveform.

The first X-ray look at SMSS J114447.77-430859.3: the most luminous quasar in the last 9 Gyr

ABSTRACT SMSS J114447.77-430859.3 (z = 0.83) has been identified in the SkyMapper Southern Survey as the most luminous quasar in the last $\sim 9\, \rm Gyr$ . In this paper, we report on the eROSITA/Spectrum–Roentgen–Gamma (SRG) observations of the source from the eROSITA All Sky Survey, along with presenting results from recent monitoring performed using Swift, XMM-Newton, and NuSTAR. The source shows a clear variability by factors of ∼10 and ∼2.7 over time-scales of a year and of a few days, respectively. When fit with an absorbed power law plus high-energy cutoff, the X-ray spectra reveal a Γ = 2.2 ± 0.2 and $E_{\rm cut}=23^{+26}_{-5}\, \rm keV$ . Assuming Comptonization, we estimate a coronal optical depth and electron temperature of $\tau =2.5-5.3\, (5.2-8)$ and $kT=8-18\, (7.5-14)\, \rm keV$ , respectively, for a slab (spherical) geometry. The broadband SED is successfully modelled by assuming either a standard accretion disc illuminated by a central X-ray source, or a thin disc with a slim disc emissivity profile. The former model results in a black hole mass estimate of the order of $10^{10}\, \mathrm{ M}_\odot$ , slightly higher than prior optical estimates; meanwhile, the latter model suggests a lower mass. Both models suggest sub-Eddington accretion when assuming a spinning black hole, and a compact ($\sim 10\, r_{\rm g}$ ) X-ray corona. The measured intrinsic column density and the Eddington ratio strongly suggest the presence of an outflow driven by radiation pressure. This is also supported by variation of absorption by an order of magnitude over the period of $\sim 900 \ \rm d$ .

The need for obscured supermassive black hole growth to explain quasar proximity zones in the epoch of reionization

ABSTRACTProximity zones of high-redshift quasars are unique probes of supermassive black hole formation, but simultaneously explaining proximity zone sizes and black hole masses has proved to be challenging. We study the robustness of some of the assumptions that are usually made to infer quasar lifetimes from proximity zone sizes. We show that small proximity zones can be readily explained by quasars that vary in brightness with a short duty cycle of fduty ∼ 0.1 and short bright periods of ton ∼ 104 yr, even for long lifetimes. We further show that reconciling this with black hole mass estimates requires the black hole to continue to grow and accrete during its obscured phase. The consequent obscured fractions of ≳ 0.7 or higher are consistent with low-redshift measurements and models of black hole accretion. Such short duty cycles and long obscured phases are also consistent with observations of large proximity zones, thus providing a simple, unified model for proximity zones of all sizes. The large dynamic range of our simulation, and its calibration to the Lyman-α forest, allows us to investigate the influence of the large-scale topology of reionization and the quasar’s host halo mass on proximity zones. We find that incomplete reionization can impede the growth of proximity zones and make them smaller up to 30 per cent, but the quasar host halo mass only affects proximity zones weakly and indirectly. Our work suggests that high-redshift proximity zones can be an effective tool to study quasar variability and black hole growth.