We utilize GIZMO, coupled with newly developed sub-grid models for Population~III (Pop~III) and Population~II (Pop~II), to study the legacy of star formation in the pre-reionization Universe. We find that the Pop~II star formation rate density (SFRD), produced in our simulation ($\sim 10^{-2}\ \Msun{\rm yr^{-1}\, Mpc^{-3}}$ at $z\simeq 10$), matches the total SFRD inferred from observations within a factor of $<2$ at $7\lesssim z \lesssim10$. The Pop~III SFRD, however, reaches a plateau at $\sim10^{-3}\ \Msun{\rm yr^{-1}\, Mpc^{-3}}$ by $z\approx10$, remaining largely unaffected by the presence of Pop~II feedback. At $z$=7.5, $\sim 20\%$ of Pop~III star formation occurs in isolated haloes which have never experienced any Pop~II star formation (i.e. primordial haloes). We predict that Pop~III-only galaxies exist at magnitudes $M_{\rm UV}\gtrsim -11$, beyond the limits for direct detection with the {\it James Webb Space Telescope (JWST)}. We assess that our stellar mass function (SMF) and UV luminosity function (UVLF) agree well with the observed low mass/faint-end behavior at $z=8$ and $10$. However, beyond the current limiting magnitudes, we find that both our SMF and UVLF demonstrate a deviation/turnover from the expected power-law slope ($M_{\rm UV,turn}= -13.4\pm1.1$ at $z$=10). This could impact observational estimates of the true SFRD by a factor of $2 (10)$ when integrating to $M_{\rm UV} = -$12 ($-$8) at $z\sim 10$, depending on integration limits. Our turnover correlates well with the transition from dark matter haloes dominated by molecular cooling to those dominated by atomic cooling, for a mass $M_{\rm halo}\approx 10^{8} \Msun$ at $z\simeq 10$.

We use cosmological hydrodynamic zoom-in simulations to study early structure formation in two dark matter (DM) cosmologies, the standard CDM model, and a thermal warm DM (WDM) model with a particle mass of $m_{\chi}c^{2}=3\ \mathrm{keV}$. We focus on DM haloes with virial masses $M\sim 10^{10}\ M_{\odot}$. We find that the first star formation activity is delayed by $\sim 200\ \mathrm{Myr}$ in the WDM model, with similar delays for metal enrichment and the formation of the second generation of stars. However, the differences between the two models in globally-averaged properties, such as star formation rate density and mean metallicity, decrease towards lower redshifts ($z\lesssim 10$). Metal enrichment in the WDM cosmology is restricted to dense environments, while low-density gas can also be significantly enriched in the CDM case. The free-free contribution from early structure formation at redshifts $z>6$ to the cosmic radio background (CRB) is $3_{-1.5}^{+13}\%$ ($8_{-3.5}^{+33}\%$) of the total signal inferred from radio experiments such as ARCADE 2, in the WDM (CDM) model. The direct detection of the $\mathrm{H_{2}}$ emission from early structure formation ($z\gtrsim 7.2$), originating from the low-mass haloes explored here, will be challenging even with the next generation of far-infrared space telescopes, unless the signal is magnified by at least a factor of 10 via gravitational lensing or shocks. However, more massive haloes with $M\gtrsim 10^{12}\ M_{\odot}$ may be observable for $z\gtrsim 10$, even without magnification, provided that our extrapolation from the scale of our simulated haloes is valid.

Using cosmological volume simulations and a custom built sub-grid model for Pop~III star formation, we examine the baseline dust extinction in the first galaxies due to Pop~III metal enrichment in the first billion years of cosmic history. We find that while the most enriched, high-density lines of sight in primordial galaxies can experience a measurable amount of extinction from Pop~III dust ($E(B-V)_{\rm max}=0.07,\ A_{\rm V,max}\approx0.28$), the average extinction is very low with $\left< E(B-V) \right> \lesssim 10^{-3}$. We derive a power-law relationship between dark matter halo mass and extinction of $E(B-V)\propto M_{\rm halo}^{0.80}$. Performing a Monte Carlo parameter study, we establish the baseline reddening of the UV spectra of dwarf galaxies at high redshift due to Pop~III enrichment only. With this method, we find $\left<\beta_{\rm UV}\right>-2.51\pm0.07$, which is both nearly halo mass and redshift independent.

We utilize the hydrodynamic and N-body code {\small GIZMO} coupled with our newly developed sub-grid Population~III (Pop~III) Legacy model, designed specifically for cosmological volume simulations, to study the baseline metal enrichment from Pop~III star formation at $z>7$. In this idealized numerical experiment, we only consider Pop~III star formation. We find that our model Pop~III star formation rate density (SFRD), which peaks at $\sim 10^{-3}\ {\rm M_\odot yr^{-1} Mpc^{-1}}$ near $z\sim10$, agrees well with previous numerical studies and is consistent with the observed estimates for Pop~II SFRDs. The mean Pop~III metallicity rises smoothly from $z=25-7$, but does not reach the critical metallicity value, $Z_{\rm crit}=10^{-4}\ Z_\odot$, required for the Pop~III to Pop~II transition in star formation mode until $z\simeq7$. This suggests that, while individual halos can suppress in-situ Pop~III star formation, the external enrichment is insufficient to globally terminate Pop~III star formation. The maximum enrichment from Pop~III star formation in star forming dark matter halos is $Z\sim10^{-2}\ Z_\odot$, whereas the minimum found in externally enriched haloes is $Z\gtrsim10^{-7}\ Z_\odot$. Finally, mock observations of our simulated IGM enriched with Pop~III metals produce equivalent widths similar to observations of an extremely metal poor damped Lyman alpha (DLA) system at $z=7.04$, which is thought to be enriched by Pop~III star formation only.

Using the cosmological smoothed particle hydrodynamical code GADGET-3 we make a realistic assessment of the technique of using constant cumulative number density as a tracer of galaxy evolution at high redshift. We find that over a redshift range of $3\leq z \leq7$ one can on average track the growth of the stellar mass of a population of galaxies selected from the same cumulative number density bin to within $\sim 0.20$ dex. Over the stellar mass range we probe ($10^{10.39}\leq M_s/M_\odot \leq 10^{10.75}$ at $z =$ 3 and $10^{8.48}\leq M_s/M_\odot \leq 10^{9.55}$ at $z =$ 7) one can reduce this bias by selecting galaxies based on an evolving cumulative number density. We find the cumulative number density evolution exhibits a trend towards higher values which can be quantified by simple linear formulations going as $-0.10\Delta z$ for descendants and $0.12\Delta z$ for progenitors. Utilizing such an evolving cumulative number density increases the accuracy of descendant/progenitor tracking by a factor of $\sim2$. This result is in excellent agreement, within $0.10$ dex, with abundance matching results over the same redshift range. However, we find that our more realistic cosmological hydrodynamic simulations produce a much larger scatter in descendant/progenitor stellar masses than previous studies, particularly when tracking progenitors. This large scatter makes the application of either the constant cumulative number density or evolving cumulative number density technique limited to average stellar masses of populations only, as the diverse mass assembly histories caused by stochastic physical processes such as gas accretion, mergers, and star formation of individual galaxies will lead to a larger scatter in other physical properties such as metallicity and star-formation rate.

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z=4-8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Year 1 Hubble Frontier Field deep parallel observations. These surveys provides an effective volume of 0.6-1.2 x 10^6 Mpc^3 over this epoch, allowing us to perform a robust search for faint (M_UV=-18) and bright (M_UV < -21) galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5<z<8.5, with >1000 galaxies at z~6-8. We measure the luminosity function using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our results agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z>6, with M* ~ -21 at z>5, different than that inferred based on previous trends at lower redshift. At z=8, a single power-law provides an equally good fit to the UV luminosity function, while at z=6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift. We measure the evolution of the cosmic star-formation rate density, correcting for dust attenuation, and find that it declines as (1+z)^(-4.3 +/- 0.5) at z>4, consistent with observations at z>9. Our observations are consistent with a reionization history that starts at z>10, completes at z>6, and reaches a midpoint (x_HII = 0.5) at 6.7<z<9.4. Finally, our observations predict that the abundance of bright z=9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z~10 galaxies. [abridged]

It has been shown observationally that star formation (SF) correlates tightly with the presence of molecular hydrogen (H$_2$). Therefore it would be important to investigate its implication on galaxy formation in a cosmological context. In the present work, we track the H$_2$ mass fraction within our cosmological smoothed particle hydrodynamics (SPH) code GADGET-3 using an equilibrium analytic model by Krumholz et al. This model allows us to regulate the star formation in our simulation by the local abundance of H$_2$ rather than the total cold gas density, and naturally introduce the dependence of star formation on metallicity. We investigate implications of the equilibrium H$_2$-based SF model on galaxy population properties, such as the stellar-to-halo mass ratio (SHMR), baryon fraction, cosmic star formation rate density (SFRD), galaxy specific SFR, galaxy stellar mass functions (GSMF), and Kennicutt-Schmidt (KS) relationship. The advantage of our work over the previous ones is having a large sample of simulated galaxies in a cosmological volume from high-redshift to $z=0$. We find that low-mass halos with M$_{DM}<10^{10.5}$ M$_\odot$ are less efficient in producing stars in the H$_2$-based SF model at $z\geq6$, which brings the simulations to a better agreement with observational estimates of SHMR and GSMF at the low-mass end. This is particularly evident by a reduction in the number of low-mass galaxies at M$_\star\leq10^{8}$ M$_\odot$ in the GSMF. The overall SFRD is also reduced at high-$z$ in the H$_2$ run, which results in slightly higher SFRD at low-redshift due to more abundant gas available for star formation at later times. This new H$_2$ model is able to reproduce the empirical KS relationship at $z=0$ naturally without the need for setting its normalization by hand, and overall it seems to have more advantages than the previous pressure-based SF model.

We present the results of a numerical study examining the effect of H_2-based star formation (SF) model on the rest-frame UV luminosity function and star formation rate function (SFRF) of z>=6 galaxies, and the implications for reionization. Using cosmological hydrodynamical simulations outfitted with an H_2-SF model, we find good agreement with our previous results (non-H_2 SF model) and observations at Muv<=-18. However at Muv>-18, we find that the LF deviates from both our previous work and current observational extrapolations, producing significantly fewer low-luminosity galaxies and exhibiting additional turnover at the faint end. We constrain the redshift evolution of this turnover point using a modified Schechter function that includes additional terms to quantify the turnover magnitude (Muv^t) and subsequent slope (Beta). We find that Muv^t evolves from Muv^t=-17.33 (at z=8) to -15.38 (z=6), while Beta becomes shallower by \Delta Beta=0.22 during the same epoch. This occurs in an Muv range which will be observable by JWST. By integrating the SFRF, we determine that even though H_2-SF model significantly reduces the number density of low luminosity galaxies at Muv>-18, it does not suppress the total SFR density enough to affect the capability of SF to maintain reionization.

We examine the duty cycle and the history of star formation (SFH) for high-redshift galaxies at z>=6 using cosmological hydrodynamic simulations. We find that, even though individual galaxies have bursty SFH, the averaged SFH between z~15 to z=6 can be characterized well by either an exponentially increasing functional form with characteristic time-scales of 70 Myr to 200 Myr for galaxies with stellar masses Ms~10^6 Msun to >10^10 Msun respectively, or by a simple power-law form which exhibits a similar mass dependent time-scales. Using the SFH of individual galaxies, we measure the duty cycle of star formation (DC_SFH); i.e., the fraction of time a galaxy of a particular mass spends above a star formation rate (SFR) threshold which would make it observable to the Hubble Space Telescope (HST) during a given epoch. We also examine the fraction of galaxies at a given redshift that are brighter than a rest-frame UV magnitude (Muv ~ -18), which is sufficient enough to make them observable (DC_Muv). We find that both DC_SFH and DC_Muv make a sharp transition from zero (for galaxies with Ms <= 10^7 Msun) to unity (for Ms > 10^9 Msun). The measured duty cycle is also manifested in the intrinsic scatter in the Ms-SFR relationship (~ 1 dex) and Ms-Muv relationship (\Delta Muv ~ +-1 mag). We provide analytic fits to the DC as a function of Ms using a sigmoid function, which can be used to correct for catalogue incompleteness. We consider the effects of duty cycle to the observational estimate of galaxy stellar mass functions (GSMF) and the star formation rate density (SFRD), and find that it results in a much shallower low-mass end slopes of the GSMF and a reduction of >~ 70% of our intrinsic SFRD, making our simulation results more compatible with observational estimates.

We present the results of a numerical study comparing photometric and physical properties of simulated z=6-9 galaxies to the observations taken by the WFC3 instrument aboard the Hubble Space Telescope. Using cosmological hydrodynamical simulations we find good agreement with observations in color-color space at all studied redshifts. We also find good agreement between observations and our Schechter luminosity function fit in the observable range, Muv<= -18, provided that a moderate dust extinction effect exists for massive galaxies. However beyond what currently can be observed, simulations predict a very large number of low-mass galaxies and evolving steep faint-end slopes from alpha_L = -2.15 at z=6 to alpha_L = -2.64 at z=9, with a dependence of |alpha_L| \propto (1+z)^0.59. During the same epoch, the normalization phi* increases and the characteristic magnitude Muv* becomes moderately brighter with decreasing redshift. We find similar trends for galaxy stellar mass function with evolving low-mass end slope from alpha_M = - 2.26 at z=6 to alpha_M = -2.87 at z=9, with a dependence of |alpha_M| \propto (1+z)^0.65. Together with our recent result on the high escape fraction of ionizing photons for low-mass galaxies, our results suggest that the low-mass galaxies are important contributor of ionizing photons for the reionisation of the Universe at z>=6.