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

© 2023 The Author(s). Artículo firmado por 23 autores. We sincerely thank the referee for the thorough reading and insightful comments that helped to improve the quality of this work. This research has been funded by grants PGC2018- 093499-B-I00 and RTI2018-096188-B-I00 funded by MCIN/ AEI/10.13039/501100011033. Á.G.A. acknowledges the support of the Universidad Complutense de Madrid through the predoctoral grant CT17/17-CT18/17. A.Y. is supported by an appointment to the NASA Postdoctoral Program (NPP) at NASA Goddard Space Flight Center, administered by Oak Ridge Associated Universities under contract with NASA. L.C. acknowledges financial support from Comunidad de Madrid under Atracción de Talento grant 2018-T2/TIC-11612. R.M. M. acknowledges support from Spanish Ministerio de Ciencia, Innovación y Universidades through grant PGC2018-093499- B-I00, from MDM-2017-0737 Unidad de Excelencia “Maria de Maeztu”-Centro de Astrobiología (INTA-CSIC) by the Spanish Ministry of Science and Innovation/State Agency of Research MCIN/AEI/ 10.13039/501100011033 and by “ERDF A way of making Europe,” and also from the Instituto Nacional de Técnica Aeroespacial SHARDS-JWST project through the PRE-SHARDSJWST/2020 grant.We use the Illustris-1 simulation to explore the capabilities of the Hubble Space Telescope (HST) and James Webb Space Telescope (JWST) 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 broadband ACS, WFC3, and NIRCam data and two-dimensional stellar population synthesis (2D-SPS) to derive the integrated star formation history (SFH) of massive (M_(*) > 10^(10) Mꙩ) simulated galaxies at 1 10^(11) Mꙩ galaxy. In particular, we explore the potential of HST and JWST data sets reaching a depth similar to those of the CANDELS and ongoing CEERS observations, respectively, and concentrate on determining the capabilities of this data set for characterizing the first episodes in the SFH of local M_(*) > 10^(11) Mꙩ 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%–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_(*) > 10^(11) Mꙩ galaxies would have started their formation by z = 16, forming the first 5% of their stellar mass present at z ∼ 1 by z = 4.5, 10% by z = 3.7, and 25% by z = 2.7.Depto. de Física de la Tierra y AstrofísicaFac. de Ciencias FísicasTRUEMinisterio de Ciencia e Innovación (MCIN)Universidad Complutense de Madrid through the predoctoral grantNASA Postdoctoral Program (NPP) at NASA Goddard Space Flight CenterComunidad de Madrid under Atracción de Talento grantUnidad de Excelencia “Maria de Maeztu”Agencia Estatal de Investigación (España)Instituto Nacional de Técnica Aeroespacial SHARDS-JWSTERDF "A way of making Europe"pu

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

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 ($\Sigma = \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.Comment: Accepted for publication in MNRA

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

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 $\sim28.2$ mag, $\sim27.9$ mag), blue (UV-slope $\sim-1.9$), star-forming (rest-frame colors $\sim0.10$ mag, $\sim0.17$ mag) galaxies. These observational characteristics are identified with young (mass-weighted age $\sim0.014$ Gyr) stellar populations subject to low attenuations ($\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.Comment: 34 pages, 16 figures, 9 tables. Accepted for publication to Ap

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

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 $110^{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$.Comment: 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 update

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

We present the ultraviolet luminosity function and an estimate of the cosmic star formation rate density at 8 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 ∼ 13 to z ∼ 8, i.e., only a small increase in the average number density of ∼0.2 dex, while the faint-end slope and absolute magnitude of the knee remain approximately constant, with values α = − 2.2 ± 0.1, and M* = − 20.8 ± 0.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 ∼108M⊙ during the first 350 Myr of the universe, z ∼ 12, with models matching better the luminosity density observational estimations ∼150 Myr later, by z ∼ 9

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

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 (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)Comment: Published in CEERS ApJL Focus Issue, ApJL 946, L1

CEERS key paper. IV. A triality in the nature of HST-dark galaxies

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, with HST-dark galaxies among them. We gather JWST data from the CEERS survey in the Extended Groth Strip, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions to estimate photometric redshifts in two dimensions and stellar population properties on a pixel-by-pixel basis for red galaxies detected by NIRCam. We select 138 galaxies with F150W - F356W > 1.5 mag and F356W 100 Gyr-1); (2) 18% are quiescent/dormant (i.e., subject to reignition/rejuvenation) galaxies (QGs) at 3 < z < 5, with log M ? / M ? ~ 10 and poststarburst mass-weighted ages (0.5-1.0 Gyr); and (3) 11% are strong young starbursts with indications of high equivalent width emission lines (typically, [O iii]+Hß) at 6 < z < 7 (XELG-z6) and log M ? / M ? ~ 9.5 . The sample is dominated by disk-like galaxies with 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 submillimeter 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 108 M ?, 60% had assembled 109 M ?, and 10% up to 1010 M ? (in situ or ex situ)

Life beyond 30: probing the −20<MUV<−17-20<M_\mathrm{UV}<-17 luminosity function at 8<z<138<z<13 with the NIRCam parallel field of the MIRI Deep Survey

International audienceWe present the ultraviolet luminosity function and an estimate of the cosmic star formation rate density at $88$ 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 $\chi^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 $\alpha=-2.3\pm0.2$ and $M^*=-20.8\pm0.2$~mag. Comparing our results with the predictions of a wide range 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 observations $\sim$150~Myr later, by $z\sim9$

Life beyond 30: Probing the −20 < MUV < −17 Luminosity Function at 8 < z < 13 with the NIRCam Parallel Field of the MIRI Deep Survey

We present the ultraviolet luminosity function and an estimate of the cosmic star formation rate density at 8 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 ∼ 13 to z ∼ 8, i.e., only a small increase in the average number density of ∼0.2 dex, while the faint-end slope and absolute magnitude of the knee remain approximately constant, with values α = − 2.2 ± 0.1, and M * = − 20.8 ± 0.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 ∼108 M ⊙ during the first 350 Myr of the universe, z ∼ 12, with models matching better the luminosity density observational estimations ∼150 Myr later, by z ∼ 9.</p