We report the discovery of two galaxy candidates at redshifts between 15.7 < z < 16.4 in James Webb Space Telescope (JWST) observations from the GLIMPSE survey. These robust sources were identified using a combination of Lyman break selection and photometric redshift estimates. The ultradeep NIRCam imaging from GLIMPSE, combined with the strong gravitational lensing of the AS1063 galaxy cluster, allows us to probe an intrinsically fainter population (down to M_UV = −17.0 mag) than previously achievable. These galaxies have absolute magnitudes ranging from M_UV = −17.0 to −17.2 mag, with blue (β ≃ −2.87) ultraviolet (UV) continuum slopes, consistent with young, dust-free stellar populations. The number density of these objects, log₁₀(ϕ/[Mpc⁻³ mag⁻¹]) = −3.47−0.10+0.13 at M_UV = −17, is in clear tension with pre-JWST theoretical predictions, extending the overabundance of galaxies from z ∼ 10 to z ∼ 17. These results, together with the scarcity of brighter galaxies in other public surveys, suggest a steep decline in the bright end of the UV luminosity function at z ∼ 16, implying efficient star formation and possibly a close connection to the halo mass function at these redshifts. Testing a variety of star formation histories suggests that these sources are plausible progenitors of the unusually UV-bright galaxies that JWST now routinely uncovers at z = 10–14. Overall, our results indicate that the luminosity distribution of the earliest star-forming galaxies could be shifting toward fainter luminosities, implying that future surveys of cosmic dawn will need to explore this faint luminosity regime.

Power spectra (PS) of high-resolution images of M51 (NGC 5194) taken with the Hubble Space Telescope and the James Webb Space Telescope (JWST) have been examined for evidence of disk thickness in the form of a change in slope between large scales, which map two-dimensional correlated structures, and small scales, which map three-dimensional correlated structures. Such a slope change is observed here in Hα, and possibly Paα, using average PS of azimuthal intensity scans that avoid bright peaks. The physical scale of the slope change occurs at ∼120 pc and ∼170 pc for these two transitions, respectively. A radial dependence in the shape of the Hα PS also suggests that the length scale drops from ∼180 pc at 5 kpc, to ∼90 pc at 2 kpc, to ∼25 pc in the central ∼kpc. We interpret these lengths as comparable to the thicknesses of the star-forming disk traced by H ii regions. The corresponding emission measure is ∼100 times larger than what is expected from the diffuse ionized gas. The PS of JWST Mid-IR Instrument images in eight passbands have more gradual changes in slope, making it difficult to determine a specific value of the thickness for this emission.

We present a detailed JWST/NIRSpec and NIRCam analysis of a gravitationally lensed galaxy (μ = 17 - 21) at a redshift of 6.14 magnified by the Hubble Frontier Field galaxy cluster MACS J0416. The target galaxy is a typical compact and UV-faint (MUV = - 17.8) Lyman-α emitter, yet the large magnification allows the detailed characterization of structures on sub-galactic scales (down to a few parsecs). Prominent optical Hα, Hβ, and [OIII]λλ4959, 5007 lines are spatially resolved with the high-spectral-resolution grating (G395H, R 2700), with large equivalent widths, EW(Hβ+OIII) ≳ 1000 Å, and elevated ionizing photon production efficiencies, log(ξion/erg- 1Hz) = 25.2-25.7. NIRCam deep imaging reveals the presence of compact rest-UV-bright regions along with individual star clusters of Reff = 3 - 8 pc in size and M ~ 2 · 105 - 5 · 106 M⊙ in mass. These clusters are characterized by steep UV slopes, βUV ≲ - 2.5, which in some cases are associated with a dearth of line emission, indicating possible leaking of the ionizing radiation, as also supported by a Lyman-α emission peaking at ~100 km s- 1 from the systemic redshift. While the entire system is characterized by low metallicity, ~0.1 Z⊙, the NIRSpec-IFU map also reveals the presence of a low-luminosity, metal-poor region with Z ≲ 2% Z⊙, which is barely detected in NIRCam imaging; this region is displaced by > 200 pc from one of the brightest structures of the system in UV, and would have been too faint to detect if not for the large magnification of the system.

The simple stellar population models produced by stellar population and spectral synthesis (SPS) codes are used as spectral templates in a variety of astrophysical contexts. In this paper, we test the predictions of four commonly used stellar population synthesis codes (YGGDRASIL, BPASS, FSPS, and a modified form of GALAXEV, which we call GALAXEVneb) by using them as spectral templates for photometric spectral energy distribution (SED) fitting with a sample of 18 young stellar clusters. All clusters have existing Hubble Space Telescope (HST) Cosmic Origins Spectrograph far-UV spectroscopy that provides constraints on their ages as well as broadband photometry from HST Advanced Camera for Surveys and Wide Field Camera 3. We use model spectra that account for both nebular and stellar emission, and additionally test four extinction curves at different values of R V . We find that for individual clusters, choice of extinction curve and SPS model can introduce significant scatter into the results of SED fitting. Model choice can introduce scatter of 34.8 Myr in age, a factor of 9.5 in mass, and 0.40 mag in extinction. Extinction curve choice can introduce scatter of up to a factor of 32.3 Myr in age, a factor of 10.4 in mass, and 0.41 mag in extinction. We caution that because of this scatter, one-to-one comparisons between the properties of individual objects derived using different SED fitting setups may not be meaningful. However, our results also suggest that SPS model and extinction curve choice do not introduce major systematic differences into SED fitting results when the entire cluster population is considered. The distribution of cluster properties for a large enough sample is relatively robust to user choice of SPS code and extinction curve.

Deuteration of hydrocarbon material, including polycyclic aromatic hydrocarbons (PAHs), has been proposed to account for the low gas-phase abundances of D in the interstellar medium (ISM). JWST spectra of four star-forming regions in M51 show an emission feature, with central wavelength ∼4.647 μm and FWHM 0.0265 μm, corresponding to the C-D stretching mode in aliphatic hydrocarbons. The emitting aliphatic material is estimated to have (D/H)aliph. ≈0.17 ± 0.02—a factor of ∼104 enrichment relative to the overall ISM. On ∼50 pc scales, deuteration levels toward four H ii regions in M51 are 2-3 times higher than in the Orion Bar photodissociation region (PDR), with implications for the processes responsible for the formation and evolution of hydrocarbon nanoparticles, including PAHs. The deuteration of the aliphatic material is found to anticorrelate with helium ionization in the associated H ii, suggesting that harsh far-UV radiation may act to lower the deuteration of aliphatics in PDRs near massive stars. No evidence is found for deuteration of aromatic material, with (D/H)arom. ≲ 0.016: deuteration of the aliphatic material exceeds that of the aromatic material by at least a factor of 10. The observed levels of deuteration may account for the depletion of D observed in the Galactic ISM. If so, the 4.65 μm feature may be detectable in absorption.

JWST observations have unveiled faint active galactic nuclei (AGNs) at high redshift that provide insights into the formation of supermassive black holes (SMBHs). However, disentangling their stellar from AGN light is challenging. Here, we use an empirical approach to infer the average stellar mass of five faint broad-line (BL) Hα emitters at z = 4–5 with BH masses ≈6 × 10⁶ M_⊙, with a method independent of their spectral energy distribution (SED). We use the deep JWST/NIRcam grism survey “All the Little Things” to measure the overdensities around BL-Hα emitters and around a spectroscopic reference sample of ∼300 galaxies. In our reference sample, we find that megaparsec-scale overdensity correlates with stellar mass. Their large-scale environments suggest that BL-Hα emitters are hosted by galaxies with stellar masses ≈5 × 10⁷ M_⊙, ≈40 times lower than those inferred from galaxy-only SED fits. Adding measurements around more luminous z ≈ 6 AGNs, we find tentative correlations between line width, BH mass, and the overdensity, suggestive of a steep BH to halo mass relation. The main implications are (1) when BH masses are taken at face value, we confirm extremely high BH to stellar mass ratios of ≈10%, (2) the galaxies of low stellar mass that host growing SMBHs are in tension with typical hydrodynamical simulations, except those without feedback, (3) a 1% duty cycle implied by the host mass hints at super-Eddington accretion, (4) the masses are at odds with an interpretation of the line broadening in terms of high stellar density, (5) our results imply a luminosity-dependent diversity of galaxy masses, environments, and SEDs among AGN samples.

We present new JWST/NIRCam observations of the starburst irregular galaxy NGC 4449, obtained in Cycle 1 as part of the Feedback in Emerging extrAgalactic Star clusTers program, which we use to investigate its resolved stellar populations and their spatial distributions. NGC 4449 near-IR color-magnitude diagrams reveal a broad range of stellar populations, spanning different evolutionary phases, from young main sequence stars, to old red giant branch stars and asymptotic giant branch (AGB) stars. The analysis of their spatial distributions shows that younger (≤10 Myr) populations form an S-shaped distribution aligned with the galaxy’s north-south axis, while stars aged 10-60 Myr show shifting concentrations from the north to the south, consistent with the possibility that external interactions or tidal effects may have triggered star formation in spatially distinct bursts. Clusters of comparable ages generally follow these distributions, suggesting that cluster and field stars form at the same pace in each galaxy region. Thanks to the unprecedented high-spatial resolution and sensitivity of the JWST data, we recover a clear gap between oxygen-rich and the carbon star branch of the AGB population, as well as the presence of a massive AGB star “finger.” The analysis of these stars can provide constraints on AGB evolution models and dust production in this galaxy. These results confirm NGC 4449's status as a compelling example of a local dwarf starburst galaxy undergoing complex and possibly externally driven star formation and underscore the power of JWST in probing the full lifecycle of stars in nearby starburst systems.

Haro 11 is a metal-poor, starburst galaxy believed to be the result of an ongoing merger, which is shaping the properties of the galaxy. In this study, we carry out a large suite of numerical simulations of a merger between two disc galaxies, to study possible origins of Haro 11 and understand under which conditions various features of the galaxy are formed. By varying galaxy parameters describing the orbital configurations, masses, and their inclination, we perform a total of ∼500 simulations. We demonstrate that a two-disc galaxy merger reproduces key, observed features of Haro 11, including its morphology, gas kinematics, star formation history, and stellar population ages and masses. In particular, we present a fiducial Haro 11 model that produces the single observed tidal tail, three stellar knots, and inner gas morphology and kinematics. The resulting orbit and galactic morphology are robust against small variations of the initial parameters. By performing mock observations, we compare with the results of observational data and discuss possible origins for various features. Furthermore, we present newly gathered observational data that confirms the presence of a stellar tidal tail with similar length and morphology as our simulations.

Haro 11 is the closest known Lyman continuum leaking galaxy and serves as an important laboratory for studying the escape of Lyman continuum radiation. The galaxy is a metal-poor, starburst galaxy believed to be undergoing a merger that might help facilitate the escape of radiation. In this study, we carry out a large suite of numerical simulations of a merger between two disc galaxies, to study possible origins of Haro 11 and understand under which conditions various features of the galaxy are formed. By varying galaxy parameters describing the orbital configurations, masses, and their inclination, we perform a total of ~500 simulations. We demonstrate that a two-disc galaxy merger is able to reproduce key, observed features of Haro 11, including its morphology, gas kinematics, star formation history, and stellar population ages and masses. We also find that small parameter variations have minimal impact on the orbits and resulting galaxy properties. In particular, we present a fiducial Haro 11 model that produces the single observed tidal tail, the presence of three stellar knots, and inner gas morphology and kinematics. By performing mock observations, we compare with the results of observational data and discuss possible origins for various features. Furthermore, we present newly gathered observational data that confirms the presence of a stellar tidal tail with similar length and direction as our simulations.

Detecting the first generation of stars, Population III (Pop III), has been a long-standing goal in astrophysics, yet they remain elusive even in the JWST era. Here we present a novel NIRCam-based selection method for Pop III galaxies, and carefully validate it through completeness and contamination simulations. We systematically search similar or equal to 500 arcmin(2) across JWST legacy fields for Pop III candidates, including GLIMPSE, which, assisted by gravitational lensing, has produced JWST's deepest NIRCam imaging thus far. We discover one promising Pop III galaxy candidate (GLIMPSE-16043) at z = 6.50(-0.24 )(+0.03) , a moderately lensed galaxy ( mu = 2.9(-0.2)(+0.1) ) with an intrinsic UV magnitude of = M-UV = -15.89 (+0.12)(-0.14) . It exhibits key Pop III features: strong H alpha emission (rest-frame EW 2810 +/- 550 & Aring;); a Balmer jump; no dust (UV slope beta = -2.34 +/- 0.36); and undetectable metal lines (e.g., [O III]; [O III]/H beta < 0.44), implying a gas-phase metallicity of Z(gas)/Z(circle dot) < 0.5%. These properties indicate the presence of a nascent, metal-deficient young stellar population (<5 Myr) with a stellar mass of similar or equal to 10(5) M-circle dot. Intriguingly, this source deviates significantly from the extrapolated UV-metallicity relation derived from recent JWST observations at z = 4-10, consistent with UV enhancement by a top-heavy Pop III initial mass function or the presence of an extremely metal-poor active galactic nucleus. We also derive the first observational constraints on the Pop III UV luminosity function at z similar or equal to 6-7. The volume density of GLIMPSE-16043 (approximate to 10(-4) cMpc(-3) ) is in excellent agreement with theoretical predictions, independently reinforcing its plausibility. This study demonstrates the power of our novel NIRCam method to finally reveal distant galaxies even more pristine than the Milky Way's most metal-poor satellites, thereby promising to bring us closer to the first generation of stars than we have ever been before