A gas-enshrouded and gas-reddened black hole at cosmic dawn

Observations and information discount

MoM-BH*-1 has been noticed with JWST by three programmes. It was imaged in cycle 1 by the PRIMER survey¹⁴ (JWST-GO-1837; principal investigator: J. Dunlop) utilizing the MIRI (5 January 2023 and 16 January 2023) and NIRCam (7 August 2023 and 9 August 2023) devices. In cycle 2 (19 December 2023), the EXCELS survey¹⁵ (JWST-GO-3543) obtained 1.5 h of NIRSpec G395M spectroscopy. In cycle 3 (15 December 2024), we focused MoM-BH*-1 as a part of the ‘Mirage or Miracle’ NIRSpec prism survey (JWST-GO-5224). We included MoM-BH*-1 as a high-priority goal second in significance solely to luminous z > 10 sources in our UDS masks, as a result of it appeared in a number of precedence goal lists—AGN/LRD candidates chosen primarily based on compact morphology and template becoming with EAZY⁴⁷, extraordinarily large galaxy candidates, sources with peculiar purple colors and the literature LRD candidates⁴⁸. 

We use the v.7.2 pictures of the PRIMER subject launched on the DAWN JWST archive (DJA) diminished utilizing the grizli software program⁴⁹. PSF-matched photometric catalogues primarily based on these pictures had been produced in ref. ⁵⁰. We use the general public v.3 NIRSpec reductions of the EXCELS grating information from the DJA derived utilizing the msaexp software program^37,51,52. The MoM information are diminished with the identical pipeline following the identical decisions.

Although we discovered no related radio or ALMA archival information, MoM-BH*-1 has been noticed with Chandra⁵³. Similar to nearly all LRDs^9,54,55, it stays undetected within the X-rays (L_x < 44.5 erg s⁻¹ (1σ) at relaxation body 5–90 keV). Key empirical properties of the supply are summarized in Extended Data Table 1.

Emission line becoming

We use a customized NIRSpec emission line becoming bundle⁵⁶ to concurrently match emission strains within the grating and prism spectra. The benefit of this strategy is that regardless of the low signal-to-noise ratio (SNR) in both mode, options could also be robustly recovered due to their prevalence at the identical wavelength throughout each dispersers.

We first match the Hβ line and [Oiii] doublet, and then use the redshift as a previous to suit Hγ (Fig. 1). We mannequin Hβ as a single emission line with three absorbers constrained to have adverse flux (one at line-centre and two on both facet of zero velocity) motivated by the symmetric absorption troughs on both facet of the central double-peak (Extended Data Fig. 2). The systemic redshift is tied to [Oiii] and the broad Hβ part, with the absorbers allowed to vary freely. The variety of absorbers is set primarily based on the maxima reached within the diminished χ², which is analogous for 3 and 4 absorbers, however we go for parsimony. This giant variety of absorbers could also be merited to account for secondary peaks at ±2,400 km s⁻¹ (Extended Data Fig. 2). Including slim Hβ at the systemic redshift results in fully unconstrained flux degenerate with absorption and no enchancment, so we neglect this part. Resulting suits, during which we pattern the posterior with the NUTS sampler applied in numpyro⁵⁷, are proven in Extended Data Fig. 1 and reported in Extended Data Table 2.

Although the formal errors on the derived fluxes and line widths of lots of the parts is important, the important thing options related to this evaluation are robustly recovered—extraordinarily broad Hβ emission, and deep, broad absorption wiping out about 25% of the emission flux, together with roughly 100% of the flux at line-centre. Another notable facet of those suits is the placement of the 2 absorbers—they’re recovered at very comparable velocities, however on both facet of line-centre (roughly ±1,500 km s⁻¹), albeit with important uncertainties ((-1,53{2}_{-113}^{+345},mathrm{km},{{rm{s}}}^{-1}) and (+1,55{6}_{-1,378}^{+232},mathrm{km},{{rm{s}}}^{-1})). We discover this symmetry, which extends not solely to the placement of the absorbers but in addition to the detailed construction of the complete emission line profile over a number of 1,000 km s⁻¹ in Extended Data Fig. 2.

Morphology

We use pysersic⁵⁸ to suit a Sersic profile to the imaging. We give attention to the F356W and F444W imaging, during which the supply is well-detected. We comply with the process described in ref. ⁵⁹—we construct an empirical PSF utilizing stars within the subject, and then use this PSF with pysersic to pattern the posterior Sersic parameters with numpyro. The supply is unresolved, and we’re capable of place a 99% higher restrict on the efficient radius of <117 computer in line with the BH-dominated interpretation.

Cloudy modelling

That extraordinarily dense gasoline would possibly blanket the LRDs was already inferred within the work that outlined this class of sources⁹—about 10% of those LRDs confirmed clear indicators of Balmer absorption, with this being a decrease restrict as a result of decision and SNR¹¹ versus the ≪0.1% in pre-JWST AGN⁶⁰, with solely a handful examples⁹. Local AGN⁶¹ don’t show Balmer breaks wherever as robust as MoM-BH*-1 (about 8), with the strongest breaks (≲2.5) occurring in quiescent galaxy AGN that attain the stellar inhabitants most proven in Fig. 2b. Inspired by Balmer absorption in LRDs, and utilizing Cloudy fashions for continuum absorption⁵, demonstrated that the very gasoline producing robust Balmer absorption possible additionally produces Balmer breaks. However, observe the breaks in these fashions are weaker than what we require to elucidate MoM-BH*-1 and are abrupt as an alternative of the sleek rollover seen on this supply. Here we construct on this work.

The grid of parameters we discover is summarized in Extended Data Table 3. This grid spans a variety to seize the acute spectral form at hand. We begin with an intrinsic AGN continuum SED that’s parametrized utilizing a collection of energy legal guidelines and with a ‘big bump’ temperature³². This AGN continuum is handed by means of gasoline surrounding the central supply that’s outlined by way of its gasoline density, column density, metallicity and turbulent velocity. The irradiation of the cloud is modulated by an ionization parameter (log(U)). A giant turbulent velocity of the absorbing gasoline 500 km s⁻¹ is motivated by the width of the absorption strains we observe (see Fig. 1). This implies a excessive Mach quantity, which is in line with current fashions of AGN disks⁶²—whether or not such excessive turbulence could be sustained throughout a big envelope stays to be seen. The turbulence is necessary in producing a easy rollover as an alternative of a pointy break⁶. Finally, we apply a uniform mud display screen⁶³ with A_V = 0−3 within the post-processing—that is the one post-processing step we apply.

To choose fashions which can be in line with the info, we require: (1) Hβ in web emission with 30 < EW [Å] < 45; (2) Hγ in absorption with −5 < EW [Å] < 0 Å; (3) a robust Balmer break and excessive optical-to-UV ratio such that ({f}_{{rm{4.5mu m}}}^{lambda }/{f}_{{rm{1.8mu m}}}^{lambda } > 6) and ({f}_{{rm{4.5mu m}}}^{lambda }/{f}_{{rm{2.8mu m}}}^{lambda } > 3); and (4) fluxes in MIRI bands inside 2σ of the observations. The few thousand fashions that fulfill these constraints are re-simulated at greater decision, retaining solely hydrogen that’s related to the important thing options for simplicity and pace. Of these, we choose the one which intently follows the detailed form of the continuum. Furthermore, we experiment with the overlaying issue and distance between the gasoline and central supply, and discover these portions to be degenerate with the AGN SED and ionization parameter. We discover the ‘net transmitted’ flux (that’s, the sum of the attenuated incident continuum and diffuse continua or strains) produces a great match to the info, not the ‘total’ flux (which additionally contains mirrored continua or strains).

Key options of fiducial Cloudy mannequin: blanketed by gasoline, not by mud

First, we emphasize that this mannequin (Fig. 3 and Extended Data Fig. 3) is a mannequin that matches numerous options of curiosity primarily based on our restricted grid. The related parameter house that we’ve got achieved our greatest to pattern is high-dimensional, degenerate and a lot stays unknown (for instance, the intrinsic SEDs of early AGN with results comparable to photon trapping). Therefore, we warning towards detailed inferences past the feasibility of the broad bodily image that we current right here (an accretion disk embedded in dense gasoline).

Notably, the chosen mannequin options an excessive column density (about 10^25.8 cm⁻²) similar to essentially the most enshrouded techniques noticed⁶⁴ and a excessive gasoline density (10¹¹ cm⁻³) conducive to Balmer absorption. The turbulent velocity of 500 km s⁻¹ is commensurate with the width of the central Hβ absorber. Metal-poor gasoline anticipated of a dwarf galaxy at z ≈ 8 is most popular. The AGN slope parameters are inside the vary of literature SEDs^6,65.

An necessary characteristic of the fiducial mannequin is that it’s nearly dust-free (A_V = 0.15 magazine). The UV-weakness arises completely due to the acute hydrogen opacity. This is a vital constraint on LRD fashions that usually invoke important quantities of mud to suppress robust UV emission from classical AGN^21,24,66 or to elucidate the weak spot of Hβ relative to Hα¹². Although important A_V ≈2–3 helps dense gasoline AGN fashions produce a easy Balmer break within the rest-optical⁶, that is dominated out by longer wavelength constraints in our supply (Extended Data Fig. 3) and extra usually by stringent infrared constraints on the LRDs^34,35. We observe that given the very low A_V ≈ 0, the precise particulars of the mud geometry (for instance, display screen in contrast with clumps) wouldn’t have a bearing on our outcomes.

Resonant Balmer scattering and insights from Lyα-like shell fashions

The Hβ profile of MoM-BH*-1 bears a exceptional resemblance to double-peaked Lyα^67,68,69. This motivates us to discover whether or not Hβ is behaving like a resonant line beneath the excessive densities in MoM-BH*-1. It will be the case that this can be a widespread phenomena within the LRDs, however {that a} symmetric double-peaked line with giant separation is less complicated to look at in MoM-BH*-1 due to negligible Hβ emission from the sub-dominant host galaxy.

Although Hβ usually has a simple cascade escape route, saturating the twop state can successfully lure Hβ photons, which can be believable by means of Lyα pumping in these optically thick environments⁷⁰. We check this concept with the COLT⁷¹ radiative switch code with which we implement a easy shell mannequin (straight analogous to Lyα shell fashions) for Hβ. The dominant parameters in these calculations are the thermal velocities of the internal and outer shells, the relative velocities between the shells and the optical depth encountered by Hβ. Strong turbulence shouldn’t be accounted for in these toy fashions.

In Extended Data Fig. 4, we present how a comparatively slim Gaussian Hβ line manifests as a double-peaked profile. We verify the fundamental tendencies seen with the Lyα line maintain up for Balmer resonant scattering and can function a strong information to interpret MoM-BH*-1. With this easy toy mannequin, we’re capable of generate a match for the rate separation and depth of the strongest peaks seen within the MoM-BH*-1 Hβ profile.

If scattering is underway, this has some necessary implications for the physics of the scenario. The noticed width of Balmer strains (particularly when fitted as an absorption system to a broad Gaussian line; that’s, the commonest baseline mannequin^9,21,72) might not hint the kinematics of the broad-line area however are a consequence of resonant scattering. Therefore, SMBH plenty in these techniques which can be primarily based on line widths of the Balmer strains could also be severely overestimated by as much as 2 dex (Extended Data Table 4). More subtle remedies of resonant scattering results in LRDs⁷³ differ in assumptions and particulars, however make the identical fundamental level.

SMBH properties and ramifications for LRD parameters

Given the distinctive situations (for instance, uncommon gasoline density) on this supply, we should train care in making use of the native scaling relations usually used to derive SMBH properties. Here we discover implications for a wide range of approaches to again out the SMBH properties which can be summarized in Extended Data Table 4. Given the numerous systematic uncertainties, these calculations have to be seen as order-of-magnitude estimates to bracket the vary of potentialities.

First, we merely apply Hβ-based native scaling relations³⁹ to the absorption-corrected and dust-corrected Hβ line. The absorption correction is often achieved by emission line becoming assuming an underlying Gaussian or Lorentzian^9,21 (as in Extended Data Fig. 1). For the mud correction, observe that the Hα and Hβ line ratio¹³ is difficult to make use of provided that their noticed fluxes mirror radiative switch within the dense gasoline. The commonplace approaches used within the literature at the second infer a big A_V of some magnitudes both from continuum slope becoming or by means of SED modeling^6,10,12,13. Going by the optical continuum slope¹⁰ would indicate an A_V ≈ 2 for MoM-BH*-1 implying a BH mass of about 10^8.3M_⊙. This is similar to the stellar mass of the host galaxy (<10^8.5M_⊙ at 95% confidence), that’s, that is an apparently ‘overmassive’ BH relative to native scaling relations between host galaxy mass and BH mass^13,41.

However, a key perception from the SED of MoM-BH*-1 is that the A_V required to elucidate the LRDs is probably going negligible because the SED is intrinsically UV-weak under the Balmer break. That is, the LRDs are purple not due to obscuration from mud, however due to opacity from gasoline. Accounting for this, with A_V = 0, the SMBH mass derived now’s about 5 × 10⁷M_⊙ with a luminosity of roughly 15% the Eddington restrict.

Next, we observe that the noticed line width (after correcting for absorbers) might not be faithfully tracing the broad-line area, violating the fundamental ansatz for scaling relations (Extended Data Figs. 4 and 5). The underlying BH emission line profile in LRDs could also be as advanced because the one proven in Extended Data Fig. 2, however is maybe tough to disentangle from robust emission from the host. The advanced construction might signify radiative processes which can be modifying the emission arising from the broad-line area. As an instance of 1 such course of, if Hβ is present process resonant scattering by means of the environment across the BH, then the intrinsic BLR width earlier than scattering may very well be as little as round 600 km s⁻¹ (Extended Data Fig. 5). This would yield an SMBH mass of roughly 10⁶M_⊙.

Finally, we’re capable of constrain the BH properties utilizing our Cloudy mannequin. The bolometric luminosity straight follows from integrating the SED proven in Extended Data Fig. 3 (about 10^44.5 erg s⁻¹). We observe that that is an order of magnitude decrease than the worth anticipated from native calibrations, additional underscoring their inapplicability. We can convert this luminosity right into a mass by observing that theoretical fashions of comparable eventualities^74,75,76 predict near Eddington or super-Eddington accretion to maintain these convective envelopes. Furthermore, the ubiquity of OI emission thought to come up from Lyβ fluorescence in LRDs^21,66,77 is considered a trademark of super-Eddington accretion⁷⁸. Inspired by these strains of proof, requiring L/L_bol ≈ 1, due to this fact, permits us to estimate the BH mass as roughly 10^6.3M_⊙. Correcting the doubtless overestimated literature BH plenty by about 10−100× would deliver the everyday M_BH/M_⋆ of JWST AGN (about 1–10%; refs. ^79,80,81) nearer to the ratio within the native Universe (0.01%; ref. ⁸²).

Host galaxy properties

To constrain the host mass, we will leverage the perception from clustering analyses that the rest-UV mild within the typical LRD (FWHM = 1,000−2,000 km s⁻¹), on common, originates nearly completely from the host⁴¹. Of course, this isn’t true for all LRDs, significantly essentially the most luminous sources with greater FWHM broad strains that show AGN signatures even within the rest-UV²¹, however to date they appear the exception. To assemble an empirical M_⋆−M_UV relation, we use the compilation of low-luminosity galaxies at z = 3−7 from the All the Little Things (ALT) survey⁴⁰ within the Abell-2744 subject. This is the biggest spectroscopic pattern of M_UV < −15 galaxies at these redshifts. The ALT stellar plenty are derived with the Prospector SED becoming code utilized to 27 bands of NIRCam + HST photometry, together with all JWST medium and broad bands^83,84. We estimate {that a} supply with M_UV between −17.9 and −18.3 has a median stellar mass of (log ({M}_{star }/{M}_{odot })=7.{4}_{-0.3}^{+0.6}) and a 95% higher restrict of (log ({M}_{star }/{M}_{odot }) < 8.5).

Hints of variability

It is of explicit curiosity to check for indicators of variability in MoM-BH*-1 as an unbiased constraint on the physics of the supply. In LRDs, the SED is a summation of the host and BH at all wavelengths to various levels, however on this case, the stark outshining of the host galaxy means an ‘undiluted’ variability sign could also be stronger and extra simply detected. One of the handful experiences of variability in an LRD to this point has been in a supply with a big Balmer break implying a excessive BH fraction³⁰.

Three units of observations overlaying 3−5 μm exist for MoM-BH*-1 separated by ≈ 60 days within the rest-frame (Extended Data Fig. 6). It is actually not preferrred to check for variability throughout three totally different observing modes with distinct systematics. Nevertheless, it’s notable that the supply seems to have brightened by 30 ± 7% between the primary epoch (NIRCam) and third epoch (NIRSpec prism, no post-processing renormalization to photometry). The prism spectrum suffers from slit losses, and but this supply seems brighter than anticipated from the NIRCam photometry. It is the one supply throughout the 136 sources with high-SNR photometry (SN > 10) and spectra (50^th percentile of SN > 5) noticed in two masks as a part of the MoM program that reveals this diploma of brightening. We additionally observe that usually, the NIRSpec/G395M flux is systematically  ≈ 10–20% decrease than the prism flux as a result of calibration uncertainties⁸⁵ – accounting for this, each the prism and NIRSpec fluxes for MoM-BH*-1 seem like in settlement and  ≈ 30% brighter than the NIRCam flux. The form of the SED (and the depth of the Balmer break) is constant throughout the NIRCam and prism information (for instance, ({f}_{{rm{F410M}}}^{nu }/{f}_{{rm{F356W}}}^{nu }=2.0pm 0.2) within the prism vs. 2.2 ± 0.2 in NIRCam). The distinction is within the absolute brightness.

Taken at face worth, these measurements indicate a ≈4σ detection of variability over a mere two months, marking MoM-BH*-1 as a superb goal for future monitoring campaigns. This could also be unbiased proof not just for the AGN nature of the supply, but in addition for the unstable surroundings that exists round it.