@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{21380,
  abstract     = {Context. Extreme emission line galaxies (EELGs) are believed to significantly contribute to the star formation activity and mass assembly in galaxies. EELGs likely also play a leading role in the cosmic re-ionization as their interstellar medium may allow a significant fraction of their ionizing photons to escape (> 5%). Finding low-redshift analogues of these high-z galaxies is therefore essential to characterizing the physical conditions in the interstellar medium of these galaxies and understanding the processes that re-ionized the Universe.

Aims. We aimed to develop a robust and efficient method for the photometric identification of EELGs using the J-PAS survey. J-PAS will cover approximately 8500 deg2 of the sky with 54 narrow-band filters in the optical range plus i-SDSS, enabling detailed studies of the physical properties of these galaxies. In this work we focused on an initial subset of the survey: a 30 square degree area with complete observations in all bands.

Methods. We combine equivalent width (EW) measurements from J-PAS narrow-band photometry with artificial intelligence techniques to identify galaxies with emission lines exceeding 300 Å. We validated our selection using spectroscopic data from DESI DR1 and characterized the selected sample through spectral energy distribution fitting with CIGALE.

Results. We identify 917 EELGs up to z = 0.8 over 30 deg2, achieving a purity of 95% and a completeness of 96% for i-SDSS < 22.5 mag. Importantly, active galactic nucleus contamination was carefully considered and is estimated to be around 5%. Furthermore, a cross-match with DESI yielded 79 counterparts; their redshifts are in excellent agreement with our photometric estimates, thereby confirming the reliability of our redshift determination. In addition, the derived emission line fluxes are in good agreement with spectroscopic measurements. Moreover, the selected sample reveals strong correlations between the ionizing photon production efficiency (ξion) and EW(Hβ), which are consistent with previous observational studies at low and high redshifts and theoretical expectations. Finally, most of the sources surpass the ionizing efficiency threshold required for re-ionization, highlighting their relevance as local analogues of early-Universe galaxies.},
  author       = {Giménez-Alcázar, A. and Amorín, R. and Vílchez, J. M. and Hernán-Caballero, A. and González-Otero, M. and Arroyo-Polonio, A. and Iglesias-Páramo, J. and Lumbreras-Calle, A. and Fernández-Ontiveros, J. A. and López-Sanjuan, C. and Bonatto, L. and González Delgado, R. M. and Kehrig, C. and Torralba Torregrosa, Alberto and Rahna, P. T. and Jiménez-Teja, Y. and Márquez, I. and Breda, I. and Álvarez-Candal, A. and Abramo, R. and Alcaniz, J. and Benitez, N. and Bonoli, S. and Carneiro, S. and Cenarro, J. and Cristóbal-Hornillos, D. and Dupke, R. and Ederoclite, A. and Hernández-Monteagudo, C. and Marín-Franch, A. and Mendes de Oliveira, C. and Moles, M. and Sodré, L. and Taylor, K. and Varela, J. and Vázquez Ramió, H.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{J-PAS: First identification, physical properties, and ionization efficiency of extreme emission line galaxies}},
  doi          = {10.1051/0004-6361/202557358},
  volume       = {706},
  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{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{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 & 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{22263,
  abstract     = {Recent studies at high redshift have revealed an enigmatic class of little red dots (LRDs) with extreme Balmer breaks, stronger than in any stellar atmosphere. However, it is unclear whether such objects exist at lower redshift, especially given the low number of LRDs reported at z ≲ 2. Here, we report the discovery of PAN-BH*-1, an LRD with an extreme Balmer break at z = 1.73, identified from JWST/NIRCam pure-parallel imaging taken by the PANORAMIC survey, and confirmed by deep VLT/X-Shooter spectroscopy. The rest-optical to near-infrared spectral energy distribution of PAN-BH*-1 is consistent with a photospheric continuum with effective temperature Teff ≈ 4800 K. The broad Hα emission line shows remarkably deep absorption, stronger than previously measured in any LRD. The absorption trough spans from −520 to +267 km s−1 with respect to the systemic redshift. The presence of blue- and red-shifted absorption suggests complex dynamics of the obscuring gas along the line of sight. We speculate that the absorption trough can be produced by a thick wind launched from a thick, rotating photospheric disk, the latter being the source of the red optical continuum. While the source is unresolved in the rest-optical JWST data (reff < 47 pc), the rest-near-UV Hubble Space Telescope imaging shows an extended morphology with (formular displayed) kpc, which we interpret as a host galaxy with a stellar mass of ∼10^8 M⊙, in line with the narrow Hα emission. The discovery of this object at cosmic noon highlights the feasibility of systematic searches for extreme LRDs with wide-area facilities such as Euclid and Roman.},
  author       = {Torralba Torregrosa, Alberto and Matthee, Jorryt J and Weibel, Andrea and Naidu, Rohan P. and Ma, Yilun and Cloonan, Aidan P. and Desai, Aayush A and De Graaff, Anna and Greene, Jenny E. and Jespersen, Christian Kragh and Kramarenko, Ivan and Mascia, Sara and Oesch, Pascal A. and Sun, Wendy Q. and Williams, Christina C.},
  issn         = {2041-8213},
  journal      = {The Astrophysical Journal Letters},
  number       = {2},
  publisher    = {IOP Publishing},
  title        = {{A black hole star at cosmic noon: Extreme Balmer break, photospheric continuum, and broad absorption by thick winds in a Little Red Dot at z = 1.7}},
  doi          = {10.3847/2041-8213/ae7bfd},
  volume       = {1005},
  year         = {2026},
}

@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{19929,
  abstract     = {Context. The observed Lyman-alpha (Lyα) line profile is a convolution of the complex Lyα radiative transfer taking place in the interstellar, circumgalactic, and intergalactic media (ISM, CGM, and IGM, respectively). Discerning the different components of the Lyα line is crucial in order to use it as a probe of galaxy formation or the evolution of the IGM.

Aims. We aim to present the second version of zELDA (redshift Estimator for Line profiles of Distant Lyman-Alpha emitters), an open-source Python module focused on modelling and fitting observed Lyα line profiles. This new version of zELDA focuses on disentangling the galactic from the IGM effects.

Methods. We built realistic Lyα line profiles that include the ISM and IGM contributions by combining the Monte Carlo radiative-transfer simulations for the so-called shell model (ISM) and IGM transmission curves generated from TNG100. We used these mock line profiles to train different artificial neural networks. These use the observed spectrum as input and the outflow parameters of the best fitting ‘shell model’ as output along with the redshift and Lyα emission IGM escape fraction of the source.

Results. We measured the accuracy of zELDA on mock Lyα line profiles. We find that zELDA is capable of reconstructing the ISM emerging Lyα line profile with high levels of accuracy (Kolmogórov-Smirnov<0.1) for 95% of the cases for HST/COS-like observations and 80% for MUSE-WIDE-like observations. zELDA is able to measure the IGM transmission with typical uncertainties below 10% for HST/COS and MUSE-WIDE data.

Conclusions. This work represents a step forward in the high-precision reconstruction of IGM-attenuated Lyα line profiles. zELDA allows the disentanglement of the galactic and IGM contribution shaping the Lyα line shape and thus allows us to use Lyα as a tool to study galaxy and ISM evolution.},
  author       = {Gurung-López, Siddhartha and Byrohl, Chris and Gronke, Max and Spinoso, Daniele and Torralba Torregrosa, Alberto and Fernández-Soto, Alberto and Arnalte-Mur, Pablo and Martínez, Vicent J.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{zELDA II: Reconstruction of galactic Lyman-alpha spectra attenuated by the intergalactic medium using neural networks}},
  doi          = {10.1051/0004-6361/202453547},
  volume       = {698},
  year         = {2025},
}

