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

Models of the decoupling of baryons and photons during the recombination epoch predict the existence of a large-scale velocity offset between baryons and dark matter at later times, the so-called streaming velocity. In this paper, we use high resolution numerical simulations to investigate the impact of this streaming velocity on the spin and shape distributions of high-redshift minihalos, the formation sites of the earliest generation of stars. We find that the presence of a streaming velocity has a negligble effect on the spin and shape of the dark matter component of the minihalos. However, it strongly affects the behaviour of the gas component. The most probable spin parameter increase from $\sim$0.03 in the absence of streaming to $\sim$0.15 for a run with a streaming velocity of three times $\sigma_{\rm rms}$. The gas within the minihalos becomes increasingly less spherical and more oblate as the streaming velocity increases, with dense clumps being found at larger distances from the halo centre. The impact of the streaming velocity is also mass-dependent: less massive objects are influenced more strongly, on account of their shallower potential wells. The number of halos in which gas cooling and runaway gravitational collapse occurs decreases substantially as the streaming velocity increases. However, the spin and shape distributions of gas that does manage to cool and collapse are insensitive to the value of the streaming velocity and uncorrelated with the same properties measured on the scale of the halo as a whole.

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.

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.

This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources.

How, when and where the first stars formed are fundamental questions regarding the epoch of Cosmic Dawn. A second order effect in the fluid equations was recently found to make a significant contribution: an offset velocity between gas and dark matter, the so-called streaming velocity. Previous simulations of a limited number of low-mass dark matter haloes suggest that this streaming velocity can delay the formation of the first stars and decrease halo gas fractions and the halo mass function in the low mass regime. However, a systematic exploration of its effects in a large sample of haloes has been lacking until now. In this paper, we present results from a set of cosmological simulations of regions of the Universe with different streaming velocities performed with the moving mesh code AREPO. Our simulations have very high mass resolution, enabling us to accurately resolve minihaloes as small as $10^5 \: {\rm M_{\odot}}$. We show that in the absence of streaming, the least massive halo that contains cold gas has a mass $M_{\rm halo, min} = 5 \times 10^{5} \: {\rm M_{\odot}}$, but that cooling only becomes efficient in a majority of haloes for halo masses greater than $M_{\rm halo,50\%} = 1.6 \times 10^6 \: {\rm M_{\odot}}$. In regions with non-zero streaming velocities, $M_{\rm halo, min}$ and $M_{\rm halo,50\%}$ both increase significantly, by around a factor of a few for each one sigma increase in the value of the local streaming velocity. As a result, in regions with streaming velocities $v_\mathrm{stream} \ge 3\,\sigma_\mathrm{rms}$, cooling of gas in minihaloes is completely suppressed, implying that the first stars in these regions form within atomic cooling haloes.

The spin and shape of galaxies at the present day have been well-studied both observationally and theoretically. At high redshifts, however, we have to rely on numerical simulations. In this study, we investigate the shape and spin of minihaloes with masses of $M \sim 10^5$--$10^7 \, {\rm M_{\odot}}$ which are of particular interest as they are the sites where the first stars in the Universe form. We analyse a large sample of these minihaloes, selected from a high resolution cosmological simulation. The first minihaloes form at $z \simeq 24$ and by the end of the simulation at $z \simeq 14$ our sample includes $\sim 9000$ minihaloes. We find that the spin parameter of the minihaloes follows a log-normal distribution with minimal dependence on redshift. Most minihaloes are prolate, but those formed at the highest redshifts are more prolate than those formed at lower redshifts. On the scale of the virial radius, there is a good correlation between the shape and spin of the gas and that of the dark matter. However, this correlation breaks down in gas which is cooling and undergoing gravitational collapse. We show, contrary to previous assumptions, that although the direction of the spin of the central dense gas correlates well with that of the halo, the magnitude of the spin of the dense gas is uncorrelated with that of the halo. Therefore, measurements of the spin of minihaloes on large scales tell us little about the angular momentum of the gas responsible for forming the first stars.

We study the influence of a high baryonic streaming velocity on the formation of direct collapse black holes (DCBHs) with the help of cosmological simulations carried out using the moving mesh code {\sc arepo}. We show that a streaming velocity that is as large as three times the root-mean-squared value is effective at suppressing the formation of H$_{2}$-cooled minihaloes, while still allowing larger atomic cooling haloes (ACHs) to form. We find that enough H$_{2}$ forms in the centre of these ACHs to effectively cool the gas, demonstrating that a high streaming velocity by itself cannot produce the conditions required for DCBH formation. However, we argue that high streaming velocity regions do provide an ideal environment for the formation of DCBHs in close pairs of ACHs (the "synchronised halo" model). Due to the absence of star formation in minihaloes, the gas remains chemically pristine until the ACHs form. If two such haloes form with only a small separation in time and space, then the one forming stars earlier can provide enough ultraviolet radiation to suppress H$_{2}$ cooling in the other, allowing it to collapse to form a DCBH. Baryonic streaming may therefore play a crucial role in the formation of the seeds of the highest redshift quasars.

Direct collapse black holes forming in pristine, atomically-cooling haloes at $z \approx 10-20$ may act as the seeds of supermassive black holes (BH) at high redshifts. In order to create a massive BH seed, the host halo needs to be prevented from forming stars. H$_2$ therefore needs to be irradiated by a large flux of Lyman-Werner (LW) UV photons in order to suppress H$_2$ cooling. A key uncertainty in this scenario is the escape fraction of LW radiation from first galaxies, the dominant source of UV photons at this epoch. To better constrain this escape fraction, we have performed radiation-hydrodynamical simulations of the growth of HII regions and their associated photodissociation regions in the first galaxies using the ZEUS-MP code. We find that the LW escape fraction crucially depends on the propagation of the ionisation front (I-front). For an R-type I-front overrunning the halo, the LW escape fraction is always larger than 95%. If the halo recombines later from the outside--in, due to a softened and weakened spectrum, the LW escape fraction in the rest-frame of the halo (the near-field) drops to zero. A detailed and careful analysis is required to analyse slowly moving, D-type I-fronts, where the escape fraction depends on the microphysics and can be as small as 3% in the near-field and 61% in the far-field or as large as 100% in both the near-field and the far-field.

Population III stars can regulate star formation in the primordial Universe in several ways. They can ionize nearby halos, and even if their ionizing photons are trapped by their own halos, their Lyman-Werner (LW) photons can still escape and destroy H$_2$ in other halos, preventing them from cooling and forming stars. LW escape fractions are thus a key parameter in cosmological simulations of early reionization and star formation but have not yet been parametrized for realistic halos by halo or stellar mass. To do so, we perform radiation hydrodynamical simulations of LW UV escape from 9--120 M$_{\odot}$ Pop III stars in $10^5$ to $10^7$ M$_{\odot}$ halos with ZEUS-MP. We find that photons in the LW lines (i.e. those responsible for destroying H$_{2}$ in nearby systems) have escape fractions ranging from 0% to 85%. No LW photons escape the most massive halo in our sample, even from the most massive star. Escape fractions for photons elsewhere in the 11.18--13.6~eV energy range, which can be redshifted into the LW lines at cosmological distances, are generally much higher, being above 60% for all but the least massive stars in the most massive halos. We find that shielding of H$_2$ by neutral hydrogen, which has been neglected in most studies to date, produces escape fractions that are up to a factor of three smaller than those predicted by H$_2$ self-shielding alone.

The stellar velocity dispersion as a function of the galactocentric radius of an early-type galaxy can generally be well approximated by a power law $\sigma \propto r^{\beta}$. However, some observed dispersion profiles show a deviation from this fit at intermediate radii, usually between one and three $R_{\mathrm{eff}}$, where the velocity dispersion remains constant with radius, showing a bump-like behavior, which we term the "$\sigma$-"bump. To understand the origin of this $\sigma$-bump, we study a set of simulated early-type galaxies formed in major mergers. We find the $\sigma$-bump in all of our simulated early-type galaxies, with the size and position of the bump slightly varying from galaxy to galaxy, leading to the assumption that the bump is a characteristic of the major merger formation scenario. The feature can be seen both in the intrinsic and projected stellar velocity dispersions. In contrast to shells that form during the merger event but evolve with time and finally disappear, the $\sigma$-bump stays nearly constant with radius and is a permanent feature that is preserved until the end of the simulation. The $\sigma$-bump is not seen in the dark matter and gas components and we therefore conclude that it is a purely stellar feature of merger remnants.