Next generation telescopes such as the James Webb Space Telescope (JWST) and the Wide Field Infrared Survey Telescope (WFIRST) will enable us to study the first billion years of our Universe in unprecedented detail. In this work we use the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS) framework, that couples galaxy formation and reionization (for a wide range of reionization feedback models), to estimate the cosmic variance expected in the UV Luminosity Function (UV LF) and the Stellar Mass Function (SMF) in JWST surveys. We find that different reionization scenarios play a minor role in the cosmic variance. Most of the cosmic variance is completely driven by the underlying density field and increases above 100% for M$_{UV} \sim$ -17.5 (-20) at z = 12 (6) for the JADES-deep survey (the deep JWST Advanced Extragalactic Survey with an area of 46 arcmin$^2$); the cosmic variance decreases with an increasing survey area roughly independently of redshift. We find that the faint-end slope of the Lyman Break Galaxies (LBGs) UV LF becomes increasingly shallower with increasing reionization feedback and show how JWST observations will be able to distinguish between different models of reionization feedback at z > 9, even accounting for cosmic variance. We also show the environments (in terms of density and ionization fields) of Lyman Break Galaxies during the EoR. Finally, we also provide a public software tool to allow interested readers to compute cosmic variance for different redshifts and survey areas.

We introduce a new self-consistent model of galaxy evolution and reionization, ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS), which couples a state-of-the-art N-body simulation with the semi-analytical galaxy evolution DELPHI and the semi-numerical reionization scheme CIFOG. ASTRAEUS includes all the key processes of galaxy formation and evolution (including accretion, mergers, supernova and radiative feedback) and follows the time and spatial evolution of the ionized regions in the intergalactic medium (IGM). Importantly, it explores different radiative feedback models that cover the physically plausible parameter space, ranging from a weak and delayed to a strong and immediate reduction of gas mass available for star formation. From our simulation suite that covers the different radiative feedback prescriptions and ionization topologies, we find that radiative feedback continuously reduces star formation in galaxies with $M_h<10^{9.5}M_{\odot}$ upon local reionization; larger mass halos are unaffected even for the strongest and immediate radiative feedback cases during reionization. For this reason, the ionization topologies of different radiative feedback scenarios differ only on scales smaller than $1-2$Mpc, and significant deviations are only found when physical parameters (e.g. the escape fraction of ionizing photons) are altered based on galactic properties. Finally, we find observables (the ultra-violet luminosity function, stellar mass function, reionization histories and ionization topologies) are hardly affected by the choice of the used stellar population synthesis models that either model single stars or binaries.

We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties and cover the physically plausible range, i.e. $f_\mathrm{esc}$ effectively decreasing, being constant, or increasing with halo mass. We find the 21cm bispectrum to be sensitive to the resulting ionization topologies that significantly differ in their size distribution of ionized and neutral regions throughout reionization. From squeezed to stretched triangles, the 21cm bispectra features a change of sign from negative to positive values, with ionized and neutral regions representing below-average and above-average concentrations contributing negatively and positively, respectively. The position of the change of sign provides a tracer of the size distribution of the ionized and neutral regions, and allows us to identify three major regimes that the 21cm bispectrum undergoes during reionization. In particular the regime during the early stages of reionization, where the 21cm bispectrum tracks the peak of the size distribution of the ionized regions, provides exciting prospects for pinning down reionization with the forthcoming Square Kilometre Array.

One of the most important parameters in characterizing the Epoch of Reionization, the escape fraction of ionizing photons, $f_\mathrm{esc}$, remains unconstrained both observationally and theoretically. With recent work highlighting the impact of galaxy-scale feedback on the instantaneous value of $f_\mathrm{esc}$, it is important to develop a model in which reionization is self-consistently coupled to galaxy evolution. In this work, we present such a model and explore how physically motivated functional forms of $f_\mathrm{esc}$ affect the evolution of ionized hydrogen within the intergalactic medium. Using the $21$cm power spectrum evolution, we investigate the likelihood of observationally distinguishing between a constant $f_\mathrm{esc}$ and other models that depend upon different forms of galaxy feedback. We find that changing the underlying connection between $f_\mathrm{esc}$ and galaxy feedback drastically alters the large-scale $21$cm power. The upcoming Square Kilometre Array Low Frequency instrument possesses the sensitivity to differentiate between our models at a fixed optical depth, requiring only $200$ hours of integration time focused on redshifts $z = 7.5-8.5$. Generalizing these results to account for a varying optical depth will require multiple $800$ hour observations spanning redshifts $z = 7-10$. This presents an exciting opportunity to observationally constrain one of the most elusive parameters during the Epoch of Reionization.

Weakly interacting cold dark matter (CDM) particles, which are otherwise extremely successful in explaining various cosmological observations, exhibit a number of problems on small scales. One possible way of solving these problems is to invoke (so-called) warm dark matter (WDM) particles with masses $m_x \sim$ keV. Since the formation of structure is delayed in such WDM models, it is natural to expect that they can be constrained using observations related to the first stars, e.g., the 21 cm signal from cosmic dawn. In this work, we use a detailed galaxy formation model, Delphi, to calculate the 21 cm signal at high-redshifts and compare this to the recent EDGES observations. We find that while CDM and 5 keV WDM models can obtain a 21 cm signal within the observed redshift range, reproducing the amplitude of the observations requires the introduction of an excess radio background. On the other hand, WDM models with $m_x \leq 3$ keV can be ruled out since they are unable to match either the redshift range or the amplitude of the EDGES signal, irrespective of the parameters used. Comparable to values obtained from the low-redshift Lyman Alpha forest, our results extend constraints on the WDM particle to an era inaccessible by any other means; additional forthcoming 21 cm data from the era of cosmic dawn will be crucial in refining such constraints.

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.

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.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

We present deep spectroscopic follow-up observations of the Bremer Deep Field (BDF) where the two $z\sim$7 bright Ly$\alpha$ emitters (LAE) BDF521 and BDF3299 were previously discovered by Vanzella et al. (2011) and where a factor of $\sim$3-4 overdensity of faint LBGs has been found by Castellano et al. (2016). We confirm a new bright Ly$\alpha$ emitter, BDF2195, at the same redshift of BDF521, $z=7.008$, and at only $\sim$90 kpc physical distance from it, confirming that the BDF area is likely an overdense, reionized region. A quantitative assessment of the Ly$\alpha$ fraction shows that the number of detected bright emitters is much higher than the average found at z$\sim$7, suggesting a high Ly$\alpha$ transmission through the inter-galactic medium (IGM). However, the line visibility from fainter galaxies is at odds with this finding, since no Ly$\alpha$ emission is found in any of the observed candidates with $M_{UV}>$-20.25. This discrepancy can be understood either if some mechanism prevents Ly$\alpha$ emission from fainter galaxies within the ionized bubbles from reaching the observer, or if faint galaxies are located outside the reionized area and bright LAEs are solely responsible for the creation of their own HII regions. A thorough assessment of the nature of the BDF region and of its sources of re-ionizing radiation will be made possible by JWST spectroscopic capabilities.

The exact role of quasars during the Epoch of Reionization remains uncertain. With consensus leaning towards quasars producing a negligible amount of ionizing photons, we pose an alternate question: Can quasars indirectly contribute to reionization by allowing ionizing photons from stars to escape more easily? Using the Semi-Analytic Galaxy Evolution model to evolve a galaxy population through cosmic time, we construct an idealized scenario in which the escape fraction of stellar ionizing photons ($f_\mathrm{esc}$) is boosted following quasar wind events, potentially for several dynamical times. We find that under this scenario, the mean value of $f_\mathrm{esc}$ as a function of galaxy stellar mass peaks for intermediate mass galaxies. This mass dependence will have consequences for the 21cm power spectrum, enhancing power at small scales and suppressing it at large scales. This hints that whilst quasars may not directly contribute to the ionizing photon budget, they could influence reionization indirectly by altering the topology of ionized regions.

Detections of the cross correlation signal between the 21cm signal during reionization and high-redshift Lyman Alpha emitters (LAEs) are subject to observational uncertainties which mainly include systematics associated with radio interferometers and LAE selection. These uncertainties can be reduced by increasing the survey volume and/or the survey luminosity limit, i.e. the faintest detectable Lyman Alpha (Ly$\alpha$) luminosity. We use our model of high-redshift LAEs and the underlying reionization state to compute the uncertainties of the 21cm-LAE cross correlation function at $z\simeq6.6$ for observations with SKA1-Low and LAE surveys with $\Delta z=0.1$ for three different values of the average IGM ionization state ($\langle\chi_\mathrm{HI}\rangle$=0.1, 0.25, 0.5). At $z\simeq6.6$, we find SILVERRUSH type surveys, with a field of view of 21 deg$^2$ and survey luminosity limits of $L_\alpha\geq7.9\times10^{42}$erg~s$^{-1}$, to be optimal to distinguish between an inter-galactic medium (IGM) that is 50%, 25% and 10% neutral, while surveys with smaller fields of view and lower survey luminosity limits, such as the 5 and 10 deg$^2$ surveys with WFIRST, can only discriminate between a 50% and 10% neutral IGM.

We have developed a modular semi-numerical code that computes the time and spatially dependent ionization of neutral hydrogen (HI), neutral (HeI) and singly ionized helium (HeII) in the intergalactic medium (IGM). The model accounts for recombinations and provides different descriptions for the photoionization rate that are used to calculate the residual HI fraction in ionized regions. We compare different semi-numerical reionization schemes to a radiative transfer (RT) simulation. We use the RT simulation as a benchmark, and find that the semi-numerical approaches produce similar HII and HeII morphologies and power spectra of the HI 21cm signal throughout reionization. As we do not track partial ionization of HeII, the extent of the double ionized helium (HeIII) regions is consistently smaller. In contrast to previous comparison projects, the ionizing emissivity in our semi-numerical scheme is not adjusted to reproduce the redshift evolution of the RT simulation, but directly derived from the RT simulation spectra. Among schemes that identify the ionized regions by the ratio of the number of ionization and absorption events on different spatial smoothing scales, we find those that mark the entire sphere as ionized when the ionization criterion is fulfilled to result in significantly accelerated reionization compared to the RT simulation. Conversely, those that flag only the central cell as ionized yield very similar but slightly delayed redshift evolution of reionization, with up to 20% ionizing photons lost. Despite the overall agreement with the RT simulation, our results suggests that constraining ionizing emissivity sensitive parameters from semi-numerical galaxy formation-reionization models are subject to photon nonconservation.

We study the signatures of reionization and ionizing properties of the early galaxies in the cross-correlations between the 21cm emission from the spin-flip transition of neutral hydrogen (HI) and the underlying galaxy population, in particular a sub-population of galaxies visible as Lyman Alpha Emitters (LAEs). With both observables simultaneously derived from a $z\simeq6.6$ hydrodynamical simulation (GADGET-2) snapshot post-processed with a radiative transfer code (pCRASH) and a dust model, we perform a parameter study and aim to constrain both the average intergalactic medium (IGM) ionization state ($1-\langle \chi_{HI} \rangle$) and the reionization topology (outside-in versus inside-out). We find that in our model LAEs occupy the densest and most-ionized regions resulting in a very strong anti-correlation between the LAEs and the 21cm emission. A 1000h SKA-LOW1 - Subaru Hyper Suprime Cam experiment can provide constraints on $\langle \chi_{HI} \rangle$, allowing us to distinguish between IGM ionization levels of 50%, 25%, 10% and fully ionized at scales $r<10$ comoving Mpc (assuming foreground avoidance for SKA). Our results support the inside-out reionization scenario where the densest knots (under-dense voids) are ionized first (last) for $\langle \chi_{HI} \rangle >= 0.1$. Further, 1000h SKA-LOW1 observations should be able to confirm the inside-out scenario by detecting a lower 21cm brightness temperature (by about 2-10 mK) in the densest regions ($> 2$ arcminute scales) hosting LAEs compared to lower-density regions devoid of them.

We present the analysis of deep HST multi-band imaging of the BDF field specifically designed to identify faint companions around two of the few Ly-alpha emitting galaxies spectroscopically confirmed at z~7 (Vanzella et al. 2011). Although separated by only 4.4 proper Mpc these galaxies cannot generate HII regions large enough to explain visibility of their Ly-alpha line, thus requiring a population of fainter ionizing sources in their vicinity. We use deep HST and VLT-Hawk-I data to select z~7 Lyman break galaxies around the emitters. We select 6 new robust z~7 LBGs at Y~26.5-27.5 whose average spectral energy distribution is consistent with the objects being at the redshift of the close-by Ly-alpha emitters. The resulting number density of z~7 LBGs in the BDF field is a factor ~3-4 higher than expected in random pointings of the same size. We compare these findings with cosmological hydrodynamic plus radiative transfer simulations of a universe with a half neutral IGM: we find that indeed Ly-alpha emitter pairs are only found in completely ionized regions characterized by significant LBG overdensities. Our findings match the theoretical prediction that the first ionization fronts are generated within significant galaxy overdensities and support a scenario where faint, "normal" star-forming galaxies are responsible for reionization.

We calculate Lyman Alpha Emitter (LAE) angular correlation functions (ACFs) at $z \simeq 6.6$ and the fraction of lifetime (for the 100 Myrs preceding $z\simeq6.6$) galaxies spend as Lyman Break Galaxies (LBGs) or as LBGs with Lyman Alpha (Ly$\alpha$) emission using a model that combines SPH cosmological simulations (GADGET-2), dust attenuation and a radiative transfer code (pCRASH). The ACFs are a powerful tool that significantly narrows the 3D parameter space allowed by LAE Ly$\alpha$ and UV luminosity functions (LFs) alone. With this work, we simultaneously constrain the escape fraction of ionizing photons $f_{esc}=0.05-0.5$, the mean fraction of neutral hydrogen in the intergalactic medium (IGM) $\langle \chi_{HI} \rangle \leq 0.01$ and the dust-dependent ratio of the escape fractions of Ly$\alpha$ and UV continuum photons $f_{\alpha}/f_c=0.6-1.2$. Our results show that reionization has the largest impact on the amplitude of the ACFs, and its imprints are clearly distinguishable from those of $f_{esc}$ and $f_\alpha/f_c$. We also show that galaxies with a critical stellar mass of $M_* = 10^{8.5} (10^{9.5})M_{\odot}$ produce enough luminosity to stay visible as LBGs (LAEs). Finally, the fraction of time during the past 100 Myrs prior to z=6.6 a galaxy spends as a LBG or as a LBG with Ly$\alpha$ emission increases with the UV magnitude (and the stellar mass $M_*$): considering observed (dust and IGM attenuated) luminosities, the fraction of time a galaxy spends as a LBG (LAE) increases from 65% to 100% (0-100%) as $M_{UV}$ decreases from $M_{UV} = -18.0$ to $-23.5$ ($M_*$ increases from $10^8-10^{10.5} M_{\odot}$). Thus in our model the brightest (most massive) LBGs most often show Ly$\alpha$ emission.

We build a physical model for high-redshift Lyman Alpha emitters (LAEs) by coupling state of the art cosmological simulations (GADGET-2) with a dust model and a radiative transfer code (pCRASH). We post-process the cosmological simulation with pCRASH using five different values of the escape fraction of hydrogen ionizing photons (f_esc=0.05,0.25,0.5,0.75,0.95) until reionization is complete, i.e. the average neutral hydrogen fraction drops to <X_HI>~10^-4. Then, the only free-parameter left to match model results to the observed Lya and UV luminosity functions of LAEs at z~6.6 is the relative escape of Lyman Alpha (Lya) and continuum photons from the galactic environment (f_alpha/f_c). We find a three-dimensional degeneracy such that the theoretical model can be reconciled with observations for an IGM Lya transmission <T_alpha>_LAE~38-50% (which translates to <X_HI>~0.5-10^-4 for Gaussian emission lines), f_esc~0.05-0.50 and f_alpha/f_c~0.6-1.8.

We consider a relativistically moving blob consisting of an isotropic electron distribution that Compton-scatters photons from an external isotropic radiation field. We compute the resulting beaming pattern, i.e. the distribution of the scattered photons, in the blob frame as well as in the observer's frame by using the full Klein-Nishina cross section and the exact incident photon distribution. In the Thomson regime the comparison of our approach with Dermer 1995 results in concurrent characteristics but different absolute number of the scattered photons by a factor of f_corr = 3.09. Additionally, our calculation yields a slightly lower boost factor which varies the more from the corresponding value in Dermer 1995 the higher the spectral index p of the electron distribution gets.