The launch of the James Webb Space Telescope will open up a new window for observations at the highest redshifts, reaching out to $z\approx 15$. However, even with this new facility, the first stars will remain out of reach, as they are born in small minihalos with luminosities too faint to be detected even by the longest exposure times. In this Letter, we investigate the basic properties of the "Ultimately Large Telescope", a facility that can detect Population III stars. Observations will take place in the near-infrared and therefore a moon-based facility is proposed. An instrument needs to reach magnitudes as faint as 39mag$_\mathrm{AB}$, corresponding to a primary mirror size of about 100m in diameter. Assuming JWST NIRCam filters, we estimate that Pop~III sources will have unique signatures in a colour-colour space and can be identified unambiguously.

We construct a theoretical framework to study Population III (Pop III) star formation in the post-reionization epoch ($z\lesssim 6$) by combining cosmological simulation data with semi-analytical models. We find that due to radiative feedback (i.e. Lyman-Werner and ionizing radiation) massive haloes ($M_{\rm halo}\gtrsim 10^{9}\ \rm M_{\odot}$) are the major ($\gtrsim 90$\%) hosts for potential Pop III star formation at $z\lesssim 6$, where dense pockets of metal-poor gas may survive to form Pop III stars, under inefficient mixing of metals released by supernovae. Metal mixing is the key process that determines not only when Pop III star formation ends, but also the total mass, $M_{\rm PopIII}$, of \textit{active} Pop III stars per host halo, which is a crucial parameter for direct detection and identification of Pop III hosts. Both aspects are still uncertain due to our limited knowledge of metal mixing during structure formation. Current predictions range from early termination at the end of reionization ($z\sim 5$) to continuous Pop III star formation extended to $z=0$ at a non-negligible rate $\sim 10^{-7}\ \rm M_{\odot}\ yr^{-1}\ Mpc^{-3}$, with $M_{\rm PopIII}\sim 10^{3}-10^{6}\ \rm M_{\odot}$. This leads to a broad range of redshift limits for direct detection of Pop III hosts, $z_{\rm PopIII}\sim 0.5-12.5$, with detection rates $\lesssim 0.1-20\ \rm arcmin^{-2}$, for current and future space telescopes (e.g. HST, WFIRST and JWST). Our model also predicts that the majority ($\gtrsim 90$\%) of the cosmic volume is occupied by metal-free gas. Measuring the volume filling fractions of this metal-free phase can constrain metal mixing parameters and Pop III star formation.

We use cosmological hydrodynamic simulations to study the gravitational wave (GW) signals from high-redshift binary black holes (BBHs) formed by dynamical capture (ex-situ formation channel). We in particular focus on BHs originating from the first generation of massive, metal-poor, so-called Population~III (Pop~III) stars. An alternative (in-situ) formation pathway arises in Pop~III binary stars, whose GW signature has been intensively studied. In our optimistic model, we predict a local GW event rate density for ex-situ BBHs (formed at z>4) of ~0.04 /(yr Gpc^3). This is comparable to or even higher than the conservative predictions of the rate density for in-situ BBHs ~0.01-0.1 /(yr Gpc^3), indicating that the ex-situ formation channel may be as important as the in-situ one for producing GW events. We also evaluate the detectability of our simulated GW events for selected planned GW instruments, such as the Einstein Telescope (ET). For instance, we find the all-sky detection rate with signal-to-noise ratios above 10 to be <100 per year for the xylophone configuration of ET. However, our results are highly sensitive to the sub-grid models for BBH identification and evolution, such that the GW event efficiency (rate) is reduced by a factor of 4 (20) in the pessimistic case. The ex-situ channel of Pop III BBHs deserves further investigation with better modeling of the environments around Pop III-seeded BHs.

We demonstrate a new approach of indirectly constraining both early star and structure formation via mature galaxy clusters at cosmic noon (z~2), using the cluster XLSSC 122 as an example. With the standard Press Schechter formalism, we infer a rapid evolution of the star formation efficiency (the ratio of stellar to halo mass) from 10^-4 to 0.01 during z~20-13, based on the age distribution of stars in post-starburst galaxies of XLSSC 122, measured by HST photometry assuming no dust extinction. Here, we consider all low-mass haloes, including minihaloes, that host the first stars and galaxies (5*10^5 Msun < M_halo < 10^10 Msun). We also place new constraints on fuzzy dark matter models of m_a < 5*10^-21 eV/c^2 for the ultra-light boson mass, from the abundance of galaxies with star formation at z > 13 in XLSSC 122. Our exploratory results are consistent with existing constraints. More comprehensive results will be obtained if our approach is extended to a large sample of clusters or field post-starburst galaxies at cosmic noon, with improved modelling of halo and stellar populations.

High-redshift Lyman-alpha blobs (LABs) are an enigmatic class of objects that have been the subject of numerous observational and theoretical investigations. It is of particular interest to determine the dominant power sources for the copious luminosity, as direct emission from HII regions, cooling gas, and fluorescence due to the presence of active galactic nuclei (AGN) can all contribute significantly. In this paper, we present the first theoretical model to consider all of these physical processes in an attempt to develop an evolutionary model for the origin of high-z LABs. This is achieved by combining a series of high-resolution cosmological zoom-in simulations with ionization and Lyman-alpha (Lya) radiative transfer models. We find that massive galaxies display a range of Lya luminosities and spatial extents (which strongly depend on the limiting surface brightness used) over the course of their lives, though regularly exhibit luminosities and sizes consistent with observed LABs. The model LABs are typically powered from a combination of recombination in star-forming galaxies, as well as cooling emission from gas associated with accretion. When AGN are included in the model, the fluorescence caused by AGN-driven ionization can be a significant contributor to the total Lya luminosity as well. We propose that the presence of an AGN may be predicted from the Gini coefficient of the blob's surface brightness. Within our modeled mass range, there are no obvious threshold physical properties that predict appearance of LABs, and only weak correlations of the luminosity with the physical properties of the host galaxy. This is because the emergent Lya luminosity from a system is a complex function of the gas temperature, ionization state, and Lya escape fraction.

We explore the possibility of detecting the first galaxies with the next generation Origins Space Telescope (OST) by applying an analytical model of primordial dust emission. By analysing source densities as a function of redshift (z), and considering deep-field exposures with the Origins Survey Spectrometer, we estimate that the redshift horizon for detecting one individual source would be above z~7 for systems with dust-to-metal ratios higher than those expected for typical primeval galaxies. On the other hand, if confusion limits could be overcome, the Far-infrared Imager and Polarimeter would enable the detection of more typical fainter systems at z>7. Given the dependence of the obtained results with the properties of primeval dust, we conclude that the OST could provide important clues to the nature of the interestellar medium in the early Universe.

We review the current status of knowledge concerning the early phases of star formation during cosmic dawn. This includes the first generations of stars forming in the lowest mass dark matter halos in which cooling and condensation of gas with primordial composition is possible at very high redshift ($z > 20$), namely metal-free Population III stars, and the first generation of massive black holes forming at such early epochs, the so-called black hole seeds. The formation of black hole seeds as end states of the collapse of Population III stars, or via direct collapse scenarios, is discussed. In particular, special emphasis is given to the physics of supermassive stars as potential precursors of direct collapse black holes, in light of recent results of stellar evolution models, and of numerical simulations of the early stages of galaxy formation. Furthermore, we discuss the role of the cosmic radiation produced by the early generation of stars and black holes at high redshift in the process of reionization.

We compare observed far-infrared (FIR) galaxy spectral energy distributions (SEDs) at $z\sim6$ and theoretical predictions for first massive population II (Pop II) galaxies. Observed FIR SEDs at $z\gtrsim5$ are broad and shifted to bluer wavelengths when compared to galaxies at $z\sim3$. By implementing an analytical model for dust emission from Pop II massive galaxies,we were able to reproduce the observed behaviour as a consequence of the high energy densities and silicate-rich dust composition of high-$z$ model galaxies. As noted in a previous work, the non-blackbody nature of galaxy SEDs at $z \sim 6$ should be taken into account when interpreting measurements of FIR luminosities to avoid underestimating star formation rates.

We explore the possibility of detecting first galaxies with a generic far-infrared/sub-millimeter telescope by applying an analytical model of primordial dust emission. As shown in previous works, galaxies at redshifts $z>7$ experience a strong negative K-correction in such a way that systems of similar masses are brighter at higher $z$. In addition, at a given mass and $z$, our model predicts that luminosity fluxes increase proportionally to the dust-to-metal ratios ($D/M$) of primeval sources. We evaluate the observability of model sources at different $z>7$ as a function of the observed survey area ($\Delta \Omega$) and sensitivity ($S$) of a generic instrument. Assuming $\Delta \Omega \sim 10~{\rm deg}^2$ and a plausible $S \sim 1~{\mu}{\rm Jy}$ for a near future survey, we could assure the detection of at least one typical source with $D/M \sim 5\times10^{-3}$ at $z>7$. For $S\gtrsim 1~{\mu}{\rm Jy}$ and $\Delta \Omega \lesssim 10~{\rm deg}^2$, higher than typical $D/M$ are required to detect at least one individual source at $z>7$. The observability of model galaxies is also affected by the size distribution of dust grains, specially towards higher $z$.

We explore the possibility of detecting the first galaxies with the next generation of space-based far infrared (FIR) telescopes by applying an analytical model of primordial dust emission. Our results indicate that FIR/sub-mm sources at $z \gtrsim 7$ will experience a strong negative K-correction. Systems of a given virial mass would exhibit larger dust luminosities at higher $z$, as a consequence of the increase in dust temperature driven by the higher temperature floor set by the cosmic microwave background. In addition, high-$z$ systems are more concentrated, which enhances the heating efficiency associated with stellar radiation. By analysing source densities as a function of $z$, and considering survey areas of 0.1 ${\rm deg}^2$ and 10 ${\rm deg}^2$, we find that the redshift horizon for detecting at least one source would be above $z\sim 7$ for instrument sensitivities $\lesssim 0.1-0.5 \ {\mu}{\rm Jy}$ and $\lesssim 0.5-3.0 \ {\mu}{\rm Jy}$, respectively, with the exact values depending on the nature of primordial dust. However, galaxy populations with higher than typical metallicities, star formation efficiencies and/or dust-to-metal ratios could relax such sensitivity requirements. In addition, the redshift horizon shows a significant dependence on the nature of primordial dust. We conclude that future FIR campaigns could play a crucial role in exploring the nature of dust and star formation in the early universe.

We derive new constraints on the non-gravitational baryon-dark-matter scattering (BDMS) by evaluating the mass thresholds of dark matter (DM) haloes in which primordial gas can cool efficiently to form Population III (Pop III) stars, based on the timing of the observed 21-cm absorption signal. We focus on the BDMS model with interaction cross-section $\sigma=\sigma_{1}[v/(1\ \mathrm{km\ s^{-1}})]^{-4}$, where $v$ is the relative velocity of the encounter. Our results rule out the region in parameter space with $\sigma_{1}\gtrsim 10^{-19}\ \mathrm{cm^{2}}$ and DM particle mass $m_{\chi}c^{2}\lesssim 3\times 10^{-2}\ \mathrm{GeV}$, where the cosmic number density of Pop III hosts at redshift $z\sim 20$ is at least three orders of magnitude smaller than in the standard Lambda cold DM ($\Lambda$CDM) case. In these BDMS models, the formation of Pop III stars is significantly suppressed for $z\gtrsim 20$, inconsistent with the timing of the observed global 21-cm absorption signal. For the fiducial BDMS model with $m_{\chi}c^{2}=0.3$ GeV and $\sigma_{1}=8\times 10^{-20}\ \mathrm{cm^{2}}$, capable of accommodating the measured absorption depth, the number density of Pop III hosts is reduced by a factor of $3-10$ at $z\sim 15-20$, when the 21-cm signal is imprinted, compared with the $\Lambda$CDM model. The confluence of future detailed cosmological simulations with improved 21-cm observations promises to probe the particle-physics nature of DM at the small-scale frontier of early structure formation.

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$.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

We investigate the absorption features associated with a gas-rich dwarf galaxy using cosmological hydrodynamics simulations. Our goal is to explore whether the progenitors of the lowest mass dwarf galaxies known to harbor neutral hydrogen today (M_star~10^6 solar mass, M_halo=4x10^9 solar mass) could possibly be detected as Damped Lyman-alpha Absorbers (DLAs) over cosmic time. We trace the evolution of a single dwarf galaxy, pre-selected to contain DLAs, from the era of the first metal-free, so-called Population~III (Pop~III), stars, down to z=0, thus allowing us to study the metal enrichment history of DLAs associated with the simulated galaxy. We find that the progenitors of the simulated dwarf are expected to be seen for most of their evolution as DLAs that are contaminated by normal, Population~II, stars. The time period during which DLAs are only metal-enriched by Pop~III stars, on the other hand, is likely very brief, confined to high redshifts, z~6. The susceptibility of the dwarfs to the external UV radiation background allows them to preserve neutral gas only at the centre (a few ~100 pc). This results in a small probability that the simulated dwarf would be observed as a DLA. This study suggests that DLAs are unlikely to be hosted in the lowest mass dwarfs that can harbor neutral gas (M_halo~ 4x10^9 solar mass), below which neutral gas is unlikely to exist. However, this study does illustrate that, when detected, absorption lines provide a powerful method for probing ISM conditions inside the smallest dwarf galaxies at intermediate to high redshifts.

Within standard $\Lambda$CDM cosmology, Population III (Pop III) star formation in minihalos of mass $M_\mathrm{halo}\gtrsim 5\times10^5$ M$_\odot$ provides the first stellar sources of Lyman$\alpha$ (Ly$\alpha$) photons. The Experiment to Detect the Global Epoch of Reionization Signature (EDGES) has measured a strong absorption signal of the redshifted 21 cm radiation from neutral hydrogen at $z\approx 17$, requiring efficient formation of massive stars before then. In this paper, we investigate whether star formation in minihalos plays a significant role in establishing the early Ly$\alpha$ background required to produce the EDGES absorption feature. We find that Pop III stars are important in providing the necessary Ly$\alpha$-flux at high redshifts, and derive a best-fitting average Pop III stellar mass of $\sim$ 750M$_\odot{}$ per minihalo, corresponding to a star formation efficiency of 0.1%. Further, it is important to include baryon-dark matter streaming velocities in the calculation, to limit the efficiency of Pop~III star formation in minihalos. Without this effect, the cosmic dawn coupling between 21 cm spin temperature and that of the gas would occur at redshifts higher than what is implied by EDGES.

The emergence of the first black holes during the first billion years of cosmic history marks a key event in cosmology. Their formation is part of the overall process of ending the cosmic dark ages, when the first stars appeared in low-mass dark matter haloes about a few 100 million years after the Big Bang. The first stars, galaxies, and black holes transformed the Universe from its simple initial state into one of ever increasing complexity. We review recent progress on our emerging theoretical picture of how the first black holes appeared on the cosmic scene, and how they impacted the subsequent history of the Universe. Our focus is on supermassive black holes, in particular assessing possible pathways to the formation of the billion-solar-mass black holes inferred to power the luminous quasars at high redshifts. We conclude with a discussion of upcoming empirical probes, such as the James Webb Space Telescope (JWST), and the Laser Interferometer Space Antenna (LISA), further ahead in time.

We examine the assembly process and the observability of a first galaxy (M_vir~10^9 solar mass at z~8) with cosmological zoom-in, hydrodynamic simulations, including the radiative, mechanical, and chemical feedback exerted by the first generations of stars. To assess the detectability of such dwarf systems with the upcoming James Webb Space Telescope (JWST), we construct the spectral energy distribution for the simulated galaxy in a post-processing fashion. We find that while the non-ionizing UV continuum emitted by the simulated galaxy is expected to be below the JWST detection limit, the galaxy might be detectable using its nebular emission, specifically in the H-alpha recombination line. This requires that the galaxy experiences an active starburst with a star formation rate of Mdot_star > 0.1 Msun/yr at z~9. Due to the bursty nature of star formation in the first galaxies, the time interval for strong nebular emission is short, less than 2-3 Myr. The probability of capturing such primordial dwarf galaxies during the observable part of their duty cycle is thus low, resulting in number densities of order one source in a single pointing with MIRI onboard the JWST, for very deep exposures. Gravitational lensing, however, will boost their observability beyond this conservative baseline. The first sources of light will thus come firmly within our reach.

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.

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

A major scientific goal of JWST is to probe the epoch of re-ionization of the Universe at z above 6, and up to 20 and beyond. At these redshifts, galaxies are just beginning to form and the observable objects are early black holes, supernovae, and cosmic infrared background. The JWST has the necessary sensitivity to observe these targets individually, but a public deep and wide science enabling survey in the wavelength range from 2-5 $\mu$m is needed to discover these black holes and supernovae and to cover the area large enough for cosmic infrared background to be reliably studied. This enabling survey will also discover a large number of other transients and enable sciences such as supernova cosmology up to z $\sim$ 5, star formation history at high redshift through supernova explosions, faint stellar objects in the Milky Way, and galaxy evolution up to z approaching 10. The results of this survey will also serve as an invaluable target feeder for the upcoming era of ELT and SKA.

Galaxy formation depends on a complex interplay between gravitational collapse, gas accretion, merging, and feedback processes. Yet, after many decades of investigation, these concepts are poorly understood. This paper presents the argument that warm H$_2$ can be used as a tool to unlock some of these mysteries. Turbulence, shocks and outflows, driven by star formation, AGN activity or inflows, may prevent the rapid buildup of star formation in galaxies. Central to our understanding of how gas is converted into stars is the process by which gas can dissipate its mechanical energy through turbulence and shocks in order to cool. H$_2$ lines provide direct quantitative measurements of kinetic energy dissipation in molecular gas in galaxies throughout the Universe. Based on the detection of very powerful H$_2$ lines from z = 2 galaxies and proto-clusters at the detection limits of {\it Spitzer}, we are confident that future far-IR and UV H$_2$ observations will provide a wealth of new information and insight into galaxy evolution to high-z. Finally, at the very earliest epoch of star and galaxy formation, warm H$_2$ may also provide a unique glimpse of molecular gas collapse at 7 $<$ z $<$ 12 in massive dark matter (DM) halos on their way to forming the very first galaxies. Such measurements are beyond the reach of existing and planned observatories.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

While the upcoming telescopes will reveal correspondingly fainter, more distant galaxies, a question will persist: what more is there that these telescopes cannot see? One answer is the source-subtracted Cosmic Infrared Background (CIB). The CIB is comprised of the collective light from all sources remaining after known, resolved sources are accounted for. Ever-more-sensitive surveys will identify the brightest of these, allowing them to be removed, and - like peeling layers off an onion - reveal deeper layers of the CIB. In this way it is possible to measure the contributions from populations not accessible to direct telescopic observation. Measurement of fluctuations in the source-subtracted CIB, i.e., the spatial power spectrum of the CIB after subtracting resolved sources, provides a robust means of characterizing its faint, and potentially new, populations. Studies over the past 15 years have revealed source-subtracted CIB fluctuations on scales out to ~100' which cannot be explained by extrapolating from known galaxy populations. Moreover, they appear highly coherent with the unresolved Cosmic X-ray Background, hinting at a significant population of accreting black holes among the CIB sources. Characterizing the source-subtracted CIB with high accuracy, and thereby constraining the nature of the new populations, is feasible with upcoming instruments and would produce critically important cosmological information in the next decade. New coextensive deep and wide-area near-infrared, X-ray, and microwave surveys will bring decisive opportunities to examine, with high fidelity, the spatial spectrum and origin of the CIB fluctuations and their cross-correlations with cosmic microwave and X-ray backgrounds, and determine the formation epochs and the nature of the new sources (stellar nucleosynthetic or accreting black holes).

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

Massive Population II galaxies undergoing the first phase of vigorous star formation after the initial Population III stage should have high energy densities and silicate-rich interstellar dust. We have modeled the resulting far-infrared spectral energy distributions (SEDs), demonstrating that they are shifted substantially to bluer (`warmer') wavelengths relative to the best fitting ones at z ~ 3, and with strong outputs in the 10 - 40 micron range. When combined with a low level of emission by carbon dust, their SEDs match that of Haro 11, a local moderately-low-metallicity galaxy undergoing a very young and vigorous starburst that is likely to approximate the relevant conditions in young Population II galaxies. We expect to see similar SEDs at high redshifts (z >= 5) given the youth of galaxies at this epoch. In fact, we find a progression with redshift in observed galaxy SEDs, from those resembling local ones at 2 < z < 4 to a closer resemblance with Haro 11 at 5 < z < 7. In addition to the insight on conditions in high redshift galaxies, this result implies that estimates of the total infrared luminosities at z ~ 6 based on measurements near lambda ~ 1 mm can vary by factors of 2 - 4, depending on the SED template used. Currently popular modified blackbodies or local templates can result in significant underestimates compared with the preferred template based on the SED of Haro 11.

Lyman-alpha (Ly{\alpha}) photons from ionizing sources and cooling radiation undergo a complex resonant scattering process that generates unique spectral signatures in high-redshift galaxies. We present a detailed Ly{\alpha} radiative transfer study of a cosmological zoom-in simulation from the Feedback In Realistic Environments (FIRE) project. We focus on the time, spatial, and angular properties of the Ly{\alpha} emission over a redshift range of z = 5-7, after escaping the galaxy and being transmitted through the intergalactic medium (IGM). Over this epoch, our target galaxy has an average stellar mass of $M_{\rm star} \approx 5 \times 10^8 {\rm M}_\odot$. We find that many of the interesting features of the Ly{\alpha} line can be understood in terms of the galaxy's star formation history. The time variability, spatial morphology, and anisotropy of Ly{\alpha} properties are consistent with current observations. For example, the rest frame equivalent width has a ${\rm EW}_{{\rm Ly}\alpha,0} > 20 {\rm \AA}$ duty cycle of 62% with a non-negligible number of sightlines with $> 100 {\rm \AA}$, associated with outflowing regions of a starburst with greater coincident UV continuum absorption, as these conditions generate redder, narrower (or single peaked) line profiles. The lowest equivalent widths correspond to cosmological filaments, which have little impact on UV continuum photons but efficiently trap Ly{\alpha} and produce bluer, broader lines with less transmission through the IGM. We also show that in dense self-shielding, low-metallicity filaments and satellites Ly{\alpha} radiation pressure can be dynamically important. Finally, despite a significant reduction in surface brightness with increasing redshift, Ly{\alpha} detections and spectroscopy of high-$z$ galaxies with the upcoming James Webb Space Telescope is feasible.

Elucidating the particle physics nature of dark matter (DM) is one of the great challenges in modern science. The current lack of any direct DM detections in the laboratory heightens the need for astrophysical constraints, extending the search to DM models beyond the popular weakly interacting massive particle (WIMP) scenario. We here apply the classical Rees-Ostriker-Silk cooling criterion for galaxy formation to models with ultralight axion DM, also known as fuzzy dark matter (FDM). The resulting constraints provide a heuristic framework for upcoming observations, and our approximate analysis motivates the need for future, self-consistent simulations of FDM structure formation. We use observational constraints for the DM hosts of ultra faint dwarf (UFD) galaxies in the Local Group, together with the redshift constraints for the onset of primordial star formation from the recent EDGES 21-cm cosmology measurement, to illustrate this approach. We find that the existing constraints are straightforward to reconcile with standard $\Lambda$CDM, but disfavour FDM axion masses below $\sim 10^{-21}\,{\rm eV}/c^2$. The future potential for harnessing astrophysical probes of DM particle physics is compelling.

We present a suite of six fully cosmological, three-dimensional simulations of the collapse of an atomic cooling halo in the early Universe. We use the moving-mesh code arepo with an improved primordial chemistry network to evolve the hydrodynamical and chemical equations. The addition of a strong Lyman-Werner background suppresses molecular hydrogen cooling and permits the gas to evolve nearly isothermally at a temperature of about 8000 K. Strong gravitational torques effectively remove angular momentum and lead to the central collapse of gas, forming a supermassive protostar at the center of the halo. We model the protostar using two methods: sink particles that grow through mergers with other sink particles, and a stiff equation of state that leads to the formation of an adiabatic core. We impose threshold densities of $10^8$, $10^{10}$, and $10^{12}\,\text{cm}^{-3}$ for the sink particle formation and the onset of the stiff equation of state to study the late, intermediate, and early stages in the evolution of the protostar, respectively. We follow its growth from masses $\simeq 10\,\text{M}_\odot$ to $\simeq 10^5\,\text{M}_\odot$, with an average accretion rate of $\langle\dot{M}_\star\rangle \simeq 2\,\text{M}_\odot\,\text{yr}^{-1}$ for sink particles, and $\simeq 0.8 - 1.4\,\text{M}_\odot\,\text{yr}^{-1}$ for the adiabatic cores. At the end of the simulations, the HII region generated by radiation from the central object has long detached from the protostellar photosphere, but the ionizing radiation remains trapped in the inner host halo, and has thus not yet escaped into the intergalactic medium. Fully coupled, radiation-hydrodynamics simulations hold the key for further progress.

The recent detection of the sky-averaged 21-cm cosmological signal indicates a stronger absorption than the maximum allowed value based on the standard model. One explanation for the required colder primordial gas is the energy transfer between the baryon and dark matter fluids due to non-gravitational scattering. Here, we explore the thermal evolution of primordial gas, collapsing to form Population III (Pop III) stars, when this energy transfer is included. Performing a series of one-zone calculations, we find that the evolution results in stars more massive than in the standard model, provided that the dark matter is described by the best-fit parameters inferred from the 21-cm observation. On the other hand, a significant part of the dark matter parameter space can be excluded by the requirement to form massive Pop III stars sufficiently early in cosmic history. Otherwise, the radiation background needed to bring about the strong Wouthuysen-Field coupling at z >~ 17, inferred to explain the 21-cm absorption feature, could not be builtup. Intriguingly, the independent constraint from the physics of first star formation at high densities points to a similarly narrow range in dark matter properties. This exploratory study has to be followed-up with self-consistent three-dimensional simulations for a more rigorous derivation.

Due to its accuracy and generality, Monte Carlo radiative transfer (MCRT) has emerged as the prevalent method for Ly$\alpha$ radiative transfer in arbitrary geometries. The standard MCRT encounters a significant efficiency barrier in the high optical depth, diffusion regime. Multiple acceleration schemes have been developed to improve the efficiency of MCRT but the noise from photon packet discretization remains a challenge. The discrete diffusion Monte Carlo (DDMC) scheme has been successfully applied in state-of-the-art radiation hydrodynamics (RHD) simulations. Still, the established framework is not optimal for resonant line transfer. Inspired by the DDMC paradigm, we present a novel extension to resonant DDMC (rDDMC) in which diffusion in space and frequency are treated on equal footing. We explore the robustness of our new method and demonstrate a level of performance that justifies incorporating the method into existing Ly$\alpha$ codes. We present computational speedups of $\sim 10^2$-$10^6$ relative to contemporary MCRT implementations with schemes that skip scattering in the core of the line profile. This is because the rDDMC runtime scales with the spatial and frequency resolution rather than the number of scatterings - the latter is typically $\propto \tau_0$ for static media, or $\propto (a \tau_0)^{2/3}$ with core-skipping. We anticipate new frontiers in which on-the-fly Ly$\alpha$ radiative transfer calculations are feasible in 3D RHD. More generally, rDDMC is transferable to any computationally demanding problem amenable to a Fokker-Planck approximation of frequency redistribution.

Pair-instability supernovae (PISNe) are very luminous explosions of massive, low metallicity stars. They can potentially be observed out to high redshifts due to their high explosion energies, thus providing a probe of the Universe prior to reionization. The near-infrared camera, NIRCam, on board the James Webb Space Telescope is ideally suited for detecting their redshifted ultraviolet emission. We calculate the photometric signature of high-redshift PISNe and derive the optimal detection strategy for identifying their prompt emission and possible afterglow. We differentiate between PISNe and other sources that could have a similar photometric signature, such as active galactic nuclei or high-redshift galaxies. We demonstrate that the optimal strategy, which maximizes the visibility time of the PISN lightcurve per invested exposure time, consists of the two wide-band filters F200W and F356W with an exposure time of 600s. For such exposures, we expect one PISN at $z \lesssim 7.5$ per at least 50,000 different field of view, which can be accomplished with parallel observations and an extensive archival search. The PISN afterglow, caused by nebular emission and reverberation, is very faint and requires unfeasibly long exposure times to be uniquely identified. However, this afterglow would be visible for several hundred years, about two orders of magnitude longer than the prompt emission, rendering PISNe promising targets for future, even more powerful telescopes.

We present a model for the evolution of supermassive protostars from their formation at $M_\star \simeq 0.1\,\text{M}_\odot$ until their growth to $M_\star \simeq 10^5\,\text{M}_\odot$. To calculate the initial properties of the object in the optically thick regime we follow two approaches: based on idealized thermodynamic considerations, and on a more detailed one-zone model. Both methods derive a similar value of $n_{\rm F} \simeq 2 \times 10^{17} \,\text{cm}^{-3}$ for the density of the object when opacity becomes important, i.e. the opacity limit. The subsequent evolution of the growing protostar is determined by the accretion of gas onto the object and can be described by a mass-radius relation of the form $R_\star \propto M_\star^{1/3}$ during the early stages, and of the form $R_\star \propto M_\star^{1/2}$ when internal luminosity becomes important. For the case of a supermassive protostar, this implies that the radius of the star grows from $R_\star \simeq 0.65 \,{\rm AU}$ to $R_\star \simeq 250 \,{\rm AU}$ during its evolution. Finally, we use this model to construct a sub-grid recipe for accreting sink particles in numerical simulations. A prime ingredient thereof is a physically motivated prescription for the accretion radius and the effective temperature of the growing protostar embedded inside it. From the latter, we can conclude that photo-ionization feedback can be neglected until very late in the assembly process of the supermassive object.

We investigate the rotation velocity of the first stars by modelling the angular momentum transfer in the primordial accretion disc.Assessing the impact of magnetic braking, we consider the transition in angular momentum transport mode at the Alfv$\acute{\rm e}$n radius, from the dynamically dominated free-fall accretion to the magnetically dominated solid-body one.The accreting protostar at the centre of the primordial star-forming cloud rotates with close to breakup speed in the case without magnetic fields.Considering a physically-motivated model for small-scale turbulent dynamo amplification, we find that stellar rotation speed quickly declines if a large fraction of the initial turbulent energy is converted to magnetic energy ($\gtrsim 0.14$). Alternatively, if the dynamo process were inefficient, for amplification due to flux-freezing, stars would become slow rotators if the pre-galactic magnetic field strength is above a critical value, $\simeq 10^{-8.2}$G, evaluated at a scale of $n_{\rm H} = 1 {\rm cm}^{-3}$, which is significantly higher than plausible cosmological seed values ($\sim 10^{-15}$G). Because of the rapid decline of the stellar rotational speed over a narrow range in model parameters, the first stars encounter a bimodal fate: rapid rotation at almost the breakup level, or the near absence of any rotation.

We complete the formulation of the standard model of first star formation by exploring the possible impact of $\mathrm{LiH}$ cooling, which has been neglected in previous simulations of non-linear collapse. Specifically, we find that at redshift $z\gtrsim 5$, the cooling by $\mathrm{LiH}$ has no effect on the thermal evolution of shocked primordial gas, and of collapsing primordial gas into minihaloes or relic HII regions, even if the primordial lithium abundance were enhanced by one order of magnitude. Adding the most important lithium species to a minimum network of primordial chemistry, we demonstrate that insufficient $\mathrm{LiH}$ is produced in all cases considered, about $[\mathrm{LiH/Li}]\sim 10^{-9}$ for $T\lesssim 100$ K. Indeed, $\mathrm{LiH}$ cooling would only be marginally significant in shocked primordial gas for the highly unlikely case that the $\mathrm{LiH}$ abundance were increased by nine orders of magnitude, implying that $all$ lithium would have to be converted into $\mathrm{LiH}$. In this study, photo-destruction processes are not considered, and the collisional disassociation rate of $\mathrm{LiH}$ is possibly underestimated, rendering our results an extreme upper limit. Therefore, the cooling by $\mathrm{LiH}$ can safely be neglected for the thermal evolution of Population~III star-forming gas.

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.

JWST was conceived and built to answer one of the most fundamental questions that humans can address empirically: "How did the Universe make its first stars?". Our First Lights At REionization (FLARE) project transforms the quest for the epoch of reionization from the static to the time domain. It targets the complementary question: "What happened to those first stars?". It will be answered by observations of the most luminous events: supernovae and accretion on to black holes formed by direct collapse from the primordial gas clouds. These transients provide direct constraints on star-formation rates (SFRs) and the truly initial Initial Mass Function (IMF), and they may identify possible stellar seeds of supermassive black holes (SMBHs). Furthermore, our knowledge of the physics of these events at ultra-low metallicity will be much expanded. JWST's unique capabilities will detect these most luminous and earliest cosmic messengers easily in fairly shallow observations. However, these events are very rare at the dawn of cosmic structure formation and so require large area coverage. Time domain astronomy can be advanced to an unprecedented depth by means of a shallow field of JWST reaching 27 mag AB in 2 and 4.4 microns over a field as large as 0.1 square degree visited multiple times each year. Such a survey may set strong constraints or detect massive Pop III SNe at redshifts beyond 10, pinpointing the redshift of the first stars, or at least their death. Based on our current knowledge of superluminous supernovae (SLSNe), such a survey will find one or more SLSNe at redshifts above 6 in five years and possibly several direct collapse black holes. Although JWST is not designed as a wide field survey telescope, we show that such a wide field survey is possible with JWST and is critical in addressing several of its key scientific goals.

We investigate the star formation history and chemical evolution of isolated analogues of Local Group (LG) ultra faint dwarf galaxies (UFDs; stellar mass range of 10^2 solar mass < M_star <10^5 solar mass) and gas rich, low mass dwarfs (Leo P analogs; stellar mass range of 10^5 solar mass < M_star <10^6 solar mass). We perform a suite of cosmological hydrodynamic zoom-in simulations to follow their evolution from the era of the first generation of stars down to z=0. We confirm that reionization, combined with supernova (SN) feedback, is primarily responsible for the truncated star formation in UFDs. Specifically, haloes with a virial mass of M_vir < 2 x 10^9 solar mass form> 90\% of stars prior to reionization. Our work further demonstrates the importance of Population~III (Pop~III) stars, with their intrinsically high $\rm [C/Fe]$ yields, and the associated external metal-enrichment, in producing low-metallicity stars ($\rm [Fe/H]\lesssim-4$) and carbon-enhanced metal-poor (CEMP) stars. We find that UFDs are composite systems, assembled from multiple progenitor haloes, some of which hosted only Population~II (Pop~II) stars formed in environments externally enriched by SNe in neighboring haloes, naturally producing, extremely low-metallicity Pop~II stars. We illustrate how the simulated chemical enrichment may be used to constrain the star formation histories (SFHs) of true observed UFDs. We find that Leo P analogs can form in haloes with M_vir ~ 4 x 10^9 solar mass (z=0). Such systems are less affected by reionization and continue to form stars until z=0, causing higher metallicity tails. Finally, we predict the existence of extremely low-metallicity stars in LG UFD galaxies that preserve the pure chemical signatures of Pop~III nucleosynthesis.

The formation of the first stars in the high-redshift Universe is a sensitive probe of the small-scale, particle physics nature of dark matter (DM). We carry out cosmological simulations of primordial star formation in ultra-light, axion-like particle DM cosmology, with masses of $10^{-22}$ and $10^{-21}\,{\rm eV}$, with de Broglie wavelengths approaching galactic scales ($\sim$kpc). The onset of star formation is delayed, and shifted to more massive host structures. For the lightest DM particle mass explored here, first stars form at $z \sim 7$ in structures with $\sim 10^9\,{\rm M}_\odot$, compared to the standard minihalo environment within the $\Lambda$ cold dark matter ($\Lambda$CDM) cosmology, where $z \sim 20 - 30$ and $\sim 10^5 - 10^6\,{\rm M}_\odot$. Despite this greatly altered DM host environment, the thermodynamic behaviour of the metal-free gas as it collapses into the DM potential well asymptotically approaches a very similar evolutionary track. Thus, the fragmentation properties are predicted to remain the same as in $\Lambda$CDM cosmology, implying a similar mass scale for the first stars. These results predict intense starbursts in the axion cosmologies, which may be amenable to observations with the {\it James Webb Space Telescope}.

We use a semi-analytic model, ${\it Delphi}$, that jointly tracks the dark matter and baryonic assembly of high-redshift ($z \simeq 4-20$) galaxies to gain insight on the number density of Direct Collapse Black Hole (DCBH) hosts in three different cosmologies: the standard Cold Dark Matter (CDM) model and two Warm Dark Matter (WDM) models with particle masses of 3.5 and 1.5 keV. Obtaining the Lyman-Werner (LW) luminosity of each galaxy from ${\it Delphi}$, we use a clustering bias analysis to identify all, pristine halos with a virial temperature $T_{vir}>=10^4$ K that are irradiated by a LW background above a critical value as, DCBH hosts. In good agreement with previous studies, we find the DCBH number density rises from $\sim10^{-6.1}$ to $\sim 10^{-3.5}\, \mathrm{cMpc^{-3}}$ from $z\simeq 17.5$ to $8$ in the CDM model using a critical LW background value of $30 J_{21}$ (where $J_{21}= 10^{-21} \, {\rm erg\, s^{-1}\, Hz^{-1} \, cm^{-2} \, sr^{-1}}$). We find that a combination of delayed structure formation and an accelerated assembly of galaxies results in a later metal-enrichment and an accelerated build-up of the LW background in the 1.5 keV WDM model, resulting in DCBH hosts persisting down to much lower redshifts ($z \simeq 5$) as compared to CDM where DCBH hosts only exist down to $z \simeq 8$. We end by showing how the expected colours in three different bands of the Near Infrared Camera (NIRCam) onboard the forthcoming James Webb Space Telescope (${\it JWST}$) can be used to hunt for potential $z \simeq 5-9$ DCBHs, allowing hints on the WDM particle mass.

We perform a post-processing radiative feedback analysis on a 3D ab initio cosmological simulation of an atomic cooling halo under the direct collapse black hole (DCBH) scenario. We maintain the spatial resolution of the simulation by incorporating native ray-tracing on unstructured mesh data, including Monte Carlo Lyman-alpha (Ly{\alpha}) radiative transfer. DCBHs are born in gas-rich, metal-poor environments with the possibility of Compton-thick conditions, $N_H \gtrsim 10^{24} {\rm cm}^{-2}$. Therefore, the surrounding gas is capable of experiencing the full impact of the bottled-up radiation pressure. In particular, we find that multiple scattering of Ly{\alpha} photons provides an important source of mechanical feedback after the gas in the sub-parsec region becomes partially ionized, avoiding the bottleneck of destruction via the two-photon emission mechanism. We provide detailed discussion of the simulation environment, expansion of the ionization front, emission and escape of Ly{\alpha} radiation, and Compton scattering. A sink particle prescription allows us to extract approximate limits on the post-formation evolution of the radiative feedback. Fully coupled Ly{\alpha} radiation hydrodynamics will be crucial to consider in future DCBH simulations.

We briefly review the historical development of the ideas regarding the first supermassive black hole seeds, the physics of their formation and radiative feedback, recent theoretical and observational progress, and our outlook for the future.

Nascent neutron stars with millisecond periods and magnetic fields in excess of $10^{16}$ Gauss can drive highly energetic and asymmetric explosions known as magnetar-powered supernovae. These exotic explosions are one theoretical interpretation for supernovae Ic-BL which are sometimes associated with long gamma-ray bursts. Twisted magnetic field lines extract the rotational energy of the neutron star and release it as a disk wind or a jet with energies greater than 10$^{52}$ erg over $\sim 20$ sec. What fractions of the energy of the central engine go into the wind and the jet remain unclear. We have performed two-dimensional hydrodynamical simulations of magnetar-powered supernovae (SNe) driven by disk winds and jets with the CASTRO code to investigate the effect of the central engine on nucleosynthetic yields, mixing, and light curves. We find that these explosions synthesize less than 0.05 Msun of Ni and that this mass is not very sensitive to central engine type. The morphology of the explosion can provide a powerful diagnostic of the properties of the central engine. In the absence of a circumstellar medium these events are not very luminous, with peak bolometric magnitudes $M_b \sim -16.5 $ due to low Ni production.

Some ancient, dim, metal-poor stars may have formed in the ashes of the first supernovae (SNe). If their chemical abundances can be reconciled with the elemental yields of specific Population III (Pop III) explosions, they could reveal the properties of primordial stars. But multidimensional simulations of such explosions are required to predict their yields because dynamical instabilities can dredge material up from deep in the ejecta that would otherwise be predicted to fall back on to the central remnant and be lost in one-dimensional (1D) models. We have performed two-dimensional (2D) numerical simulations of two low-energy Pop III SNe, a 12.4 Msun explosion and a 60 Msun explosion, and find that they produce elemental yields that are a good fit to those measured in the most iron-poor star discovered to date, SMSS J031300.36-670839.3 (J031300). Fallback on to the compact remnant in these weak explosions accounts for the lack of measurable iron in J031300 and its low iron-group abundances in general. Our 2D explosions produce higher abundances of heavy elements (atomic number Z > 20) than their 1D counterparts due to dredge-up by fluid instabilities. Since almost no Ni is ejected by these weak SNe, their low luminosities will prevent their detection in the near-infrared with the James Webb Space Telescope and future 30-m telescopes on the ground. The only evidence that they ever occurred will be in the fossil abundance record.

The initial mass function of the first, Population III (Pop III), stars plays a vital role in shaping galaxy formation and evolution in the early Universe. One key remaining issue is the final fate of secondary protostars formed in the accretion disc, specifically whether they merge or survive. We perform a suite of hydrodynamic simulations of the complex interplay between fragmentation, protostellar accretion, and merging inside dark matter minihaloes. Instead of the traditional sink particle method, we employ a stiff equation of state approach, so that we can more robustly ascertain the viscous transport inside the disc. The simulations show inside-out fragmentation because the gas collapses faster in the central region. Fragments migrate on the viscous timescale, over which angular momentum is lost, enabling them to move towards the disc centre, where merging with the primary protostar can occur. This process depends on the fragmentation scale, such that there is a maximum scale of $(1 - 5) \times 10^4$ au, inside which fragments can migrate to the primary protostar. Viscous transport is active until radiative feedback from the primary protostar destroys the accretion disc. The final mass spectrum and multiplicity thus crucially depends on the effect of viscosity in the disc. The entire disc is subjected to efficient viscous transport in the primordial case with viscous parameter $\alpha \le 1$. An important aspect of this question is the survival probability of Pop III binary systems, possible gravitational wave sources to be probed with the Advanced LIGO detectors.

We compare model results from a semi-analytic (merger-tree based) framework for high-redshift (z ~ 5-20) galaxy formation against reionization indicators, including the Planck electron scattering optical depth and the ionizing photon emissivity, to shed light on the reionization history and sources in Cold (CDM) and Warm Dark Matter (WDM; particle masses of $m_x = 1.5,3$ and 5 keV) cosmologies. This model includes all the key processes of star formation, supernova feedback, the merger/accretion/ejection driven evolution of gas and stellar mass and the effect of the ultra-violet background (UVB) in photo-evaporating the gas content of low-mass galaxies. We find that the delay in the start of reionization in light (1.5 keV) WDM models can be compensated by a steeper redshift evolution of the ionizing photon escape fraction and a faster mass assembly, resulting in reionization ending at comparable redshifts (z~5.5) in all the DM models considered. We find the bulk of the reionization photons come from galaxies with a halo mass $M_h < 10^9 M_\odot$ and a UV magnitude $ -15 < M_{UV} < -10$ in CDM. The progressive suppression of low-mass halos with decreasing $m_x$ leads to a shift in the reionization population to larger halo masses of $M_h > 10^9 M_\odot$ and $ -17 < M_{UV} < -13$ for 1.5 keV WDM. We find that current observations of the electron scattering optical depth and the Ultra-violet luminosity function are equally compatible with all the (cold and warm) DM models considered in this work. We propose that global indicators including the redshift evolution of the stellar mass density and the stellar mass-halo mass relation, observable with the James Webb Space Telescope, can be used to distinguish between CDM and WDM (1.5 keV) cosmologies.

We study the contribution of the first galaxies to the far-infrared/sub-millimeter (FIR/sub-mm) extragalactic background light (EBL) by implementing an analytical model for dust emission. We explore different dust models, assuming different grain size distributions and chemical compositions. According to our findings, observed re-radiated emission from dust in dwarf-size galaxies at $z \sim 10$ would peak at a wavelength of $\sim 500 \mu {\rm m}$ with observed fluxes of $\sim 10^{-3} - 10^{-2}$ nJy, which is below the capabilities of current observatories. In order to be detectable, model sources at these high redshifts should exhibit luminosities of $\gtrsim 10^{12} L_{\odot}$, comparable to that of local ultra-luminous systems. The FIR/sub-mm EBL generated by primeval galaxies peaks at $\sim 500 \mu {\rm m}$, with an intensity ranging from $\sim 10^{-4}$ to $10^{-3} {\rm nW \ m^{-2} \ sr^{-1}}$, depending on dust properties. These values are $\sim 3 - 4$ orders of magnitude below the absolute measured cosmic background level, suggesting that the first galaxies would not contribute significantly to the observed FIR/sub-mm EBL. Our model EBL exhibits a strong correlation with the dust-to-metal ratio, where we assume a fiducial value of $D = 0.005$, increasing almost proportionally to it. Thus, measurements of the FIR/sub-mm EBL could provide constraints on the amount of dust in the early Universe. Even if the absolute signal from primeval dust emission may be undetectable, it might still be possible to obtain information about it by exploring angular fluctuations at $\sim 500 \mu {\rm m}$, close to the peak of dust emission from the first galaxies.

The dynamical impact of Lyman-alpha (Ly{\alpha}) radiation pressure on galaxy formation depends on the rate and duration of momentum transfer between Ly{\alpha} photons and neutral hydrogen gas. Although photon trapping has the potential to multiply the effective force, ionizing radiation from stellar sources may relieve the Ly{\alpha} pressure before appreciably affecting the kinematics of the host galaxy or efficiently coupling Ly{\alpha} photons to the outflow. We present self-consistent Ly{\alpha} radiation-hydrodynamics simulations of high-$z$ galaxy environments by coupling the Cosmic Ly{\alpha} Transfer code (COLT) with spherically symmetric Lagrangian frame hydrodynamics. The accurate but computationally expensive Monte-Carlo radiative transfer calculations are feasible under the one-dimensional approximation. The initial starburst drives an expanding shell of gas from the centre and in certain cases Ly{\alpha} feedback significantly enhances the shell velocity. Radiative feedback alone is capable of ejecting baryons into the intergalactic medium (IGM) for protogalaxies with a virial mass of $M_{\rm vir} \lesssim 10^8~{\rm M}_\odot$. We compare the Ly{\alpha} signatures of Population III stars with $10^5$ K blackbody emission to that of direct collapse black holes with a nonthermal Compton-thick spectrum and find substantial differences if the Ly{\alpha} spectra are shaped by gas pushed by Ly{\alpha} radiation-driven winds. For both sources, the flux emerging from the galaxy is reprocessed by the IGM such that the observed Ly{\alpha} luminosity is reduced significantly and the time-averaged velocity offset of the Ly{\alpha} peak is shifted redward.

Extremely metal poor stars have been the focus of much recent attention owing to the expectation that their chemical abundances can shed light on the metal and dust yields of the earliest supernovae. We present our most realistic simulation to date of the astrophysical pathway to the first metal enriched stars. We simulate the radiative and supernova hydrodynamic feedback of a $60\,M_\odot$ Population III star starting from cosmological initial conditions realizing Gaussian density fluctuations. We follow the gravitational hydrodynamics of the supernova remnant at high spatial resolution through its freely-expanding, adiabatic, and radiative phases, until gas, now metal-enriched, has resumed runaway gravitational collapse. Our findings are surprising: while the Population III progenitor exploded with a low energy of $10^{51}\,\text{erg}$ and injected an ample metal mass of $6\,M_\odot$, the first cloud to collapse after the supernova explosion is a dense surviving primordial cloud on which the supernova blastwave deposited metals only superficially, in a thin, unresolved layer. The first metal-enriched stars can form at a very low metallicity, of only $2-5\times10^{-4}\,Z_\odot$, and can inherit the parent cloud's highly elliptical, radially extended orbit in the dark matter gravitational potential.

CR7 is the brightest Lyman-$\alpha$ emitter observed at $z>6$, which shows very strong Lyman-$\alpha$ and HeII 1640\AA\ line luminosities, but no metal line emission. Previous studies suggest that CR7 hosts either young primordial stars with a total stellar mass of $\sim 10^7\,\mathrm{M}_\odot$ or a black hole of $\gtrsim 10^6\,\mathrm{M}_\odot$. Here, we explore different formation scenarios for CR7 with a semianalytical model, based on the random sampling of dark matter merger trees. We are unable to reproduce the observational constraints with a primordial stellar source, given our model assumptions, due to the short stellar lifetimes and the early metal enrichment. Black holes that are the remnants of the first stars are either not massive enough, or reside in metal-polluted haloes, ruling out this possible explanation of CR7. Our models instead suggest that direct collapse black holes, which form in metal-free haloes exposed to large Lyman-Werner fluxes, are more likely the origin of CR7. However, this result is derived under optimistic assumptions and future observations are necessary to further constrain the nature of CR7.

We simulate the growth of a Population III stellar system, starting from cosmological initial conditions at z=100. We follow the formation of a minihalo and the subsequent collapse of its central gas to high densities, resolving scales as small as ~ 1 AU. Using sink particles to represent the growing protostars, we model the growth of the photodissociating and ionizing region around the first sink, continuing the simulation for ~ 5000 yr after initial protostar formation. Along with the first-forming sink, several tens of secondary sinks form before an ionization front develops around the most massive star. The resulting cluster has high rates of sink formation, ejections from the stellar disc, and sink mergers during the first ~ 2000 yr, before the onset of radiative feedback. By this time a warm ~ 5000 K phase of neutral gas has expanded to roughly the disc radius of 2000 AU, slowing mass flow onto the disc and sinks. By 5000 yr the most massive star grows to 20 M_sol, while the total stellar mass approaches 75 M_sol. Out of the ~ 40 sinks, approximately 30 are low-mass (M_* < 1 M_sol), and if the simulation had resolved smaller scales an even greater number of sinks might have formed. Thus, protostellar radiative feedback is insufficient to prevent rapid disc fragmentation and the formation of a high-member Pop III cluster before an ionization front emerges. Throughout the simulation, the majority of stellar mass is contained within the most massive stars, further implying that the Pop III initial mass function is top-heavy.

Gravitational waves (GWs) provide a revolutionary tool to investigate yet unobserved astrophysical objects. Especially the first stars, which are believed to be more massive than present-day stars, might be indirectly observable via the merger of their compact remnants. We develop a self-consistent, cosmologically representative, semi-analytical model to simulate the formation of the first stars. By extrapolating binary stellar-evolution models at 10% solar metallicity to metal-free stars, we track the individual systems until the coalescence of the compact remnants. We estimate the contribution of primordial stars to the merger rate density and to the detection rate of the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO). Owing to their higher masses, the remnants of primordial stars produce strong GW signals, even if their contribution in number is relatively small. We find a probability of $\gtrsim1\%$ that the current detection GW150914 is of primordial origin. We estimate that aLIGO will detect roughly 1 primordial BH-BH merger per year for the final design sensitivity, although this rate depends sensitively on the primordial initial mass function (IMF). Turning this around, the detection of black hole mergers with a total binary mass of $\sim 300\,\mathrm{M}_\odot$ would enable us to constrain the primordial IMF.

Throughout the epoch of reionization the most luminous Ly{\alpha} emitters are capable of ionizing their own local bubbles. The CR7 galaxy at $z \approx 6.6$ stands out for its combination of exceptionally bright Ly{\alpha} and HeII 1640 Angstrom line emission but absence of metal lines. As a result CR7 may be the first viable candidate host of a young primordial starburst or direct collapse black hole. High-resolution spectroscopy reveals a +160 km s$^{-1}$ velocity offset between the Ly{\alpha} and HeII line peaks while the spatial extent of the Ly{\alpha} emitting region is $\sim 16$ kpc. The observables are indicative of an outflow signature produced by a strong central source. We present one-dimensional radiation-hydrodynamics simulations incorporating accurate Ly{\alpha} feedback and ionizing radiation to investigate the nature of the CR7 source. We find that a Population III star cluster with $10^5$ K blackbody emission ionizes its environment too efficiently to generate a significant velocity offset. However, a massive black hole with a nonthermal Compton-thick spectrum is able to reproduce the Ly{\alpha} signatures as a result of higher photon trapping and longer potential lifetime. For both sources, Ly{\alpha} radiation pressure turns out to be dynamically important.

Gamma-ray bursts (GRBs) are ideal probes of the epoch of the first stars and galaxies. We review the recent theoretical understanding of the formation and evolution of the first (so-called Population III) stars, in light of their viability of providing GRB progenitors. We proceed to discuss possible unique observational signatures of such bursts, based on the current formation scenario of long GRBs. These include signatures related to the prompt emission mechanism, as well as to the afterglow radiation, where the surrounding intergalactic medium might imprint a telltale absorption spectrum. We emphasize important remaining uncertainties in our emerging theoretical framework.

We explore the impact of a cosmic ray (CR) background generated by supernova explosions from the first stars on star-forming metal-free gas in a minihalo at $z\sim25$. Starting from cosmological initial conditions, we use the smoothed particle hydrodynamics code GADGET-2 to follow gas collapsing under the influence of a CR background up to densities of $n=10^{12}\,{\rm cm}^{-3}$, at which point we form sink particles. Using a suite of simulations with two sets of initial conditions and employing a range of CR background models, we follow each simulation for $5000\,$yr after the first sink forms. CRs both heat and ionise the gas, boosting ${\rm H}_2$ formation. Additional ${\rm H}_2$ enhances the cooling efficiency of the gas, allowing it to fulfil the Rees-Ostriker criterion sooner and expediting the collapse, such that each simulation reaches high densities at a different epoch. As it exits the loitering phase, the thermodynamic path of the collapsing gas converges to that seen in the absence of any CR background. By the time the gas approaches sink formation densities, the thermodynamic state of the gas is thus remarkably similar across all simulations. This leads to a robust characteristic mass that is largely independent of the CR background, of order $\sim$ a few $\times10\,{\rm M}_{\odot}$ even as the CR background strength varies by five orders of magnitude.

We model early star forming regions and their chemical enrichment by Population III (Pop III) supernovae with nucleosynthetic yields featuring high [C/Fe] ratios and pair-instability supernova (PISN) signatures. We aim to test how well these chemical abundance signatures are preserved in the gas prior to forming the first long-lived low-mass stars (or second-generation stars). Our results show that second-generation stars can retain the nucleosynthetic signature of their Pop III progenitors, even in the presence of nucleosynthetically normal Pop III core-collapse supernovae. We find that carbon-enhanced metal-poor stars are likely second-generation stars that form in minihaloes. Furthermore, it is likely that the majority of Pop III supernovae produce high [C/Fe] yields. In contrast, metals ejected by a PISN are not concentrated in the first star forming haloes, which may explain the absence of observed PISN signatures in metal-poor stars. We also find that unique Pop III abundance signatures in the gas are quickly wiped out by the emergence of Pop II supernovae. We caution that the observed fractions of stars with Pop III signatures cannot be directly interpreted as the fraction of Pop III stars producing that signature. Such interpretations require modelling the metal enrichment process prior to the second-generation stars' formation, including results from simulations of metal mixing. The full potential of stellar archaeology can likely be reached in ultra-faint dwarf galaxies, where the simple formation history may allow for straightforward identification of second-generation stars.

The study of the Cosmic Near-Infrared Background (CIB) light after subtraction of resolved sources can push the limits of current observations and infer the level of galaxy and black hole activity in the early universe. However, disentangling the relative contribution from low- and high-redshift sources is not trivial. Spatial fluctuations of the CIB exhibit a clustering excess at angular scales $\sim 1^\circ$ whose origin has not been conclusively identified. We explore the likelihood that this signal is dominated by emission from galaxies and accreting black holes in the early Universe. We find that, if the first small mass galaxies have a normal IMF, the light of their ageing stars (fossils) integrated over cosmic time contributes a comparable amount to the CIB as their pre-reionization progenitors. However, the measured fluctuation signal is too large to be produced by galaxies at redshifts $z>8$ unless their star formation efficiencies are much larger than those inferred from the observed Lyman-dropout population. In order to produce the observed level of CIB fluctuation without violating constraints from galaxy counts and the electron optical depth of the IGM, minihalos at $z>12$ must form stars with efficiency $f_\star \gtrsim 0.1$ and, although a top-heavy IMF is preferred, have a very low escape fraction of ionizing radiation, $f_{\rm esc}<0.01$. If instead the CIB fluctuations are produced by high-$z$ black holes, one requires vigorous accretion in the early universe reaching $\rho_{\rm acc} \gtrsim 10^5M_\odot{\rm Mpc^{-3}}$ by $z\simeq 10$. This growth must stop by $z \sim 6$ and be significantly obscured not to overproduce the soft cosmic X-ray background (CXB) and its observed coherence with the CIB. We therefore find the range of suitable possibilities at high-$z$ to be narrow, but could possibly be widened by including additional physics and evolution at those epochs.

We investigate the impact of a cosmic X-ray background (CXB) on Population III stars forming in a minihalo at $z\simeq25$. Using the smoothed particle hydrodynamics code GADGET-2, we attain sufficient numerical resolution to follow gas collapsing into the centre of the minihalo from cosmological initial conditions up to densities of $10^{12}\,{\rm cm}^{-3}$, at which point we form sink particles. This allows us to study how the presence of a CXB affects the formation of H$_2$ and HD in the gas prior to becoming fully molecular. Using a suite of simulations for a range of possible CXB models, we follow each simulation for 5000\yr after the first sink particle forms. The CXB provides two competing effects, with X-rays both heating the gas and increasing the free electron fraction, allowing more H$_2$ to form. X-ray heating dominates below $n\sim1\,{\rm cm}^{-3}$, while the additional H$_2$ cooling becomes more important above $n\sim10^2\,{\rm cm}^{-3}$. The gas becomes optically thick to X-rays as it exits the quasi-hydrostatic `loitering phase,' such that the primary impact of the CXB is to cool the gas at intermediate densities, resulting in an earlier onset of baryonic collapse into the dark matter halo. At the highest densities, self-shielding results in similar thermodynamic behaviour across a wide range of CXB strengths. Consequently, we find that star formation is relatively insensitive to the presence of a CXB; both the number and the characteristic mass of the stars formed remains quite similar even as the strength of the CXB varies by several orders of magnitude.

We investigate the formation of a galaxy reaching a virial mass of $~ 10^8$ solar mass at $z=10$ by carrying out a zoomed radiation-hydrodynamical cosmological simulation. This simulation traces Population~III (Pop~III) star formation, characterized by a modestly top-heavy initial mass function (IMF), and considers stellar feedback such as photoionization heating from Pop III and Population~II (Pop~II) stars, mechanical and chemical feedback from supernovae (SNe), and X-ray feedback from accreting black holes (BHs) and high-mass X-ray binaries (HMXBs). We self-consistently impose a transition in star formation mode from top-heavy Pop III to low-mass Pop~II, and find that the star formation rate in the computational box is dominated by Pop~III until $z=13$, and by Pop~II thereafter. The simulated galaxy experiences bursty star formation, with a substantially reduced gas content due to photoionization heating from Pop~III and Pop~II stars, together with SN feedback. All the gas within the simulated galaxy is metal-enriched above $10^{-5}$ solar, such that there are no remaining pockets of primordial gas. The simulated galaxy has an estimated observed flux of $~10^{-3} nJy$, which is too low to be detected by the James Webb Space Telescope (JWST) without strong lensing amplification. We also show that our simulated galaxy is similar in terms of stellar mass to Segue 2, the least luminous dwarf known in the Local Group.

We present a simulation of the long-term evolution of a Population III supernova remnant in a cosmological minihalo. Employing passive Lagrangian tracer particles, we investigate how chemical stratification and anisotropy in the explosion can affect the abundances of the first low-mass, metal-enriched stars. We find that reverse shock heating can leave the inner mass shells at entropies too high to cool, leading to carbon-enhancement in the re-collapsing gas. This hydrodynamic selection effect could explain the observed incidence of carbon-enhanced metal-poor (CEMP) stars at low metallicity. We further explore how anisotropic ejecta distributions, recently seen in direct numerical simulations of core-collapse explosions, may translate to abundances in metal-poor stars. We find that some of the observed scatter in the Population II abundance ratios can be explained by an incomplete mixing of supernova ejecta, even in the case of only one contributing enrichment event. We demonstrate that the customary hypothesis of fully-mixed ejecta clearly fails if post-explosion hydrodynamics prefers the recycling of some nucleosynthetic products over others. Furthermore, to fully exploit the stellar-archaeological program of constraining the Pop III initial mass function from the observed Pop II abundances, considering these hydrodynamical transport effects is crucial. We discuss applications to the rich chemical structure of ultra-faint dwarf satellite galaxies, to be probed in unprecedented detail with upcoming spectroscopic surveys.

We present a new near-field cosmological probe of the initial mass function (IMF) of the first stars. Specifically, we constrain the lower-mass limit of the Population III (Pop III) IMF with the total number of stars in large, unbiased surveys of the Milky Way. We model the early star formation history in a Milky Way-like halo with a semi-analytic approach, based on Monte-Carlo sampling of dark matter merger trees, combined with a treatment of the most important feedback mechanisms. Assuming a logarithmically flat Pop III IMF and varying its low mass limit, we derive the number of expected survivors of these first stars, using them to estimate the probability to detect any such Pop III fossil in stellar archaeological surveys. Following our analysis, the most promising region to find possible Pop III survivors is the stellar halo of the Milky Way, which is the best target for future surveys. We find that if no genuine Pop III survivor is detected in a sample size of $4 \times 10^6$ ($2 \times 10^7$) halo stars with well-controlled selection effects, then we can exclude the hypothesis that the primordial IMF extended down below $0.8 M_\odot$ at a confidence level of 68% (99%). With the sample size of the Hamburg/ESO survey, we can tentatively exclude Pop III stars with masses below $0.65 M_\odot$ with a confidence level of 95%, although this is subject to significant uncertainties. To fully harness the potential of our approach, future large surveys are needed that employ uniform, unbiased selection strategies for high-resolution spectroscopic follow-up.

It is widely recognized that nucleosynthetic output of the first, Population III supernovae was a catalyst defining the character of subsequent stellar generations. Most of the work on the earliest enrichment was carried out assuming that the first stars were extremely massive and that the associated supernovae were unusually energetic, enough to completely unbind the baryons in the host cosmic minihalo and disperse the synthesized metals into the intergalactic medium. Recent work, however, suggests that the first stars may in fact have been somewhat less massive, with a characteristic mass scale of a few tens of solar masses. We present a cosmological simulation following the transport of the metals synthesized in a Population III supernova assuming that it had an energy of 1e51 ergs, compatible with standard Type II supernovae. A young supernova remnant is inserted in the first star's relic HII region in the free expansion phase and is followed for 40 Myr employing adaptive mesh refinement and Lagrangian tracer particle techniques. The supernova remnant remains partially trapped within the minihalo and the thin snowplow shell develops pronounced instability and fingering. Roughly half of the ejecta turn around and fall back toward the center of the halo, with 1% of the ejecta reaching the center in 30 kyr and 10% in 10 Myr. The average metallicity of the combined returning ejecta and the pristine filaments feeding into the halo center from the cosmic web is 0.001 - 0.01 Z_sun, but the two remain unmixed until accreting onto the central hydrostatic core that is unresolved at the end of the simulation. We conclude that if Population III stars had less extreme masses, they promptly enriched the host minihalos with metals and triggered Population II star formation.

To constrain the properties of the first stars with the chemical abundance patterns observed in metal-poor stars, one must identify any non-trivial effects that the hydrodynamics of metal dispersal can imprint on the abundances. We use realistic cosmological hydrodynamic simulations to quantify the distribution of metals resulting from one Population III supernova and from a small number of such supernovae exploding in close succession. Overall, supernova ejecta are highly inhomogeneously dispersed throughout the simulations. When the supernova bubbles collapse, quasi-virialized metal-enriched clouds, fed by fallback from the bubbles and by streaming of metal-free gas from the cosmic web, grow in the centers of the dark matter halos. Partial turbulent homogenization on scales resolved in the simulation is observed only in the densest clouds where the vortical time scales are short enough to ensure true homogenization on subgrid scales. However, the abundances in the clouds differ from the gross yields of the supernovae. Continuing the simulations until the cloud have gone into gravitational collapse, we predict that the abundances in second-generation stars will be deficient in the innermost mass shells of the supernova (if only one has exploded) or in the ejecta of the latest supernovae (when multiple have exploded). This indicates that hydrodynamics gives rise to biases complicating the identification of nucleosynthetic sources in the chemical abundance spaces of the surviving stars.

We present new methodology to use cosmic infrared background (CIB) fluctuations to probe sources at 10<z<30 from a JWST/NIRCam configuration that will isolate known galaxies to 28 AB mag at 0.5--5 micron. At present significant mutually consistent source-subtracted CIB fluctuations have been identified in the Spitzer and Akari data at 2--5 micron, but we demonstrate internal inconsistencies at shorter wavelengths in the recent CIBER data. We evaluate CIB contributions from remaining galaxies and show that the bulk of the high-z sources will be in the confusion noise of the NIRCam beam, requiring CIB studies. The accurate measurement of the angular spectrum of the fluctuations and probing the dependence of its clustering component on the remaining shot noise power would discriminate between the various currently proposed models for their origin and probe the flux distribution of its sources. We show that the contribution to CIB fluctuations from remaining galaxies is large at visible wavelengths for the current instruments precluding probing the putative Lyman-break of the CIB fluctuations. We demonstrate that with the proposed JWST configuration such measurements will enable probing the Lyman break. We develop a Lyman-break tomography method to use the NIRCam wavelength coverage to identify or constrain, via the adjacent two-band subtraction, the history of emissions over 10<z<30 as the Universe comes out of the 'Dark Ages'. We apply the proposed tomography to the current Spitzer/IRAC measurements at 3.6 and 4.5 micron, to find that it already leads to interestingly low upper limit on emissions at z>30.

We present the Cosmic Lyman-$\alpha$ Transfer code (COLT), a massively parallel Monte-Carlo radiative transfer code, to simulate Lyman-$\alpha$ (Ly$\alpha$) resonant scattering through neutral hydrogen as a probe of the first galaxies. We explore the interaction of centrally produced Ly$\alpha$ radiation with the host galactic environment. Ly$\alpha$ photons emitted from the luminous starburst region escape with characteristic features in the line profile depending on the density distribution, ionization structure, and bulk velocity fields. For example, anisotropic ionization exhibits a tall peak close to line centre with a skewed tail that drops off gradually. Idealized models of first galaxies explore the effect of mass, anisotropic H II regions, and radiation pressure driven winds on Ly$\alpha$ observables. We employ mesh refinement to resolve critical structures. We also post-process an ab initio cosmological simulation and examine images captured at various escape distances within the 1 Mpc$^3$ comoving volume. Finally, we discuss the emergent spectra and surface brightness profiles of these objects in the context of high-$z$ observations. The first galaxies will likely be observed through the red damping wing of the Ly$\alpha$ line. Observations will be biased toward galaxies with an intrinsic red peak located far from line centre that reside in extensive H II super bubbles, which allows Hubble flow to sufficiently redshift photons away from line centre and facilitate transmission through the intergalactic medium (IGM). Even with gravitational lensing to boost the luminosity this preliminary work indicates that Ly$\alpha$ emission from stellar clusters within haloes of $M_{\rm vir}<10^9~{\rm M}_\odot$ is generally too faint to be detected by the James Webb Space Telescope (JWST).

We present the results of the stellar feedback from Pop III binaries by employing improved, more realistic Pop III evolutionary stellar models. To facilitate a meaningful comparison, we consider a fixed mass of 60 solar masses (Msun) incorporated in Pop III stars, either contained in a single star, or split up in binary stars of 30 Msun each or an asymmetric case of one 45 Msun and one 15 Msun star. Whereas the sizes of the resulting HII regions are comparable across all cases, the HeIII regions around binary stars are significantly smaller than that of the single star. Consequently, the He$^{+}$ 1640 angstrom recombination line is expected to become much weaker. Supernova feedback exhibits great variety due to the uncertainty in possible explosion pathways. If at least one of the component stars dies as a hypernova about ten times more energetic than conventional core-collapse supernovae, the gas inside the host minihalo is effectively blown out, chemically enriching the intergalactic medium (IGM) to an average metallicity of $10^{-4}-10^{-3}$ solar metallicity (Zsun), out to $\sim 2$ kpc. The single star, however, is more likely to collapse into a black hole, accompanied by at most very weak explosions. The effectiveness of early chemical enrichment would thus be significantly reduced, in difference from the lower mass binary stars, where at least one component is likely to contribute to heavy element production and dispersal. Important new feedback physics is also introduced if close binaries can form high-mass x-ray binaries, leading to the pre-heating and -ionization of the IGM beyond the extent of the stellar HII regions.

We simulate the formation of a low metallicity (0.01 Zsun) stellar cluster in a dwarf galaxy at redshift z~14. Beginning with cosmological initial conditions, the simulation utilizes adaptive mesh refinement and sink particles to follow the collapse and evolution of gas past the opacity limit for fragmentation, thus resolving the formation of individual protostellar cores. A time- and location-dependent protostellar radiation field, which heats the gas by absorption on dust, is computed by integration of protostellar evolutionary tracks with the MESA code. The simulation also includes a robust non-equilibrium chemical network that self-consistently treats gas thermodynamics and dust-gas coupling. The system is evolved for 18 kyr after the first protostellar source has formed. In this time span, 30 sink particles representing protostellar cores form with a total mass of 81 Msun. Their masses range from ~0.1 Msun to 14.4 Msun with a median mass ~0.5-1 Msun. Massive protostars grow by competitive accretion while lower-mass protostars are stunted in growth by close encounters and many-body ejections. In the regime explored here, the characteristic mass scale is determined by the temperature floor set by the cosmic microwave background and by the onset of efficient dust-gas coupling. It seems unlikely that host galaxies of the first bursts of metal-enriched star formation will be detectable with the James Webb Space Telescope or other next-generation infrared observatories. Instead, the most promising access route to the dawn of cosmic star formation may lie in the scrutiny of metal-poor, ancient stellar populations in the Galactic neighborhood. The observable targets that correspond to the system simulated here are ultra-faint dwarf satellite galaxies such as Bootes II, Segue I and II, and Willman I.

We use cosmological simulations to assess how the explosion of the first stars in supernovae (SNe) influences early cosmic history. Specifically, we investigate the impact by SNe on the host systems for Population~III (Pop~III) star formation and explore its dependence on halo environment and Pop~III progenitor mass. We then trace the evolution of the enriched gas until conditions are met to trigger second-generation star formation. To this extent, we quantify the recovery timescale, which measures the time delay between a Pop~III SN explosion and the appearance of cold, dense gas, out of which second-generation stars can form. We find that this timescale is highly sensitive to the Pop~III progenitor mass, and less so to the halo environment. For more massive progenitors, including those exploding in pair instability SNe, second-generation star formation is delayed significantly, for up to a Hubble time. The dependence of the recovery time on the mass of the SN progenitor is mainly due to the ionizing impact of the progenitor star. Photoionization heating increases the gas pressure and initiates a hydrodynamical response that reduces the central gas density, an effect that is stronger in more massive. The gas around lower mass Pop~III stars remains denser and hence the SN remnants cool more rapidly, facilitating the subsequent re-condensation of the gas and formation of a second generation of stars. In most cases, the second-generation stars are already metal-enriched to ~2-5 X 10^{-4}\zsun, thus belonging to Population~II. The recovery timescale is a key quantity governing the nature of the first galaxies, able to host low-mass, long-lived stellar systems. These in turn are the target of future deep-field campaigns with the James Webb Space Telescope.

We use cosmological simulations of high-redshift minihalos to investigate the effect of dark matter annihilation (DMA) on the collapse of primordial gas. We numerically investigate the evolution of the gas as it assembles in a Population III stellar disk. We find that when DMA effects are neglected, the disk undergoes multiple fragmentation events beginning at ~ 500 yr after the appearance of the first protostar. On the other hand, DMA heating and ionization of the gas speeds the initial collapse of gas to protostellar densities and also affects the stability of the developing disk against fragmentation, depending on the DM distribution. We compare the evolution when we model the DM density with an analytical DM profile which remains centrally peaked, and when we simulate the DM profile using N-body particles (the 'live' DM halo). When utilizing the analytical DM profile, DMA suppresses disk fragmentation for ~ 3500 yr after the first protostar forms, in agreement with earlier work. However, when using a 'live' DM halo, the central DM density peak is gradually flattened due to the mutual interaction between the DM and the rotating gaseous disk, reducing the effects of DMA on the gas, and enabling secondary protostars of mass ~ 1 M_sol to be formed within ~ 900 yr. These simulations demonstrate that DMA is ineffective in suppressing gas collapse and subsequent fragmentation, rendering the formation of long-lived dark stars unlikely. However, DMA effects may still be significant in the early collapse and disk formation phase of primordial gas evolution.

Recent work suggests that the first generation of stars, the so-called Population III (Pop III), could have formed primarily in binaries or as members of small multiple systems. Here we investigate the impact of X-ray feedback from High-Mass X-ray Binaries (HMXBs) left behind in stellar binary systems after the primary forms a black hole (BH), accreting gas at a high rate from the companion, a process that is thought to be favored at the low metallicities characteristic of high-redshift gas. Thanks to their large mean free path, X-rays are capable of preionizing and preheating the gas in the intergalactic medium (IGM) and in haloes long before the reionization of the Universe is complete, and thus could have strongly affected the formation of subsequent generations of stars as well as reionization. We have carried out zoomed hydrodynamical cosmological simulations of minihaloes, accounting for the formation of Pop III stars and their collapse into BHs and HMXBs, and the associated radiation-hydrodynamic feedback from UV and X-ray photons. We find no strong net feedback from HMXBs on the simulated star formation history. On the other hand, the preheating of the IGM by HMXBs leads to a strong suppression of small-scale structures and significantly lowers the recombination rate in the IGM, thus yielding a net positive feedback on reionization. We further show that X-ray feedback from HMXBs can augment the ionizing feedback from the Pop III progenitor stars to suppress gas accretion onto the first BHs, limiting their growth into supermassive BHs. Finally, we show that X-ray ionization by HMXBs leaves distinct signatures in the properties of the high-redshift hydrogen that may be probed in upcoming observations of the redshifted 21cm spin-flip line.

We perform numerical simulations of the growth of a Population III stellar system under photodissociating feedback. We start from cosmological initial conditions at z = 100, self-consistently following the formation of a minihalo at z = 15 and the subsequent collapse of its central gas to high densities. The simulations resolve scales as small as ~ 1 AU, corresponding to gas densities of 10^16 cm^-3. Using sink particles to represent the growing protostars, we evolve the stellar system for the next 5000 years. We find that this emerging stellar group accretes at an unusually low rate compared with minihalos which form at earlier times (z = 20 - 30), or with lower baryonic angular momentum. The stars in this unusual system will likely reach masses ranging from < 1 M_sun to 5 M_sun by the end of their main-sequence lifetimes, placing them in the mass range for which stars will undergo an asymptotic giant branch (AGB) phase. Based upon the simulation, we predict the rare existence of Population III stars that have survived to the present day and have been enriched by mass overflow from a previous AGB companion.

Dwarf spheroidal satellite galaxies have emerged a powerful probe of small-scale dark matter clustering and of cosmic reionization. They exhibit structural and chemical continuity with dwarf irregular galaxies in the field and with spheroidal galaxies in high-density environments. By combining empirical constraints derived for star formation at low gas column densities and metallicities in the local universe with a model for dark matter and baryonic mass assembly, we provide an analytical description of how the dwarf spheroidals acquired their stellar content. Their progenitors formed stars until the gas content, initially reduced from the cosmic average by the thermal pressure of the reionized intergalactic medium, was finally ram pressure stripped during the progenitors' accretion on to the host galaxy. Dwarf spheroidal satellites of differing luminosities seem to share very similar most massive progenitor histories that reach thresholds for gas cooling by atomic line emission at epochs at which the Lagrangian volume of the Local Group should have been reionized. We hypothesize that dwarf spheroidals formed the bulk of their stars in partially rotationally supported HI disks in a reionized universe. This model provides an explanation for the "common mass scale" relation and reproduces the empirical luminosity-size and luminosity-metallicity relations. Explosive feedback phenomena, such as outflows driven by the concerted action of supernovae, need not have been significant in the dwarf spheroidals' formation. We further speculate that the true pre-reionization fossils should exhibit a structure distinct from that of the dwarf spheroidals, e.g., in the form of dense isolated or nuclear star clusters.

We simulate the formation of a metal-poor (10^-2 Zsun) stellar cluster in one of the first galaxies to form in the early Universe, specifically a high-redshift atomic cooling halo (z~14). This is the first calculation that resolves the formation of individual metal-enriched stars in simulations starting from realistic cosmological initial conditions. We follow the evolution of a single dense clump among several in the parent halo. The clump forms a cluster of ~40 stars and sub-stellar objects within 7000 years and could continue forming stars ~5 times longer. Protostellar dust heating has a negligible effect on the star formation efficiency, at least during the early evolutionary stages, but it moderately suppresses gaseous fragmentation and brown dwarf formation. We observe fragmentation in thin gaseous filaments and sustained accretion in larger, rotating structures as well as ejections by binary interactions. The stellar initial mass function above 0.1 Msun, evaluated after ~10^4 years of fragmentation and accretion, seems in agreement with the recent measurement in ultra-faint dwarf spheroidal Galactic satellites of Geha et al. (2013).

We investigate the impact of dust-induced gas fragmentation on the formation of the first low-mass, metal-poor stars (< 1Msun) in the early universe. Previous work has shown the existence of a critical dust-to-gas ratio, below which dust thermal cooling cannot cause gas fragmentation. Assuming the first dust is silicon-based, we compute critical dust-to-gas ratios and associated critical silicon abundances ([Si/H]crit). At the density and temperature associated with protostellar disks, we find that a standard Milky Way grain size distribution gives [Si/H]crit = -4.5 +/- 0.1, while smaller grain sizes created in a supernova reverse shock give [Si/H]crit = -5.3 +/- 0.1. Other environments are not dense enough to be influenced by dust cooling. We test the silicate dust cooling theory by comparing to silicon abundances observed in the most iron-poor stars ([Fe/H] < -4.0). Several stars have silicon abundances low enough to rule out dust-induced gas fragmentation with a standard grain size distribution. Moreover, two of these stars have such low silicon abundances that even dust with a shocked grain size distribution cannot explain their formation. Adding small amounts of carbon dust does not significantly change these conclusions. Additionally, we find that these stars exhibit either high carbon with low silicon abundances or the reverse. A silicate dust scenario thus suggests that the earliest low-mass star formation in the most metal-poor regime may have proceeded through two distinct cooling pathways: fine structure line cooling and dust cooling. This naturally explains both the carbon-rich and carbon-normal stars at extremely low [Fe/H].

Events such as GRB130606A at z=5.91, offer an exciting new window into pre-galactic metal enrichment in these very high redshift host galaxies. We study the environment and host galaxy of GRB 130606A, a high-z event, in the context of a high redshift population of GRBs. We have obtained multiwavelength observations from radio to gamma-ray, concentrating particularly on the X-ray evolution as well as the optical photometric and spectroscopic data analysis. With an initial Lorentz bulk factor in the range \Gamma_0 ~ 65-220, the X-ray afterglow evolution can be explained by a time-dependent photoionization of the local circumburst medium, within a compact and dense environment. The host galaxy is a sub-DLA (log N (HI) = 19.85+/-0.15), with a metallicity content in the range from ~1/7 to ~1/60 of solar. Highly ionized species (N V and Si IV) are also detected. This is the second highest redshift burst with a measured GRB-DLA metallicity and only the third GRB absorber with sub-DLA HI column density. GRB ' lighthouses' therefore offer enormous potential as backlighting sources to probe the ionization and metal enrichment state of the IGM at very high redshifts for the chemical signature of the first generation of stars.

Population III stars are believed to have been more massive than typical stars today and to have formed in relative isolation. The thermodynamic impact of metals is expected to induce a transition leading to clustered, low-mass Population II star formation. In this work, we present results from three cosmological simulations, only differing in gas metallicity, that focus on the impact of metal fine-structure line cooling on the formation of stellar clusters in a high-redshift atomic cooling halo. Introduction of sink particles allows us to follow the process of gas hydrodynamics and accretion onto cluster stars for 4 Myr corresponding to multiple local free-fall times. At metallicities at least $10^{-3}\, Z_{\odot}$, gas is able to reach the CMB temperature floor and fragment pervasively resulting in a stellar cluster of size $\sim1$ pc and total mass $\sim1000\, M_{\odot}$. The masses of individual sink particles vary, but are typically $\sim100\, M_{\odot}$, consistent with the Jeans mass when gas cools to the CMB temperature, though some solar mass fragments are also produced. At the low metallicity of $10^{-4}\, Z_{\odot}$, fragmentation is completely suppressed on scales greater than 0.01 pc and total stellar mass is lower by a factor of 3 than in the higher metallicity simulations. The sink particle accretion rates, and thus their masses, are determined by the mass of the gravitationally unstable gas cloud and the prolonged gas accretion over many Myr. The simulations thus exhibit features of both monolithic collapse and competitive accretion. Even considering possible dust induced fragmentation that would occur at higher densities, the formation of a bona fide stellar cluster seems to require metal line cooling and metallicities of at least $10^{-3}\, Z_{\odot}$.

Understanding the formation of the first stars is one of the frontier topics in modern astrophysics and cosmology. Their emergence signaled the end of the cosmic dark ages, a few hundred million years after the Big Bang, leading to a fundamental transformation of the early Universe through the production of ionizing photons and the initial enrichment with heavy chemical elements. We here review the state of our knowledge, separating the well understood elements of our emerging picture from those where more work is required. Primordial star formation is unique in that its initial conditions can be directly inferred from the Lambda Cold Dark Matter (LCDM) model of cosmological structure formation. Combined with gas cooling that is mediated via molecular hydrogen, one can robustly identify the regions of primordial star formation, the so-called minihalos, having total masses of ~10^6 M_sun and collapsing at redshifts z~20-30. Within this framework, a number of studies have defined a preliminary standard model, with the main result that the first stars were predominantly massive. This model has recently been modified to include a ubiquitous mode of fragmentation in the protostellar disks, such that the typical outcome of primordial star formation may be the formation of a binary or small multiple stellar system. We will also discuss extensions to this standard picture due to the presence of dynamically significant magnetic fields, of heating from self-annihalating WIMP dark matter, or cosmic rays. We conclude by discussing possible strategies to empirically test our theoretical models.

We investigate the operation of the chemothermal instability in primordial star-forming clouds with a suite of three-dimensional, moving-mesh simulations. In line with previous studies, we find that the gas at the centre of high-redshift minihaloes becomes chemothermally unstable as three-body reactions convert the atomic hydrogen into a fully molecular gas. The competition between the increasing rate at which the gas cools and the increasing optical depth to H2 line emission creates a characteristic dip in the cooling time over the free-fall time on a scale of 100 au. As a result, the free-fall time decreases to below the sound-crossing time, and the cloud may become gravitationally unstable and fragment on a scale of a few tens of au during the initial free-fall phase. In three of the nine haloes investigated, secondary clumps condense out of the parent cloud, which will likely collapse in their own right before they are accreted by the primary clump. In the other haloes, fragmentation at such an early stage is less likely. However, given that previous simulations have shown that the infall velocity decreases substantially once the gas becomes rotationally supported, the amount of time available for perturbations to develop may be much greater than is evident from the limited period of time simulated here.

We perform a cosmological simulation in order to model the growth and evolution of Population III (Pop III) stellar systems in a range of host minihalo environments. A Pop III multiple system forms in each of the ten minihaloes, and the overall mass function is top-heavy compared to the currently observed initial mass function in the Milky Way. Using a sink particle to represent each growing protostar, we examine the binary characteristics of the multiple systems, resolving orbits on scales as small as 20 AU. We find a binary fraction of ~36%, with semi-major axes as large as 3000 AU. The distribution of orbital periods is slightly peaked at < 900 yr, while the distribution of mass ratios is relatively flat. Of all sink particles formed within the ten minihaloes, ~50% are lost to mergers with larger sinks, and ~50% of the remaining sinks are ejected from their star-forming disks. The large binary fraction may have important implications for Pop III evolution and nucleosynthesis, as well as the final fate of the first stars.

We investigate how radiative feedback from the first stars affects the assembly of the first dwarf galaxies. We perform cosmological zoomed SPH simulations of a dwarf galaxy assembling inside a halo of virial mass 10^9 solar at z = 10. The simulations follow the non-equilibrium chemistry/cooling of primordial gas and the conversion of the gas into metal-free stars. To quantify the radiative feedback, we compare a simulation in which stars emit both molecular hydrogen dissociating and hydrogen/helium ionizing radiation with a simulation in which stars emit only dissociating radiation, and with a simulation in which stars remain dark. Photodissociation and -ionization exert a strong negative feedback on the assembly of the simulated galaxy. Gas condensation is strongly impeded, and star formation is strongly suppressed in comparison with the simulation in which stars remain dark. The feedback on the gas implies a suppression of the central dark matter densities in the minihalo progenitor by factors of up to a few, which is a significant deviation from the singular isothermal density profile characterizing the dark matter distribution in the absence of radiative feedback. The evolution of gas densities, star formation rates, and the distribution of dark matter becomes insensitive to the inclusion of dissociating radiation in the late stages of the minihalo assembly, and it becomes insensitive to the inclusion of ionizing radiation once the minihalo turns into an atomically cooling galaxy. The formation of an extended disk inside the dwarf galaxy is a robust outcome not affected by the inclusion of radiation. We estimate that dwarf galaxies such as simulated here will be among the faintest galaxies the upcoming James Webb Space Telescope will detect. Our conclusions are subject to our neglect of feedback from supernovae and chemical enrichment as well as to cosmic variance. [abridged]

We analyze the cosmological simulations performed in the recent work of Greif et al. (2012), which followed the early growth and merger history of Pop III stars while resolving scales as small as 0.05 R_sol. This is the first set of cosmological simulations to self-consistently resolve the rotation and internal structure of Pop III protostars. We find that Pop III stars form under significant rotational support which is maintained for the duration of the simulations. The protostellar surfaces spin from ~50% to nearly 100% of Keplerian rotational velocity. These rotation rates persist after experiencing multiple stellar merger events. In the brief time period simulated (~ 10 yr), the protostars show little indication of convective instability, and their properties furthermore show little correlation with the properties of their host minihaloes. If Pop III protostars within this range of environments generally form with high degrees of rotational support, and if this rotational support is maintained for a sufficient amount of time, this has a number of crucial implications for Pop III evolution and nucleosynthesis, as well as the possibility for Pop III pair-instability supernovae, and the question of whether the first stars produced gamma-ray bursts.

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Pop III GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly-resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature >10^4 K. We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1+z) with values of about a few cm^{-3}. In more massive halos, corresponding to the first galaxies, the density may be larger, n>100 cm^{-3}. The resulting afterglow fluxes may be detectable with the JWST and VLA in the near-IR and radio wavebands, respectively, out to redshift z>20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages. We show that if the afterglow can be observed sufficiently early, analysis of the metal lines can distinguish whether the first heavy elements were produced in a PISN, or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

The emergence of the first sources of light at redshifts of z ~ 10-30 signaled the transition from the simple initial state of the Universe to one of increasing complexity. We review recent progress in our understanding of the formation of the first stars and galaxies, starting with cosmological initial conditions, primordial gas cooling, and subsequent collapse and fragmentation. We emphasize the important open question of how the pristine gas was enriched with heavy chemical elements in the wake of the first supernovae. We conclude by discussing how the chemical abundance patterns conceivably allow us to probe the properties of the first stars and subsequent stellar generations, and allow us to test models of early metal enrichment.

We utilize metal-poor stars in the local, ultra-faint dwarf galaxies (UFDs; L_tot < 10^5 L_sun) to empirically constrain the formation process of the first galaxies. Since UFDs have much simpler star formation histories than the halo of the Milky Way, their stellar populations should preserve the fossil record of the first supernova (SN) explosions in their long-lived, low-mass stars. Guided by recent hydrodynamical simulations of first galaxy formation, we develop a set of stellar abundance signatures that characterize the nucleosynthetic history of such an early system if it was observed in the present-day universe. Specifically, we argue that the first galaxies are the product of chemical "one-shot" events, where only one (long-lived) stellar generation forms after the first, Population III, SN explosions. Our abundance criteria thus constrain the strength of negative feedback effects inside the first galaxies. We compare the stellar content of UFDs with these one-shot criteria. Several systems (Ursa Major II, and also Coma Berenices, Bootes I, Leo IV, Segue 1) largely fulfill the requirements, indicating that their high-redshift predecessors did experience strong feedback effects that shut off star formation. We term the study of the entire stellar population of a dwarf galaxy for the purpose of inferring details about the nature and origin of the first galaxies "dwarf galaxy archaeology". This will provide clues to the connection of the first galaxies, the surviving, metal-poor dwarf galaxies, and the building blocks of the Milky Way.

We investigate the process of metal-free star formation in the first galaxies with a high-resolution cosmological simulation. We consider the cosmologically motivated scenario in which a strong molecule-destroying Lyman-Werner (LW) background inhibits effective cooling in low-mass haloes, delaying star formation until the collapse or more massive haloes. Only when molecular hydrogen (H2) can self-shield from LW radiation, which requires a halo capable of cooling by atomic line emission, will star formation be possible. To follow the formation of multiple gravitationally bound objects, at high gas densities we introduce sink particles which accrete gas directly from the computational grid. We find that in a 1 Mpc^3 (comoving) box, runaway collapse first occurs in a 3x10^7 M_sun dark matter halo at z~12 assuming a background intensity of J21=100. Due to a runaway increase in the H2 abundance and cooling rate, a self-shielding, supersonically turbulent core develops abruptly with ~10^4 M_sun in cold gas available for star formation. We analyze the formation of this self-shielding core, the character of turbulence, and the prospects for star formation. Due to a lack of fragmentation on scales we resolve, we argue that LW-delayed metal-free star formation in atomic cooling haloes is very similar to star formation in primordial minihaloes, although in making this conclusion we ignore internal stellar feedback. Finally, we briefly discuss the detectability of metal-free stellar clusters with the James Webb Space Telescope.

Theoretical models predict that some of the first stars ended their lives as extremely energetic pair-instability supernovae (PISNe). With energies approaching 10^53 ergs, these supernovae are expected to be within the detection limits of the upcoming James Webb Space telescope (JWST), allowing observational constraints to be placed on the properties of the first stars. We estimate the source density of PISNe using a semi-analytic halo mass function based approach, accounting for the effects of feedback from star formation on the PISN rates using cosmological simulations. We estimate an upper limit of ~0.2 PISNe per JWST field of view at any given time. Feedback can reduce this rate significantly, e.g., lowering it to as little as one PISN per 4000 JWST fields of view for the most pessimistic explosion models. We also find that the main obstacle to observing PISNe from the first stars is their scarcity, not their faintness; exposures longer than a few times 10^4s will do little to increase the number of PISNe found. Given this we suggest a mosaic style search strategy for detecting PISNe from the first stars. Even rather high redshift PISNe are unlikely to be missed by moderate exposures, and a large number of pointings will be required to ensure a detection.

We investigate the formation of the first stars at the end of the cosmic dark ages with a suite of three-dimensional, moving mesh simulations that directly resolve the collapse of the gas beyond the formation of the first protostar at the centre of a dark matter minihalo. The simulations cover more than 25 orders of magnitude in density and have a maximum spatial resolution of 0.05 R_sun, which extends well below the radius of individual protostars and captures their interaction with the surrounding gas. In analogy to previous studies that employed sink particles, we find that the Keplerian disc around the primary protostar fragments into a number of secondary protostars, which is facilitated by H2 collisional dissociation cooling and collision-induced emission. The further evolution of the protostellar system is characterized by strong gravitational torques that transfer angular momentum between the secondary protostars formed in the disc and the surrounding gas. This leads to the migration of about half of the secondary protostars to the centre of the cloud in a free-fall time, where they merge with the primary protostar and enhance its growth to about five times the mass of the second most massive protostar. By the same token, a fraction of the protostars obtain angular momentum from other protostars via N-body interactions and migrate to higher orbits. On average, only every third protostar survives until the end of the simulation. However, the number of protostars present at any given time increases monotonically, suggesting that the system will continue to grow beyond the limited period of time simulated here.

We study how the first galaxies were assembled under feedback from the accretion onto a central black hole (BH) that is left behind by the first generation of metal-free stars through self-consistent, cosmological simulations. X-ray radiation from the accretion of gas onto BH remnants of Population III (Pop III) stars, or from high-mass X-ray binaries (HMXBs), again involving Pop III stars, influences the mode of second generation star formation. We track the evolution of the black hole accretion rate and the associated X-ray feedback starting with the death of the Pop III progenitor star inside a minihalo and following the subsequent evolution of the black hole as the minihalo grows to become an atomically cooling galaxy. We find that X-ray photoionization heating from a stellar-mass BH is able to quench further star formation in the host halo at all times before the halo enters the atomic cooling phase. X-ray radiation from a HMXB, assuming a luminosity close to the Eddington value, exerts an even stronger, and more diverse, feedback on star formation. It photoheats the gas inside the host halo, but also promotes the formation of molecular hydrogen and cooling of gas in the intergalactic medium and in nearby minihalos, leading to a net increase in the number of stars formed at early times. Our simulations further show that the radiative feedback from the first BHs may strongly suppress early BH growth, thus constraining models for the formation of supermassive BHs.

Understanding the formation of the first stars and galaxies is a key problem in modern cosmology. In these lecture notes, we will derive some of the basic physical principles underlying this emerging field. We will consider the basic cosmological context, the cooling and chemistry in primordial gas, the physics of gravitational instability, and the main properties of the first stars. We will conclude with a discussion of the observational signature of the first sources of light, to be probed with future telescopes, such as the James Webb Space Telescope.

We perform three-dimensional cosmological simulations to examine the growth of metal-free, Population III (Pop III) stars under radiative feedback. We begin our simulation at z=100 and trace the evolution of gas and dark matter until the formation of the first minihalo. We then follow the collapse of the gas within the minihalo up to densities of n = 10^12 cm^-3, at which point we replace the high-density particles with a sink particle to represent the growing protostar. We model the effect of Lyman-Werner (LW) radiation emitted by the protostar, and employ a ray-tracing scheme to follow the growth of the surrounding H II region over the next 5000 yr. We find that a disk assembles around the first protostar, and that radiative feedback will not prevent further fragmentation of the disk to form multiple Pop III stars. Ionization of neutral hydrogen and photodissociation of H_2 by LW radiation leads to heating of the dense gas to several thousand Kelvin, and this warm region expands outward at the gas sound speed. Once the extent of this warm region becomes equivalent to the size of the disk, the disk mass declines while the accretion rate onto the protostars is reduced by an order of magnitude. This occurs when the largest sink has grown to ~ 20 M_sol while the second sink has grown to 7 M_sol, and we estimate the main sink will approach an asymptotic value of ~ 30 M_sol by the time it reaches the main sequence. Our simulation thus indicates that the most likely outcome is a massive Pop III binary. However, we simulate only one minihalo, and the statistical variation between minihaloes may be substantial. If Pop III stars were typically unable to grow to more than a few tens of solar masses, this would have important consequences for the occurence of pair-instability supernovae in the early Universe as well as the Pop III chemical signature in the oldest stars observable today.

We numerically study the mutual interaction between dark matter (DM) and Population III (Pop III) stellar systems in order to explore the possibility of Pop III dark stars within this physical scenario. We perform a cosmological simulation, initialized at z ~ 100, which follows the evolution of gas and DM. We analyze the formation of the first minihalo at z ~ 20 and the subsequent collapse of the gas to densities of 10^12 cm^-3. We then use this simulation to initialize a set of smaller-scale `cut-out' simulations in which we further refine the DM to have spatial resolution similar to that of the gas. We test multiple DM density profiles, and we employ the sink particle method to represent the accreting star-forming region. We find that, for a range of DM configurations, the motion of the Pop III star-disk system serves to separate the positions of the protostars with respect to the DM density peak, such that there is insufficient DM to influence the formation and evolution of the protostars for more than ~ 5000 years. In addition, the star-disk system causes gravitational scattering of the central DM to lower densities, further decreasing the influence of DM over time. Any DM-powered phase of Pop III stars will thus be very short-lived for the typical multiple system, and DM will not serve to significantly prolong the life of Pop III stars.

The cosmic dark ages ended a few hundred million years after the Big Bang, when the first stars began to fill the universe with new light. It has generally been argued that these stars formed in isolation and were extremely massive - perhaps 100 times as massive as the Sun. In a recent study, Clark and collaborators showed that this picture requires revision. They demonstrated that the accretion disks that build up around Population III stars are strongly susceptible to fragmentation and that the first stars should therefore form in clusters rather than in isolation. We here use a series of high-resolution hydrodynamical simulations performed with the moving mesh code AREPO to follow up on this proposal and to study the influence of environmental parameters on the level of fragmentation. We model the collapse of five independent minihalos from cosmological initial conditions, through the runaway condensation of their central gas clouds, to the formation of the first protostar, and beyond for a further 1000 years. During this latter accretion phase, we represent the optically thick regions of protostars by sink particles. Gas accumulates rapidly in the circumstellar disk around the first protostar, fragmenting vigorously to produce a small group of protostars. After an initial burst, gravitational instability recurs periodically, forming additional protostars with masses ranging from ~ 0.1 to 10 M_sun. Although the shape, multiplicity, and normalization of the protostellar mass function depend on the details of the sink-particle algorithm, fragmentation into protostars with diverse masses occurs in all cases, confirming earlier reports of Population III stars forming in clusters. Depending on the efficiency of later accretion and merging, Population III stars may enter the main sequence in clusters and with much more diverse masses than are commonly assumed.

We review our current understanding of how the first galaxies formed at the end of the cosmic dark ages, a few 100 million years after the Big Bang. Modern large telescopes discovered galaxies at redshifts greater than seven, whereas theoretical studies have just reached the degree of sophistication necessary to make meaningful predictions. A crucial ingredient is the feedback exerted by the first generation of stars, through UV radiation, supernova blast waves, and chemical enrichment. The key goal is to derive the signature of the first galaxies to be observed with upcoming or planned next-generation facilities, such as the James Webb Space Telescope or Atacama Large Millimeter Array. From the observational side, ongoing deep-field searches for very high-redshift galaxies begin to provide us with empirical constraints on the nature of the first galaxies.

The James Webb Space Telescope (JWST) will enable observations of galaxies at redshifts z > 10 and hence allow to test our current understanding of structure formation at very early times. Previous work has shown that the very first galaxies inside halos with virial temperatures T < 10^4 K and masses M < 10^8 M_sun at z > 10 are probably too faint, by at least one order of magnitude, to be detected even in deep exposures with JWST. The light collected with JWST may therefore be dominated by radiation from galaxies inside ten times more massive halos. We use cosmological zoomed smoothed particle hydrodynamics simulations to investigate the assembly of such galaxies and assess their observability with JWST. We compare two simulations that are identical except for the inclusion of non-equilibrium H/D chemistry and radiative cooling by molecular hydrogen. In both simulations a large fraction of the halo gas settles in two nested, extended gas disks which surround a compact massive gas core. The presence of molecular hydrogen allows the disk gas to reach low temperatures and to develop marked spiral structure but does not qualitatively change its stability against fragmentation. We post-process the simulated galaxies by combining idealized models for star formation with stellar population synthesis models to estimate the luminosities in nebular recombination lines as well as in the ultraviolet continuum. We demonstrate that JWST will be able to constrain the nature of the stellar populations in galaxies such as simulated here based on the detection of the He1640 recombination line. Extrapolation of our results to halos with masses both lower and higher than those simulated shows that JWST may find up to a thousand star-bursting galaxies in future deep exposures of the z > 10 universe.

We evaluate the effect of a supersonic relative velocity between the baryons and dark matter on the thermal and density evolution of the first gas clouds at z < 50. Through a series of cosmological simulations, initialized at z=100 with a range of relative streaming velocities and minihalo formation redshifts, we find that the typical streaming velocities will have little effect on the gas evolution. Once the collapse begins, the subsequent evolution of the gas will be nearly indistinguishable from the case of no streaming, and star formation will still proceed in the same way, with no change in the characteristic Pop III stellar masses. Reionization is expected to be dominated by halo masses of > 10^8 M_sun, for which the expected effect of streaming is negligible.

The very first stars to form in the Universe heralded an end to the cosmic dark ages and introduced new physical processes that shaped early cosmic evolution. Until now, it was thought that these stars lived short, solitary lives, with only one extremely massive star, or possibly a very wide binary system, forming in each dark matter minihalo. Here we describe numerical simulations that show that these stars were, to the contrary, often members of tight multiple systems. Our results show that the disks that formed around the first young stars were unstable to gravitational fragmentation, possibly producing small binary and higher-order systems that had separations as small as the distance between the Earth and the Sun.

We report results from numerical simulations of star formation in the early universe that focus on the dynamical behavior of metal-free gas under different initial and environmental conditions. In particular we investigate the role of turbulence, which is thought to ubiquitously accompany the collapse of high-redshift halos. We distinguish between two main cases: the birth of Population III.1 stars - those which form in the pristine halos unaffected by prior star formation - and the formation of Population III.2 stars - those forming in halos where the gas is still metal free but has an increased ionization fraction. This latter case can arise either from exposure to the intense UV radiation of stellar sources in neighboring halos, or from the high virial temperatures associated with the formation of massive halos, that is, those with masses greater than 1e8 solar masses. We find that turbulent primordial gas is highly susceptible to fragmentation in both cases, even for turbulence in the subsonic regime, i.e. for rms velocity dispersions as low as 20 % of the sound speed. Contrary to our original expectations, fragmentation is more vigorous and more widespread in pristine halos compared to pre-ionized ones. We therefore predict Pop III.1 stars to be on average of somewhat lower mass, and form in larger groups, than Pop III.2 stars. We find that fragment masses cover over two orders of magnitude, indicating that the resulting Population III initial mass function was significantly extended in mass as well. This prompts the need for a large, high-resolution study of the formation of dark matter minihalos that is capable of resolving the turbulent flows in the gas at the moment when the baryons become self-gravitating. This would help determine which, if any, of the initial conditions presented in our study are realized in nature.

We estimate the rotation speed of Population III (Pop III) stars within a minihalo at z ~ 20 using a smoothed particle hydrodynamics (SPH) simulation, beginning from cosmological initial conditions. We follow the evolution of the primordial gas up to densities of 10^12 cm^-3. Representing the growing hydrostatic cores with accreting sink particles, we measure the velocities and angular momenta of all particles that fall onto these protostellar regions. This allows us to record the angular momentum of the sinks and estimate the rotational velocity of the Pop III stars expected to form within them. The rotation rate has important implications for the evolution of the star, the fate encountered at the end of its life, and the potential for triggering a gamma-ray burst (GRB). We find that there is sufficient angular momentum to yield rapidly rotating stars (> 1000 km s^-1, or near break-up speeds). This indicates that Pop III stars likely experienced strong rotational mixing, impacting their structure and nucleosynthetic yields. A subset of them was also likely to result in hypernova explosions, and possibly GRBs.

The transformation of atomic hydrogen to molecular hydrogen through three-body reactions is a crucial stage in the collapse of primordial, metal-free halos, where the first generation of stars (Population III stars) in the Universe are formed. However, in the published literature, the rate coefficient for this reaction is uncertain by nearly an order of magnitude. We report on the results of both adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of the collapse of metal-free halos as a function of the value of this rate coefficient. For each simulation method, we have simulated a single halo three times, using three different values of the rate coefficient. We find that while variation between halo realizations may be greater than that caused by the three-body rate coefficient being used, both the accretion physics onto Population III protostars as well as the long-term stability of the disk and any potential fragmentation may depend strongly on this rate coefficient.

We study the gravitational fragmentation of cold accretion streams flowing into a typical first galaxy. We use a one-zone hydrodynamical model to examine the thermal evolution of the gas entering a 10^8 M_sun DM halo at z=10. The goal is to find the expected fragmentation mass scale and thus a characteristic mass of the first population of stars to form by shock fragmentation at high redshift. Our model accurately describes the chemical and thermal evolution of the gas as we are specifically concerned with how the cooling of the gas alters its fragmentation properties. We find there to be a sharp drop in the fragmentation mass at a metallicity of ~10^-4 Z_sun when a strong molecule destroying, LW background is present. However, If molecules can efficiently form, they dominate the cooling at T < 10^4 K, demonstrating no 'critical metallicity'. Dust grains are not included in our chemical model, but we argue their inclusion would not significantly the results. We also find that this physical scenario allows for the formation of a cluster of solar mass fragments, or a single 10^4 M_sun fragment, possibly the precursors to primeval clusters and SMBHs. Lastly, we conclude that the usual assumption of isobaricity for galactic shocks breaks down in gas of sufficiently high metallicity, suggesting that metal cooling may lead to thermal instabilities.

We report results from numerical simulations of star formation in the early universe that focus on the role of subsonic turbulence, and investigate whether it can induce fragmentation of the gas. We find that dense primordial gas is highly susceptible to fragmentation, even for rms turbulent velocity dispersions as low as 20% of the initial sound speed. The resulting fragments cover over two orders of magnitude in mass, ranging from 0.1 to 40 solar masses. However, our results suggest that the details of the fragmentation depend on the local properties of the turbulent velocity field and hence we expect considerable variations in the resulting stellar mass spectrum in different halos.

Using three-dimensional cosmological simulations, we study the assembly process of one of the first galaxies, with a total mass of 10^8 M_sun, collapsing at z = 10. Our main goal is to trace the transport of the heavy chemical elements produced and dispersed by a pair-instability supernova exploding in one of the minihalo progenitors. To this extent, we incorporate an efficient algorithm into our smoothed particle hydrodynamics code which approximately models turbulent mixing as a diffusion process. We study this mixing with and without the radiative feedback from Population III stars that subsequently form in neighboring minihalos. Our simulations allow us to constrain the initial conditions for second-generation star formation, within the first galaxy itself, and inside of minihalos that virialize after the supernova explosion. We find that most minihalos remain unscathed by ionizing radiation or the supernova remnant, while some are substantially photoheated and enriched to supercritical levels, likely resulting in the formation of low-mass Population III or even Population II stars. At the center of the newly formed galaxy, 10^5 M_sun of cold, dense gas uniformly enriched to 10^-3 Z_sun are in a state of collapse, suggesting that a cluster of Population II stars will form. The first galaxies, as may be detected by the James Webb Space Telescope, would therefore already contain stellar populations familiar from lower redshifts.

We investigate the formation of metal-free, Population III (Pop III), stars within a minihalo at z ~ 20 with a smoothed particle hydrodynamics (SPH) simulation, starting from cosmological initial conditions. Employing a hierarchical, zoom-in procedure, we achieve sufficient numerical resolution to follow the collapsing gas in the center of the minihalo up to number densities of 10^12 cm^-3. This allows us to study the protostellar accretion onto the initial hydrostatic core, which we represent as a growing sink particle, in improved physical detail. The accretion process, and in particular its termination, governs the final masses that were reached by the first stars. The primordial initial mass function (IMF), in turn, played an important role in determining to what extent the first stars drove early cosmic evolution. We continue our simulation for 5000 yr after the first sink particle has formed. During this time period, a disk-like configuration is assembled around the first protostar. The disk is gravitationally unstable, develops a pronounced spiral structure, and fragments into several other protostellar seeds. At the end of the simulation, a small multiple system has formed, dominated by a binary with masses ~ 40 M_Sun and ~ 10 M_Sun. If Pop III stars were to form typically in binaries or small multiples, the standard model of primordial star formation, where single, isolated stars are predicted to form in minihaloes, would have to be modified. This would have crucial consequences for the observational signature of the first stars, such as their nucleosynthetic pattern, and the gravitational-wave emission from possible Pop III black-hole binaries.

We use three-dimensional smoothed particle hydrodynamics simulations together with a dynamical ray-tracing scheme to investigate the build-up of the first H ii regions around massive Population III stars in minihaloes. We trace the highly anisotropic breakout of the ionising radiation into the intergalactic medium, allowing us to predict the resulting recombination radiation with greatly increased realism. Our simulations, together with Press-Schechter type arguments, allow us to predict the Population III contribution to the radio background at 100 MHz via bremsstrahlung and 21 cm emission. We find a global bremsstrahlung signal of around 1 mK, and a combined 21 cm signature which is an order of magnitude larger. Both might be within reach of the planned Square Kilometer Array experiment, although detection of the free-free emission is only marginal. The imprint of the first stars on the cosmic radio background might provide us with one of the few diagnostics to test the otherwise elusive minihalo star formation site.

The formation of the first galaxies at redshifts z~10-15 signaled the transition from the simple initial state of the universe to one of ever increasing complexity. We here review recent progress in understanding their assembly process with numerical simulations, starting with cosmological initial conditions and modelling the detailed physics of star formation. In particular, we study the role of HD cooling in ionized primordial gas, the impact of UV radiation produced by the first stars, and the propagation of the supernova blast waves triggered at the end of their brief lives. We conclude by discussing promising observational diagnostics that will allow us to probe the properties of the first galaxies, such as their contribution to reionization and the chemical abundance pattern observed in extremely low-metallicity stars.

We investigate the properties of the first galaxies at z > 10 with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds 10^4 K, and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a 5*10^7 M_sun system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures > 10^4 K catalyses the formation of H2 and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock-heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation.

The formation of the first galaxies at redshifts z ~ 10-15 signaled the transition from the simple initial state of the universe to one of ever increasing complexity. We here review recent progress in understanding their assembly process with numerical simulations, starting with cosmological initial conditions and modelling the detailed physics of star formation. In this context we emphasize the importance and influence of selecting appropriate initial conditions for the star formation process. We revisit the notion of a critical metallicity resulting in the transition from primordial to present-day initial mass functions and highlight its dependence on additional cooling mechanisms and the exact initial conditions. We also review recent work on the ability of dust cooling to provide the transition to present-day low-mass star formation. In particular, we highlight the extreme conditions under which this transition mechanism occurs, with violent fragmentation in dense gas resulting in tightly packed clusters.

We perform three-dimensional smoothed particle hydrodynamics simulations in a realistic cosmological setting to investigate the expansion, feedback, and chemical enrichment properties of a 200 M_sun pair-instability supernova in the high-redshift universe. We find that the SN remnant propagates for a Hubble time at z = 20 to a final mass-weighted mean shock radius of 2.5 kpc (proper), roughly half the size of the HII region, and in this process sweeps up a total gas mass of 2.5*10^5 M_sun. The morphology of the shock becomes highly anisotropic once it leaves the host halo and encounters filaments and neighboring minihalos, while the bulk of the shock propagates into the voids of the intergalactic medium. The SN entirely disrupts the host halo and terminates further star formation for at least 200 Myr, while in our specific case it exerts positive mechanical feedback on neighboring minihalos by shock-compressing their cores. In contrast, we do not observe secondary star formation in the dense shell via gravitational fragmentation, due to the previous photoheating by the progenitor star. We find that cooling by metal lines is unimportant for the entire evolution of the SN remnant, while the metal-enriched, interior bubble expands adiabatically into the cavities created by the shock, and ultimately into the voids with a maximum extent similar to the final mass-weighted mean shock radius. Finally, we conclude that dark matter halos of at least M_vir > 10^8 M_sun must be assembled to recollect all components of the swept-up gas.

We construct star formation histories at redshifts z > 5 for two physically distinct populations of primordial, metal-free stars, motivated by theoretical and observational arguments that have hinted towards the existence of an intermediate stellar generation between Population III and Population I/II. Taking into account the cosmological parameters as recently revised by WMAP after three years of operation, we determine self-consistent reionization histories and discuss the resulting chemical enrichment from these early stellar generations. We find that the bulk of ionizing photons and heavy elements produced at high redshifts must have originated in Population II.5 stars, which formed out of primordial gas in haloes with virial temperatures >= 10^4 K, and had typical masses >= 10 M_sun. Classical Population III stars, formed in minihaloes and having masses >= 100 M_sun, on the other hand, had only a minor impact on reionization and early metal enrichment. Specifically, we find that only 10 per cent by mass of metal-free star formation went into Population III.

Detection of the radiation emitted from the first galaxies at z > 10 will be made possible in the next decade, with the launch of the James Webb Space Telescope (JWST). We carry out cosmological radiation hydrodynamics simulations of Population III (Pop III) starbursts in a 10^8 M_Sun dwarf galaxy at z = 12.5. For different star formation efficiencies and stellar initial mass functions (IMFs), we calculate the luminosities and equivalent widths (EWs) of the recombination lines H_alpha, Ly_alpha, and He II 1640, under the simple assumption that the stellar population does not evolve over the first ~3 Myr of the starburst. Although only < 40 percent of the gas in the central 100 pc of the galaxy is photoionized, we find that photoheating by massive stars causes a strong dynamical response, which results in a weak correlation between luminosity emitted in hydrogen recombination lines and the total mass in stars. However, owing to the low escape fraction of He II-ionizing photons, the luminosity emitted in He II 1640 is much more strongly correlated with the total stellar mass. The ratio of the luminosity in He II 1640 to that in Ly_alpha or H_alpha is found to be a good indicator of the IMF in many cases. The ratio of observable fluxes is F_1640/F_Halpha ~ 1 for clusters of 100 M_Sun Pop III stars and F_1640/F_Halpha ~ 0.1 for clusters of 25 M_Sun Pop III stars. The EW of the He II 1640 emission line is the most reliable IMF indicator, its value varying between ~ 20 and ~ 200 angstrom for a massive and very massive Pop III IMF, respectively. Even the bright, initial stages of Pop III starbursts in the first dwarf galaxies will likely not be directly detectable by the JWST. Instead, the JWST may discover only more massive, and hence more chemically evolved, galaxies which host normal, Pop I/II, star formation.

Observations made using large ground-based and space-borne telescopes have probed cosmic history all the way from the present-day to a time when the Universe was less than a tenth of its present age. Earlier on lies the remaining frontier, where the first stars, galaxies, and massive black holes formed. They fundamentally transformed the early Universe by endowing it with the first sources of light and chemical elements beyond the primordial hydrogen and helium produced in the Big Bang. The interplay of theory and upcoming observations promises to answer the key open questions in this emerging field.

Probing the growth of structure from the epoch of hydrogen recombination to the formation of the first stars and galaxies is one of the most important uncharted areas of observational cosmology. Far-IR spectroscopy covering $\lambda$ 100-500 microns from space, and narrow partial transmission atmospheric bands available from the ground, opens up the possibility of probing the molecular hydrogen and metal fine-structure lines from primordial clouds from which the first stars and galaxies formed at 6 < z $<$ 15. Building on Spitzer observations of unexpectedly powerful H2 emission from shocks, we argue that next-generation far-IR space telescopes may open a new window into the main cloud cooling processes and feedback effects which characterized this vital, but unexplored epoch. Without this window, we are essential blind to the dominant cloud cooling which inevitably led to star formation and cosmic reionization.

The cosmic dark ages are the mysterious epoch during which the pristine gas began to condense and ultimately form the first stars. Although these beginnings have long been a topic of theoretical interest, technology has only recently allowed the beginnings of observational insight into this epoch. Many questions surround the formation of stars in metal-free gas and the history of the build-up of metals in the intergalactic medium: (1) What were the properties of the first stellar and galactic sources to form in pristine (metal-free) gas? (2) When did the epoch of Population III (metal-free) star formation take place and how long did it last? (3) Was the stellar initial mass function dramatically different for the first stars and galaxies? These questions are all active areas of theoretical research. However, new observational constraints via the direct detection of Population III star formation are vital to making progress in answering the broader questions surrounding how galaxies formed and how the cosmological properties of the universe have affected the objects it contains.

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

We introduce a simple and efficient algorithm for diffusion in smoothed particle hydrodynamics (SPH) simulations and apply it to the problem of chemical mixing. Based on the concept of turbulent diffusion, we link the diffusivity of a pollutant to the local physical conditions and can thus resolve mixing in space and time. We apply our prescription to the evolution of an idealized supernova remnant and find that we can model the distribution of heavy elements without having to explicitly resolve hydrodynamic instabilities in the post-shock gas. Instead, the dispersal of the pollutant is implicitly modelled through its dependence on the local velocity dispersion. Our method can thus be used in any SPH simulation that investigates chemical mixing but lacks the necessary resolution on small scales. Potential applications include the enrichment of the interstellar medium in present-day galaxies, as well as the intergalactic medium at high redshifts.

The validity of the hypothesis that the massive black holes in high redshift quasars grew from stellar-sized "seeds" is contingent on a seed's ability to double its mass every few ten million years. This requires that the seed accrete at approximately the Eddington-limited rate. In the specific case of radiatively efficient quasiradial accretion in a metal-poor protogalactic medium, for which the Bondi accretion rate is often prescribed in cosmological simulations of massive black hole formation, we examine the effects of the radiation emitted near the black hole's event horizon on the structure of the surrounding gas flow. We find that the radiation pressure from photoionization significantly reduces the steady-state accretion rate and renders the quasiradial accretion flow unsteady and inefficient. The time-averaged accretion rates are a small fraction of the Eddington-limited accretion rate for Thomson scattering. The pressure of Ly-alpha photons trapped near the HII region surrounding the black hole may further attenuate the inflow. These results suggest that an alternative to quasiradial, radiatively efficient Bondi-like accretion should be sought to explain the rapid growth of quasar-progenitor seed black holes.

We present the first results from two-dimensional simulations of radiatively-efficient accretion of metal-free gas onto intermediate-mass black holes. We fix the shape of the spectral energy distribution of the radiation produced near the event horizon and study the structure of the irradiated low-angular-momentum accretion flow over three orders of magnitude in radius from the black hole, 10^{14}-10^{17} cm for a 100 M_sun black hole. The luminosity of the central source is made to be proportional to the rate at which gas accretes across the inner boundary, which we set just inside the sonic radius. We find that photoionization heating and radiation pressure modify the structure of the flow. When the ambient gas density is 10^7 cm^{-3}, accretion is intermittent and on average reduced to 32% of the Eddington-limited rate, two orders of magnitude below the "Bondi" rate evaluated ignoring radiation, in agreement with simplified theoretical models. Even if the vicinity of the black hole is supplied with high density gas, accretion is rendered inefficient through heating and radiation pressure.

The first stars fundamentally transformed the early Universe through their production of energetic radiation and the first heavy chemical elements. The impact on cosmic evolution sensitively depends on their initial mass function (IMF), which can be empirically constrained through detailed studies of ancient, metal-poor halo stars in our Galaxy. We compare the lowest magnesium and iron abundances measured in Galactic halo stars with theoretical predictions for the minimum stellar enrichment provided by Population III stars under the assumption of a top-heavy IMF. To demonstrate that abundances measured in metal-poor stars reflect the chemical conditions at their formation, and that they can thus be used to derive constraints on the primordial IMF, we carry out a detailed kinematic analysis of a large sample of metal-poor stars drawn from the SDSS survey. We assess whether interstellar accretion has altered their surface abundances. We find that accretion is generally negligible, even at the extremely low levels where the primordial IMF can be tested. We conclude that the majority of the first stars were very massive, but had likely masses below ~140 M.

The first stars were key drivers of early cosmic evolution. We review the main physical elements of the current consensus view, positing that the first stars were predominantly very massive. We continue with a discussion of important open questions that confront the standard model. Among them are uncertainties in the atomic and molecular physics of the hydrogen and helium gas, the multiplicity of stars that form in minihalos, and the possible existence of two separate modes of metal-free star formation.

The character of the first galaxies at redshifts z > 10 strongly depends on their level of pre-enrichment, which is in turn determined by the rate of primordial star formation prior to their assembly. In order for the first galaxies to remain metal-free, star formation in minihaloes must be highly suppressed, most likely by H2-dissociating Lyman-Werner (LW) radiation. We show that the build-up of such a strong LW background is hindered by two effects. Firstly, the level of the LW background is self-regulated, being produced by the Population III (Pop III) star formation which it, in turn, suppresses. Secondly, the high opacity to LW photons which is built up in the relic H II regions left by the first stars acts to diminish the global LW background. Accounting for a self-regulated LW background, we estimate a lower limit for the rate of Pop III star formation in minihaloes at z > 15. Further, we simulate the formation of a 'first galaxy' with virial temperature T > 10^4 K and total mass > 10^8 M_Sun at z > 10, and find that complete suppression of previous Pop III star formation is unlikely, with stars of > 100 M_Sun (Pop III.1) and > 10 M_Sun (Pop III.2) likely forming. Finally, we discuss the implications of these results for the nature of the first galaxies, which may be observed by future missions such as the James Webb Space Telescope.

The relationship between galaxies and supermassive black holes (SMBH) found in their cores plays a key role in the formation and evolution of both of these major constituents of the universe, as well as the evolution of the intergalactic medium. Neither can be fully understood on their own, and studies of galaxy and SMBH co-formation and co-evolution are now among the central topics of research in cosmology. Yet the very origins, and the early growth phases of the SMBH are still not firmly established. We review our current understanding of the relevant processes and their astrophysical and cosmological context, with an emphasis on the observability of the SMBH growth mechanisms at high redshifts, and their leftover progeny at low redshifts.

The formation of the first generations of stars at redshifts z > 15-20 signaled the transition from the simple initial state of the universe to one of increasing complexity. We here review recent progress in understanding the assembly process of the first galaxies, starting with cosmological initial conditions and modelling the detailed physics of star formation. In particular, we study the role of HD cooling in ionized primordial gas, the impact of UV radiation produced by the first stars, and the propagation of the supernova blast waves triggered at the end of their brief lives. We conclude by discussing how the chemical abundance patterns observed in extremely low-metallicity stars allow us to probe the properties of the first stars.

The chemical abundance patterns observed in metal-poor Galactic halo stars contain the signature of the first supernovae, and thus allows us to probe the first stars that formed in the universe. We construct a theoretical model for the early chemical enrichment history of the Milky Way, aiming in particular at the contribution from pair-instability supernovae (PISNe). These are a natural consequence of current theoretical models for primordial star formation at the highest masses. However, no metal-poor star displaying the distinct PISN signature has yet been observed. We here argue that this apparent absence of any PISN signature is due to an observational selection effect. Whereas most surveys traditionally focus on the most metal-poor stars, we predict that early PISN enrichment tends to `overshoot', reaching enrichment levels of [Ca/H] ~= -2.5 that would be missed by current searches. We utilize existing observational data to place constraints on the primordial initial mass function (IMF). The number fraction of PISNe in the primordial stellar population is estimated to be <0.07, or ~< 40% by mass, assuming that metal-free stars have masses in excess of 10 M_sun. We further predict, based on theoretical estimates for the relative number of PISNe, that the expected fraction of second generation stars below [Ca/H]=-2 with a dominant (i.e., >90%) contribution from PISNe is merely ~1e-4 - 5e-4. The corresponding fraction of stars formed from gas exclusively enriched by PISNe is a factor of ~4 smaller. With the advent of next generation telescopes and new, deeper surveys, we should be able to test these predictions.

We describe our visualization process for a particle-based simulation of the formation of the first stars and their impact on cosmic history. The dataset consists of several hundred time-steps of point simulation data, with each time-step containing approximately two million point particles. For each time-step, we interpolate the point data onto a regular grid using a method taken from the radiance estimate of photon mapping. We import the resulting regular grid representation into ParaView, with which we extract isosurfaces across multiple variables. Our images provide insights into the evolution of the early universe, tracing the cosmic transition from an initially homogeneous state to one of increasing complexity. Specifically, our visualizations capture the build-up of regions of ionized gas around the first stars, their evolution, and their complex interactions with the surrounding matter. These observations will guide the upcoming James Webb Space Telescope, the key astronomy mission of the next decade.

Current observations are about to open up a direct window into the final frontier of cosmology: the first billion years in cosmic history when the first stars and galaxies formed. Even before the launch of the James Webb Space Telescope, it might be possible to utilize Gamma-ray Bursts (GRBs) as unique probes of cosmic star formation and the state of the intergalactic medium (IGM) up to redshifts of several tens, when the first (Population III) stars had formed. The Swift mission, or future satellites such as EXIST, might be the first observatories to detect individual Population III stars, provided that massive metal-free stars were able to trigger GRBs. Spectroscopic follow-up observations of the GRB afterglow emission would allow to probe the ionization state and metal enrichment of the IGM as a function of redshift.

We explore the implications of a possible cosmic ray (CR) background generated during the first supernova explosions that end the brief lives of massive Population III stars. We show that such a CR background could have significantly influenced the cooling and collapse of primordial gas clouds in minihaloes around redshifts of z ~ 15 - 20, provided the CR flux was sufficient to yield an ionization rate greater than about 10^-19 s^-1 near the center of the minihalo. The presence of CRs with energies less than approximately 10^7 eV would indirectly enhance the molecular cooling in these regions, and we estimate that the resulting lower temperatures in these minihaloes would yield a characteristic stellar mass as low as ~ 10 M_sun. CRs have a less pronounced effect on the cooling and collapse of primordial gas clouds inside more massive dark matter haloes with virial masses greater than approximately 10^8 M_sun at the later stages of cosmological structure formation around z ~ 10 - 15. In these clouds, even without CR flux the molecular abundance is already sufficient to allow cooling to the floor set by the temperature of the cosmic microwave background.

We investigate the conditions under which the first low-mass stars formed in the universe by confronting theoretical predictions governing the transition from massive Population III to low-mass Population II stars with recent observational C and/or O abundance data of metal-poor Galactic stars. Based on an existing theory, we introduce a new "observer-friendly" function, the transition discriminant D_trans, which provides empirical constraints as well as a powerful comparison between the currently available data of metal-poor halo stars and theoretical predictions of the formation of the first low-mass stars (<1 M_sun). Specifically, we compare the empirical stellar results with the theory that fine-structure lines of C and O dominate the transition from Pop III to Pop II in the early universe. We find the observational data for halo objects as well as for dSph galaxies and globular clusters to be consistent with this theory. An explanation for the observed lack of metal-poor stars in dSph galaxies and globular clusters is also suggested. Finally, we predict that any star to be found with [Fe/H]<-4 should have enhanced C and/or O abundances. The high C and O abundances of the two most iron-poor stars are in line with our prediction.

The formation of the first galaxies is influenced by the radiative feedback from the first generations of stars. This feedback is manisfested by the heating and ionization of the gas which lies within the H II regions surrounding the first stars, as well as by the photodissociation of hydrogen molecules within the larger Lyman-Werner (LW) bubbles that surround these sources. Using a ray-tracing method in three-dimensional cosmological simulations, we self-consistently track the formation of, and radiative feedback from, individual stars in the course of the formation of a protogalaxy. We compute in detail the H II regions of each of these sources, as well as the regions affected by their molecule-dissociating radiation. We follow the thermal, chemical, and dynamical evolution of the primordial gas, as it becomes incorporated into the protogalaxy. While the IGM is, in general, optically thick to LW photons only over physical distances of > 30 kpc at redshifts z < 20, the high molecule fraction that is built up in relic H II regions and their increasing volume-filling fraction renders even the local IGM optically thick to LW photons over physical distances of the order of a few kiloparsecs. We find that efficient accretion onto Population III relic black holes may occur after ~ 60 Myr from the time of their formation, by which time the photo-heated relic H II region gas can cool and re-collapse into the 10^6 M_solar minihalo which hosts the black hole. Also, Pop II.5 stars, postulated to have masses of the order of 10 M_solar, can likely form from this re-collapsing relic H II region gas. Overall, we find that the local radiative feedback from the first generations of stars suppresses the star formation rate by only a factor of, at most, a few.

We investigate the evolution of the primordial gas surrounding the first massive black holes formed by the collapse of Population III stars at redshifts z > 20. Carrying out three-dimensional hydrodynamical simulations using GADGET, we study the dynamical, thermal and chemical evolution of the first relic H II regions. We also carry out simulations of the mergers of relic H II regions with neighboring neutral minihaloes, which contain high density primordial gas that can accrete onto a Pop III remnant black hole. We find that there may have been a significant time delay, of order ~10^8 yr, between black hole formation and the onset of efficient accretion. The build-up of supermassive black holes, believed to power the z > 6 quasars observed in the Sloan Digital Sky Survey, therefore faces a crucial early bottleneck. More massive seed black holes may thus be required, such as those formed by the direct collapse of a primordial gas cloud facilitated by atomic line cooling. The high optical depth to Lyman-Werner (LW) photons that results from the high fraction of H_2 molecules that form in relic H II regions, combined with the continued formation of H_2 inside the dynamically expanding relic H II region, leads to shielding of the molecules inside these regions at least until a critical background LW flux of \~10^{-24} ergs s^{-1} cm^{-2} Hz^{-1} sr^{-1}, is established. Furthermore, we find that a high fraction of HD molecules, X_{HD} > 10^{-7}, is formed, potentially enabling the formation of Pop II.5 stars during later stages of structure formation when the relic H II region gas is assembled into a sufficiently deep potential well to gravitationally confine the gas again.

Gamma-Ray Bursts (GRBs) are unique probes of the cosmic star formation history and the state of the intergalactic medium up to the redshifts of the first stars. In particular, the ongoing {\it Swift} mission might be the first observatory to detect individual Population~III stars, provided that the massive, metal-free stars were able to trigger GRBs. {\it Swift} will empirically constrain the redshift at which Population~III star formation was terminated, thus providing crucial input to models of cosmic reionization and metal enrichment.

Detection of gamma-ray bursts (GRBs) from redshifts z > 7 would open a new window into the earliest epoch of cosmic star formation. We construct separate star formation histories at high redshifts for normal (Pop I and II) stars, and for predominantly massive (Pop III) stars. Based on these separate histories, we predict the GRB redshift distribution to be observed by the {\it Swift} mission. Regardless of whether Pop III progenitors are able to trigger GRBs, we find that a fraction 10% of all bursts detected by {\it Swift} will originate at z > 5. This baseline contribution is due to Pop I/II star formation which must have extended out to high redshifts in rare massive galaxies that were enriched by heavy elements earlier than the typical galaxies. In addition, we consider the possible contribution of Pop III progenitors to the observable GRB rate. Pop III stars are viable progenitors for long-duration GRBs which are triggered by the collapsar mechanism, as long as they can lose their outer envelope through mass transfer to a companion star in a close binary. We find that the likelihood of Pop III binaries to satisfy the conditions required by the collapsar mechanism could be enhanced significantly relative to Pop I/II binaries. If Pop III binaries are common, {\it Swift} will be first observatory to probe Pop III star formation at redshifts z > 7.

Simulations predict that the first stars in an LCDM universe formed at redshifts z>20 in minihalos with masses of about 10^6 M_sun. We have studied their radiative feedback by simulating the propagation of ionization fronts (I-fronts) created by these first Population III stars (M_* = 15-500 M_sun) at z=20, within the density field of a cosmological simulation of primordial star formation, outward thru the host minihalo and into the surrounding gas. A three-dimensional ray-tracing calculation tracks the I-front once the H II region evolves a "champagne flow'' inside the minihalo, after the early D-type I-front detaches from the shock and runs ahead, becoming R-type. We take account of the hydrodynamical back-reaction by an approximate model of the central wind. We find that the escape fraction of ionizing radiation from the host halo increases with stellar mass, with 0.7 < f_esc <0.9 for 80 < M_*/M_sun < 500. To quantify the ionizing efficiency of these stars as they begin cosmic reionization, we find that, for M_* > 80 M_sun, the ratio of gas mass ionized to stellar mass is ~60,000, roughly half the number of ionizing photons released per stellar baryon. Nearby minihalos are shown to trap the I-front, so their centers remain neutral. This is contrary to the recent suggestion that these stars would trigger formation of a second generation by fully ionizing neighboring minihalos, stimulating H_2 formation in their cores. Finally, we discuss how the evacuation of gas from the host halo reduces the growth and luminosity of "miniquasars'' that may form from black hole remnants of the first stars.

We find that at redshifts z > 10, HD line cooling allows strongly-shocked primordial gas to cool to the temperature of the cosmic microwave background (CMB). This temperature is the minimum value attainable via radiative cooling. Provided that the abundance of HD, normalized to the total number density, exceeds a critical level of ~ 10^{-8}, the CMB temperature floor is reached in a time which is short compared to the Hubble time. We estimate the characteristic masses of stars formed out of shocked primordial gas in the wake of the first supernovae, and resulting from the mergers of dark matter haloes during hierarchical structure formation to be ~ 10 M_{solar}. In addition, we show that cooling by HD enables the primordial gas in relic H II regions to cool to temperatures considerably lower than those reached via H_2 cooling alone. We confirm that HD cooling is unimportant in cases where the primordial gas does not go through an ionized phase, as in the formation process of the very first stars in z ~ 20 minihaloes of mass ~ 10^{6} M_{solar}.

How and when did the first generation of stars form at the end of the cosmic dark ages? Quite generically, within variants of the cold dark matter model of cosmological structure formation, the first sources of light are expected to form in ~ 10^{6} M_sun dark matter potential wells at redshifts z > 20. I discuss the physical processes that govern the formation of the first stars. These so-called Population~III stars are predicted to be predominantly very massive, and to have contributed significantly to the early reionization of the intergalactic medium. Such an early reionization epoch is inferred from the recent measurement of the Thomson optical depth by the {\it WMAP} satellite. I address the importance of heavy elements in bringing about the transition from an early star formation mode dominated by massive stars, to the familiar mode dominated by low mass stars, at later times, and present possible observational probes. This transition could have been gradual, giving rise to an intermediate-mass population of still virtually metal-free stars (``Population II.5''). These stars could have given rise to the peculiar class of black-hole forming supernovae inferred from the abundance pattern of extremely iron-poor stars.

We present a three-dimensional numerical simulation that resolves the formation process of a Population III star down to a scale of ~100 AU. The simulation is initialized on the scale of a dark matter halo of mass ~10^6 M_sun that virializes at z~20. It then follows the formation of a fully-molecular central core, and traces the accretion from the diffuse dust-free cloud onto the protostellar core for as long as ~10^4 yr, at which time the core has grown to ~50 M_sun. We find that the accretion rate starts very high, ~0.1 M_sun yr^-1, and declines rapidly thereafter approaching a power-law temporal scaling. Asymptotically, at times t > 10^3 yr after core formation, the stellar mass grows approximately as M_star=20 M_sun (t/10^3 yr)^0.4. Earlier on, accretion is faster with M_star \propto t^0.75. By extrapolating this growth over the full lifetime of very massive stars, t~3x10^6 yr, we obtain the conservative upper limit M_star < 500 M_sun. The actual stellar mass is, however, likely to be significantly smaller than this mass limit due to radiative and mechanical feedback from the protostar.

The birth and death of the first generation of stars have important implications for the thermal state and chemical properties of the intergalactic medium (IGM) in the early universe. Sometime after recombination, the neutral, chemically pristine gas was reionized by ultraviolet photons emitted from the first stars, but also enriched with heavy elements when these stars ended their lives as energetic supernovae. Using the results from previous high-resolution cosmological simulations of early structure formation that include radiative transfer, we show that a significant volume fraction of the IGM can be metal-polluted, as well as ionized, by massive Population III stars formed in small-mass (10^6-10^7 Msun) halos early on. If most of the early generation stars die as pair-instability supernovae with energies up to 10^{53} ergs, the volume-averaged mean metallicity will quickly reach Z ~ 10^{-4}Zsun by a redshift of 15-20, possibly causing a prompt transition to the formation of a stellar population that is dominated by low-mass stars. In this scenario, the early chemical enrichment history should closely trace the reionization history of the IGM, and the end of the Population III era is marked by the completion of reionization and pre-enrichment by z=15. We conclude that, while the pre-enrichment may partially account for the ``metallicity-floor'' in high-redshift Lyman-alpha clouds, it does not significantly affect the elemental abundance in the intracluster medium.

I review recent progress in understanding the formation of the first stars and quasars. The initial conditions for their emergence are given by the now firmly established model of cosmological structure formation. Numerical simulations of the collapse and fragmentation of primordial gas indicate that the first stars formed at redshifts z ~ 20 - 30, and that they were predominantly very massive, with M_* > 100 M_sun. Important uncertainties, however, remain. Paramount among them is the accretion process, which builds up the final stellar mass by incorporating part of the diffuse, dust-free envelope into the central protostellar core. The first quasars, on the other hand, are predicted to have formed later on, at z ~ 10, in more massive dark matter halos, with total masses, ~ 10^8 M_sun, characteristic of dwarf galaxies.

We review recent theoretical results on the formation of the first stars in the universe, and emphasize related open questions. In particular, we discuss the initial conditions for Population III star formation, as given by variants of the cold dark matter cosmology. Numerical simulations have investigated the collapse and the fragmentation of metal-free gas, showing that the first stars were predominantly very massive. The exact determination of the stellar masses, and the precise form of the primordial initial mass function, is still hampered by our limited understanding of the accretion physics and the protostellar feedback effects. We address the importance of heavy elements in bringing about the transition from an early star formation mode dominated by massive stars, to the familiar mode dominated by low mass stars, at later times. We show how complementary observations, both at high redshifts and in our local cosmic neighborhood, can be utilized to probe the first epoch of star formation.

The first stars in the Universe are predicted to have been much more massive than the Sun. Gravitational condensation accompanied by cooling of the primordial gas due to molecular hydrogen, yields a minimum fragmentation scale of a few hundred solar masses. Numerical simulations indicate that once a gas clump acquires this mass, it undergoes a slow, quasi-hydrostatic contraction without further fragmentation. Here we show that as soon as the primordial gas - left over from the Big Bang - is enriched by supernovae to a carbon or oxygen abundance as small as ~0.01-0.1% of that found in the Sun, cooling by singly-ionized carbon or neutral oxygen can lead to the formation of low-mass stars. This mechanism naturally accommodates the discovery of solar mass stars with unusually low (10^{-5.3} of the solar value) iron abundance but with a high (10^{-1.3} solar) carbon abundance. The minimum stellar mass at early epochs is partially regulated by the temperature of the cosmic microwave background. The derived critical abundances can be used to identify those metal-poor stars in our Milky Way galaxy with elemental patterns imprinted by the first supernovae.

We investigate the supernova explosions that end the lives of massive Population III stars in low-mass minihalos (M~10^6 M_sun) at redshifts z~20. Employing the smoothed particle hydrodynamics method, we carry out numerical simulations in a cosmological set-up of pair-instability supernovae with explosion energies of E_SN=10^51 and 10^53 ergs. We find that the more energetic explosion leads to the complete disruption of the gas in the minihalo, whereas the lower explosion energy leaves much of the halo intact. The higher energy supernova expels > 90% of the stellar metals into a region ~1 kpc across over a timescale of 3-5 Myr. Due to this burst-like initial star formation episode, a large fraction of the universe could have been endowed with a metallicity floor, Z_min>10^-4 Z_sun, already at z>15.

We consider the physical conditions under which supermassive black holes could have formed inside the first galaxies. Our SPH simulations indicate that metal-free galaxies with a virial temperature ~10^4 K and with suppressed H2 formation (due to an intergalactic UV background) tend to form a binary black hole system which contains a substantial fraction (>10%) of the total baryonic mass of the host galaxy. Fragmentation into stars is suppressed without substantial H2 cooling. Our simulations follow the condensation of ~5x10^6 M_sun around the two centers of the binary down to a scale of < 0.1pc. Low-spin galaxies form a single black hole instead. These early black holes lead to quasar activity before the epoch of reionization. Primordial black hole binaries lead to the emission of gravitational radiation at redshifts z>10 that would be detectable by LISA.

We present a model for the star formation process during the initial collapse of dark matter haloes at redshifts z=0-30. We derive a simple expression for the characteristic stellar mass scale during this initial burst of star formation. In our picture, this characteristic scale reflects both the minimum temperature to which the gas can cool (determined by the metallicity and the temperature of the cosmic microwave background) and the pressure of overlying baryons in the collapsing halo. This prescription reproduces both the large mass scales found in simulations of Population III star formation and the near solar values observed for star formation at low redshift.

The earliest stars to form in the Universe were the first sources of light, heat and metals after the Big Bang. The products of their evolution will have had a profound impact on subsequent generations of stars. Recent studies of primordial star formation have shown that, in the absence of metals (elements heavier than helium), the formation of stars with masses 100 times that of the Sun would have been strongly favoured, and that low-mass stars could not have formed before a minimum level of metal enrichment had been reached. The value of this minimum level is very uncertain, but is likely to be between 10^{-6} and 10^{-4} that of the Sun. Here we show that the recent discovery of the most iron-poor star known indicates the presence of dust in extremely low-metallicity gas, and that this dust is crucial for the formation of lower-mass second-generation stars that could survive until today. The dust provides a pathway for cooling the gas that leads to fragmentation of the precursor molecular cloud into smaller clumps, which become the lower-mass stars.

We review recent theoretical results on the formation of the first stars and quasars in the universe, and emphasize related open questions. In particular, we list important differences between the star formation process at high redshifts and in the present-day universe. We address the importance of heavy elements in bringing about the transition from an early star formation mode dominated by massive stars, to the familiar mode dominated by low mass stars, at later times. We show how gamma-ray bursts can be utilized to probe the first epoch of star formation. Finally, we discuss how the first supermassive black holes could have formed through the direct collapse of primordial gas clouds.

We investigate the impact of an early population of massive stars on their surroundings. Dissociation of molecular hydrogen by strong UV emission from such stars is expected to produce a global transition in the cooling mechanism of minihalos at a redshift of approximately 30, strongly inhibiting star formation until more massive halos can collapse. Furthermore, chemical enrichment from Pop III supernovae produces a second transition at z ~ 15-20, when the mean metallicity of the universe exceeds a critical threshold and Pop III star formation gives way to Pop II. We show that observations of high redshift supernovae with the Next Generation Space Telescope (NGST) have the potential to trace the cosmic star formation rate out to z >~ 20, provided that Pop III supernovae are at least as bright as, and ideally brighter than, type Ia supernovae. We also propose a mechanism for the formation of a novel population of extremely low metallicity stars of intermediate mass at very high redshifts, which we term Pop II.5. In our model shock compression, heating, and subsequent cooling to high density reduces the fragment mass in primordial gas to ~10 M_sun, allowing low mass stars to form. We predict the number density of relic Pop II.5 stars in the Milky Way halo today and find that, with certain assumptions, there should be ~ 10 kpc^{-3} in the solar neighborhood.

We present results from the first hydrodynamical star formation calculation to demonstrate that close binary stellar systems (separations $\lsim 10$ AU) need not be formed directly by fragmentation. Instead, a high frequency of close binaries can be produced through a combination of dynamical interactions in unstable multiple systems and the orbital decay of initially wider binaries. Orbital decay may occur due to gas accretion and/or the interaction of a binary with its circumbinary disc. These three mechanisms avoid the problems associated with the fragmentation of optically-thick gas to form close systems directly. They also result in a preference for close binaries to have roughly equal-mass components because dynamical exchange interactions and the accretion of gas with high specific angular momentum drive mass ratios towards unity. Furthermore, due to the importance of dynamical interactions, we find that stars with greater masses ought to have a higher frequency of close companions, and that many close binaries ought to have wide companions. These properties are in good agreement with the results of observational surveys.

We present results from the largest numerical simulation of star formation to resolve the fragmentation process down to the opacity limit. The simulation follows the collapse and fragmentation of a large-scale turbulent molecular cloud to form a stellar cluster and, simultaneously, the formation of circumstellar discs and binary stars. This large range of scales enables us to predict a wide variety of stellar properties for comparison with observations. The calculation clearly demonstrates that star formation is a highly-dynamic and chaotic process. Star-disc encounters form binaries and truncate discs. Stellar encounters disrupt bound multiple systems. The cloud produces roughly equal numbers of stars and brown dwarfs, with masses down to the opacity limit for fragmentation (~5 Jupiter masses). The initial mass function is consistent with a Salpeter slope (Gamma=-1.35) above 0.5 Msun, a roughly flat distribution (Gamma=0) in the range 0.006-0.5 Msun, and a sharp cutoff below ~0.005 Msun. This is consistent with recent observational surveys. The brown dwarfs form by the dynamical ejection of low-mass fragments from dynamically unstable multiple systems before the fragments have been able to accrete to stellar masses. Close binary systems (with separations <10 AU) are not formed by fragmentation in situ. Rather, they are produced by hardening of initially wider multiple systems through a combination of dynamical encounters, gas accretion, and/or the interaction with circumbinary and circumtriple discs. Finally, we find that the majority of circumstellar discs have radii less than 20 AU due to truncation in dynamical encounters. This is consistent with observations of the Orion Trapezium Cluster and implies that most stars and brown dwarfs do not form large planetary systems.

We calculate the contribution to the cosmic infrared background from very massive metal-free stars at high redshift. We explore two plausible star-formation models and two limiting cases for the reprocessing of the ionizing stellar emission. We find that Population III stars may contribute significantly to the cosmic near-infrared background if the following conditions are met: (i) The first stars were massive, with M > ~100 M_sun. (ii) Molecular hydrogen can cool baryons in low-mass haloes. (iii) Pop III star formation is ongoing, and not shut off through negative feedback effects. (iv) Virialized haloes form stars at about 40 per cent efficiency up to the redshift of reionization, z~7. (v) The escape fraction of the ionizing radiation into the intergalactic medium is small. (vi) Nearly all of the stars end up in massive black holes without contributing to the metal enrichment of the Universe.

We present results from the first hydrodynamical star formation calculation to demonstrate that brown dwarfs are a natural and frequent product of the collapse and fragmentation of a turbulent molecular cloud. The brown dwarfs form via the fragmentation of dense molecular gas in unstable multiple systems and are ejected from the dense gas before they have been able to accrete to stellar masses. Thus, they can be viewed as `failed stars'. Approximately three quarters of the brown dwarfs form in gravitationally-unstable circumstellar discs while the remainder form in collapsing filaments of molecular gas. These formation mechanisms are very efficient, producing roughly the same number of brown dwarfs as stars, in agreement with recent observations. However, because close dynamical interactions are involved in their formation, we find a very low frequency of binary brown dwarf systems ($\lsim 5$%) and that those binary brown dwarf systems that do exist must be close $\lsim 10$ AU. Similarly, we find that young brown dwarfs with large circumstellar discs (radii $\gsim 10$ AU) are rare ($\approx 5$%).

We predict the redshift distribution of Gamma-Ray Bursts (GRBs) assuming that they trace the cosmic star formation history. We find that a fraction >~50% of all GRBs on the sky originate at a redshift z >~5, even though the fraction of the total stellar mass formed by z~5 is only ~15%. These two fractions are significantly different because they involve different cosmological factors when integrating the star formation rate over redshift. Hence, deep observations of transient events, such as GRB afterglows or supernovae, provide an ideal strategy for probing the high-redshift universe. We also show that the intrinsic distribution of GRB durations is bimodal but significantly narrower and shifted towards shorter durations than the observed distribution.

We explore a mechanism for the formation of the first globular clusters, operating during the assembly of dwarf galaxies at high redshifts, z > 10. The substructure in the dark matter and the corresponding potential wells are responsible for setting the cluster scale of ~10^5 M_sun. The second mass scale in the formation problem, the stellar scale of ~1 M_sun, is determined in turn by the processes that cool the gas. We address the origin of the first, cluster scale by means of three-dimensional numerical simulations of the collapsing dark matter and gaseous components. We find that the gas falls into the deepest dark subhalos, resulting in a system of ~ 5 proto-globular clouds. The incipient globular clusters lose their individual dark halos in the process of violent relaxation, leading to the build-up of the general dark halo around the dwarf galaxy.

Paramount among the processes that ended the cosmic `dark ages' must have been the formation of the first generation of stars (the so-called Population III). We summarize recent progress in constraining its nature, and we discuss the basic physical reasons for why the first stars might have been quite massive, possibly even very massive with M > 100 M_sun.

Recent theoretical investigations have suggested that the formation of the very first stars, forming out of metal-free gas, was fundamentally different from the present-day case. In this paper, we study the effect of metallicity on the evolution of the gas in a collapsing dark matter mini-halo. We model such a system as an isolated 3\sigma peak of mass 2x10^6 M_sun that collapses at z_coll=30, using smoothed particle hydrodynamics. The gas has a supposed level of pre-enrichment of either 10^-4 Z_sun or 10^-3 Z_sun. We find that the evolution proceeds very differently for the two cases. The gas in the lower metallicity simulation fails to undergo continued collapse and fragmentation, whereas the gas in the higher metallicity case dissipatively settles into the center of the dark matter halo. The central gas, characterized by densities n > 10^4 cm^-3, and a temperature, T \sim 90 K, which closely follows that of the CMB, is gravitationally unstable and undergoes vigorous fragmentation. We discuss the physical reason for the existence of a critical metallicity, Z_crit \sim 5x10^-4 Z_sun, and its possible dependence on redshift. Compared to the pure H/He case, the fragmentation of the 10^-3 Z_sun gas leads to a larger relative number of low-mass clumps.

To constrain the nature of the very first stars, we investigate the collapse and fragmentation of primordial, metal-free gas clouds. We explore the physics of primordial star formation by means of three-dimensional simulations of the dark matter and gas components, using smoothed particle hydrodynamics, under a wide range of initial conditions, including the initial spin, the total mass of the halo, the redshift of virialization, the power spectrum of the DM fluctuations, the presence of HD cooling, and the number of particles employed in the simulation. We find characteristic values for the temperature, T ~ a few 100 K, and the density, n ~ 10^3-10^4 cm^-3, characterising the gas at the end of the initial free-fall phase. These values are rather insensitive to the initial conditions. The corresponding Jeans mass is M_J ~ 10^3 M_sun. The existence of these characteristic values has a robust explanation in the microphysics of H2 cooling, connected to the minimum temperature that can be reached with the H2 coolant, and to the critical density at which the transition takes place betweeb levels being populated according to NLTE, and according to LTE. In all cases, the gas dissipatively settles into an irregular, central configuration which has a filamentary and knotty appearance. The fluid regions with the highest densities are the first to undergo runaway collapse due to gravitational instability, and to form clumps with initial masses ~ 10^3 M_sun, close to the characteristic Jeans scale. These results suggest that the first stars might have been quite massive, possibly even very massive with M_star > 100 M_sun.

We calculate the generic spectral signature of an early population of massive stars at high redshifts. For metal-free stars in the mass range of 100-1000 solar masses, we find that the combined spectral luminosity per unit stellar mass is almost independent of the mass distribution of these stars. To zeroth order, the generic spectrum resembles a black-body with an effective temperature of ~10^5 K, making these stars highly efficient at ionizing hydrogen and helium. The production rate of ionizing radiation per stellar mass by stars more massive than 100 solar masses is larger by about an order of magnitude for hydrogen and He I, and by about two orders of magnitude for He II, than the emission from a standard initial mass function. This would result in unusually strong hydrogen and helium recombination lines from the surrounding interstellar medium. It could also alleviate the current difficulty of ionizing the intergalactic medium at z>6 with the cosmic star formation rate inferred at somewhat lower redshifts.

We briefly review the motivations for studying the formation of the first ``Population III'' stars and present recent results from our numerical simulations in this area. We discuss the new questions raised as a result of the simulations presented by us and others at this meeting.

In order to constrain the initial mass function (IMF) of the first generation of stars (Population III), we investigate the fragmentation properties of metal-free gas in the context of a hierarchical model of structure formation. We investigate the evolution of an isolated 3-sigma peak of mass 2x10^6 M_solar which collapses at z_coll=30 using Smoothed Particle Hydrodynamics. We find that the gas dissipatively settles into a rotationally supported disk which has a very filamentary morphology. The gas in these filaments is Jeans unstable with M_J~10^3 M_solar. Fragmentation leads to the formation of high density (n>10^8 cm^-3) clumps which subsequently grow in mass by accreting surrounding gas and by merging with other clumps up to masses of ~10^4 M_solar. This suggests that the very first stars were rather massive. We explore the complex dynamics of the merging and tidal disruption of these clumps by following their evolution over a few dynamical times.

We have constructed a broad-band spectrum for GRB 930131 (the ``Superbowl Burst''), ranging from 20 keV to 200 MeV, by combining spectral information from the Gamma Ray Observatory's BATSE, COMPTEL and EGRET instruments. We present general methods for combining spectra from different time intervals obtained by the same instrument as well as for combining spectra from the same time interval taken by different instruments. The resulting spectrum is remarkably flat (in nu F_nu-space) up to high energies. We find that the spectral shape can be successfully fitted by the shocked synchrotron emission model of Tavani. We present evidence that the flatness of the spectrum at high energies is not due to spectral time-variability.

In order to constrain the initial mass function of the first generation of stars (Population III), one has to study the fragmentation properties of primordial gas. We present results from 3D simulations, based on Smoothed Particle Hydrodynamics, which explore the idea that small-scale fluctuations in the (cold) dark matter recreate a filamentary and clumpy structure in the gas component on scales smaller than the initial Jeans mass, where all primordial fluctuations would have been wiped out.