Search

We present the results of a James Webb Space Telescope NIRCam and NIRSpec investigation into the young massive star cluster (YMC) population of NGC 3256, the most cluster-rich luminous infrared galaxy in the Great Observatories All Sky LIRG Survey. We detect 3061 compact YMC candidates with a signal-to-noise ratio ≥3 at F150W, F200W, and F335M. Based on yggdrasil stellar population models, we identify 116/3061 sources with F150W – F200W > 0.47 and F200W – F355M > −1.37 colors, suggesting that they are young ( t ≤ 5 Myr), dusty ( A _V = 5−15), and massive ( M _⊙ > 10 ^5 ). This increases the sample of dust-enshrouded YMCs detected in this system by an order of magnitude relative to previous Hubble Space Telescope studies. With NIRSpec integral field unit pointings centered on the northern and southern nucleus, we extract the Pa α and 3.3 μ m polycyclic aromatic hydrocarbon (PAH) equivalent widths for eight bright and isolated YMCs. Variations in both the F200W – F335M color and 3.3 μ m PAH emission with the Pa α line strength suggest a rapid dust clearing ( 0) correspond to sources with E ( B − V ) > 3, which are typically missed in UV-optical studies. This underscores the importance of deep near-infrared imaging for finding and characterizing these very young and dust-embedded sources.

The unprecedented infrared spectroscopic capabilities of JWST have provided high-quality interstellar medium metallicity measurements and enabled characterization of the gas-phase mass–metallicity relation (MZR) for galaxies at z ≳ 5 for the first time. We analyze the gas-phase MZR and its evolution in a high-redshift suite of FIRE-2 cosmological zoom-in simulations at z = 5–12 and for stellar masses M _* ∼ 10 ^6 –10 ^10 M _⊙ . These simulations implement a multichannel stellar feedback model and produce broadly realistic galaxy properties, including when evolved to z = 0. The simulations predict very weak redshift evolution of the MZR over the redshift range studied, with the normalization of the MZR increasing by less than 0.01 dex as redshift decreases from z = 12 to z = 5. The median MZR in the simulations is well approximated as a constant power-law relation across this redshift range given by $\mathrm{log}(Z/{Z}_{\odot })=0.37\mathrm{log}({M}_{* }/{{\rm{M}}}_{\odot })-4.3$ . We find good agreement between our best-fit model and recent observations made by JWST at high redshift. The weak evolution of the MZR at z > 5 contrasts with the evolution at z ≲ 3, where increasing normalization of the MZR with decreasing redshift is observed and predicted by most models. The FIRE-2 simulations predict increasing scatter in the gas-phase MZR with decreasing stellar mass, in qualitative agreement with some observations.

The James Webb Space Telescope is unveiling a surprising lack of evolution in the number densities of ultraviolet (UV)-selected galaxies at redshift z ≳ 10. At the same time, observations and simulations are providing evidence for highly bursty star formation in high- z galaxies, resulting in significant scatter in their UV luminosities. Galaxies in low-mass dark matter halos are expected to experience most stochasticity due to their shallow potential wells. Here, we explore the impact of a mass-dependent stochasticity using a simple analytical model. We assume that scatter in the M _UV – M _h relation increases toward lower halo masses, following the decrease in halo escape velocity, ${\sigma }_{\mathrm{UV}}\sim {M}_{h}^{-1/3}$ , independent of redshift. Since low-mass halos are more dominant in the early universe, this model naturally predicts an increase in UV luminosity functions (LFs) at high redshifts compared to models without scatter. We make predictions for additional observables, which would be affected by stochasticity and could be used to constrain its amplitude, finding (i) galaxies are less clustered compared to the no-scatter scenario, with the difference increasing at higher- z ; (ii) assuming that star-bursting galaxies dominate the ionizing photon budget implies reionization starts earlier and is more gradual compared to the no-scatter case; (iii) at fixed UV magnitude, galaxies should exhibit wide ranges of UV slopes, nebular emission line strengths, and Balmer breaks. Comparing to observations, the mass-dependent stochasticity model successfully reproduces the observed LFs up to z ∼ 12. However, the model cannot match the observed z ∼ 14 LFs, implying additional physical processes enhance star formation efficiency in the earliest galaxies.

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.

We present observations and analysis of the starburst PACS-819 at z = 1.45 ( M _* = 10 ^10.7 M _⊙ ), using high-resolution (0.″1; 0.8 kpc) Atacama Large Millimeter/submillimeter Array (ALMA) and multiwavelength JWST images from the COSMOS-Web program. Dissimilar to Hubble Space Telescope (HST) ACS images in the rest-frame UV, the redder NIRCam and MIRI images reveal a smooth central mass concentration and spiral-like features, atypical for such an intense starburst. Through dynamical modeling of the CO ( J = 5–4) emission with ALMA, PACS-819 is rotation dominated and thus consistent with having a disk-like nature. However, kinematic anomalies in CO and asymmetric features in the bluer JWST bands (e.g., F150W) support a more disturbed nature likely due to interactions. The JWST imaging further enables us to map the distribution of stellar mass and dust attenuation, thus clarifying the relationships between different structural components not discernible in the previous HST images. The CO ( J = 5–4) and far-infrared dust continuum emission are cospatial with a heavily obscured starbursting core (

We present JWST and Atacama Large Millimeter/submillimeter Array (ALMA) imaging for the lensing system SPT0418−47, which includes a strongly lensed, dusty, star-forming galaxy at redshift z = 4.225 and an associated multiply imaged companion. The JWST NIRCam and MIRI imaging observations presented in this paper were acquired as part of the Early Release Science program Targeting Extremely Magnified Panchromatic Lensed Arcs and Their Extended Star formation (TEMPLATES). This data set provides robust mutiwavelength detections of stellar light in both the main (SPT0418A) and companion (SPT0418B) galaxies, while the ALMA detection of [C ii ] emission confirms that SPT0418B lies at the same redshift as SPT0418A. We infer that the projected physical separation of the two galaxies is 4.42 ± 0.05 kpc. We derive total magnifications of μ = 29 ± 1 and μ = 4.1 ± 0.7 for SPT0418A and SPT0418B, respectively. We use both prospector and cigale to derive stellar masses. We find that SPT0418A has a stellar mass of ${M}_{* }={3.4}_{-0.6}^{+1.1}\times {10}^{10}{M}_{\odot }$ from prospector or M _* = 1.5 ± 0.3 × 10 ^10 M _⊙ from cigale . The stellar mass ratio of SPT0418A and SPT0418B is roughly between 4 and 7 ( ${4.2}_{-1.6}^{+1.9}$ for prospector and 7.5 ± 3.7 for cigale ). We see evidence of extended structure associated with SPT0418A that is suggestive of a tidal feature. These features, along with the close projected proximity, imply that the system is interacting. Interestingly, the star formation rates and stellar masses of both galaxies are consistent with the main sequence of star-forming galaxies at this epoch, indicating that this ongoing interaction has not noticeably elevated the star formation levels.

The enormous increase in mid-IR sensitivity and spatial and spectral resolution provided by the JWST spectrographs enables, for the first time, detailed extragalactic studies of molecular vibrational bands. This opens an entirely new window for the study of the molecular interstellar medium in luminous infrared galaxies (LIRGs). We present a detailed analysis of rovibrational bands of gas-phase CO, H _2 O, C _2 H _2 , and HCN toward the heavily obscured eastern nucleus of the LIRG VV 114, as observed by NIRSpec and the medium resolution spectrograph on the Mid-InfraRed Instrument (MIRI MRS). Spectra extracted from apertures of 130 pc in radius show a clear dichotomy between the obscured active galactic nucleus (AGN) and two intense starburst regions. We detect the 2.3 μ m CO bandheads, characteristic of cool stellar atmospheres, in the star-forming regions, but not toward the AGN. Surprisingly, at 4.7 μ m, we find highly excited CO ( T _ex ≈ 700–800 K out to at least rotational level J = 27) toward the star-forming regions, but only cooler gas ( T _ex ≈ 200 K) toward the AGN. We conclude that only mid-infrared pumping through the rovibrational lines can account for the equilibrium conditions found for CO and H _2 O in the deeply embedded starbursts. Here, the CO bands probe regions with an intense local radiation field inside dusty young massive star clusters or near the most massive young stars. The lack of high-excitation molecular gas toward the AGN is attributed to geometric dilution of the intense radiation from the bright point source. An overview of the relevant excitation and radiative transfer physics is provided in an appendix.

We present new JWST NIRSpec integral field spectroscopy (IFS) data for the luminous infrared galaxy NGC 7469, a nearby (70.6 Mpc) active galaxy with a Seyfert 1.5 nucleus that drives a highly ionized gas outflow and a prominent nuclear star-forming ring. Using the superb sensitivity and high spatial resolution of the JWST instrument NIRSpec IFS, we investigate the role of the Seyfert nucleus in the excitation and dynamics of the circumnuclear gas. Our analysis focuses on the [Fe ii ], H _2 , and hydrogen recombination lines that trace the radiation/shocked-excited molecular and ionized interstellar medium around the active galactic nucleus (AGN). We investigate gas excitation through H _2 /Br γ and [Fe ii ]/Pa β emission line ratios and find that photoionization by the AGN dominates within the central 300 pc of the galaxy except in a small region that shows signatures of shock-heated gas; these shock-heated regions are likely associated with a compact radio jet. In addition, the velocity field and velocity dispersion maps reveal complex gas kinematics. Rotation is the dominant feature, but we also identify noncircular motions consistent with gas inflows as traced by the velocity residuals and the spiral pattern in the Pa α velocity dispersion map. The inflow is 2 orders of magnitude higher than the AGN accretion rate. The compact nuclear radio jet has enough power to drive the highly ionized outflow. This scenario suggests that the inflow and outflow are in a self-regulating feeding–feedback process, with a contribution from the radio jet helping to drive the outflow.

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 have observed the z = 4.3 protocluster SPT2349−56 with the Australia Telescope Compact Array (ATCA) with the aim of detecting radio-loud active galactic nuclei (AGNs) among the ∼30 submillimeter (submm) galaxies (SMGs) identified in the structure. We detect the central complex of submm sources at 2.2 GHz with a luminosity of L _2.2 = (4.42 ± 0.56) × 10 ^25 W Hz ^−1 . MeerKAT and the Australian Square Kilometre Array Pathfinder also detect the source at 816 MHz and 888 MHz, respectively, constraining the radio spectral index to α = −1.45 ± 0.16, implying L _1.4,rest = (2.2 ± 0.2) × 10 ^26 W Hz ^−1 . The radio observations do not have sufficient spatial resolution to uniquely identify one of the three Atacama Large Millimeter/submillimeter Array (ALMA) galaxies as the AGN, however the ALMA source properties themselves suggest a likely host. This radio luminosity is ∼100× higher than expected from star formation, assuming the usual far-infrared–radio correlation, indicating an AGN driven by a forming brightest cluster galaxy. None of the SMGs in SPT2349−56 show signs of AGNs in any other diagnostics available to us, highlighting the radio continuum as a powerful probe of obscured AGNs. We compare these results to field samples of radio sources and SMGs, along with the 22 gravitationally lensed SPT-SMGs also observed in the ATCA program, as well as powerful radio galaxies at high redshifts. The (3.3 ± 0.7) × 10 ^38 W of power from the radio-loud AGN sustained over 100 Myr is comparable to the binding energy of the gas mass of the central halo, and similar to the instantaneous energy injection from supernova feedback from the SMGs in the core region. The SPT2349−56 radio-loud AGNs may be providing strong feedback on a nascent intracluster medium.

Dusty star-forming galaxies (DSFGs) are among the most massive and active star-forming galaxies during the cosmic noon. Theoretical studies have proposed various formation mechanisms of DSFGs, including major merger-driven starbursts and secular star-forming disks. Here, we report J0107a, a bright (∼8 mJy at observed-frame 888 μ m) DSFG at z = 2.467 that appears to be a gas-rich massive disk and might be an extreme case of the secular disk scenario. J0107a has a stellar mass M _⋆ ∼ 5 × 10 ^11 M _⊙ , molecular gas mass M _mol ≳ 10 ^11 M _⊙ , and a star formation rate of ∼500 M _⊙ yr ^−1 . J0107a does not have a gas-rich companion. The rest-frame 1.28 μ m JWST NIRCam image of J0107a shows a grand-design spiral with a prominent stellar bar extending ∼15 kpc. The Atacama Large Millimeter/submillimeter Array Band 7 continuum map reveals that the dust emission originates from both the central starburst and the stellar bar. 3D disk modeling of the CO(4–3) emission line indicates a dynamically cold disk with rotation-to-dispersion ratio ${V}_{\max }/\sigma \sim 8$ . The results suggest a bright DSFG may have a nonmerger origin, and its vigorous star formation may be triggered by the bar and/or rapid gas inflow.

We present JWST Mid-InfraRed Instrument/Medium-Resolution Spectrometer (MIRI/MRS) observations of an infrared luminous disk galaxy, FLS1, at z ∼ 0.54. With a lookback time of 5 Gyr, FLS1 is chronologically at the midpoint between the peak epoch of star formation and the present day. The MRS data provide maps of the atomic fine structure lines [Ar ii ]6.99, [Ar iii ]8.99, [Ne ii ]12.81, and [Ne iii ]15.55 μ m, polycyclic aromatic hydrocarbon (PAH) features at 3.3, 6.2, and 11.3 μ m, and the warm molecular gas indicators H _2 S(5) and H _2 S(3); all these emission features are spatially resolved. We find that the PAH emission is more extended along the northern side of the galaxy when compared to the well-studied star formation tracer [Ne ii ]. The H _2 rotational lines, which are shock indicators, are strongest and most extended on the southern side of the galaxy. [Ar ii ] is the second brightest fine structure line detected in FLS1 and we show that it is a useful kinematic probe that can be detected with JWST out to z ∼ 3. Velocity maps of [Ar ii ] show a rotating disk with signs of turbulence. Our results provide an example of how spatially resolved mid-infrared spectroscopy can allow us to better understand the star formation and interstellar medium conditions in a galaxy halfway back to the peak epoch of galaxy evolution.

We present James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec) integral field spectroscopy of the nearby luminous infrared galaxy NGC 7469. We take advantage of the high spatial/spectral resolution and wavelength coverage of JWST/NIRSpec to study the 3.3 μ m neutral polycyclic aromatic hydrocarbon (PAH) grain emission on ∼200 pc scales. A clear change in the average grain properties between the star-forming ring and the central AGN is found. Regions in the vicinity of the AGN, with [Ne iii ]/[Ne ii ] > 0.25, tend to have larger grain sizes and lower aliphatic-to-aromatic (3.4/3.3) ratios, indicating that smaller grains are preferentially removed by photodestruction in the vicinity of the AGN. PAH emission at the nucleus is weak and shows a low 11.3/3.3 PAH ratio. We find an overall suppression of the total PAH emission relative to the ionized gas in the central 1 kpc region of the AGN in NGC 7469 compared to what has been observed with Spitzer on 3 kpc scales. However, the fractional 3.3 μ m–to–total PAH power is enhanced in the starburst ring, possibly due to a variety of physical effects on subkiloparsec scales, including recurrent fluorescence of small grains or multiple photon absorption by large grains. Finally, the IFU data show that while the 3.3 μ m PAH-derived star formation rate (SFR) in the ring is 27% higher than that inferred from the [Ne ii ] and [Ne iii ] emission lines, the integrated SFR derived from the 3.3 μ m feature would be underestimated by a factor of 2 due to the deficit of PAHs around the AGN, as might occur if a composite system like NGC 7469 were to be observed at high redshift.

Recent discoveries of a significant population of bright galaxies at cosmic dawn $\left(z\gtrsim 10\right)$ have enabled critical tests of cosmological galaxy formation models. In particular, the bright end of the galaxys’ UV luminosity functions (UVLFs) appear higher than predicted by many models. Using approximately 25,000 galaxy snapshots at 8 ≤ z ≤ 12 in a suite of FIRE-2 cosmological “zoom-in” simulations from the Feedback in Realistic Environments (FIRE) project, we show that the observed abundance of UV-bright galaxies at cosmic dawn is reproduced in these simulations with a multichannel implementation of standard stellar feedback processes, without any fine-tuning. Notably, we find no need to invoke previously suggested modifications, such as a nonstandard cosmology, a top-heavy stellar initial mass function, or a strongly enhanced star formation efficiency. We contrast the UVLFs predicted by bursty star formation in these original simulations to those derived from star formation histories (SFHs) smoothed over prescribed timescales (e.g., 100 Myr). The comparison demonstrates that the strongly time-variable SFHs predicted by the FIRE simulations play a key role in correctly reproducing the observed, bright-end UVLFs at cosmic dawn: the bursty SFHs induce order-or-magnitude changes in the abundance of UV-bright ( M _UV ≲ −20) galaxies at z ≳ 10. The predicted bright-end UVLFs are consistent with both the spectroscopically confirmed population and the photometrically selected candidates. We also find good agreement between the predicted and observationally inferred integrated UV luminosity densities, which evolve more weakly with redshift in FIRE than suggested by some other models.

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 James Webb Space Telescope (JWST) imaging of NGC 7469 with the Near-Infrared Camera and the Mid-InfraRed Instrument. NGC 7469 is a nearby, z = 0.01627, luminous infrared galaxy that hosts both a Seyfert Type-1.5 nucleus and a circumnuclear starburst ring with a radius of ∼0.5 kpc. The new near-infrared (NIR) JWST imaging reveals 66 star-forming regions, 37 of which were not detected by Hubble Space Telescope (HST) observations. Twenty-eight of the 37 sources have very red NIR colors that indicate obscurations up to A _v ∼ 7 and a contribution of at least 25% from hot dust emission to the 4.4 μ m band. Their NIR colors are also consistent with young (

We present the results of a James Webb Space Telescope NIRCam investigation into the young massive star cluster (YMC) population in the luminous infrared galaxy VV 114. We identify 374 compact YMC candidates with signal-to-noise ratios ≥ 3, 5, and 5 at F150W, F200W, and F356W, respectively. A direct comparison with our HST cluster catalog reveals that ∼20% of these sources are undetected at optical wavelengths. Based on yggdrasil stellar population models, we identify 17 YMC candidates in our JWST imaging alone with F150W – F200W and F200W – F356W colors suggesting they are all very young, dusty ( A _V = 5–15), and massive (10 ^5.8 < M _⊙ < 10 ^6.1 ). The discovery of these “hidden” sources, many of which are found in the “overlap” region between the two nuclei, quadruples the number of t < 3 Myr clusters and nearly doubles the number of t < 6 Myr clusters detected in VV 114. Now extending the cluster age distribution ( ${dN}/d\tau \propto {\tau }^{\gamma }$ ) to the youngest ages, we find a slope of γ = −1.30 ± 0.39 for 10 ^6 < τ (yr) < 10 ^7 , which is consistent with the previously determined value from 10 ^7 < τ (yr) < 10 ^8.5 , and confirms that VV 114 has a steep age distribution slope for all massive star clusters across the entire range of cluster ages observed. Finally, the consistency between our JWST- and HST-derived age distribution slopes indicates that the balance between cluster formation and destruction has not been significantly altered in VV 114 over the last 0.5 Gyr.

James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) images of the luminous infrared (IR) galaxy VV 114 are presented. This redshift ∼0.020 merger has a western component (VV 114W) rich in optical star clusters and an eastern component (VV 114E) hosting a luminous mid-IR nucleus hidden at UV and optical wavelengths by dust lanes. With MIRI, the VV 114E nucleus resolves primarily into bright NE and SW cores separated by 630 pc. This nucleus comprises 45% of the 15 μ m light of VV 114, with the NE and SW cores having IR luminosities, L _IR (8 − 1000 μ m) ∼ 8 ± 0.8 × 10 ^10 L _⊙ and ∼ 5 ± 0.5 × 10 ^10 L _⊙ , respectively, and IR densities, Σ _IR ≳ 2 ± 0.2 × 10 ^13 L _⊙ kpc ^−2 and ≳ 7 ± 0.7 × 10 ^12 L _⊙ kpc ^−2 , respectively—in the range of Σ _IR for the Orion star-forming core and the nuclei of Arp 220. The NE core, previously speculated to have an active galactic nucleus (AGN), has starburst-like mid-IR colors. In contrast, the VV 114E SW core has AGN-like colors. Approximately 40 star-forming knots with L _IR ∼ 0.02–5 × 10 ^10 L _⊙ are identified, 28% of which have no optical counterpart. Finally, diffuse emission accounts for 40%–60% of the mid-IR emission. Mostly notably, filamentary polycyclic aromatic hydrocarbon (PAH) emission stochastically excited by UV and optical photons accounts for half of the 7.7 μ m light of VV 114. This study illustrates the ability of JWST to detect obscured compact activity and distributed PAH emission in the most extreme starburst galaxies in the local universe.

We have used the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST) to obtain the first spatially resolved, mid-infrared images of IIZw096, a merging luminous infrared galaxy (LIRG) at z = 0.036. Previous observations with the Spitzer Space Telescope suggested that the vast majority of the total IR luminosity ( L _IR ) of the system originated from a small region outside of the two merging nuclei. New observations with JWST/MIRI now allow an accurate measurement of the location and luminosity density of the source that is responsible for the bulk of the IR emission. We estimate that 40%–70% of the IR bolometric luminosity, or 3–5 × 10 ^11 L _⊙ , arises from a source no larger than 175 pc in radius, suggesting a luminosity density of at least 3–5 × 10 ^12 L _⊙ kpc ^−2 . In addition, we detect 11 other star-forming sources, five of which were previously unknown. The MIRI F1500W/F560W colors of most of these sources, including the source responsible for the bulk of the far-IR emission, are much redder than the nuclei of local LIRGs. These observations reveal the power of JWST to disentangle the complex regions at the hearts of merging, dusty galaxies.

The nearby, luminous infrared galaxy NGC 7469 hosts a Seyfert nucleus with a circumnuclear star-forming ring and is thus the ideal local laboratory for investigating the starburst–AGN (active galactic nucleus) connection in detail. We present integral-field observations of the central 1.3 kpc region in NGC 7469 obtained with the JWST Mid-InfraRed Instrument. Molecular and ionized gas distributions and kinematics at a resolution of ∼100 pc over the 4.9–7.6 μ m region are examined to study the gas dynamics influenced by the central AGN. The low-ionization [Fe ii ] λ 5.34 μ m and [Ar ii ] λ 6.99 μ m lines are bright on the nucleus and in the starburst ring, as opposed to H _2 S(5) λ 6.91 μ m, which is strongly peaked at the center and surrounding ISM. The high-ionization [Mg v ] line is resolved and shows a broad, blueshifted component associated with the outflow. It has a nearly face-on geometry that is strongly peaked on the nucleus, where it reaches a maximum velocity of −650 km s ^−1 , and extends about 400 pc to the east. Regions of enhanced velocity dispersion in H _2 and [Fe ii ] ∼ 180 pc from the AGN that also show high L (H _2 )/ L (PAH) and L ([Fe ii ])/ L (Pf α ) ratios to the W and N of the nucleus pinpoint regions where the ionized outflow is depositing energy, via shocks, into the dense interstellar medium between the nucleus and the starburst ring. These resolved mid-infrared observations of the nuclear gas dynamics demonstrate the power of JWST and its high-sensitivity integral-field spectroscopic capability to resolve feedback processes around supermassive black holes in the dusty cores of nearby luminous infrared galaxies.

We present James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) integral-field spectroscopy of the nearby merging, luminous infrared galaxy, NGC 7469. This galaxy hosts a Seyfert type-1.5 nucleus, a highly ionized outflow, and a bright, circumnuclear star-forming ring, making it an ideal target to study active galactic nucleus (AGN) feedback in the local universe. We take advantage of the high spatial/spectral resolution of JWST/MIRI to isolate the star-forming regions surrounding the central active nucleus and study the properties of the dust and warm molecular gas on ∼100 pc scales. The starburst ring exhibits prominent polycyclic aromatic hydrocarbon (PAH) emission, with grain sizes and ionization states varying by only ∼30%, and a total star formation rate of 10–30 M _⊙ yr ^−1 derived from fine structure and recombination emission lines. Using pure rotational lines of H _2 we detect 1.2 × 10 ^7 M _⊙ of warm molecular gas at a temperature higher than 200 K in the ring. All PAH bands get significantly weaker toward the central source, where larger and possibly more ionized grains dominate the emission, likely the result of the ionizing radiation and/or the fast wind emerging from the AGN. The small grains and warm molecular gas in the bright regions of the ring however display properties consistent with normal star-forming regions. These observations highlight the power of JWST to probe the inner regions of dusty, rapidly evolving galaxies for signatures of feedback and inform models that seek to explain the coevolution of supermassive black holes and their hosts.

We present the first spatially resolved maps of gas-phase metallicity for two dust-obscured star-forming galaxies at z ∼ 4, from the JWST TEMPLATES Early Release Science program, derived from NIRSpec integral field unit spectroscopy of the H α and [N ii ] emission lines. Empirical optical line calibrations are used to determine that the sources are globally enriched to near-solar levels. While one source shows elevated [N ii ]/H α ratios and broad H α emission consistent with the presence of an active galactic nucleus in a ≳1 kpc region, we argue that both systems have already undergone significant metal enrichment as a result of their extremely high star formation rates. Utilizing Atacama Large Millimeter/submillimeter Array rest-frame 380 μ m continuum and [C i ]( ^3 P _2 – ^3 P _1 ) line maps we compare the spatial variation of the metallicity and gas-to-dust ratio in the two galaxies, finding the two properties to be anticorrelated on highly resolved spatial scales, consistent with various literature studies of z ∼ 0 galaxies. The data are indicative of the enormous potential of JWST to probe the enrichment of the interstellar medium on ∼kpc scales in extremely dust-obscured systems at z ∼ 4 and beyond.

A growing number of far-infrared (FIR) bright sources completely invisible in deep extragalactic optical surveys hint at an elusive population of z > 4 dusty, star-forming galaxies. Cycle 1 JWST surveys are now detecting their rest-frame optical light, which provides key insight into their stellar properties and statistical constraints on the population as a whole. This work presents the JWST Near Infrared Camera (NIRCam) counterpart from the COSMOS-Web survey to an FIR SCUBA-2 and Atacama Large Millimeter/submillimeter Array (ALMA) source, AzTECC71, which was previously undetected at wavelengths shorter than 850 μ m. AzTECC71, among the reddest galaxies in COSMOS-Web with F277W − F444W ∼ 0.9, is undetected in NIRCam/F150W and F115W and fainter in F444W than other submillimeter galaxies identified in COSMOS-Web by 2–4 magnitudes. This is consistent with the system having both a lower stellar mass and higher redshift than the median dusty, star-forming galaxy. With deep ground- and space-based upper limits combined with detections in F277W, F444W, and the FIR including ALMA Band 6, we find a high probability (99%) that AzTECC71 is at z > 4 with ${z}_{\mathrm{phot}}={5.7}_{-0.7}^{+0.8}$ . This galaxy is massive ( $\mathrm{log}\,{M}_{* }/{M}_{\odot }\sim 10.7$ ) and infrared-luminous ( $\mathrm{log}\,{L}_{\mathrm{IR}}/{L}_{\odot }\sim 12.7$ ), comparable to other optically undetected but FIR-bright dusty, star-forming galaxies at z > 4. This population of luminous, infrared galaxies at z > 4 is largely unconstrained but comprises an important bridge between the most extreme dust-obscured galaxies and more typical high-redshift star-forming galaxies. If further FIR-selected galaxies that drop out of the F150W filter in COSMOS-Web have redshifts z > 4 like AzTECC71, then the volume density of such sources may be ∼3–10 × greater than previously estimated.

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.

The surface densities of gas, dust, and stars provide a window into the physics of star formation that, until the advent of high-resolution far-IR/submillimeter observations, has been historically difficult to assess among dusty galaxies. To study the link between IR surface densities and dust properties, we leverage the Atacama Large Millimetre/Submillimetre Array archive to measure the extent of cold dust emission in 15 z ∼ 2 IR-selected galaxies selected on the basis of having available mid-IR spectroscopy from Spitzer. We use the mid-IR spectra to constrain the relative balance between dust heating from star formation and active galactic nuclei (AGNs), and to measure emission from polycylic aromatic hydrocarbons (PAHs), small dust grains that play a key role in the photoelectric heating of gas. In general, we find that dust-obscured star formation at high IR surface densities exhibits similar properties at low and high redshift, namely, local luminous IR galaxies (LIRGs) have comparable PAH luminosity to total dust mass ratios as high- z galaxies, and star formation at z ∼ 0–2 is more efficient at high IR surface densities despite the fact that our sample of high- z galaxies is closer to the main sequence than local LIRGs. High star formation efficiencies are coincident with a decline in the PAH-to-IR luminosity ratio reminiscent of the deficit observed in far-IR fine-structure lines. Changes in the gas and dust conditions arising from high star formation surface densities might help drive the star formation efficiency up. This could help explain the high efficiencies needed to reconcile star formation and gas volume densities in dusty galaxies at cosmic noon.

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.

We employ the Feedback In Realistic Environments (FIRE-2) physics model to study how the properties of giant molecular clouds (GMCs) evolve during galaxy mergers. We conduct a pixel-by-pixel analysis of molecular gas properties in both the simulated control galaxies and galaxy major mergers. The simulated GMC pixels in the control galaxies follow a similar trend in a diagram of velocity dispersion ( σ _v ) versus gas surface density (Σ _mol ) to the one observed in local spiral galaxies in the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) survey. For GMC pixels in simulated mergers, we see a significant increase of a factor of 5–10 in both Σ _mol and σ _v , which puts these pixels above the trend of PHANGS galaxies in the σ _v versus Σ _mol diagram. This deviation may indicate that GMCs in the simulated mergers are much less gravitationally bound compared with simulated control galaxies with virial parameters ( α _vir ) reaching 10–100. Furthermore, we find that the increase in α _vir happens at the same time as the increase in global star formation rate, which suggests that stellar feedback is responsible for dispersing the gas. We also find that the gas depletion time is significantly lower for high- α _vir GMCs during a starburst event. This is in contrast to the simple physical picture that low- α _vir GMCs are easier to collapse and form stars on shorter depletion times. This might suggest that some other physical mechanisms besides self-gravity are helping the GMCs in starbursting mergers collapse and form stars.

SPT0311-58, a system of two interacting galaxies in the Epoch of Reionization, exists in one of the rarest, most massive dark matter halos theoretically possible in that era. Studying the interstellar medium (ISM) in these galaxies can illuminate the process of galaxy formation in the early Universe. In this work, we explore the multiphase ISM in this system, using ALMA observations of the [C ii ] 158, [O i ] 146, [N ii ] 122, and [O iii ] 88 fine-structure lines and dust continuum. We find wide variations in line ratios between the eastern and western galaxies, as well as across the western galaxy. Continuum colors indicate that SPT0311-58 E has a higher ionization parameter ( $\mathrm{log}U\approx -2.8$ ) than SPT0311-58 W ( $\mathrm{log}U\approx -3.1$ ). The ratio of [O iii ] 88–[N ii ] 122 and the ionization parameter constraints combine to demonstrate near-solar metallicity in these objects just 800 Myr after the Big Bang.

The driving of turbulence in galaxies is deeply connected with the physics of feedback, star formation, outflows, accretion, and radial transport in disks. The velocity dispersion of gas in galaxies therefore offers a promising observational window into these processes. However, the relative importance of each of these mechanisms remains controversial. In this work we revisit the possibility that turbulence on galactic scales is driven by the direct impact of accreting gaseous material on the disk. We measure this effect in a disk-like star-forming galaxy in IllustrisTNG, using the high-resolution cosmological magnetohydrodynamical simulation TNG50. We employ Lagrangian tracer particles with a high time cadence of only a few million years to identify accretion and other events. The energies of particles are measured by stacking the events in bins of time around the event. The average effect of each event is measured by fitting explicit models for the kinetic and turbulent energies as a function of time. These measurements are corroborated by cross-correlating the turbulent energy with other time series and searching for signals of causality, i.e., asymmetries across zero time lag. We find that accretion contributes to the large-scale turbulent kinetic energy even if it does not dominate in this ∼5 × 10 ^9 M _⊙ stellar mass galaxy. Extrapolating this finding to a range of galaxy masses, we find that there are regimes where energy from direct accretion may dominate the turbulent energy budget, particularly in disk outskirts, galaxies less massive than the Milky Way, and at redshift ∼2.

We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation (available at http://flathub.flatironinstitute.org/fire ) from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multiphase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Type Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from three suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way (MW)–mass galaxies, five SMC/LMC-mass galaxies, and four lower-mass galaxies including one ultrafaint; we release 39 snapshots across z = 0–10. The second comprises four massive galaxies, with 19 snapshots across z = 1–10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5–10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all MW-mass galaxies, we release the formation coordinates and an “ex situ” flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations.

We present the characteristics of 2 mm selected sources from the largest Atacama Large Millimeter/submillimeter Array (ALMA) blank-field contiguous survey conducted to date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184 arcmin ^2 at 2 mm. Twelve of 13 detections above 5 σ are attributed to emission from galaxies, 11 of which are dominated by cold dust emission. These sources have a median redshift of $\langle {z}_{2\,\mathrm{mm}}\rangle ={3.6}_{-0.3}^{+0.4}$ primarily based on optical/near-infrared photometric redshifts with some spectroscopic redshifts, with 77% ± 11% of sources at z > 3 and 38% ± 12% of sources at z > 4. This implies that 2 mm selection is an efficient method for identifying the highest-redshift dusty star-forming galaxies (DSFGs). Lower-redshift DSFGs ( z < 3) are far more numerous than those at z > 3 yet are likely to drop out at 2 mm. MORA shows that DSFGs with star formation rates in excess of 300 M _⊙ yr ^−1 and a relative rarity of ∼10 ^−5 Mpc ^−3 contribute ∼30% to the integrated star formation rate density at 3 < z < 6. The volume density of 2 mm selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at z > 2. Analysis of MORA sources’ spectral energy distributions hint at steeper empirically measured dust emissivity indices than reported in typical literature studies, with $\langle \beta \rangle ={2.2}_{-0.4}^{+0.5}$ . The MORA survey represents an important step in taking census of obscured star formation in the universe’s first few billion years, but larger area 2 mm surveys are needed to more fully characterize this rare population and push to the detection of the universe’s first dusty galaxies.