@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{21727, abstract = {We present a comprehensive analysis of the MIRI Extremely Red Object Virgil, a Lyα emitter at zspec = 6.6379 ± 0.0035 with the photometric properties of a Little Red Dot. Leveraging new JWST/MIRI imaging from the MIDIS and PAHSPECS programs, we confirm Virgil’s extraordinary nature among galaxies in JADES/GOODS-South, exhibiting a strikingly red NIRCam-to-MIRI color (F444W–F1500W = 2.84 ± 0.04 mag). Deep NIRSpec/PRISM spectroscopy from the OASIS program offers key insights into the host galaxy, revealing properties of an average star-forming galaxy during Cosmic Reionization, such as a subsolar metallicity, low-to-moderate dust content, and a relatively high ionization parameter and electron temperature. By estimating the star formation rate of Virgil from UV and Hα, we find evidence that the galaxy is either entering or fading out of a bursty episode. Although line-ratio diagnostics employed at high z would classify Virgil as an active galactic nucleus (AGN), this classification becomes ambiguous once redshift evolution is considered. Nonetheless, Virgil occupies the same parameter space as recently confirmed AGNs at similar redshifts. The new deep MIRI data at 15 μm reinforce the AGN nature of Virgil, as inferred from multiple spectral energy distribution (SED) fitting codes. Virgil’s rising infrared SED and UV excess resemble those of Dust-Obscured Galaxies (DOGs) studied with Spitzer at Cosmic Noon, particularly blue-excess HotDOGs. Our results highlight the need for a multiwavelength approach incorporating MIRI to uncover such extreme sources at z ≳ 6 and to shed light on the interplay between galaxy evolution and early black hole growth during Cosmic Reionization.}, author = {Rinaldi, Pierluigi and Pérez-González, Pablo G. and Rieke, George H. and Lyu, Jianwei and D’Eugenio, Francesco and Wu, Zihao and Carniani, Stefano and Looser, Tobias J. and Shivaei, Irene and Boogaard, Leindert A. and Diaz-Santos, Tanio and Colina, Luis and Östlin, Göran and Alberts, Stacey and Álvarez-Márquez, Javier and Annuziatella, Marianna and Aravena, Manuel and Bhatawdekar, Rachana and Bunker, Andrew J. and Caputi, Karina I. and Charlot, Stéphane and Crespo Gómez, Alejandro and Curti, Mirko and Eckart, Andreas and Gillman, Steven and Hainline, Kevin and Kumari, Nimisha and Hjorth, Jens and Iani, Edoardo and Inami, Hanae and Ji, Zhiyuan and Johnson, Benjamin D. and Jones, Gareth C. and Labiano, Álvaro and Maiolino, Roberto and Melinder, Jens and Moutard, Thibaud and Peissker, Florian and Rieke, Marcia and Robertson, Brant and Scholtz, Jan and Tacchella, Sandro and Van Der Werf, Paul P. and Walter, Fabian and Williams, Christina C. and Willott, Chris and Witstok, Joris and Übler, Hannah and Zhu, Yongda}, issn = {1538-4357}, journal = {The Astrophysical Journal}, number = {1}, publisher = {IOP Publishing}, title = {{Deciphering the nature of Virgil: An obscured active galactic nucleus lurking within an apparently normal Lyα emitter during cosmic reionization}}, doi = {10.3847/1538-4357/ae089c}, volume = {994}, year = {2025}, } @article{17889, abstract = {The increasingly neutral intergalactic gas at z > 6 impacts the Lyman-α (Lyα) flux observed from galaxies. One luminous galaxy, COLA1, stands out because of its unique double-peaked Lyα line at z = 6.6, unseen in any simulation of reionization. Here, we present JWST/NIRCam wide-field slitless spectroscopy in a 21 arcmin2 field centered on COLA1. We find 141 galaxies spectroscopically selected through the [O III] doublet at 5.35 < z < 6.95, with 40 of these sources showing Hβ. For COLA1, we additionally detect [O III]4363 as well as Hγ. We measure a systemic redshift of z = 6.5917 for COLA1, confirming the classical double-peak nature of the Lyα profile. This implies that it resides in a highly ionized bubble and that it is leaking ionizing photons with a high escape fraction of fesc(LyC) = 20–50%, making it a prime laboratory to study Lyman continuum escape in the Epoch of Reionization. COLA1 shows all the signs of a prolific ionizer with a Lyα escape fraction of 81 ± 5%, Balmer decrement indicating no dust, a steep UV slope (βUV = −3.2 ± 0.4), and a star-formation surface density ≳10× that of typical galaxies at similar redshift. We detect five galaxies in COLA1’s close environment (Δz < 0.02). Exploiting the high spectroscopic completeness inherent to grism surveys, and using mock simulations that fully mimic the selection function, we show that the number of detected companions is very typical for a normal similarly UV-bright (MUV ∼ −21.3) galaxy – that is, the ionized bubble around COLA1 is unlikely to be due to an excessively large over-density. Instead, the measured ionizing properties suggest that COLA1 by itself might be powering the bubble required to explain its double-peaked Lyα profile (Rion ≈ 0.7 pMpc), with only minor contributions from detected neighbors (−19.5 ≲ MUV ≲ −17.5).}, author = {Torralba-Torregrosa, Alberto and Matthee, Jorryt J and Naidu, Rohan P. and Mackenzie, Ruari and Pezzulli, Gabriele and Hutter, Anne and Arnalte-Mur, Pablo and Gurung-López, Siddhartha and Tacchella, Sandro and Oesch, Pascal and Kashino, Daichi and Conroy, Charlie and Sobral, David}, issn = {1432-0746}, journal = {Astronomy and Astrophysics}, publisher = {EDP Sciences}, title = {{Anatomy of an ionized bubble: NIRCam grism spectroscopy of the z = 6.6 double-peaked Lyman- α emitter COLA1 and its environment}}, doi = {10.1051/0004-6361/202450318}, volume = {689}, year = {2024}, } @article{18069, abstract = {We present results from the JWST First Reionization Epoch Spectroscopically Complete Observations survey on the star-forming sequence (SFS) of galaxies at 1.0 < z < 1.7, around the peak of the cosmic star formation history. Star formation rates (SFRs) are measured from the redshifted, relatively dust-insensitive Paschen-α emission line, and stellar mass measurements include the F444W (4.4 μm; rest-frame H) band. We find SFRs of galaxies with log(M*/M⊙) > 9.5 that are lower than found in many earlier studies by up to 0.6 dex, but in good agreement with recent results obtained with the Prospector fitting framework. The difference (log(SFR(Paα)-SFR(Prospector)) is −0.09 ± 0.04 dex at 1010−11M⊙. We also measure the empirical relation between Paschen-α luminosity and rest-frame H-band magnitude and find that the scatter is only 0.04 dex lower than that of the SFR–M* relation and is much lower than the systematic differences among relations in the literature due to various methods of converting observed measurements to physical properties. We additionally identify examples of sources—that, with standard cutoffs via the UVJ diagram, would be deemed quiescent—with significant (log(sSFR)> −11 yr−1), typically extended, Paschen-α emission. Our results may be indicative of the potential unification of methods used to derive the SFS with careful selection of star-forming galaxies and independent SFR and stellar mass indicators.}, author = {Neufeld, Chloe and Van Dokkum, Pieter and Asali, Yasmeen and Covelo-Paz, Alba and Leja, Joel and Lin, Jamie and Matthee, Jorryt J and Oesch, Pascal A. and Reddy, Naveen A. and Shivaei, Irene and Whitaker, Katherine E. and Wuyts, Stijn and Brammer, Gabriel and Marchesini, Danilo and Maseda, Michael V. and Naidu, Rohan P. and Nelson, Erica J. and Velichko, Anna and Weibel, Andrea and Xiao, Mengyuan}, issn = {1538-4357}, journal = {Astrophysical Journal}, number = {2}, publisher = {IOP Publishing}, title = {{FRESCO: The Paschen-α star-forming sequence at cosmic noon}}, doi = {10.3847/1538-4357/ad6158}, volume = {972}, year = {2024}, } @article{18447, abstract = {We present the first results on the spatial distribution of star formation in 454 star-forming galaxies just after the epoch of reionisation (4.8 < z < 6.5) using Hα emission-line maps and F444W imaging that traces the stellar continuum from the JWST FRESCO NIRCam Slitless Spectroscopy Survey. The Hα equivalent width profiles of star-forming galaxies across the main sequence at z ∼ 5.3 with stellar masses 6.8≤ log(M*/M⊙) < 11.1 increase with radius, which provides direct evidence for the inside-out growth of star-forming galaxies just after the epoch of reionisation. GALFIT was used to calculate half-light radii, Reff, and central surface densities within 1 kiloparsec, Σ1kpc of Hα and the continuum. At a fixed stellar mass of Log(M*/M⊙) = 9.5, Σ1kpc, Hα is 1.04 ± 0.05 times higher than Σ1kpc, C, Reff, Hα is 1.18 ± 0.03 times larger than Reff, C and both Reff measurements are smaller than 1 kiloparsec. These measurements suggest the rapid build-up of compact bulges via star formation just after the epoch of reionisation. By comparison to analogous work done at lower redshifts with Hubble Space Telescope WFC3 slitless spectroscopy as part of the 3D-HST (z ∼ 1) and CLEAR (z ∼ 0.5) surveys, we find that Reff(z) evolves at the same pace for Hα and the continuum, but Σ1kpc(z) evolves faster for Hα than the stellar continuum. As a function of the Hubble parameter, Reff, Hα/Reff,C = 1.1h(z) and Σ1 kpc, Hα/Σ1 kpc,C = h(z)1.3. These parametrisations suggest that the inside-out growth of the disk starts to dominate the inside-out growth of the bulge towards lower redshifts. This is supported by the redshift evolution in the EW(Hα) profiles from FRESCO, 3D-HST, and CLEAR at fixed stellar mass and when star-forming progenitors are traced, in which in EW(Hα) rapidly increases with radius within the half-light radius at z ∼ 5.3, but EW(Hα) increases only significantly with radius in the outer disk at z ∼ 0.5.}, author = {Matharu, Jasleen and Nelson, Erica J. and Brammer, Gabriel and Oesch, Pascal A. and Allen, Natalie and Shivaei, Irene and Naidu, Rohan P. and Chisholm, John and Covelo-Paz, Alba and Fudamoto, Yoshinobu and Giovinazzo, Emma and Herard-Demanche, Thomas and Kerutt, Josephine and Kramarenko, Ivan and Marchesini, Danilo and Meyer, Romain A. and Prieto-Lyon, Gonzalo and Reddy, Naveen and Shuntov, Marko and Weibel, Andrea and Wuyts, Stijn and Xiao, Mengyuan}, issn = {1432-0746}, journal = {Astronomy and Astrophysics}, publisher = {EDP Sciences}, title = {{A first look at spatially resolved star formation at 4.8 < z < 6.5 with JWST FRESCO NIRCam slitless spectroscopy}}, doi = {10.1051/0004-6361/202450522}, volume = {690}, year = {2024}, } @article{18448, abstract = {Aims. This paper utilises the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) to extend the observational studies of dust and polycyclic aromatic hydrocarbon (PAH) emission to a new mass and star formation rate (SFR) parameter space beyond our local Universe. The combination of fully sampled spectral energy distributions (SEDs) with multiple mid-infrared (mid-IR) bands and the unprecedented sensitivity of MIRI allows us to investigate dust obscuration and PAH behaviour from z = 0.7 up to z = 2 in typical main-sequence galaxies. Our focus is on constraining the evolution of PAH strength and the dust-obscured luminosity fraction before and during cosmic noon, the epoch of peak star formation activity in the Universe. Methods. In this study, we utilise MIRI multi-band imaging data from the SMILES survey (5 to 25 μm), complemented with NIRCam photometry from the JADES survey (1 to 5 μm), available HST photometry (0.4 to 0.9 μm), and spectroscopic redshifts from the FRESCO and JADES surveys in GOODS-S for 443 star-forming (without dominant active galactic nucleus (AGN)) galaxies at z = 0.7 − 2.0. This redshift range was chosen to ensure that the MIRI data cover mid-IR dust emission. Our methodology involved employing ultraviolet (UV) to IR energy balance SED fitting to robustly constrain the fraction of dust mass in PAHs and dust-obscured luminosity. Additionally, we inferred dust sizes from MIRI 15 μm imaging data, enhancing our understanding of the physical characteristics of dust within these galaxies. Results. We find a strong correlation between the fraction of dust in PAHs (PAH fraction, qPAH) with stellar mass. Moreover, the sub-sample with robust qPAH measurements (N = 216) shows a similar behaviour between qPAH and gas-phase metallicity to that at z ∼ 0, suggesting a universal relation: qPAH is constant (∼3.4%) above a metallicity of Z ∼ 0.5 Z⊙ and decreases to < 1% at metallicities ≲0.3 Z⊙. This indicates that metallicity is a good indicator of the interstellar medium properties that affect the balance between the formation and destruction of PAHs. The lack of a redshift evolution from z ∼ 0 − 2 also implies that above Z ∼ 0.5 Z⊙ the PAH emission effectively traces obscured luminosity and the previous locally calibrated PAH-SFR calibrations remain applicable in this metallicity regime. We observe a strong correlation between the obscured UV luminosity fraction (ratio of obscured to total luminosity) and stellar mass. Above the stellar mass of M* > 5 × 109 M⊙, on average, more than half of the emitted luminosity is obscured, while there exists a non-negligible population of lower-mass galaxies with > 50% obscured fractions. At a fixed mass, the obscured fraction correlates with SFR surface density. This is a result of higher dust covering fractions in galaxies with more compact star-forming regions. Similarly, galaxies with high IRX (IR to UV luminosity) at a given mass or UV continuum slope (β) tend to have higher ΣSFR and shallower attenuation curves, owing to their higher effective dust optical depths and more compact star-forming regions.}, author = {Shivaei, Irene and Alberts, Stacey and Florian, Michael and Rieke, George and Wuyts, Stijn and Bodansky, Sarah and Bunker, Andrew J. and Cameron, Alex J. and Curti, Mirko and Da'Eugenio, Francesco and Dudzevičiūte, Ugne and Ji, Zhiyuan and Johnson, Benjamin D. and Kramarenko, Ivan and Lyu, Jianwei and Matthee, Jorryt J and Morrison, Jane and Naidu, Rohan and Pérez-González, Pablo G. and Reddy, Naveen and Robertson, Brant and Sun, Yang and Tacchella, Sandro and Whitaker, Katherine and Williams, Christina C. and Willmer, Christopher N.A. and Witstok, Joris and Xiao, Mengyuan and Zhu, Yongda}, issn = {1432-0746}, journal = {Astronomy and Astrophysics}, publisher = {EDP Sciences}, title = {{A new census of dust and polycyclic aromatic hydrocarbons at z = 0.7-2 with JWST MIRI}}, doi = {10.1051/0004-6361/202449579}, volume = {690}, year = {2024}, } @article{18493, abstract = {Context. The escape of Lyman-α photons at redshifts greater than two is an ongoing subject of study and an important quantity to further understanding of Lyman-α emitters (LAEs), the transmission of Lyman-α photons through the interstellar medium and intergalactic medium, and the impact these LAEs have on cosmic reionisation. Aims. This study aims to assess the Lyman-α escape fraction, fesc, Lyα, over the redshift range 2.9 < z < 6.7, focusing on Very Large Telescope/Multi Unit Spectroscopic Explorer (VLT/MUSE) selected, gravitationally lensed, intrinsically faint LAEs. These galaxies are of particular interest as the potential drivers of cosmic reionisation. Methods. We assessed fesc, Lyα in two ways: through an individual study of 96 LAEs behind the A2744 lensing cluster, with James Webb Space Telescope/Near-Infrared Camera (JWST/NIRCam) and HST data, and through a study of the global evolution of fesc, Lyα using the state-of-the-art luminosity functions for LAEs and the UV-selected ‘parent’ population (dust-corrected). We compared these studies to those in the literature based on brighter samples. Results. We find a negligible redshift evolution of fesc, Lyα for our individual galaxies; it is likely that it was washed out by significant intrinsic scatter. We observed a more significant evolution towards higher escape fractions with decreasing UV magnitude and fit this relation. When comparing the two luminosity functions to derive fesc, Lyα in a global sense, we saw agreement with previous literature when integrating the luminosity functions to a bright limit. However, when integrating using a faint limit equivalent to the observational limits of our samples, we observed enhanced values of fesc, Lyα, particularly around z ∼ 6, where fesc, Lyα becomes consistent with 100% escape. This indicates for the faint regimes we sampled that galaxies towards reionisation tend to allow very large fractions of Lyman-α photons to escape. We interpret this as evidence of a lack of any significant dust in these populations; our sample is likely dominated by young, highly star-forming chemically unevolved galaxies. Finally, we assessed the contribution of the LAE population to reionisation using our latest values for fesc, Lyα and the LAE luminosity density. The dependence on the escape fraction of Lyman continuum photons is strong, but for values similar to those observed recently in z ∼ 3 LAEs and high-redshift analogues, LAEs could provide all the ionising emissivity necessary for reionisation.}, author = {Goovaerts, I. and Thai, T. T. and Pello, R. and Tuan-Anh, P. and Laporte, N. and Matthee, Jorryt J and Nanayakkara, T. and Pharo, J.}, issn = {1432-0746}, journal = {Astronomy and Astrophysics}, publisher = {EDP Sciences}, title = {{Charting the Lyman-α escape fraction in the range 2.9 < z < 6.7 and consequences for the LAE reionisation contribution}}, doi = {10.1051/0004-6361/202451432}, volume = {690}, 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}, } @article{18523, abstract = {Recent observations from the EIGER JWST program have measured for the first time the quasar–galaxy cross-correlation function at z ≈ 6. The autocorrelation function of faint z ≈ 6 quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter haloes. In this work, we interpret these data building upon an empirical quasar population model that has been applied successfully to quasar clustering and demographic measurements at z ≈ 2–4. We use a new, large-volume N-body simulation with more than a trillion particles, FLAMINGO-10k, to model quasars and galaxies simultaneously. We successfully reproduce observations of z ≈ 6 quasars and galaxies (i.e. their clustering properties and luminosity functions), and infer key quantities such as their luminosity–halo mass relation, the mass function of their host haloes, and their duty cycle/occupation fraction. Our key findings are (i) quasars reside on average in ≈ 1012.5 M haloes (corresponding to ≈ 5σ fluctuations in the initial conditions of the linear density field), but the distribution of host halo masses is quite broad; (ii) the duty cycle of (UV-bright) quasar activity is relatively low (≈ 1 per cent); (iii) galaxies (that are bright in [O III]) live in much smaller haloes (≈ 1010.9 M) and have a larger duty cycle (occupation fraction) of ≈ 13 per cent. Finally, we focus on the inferred properties of quasars and present a homogeneous analysis of their evolution with redshift. The picture that emerges reveals a strong evolution of the host halo mass and duty cycle of quasars at z ≈ 2–6, and calls for new investigations of the role of quasar activity across cosmic time.}, author = {Pizzati, Elia and Hennawi, Joseph F. and Schaye, Joop and Schaller, Matthieu and Eilers, Anna Christina and Wang, Feige and Frenk, Carlos S. and Elbers, Willem and Helly, John C. and Mackenzie, Ruari and Matthee, Jorryt J and Bordoloi, Rongmon and Kashino, Daichi and Naidu, Rohan P. and Yue, Minghao}, issn = {1365-2966}, journal = {Monthly Notices of the Royal Astronomical Society}, number = {4}, pages = {3155--3175}, publisher = {Oxford University Press}, title = {{A unified model for the clustering of quasars and galaxies at z ≈ 6}}, doi = {10.1093/mnras/stae2307}, volume = {534}, year = {2024}, } @article{18527, abstract = {Context. Galaxies evolve through a dynamic exchange of material with their immediate surrounding environment, the so-called circumgalactic medium (CGM). Understanding the physics of gas flows and the nature of the CGM is fundamental to studying galaxy evolution, especially at 4 ≤ z ≤ 6 (i.e., after the Epoch of Reionization) when galaxies rapidly assembled their masses and reached their chemical maturity. Galactic outflows are predicted to enrich the CGM with metals, although it has also been suggested that gas stripping in systems undergoing a major merger may play a role. Aims. In this work, we explore the metal enrichment of the medium around merging galaxies at z ∼ 4.5, observed by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). To do so, we study the nature of the [CII] 158 μm emission in the CGM around these systems, using simulations to help disentangle the mechanisms contributing to the CGM metal pollution. Methods. By adopting an updated classification of major merger systems in the ALPINE survey, we selected and analyzed merging galaxies whose components can be spatially and/or spectrally resolved in a robust way. This makes it possible to distinguish between the [CII] emission coming from the single components of the system and that coming from the system as a whole. We also made use of the dustyGadget cosmological simulation to select synthetic analogs of observed galaxies and guide the interpretation of the observational results. Results. We find a large diffuse [CII] envelope (≳20 kpc) embedding all the merging systems, with at least 25% of the total [CII] emission coming from the medium between the galaxies. Using predictions from dustyGadget, we suggest that this emission has a multi-fold nature, with dynamical interactions between galaxies playing a major role in stripping the gas and enriching the medium with heavy elements.}, author = {Di Cesare, Claudia and Ginolfi, M. and Graziani, L. and Schneider, R. and Romano, M. and Popping, G.}, issn = {1432-0746}, journal = {Astronomy and Astrophysics}, publisher = {EDP Sciences}, title = {{Carbon envelopes around merging galaxies at z ~ 4.5}}, doi = {10.1051/0004-6361/202449164}, volume = {690}, year = {2024}, } @article{18584, abstract = {In this paper, we describe the "Medium Bands, Mega Science" JWST Cycle 2 survey (JWST-GO-4111) and demonstrate the power of these data to reveal both the spatially integrated and spatially resolved properties of galaxies from the local Universe to the era of cosmic dawn. Executed in 2023 November, MegaScience obtained ∼30 arcmin2 of deep multiband NIRCam imaging centered on the z ∼ 0.3 A2744 cluster, including 11 medium-band filters and the two shortest-wavelength broadband filters, F070W and F090W. Together, MegaScience and the UNCOVER Cycle 1 treasury program provide a complete set of deep (∼28–30 magAB) images in all NIRCam medium- and broadband filters. This unique data set allows us to precisely constrain photometric redshifts, map stellar populations and dust attenuation for large samples of distant galaxies, and examine the connection between galaxy structures and formation histories. MegaScience also includes ∼17 arcmin2 of NIRISS parallel imaging in two broadband and four medium-band filters from 0.9 to 4.8 μm, expanding the footprint where robust spectral energy distribution (SED) fitting is possible. We provide example SEDs and multiband cutouts at a variety of redshifts, and use a catalog of JWST spectroscopic redshifts to show that MegaScience improves both the scatter and catastrophic outlier rate of photometric redshifts by factors of 2–3. Additionally, we demonstrate the spatially resolved science enabled by MegaScience by presenting maps of the [O iii] line emission and continuum emission in three spectroscopically confirmed z > 6 galaxies. We show that line emission in reionization-era galaxies can be clumpy, extended, and spatially offset from continuum emission, implying that galaxy assembly histories are complex even at these early epochs. We publicly release fully reduced mosaics and photometric catalogs for both the NIRCam primary and NIRISS parallel fields (jwst-uncover.github.io/megascience).}, author = {Suess, Katherine A. and Weaver, John R. and Price, Sedona H. and Pan, Richard and Wang, Bingjie and Bezanson, Rachel and Brammer, Gabriel and Cutler, Sam E. and Labbé, Ivo and Leja, Joel and Williams, Christina C. and Whitaker, Katherine E. and Atek, Hakim and Dayal, Pratika and De Graaff, Anna and Feldmann, Robert and Franx, Marijn and Fudamoto, Yoshinobu and Fujimoto, Seiji and Furtak, Lukas J. and Goulding, Andy D. and Greene, Jenny E. and Khullar, Gourav and Kokorev, Vasily and Kriek, Mariska and Lorenz, Brian and Marchesini, Danilo and Maseda, Michael V. and Matthee, Jorryt J and Miller, Tim B. and Mitsuhashi, Ikki and Mowla, Lamiya A. and Muzzin, Adam and Naidu, Rohan P. and Nanayakkara, Themiya and Nelson, Erica J. and Oesch, Pascal A. and Setton, David J. and Shipley, Heath and Smit, Renske and Spilker, Justin S. and Van Dokkum, Pieter and Zitrin, Adi}, issn = {1538-4357}, journal = {Astrophysical Journal}, number = {1}, publisher = {IOP Publishing}, title = {{Medium bands, mega science: A JWST/NIRCam medium-band imaging survey of A2744}}, doi = {10.3847/1538-4357/ad75fe}, volume = {976}, year = {2024}, } @article{18585, abstract = {We present the census of Hβ+[OIII] 4960,5008 Åemitters at 6.87 galaxies. We report a rapid decline of the [OIII] luminosity density at z≳6−7 which cannot be explained by the evolution of the cosmic star-formation rate density. Finally we find that FRESCO detects in only 2h galaxies likely accounting for ∼10−20% of the ionising budget at z=7−8 (assuming an escape fraction of 10%), raising the prospect of directly detecting a significant fraction of the sources of reionisation with JWST.}, author = {Meyer, R. A. and Oesch, P. A. and Giovinazzo, E. and Weibel, A. and Brammer, G. and Matthee, Jorryt J and Naidu, R. P. and Bouwens, R. J. and Chisholm, J. and Covelo-Paz, A. and Fudamoto, Y. and Maseda, M. and Nelson, E. and Shivaei, I. and Xiao, M. and Herard-Demanche, T. and Illingworth, G. D. and Kerutt, J. and Kramarenko, Ivan and Labbe, I. and Leonova, E. and Magee, D. and Matharu, J. and Prieto Lyon, G. and Reddy, N. and Schaerer, D. and Shapley, A. and Stefanon, M. and Wozniak, M. A. and Wuyts, S.}, issn = {1365-2966}, journal = {Monthly Notices of the Royal Astronomical Society}, number = {1}, pages = {1067--1094}, publisher = {Oxford University Press}, title = {{JWST FRESCO: A comprehensive census of H β + [O iii] emitters at 6.8 < z < 9.0 in the GOODS fields}}, doi = {10.1093/mnras/stae2353}, volume = {535}, year = {2024}, } @article{18760, abstract = {With the remarkable sensitivity and resolution of JWST in the infrared, measuring rest-optical kinematics of galaxies at z > 5 has become possible for the first time. This study pilots a new method for measuring galaxy dynamics for highly multiplexed, unbiased samples by combining FRESCO NIRCam grism spectroscopy and JADES medium-band imaging. Here we present one of the first JWST kinematic measurements for a galaxy at z > 5. We find a significant velocity gradient, which, if interpreted as rotation, yields Vrot = 305 ± 70 km s−1, and we hence refer to this galaxy as Twister-z5. With a rest-frame optical effective radius of re = 2.25 kpc, the high rotation velocity in this galaxy is not due to a compact size, as may be expected in the early Universe, but rather to a high total mass, (math formula). This is a factor of roughly 10× higher than the stellar mass within re. We also observe that the radial Hα equivalent width profile and the specific star formation rate map from resolved stellar population modeling are centrally depressed by a factor of ∼1.5 from the center to re. Combined with the morphology of the line-emitting gas in comparison to the continuum, this centrally suppressed star formation is consistent with a star-forming disk surrounding a bulge growing inside out. While large, rapidly rotating disks are common to z ∼ 2, the existence of one after only 1 Gyr of cosmic time, shown for the first time in ionized gas, adds to the growing evidence that some galaxies matured earlier than expected in the history of the Universe.}, author = {Nelson, Erica and Brammer, Gabriel and Giménez-Arteaga, Clara and Oesch, Pascal A. and Naidu, Rohan P. and Übler, Hannah and Matharu, Jasleen and Shapley, Alice E. and Whitaker, Katherine E. and Wisnioski, Emily and Förster Schreiber, Natascha M. and Smit, Renske and Van Dokkum, Pieter and Chisholm, John and Endsley, Ryan and Hartley, Abigail I. and Gibson, Justus and Giovinazzo, Emma and Illingworth, Garth and Labbe, Ivo and Maseda, Michael V. and Matthee, Jorryt J and Covelo Paz, Alba and Price, Sedona H. and Reddy, Naveen A. and Shivaei, Irene and Weibel, Andrea and Wuyts, Stijn and Xiao, Mengyuan and Alberts, Stacey and Baker, William M. and Bunker, Andrew J. and Cameron, Alex J. and Charlot, Stephane and Eisenstein, Daniel J. and De Graaff, Anna and Ji, Zhiyuan and Johnson, Benjamin D. and Jones, Gareth C. and Maiolino, Roberto and Robertson, Brant and Sandles, Lester and Suess, Katherine A. and Tacchella, Sandro and Williams, Christina C. and Witstok, Joris}, issn = {2041-8213}, journal = {Astrophysical Journal Letters}, number = {2}, publisher = {IOP Publishing}, title = {{Ionized gas kinematics with FRESCO: An extended, massive, rapidly rotating galaxy at z = 5.4}}, doi = {10.3847/2041-8213/ad7b17}, volume = {976}, year = {2024}, } @article{14852, abstract = {The physical conditions giving rise to high escape fractions of ionizing radiation (LyC fesc) in star-forming galaxies – most likely protagonists of cosmic reionization – are not yet fully understood. Using the VLT/MUSE observations of ∼1400 Ly α emitters at 2.9 &lt; z &lt; 6.7, we compare stacked rest-frame UV spectra of candidates for LyC leakers and non-leakers selected based on their Ly α profiles. We find that the stacks of potential LyC leakers, i.e. galaxies with narrow, symmetric Ly α profiles with small peak separation, generally show (i) strong nebular O iii]λ1666, [Si iii]λ1883, and [C iii]λ1907 +C iii]λ1909 emission, indicating a high-ionization state of the interstellar medium (ISM); (ii) high equivalent widths of He iiλ1640 (∼1 − 3 Å), suggesting the presence of hard ionizing radiation fields; (iii) Si ii*λ1533 emission, revealing substantial amounts of neutral hydrogen off the line of sight; (iv) high C ivλλ1548,1550 to [C iii]λ1907 +C iii]λ1909 ratios (C iv/C iii] ≳0.75) , signalling the presence of low column density channels in the ISM. In contrast, the stacks with broad, asymmetric Ly α profiles with large peak separation show weak nebular emission lines, low He iiλ1640 equivalent widths (≲1 Å), and low C iv/C iii] (≲0.25), implying low-ionization states and high-neutral hydrogen column densities. Our results suggest that C iv/C iii] might be sensitive to the physical conditions that govern LyC photon escape, providing a promising tool for identification of ionizing sources among star-forming galaxies in the epoch of reionization.}, author = {Kramarenko, Ivan and Kerutt, J and Verhamme, A and Oesch, P A and Barrufet, L and Matthee, Jorryt J and Kusakabe, H and Goovaerts, I and Thai, T T}, issn = {1365-2966}, journal = {Monthly Notices of the Royal Astronomical Society}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {4}, pages = {9853--9871}, publisher = {Oxford University Press}, title = {{Linking UV spectral properties of MUSE Ly α emitters at z ≳ 3 to Lyman continuum escape}}, doi = {10.1093/mnras/stad3853}, volume = {527}, year = {2024}, } @article{15170, abstract = {The James Webb Space Telescope is revealing a new population of dust-reddened broad-line active galactic nuclei (AGN) at redshifts z ≳ 5. Here we present deep NIRSpec/Prism spectroscopy from the Cycle 1 Treasury program Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) of 15 AGN candidates selected to be compact, with red continua in the rest-frame optical but with blue slopes in the UV. From NIRCam photometry alone, they could have been dominated by dusty star formation or an AGN. Here we show that the majority of the compact red sources in UNCOVER are dust-reddened AGN: 60% show definitive evidence for broad-line Hα with a FWHM > 2000 km s −1, 20% of the current data are inconclusive, and 20% are brown dwarf stars. We propose an updated photometric criterion to select red z > 5 AGN that excludes brown dwarfs and is expected to yield >80% AGN. Remarkably, among all zphot > 5 galaxies with F277W – F444W > 1 in UNCOVER at least 33% are AGN regardless of compactness, climbing to at least 80% AGN for sources with F277W – F444W > 1.6. The confirmed AGN have black hole masses of 107–109M⊙. While their UV luminosities (−16 > MUV > −20 AB mag) are low compared to UV-selected AGN at these epochs, consistent with percent-level scattered AGN light or low levels of unobscured star formation, the inferred bolometric luminosities are typical of 107–109M⊙ black holes radiating at ∼10%–40% the Eddington limit. The number densities are surprisingly high at ∼10−5 Mpc−3 mag−1, 100 times more common than the faintest UV-selected quasars, while accounting for ∼1% of the UV-selected galaxies. While their UV faintness suggests they may not contribute strongly to reionization, their ubiquity poses challenges to models of black hole growth.}, author = {Greene, Jenny E. and Labbe, Ivo and Goulding, Andy D. and Furtak, Lukas J. and Chemerynska, Iryna and Kokorev, Vasily and Dayal, Pratika and Volonteri, Marta and Williams, Christina C. and Wang, Bingjie and Setton, David J. and Burgasser, Adam J. and Bezanson, Rachel and Atek, Hakim and Brammer, Gabriel and Cutler, Sam E. and Feldmann, Robert and Fujimoto, Seiji and Glazebrook, Karl and De Graaff, Anna and Khullar, Gourav and Leja, Joel and Marchesini, Danilo and Maseda, Michael V. and Matthee, Jorryt J and Miller, Tim B. and Naidu, Rohan P. and Nanayakkara, Themiya and Oesch, Pascal A. and Pan, Richard and Papovich, Casey and Price, Sedona H. and Van Dokkum, Pieter and Weaver, John R. and Whitaker, Katherine E. and Zitrin, Adi}, issn = {1538-4357}, journal = {Astrophysical Journal}, publisher = {IOP Publishing}, title = {{UNCOVER spectroscopy confirms the surprising ubiquity of active galactic nuclei in red sources at z > 5}}, doi = {10.3847/1538-4357/ad1e5f}, volume = {964}, year = {2024}, } @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{15336, abstract = {Submillimeter surveys toward overdense regions in the early Universe are essential for uncovering the obscured star formation and the cold gas content of assembling galaxies within massive dark matter halos. In this work, we present deep ALMA mosaic observations covering an area of ∼2′×2′ around MUSE Quasar Nebula 01 (MQN01), one of the largest and brightest Ly-α emitting nebulae discovered thus far; it surrounds a radio-quiet quasar at z ≃ 3.25. Our observations target the 1.2 and the 3 mm dust continuum as well as the carbon monoxide CO(4–3) transition in galaxies in the vicinity of the quasar. We identify a robust sample of 11 CO-line-emitting galaxies (including a closely separated quasar companion) that lie within ±4000 km s−1 of the quasar systemic redshift. A fraction of these objects were missed in previous deep rest-frame optical/UV surveys, which highlights the critical role of (sub)millimeter imaging. We also detect a total of 11 sources revealed in the dust continuum at 1.2 mm; six of them have either high-fidelity spectroscopic redshift information from rest-frame UV metal absorptions or the CO(4–3) line that places them in the same narrow redshift range. A comparison of the CO luminosity function and 1.2 mm number count density with those of the general fields points to a galaxy overdensity of δ > 10. We find evidence of a systematic flattening at the bright end of the CO luminosity function with respect to the trend measured in blank fields. Our findings reveal that galaxies in dense regions at z ∼ 3 are more massive and significantly richer in molecular gas than galaxies in fields, which enables a faster and accelerated assembly. This is the first in a series of studies aimed at characterizing one of the densest regions of the Universe found so far at z > 3.}, author = {Pensabene, A. and Cantalupo, S. and Cicone, C. and Decarli, R. and Galbiati, M. and Ginolfi, M. and De Beer, S. and Fossati, M. and Fumagalli, M. and Lazeyras, T. and Pezzulli, G. and Travascio, A. and Wang, W. and Matthee, Jorryt J and Maseda, M. V.}, issn = {1432-0746}, journal = {Astronomy and Astrophysics}, publisher = {EDP Sciences}, title = {{ALMA survey of a massive node of the Cosmic Web at z ∼ 3: I. Discovery of a large overdensity of CO emitters}}, doi = {10.1051/0004-6361/202348659}, volume = {684}, year = {2024}, } @article{15405, abstract = {We report JWST/NIRCam measurements of quasar host galaxy emissions and supermassive black hole (SMBH) masses for six quasars at 5.9 < z < 7.1 in the Emission-line galaxies and Intergalactic Gas in the Epoch of Reionization (EIGER) project. We obtain deep NIRCam imaging in the F115W, F200W, and F356W bands, as well as F356W grism spectroscopy of the quasars. We use bright unsaturated stars to construct models of the point-spread functions (PSFs) and estimate the errors of these PSFs. We then measure or constrain the fluxes and morphology of the quasar host galaxies by fitting the quasar images as a point source plus an exponential disk. We successfully detect the host galaxies of three quasars, which have host-to-quasar-flux ratios of ∼1%–5%. Spectral energy distribution fitting suggests that these quasar host galaxies have stellar masses of M* ≳ 1010M⊙. For quasars with host galaxy nondetections, we estimate the upper limits of their stellar masses. We use the grism spectra to measure the Hβ line profile and the continuum luminosity, then estimate the SMBH masses for the quasars. Our results indicate that the positive relation between SMBH masses and host galaxy stellar masses already exists at redshift z ≳ 6. The quasars in our sample show a high BH-to-stellar-mass ratio of MBH/M* ∼ 0.15, which is about ∼2 dex higher than local relations. We find that selection effects only contribute partially to the high MBH/M* ratios of high-redshift quasars. This result hints at a possible redshift evolution of the MBH–M* relation.}, author = {Yue, Minghao and Eilers, Anna Christina and Simcoe, Robert A. and Mackenzie, Ruari and Matthee, Jorryt J and Kashino, Daichi and Bordoloi, Rongmon and Lilly, Simon J. and Naidu, Rohan P.}, issn = {1538-4357}, journal = {Astrophysical Journal}, number = {2}, publisher = {IOP Publishing}, title = {{EIGER. V. Characterizing the host galaxies of luminous quasars at z ≳ 6}}, doi = {10.3847/1538-4357/ad3914}, volume = {966}, year = {2024}, } @article{21064, abstract = {The James Webb Space Telescope (JWST) is revolutionizing our knowledge of z > 5 galaxies and their actively accreting black holes. Using the JWST Cycle 1 Treasury program Ultradeep NIRSpec and NIRCam Observations before the Epoch of Reionization (UNCOVER) in the lensing field A2744, we report the identification of a sample of little red dots at 3 < zphot < 7 that likely contain highly reddened accreting supermassive black holes. Using a NIRCam-only selection to F444W < 27.7 mag, we find 26 sources over the ∼45 arcmin^2 field that are blue in F115W − F200W ∼ 0 (or βUV ∼ –2.0 for fλ ∝ λ^β), red in F200W − F444W = 1−4 (βopt ∼ +2.0), and are dominated by a point-source-like central component. Of the 20 sources with deep Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm coverage, none are detected individually or in a stack. For the majority of the sample, spectral energy distribution fits to the JWST+ALMA observations prefer models with hot dust rather than obscured star formation to reproduce the red NIRCam colors and ALMA 1.2 mm nondetections. While compact dusty star formation cannot be ruled out, the combination of extremely small sizes (〈re〉 ≈ 50 pc after correction for magnification), red rest-frame optical slopes, and hot dust can be explained by reddened broad-line active galactic nuclei (AGNs). Our targets have faint M1450 ≈ −14 to −18 mag but inferred bolometric luminosities of Lbol = 10^43–10^46 erg s^−1, reflecting their obscured nature. If the candidates are confirmed as AGNs with upcoming UNCOVER spectroscopy, then we have found an abundant population of reddened luminous AGNs that are at least ten times more numerous than UV-luminous AGNs at the same intrinsic bolometric luminosity.}, author = {Labbe, Ivo and Greene, Jenny E. and Bezanson, Rachel and Fujimoto, Seiji and Furtak, Lukas J. and Goulding, Andy D. and Matthee, Jorryt J and Naidu, Rohan P. and Oesch, Pascal A. and Atek, Hakim and Brammer, Gabriel and Chemerynska, Iryna and Coe, Dan and Cutler, Sam E. and Dayal, Pratika and Feldmann, Robert and Franx, Marijn and Glazebrook, Karl and Leja, Joel and Maseda, Michael and Marchesini, Danilo and Nanayakkara, Themiya and Nelson, Erica J. and Pan, Richard and Papovich, Casey and Price, Sedona H. and Suess, Katherine A. 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}, publisher = {IOP Publishing}, title = {{UNCOVER: Candidate red active galactic nuclei at 3 < z < 7 with JWST and ALMA}}, doi = {10.3847/1538-4357/ad3551}, volume = {978}, year = {2024}, }