Inside the cocoon: a comprehensive explanation of the spectra of Little Red Dots

TL;DR

This study tackles the puzzle of Little Red Dots (LRDs) in the early Universe by constructing a self-consistent, first-principles model in which a supermassive black hole accretes from a dense, ionized gas cocoon. Using the Sirocco Monte Carlo radiative-transfer framework, the authors reproduce the distinctive Balmer breaks, broad exponential hydrogen lines, and Balmer/He I absorptions observed in $z>3$ LRD-like objects, accounting for non-spherical, inflow/outflow kinematics and electron scattering. They derive physical properties of the cocoon and SMBHs, finding $oxed{\log_{10}(M/M_{ootnote})=5.7\pm0.8}$ and $oxed{\log_{10}(L/L_{ m edd})=-0.2^{+1.2}_{-0.9}}$, with cocoon masses $\sim 1$–$100\,M_{ootnote}$ and electron columns $N_e\sim10^{24}\,\text{cm}^{-2}$, while dust is negligible and a cold-gas reservoir of $\sim10^{6}\,M_{ootnote}$ may sustain accretion for $\gtrsim$Myr. The framework also predicts correlations between Balmer break strength and line widths and demonstrates that a non-spherical, balanced inflow/outflow geometry is required to match the near-symmetric scattering wings, offering new diagnostics for LRDs and implications for SMBH growth in the early Universe.

Abstract

JWST has revealed a population of compact galaxies in the early Universe with broad emission lines and strong Balmer breaks; among them the so-called ''little red dots'' (LRDs). Their nature remains uncertain with hypotheses including exotic phenomena. We assemble a sample of LRD-like objects at $z>3$ and use self-consistent radiative-transfer calculations to show that a supermassive black hole accreting from a dense gas cocoon accurately reproduces the detailed spectra. We show that the cocoons must be non-spherical, with comparable amounts of inflowing and outflowing material. And we predict correlations between Balmer break strength, Balmer line-absorption and scattering line width, which we confirm in our observed sample. We reproduce all LRD-like properties without requiring star-like atmospheres and we determine the typical black hole in our sample to be of order a million solar masses, with ionized cocoon masses of tens of solar masses potentially supplied from a much larger cold-gas reservoir.

Figures (23)

• Figure 1: Properties of the LRD sample ($3.1<z<7$). Left:JWST NIRSpec/PRISM spectra highlighting the Balmer break and the 'v-shape' SED. Middle:JWST imaging of a subset of the sample showing the compact nature of the sources (see imaging information for the full sample in App. \ref{['app:obs_sample']}). Right: Observed H$\alpha$ lines on a semi-logarithmic plot showing the exponential wings (a straight line on this plot) and the red/blue symmetry of the lines. The ratio of the exponential slopes on the blue and red sides, $v_{e,blue}/v_{e,red}$, is indicated over each line, showing a high degree of symmetry; we note a $>2\sigma$ asymmetry in the bottom object, ID 24 ('the Cliff').
• Figure 2: Overview of the physical structure of an LRD-like object, showing the gas flows and the spectral-formation regions in the ionised gas cocoon surrounding the central SMBH. Left: Inputs to the model are represented by the vector field, with arrows denoting the velocity (arrow length) and density (arrow density) of the gas. The data require both inflow and outflow (see Fig. \ref{['fig:sphericity']} and §\ref{['sec:non_spherical']}). Outputs of the model are shown in the different quadrants: colours indicate the regions producing the Balmer break (red/orange), the electron density (green), and the hydrogen-line emissivity (sum of Balmer and Paschen lines, blue). Top right: radial distributions for the emission of escaped photons in H$\alpha$ (blue), H$\beta$ (orange), and Pa$\beta$ (green). The different recombination lines are not released from identical radii. Bottom right: the local monochromatic mean intensity at various radii through the cocoon (with a greyscale contour plot showing the projection). This illustrates that recombination lines form in the inner ionised gas, whereas the Balmer break is formed further out, where a sufficiently large $n=2$ population is present.
• Figure 3: Comparison of the spectra and line shapes of observed LRD-like objects (left) and Sirocco models (right). The spectra are ordered by Balmer break strength, with the corresponding H$\alpha$ line profiles shown in the same ordering. Top left: Sequence of observed sources with increasingly prominent Balmer breaks. Top right: Spectra from a sequence of Sirocco models with varying column density in a homologous wind (keeping all other properties constant). Significant electron-scattering wings require a lower column density ($N_e\gtrsim 10^{24}\,\mathrm{cm^{-2}}$) compared to the Balmer breaks and absorption lines ($N_H\gtrsim 10^{25}\,\mathrm{cm^{-2}}$, Taylor2025). Bottom: H$\alpha$ spectra from the same observational sample (left) and Sirocco density-sequence (right) scaled to the luminosity distance at $z=4$. For both models and observations, the most luminous H$\alpha$ lines and the strongest Balmer break systems (i.e. the most massive cocoons) typically display broader wings and are more likely to have Balmer-absorption line profiles.
• Figure 4: Comparison of the H$\alpha$ broad line width and the Balmer break strength for the observed sample. The curve shows an illustrative Sirocco model density-sequence. To account for reddening, we apply a dust screen with $E(B-V)=0.8\pm0.2$; the coloured curve shows the mean sequence, while the grey curves shift the models horizontally in Balmer-break strength by varying $E(B-V)$ within the quoted range. The initial positive correlation between line width and break strength arises because both the free-electron column density and the $n=2$ population increase with total gas column density. At sufficiently large column densities, however, the extreme scattered wings (requiring multiple scatterings) are decreased due to continuum absorption opacity. The observed data appear to follow this predicted pattern.
• Figure 5: Symmetry of the H$\alpha$ scattering wings as a function of the outflow velocity. The blue curves show the line asymmetry predicted for spherical kinematics. Sample sources with significant blueshifts are plotted. All objects are strongly discrepant with a spherical outflow model, implying that both inflow and outflow must be present. The symmetry is measured by taking the ratio of the exponential slope on the blue side of the line to that on the red side. The outflow velocity is derived from the H$\alpha$ absorption line and is plotted in units of the thermal electron velocity. The models do not predict perfectly symmetric profiles, and indeed, several objects show small, but measurable asymmetry.