Showing 1–50 of 72 results for author: Annunziatella, M

MIDIS: Strong H$β$+[OIII] Line Emitters at $z \geq 9$

Authors: Thomas R. Greve, Steven Gillman, Pierluigi Rinaldi, Iris Jermann, Jens Melinder, Göran Östlin, Pablo G. Pérez-González, Luis Colina, Fabian Walter, Javier Álvarez-Márquez, Martin J. Ward, Alejandro Crespo Gómez, John P. Pye, Tuomo V. Tikkanen, Edoardo Iani, Roman A. Meyer, Leindert A. Boogaard, Jens Hjorth, Danial Langeroodi, Paul van der Werf, Sarah E. I. Bosman, Karina I. Caputi, Luca Costantin, Marianna Annunziatella, Arjan Bik , et al. (2 additional authors not shown)

Abstract: We present a search for strong H$β$+[O III] emitters at $z=9.4-11.3$ in the HUDF using ultra-deep JWST/MIRI F560W imaging from the MIDIS survey. Three galaxies are identified via pronounced F560W flux excesses, consistent with strong rest-frame optical line emission. SED modelling yields rest-frame H$β$+[O III] equivalent widths of $\sim 600-1300$AA (median $\simeq 1260$AA), placing these sources… ▽ More We present a search for strong H$β$+[O III] emitters at $z=9.4-11.3$ in the HUDF using ultra-deep JWST/MIRI F560W imaging from the MIDIS survey. Three galaxies are identified via pronounced F560W flux excesses, consistent with strong rest-frame optical line emission. SED modelling yields rest-frame H$β$+[O III] equivalent widths of $\sim 600-1300$AA (median $\simeq 1260$AA), placing these sources among the most extreme known at these epochs. Combining these with a literature sample of 16 spectroscopically confirmed galaxies at $z\geq 9$, we find a median ${\rm EW}^{\rm Hβ+[O III]}_{\rm rest}\simeq 1300$AA, similar to values at $z\sim6-9$. We find no evidence for either a strong increase or decline in EW beyond $z\sim9$. A tentative trend of higher EW with increasing UV luminosity is observed, while no statistically significant anti-correlation with stellar mass is found. We place a first constraint on the H$β$+[O III] luminosity function at $z\simeq9-11$ ($Φ\sim10^{-3.4}\,{\rm Mpc^{-3}\,dex^{-1}}$ at $\log( L_{\rm Hβ+[OIII]}/{\rm erg\,s^{-1}})=42.5$), consistent with a decline relative to $z\sim7-8$. The MIDIS sources have $\log(ξ_{\rm ion}/{\rm Hz\,erg^{-1}})=25.1-25.4$. We find significant correlations between $ξ_{\rm ion}$ and EW and $β$, but not with UV luminosity, consistent with trends at lower redshift. These results suggest that the physical conditions governing nebular emission and ionising efficiency are already in place at $z\sim9-11$, extending trends established at $z\sim6-9$. △ Less

Submitted 5 April, 2026; originally announced April 2026.

Comments: 19 pages, 12 figures. Submitted to MNRAS

JWST Reveals Two Overmassive Black Hole Candidates in Dwarf Galaxies at z $\approx$ 0.7: Pushing Black Hole Searches into the Dwarf-Galaxy Regime

Authors: E. Iani, P. Rinaldi, A. Torralba, J. Lyu, R. Navarro-Carrera, G. H. Rieke, F. Sun, C. Willott, Y. Zhu, A. Alonso-Herrero, M. Annunziatella, P. Bergamini, K. Caputi, M. Catone, L. Colina, R. Cooper, L. Costantin, A. Crespo Gómez, G. Desprez, C. Di Cesare, M. J. Hayes, I. Jermann, G. Kotiwale, I. Kramarenko, D. Langeroodi , et al. (13 additional authors not shown)

Abstract: We report the discovery and characterization of two compact galaxies, Pelias and Neleus, at z ~ 0.71 and z ~ 0.75, identified in MACS J0416.1-2403 and GOODS-North. Both exhibit unusual spectral energy distributions (SEDs), with very blue rest-frame UV-optical emission and a steep rise toward near- and mid-infrared wavelengths. JWST/NIRISS and JWST/NIRSpec spectroscopy show strong rest-frame optica… ▽ More We report the discovery and characterization of two compact galaxies, Pelias and Neleus, at z ~ 0.71 and z ~ 0.75, identified in MACS J0416.1-2403 and GOODS-North. Both exhibit unusual spectral energy distributions (SEDs), with very blue rest-frame UV-optical emission and a steep rise toward near- and mid-infrared wavelengths. JWST/NIRISS and JWST/NIRSpec spectroscopy show strong rest-frame optical lines ([O III] 4959,5007 and Halpha) with extreme equivalent widths (>= 1000 Angstrom), indicating young burst-dominated populations with low metallicities (Z ~ 0.1-0.4 Zsun), low dust attenuation (Av ~ 0.2 mag), and stellar masses of Mstar ~ 10^7 Msun. Nonetheless, JWST/MIRI photometry reveals a strong mid-infrared excess that cannot be explained by stellar populations or star-formation-heated dust alone, requiring a hot-dust component most naturally associated with a deeply embedded active galactic nucleus (AGN). SED modelling yields log10(Lbol [erg/s]) ~ 43.7-44.0, implying black hole masses of log10(MBH [Msun]) ~ 5.7-6.7 under the assumption of Eddington-limited accretion. Given the very low stellar masses of the hosts, this corresponds to black-hole-to-stellar mass ratios of about 6-60%, well above the extrapolation of local scaling relations. The lack of X-ray detections suggests that the accretion may be either heavily obscured or intrinsically X-ray weak. Their SEDs also resemble those of Blue Excess Hot Dust Obscured Galaxies and show the characteristic V-shaped continuum seen in Little Red Dots, although with the inflection occurring at redder wavelengths. △ Less

Submitted 18 March, 2026; originally announced March 2026.

Comments: 16 pages (main text), 9 figures (+2 in the appendix), 2 tables (+6 in the appendix)

Radio selection of heavily obscured AGN in the J1030 field: unraveling a missing Compton-thick population

Authors: Giovanni Mazzolari, Roberto Gilli, Marco Mignoli, Marcella Brusa, Isabella Prandoni, Fabio Vito, Ivan Delvecchio, Giorgio Lanzuisi, Alessandro Peca, Andrea Comastri, Stefano Marchesi, Marco Chiaberge, Marisa Brienza, Cristian Vignali, Matilde Signorini, Quirino D'Amato, Fabrizio Gentile, Kazushi Iwasawa, Colin Norman, Alberto Traina, Federica Loiacono, Pietro Baldini, Marianna Annunziatella, Roberto Decarli

Abstract: We tested the effectiveness of radio selection to discover heavily obscured AGNs, particularly at high-z, and we measured their abundance for the first time from a radio perspective. We consider the radio sources detected in the J1030 field, which is one of the fields with the deepest combination of 1.4 GHz radio and X-ray observations. We defined a radio excess parameter as the ratio between the… ▽ More We tested the effectiveness of radio selection to discover heavily obscured AGNs, particularly at high-z, and we measured their abundance for the first time from a radio perspective. We consider the radio sources detected in the J1030 field, which is one of the fields with the deepest combination of 1.4 GHz radio and X-ray observations. We defined a radio excess parameter as the ratio between the star formation rate (SFR) that would correspond to the observed radio luminosity and the one directly derived from the spectral energy distribution (SED) fitting, $\rm REX=SFR_{1.4GHz}/SFR^{corr}_{SED}$. We then select as radio excess AGN those sources with $\rm REX>8.5$, corresponding to a $3σ$ excess above the median value. In this way, we find 145 radio-excess sources falling into the \textit{Chandra} X-ray image footprint but without X-ray detection. From the deep X-ray upper limits, we estimated a lower limit to the obscuration of each radio-excess AGN, finding on average $\log (N_H/\rm{cm^{-2}})>23.7$. A CTK AGN scenario is also supported by the results of the X-ray stacking analysis performed on sources at $z>1.5$, which revealed X-ray luminosities and hardness ratios compatible with very highly obscured AGN. Finally, we computed the number density of these radio-selected CTK AGN. While at $z\sim 2$ the radio number density agrees well with the CTK AGN predictions of different population synthesis models, at $z\sim3$ the radio selection returns a CTK AGN number density $\sim 2-3$ times larger than what is predicted by the X-ray models and observations. This result supports the effectiveness of radio emission in selecting the most obscured sources, unraveling a population of AGN potentially missed by X-rays surveys at $z>3$, paving the way to a synergistic use of the future radio and X-ray facilities such as the \textit{SKAO} and \textit{NewAthena}. △ Less

Submitted 10 March, 2026; originally announced March 2026.

Comments: Accepted for publication on A&A. 15 pages + appendix. Main figures: Fig. 3 (radio-excess selection), Fig. 5 (column density distribution of radio-excess sources non detected in the X-rays), Fig. 8 (radio-selected Compton-thick AGN number density)

There Is More to Outshining: 2D Dust Effects on Stellar Mass Estimates at $3 \leq z < 9$ with JWST in the JADES Field

Authors: M. Hamed, P. G. Pérez-González, M. Annunziatella, L. Colina, I. Shivaei, M. Perna, A. J. Bunker, K. Małek, S. Arribas, J. Álvarez-Márquez, C. N. A. Willmer, H. Übler, R. Bhatawdekar, J. Chevallard, E. Curtis-Lake, Z. Ji, P. Rinaldi, C. C. Williams

Abstract: Dust attenuation modifies the observed spectral energy distribution (SED), leading to biases in the properties inferred from integrated SED fitting. As spatially resolved SED modeling becomes feasible for large high-redshift samples, it is increasingly important to assess how dust attenuation affects resolved mass estimates. We evaluate the impact of dust attenuation on stellar mass estimates deri… ▽ More Dust attenuation modifies the observed spectral energy distribution (SED), leading to biases in the properties inferred from integrated SED fitting. As spatially resolved SED modeling becomes feasible for large high-redshift samples, it is increasingly important to assess how dust attenuation affects resolved mass estimates. We evaluate the impact of dust attenuation on stellar mass estimates derived from integrating spatially resolved SED fitting results. We perform spatially resolved and integrated SED fitting on a sample of 3408 galaxies at $3 \leq z < 9$ from the GOODS South field, combining deep NIRCam from the JWST Advanced Deep Extragalactic Survey (JADES) and HST/ACS imaging from GOODS and CANDELS. We compare galaxy-integrated properties derived from fitting the summed SED with those obtained from spatially resolved SED modeling. Using a two-component dust attenuation model with a variable slope, we investigate how the dust attenuation slope, A(V), and stellar population properties contribute to discrepancies in the resulting stellar mass estimates. Resolved stellar masses are systematically higher than integrated estimates, with a median offset of +0.24 dex. Resolved analyses recover higher dust attenuations ($ΔA(V)\approx +0.08$ mag), lower birth cloud fractions ($Δμ\approx -0.28$), and grayer attenuation curves ($Δδ_{\mathrm{ISM}} = +0.08$), arising from preferential sampling of compact star-forming regions. Integrated fits underestimate stellar ages by $\sim23\%$ at $z < 5$ and 31$\%$ at $z \gtrsim 5$. The stellar mass offset correlates strongly with the age difference and the attenuation slope difference, indicating that age-dependent outshining and spatially varying dust geometry are primary drivers of the discrepancy between resolved and integrated stellar masses. △ Less

Submitted 22 January, 2026; originally announced January 2026.

Comments: Submitted

MIRI spectrophotometry of GN-z11: Detection and nature of an optical red continuum component

Authors: A. Crespo Gómez, L. Colina, P. G. Pérez-González, J. Álvarez-Márquez, M. García-Marín, A. Alonso-Herrero, M. Annunziatella, A. Bik, S. Bosman, A. J. Bunker, A. Labiano, D. Langeroodi, P. Rinaldi, G. Östlin, L. Boogaard, S. Gillman, G. Barro, S. L. Finkelstein, G. C. K. Leung

Abstract: We present new MIRI F560W, F770W and F1000W imaging of GN-z11, extending the previous rest-frame coverage from 0.38 to 0.86$μ$m. We report significant detections (14$σ$) in the F560W and F770W images, and a marginal detection (3.2$σ$) in F1000W. Here, we analyse its SED combining new MIRI imaging data with archival NIRSpec/Prism and MRS spectroscopy, and NIRCam imaging. The continuum emission show… ▽ More We present new MIRI F560W, F770W and F1000W imaging of GN-z11, extending the previous rest-frame coverage from 0.38 to 0.86$μ$m. We report significant detections (14$σ$) in the F560W and F770W images, and a marginal detection (3.2$σ$) in F1000W. Here, we analyse its SED combining new MIRI imaging data with archival NIRSpec/Prism and MRS spectroscopy, and NIRCam imaging. The continuum emission shows a flat energy distribution, in f$_ν$, up to 0.5$μ$m, compatible with the presence of a mixed stellar population of young (4$\pm$1 Myr) and mature (63$\pm$23 Myr) stars that also account for the [O III], H$β$ and H$α$ emission lines. The continuum at rest-frame 0.66$μ$m shows a 36$\pm$3% flux excess above the predicted flux for a mixed stellar population, pointing to the presence of an additional source contributing at these wavelengths. This excess increases to 91$\pm$28% at rest-frame 0.86$μ$m, although with a large uncertainty due to the marginal detection in the F1000W filter. We consider that hot dust emission in the dusty torus around a type 2 AGN could be responsible for the observed excess. Alternatively, this excess could be due to hot dust emission or to a photoluminiscence dust process (Extended Red Emission, ERE) under the extreme UV radiation field, as observed in some local metal-poor galaxies and in young compact starbursts. The presence of a type 1 AGN is not supported by the observed SED since high-z QSOs contribute at wavelengths above rest-frame 1$μ$m, and an additional ad-hoc red source would be required to explain the observed flux excess at 0.66 and 0.86$μ$m. Additional deep MIRI imaging covering the rest-frame near-IR are needed to confirm the flux detection at 10$μ$m with higher significance, and to discriminate between the different hot dust emission in the extreme starburst and AGN scenarios with MIRI imaging at longer wavelengths. △ Less

Submitted 10 February, 2026; v1 submitted 2 December, 2025; originally announced December 2025.

Comments: 14 pages, 5 figures

Journal ref: A&A 706, A46 (2026)

No Luminous Little Red Dots: A Sharp Cutoff in Their Luminosity Function

Authors: Yilun Ma, Jenny E. Greene, Marta Volonteri, Andy D. Goulding, David J. Setton, Marianna Annunziatella, Eiichi Egami, Xiaohui Fan, Vasily Kokorev, Ivo Labbe, Xiaojing Lin, Danilo Marchesini, Jorryt Matthee, Themiya Nanayakkara, Luke Robbins, Anna Sajina, Marcin Sawicki

Abstract: One of the most surprising results of early James Webb Space Telescope (JWST) observations is the discovery of an abundance of red, compact, broad-line objects dubbed "little red dots" (LRDs) at $z>4$. Their spatial density ($\sim10^{-4}$-$10^{-5}\,\mathrm{cMpc^{-3}}$) is 100 times more abundant than UV-selected quasars at those redshift if one extrapolates the quasar luminosity function (QLF) dow… ▽ More One of the most surprising results of early James Webb Space Telescope (JWST) observations is the discovery of an abundance of red, compact, broad-line objects dubbed "little red dots" (LRDs) at $z>4$. Their spatial density ($\sim10^{-4}$-$10^{-5}\,\mathrm{cMpc^{-3}}$) is 100 times more abundant than UV-selected quasars at those redshift if one extrapolates the quasar luminosity function (QLF) down to the LRD regime. However, whether LRDs dominate black hole accretion at quasar-like luminosities ($L_\mathrm{bol}\gtrsim 10^{45-46}\,\mathrm{erg\,s^{-1}}$) remains unanswered, as probing the bright end of the LRD luminosity function requires a much larger area than those able to be surveyed by JWST. In this work, we present our search for the brightest LRDs ($K<23.7$) at $4.5<z<4.9$ using wide-area multiwavelength imaging surveys from the near-UV to the infrared bands. With over 15 square degrees of sky coverage, we only identify one single LRD candidate at $z_\mathrm{phot}\approx4.6$, which translates into a spatial density of $n(M_{5100}<-23.5)\approx10^{-8}\,\mathrm{cMpc^{-3}}$ -- this is nearly 10 times less abundant than the UV-selected quasars at similar optical luminosity. When combined with the LRD sample identified by JWST at the same redshift range, we find a sharp cutoff in the optical luminosity function at $λL_{5100}\approx2.5\times10^{44}\,\mathrm{erg\,s^{-1}}$, while the QLF turnover occurs at $\gtrsim20$ times higher luminosity. We therefore confirm the exclusively low-luminosity nature of LRDs, ruling out that LRDs are the counter parts of quasars. Furthermore, we speculate that, if the shape of the luminosity function holds up, it points to LRDs being powered by low-mass black holes with a narrow range of Eddington-level accretion rates. △ Less

Submitted 2 September, 2025; originally announced September 2025.

Comments: 14 pages, 6 figures, 3 tables. Submitted to ApJL. Comments are welcome!!

MIDIS: Unveiling the Star Formation History in massive galaxies at $1<z<4.5$ with spectro-photometric analysis

Authors: M. Annunziatella, P. G. Pérez-González, J. Álvarez-Márquez, L. Costantin, E. Iani, Á. Labiano, P. Rinaldi, L. Boogaard, R. A. Meyer, G. Östlin, L. Colina, J. Melinder, I. Jermann, S. Gillman, D. Langeroodi, J. Hjorth, A. Alonso-Herrero, A. Eckart, F. Walter, P. P. van der Werf, A. Bik, F. Peißker, K. I. Caputi, M. García-Marín, G. Wright , et al. (1 additional authors not shown)

Abstract: We investigate the star formation histories (SFHs) of a sample of massive galaxies ($M_\star \geq 10^{10} \, M_\odot$) in the redshift range $1 < z < 4.5$. We analyze spectro-photometric data combining broadband photometry from HST and JWST with low-resolution grism spectroscopy from JWST/NIRISS, obtained as part of the MIDIS (MIRI Deep Imaging Survey) program. SFHs are derived through spectral en… ▽ More We investigate the star formation histories (SFHs) of a sample of massive galaxies ($M_\star \geq 10^{10} \, M_\odot$) in the redshift range $1 < z < 4.5$. We analyze spectro-photometric data combining broadband photometry from HST and JWST with low-resolution grism spectroscopy from JWST/NIRISS, obtained as part of the MIDIS (MIRI Deep Imaging Survey) program. SFHs are derived through spectral energy distribution (SED) fitting using two independent codes, BAGPIPES and Synthesizer, under various SFH assumptions. This approach enables a comprehensive assessment of the biases introduced by different modeling choices. The inclusion of NIRISS spectroscopy, even with its low resolution, significantly improves constraints on key physical parameters, such as the mass-weighted stellar age ($t_M$) and formation redshift ($z_{\mathrm{form}}$), by narrowing their posterior distributions. The massive galaxies in our sample exhibit rapid stellar mass assembly, forming 50\% of their mass between $3 \leq z \leq 9$. The highest inferred formation redshifts are compatible with elevated star formation efficiencies ($ε$) at early epochs. Non-parametric SFHs generally imply an earlier and slower mass assembly compared to parametric forms, highlighting the sensitivity of inferred formation timescales to the chosen SFH model, particularly for galaxies at $z < 2$. Quiescent galaxies are, on average, older ($t_M \sim 1.1$ Gyr) and assembled more rapidly at earlier times than their star-forming counterparts. These findings support the ``downsizing'' scenario, in which more massive and passive systems form earlier and more efficiently. △ Less

Submitted 23 August, 2025; originally announced August 2025.

Comments: accepted for publication in A&A

Journal ref: A&A 702, A224 (2025)

CLASH-VLT: The variance in the velocity anisotropy profiles of galaxy clusters

Authors: A. Biviano, E. A. Maraboli, L. Pizzuti, P. Rosati, A. Mercurio, G. De Lucia, C. Ragone-Figueroa, C. Grillo, G. L. Granato, M. Girardi, B. Sartoris, M. Annunziatella

Abstract: The velocity anisotropy profiles, $β(r)$, of galaxy clusters are directly related to the shape of the orbits of their member galaxies. Knowledge of $β(r)$ is important to understand the assembly process of clusters and the evolutionary processes of their galaxies, and to improve the determination of cluster masses based on cluster kinematics. We determined the $β(r)$ of nine massive clusters at re… ▽ More The velocity anisotropy profiles, $β(r)$, of galaxy clusters are directly related to the shape of the orbits of their member galaxies. Knowledge of $β(r)$ is important to understand the assembly process of clusters and the evolutionary processes of their galaxies, and to improve the determination of cluster masses based on cluster kinematics. We determined the $β(r)$ of nine massive clusters at redshift $0.19 \leq z \leq 0.45$ from the CLASH-VLT data set, with 150 to 950 spectroscopic members each. We selected spectroscopic cluster members with the CLUMPS algorithm calibrated on cosmological simulations. We applied the MAMPOSSt code to the distribution of cluster members in projected phase-space to constrain the cluster mass profile, $M(r)$, using priors derived from a previous gravitational lensing analysis. Given the MAMPOSSt best-fit solution for $M(r)$, we then solved the inversion of the Jeans equation to determine $β(r)$ without assumptions of its functional form. We also ran the DS+ code to identify subclusters and characterize the dynamical status of our clusters. The average $β(r)$ is slightly radial, with the anisotropy increasing from $β\simeq 0.2$ at the cluster center, to $β\simeq 0.5$ at the virial radius. There is substantial variance in the $β(r)$ of the individual clusters, that cannot be entirely accounted for by the observational uncertainties. Clusters of lower mass and with a low concentration per given mass have more tangential $β(r)$. A comparison with cluster-sized halos in cosmological hydrodynamical simulations indicates a very good agreement for the average $β(r)$, but a smaller variance in the profiles than observed. We conclude that massive clusters cannot be characterized by a unique universal $β(r)$ and that the orbital distribution of cluster galaxies carries information on the merging history of the cluster. △ Less

Submitted 6 February, 2026; v1 submitted 7 August, 2025; originally announced August 2025.

Comments: Accepted for publication in A&A on Jan 31, 2026

Journal ref: A&A 707, A153 (2026)

Galaxies at the edges: a complete census of MACS J0416.1-2403 cluster

Authors: Rossella Ragusa, Maurizio D'Addona, Amata Mercurio, Marcella Longhetti, Marisa Girardi, Marianna Annunziatella, Nicolas Estrada, Claudio Grillo, Angela Iovino, Giula Rodighiero, Piero Rosati, Benedetta Vulcani, Giuseppe Angora, Hans Bohringer, Massimo Brescia, Gabriel Bartosch Caminha, Gayoung Chon, Fedor Getman, Aniello Grado, Marco Gullieuszik, Luca Limatola, Alessia Moretti, Luisa Pecoraro

Abstract: Numerous studies have established that the environment influences the physical properties of a galaxy. While gas inflows supply the fuel for SF, high density and temperature conditions suppress SF activity through various quenching processes. Investigations into large scale structures, such as filaments and overdense regions in the cluster outskirts, have focused on the low z. To move to intermedi… ▽ More Numerous studies have established that the environment influences the physical properties of a galaxy. While gas inflows supply the fuel for SF, high density and temperature conditions suppress SF activity through various quenching processes. Investigations into large scale structures, such as filaments and overdense regions in the cluster outskirts, have focused on the low z. To move to intermediate z and explore galaxy pathways combined with environmental effects, it is crucial to join wide field spectroscopy and deep photometry. Our primary objective is to spectroscopically analyze the photometric overdensities observed by Estrada et al.(2023) in the outskirts of massive cluster MACS J0416.1-2403 (z=0.397), interpreted as evidence of ongoing group infall, i.e. the pre processing scenario, and to investigate the behavior of galaxies in the outskirts about their g-r color, Mstar, and local density, emphasizing the influence of the environment on galaxy evolution. We conducted a spectroscopic analysis out to 5.5R200, using the AAOmega spectrograph. The large FoV and depth allowed us to explore galaxies up to the cluster periphery and across a wide Mstar range, reaching down to the limit of dwarf galaxies. Redshifts were obtained through independent but comparable methods: Redrock, EZ, and Redmost. We identified 148 new spectroscopic cluster members from a sample of 1236 objects. We found 81 galaxies located in filamentary and overdense regions, supporting the role of filamentary infall in the cluster mass assembly history. Our analysis revealed that galaxies in high density regions are more massive, redder, and more passive, compared to galaxies in low density regions that appear to be bluer, less massive, and more SF. These findings underscore the significance of environmental effects and the role of pre processing in shaping galaxy properties before cluster infall. △ Less

Submitted 12 June, 2025; originally announced June 2025.

Comments: Accepted for publication in the A&A journal. 16 pages, 9 figures, and 1 table

Journal ref: A&A 701, A38 (2025)

Counting Little Red Dots at $z<4$ with Ground-based Surveys and Spectroscopic Follow-up

Authors: Yilun Ma, Jenny E. Greene, David J. Setton, Andy D. Goulding, Marianna Annunziatella, Xiaohui Fan, Vasily Kokorev, Ivo Labbe, Jiaxuan Li, Xiaojing Lin, Danilo Marchesini, Jorryt Matthee, Luke Robbins, Anna Sajina, Marcin Sawicki, O. Grace Telford

Abstract: Little red dots (LRDs) are a population of red, compact objects discovered by JWST at $z>4$. At $4<z<8$, they are roughly 100 times more abundant than UV-selected quasars. However, their number density is uncertain at $z<4$ due to the small sky coverage and limited blue wavelength coverage of JWST. We present our ground-based search for LRDs at $2\lesssim z\lesssim4$, combining ultra-deep Hyper Su… ▽ More Little red dots (LRDs) are a population of red, compact objects discovered by JWST at $z>4$. At $4<z<8$, they are roughly 100 times more abundant than UV-selected quasars. However, their number density is uncertain at $z<4$ due to the small sky coverage and limited blue wavelength coverage of JWST. We present our ground-based search for LRDs at $2\lesssim z\lesssim4$, combining ultra-deep Hyper Suprime-Cam photometry and various (near-)infrared surveys within a total area of $\sim3.1\,\mathrm{deg^{2}}$. We find that for LRDs with $M_{5500}<-22.5$, their number density declines from $\sim10^{-4.5}\,\mathrm{cMpc^{-3}}$ at $z>4$ to $\sim10^{-5.3}\,\mathrm{cMpc^{-3}}$ at $2.7<z<3.7$ and $\sim10^{-5.7}\,\mathrm{cMpc^{-3}}$ at $1.7<z<2.7$. We also present the Magellan/FIRE spectrum of our first followed-up candidate, DEEP23-z2LRD1 at $z_\mathrm{spec}=2.26$, as a proof of concept for our sample selection. Similar to high-redshift LRDs, the spectrum of DEEP23-z2LRD1 exhibits broad H$α$ emission with $\mathrm{FWHM}\approx2400\,\mathrm{km\,s^{-1}}$ and with nearly symmetric narrow H$α$ absorption. Additionally, DEEP23-z2LRD1 has extremely narrow [OIII] lines with $\mathrm{FWHM}\approx140\,\mathrm{km\,s^{-1}}$, suggesting the presence of an accreting black hole in a low-mass host galaxy. Limited by the angular resolution of ground-based surveys, we emphasize that spectroscopic follow-ups are required to characterize the contamination fraction of this sample and pin down LRD number density at $z<4$. △ Less

Submitted 10 April, 2025; originally announced April 2025.

Comments: 15 pages, 8 figures, submitted to ApJL. Comments are welcome

Effects of Environment on the Size Evolution of Quiescent Galaxies: Comparing Galaxies in Clusters and in the Field at Two Rest-frame Wavelengths

Authors: Angelo George, Ivana Damjanov, Marcin Sawicki, Devin J. Williams, Lingjian Chen, Guillaume Desprez, Marianna Annunziatella, Stéphane Arnouts, Stephen Gwyn, Danilo Marchesini, Thibaud Moutard, Anna Sajina

Abstract: We investigate the impact of environment on quiescent galaxy (QG) size evolution using the CLAUDS+HSC imaging covering 18.6~deg$^2$ in five broad filters ($Ugriz$) and the effective radius of a single-Sérsic fit as a proxy for galaxy size. We estimate sizes in two rest-frame wavelengths -- 3000Å (UV) and 5000Å (optical) -- for $\sim86,000$ massive ($M_*>10^{9.5}$M$_\odot$) field QGs and for… ▽ More We investigate the impact of environment on quiescent galaxy (QG) size evolution using the CLAUDS+HSC imaging covering 18.6~deg$^2$ in five broad filters ($Ugriz$) and the effective radius of a single-Sérsic fit as a proxy for galaxy size. We estimate sizes in two rest-frame wavelengths -- 3000Å (UV) and 5000Å (optical) -- for $\sim86,000$ massive ($M_*>10^{9.5}$M$_\odot$) field QGs and for $1,000$ of their similarly massive counterparts from 47 clusters at $0.1<z<0.85$. We fit the size-mass relation (SMR) for field and cluster QGs in five $Δz=0.15$ redshift bins and use the characteristic size of $M_*=5\times10^{10}$M$_\odot$ QGs (SMR's zero point) to trace the change in galaxy size over cosmic time and in two types of environment. Sizes of QGs are larger in the rest-frame UV than in the rest-frame optical in both clusters and the field, and this difference is more prominent in the field sample. However, QGs in clusters are systematically smaller than the field QGs, and this difference is significantly more pronounced if measured in the rest-frame UV light. Modeling of the redshift evolution in the characteristic QG size as $R_e\varpropto(1+z)^β$ shows that the cluster QGs ($β=-1.02$ in UV and $β= -1.00$ in optical) grow in size as fast as the field QGs ($β=-0.95$ in UV and $-1.22$ in optical). This fast growth of cluster QGs is consistent with size increase driven by the accretion of two subpopulations onto clusters: a) field QGs that are larger than their quiescent counterparts in clusters, and b) environmentally quenched galaxies (newcomers) that are larger than the rest of the quiescent population. △ Less

Submitted 27 June, 2025; v1 submitted 28 March, 2025; originally announced March 2025.

Comments: 24 pages, 11 figures, published in the Astrophysical Journal

Journal ref: pJ, 987, 45, 19pp (2025)

The rise of the galactic empire: luminosity functions at $z\sim17$ and $z\sim25$ estimated with the MIDIS$+$NGDEEP ultra-deep JWST/NIRCam dataset

Authors: Pablo G. Pérez-González, Göran Östlin, Luca Costantin, Jens Melinder, Steven L. Finkelstein, Rachel S. Somerville, Marianna Annunziatella, Javier Álvarez-Márquez, Luis Colina, Avishai Dekel, Mark Dickinson, Henry C. Ferguson, Zhaozhou Li, L. Y. Aaron Yung, Mic B. Bagley, Leindert A. Boogaard, Denis Burgarella, Antonello Calabrò, Karina I. Caputi, Yingjie Cheng, Andreas Eckart, Mauro Giavalisco, Steven Gillman, Thomas R. Greve, Mahmoud Hamed , et al. (17 additional authors not shown)

Abstract: We present a sample of six F200W and three F277W dropout sources identified as $16<z<25$ galaxy candidates using the deepest JWST/NIRCam data to date (5$σ$ depths $\sim31.5$ mag at $\geq2$ $μ$m), provided by the MIRI Deep Imaging Survey (MIDIS) and the Next Generation Deep Extragalactic Exploratory Public survey (NGDEEP). We estimate ultraviolet (UV) luminosity functions and densities at… ▽ More We present a sample of six F200W and three F277W dropout sources identified as $16<z<25$ galaxy candidates using the deepest JWST/NIRCam data to date (5$σ$ depths $\sim31.5$ mag at $\geq2$ $μ$m), provided by the MIRI Deep Imaging Survey (MIDIS) and the Next Generation Deep Extragalactic Exploratory Public survey (NGDEEP). We estimate ultraviolet (UV) luminosity functions and densities at $z\sim17$ and $z\sim25$. The number density of galaxies with absolute magnitudes $-19<M_\mathrm{UV}<-18$ at $z\sim17$ ($z\sim25$) is a factor of 4 (25) smaller than at $z\sim12$; the luminosity density presents a similar evolution. Compared to state-of-the-art galaxy simulations, we find the need for an enhanced UV-photon production at $z=17-25$ in $\mathrm{M}_\mathrm{DM}=10^{8.5-9.5}$ M$_\odot$ dark matter halos, provided by an increase in the star formation efficiency at early times and/or by intense compact starbursts with enhanced emissivity linked to strong burstiness, low or primordial gas metallicities, and/or a top-heavy initial mass function. There are few robust theoretical predictions for the evolution of galaxies above $z\sim20$ in the literature, however, the continuing rapid drop in the halo mass function would predict a more rapid evolution than we observe if photon production efficiencies remained constant. Our $z>16$ candidates present mass-weighted ages around 30 Myr, and attenuations $\mathrm{A(V)}<0.1$ mag. Their average stellar mass is $\mathrm{M}_\bigstar\sim10^{7}\,\mathrm{M}_\odot$, implying a stellar-to-baryon mass fraction around 10% if the emissivity increases with redshift, or significantly higher otherwise. Three candidates present very blue UV spectral slopes ($β\sim-3$) compatible with Pop III young ($\lesssim10$ Myr) stars and/or high escape fractions of ionizing photons; the rest have $β\sim-2.5$ similar to $z=10-12$ samples. △ Less

Submitted 30 September, 2025; v1 submitted 19 March, 2025; originally announced March 2025.

Comments: Published in ApJ

Is there a tilt in the fundamental (hyper)plane?

Authors: M. D'Addona, A. Mercurio, C. Grillo, P. Rosati, G. Granata, G. Angora, M. Annunziatella, P. Bergamini, V. Bozza, G. B. Caminha, A. Gargiulo, F. Getman, M. Girardi, A. Grado, L. Limatola, M. Lombardi, M. Meneghetti, L. Pecoraro, R. Ragusa, L. Tortorelli, E. Vanzella

Abstract: We investigate the fundamental plane (FP) of selected early-type (ETG) member galaxies of the galaxy cluster PLCK G287.0+32.9 ($ z_c = 0.3833 $), exploring also four-dimensional hyperplane extensions. We measure ETGs structural parameters and photometry from Hubble Space Telescope (HST) observations. We use high-quality spectroscopic data from the Multi Unit Spectroscopic Explorer (MUSE) to measur… ▽ More We investigate the fundamental plane (FP) of selected early-type (ETG) member galaxies of the galaxy cluster PLCK G287.0+32.9 ($ z_c = 0.3833 $), exploring also four-dimensional hyperplane extensions. We measure ETGs structural parameters and photometry from Hubble Space Telescope (HST) observations. We use high-quality spectroscopic data from the Multi Unit Spectroscopic Explorer (MUSE) to measure the galaxy central stellar velocity dispersions and stellar population properties. With this data, we construct the FP through a robust fitting procedure and analyze its tilt and scatter. We then introduce two hyperplane extensions, one including the stellar mass ($M^\star$-HP) and another including the stellar over total mass fraction ($f_{\mathrm{e}}^\star$-HP), and compare their coefficients and scatter to those of the FP. The FP of PLCK G287.0+32.9 is found to have best-fit parameter values consistent with those in the literature ($α= 1.2 \pm 0.1$ and $β= -0.75 \pm 0.05$), with a scatter of $0.09$ dex. The ($f_{\mathrm{e}}^\star$-HP) shows no tilt compared to the theoretical plane ($α= 2.1 \pm 0.2$ and $β= -1.12 \pm 0.07$), with a scatter of $0.042$ dex, and the ($M^\star$-HP) reveals an even tighter relation, with a scatter of only $0.023$. Our findings support the idea that the FP is a lower-dimensional projection of a more complex hyperplane and confirm that the variations in the dark matter content contribute significantly to the tilt of the FP. Future studies incorporating larger samples of galaxies and additional physical parameters may further refine our understanding of the FP and its higher-dimensional extensions. △ Less

Submitted 5 September, 2025; v1 submitted 17 March, 2025; originally announced March 2025.

Comments: Published in Astronomy and Astrophysics on 15 August 2025

Journal ref: A&A 700, A179 (2025)

The Growth of Galaxy Stellar Haloes Over $0.2 \leq z \leq 1.1$

Authors: Devin J. Williams, Ivana Damjanov, Marcin Sawicki, Harrison Souchereau, Lingjian Chen, Guillaume Desprez, Angelo George, Marianna Annunziatella, Stéphane Arnouts, Stephen Gwyn, Danilo Marchesini, Anna Sajina

Abstract: Galaxies are predicted to assemble their stellar haloes through the accretion of stellar material from interactions with their cosmic environment. Observations that trace stellar halo buildup probe the processes that drive galaxy size and stellar mass growth. We investigate stellar halo assembly over $0.2 \leq z \leq 1.1$ in a mass-complete ($M_{\star} \geq 10^{9.5}M_{\odot}$) sample of 242,456 st… ▽ More Galaxies are predicted to assemble their stellar haloes through the accretion of stellar material from interactions with their cosmic environment. Observations that trace stellar halo buildup probe the processes that drive galaxy size and stellar mass growth. We investigate stellar halo assembly over $0.2 \leq z \leq 1.1$ in a mass-complete ($M_{\star} \geq 10^{9.5}M_{\odot}$) sample of 242,456 star-forming and 88,421 quiescent galaxies (SFGs and QGs) from the CLAUDS and HSC-SSP surveys. We extract galaxy rest-frame $g$-band surface brightness ($μ_g$) profiles to study faint, extended emission in galaxy outskirts. We examine trends in galaxy assembly by analyzing the median $μ_g$ profiles in different SFG and QG \msS ranges with decreasing redshift and connecting evolution in galaxy $μ_g$ profiles with the underlying stellar mass growth in galaxies. Since $z=1.1$, the majority of evolution in the median $μ_g$ profiles of galaxies ($\sim$64$\%$ in SFGs and $\sim$71$\%$ in QGs) occurs throughout their stellar halo regions (2-10$R_e$). More massive galaxies assemble stellar halo material more rapidly at $0.2 \leq z \leq 1.1$. Over this period, QGs grow a larger fraction of their stellar haloes than SFGs at fixed $M_{\star}$ (factor of $\sim$1.2). Although star formation can account for the stellar halo growth observed in low-mass SFGs ($10^{9.5}M_\odot \leq M_\star < 10^{10.5}M_\odot$), high-mass SFGs ($M_\star \geq 10^{10.5}M_\odot$) and both low- and high-mass QGs require an additional assembly mechanism. Our results suggest accretion via minor mergers drives additional stellar halo growth in these galaxies. The contribution from accretion is larger in more massive galaxies (over $M_{\star} \geq 10^{9.5}M_{\odot}$), and QGs exhibit larger fractional increases to their ex-situ fractions over $0.2 \leq z \leq 1.1$ than SFGs at fixed $M_{\star}$. △ Less

Submitted 29 May, 2025; v1 submitted 4 December, 2024; originally announced December 2024.

Comments: 31 pages, 17 figures, 2 tables. Submitted for publication in ApJ

MIRI Deep Imaging Survey (MIDIS) of the Hubble Ultra Deep Field

Authors: Göran Östlin, Pablo G. Pérez-González, Jens Melinder, Steven Gillman, Edoardo Iani, Luca Costantin, Leindert A. Boogaard, Pierluigi Rinaldi, Luis Colina, Hans Ulrik Nørgaard-Nielsen, Daniel Dicken, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Javier Alvarez-Marquez, Marianna Annunziatella, Arjan Bik, Sarah E. I. Bosman, Karina I. Caputi, Alejandro Crespo Gomez, Andreas Eckart, Macarena Garcia-Marin, Jens Hjorth, Olivier Ilbert, Iris Jermann , et al. (19 additional authors not shown)

Abstract: The recently launched James Webb Space Telescope (JWST) is opening new observing windows on the distant universe. Among JWST's instruments, the Mid Infrared Instrument (MIRI) offers the unique capability of imaging observations at wavelengths $λ> 5μ$m. This enables unique access to the rest frame near infra-red (NIR, $λ\ge 1$\mum) emission from galaxies at redshifts $z>4$ and the visual (… ▽ More The recently launched James Webb Space Telescope (JWST) is opening new observing windows on the distant universe. Among JWST's instruments, the Mid Infrared Instrument (MIRI) offers the unique capability of imaging observations at wavelengths $λ> 5μ$m. This enables unique access to the rest frame near infra-red (NIR, $λ\ge 1$\mum) emission from galaxies at redshifts $z>4$ and the visual ($λ\gtrsim 5000$Å) rest frame for $z>9$. We here report on the guaranteed time observations (GTO) from the MIRI European Consortium, of the Hubble Ultra Deep Field (HUDF), forming the MIRI Deep Imaging Survey (MIDIS), consisting of an on source integration time of $\sim41$ hours in the MIRI/F560W (5.6 $μ$m) filter. To our knowledge, this constitutes the longest single filter exposure obtained with JWST of an extragalactic field as yet. △ Less

Submitted 29 November, 2024; originally announced November 2024.

Comments: submitted to A&A on July 30, 2024

Journal ref: A&A 696, A57 (2025)

MAGAZ3NE: Evidence for Galactic Conformity in $z\gtrsim3$ Protoclusters

Authors: Ian McConachie, Gillian Wilson, Ben Forrest, Z. Cemile Marsan, Adam Muzzin, M. C. Cooper, Marianna Annunziatella, Danilo Marchesini, Percy Gomez, Wenjun Chang, Stephanie M. Urbano Stawinski, Michael McDonald, Tracy Webb, Allison Noble, Brian C. Lemaux, Ekta A. Shah, Priti Staab, Lori M. Lubin, Roy R. Gal

Abstract: We examine the quiescent fractions of massive galaxies in six $z\gtrsim3$ spectroscopically-confirmed protoclusters in the COSMOS field, one of which is newly confirmed and presented here. We report the spectroscopic confirmation of MAGAZ3NE~J100143+023021 at $z=3.122^{+0.007}_{-0.004}$ by the Massive Ancient Galaxies At $z>3$ NEar-infrared (MAGAZ3NE) survey. MAGAZ3NE~J100143+023021 contains a tot… ▽ More We examine the quiescent fractions of massive galaxies in six $z\gtrsim3$ spectroscopically-confirmed protoclusters in the COSMOS field, one of which is newly confirmed and presented here. We report the spectroscopic confirmation of MAGAZ3NE~J100143+023021 at $z=3.122^{+0.007}_{-0.004}$ by the Massive Ancient Galaxies At $z>3$ NEar-infrared (MAGAZ3NE) survey. MAGAZ3NE~J100143+023021 contains a total of 79 protocluster members (28 spectroscopic and 51 photometric). Three spectroscopically-confirmed members are star-forming ultra-massive galaxies ($\log(M_{\star}/{\rm M}_\odot)>11$; UMGs), the most massive of which has $\log(M_{\star}/{\rm M}_\odot)=11.15^{+0.05}_{-0.06}$. Combining Keck/MOSFIRE spectroscopy and the COSMOS2020 photometric catalog, we use a weighted Gaussian kernel density estimator to map the protocluster and measure its total mass $2.25^{+1.55}_{-0.65}\times10^{14}~{\rm M}_{\odot}$ in the dense ``core'' region. For each of the six COSMOS protoclusters, we compare the quiescent fraction to the status of the central UMG as star-forming or quiescent. We observe that galaxies in these protoclusters appear to obey galactic conformity: elevated quiescent fractions are found in protoclusters with $UVJ$ quiescent UMGs and low quiescent fractions are found in protoclusters containing $UVJ$ star-forming UMGs. This correlation of star-formation/quiescence in UMGs and the massive galaxies nearby in these protoclusters is the first evidence for the existence of galactic conformity at $z>3$. Despite disagreements over mechanisms behind conformity at low redshifts, its presence at these early cosmic times would provide strong constraints on the physics proposed to drive galactic conformity. △ Less

Submitted 21 November, 2024; originally announced November 2024.

Comments: 6 figures, 5 tables. Accepted for publication in the Astrophysical Journal

Extreme Ionizing Properties of Metal-Poor, Muv ~ -12 Star Complex in the first Gyr

Authors: E. Vanzella, F. Loiacono, M. Messa, M. Castellano, P. Bergamini, A. Zanella, F. Annibali, B. Sun, M. Dickinson, A. Adamo, F. Calura, M. Ricotti, P. Rosati, M. Meneghetti, C. Grillo, M. Bradac, C. J. Conselice, H. Yan, A. Bolamperti, U. Mestric, R. Gilli, M. Gronke, C. Willott, E. Sani, A. Acebron , et al. (16 additional authors not shown)

Abstract: We report the serendipitous discovery of a faint (M_UV > -12.2), low-metallicity (Z ~ 0.02 Zsun), ionizing source (dubbed T2c) with a spectroscopic redshift of z=6.146. T2c is part of a larger structure amplified by the Hubble Frontier Field galaxy cluster MACSJ0416, and was observed with JWST/NIRSpec IFU. Stacking the short-wavelength NIRCam data reveals no stellar continuum detection down to a m… ▽ More We report the serendipitous discovery of a faint (M_UV > -12.2), low-metallicity (Z ~ 0.02 Zsun), ionizing source (dubbed T2c) with a spectroscopic redshift of z=6.146. T2c is part of a larger structure amplified by the Hubble Frontier Field galaxy cluster MACSJ0416, and was observed with JWST/NIRSpec IFU. Stacking the short-wavelength NIRCam data reveals no stellar continuum detection down to a magnitude limit of m_UV ~ 31.0 (3 sigma). However, prominent Hb, [OIII]4959,5007, and Ha emissions are detected, with equivalent widths exceeding 200A, 800A, and 1300A (3 sigma), respectively. The corresponding intrinsic (magnification-corrected x23 +/- 3) ultraviolet and optical rest-frame magnitudes exceed 34.4 and 33.9 (corresponding to M_uv and M_opt fainter than -12.2 and -12.8, at lambda_rest ~ 2000A and ~5000A, respectively), suggesting a stellar mass lower than a few 10^4 Msun under an instantaneous burst scenario. The inferred ionizing photon production efficiency (xi_ion) is high, xi_ion >~ 26.08(25.86) 3(5)sigma, assuming no dust attenuation and no Lyman continuum leakage, indicating the presence of massive stars despite the low mass of the object. The very poor sampling of the initial mass function at such low mass star-forming complex suggests that the formation of very massive stars might be favored in very low metallicity environments. T2c is surrounded by Balmer and weak oxygen emission on a spatial scale of a few hundred parsecs after correcting for lensing effects. This system resembles an HII region potentially powered by currently undetected, extremely efficient, low-metallicity star complexes or clusters. We propose that massive O-type stars populate this low-mass and metallicity high-redshift satellites, likely caught in an early and short formation phase, contributing to the ionization of the surrounding medium. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: 9 pages, 5 figures, 1 table. Submitted to A&A. Comments are welcome

Journal ref: A&A 691, A251 (2024)

MIDIS. Near-infrared rest-frame morphology of massive galaxies at $3<z<5$ in the Hubble eXtreme Deep Field

Authors: L. Costantin, S. Gillman, L. A. Boogaard, P. G. Pérez-González, E. Iani, P. Rinaldi, J. Melinder, A. Crespo Gómez, L. Colina, T. R. Greve, G. Östlin, G. Wright, A. Alonso-Herrero, J. Álvarez-Márquez, M. Annunziatella, A. Bik., K. I. Caputi, D. Dicken, A. Eckart, J. Hjorth, O. Ilbert, I. Jermann, A. Labiano, D. Langeroodi, F. Peißker , et al. (7 additional authors not shown)

Abstract: Thanks to decades of observations using the Hubble Space Telescope (HST), the structure of galaxies at redshift $z>2$ has been widely studied in the rest-frame ultraviolet regime, which traces recent star formation from young stellar populations. But, we still have little information about the spatial distribution of the older, more evolved stellar populations, constrained by the rest-frame infrar… ▽ More Thanks to decades of observations using the Hubble Space Telescope (HST), the structure of galaxies at redshift $z>2$ has been widely studied in the rest-frame ultraviolet regime, which traces recent star formation from young stellar populations. But, we still have little information about the spatial distribution of the older, more evolved stellar populations, constrained by the rest-frame infrared portion of the galaxies' spectral energy distribution. We present the morphological characterization of a sample of 49 massive galaxies ($\log(M_{\star}/M_{\odot})>9$) at redshift $3<z<5$. The MIRI 5.6~$μ$m imaging allows us to characterize the rest-frame near-infrared structure of galaxies beyond cosmic noon, at higher redshifts than possible with NIRCam, tracing their older and dust-insensitive stellar populations. We derive the non-parametric morphology of galaxies and model the light distribution of galaxies with a single Sérsic component and derive their parametric morphology. We find that at $z>3$ massive galaxies show a smooth distribution of their rest-infrared light, strongly supporting the increasing number of regular disk galaxies already in place at early epochs. On the contrary, the ultraviolet structure obtained from HST/WFC3 and JWST/NIRCam observations at $\sim1.5~μ$m is generally more irregular, catching the most recent episodes of star formation. Importantly, we find a segregation of morphologies across cosmic time, where galaxies at redshift $z>3.75$ show later-type morphologies compared to $z\sim3$ galaxies. These findings suggest a transition phase in galaxy assembly and central mass build-up, which is already taking place at $z\sim3-4$. The combined analysis of NIRCam and MIRI imaging datasets allows us to prove that the rest-frame near-infrared morphology of massive galaxies at cosmic noon is typical of compact disk galaxies with a smooth mass distribution. △ Less

Submitted 9 May, 2025; v1 submitted 28 June, 2024; originally announced July 2024.

Comments: Accepted for publication in A&A. 10 pages, 6+2 figures

Journal ref: A&A 699, A360 (2025)

MIDIS: MIRI uncovers Virgil, the first Little Red Dot with clear detection of its host galaxy at z ~ 6.6

Authors: Edoardo Iani, Pierluigi Rinaldi, Karina I. Caputi, Marianna Annunziatella, Danial Langeroodi, Jens Melinder, Pablo G. Pérez-González, Javier Álvarez-Márquez, Leindert A. Boogaard, Sarah E. I. Bosman, Luca Costantin, Thibaud Moutard, Luis Colina, Göran Östlin, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Arjan Bik, Steven Gillman, Alejandro Crespo Gómez, Jens Hjorth, Sarah Kendrew, Alvaro Labiano, John P. Pye, Tuomo V. Tikkanen , et al. (2 additional authors not shown)

Abstract: We present Virgil, a MIRI extremely red object (MERO) detected with the F1000W filter as part of the MIRI Deep Imaging Survey (MIDIS) observations of the Hubble Ultra Deep Field (HUDF). Virgil is a Lyman-$α$ emitter (LAE) at $z_{spec} = 6.6312\pm 0.0019$ (from VLT/MUSE) with a rest-frame UV-to-optical spectral energy distribution (SED) typical of LAEs at similar redshifts. However, MIRI observatio… ▽ More We present Virgil, a MIRI extremely red object (MERO) detected with the F1000W filter as part of the MIRI Deep Imaging Survey (MIDIS) observations of the Hubble Ultra Deep Field (HUDF). Virgil is a Lyman-$α$ emitter (LAE) at $z_{spec} = 6.6312\pm 0.0019$ (from VLT/MUSE) with a rest-frame UV-to-optical spectral energy distribution (SED) typical of LAEs at similar redshifts. However, MIRI observations reveal an unexpected extremely red color at rest-frame near-infrared wavelengths, $\rm F444W - F1000W = 2.33 \pm 0.06$. Such steep rise in the near-infrared, completely missed without MIRI imaging, is poorly reproduced by models including only stellar populations and hints towards the presence of an Active Galactic Nucleus (AGN). According to the shape of its overall SED, Virgil belongs to the recently discovered population of Little Red Dots (LRDs) but displays an extended rest-frame UV-optical wavelengths morphology following a 2D-Sérsic profile with average index $n = 0.93^{+0.85}_{-0.31}$ and $r_e = 0.49^{+0.05}_{-0.11}$ pkpc. Only at MIRI wavelengths Virgil is unresolved due to the coarser PSF. We also estimate a bolometric luminosity $L_{\rm bol, AGN} = (8.9-11)\times 10^{44}\rm~erg~s^{-1}$ and a supermassive black hole mass $M_{\rm BH} = (7-9)\times 10^6\rm ~ M_\odot$ (if $λ_{\rm Edd} = 1$) in agreement with recently reported values for LRDs. This discovery demonstrates the crucial importance of deep MIRI surveys to find AGN amongst high-$z$ galaxies that otherwise would be completely missed and raises the question of how common Virgil-like objects could be in the early Universe. △ Less

Submitted 21 February, 2025; v1 submitted 26 June, 2024; originally announced June 2024.

Comments: 19 pages, 11 figures, 3 tables. Submitted to ApJ

The emergence of the Star Formation Main Sequence with redshift unfolded by JWST

Authors: P. Rinaldi, R. Navarro-Carrera, K. I. Caputi, E. Iani, G. Ostlin, L. Colina, S. Alberts, J. Alvarez-Marquez, M. Annunziatella, L. Boogaard, L. Costantin, J. Hjorth, D. Langeroodi, J. Melinder, T. Moutard, F. Walter

Abstract: We investigate the correlation between stellar mass (M*) and star formation rate (SFR) across the stellar mass range log10(M*/Msun)~6-11. We consider almost 50,000 star-forming galaxies at z~3-7, leveraging data from COSMOS/SMUVS, JADES/GOODS-SOUTH, and MIDIS/XDF. This is the first study spanning such a wide stellar mass range without relying on gravitational lensing effects. We locate our galaxie… ▽ More We investigate the correlation between stellar mass (M*) and star formation rate (SFR) across the stellar mass range log10(M*/Msun)~6-11. We consider almost 50,000 star-forming galaxies at z~3-7, leveraging data from COSMOS/SMUVS, JADES/GOODS-SOUTH, and MIDIS/XDF. This is the first study spanning such a wide stellar mass range without relying on gravitational lensing effects. We locate our galaxies on the SFR-M* plane to assess how the location of galaxies in the star-formation main sequence (MS) and starburst (SB) region evolves with stellar mass and redshift. We find that the two star-forming modes tend to converge at log10(M*/Msun) < 7, with all galaxies found in the SB mode. However, deeper observations will be instrumental for reaching lower SFRs and Msun to further validate this scenario. By dissecting our galaxy sample in stellar mass and redshift, we show that the emergence of the star-formation MS is stellar-mass dependent: while in galaxies with log10(M*/Msun) > 9 the MS is already well in place at z = 5-7, for galaxies with log10(M*/Msun)~7-8 it only becomes significant at z<4. Overall, our results are in line with previous findings that the SB mode dominates amongst low stellar-mass galaxies. The earlier emergence of the MS for massive galaxies is consistent with galaxy downsizing. △ Less

Submitted 1 February, 2025; v1 submitted 19 June, 2024; originally announced June 2024.

Comments: 16 pages, 5 figures, 1 table. Accepted for publication in ApJ

MAGAZ3NE: Massive, Extremely Dusty Galaxies at $z\sim2$ Lead to Photometric Overestimation of Number Densities of the Most Massive Galaxies at $3<z<4$

Authors: Ben Forrest, M. C. Cooper, Adam Muzzin, Gillian Wilson, Danilo Marchesini, Ian McConachie, Percy Gomez, Marianna Annunziatella, Z. Cemile Marsan, Joey Braspenning, Wenjun Chang, Gabriella de Lucia, Fabio Fontanot, Michaela Hirschmann, Dylan Nelson, Annalisa Pillepich, Joop Schaye, Stephanie M. Urbano Stawinski, Mauro Stefanon, Lizhi Xie

Abstract: We present rest-frame optical spectra from Keck/MOSFIRE and Keck/NIRES of 16 candidate ultramassive galaxies targeted as part of the Massive Ancient Galaxies at $z>3$ Near-Infrared (MAGAZ3NE) Survey. These candidates were selected to have photometric redshifts $3\lesssim z_{\rm phot}<4$, photometric stellar masses log($M$/M$_\odot$)$>11.7$, and well-sampled photometric spectral energy distribution… ▽ More We present rest-frame optical spectra from Keck/MOSFIRE and Keck/NIRES of 16 candidate ultramassive galaxies targeted as part of the Massive Ancient Galaxies at $z>3$ Near-Infrared (MAGAZ3NE) Survey. These candidates were selected to have photometric redshifts $3\lesssim z_{\rm phot}<4$, photometric stellar masses log($M$/M$_\odot$)$>11.7$, and well-sampled photometric spectral energy distributions (SEDs) from the UltraVISTA and VIDEO surveys. In contrast to previous spectroscopic observations of blue star-forming and post-starburst ultramassive galaxies, candidates in this sample have very red SEDs implying significant dust attenuation, old stellar ages, and/or active galactic nuclei (AGN). Of these galaxies, eight are revealed to be heavily dust-obscured $2.0<z<2.7$ galaxies with strong emission lines, some showing broad features indicative of AGN, three are Type I AGN hosts at $z>3$, one is a $z\sim1.2$ dusty galaxy, and four galaxies do not have a confirmed spectroscopic redshift. In fact, none of the sample has |$z_{\rm spec}-z_{\rm phot}$|$<0.5$, suggesting difficulties for photometric redshift programs in fitting similarly red SEDs. The prevalence of these red interloper galaxies suggests that the number densities of high-mass galaxies are overestimated at $z\gtrsim3$ in large photometric surveys, helping to resolve the `impossibly early galaxy problem' and leading to much better agreement with cosmological galaxy simulations. A more complete spectroscopic survey of ultramassive galaxies is required to pin down the uncertainties on their number densities in the early universe. △ Less

Submitted 22 October, 2024; v1 submitted 29 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ

A NIRCam-dark galaxy detected with the MIRI/F1000W filter in the MIDIS/JADES Hubble Ultra Deep Field

Authors: Pablo G. Pérez-González, Pierluigi Rinaldi, Karina I. Caputi, Javier Álvarez-Márquez, Marianna Annunziatella, Danial Langeroodi, Thibaud Moutard, Leindert Boogaard, Edoardo Iani, Jens Melinder, Luca Costantin, Goran Östlin, Luis Colina, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Arjan Bik, Sarah E. I. Bosman, Alejandro Crespo Gómez, Daniel Dicken, Andreas Eckart, Macarena García-Marín, Steven Gillman, Manuel Güdel, Thomas Henning , et al. (10 additional authors not shown)

Abstract: We report the discovery of Cerberus, an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at $S/N\sim6$, with $\mathrm{F1000W}\sim27$ mag, and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey, JADES, fainter than the 30.0-30.5 mag $5σ$ detection limits… ▽ More We report the discovery of Cerberus, an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at $S/N\sim6$, with $\mathrm{F1000W}\sim27$ mag, and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey, JADES, fainter than the 30.0-30.5 mag $5σ$ detection limits in individual bands, as well as in the MIDIS F560W ultra-deep data ($\sim$29 mag, $5σ$). Analyzing the spectral energy distribution built with low-$S/N$ ($<5$) measurements in individual optical-to-mid-infrared filters and higher $S/N$ ($\gtrsim5$) in stacked NIRCam data, we discuss the possible nature of this red NIRCam-dark source using a battery of codes. We discard the possibility of Cerberus being a Solar System body based on the $<0.016$" proper motion in the 1-year apart JADES and MIDIS observations. A sub-stellar Galactic nature is deemed unlikely, given that the Cerberus' relatively flat NIRCam-to-NIRCam and very red NIRCam-to-MIRI flux ratios are not consistent with any brown dwarf model. The extragalactic nature of Cerberus offers 3 possibilities: (1) A $z\sim0.4$ galaxy with strong emission from polycyclic aromatic hydrocarbons; the very low inferred stellar mass, $\mathrm{M}_\star=10^{5-6}$ M$_\odot$, makes this possibility highly improbable. (2) A dusty galaxy at $z\sim4$ with an inferred stellar mass $\mathrm{M}_\star\sim10^{8}$ M$_\odot$. (3) A galaxy with observational properties similar to those of the reddest little red dots discovered around $z\sim7$, but Cerberus lying at $z\sim15$, with the rest-frame optical dominated by emission from a dusty torus or a dusty starburst. △ Less

Submitted 26 May, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL, results unchanged, manuscript shortened to publish as a letter

The powerful lens galaxy cluster PLCK G287.0+32.9 ($θ_E \sim 43''$)

Authors: Maurizio D'Addona, Amata Mercurio, Piero Rosati, Claudio Grillo, Gabriel Caminha, Ana Acebron, Giuseppe Angora, Pietro Bergamini, Valerio Bozza, Giovanni Granata, Marianna Annunziatella, Adriana Gargiulo, Raphael Gobat, Paolo Tozzi, Marisa Girardi, Marco Lombardi, Massimo Meneghetti, Pietro Schipani, Luca Tortorelli, Eros Vanzella

Abstract: We present a new high-precision strong lensing model of PLCK G287.0$+$32.9, a massive lens galaxy cluster at $z=0.383$, with the aim to get an accurate estimation of its effective Einstein radius and total mass distribution. We also present a spectroscopic catalog containing accurate redshift measurements for 490 objects, including multiply-lensed sources and cluster member galaxies. We exploit hi… ▽ More We present a new high-precision strong lensing model of PLCK G287.0$+$32.9, a massive lens galaxy cluster at $z=0.383$, with the aim to get an accurate estimation of its effective Einstein radius and total mass distribution. We also present a spectroscopic catalog containing accurate redshift measurements for 490 objects, including multiply-lensed sources and cluster member galaxies. We exploit high-quality spectroscopic data from three pointings of the VLT Multi Unit Spectroscopic Explorer, covering a central $3~\rm{arcmin}^2$ region of the cluster. We complete the spectroscopic catalog by including redshift measurements from VLT-VIMOS and KECK-DEIMOS. We identify 129 spectroscopic cluster member galaxies, with redshift values $0.360 \leq z \leq 0.405$ and $m_{\rm{F160W}} \leq 21$, and 24 photometric ones identified with a Convolutional Neural Network from ancillary HST imaging. We also identify 114 multiple images from 28 background sources, of which 84 images from 16 sources are new and the remaining ones were identified in previous work. The best-fitting lens model shows a root mean square separation value between the predicted and observed positions of the multiple images of $0.75''$, corresponding to an improvement in reconstructing the observed positions of the multiple images of a factor of $2.5$ with respect to previous models. Using the predictive power of our new lens model we find 3 new multiple images and we confirm the configuration of three systems of multiple images that were not used for the optimization of the model. The derived total mass distribution confirms this cluster to be a very prominent gravitational lens with an effective Einstein $θ_{E} = 43.4'' \pm 0.1''$, that is in agreement with previous estimates and corresponds to a total mass enclosed in the critical curve of $M_E = {3.33}_{-0.07}^{+0.02} \times{ 10^{14} M_\odot}$. △ Less

Submitted 23 February, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: submitted to "Astronomy & Astrophysics" and accepted the on January 29th 2024

Journal ref: A&A Volume 686, June 2024, A4

What is the nature of Little Red Dots and what is not, MIRI SMILES edition

Authors: Pablo G. Pérez-González, Guillermo Barro, George H. Rieke, Jianwei Lyu, Marcia Rieke, Stacey Alberts, Christina Williams, Kevin Hainline, Fengwu Sun, David Puskas, Marianna Annunziatella, William M. Baker, Andrew J. Bunker, Eiichi Egami, Zhiyuan Ji, Benjamin D. Johnson, Brant Robertson, Bruno Rodriguez Del Pino, Wiphu Rujopakarn, Irene Shivaei, Sandro Tacchella, Christopher N. A. Willmer, Chris Willott

Abstract: We study little red dots (LRD) detected by JADES and covered by the SMILES MIRI survey. Our sample contains 31 sources, $\sim70$% detected in the two bluest MIRI bands, 40% in redder filters. The median/quartiles redshifts are $z=6.9_{5.9}^{7.7}$ (55% spectroscopic). We analyze the rest-frame ultraviolet through near/mid-infrared spectral energy distributions of LRDs combining NIRCam and MIRI obse… ▽ More We study little red dots (LRD) detected by JADES and covered by the SMILES MIRI survey. Our sample contains 31 sources, $\sim70$% detected in the two bluest MIRI bands, 40% in redder filters. The median/quartiles redshifts are $z=6.9_{5.9}^{7.7}$ (55% spectroscopic). We analyze the rest-frame ultraviolet through near/mid-infrared spectral energy distributions of LRDs combining NIRCam and MIRI observations, using a variety of modeling techniques that include emission from stars, dust, and (un)obscured active galactic nuclei (AGN). The NIRCam$-$MIRI colors, for $\geq10$ $μ$m, are bluer than direct pure emission from AGN tori; the spectral slope flattens in the rest-frame near-infrared, consistent with a 1.6 $μ$m stellar bump. Both observations imply that stellar emission makes the dominant contribution at these wavelengths, expediting a stellar mass estimation: the median/quartiles are $\log \mathrm{M_\star/M_\odot}=9.4_{9.1}^{9.7}$. The number density of LRDs is $10^{-4.0\pm0.1}$ Mpc$^{-3}$, accounting for $14\pm3$% of the global population of galaxies with similar redshifts and masses. The flat ultraviolet spectral range is dominated by young stars. The rest-frame near/mid-infrared (2-4 $μ$m) spectral slope reveals significant amounts of dust (bolometric stellar attenuation $\sim3-4$ mag) heated by strong radiation fields arising from highly embedded compact sources. Our models imply $<0.4$ kpc heating knots, containing dust-enshrouded OB stars or an AGN producing a similar radiation field, obscured by $\mathrm{A(V)}>10$ mag. We conclude that LRDs are extremely intense and compact starburst galaxies with mass-weighted ages 5-10 Myr, very efficient in producing dust, their global energy output dominated by the direct and dust-recycled emission from OB stars, with some contribution from obscured AGN in the mid-infrared. △ Less

Submitted 26 March, 2024; v1 submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted for publication in ApJ, LRD average models provided in https://tinyurl.com/lrdmodels

The FENIKS Survey: Multi-wavelength Photometric Catalog in the UDS Field, and Catalogs of Photometric Redshifts and Stellar Population Properties

Authors: Kumail Zaidi, Danilo Marchesini, Casey Papovich, Jacqueline Antwi-Danso, Mario Nonino, Marianna Annunziatella, Gabriel Brammer, James Esdaile, Karl Glazebrook, Kartheik Iyer, Ivo Labbé, Z. Cemile Marsan, Adam Muzzin, David A. Wake

Abstract: We present the construction of a deep multi-wavelength PSF-matched photometric catalog in the UDS field following the final UKIDSS UDS DR11 release. The catalog includes photometry in 24 filters, from the MegaCam-uS (0.38 microns) band to the Spitzer-IRAC (8 microns) band, over ~ 0.9 sq. deg. and with a 5-sigma depth of 25.3 AB in the K-band detection image. The catalog, containing approximately 1… ▽ More We present the construction of a deep multi-wavelength PSF-matched photometric catalog in the UDS field following the final UKIDSS UDS DR11 release. The catalog includes photometry in 24 filters, from the MegaCam-uS (0.38 microns) band to the Spitzer-IRAC (8 microns) band, over ~ 0.9 sq. deg. and with a 5-sigma depth of 25.3 AB in the K-band detection image. The catalog, containing approximately 188,564 (136,235) galaxies at 0.2 < z < 8.0 with stellar mass > 10$^{8}$ solar masses and K-band total magnitude K < 25.2 (24.3) AB, enables a range of extragalactic studies. We also provide photometric redshifts, corresponding redshift probability distributions, and rest-frame absolute magnitudes and colors derived using the template-fitting code eazy-py. Photometric redshift errors are less than 3 to 4 percent at z < 4 across the full brightness range in K-band and stellar mass range 10$^{8}$-10$^{12}$ solar masses. Stellar population properties (e.g., stellar mass, star-formation rate, dust extinction) are derived from the modeling of the spectral energy distributions (SEDs) using the codes FAST and Dense Basis. △ Less

Submitted 21 May, 2024; v1 submitted 5 January, 2024; originally announced January 2024.

Comments: 26 pages, 18 figures, accepted in ApJ. Photometric Catalog, catalogs of stellar population properties derived using Eazy and Dense Basis and other high-level products available for public use at the following webpage - https://www.zaidikumail.com/feniks-uds-catalogs (DOI:10.5281/zenodo.11002298)

MIDIS: Unveiling the Role of Strong Ha-emitters during the Epoch of Reionization with JWST

Authors: P. Rinaldi, K. I. Caputi, E. Iani, L. Costantin, S. Gillman, P. G. Perez-Gonzalez, G. Ostlin, L. Colina, T. R. Greve, H. U. Noorgard-Nielsen, G. S. Wright, J. Alvarez-Marquez, A. Eckart, M. Garcia-Marin, J. Hjorth, O. Ilbert, S. Kendrew, A. Labiano, O. Le Fevre, J. Pye, T. Tikkanen, F. Walter, P. van der Werf, M. Ward, M. Annunziatella , et al. (18 additional authors not shown)

Abstract: By using the ultra-deep \textit{JWST}/MIRI image at 5.6 $μm$ in the Hubble eXtreme Deep Field, we constrain the role of strong H$α$-emitters (HAEs) during Cosmic Reionization at $z\simeq7-8$. Our sample of HAEs is comprised of young ($<35\;\rm Myr$) galaxies, except for one single galaxy ($\approx 300\;\rm Myr$), with low stellar masses ($\lesssim 10^{9}\;\rm M_{\odot}$). These HAEs show a wide ra… ▽ More By using the ultra-deep \textit{JWST}/MIRI image at 5.6 $μm$ in the Hubble eXtreme Deep Field, we constrain the role of strong H$α$-emitters (HAEs) during Cosmic Reionization at $z\simeq7-8$. Our sample of HAEs is comprised of young ($<35\;\rm Myr$) galaxies, except for one single galaxy ($\approx 300\;\rm Myr$), with low stellar masses ($\lesssim 10^{9}\;\rm M_{\odot}$). These HAEs show a wide range of UV-$β$ slopes, with a median value of $β= -2.15\pm0.21$ which broadly correlates with stellar mass. We estimate the ionizing photon production efficiency ($ξ_{ion,0}$) of these sources (assuming $f_{esc,LyC} = 0\%$), which yields a median value $\rm log_{10}(ξ_{ion,0}/(Hz\;erg^{-1})) = 25.50^{+0.10}_{-0.12}$. We show that $ξ_{ion,0}$ positively correlates with EW$_{0}$(H$α$) and specific star formation rate (sSFR). Instead $ξ_{ion,0}$ weakly anti-correlates with stellar mass and $β$. Based on the $β$ values, we predict $f_{esc, LyC}=4\%^{+3}_{-2}$, which results in $\rm log_{10}(ξ_{ion}/(Hz\;erg^{-1})) = 25.55^{+0.11}_{-0.13}$. Considering this and related findings from the literature, we find a mild evolution of $ξ_{ion}$with redshift. Additionally, our results suggest that these HAEs require only modest escape fractions ($f_{esc, rel}$) of 6$-$15\% to reionize their surrounding intergalactic medium. By only considering the contribution of these HAEs, we estimated their total ionizing emissivity ($\dot{N}_{ion}$) as $\dot{N}_{ion} = 10^{50.53 \pm 0.45}; \text{s}^{-1}\text{Mpc}^{-3}$. When comparing their $\dot{N}_{ion}$ with "non-H$α$ emitter" galaxies across the same redshift range, we find that that strong, young, and low-mass emitters may have played an important role during Cosmic Reionization. △ Less

Submitted 13 June, 2024; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 21 pages, 11 Figures, 1 table. Accepted for publication in ApJ. Updated version 13/06/2024 -- found a mistake in the reported eq. 2 (now corrected)

MIDIS: JWST NIRCam and MIRI unveil the stellar population properties of Ly$α$-emitters and Lyman-Break galaxies at z ~ 3-7

Authors: Edoardo Iani, Karina I. Caputi, Pierluigi Rinaldi, Marianna Annunziatella, Leindert A. Boogaard, Göran Östlin, Luca Costantin, Steven Gillman, Pablo G. Pérez-González, Luis Colina, Gillian Wright, Almudena Alonso-Herrero, Javier Álvarez-Márquez, Arjan Bik, Sarah E. I. Bosman, Alejandro Crespo-Gómez, Andreas Eckart, Thomas R. Greve, Thomas K. Henning, Jens Hjorth, Iris Jermann, Alvaro Labiano, Danial Langeroodi, Jens Melinder, Thibaud Moutard , et al. (7 additional authors not shown)

Abstract: We study the stellar population properties of 182 spectroscopically-confirmed (MUSE/VLT) Lyman-$α$ emitters (LAEs) and 450 photometrically-selected Lyman-Break galaxies (LBGs) at z = 2.8 - 6.7 in the Hubble eXtreme Deep Field (XDF). Leveraging the combined power of HST and JWST NIRCam and MIRI observations, we analyse their rest-frame UV-through-near-IR spectral energy distributions (SEDs) with MI… ▽ More We study the stellar population properties of 182 spectroscopically-confirmed (MUSE/VLT) Lyman-$α$ emitters (LAEs) and 450 photometrically-selected Lyman-Break galaxies (LBGs) at z = 2.8 - 6.7 in the Hubble eXtreme Deep Field (XDF). Leveraging the combined power of HST and JWST NIRCam and MIRI observations, we analyse their rest-frame UV-through-near-IR spectral energy distributions (SEDs) with MIRI playing a crucial role in robustly assessing the LAE's stellar mass and ages. Our LAEs are low-mass objects (log$_{10}$(M$_\star$[M$_\odot$]) ~ 7.5), with little or no dust extinction (E(B - V) ~ 0.1) and a blue UV continuum slope ($β$ ~ -2.2). While 75% of our LAEs are young (< 100 Myr), the remaining 25% have significantly older stellar populations (> 100 Myr). These old LAEs are statistically more massive, less extinct and have lower specific star formation rate (sSFR) compared to young LAEs. Besides, they populate the M$_\star$ - SFR plane along the main-sequence (MS) of star-forming galaxies, while young LAEs populate the starburst region. The comparison between the LAEs properties to those of a stellar-mass matched sample of LBGs shows no statistical difference between these objects, except for the LBGs redder UV continuum slope and marginally larger E(B - V) values. Interestingly, 48% of the LBGs have ages < 10 Myr and are classified as starbursts, but lack detectable Ly$α$ emission. This is likely due to HI resonant scattering and/or selective dust extinction. Overall, we find that JWST observations are crucial in determining the properties of LAEs and shedding light on the properties and similarities between LAEs and LBGs. △ Less

Submitted 25 January, 2024; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: 22 pages, 15 figures, 5 tables. Accepted for publication in APJ

Spatially-resolved H$α$ and ionizing photon production efficiency in the lensed galaxy MACS1149-JD1 at a redshift of 9.11

Authors: J. Álvarez-Márquez, L. Colina, A. Crespo Gómez, P. Rinaldi, J. Melinder, G. Östlin, M. Annunziatella, A. Labiano, A. Bik, S. Bosman, T. R. Greve, G. Wright, A. Alonso-Herrero, L. Boogaard, R. Azollini, K. I. Caputi, L. Costantin, A. Eckart, M. GarcÍa-MarÍn, S. Gillman, J. Hjorth, E. Iani, O. Ilbert, I. Jermann, D. Langeroodi , et al. (10 additional authors not shown)

Abstract: We present MIRI/JWST medium-resolution spectroscopy (MRS) and imaging (MIRIM) of the lensed galaxy MACS1149-JD1 at a redshift of $z$=9.1092$\pm$0.0002 (Universe age about 530 Myr). We detect, for the first time, spatially resolved H$α$ emission in a galaxy at a redshift above nine. The structure of the H$α$ emitting gas consists of two clumps, S and N. The total H$α$ luminosity implies an instanta… ▽ More We present MIRI/JWST medium-resolution spectroscopy (MRS) and imaging (MIRIM) of the lensed galaxy MACS1149-JD1 at a redshift of $z$=9.1092$\pm$0.0002 (Universe age about 530 Myr). We detect, for the first time, spatially resolved H$α$ emission in a galaxy at a redshift above nine. The structure of the H$α$ emitting gas consists of two clumps, S and N. The total H$α$ luminosity implies an instantaneous star-formation of 5.3$\pm$0.4 $M_{\odot}$ yr$^{-1}$ for solar metallicities. The ionizing photon production efficiency, $\log(ζ_\mathrm{ion})$, shows a spatially resolved structure with values of 25.55$\pm$0.03, 25.47$\pm$0.03, and 25.91$\pm$0.09 Hz erg$^{-1}$ for the integrated galaxy, and clumps S and N, respectively. The H$α$ rest-frame equivalent width, EW$_{0}$(H$α$), is 726$^{+660}_{-182}$ Ángstrom for the integrated galaxy, but presents extreme values of 531$^{+300}_{-96}$ Ángstrom and $\geq$1951 Ángstrom for clumps S and N, respectively. The spatially resolved ionizing photon production efficiency is within the range of values measured in galaxies at redshift above six, and well above the canonical value (25.2$\pm$0.1 Hz erg$^{-1}$). The extreme difference of EW$_{0}$(H$α$) for Clumps S and N indicates the presence of a recent (<5 Myrs) burst in clump N and a star formation over a larger period of time (e.g., $\sim$50 Myr) in clump S. Finally, clump S and N show very different H$α$ kinematics with velocity dispersions of 56$\pm$4 km s$^{-1}$ and 113$\pm$33 km s$^{-1}$, likely indicating the presence of outflows or increased turbulence in the clump N. The dynamical mass, $M_\mathrm{dyn}$= (2.4$\pm$0.5)$\times$10$^{9}$ $M_{\odot}$, is within the range measured with the spatially resolved [OIII]88$μ$m line. △ Less

Submitted 24 March, 2024; v1 submitted 12 September, 2023; originally announced September 2023.

Comments: 11 pages, 6 figures, published version in A&A

Journal ref: A&A 686, A85 (2024)

MIDIS: JWST/MIRI reveals the Stellar Structure of ALMA-selected Galaxies in the Hubble-UDF at Cosmic Noon

Authors: Leindert A. Boogaard, Steven Gillman, Jens Melinder, Fabian Walter, Luis Colina, Göran Östlin, Karina I. Caputi, Edoardo Iani, Pablo Pérez-González, Paul van der Werf, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Javier Álvarez-Márquez, Marianna Annunziatella, Arjan Bik, Sarah Bosman, Luca Costantin, Alejandro Crespo Gómez, Dan Dicken, Andreas Eckart, Jens Hjorth, Iris Jermann, Alvaro Labiano, Danial Langeroodi , et al. (8 additional authors not shown)

Abstract: We present deep James Webb Space Telescope (JWST)/MIRI F560W observations of a flux-limited, ALMA-selected sample of 28 galaxies at z=0.5-3.6 in the Hubble Ultra Deep Field (HUDF). The data from the MIRI Deep Imaging Survey (MIDIS) reveal the stellar structure of the HUDF galaxies at rest-wavelengths of >1 micron for the first time. We revise the stellar mass estimates using new JWST photometry an… ▽ More We present deep James Webb Space Telescope (JWST)/MIRI F560W observations of a flux-limited, ALMA-selected sample of 28 galaxies at z=0.5-3.6 in the Hubble Ultra Deep Field (HUDF). The data from the MIRI Deep Imaging Survey (MIDIS) reveal the stellar structure of the HUDF galaxies at rest-wavelengths of >1 micron for the first time. We revise the stellar mass estimates using new JWST photometry and find good agreement with pre-JWST analysis; the few discrepancies can be explained by blending issues in the earlier lower-resolution Spitzer data. At z~2.5, the resolved rest-frame near-infrared (1.6 micron) structure of the galaxies is significantly more smooth and centrally concentrated than seen by HST at rest-frame 450 nm (F160W), with effective radii of Re(F560W)=1-5 kpc and Sérsic indices mostly close to an exponential (disk-like) profile (n~1), up to n~5 (excluding AGN). We find an average size ratio of Re(F560W)/Re(F160W)~0.7 that decreases with stellar mass. The stellar structure of the ALMA-selected galaxies is indistinguishable from a HUDF reference sample of galaxies with comparable MIRI flux density. We supplement our analysis with custom-made, position-dependent, empirical PSF models for the F560W observations. The results imply that a smoother stellar structure is in place in massive gas-rich, star-forming galaxies at Cosmic Noon, despite a more clumpy rest-frame optical appearance, placing additional constraints on galaxy formation simulations. As a next step, matched-resolution, resolved ALMA observations will be crucial to further link the mass- and light-weighted galaxy structures to the dusty interstellar medium. △ Less

Submitted 26 April, 2024; v1 submitted 31 August, 2023; originally announced August 2023.

Comments: 19 pages, 12 figures, 1 table, Accepted for publication in ApJ

An extremely metal poor star complex in the reionization era: Approaching Population III stars with JWST

Authors: E. Vanzella, F. Loiacono, P. Bergamini, U. Mestric, M. Castellano, P. Rosati, M. Meneghetti, C. Grillo, F. Calura, M. Mignoli, M. Bradac, A. Adamo, G. Rihtarsic, M. Dickinson, M. Gronke, A. Zanella, F. Annibali, C. Willott, M. Messa, E. Sani, A. Acebron, A. Bolamperti, A. Comastri, R. Gilli, K. I. Caputi , et al. (9 additional authors not shown)

Abstract: We present JWST/NIRSpec integral field spectroscopy (IFS) of a lensed Population III candidate stellar complex (dubbed Lensed And Pristine 1, LAP1), with a lensing-corrected stellar mass ~<10^4 Msun, absolute luminosity M_UV > -11.2 (m_UV > 35.6), confirmed at redshift 6.639 +/- 0.004. The system is strongly amplified (μ>~ 100) by straddling a critical line of the Hubble Frontier Field galaxy clus… ▽ More We present JWST/NIRSpec integral field spectroscopy (IFS) of a lensed Population III candidate stellar complex (dubbed Lensed And Pristine 1, LAP1), with a lensing-corrected stellar mass ~<10^4 Msun, absolute luminosity M_UV > -11.2 (m_UV > 35.6), confirmed at redshift 6.639 +/- 0.004. The system is strongly amplified (μ>~ 100) by straddling a critical line of the Hubble Frontier Field galaxy cluster MACS J0416. Despite the stellar continuum is currently not detected in the Hubble and JWST/NIRCam and NIRISS imaging, arclet-like shapes of Lyman and Balmer lines, Lya, Hg, Hb and Ha are detected with NIRSpec IFS with signal-to-noise ratios SNR=5-13 and large equivalent widths (>300-2000A), along with a remarkably weak [OIII]4959-5007 at SNR ~ 4. LAP1 shows a large ionizing photon production efficiency, log(ξ_{ion}[erg~Hz^{-1}])>26. From the metallicity indexes R23 = ([OIII]4959-5007 + [OII]3727) / Hb ~< 0.74 and R3 = ([OIII]5007 / Hb) = 0.55 +/- 0.14, we derive an oxygen abundance 12+log(O/H) ~< 6.3. Intriguingly, the Ha emission is also measured in mirrored sub-components where no [OIII] is detected, providing even more stringent upper limits on the metallicity if in-situ star formation is ongoing in this region (12+log(O/H) < 6, or Z < 0.002 Zsun). The formal stellar mass limit of the sub-components would correspond to ~10^{3} Msun or M_UV fainter than -10. Alternatively, such a metal-free pure line emitting region could be the first case of a fluorescing HI gas region, induced by transverse escaping ionizing radiation from a nearby star-complex. The presence of large equivalent-width hydrogen lines and the deficiency of metal lines in such a small region, make LAP1 the most metal poor star-forming region currently known in the reionization era and a promising site that may host isolated, pristine stars. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 15 pages, 10 figures, 1 table. Submitted to A&A. Comments welcome

Journal ref: A&A 678, A173 (2023)

The Spitzer Coverage of HSC-Deep with IRAC for Z studies (SHIRAZ) I: IRAC mosaics

Authors: Marianna Annunziatella, Anna Sajina, Mauro Stefanon, Danilo Marchesini, Mark Lacy, Ivo Labbe, Lilianna Houston, Rachel Bezanson, Eiichi Egami, Xiaohui Fan, Duncan Farrah, Jenny Greene, Andy Goulding, Yen-Ting Lin, Xin Liu, Thibaud Moutard, Yoshiaki Ono, Masami Ouchi, Marcin Sawicki, Jason Surace, Katherine Whitaker

Abstract: We present new Spitzer Infrared Array Camera (IRAC) 3.6 and 4.5μm mosaics of three fields, E-COSMOS, DEEP2-F3, and ELAIS-N1. Our mosaics include both new IRAC observations as well as re-processed archival data in these fields. These fields are part of the HSC-Deep grizy survey and have a wealth of additional ancillary data. The addition of these new IRAC mosaics is critical in allowing for improve… ▽ More We present new Spitzer Infrared Array Camera (IRAC) 3.6 and 4.5μm mosaics of three fields, E-COSMOS, DEEP2-F3, and ELAIS-N1. Our mosaics include both new IRAC observations as well as re-processed archival data in these fields. These fields are part of the HSC-Deep grizy survey and have a wealth of additional ancillary data. The addition of these new IRAC mosaics is critical in allowing for improved photometric redshifts and stellar population parameters at cosmic noon and earlier epochs. The total area mapped by this work is {\sim} 17 deg2 with a mean integration time of {\sim}1200s, providing a median 5σ depth of 23.7(23.3) at 3.6(4.5)μm in AB. We perform SExtractor photometry both on the combined mosaics as well as the single-epoch mosaics taken {\sim}6 months apart. The resultant IRAC number counts show good agreement with previous studies. In combination with the wealth of existing and upcoming spectro-photometric data in these fields, our IRAC mosaics will enable a wide range of galactic evolution and AGN studies. With that goal in mind, we make the combined IRAC mosaics and coverage maps of these three fields publicly available. counts show good agreement with previous studies. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: accepted for publication in AJ

Probing the Star Formation Main Sequence down to $10^{8}$ M$_\odot$ at $1.0<z<3.0$

Authors: Rosa M. Mérida, Pablo G. Pérez-González, Patricia Sánchez-Blázquez, Ángela García-Argumánez, Marianna Annunziatella, Luca Costantin, Alejandro Lumbreras-Calle, Belén Alcalde-Pampliega, Guillermo Barro, Néstor Espino-Briones, Anton M. Koekemoer

Abstract: We investigate the star formation main sequence (MS) (SFR-M$_{\star}$) down to 10$^{8-9}\mathrm{M}_\odot$ using a sample of 34,061 newly-discovered ultra-faint ($27\lesssim i \lesssim 30$ mag) galaxies at $1<z<3$ detected in the GOODS-N field. Virtually these galaxies are not contained in previous public catalogs, effectively doubling the number of known sources in the field. The sample was constr… ▽ More We investigate the star formation main sequence (MS) (SFR-M$_{\star}$) down to 10$^{8-9}\mathrm{M}_\odot$ using a sample of 34,061 newly-discovered ultra-faint ($27\lesssim i \lesssim 30$ mag) galaxies at $1<z<3$ detected in the GOODS-N field. Virtually these galaxies are not contained in previous public catalogs, effectively doubling the number of known sources in the field. The sample was constructed by stacking the optical broad-band observations taken by the HST/GOODS-CANDELS surveys as well as the 25 ultra-deep medium-band images gathered by the GTC/SHARDS project. Our sources are faint (average observed magnitudes $<i>\sim28.2$ mag, $<H>\sim27.9$ mag), blue (UV-slope $<β>\sim-1.9$), star-forming (rest-frame colors $<U-V>\sim0.10$ mag, $<V-J>\sim0.17$ mag) galaxies. These observational characteristics are identified with young (mass-weighted age $<\mathrm{t_{M-w}}>\sim0.014$ Gyr) stellar populations subject to low attenuations ($<\mathrm{A(V)}>\sim0.30$ mag). Our sample allows us to probe the MS down to $10^{8.0}\,\mathrm{M}_\odot$ at $z=1$ and $10^{8.5}\,\mathrm{M}_\odot$ at $z=3$, around 0.6 dex deeper than previous analysis. In the low-mass galaxy regime, we find an average value for the slope of 0.97 at $1<z<2$ and 1.12 at $2<z<3$. Nearly $\sim$60% of our sample presents stellar masses in the range $10^{6-8}$ M$_\odot$ between $1<z<3$. If the slope of the MS remained constant in this regime, the sources populating the low-mass tail of our sample would qualify as starburst galaxies. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: 34 pages, 16 figures, 9 tables. Accepted for publication to ApJ

Life beyond 30: probing the -20<M_UV<-17 luminosity function at 8<z<13 with the NIRCam parallel field of the MIRI Deep Survey

Authors: Pablo G. Pérez-González, Luca Costantin, Danial Langeroodi, Pierluigi Rinaldi, Marianna Annunziatella, Olivier Ilbert, Luis Colina, Hans Ulrik Noorgaard-Nielsen, Thomas Greve, Göran Ostlin, Gillian Wright, Almudena Alonso-Herrero, Javier Álvarez-Márquez, Karina I. Caputi, Andreas Eckart, Olivier Le Fèvre, Álvaro Labiano, Macarena García-Marín, Jens Hjorth, Sarah Kendrew, John P. Pye, Tuomo Tikkanen, Paul van der Werf, Fabian Walter, Martin Ward , et al. (19 additional authors not shown)

Abstract: We present the ultraviolet luminosity function and an estimate of the cosmic star formation rate density at $8<z<13$ derived from deep NIRCam observations taken in parallel with the MIRI Deep Survey (MDS) of the Hubble Ultra Deep Field (HUDF), NIRCam covering the parallel field 2 (HUDF-P2). Our deep (40 hours) NIRCam observations reach a F277W magnitude of 30.8 ($5σ$), more than 2 magnitudes deepe… ▽ More We present the ultraviolet luminosity function and an estimate of the cosmic star formation rate density at $8<z<13$ derived from deep NIRCam observations taken in parallel with the MIRI Deep Survey (MDS) of the Hubble Ultra Deep Field (HUDF), NIRCam covering the parallel field 2 (HUDF-P2). Our deep (40 hours) NIRCam observations reach a F277W magnitude of 30.8 ($5σ$), more than 2 magnitudes deeper than JWST public datasets already analyzed to find high redshift galaxies. We select a sample of 44 $z>8$ galaxy candidates based on their dropout nature in the F115W and/or F150W filters, a high probability for their photometric redshifts, estimated with three different codes, being at $z>8$, good fits based on $χ^2$ calculations, and predominant solutions compared to $z<8$ alternatives. We find mild evolution in the luminosity function from $z\sim13$ to $z\sim8$, i.e., only a small increase in the average number density of $\sim$0.2 dex, while the faint-end slope and absolute magnitude of the knee remain approximately constant, with values $α=-2.2\pm0.1$ and $M^*=-20.8\pm0.2$ mag. Comparing our results with the predictions of state-of-the-art galaxy evolution models, we find two main results: (1) a slower increase with time in the cosmic star formation rate density compared to a steeper rise predicted by models; (2) nearly a factor of 10 higher star formation activity concentrated in scales around 2 kpc in galaxies with stellar masses $\sim10^8$ M$_\odot$ during the first 350 Myr of the Universe, $z\sim12$, with models matching better the luminosity density observational estimations $\sim$150 Myr later, by $z\sim9$. △ Less

Submitted 22 May, 2023; v1 submitted 5 February, 2023; originally announced February 2023.

Comments: Accepted for publication in ApJL, in memoriam Hans Ulrik Noorgaard-Nielsen and Olivier Le Fèvre

MIDIS: Strong (Hb + [OIII]) and Ha emitters at redshift $z \simeq 7-8$ unveiled with JWST/NIRCam and MIRI imaging in the Hubble eXtreme Deep Field (XDF)

Authors: P. Rinaldi, K. I. Caputi, L. Costantin, S. Gillman, E. Iani, P. G. Perez Gonzalez, G. Oestlin, L. Colina, T. Greve, H. U. Noorgard-Nielsen, G. S. Wright, A. Alonso-Herrero, J. Alvarez-Marquez, A. Eckart, M. Garcia-Marin, J. Hjorth, O. Ilbert, S. Kendrew, A. Labiano, O. Le Fevre, J. Pye, T. Tikkanen, F. Walter, P. van der Werf, M. Ward , et al. (22 additional authors not shown)

Abstract: We make use of \textit{JWST} medium and broad-band NIRCam imaging, along with ultra-deep MIRI $5.6 \rm μm$ imaging, in the Hubble eXtreme Deep Field (XDF) to identify prominent line emitters at $z\simeq 7-8$. Out of a total of 58 galaxies at $z\simeq 7-8$, we find 18 robust candidates ($\simeq$31\%) for (H$β$ + [OIII]) emitters, based on their enhanced fluxes in the F430M and F444W filters, with E… ▽ More We make use of \textit{JWST} medium and broad-band NIRCam imaging, along with ultra-deep MIRI $5.6 \rm μm$ imaging, in the Hubble eXtreme Deep Field (XDF) to identify prominent line emitters at $z\simeq 7-8$. Out of a total of 58 galaxies at $z\simeq 7-8$, we find 18 robust candidates ($\simeq$31\%) for (H$β$ + [OIII]) emitters, based on their enhanced fluxes in the F430M and F444W filters, with EW$_{0}$(H$β$ +[OIII]) $\simeq 87 - 2100$ Å. Among these emitters, 16 lie in the MIRI coverage area and 12 exhibit a clear flux excess at $5.6 \, \rm μm$, indicating the simultaneous presence of a prominent H$α$ emission line with EW$_{0}$(H$α$) $\simeq 200-3000$ Å. This is the first time that H$α$ emission can be detected in individual galaxies at $z>7$. The H$α$ line, when present, allows us to separate the contributions of H$β$ and [OIII] to the (H$β$ +[OIII]) complex, and derive H$α$-based star formation rates (SFRs). We find that in most cases [OIII]/H$β> 1$. Instead, two galaxies have [OIII]/H$β< 1$, indicating that the NIRCam flux excess is mainly driven by H$β$. This could potentially imply extremely low metallicities. Most prominent line emitters are very young starbursts or galaxies on their way to/from the starburst cloud. They make for a cosmic SFR density $\rm log_{10}(ρ_{SFR_{Hα}}) \simeq -2.35$, which is about a quarter of the total value ($\rm log_{10}(ρ_{SFR_{tot}}) \simeq -1.76$) at $z\simeq 7-8$. Therefore, the strong H$α$ emitters likely had a significant role in reionization. △ Less

Submitted 5 June, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: 18 pages, 11 figures. The paper has been accepted for publication in ApJ

VST-GAME: Galaxy Assembly as a function of Mass and Environment with VST. Photometric assessment and density field of MACSJ0416

Authors: Nicolas Estrada, Amata Mercurio, Benedetta Vulcani, Giulia Rodighiero, Mario Nonino, Marianna Annunziatella, Piero Rosati, Claudio Grillo, Gabriel Bartosch Caminha, Giuseppe Angora, Andrea Biviano, Massimo Brescia, Gabriella De Lucia, Ricardo Demarco, Marisa Girardi, Raphael Gobat, Brian C. Lemaux

Abstract: Observational studies have widely demonstrated that galaxy physical properties are strongly affected by the surrounding environment. On one side, gas inflows provide galaxies with new fuel for star formation. On the other side, the high temperatures and densities of the medium are expected to induce quenching in the star formation. Observations of large structures, in particular filaments at the c… ▽ More Observational studies have widely demonstrated that galaxy physical properties are strongly affected by the surrounding environment. On one side, gas inflows provide galaxies with new fuel for star formation. On the other side, the high temperatures and densities of the medium are expected to induce quenching in the star formation. Observations of large structures, in particular filaments at the cluster outskirts (r>2r$_{200}$), are currently limited to the low redshift Universe. We present a multi-band dataset for the cluster MACS J0416.1-2403 (z=0.397), observed in the context of the Galaxy Assembly as a function of Mass and Environment with VST (VST-GAME) survey. The project aims at gathering deep ($r$<24.4) and wide (20x20Mpc$^2$) observations at optical wavelengths for six massive galaxy clusters at 0.2<z<0.6, complemented with near infrared data. This work describes the photometric analysis of the cluster, defines a density field and studies galaxy properties in the cluster outskirts. We extract sources paying particular attention to recover the faintest ones. We combine all the extractions in a multi-band catalog and compute photometric redshifts. We then define cluster memberships up to 5r$_{200}$ from the cluster core and measure the density field, comparing galaxy properties in different environments. We found that the $g-r$ colors show bimodal behaviours in all the environments, but the peak of the distribution of red galaxies shifts toward redder colors with increasing density and the fraction of galaxies in the blue cloud increases with decreasing density. We also found 3 overdense regions in the cluster outskirts at r$\sim$5r$_{200}$. The color of galaxies suggests the presence of evolved galaxy populations, an insight for pre-processing phenomena over these substructures. We release the multi-band catalog, down to the completeness limit $r$<24.4 mag. △ Less

Submitted 23 March, 2023; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: 19 pages, 15 figures and 5 tables

Journal ref: A&A 671, A146 (2023)

CEERS Key Paper V: A triality on the nature of HST-dark galaxies

Authors: Pablo G. Pérez-González, Guillermo Barro, Marianna Annunziatella, Luca Costantin, Ángela García-Argumánez, Elizabeth J. McGrath, Rosa M. Mérida, Jorge A. Zavala, Pablo Arrabal Haro, Micaela B. Bagley, Bren E. Backhaus, Peter Behroozi, Eric F. Bell, Laura Bisigello, Véronique Buat, Antonello Calabrò, Caitlin M. Casey, Nikko J. Cleri, Rosemary T. Coogan, M. C. Cooper, Asantha R. Cooray, Avishai Dekel, Mark Dickinson, David Elbaz, Henry C. Ferguson , et al. (30 additional authors not shown)

Abstract: The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometr… ▽ More The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometric redshifts in 2 dimensions and stellar populations properties in a pixel-by-pixel basis. We select 138 galaxies with F150W-F356W>1.5 mag, F356W<27.5 mag. The nature of these sources is threefold: (1) 71% are dusty star-forming galaxies at 2<z<6 with masses 9<log M/M_sun<11 and a variety of specific SFRs (<1 to >100 Gyr^-1); (2) 18% are quiescent/dormant (i.e., subject to reignition and rejuvenation) galaxies at 3<z<5, masses log M/M_sun~10 and post-starburst stellar mass-weighted ages (0.5-1 Gyr); and (3) 11% are strong young starbursts with indications of high-EW emission lines (typically, [OIII]+Hbeta) at 6<z<7 and log M/M_sun~9.5. The sample is dominated by disk-like galaxies with a remarkable compactness for XELG-z6 (effective radii smaller than 0.4 kpc). Large attenuations in SFGs, 2<A(V)<5 mag, are found within 1.5 times the effective radius, approximately 2 kpc, while QGs present A(V)~0.2 mag. Our SED-fitting technique reproduces the expected dust emission luminosities of IR-bright and sub-millimeter galaxies. This study implies high levels of star formation activity between z~20 and z~10, where virtually 100% of our galaxies had already formed 10^8 M_sun of their stellar content, 60% of them had assembled 10^9 M_sun, and 10% up to 10^10 M_sun (in situ or ex situ). (abridged) △ Less

Submitted 3 April, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: Published in CEERS ApJL Focus Issue, ApJL 946, L16

Journal ref: ApJL 946, L16 https://iopscience.iop.org/article/10.3847/2041-8213/acb3a5

Lack of influence of the environment in the earliest stages of massive galaxy formation

Authors: Marianna Annunziatella, PabloG. Pérez-González, Ángela García Argumánez, Guillermo Barro, Belén Alcalde Pampliega, Luca Costantin, Anton M. Koekemoer, Rosa M. Mérida

Abstract: We investigate how the environment affects the assembly history of massive galaxies. For that purpose, we make use of SHARDS and HST spectro-photometric data, whose depth, spectral resolution, and wavelength coverage allow to perform a detailed analysis of the stellar emission as well as obtaining unprecedentedly accurate photometric redshifts. This expedites a sufficiently accurate estimate of th… ▽ More We investigate how the environment affects the assembly history of massive galaxies. For that purpose, we make use of SHARDS and HST spectro-photometric data, whose depth, spectral resolution, and wavelength coverage allow to perform a detailed analysis of the stellar emission as well as obtaining unprecedentedly accurate photometric redshifts. This expedites a sufficiently accurate estimate of the local environment and a robust derivation of the star formation histories of a complete sample of 332 massive galaxies ($\mathrm{>10^{10}M_{\odot}}$) at redshift $1\leq z \leq 1.5$ in the GOODS-N field. We find that massive galaxies in this redshift range avoid the lowest density environments. Moreover, we observed that the oldest galaxies in our sample with with mass-weighted formation redshift $\mathrm{\overline{z}_{M-w} \geq 2.5}$, avoid the highest density regions, preferring intermediate environments. Younger galaxies, including those with active star formation, tend to live in denser environments ($Σ= \mathrm{5.0_{1.1}^{24.8}\times 10^{10}M_{\odot}Mpc^{-2}}$). This behavior could be expected if those massive galaxies starting their formation first would merge with neighbors and sweep their environment earlier. On the other hand, galaxies formed more recently ($\overline{z}_{M-w} < 2.5$) are accreted into large scale structures at later times and we are observing them before sweeping their environment or, alternatively, they are less likely to affect their environment. However, given that both number and mass surface densities of neighbor galaxies is relatively low for the oldest galaxies, our results reveal a very weak correlation between environment and the first formation stages of the earliest massive galaxies. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: Accepted for publication in MNRAS

MAGAZ3NE: High Stellar Velocity Dispersions for Ultra-Massive Quiescent Galaxies at $z\gtrsim3$

Authors: Ben Forrest, Gillian Wilson, Adam Muzzin, Danilo MArchesini, M. C. Cooper, Z. Cemile Marsan, Marianna Annunziatella, Ian McConachie, Kumail Zaidi, Percy Gomez, Stephanie M. Urbano Stawinski, Wenjun Chang, Gabriella de Lucia, Francesco La Barbera, Lori Lubin, Julie Nantais, Theodore Peña, Paolo Saracco, Jason Surace, Mauro Stefanon

Abstract: In this work we publish stellar velocity dispersions, sizes, and dynamical masses for 8 ultra-massive galaxies (UMGs; log($M$/M$_\odot>11$, $z\gtrsim3$) from the Massive Ancient Galaxies At $z>3$ NEar-infrared (MAGAZ3NE) Survey, more than doubling the number of such galaxies with velocity dispersion measurements at this epoch. Using the deep Keck/MOSFIRE and Keck/NIRES spectroscopy of these object… ▽ More In this work we publish stellar velocity dispersions, sizes, and dynamical masses for 8 ultra-massive galaxies (UMGs; log($M$/M$_\odot>11$, $z\gtrsim3$) from the Massive Ancient Galaxies At $z>3$ NEar-infrared (MAGAZ3NE) Survey, more than doubling the number of such galaxies with velocity dispersion measurements at this epoch. Using the deep Keck/MOSFIRE and Keck/NIRES spectroscopy of these objects in the $H$- and $K$-bandpasses, we obtain large velocity dispersions of $\sim400$ km s$^{-1}$ for most of the objects, which are some of the highest stellar velocity dispersions measured, and $\sim40$\% larger than those measured for galaxies of similar mass at $z\sim1.7$. The sizes of these objects are also smaller by a factor of 1.5-3 compared to this same $z\sim1.7$ sample. We combine these large velocity dispersions and small sizes to obtain dynamical masses. The dynamical masses are similar to the stellar masses of these galaxies, consistent with a Chabrier initial mass function (IMF). Considered alongside previous studies of massive quiescent galaxies across $0.2<z<4.0$, there is evidence for an evolution in the relation between the dynamical mass - stellar mass ratio and velocity dispersion as a function of redshift. This implies an IMF with fewer low mass stars (e.g., Chabrier IMF) for massive quiescent galaxies at higher redshifts in conflict with the bottom-heavy IMF (e.g., Salpeter IMF) found in their likely $z\sim0$ descendants, though a number of alternative explanations such as a different dynamical structure or significant rotation are not ruled out. Similar to data at lower redshifts, we see evidence for an increase of IMF normalization with velocity dispersion, though the $z\gtrsim3$ trend is steeper than that for $z\sim0.2$ early-type galaxies and offset to lower dynamical-to-stellar mass ratios. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: 17 pages and 8 figures in the main text, plus references and Appendices for a total of 25 pages and 14 figures. Accepted for publication in The Astrophysical Journal

Dusty Starbursts Masquerading as Ultra-high Redshift Galaxies in JWST CEERS Observations

Authors: Jorge A. Zavala, Veronique Buat, Caitlin M. Casey, Denis Burgarella, Steven L. Finkelstein, Micaela B. Bagley, Laure Ciesla, Emanuele Daddi, Mark Dickinson, Henry C. Ferguson, Maximilien Franco, E. F. Jim'enez-Andrade, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Aurélien Le Bail, E. J. Murphy, Casey Papovich, Sandro Tacchella, Stephen M. Wilkins, Itziar Aretxaga, Peter Behroozi, Jaclyn B. Champagne, Adriano Fontana, Mauro Giavalisco, Andrea Grazian , et al. (99 additional authors not shown)

Abstract: Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may als… ▽ More Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations. △ Less

Submitted 30 January, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: Published in The Astrophysical Journal Letters (updated to match the published version)

Journal ref: ApJL 943 L9 (2023)

Probing the earliest phases in the formation of massive galaxies with simulated HST+JWST imaging data from Illustris

Authors: Ángela García-Argumánez, Pablo G. Pérez-González, Armando Gil de Paz, Gregory F. Snyder, Pablo Arrabal Haro, Micaela B. Bagley, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Anton Koekemoer, Casey Papovich, Nor Pirzkal, Harry C. Ferguson, L. Y. Aaron Yung, Marianna Annunziatella, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Benne W. Holwerda, Rosa María Mérida González, Caitlin Rose, Mauro Giavalisco, Norman A. Grogin, Dale D. Kocevski

Abstract: We use the Illustris-1 simulation to explore the capabilities of the $\textit{Hubble}$ and $\textit{James Webb Space Telescope}$ data to analyze the stellar populations in high-redshift galaxies, taking advantage of the combined depth, spatial resolution, and wavelength coverage. For that purpose, we use simulated broad-band ACS, WFC3 and NIRCam data and 2-dimensional stellar population synthesis… ▽ More We use the Illustris-1 simulation to explore the capabilities of the $\textit{Hubble}$ and $\textit{James Webb Space Telescope}$ data to analyze the stellar populations in high-redshift galaxies, taking advantage of the combined depth, spatial resolution, and wavelength coverage. For that purpose, we use simulated broad-band ACS, WFC3 and NIRCam data and 2-dimensional stellar population synthesis (2D-SPS) to derive the integrated star formation history (SFH) of massive (M$_{\ast}>10^{10}\,$M$_{\odot}$) simulated galaxies at $1<z<4$ that evolve into a local M$_{\ast}>10^{11}\,$M$_{\odot}$ galaxy. In particular, we explore the potential of HST and JWST datasets reaching a depth similar to those of the CANDELS and ongoing CEERS observations, respectively, and concentrate on determining the capabilities of this dataset for characterizing the first episodes in the SFH of local M$_{\ast}>10^{11}\,$M$_{\odot}$ galaxies by studying their progenitors at $z>1$. The 2D-SPS method presented in this paper has been calibrated to robustly recover the cosmic times when the first star formation episodes occurred in massive galaxies, i.e., the first stages in their integrated SFHs. In particular, we discuss the times when the first 1% to 50% of their total stellar mass formed in the simulation. We demonstrate that we can recover these ages with typical median systematic offset of less than 5% and scatter around 20%-30%. According to our measurements on Illustris data, we are able to recover that local M$_{\ast}>10^{11}\,$M$_{\odot}$ galaxies would have started their formation by $z=16$, forming the first 5% of their stellar mass present at $z \sim 1$ by $z=4.5$, 10% by $z=3.7$, and 25% by $z=2.7$. △ Less

Submitted 19 December, 2022; v1 submitted 28 July, 2022; originally announced July 2022.

Comments: 28 pages, 13 figures, 4 tables. ApJ in press. Summary of changes from original submission: the major change is that we now include in Sec. 6 the comparison of the results obtained for our sample of massive 1 < z < 4 progenitors with those obtained by considering all massive galaxies at 1 < z < 4 in the simulated images. Several figures and sections have been updated

A Long Time Ago in a Galaxy Far, Far Away: A Candidate z ~ 12 Galaxy in Early JWST CEERS Imaging

Authors: Steven L. Finkelstein, Micaela B. Bagley, Pablo Arrabal Haro, Mark Dickinson, Henry C. Ferguson, Jeyhan S. Kartaltepe, Casey Papovich, Denis Burgarella, Dale D. Kocevski, Marc Huertas-Company, Kartheik G. Iyer, Rebecca L. Larson, Pablo G. Pérez-González, Caitlin Rose, Sandro Tacchella, Stephen M. Wilkins, Katherine Chworowsky, Aubrey Medrano, Alexa M. Morales, Rachel S. Somerville, L. Y. Aaron Yung, Adriano Fontana, Mauro Giavalisco, Andrea Grazian, Norman A. Grogin , et al. (95 additional authors not shown)

Abstract: We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging f… ▽ More We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang. △ Less

Submitted 7 September, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

Redshift identification of X-ray selected active galactic nuclei in the J1030 field: searching for large-scale structures and high-redshift sources

Authors: Stefano Marchesi, Marco Mignoli, Roberto Gilli, Alessandro Peca, Micol Bolzonella, Riccardo Nanni, Marianna Annunziatella, Barbara Balmaverde, Marcella Brusa, Francesco Calura, Letizia P. Cassarà, Marco Chiaberge, Andrea Comastri, Felice Cusano, Quirino D'Amato, Kazushi Iwasawa, Giorgio Lanzuisi, Danilo Marchesini, Takahiro Morishita, Isabella Prandoni, Andrea Rossi, Paolo Tozzi, Cristian Vignali, Fabio Vito, Giovanni Zamorani , et al. (1 additional authors not shown)

Abstract: We publicly release the spectroscopic and photometric redshift catalog of the sources detected with Chandra in the field of the $z$=6.3 quasar SDSS J1030+0525. This is currently the fifth deepest X-ray field, and reaches a 0.5-2 keV flux limit $f_{\rm 0.5-2}$=6$\times$10$^{-17}$ erg s$^{-1}$ cm$^{-2}$. By using two independent methods, we measure a photometric redshift for 243 objects, while 123 (… ▽ More We publicly release the spectroscopic and photometric redshift catalog of the sources detected with Chandra in the field of the $z$=6.3 quasar SDSS J1030+0525. This is currently the fifth deepest X-ray field, and reaches a 0.5-2 keV flux limit $f_{\rm 0.5-2}$=6$\times$10$^{-17}$ erg s$^{-1}$ cm$^{-2}$. By using two independent methods, we measure a photometric redshift for 243 objects, while 123 (51%) sources also have a spectroscopic redshift, 110 of which coming from an INAF-Large Binocular Telescope (LBT) Strategic Program. We use the spectroscopic redshifts to determine the quality of the photometric ones, and find it in agreement with that of other X-ray surveys which used a similar number of photometric data-points. In particular, we measure a sample normalized median absolute deviation $σ_{NMAD}$=1.48||$z_{phot}$-$z_{spec}$||/(1+$z_{spec}$)=0.065. We use these new spectroscopic and photometric redshifts to study the properties of the Chandra J1030 field. We observe several peaks in our spectroscopic redshift distribution between $z$=0.15 and $z$=1.5, and find that the sources in each peak are often distributed across the whole Chandra field of view. This evidence confirms that X-ray selected AGN can efficiently track large-scale structures over physical scales of several Mpc. Finally, we computed the Chandra J1030 $z>$3 number counts: while the spectroscopic completeness at high-redshift of our sample is limited, our results point towards a potential source excess at $z\geq$4, which we plan to either confirm or reject in the near future with dedicated spectroscopic campaigns. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Comments: 21 pages, 15 figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 656, A117 (2021)

Spectroscopic Confirmation of a Protocluster at $z=3.37$ with a High Fraction of Quiescent Galaxies

Authors: Ian McConachie, Gillian Wilson, Ben Forrest, Z. Cemile Marsan, Adam Muzzin, M. C. Cooper, Marianna Annunziatella, Danilo Marchesini, Jeffrey C. C. Chan, Percy Gomez, Mohamed H. Abdullah, Paolo Saracco, Julie Nantais

Abstract: We report the discovery of MAGAZ3NE J095924+022537, a spectroscopically-confirmed protocluster at $z = 3.3665^{+0.0009}_{-0.0012}$ around a spectroscopically-confirmed $UVJ$-quiescent ultra-massive galaxy (UMG; $M_{\star}=2.34^{+0.23}_{-0.34}\times10^{11} {\rm M}_\odot$) in the COSMOS UltraVISTA field. We present a total of 38 protocluster members (14 spectroscopic and 24 photometric), including t… ▽ More We report the discovery of MAGAZ3NE J095924+022537, a spectroscopically-confirmed protocluster at $z = 3.3665^{+0.0009}_{-0.0012}$ around a spectroscopically-confirmed $UVJ$-quiescent ultra-massive galaxy (UMG; $M_{\star}=2.34^{+0.23}_{-0.34}\times10^{11} {\rm M}_\odot$) in the COSMOS UltraVISTA field. We present a total of 38 protocluster members (14 spectroscopic and 24 photometric), including the UMG. Notably, and in marked contrast to protoclusters previously reported at this epoch which have been found to contain predominantly star-forming members, we measure an elevated fraction of quiescent galaxies relative to the coeval field ($73.3^{+26.7}_{-16.9}\%$ versus $11.6^{+7.1}_{-4.9}\%$ for galaxies with stellar mass $M_{\star} \geq 10^{11} {\rm M}_\odot$). This high quenched fraction provides a striking and important counterexample to the seeming ubiquitousness of star-forming galaxies in protoclusters at $z>2$ and suggests, rather, that protoclusters exist in a diversity of evolutionary states in the early Universe. We discuss the possibility that we might be observing either "early mass quenching" or non-classical "environmental quenching." We also present the discovery of MAGAZ3NE J100028+023349, a second spectroscopically-confirmed protocluster, at a very similar redshift of $z = 3.3801^{+0.0213}_{-0.0281}$. We present a total of 20 protocluster members, 12 of which are photometric and 8 spectroscopic including a post-starburst UMG ($M_{\star}=2.95^{+0.21}_{-0.20}\times10^{11} {\rm M}_\odot$). Protoclusters MAGAZ3NE J0959 and MAGAZ3NE J1000 are separated by 18 arcminutes on the sky (35 comoving Mpc), in good agreement with predictions from simulations for the size of "Coma"-type cluster progenitors at this epoch. It is highly likely that the two UMGs are the progenitors of Brightest Cluster Galaxies (BCGs) seen in massive virialized clusters at lower redshift. △ Less

Submitted 7 October, 2021; v1 submitted 15 September, 2021; originally announced September 2021.

Comments: 19 pages, 6 figures, 3 tables. Accepted for publication in ApJ

CLASH-VLT: Abell~S1063. Cluster assembly history and spectroscopic catalogue

Authors: A. Mercurio, P. Rosati, A. Biviano, M. Annunziatella, M. Girardi, B. Sartoris, M. Nonino, M. Brescia, G. Riccio, C. Grillo, I. Balestra, G. B. Caminha, G. De Lucia, R. Gobat, S. Seitz, P. Tozzi, M. Scodeggio, E. Vanzella, G. Angora, P. Bergamini, S. Borgani, R. Demarco, M. Meneghetti, V. Strazzullo, L. Tortorelli , et al. (9 additional authors not shown)

Abstract: Using the CLASH-VLT survey, we assembled an unprecedented sample of 1234 spectroscopically confirmed members in Abell~S1063, finding a dynamically complex structure at z_cl=0.3457 with a velocity dispersion σ_v=1380 -32 +26 km s^-1. We investigate cluster environmental and dynamical effects by analysing the projected phase-space diagram and the orbits as a function of galaxy spectral properties. W… ▽ More Using the CLASH-VLT survey, we assembled an unprecedented sample of 1234 spectroscopically confirmed members in Abell~S1063, finding a dynamically complex structure at z_cl=0.3457 with a velocity dispersion σ_v=1380 -32 +26 km s^-1. We investigate cluster environmental and dynamical effects by analysing the projected phase-space diagram and the orbits as a function of galaxy spectral properties. We classify cluster galaxies according to the presence and strength of the [OII] emission line, the strength of the H$δ$ absorption line, and colours. We investigate the relationship between the spectral classes of galaxies and their position in the projected phase-space diagram. We analyse separately red and blue galaxy orbits. By correlating the observed positions and velocities with the projected phase-space constructed from simulations, we constrain the accretion redshift of galaxies with different spectral types. Passive galaxies are mainly located in the virialised region, while emission-line galaxies are outside r_200, and are accreted later into the cluster. Emission-lines and post-starbursts show an asymmetric distribution in projected phase-space within r_200, with the first being prominent at Delta_v/sigma <~-1.5$, and the second at Delta_v/ sigma >~ 1.5, suggesting that backsplash galaxies lie at large positive velocities. We find that low-mass passive galaxies are accreted in the cluster before the high-mass ones. This suggests that we observe as passives only the low-mass galaxies accreted early in the cluster as blue galaxies, that had the time to quench their star formation. We also find that red galaxies move on more radial orbits than blue galaxies. This can be explained if infalling galaxies can remain blue moving on tangential orbits. △ Less

Submitted 3 November, 2021; v1 submitted 7 September, 2021; originally announced September 2021.

Comments: 24 pages, 28 figures, 3 tables, Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 656, A147 (2021)

The Rapid Build-up of Massive Early-type Galaxies. Supersolar Metallicity, High Velocity Dispersion and Young Age for an ETG at z=3.35

Authors: Paolo Saracco, Danilo Marchesini, Francesco La Barbera, Adriana Gargiulo, Marianna Annunziatella, Ben Forrest, Daniel J. Lange Vagle, Z. Cemile Marsan, Adam Muzzin, Mauro Stefanon, Gillian Wilson

Abstract: Thanks to very deep spectroscopic observations carried out at the Large Binocular Telescope, we measured simultaneously stellar age, metallicity and velocity dispersion for C1-23152, an ETG at redshift $z$=3.352, corresponding to an epoch when the Universe was $\sim$1.8 Gyr old. The analysis of its spectrum shows that this galaxy, hosting an AGN, formed and assembled $\sim$2$\times$10$^{11}$ M… ▽ More Thanks to very deep spectroscopic observations carried out at the Large Binocular Telescope, we measured simultaneously stellar age, metallicity and velocity dispersion for C1-23152, an ETG at redshift $z$=3.352, corresponding to an epoch when the Universe was $\sim$1.8 Gyr old. The analysis of its spectrum shows that this galaxy, hosting an AGN, formed and assembled $\sim$2$\times$10$^{11}$ M$_\odot$ shaping its morphology within the $\sim$600 Myr preceding the observations, since $z$$\sim$4.6. The stellar population has a mean mass-weighted age 400$^{+30}_{-70}$ Myr and it is formed between $\sim$600 Myr and $\sim$150 Myr before the observed epoch, this latter being the time since quenching. Its high stellar velocity dispersion, $σ_e$=409$\pm$60 km s$^{-1}$, confirms the high mass (M$_{dyn}$=$2.2(\pm0.4)$$\times$10$^{11}$ M$_\odot$) and the high mass density ($Σ_e^{M^*}$=$Σ_{1kpc}=3.2(\pm0.7)\times10^{10}$ M$_\odot$ kpc$^{-2}$), suggesting a fast dissipative process at its origin. The analysis points toward a supersolar metallicity, [Z/H]=0.25$^{+0.006}_{-0.10}$, in agreement with the above picture, suggesting a star formation efficiency much higher than the replenishment time. However, sub-solar metallicity values cannot be firmly ruled out by our analysis. Quenching must have been extremely efficient to reduce the star formation to SFR$<$6.5 M$_\odot$ yr$^{-1}$ in less than 150 Myr. This could be explained by the presence of the AGN, even if a causal relation cannot be established from the data. C1-23152 has the same stellar and physical properties of the densest ETGs in the local Universe of comparable mass, suggesting that they are C1-23152-like galaxies which evolved to $z=0$ unperturbed. △ Less

Submitted 10 December, 2020; v1 submitted 9 November, 2020; originally announced November 2020.

Comments: 20 pages, 12 figures. Accepted for publication in ApJ (revised to match the ApJ version)

Journal ref: 2020, ApJ, 905, 40

The Number Densities and Stellar Populations of Massive Galaxies at 3 < z < 6: A Diverse, Rapidly Forming Population in the Early Universe

Authors: Z. Cemile Marsan, Adam Muzzin, Danilo Marchesini, Mauro Stefanon, Nicholas Martis, Marianna Annunziatella, Jeffrey C. C. Chan, Michael C. Cooper, Ben Forrest, Percy Gomez, Ian McConachie, Gillian Wilson

Abstract: We present the census of massive (log(M$_{*}$/M$_{\odot}$)$\geq 11$) galaxies at $3<z<6$ identified over the COSMOS/UltraVISTA Ultra-Deep field stripes: consisting of $\approx100$ and $\approx20$ high-confidence candidates at $3<z<4$ and $4<z<6$, respectively. The $3<z<4$ population is comprised of post-starburst, UV star-forming and dusty-star forming galaxies in roughly equal fractions, while UV… ▽ More We present the census of massive (log(M$_{*}$/M$_{\odot}$)$\geq 11$) galaxies at $3<z<6$ identified over the COSMOS/UltraVISTA Ultra-Deep field stripes: consisting of $\approx100$ and $\approx20$ high-confidence candidates at $3<z<4$ and $4<z<6$, respectively. The $3<z<4$ population is comprised of post-starburst, UV star-forming and dusty-star forming galaxies in roughly equal fractions, while UV-star-forming galaxies dominate at $4<z<6$ . We account for various sources of biases in SED modelling, finding that the treatment of emission line contamination is essential for understanding the number densities and mass growth histories of massive galaxies at $z>3$. The significant increase in observed number densities at $z\sim4$ ($>\times$ 5 in $\lesssim600$ Myrs) implies that this is the epoch at which log(M$_{*}$/M$_{\odot}$)$\geq 11$ galaxies emerge in significant numbers, with stellar ages ($\approx500-900$ Myrs) indicating rapid formation epochs as early as $z\sim7$. Leveraging ancillary multi-wavelength datasets, we perform panchromatic SED modelling to constrain the total star-formation activity of the sample. The star-formation activity of the sample is generally consistent with being on the star-formation main sequence at the considered redshifts, with $\approx15-25\%$ of the population showing evidence of suppressed star-formation rates, indicating that quenching mechanisms are already at play by $z\sim4$. We stack available HST imaging, confirming their compact nature ($r_{e}\lesssim2.2$ kpc), consistent with expected sizes of high-$z$ star-forming galaxies. Finally, we discuss how our results are in-line with the early formation epochs and short formation timescales inferred from the fossil records of the most massive galaxies in the Universe. △ Less

Submitted 9 October, 2020; originally announced October 2020.

Comments: 22 pages, 13 figures, 3 tables (31 pages, 16 figures, 4 tables with Appendix). ApJ submitted. Comments welcome!

The Massive Ancient Galaxies At $z>3$ NEar-infrared (MAGAZ3NE) Survey: Confirmation of Extremely Rapid Star-Formation and Quenching Timescales for Massive Galaxies in the Early Universe

Authors: Ben Forrest, Z. Cemile Marsan, Marianna Annunziatella, Gillian Wilson, Adam Muzzin, Danilo Marchesini, M. C. Cooper, Jeffrey C. C. Chan, Ian McConachie, Percy Gomez, Erin Kado-Fong, Francesco La Barbera, Daniel Lange-Vagle, Julie Nantais, Mario Nonino, Paolo Saracco, Mauro Stefanon, Remco F. J. van der Burg

Abstract: We present near-infrared spectroscopic confirmations of a sample of 16 photometrically-selected galaxies with stellar masses log(M_star/M_sun) > 11 at redshift z > 3 from the XMM-VIDEO and COSMOS-UltraVISTA fields using Keck/MOSFIRE as part of the MAGAZ3NE survey. Eight of the ultra-massive galaxies (UMGs) have specific star formation rates (sSFR) < 0.03 Gyr-1, with negligible emission lines. Anot… ▽ More We present near-infrared spectroscopic confirmations of a sample of 16 photometrically-selected galaxies with stellar masses log(M_star/M_sun) > 11 at redshift z > 3 from the XMM-VIDEO and COSMOS-UltraVISTA fields using Keck/MOSFIRE as part of the MAGAZ3NE survey. Eight of the ultra-massive galaxies (UMGs) have specific star formation rates (sSFR) < 0.03 Gyr-1, with negligible emission lines. Another seven UMGs show emission lines consistent with active galactic nuclei and/or star formation, while only one UMG has sSFR > 1 Gyr-1. Model star formation histories of these galaxies describe systems that formed the majority of their stars in vigorous bursts of several hundred Myr duration around 4 < z < 6during which hundreds to thousands of solar masses were formed per year. These formation ages of < 1 Gyr prior to observation are consistent with ages derived from measurements of Dn(4000) and EW0(Hδ). Rapid quenching followed these bursty star-forming periods, generally occurring less than 350 Myr before observation, resulting in post-starburst SEDs and spectra for half the sample. The rapid formation timescales are consistent with the extreme star formation rates observed in 4 < z < 7 dusty starbursts observed with ALMA, suggesting that such dusty galaxies are progenitors of these UMGs. While such formation histories have been suggested in previous studies, the large sample introduced here presents the most compelling evidence yet that vigorous star formation followed by rapid quenching is almost certainly the norm for high mass galaxies in the early universe. The UMGs presented here were selected to be brighter than Ks = 21.7 raising the intriguing possibility that even (fainter) older quiescent UMGs could exist at this epoch. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJ. 18 pages, 11 figures in main text

The deep Chandra survey in the SDSS J1030+0524 field

Authors: R. Nanni, R. Gilli, C. Vignali, M. Mignoli, A. Peca, S. Marchesi, M. Annunziatella, M. Brusa, F. Calura, N. Cappelluti, M. Chiaberge, A. Comastri, K. Iwasawa, G. Lanzuisi, E. Liuzzo, D. Marchesini, I. Prandoni, P. Tozzi, F. Vito, G. Zamorani, C. Norman

Abstract: We present the X-ray source catalog for the 479 ks Chandra exposure of the SDSS J1030+0524 field, that is centered on a region that shows the best evidence to date of an overdensity around a z > 6 quasar, and also includes a galaxy overdensity around a Compton-thick Fanaroff-Riley type II radio galaxy at z = 1.7. Using wavdetect for initial source detection and ACIS Extract for source photometry a… ▽ More We present the X-ray source catalog for the 479 ks Chandra exposure of the SDSS J1030+0524 field, that is centered on a region that shows the best evidence to date of an overdensity around a z > 6 quasar, and also includes a galaxy overdensity around a Compton-thick Fanaroff-Riley type II radio galaxy at z = 1.7. Using wavdetect for initial source detection and ACIS Extract for source photometry and significance assessment, we create preliminary catalogs of sources that are detected in the full, soft, and hard bands, respectively. We produce X-ray simulations that mirror our Chandra observation to filter our preliminary catalogs and get a completeness level of > 91% and a reliability level of 95% in each band. The catalogs in the three bands are then matched into a final main catalog of 256 unique sources. Among them, 244, 193, and 208 are detected in the full, soft, and hard bands, respectively. This makes J1030 field the fifth deepest extragalactic X-ray survey to date. The field is part of the Multiwavelength Survey by Yale-Chile (MUSYC), and is also covered by optical imaging data from the Large Binocular Camera (LBC) at the Large Binocular Telescope, near-IR imaging data from the Canada France Hawaii Telescope WIRCam, and Spitzer IRAC. Thanks to its dense multi-wavelength coverage, J1030 represents a legacy field for the study of large-scale structures around distant accreting supermassive black holes. Using a likelihood ratio analysis, we associate multi-band counterparts for 252 (98.4%) of the 256 Chandra sources, with an estimated reliability of 95%. Finally, we compute the cumulative number of sources in each X-ray band, finding that they are in general agreement with the results from the Chandra Deep Fields. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: A&A in press, 20 pages, 15 figures. The figures A.2 and and A.3 (postage-stamp images of the catalog) can be found online at: http://www.oabo.inaf.it/~LBTz6/1030/chandra_1030

The Fundamental Plane of cluster spheroidal galaxies at z$\sim1.3$. Evidence for mass-dependent evolution

Authors: P. Saracco, A. Gargiulo, F. La Barbera, M. Annunziatella, D. Marchesini

Abstract: We present spectroscopic observations obtained at the {\it Large Binocular Telescope} in the field of the cluster XLSSJ0223-0436 at $z=1.22$. We confirm 12 spheroids cluster members and determine stellar velocity dispersion for 7 of them. We combine these data with those in the literature for clusters RXJ0848+4453 at $z=1.27$ (8 galaxies) and XMMJ2235-2557 at $z=1.39$ (7 galaxies) to determine the… ▽ More We present spectroscopic observations obtained at the {\it Large Binocular Telescope} in the field of the cluster XLSSJ0223-0436 at $z=1.22$. We confirm 12 spheroids cluster members and determine stellar velocity dispersion for 7 of them. We combine these data with those in the literature for clusters RXJ0848+4453 at $z=1.27$ (8 galaxies) and XMMJ2235-2557 at $z=1.39$ (7 galaxies) to determine the Fundamental Plane of cluster spheroids. We find that the FP at $z\sim1.3$ is offset and { rotated ($\sim3σ$)} with respect to the local FP. The offset corresponds to a mean evolution $Δ$\rm{log}(M$_{dyn}$/L$_B$)=(-0.5$\pm$0.1)$z$. High-redshift galaxies follow a steeper mass-dependent M$_{dyn}$/L$_B$-M$_{dyn}$ relation than local ones. Assuming $Δ$ log$(M_{dyn}/L_B)$=$Δ$ log$(M^*/L_B)$, higher-mass galaxies (log(M$_{dyn}$/M$_\odot$)$\geq$11.5) have a higher-formation redshift ($z_f\geq$6.5) than lower-mass ones ($z_f\leq$2 for log(M$_{dyn}$/M$_\odot$$\leq$10)), with a median $z_f\simeq2.5$ for the whole sample. Also, galaxies with higher stellar mass density host stellar populations formed earlier than those in lower density galaxies. At fixed IMF, M$_{dyn}$/M$^*$ varies systematically with mass and mass density. It follows that the evolution of the stellar populations (M$^*/L_B$) accounts for the observed evolution of M$_{dyn}/L_B$ for M$_{dyn}$$>10^{11}$ M$_\odot$ galaxies, while accounts for $\sim$85\% of the evolution at M$_{dyn}$$<10^{11}$ M$_\odot$. We find no evidence in favour of structural evolution of individual galaxies, while we find evidences that spheroids later added to the population account for the observed discrepancy at masses $<10^{11}$ M$_\odot$. [Abridged] △ Less

Submitted 4 November, 2019; originally announced November 2019.

Comments: 19 pages (including appendices), 16 figures, 5 tables, accepted for publication in MNRAS

An Extremely Massive Quiescent Galaxy at $z=3.493$: Evidence of Insufficiently Rapid Quenching Mechanisms in Theoretical Models

Authors: Ben Forrest, Marianna Annunziatella, Gillian Wilson, Danilo Marchesini, Adam Muzzin, M. C. Cooper, Z. Cemile Marsan, Ian McConachie, Jeffrey C. C. Chan, Percy Gomez, Erin Kado-Fong, Francesco La Barbera, Ivo Labbé, Daniel Lange-Vagle, Julie Nantais, Mario Nonino, Theodore Peña, Paolo Saracco, Mauro Stefanon, Remco F. J. van der Burg

Abstract: We present spectra of the most massive quiescent galaxy yet discovered at $z>3$, spectroscopically confirmed via the detection of Balmer absorption features in the $H-$ and $K-$bands of Keck/MOSFIRE. The spectra confirm a galaxy with no significant ongoing star formation, consistent with the lack of rest-frame UV flux and overall photometric spectral energy distribution. With a stellar mass of… ▽ More We present spectra of the most massive quiescent galaxy yet discovered at $z>3$, spectroscopically confirmed via the detection of Balmer absorption features in the $H-$ and $K-$bands of Keck/MOSFIRE. The spectra confirm a galaxy with no significant ongoing star formation, consistent with the lack of rest-frame UV flux and overall photometric spectral energy distribution. With a stellar mass of $3.1^{+0.1}_{-0.2} \times 10^{11} ~\rm{M}_\odot$ at $z = 3.493$, this galaxy is nearly three times more massive than the highest redshift spectroscopically confirmed absorption-line identified galaxy known. The star-formation history of this quiescent galaxy implies that it formed $>1000 ~\rm{M}_\odot$/yr for almost 0.5 Gyr beginning at $z\sim7.2$, strongly suggestive that it is the descendant of massive dusty star-forming galaxies at $5<z<7$ recently observed with ALMA. While galaxies with similarly extreme stellar masses are reproduced in some simulations at early times, such a lack of ongoing star formation is not seen there. This suggests the need for a more rapid quenching process than is currently prescribed, challenging our current understanding of how ultra-massive galaxies form and evolve in the early Universe. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 6 pages of text + 5 figures + references. Submitted to ApJL