The warm outer layer of a Little Red Dot as the source of [Fe II] and collisional Balmer lines with scattering wings

TL;DR

This study uses deep JWST/NIRSpec observations of GN-9771 at $z=5.535$ to reveal a dense, warm gas envelope that shapes the optical/UV spectrum of a Little Red Dot. Through a combination of high-resolution line measurements and Cloudy-based photoionization modeling (BH*), the authors show that $n_{\rm H}\sim10^{9-10}$ cm$^{-3}$ and $T_{\rm e}\sim6000$–$7500$ K in the envelope reproduce a forest of [Fe II] lines and explain the unusually strong Balmer break and Balmer-line profiles via collisional excitation and electron scattering. A low-mass host galaxy is inferred from narrow lines such as [O III], with a star formation rate of ~5 $M_\odot$ yr$^{-1}$, while the central BH mass inferred from traditional virial methods would be biased if envelope effects are ignored. The results imply that Balmer lines in LRDs trace envelope physics rather than the SMBH vicinity and suggest LRDs represent a phase of rapid BH growth in low-mass galaxies, potentially powered by super-Eddington accretion in an extended envelope. Overall, the work provides a unified physical picture for the spectral peculiarities of LRDs and their role in early black hole and galaxy co-evolution.

Abstract

The population of the Little Red Dots (LRDs) may represent a key phase of supermassive black hole (SMBH) growth. A cocoon of dense excited gas is emerging as key component to explain the most striking properties of LRDs, such as strong Balmer breaks and Balmer absorption, as well as the weak IR emission. To dissect the structure of LRDs, we analyze new deep JWST/NIRSpec PRISM and G395H spectra of FRESCO-GN-9771, one of the most luminous known LRDs at $z=5.5$. These reveal a strong Balmer break, broad Balmer lines and very narrow [O III] emission. We unveil a forest of optical [Fe II] lines, which we argue is emerging from a dense ($n_{\rm H}=10^{9-10}$ cm$^{-3}$) warm layer with electron temperature $T_{\rm e}\approx7000$ K. The broad wings of H$α$ and H$β$ have an exponential profile due to electron scattering in this same layer. The high $\rm Hα:Hβ:Hγ$ flux ratio of $\approx10.4:1:0.14$ is an indicator of collisional excitation and resonant scattering dominating the Balmer line emission. A narrow H$γ$ component, unseen in the other two Balmer lines due to outshining by the broad components, could trace the ISM of a normal host galaxy with a star formation rate $\sim5$ M$_{\odot}$ yr$^{-1}$. The warm layer is mostly opaque to Balmer transitions, producing a characteristic P-Cygni profile in the line centers suggesting outflowing motions. This same layer is responsible for shaping the Balmer break. The broad-band spectrum can be reasonably matched by a simple photoionized slab model that dominates the $λ>1500$ Å continuum and a low mass ($\sim10^8$ M$_{\odot}$) galaxy that could explain the narrow [O III], with only subdominant contribution to the UV continuum. Our findings indicate that Balmer lines are not directly tracing gas kinematics near the SMBH and that the BH mass scale is likely much lower than virial indicators suggest.

Figures (20)

• Figure 1: Broad-band spectrum of GN-9771 in comparison to other LRDs. We have highlighted the main emission lines on the top axis. We compare to the PRISM spectrum of the other LRDs in our IFU program (Matthee et al., in prep) and of A2744-45924 labbe2024, after normalizing their flux to match GN-9771 in the 5400--5700 Å band. The spectrum of GN-9771 is very similar to the spectrum of A2744-45924, with only a slightly weaker Balmer break in the case of GN-9771.
• Figure 2: Triangular P-Cygni H$\pmb\alpha$ spectrum of GN-9771. The continuum-subtracted H$\alpha$ spectrum based on the G395H data is shown in black, whereas the red line shows the best-fit combined model ($\rm BIC = 932$). Residuals to the model are shown in the middle panel. The fiducial fitting model is described in Sect. \ref{['sec:balmer_O3']}. We indicate the masked regions based on the locations of possible narrow [Feii] emission in veron-cetty2004 (teal) and Hei $\lambda$6680 (purple). The [Nii] component is not shown due to its relative flux being negligible. The bottom panel shows the H$\alpha$ spectrum and best-fit after subtracting the exponential component to highlight the P-Cygni profile.
• Figure 3: The H$\pmb\beta$ spectrum of GN-9771. A similar model set-up was used as for H$\alpha$ described in Fig. \ref{['fig:Ha_line_fit']} ($\rm BIC = 74$). Teal regions highlight that are masked due to possible Feii emission. The orange region indicates the masked [Oiii] wavelengths. The bottom panel shows the H$\beta$ spectrum and best fit (red) after subtracting the exponential component to highlight the P-Cygni profile. Figure \ref{['fig:Hb_fit_fixed_exp']} shows a version of the same fit with the exponential scale fixed to that fitted for H$\alpha$.
• Figure 4: The [Feii] forest detected in the G395H data over the $\bf\sim5000$ Å region of GN-9771. We show the fitted [Feii] spectrum as described in Sect. \ref{['sec:forb_feii_model']} (green), the Hei lines (pink), H$\gamma$ (blue dot-dashed), the H$\beta$ model (blue dashed), [Oiii] $\lambda\lambda 4960,5008$ (orange dot-dashed), and [Oiii] $\lambda4364$ (brown). The residuals of the fit as well as the observational uncertainties on the spectrum are shown in the bottom panel.
• Figure 5: Same as Fig. \ref{['fig:feii_fit_1']}, zooming-in the wavelength range containing H$\gamma$ and [Oiii] $\lambda 4364$.