We investigate the effects of Population III (Pop III) stars and their remnants on the cosmological 21-cm global signal. By combining a semi-analytic model of Pop III star formation with a global 21-cm simulation code, we investigate how X-ray and radio emission from accreting Pop III black holes may affect both the timing and depth of the 21-cm absorption feature that follows the initial onset of star formation during the Cosmic Dawn. We compare our results to the findings of the EDGES experiment, which has reported the first detection of a cosmic 21-cm signal. In general, we find that our fiducial Pop III models, which have peak star formation rate densities of $\sim 10^{-4}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ between $z \sim 10$ and $z \sim 15$, are able to match the timing of the EDGES signal quite well, in contrast to models that ignore Pop III stars. To match the unexpectedly large depth of the EDGES signal without recourse to exotic physics, we vary the parameters of emission from accreting black holes (formed as Pop III remnants) including the intrinsic strength of X-ray and radio emission as well as the local column density of neutral gas. We find that models with strong radio emission and relatively weak X-ray emission can self-consistently match the EDGES signal, though this solution requires fine-tuning. We are only able to produce signals with sharp features similar to the EDGES signal if the Pop~III IMF is peaked narrowly around $140 \, M_\odot$.

We present a semi-analytic model of star formation in the early universe, beginning with the first metal-free stars. By employing a completely feedback-limited star formation prescription, stars form at maximum efficiency until the self-consistently calculated feedback processes halt formation. We account for a number of feedback processes including a meta-galactic Lyman-Werner background, supernovae, photoionization, and chemical feedback. Halos are evolved combining mass accretion rates found through abundance matching with our feedback-limited star formation prescription, allowing for a variety of Population III (Pop III) initial mass functions (IMFs). We find that, for a number of models, massive Pop III star formation can continue on until at least $z \sim 20$ and potentially past $z \sim 6$ at rates of around $10^{-4}$ to $10^{-5}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$, assuming these stars form in isolation. At this point Lyman-Werner feedback pushes the minimum halo mass for star formation above the atomic cooling threshold, cutting off the formation of massive Pop III stars. We find that, in most models, Pop II and Pop III star formation co-exist over cosmological time-scales, with the total star formation rate density and resulting radiation background strongly dominated by the former before Pop III star formation finally ends. These halos form at most $\sim 10^3$ M$_\odot$ of massive Pop III stars during this phase and typically have absolute magnitudes in the range of $M_\text{AB} = -5 $ to $ -10$. We also briefly discuss how future observations from telescopes such as JWST or WFIRST and 21-cm experiments may be able to constrain unknown parameters in our model such as the IMF, star formation prescription, or the physics of massive Pop III stars.

We investigate the effects of Population III stars on the sky-averaged 21-cm background radiation, which traces the collective emission from all sources of ultraviolet and X-ray photons before reionization is complete. While UV photons from PopIII stars can in principle shift the onset of radiative coupling of the 21-cm transition -- and potentially reionization -- to early times, we find that the remnants of PopIII stars are likely to have a more discernible impact on the 21-cm signal than PopIII stars themselves. The X-rays from such sources preferentially heat the IGM at early times, which elongates the epoch of reheating and results in a more gradual transition from an absorption signal to emission. This gradual heating gives rise to broad, asymmetric wings in the absorption signal, which stand in contrast to the relatively sharp, symmetric signals that arise in models treating PopII sources only. A stronger signature of PopIII, in which the position of the absorption minimum becomes inconsistent with PopII-only models, requires extreme star-forming events that may not be physically plausible, lending further credence to predictions of relatively high frequency absorption troughs, $\nu_{\min} \sim 100$ MHz. As a result, though the trough location alone may not be enough to indicate the presence of PopIII, the asymmetric wings should arise even if only a few PopIII stars form in each halo before the transition to PopII star formation occurs, provided that the PopIII IMF is sufficiently top-heavy and at least some PopIII stars form in binaries.

Near-infrared surveys have now determined the luminosity functions of galaxies at 6<z<9 to impressive precision and identified a number of candidates at even earlier times. Here we develop a simple analytic model to describe these populations that allows physically-motivated extrapolation to earlier times and fainter luminosities. We assume that galaxies grow through accretion onto dark matter halos, which we model by matching halos at fixed number density across redshift, and that stellar feedback limits the star formation rate. We allow for a variety of feedback mechanisms, including regulation through supernova energy and momentum from radiation pressure. We show that reasonable choices for the feedback parameters can fit the available galaxy data, which in turn substantially limits the range of plausible extrapolations of the luminosity function to earlier times and fainter luminosities: for example, the global star formation rate declines rapidly at z>10, but the bright galaxies accessible to observations decline much faster than the total. Deviations from our predictions would provide evidence for new astrophysics within the first generations of galaxies. We also provide predictions for galaxy measurements by future facilities, including JWST and WFIRST.