@article{21038,
  abstract     = {Little Red Dots (LRDs) are compact sources at z > 5 discovered through James Webb Space Telescope spectroscopy. Their spectra exhibit broad Balmer emission lines (>~1000 km s^−1), alongside absorption features and a pronounced Balmer break – evidence for a dense, neutral hydrogen medium, in which the n = 2 state is significantly populated. When interpreted as arising
from active galactic nucleus broad-line regions, inferred black hole masses from local scaling relations exceed expectations given their stellar masses, challenging models of early black hole–galaxy co-evolution. However, radiative transfer effects in dense media may also impact the formation of hydrogen emission lines. We model three scattering processes shaping hydrogen
line profiles: resonance scattering by hydrogen in the n = 2 state, Raman scattering of ultraviolet (UV) radiation by ground-state hydrogen, and Thomson scattering by free electrons. Using 3D Monte Carlo radiative transfer simulations, we examine their imprint on line shapes and ratios. Resonance scattering produces strong deviations from Case B flux ratios, clear differences
between Hα and Hβ, and encodes gas kinematics in line profiles but cannot broaden Hβ due to conversion to Paα. While Raman scattering can yield broad wings, scattering of the UV continuum is disfavoured given the absence of strong full width at half-maximum variations across transitions. Raman scattering of higher Lyman-series emission can produce Hα/Hβ wing
width ratios of  >~1.28, agreeing with observations. Thomson scattering can reproduce the observed >~ 1000 km s^−1 wings under plausible conditions – e.g. Te ∼ 10^4 K and Ne ∼ 10^24 cm^−2 – and lead to black hole mass overestimates by factors  10. Our results provide a framework for interpreting hydrogen lines in LRDs and similar systems.},
  author       = {Chang, Seok Jun and Gronke, Max and Matthee, Jorryt J and Mason, Charlotte},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  number       = {4},
  publisher    = {Oxford University Press},
  title        = {{Impact of resonance, Raman, and Thomson scattering on hydrogen line formation in Little Red Dots}},
  doi          = {10.1093/mnras/staf2131},
  volume       = {545},
  year         = {2026},
}

@article{21045,
  abstract     = {The abundant population of little red dots (LRDs), compact objects with red UV to optical colors and broad Balmer lines at high redshift, is revealing new insights into the properties of early active galactic nuclei (AGN). Perhaps the most surprising features of this population are the presence of Balmer absorption and ubiquitous strong Balmer breaks. Recent models link these features to an active supermassive black hole (SMBH) cocooned in very dense gas (NH ∼ 1024 cm−2). We present a stringent test of such models using VLT/MUSE observations of A2744-45924, the most luminous LRD known to date (LHα ≈ 1044 erg s−1), located behind the Abell-2744 lensing cluster at z = 4.464 (μ = 1.8). We detect a moderately extended Lyα nebula (h ≈ 5.7 pkpc), spatially offset from the point-like Hα seen by JWST by ≈1.6 pkpc. The Lyα emission is narrow (FWHM = 270 ± 15 km s−1), and faint (Lyα = 0.07Hα) compared to Lyα nebulae typically observed around quasars of similar luminosity. We detect compact N IV]λ1486 emission, spatially aligned with Hα, and a spatial shift in the far-UV continuum matching the Lyα offset. We discuss that Hα and Lyα have distinct physical origins: Hα originates from the AGN, while Lyα is powered by star formation. In the environment of A2744-45924, we identified four extended Lyα halos (Δz < 0.02, Δr < 100 pkpc). Their Lyα luminosities match the expectations based on Hα emission, and show no evidence for radiation from A2744-45924 affecting its surroundings. The lack of strong, compact, and broad Lyα and the absence of a luminous extended halo, suggest that the UV AGN light is obscured by dense gas cloaking the SMBH with a covering factor close to unity.},
  author       = {Torralba Torregrosa, Alberto and Matthee, Jorryt J and Pezzulli, Gabriele and Urrutia, Tanya and Gronke, Max and Mascia, Sara and D’Eugenio, Francesco and Di Cesare, Claudia and Eilers, Anna Christina and Greene, Jenny E. and Iani, Edoardo and Ishikawa, Yuzo and Mackenzie, Ruari and Naidu, Rohan P. and Navarrete, Benjamín and Kotiwale, Gauri},
  issn         = {1432-0746},
  journal      = {Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{A weak Ly α halo for an extremely bright little red dot. Indications of enshrouded supermassive black hole growth}},
  doi          = {10.1051/0004-6361/202555596},
  volume       = {705},
  year         = {2026},
}

@article{21341,
  abstract     = {We aim to characterise the mass-metallicity relation (MZR) and the 3D correlation between the stellar mass, metallicity, and star formation rate (SFR) known as the fundamental metallicity relation (FMR) for galaxies at 5 < z < 7. Using ∼800 [O III] selected galaxies from deep NIRCam grism surveys, we present our stacked measurements of direct-Te metallicities, which we used to test recent strong-line metallicity calibrations. Our measured direct-Te metallicities (0.1–0.2 Z⊙ for M★ ≈ 5 × 107 − 9 M⊙, respectively) match recent JWST/NIRSpec-based results. However, there are significant inconsistencies between observations and hydrodynamical simulations. We observe a flatter MZR slope than the SPHINX20 and FLARES simulations, which cannot be attributed to selection effects. With simple models, we show that the effect of an [O III] flux-limited sample on the observed shape of the MZR is strongly dependent on the FMR. If the FMR is similar to the one in the local Universe, the intrinsic high-redshift MZR should be even flatter than is observed. In turn, a 3D relation where SFR correlates positively with metallicity at fixed mass would imply an intrinsically steeper MZR. Our measurements indicate that metallicity variations at fixed mass show little dependence on the SFR, suggesting a flat intrinsic MZR. This could indicate that the low-mass galaxies at these redshifts are out of equilibrium and that metal enrichment occurs rapidly in low-mass galaxies. However, being limited by our stacking analysis, we are yet to probe the scatter in the MZR and its dependence on SFR. Large carefully selected samples of galaxies with robust metallicity measurements can put tight constraints on the high-redshift FMR and help us to understand the interplay between gas flows, star formation, and feedback in early galaxies.},
  author       = {Kotiwale, Gauri and Matthee, Jorryt J and Kashino, Daichi and Vijayan, Aswin P. and Torralba Torregrosa, Alberto and Di Cesare, Claudia and Iani, Edoardo and Bordoloi, Rongmon and Leja, Joel and Maseda, Michael V. and Tacchella, Sandro and Shivaei, Irene and Heintz, Kasper E. and Danhaive, A. Lola and Mascia, Sara and Kramarenko, Ivan and Navarrete, Benjamín and Mackenzie, Ruari and Naidu, Rohan P. and Sobral, David},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{Rapid, out-of-equilibrium metal enrichment indicated by a flat mass-metallicity relation at z ∼ 6 from NIRCam grism spectroscopy}},
  doi          = {10.1051/0004-6361/202556597},
  volume       = {706},
  year         = {2026},
}

@article{21342,
  abstract     = {JWST has revealed a stunning population of bright galaxies at surprisingly early epochs, z > 10,
where few such sources were expected. Here we present the most distant example of this class yet – MoM-z14, a luminous (MUV = −20.2) source in the COSMOS legacy field at zspec = 14.44+0.02−0.02 that expands the observational frontier to a mere 280 million years after the Big Bang. The redshift is confirmed with NIRSpec/prism spectroscopy through a sharp Lyman-α break and ≈ 3σ detections of five rest-UV emission lines. The number density of bright zspec ≈ 14 − 15 sources implied by our “Mirage or Miracle” survey spanning ≈ 350 arcmin2 is > 100× larger (182+329 −105×) than pre-JWST consensus models. The high EWs of UV lines (≈15−35˚A) signal a rising star-formation history, with a ≈10× increase in the last 5 Myr (SFR5Myr/SFR50Myr = 9.9 +3.0 −5.8). The source is extremely compact (circularized re = 74+15
−12 pc), and yet elongated (b/a = 0.25+0.11−0.06), suggesting an AGN is not the dominant source of UV light. The steep UV slope (β = −2.5 +0.2 −0.2) implies negligible dust attenuation
and a young stellar population. The absence of a strong damping wing provides tentative evidence that the immediate surroundings of MoM-z14 may be partially ionized at a redshift where virtually every reionization model predicts a ≈ 100% neutral fraction. The nitrogen emission and highly supersolar [N/C]> 1 hint at an abundance pattern similar to local globular clusters that may have once hosted luminous supermassive stars. Since this abundance pattern is also common among the most ancient stars born in the Milky Way, we may be directly witnessing the formation of such stars in dense clusters, connecting galaxy evolution across the entire sweep of cosmic time. },
  author       = {Naidu, Rohan P. and Oesch, Pascal A. and Brammer, Gabriel and Weibel, Andrea and Li, Yijia and Matthee, Jorryt J and Chisholm, John and Pollock, Clara L. and Heintz, Kasper E. and Johnson, Benjamin D. and Shen, Xuejian and Hviding, Raphael E. and Leja, Joel and Tacchella, Sandro and Ganguly, Arpita and Witten, Callum and Atek, Hakim and Belli, Sirio and Bose, Sownak and Bouwens, Rychard and Dayal, Pratika and Decarli, Roberto and De Graaff, Anna and Fudamoto, Yoshinobu and Giovinazzo, Emma and Greene, Jenny E. and Illingworth, Garth and Inoue, Akio K. and Kane, Sarah G. and Labbe, Ivo and Leonova, Ecaterina and Marques-Chaves, Rui and Meyer, Romain A. and Nelson, Erica J. and Roberts-Borsani, Guido and Schaerer, Daniel and Simcoe, Robert A. and Stefanon, Mauro and Sugahara, Yuma and Toft, Sune and Van Der Wel, Arjen and Van Dokkum, Pieter and Walter, Fabian and Watson, Darach and Weaver, John R. and Whitaker, Katherine E.},
  issn         = {2565-6120},
  journal      = {The Open Journal of Astrophysics},
  publisher    = {Maynooth Academic Publishing},
  title        = {{A cosmic miracle: A remarkably luminous galaxy at zspec = 14.44 confirmed with JWST}},
  doi          = {10.33232/001c.156033},
  volume       = {9},
  year         = {2026},
}

@article{21371,
  abstract     = {There may be a newly identified early phase of supermassive black hole growth},
  author       = {Matthee, Jorryt J},
  issn         = {1095-9203},
  journal      = {Science},
  number       = {6787},
  pages        = {767--768},
  publisher    = {AAAS},
  title        = {{Black holes disguised as little red dots}},
  doi          = {10.1126/science.adz8603},
  volume       = {391},
  year         = {2026},
}

@article{21451,
  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 a 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 analyzed new deep JWST/NIRSpec PRISM and G395H spectra of FRESCO-GN-9771, one of the most luminous known LRDs at z = 5.5. These spectra reveal a strong Balmer break, broad Balmer lines, and very narrow [O III] emission. We revealed a forest of optical [Fe II] lines, which we argue are emerging from a dense (nH = 109 − 10 cm−3) warm layer with electron temperature Te ≈ 7000 K. The broad wings of Hα and Hβ have an exponential profile due to electron scattering in this same layer. The high Hα : Hβ : Hγ flux ratio of ≈10.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 of ∼5 M⊙ 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 broadband spectrum can be reasonably matched by a simple photoionized slab model that dominates the λ > 1500 Å continuum and a low-mass (∼108 M⊙) galaxy that could explain the narrow [O III], with only a subdominant contribution to the UV continuum. Our findings indicate that Balmer lines are not directly tracing the gas kinematics near the SMBH and that the BH mass scale is likely much lower than virial indicators suggest.},
  author       = {Torralba Torregrosa, Alberto and Matthee, Jorryt J and Pezzulli, Gabriele and Naidu, Rohan P. and Ishikawa, Yuzo and Brammer, Gabriel B. and Chang, Seok Jun and Chisholm, John and De Graaff, Anna and D’Eugenio, Francesco and Di Cesare, Claudia and Eilers, Anna Christina and Greene, Jenny E. and Gronke, Max and Iani, Edoardo and Kokorev, Vasily and Kotiwale, Gauri and Kramarenko, Ivan and Ma, Yilun and Mascia, Sara and Navarrete, Benjamín and Nelson, Erica and Oesch, Pascal and Simcoe, Robert A. and Wuyts, Stijn},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{The warm outer layer of a little red dot as the source of [Fe ii] and collisional Balmer lines with scattering wings}},
  doi          = {10.1051/0004-6361/202557537},
  volume       = {707},
  year         = {2026},
}

@article{21452,
  abstract     = {Galaxies exhibit a tight correlation between their star formation rate (SFR) and stellar mass over a wide redshift range known as the star-forming main sequence (SFMS). With JWST, the SFMS can now be investigated at high redshifts down to masses of ∼106 M⊙, using sensitive star formation rate tracers such as the Hα emission, which allow us to probe the variability in the star formation histories. We present inferences of the SFMS based on 316 Hα-selected galaxies at z ∼ 4 − 5 with log(M★/M⊙) = 6.4 − 10.6. These galaxies were identified behind the Abell 2744 lensing cluster with NIRCam grism spectroscopy from the survey All the Little Things (ALT). At face value, our data suggest a shallow slope in the SFMS (SFR ∝ M★α, with α = 0.45). After we corrected this for the Hα-flux limited nature of our survey using a Bayesian framework, the slope steepened to α = 0.59+0.10−0.09, whereas current data on their own are inconclusive on the mass dependence of the scatter. These slopes differ significantly from the slope of ∼1 that is expected from the observed evolution of the galaxy stellar mass function and from simulations. When we fixed the slope to α = 1, we found evidence for a decreasing intrinsic scatter with stellar mass (from ∼0.5 dex at M★ = 108 M⊙ to 0.4 dex at M★ = 1010 M⊙). This difference might be explained by a (combination of) luminosity-dependent SFR(Hα) calibration, a population of (mini)-quenched low-mass galaxies, or underestimated dust attenuation in high-mass galaxies. Future deep observations with different facilities can quantify these processes, which will enable us to achieve better insights into the variability of the star formation histories.},
  author       = {Di Cesare, Claudia and Matthee, Jorryt J and Naidu, Rohan P. and Torralba, Alberto and Kotiwale, Gauri and Kramarenko, Ivan and Blaizot, Jeremy and Rosdahl, Joakim and Leja, Joel and Iani, Edoardo and Adamo, Angela and Covelo-Paz, Alba and Furtak, Lukas J. and Heintz, Kasper E. and Mascia, Sara and Navarrete, Benjamín and Oesch, Pascal A. and Romano, Michael and Shivaei, Irene and Tacchella, Sandro},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{The slope and scatter of the star-forming main sequence at z ∼ 5: Reconciling observations with simulations}},
  doi          = {10.1051/0004-6361/202557790},
  volume       = {707},
  year         = {2026},
}

@article{21481,
  abstract     = {The Hα emission line in galaxies is a powerful tracer of their recent star formation activity. With the advent of JWST, we are now able to routinely observe Hα in galaxies at high redshift (z ≳ 3) and thus measure their star formation rates (SFRs). However, using classical SFR(Hα) calibrations to derive the SFRs leads to biased results because high-redshift galaxies are commonly characterized by low metallicities and bursty star formation histories, affecting the conversion factor between the Hα luminosity (LHα) and the SFR. We developed a set of new SFR(Hα) calibrations that allowed us to predict the SFRs of Hα-emitters at z ≳ 3 with very little error. We used the SPHINX cosmological simulations to select a sample of star-forming galaxies representative of the Hα-emitter population observed with JWST. We then derived linear corrections to the classical SFR(Hα) calibrations that took variations in the physical properties (e.g., stellar metallicities) among individual galaxies into account. We obtained two new SFR(Hα) calibrations that compared to the classical calibrations reduce the root mean squared error (RMSE) in the predicted SFRs by ΔRMSE ≈ 0.04 dex and ΔRMSE ≈ 0.06 dex, respectively. Using the recent JWST NIRCam/grism observations of Hα-emitters at z ∼ 6, we show that the new calibrations affect the high-redshift galaxy population statistics: (i) the estimated cosmic SFR density decreases by ΔρSFR ≈ 12%, and (ii) the observed slope of the star formation main sequence increases by Δ∂logSFR/∂logM★ = 0.08 ± 0.02.},
  author       = {Kramarenko, Ivan and Rosdahl, J. and Blaizot, J. and Matthee, Jorryt J and Katz, H. and Di Cesare, Claudia},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{H α as a tracer of star formation in the SPHINX cosmological simulations}},
  doi          = {10.1051/0004-6361/202557114},
  volume       = {707},
  year         = {2026},
}

@article{21715,
  abstract     = {New populations of red active galactic nuclei (known as “little red dots”) discovered by JWST exhibit remarkable spectral energy distributions. Leveraging X-ray through far-infrared observations of two of the most luminous known little red dots, we directly measure their bolometric luminosities. We find evidence that more than half of the bolometric luminosity likely emerges in the rest-frame optical, with Lbol/L5100 = 5, roughly half the value for “standard” active galactic nuclei. Meanwhile, the X-ray emitting corona, UV-emitting blackbody, and reprocessed mid to far-infrared emission are all considerably subdominant, assuming that the far-infrared luminosity is well below current measured limits. We present new bolometric corrections that dramatically lower inferred bolometric luminosities by a factor of 10 compared to published values in the literature. These bolometric corrections are in accord with expectations from models in which gas absorption and reprocessing are responsible for the red rest-frame optical colors of little red dots. We discuss how this lowered luminosity scale suggests a lower mass scale for the population by at least an order of magnitude (e.g., ∼105–107 M⊙ black holes, and ∼108 M⊙ galaxies), alleviating tensions with clustering, overmassive black holes, and the integrated black hole mass density in the Universe.},
  author       = {Greene, Jenny E. and Setton, David J. and Furtak, Lukas J. and Naidu, Rohan P. and Volonteri, Marta and Dayal, Pratika and Labbe, Ivo and Van Dokkum, Pieter and Bezanson, Rachel and Brammer, Gabriel and Cutler, Sam E. and Glazebrook, Karl and De Graaff, Anna and Hirschmann, Michaela and Hviding, Raphael E. and Kokorev, Vasily and Leja, Joel and Liu, Hanpu and Ma, Yilun and Matthee, Jorryt J and Nanayakkara, Themiya and Oesch, Pascal A. and Pan, Richard and Price, Sedona H. and Spilker, Justin S. and Wang, Bingjie and Weaver, John R. and Whitaker, Katherine E. and Williams, Christina C. and Zitrin, Adi},
  issn         = {1538-4357},
  journal      = {The Astrophysical Journal},
  number       = {2},
  publisher    = {IOP Publishing},
  title        = {{What you see is what you get: Empirically measured bolometric luminosities of Little Red Dots}},
  doi          = {10.3847/1538-4357/ae1836},
  volume       = {996},
  year         = {2026},
}

@article{21846,
  abstract     = {We compile a sample of 83 little red dots (LRDs) with JWST imaging and find that a substantial fraction (∼43%, rising to ≳80% for the most luminous LRDs) host one or more spatially offset, UV-bright companions at projected separations of 0.5 kpc ≲ d ≲ 5 kpc, with median 〈d〉 = 1.0 kpc. This fraction is even higher when smaller spatial scales are probed at high signal-to-noise ratio: the two most strongly lensed LRDs, A383-LRD1 and the newly discovered A68-LRD1, both have UV-bright companions at separations of only d ∼ 0.3 kpc, below the resolution limit of most unlensed JWST samples. We explore whether these ubiquitous red/blue configurations may be physically linked to the formation of LRDs, in analogy with the “synchronized pair” scenario originally proposed for direct-collapse black hole formation. In this picture, UV radiation from the companions, with typically modest stellar masses (M∗ ∼ 108−109 M⊙), suppresses molecular hydrogen cooling in nearby gas, allowing nearly isothermal collapse and the formation of extremely compact objects, such as massive black holes, supermassive stars, or quasi-stars. Using component-resolved photometry and spectral energy distribution modeling, we infer Lyman–Werner radiation fields of J21,LW ∼ 102.5–105 at the locations of the red components, comparable to those required in direct-collapse models, suggesting that the necessary photodissociation conditions are realized in many LRD systems. This framework provides a simple and self-consistent explanation for the extreme compactness and distinctive spectral properties of LRDs and links long-standing theoretical models for early compact object formation directly to a population now observed with JWST in the early Universe.},
  author       = {Baggen, Josephine F.W. and Scoggins, Matthew T. and Van Dokkum, Pieter and Haiman, Zoltán and Torralba Torregrosa, Alberto and Matthee, Jorryt J},
  issn         = {2041-8213},
  journal      = {The Astrophysical Journal Letters},
  number       = {1},
  publisher    = {IOP Publishing},
  title        = {{Connecting the dots: UV-bright companions of Little Red Dots as Lyman–Werner sources enabling direct-collapse Black Hole formation}},
  doi          = {10.3847/2041-8213/ae58a5},
  volume       = {1002},
  year         = {2026},
}

@article{21951,
  abstract     = {The central engines of Little Red Dots (LRDs) may be “black hole stars” (BH*s), early stages of
black hole growth characterized by dense gas envelopes. So far, the most direct evidence for BH*s
comes from a handful of sources where the host galaxy is completely outshone as suggested by their
remarkably steep Balmer breaks. Here we present a novel scheme to disentangle BH*s from their
host galaxies assuming that the [O III]5008˚A line arises exclusively from the host. Using a sample
of 98 LRDs (z ≈ 2 − 9) with high quality NIRSpec/PRISM spectra, we demonstrate that the hostsubtracted median stack displays a Balmer break > 2× stronger than massive quiescent galaxies,
with the rest-optical continuum resembling a blackbody-like SED (Teff ≈ 4050 K, log(Lbol) ≈ 43.9
erg s−1
, Reff ≈ 1300 au). We measure a steep Balmer decrement (Hα/Hβ > 10) and numerous
density-sensitive features (e.g., Fe II, He I, O I). These are hallmark signatures of dense gas envelopes,
providing population-level evidence that BH*s indeed power LRDs. In the median LRD, BH*s account
for ∼ 20% of the UV emission, ∼ 50% at the Balmer break, and ∼ 90% at wavelengths longer
than Hα with the remainder arising from the host. BH*s preferentially reside in low-mass galaxies
(M⋆ ≈ 108 M⊙) undergoing recent starbursts, as evidenced by extreme emission line EWs (e.g.,
[O III]5008˚A≈ 1100˚A, C III]≈ 12˚A), thereby favoring BH* origins linked to star-formation. We show
V-shaped LRD selections are biased to high BH*/host fractions (≳ 60% at 5500˚A) – less dominant
BH*s may be powering JWST’s blue broad-line AGN. We find BH*s are so commonplace and transient
(duty cycle ∼ 1%, lifetime ∼ 10 Myrs) that every massive black hole may have once shone as a BH*.
},
  author       = {Sun, Wendy Q. and Naidu, Rohan P. and Matthee, Jorryt J and De Graaff, Anna and Chisholm, John and Greene, Jenny E. and Oesch, Pascal A. and Torralba Torregrosa, Alberto and Hviding, Raphael E. and Brammer, Gabriel and Simcoe, Robert A. and Bose, Sownak and Bouwens, Rychard and Dayal, Pratika and Eilers, Anna Christina and Fei, Qinyue and Furtak, Lukas J. and Gottumukkala, Rashmi and Goulding, Andy and Heintz, Kasper E. and Hirschmann, Michaela and Kokorev, Vasily and Leja, Joel and Liu, Zhaoran and Natarajan, Priyamvada and Santarelli, Andrew D. and Setton, David J. and Smith, Aaron and Tacchella, Sandro and Volonteri, Marta and Walter, Fabian and Weibel, Andrea and Williams, Christina C.},
  issn         = {2565-6120},
  journal      = {The Open Journal of Astrophysics},
  publisher    = {Maynooth Academic Publishing},
  title        = {{Little Red Dot - Host Galaxy = Black Hole Star: A gas-enshrouded heart at the center of every Little Red Dot}},
  doi          = {10.33232/001c.162505},
  volume       = {9},
  year         = {2026},
}

@article{20864,
  abstract     = {Studies of the distant Universe are providing key insights into our understanding of the formation of galaxies. The advent of the James Webb Space Telescope (JWST) has significantly enhanced our observational capabilities, leading to an expanded redshift frontier, providing unprecedented detail in the characterisation of early galaxies and enabling the discovery of new populations of accreting black holes. This review aims to provide an introduction to the basic processes and components that shape the observed spectra of galaxies, with a focus on their relevance to techniques with which high-redshift galaxies are selected. The review further introduces specific topics that have attracted significant attention in recent literature, including the discovery of highly efficient galaxy formation in the early Universe, the relation between galaxies and the process of reionization, new insights into the formation of the first stars and the enrichment of interstellar gas with heavy elements, and breakthroughs in our understanding of the origins of supermassive black holes.},
  author       = {Matthee, Jorryt J},
  issn         = {1366-5812},
  journal      = {Contemporary Physics},
  number       = {1-4},
  pages        = {116--151},
  publisher    = {Taylor & Francis},
  title        = {{JWST provides a new view of cosmic dawn: Latest developments in studies of early galaxies}},
  doi          = {10.1080/00107514.2025.2586370},
  volume       = {66},
  year         = {2025},
}

@article{21060,
  abstract     = {Compact, star-forming galaxies with high star formation rate surface densities (ΣSFR) are often efficient Lyman continuum (LyC) emitters at z ≤ 4.5, likely because intense stellar feedback creates low-density channels that allow photons to escape. Irregular or disturbed morphologies, such as those resulting from mergers, can also facilitate LyC escape by creating anisotropic gas distributions. We investigated the influence of galaxy morphology on LyC production and escape at redshifts 5 ≤ z ≤ 7 using observations from various James Webb Space Telescope (JWST) surveys. Our sample consists of 436 sources, which are predominantly low-mass (∼10^8.15 M), star-forming galaxies with ionizing photon efficiency (ξion) values consistent with canonical expectations. Since direct measurements of fesc are not possible during the Epoch of  Reionization (EoR), we predicted fesc for high-redshift galaxies by applying survival analysis to a subsample of LyC emitters from the Low-Redshift Lyman Continuum Survey (LzLCS), selected to be direct analogs of reionization-era galaxies. We find that these galaxies exhibit, on average, modest predicted escape fractions (∼0.04). In addition, we evaluated the correlation between morphological features and LyC emission. Our findings indicate that neither ξion nor the predicted fesc values show a significant correlation with the presence of merger signatures. This suggests that in low-mass galaxies at z ≥ 5, strong morphological disturbances are not the primary mechanism driving LyC emission and leakage. Instead, compactness and star formation activity likely play a more pivotal role in regulating LyC escape. },
  author       = {Mascia, Sara and Pentericci, L. and Llerena, M. and Calabrò, A. and Matthee, Jorryt J and Flury, S. and Pacucci, F. and Jaskot, A. and Amorín, R. O. and Bhatawdekar, R. and Castellano, M. and Cleri, N. and Costantin, L. and Davis, K. and Di Cesare, Claudia and Dickinson, M. and Fontana, A. and Guo, Y. and Giavalisco, M. and Holwerda, B. W. and Hu, W. and Huertas-Company, M. and Jung, Intae and Kartaltepe, J. and Kashino, D. and Koekemoer, A. M. and Lucas, R. A. and Lotz, J. and Napolitano, L. and Jogee, S. and Wilkins, S.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{Little impact of mergers and galaxy morphology on the production and escape of ionizing photons in the early Universe}},
  doi          = {10.1051/0004-6361/202553760},
  volume       = {701},
  year         = {2025},
}

@article{21062,
  abstract     = {JWST observations have unveiled faint active galactic nuclei (AGNs) at high redshift that provide insights into the formation of supermassive black holes (SMBHs). However, disentangling their stellar from AGN light is challenging. Here, we use an empirical approach to infer the average stellar mass of five faint broad-line (BL) Hα emitters at z = 4–5 with BH masses ≈6 × 10^6 M⊙, with a method independent of their spectral energy distribution (SED). We use the deep JWST/NIRcam grism survey “All the Little Things” to measure the overdensities around BL-Hα emitters and around a spectroscopic reference sample of ∼300 galaxies. In our reference sample, we find that megaparsec-scale overdensity correlates with stellar mass. Their large-scale environments suggest that BL-Hα emitters are hosted by galaxies with stellar masses ≈5 × 10^7 M⊙, ≈40 times lower than those inferred from galaxy-only SED fits. Adding measurements around more luminous z ≈ 6 AGNs, we find tentative correlations between line width, BH mass, and the overdensity, suggestive of a steep BH to halo mass relation. The main implications are (1) when BH masses are taken at face value, we confirm extremely high BH to stellar mass ratios of ≈10%, (2) the galaxies of low stellar mass that host growing SMBHs are in tension with typical hydrodynamical simulations, except those without feedback, (3) a 1% duty cycle implied by the host mass hints at super-Eddington accretion, (4) the masses are at odds with an interpretation of the line broadening in terms of high stellar density, (5) our results imply a luminosity-dependent diversity of galaxy masses, environments, and SEDs among AGN samples.},
  author       = {Matthee, Jorryt J and Naidu, Rohan P. and Kotiwale, Gauri and Furtak, Lukas J. and Kramarenko, Ivan and Mackenzie, Ruari and Greene, Jenny and Adamo, Angela and Bouwens, Rychard J. and Di Cesare, Claudia and Eilers, Anna-Christina and de Graaff, Anna and Heintz, Kasper E. and Kashino, Daichi and Maseda, Michael V. and Tacchella, Sandro and Torralba Torregrosa, Alberto},
  issn         = {1538-4357},
  journal      = {The Astrophysical Journal},
  number       = {2},
  publisher    = {IOP Publishing},
  title        = {{Environmental evidence for overly massive Black Holes in low-mass galaxies and a Black Hole–Halo mass relation at z ∼ 5}},
  doi          = {10.3847/1538-4357/ade886},
  volume       = {988},
  year         = {2025},
}

@article{19066,
  abstract     = {We present a sample of 1956 individual stellar clumps at redshift 0.7 < z < 10, detected with JWST/NIRCam in 476 galaxies lensed by the galaxy cluster Abell2744. The lensed clumps present magnifications ranging between μ = 1.8 and μ = 300. We perform simultaneous size-photometry estimates in 20 JWST/NIRCam median and broad-band filters from 0.7 to 5 μm.
Spectral energy distribution (SED) fitting analyses enable us to recover the physical properties of the clumps. The majority of the clumps are spatially resolved and have effective radii in the range Reff = 10–700 pc. We restrict this first study to the 1751 post-reionization era clumps with redshift < 5.5. We find a significant evolution of the average clump ages, star formation rates (SFRs), SFR surface densities, and metallicity with increasing redshift, while median stellar mass and stellar mass surface densities are similar in the probed redshift range. We observe a strong correlation between the clump properties and the properties of their host galaxies, with more massive galaxies hosting more massive and older clumps. We find that clumps closer to their host galactic centre are on average more massive, while their ages do not show clear sign of migration. We find that clumps at cosmic noon sample the upper-mass end of the mass function to higher masses than at z > 3, reflecting the rapid increase towards the peak of the cosmic star formation history. We conclude that the results achieved over the studied redshift range are in agreement with expectation of in situ clump formation scenario from large-scale disc fragmentation. },
  author       = {Claeyssens, Adélaïde and Adamo, Angela and Messa, Matteo and Dessauges-Zavadsky, Miroslava and Richard, Johan and Kramarenko, Ivan and Matthee, Jorryt J and Naidu, Rohan P.},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  number       = {3},
  pages        = {2535--2558},
  publisher    = {Oxford University Press},
  title        = {{Tracing star formation across cosmic time at tens of parsec-scales in the lensing cluster field Abell 2744}},
  doi          = {10.1093/mnras/staf058},
  volume       = {537},
  year         = {2025},
}

@article{20425,
  abstract     = {Ultraviolet (UV) radiation from accreting black holes ionizes the intergalactic gas around early quasars, carving out highly ionized bubbles in their surroundings. Any changes in a quasar’s luminosity are therefore predicted to produce outward-propagating ionization gradients, affecting the Lyα absorption opacity near the quasar’s systemic redshift. This “proximity effect” is well-documented in rest-UV quasar spectra but only provides a one-dimensional probe along our line of sight. Here we present deep spectroscopic observations with the James Webb Space Telescope (JWST) of galaxies in the background of a superluminous quasar at zQSO ≈ 6.3, which reveal the quasar’s “light echo” with Lyα tomography in the transverse direction. This transverse proximity effect is detected for the first time toward multiple galaxy sightlines, allowing us to map the extent and geometry of the quasar’s ionization cone. We obtain constraints on the orientation and inclination of the cone, as well as an upper limit on the obscured solid angle fraction of fobsc < 91%. Additionally, we find a timescale of the quasar’s UV radiation of tqso = 10^5.6+0.1-0.3 yr, which is significantly shorter than would be required to build up the central supermassive black hole (SMBH) with conventional growth models, but is consistent with independent measurements of the quasars’ duty cycle. Our inferred obscured fraction disfavors a scenario where short quasar lifetimes can be explained exclusively by geometric obscuration, and instead supports the idea that radiatively inefficient accretion or growth in initially heavily enshrouded cocoons plays a pivotal role in early SMBH growth. Our results pave the way for novel studies of quasars’ ionizing geometries and radiative histories at early cosmic times.},
  author       = {Eilers, Anna Christina and Yue, Minghao and Matthee, Jorryt J and Hennawi, Joseph F. and Davies, Frederick B. and Simcoe, Robert A. and Teague, Richard and Bordoloi, Rongmon and Brammer, Gabriel and Kang, Yi and Kashino, Daichi and Mackenzie, Ruari and Naidu, Rohan P. and Navarrete, Benjamín},
  issn         = {2041-8213},
  journal      = {The Astrophysical Journal Letters},
  number       = {2},
  publisher    = {IOP Publishing},
  title        = {{The light echo of a high-redshift quasar mapped with Lyα tomography}},
  doi          = {10.3847/2041-8213/ae057a},
  volume       = {991},
  year         = {2025},
}

@article{15180,
  abstract     = {Characterizing the prevalence and properties of faint active galactic nuclei (AGNs) in the early Universe is key for understanding the formation of supermassive black holes (SMBHs) and determining their role in cosmic reionization. We perform a spectroscopic search for broad Hα emitters at z ≈ 4–6 using deep JWST/NIRCam imaging and wide field slitless spectroscopy from the EIGER and FRESCO surveys. We identify 20 Hα lines at z = 4.2–5.5 that have broad components with line widths from ∼1200–3700 km s−1, contributing ∼30%–90% of the total line flux. We interpret these broad components as being powered by accretion onto SMBHs with implied masses ∼107–8M⊙. In the UV luminosity range MUV,AGN+host = −21 to −18, we measure number densities of ≈10−5 cMpc−3. This is an order of magnitude higher than expected from extrapolating quasar UV luminosity functions (LFs). Yet, such AGN are found in only <1% of star-forming galaxies at z ∼ 5. The number density discrepancy is much lower when compared to the broad Hα LF. The SMBH mass function agrees with large cosmological simulations. In two objects, we detect complex Hα profiles that we tentatively interpret as caused by absorption signatures from dense gas fueling SMBH growth and outflows. We may be witnessing early AGN feedback that will clear dust-free pathways through which more massive blue quasars are seen. We uncover a strong correlation between reddening and the fraction of total galaxy luminosity arising from faint AGN. This implies that early SMBH growth is highly obscured and that faint AGN are only minor contributors to cosmic reionization.},
  author       = {Matthee, Jorryt J and Naidu, Rohan P. and Brammer, Gabriel and Chisholm, John and Eilers, Anna-Christina and Goulding, Andy and Greene, Jenny and Kashino, Daichi and Labbe, Ivo and Lilly, Simon J. and Mackenzie, Ruari and Oesch, Pascal A. and Weibel, Andrea and Wuyts, Stijn and Xiao, Mengyuan and Bordoloi, Rongmon and Bouwens, Rychard and van Dokkum, Pieter and Illingworth, Garth and Kramarenko, Ivan and Maseda, Michael V. and Mason, Charlotte and Meyer, Romain A. and Nelson, Erica J. and Reddy, Naveen A. and Shivaei, Irene and Simcoe, Robert A. and Yue, Minghao},
  issn         = {1538-4357},
  journal      = {The Astrophysical Journal},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  number       = {2},
  publisher    = {American Astronomical Society},
  title        = {{Little Red Dots: An abundant population of faint active galactic nuclei at z ∼ 5 revealed by the EIGER and FRESCO JWST surveys}},
  doi          = {10.3847/1538-4357/ad2345},
  volume       = {963},
  year         = {2024},
}

@article{15249,
  abstract     = {Observationally mapping the relation between galaxies and the intergalactic medium (IGM) is of key interest for studies of cosmic reionization. Diffuse hydrogen gas has typically been observed in H I Lyman-α (Lyα) absorption in the spectra of bright background quasars. However, it is important to extend these measurements to background galaxies as quasars become increasingly rare at high redshift and rarely probe closely separated sight lines. Here, we use deep integral field spectroscopy in the MUSE eXtremely Deep Field to demonstrate the measurement of the Lyα transmission at z ≈ 4 in absorption to a background galaxy at z = 4.77. The H I transmission is consistent with independent quasar sight lines at similar redshifts. Exploiting the high number of spectroscopic redshifts of faint galaxies (500 between z = 4.0–4.7 within a radius of 8 arcmin) that are tracers of the density field, we show that Lyα transmission is inversely correlated with galaxy density, i.e. transparent regions in the Lyα forest mark underdense regions at z ≈ 4. Due to large-scale clustering, galaxies are surrounded by excess H I absorption over the cosmic mean out to 4 cMpc/h70. We also find that redshifts from the peak of the Lyα line are typically offset from the systemic redshift by +170 km s−1. This work extends results from z ≈ 2–3 to higher redshifts and demonstrates the power of deep integral field spectroscopy to simultaneously measure the ionization structure of the IGM and the large-scale density field in the early Universe.},
  author       = {Matthee, Jorryt J and Golling, Christopher and Mackenzie, Ruari and Pezzulli, Gabriele and Lilly, Simon and Schaye, Joop and Bacon, Roland and Kusakabe, Haruka and Urrutia, Tanya and Boogaard, Leindert and Brinchmann, Jarle and Maseda, Michael V. and Garel, Thibault and Bouché, Nicolas F. and Wisotzki, Lutz},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  number       = {3},
  pages        = {2794--2806},
  publisher    = {Oxford University Press},
  title        = {{Large-scale excess H I absorption around z ≈ 4 galaxies detected in a background galaxy spectrum in the MUSE eXtremely deep field}},
  doi          = {10.1093/mnras/stae673},
  volume       = {529},
  year         = {2024},
}

@article{18494,
  abstract     = {We expect luminous (M 1450 ≲ −26.5) high-redshift quasars to trace the highest-density peaks in the early Universe. Here, we present observations of four z ≳ 6 quasar fields using JWST/NIRCam in the imaging and wide-field slitless spectroscopy mode and report a wide range in the number of detected [O iii]-emitting galaxies in the quasars’ environments, ranging between a density enhancement of δ ≈ 65 within a 2 cMpc radius—one of the largest protoclusters during the Epoch of Reionization discovered to date—to a density contrast consistent with zero, indicating the presence of a UV-luminous quasar in a region comparable to the average density of the Universe. By measuring the two-point cross-correlation function of quasars and their surrounding galaxies, as well as the galaxy autocorrelation function, we infer a correlation length of quasars at 〈z〉 = 6.25 of r 0 QQ = 22.0 − 2.9 + 3.0 cMpc h − 1 , while we obtain a correlation length of the [O iii]-emitting galaxies of r 0 GG = 4.1 ± 0.3 cMpc h − 1 . By comparing the correlation functions to dark-matter-only simulations we estimate the minimum mass of the quasars’ host dark matter halos to be log 10 ( M halo , min / M ⊙ ) = 12.43 − 0.15 + 0.13 (and log 10 ( M halo , min [ OIII ] / M ⊙ ) = 10.56 − 0.03 + 0.05 for the [O iii] emitters), indicating that (a) luminous quasars do not necessarily reside within the most overdense regions in the early Universe, and that (b) the UV-luminous duty cycle of quasar activity at these redshifts is f duty ≪ 1. Such short quasar activity timescales challenge our understanding of early supermassive black hole growth and provide evidence for highly dust-obscured growth phases or episodic, radiatively inefficient accretion rates.},
  author       = {Eilers, Anna Christina and Mackenzie, Ruari and Pizzati, Elia and Matthee, Jorryt J and Hennawi, Joseph F. and Zhang, Haowen and Bordoloi, Rongmon and Kashino, Daichi and Lilly, Simon J. and Naidu, Rohan P. and Simcoe, Robert A. and Yue, Minghao and Frenk, Carlos S. and Helly, John C. and Schaller, Matthieu and Schaye, Joop},
  issn         = {1538-4357},
  journal      = {Astrophysical Journal},
  number       = {2},
  publisher    = {IOP Publishing},
  title        = {{EIGER. VI. The correlation function, host halo mass, and duty cycle of luminous quasars at z ≳ 6}},
  doi          = {10.3847/1538-4357/ad778b},
  volume       = {974},
  year         = {2024},
}

