Search

We present a sample of 88 candidate z ∼ 8.5–14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science survey. These data cover ∼90 arcmin ^2 (10 NIRCam pointings) in six broadband imaging filters and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than predicted by most theoretical models. We construct the rest-frame ultraviolet luminosity functions at z ∼ 9, 11, and 14 and show that the space density of bright ( M _UV = −20) galaxies changes only modestly from z ∼ 14 to z ∼ 9, compared to a steeper increase from z ∼ 8 to z ∼ 4. While our candidates are photometrically selected, spectroscopic follow-up has now confirmed 13 of them, with only one significant interloper, implying that the fidelity of this sample is high. Successfully explaining the evidence for a flatter evolution in the number densities of UV-bright z > 10 galaxies may thus require changes to the dominant physical processes regulating star formation. While our results indicate that significant variations of dust attenuation with redshift are unlikely to be the dominant factor at these high redshifts, they are consistent with predictions from models that naturally have enhanced star formation efficiency and/or stochasticity. An evolving stellar initial mass function could also bring model predictions into better agreement with our results. Deep spectroscopic follow-up of a large sample of early galaxies can distinguish between these competing scenarios.

We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy ( ${f}_{\mathrm{lim},\mathrm{line},5\sigma }\approx 1.2\,\times \,$ 10 ^−18 erg s ^−1 cm ^−2 ) to measure metallicities and star formation rates (SFRs) for low-mass galaxies through the peak of the cosmic SFR density (0.5 < z < 4). In parallel, NGDEEP targets the HUDF-Par2 parallel field with NIRCam ( ${m}_{\mathrm{lim},5\sigma }=30.6-30.9$ ) to discover galaxies to z > 12, constraining the slope of the faint end of the rest-ultraviolet luminosity function. NGDEEP overlaps with the deepest HST Advanced Camera for Surveys optical imaging in the sky, F435W in the HUDF ( ${m}_{\mathrm{lim},{\rm{F}}435{\rm{W}}}=29.6$ ) and F814W in HUDF-Par2 ( ${m}_{\mathrm{lim},{\rm{F}}814{\rm{W}}}=30$ ), making this a premier HST+JWST deep field. As a treasury survey, NGDEEP data are public immediately, and we will rapidly release data products and catalogs in the spirit of previous deep-field initiatives. In this paper we present the NGDEEP survey design, summarize the science goals, and detail plans for the public release of NGDEEP reduced data products.

We present measurements of the rest-frame UV spectral slope, β , for a sample of 36 faint star-forming galaxies at z ∼ 9–16 discovered in one of the deepest JWST NIRCam surveys to date, the Next Generation Deep Extragalactic Exploratory Public Survey. We use robust photometric measurements for UV-faint galaxies (down to M _UV ∼ −16), originally published in Leung et al., and measure values of the UV spectral slope via photometric power-law fitting to both the observed photometry and stellar population models obtained through spectral energy distribution (SED) fitting with Bagpipes . We obtain a median and 68% confidence interval for β from photometric power-law fitting of ${\beta }_{\mathrm{PL}}=-{2.7}_{-0.5}^{+0.5}$ and from SED fitting, ${\beta }_{\mathrm{SED}}=-{2.3}_{-0.1}^{+0.2}$ for the full sample. We show that when only two to three photometric detections are available, SED fitting has a lower scatter and reduced biases than photometric power-law fitting. We quantify this bias and find that after correction the median ${\beta }_{\mathrm{SED},\mathrm{corr}}=-{2.5}_{-0.2}^{+0.2}$ . We measure physical properties for our galaxies with Bagpipes and find that our faint ( ${M}_{\mathrm{UV}}=-{18.1}_{-0.9}^{+0.7}$ ) sample is low in mass ( $\mathrm{log}[{M}_{* }/{M}_{\odot }]={7.7}_{-0.5}^{+0.5}$ ), fairly dust-poor ( ${A}_{{\rm{v}}}={0.1}_{-0.1}^{+0.2}$ mag), and modestly young ( $\mathrm{log}[\mathrm{age}]={7.8}_{-0.8}^{+0.2}$ yr) with a median star formation rate of $\mathrm{log}(\mathrm{SFR})=-{0.3}_{-0.4}^{+0.4}{M}_{\odot }\,{\mathrm{yr}}^{-1}$ . We find no strong evidence for ultrablue UV spectral slopes ( β ∼ −3) within our sample, as would be expected for exotically metal-poor ( Z / Z _⊙ < 10 ^−3 ) stellar populations with very high Lyman continuum escape fractions. Our observations are consistent with model predictions that galaxies of these stellar masses at z ∼ 9–16 should have only modestly low metallicities ( Z / Z _⊙ ∼ 0.1–0.2).

We report the identification of 15 galaxy candidates at z ≥ 9 using the initial COSMOS-Web JWST observations over 77 arcmin ^2 through four Near Infrared Camera filters (F115W, F150W, F277W, and F444W) with an overlap with the Mid-Infrared Imager (F770W) of 8.7 arcmin ^2 . We fit the sample using several publicly available spectral energy distribution (SED) fitting and photometric redshift codes and determine their redshifts between z = 9.3 and z = 10.9 (〈 z 〉 = 10.0), UV magnitudes between M _UV = −21.2 and −19.5 (with 〈 M _UV 〉 = −20.2), and rest-frame UV slopes (〈 β 〉 = −2.4). These galaxies are, on average, more luminous than most z ≥ 9 candidates discovered by JWST so far in the literature, while exhibiting similar blue colors in their rest-frame UV. The rest-frame UV slopes derived from SED fitting are blue ( β ∼ [−2.0, −2.7]) without reaching extremely blue values as reported in other recent studies at these redshifts. The blue color is consistent with models that suggest the underlying stellar population is not yet fully enriched in metals like similarly luminous galaxies in the lower-redshift Universe. The derived stellar masses with $\langle {\mathrm{log}}_{10}($ M _⋆ / M _⊙ )〉 ≈ 8–9 are not in tension with the standard Lambda cold dark matter (ΛCDM) model, and our measurement of the volume density of such UV-luminous galaxies aligns well with previously measured values presented in the literature at z ∼ 9–10. Our sample of galaxies, although compact, is significantly resolved.

We analyze the evolution of massive (log _10 [ M _⋆ / M _⊙ ] > 10) galaxies at z ∼ 1–4 selected from JWST Cosmic Evolution Early Release Survey (CEERS). We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting with dense basis to select a sample of high-redshift massive galaxies. Where available we include constraints from additional CEERS observing modes, including 18 sources with MIRI photometric coverage, and 28 sources with spectroscopic confirmations from NIRSpec or NIRCam WFSS. We sample the recovered posteriors in stellar mass from SED fitting to infer the volume densities of massive galaxies across cosmic time, taking into consideration the potential for sample contamination by active galactic nuclei. We find that the evolving abundance of massive galaxies tracks expectations based on a constant baryon conversion efficiency in dark matter halos for z ∼ 1–4. At higher redshifts, we observe an excess abundance of massive galaxies relative to this simple model, resulting in a shallower decline of observed volume densities of massive galaxies. These higher abundances can be explained by modest changes to star formation physics and/or the efficiencies with which star formation occurs in massive dark matter halos, and are not in tension with modern cosmology.

Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at z ≥ 3. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to be robust to noise and to correlate well with the physical properties of the simulated galaxies, including their 3D structure. We apply the method simultaneously to F200W and F356W galaxy images of a mass-complete sample ( M _* / M _⊙ > 10 ^9 ) at 3 ≤ z ≤ 6 from the first JWST/NIRCam CEERS data release. We find that the simulated and observed galaxies do not exactly populate the same manifold in the representation space from contrastive learning. We also find that half the galaxies classified as disks—either convolutional neural network-based or visually—populate a similar region of the representation space as TNG50 galaxies with low stellar specific angular momentum and nonoblate structure. Although our data-driven study does not allow us to firmly conclude on the true nature of these galaxies, it suggests that the disk fraction at z ≥ 3 remains uncertain and possibly overestimated by traditional supervised classifications. Deeper imaging and spectroscopic follow-ups as well as comparisons with other simulations will help to unambiguously determine the true nature of these galaxies, and establish more robust constraints on the emergence of disks at very high redshift.

We present a robust sample of very high redshift galaxy candidates from the first epoch of JWST/NIRCam imaging from the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey. The NGDEEP NIRCam imaging, spanning 9.7 arcmin ^2 in the Hubble Ultra Deep Field Parallel Field 2, reaches m = 30.4 (5 σ , point-source, 2″ diameter apertures corrected to total) in F277W, making it the deepest public JWST GO imaging data set to date. We describe our detailed data reduction process of the six-filter broadband JWST/NIRCam imaging, incorporating custom corrections for systematic effects to produce high-quality calibrated images. Using robust photometric redshift selection criteria, we identify a sample of 38 z ≳ 9 galaxy candidates. These objects span a redshift range of z = 8.5–15.8 and apparent magnitudes of m _F277W = 27–30.5 AB mag, reaching ∼1.5 mag deeper than previous public JWST imaging surveys. We calculate the rest-frame ultraviolet luminosity function at z ∼ 9 and 11 and present a new measurement of the luminosity function faint-end slope at z ∼ 11. We find a faint-end slope of α = −2.5 ± 0.4 and −2.2 ± 0.2 at z ∼ 9 and 11, respectively. This is consistent with no significant evolution in the faint-end slope and number density from z = 9 to 11. Comparing our results with theoretical predictions, we find that some models produce better agreement at the faint end than the bright end. These results will help to constrain how stellar feedback impacts star formation at these early epochs.

We present JWST/NIRSpec prism spectroscopy of seven galaxies selected from Cosmic Evolution Early Release Science (CEERS) survey NIRCam imaging with photometric redshifts z _phot > 8. We measure emission line redshifts of z = 7.65 and 8.64 for two galaxies. For two other sources without securely detected emission lines we measure $z={9.77}_{-0.29}^{+0.37}$ and ${10.01}_{-0.19}^{+0.14}$ by fitting model spectral templates to the prism data, from which we detect continuum breaks consistent with Ly α opacity from a mostly neutral intergalactic medium. The presence of strong breaks and the absence of strong emission lines give high confidence that these two galaxies have redshifts z > 9.6, but the redshift values derived from the breaks alone have large uncertainties given the low spectral resolution and relatively low S/N of the CEERS NIRSpec prism data. The two z ∼ 10 sources observed are relatively luminous ( M _UV < −20), with blue continua (−2.3 ≲ β ≲ −1.9) and low dust attenuation ( ${A}_{V}\simeq {0.15}_{-0.1}^{+0.3}$ ); and at least one of them has a high stellar mass for a galaxy at that redshift ( $\mathrm{log}({M}_{\star }/{M}_{\odot })\simeq {9.3}_{-0.3}^{+0.2}$ ). Considered together with spectroscopic observations of other CEERS NIRCam-selected high- z galaxy candidates in the literature, we find a high rate of redshift confirmation and low rate of confirmed interlopers (8%). Ten out of 35 z > 8 candidates with CEERS NIRSpec spectroscopy do not have secure redshifts, but the absence of emission lines in their spectra is consistent with redshifts z > 9.6. We find that z > 8 photometric redshifts are generally in agreement (within their uncertainties) with the spectroscopic values, but also that the photometric redshifts tend to be slightly overestimated (〈Δ z 〉 = 0.45 ± 0.11), suggesting that current templates do not fully describe the spectra of very-high- z sources. Overall, the spectroscopy solidifies photometric redshift evidence for a high spatial density of bright galaxies at z > 8 compared to theoretical model predictions, and further disfavors an accelerated decline in the integrated UV luminosity density at z > 8.

We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of z ≃ 9 − 13 and M _UV ∈ [ −21, −18] newly identified in NIRCam images in the Cosmic Evolution Early Release Science Survey. We confirm emission line redshifts for 7 galaxies at z = 7.762–8.998 using spectra at ∼1–5 μ m either with the NIRSpec prism or its three medium-resolution ( R ∼ 1000) gratings. For z ≃ 9 photometric candidates, we achieve a high confirmation rate of ≃90%, which validates the classical dropout selection from NIRCam photometry. No robust emission lines are identified in three galaxy candidates at z > 10, where the strong [O iii ] and H β lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Ly α continuum break is not detected with the sensitivity of the current data. Compared with Hubble Space Telescope-selected bright galaxies ( M _UV ≃ −22) that are similarly spectroscopically confirmed at z ≃ 8 − 9, these NIRCam-selected galaxies are characterized by lower star formation rates (SFRs; SFR ≃ 4 M _⊙ yr ^−1 ) and lower stellar masses (≃10 ^8 M _⊙ ), but with higher specific SFR (≃40 Gyr ^−1 ), higher [O iii ]+H β equivalent widths (≃1100 Å), and elevated production efficiency of ionizing photons ( $\mathrm{log}({\xi }_{\mathrm{ion}}/\mathrm{Hz}\,{\mathrm{erg}}^{-1})\simeq 25.8$ ) induced by young stellar populations (

We present results from the Cosmic Evolution Early Release Survey on the stellar population parameters for 28 galaxies with redshifts 4 < z < 9 using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 μ m data extend the coverage of the rest-frame spectral energy distribution to nearly 1 μ m for galaxies in this redshift range. By modeling the galaxies’ SEDs the MIRI data show that the galaxies have, on average, rest-frame UV (1600 Å)— I -band colors 0.4 mag bluer than derived when using photometry that lacks MIRI. Therefore, the galaxies have lower ratios of stellar mass to light. The MIRI data reduce the stellar masses by $\langle {\rm{\Delta }}\mathrm{log}{M}_{* }\rangle =0.25$ dex at 4 < z < 6 and 0.37 dex at 6 < z < 9. This also reduces the star formation rates (SFRs) by 〈ΔlogSFR〉 = 0.14 dex at 4 < z < 6 and 0.27 dex at 6 < z < 9. The MIRI data also improve constraints on the allowable stellar mass formed in early star formation. We model this using a star formation history that includes both a “burst” at z _f = 100 and a slowly varying (“delayed- τ ”) model. The MIRI data reduce the allowable stellar mass by 0.6 dex at 4 < z < 6 and by ≈1 dex at 6 < z < 9. Applying these results globally, this reduces the cosmic stellar-mass density by an order of magnitude in the early Universe ( z ≈ 9). Therefore, observations of rest-frame ≳1 μ m are paramount for constraining the stellar-mass buildup in galaxies at very high redshifts.

We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z ∼ 0.5 to >10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 arcmin ^2 , to search for candidate galaxies at z > 9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multiband photometry across seven NIRCam broad- and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z ∼ 9–16. These objects are compact with a median half-light radius of ∼0.5 kpc. We present an early estimate of the z ∼ 11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M _UV ∼ −20 appears to evolve very little from z ∼ 9 to 11. We also find that the abundance (surface density [arcmin ^−2 ]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z > 10, star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultrahigh-redshift galaxies with which to further explore these conclusions.

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 < 27.5 mag. The nature of these sources is threefold: (1) 71% are dusty star-forming galaxies (SFGs) at 2 < z < 6 with $9\lt \mathrm{log}{M}_{\star }/{M}_{\odot }\lt 11$ and a variety of specific SFRs (100 Gyr ^−1 ); (2) 18% are quiescent/dormant (i.e., subject to reignition/rejuvenation) galaxies (QGs) at 3 < z < 5, with $\mathrm{log}{M}_{\star }/{M}_{\odot }\sim 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- z 6) and $\mathrm{log}{M}_{\star }/{M}_{\odot }\sim 9.5$ . The sample is dominated by disk-like galaxies with remarkable compactness for XELG- z 6 (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 10 ^8 M _⊙ , 60% had assembled 10 ^9 M _⊙ , and 10% up to 10 ^10 M _⊙ (in situ or ex situ).

We report the discovery of a candidate galaxy with a photo- z of z ∼ 12 in the first epoch of the James Webb Space Telescope (JWST) Cosmic Evolution Early Release Science Survey. Following conservative selection criteria, we identify a source with a robust z _phot = ${11.8}_{-0.2}^{+0.3}$ (1 σ uncertainty) with m _F200W = 27.3 and ≳7 σ detections in five filters. The source is not detected at λ < 1.4 μ m in deep imaging from both Hubble Space Telescope (HST) and JWST and has faint ∼3 σ detections in JWST F150W and HST F160W, which signal a Ly α break near the red edge of both filters, implying z ∼ 12. This object (Maisie’s Galaxy) exhibits F115W − F200W > 1.9 mag (2 σ lower limit) with a blue continuum slope, resulting in 99.6% of the photo- z probability distribution function 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 _* / M _⊙ ∼ 8.5 and is highly star-forming (log sSFR ∼ −8.2 yr ^−1 ), with a blue rest-UV color ( β ∼ −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 that smoothly decline with increasing redshift. Should follow-up spectroscopy validate this redshift, our universe was already aglow with galaxies less than 400 Myr after the Big Bang.

We present a search for extremely red, dust-obscured, z > 7 galaxies with JWST/NIRCam+MIRI imaging over the first 20 arcmin ^2 of publicly available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W−F444W (∼2.5 mag) and detection in MIRI/F770W (∼25 mag), we identify two galaxies, COS-z8M1 and CEERS-z7M1, that have best-fit photometric redshifts of $z={8.4}_{-0.4}^{+0.3}$ and ${7.6}_{-0.1}^{+0.1}$ , respectively. We perform spectral energy distribution fitting with a variety of codes (including bagpipes , prospector , beagle , and cigale ) and find a >95% probability that these indeed lie at z > 7. Both sources are compact ( R _eff ≲ 200 pc) and highly obscured ( A _V ∼ 1.5–2.5) and, at our best-fit redshift estimates, likely have strong [O iii ]+H β emission contributing to their 4.4 μ m photometry. We estimate stellar masses of ∼10 ^10 M _⊙ for both sources; by virtue of detection in MIRI at 7.7 μ m, these measurements are robust to the inclusion of bright emission lines, for example, from an active galactic nucleus. We identify a marginal (2.9 σ ) Atacama Large Millimeter/submillimeter Array detection at 2 mm within 0.″5 of COS-z8M1, which, if real, would suggest a remarkably high IR luminosity of ∼10 ^12 L _⊙ . These two galaxies, if confirmed at z ∼ 8, would be extreme in their stellar and dust masses and may be representative of a substantial population of highly dust-obscured galaxies at cosmic dawn.

We present Atacama Large Millimeter/submillimeter Array (ALMA) Band 7 observations of a remarkably bright galaxy candidate at ${z}_{\mathrm{phot}}={16.7}_{-0.3}^{+1.9}$ ( M _UV = −21.6), S5-z17-1, identified in James Webb Space Telescope (JWST) Early Release Observation data of Stephen’s Quintet. We do not detect the dust continuum at 866 μ m, ruling out the possibility that S5-z17-1 is a low- z dusty starburst with a star formation rate of ≳30 M _⊙ yr ^−1 . We detect a 5.1 σ line feature at 338.726 ± 0.007 GHz exactly coinciding with the JWST source position, with a 2% likelihood of the signal being spurious. The most likely line identification would be [O iii ]52 μ m at z = 16.01 or [C ii ]158 μ m at z = 4.61, whose line luminosities do not violate the nondetection of the dust continuum in both cases. Together with three other z ≳ 11–13 candidate galaxies recently observed with ALMA, we conduct a joint ALMA and JWST spectral energy distribution (SED) analysis and find that the high- z solution at z ∼ 11–17 is favored in every candidate as a very blue (UV continuum slope of ≃−2.3) and luminous ( M _UV ≃ [ − 24:−21]) system. Still, we find in several candidates that reasonable SED fits (Δ χ ^2 ≲ 4) are reproduced by type II quasar and/or quiescent galaxy templates with strong emission lines at z ∼ 3–5, where such populations predicted from their luminosity functions and EW([O iii ]+H β ) distributions are abundant in survey volumes used for the identification of the z ∼ 11–17 candidates. While these recent ALMA observation results have strengthened the likelihood of the high- z solutions, lower- z possibilities are not completely ruled out in several of the z ∼ 11–17 candidates, indicating the need to consider the relative surface densities of the lower- z contaminants in the ultra-high- z galaxy search.

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 < z < 4 that evolve into a local M _* > 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.