This paper presents a deep MIRI/JWST medium-resolution spectroscopy (MRS) covering the rest-frame optical spectrum of the GN-z11 galaxy. The [O III] 5008 Å and Hα emission lines are detected and spectroscopically resolved. The line profiles are well modeled by a narrow Gaussian component with intrinsic full widths at half maximum of 189 ± 25 and 231 ± 52 km s-1, respectively. We do not find any evidence of a dominant broad Hα emission line component tracing a broad-line region in a type 1 active galactic nucleus (AGN). The existence of an accreting black hole dominating the optical continuum and emission lines of GN-z11 is not compatible with the measured Hα and [O III] 5008 Å luminosities. If the well-established relations for low-z AGNs apply in GN-z11, the [O III] 5008 Å and Hα luminosities would imply extremely high super-Eddington ratios (λE > 290), and bolometric luminosities ∼20 times those derived from the UV/optical continuum. However, a broad (∼430-470 km s-1) and weak (< 20-30%) Hα line component, tracing a minor AGN contribution in the optical, cannot be completely ruled out with the sensitivity of the current data. The physical and excitation properties of the ionized gas are consistent with a low-metallicity starburst with a star formation rate of 24 ± 3 M⊙ yr-1. The electron temperature of the ionized gas is Te (O++) = 14 000 ± 2100 K, while the direct-Te gas-phase metallicity is 12 + log(O/H) = 7.91 ± 0.07 (Z = 0.17 ± 0.03 Z⊙). The optical line ratios locate GN-z11 in the starburst or AGN region, but they are more consistent with those of local low-metallicity starbursts and high-z luminous galaxies detected at redshifts similar to GN-z11. We conclude that the MRS optical spectrum of GN-z11 is consistent with that of a massive, compact, and low-metallicity starburst galaxy. Its high star formation and stellar mass surface densities are close to those of the densest stellar clusters, and we therefore speculate that GN-z11 might undergo a feedback-free, highly efficient starburst phase. Additional JWST data are needed to validate this scenario and other recently proposed alternatives to explain the existence of bright compact galaxies in the early Universe.

Context. The recently launched James Webb Space Telescope (JWST) is opening new observing windows on the distant Universe. Among JWST’s instruments, the Mid Infrared Instrument (MIRI) offers the unique capability of imaging observations at wavelengths of λ > 5 μm. This enables unique access to the rest frame near-infrared (NIR, λ ≥ 1 μm) emission from galaxies at redshifts of z > 4 and the visual (λ ≳ 5000 Å) rest frame for z > 9. We report here on the guaranteed time observations (GTO), from the MIRI European Consortium, of the Hubble Ultra Deep Field (HUDF), forming the MIRI Deep Imaging Survey (MIDIS), consisting of an on source integration time of ∼41 hours in the MIRI/F560W (5.6 μm) filter. The F560W filter was selected since it would produce the deepest data in terms of AB magnitudes in a given time. To our knowledge, this constitutes the longest single filter exposure obtained with JWST of an extragalactic field as of yet.

Aims. The HUDF is one of the most observed extragalactic fields, with extensive multi-wavelength coverage, where (before JWST) galaxies up to z ∼ 7 have been confirmed, and at z > 10 suggested, from HST photometry. We aim to characterise the galaxy population in HUDF at 5.6 μm, enabling studies such as: the rest frame NIR morphologies for galaxies at z ≲ 4.6, probing mature stellar populations and emission lines in z > 6 sources, intrinsically red and dusty galaxies, and active galactic nuclei (AGNs) and their host galaxies at intermediate redshifts.

Methods. We reduced the MIRI data using the official JWST pipeline, augmented by in-house custom scripts. We measured the noise characteristics of the resulting image. Galaxy photometry was obtained, and photometric redshifts were estimated for sources with available multi-wavelength photometry (and compared to spectroscopic redshifts when available).

Results. Over the deepest part of our image, the 5σ point source limit is 28.65 mag AB (12.6 nJy), ∼0.35 mag better than predicted by the JWST exposure time calculator. We find ∼2500 sources, the overwhelming majority of which are distant galaxies, but we note that spurious sources likely remain at faint magnitudes due to imperfect cosmic ray rejection in the JWST pipeline. More than 500 galaxies with available spectroscopic redshifts, up to z ≈ 11, have been identified, the majority of which are at z < 6. More than 1000 galaxies have reliable photometric redshift estimates, of which ∼25 are at 6 < z < 12. The point spread function in the F560W filter has a full width at half maximum (FWHM) of ≈0.2″ (corresponding to 1.4 kpc at z = 4), allowing the NIR rest frame morphologies and stellar mass distributions to be resolved for z < 4.5. Moreover, > 100 objects with very red NIRCam vs MIRI (3.6–5.6 μm > 1 mag) colours have been found, suggestive of dusty or old stellar populations at high redshifts.

Conclusions. We conclude that MIDIS surpasses preflight expectations and that deep MIRI imaging has great potential to characterise the galaxy population from cosmic noon to dawn.

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

Gas kinematics affect the radiative transfer and escape of hydrogen Lyman-α (Lyα) emission from galaxies. We investigate this interplay empirically by relating the ionised gas kinematics of 42 galaxies in the extended Lyα Reference Sample (eLARS) with their Lyα escape fractions, fescLyα, Lyα equivalent widths, EWLyα, and Lyα luminosities, LLyα. To this aim we use PMAS integral-field spectroscopic observations of the Balmer-α line. Our sample contains 18 rotating discs, 13 perturbed rotators, and 13 galaxies with more complex kinematics. The distributions of fescLyα, EWLyα, and LLyα do not differ significantly between these kinematical classes, but the largest Lyα observables are found amongst the kinematically complex systems. We find no trends between either fescLyα or EWLyα and kinematic or photometric inclinations. We calculate shearing velocities, υshear, and intrinsic velocity dispersions, σobs0 (empirically corrected for beam-smearing effects), as global kinematical measures for each galaxy. The sample is characterised by highly turbulent motions (30 km s−1 ≲ σobs0 ≲ 80 km s−1) and more than half of the sources show dispersion-dominated kinematics. We uncover clear trends between Lyα observables and global kinematical statistics: EWLyα and LLyα correlate with σobs0 , while fescLyα anti-correlates with υshear and υshear/σobs0 . Moreover, we find, that galaxies with EWLyα ≥ 20 Å are characterised by higher σ0 and lower υshear/σobs0 than galaxies below this threshold. We discuss the statistical importance of υshear, σobs0 , and υshear/σobs0 for regulating the Lyα observables in comparison to other galaxy parameters. It emerges that σobs0 is the dominating parameter for regulating EWLyα and that is as important as nebular extinction, gas covering fraction, and ionising photon production efficiency in regulating fescLyα. A simple scenario where the starburst age is simultaneously regulating turbulence, EWLyα, and fescesc is not supported by our observations. However, we show that the small-scale distribution of dust appears to be influenced by turbulence in some galaxies. In support of our observational result, we discuss how turbulence is theoretically expected to play a significant role in modulating fescLyα

Context. The Lyman-α (Lyα) line of hydrogen is a well-known tracer of galaxies at high redshift. However, the connection between Lyα observables and galaxy properties has not been fully established, limiting the use of the line to probe the physics of galaxies. Aims. Here, we derive the global neutral hydrogen gas (Hi) properties of nearby Lyα-emitting galaxies to assess the impact of neutral gas on the Lyα output of galaxies. Methods. We observed the 21 cm line emission using the Karl G. Jansky Very Large Array in D-array configuration (∼5500 resolution, ∼38 kpc) for 37 star-forming galaxies with available Hubble Space Telescope (HST) Lyα imaging from the Lyman Alpha Reference Samples. Results. We detected 21 cm emission for 33 out of the 37 galaxies observed. We found no significant correlation of global Hi properties (including Hi mass, column density, gas fraction, depletion time, line width, or velocity shift between Hi and Lyα), with the Lyα luminosity, escape fractions, or equivalent widths (EW) derived with HST photometry. Additionally, both Lyα-emitters and weak or non-emitters are distributed evenly along the Hi parameter space of optically selected z = 0 galaxies. Around 74% of the sample is undergoing galaxy interaction, this fraction is higher for Lyα-emitters (83% for galaxies with EW ≥ 20 Å) than for weak or non-emitters (70%). Nevertheless, galaxies identified as interacting have Lyα and Hi properties statistically consistent with those of non-interacting galaxies. Conclusions. Our results show that global Hi properties (on scales >30kpc) have little direct impact on the Lyα output from galaxies. Instead, Hi likely regulates Lyα emission on small scales: statistical comparisons of Lyα and high angular resolution 21 cm observations are required to fully assess the role of Hi in Lyα radiative transfer. While our study indicates that major and minor galaxy mergers could play a role in the emission of Lyα photons in the Local Universe, especially for galaxies with high Hi fractions, the line of sight that a system is observed through ultimately determines the Lyα observables.

We report the discovery of Cerberus, an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at signal-to-noise ratio (S/N) ∼ 6, with F1000W ∼ 27 mag, and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey (JADES), fainter than the 30.0-30.5 mag 5σ detection limits in individual bands, as well as in the MIDIS F560W ultradeep data (∼29 mag, 5σ). Analyzing the spectral energy distribution built with low-S/N (<5) measurements in individual optical-to-mid-infrared filters and higher-S/N (≳5) measurements in stacked NIRCam data, we discuss the possible nature of this red NIRCam-dark source using a battery of codes. We discard the possibility of Cerberus being a solar system body based on the <0.″016 proper motion in the 1 yr apart JADES and MIDIS observations. A substellar Galactic nature is deemed unlikely, given that the Cerberus’s relatively flat NIRCam-to-NIRCam and very red NIRCam-to-MIRI flux ratios are not consistent with any brown dwarf model. The extragalactic nature of Cerberus offers three possibilities: (1) a z ∼ 0.4 galaxy with strong emission from polycyclic aromatic hydrocarbons—the very low inferred stellar mass, M ⋆ = 105-106 M ⊙, makes this possibility highly improbable; (2) a dusty galaxy at z ∼ 4 with an inferred stellar mass M ⋆ ∼ 108 M ⊙; and (3) a galaxy with observational properties similar to those of the reddest little red dots discovered around z ∼ 7, but Cerberus lying at z ∼ 15, with the rest-frame optical dominated by emission from a dusty torus or a dusty starburst.

Dusty star-forming galaxies emit most of their light at far-infrared to millimeter wavelengths as their star formation is highly obscured. Far-infrared and millimeter observations have revealed their dust, neutral and molecular gas properties. The sensitivity of JWST at rest-frame optical and near-infrared wavelengths now allows the study of the stellar and ionized gas content. We investigate the spatially resolved distribution and kinematics of the ionized gas in GN20, a dusty star-forming galaxy at z= 4.0548. We present deep MIRI/MRS integral field spectroscopy of the near-infrared rest-frame emission of GN20. We detect spatially resolved Paα, out to a radius of 6 kpc, distributed in a clumpy morphology. The star formation rate derived from Paα (144 ± 9 M yr1) is only 7.7 ± 0.5% of the infrared star formation rate (1860 ± 90 M yr1). We attribute this to very high extinction (AV = 17.2 ± 0.4 mag, or AV, mixed = 44 ± 3 mag), especially in the nucleus of GN20, where only faint Paα is detected, suggesting a deeply buried starburst. We identify four, spatially unresolved, clumps in the Paα emission. Based on the double peaked Paα profile, we find that each clump consists of at least two sub-clumps. We find mass upper limits consistent with them being formed in a gravitationally unstable gaseous disk. The ultraviolet bright region of GN20 does not have any detected Paα emission, suggesting an age of more than 10 Myr for this region of the galaxy. From the rotation profile of Paα, we conclude that the gas kinematics are rotationally dominated and the vrot/I m = 3.8 ± 1.4 is similar to low-redshift luminous infrared galaxies. From the Paα kinematics, we cannot distinguish between a rotational profile of a large disk and a late stage merger mimicking a disk. We speculate that GN20 is in the late stage of a major merger, where the clumps in a large gas-rich disk are created by the major merger, while the central starburst is driven by the merger event.

We investigate the emergence phase of young star clusters in the nearby spiral galaxy NGC 628. We use JWST NIRCam and MIRI observations to create spatially resolved maps of the Paα 1.87 μm and Brα 4.05 μm hydrogen recombination lines, as well as 3.3 and 7.7 μm emission from polycyclic aromatic hydrocarbons (PAHs). We extract 953 compact H ii regions and analyze the PAH emission and morphology at ∼10 pc scales in the associated photodissociation regions. While H ii regions remain compact, radial profiles help us to define three PAH morphological classes: compact (∼42%), extended (∼34%), and open (∼24%). The majority of compact and extended PAH morphologies are associated with very young star clusters (<5 Myr), while open PAH morphologies are mainly associated with star clusters older than 3 Myr. We observe a general decrease in the 3.3 and 7.7 μm PAH band emission as a function of cluster age, while their ratio remains constant with age out to 10 Myr and morphological class. The recovered PAH3.3μm/PAH7.7μm ratio is lower than values reported in the literature for reference models that consider neutral and ionized PAH populations and analyses conducted at galactic physical scales. The 3.3 and 7.7 μm bands are typically associated with neutral and ionized PAHs, respectively. While we expected neutral PAHs to be suppressed in proximity to an ionizing source, the constant PAH3.3μm/PAH7.7μm ratio would indicate that both families of molecules disrupt at similar rates in proximity to H ii regions.

We present maps of ionized gas (traced by Paα and Brα) and 3.3 μm polycyclic aromatic hydrocarbon (PAH) emission in the nearby spiral galaxy NGC 628, derived from new JWST/NIRCam data from the Feedback in Emerging extrAgalactic Star clusTers (FEAST) survey. With this data, we investigate and calibrate the relation between 3.3 μm PAH emission and star formation rate (SFR) in and around emerging young star clusters (eYSCs) on a scale of ∼40 pc. We find a tight (correlation coefficient ρ ∼ 0.9) sublinear (power-law exponent α ∼ 0.75) relation between the 3.3 μm PAH luminosity surface density and SFR traced by Brα for compact, cospatial (within 0.″16 or ∼7 pc) peaks in Paα, Brα, and 3.3 μm (eYSC-I). The scatter in the relationship does not correlate well with variations in local interstellar medium metallicity, due to a radial metallicity gradient, but rather is likely due to stochastic sampling of the stellar initial mass function (IMF) and variations in the PAH heating and age of our sources. The deviation from a linear relation may be explained by PAH destruction in more intense ionizing environments, variations in age, and IMF stochasticity at intermediate to low luminosities. We test our results with various continuum subtraction techniques using combinations of NIRCam bands and find that they remain robust with only minor differences in the derived slope and intercept. An unexpected discrepancy is identified between the relations of hydrogen recombination lines (Paα versus Brα; Hα versus Brα).

We report the first results from a deep near-infrared campaign with the Hubble Space Telescope to obtain late-epoch images of the Hubble Ultra Deep Field, 10-15 yr after the first epoch data were obtained. The main objectives are to search for faint active galactic nuclei (AGN) at high redshifts by virtue of their photometric variability and measure (or constrain) the comoving number density of supermassive black holes (SMBHs), n SMBH, at early times. In this Letter, we present an overview of the program and preliminary results concerning eight objects. Three variables are supernovae, two of which are apparently hostless with indeterminable redshifts, although one has previously been recorded as a z ≈ 6 object precisely because of its transient nature. Two further objects are clear AGN at z = 2.0 and 3.2, based on morphology and/or infrared spectroscopy from JWST. Three variable targets are identified at z = 6-7 that are also likely AGN candidates. These sources provide a first measure of n SMBH in the reionization epoch by photometric variability, which places a firm lower limit of 3 × 10−4 cMpc−3. After accounting for variability and luminosity incompleteness, we estimate n SMBH ≳ 8 × 10−3 cMpc−3, which is the largest value so far reported at these redshifts. This SMBH abundance is also strikingly similar to estimates of n SMBH in the local Universe. We discuss how these results test various theories for SMBH formation.