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.
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
ABSTRACT -
We use FIRE-2 zoom simulations of Milky Way size disk galaxies to derive easy-to-use relationships between the observed circular speed of the Galaxy at the Solar location,more » « less
v c, and dark matter properties of relevance for direct detection experiments: the dark matter density, the dark matter velocity dispersion, and the speed distribution of dark matter particles near the Solar location. We find that both the local dark matter density and 3D velocity dispersion follow tight power laws withv c. Using this relation together with the observed circular speed of the Milky Way at the Solar radius, we infer the local dark matter density and velocity dispersion near the Sun to beρ = 0.42±0.06 GeV cm-3andσ3D= 280+19-18km s-1. We also find that the distribution of dark matter particle speeds is well-described by a modified Maxwellian with two shape parameters, both of which correlate with the observedv c. We use that modified Maxwellian to predict the speed distribution of dark matter near the Sun and find that it peaks at a most probable speed of 257 km s-1and begins to truncate sharply above 470 km s-1. This peak speed is somewhat higher than expected from the standard halo model, and the truncation occurs well below the formal escape speed to infinity, with fewer very-high-speed particles than assumed in the standard halo model.Free, publicly-accessible full text available August 1, 2025 -
ABSTRACT 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.
-
Abstract 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 *∼ 105–109M ⊙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 cm2g−1. On the contrary, models with (σ /m ) = 10 cm2g−1produce 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 cm2g−1are enhanced by about a factor of 2, which are still consistent with the kinematic measurements of Local Group dwarfs but in tension with the Hi rotation 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 cm2g−1,f diss= 0.75 are in tension in massive dwarfs (M halo∼ 1011M ⊙) 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.Free, publicly-accessible full text available April 30, 2025 -
ABSTRACT 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.
-
Abstract 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 , , and ∇[Fe/H]= −0.23 ± 0.15 dex . 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 , , and ∇[Fe/H]= −0.46 ± 0.10 dex , 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 ≲ 106L ⊙) 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. -
ABSTRACT 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.
-
ABSTRACT 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.
-
ABSTRACT 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.
-
Abstract Early data from the James Webb Space Telescope (JWST) have revealed a bevy of high-redshift galaxy candidates with unexpectedly high stellar masses. An immediate concern is the consistency of these candidates with galaxy formation in the standard
Λ CDM cosmological model, wherein the stellar mass (M ⋆) of a galaxy is limited by the available baryonic reservoir of its host dark matter halo. The mass function of dark matter haloes therefore imposes an absolute upper limit on the number densityn (>M ⋆,z ) and stellar mass densityρ ⋆(>M ⋆,z ) of galaxies more massive thanM ⋆at any epochz . Here I show that the most massive galaxy candidates in JWST observations atz ≈ 7–10 lie at the very edge of these limits, indicating an important unresolved issue with the properties of galaxies derived from the observations, how galaxies form at early times inΛ CDM or within this standard cosmology itself.