Search for: All records

Recent observations with JWST have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the big bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox$^{\it HR}$, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments (FIRE) project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox$^{\it HR}$ re-simulates the cosmic volume ($L=22.1$ cMpc) of the original FIREbox run with eight times higher mass resolution ($m_{\rm b}\sim {}7800\, M_\odot$), but with identical physics, down to $z\sim {}6$. FIREbox$^{\it HR}$ predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at $z\sim {}6{\!-\!}14$, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox$^{\it HR}$, the SFE–halo mass relation for intermediate mass haloes ($M_{\rm halo}\sim {}10^9{\!-\!}10^{11}\, {\rm M}_\odot$) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE–halo mass relation lead to a larger contribution from lower mass haloes at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE–halo mass relation inferred from FIREbox$^{\it HR}$ allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at $z\gt 12$ will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn.

Fuzzy dark matter (FDM), comprised of ultralight ($m \sim 10^{-22}\,{\rm eV}$) boson particles, has received significant attention as a viable alternative to cold dark matter (CDM), as it approximates CDM on large scales (${\gtrsim}1$ Mpc) while potentially resolving some of its small-scale problems via kiloparsec-scale quantum interference. However, the most basic FDM model, with one free parameter (the boson mass), is subject to a tension: small boson masses yield the desired cores of dwarf galaxies but underpredict structure in the Lyman-α forest, while large boson masses render FDM effectively identical to CDM. This Catch-22 problem may be alleviated by considering an axion-like particle with attractive particle self-interactions. We simulate an idealized FDM halo with self-interactions parametrized by an energy decay constant $f \sim 10^{15}~\rm {GeV}$ related to the axion symmetry-breaking conjectured to solve the strong-CP problem in particle physics. We observe solitons, a hallmark of FDM, condensing within a broader halo envelope, and find that the density profile and soliton mass depend on self-interaction strength. We propose generalized formulae to extend those from previous works to include self-interactions. We also investigate a critical mass threshold predicted for strong interactions at which the soliton collapses into a compact, unresolved state. We find that the collapse happens quickly, and its effects are initially contained to the central region of the halo.

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.

We analyze the evolution of massive (log₁₀[M_⋆/M_⊙] > 10) galaxies atz∼ 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 withdense basisto 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 forz∼ 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 ∼300 stellar metallicity measurements in two faint M31 dwarf galaxies, Andromeda XVI (M_V= −7.5) and Andromeda XXVIII (M_V= –8.8), derived using metallicity-sensitive calcium H and K narrowband Hubble Space Telescope imaging. These are the first individual stellar metallicities in And XVI (95 stars). Our And XXVIII sample (191 stars) is a factor of ∼15 increase over literature metallicities. For And XVI, we measure$〈\mathrm{[Fe/H]}〉=-{2.17}_{-0.05}^{+0.05}$,${\sigma }_{\mathrm{[Fe/H]}}={0.33}_{-0.07}^{+0.07}$, and ∇_[Fe/H]= −0.23 ± 0.15 dex$R_e^{-1}$. We find that And XVI is more metal-rich than Milky Way ultrafaint dwarf galaxies of similar luminosity, which may be a result of its unusually extended star formation history. For And XXVIII, we measure$〈\mathrm{[Fe/H]}〉=-{1.95}_{-0.04}^{+0.04}$,${\sigma }_{\mathrm{[Fe/H]}}={0.34}_{-0.05}^{+0.05}$, and ∇_[Fe/H]= −0.46 ± 0.10 dex$R_e^{-1}$, placing it on the dwarf galaxy mass–metallicity relation. Neither galaxy has a metallicity distribution function (MDF) with an abrupt metal-rich truncation, suggesting that star formation fell off gradually. The stellar metallicity gradient measurements are among the first for faint (L≲ 10⁶L_⊙) galaxies outside the Milky Way halo. Both galaxies’ gradients are consistent with predictions from the FIRE simulations, where an age–gradient strength relationship is the observational consequence of stellar feedback that produces dark matter cores. We include a catalog for community spectroscopic follow-up, including 19 extremely metal-poor ([Fe/H] < –3.0) star candidates, which make up 7% of And XVI’s MDF and 6% of And XXVIII’s.

We analyze the first cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modeled as a constant fractional energy dissipation (f_diss= 0.75). In this paper, we examine the properties of dwarf galaxies withM_*∼ 10⁵–10⁹M_⊙in both isolation and within Milky Way–mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of disk formation in dSIDM with (σ/m) ≤ 1 cm²g⁻¹. On the contrary, models with (σ/m) = 10 cm²g⁻¹produce puffier stellar distributions that are in tension with the observed size–mass relation. In addition, owing to the steeper central density profiles, the subkiloparsec circular velocities of isolated dwarfs when (σ/m) ≥ 0.1 cm²g⁻¹are enhanced by about a factor of 2, which are still consistent with the kinematic measurements of Local Group dwarfs but in tension with the Hirotation curves of more massive field dwarfs. Meanwhile, for satellites of Milky Way–mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help promote the structural diversity of dwarf satellites. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with (σ/m) ≳ 0.1 cm²g⁻¹,f_diss= 0.75 are in tension in massive dwarfs (M_halo∼ 10¹¹M_⊙) due to circular velocity constraints. However, models with lower effective cross sections (at this halo mass/velocity scale) are still viable and can produce nontrivial observable signatures.

We introduce an analytic surface density profile for dark matter haloes that accurately reproduces the structure of simulated haloes of mass Mvir = 107–1011 M⊙, making it useful for modelling line-of-sight (LOS) perturbers in strong gravitational lensing models. The two-parameter function has an analytic deflection potential and is more accurate than the projected Navarro, Frenk, and White profile commonly adopted at this mass scale for perturbers, especially at the small radii of most relevant for lensing perturbations. Using a characteristic radius, R−1, where the log slope of surface density is equal to −1, and an associated surface density, Σ−1, we can represent the expected lensing signal from LOS haloes statistically, for an ensemble of halo orientations, using a distribution of projected concentration parameters, $\mathcal {C}_{\rm vir} := r_{\rm vir}/ R_{-1}$. Though an individual halo can have a projected concentration that varies with orientation with respect to the observer, the range of projected concentrations correlates with the usual three-dimensional halo concentration in a way that enables ease of use.

The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at $z \gtrsim 6$ is an open question for understanding which galaxies primarily drove cosmic reionization. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil record provides a powerful ‘near-far’ technique for studying the reionization-era SMF/UVLF, probing orders of magnitude lower in mass than direct HST/JWST observations. Using 882 low-mass ($M_{\rm star}\lesssim 10^{9}\, \rm {M_\odot }$) galaxies across 11 Milky Way (MW)- and Local Group-analogue environments from the FIRE-2 cosmological baryonic zoom-in simulations, we characterize their progenitors at $z=6\!-\!9$, the mergers/disruption of those progenitors over time, and how well their present-day fossil record traces the high-redshift SMF. A present-day galaxy with $M_{\rm star}\sim 10^5\, \rm {M_\odot }$ ($\sim 10^9\, \rm {M_\odot }$) had $\approx 1$ ($\approx 30$) progenitors at $z\approx 7$, and its main progenitor comprised $\approx 100~{{\ \rm per\ cent}}$ ($\approx 10~{{\ \rm per\ cent}}$) of the total stellar mass of all its progenitors at $z\approx 7$. We show that although only $\sim 15~{{\ \rm per\ cent}}$ of the early population of low-mass galaxies survives to present day, the fossil record of surviving Local Group galaxies accurately traces the low-mass slope of the SMF at $z \sim 6 \!-\! 9$. We find no obvious mass dependence to the mergers and accretion, and show that applying this reconstruction technique to just low-mass galaxies at $z = 0$ and not the MW/M31 hosts correctly recovers the slope of the SMF down to $M_{\rm star} \sim 10^{4.5}\, \rm {{\rm M}_{\odot }}$ at $z \gtrsim 6$. Thus, we validate the ‘near-far’ approach as an unbiased tool for probing low-mass reionization-era galaxies.

The radial acceleration relation (RAR) connects the total gravitational acceleration of a galaxy at a given radius, atot(r), with that accounted for by baryons at the same radius, abar(r). The shape and tightness of the RAR for rotationally-supported galaxies have characteristics in line with MOdified Newtonian Dynamics (MOND) and can also arise within the cosmological constant + cold dark matter (ΛCDM) paradigm. We use zoom simulations of 20 galaxies with stellar masses of M⋆ ≃ 107–11 M⊙ to study the RAR in the FIRE-2 simulations. We highlight the existence of simulated galaxies with non-monotonic RAR tracks that ‘hook’ down from the average relation. These hooks are challenging to explain in Modified Inertia theories of MOND, but naturally arise in all of our ΛCDM-simulated galaxies that are dark-matter dominated at small radii and have feedback-induced cores in their dark matter haloes. We show, analytically and numerically, that downward hooks are expected in such cored haloes because they have non-monotonic acceleration profiles. We also extend the relation to accelerations below those traced by disc galaxy rotation curves. In this regime, our simulations exhibit ‘bends’ off of the MOND-inspired extrapolation of the RAR, which, at large radii, approach atot ≈ abar/fb, where fb is the cosmic baryon fraction. Future efforts to search for these hooks and bends in real galaxies will provide interesting tests for MOND and ΛCDM.

We report the discovery of 15 exceptionally luminous 10 ≲z≲ 14 candidate galaxies discovered in the first 0.28 deg²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 luminousz≳ 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 genuinez≳ 10 sources and 3/15 (20%) likely low-redshift contaminants. Three of ourz∼ 12 candidates push the limits of early stellar mass assembly: they have estimated stellar masses ∼ 5 × 10⁹M_⊙, implying an effective stellar baryon fraction ofϵ_⋆∼ 0.2−0.5, whereϵ_⋆≡M_⋆/(f_bM_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 forM_⋆∼ 10¹⁰M_⊙galaxies relative toM_⋆∼ 10⁹M_⊙—both about 10⁻⁶Mpc⁻³—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 atz≈ 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.