We present the SPHINX suite of cosmological adaptive mesh refinement simulations, the first radiation-hydrodynamical simulations to simultaneously capture large-scale reionization and the escape of ...ionizing radiation from thousands of resolved galaxies. Our 5 and 10 co-moving Mpc volumes resolve haloes down to the atomic cooling limit and model the interstellar medium with better than ≈10 pc resolution. The project has numerous goals in improving our understanding of reionization and making predictions for future observations. In this first paper, we study how the inclusion of binary stars in computing stellar luminosities impacts reionization, compared to a model that includes only single stars. Owing to the suppression of galaxy growth via strong feedback, our galaxies are in good agreement with observational estimates of the galaxy luminosity function. We find that binaries have a significant impact on the timing of reionization: with binaries, our boxes are 99.9 percent ionized by volume at z ≈ 7, while without them our volumes fail to reionize by z = 6. These results are robust to changes in volume size, resolution, and feedback efficiency. The escape of ionizing radiation from individual galaxies varies strongly and frequently. On average, binaries lead to escape fractions of ≈7-10 percent, about three times higher than with single stars only. The higher escape fraction is a result of a shallower decline in ionizing luminosity with age, and is the primary reason for earlier reionization, although the higher integrated luminosity with binaries also plays a subdominant role.
Abstract
To understand the impact of radiation feedback during the formation of a globular cluster (GC), we simulate a head-on collision of two turbulent giant molecular clouds (GMCs). A series of ...idealized radiation-hydrodynamic simulations is performed, with and without stellar radiation or Type II supernovae. We find that a gravitationally bound, compact star cluster of mass
M
GC
∼ 10
5
M
⊙
forms within ≈3 Myr when two GMCs with mass
M
GMC
= 3.6 × 10
5
M
⊙
collide. The GC candidate does not form during a single collapsing event but emerges due to the mergers of local dense gas clumps and gas accretion. The momentum transfer due to the absorption of the ionizing radiation is the dominant feedback process that suppresses the gas collapse, and photoionization becomes efficient once a sufficient number of stars form. The cluster mass is larger by a factor of ∼2 when the radiation feedback is neglected, and the difference is slightly more pronounced (16%) when extreme Ly
α
feedback is considered in the fiducial run. In the simulations with radiation feedback, supernovae explode after the star-forming clouds are dispersed, and their metal ejecta are not instantaneously recycled to form stars.
ABSTRACT We investigate the formation of metal-poor globular clusters (GCs) at the center of two dark matter halos with at using cosmological radiation-hydrodynamics simulations. We find that very ...compact ( 1 pc) and massive ( ) clusters form rapidly when pristine gas collapses isothermally with the aid of efficient Ly emission during the transition from molecular-cooling halos to atomic-cooling halos. Because the local free-fall time of dense star-forming gas is very short ( ), a large fraction of the collapsed gas is turned into stars before stellar feedback processes blow out the gas and shut down star formation. Although the early stage of star formation is limited to a small region of the central star-forming disk, we find that the disk quickly fragments due to metal enrichment from supernovae. Sub-clusters formed in the fragmented clouds eventually merge with the main cluster at the center. The simulated clusters closely resemble the local GCs in mass and size but show a metallicity spread that is much wider than found in the local GCs. We discuss a role of pre-enrichment by Pop III and II stars as a potential solution to the latter issue. Although not without shortcomings, it is encouraging that a naive blind (not tuned) cosmological simulation presents a possible channel for the formation of at least some massive GCs.
Supermassive black holes (BH) are powerful sources of energy that are already in place at very early epochs of the Universe (by z = 6). Using hydrodynamical simulations of the formation of a massive ...M
vir = 5 × 1011 M halo by z = 6 (the most massive progenitor of a cluster of M
vir = 2 × 1015 M at z = 0), we evaluate the impact of active galactic nuclei (AGN) on galaxy mass content, BH self-regulation and gas distribution inside this massive halo. We find that supernova feedback has a marginal influence on the stellar structure, and no influence on the mass distribution on large scales. In contrast, AGN feedback alone is able to significantly alter the stellar-bulge mass content by quenching star formation when the BH is self-regulating, and by depleting the cold gas reservoir in the centre of the galaxy. The growth of the BH proceeds first by a rapid Eddington-limited period fed by direct cold filamentary infall. When the energy delivered by the AGN is sufficiently large to unbind the cold gas of the bulge, the accretion of gas on to the BH is maintained by both smooth gas inflow and clump migration through the galactic disc triggered by merger-induced torques. The feedback from the AGN has also a severe consequence on the baryon mass content within the halo, producing large-scale hot superwinds, able to blow away some of the cold filamentary material from the centre and reduce the baryon fraction by more than 30 per cent within the halo's virial radius. Thus, in the very young universe, AGN feedback is likely to be a key process, shaping the properties of the most massive galaxies.
ABSTRACT
We measure escape fractions, fesc, of ionizing radiation from galaxies in the sphinx suite of cosmological radiation-hydrodynamical simulations of reionization, resolving haloes with $M_{\rm ...vir}\gtrsim 7.5 \times 10^7 \ {\rm {M}_{\odot }}$ with a minimum cell width of ≈10 pc. Our new and largest 20 co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with halo masses up to a few times 1011 M⊙. The simulated galaxies agree well with observational constraints of the ultraviolet (UV) luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over time-scales of a few Myr, due to its regulation by supernova and radiation feedback, and at any given time a tiny fraction of star-forming galaxies emits a large fraction of the ionizing radiation escaping into the intergalactic medium. Statistically, fesc peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies (M* ≈ 107 M⊙, M1500 ≈ −17, Z ≲ 5 × 10−3 Z⊙), dropping strongly for lower and higher masses, brighter and dimmer galaxies, and more metal-rich galaxies. The escape fraction correlates positively with both the short-term and long-term specific star formation rate. According to sphinx, galaxies too dim to be yet observed, with ${M_{1500}}\gtrsim -17$, provide about 55 per cent of the photons contributing to reionization. The global averaged fesc naturally decreases with decreasing redshift, as predicted by UV background models and low-redshift observations. This evolution is driven by decreasing specific star formation rates over cosmic time.
Abstract
Understanding the origin of strong galactic outflows and the suppression of star formation in dwarf galaxies is a key problem in galaxy formation. Using a set of radiation-hydrodynamic ...simulations of an isolated dwarf galaxy embedded in a 1010 M⊙ halo, we show that the momentum transferred from resonantly scattered Lyman-α (Lyα) photons is an important source of stellar feedback which can shape the evolution of galaxies. We find that Lyα feedback suppresses star formation by a factor of two in metal-poor galaxies by regulating the dynamics of star-forming clouds before the onset of supernova explosions (SNe). This is possible because each Lyα photon resonantly scatters and imparts ∼10–300 times greater momentum than in the single scattering limit. Consequently, the number of star clusters predicted in the simulations is reduced by a factor of ∼5, compared to the model without the early feedback. More importantly, we find that galactic outflows become weaker in the presence of strong Lyα radiation feedback, as star formation and associated SNe become less bursty. We also examine a model in which radiation field is arbitrarily enhanced by a factor of up to 10, and reach the same conclusion. The typical mass-loading factors in our metal-poor dwarf system are estimated to be ∼5–10 near the mid-plane, while it is reduced to ∼1 at larger radii. Finally, we find that the escape of ionizing radiation and hence the reionization history of the Universe is unlikely to be strongly affected by Lyα feedback.
We present the O
BELISK
project, a cosmological radiation-hydrodynamics simulation that follows the assembly and reionization of a protocluster progenitor during the first two billion years after the ...big bang, down to
z
= 3.5. The simulation resolves haloes down to the atomic cooling limit and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The O
BELISK
project is specifically designed to study the coevolution of high-redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high-redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to those of high-redshift protoclusters. We find that hydrogen reionization happens inside-out, is completed by
z
∼ 6 in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at
z
∼ 4.
ABSTRACT
We present the first results from SPHINX-MHD, a suite of cosmological radiation-magnetohydrodynamics simulations designed to study the impact of primordial magnetic fields (PMFs) on galaxy ...formation and the evolution of the intergalactic medium (IGM) during the epoch of reionization. The simulations are among the first to employ multifrequency, on-the-fly radiation transfer and constrained transport ideal MHD in a cosmological context to simultaneously model the inhomogeneous process of reionization as well as the growth of primordial magnetic fields. We run a series of $(5\, \text{cMpc})^3$ cosmological volumes, varying both the strength of the seed magnetic field and its spectral index. We find that PMFs with a spectral index (nB) and a comoving amplitude (B0) that have $\scriptstyle{n_B \gt -0.562\log _{10}\left(\frac{B_0}{1{\rm n}G}\right) - 3.35}$ produce electron optical depths (τe) that are inconsistent with CMB constraints due to the unrealistically early collapse of low-mass dwarf galaxies. For nB ≥ −2.9, our constraints are considerably tighter than the ∼nG constraints from Planck. PMFs that do not satisfy our constraints have little impact on the reionization history or the shape of the UV luminosity function. Likewise, detecting changes in the Ly α forest due to PMFs will be challenging because photoionization and photoheating efficiently smooth the density field. However, we find that the first absorption feature in the global 21-cm signal is a particularly sensitive indicator of the properties of the PMFs, even for those that satisfy our τe constraint. Furthermore, strong PMFs can marginally increase the escape of LyC photons by up to 25 per cent and shrink the effective radii of galaxies by $\sim 44{{\ \rm per\ cent}}$ which could increase the completeness fraction of galaxy surveys. Finally, our simulations show that surveys with a magnitude limit of ${\rm \mathit{ M}_{UV,1500\,\mathring{\rm A}}=-13}$ can probe the sources that provide the majority of photons for reionization out to z = 12.
The environmental effect is commonly used to explain the excess of gas-poor galaxies in galaxy clusters. Meanwhile, the presence of gas-poor galaxies at cluster outskirts, where galaxies have not ...spent enough time to feel the cluster environmental effect, hints at the presence of preprocessing. Using cosmological hydrodynamic simulations on 16 clusters, we investigate the mechanisms of gas depletion of galaxies found inside clusters. The gas-depletion mechanisms can be categorized into three channels based on where and when they took place. First, 34% of our galaxies are gas poor before entering clusters ("preprocessing"). They are mainly satellites that have undergone the environmental effect inside group halos. Second, 43% of the sample quickly became gas deficient in clusters before the first pericentric pass ("fast cluster processing"). Some of them were group satellites that are low in gas at the time of cluster entry compared to the galaxies directly coming from the field. Even the galaxies with large gas fractions take this channel if they fall into massive clusters ( 1014.5 M ) or approach cluster centers through radial orbits. Third, 24% of our sample retain gas even after their first pericentric pass ("slow cluster processing") as they fall into the less massive clusters or have circular orbits. The relative importance of each channel varies with a cluster's mass, while the exact degree of significance is subject to large uncertainties. Group preprocessing accounts for one-third of the total gas depletion, but it also determines the gas fraction of galaxies at their cluster entry, which in turn determines whether a galaxy should take the fast or slow cluster processing.