SCORCH (Simulations and Constructions of the Reionization of Cosmic Hydrogen) is a new project to study the Epoch of Reionization (EoR). In this first paper, we probe the connection between observed high-redshift galaxies and simulated dark matter ha
los to better understand the abundance and evolution of the primary source of ionizing radiation. High-resolution N-body simulations are run to quantify the abundance of dark matter halos as a function of mass $M$, accretion rate $dot{M}$, and redshift $z$. A new fit for the halo mass function $dn/dM$ is $approx 20%$ more accurate at the high-mass end where bright galaxies are expected to reside. A novel approach is used to fit the halo accretion rate function $dn/ddot{M}$ in terms of the halo mass function. Abundance matching against the observed galaxy luminosity function is used to estimate the luminosity-mass relation and the luminosity-accretion-rate relation. The inferred star formation efficiency is not monotonic with $M$ nor $dot{M}$, but reaches a maximum value at a characteristic mass $sim 2 times 10^{11} M_odot$ and a characteristic accretion rate $sim 6 times 10^2 M_odot/{rm yr}$ at $z approx 6$. We find a universal EoR luminosity-accretion-rate relation and construct a fiducial model for the galaxy luminosity function. The Schechter parameters evolve such that $phi_star$ decreases, $M_star$ is more positive (fainter), and $alpha$ is more negative (steeper) at higher redshifts. We forecast for the upcoming James Webb Space Telescope and show that with apparent magnitude limit $m_{rm AB} approx 31 (32)$, it can observe $gtrsim 11 (24)$ unlensed galaxies per square degree per unit redshift at least down to $M_star$ at $z lesssim 13 (14)$.
We reason that, without physical fine-tuning, neither the supermassive black holes (SMBHs) nor the stellar bulges can self-regulate or inter-regulate by driving away already fallen cold gas to produce the observed correlation between them. We suggest
an alternative scenario where the observed mass ratios of the SMBHs to bulges reflect the angular momentum distribution of infallen gas such that the mass reaching the stable accretion disc is a small fraction of that reaching the bulge region, averaged over the cosmological time scales. We test this scenario using high resolution, large-scale cosmological hydrodynamic simulations (without AGN feedback), assuming the angular momentum distribution of gas landing in the bulge region to yield a Mestel disc that is supported by independent simulations resolving the Bondi radii of SMBHs. A mass ratio of $0.1-0.3%$ between the very low angular momentum gas that free-falls to the sub-parsec region to accrete to the SMBH and the overall star formation rate is found. This ratio is found to increase with increasing redshift to within a factor of $sim 2$, suggesting that the SMBH to bulge ratio is nearly redshift independent, with a modest increase with redshift, a testable prediction. Furthermore, the duty cycle of active galactic nuclei (AGN) with high Eddington ratios is expected to increase significantly with redshift. Finally, while SMBHs and bulges are found to coevolve on $sim 30-150$Myr time scales or longer, there is indication that, on shorer time scales, the SMBH accretion rate and star formation may be less correlated.
Simulations have indicated that most of the escaped Lyman continuum photons escape through a minority of solid angles with near complete transparency, with the remaining majority of the solid angles largely opaque, resulting in a very broad and skewe
d probability distribution function (PDF) of the escape fraction when viewed at different angles. Thus, the escape fraction of Lyman continuum photons of a galaxy observed along a line of sight merely represents the properties of the interstellar medium along that line of sight, which may be an ill-representation of true escape fraction of the galaxy averaged over its full sky. Here we study how Lyman continuum photons escape from galaxies at $z=4-6$, utilizing high-resolution large-scale cosmological radiation-hydrodynamic simulations. We compute the PDF of the mean escape fraction ($left<f_{rm esc,1D}right>$) averaged over mock observational samples, as a function of the sample size, compared to the true mean (had you an infinite sample size). We find that, when the sample size is small, the apparent mean skews to the low end. For example, for a true mean of 6.7%, an observational sample of (2,10,50) galaxies at $z=4$ would have have 2.5% probability of obtaining the sample mean lower than $left<f_{rm esc,1D}right>=$(0.007%, 1.8%, 4.1%) and 2.5% probability of obtaining the sample mean being greater than (43%, 18%, 11%). Our simulations suggest that at least $sim$ 100 galaxies should be stacked in order to constrain the true escape fraction within 20% uncertainty.
Ly$alpha$ photons scattered by neutral hydrogen atoms in the circumgalactic media or produced in the halos of star-forming galaxies are expected to lead to extended Ly$alpha$ emission around galaxies. Such low surface brightness Ly$alpha$ halos (LAHs
) have been detected by stacking Ly$alpha$ images of high-redshift star-forming galaxies. We study the origin of LAHs by performing radiative transfer modeling of nine $z=3.1$ Lyman-Alpha Emitters (LAEs) in a high resolution hydrodynamic cosmological galaxy formation simulation. We develop a method of computing the mean Ly$alpha$ surface brightness profile of each LAE by effectively integrating over many different observing directions. Without adjusting any parameters, our model yields an average Ly$alpha$ surface brightness profile in remarkable agreement with observations. We find that observed LAHs cannot be accounted for solely by photons originating from the central LAE and scattered to large radii by hydrogen atoms in the circumgalactic gas. Instead, Ly$alpha$ emission from regions in the outer halo is primarily responsible for producing the extended LAHs seen in observations, which potentially includes both star-forming and cooling radiation. With the limit on the star formation contribution set by the ultra-violet (UV) halo measurement, we find that cooling radiation can play an important role in forming the extended LAHs. We discuss the implications and caveats of such a picture.
Recent large surveys have found a reversal of the star formation rate (SFR)-density relation at z=1 from that at z=0 (e.g. Elbaz et al.; Cooper et al.), while the sign of the slope of the color-density relation remains unchanged (e.g. Cucciati et al.
; Quadri et al.). We use state-of-the-art adaptive mesh refinement cosmological hydrodynamic simulations of a 21x24x20 (Mpc/h)$^3$ region centered on a cluster to examine the SFR-density and color-density relations of galaxies at z=0 and z=1. The local environmental density is defined by the dark matter mass in spheres of radius 1 Mpc/h, and we probe two decades of environmental densities. Our simulations produce a large increase of SFR with density at z=1, as in the observations of Elbaz et al. We also find a significant evolution to z=0, where the SFR-density relation is much flatter. The color-density relation in our simulations is consistent from z=1 to z=0, in agreement with observations. We find that the increase in the median SFR with local density at z=1 is due to a growing population of star-forming galaxies in higher-density environments. At z=0 and z=1 both the SFR and cold gas mass are tightly correlated with the galaxy halo mass, and therefore the correlation between median halo mass and local density is an important cause of the SFR-density relation at both redshifts. We also show that the local density on 1 Mpc/h scales affects galaxy SFRs as much as halo mass at z=0. Finally, we find indications that the role of the 1 Mpc/h scale environment reverses from z=0 to z=1: at z=0 high-density environments depress galaxy SFRs, while at z=1 high-density environments tend to increase SFRs.
Halo gas in low-z (z<0.5) >0.1L* galaxies in high-resolution, large-scale cosmological hydrodynamic simulations is examined with respect to three components: (cold, warm, hot) with temperatures equal to (<10^5, 10^{5-6}, >10^6)K, respectively. The wa
rm component is compared, utilizing O VI lambdalambda 1032, 1038 absorption lines, to observations and agreement is found with respect to the galaxy-O VI line correlation, the ratio of O VI line incidence rate in blue to red galaxies and the amount of O VI mass in star-forming galaxies. A detailed account of the sources of warm halo gas (stellar feedback heating, gravitational shock heating and accretion from the intergalactic medium), inflowing and outflowing warm halo gas metallicity disparities and their dependencies on galaxy types and environment is also presented. Having the warm component securely anchored, our simulations make the following additional predictions. First, cold gas is the primary component in inner regions, with its mass comprising 50% of all gas within galacto-centric radius r=(30,150)kpc in (red, blue) galaxies. Second, at r>(30,200)kpc in (red, blue) galaxies the hot component becomes the majority. Third, the warm component is a perpetual minority, with its contribution peaking at ~30% at r=100-300kpc in blue galaxies and never exceeding 5% in red galaxies. The significant amount of cold gas in low-z early-type galaxies found in simulations, in agreement with recent observations (Thom et al.), is intriguing, so is the dominance of hot gas at large radii in blue galaxies.
Lyman-alpha (Lya) photons that escape the interstellar medium of star-forming galaxies may be resonantly scattered by neutral hydrogen atoms in the circumgalactic and intergalactic media, thereby increasing the angular extent of the galaxys Lya emiss
ion. We present predictions of this extended, low surface brightness Lya emission based on radiative transfer modeling in a cosmological reionization simulation. The extended emission can be detected from stacked narrowband images of Lya emitters (LAEs) or of Lyman break galaxies (LBGs). Its average surface brightness profile has a central cusp, then flattens to an approximate plateau beginning at an inner characteristic scale below ~0.2 Mpc (comoving), then steepens again beyond an outer characteristic scale of ~1 Mpc. The inner scale marks the transition from scattered light of the central source to emission from clustered sources, while the outer scale marks the spatial extent of scattered emission from these clustered sources. Both scales tend to increase with halo mass, UV luminosity, and observed Lya luminosity. The extended emission predicted by our simulation is already within reach of deep narrowband photometry using large ground-based telescopes. Such observations would test radiative transfer models of emission from LAEs and LBGs, and they would open a new window on the circumgalactic environment of high-redshift star-forming galaxies.
In the context stellar reionization in the standard cold dark matter model, we analyze observations at z~6 and are able to draw three significant conclusions with respect to star formation and the state of the intergalactic medium (IGM) at z~6. (1) A
n initial stellar mass function (IMF) more efficient, by a factor of 10-20, in producing ionizing photons than the standard Salpeter IMF is required at z~6. This may be achieved by having either (A) a metal-enriched IMF with and a lower mass cutoff of >= 30Msun or (B) 2-4% of stellar mass being Population III massive metal-free stars at z~6. While there is no compelling physical reason or observational evidence to support (A), (B) could be fulfilled plausibly by continued existence of some pockets of uncontaminated, metal-free gas for star formation. (2) The volume-weighted neutral fraction of the IGM of <f_HI>_V~ 10^-4 at z=5.8 inferred from the SDSS observations of QSO absorption spectra provides enough information to ascertain that reionization is basically complete with at most ~0.1-1% of IGM that is un-ionized at z=5.8. (3) Barring some extreme evolution of the IMF, the neutral fraction of the IGM is expected to rise quickly toward high redshift from the point of HII bubble percolation, with the mean neutral fraction of the IGM expected to reach 6-12% at z=6.5, 13-27% at z=7.7 and 22-38% at z=8.8.
We study the clustering properties of z~5.7 Lyman-alpha emitters (LAEs) in a cosmological reionization simulation with a full Lya radiative transfer calculation. Lya radiative transfer substantially modifies the intrinsic Lya emission properties, com
pared to observed ones, depending on the density and velocity structure environment around the Lya emitting galaxy. This environment-dependent Lya selection introduces new features in LAE clustering, suppressing (enhancing) the line-of-sight (transverse) density fluctuations and giving rise to scale-dependent galaxy bias. In real space, the contours of the three-dimensional two-point correlation function of LAEs appear to be prominently elongated along the line of sight on large scales, an effect that is opposite to and much stronger than the linear redshift-space distortion effect. The projected two-point correlation function is greatly enhanced in amplitude by a factor of up to a few, compared to the case without the environment dependent selection effect. The new features in LAE clustering can be understood with a simple, physically motivated model, where Lya selection depends on matter density, velocity, and their gradients. We discuss the implications and consequences of the effects on galaxy clustering from Lya selection in interpreting clustering measurements and in constraining cosmology and reionization from LAEs.
We combine a cosmological reionization simulation with box size of 100Mpc/h on a side and a Monte Carlo Lyman-alpha (Lya) radiative transfer code to model Lyman Alpha Emitters (LAEs) at z~5.7. The model introduces Lya radiative transfer as the single
factor for transforming the intrinsic Lya emission properties into the observed ones. Spatial diffusion of Lya photons from radiative transfer results in extended Lya emission and only the central part with high surface brightness can be observed. Because of radiative transfer, the appearance of LAEs depends on density and velocity structures in circumgalactic and intergalactic media as well as the viewing angle, which leads to a broad distribution of apparent (observed) Lya luminosity for a given intrinsic Lya luminosity. Radiative transfer also causes frequency diffusion of Lya photons. The resultant Lya line is asymmetric with a red tail. The peak of the Lya line shifts towards longer wavelength and the shift is anti-correlated with the apparent to intrinsic Lya luminosity ratio. The simple radiative transfer model provides a new framework for studying LAEs. It is able to explain an array of observed properties of z~5.7 LAEs in Ouchi et al. (2008), producing Lya spectra, morphology, and apparent Lya luminosity function (LF) similar to those seen in observation. The broad distribution of apparent Lya luminosity at fixed UV luminosity provides a natural explanation for the observed UV LF, especially the turnover towards the low luminosity end. The model also reproduces the observed distribution of Lya equivalent width (EW) and explains the deficit of UV bright, high EW sources. Because of the broad distribution of the apparent to intrinsic Lya luminosity ratio, the model predicts effective duty cycles and Lya escape fractions for LAEs.
