An official website of the United States government
Abstract A key obstacle to accurate models of the first stars and galaxies is the vast range of distance scales that must be considered. While star formation occurs on sub-parsec scales within dark matter (DM) minihalos, it is influenced by large-scale baryon-dark matter streaming velocities (vbc) and Lyman-Werner (LW) radiative feedback which vary significantly on scales of ∼100 Mpc. We present a novel approach to this issue in which we utilize artificial neural networks (NNs) to emulate the Population III (PopIII) and Population II (PopII) star formation histories of many small-scale cells given by a more complex semi-analytic framework based on DM halo merger trees. Within each simulation cell, the NN takes a set of input parameters that depend on the surrounding large-scale environment, such as the cosmic overdensity,δ(x⃗), andvbcof the cell, then outputs the resulting star formation far more efficiently than is possible with the semi-analytic model. This rapid emulation allows us to self-consistently determine the LW background intensity on ∼100 Mpc scales, while simultaneously including the detailed merger histories (and corresponding star formation histories) of the low-mass minihalos that host the first stars. Comparing with the full semi-analytic framework utilizing DM halo merger trees, our NN emulators yield star formation histories with redshift-averaged errors of ∼7.3% and ∼5.2% for PopII and PopIII, respectively. When compared to a simpler sub-grid star formation prescription reliant on halo mass function integration, we find that the diversity of halo merger histories in our simulation leads to enhanced spatial fluctuations, an earlier transition from PopIII to PopII dominated star formation, and more scatter in star formation histories overall.
more »
« less
A Global Semianalytic Model of the First Stars and Galaxies Including Dark Matter Halo Merger Histories
Abstract We present a new self-consistent semianalytic model of the first stars and galaxies to explore the high-redshift (z≥ 15) Population III (PopIII) and metal-enriched star formation histories. Our model includes the detailed merger history of dark matter halos generated with Monte Carlo merger trees. We calibrate the minimum halo mass for PopIII star formation from recent hydrodynamical cosmological simulations that simultaneously include the baryon–dark matter streaming velocity, Lyman–Werner (LW) feedback, and molecular hydrogen self-shielding. We find an overall increase in the resulting star formation rate density (SFRD) compared to calibrations based on previous simulations (e.g., the PopIII SFRD is over an order of magnitude higher atz= 35−15). We evaluate the effect of the halo-to-halo scatter in this critical mass and find that it increases the PopIII stellar mass density by a factor ∼1.5 atz≥ 15. Additionally, we assess the impact of various semianalytic/analytic prescriptions for halo assembly and star formation previously adopted in the literature. For example, we find that models assuming smooth halo growth computed via abundance matching predict SFRDs similar to the merger tree model for our fiducial model parameters, but that they may underestimate the PopIII SFRD in cases of strong LW feedback. Finally, we simulate subvolumes of the Universe with our model both to quantify the reduction in total star formation in numerical simulations due to a lack of density fluctuations on spatial scales larger than the simulation box, and to determine spatial fluctuations in SFRD due to the diversity in halo abundances and merger histories.
more »
« less
- Award ID(s):
- 2009309
- PAR ID:
- 10489881
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 962
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 62
- Size(s):
- Article No. 62
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)ABSTRACT Cosmic variance is the intrinsic scatter in the number density of galaxies due to fluctuations in the large-scale dark matter density field. In this work, we present a simple analytic model of cosmic variance in the high-redshift Universe (z ∼ 5–15). We assume that galaxies grow according to the evolution of the halo mass function, which we allow to vary with large-scale environment. Our model produces a reasonable match to the observed ultraviolet (UV) luminosity functions in this era by regulating star formation through stellar feedback and assuming that the UV luminosity function is dominated by recent star formation. We find that cosmic variance in the UV luminosity function is dominated by the variance in the underlying dark matter halo population, and not by differences in halo accretion or the specifics of our stellar feedback model. We also find that cosmic variance dominates over Poisson noise for future high-z surveys except for the brightest sources or at very high redshifts (z ≳ 12). We provide a linear approximation of cosmic variance for a variety of redshifts, magnitudes, and survey areas through the public python package galcv. Finally, we introduce a new method for incorporating priors on cosmic variance into estimates of the galaxy luminosity function and demonstrate that it significantly improves constraints on that important observable.more » « less
-
Abstract Globular clusters (GCs) provide valuable insight into the properties of their host galaxies’ dark matter halos. UsingN-body simulations incorporating semianalytic dynamical friction and GC−GC merger prescriptions, we study the evolution of GC radial distributions and mass functions in cuspy and cored dark matter halos. Modeling the dynamics of the GC-rich system in the dwarf galaxy UGC 7369, we find that friction-induced inspiral and subsequent mergers of massive GCs can naturally and robustly explain the mass segregation of the GCs and the existence of a nuclear star cluster (NSC). However, the multiple mergers required to form the NSC only take place when the dark matter halo is cuspy. In a cored halo, stalling of the dynamical friction within the core halts the inspiral of the GCs, and so the GC merger rate falls significantly, precluding the formation of an NSC. We therefore argue that the presence of an NSC requires a cusp in UGC 7369. More generally, we propose that the presence of an NSC and the corresponding alteration of the GC mass function due to mergers may be used as an indicator of a cuspy halo for galaxies in which we expect NSC formation to be merger dominated. These observables represent a simple, powerful complement to other inner halo density profile constraint techniques and should allow for straightforward extension to larger samples.more » « less
-
Abstract We present a method that calibrates a semianalytic model to the Renaissance Simulations, a suite of cosmological hydrodynamical simulations with high-redshift galaxy formation. This approach combines the strengths of semianalytic techniques and hydrodynamical simulations, enabling the extension to larger volumes and lower redshifts that are inaccessible to simulations due to computational expense. Using a sample of Renaissance star formation histories from an average density region of the Universe, we construct a four-parameter prescription for metal-enriched star formation characterized by an initial bursty stage followed by a steady stage where stars are formed at constant efficiencies. Our model also includes a treatment of Pop III star formation where a minimum halo mass and log-normal distribution of stellar mass are adopted to match the numerical simulations. Star formation is generally well reproduced for halos with masses ≲109M⊙. Between 11 <z< 25 our model produces metal-enriched star formation rate densities (SFRDs) that typically agree with Renaissance within a factor of ∼2 for the average density region. Additionally, the total metal-enriched stellar mass only differs from Renaissance by about 10% atz∼ 11. For regions that are either more overdense or rarefied and not included in the calibration, we produce metal-enriched SFRDs that agree with Renaissance within a factor of ∼2 at high-zbut eventually differ by higher factors for later times. This is likely due to environmental dependencies not included in the model. Our star formation prescriptions can easily be adopted in other analytic or semianalytic works to match our calibration to Renaissance.more » « less
-
ABSTRACT The mass assembly history (MAH) of dark matter haloes plays a crucial role in shaping the formation and evolution of galaxies. MAHs are used extensively in semi-analytic and empirical models of galaxy formation, yet current analytic methods to generate them are inaccurate and unable to capture their relationship with the halo internal structure and large-scale environment. This paper introduces florah (FLOw-based Recurrent model for Assembly Histories), a machine-learning framework for generating assembly histories of ensembles of dark matter haloes. We train florah on the assembly histories from the Gadget at Ultra-high Redshift with Extra Fine Time-steps and vsmdplN-body simulations and demonstrate its ability to recover key properties such as the time evolution of mass and concentration. We obtain similar results for the galaxy stellar mass versus halo mass relation and its residuals when we run the Santa Cruz semi-analytic model on florah-generated assembly histories and halo formation histories extracted from an N-body simulation. We further show that florah also reproduces the dependence of clustering on properties other than mass (assembly bias), which is not captured by other analytic methods. By combining multiple networks trained on a suite of simulations with different redshift ranges and mass resolutions, we are able to construct accurate main progenitor branches with a wide dynamic mass range from $z=0$ up to an ultra-high redshift $$z \approx 20$$, currently far beyond that of a single N-body simulation. florah is the first step towards a machine learning-based framework for planting full merger trees; this will enable the exploration of different galaxy formation scenarios with great computational efficiency at unprecedented accuracy.more » « less
