Renaissance Simulations Research Articles

Predicting the UV Escape Fraction of the First Galaxies in the Renaissance Simulations with Machine Learning

Cosmic reionization is likely driven by UV starlight emanating from the first generations of galaxies. A galaxy’s UV escape fraction, or the fraction of photons escaping from the galaxy, is useful to quantify its contribution to reionization. However, the UV escape fraction is notoriously difficult to predict due to local environment dependency and variability over time. Using data from the Renaissance Simulations, we attempt to make predictions about the impact of the first stars and galaxies on their environments. We present a time-independent classification model using a general artificial neural network architecture to predict the UV escape fraction given other galaxy properties—namely halo mass, stellar mass, redshift, star formation rate, lookback time, and gas fraction. We find our validation accuracy to be approximately 50%–65%, depending on the data set size from each zoom-in region of the Renaissance Simulations.

How long do high redshift massive black hole seeds remain outliers in black hole versus host galaxy relations?

ABSTRACT The existence of 109 M⊙ supermassive black holes (SMBHs) within the first billion years of the Universe remains a puzzle in our conventional understanding of black hole formation and growth. Several suggested formation pathways for these SMBHs lead to a heavy seed, with an initial black hole mass of 104–106 M⊙. This can lead to an overly massive BH galaxy (OMBG), whose nuclear black hole’s mass is comparable to or even greater than the surrounding stellar mass: the black hole to stellar mass ratio is Mbh/M* ≫ 10−3, well in excess of the typical values at lower redshift. We investigate how long these newborn BHs remain outliers in the Mbh − M* relation, by exploring the subsequent evolution of two OMBGs previously identified in the Renaissance simulations. We find that both OMBGs have Mbh/M* > 1 during their entire life, from their birth at z ≈ 15 until they merge with much more massive haloes at z ≈ 8. We find that the OMBGs are spatially resolvable from their more massive, 1011 M⊙, neighbouring haloes until their mergers are complete at z ≈ 8. This affords a window for future observations with JWST and sensitive X-ray telescopes to diagnose the heavy-seed scenario, by detecting similar OMBGs and establishing their uniquely high black hole-to-stellar mass ratio.

Massive Star Formation in Metal-Enriched Haloes at High Redshift

The formation of supermassive stars has generally been studied under the assumption of rapid accretion of pristine metal-free gas. Recently it was found, however, that gas enriched to metallicities up to $Z \sim 10^{-3}$ Z$_{\odot}$ can also facilitate supermassive star formation, as long as the total mass infall rate onto the protostar remains sufficiently high. We extend the analysis further by examining how the abundance of supermassive star candidate haloes would be affected if all haloes with super-critical infall rates, regardless of metallicity were included. We investigate this scenario by identifying all atomic cooling haloes in the Renaissance simulations with central mass infall rates exceeding a fixed threshold. We find that among these haloes with central mass infall rates above 0.1 M$_{\odot}$ yr$^{-1}$ approximately two-thirds of these haloes have metallicities of $Z > 10^{-3}$ Z$_{\odot}$. If metal mixing within these haloes is inefficient early in their assembly and pockets of metal-poor gas can remain then the number of haloes hosting supermassive stars can be increased by at least a factor of four. Additionally the centres of these high infall-rate haloes provide ideal environments in which to grow pre-existing black holes. Further research into the (supermassive) star formation dynamics of rapidly collapsing haloes, with inhomogeneous metal distributions, is required to gain more insight into both supermassive star formation in early galaxies as well as early black hole growth.

Fully Coupled Simulation of Cosmic Reionization. III. Stochastic Early Reionization by the Smallest Galaxies

Abstract Previously we identified a new class of early galaxy that we estimate contributes up to 30% of the ionizing photons responsible for reionization. These are low-mass halos in the range = 106.5–108 that have been chemically enriched by supernova ejecta from prior Population III star formation. Despite their low star formation rates, these metal cooling halos (MCs) are significant sources of ionizing radiation, especially at the onset of reionization, owing to their high number density and ionizing escape fractions. Here we present a fully coupled radiation hydrodynamic simulation of reionization that includes these MCs, as well the more massive hydrogen atomic line cooling halos. Our method is novel: we perform halo finding in line with the radiation hydrodynamical simulation and assign escaping ionizing fluxes to halos using a probability distribution function (pdf) measured from the galaxy-resolving Renaissance Simulations. The pdf captures the mass dependence of the ionizing escape fraction, as well as the probability that a halo is actively forming stars. With MCs, reionization starts earlier than if only halos of 108 and above are included; however, the redshift when reionization completes is only marginally affected, as this is driven by more massive galaxies. Because star formation is intermittent in MCs, the earliest phase of reionization exhibits a stochastic nature, with small H ii regions forming and recombining. Only later, once halos of mass ∼109 and above begin to dominate the ionizing emissivity, does reionization proceed smoothly in the usual manner deduced from previous studies. This occurs at z ≈ 10 in our simulation.

The growth of black holes from Population III remnants in the Renaissance simulations

The formation of stellar mass black holes from the remnants of Population III stars provides a source of initial black hole seeds with the potential to grow into intermediate or, in rare cases, possibly supermassive black holes. We use the Renaissance simulation suite to follow the growth of over 15,000 black holes born into mini-haloes in the early Universe. We compute the evolution of the black holes by post-processing individual remnant Population III star particles in the Renaissance simulation snapshots. The black holes populate haloes from 10$^{6}$ M$_{\odot}$ up to 10$^{9}$ M$_{\odot}$. We find that all of the black holes display very inefficient growth. On average the black holes increase their initial mass by a factor 10$^{-5}$, with the most active black holes increasing their mass by approximately 10%. Only a single black hole experiences any period of super-Eddington accretion, but the duration is very short and not repeated. Furthermore, we find no correlation of black hole accretion with halo mass in the mass range sampled. Within most haloes, we identify clumps of cool, dense gas for which accretion rates would be high, but instances of black holes encountering these clumps are rare and short-lived. Star formation competes with black hole growth by consuming available gas and driving down accretion rates through feedback. We conclude that the black holes born from Population III remnants do not form a significant population of intermediate mass black holes in the early Universe and will need to wait until later times to undergo significant accretion, if at all.

First light – II. Emission line extinction, population III stars, and X-ray binaries

We produce synthetic spectra and observations for metal-free stellar populations and high mass X-ray binaries in the Renaissance Simulations at a redshift of 15. We extend our methodology from the first paper in the series by modelling the production and extinction of emission lines throughout a dusty and metal-enriched interstellar and circum-galactic media extracted from the simulation, using a Monte Carlo calculation. To capture the impact of high-energy photons, we include all frequencies from hard X-ray to far infrared with enough frequency resolution to discern line emission and absorption profiles. The most common lines in our sample in order of their rate of occurrence are Ly$\alpha$, the C IV $\lambda\lambda1548,1551$ doublet, H-$\alpha$, and the Ca II $\lambda\lambda\lambda8498,8542,8662$ triplet. The best scenario for a direct observation of a metal-free stellar population is a merger between two Population III galaxies. In mergers between metal-enriched and metal-free stellar populations, some characteristics may be inferred indirectly. Single Population III galaxies are too dim to be observed photometrically at $z = 15$. Ly$\alpha$ emission is discernible by $JWST$ as an increase in $\rm{J_{200w} - J_{277w}}$ colour off the intrinsic stellar tracks. Observations of metal-free stars will be difficult, though not impossible, with the next generation of space telescopes.

First light: exploring the spectra of high-redshift galaxies in the Renaissance Simulations

We present synthetic observations for the first generations of galaxies in the Universe and make predictions for future deep field observations for redshifts greater than 6. Due to the strong impact of nebular emission lines and the relatively compact scale of HII regions, high resolution cosmological simulations and a robust suite of analysis tools are required to properly simulate spectra. We created a software pipeline consisting of FSPS, Hyperion, Cloudy and our own tools to generate synthetic IR observations from a fully three-dimensional arrangement of gas, dust, and stars. Our prescription allows us to include emission lines for a complete chemical network and tackle the effect of dust extinction and scattering in the various lines of sight. We provide spectra, 2-D binned photon imagery for both HST and JWST IR filters, luminosity relationships, and emission line strengths for a large sample of high redshift galaxies in the Renaissance Simulations. Our resulting synthetic spectra show high variability between galactic halos with a strong dependence on stellar mass, metallicity, gas mass fraction, and formation history. Halos with the lowest stellar mass have the greatest variability in [OIII]/H$\beta$, [OIII] and CIII] while halos with higher masses are seen to show consistency in their spectra and [OIII] equivalent widths (EW) between 1\AA\ and 10\AA. Viewing angle accounted for three-fold difference in flux due to the presence of ionized gas channels in a halo. Furthermore, JWST color plots show a discernible relationship between redshift, color, and mean stellar age.

GALAXY PROPERTIES AND UV ESCAPE FRACTIONS DURING THE EPOCH OF REIONIZATION: RESULTS FROM THE RENAISSANCE SIMULATIONS

ABSTRACT Cosmic reionization is thought to be primarily fueled by the first generations of galaxies. We examine their stellar and gaseous properties, focusing on the star formation rates and the escape of ionizing photons, as a function of halo mass, redshift, and environment using the full suite of the Renaissance Simulations with an eye to provide better inputs to global reionization simulations. This suite probes overdense, average, and underdense regions of the universe of several hundred comoving Mpc3, each yielding a sample of over 3000 halos in the mass range of 107–109.5 at their final redshifts of 15, 12.5, and 8, respectively. In the process, we simulate the effects of radiative and supernova feedback from 5000 to 10,000 Population III stars in each simulation. We find that halos as small as 107 are able to host bursty star formation due to metal-line cooling from earlier enrichment by massive Population III stars. Using our large sample, we find that the galaxy-halo occupation fraction drops from unity at virial masses above 108.5 to ∼50% at 108 and ∼10% at 107 , quite independent of redshift and region. Their average ionizing escape fraction is ∼5% in the mass range of 108–109 and increases with decreasing halo mass below this range, reaching 40%–60% at 107 . Interestingly, we find that the escape fraction varies between 10%–20% in halos with virial masses of ∼3 × 109 . Taken together, our results confirm the importance of the smallest galaxies as sources of ionizing radiation contributing to the reionization of the universe.

X-RAY BACKGROUND AT HIGH REDSHIFTS FROM POP III REMNANTS: RESULTS FROM POP III STAR FORMATION RATES IN THE RENAISSANCE SIMULATIONS

Due to their long mean free paths, X-rays are expected to have many significant impacts globally on the properties of the intergalactic medium (IGM) by their heating and ionizing processes on large scales. At high redshifts, X-rays from Population (Pop) III binaries might have important effects on cosmic reionization and the Lyman alpha forest. As a continuation of our previous work on Pop III binary X-rays (Xu et al. 2014), we use the Pop III distribution and evolution from the Renaissance Simulations, a suite of self-consistent cosmological radiation hydrodynamics simulations of the formation of the first galaxies, to calculate the X-ray luminosity density and background over the redshift range 20 > z > 7.6. As we find that Pop III star formation continues at a low, nearly constant rate to the end of reionization, X-rays are being continuously produced at significant rates compared to other possible X-ray sources, such as AGNs and normal X-ray binaries during the same period of time. We estimate that Pop III binaries produce approximately 6 eV of energy in the X-rays per hydrogen atom. We calculate the X-ray background for different monochromatic photon energies. KeV X-rays redshift and accumulate to produce a strong X-ray background spectrum extending to roughly 500 eV. The X-ray background is strong enough to heat the IGM to ~ 1000 K and to ionize a few percent of the neutral hydrogen. These effects are important for an understanding of the neutral hydrogen hyperfine transition 21-cm line signatures, the Ly alpha forest, and optical depth of the CMB to Thomson scattering.

LATE POP III STAR FORMATION DURING THE EPOCH OF REIONIZATION: RESULTS FROM THE RENAISSANCE SIMULATIONS

ABSTRACT We present results on the formation of Population III (Pop III) stars at redshift 7.6 from the Renaissance Simulations, a suite of extremely high-resolution and physics-rich radiation transport hydrodynamics cosmological adaptive-mesh refinement simulations of high-redshift galaxy formation performed on the Blue Waters supercomputer. In a survey volume of about 220 comoving Mpc3, we found 14 Pop III galaxies with recent star formation. The surprisingly late formation of Pop III stars is possible due to two factors: (i) the metal enrichment process is local and slow, leaving plenty of pristine gas to exist in the vast volume; and (ii) strong Lyman–Werner radiation from vigorous metal-enriched star formation in early galaxies suppresses Pop III formation in (“not so”) small primordial halos with mass less than ∼3 × 107 M ⊙. We quantify the properties of these Pop III galaxies and their Pop III star formation environments. We look for analogs to the recently discovered luminous Ly α emitter CR7, which has been interpreted as a Pop III star cluster within or near a metal-enriched star-forming galaxy. We find and discuss a system similar to this in some respects, however, the Pop III star cluster is far less massive and luminous than CR7 is inferred to be.

PROBING THE ULTRAVIOLET LUMINOSITY FUNCTION OF THE EARLIEST GALAXIES WITH THE RENAISSANCE SIMULATIONS

In this paper, we present the first results from the Renaissance Simulations, a suite of extremely high-resolution and physics-rich AMR calculations of high-redshift galaxy formation performed on the Blue Waters supercomputer. These simulations contain hundreds of well-resolved galaxies at , and make several novel, testable predictions. Most critically, we show that the ultraviolet luminosity function of our simulated galaxies is consistent with observations of high-z galaxy populations at the bright end of the luminosity function , but at lower luminosities is essentially flat rather than rising steeply, as has been inferred by Schechter function fits to high-z observations, and has a clearly defined lower limit in UV luminosity. This behavior of the luminosity function is due to two factors: (i) the strong dependence of the star formation rate (SFR) on halo virial mass in our simulated galaxy population, with lower-mass halos having systematically lower SFRs and thus lower UV luminosities; and (ii) the fact that halos with virial masses below do not universally contain stars, with the fraction of halos containing stars dropping to zero at . Finally, we show that the brightest of our simulated galaxies may be visible to current and future ultra-deep space-based surveys, particularly if lensed regions are chosen for observation.