We introduce the Phoenix Simulations, a suite of highly resolved cosmological simulations featuring hydrodynamics, primordial gas chemistry, primordial and enriched star formation and feedback, UV radiative transfer, and saved outputs with Δkpc. 7 and 17% of regions have
The recently launched James Webb Space Telescope promises unparalleled advances in our understanding of the first stars and galaxies, but realizing this potential requires cosmological simulations that capture the key physical processes that affected these objects. Here, we show that radiative transfer and subgrid turbulent mixing are two such processes. By comparing simulations with and without radiative transfer but with exactly the same physical parameters and subgrid turbulent mixing model, we show that tracking radiative transfer suppresses the Population III star formation density by a factor ≈4. In both simulations, ≳90% of Population III stars are found in the unresolved pristine regions tracked by our subgrid model, which does a better job at modeling the regions surrounding proto-galaxy cores where metals from supernovae take tens of megayears to mix thoroughly. At the same time, radiative transfer suppresses Population III star formation, via the development of ionized bubbles that slow gas accretion in these regions, and it results in compact high-redshift galaxies that are surrounded by isolated low-mass satellites. Thus, turbulent mixing and radiative transfer are both essential processes that must be included to accurately model the morphology, composition, and growth of primordial galaxies.
more » « less- Award ID(s):
- 1715876
- PAR ID:
- 10370055
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 935
- Issue:
- 2
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 174
- Size(s):
- Article No. 174
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract t = 200 kyr. We observe 73,523 individual primordial stars within 3313 distinct regions forming 2110 second-generation enriched star clusters byz ≥ 12 within a combined 177.25 Mpc3volume across three simulations. The regions that lead to enriched star formation can contain ≳150 primordial stars, with 80% of regions having experienced combinations of primordial Type II, hypernovae, and/or pair-instability supernovae. Primordial supernovae enriched 0.8% of the volume, with 2% of enriched gas enriched by later-generation stars. We determine the extent of a primordial stellar region by its metal-rich or ionized hydrogen surrounding cloud; the metal-rich and ionized regions have time-dependent average radiir ≲ 3r > 7 kpc for metal-rich and ionized radii, respectively. We find that the metallicity distribution function of second-generation stars overlaps that of subsequent Population II star formation, spanning metal-deficient (∼7.94 × 10−8Z ⊙) to supersolar (∼3.71Z ⊙), and that 30.5% of second-generation stars haveZ > 10−2Z ⊙. We find that the metallicity of second-generation stars depends on progenitor configuration, with metals from pair-instability supernovae contributing to the most metal-rich clusters; these clusters form promptly after the supernova event. Finally, we create an interpretable regression model to predict the radius of the metal-rich influence of Population III star systems within the first 7–18 Myr after the first Population III stars form in the region. -
more » « less
ABSTRACT Galaxies comprise intricate networks of interdependent processes which together govern their evolution. Central among these are the multiplicity of feedback channels, which remain incompletely understood. One outstanding problem is the understanding and modelling of the multiphase nature of galactic winds, which play a crucial role in galaxy formation and evolution. We present the results of three-dimensional magnetohydrodynamical simulations of tall–box interstellar medium (ISM) patches with clustered supernova-driven outflows. Dynamical fragmentation of the ISM during superbubble breakout seeds the resulting hot outflow with a population of cool clouds. We focus on analyzing and modelling the origin and properties of these clouds. Their presence induces large-scale turbulence, which, in turn, leads to complex cloud morphologies. Cloud sizes are well described by a power-law distribution and mass growth rates can be modelled using turbulent radiative mixing layer theory. Turbulence provides significant pressure support in the clouds, while magnetic fields only play a minor role. We conclude that many of the physical insights and analytic scalings derived from idealized small-scale simulations of turbulent radiative mixing layers and cloud–wind interactions are directly translatable and applicable to these larger scale cloud populations. This opens the door to developing effective subgrid recipes for their inclusion in global-scale galaxy models where they are unresolved.
-
ABSTRACT The cosmic near-infrared background (NIRB) offers a powerful integral probe of radiative processes at different cosmic epochs, including the pre-reionization era when metal-free, Population III (Pop III) stars first formed. While the radiation from metal-enriched, Population II (Pop II) stars likely dominates the contribution to the observed NIRB from the reionization era, Pop III stars – if formed efficiently – might leave characteristic imprints on the NIRB, thanks to their strong Lyα emission. Using a physically motivated model of first star formation, we provide an analysis of the NIRB mean spectrum and anisotropy contributed by stellar populations at z > 5. We find that in circumstances where massive Pop III stars persistently form in molecular cooling haloes at a rate of a few times $10^{-3}\, \mathrm{ M}_\odot \ \mathrm{yr}^{-1}$, before being suppressed towards the epoch of reionization (EoR) by the accumulated Lyman–Werner background, a unique spectral signature shows up redward of $1\, \mu$m in the observed NIRB spectrum sourced by galaxies at z > 5. While the detailed shape and amplitude of the spectral signature depend on various factors including the star formation histories, initial mass function, LyC escape fraction and so forth, the most interesting scenarios with efficient Pop III star formation are within the reach of forthcoming facilities, such as the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer. As a result, new constraints on the abundance and formation history of Pop III stars at high redshifts will be available through precise measurements of the NIRB in the next few years.more » « less
-
null (Ed.)The spatial decorrelation of dense molecular gas and young stars observed on ≲ 1 kiloparsec scales in nearby galaxies indicates rapid dispersal of star-forming regions by stellar feedback. We explore the sensitivity of this decorrelation to different processes controlling the structure of the interstellar medium, the abundance of molecular gas, star formation, and feedback in a suite of simulations of an isolated dwarf galaxy with structural properties similar to NGC300 that self-consistently model radiative transfer and molecular chemistry. Our fiducial simulation reproduces the magnitude of decorrelation and its scale dependence measured in NGC300, and we show that this agreement is due to different aspects of feedback, including H2 dissociation, gas heating by the locally variable UV field, early mechanical feedback, and supernovae. In particular, early radiative and mechanical feedback affect the correlation on ≲100 pc scales, while supernovae play a significant role on ≳100 pc scales. The correlation is also sensitive to the choice of the local star formation efficiency per freefall time, eps_ff, which provides a strong observational constraint on eps_ff when the global star formation rate is independent of its value. Finally, we explicitly show that the degree of correlation between the peaks of molecular gas and star formation density is directly related to the distribution of the lifetimes of star-forming regions.more » « less
-
more » « less
Abstract Understanding the connections between galaxy stellar mass, star formation rate, and dark matter halo mass represents a key goal of the theory of galaxy formation. Cosmological simulations that include hydrodynamics, physical treatments of star formation, feedback from supernovae, and the radiative transfer of ionizing photons can capture the processes relevant for establishing these connections. The complexity of these physics can prove difficult to disentangle and obfuscate how mass-dependent trends in the galaxy population originate. Here, we train a machine-learning method called Explainable Boosting Machines (EBMs) to infer how the stellar mass and star formation rate of nearly 6 million galaxies simulated by the Cosmic Reionization on Computers project depend on the physical properties of halo mass, the peak circular velocity of the galaxy during its formation history
v peak, cosmic environment, and redshift. The resulting EBM models reveal the relative importance of these properties in setting galaxy stellar mass and star formation rate, withv peakproviding the most dominant contribution. Environmental properties provide substantial improvements for modeling the stellar mass and star formation rate in only ≲10% of the simulated galaxies. We also provide alternative formulations of EBM models that enable low-resolution simulations, which cannot track the interior structure of dark matter halos, to predict the stellar mass and star formation rate of galaxies computed by high-resolution simulations with detailed baryonic physics.
