Search

We present the properties of two candidate massive ( M _⋆ ∼ 10 ^11 M _⊙ ) and dusty ( A _v > 2.5 mag) galaxies at z = 5–7 in the first 0.28 deg ^2 of the COSMOS-Web survey. One object is spectroscopically confirmed at z _spec = 5.051, while the other has a robust z _phot = 6.7 ± 0.3. Thanks to their extremely red colors (F277W–F444W ∼ 1.7 mag), these galaxies satisfy the nominal color selection for the widely studied “little red dot” (LRD) population with the exception of their spatially resolved morphologies. The morphology of our targets allows us to conclude that their red continuum is dominated by highly obscured stellar emission and not by reddened nuclear activity. Using a variety of spectral energy distribution fitting tools and star formation histories, we estimate the stellar masses to be $\mathrm{log}({M}_{\star })={11.32}_{-0.15}^{+0.07}\,{M}_{\odot }$ and $\mathrm{log}({M}_{\star })={11.2}_{-0.2}^{+0.1}\,{M}_{\odot }$ , respectively, with a red continuum emission dominated by a recent episode of star formation. We then compare their number density to the halo mass function to infer stellar baryon fractions of ϵ _⋆ ∼ 0.25 and ϵ _⋆ ∼ 0.5. Both are significantly higher than what is commonly observed in lower- z galaxies or more dust-obscured galaxies at similar redshifts. With very bright ultra-high- z Lyman-Break Galaxies and some non-AGN-dominated LRDs, such “extended” LRDs represent another population that may require very efficient star formation at early times.

The identification of red, apparently massive galaxies at z > 7 in early James Webb Space Telescope (JWST) photometry suggests a strongly accelerated time line compared to standard models of galaxy growth. A major uncertainty in the interpretation is whether the red colors are caused by evolved stellar populations, dust, or other effects such as emission lines or active galactic nuclei (AGNs). Here we show that three of the massive galaxy candidates at z = 6.7–8.4 have prominent Balmer breaks in JWST/NIRSpec spectroscopy from the RUBIES program. The Balmer breaks demonstrate unambiguously that stellar emission dominates at λ _rest = 0.4 μ m and require formation histories extending hundreds of millions of years into the past in galaxies only 600–800 Myr after the big bang. Two of the three galaxies also show broad Balmer lines, with H β FWHM > 2500 km s ^−1 , suggesting that dust-reddened AGNs contribute to, or even dominate, the spectral energy distributions of these galaxies at λ _rest ≳ 0.6 μ m. All three galaxies have relatively narrow [O iii ] lines, seemingly ruling out a high-mass interpretation if the lines arise in dynamically relaxed, inclined disks. Yet the inferred masses also remain highly uncertain. We model the high-quality spectra using Prospector to decompose the continuum into stellar and AGN components and explore limiting cases in stellar/AGN contribution. This produces a wide range of possible stellar masses, spanning M _⋆ ∼ 10 ^9 −10 ^11 M _⊙ . Nevertheless, all fits suggest a very early and rapid formation, most of which follow with a truncation in star formation. Potential origins and evolutionary tracks for these objects are discussed, from the cores of massive galaxies to low-mass galaxies with overmassive black holes. Intriguingly, we find all of these explanations to be incomplete; deeper and redder data are needed to understand the physics of these systems.

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.

Many quiescent galaxies discovered in the early Universe by JWST raise fundamental questions on when and how these galaxies became and stayed quenched. Making use of the latest version of the semianalytic model GAEA that provides good agreement with the observed quenched fractions up to z ∼ 3, we make predictions for the expected fractions of quiescent galaxies up to z ∼ 7 and analyze the main quenching mechanism. We find that in a simulated box of 685 Mpc on a side, the first quenched massive ( M _⋆ ∼ 10 ^11 M _⊙ ), Milky Way–mass, and low-mass ( M _⋆ ∼ 10 ^9.5 M _⊙ ) galaxies appear at z ∼ 4.5, z ∼ 6.2, and before z = 7, respectively. Most quenched galaxies identified at early redshifts remain quenched for more than 1 Gyr. Independently of galaxy stellar mass, the dominant quenching mechanism at high redshift is accretion disk feedback (quasar winds) from a central massive black hole, which is triggered by mergers in massive and Milky Way–mass galaxies and by disk instabilities in low-mass galaxies. Environmental stripping becomes increasingly more important at lower redshift.

In the past, high- z active galactic nuclei (AGNs) were given a minor role as possible drivers of reionization, despite initial evidence in favor of their large space densities at low luminosities by Chandra and the Hubble Space Telescope. Recent observations from JWST are finding relatively large numbers of faint AGNs at z > 4, convincingly confirming these early results. We present a sample of z ∼ 5 AGNs, both from wide, shallow ground-based surveys and from deep, pencil-beam observations from JWST, allowing us to estimate their space densities with unprecedented accuracy. The bright end ( M _1450 < −26) of the z ∼ 5 AGN luminosity function is well constrained, with a rather steep slope. The faint end ( M _1450 ≥ −22) indicates a high space density, the scatter is significant, and the knee ( M _1450 ∼ −24) is mostly undetermined. Comparisons with state-of-the-art models find reasonable agreement with the observed AGN luminosity function at z = 5, while the predicted space density evolution at higher redshifts appears to be too fast with respect to observational constraints. Given the large variance at the faint end, we consider different options in fitting the luminosity functions and deriving the ionizing emissivity. Even in the most conservative scenario, the photoionization rate produced by z ∼ 5 AGNs is consistent with the ultraviolet background measurements. A slow evolution of the space density of faint AGNs is observed, indicating that active SMBHs are probably producing large amounts of ionizing photons at z > 6, well into the Epoch of Reionization. This is an important indication that high- z AGNs could be major contributors to the reionization of the Universe.

We present a sample of 1165 extreme emission-line galaxies (EELGs) at 4 < z < 9 selected using James Webb Space Telescope (JWST) NIRCam photometry in the Cosmic Evolution Early Release Science (CEERS) program. We use a simple method to photometrically identify EELGs with H β + [O iii ] (combined) or H α emission of observed-frame equivalent width (EW) > 5000 Å. JWST/NIRSpec spectroscopic observations of a subset (34) of the photometrically selected EELGs validate our selection method: All spectroscopically observed EELGs confirm our photometric identification of extreme emission, including some cases where the spectral-energy-distribution-derived photometric redshifts are incorrect. We find that the medium-band F410M filter in CEERS is particularly efficient at identifying EELGs, both in terms of including emission lines in the filter and in correctly identifying the continuum between H β + [O iii ] and H α in the neighboring broadband filters. We present examples of EELGs that could be incorrectly classified as ultrahigh redshift ( z > 12) as a result of extreme H β + [O iii ] emission blended across the reddest photometric filters. We compare the EELGs to the broader (subextreme) galaxy population in the same redshift range and find that they are consistent with being the bluer, high-EW tail of a broader population of emission-line galaxies. The highest-EW EELGs tend to have more compact emission-line sizes than continuum sizes, suggesting that active galactic nuclei are responsible for at least some of the most extreme EELGs. The photometrically inferred emission-line ratios are consistent with interstellar medium conditions with high ionization and moderately low metallicity, consistent with previous spectroscopic studies.

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.

Ultraviolet (UV; rest-frame ∼1200–2000 Å) spectra provide a wealth of diagnostics to characterize fundamental galaxy properties, such as their chemical enrichment, the nature of their stellar populations, and their amount of Lyman-continuum (LyC) radiation. In this work, we leverage publicly released JWST data to construct the rest-frame UV-to-optical composite spectrum of a sample of 63 galaxies at 5.6 < z < 9, spanning a wavelength range from 1500 to 5200 Å. Based on the composite spectrum, we derive an average dust attenuation $E{(B-V)}_{\mathrm{gas}}\,=\,{0.10}_{-0.11}^{+0.10}$ from H β /H γ , an electron density ${n}_{e}={570}_{-290}^{+510}$ cm ^−3 from the [O ii ] doublet ratio, an electron temperature ${T}_{e}={16700}_{-1500}^{+1500}$ K from the [O iii ] λ 4363/[O iii ] λ 5007 ratio, and an ionization parameter $\mathrm{log}(U)=-{2.15}_{-0.03}^{+0.03}$ from the [O iii ]/[O ii ] ratio. Using a direct T _e method, we calculate an oxygen abundance $12+\mathrm{log}({\rm{O}}/{\rm{H}})\,=\,7.67\pm 0.08$ and a carbon-to-oxygen (C/O) abundance ratio $\mathrm{log}({\rm{C}}/{\rm{O}})\,=\,-{0.86}_{-0.10}^{+0.13}$ . This C/O ratio is smaller than compared to z = 0 and z = 2–4 star-forming galaxies, albeit with moderate significance. This indicates the reionization-era galaxies might be undergoing a rapid buildup of stellar mass with high specific star formation rates. A UV diagnostic based on the ratios of C iii ] λ λ 1907, 1909/He ii λ 1640 versus O iii ] λ 1666/He ii λ 1640 suggests that the star formation is the dominant source of ionization, similar to the local extreme dwarf galaxies and z ∼ 2–4 He ii –detected galaxies. The [O iii ]/[O ii ] and C iv /C iii ] ratios of the composite spectrum are marginally larger than the criteria used to select galaxies as LyC leakers, suggesting that some of the galaxies in our sample are strong contributors to the reionizing radiation.

Substantial populations of massive quiescent galaxies at z ≥ 3 challenge our understanding of rapid galaxy growth and quenching over short timescales. In order to piece together this evolutionary puzzle, more statistical samples of these objects are required. Established techniques for identifying massive quiescent galaxies are increasingly inefficient and unconstrained at z > 3. As a result, studies report that as much as 70% of quiescent galaxies at z > 3 may be missed from existing surveys. In this work, we propose a new empirical color selection technique designed to select massive quiescent galaxies at 3 ≲ z ≲ 6 using JWST NIRCam imaging data. We use empirically constrained galaxy spectral energy distribution (SED) templates to define a region in the F277W − F444W versus F150W − F277W color plane that captures quiescent galaxies at z > 3. We apply these color selection criteria to the Cosmic Evolution Early Release Science (CEERS) Survey and use SED fitting on sources in the region to identify 44 candidate z ≳ 3 quiescent galaxies. Over half of these sources are newly discovered and, on average, exhibit specific star formation rates of poststarburst galaxies. Most of these sources would not be discovered using canonical UVJ diagrams. We derive volume density estimates of n ∼ 1–4 × 10 ^−5 Mpc ^−3 at 3 < z < 5, finding excellent agreement with existing reports on similar populations in the CEERS field. Thanks to NIRCam’s wavelength coverage and sensitivity, this technique provides an efficient tool to search for large samples of these rare galaxies.

We analyze rest-frame ultraviolet to optical spectra of three z ≃ 7.47–7.75 galaxies whose Ly α emission lines were previously detected with Keck/MOSFIRE observations, using the JWST/NIRSpec observations from the Cosmic Evolution Early Release Science survey. From NIRSpec data, we confirm the systemic redshifts of these Ly α emitters, and emission-line ratio diagnostics indicate these galaxies were highly ionized and metal-poor. We investigate Ly α line properties, including the line flux, velocity offset, and spatial extent. For the one galaxy where we have both NIRSpec and MOSFIRE measurements, we find a significant offset in their flux measurements (∼1.3–5× greater in MOSFIRE) and a marginal difference in the velocity shifts. The simplest interpretation is that the Ly α emission is extended and not entirely encompassed by the NIRSpec slit. The cross-dispersion profiles in NIRSpec reveal that Ly α in one galaxy is significantly more extended than the nonresonant emission lines. We also compute the expected sizes of ionized bubbles that can be generated by the Ly α sources and discuss viable scenarios for the creation of sizable ionized bubbles (>1 physical Mpc). The source with the highest-ionization condition is possibly capable of ionizing its own bubble, while the other two do not appear to be capable of ionizing such a large region, but require additional sources of ionizing photons. Therefore, the fact that we detect Ly α from these galaxies suggests diverse scenarios for the escape of Ly α during the epoch of reionization. High-spectral-resolution spectra with JWST/NIRSpec will be extremely useful for constraining the physics of patchy reionization.

We report the discovery of 15 exceptionally luminous 10 ≲ z ≲ 14 candidate galaxies discovered in the first 0.28 deg ^2 of JWST/NIRCam imaging from the COSMOS-Web survey. These sources span rest-frame UV magnitudes of −20.5 > M _UV > −22, and thus constitute the most intrinsically luminous z ≳ 10 candidates identified by JWST to date. Selected via NIRCam imaging, deep ground-based observations corroborate their detection and help significantly constrain their photometric redshifts. We analyze their spectral energy distributions using multiple open-source codes and evaluate the probability of low-redshift solutions; we conclude that 12/15 (80%) are likely genuine z ≳ 10 sources and 3/15 (20%) likely low-redshift contaminants. Three of our z ∼ 12 candidates push the limits of early stellar mass assembly: they have estimated stellar masses ∼ 5 × 10 ^9 M _⊙ , implying an effective stellar baryon fraction of ϵ _⋆ ∼ 0.2−0.5, where ϵ _⋆ ≡ M _⋆ /( f _b M _halo ). The assembly of such stellar reservoirs is made possible due to rapid, burst-driven star formation on timescales < 100 Myr where the star formation rate may far outpace the growth of the underlying dark matter halos. This is supported by the similar volume densities inferred for M _⋆ ∼ 10 ^10 M _⊙ galaxies relative to M _⋆ ∼ 10 ^9 M _⊙ —both about 10 ^−6 Mpc ^−3 —implying they live in halos of comparable mass. At such high redshifts, the duty cycle for starbursts would be of order unity, which could cause the observed change in the shape of the UV luminosity function from a double power law to a Schechter function at z ≈ 8. Spectroscopic redshift confirmation and ensuing constraints of their masses will be critical to understand how, and if, such early massive galaxies push the limits of galaxy formation in the Lambda cold dark matter paradigm.

We investigate the fueling mechanisms of supermassive black holes (SMBHs) by analyzing 10 zoom-in cosmological simulations of massive galaxies, with stellar masses 10 ^11–12 M _⊙ and SMBH masses 10 ^8.9–9.7 M _⊙ at z = 0, featuring various major and minor merger events. By tracing the gas history in these simulations, we categorize the gas accreted by the central SMBHs based on its origin. Gas that belonged to a different galaxy before accretion onto the BH is labeled as (i) “external,” while smoothly accreted cosmic gas is classified as (ii) “smooth.” Gas produced within the primary halo through stellar evolution and subsequently accreted by the SMBH is classified as (iii) “recycled.” Our analysis, which includes stellar feedback, reveals that the primary fuel source for SMBHs is the recycled gas from dying stars. This recycled gas from stars in the inner region of the galaxy readily collapses toward the center, triggering starbursts and simultaneously fueling the SMBH. Galaxy mergers also play a crucial role in fueling SMBHs in massive galaxies, as SMBHs in massive halos tend to accrete a higher fraction of external gas from mergers compared to smoothly accreted gas. However, on average, it takes approximately 1.85 Gyr for external gas to enter the main galaxy and accrete onto the SMBH. Considering the presence of various other gas triggers for active galactic nucleus (AGN) activity alongside this time delay, the association between AGNs and mergers may not always be obvious.

We study a new population of extremely red objects (EROs) recently discovered by the James Webb Space Telescope (JWST) based on their NIRCam colors F277W − F444W > 1.5 mag. We find 37 EROs in the Cosmic Evolution Early Release Science Survey (CEERS) field with F444W < 28 mag and photometric redshifts between 5 < z < 7, with median $z={6.9}_{-1.6}^{+1.0}$ . Surprisingly, despite their red long-wavelength colors, these EROs have blue short-wavelength colors (F150W − F200W ∼ 0 mag) indicative of bimodal spectral energy distributions (SEDs) with a red, steep slope in the rest-frame optical, and a blue, flat slope in the rest-frame UV. Moreover, all these EROs are unresolved, point-like sources in all NIRCam bands. We analyze the SEDs of eight of them with MIRI and NIRSpec observations using stellar population models and active galactic nucleus (AGN) templates. We find that dusty galaxies or obscured AGNs provide similarly good SED fits but different stellar properties: massive and dusty, $\mathrm{log}{M}_{\star }/{M}_{\odot }$ ∼ 10 and A _V ≳ 3 mag, or low mass and obscured, $\mathrm{log}{M}_{\star }/{M}_{\odot }$ ∼ 7.5 and A _V ∼ 0 mag, hosting an obscured quasi-stellar object (QSO). SED modeling does not favor either scenario, but their unresolved sizes are more suggestive of AGNs. If any EROs are confirmed to have $\mathrm{log}{M}_{\star }/{M}_{\odot }$ ≳ 10.5, it would increase the pre-JWST number density at z > 7 by up to a factor ∼60. Similarly, if they are QSOs with luminosities in the L _bol > 10 ^45–46 erg s ^−1 range, their number would exceed that of bright blue QSOs by more than three orders of magnitude. Additional photometry at mid-infrared wavelengths will reveal the true nature of the red continuum emission in these EROs and will place this puzzling population in the right context of galaxy evolution.

We use James Webb Space Telescope Near-Infrared Camera Wide Field Slitless Spectroscopy (NIRCam WFSS) and the Near-Infrared spectrograph (NIRSpec) in the Cosmic Evolution Early Release survey to measure rest-frame optical emission-line ratios of 155 galaxies at z > 2. The blind NIRCam grism observations include a sample of galaxies with bright emission lines that were not observed on the NIRSpec masks. We study the changes of the H α , [O III ]/H β , and [Ne III ]/[O II ] emission lines in terms of redshift by comparing to lower-redshift SDSS, CLEAR, and MOSDEF samples. We find a significant (>3 σ ) correlation between [O III ]/H β with redshift, while [Ne III ]/[O II ] has a marginal (2 σ ) correlation with redshift. We compare [O III ]/H β and [Ne III ]/[O II ] to stellar mass and H β SFR. We find that both emission-line ratios have a correlation with H β SFR and an anticorrelation with stellar mass across the redshifts 0 < z < 9. Comparison with MAPPINGS V models indicates that these trends are consistent with lower metallicity and higher ionization in low-mass and high-SFR galaxies. We additionally compare to IllustrisTNG predictions and find that they effectively describe the highest [O III ]/H β ratios observed in our sample, without the need to invoke MAPPINGS models with significant shock ionization components.

We combine deep imaging data from the CEERS early release JWST survey and Hubble Space Telescope imaging from CANDELS to examine the size–mass relation of star-forming galaxies and the morphology–quenching relation at stellar masses M _⋆ ≥ 10 ^9.5 M _⊙ over the redshift range 0.5 < z < 5.5. In this study with a sample of 2450 galaxies, we separate star-forming and quiescent galaxies based on their star formation activity and confirm that star-forming and quiescent galaxies have different morphologies out to z = 5.5, extending the results of earlier studies out to higher redshifts. We find that star-forming and quiescent galaxies have typical Sérsic indices of n ∼ 1.3 and n ∼ 4.3, respectively. Focusing on star-forming galaxies, we find that the slope of the size–mass relation is nearly constant with redshift, as was found previously, but shows a modest increase at z ∼ 4.2. The intercept in the size–mass relation declines out to z = 5.5 at rates that are similar to what earlier studies found. The intrinsic scatter in the size–mass relation is relatively constant out to z = 5.5.

The Cosmic Evolution Early Release Science program observed the Extended Groth Strip (EGS) with the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) in 2022. In this paper, we discuss the four MIRI pointings that observed with longer-wavelength filters, including F770W, F1000W, F1280W, F1500W, F1800W, and F2100W. We compare the MIRI galaxies with the Spitzer/MIPS 24 μ m population in the EGS field. We find that MIRI can observe an order of magnitude deeper than MIPS in significantly shorter integration times, attributable to JWST's much larger aperture and MIRI’s improved sensitivity. MIRI is exceptionally good at finding faint ( L _IR < 10 ^10 L _⊙ ) galaxies at z ∼ 1–2. We find that a significant portion of MIRI galaxies are “mid-IR weak”—they have strong near-IR emission and relatively weaker mid-IR emission, and most of the star formation is unobscured. We present new IR templates that capture how the mid-to-near-IR emission changes with increasing infrared luminosity. We present two color–color diagrams to separate mid-IR weak galaxies and active galactic nuclei (AGN) from dusty star-forming galaxies and find that these color diagrams are most effective when used in conjunction with each other. We present the first number counts of 10 μ m sources and find that there are ≲10 IR AGN per MIRI pointing, possibly due to the difficulty of distinguishing AGN from intrinsically mid-IR weak galaxies (due to low metallicities or dust content). We conclude that MIRI is most effective at observing moderate-luminosity ( L _IR = 10 ^9 –10 ^10 L _⊙ ) galaxies at z = 1–2, and that photometry alone is not effective at identifying AGN within this faint population.

Extreme emission line galaxies (EELGs), where nebular emissions contribute 30%–40% of the flux in certain photometric bands, are ubiquitous in the early Universe ( z > 6). We utilize deep NIRCam imaging from the JWST Advanced Deep Extragalactic Survey ( JADES ) to investigate the properties of companion galaxies (projected distance 6.

We report on the discovery of two low-luminosity, broad-line active galactic nuclei (AGNs) at z > 5 identified using JWST NIRSpec spectroscopy from the Cosmic Evolution Early Release Science (CEERS) survey. We detect broad H α emission in the spectra of both sources, with FWHM of 2060 ± 290 km s ^−1 and 1800 ± 200 km s ^−1 , resulting in virial black hole (BH) masses that are 1–2 dex below those of existing samples of luminous quasars at z > 5. The first source, CEERS 2782 at z = 5.242, is 2–3 dex fainter than known quasars at similar redshifts and was previously identified as a candidate low-luminosity AGN based on its morphology and rest-frame optical spectral energy distribution (SED). We measure a BH mass of M _BH = (1.3 ± 0.4) × 10 ^7 M _⊙ , confirming that this AGN is powered by the least massive BH known in the Universe at the end of cosmic reionization. The second source, CEERS 746 at z = 5.624, is inferred to be a heavily obscured, broad-line AGN caught in a transition phase between a dust-obscured starburst and an unobscured quasar. We estimate its BH mass to be in the range of M _BH ≃ (0.9–4.7) × 10 ^7 M _⊙ , depending on the level of dust obscuration assumed. We perform SED fitting to derive host stellar masses, M _⋆ , allowing us to place constraints on the BH–galaxy mass relationship in the lowest mass range yet probed in the early Universe. The M _BH / M _⋆ ratio for CEERS 2782, in particular, is consistent with or higher than the empirical relationship seen in massive galaxies at z = 0. We examine the narrow emission line ratios of both sources and find that their location on the BPT and OHNO diagrams is consistent with model predictions for moderately low metallicity AGNs with Z / Z _⊙ ≃ 0.2–0.4. The spectroscopic identification of low-luminosity, broad-line AGNs at z > 5 with M _BH ≃ 10 ^7 M _⊙ demonstrates the capability of JWST to push BH masses closer to the range predicted for the BH seed population and provides a unique opportunity to study the early stages of BH–galaxy assembly.

We report the discovery of an accreting supermassive black hole at z = 8.679. This galaxy, denoted here as CEERS_1019, was previously discovered as a Ly α -break galaxy by Hubble with a Ly α redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we have observed this source with JWST/NIRSpec, MIRI, NIRCam, and NIRCam/WFSS and uncovered a plethora of emission lines. The H β line is best fit by a narrow plus a broad component, where the latter is measured at 2.5 σ with an FWHM ∼1200 km s ^−1 . We conclude this originates in the broadline region of an active galactic nucleus (AGN). This is supported by the presence of weak high-ionization lines (N V, N IV], and C III]), as well as a spatial point-source component. The implied mass of the black hole (BH) is log ( M _BH / M _⊙ ) = 6.95 ± 0.37, and we estimate that it is accreting at 1.2 ± 0.5 times the Eddington limit. The 1–8 μ m photometric spectral energy distribution shows a continuum dominated by starlight and constrains the host galaxy to be massive (log M/M _⊙ ∼9.5) and highly star-forming (star formation rate, or SFR ∼ 30 M _⊙ yr ^−1 ; log sSFR ∼ − 7.9 yr ^−1 ). The line ratios show that the gas is metal-poor ( Z / Z _⊙ ∼ 0.1), dense ( n _e ∼ 10 ^3 cm ^−3 ), and highly ionized (log U ∼ − 2.1). We use this present highest-redshift AGN discovery to place constraints on BH seeding models and find that a combination of either super-Eddington accretion from stellar seeds or Eddington accretion from very massive BH seeds is required to form this object.

A large fraction of the accreting supermassive black hole population is shrouded by copious amounts of gas and dust, particularly in the distant ( z ≳ 1) universe. While much of the obscuration is attributed to a parsec-scale torus, there is a known contribution from the larger-scale host galaxy. Using JWST/NIRCam imaging from the COSMOS-Web survey, we probe the galaxy-wide dust distribution in X-ray selected active galactic nuclei (AGNs) up to z ∼ 2. Here, we focus on a sample of three AGNs with their host galaxies exhibiting prominent dust lanes, potentially due to their edge-on alignment. These represent 27% (3 out of 11 with early NIRCam data) of the heavily obscured ( N _H > 10 ^23 cm ^−2 ) AGN population. With limited signs of a central AGN in the optical and near-infrared, the NIRCam images are used to produce reddening maps E ( B − V ) of the host galaxies. We compare the mean central value of E ( B − V ) to the X-ray obscuring column density along the line of sight to the AGN ( N _H ∼ 10 ^23−23.5 cm ^−2 ). We find that the extinction due to the host galaxy is present (0.6 ≲ E ( B − V ) ≲ 0.9; 1.9 ≲ A _V ≲ 2.8) and significantly contributes to the X-ray obscuration at a level of N _H ∼ 10 ^22.5 cm ^−2 assuming an SMC gas-to-dust ratio that amounts to ≲30% of the total obscuring column density. These early results, including three additional cases from CEERS, demonstrate the ability to resolve such dust structures with JWST and separate the different circumnuclear and galaxy-scale obscuring structures.

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 (

Lyman-break galaxy (LBG) candidates at z ≳ 10 are rapidly being identified in James Webb Space Telescope (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.6 σ SCUBA-2 detection at 850 μ m 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.

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.

We present the data release and data reduction process for the Epoch 1 NIRCam observations for the Cosmic Evolution Early Release Science Survey (CEERS). These data consist of NIRCam imaging in six broadband filters (F115W, F150W, F200W, F277W, F356W and F444W) and one medium-band filter (F410M) over four pointings, obtained in parallel with primary CEERS MIRI observations. We reduced the NIRCam imaging with the JWST Calibration Pipeline, with custom modifications and reduction steps designed to address additional features and challenges with the data. Here we provide a detailed description of each step in our reduction and a discussion of future expected improvements. Our reduction process includes corrections for known prelaunch issues such as 1/ f noise, as well as in-flight issues including snowballs, wisps, and astrometric alignment. Many of our custom reduction processes were first developed with prelaunch simulated NIRCam imaging over the full 10 CEERS NIRCam pointings. We present a description of the creation and reduction of this simulated data set in the Appendix. We provide mosaics of the real images in a public release, as well as our reduction scripts with detailed explanations to allow users to reproduce our final data products. These represent one of the first official public data sets released from the Directors Discretionary Early Release Science (DD-ERS) program.

We report on the host properties of five X-ray-luminous active galactic nuclei (AGN) identified at 3 < z < 5 in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey. Each galaxy has been imaged with the JWST Near-Infrared Camera, which provides rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host by fitting its spectral energy distribution using a combination of galaxy and AGN templates. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk, and one is dominated by pointlike emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. When compared to a sample of mass-matched inactive galaxies, we find that the AGN hosts have morphologies that are less disturbed and more bulge-dominated. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or poststarburst stellar population. The presence of AGN in passively evolving galaxies at z > 3 is significant because a rapid feedback mechanism is required in most semianalytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings show that AGN can continue to inject energy into these systems after their star formation is curtailed, potentially heating their halos and preventing renewed star formation. Additional observations will be needed to determine what role this feedback may play in helping to quench these systems and/or maintain their quiescent state.

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.

The selection of high-redshift galaxies often involves spectral energy distribution (SED) fitting to photometric data, an expectation for contamination levels, and measurement of sample completeness—all vetted through comparison to spectroscopic redshift measurements of a sub-sample. The first JWST data are now being taken over several extragalactic fields to different depths and across various areas, which will be ideal for the discovery and classification of galaxies out to distances previously uncharted. As spectroscopic redshift measurements for sources in this epoch will not be initially available to compare with the first photometric measurements of z > 8 galaxies, robust photometric redshifts are of the utmost importance. Galaxies at z > 8 are expected to have bluer rest-frame ultraviolet (UV) colors than typically used model SED templates, which could lead to catastrophic photometric redshift failures. We use a combination of BPASS and Cloudy models to create a supporting set of templates that match the predicted rest-UV colors of z > 8 simulated galaxies. We test these new templates by fitting simulated galaxies in a mock catalog, Yung et al., which mimic expected field depths and areas of the JWST Cosmic Evolution Early Release Science Survey ( m _5 _σ ∼ 28.6 over ∼100 arcmin ^2 ). We use EAZY to highlight the improvements in redshift recovery with the inclusion of our new template set and suggest criteria for selecting galaxies at 8 < z < 10 with the JWST, providing an important test case for observers venturing into this new era of astronomy.

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 the survey design, implementation, and outlook for COSMOS-Web, a 255 hr treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg ^2 NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5 σ point-source depths ranging ∼27.5–28.2 mag. In parallel, we will obtain 0.19 deg ^2 of MIRI imaging in one filter (F770W) reaching 5 σ point-source depths of ∼25.3–26.0 mag. COSMOS-Web will build on the rich heritage of multiwavelength observations and data products available in the COSMOS field. The design of COSMOS-Web is motivated by three primary science goals: (1) to discover thousands of galaxies in the Epoch of Reionization (6 ≲ z ≲ 11) and map reionization’s spatial distribution, environments, and drivers on scales sufficiently large to mitigate cosmic variance, (2) to identify hundreds of rare quiescent galaxies at z > 4 and place constraints on the formation of the universe’s most-massive galaxies ( M _⋆ > 10 ^10 M _⊙ ), and (3) directly measure the evolution of the stellar-mass-to-halo-mass relation using weak gravitational lensing out to z ∼ 2.5 and measure its variance with galaxies’ star formation histories and morphologies. In addition, we anticipate COSMOS-Web’s legacy value to reach far beyond these scientific goals, touching many other areas of astrophysics, such as the identification of the first direct collapse black hole candidates, ultracool subdwarf stars in the Galactic halo, and possibly the identification of z > 10 pair-instability supernovae. In this paper we provide an overview of the survey’s key measurements, specifications, goals, and prospects for new discovery.

Spectroscopic studies of extreme-ionization galaxies (EIGs) are critical to our understanding of exotic systems throughout cosmic time. These EIGs exhibit spectral features requiring >54.42 eV photons: the energy needed to ionize helium into He ^2+ fully and emit He ii recombination lines. Spectroscopic studies of EIGs can probe exotic stellar populations or accretion onto intermediate-mass black holes (∼10 ^2 –10 ^5 M _⊙ ), which are the possibly key contributors to the reionization of the Universe. To facilitate the use of EIGs as probes of high-ionization systems, we focus on ratios constructed from several rest-frame UV/optical emission lines: [O iii ] λ 5008, H β , [Ne iii ] λ 3870, [O ii ] λ λ 3727, 3729, and [Ne v ] λ 3427. These lines probe the relative intensity at energies of 35.12, 13.62, 40.96, 13.62, and 97.12 eV, respectively, covering a wider range of ionization than traced by other common rest-frame UV/optical techniques. We use the ratios of these lines ([Ne v ]/[Ne iii ] ≡ Ne53, [O iii ]/H β , and [Ne iii ]/[O ii ]), which are nearby in wavelength, mitigating the effects of dust attenuation and uncertainties in flux calibration. We make predictions from photoionization models constructed from Cloudy that use a broad range of stellar populations and black hole accretion models to explore the sensitivity of these line ratios to changes in the ionizing spectrum. We compare our models to observations from the Hubble Space Telescope and JWST of galaxies with strong high-ionization emission lines at z ∼ 0, z ∼ 2, and 5 < z < 8.5. We show that the Ne53 ratio can separate galaxies with ionization from “normal” stellar populations from those with active galactic nuclei and even “exotic” Population III models. We introduce new selection methods to identify galaxies with photoionization driven by Population III stars or intermediate-mass black hole accretion disks that could be identified in upcoming high-redshift spectroscopic surveys.

We present rest-frame optical emission-line flux ratio measurements for five z > 5 galaxies observed by the James Webb Space Telescope Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wavelength, despite the uncertain absolute spectrophotometry of the current version of the reductions. Compared to z ∼ 3 galaxies in the literature, the z > 5 galaxies have similar [O iii ] λ 5008/H β ratios, similar [O iii ] λ 4364/H γ ratios, and higher (∼0.5 dex) [Ne III ] λ 3870/[O II ] λ 3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [Ne III ] λ 3870/[O II ] λ 3728, [O iii ] λ 4364/H γ , and [O iii ] λ 5008/H β emission-line ratios are consistent with an interstellar medium (ISM) that has very high ionization ( $\mathrm{log}(Q)\simeq 8-9$ , units of cm s ^−1 ), low metallicity ( Z / Z _⊙ ≲ 0.2), and very high pressure ( $\mathrm{log}(P/k)\simeq 8-9$ , units of cm ^−3 ). The combination of [O iii ] λ 4364/H γ and [O iii ] λ (4960 + 5008)/H β line ratios indicate very high electron temperatures of $4.1\lt \mathrm{log}({T}_{e}/{\rm{K}})\lt 4.4$ , further implying metallicities of Z / Z _⊙ ≲ 0.2 with the application of low-redshift calibrations for “ T _e -based” metallicities. These observations represent a tantalizing new view of the physical conditions of the ISM in galaxies at cosmic dawn.

We present the results of a systematic search for candidate quiescent galaxies in the distant universe in 11 JWST fields with publicly available observations collected during the first 3 months of operations and covering an effective sky area of ∼145 arcmin ^2 . We homogeneously reduce the new JWST data and combine them with existing observations from the Hubble Space Telescope. We select a robust sample of ∼80 candidate quiescent and quenching galaxies at 3 < z < 5 using two methods: (1) based on their rest-frame UVJ colors, and (2) a novel quantitative approach based on Gaussian mixture modeling of the near-UV − U , U − V , and V − J rest-frame color space, which is more sensitive to recently quenched objects. We measure comoving number densities of massive ( M _⋆ ≥ 10 ^10.6 M _⊙ ) quiescent galaxies consistent with previous estimates relying on ground-based observations, after homogenizing the results in the literature with our mass and redshift intervals. However, we find significant field-to-field variations of the number densities up to a factor of 2–3, highlighting the effect of cosmic variance and suggesting the presence of overdensities of red quiescent galaxies at z > 3, as could be expected for highly clustered massive systems. Importantly, JWST enables the robust identification of quenching/quiescent galaxy candidates at lower masses and higher redshifts than before, challenging standard formation scenarios. All data products, including the literature compilation, are made publicly available.

Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of “Walking the Line”, a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodology. Emission lines across the spectrum from the millimeter to the UV were discussed, with most of the focus on the interstellar medium, but also some topics on the circumgalactic medium. The most important quality of a useful model is a good synergy with observations and experiments. Challenges in simulating line emission are identified, some of which are already being worked upon, and others that must be addressed in the future for models to agree with observations. Recent advances in several areas aiming at achieving that synergy are summarized here, from micro-physical to galactic and circum-galactic scale.