Search for: All records

ABSTRACT JWST has revealed a large population of UV-bright galaxies at $$z\gtrsim 10$$ and possibly overly massive galaxies at $$z\gtrsim 7$$, challenging standard galaxy formation models in the ΛCDM cosmology. We use an empirical galaxy formation model to explore the potential of alleviating these tensions through an Early Dark Energy (EDE) model, originally proposed to solve the Hubble tension. Our benchmark model demonstrates excellent agreement with the UV luminosity functions (UVLFs) at $$4\lesssim z \lesssim 10$$ in both ΛCDM and EDE cosmologies. In the EDE cosmology, the UVLF measurements at $$z\simeq 12$$ based on spectroscopically confirmed galaxies (eight galaxies at $$z\simeq 11\!-\!13.5$$) exhibit no tension with the benchmark model. Photometric constraints at $$12 \lesssim z\lesssim 16$$ can be fully explained within EDE via either moderately increased star-formation efficiencies ($$\epsilon _{\ast}\sim 3\!-\!10\ \hbox{per cent}$$ at $$M_{\rm halo}\sim 10^{10.5}{\, \rm M_\odot }$$) or enhanced UV variabilities ($$\sigma _{\rm UV}\sim 0.8\!-\!1.3$$ mag at $$M_{\rm halo}\sim 10^{10.5}{\, \rm M_\odot }$$) that are within the scatter of hydrodynamical simulation predictions. A similar agreement is difficult to achieve in $$\Lambda$$CDM, especially at $$z\gtrsim 14$$, where the required $$\sigma _{\rm UV}$$ exceeds the maximum value seen in simulations. Furthermore, the implausibly large cosmic stellar mass densities inferred from some JWST observations are no longer in tension with cosmology when the EDE is considered. Our findings highlight EDE as an intriguing unified solution to a fundamental problem in cosmology and the recent tensions raised by JWST observations. Data at the highest redshifts reached by JWST will be crucial for differentiating modified galaxy formation physics from new cosmological physics. 

ABSTRACT JWST observations have revealed a population of galaxies bright enough that potentially challenge standard galaxy formation models in the Λ cold dark matter (ΛCDM) cosmology. Using a minimal empirical framework, we investigate the influence of variability on the rest-frame ultra-violet (UV) luminosity function of galaxies at z ≥ 9. Our study differentiates between the median UV radiation yield and the variability of UV luminosities of galaxies at a fixed dark matter halo mass. We primarily focus on the latter effect, which depends on halo assembly and galaxy formation processes and can significantly increase the abundance of UV-bright galaxies due to the upscatter of galaxies in lower-mass haloes. We find that a relatively low level of variability, σUV ≈ 0.75 mag, matches the observational constraints at z ≈ 9. However, increasingly larger σUV is necessary when moving to higher redshifts, reaching $$\sigma _{\rm UV} \approx 2.0\, (2.5)\, {\rm mag}$$ at z ≈ 12 (16). This implied variability is consistent with expectations of physical processes in high-redshift galaxies such as bursty star formation and dust clearance during strong feedback cycles. Photometric constraints from JWST at z ≳ 9 therefore can be reconciled with a standard ΛCDM-based galaxy formation model calibrated at lower redshifts without the need for adjustments to the median UV radiation yield. 

ABSTRACT Using high-resolution cosmological radiation-hydrodynamic (RHD) simulations (thesan-hr), we explore the impact of alternative dark matter (altDM) models on galaxies during the Epoch of Reionization. The simulations adopt the IllustrisTNG galaxy formation model. We focus on altDM models that exhibit small-scale suppression of the matter power spectrum, namely warm dark matter (WDM), fuzzy dark matter (FDM), and interacting dark matter (IDM) with strong dark acoustic oscillations (sDAO). In altDM scenarios, both the halo mass functions and the ultraviolet luminosity functions at z ≳ 6 are suppressed at the low-mass/faint end, leading to delayed global star formation and reionization histories. However, strong non-linear effects enable altDM models to ‘catch up’ with cold dark matter (CDM) in terms of star formation and reionization. The specific star formation rates are enhanced in halos below the half-power mass in altDM models. This enhancement coincides with increased gas abundance, reduced gas depletion times, more compact galaxy sizes, and steeper metallicity gradients at the outskirts of the galaxies. These changes in galaxy properties can help disentangle altDM signatures from a range of astrophysical uncertainties. Meanwhile, it is the first time that altDM models have been studied in RHD simulations of galaxy formation. We uncover significant systematic uncertainties in reionization assumptions on the faint-end luminosity function. This underscores the necessity of accurately modeling the small-scale morphology of reionization in making predictions for the low-mass galaxy population. Upcoming James Webb Space Telescope imaging surveys of deep lensed fields hold potential for uncovering the faint low-mass galaxy population, which could provide constraints on altDM models. 

ABSTRACT JWST has recently sparked a new era of Lyα spectroscopy, delivering the first measurements of the Lyα escape fraction and velocity profile in typical galaxies at z ≃ 6−10. These observations offer new prospects for insight into the earliest stages of reionization. But to realize this potential, we need robust models of Lyα properties in galaxies at z ≃ 5−6 when the IGM is mostly ionized. Here, we use new JWST observations from the JADES and FRESCO surveys combined with VLT/MUSE and Keck/DEIMOS data to characterize statistical distributions of Lyα velocity offsets, escape fractions, and EWs in z ≃ 5−6 galaxies. We find that galaxies with large Lyα escape fractions (>0.2) are common at z ≃ 5−6, comprising 30 per cent of Lyman break selected samples. Comparing to literature studies, our census suggests that Lyα becomes more prevalent in the galaxy population towards higher redshift from z ∼ 3 to z ∼ 6, although we find that this evolution slows considerably between z ∼ 5 and z ∼ 6, consistent with modest attenuation from residual H i in the mostly ionized IGM at z ≃ 5−6. We find significant evolution in Lyα velocity profiles between z ≃ 2−3 and z ≃ 5−6, likely reflecting the influence of resonant scattering from residual intergalactic H i on the escape of Lyα emission near line centre. This effect will make it challenging to use Lyα peak offsets as a probe of Lyman continuum leakage at z ≃ 5−6. We use our z ≃ 5−6 Lyα distributions to make predictions for typical Lyα properties at z ≳ 8 and discuss implications of a recently discovered Lyα emitter at z ≃ 8.5 with a small peak velocity offset (156 km s−1). 

We present ten novel [OIII]λ4363 auroral line detections up toz ∼ 9.5 measured from ultra-deep JWST/NIRSpec MSA spectroscopy from the JWST Advanced Deep Extragalactic Survey (JADES). We leverage the deepest spectroscopic observations taken thus far with NIRSpec to determine electron temperatures and oxygen abundances using the directT_emethod. We directly compare these results against a suite of locally calibrated strong-line diagnostics and recent high-zcalibrations. We find the calibrations fail to simultaneously match our JADES sample, thus warranting a self-consistent revision of these calibrations for the high-zUniverse. We find a weak dependence between R2 and O3O2 with metallicity, thus suggesting these line ratios are inefficient in the high-zUniverse as metallicity diagnostics and degeneracy breakers. We find R3 and R23 are still correlated with metallicity, but we find a tentative flattening of these diagnostics, thus suggesting future difficulties when applying these strong line ratios as metallicity indicators in the high-zUniverse. We also propose and test an alternative diagnostic based on a different combination of R3 and R2 with a higher dynamic range. We find a reasonably good agreement (median offset of 0.002 dex, median absolute offset of 0.13 dex) with the JWST sample at low metallicity, but future investigations are required on larger samples to probe past the turnover point. At a given metallicity, our sample demonstrates higher ionization and excitation ratios than local galaxies with rest-frame EWs(Hβ) ≈200 − 300 Å. However, we find the median rest-frame EWs(Hβ) of our sample to be ∼2× less than the galaxies used for the local calibrations. This EW discrepancy combined with the high ionization of our galaxies does not offer a clear description of [OIII]λ4363 production in the high-zUniverse, thus warranting a much deeper examination into the factors influencing these processes. 

ABSTRACT Hydrogen emission lines can provide extensive information about star-forming galaxies in both the local and high-redshift Universe. We present a detailed Lyman continuum (LyC), Lyman-α (Lyα), and Balmer line (Hα and Hβ) radiative transfer study of a high-resolution isolated Milky Way simulation using the state-of-the-art Arepo-RT radiation hydrodynamics code with the SMUGGLE galaxy formation model. The realistic framework includes stellar feedback, non-equilibrium thermochemistry accounting for molecular hydrogen, and dust grain evolution in the interstellar medium (ISM). We extend our publicly available Cosmic Lyα Transfer (COLT) code with photoionization equilibrium Monte Carlo radiative transfer and various methodology improvements for self-consistent end-to-end (non-)resonant line predictions. Accurate LyC reprocessing to recombination emission requires modelling pre-absorption by dust ($$f_\text{abs} \approx 27.5\,\rm{per\,\,cent}$$), helium ionization ($$f_\text{He} \approx 8.7\,\rm{per\,\,cent}$$), and anisotropic escape fractions ($$f_\text{esc} \approx 7.9\,\rm{per\,\,cent}$$), as these reduce the available budget for hydrogen line emission ($$f_\text{H} \approx 55.9\,\rm{per\,\,cent}$$). We investigate the role of the multiphase dusty ISM, disc geometry, gas kinematics, and star formation activity in governing the physics of emission and escape, focusing on the time variability, gas-phase structure, and spatial spectral, and viewing angle dependence of the emergent photons. Isolated disc simulations are well-suited for comprehensive observational comparisons with local Hα surveys, but would require a proper cosmological circumgalactic medium (CGM) environment as well as less dust absorption and rotational broadening to serve as analogs for high-redshift Lyα emitting galaxies. Future applications of our framework to next-generation cosmological simulations of galaxy formation including radiation-hydrodynamics that resolve ≲10 pc multiphase ISM and ≲1 kpc CGM structures will provide crucial insights and predictions for current and upcoming Lyα observations. 

Abstract Several lines of evidence suggest that the Milky Way underwent a major merger at z ∼ 2 with the Gaia-Sausage-Enceladus (GSE) galaxy. Here we use H3 Survey data to argue that GSE entered the Galaxy on a retrograde orbit based on a population of highly retrograde stars with chemistry similar to the largely radial GSE debris. We present the first tailored N -body simulations of the merger. From a grid of ≈500 simulations we find that a GSE with M ⋆ = 5 × 10 8 M ⊙ , M DM = 2 × 10 11 M ⊙ best matches the H3 data. This simulation shows that the retrograde stars are stripped from GSE’s outer disk early in the merger. Despite being selected purely on angular momenta and radial distributions, this simulation reproduces and explains the following phenomena: (i) the triaxial shape of the inner halo, whose major axis is at ≈35° to the plane and connects GSE’s apocenters; (ii) the Hercules-Aquila Cloud and the Virgo Overdensity, which arise due to apocenter pileup; and (iii) the 2 Gyr lag between the quenching of GSE and the truncation of the age distribution of the in situ halo, which tracks the lag between the first and final GSE pericenters. We make the following predictions: (i) the inner halo has a “double-break” density profile with breaks at both ≈15–18 kpc and 30 kpc, coincident with the GSE apocenters; and (ii) the outer halo has retrograde streams awaiting discovery at >30 kpc that contain ≈10% of GSE’s stars. The retrograde (radial) GSE debris originates from its outer (inner) disk—exploiting this trend, we reconstruct the stellar metallicity gradient of GSE (−0.04 ± 0.01 dex r 50 − 1 ). These simulations imply that GSE delivered ≈20% of the Milky Way’s present-day dark matter and ≈50% of its stellar halo. 

ABSTRACT The nebular recombination line H α is widely used as a star formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H α radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H α emission, and its connection to the underlying gas and star formation properties. The H α and H β radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H α emission from collisional excitation amounts to fcol ∼ 5–$$10{{\ \rm per\ cent}}$$, only weakly dependent on radius and vertical height, and that scattering boosts the H α luminosity by $$\sim 40{{\ \rm per\ cent}}$$. The dust correction via the Balmer decrement works well (intrinsic H α emission recoverable within 25 per cent), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H α–SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $$f_{\rm abs}\approx 28{{\ \rm per\ cent}}$$ and $$f_{\rm He}\approx 9{{\ \rm per\ cent}}$$, respectively. Together with an escape fraction of $$f_{\rm esc}\approx 6{{\ \rm per\ cent}}$$, this reduces the available budget for hydrogen line emission by nearly half ($$f_{\rm H}\approx 57{{\ \rm per\ cent}}$$). We discuss the impact of the diffuse ionized gas, showing – among other things – that the extraplanar H α emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope. 

ABSTRACT We post-process galaxies in the IllustrisTNG simulations with skirt radiative transfer calculations to make predictions for the rest-frame near-infrared (NIR) and far-infrared (FIR) properties of galaxies at z ≥ 4. The rest-frame K- and z-band galaxy luminosity functions from TNG are overall consistent with observations, despite $${\sim}0.5\, \mathrm{dex}$$ underprediction at z = 4 for MK ≲ −25 and Mz ≲ −24. Predictions for the JWST MIRI observed galaxy luminosity functions and number counts are given. Based on theoretical estimations, we show that the next-generation survey conducted by JWST can detect 500 (30) galaxies in F1000W in a survey area of $$500\, {\rm arcmin}^{2}$$ at z = 6 (z = 8). As opposed to the consistency in the UV, optical, and NIR, we find that TNG, combined with our dust modelling choices, significantly underpredicts the abundance of most dust-obscured and thus most luminous FIR galaxies. As a result, the obscured cosmic star formation rate density (SFRD) and the SFRD contributed by optical/NIR dark objects are underpredicted. The discrepancies discovered here could provide new constraints on the sub-grid feedback models, or the dust contents, of simulations. Meanwhile, although the TNG predicted dust temperature and its relations with IR luminosity and redshift are qualitatively consistent with observations, the peak dust temperature of z ≥ 6 galaxies are overestimated by about $$20\, {\rm K}$$. This could be related to the limited mass resolution of our simulations to fully resolve the porosity of the interstellar medium (or specifically its dust content) at these redshifts. 

Abstract The astrophysical origins of r -process elements remain elusive. Neutron star mergers (NSMs) and special classes of core-collapse supernovae (rCCSNe) are leading candidates. Due to these channels’ distinct characteristic timescales (rCCSNe: prompt, NSMs: delayed), measuring r -process enrichment in galaxies of similar mass but differing star formation durations might prove informative. Two recently discovered disrupted dwarfs in the Milky Way’s stellar halo, Kraken and Gaia-Sausage Enceladus (GSE), afford precisely this opportunity: Both have M ⋆ ≈ 10 8 M ⊙ but differing star formation durations of ≈2 Gyr and ≈3.6 Gyr. Here we present R ≈ 50,000 Magellan/MIKE spectroscopy for 31 stars from these systems, detecting the r -process element Eu in all stars. Stars from both systems have similar [Mg/H] ≈ −1, but Kraken has a median [Eu/Mg] ≈ −0.1 while GSE has an elevated [Eu/Mg] ≈ 0.2. With simple models, we argue NSM enrichment must be delayed by 500–1000 Myr to produce this difference. rCCSNe must also contribute, especially at early epochs, otherwise stars formed during the delay period would be Eu free. In this picture, rCCSNe account for ≈50% of the Eu in Kraken, ≈25% in GSE, and ≈15% in dwarfs with extended star formation durations like Sagittarius. The inferred delay time for NSM enrichment is 10×–100× longer than merger delay times from stellar population synthesis—this is not necessarily surprising because the enrichment delay includes time taken for NSM ejecta to be incorporated into subsequent generations of stars. For example, this may be due to natal kicks that result in r -enriched material deposited far from star-forming gas, which then takes ≈10 8 –10 9 yr to cool in these galaxies.