Many works have concentrated on the observable signatures of the dark matter being an ultralight axion-like particle (ALP). We concentrate on a particularly dramatic signature in the late-time cosmological matter power spectrum that occurs if the symmetry breaking that establishes the ALP happens after inflation -- white-noise density fluctuations that dominate at small scales over the adiabatic fluctuations from inflation. These fluctuations alter the early history of nonlinear structure formation. We find that for symmetry breaking scales of $f_A \sim 10^{13}-10^{15}$GeV, which requires a high effective maximum temperature after inflation, ALP dark matter with particle mass of $m_A \sim 10^{-13}-10^{-20}$eV could significantly change the number of high-redshift dwarf galaxies, the reionization history, and the Ly$\alpha$ forest. We consider all three observables. We find that the Ly$\alpha$ forest is the most constraining of current observables, excluding $f_A \gtrsim 10^{15}$GeV ($m_A \lesssim 10^{-17}$eV) in the simplest model for the ALP and considerably lower values in models coupled to a hidden strongly interacting sector ($f_A \gtrsim 10^{13}$GeV and $m_A \lesssim 10^{-13}$eV). Observations that constrain the extremely high-redshift tail of reionization may disfavor similar levels of isocurvature fluctuations as the forest. Future $z\sim 20-30$ 21cm observations have the potential to improve these constraints further using that the supersonic motions of the isocurvature-enhanced abundance of $\sim10^4M_\odot$ halos would shock heat the baryons, sourcing large BAO features.

More than three quarters of the baryonic content of the Universe resides in a highly diffuse state that is difficult to observe, with only a small fraction directly observed in galaxies and galaxy clusters. Censuses of the nearby Universe have used absorption line spectroscopy to observe these invisible baryons, but these measurements rely on large and uncertain corrections and are insensitive to the majority of the volume, and likely mass. Specifically, quasar spectroscopy is sensitive either to only the very trace amounts of Hydrogen that exists in the atomic state, or highly ionized and enriched gas in denser regions near galaxies. Sunyaev-Zel'dovich analyses provide evidence of some of the gas in filamentary structures and studies of X-ray emission are most sensitive to gas near galaxy clusters. Here we report the direct measurement of the baryon content of the Universe using the dispersion of a sample of localized fast radio bursts (FRBs), thus utilizing an effect that measures the electron column density along each sight line and accounts for every ionised baryon. We augment the sample of published arcsecond-localized FRBs with a further four new localizations to host galaxies which have measured redshifts of 0.291, 0.118, 0.378 and 0.522, completing a sample sufficiently large to account for dispersion variations along the line of sight and in the host galaxy environment to derive a cosmic baryon density of $\Omega_{b} = 0.051_{-0.025}^{+0.021} \, h_{70}^{-1}$ (95% confidence). This independent measurement is consistent with Cosmic Microwave Background and Big Bang Nucleosynthesis values.

We examine the impact of the baryon-dark matter relative velocity on intergalactic small-scale structure and the 21 cm signal during the epoch of reionization (EoR). Streaming velocities reduced the clumping of the intergalactic medium (IGM) on mass scales of $\sim 10^4 - 10^8$ M$_{\odot}$. This effect produced a distinct baryon acoustic oscillation (BAO) feature in the EoR 21 cm power spectrum at wave numbers $k\sim 0.1$ h/Mpc, near which forthcoming surveys will be most sensitive. In contrast to the highly uncertain impact of streaming velocities on star formation, the effect on clumping can be calculated more precisely because it is set mainly by cosmology and straightforward gas dynamics. We quantify the latter using fully coupled radiation-hydrodynamics simulations that capture the Jeans scale of the pre-EoR gas. The clumping factor of ionized gas is reduced by 5-10\% in regions that had typical streaming velocities of 30 km/s at recombination. The suppression peaks $\approx 5$ Myr after a region is reionized, but washes out within 200 Myr due to pressure smoothing of the gas. Using these results, we model the corresponding impact on the EoR 21 cm power spectrum and find that the BAO feature is most likely to appear early in reionization ($\approx$ 10 \% ionization).During this phase, the signal may appear at the 1 \% (5 \%) level at $k \sim 0.1 (0.06)$ h/\Mpc with an amplitude that varies by less than a factor of 10 across a range of reionization histories. We also provide a simple model for the signal originating from streaming velocity's impact on ionizing sources, which can vary by 4 orders of magnitude depending on highly uncertain source properties. We find that the clumping-driven signal is likely to dominate the source-driven one unless population III star formation in halos of masses $10^6 - 10^8$ M$_{\odot}$ was efficient enough to drive the first 10\% of reionization.

The intergalactic medium is expected to clump on scales down to $10^4-10^8$ M$_{\odot}$ before the onset of reionization. The impact of these small-scale structures on reionization is poorly understood despite the modern understanding that gas clumpiness limits the growth of H II regions. We use a suite of radiation-hydrodynamics simulations that capture the $\sim 10^4$ M$_{\odot}$ Jeans mass of unheated gas to study density fluctuations during reionization. Our simulations track the complex ionization and hydrodynamical response of gas in the wake of ionization fronts. The clumping factor of ionized gas (proportional to the recombination rate) rises to a peak value of $5-20$ approximately $\Delta t = 10$ Myr after ionization front passage, depending on the incident intensity, redshift, and degree to which the gas had been pre-heated by the first X-ray sources. The clumping factor reaches its relaxed value of $\approx 3$ by $\Delta t = 300$ Myr. The mean free path of Lyman-limit photons evolves in unison, being up to several times shorter in un-relaxed, recently reionized regions compared to those that were reionized much earlier. Assessing the impact of this response on the global reionizaton process, we find that un-relaxed gaseous structures boost the total number of recombinations by $\approx 50$ % and lead to spatial fluctuations in the mean free path that persist appreciably for several hundred million years after the completion of reionization.

Massive stars are often found in binary systems and it has been argued that binary products boost the ionizing radiation of stellar populations. Accurate predictions for binary products are needed to understand and quantify their contribution to Cosmic Reionization. We investigate the contribution of stars stripped in binaries since (1) they are, arguably, the best-understood products of binary evolution, (2) we recently produced the first non-LTE radiative transfer calculations for the atmospheres of these stripped stars that predict their ionizing spectra, and (3) they are very promising sources since they boost the ionizing emission of stellar populations at late times. This allows stellar feedback to clear the surroundings such that a higher fraction of their photons can escape and ionize the intergalactic medium. Combining our detailed predictions for the ionizing spectra with a simple cosmic reionization model, we estimate that stripped stars contributed tens of percent of the photons that caused cosmic reionization of hydrogen, depending on the assumed escape fractions. More importantly, stripped stars harden the ionizing emission. We estimate that the spectral index for the ionizing part of the spectrum can increase to -1 compared to <-2 for single stars. At high redshift, stripped stars and massive single stars combined dominate the HeII-ionizing emission, but we expect active galactic nuclei drive cosmic helium reionization. Further observational consequences we expect are (1) high ionization states for the intergalactic gas surrounding stellar systems, such as CIV and SiIV and (2) additional heating of the intergalactic medium of up to a few thousand Kelvin. Quantifying these warrants the inclusion of accurate models for stripped stars and other binary products in full cosmological simulations.

Spatial fluctuations in ultraviolet backgrounds can subtly modulate the distribution of extragalactic sources, a potential signal and systematic for large-scale structure surveys. While this modulation has been shown to be significant for 3D Ly alpha forest surveys, its relevance for other large-scale structure probes has been hardly explored, despite being the only astrophysical process that likely can affect clustering measurements on the scales of >~Mpc. We estimate that the background fluctuations, modulating the amount of H I, have a fractional effect of (0.03-0.3) x (k/[10^-2 Mpc^-1])^-1 on the power spectrum of 21 cm intensity maps at z = 1-3. We find a smaller effect for H alpha and Ly alpha intensity mapping surveys of (0.001-0.1) x (k/[10^-2 Mpc^-1])-1 and even smaller effect for more traditional surveys that correlate the positions of individual H alpha or Ly alpha emitters. We also estimate the effect of backgrounds on low-redshift galaxy surveys in general based on a simple model in which background fluctuations modulate the rate halo gas cools, modulating star formation: We estimate a maximum fractional effect on the power of ~0.01 (k/[10-2 Mpc-1])-1 at z = 1. We compare sizes of these imprints to cosmological parameter benchmarks for the next generation of redshift surveys: We find that ionizing backgrounds could result in a bias on the squeezed triangle non-Gaussianity parameter fNL that can be larger than unity for power spectrum measurements with a SPHEREx-like galaxy survey, and typical values of intensity bias. Marginalizing over a shape of the form k-1PL, where PL is the linear matter power spectrum, removes much of this bias at the cost of ~ 40% larger statistical errors.

Reionization leads to large spatial fluctuations in the intergalactic temperature that can persist well after its completion. We study the imprints of such fluctuations on the $z\sim5$ Ly$\alpha$ forest flux power spectrum using a set of radiation-hydrodynamic simulations that model different reionization scenarios. We find that large-scale coherent temperature fluctuations bring $\sim20-60\%$ extra power at $k\sim0.002$ s/km, with the largest enhancements in the models where reionization is extended or ends the latest. On smaller scales ($k\gtrsim0.1$ s/km), we find that temperature fluctuations suppress power by $\lesssim10\%$. We find that the shape of the power spectrum is mostly sensitive to the reionization midpoint rather than temperature fluctuations from reionization's patchiness. However, for all of our models with reionization midpoints of $z\le 8$ ($z\le 12$) the shape differences are $\lesssim20\%$ ($\lesssim40\%$) because of a surprisingly well-matched cancellation between thermal broadening and pressure smoothing that occurs for realistic thermal histories. We also consider fluctuations in the ultraviolet background, finding their impact on the power spectrum to be much smaller than temperature fluctuations at $k\gtrsim0.01$ s/km. Furthermore, we compare our models to power spectrum measurements, finding that none of our models with reionization midpoints of $z<8$ is strongly preferred over another and that all of our models with midpoints of $z\geq8$ are excluded at $2.5\sigma$. Future measurements may be able to distinguish between viable reionization models if they can be performed at lower $k$ or, alternatively, if the error bars on the high-$k$ power can be reduced by a factor of $1.5$.

Present-day galaxies are surrounded by cool and enriched halo gas extending to hundreds of kiloparsecs. This halo gas is thought to be the dominant reservoir of material available to fuel future star formation, but direct constraints on its mass and physical properties have been difficult to obtain. We report the detection of a fast radio burst (FRB 181112) with arcsecond precision, which passes through the halo of a foreground galaxy. Analysis of the burst shows the halo gas has low net magnetization and turbulence. Our results imply predominantly diffuse gas in massive galactic halos, even those hosting active supermassive black holes, contrary to some previous results.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

Fast Radio Bursts (FRBs) are brief radio emissions from distant astronomical sources. Some are known to repeat, but most are single bursts. Non-repeating FRB observations have had insufficient positional accuracy to localize them to an individual host galaxy. We report the interferometric localization of the single pulse FRB 180924 to a position 4 kpc from the center of a luminous galaxy at redshift 0.3214. The burst has not been observed to repeat. The properties of the burst and its host are markedly different from the only other accurately localized FRB source. The integrated electron column density along the line of sight closely matches models of the intergalactic medium, indicating that some FRBs are clean probes of the baryonic component of the cosmic web.

During reionization, the intergalactic medium is heated impulsively by supersonic ionization fronts (I-fronts). The peak gas temperatures behind the I-fronts, $T_\mathrm{reion}$, are a key uncertainty in models of the thermal history after reionization. Here we use high-resolution radiative transfer simulations to study the parameter space of $T_\mathrm{reion}$. We show that $T_\mathrm{reion}$ is only mildly sensitive to the spectrum of incident radiation over most of the parameter space, with temperatures set primarily by I-front speeds. We also explore what current models of reionization predict for $T_\mathrm{reion}$ by measuring I-front speeds in cosmological radiative transfer simulations. We find that the post-I-front temperatures evolve toward hotter values as reionization progresses. Temperatures of $T_\mathrm{reion} = 17,000-22,000$ K are typical during the first half of reionization, but $T_\mathrm{reion} = 25,000 - 30,000$ K may be achieved near the end of this process if I-front speeds reach $\sim10^4$ km/s as found in our simulations. Shorter reionization epochs lead to hotter $T_\mathrm{reion}$. We discuss implications for $z>5$ Ly$\alpha$ forest observations, which potentially include sight lines through hot, recently reionized patches of the Universe. Interpolation tables from our parameter space study are made publicly available, along with a simple fit for the dependence of $T_\mathrm{reion}$ on the I-front speed.

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.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

The discovery of Fast Radio Bursts (FRBs) at cosmological distances has opened a powerful window on otherwise unseen matter in the Universe. In the 2020s, observations of $>10^{4}$ FRBs will assess the baryon contents and physical conditions in the hot/diffuse circumgalactic, intracluster, and intergalactic medium, and test extant compact-object dark matter models.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

Achieving a precise understanding of galaxy formation in a cosmological context is one of the great challenges in theoretical astrophysics, due to the vast range of spatial scales involved in the relevant physical processes. Observations in the millimeter bands, particularly those using the cosmic microwave background (CMB) radiation as a "backlight", provide a unique probe of the thermodynamics of these processes, with the capability to directly measure the density, pressure, and temperature of ionized gas. Moreover, these observations have uniquely high sensitivity into the outskirts of the halos of galaxies and clusters, including systems at high redshift. In the next decade, the combination of large spectroscopic and photometric optical galaxy surveys and wide-field, low-noise CMB surveys will transform our understanding of galaxy formation via these probes.

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.

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.

Radiation in the extreme ultraviolet (EUV) and soft X-ray holds clues to the location of the missing baryons, the energetics in stellar feedback processes, and the cosmic enrichment history. Additionally, EUV and soft X-ray photons help determine the ionization state of most intergalactic and circumgalactic metals, shaping the rate at which cosmic gas cools. Unfortunately, this band is extremely difficult to probe observationally due to absorption from the Galaxy. In this paper, we model the contributions of various sources to the cosmic EUV and soft X-ray backgrounds. We bracket the contribution from (1) quasars, (2) X-ray binaries, (3) hot interstellar gas, (4) circumgalactic gas, (5) virialized gas, and (6) supersoft sources, developing models that extrapolate into these bands using both empirical and theoretical inputs. While quasars are traditionally assumed to dominate these backgrounds, we discuss the substantial uncertainty in their contribution. Furthermore, we find that hot intrahalo gases likely emit an O(1) fraction of this radiation at low redshifts, and that interstellar and circumgalactic emission potentially contribute tens of percent to these backgrounds at all redshifts. We estimate that uncertainties in the angular-averaged background intensity impact the ionization corrections for common circumgalactic and intergalactic metal absorption lines by ~0.3-1 dex, and we show that local emissions are comparable to the cosmic background only at r_prox = 10-100 kpc from Milky Way-like galaxies.

We develop an effective perturbation theory (and, equivalently, a bias expansion) for the inhomogeneous 21cm radiation field from reionization. Using large-scale simulations of cosmological reionization, we find that this expansion describes the modes in the simulated 21cm signal over much of the wavenumber range probed by upcoming 21cm arrays. This result provides an understanding of the potential signal shapes that compliments the nonlinear numerical modeling that has been the focus of most previous work. We find that the observable signal often can be described with 2-3 bias coefficients that can be interpreted in terms of the source biases, the average neutral fraction, the characteristic size of ionized regions, and the patchiness of reionization. The 21cm signal serves as a wacky example of a bias expansion in cosmology, with our approach synthesizing key results.

We present a formalism for computing time-dependent fluctuations in the cosmological photoionizing radiation background, extending background fluctuations models beyond the steady-state approximation. We apply this formalism to estimate fluctuations in the HI LyA flux redshift space power spectrum and its spatial correlation function at redshifts 2<z<4, assuming the photoionization background is dominated by Quasi-stellar Objects (QSOs) and/or galaxies. We show the shot noise in the power spectrum due to discrete sources is strongly suppressed relative to the steady-state value at low wavenumbers by a factor proportional to the lifetime of the sources, and that this suppression may be used to constrain QSO lifetimes. The total HI LyA power spectrum including shot noise is affected at tens of percent on short scales, and by as much as an order of magnitude or more on scales exceeding the mean free path. The spatial correlation function is similarly found to be sensitive to the shot noise, although moderately insensitive to the effects of time-dependence on the non-shotnoise contribution. Photoionization rate fluctuations substantially modify the shape of the Baryonic Acoustic Oscillation peak in the correlation function, including a small increase in its position that must be accounted for to avoid biasing estimates of cosmological parameters based on the peak position. We briefly investigate solving the full frequency dependent equation, finding that it agrees with the frequency-independent to better than percent accuracy. Simple formulas are provided for the power spectrum of fluctuations in the photoionization rate that approximate the full computations.

Cosmological simulations of the low-density intergalactic medium exhibit a strikingly tight power-law relation between temperature and density that holds over two decades in density. It is found that this relation should roughly apply Delta z ~ 1-2 after a reionization event, and this limiting behavior has motivated the power-law parameterizations used in most analyses of the Ly-alpha forest. This relation has been explained by using equations linearized in the baryonic overdensity (which does not address why a tight power-law relation holds over two decades in density) or by equating the photoheating rate with the cooling rate from cosmological expansion (which we show is incorrect). Previous explanations also did not address why recombination cooling and Compton cooling off of the cosmic microwave background, which are never negligible, do not alter the character of this relation. We provide an understanding for why a tight power-law relation arises for unshocked gas at all densities for which collisional cooling is unimportant. We also use our results to comment on (1) how quickly fluctuations in temperature redshift away after reionization processes, (2) how much shock heating occurs in the low-density intergalactic medium, and (3) how the temperatures of collapsing gas parcels evolve.

We present the first simulations of the high-redshift Ly$\alpha$ intensity field that account for scattering in the intergalactic medium (IGM). Using a 3D Monte Carlo radiative transfer code, we find that Ly$\alpha$ scattering smooths spatial fluctuations in the Ly$\alpha$ intensity on small scales and that the spatial dependence of this smoothing depends strongly on the mean neutral fraction of the IGM. Our simulations find a strong effect of reionization on $k=0.1-1~{\rm Mpc^{-1}}$, with $P_{\rm Ly\alpha}\propto k^{-1.75}$ for $\bar{x}_{\rm HI} = 0.63$ and $P_{\rm Ly\alpha} \propto k^{-2.2}$ for $\bar{x}_{\rm HI} = 0.86$ in contrast to $P_{\rm Ly\alpha}\propto k^{-1.5}$ after reionization. At wavenumbers of $k>1 ~ {\rm Mpc^{-1}}$, we find that the signal is sensitive to the emergent Ly$\alpha$ line profiles from galaxies. We also demonstrate that the cross-correlation between a Ly$\alpha$ intensity map and a future galaxy redshift survey could be detected on large scales by an instrument similar to SPHEREx, and over a wide range of scales by a hypothetical intensity mapping instrument in the vein of CDIM.

A component of the dark matter could consist of two darkly charged particles with a large mass ratio and a massless force carrier. This `atomic' dark sector could behave much like the baryonic sector, cooling and fragmenting down to stellar-mass or smaller scales. Past studies have shown that cosmic microwave background and large-scale structure constraints rule out $\gtrsim 5\%$ of the dark matter to behave in this manner. However, we show that, even with percent level mass fractions, a dark atomic sector could affect some extragalactic and galactic observables. We track the cooling and merger history of an atomic dark component for much of the interesting parameter space. Unlike the baryons, where stellar feedback (driven by nuclear physics) delays the formation and growth of galaxies, cooling dark atomic gas typically results in disks forming earlier, leaving more time for their destruction via mergers. Rather than disks in Milky Way sized halos, we find the end product is typically spheroidal structures on galactic scales or dark atom fragments distributed on halo scales. This result contrasts with previous studies, which had assumed that the dark atoms would result in dark disks. Furthermore the dark atoms condense into dense clumps, analogous to how the baryons fragment on solar-mass scales. We estimate the size of these dark clumps, and use these estimates to show that viable atomic dark matter parameter space is ruled out by stellar microlensing, by the half-light radii of ultra-faint dwarf galaxies, and by Milky Way mass-to-light inferences.

We present a semi-analytic model for the Lyman-$\alpha$ forest that is inspired by the Halo Model. This model is built on the absorption line decomposition of the forest. Flux correlations are decomposed into those within each absorption line (the 1-absorber term) and those between separate lines (the 2-absorber term), treating the lines as biased tracers of the underlying matter fluctuations. While the nonlinear exponential mapping between optical depth and flux requires an infinite series of moments to calculate any statistic, we show that this series can be resumed (truncating at the desired order in the linear matter overdensity). We focus on the $z=2-3$ line-of-sight power spectrum. Our model finds that 1-absorber term dominates the power on all scales, with most of its contribution coming from HI columns of $10^{14}-10^{15}\;\mathrm{cm^{-2}}$, while the smaller 2-absorber contribution comes from lower columns that trace overdensities of a few. The prominence of the 1-absorber correlations indicates that the line-of-sight power spectrum is shaped principally by the lines' number densities and their absorption profiles, with correlations between lines contributing to a lesser extent. We present intuitive formulae for the effective optical depth as well as the large-scale limits of 1-absorber and 2-absorber terms, which simplify to integrals over the HI column density distribution with different equivalent-width weightings. With minimalist models for the bias of absorption systems and their peculiar velocity broadening, our model predicts values for the density bias and velocity gradient bias that are consistent with those found in simulations.

Motivated by the claimed detection of a large population of faint active galactic nuclei (AGN) at high redshift, recent studies have proposed models in which AGN contribute significantly to the z > 4 H I ionizing background. In some models, AGN are even the chief sources of reionization. If correct, these models would make necessary a complete revision to the standard view that galaxies dominated the high-redshift ionizing background. It has been suggested that AGN-dominated models can better account for two recent observations that appear to be in conflict with the standard view: (1) large opacity variations in the z ~ 5.5 H I Lyman-alpha forest, and (2) slow evolution in the mean opacity of the He II Lyman-alpha forest. Large spatial fluctuations in the ionizing background from the brightness and rarity of AGN may account for the former, while the earlier onset of He II reionization in these models may account for the latter. Here we show that models in which AGN emissions source >~ 50 % of the ionizing background generally provide a better fit to the observed H I Lyman-alpha forest opacity variations compared to standard galaxy-dominated models. However, we argue that these AGN-dominated models are in tension with constraints on the thermal history of the intergalactic medium (IGM). Under standard assumptions about the spectra of AGN, we show that the earlier onset of He II reionization heats up the IGM well above recent temperature measurements. We further argue that the slower evolution of the mean opacity of the He II Lyman-alpha forest relative to simulations may reflect deficiencies in current simulations rather than favor AGN-dominated models as has been suggested.

Recent observations have shown that the scatter in opacities among coeval segments of the Lyman-alpha forest increases rapidly at z > 5. In this paper, we assess whether the large scatter can be explained by fluctuations in the ionizing background in the post-reionization intergalactic medium. We find that matching the observed scatter at z ~ 5.5 requires a short spatially averaged mean free path of < 15 comoving Mpc/h, a factor of > 3 shorter than direct measurements at z ~ 5.2. We argue that such rapid evolution in the mean free path is difficult to reconcile with our measurements of the global H I photoionization rate, which stay approximately constant over the interval z ~ 4.8 - 5.5. However, we also show that measurements of the mean free path at z > 5 are likely biased towards higher values by the quasar proximity effect. This bias can reconcile the short values of the mean free path that are required to explain the large scatter in opacities. We discuss the implications of this scenario for cosmological reionization. Finally, we investigate whether other statistics applied to the z > 5 Lyman-alpha forest can shed light on the origin of the scatter. Compared to a model with a uniform ionizing background, models that successfully account for the scatter lead to enhanced power in the line-of-sight flux power spectrum on scales k < 0.1 h/Mpc. We find tentative evidence for this enhancement in observations of the high-redshift Lyman-alpha forest.

Observations reveal massive amounts of OVI around star-forming $L_*$ galaxies, with covering fractions of near unity extending to the host halo's virial radius. This OVI absorption is typically kinematically centered upon photoionized gas, with line widths that are suprathermal and kinematically offset from the galaxy. We discuss various scenarios and whether they could result in the observed phenomenology (cooling gas flows, boundary layers, shocks, virialized gas, photoionized clouds in thermal equilibrium). If predominantly collisionally ionized, as we argue is most probable, the OVI observations require that the circumgalactic medium (CGM) of $L_*$ galaxies holds nearly all the associated baryons within a virial radius ($\sim 10^{11}M_\odot$) and hosts massive flows of cooling gas with $\approx30[nT/30{\rm~cm^{-3}K}]~M_\odot~$yr$^{-1}$, which must be largely prevented from accreting onto the host galaxy. Cooling and feedback energetics considerations require $10 <\langle nT\rangle<100{\rm~cm^{-3}K}$ for the warm and hot halo gases. We argue that virialized gas, boundary layers, hot winds, and shocks are unlikely to directly account for the bulk of the OVI. Furthermore, we show that there is a robust constraint on the number density of many of the photoionized $\sim10^4$K absorption systems that yields upper bounds in the range $n<(0.1-3)\times10^{-3}(Z/0.3)$cm$^{-3}$, where $Z$ is the metallicity, suggestive that the dominant pressure in some photoionized clouds is nonthermal. This constraint, which requires minimal ionization modeling, is in accord with the low densities inferred from more complex photoionization modeling. The large amount of cooling gas that is inferred could re-form these clouds in a fraction of the halo dynamical time, as some arguments require, and it requires much of the feedback energy available from supernovae and stellar winds to be dissipated in the CGM.

Despite decades of effort, the timing and duration of He II reionization, as well as its morphology and the properties of the quasars believed to drive it, are still not well constrained. In this paper we present a new method to study both He II reionization and quasars via the thermal proximity effect -- the photoelectric heating of the intergalactic medium around quasars when their hard radiation doubly ionizes helium. We post-process a SPH simulation with 1D radiative calculations, and study how the thermal proximity effect depends on the amount of singly ionized helium, $x_{\rm HeII,0}$, which prevailed in the IGM before the quasar turned on, and the characteristic lifetime $t_{\rm Q}$ for which quasars shine. We find that the amplitude of the temperature boost in the quasar environment depends on $x_{\rm HeII,0}$, with a characteristic value of $\Delta T \simeq 10^4\,{\rm K}$ for an initially singly ionized IGM ($x_{\rm HeII,0} = 1.0$), whereas the size of the thermal proximity zone is sensitive to quasar lifetime $t_{\rm Q}$, with typical sizes of ~100 cMpc for luminous quasars shining for $t_{\rm Q}=10^8$ yr. This temperature boost is manifest as a change in the thermal broadening of H I absorption lines near the quasar. We introduce a new method based on measuring the Ly$\alpha$ forest power spectrum as a function of distance from the quasar, and conduct Bayesian MCMC analysis to demonstrate that the thermal proximity effect should be easily detectable. For a mock dataset of 50 quasars at z~4, we predict that one can measure $x_{\rm HeII,0}$ to a precision $\approx 0.04$, and $t_{\rm Q}$ to a precision of $\approx 0.1$ dex. By applying our formalism to existing high-resolution Ly$\alpha$ forest spectra of quasars at $3.1 \lesssim z \lesssim 5.0$, one should be able to detect the thermal proximity effect, and reconstruct the full reionization history of He II.

The total contribution of diffuse halo gas to the galaxy baryon budget strongly depends on its dominant ionization state. In this paper, we address the physical conditions in the highly ionized circumgalactic medium (CGM) traced by OVI absorption lines observed in COS-Halos spectra. We analyze the observed ionic column densities, absorption-line widths and relative velocities, along with the ratios of NV/OVI for 39 fitted Voigt profile components of OVI. We compare these quantities with the predictions given by a wide range of ionization models. Photoionization models that include only extragalactic UV background radiation are ruled out; conservatively, the upper limits to NV/OVI and measurements of N$_{\rm OVI}$ imply unphysically large path lengths > 100 kpc. Furthermore, very broad OVI absorption (b > 40 km s$^{-1}$) is a defining characteristic of the CGM of star-forming L* galaxies. We highlight two possible origins for the bulk of the observed OVI: (1) highly structured gas clouds photoionized primarily by local high energy sources or (2) gas radiatively cooling on large scales behind a supersonic wind. Approximately 20% of circumgalactic OVI does not align with any low-ionization state gas within $\pm$50 km s$^{-1}$ and is found only in halos with M$_{\rm halo}$ < 10$^{12}$ M$_{\odot}$. We suggest that this type of unmatched OVI absorption traces the hot corona itself at a characteristic temperature of 10$^{5.5}$ K. We discuss the implications of these very distinct physical origins for the dynamical state, gas cooling rates, and total baryonic content of L* gaseous halos.

We investigate how well the redshift distribution of a population of extragalactic objects can be reconstructed using angular cross-correlations with a sample whose redshifts are known. We derive the minimum variance quadratic estimator, which has simple analytic representations in very applicable limits and is significantly more sensitive than earlier proposed estimation procedures. This estimator is straightforward to apply to observations, it robustly finds the likelihood maximum, and it conveniently selects angular scales at which fluctuations are well approximated as independent between redshift bins and at which linear theory applies. We find that the linear bias times number of objects in a redshift bin generally can be constrained with cross-correlations to fractional error (10^2 n/N)^1/2, where N is the total number of spectra per dz and n is the number of redshift bins spanned by the bulk of the unknown population. The error is often independent of the sky area and sampling fraction. Furthermore, we find that sub-percent measurements of the angular source density per unit redshift, dN/dz, are in principle possible, although cosmic magnification needs to be accounted for at fractional errors of <~ 10 per cent. We discuss how the sensitivity to dN/dz changes as a function of photometric and spectroscopic depth and how to optimize the survey strategy to constrain dN/dz. We also quantify how well cross-correlations of photometric redshift bins can be used to self-calibrate a photometric redshift sample. Simple formulae that can be quickly applied to gauge the utility of cross correlating different samples are given.

The bulk of cosmic matter resides in a dilute reservoir that fills the space between galaxies, the intergalactic medium (IGM). The history of this reservoir is intimately tied to the cosmic histories of structure formation, star formation, and supermassive black hole accretion. Our models for the IGM at intermediate redshifts (2<z<5) are a tremendous success, quantitatively explaining the statistics of Lyman-alpha absorption of intergalactic hydrogen. However, at both lower and higher redshifts (and around galaxies) much is still unknown about the IGM. We review the theoretical models and measurements that form the basis for the modern understanding of the IGM, and we discuss unsolved puzzles (ranging from the largely unconstrained process of reionization at high-z to the missing baryon problem at low-z), highlighting the efforts that have the potential to solve them.

Recent years have brought more precise temperature measurements of the low-density intergalactic medium (IGM). These new measurements constrain the processes that heated the IGM, such as the reionization of H I and of He II. We present a semi-analytical model for the thermal history of the IGM that follows the photoheating history of primordial gas. Our model adopts a multizone approach that, compared to previous models, more accurately captures the inhomogeneous heating and cooling of the IGM during patchy reionization. We compare our model with recent temperature measurements spanning z = 1.6-4.8, finding that these measurements are consistent with scenarios in which the He II was reionized at z = 3-4 by quasars. Significantly longer duration or higher redshift He II reionization scenarios are ruled out by the measurements. For hydrogen reionization, we find that only low-redshift and high-temperature scenarios are excluded. For example, a model in which the IGM was heated to 30 000 K when an ionization front passed, and with hydrogen reionization occurring over 6 < z < 9, is ruled out. Finally, we place constraints on how much heating could owe to TeV blazars, cosmic rays, and other non-standard mechanisms. We find that by z = 2, a maximum of 1 eV of additional heat could be injected per baryon over standard photoheating-only models, with this limit becoming <~0.5 eV at z > 3.

The lifetime of quasars is fundamental for understanding the growth of supermassive black holes, and is an important ingredient in models of the reionization of the intergalactic medium. However, despite various attempts to determine quasar lifetimes, current estimates from a variety of methods are uncertain by orders of magnitude. This work combines cosmological hydrodynamical simulations and 1D radiative transfer to investigate the structure and evolution of the He II Ly$\alpha$ proximity zones around quasars at $z \simeq 3-4$. We show that the time evolution in the proximity zone can be described by a simple analytical model for the approach of the He II fraction $x_{\rm HeII}\left( t \right)$ to ionization equilibrium, and use this picture to illustrate how the transmission profile depends on the quasar lifetime, quasar UV luminosity, and the ionization state of helium in the ambient IGM (i.e. the average He II fraction, or equivalently the metagalactic He II ionizing background). A significant degeneracy exists between the lifetime and the average He II fraction, however the latter can be determined from measurements of the He II Ly$\alpha$ optical depth far from quasars, allowing the lifetime to be measured. We advocate stacking existing He II quasar spectra at $z\sim 3$, and show that the shape of this average proximity zone profile is sensitive to lifetimes as long as $\sim 30$ Myr. At higher redshift $z\sim 4$ where the He II fraction is poorly constrained, degeneracies will make it challenging to determine these parameters independently. Our analytical model for He II proximity zones should also provide a useful description of the properties of H I proximity zones around quasars at $z \simeq 6-7$.

We measure the effective optical depth of HeII Ly\alpha\ absorption \tau$_\mathsf{eff,HeII}$ at 2.3<z<3.5 in 17 UV-transmitting quasars observed with UV spectrographs on the Hubble Space Telescope (HST). The median \tau$_\mathsf{eff,HeII}$ values increase gradually from 1.95 at z=2.7 to 5.17 at z=3.4, but with a strong sightline-to-sightline variance. Many $\simeq$35 comoving Mpc regions of the z>3 intergalactic medium (IGM) remain transmissive (\tau$_\mathsf{eff,HeII}$<4), and the gradual trend with redshift appears consistent with density evolution of a fully reionized IGM. These modest optical depths imply average HeII fractions of x$_\mathsf{HeII}$<0.01 and HeII ionizing photon mean free paths of $\simeq$50 comoving Mpc at z$\simeq$3.4, thus requiring that a substantial volume of the helium in the Universe was already doubly ionized at early times; this stands in conflict with current models of HeII reionization driven by luminous quasars. Along 10 sightlines we measure the coeval HI Ly\alpha\ effective optical depths, allowing us to study the density dependence of \tau$_\mathsf{eff,HeII}$ at z$\sim$3. We establish that the dependence of \tau$_\mathsf{eff,HeII}$ on increasing \tau$_\mathsf{eff,HI}$ is significantly shallower than expected from simple models of an IGM reionized in HeII. This requires higher HeII photoionization rates in overdense regions or underdense regions being not in photoionization equilibrium. Moreover, there are very large fluctuations in \tau$_\mathsf{eff,HeII}$ at all $\tau_\mathsf{eff,HI}$, which greatly exceed the expectations from these simple models. These data present a distinct challenge to scenarios of HeII reionization - an IGM where HeII appears to be predominantly ionized at z$\simeq$3.4, and with a radiation field strength that may be correlated with the density field, but exhibits large fluctuations at all densities.

Many recent studies have highlighted certain failures of the standard Eulerian-space cosmological perturbation theory (SPT). Its problems include (1) not capturing large-scale bulk flows [leading to an O(1) error in the 1-loop SPT prediction for the baryon acoustic peak in the correlation function], (2) assuming that the Universe behaves as a pressureless, inviscid fluid, and (3) treating fluctuations on scales that are non-perturbative as if they were. Recent studies have highlighted the successes of perturbation theory in Lagrangian space or theories that solve equations for the effective dynamics of smoothed fields. Both approaches mitigate some or all of the aforementioned issues with SPT. We discuss these physical developments by specializing to the simplified 1D case of gravitationally interacting sheets, which allows us to substantially reduces the analytic overhead and still (as we show) maintain many of the same behaviors as in 3D. In 1D, linear-order Lagrangian perturbation theory ("the Zeldovich approximation") is exact up to shell crossing, and we prove that n^{th}-order Eulerian perturbation theory converges to the Zeldovich approximation as n goes to infinity. In no 1D cosmology that we consider (including a CDM-like case and power-law models) do these theories describe accurately the matter power spectrum on any mildly nonlinear scale. We find that theories based on effective equations are much more successful at describing the dynamics. Finally, we discuss many topics that have recently appeared in the perturbation theory literature such as beat coupling, the shift and smearing of the baryon acoustic oscillation feature, and the advantages of Fourier versus configuration space. Our simplified 1D case serves as an intuitive review of these perturbation theory results.

Recent observations of the Lyman-alpha forest show large-scale spatial variations in the intergalactic Lyman-alpha opacity that grow rapidly with redshift at z>5, far in excess of expectations from empirically motivated models. Previous studies have attempted to explain this excess with spatial fluctuations in the ionizing background, but found that this required either extremely rare sources or problematically low values for the mean free path of ionizing photons. Here we report that much -- or potentially all -- of the observed excess likely arises from residual spatial variations in temperature that are an inevitable byproduct of a patchy and extended reionization process. The amplitude of opacity fluctuations generated in this way depends on the timing and duration of reionization. If the entire excess is due to temperature variations alone, the observed fluctuation amplitude favors a late-ending but extended reionization process that was roughly half complete by z~9 and that ended at z~6. In this scenario, the highest opacities occur in regions that reionized earliest, since they have had the most time to cool, while the lowest opacities occur in the warmer regions that reionized most recently. This correspondence potentially opens a new observational window into patchy reionization.

We present constraints on both the kinetic temperature of the intergalactic medium (IGM) at z=8.4, and on models for heating the IGM at high-redshift with X-ray emission from the first collapsed objects. These constraints are derived using a semi-analytic method to explore the new measurements of the 21 cm power spectrum from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), which were presented in a companion paper, Ali et al. (2015). Twenty-one cm power spectra with amplitudes of hundreds of mK^2 can be generically produced if the kinetic temperature of the IGM is significantly below the temperature of the Cosmic Microwave Background (CMB); as such, the new results from PAPER place lower limits on the IGM temperature at z=8.4. Allowing for the unknown ionization state of the IGM, our measurements find the IGM temperature to be above ~5 K for neutral fractions between 10% and 85%, above ~7 K for neutral fractions between 15% and 80%, or above ~10 K for neutral fractions between 30% and 70%. We also calculate the heating of the IGM that would be provided by the observed high redshift galaxy population, and find that for most models, these galaxies are sufficient to bring the IGM temperature above our lower limits. However, there are significant ranges of parameter space that could produce a signal ruled out by the PAPER measurements; models with a steep drop-off in the star formation rate density at high redshifts or with relatively low values for the X-ray to star formation rate efficiency of high redshift galaxies are generally disfavored. The PAPER measurements are consistent with (but do not constrain) a hydrogen spin temperature above the CMB temperature, a situation which we find to be generally predicted if galaxies fainter than the current detection limits of optical/NIR surveys are included in calculations of X-ray heating.

Poorly understood feedback processes associated with highly-luminous black hole accretion in quasars may dramatically affect the properties of their host galaxies. We search for the effect of quasar feedback on surrounding gas using Planck maps of the thermal Sunyaev-Zel'dovich effect (tSZ). By stacking tSZ Compton-y maps centered on the locations of 26,686 spectroscopic quasars from the Sloan Digital Sky Survey, we detect a strong but unresolved tSZ Compton-y signal at >5 sigma significance that likely originates from a combination of virialized halo atmosphere gas and quasar feedback effects. We show that the feedback contribution to our detected quasar tSZ signal is likely to dominate over virialized halo gas by isolating the feedback tSZ component for high- and low-redshift quasars. We find that this quasar tSZ signal also scales with black hole mass and bolometric luminosity, all consistent with general expectations of quasar feedback. We estimate the mean angularly-integrated Compton-y of quasars at z~1.5 to be 3.5x10^-6 Mpc^2, corresponding to mean total thermal energies in feedback and virialized halo gas of 1.1(+/- 0.2) x 10^62 erg, and discuss the implications for quasar feedback. If confirmed, the large total thermal feedback energetics we estimate of 5% (+/-1% statistical uncertainty) of the black hole mass will have important implications for the effects of quasar feedback on the host galaxy, as well as the surrounding intergalactic medium.

Impulsive radio bursts that are detectable across cosmological distances constitute extremely powerful probes of the ionized Inter-Galactic Medium (IGM), intergalactic magnetic fields, and the properties of space-time itself. Their dispersion measures (DMs) will enable us to detect the "missing" baryons in the low-redshift Universe and make the first measurements of the mean galaxy halo profile, a key parameter in models of galaxy formation and feedback. Impulsive bursts can be used as cosmic rulers at redshifts exceeding 2, and constrain the dark energy equation-of-state parameter, $w(z)$ at redshifts beyond those readily accessible by Type Ia SNe. Both of these goals are realisable with a sample of $\sim 10^4$ fast radio bursts (FRBs) whose positions are localized to within one arcsecond, sufficient to obtain host galaxy redshifts via optical follow-up. It is also hypothesised that gravitational wave events may emit coherent emission at frequencies probed by SKA1-LOW, and the localization of such events at cosmological distances would enable their use as cosmological standard sirens. To perform this science, such bursts must be localized to their specific host galaxies so that their redshifts may be obtained and compared against their dispersion measures, rotation measures, and scattering properties. The SKA can achieve this with a design that has a wide field-of-view, a substantial fraction of its collecting area in a compact configuration (80\% within a 3\,km radius), and a capacity to attach high-time-resolution instrumentation to its signal path.

The next frontier of Lyman-alpha forest studies is the reconstruction of 3D correlations from a dense sample of background sources. The measurement of 3D correlations has the potential to improve constraints on fundamental cosmological parameters, ionizing background models, and the reionization history. This study addresses the sensitivity of spectroscopic surveys to 3D correlations in the Lyman-alpha forest. We show that the sensitivity of a survey to this signal can be quantified by just a single number, a noise-weighted number density of sources on the sky. We investigate how the sensitivity of a spectroscopic quasar (or galaxy) survey scales as a function of its depth, area, and redshift. We propose a simple method for weighting sightlines with varying S/N levels to estimate the correlation function, and we show that this estimator generally performs nearly as well as the minimum variance quadratic estimator. In addition, we show that the sensitivity of a quasar survey to the flux correlation function is generally maximized if it observes each field just long enough to achieve S/N ~ 2 in a 1 A pixel on an L_* quasar while acquiring spectra for all quasars with L > L_*: Little is gained by integrating longer on the same targets or by including fainter quasars. We quantify how these considerations relate to constraints on the angular diameter distance, the curvature of space-time, and the reionization history.

The most promising near-term observable of the cosmic dark age prior to widespread reionization (z~15-200) is the sky-averaged \lambda 21 cm background arising from hydrogen in the intergalactic medium. Though an individual antenna could in principle detect the line signature, data analysis must separate foregrounds that are orders of magnitude brighter than the \lambda 21 cm background (but that are anticipated to vary monotonically and gradually with frequency). Using more physically motivated models for foregrounds than in previous studies, we show that the intrinsic "spectral smoothness" of the foregrounds is likely not a concern, and that data analysis for an ideal antenna should be able to detect the \lambda 21 cm signal after deprojecting a ~5th order polynomial in log(\nu). However, we find that the foreground signal is corrupted by the frequency-dependent response of a real antenna. The frequency dependence complicates modeling of foregrounds commonly based on the assumption of spectral smoothness. Much of our study focuses on the Large-aperture Experiment to detect the Dark Age (LEDA), which combines both radiometric and interferometric measurements. We show that statistical uncertainty remaining after fitting antenna gain patterns to interferometric measurements does not compromise extraction of the \lambda 21 cm signal for a range of cosmological models after fitting a 7th order polynomial to radiometric data. Our results generalize to most efforts to measure the sky-averaged spectrum.

A number of experiments are currently working towards a measurement of the 21 cm signal from the Epoch of Reionization. Whether or not these experiments deliver a detection of cosmological emission, their limited sensitivity will prevent them from providing detailed information about the astrophysics of reionization. In this work, we consider what types of measurements will be enabled by a next-generation of larger 21 cm EoR telescopes. To calculate the type of constraints that will be possible with such arrays, we use simple models for the instrument, foreground emission, and the reionization history. We focus primarily on an instrument modeled after the $\sim 0.1~\rm{km}^2$ collecting area Hydrogen Epoch of Reionization Array (HERA) concept design, and parameterize the uncertainties with regard to foreground emission by considering different limits to the recently described "wedge" footprint in k-space. Uncertainties in the reionization history are accounted for using a series of simulations which vary the ionizing efficiency and minimum virial temperature of the galaxies responsible for reionization, as well as the mean free path of ionizing photons through the IGM. Given various combinations of models, we consider the significance of the possible power spectrum detections, the ability to trace the power spectrum evolution versus redshift, the detectability of salient power spectrum features, and the achievable level of quantitative constraints on astrophysical parameters. Ultimately, we find that $0.1~\rm{km}^2$ of collecting area is enough to ensure a very high significance ($\gtrsim30\sigma$) detection of the reionization power spectrum in even the most pessimistic scenarios. This sensitivity should allow for meaningful constraints on the reionization history and astrophysical parameters, especially if foreground subtraction techniques can be improved and successfully implemented.

The conditions within the ionization front of a quasar during reionization (T ~ 30,000 K, neutral hydrogen fraction ~ 0.5) are ideal for producing Lyman-alpha emission via collisional excitation of hydrogen atoms. Observations of this emission, which could subtend >10 arcmin$^2$ on the sky, would definitively demonstrate the presence of a neutral intergalactic medium at the observed epoch, placing valuable constraints on the progress of reionization. We find that the expected Lyman-alpha surface brightness is significantly weaker than previously determined and may be impossible to observe with current and near-future instruments. Past work calculated the Lyman-alpha emission from a quasar ionization front in a homogeneous medium with a clumping factor approximation to account for inhomogeneities. We find using 1D radiative transfer calculations that this approximation overestimates the emission by a factor of >3. Our calculations model the propagation of ionizing photons and compute the Lyman-alpha emission from quasar ionization fronts on sightlines from a hydrodynamic cosmological simulation at z = 7.1. To better understand the physical properties of the emission, we also develop an analytic model that accurately describes the results of the full radiative transfer calculation.

We develop and test with cosmological simulations a physically motivated theory for how the interplay between gravity, pressure, cooling, and self-shielding set the redshift--dependent mass scale at which halos can accrete intergalactic gas. This theory provides a physical explanation for the halo mass scale that can accrete unshocked intergalactic gas, which has been explained with ad hoc criteria tuned to reproduce the results of a few simulations. Furthermore, it provides an intuitive explanation for how this mass scale depends on the reionization redshift, the amplitude of the ionizing background, and the redshift. We show that accretion is inhibited onto more massive halos than had been thought because previous studies had focused on the gas fraction of halos rather than the instantaneous mass that can accrete gas. A halo as massive as 10^11 Msun cannot accrete intergalactic gas at z=0, even though typically its progenitors were able to accrete gas at higher redshifts. We describe a simple algorithm that can be implemented in semi-analytic models, and we compare the predictions of this algorithm to numerical simulations.

Previous studies of the 2.2 < z< 2.7 HeII Lyman-alpha forest measured much larger ionizing background fluctuations than are anticipated theoretically. We re-analyze recent Hubble Space Telescope data from the two HeII sightlines that have been used to make these measurements, HE2347-4342 and HS1700+6416, and find that the vast majority of the absorption is actually consistent with a single HeII photoionization rate. We show that the data constrains the RMS fractional fluctuation level smoothed at 1 Mpc to be < 2 and discuss why other studies had found larger fluctuations. Our measurement is consistent with models in which quasars dominate the z=2.5 metagalactic HeII-ionizing background (but it can accommodate less compelling models), and it suggests that quasars (rather than stars) are the dominant contributor to the HI-ionizing background. We detect a HeII transverse proximity effect that is slightly offset in redshift from a known quasar. Its profile and offset may indicate that the quasar turned on 10 Myr ago.

Recently, Thornton and coworkers (2013) confirmed a class of millisecond radio bursts likely of extragalactic origin that is well-suited for estimating dispersion measures (DMs). We calculate the probability distribution of DM(z) in different models for how the cosmic baryons are distributed (both analytically and with cosmological simulations). We show that the distribution of DM is quite sensitive to whether the "missing" baryons lie around the virial radius of 10^11-10^13 Msun halos or further out, which is not easily constrained with other observational techniques. The intrinsic contribution to DM from each source could complicate studies of the extragalactic contribution. This difficulty is avoided by stacking based on the impact parameter to foreground galaxies. We show that a stacking analysis using a sample of ~100 DM measurements from arcminute-localized, z >~ 0.5 sources would place interesting constraints at 0.2-2 halo virial radii on the baryonic mass profile surrounding different galaxy types. Conveniently for intergalactic studies, sightlines that intersect intervening galactic disks should be easily identified owing to scattering. A detectable level of scattering may also result from turbulence in the circumgalactic medium.

Several recent papers have proposed schemes by which a dispersion measure, and hence electron column, could be obtained from a time-steady, incoherent radio source at a cosmological distance (such as an active galactic nucleus). If correct, this would open a new window on the distribution of intergalactic baryons. These schemes are based on the statistical properties of the received radiation, such as the 2- or 4-point correlation function of the received electric field, and in one case on the quantum nature of the electromagnetic field. We show, on the basis of general principles, that these schemes are not sensitive to dispersion measure (or have an extremely small signal-to-noise ratio), because (i) the classical 2-point correlation function is unaffected by dispersion; (ii) for a source with a large number of incoherently emitting electrons, the central limit theorem obliterates additional information in higher-order functions; and (iii) such an emitter produces a radiation density matrix that is equivalent to a statistical distribution of coherent states, which contains no information that is not already in the statistics of the classical waveforms. Why the proposed observables do not depend on dispersion measure (or have extremely tiny dependences) is discussed in detail.

This work describes a new instrument optimized for a detection of the neutral hydrogen 21cm power spectrum between redshifts of 0.5-1.5: the Baryon Acoustic Oscillation Broadband and Broad-beam (BAOBAB) Array. BAOBAB will build on the efforts of a first generation of 21cm experiments which are targeting a detection of the signal from the Epoch of Reionization at z ~ 10. At z ~ 1, the emission from neutral hydrogen in self-shielded overdense halos also presents an accessible signal, since the dominant, synchrotron foreground emission is considerably fainter than at redshift 10. The principle science driver for these observations are Baryon Acoustic Oscillations in the matter power spectrum which have the potential to act as a standard ruler and constrain the nature of dark energy. BAOBAB will fully correlate dual-polarization antenna tiles over the 600-900MHz band with a frequency resolution of 300 kHz and a system temperature of 50K. The number of antennas will grow in staged deployments, and reconfigurations of the array will allow for both traditional imaging and high power spectrum sensitivity operations. We present calculations of the power spectrum sensitivity for various array sizes, with a 35-element array measuring the cosmic neutral hydrogen fraction as a function of redshift, and a 132-element system detecting the BAO features in the power spectrum, yielding a 1.8% error on the z ~ 1 distance scale, and, in turn, significant improvements to constraints on the dark energy equation of state over an unprecedented range of redshifts from ~0.5-1.5.

Recently, Tseliakhovich and Hirata (2010) showed that during the cosmic Dark Ages the baryons were typically moving supersonically with respect to the dark matter with a spatially variable Mach number. Such supersonic motion may source shocks that heat the Universe. This motion may also suppress star formation in the first halos. Even a small amount of coupling of the 21cm signal to this motion has the potential to vastly enhance the 21cm brightness temperature fluctuations at 15<z<40 as well as to imprint acoustic oscillations in this signal. We present estimates for the size of this coupling, which we calibrate with a suite of cosmological simulations. Our simulations, discussed in detail in a companion paper, are initialized to self-consistently account for gas pressure and the dark matter-baryon relative velocity, v_bc (in contrast to prior simulations). We find that the supersonic velocity difference dramatically suppresses structure formation at 10-100 comoving kpc scales, it sources shocks throughout the Universe, and it impacts the accretion of gas onto the first star-forming minihalos (even for halo masses as large as ~10^7 Msun). However, we find that the v_bc-sourced temperature fluctuations can contribute only as much as ~10% of the fluctuations in the 21cm signal. We do find that v_bc could source an O(1) component in the power spectrum of the 21cm signal via the X-ray (but not ultraviolet) backgrounds produced once the first stars formed. In a scenario in which ~10^6 Msun minihalos reheated the Universe via their X-ray backgrounds, we find that the pre-reionization 21cm signal would be larger than previously anticipated and exhibit significant acoustic features. We show that structure formation shocks are unable to heat the Universe sufficiently to erase a strong 21cm absorption trough at z ~ 20 that is found in most models of the sky-averaged 21cm intensity.

We perform a suite of cosmological simulations in the LCDM paradigm of the formation of the first structures in the Universe prior to astrophysical reheating and reionization (15<~z<200). These are the first simulations initialized in a manner that self consistently account for the impact of pressure on the rate of growth of modes, temperature fluctuations in the gas, and the dark matter-baryon supersonic velocity difference. Even with these improvements, these are still difficult times to simulate accurately as the Jeans length of the cold intergalactic gas must be resolved while also capturing a representative sample of the Universe. Our simulations support the finding of recent studies that the dark matter-baryon velocity difference has a surprisingly large impact on the accretion of gas onto the first star-forming minihalos (with masses of ~10^6 Msun). In fact, the halo gas is often significantly downwind of such halos and with lower densities, which delays the formation of the first stars in most locations in the Universe. We also show that dynamical friction plays an important role in the nonlinear evolution of the differential velocity, acting to erase this velocity difference quickly in overdense gas as well as sourcing visually-apparent bow shocks and Mach cones throughout the Universe. We use simulations with both the GADGET and Enzo cosmological codes to test the robustness of these conclusions. We find that particle coupling in GADGET between the gas and dark matter particles can result in spurious growth that mimics nonlinear growth in the matter power spectrum. In a companion paper, we use the simulations presented here to make detailed estimates for the impact of the dark matter--baryon velocity differential on redshifted 21cm radiation. The initial conditions generator used in this study CICsASS can be publicly downloaded.

We examine the constraints on soft X-ray emissions from the reionization era. It has generally been assumed that the Universe was reionized by ultraviolet photons from massive stars. However, it has been argued that X-ray photons associated with the death of these stars would have contributed ~10% to the total ionizations via several channels. The parameter space for a significant component of cosmological reionization to be sourced by X-rays is limited by a few observations. We revisit the unresolved soft X-ray background constraint and show that it significantly limits the contribution to reionization from several potential sources: X-rays from X-ray binaries, from Compton scattering off supernovae-accelerated electrons, and from the annihilation of dark matter particles. We discuss the additional limits on high-redshift X-ray production from (1) z~3 measurements of metal absorption lines, (2) the consensus that helium reionization was ending at z~3, and (3) measurements of the intergalactic medium's thermal history. We show that observations of z~3 metal lines allow little room for extra coeval X-ray emission from nonstandard sources. In addition, we show that the late reionization of helium makes it difficult to also ionize the hydrogen at z>6 with a single source population (such as quasars) and that it likely requires the spectrum of ionizing emissions to soften with increasing redshift. However, it is difficult to constrain an X-ray contribution to reionization from the intergalactic temperature history. We show that the gas would have been heated to a narrower range of temperatures than is typically assumed at reionization, 2-3 x10^4 K.

Telescopes aiming to measure 21cm emission from the Epoch of Reionization must toe a careful line, balancing the need for raw sensitivity against the stringent calibration requirements for removing bright foregrounds. It is unclear what the optimal design is for achieving both of these goals. Via a pedagogical derivation of an interferometer's response to the power spectrum of 21cm reionization fluctuations, we show that even under optimistic scenarios, first-generation arrays will yield low-SNR detections, and that different compact array configurations can substantially alter sensitivity. We explore the sensitivity gains of array configurations that yield high redundancy in the uv-plane -- configurations that have been largely ignored since the advent of self-calibration for high-dynamic-range imaging. We first introduce a mathematical framework to generate optimal minimum-redundancy configurations for imaging. We contrast the sensitivity of such configurations with high-redundancy configurations, finding that high-redundancy configurations can improve power-spectrum sensitivity by more than an order of magnitude. We explore how high-redundancy array configurations can be tuned to various angular scales, enabling array sensitivity to be directed away from regions of the uv-plane (such as the origin) where foregrounds are brighter and where instrumental systematics are more problematic. We demonstrate that a 132-antenna deployment of the Precision Array for Probing the Epoch of Reionization (PAPER) observing for 120 days in a high-redundancy configuration will, under ideal conditions, have the requisite sensitivity to detect the power spectrum of the 21cm signal from reionization at a 3\sigma level at k<0.25h Mpc^{-1} in a bin of \Delta ln k=1. We discuss the tradeoffs of low- versus high-redundancy configurations.

[ABRIDGED] Inhomogeneous reionization acts as a source of arcminute-scale anisotropies in the cosmic microwave background (CMB), the most important of which is the kinetic Sunyaev-Zel'dovich (kSZ) effect. Observational efforts with the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT) are poised to detect this signal for the first time. Indeed, recent SPT measurements place a bound on the dimensionless kSZ power spectrum around a multipole of l=3000 of P_tot < 2.8 (6) micro K^2 at 95% C.L., by ignoring (allowing) correlations between the thermal Sunyaev-Zel'dovich (tSZ) effect and the cosmic infrared background (CIB). To interpret these and upcoming observations, we compute the kSZ signal from a suite of ~ 100 reionization models using the publicly available code 21cmFAST. Our physically motivated reionization models are parameterized by the ionizing efficiency of high-redshift galaxies, the minimum virial temperature of halos capable of hosting stars, and the ionizing photon mean free path. We predict the contribution of patchy reionization to be P_patchy = 1.5-3.5 micro K^2. Therefore, even when conservatively adopting a low estimate of the post-reionization signal, P_OV ~ 2 micro K^2, none of our models are consistent with the aggressive 2sigma SPT bound that does not include correlations. This implies that either: (i) the early stages of reionization occurred in a much more homogeneous manner than suggested by the stellar-driven scenarios we explore, such as would be the case if, e.g., very high energy X-rays or exotic particles contributed significantly; and/or (ii) that there is a significant correlation between the CIB and the tSZ. On the other hand, the conservative SPT bound is compatible with all of our models, and is on the boarder of constraining reionization.

We study the properties of self-shielding intergalactic absorption systems and their implications for the ionizing background. We find that cosmological simulations post-processed with detailed radiative transfer calculations generally are able to reproduce the observed abundance of Lyman-limit systems, and we highlight possible discrepancies between the observations and simulations. This comparison tests cosmological simulations at overdensities of ~100. Furthermore, we show that the properties of Lyman-limit systems in these simulations, in simple semi-analytic arguments, and as suggested by recent observations indicate that a small change in the ionizing emissivity of the sources would have resulted in a much larger change in the amplitude of the intergalactic HI-ionizing background (with this scaling strengthening with increasing redshift). This strong scaling could explain the rapid evolution in the Lyman-alpha forest transmission observed at z = 6. Our calculations agree with the suggestion of simpler models that the comoving ionizing emissivity was constant or even increasing from z = 3 to 6. Our calculations also provide a more rigorous estimate than in previous studies for the clumping factor of intergalactic gas after reionization, which we estimate was 2 - 3 at z = 6.

It appears inevitable that reionization processes would have produced large-scale temperature fluctuations in the intergalactic medium. Using toy temperature models and detailed heating histories from cosmological simulations of HeII reionization, we study the consequences of inhomogeneous heating for the Ly-alpha forest. The impact of temperature fluctuations in physically well-motivated models can be surprisingly subtle. In fact, we show that temperature fluctuations at the level predicted by our reionization simulations do not give rise to detectable signatures in the types of statistics that have been employed previously. However, because of the aliasing of small-scale density power to larger scale modes in the line-of-sight Ly-alpha forest power spectrum, earlier analyses were not sensitive to 3D modes with >~ 30 comoving Mpc wavelengths -- scales where temperature fluctuations are likely to be relatively largest. The ongoing Baryon Oscillation Spectroscopic Survey (BOSS) aims to measure the 3D power spectrum of the Ly-alpha forest, P_F, from a large sample of quasars in order to avoid this aliasing. We find that physically motivated temperature models can alter P_F at an order unity level at k <~ 0.1 comoving Mpc^{-1}, a magnitude that should be easily detectable with BOSS. Fluctuations in the intensity of the ultraviolet background can also alter P_F significantly. These signatures will make it possible for BOSS to study the thermal impact of HeII reionization at 2 < z < 3 and to constrain models for the sources of the ionizing background. Future spectroscopic surveys could extend this measurement to even higher redshifts, potentially detecting the thermal imprint of hydrogen reionization.

Inhomogeneous reionization gives rise to angular fluctuations in the Cosmic Microwave Background (CMB) optical depth tau(n) to the last scattering surface, correlating different spherical harmonic modes and imprinting characteristic non-Gaussianity on CMB maps. Recently the minimum variance quadratic estimator has been derived using this mode-coupling signal, and found that the optical depth fluctuations could be detected with (S/N)^2 ~ 100 in futuristic experiments like CMBPol. We first demonstrate that the non-Gaussian signal from gravitational lensing of CMB is the dominant source of contamination for reconstructing inhomogeneous reionization signals, even with 98% of its contribution removed by delensing. We then construct unbiased estimators that simultaneously reconstruct inhomogeneous reionization signals tau(n) and gravitational lensing potential phi(n). We apply our new unbiased estimators to future CMB experiment to assess the detectability of inhomogeneous reionization signals. With more physically motivated simulations of inhomogenuous reionizations that predict an order of magnitude smaller optical depth power spectrum than previous studies, we show that a CMBPol-like experiment could achieve a marginal detection of inhomogeneous reionization,(S/N)^2 ~ O(1) with this quadratic estimator to ~O(10) with the analogous maximum likelihood estimator.

Weak gravitational lensing has become a powerful probe of large-scale structure and cosmological parameters. Precision weak lensing measurements require an understanding of the intrinsic alignment of galaxy ellipticities, which can in turn inform models of galaxy formation. It is hypothesized that elliptical galaxies align with the background tidal field and that this alignment mechanism dominates the correlation between ellipticities on cosmological scales (in the absence of lensing). We use recent large-scale structure measurements from the Sloan Digital Sky Survey to test this picture with several statistics: (1) the correlation between ellipticity and galaxy overdensity, w_{g+}; (2) the intrinsic alignment auto-correlation functions; (3) the correlation functions of curl-free, E, and divergence-free, B, modes (the latter of which is zero in the linear tidal alignment theory); (4) the alignment correlation function, w_g(r_p,theta), a recently developed statistic that generalizes the galaxy correlation function to account for the angle between the galaxy separation vector and the principle axis of ellipticity. We show that recent measurements are largely consistent with the tidal alignment model and discuss dependence on galaxy luminosity. In addition, we show that at linear order the tidal alignment model predicts that the angular dependence of w_g(r_p,theta) is simply w_{g+}*cos(2*theta) and that this dependence is consistent with recent measurements. We also study how stochastic nonlinear contributions to galaxy ellipticity impact these statistics. We find that a significant fraction of the observed LRG ellipticity can be explained by alignment with the tidal field on scales >~10 h^-1 Mpc. These considerations are relevant to galaxy formation and evolution.

We report on the detection of strongly varying intergalactic HeII absorption in HST/COS spectra of two z~3 quasars. From our homogeneous analysis of the HeII absorption in these and three archival sightlines, we find a marked increase in the mean HeII effective optical depth from tau~1 at z~2.3 to tau>5 at z~3.2, but with a large scatter of 2< tau <5 at 2.7< z <3 on scales of ~10 proper Mpc. This scatter is primarily due to fluctuations in the HeII fraction and the HeII-ionizing background, rather than density variations that are probed by the co-eval HI forest. Semianalytic models of HeII absorption require a strong decrease in the HeII-ionizing background to explain the strong increase of the absorption at z>2.7, probably indicating HeII reionization was incomplete at z>2.7. Likewise, recent three-dimensional numerical simulations of HeII reionization qualitatively agree with the observed trend only if HeII reionization completes at z=2.7 or even below, as suggested by a large tau>3 in two of our five sightlines at z<2.8. By doubling the sample size at 2.7< z <3, our newly discovered HeII sightlines for the first time probe the diversity of the second epoch of reionization when helium became fully ionized.

Analyses have found a "haze" of anomalous microwave emission surrounding the Galactic Center in the WMAP sky maps. A recent study using Fermi data detected a similar haze in the gamma-ray. Several studies have modeled these hazes as radiation from the leptonic byproducts of dark matter annihilations, and arguably no convincing astrophysical alternative has been suggested. We discuss the characteristics of astrophysical cosmic ray sources that could potentially explain this microwave and gamma-ray emission. The most promising astrophysical scenarios involve cosmic ray sources that are clustered such that many fall within ~1 kpc of the Galactic Center. For example, we show that several hundred Galactic Center supernovae in the last million years plus a diffusion-hardened electron spectrum may be consistent with present constraints on this emission. Alternatively, it could be due to a burst of activity probably associated with Sagittarius A* occurring ~1 Myr ago and producing >10^51 erg in cosmic ray electrons. Different models predict different trends for the spectral index of the microwave and gamma-ray spectrum as a function of angle from the Galactic Center that should be robust to cosmic ray propagation uncertainties. In particular, if the haze is from dark matter annihilations, it should have a very hard microwave and gamma-ray spectrum for which the spectral shape does not change significantly with angle, which we argue would be difficult to achieve with any astrophysical mechanism. Observations with the Planck and Fermi satellites can distinguish between viable haze models using these signatures.

It has been several years since the first detection of Gunn-Peterson troughs in the z > 6 Ly-alpha forest and since the first measurement of the Thomson scattering optical depth through reionization from the cosmic microwave background (CMB). Present day CMB measurements provide a significant constraint on the mean redshift of reionization, and the Ly-alpha forest provides a lower bound on the redshift at which reionization ended. However, no observation has provided definitive information on the duration and morphology of this process. This article is intended as a short review on the most promising observational methods that aim to detect and study this cosmic phase transition, focusing on CMB anisotropies, gamma ray burst afterglows, Ly-alpha emitting galaxies, and redshifted 21cm emission.

We discuss the potential of using the HeI 584 A forest to detect and study HeII reionization. Significant 584 A absorption is expected from intergalactic HeII regions, whereas there should be no detectable absorption from low density gas in HeIII regions. Unlike HeII Ly-alpha absorption (the subject of much recent study), the difficulty with using this transition to study HeII reionization is not saturation but rather that the absorption is weak. The Gunn-Peterson optical depth for this transition is tau ~ 0.1 x_{HeII} Delta^2 [(1+z)/5]^{9/2}, where x_{HeII} is the fraction of helium in HeII and Delta is the density in units of the cosmic mean. In addition, HeI 584 A absorption is contaminated by lower redshift HI Ly-alpha absorption with a comparable flux decrement. We estimate the requirements for a definitive detection of redshifted HeI absorption from low density gas (Delta ~ 1), which would indicate that HeII reionization was occurring. We find that this objective can be accomplished (using coeval HI Ly-alpha absorption to mask dense regions and in cross correlation) with a spectral resolution of 10^4 and a signal-to-noise ratio per resolution element of ~ 10. Such specifications may be achievable on a few known z ~ 3.5 quasar sightlines with the Cosmic Origins Spectrograph on the Hubble Space Telescope. We also discuss how HeI absorption can be used to measure the hardness of the ionizing background above 13.6 eV.

We compare the predictions of four different algorithms for the distribution of ionized gas during the Epoch of Reionization. These algorithms are all used to run a 100 Mpc/h simulation of reionization with the same initial conditions. Two of the algorithms are state-of-the-art ray-tracing radiative transfer codes that use disparate methods to calculate the ionization history. The other two algorithms are fast but more approximate schemes based on iterative application of a smoothing filter to the underlying source and density fields. We compare these algorithms' resulting ionization and 21 cm fields using several different statistical measures. The two radiative transfer schemes are in excellent agreement with each other (with the cross-correlation coefficient of the ionization fields >0.8 for k < 10 h/Mpc and in good agreement with the analytic schemes (>0.6 for k < 1 h/Mpc). When used to predict the 21cm power spectrum at different times during reionization, all ionization algorithms agree with one another at the 10s of percent level. This agreement suggests that the different approximations involved in the ray tracing algorithms are sensible and that semi-numerical schemes provide a numerically-inexpensive, yet fairly accurate, description of the reionization process.

We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the intergalactic medium, to detail the growth of HII regions around high redshift galaxies. Our simulation tracks 1024^3 dark matter particles, in a box of co-moving side length 65.6 Mpc/h. This large volume allows us to accurately characterize the size distribution of HII regions throughout most of the reionization process. At the same time, our simulation resolves many of the small galaxies likely responsible for reionization. It confirms a picture anticipated by analytic models: HII regions grow collectively around highly-clustered sources, and have a well-defined characteristic size, which evolves from a sub-Mpc scale at the beginning of reionization to R > 10 Mpc towards the end. We present a detailed statistical description of our results, and compare them with a numerical scheme based on the analytic model by Furlanetto, Zaldarriaga, and Hernquist. We find that the analytic calculation reproduces the size distribution of HII regions and the 21 cm power spectrum of the radiative transfer simulation remarkably well. The ionization field from the simulation, however, has more small scale structure than the analytic calculation, owing to Poisson scatter in the simulated abundance of galaxies on small scales. We propose and validate a simple scheme to incorporate this scatter into our calculations. Our results suggest that analytic calculations are sufficiently accurate to aid in predicting and interpreting the results of future 21 cm surveys. In particular, our fast numerical scheme is useful for forecasting constraints from future 21 cm surveys, and in constructing mock surveys to test data analysis procedures.

Motivated by recent interest in redshifted 21 cm emission of intergalactic hydrogen, we investigate the 8.7 GHz ^2S_{1/2} F=0-1 hyperfine transition of ^3He^+. While the primordial abundance of 3He relative to hydrogen is 10^-5, the hyperfine spontaneous decay rate is 680 times larger. Furthermore, the antenna temperature is much lower at the frequencies relevant for the ^3He^+ transition compared to that of z>6 21 cm emission. We find that the spin temperature of this 8.7 GHz line in the intergalactic medium is approximately the cosmic microwave background temperature, such that this transition is best observed in absorption against high-redshift, radio-bright quasars. We show that intergalactic 8.7 GHz absorption is a promising, unsaturated observable of the ionization history of intergalactic helium (for which HeII->HeIII reionization is believed to be complete at z~3) and of the primordial 3He abundance. Instruments must reach ~1 \mu Jy RMS noise in bands of 1 MHz on a 1 Jy source to directly resolve this absorption. However, in combination with HI Ly\alpha forest measurements, an instrument can statistically detect this absorption from z > 3 with 30 \mu Jy RMS noise in 0.1 MHz spectral bands over 100 MHz, which may be within the reach of present instruments.

The amplitude of fluctuations in the Ly-a forest on small spatial scales is sensitive to the temperature of the IGM and its spatial fluctuations. The temperature of the IGM and its spatial variations contain important information about hydrogen and helium reionization. We present a new measurement of the small-scale structure in the Ly-a forest from 40 high resolution, high signal-to-noise, VLT spectra at z=2.2-4.2. We convolve each Ly-a forest spectrum with a suitably chosen wavelet filter, which allows us to extract the amount of small-scale structure in the forest as a function of position across each spectrum. We compare these measurements with high resolution hydrodynamic simulations of the Ly-a forest which track more than 2 billion particles. This comparison suggests that the IGM temperature close to the cosmic mean density (T_0) peaks near z=3.4, at which point it is greater than 20,000 K at 2-sigma confidence. The temperature at lower redshift is consistent with the fall-off expected from adiabatic cooling ($T_0 \propto (1+z)^2$), after the peak temperature is reached near z=3.4. At z=4.2 our results favor a temperature of T_0 = 15-20,000 K. However, owing mostly to uncertainties in the mean transmitted flux at this redshift, a cooler IGM model with T_0 = 10,000 K is only disfavored at the 2-sigma level here, although such cool IGM models are strongly discrepant with the z ~ 3-3.4 measurement. We do not detect large spatial fluctuations in the IGM temperature at any redshift covered by our data set. The simplest interpretation of our measurements is that HeII reionization completes sometime near z ~ 3.4, although statistical uncertainties are still large [Abridged].

Many experts believe that the z ~ 3 HeII Ly-alpha forest will suffer from the same saturation issues as the z ~ 6 HI Ly-alpha forest and, therefore, will not be a sensitive probe of HeII reionization. However, there are several factors that make HeII Ly-alpha absorption more sensitive than HI Ly-alpha. We show that observations of HeII Ly-alpha and Ly-beta Gunn-Peterson troughs can provide a relatively model-independent constraint on the volume-averaged HeII fraction of x_HeII >~ 0.1. This bound derives from first using the most underdense regions in the HeII forest to constrain the local HeII fraction and, then, assuming photoionization equilibrium with the maximum allowed photoionization rate to calculate the ionization state of nearby gas. It is possible to evade this constraint by a factor of ~2, but only if the HeII were reionized recently. We argue that HeII Ly-alpha Gunn-Peterson troughs observed in the spectra of Q0302-003 and HE2347-4342 signify the presence of >~ 10 comoving Mpc patches in which x_HeII > 0.03. This is a factor of 20 improvement over previous constraints from these spectra and 100 times stronger than the tightest constraint on the HI volume-filling fraction from the z>6 HI Lyman forest.

The first structures in the Universe formed at z>7, at higher redshift than all currently known galaxies. Since GRBs are brighter than other cosmological sources at high redshift and exhibit simple power-law afterglow spectra that is ideal for absorption studies, they serve as powerful tools for studying the early universe. New facilities planned for the coming decade will be able to obtain a large sample of high-redshift GRBs. Such a sample would constrain the nature of the first stars, galaxies, and the reionization history of the Universe.

One of the next decade's most exciting prospects is to explore the cosmic "dark ages," during which the first stars in the Universe formed, with the 21 cm line of neutral hydrogen. At z>6, this light redshifts into the low-frequency radio regime and can be observed with new generations of wide-field arrays. These experiments have the potential to observe enormous volumes and to improve both cosmological and astrophysical constraints dramatically. In this white paper, we describe how the next decade will see the first steps toward unlocking this potential, allowing us to answer two key questions about the fundamental physics of our Universe and the intergalactic medium: Does the standard cosmological model describe the Universe during the "dark ages?" How did the intergalactic medium evolve during this important time, ending with the reionization of hydrogen?

The cosmic dark ages and the epoch of reionization, during which the first generations of stars and galaxies formed, are among the most compelling frontiers of extragalactic astrophysics and cosmology. Here we describe an exciting new probe of these eras: the 21 cm line of neutral hydrogen, which will allow us to map the neutral intergalactic medium throughout the era of first galaxy formation. In this white paper, we describe how these studies can answer two key questions about galaxy formation: What were the properties of high-redshift galaxies? How did they affect the Universe around them?

Observations of the intergalactic medium (IGM) suggest that quasars reionize HeII in the IGM at z ~ 3. We have run a set of 190 and 430 comoving Mpc simulations of HeII being reionized by quasars to develop an understanding of the nature of HeII reionization and its potential impact on observables. We find that HeII reionization heats regions in the IGM by as much as 25,000 K above the temperature that is expected otherwise, with the volume-averaged temperature increasing by ~ 12,000 K and with large temperature fluctuations on ~ 50 Mpc scales. Much of the heating occurs far from QSOs by hard photons. We find a temperature-density equation of state of gamma -1 ~ 0.3 during HeII reionization, but with a wide dispersion in this relation having sigma ~ 10^4 K. HeII reionization by the observed population of quasars cannot produce an inverted relation (gamma - 1 < 0). Our simulations are consistent with the observed evolution in the mean transmission of the HeII Ly-alpha forest. We argue that the heat input due to HeII reionization is unable to cause the observed depression at z = 3.2 in the HI Ly-alpha forest opacity as has been suggested. We investigate how uncertainties in the properties of QSOs and of HeII Lyman-limit systems influence our predictions.

There is growing interest in using 3-dimensional neutral hydrogen mapping with the redshifted 21 cm line as a cosmological probe, as it has been argued to have a greater long-term potential than the cosmic microwave background. However, its utility depends on many assumptions. To aid experimental planning and design, we quantify how the precision with which cosmological parameters can be measured depends on a broad range of assumptions. We cover assumptions related to modeling of the ionization power spectrum and associated nonlinearity, experimental specifications like array layout and noise, cosmological assumptions about reionization history and inter-galactic medium (IGM) evolution, and assumptions about astrophysical foregrounds. We derive simple analytic approximations for how various assumptions affect the results, and find that ionization power modeling is most important, followed by array layout (crudely, the more compact, the better). We also present an accurate yet robust method for measuring cosmological parameters in practice, separating the physics from the astrophysics by exploiting both gravitationally induced clustering anisotropy and the fact that the ionization power spectra are rather smooth functions that can be accurately fit by 7 phenomenological parameters. For example, a future square kilometer array optimized for 21 cm tomography could improve the sensitivity of the Planck CMB satellite to spatial curvature and neutrino masses by up to two orders of magnitude, to Delta-Omega_k ~ 0.0002 and Delta m_nu ~ 0.007 eV, and give a 4 sigma detection of the spectral index running predicted by the simpliest inflation models.

It is possible that the properties of HII regions during reionization depend sensitively on many poorly constrained quantities (the nature of the ionizing sources, the clumpiness of the gas in the IGM, the degree to which photo-ionizing feedback suppresses the abundance of low mass galaxies, etc.), making it extremely difficult to interpret upcoming observations of this epoch. We demonstrate that the actual situation is more encouraging, using a suite of radiative transfer simulations, post-processed on outputs from a 1024^3, 94 Mpc N-body simulation. Analytic prescriptions are used to incorporate small-scale structures that affect reionization, yet remain unresolved in the N-body simulation. We show that the morphology of the HII regions is most dependent on the global ionization fraction x_i. This is not to say that the bubble morphology is completely independent of all parameters besides x_i. The next most important dependence is that of the nature of the ionizing sources. The rarer the sources, the larger and more spherical the HII regions become. The typical bubble size can vary by as much as a factor of 4 at fixed x_i between different possible source prescriptions. The final relevant factor is the abundance of minihalos or of Lyman-limit systems. These systems suppress the largest bubbles from growing, and the magnitude of this suppression depends on the thermal history of the gas as well as the rate at which minihalos are photo-evaporated. We find that neither source suppression owing to photo-heating nor gas clumping significantly affect the large-scale structure of the HII regions. We discuss how observations of the 21cm line with MWA and LOFAR can constrain properties of reionization, and we study the effect patchy reionization has on the statistics of Lyman-alpha emitting galaxies. [abridged]

The cross-correlation between high redshift galaxies and 21 cm emission from the high redshift intergalactic medium (IGM) promises to be an excellent probe of the Epoch of Reionization (EoR). On large scales, the 21 cm and galaxy fields are anti-correlated during most of the reionization epoch. However, on scales smaller than the size of the H II regions around detectable galaxies, the two fields become roughly uncorrelated. Consequently, the 21 cm-galaxy cross power spectrum provides a tracer of bubble growth during reionization, with the signal turning over on progressively larger scales as reionization proceeds. The precise turnover scale depends on the minimum host mass of the detectable galaxies, and the galaxy selection technique. Measuring the turnover scale as a function of galaxy luminosity constrains the characteristic bubble size around galaxies of different luminosities. The cross spectrum becomes positive on small scales if ionizing photons fail to escape from low mass galaxies, and these galaxies are detectable longward of the hydrogen ionization edge, because in this case some identifiable galaxies lie outside of ionized regions. LOFAR can potentially measure the 21 cm-galaxy cross spectrum in conjunction with mild extensions to the existing Subaru survey for $z=6.6$ Lyman-alpha emitters, while the MWA is slightly less sensitive for detecting the cross spectrum. A futuristic galaxy survey covering a sizable fraction of the MWA field of view ($\sim 800$ deg$^2$) can probe the scale dependence of the cross spectrum, constraining the filling factor of H II regions at different redshifts during reionization, and providing other valuable constraints on reionization models.

We show that near-infrared observations of the red side of the Ly-alpha line from a single gamma ray burst (GRB) afterglow cannot be used to constrain the global neutral fraction of the intergalactic medium (IGM), x_H, at the GRB's redshift to better than ~0.3. Some GRB sight-lines will encounter more neutral hydrogen than others at fixed x_H owing to the patchiness of reionisation. GRBs during the epoch of reionization will often bear no discernible signature of a neutral IGM in their afterglow spectra. We discuss the constraints on x_H from the z = 6.3 burst, GRB050904, and quantify the probability of detecting a neutral IGM using future spectroscopic observations of high-redshift, near-infrared GRB afterglows. Assuming an observation with signal-to-noise similar to the Subaru FOCAS spectrum of GRB050904 and that the column density distribution of damped Ly-alpha absorbers is the same as measured at lower redshifts, a GRB from an epoch when x_H = 0.5 can be used to detect a partly neutral IGM at 98% confidence level 10% of the time (and, for an observation with three times the sensitivity, 30% of the time).

We forecast the sensitivity with which the Murchison Widefield Array (MWA) can measure the 21 cm power spectrum of cosmic hydrogen, using radiative transfer simulations to model reionization and the 21 cm signal. The MWA is sensitive to roughly a decade in scale (wavenumbers of k ~ 0.1 - 1 h Mpc^{-1}), with foreground contamination precluding measurements on larger scales, and thermal detector noise limiting the small scale sensitivity. This amounts primarily to constraints on two numbers: the amplitude and slope of the 21 cm power spectrum on the scales probed. We find, however, that the redshift evolution in these quantities can yield important information about reionization. Although the power spectrum differs substantially across plausible models, a generic prediction is that the amplitude of the 21 cm power spectrum on MWA scales peaks near the epoch when the intergalactic medium (IGM) is ~ 50% ionized. Moreover, the slope of the 21 cm power spectrum on MWA scales flattens as the ionization fraction increases and the sizes of the HII regions grow. Considering detection sensitivity, we show that the optimal MWA antenna configuration for power spectrum measurements would pack all 500 antenna tiles as close as possible in a compact core. The MWA is sensitive enough in its optimal configuration to measure redshift evolution in the slope and amplitude of the 21 cm power spectrum. Detecting the characteristic redshift evolution of our models will confirm that observed 21 cm fluctuations originate from the IGM, and not from foregrounds, and provide an indirect constraint on the volume-filling factor of HII regions during reionization. After two years of observations under favorable conditions, the MWA can constrain the filling factor at an epoch when <x_i> ~ 0.5 to within roughly +/- 0.1 at 2-sigma.

We show that observations of high-redshift Ly-alpha emitters (LAEs) have the potential to provide definitive evidence for reionization in the near future. Using 200 Mpc radiative transfer simulations, we calculate the effect that patchy reionization has on the line profile, on the luminosity function, and, most interestingly, on the clustering of emitters for several realistic models of reionization. Reionization increases the measured clustering of emitters, and we show that this enhancement would be essentially impossible to attribute to anything other than reionization. Our results motivate looking for the signature of reionization in existing LAE data. We find that for stellar reionization scenarios the angular correlation function of the 58 LAEs in the Subaru Deep Field z = 6.6 photometric sample is more consistent with a fully ionized universe (mean volume ionized fraction x_i = 1) than a universe with x_i < 0.5 at >2-sigma confidence level. Measurements in the next year on Subaru will increase their z = 6.6 LAE sample by a factor of five and tighten these limits. If the clustering signature of reionization is detected in a LAE survey, a comparison with a Lyman-break or a H-alpha survey in the same field would confirm the reionization hypothesis. We discuss the optimal LAE survey specifications for detecting reionization, with reference to upcoming programs.

Lyman-alpha forest absorption spectra towards quasars at z ~ 6 show regions of enhanced transmission close to their source. Several authors have argued that the apparently small sizes of these regions indicate that quasar ionization fronts at z >~ 6 expand into a largely or partly neutral intergalactic medium (IGM). Assuming that the typical region in the IGM is reionized by z <= 6, as is suggested by Ly-a forest observations, we argue that at {\em least} 50% of the volume of the IGM was reionized before the highest redshift quasars turned on. Further, even if the IGM is as much as 50% neutral at quasar turn-on, the quasars are likely born into large galaxy-generated HII regions. The HII regions during reionization are themselves clustered, and using radiative transfer simulations, we find that long skewers through the IGM towards quasar progenitor halos pass entirely through ionized bubbles, even when the IGM is half neutral. These effects have been neglected in most previous analyses of quasar proximity zones, which assumed a spatially {\em uniform} neutral fraction. We model the subsequent ionization from a quasar, and construct mock Ly-a forest spectra. Our mock absorption spectra are more sensitive to the level of small-scale structure in the IGM than to the volume-averaged neutral fraction, and suggest that existing proximity-zone size measurements are compatible with a fully ionized IGM. However, we mention several improvements in our modeling that are necessary to make more definitive conclusions.

A number of radio interferometers are currently being planned or constructed to observe 21 cm emission from reionization. Not only will such measurements provide a detailed view of that epoch, but, since the 21 cm emission also traces the distribution of matter in the Universe, this signal can be used to constrain cosmological parameters at 6 < z < 20. The sensitivity of an interferometer to the cosmological information in the signal may depend on how precisely the angular dependence of the 21 cm 3-D power spectrum can be measured. Utilizing an analytic model for reionization, we quantify all the effects that break the spherical symmetry of the 3-D 21 cm power spectrum and produce physically motivated predictions for this power spectrum. We find that upcoming observatories will be sensitive to the 21 cm signal over a wide range of scales, from larger than 100 to as small as 1 comoving Mpc. We consider three methods to measure cosmological parameters from the signal: (1) direct fitting of the density power spectrum to the signal, (2) using only the velocity field fluctuations in the signal, (3) looking at the signal at large enough scales such that all fluctuations trace the density field. With the foremost method, the first generation of 21 cm observations should moderately improve existing constraints on cosmological parameters for certain low-redshift reionization scenarios, and a two year observation with the second generation interferometer MWA5000 can improve constraints on Omega_w, Omega_m h^2, Omega_b h^2, Omega_nu, n_s, and alpha_s. If the Universe is substantially ionized by z = 12 or if spin temperature fluctuations are important, we show that it will be difficult to place competitive constraints on cosmological parameters with any of the considered methods.

We consider the contribution of 3rd and 4th order terms to the power spectrum of 21 cm brightness temperature fluctuations during the epoch of reionization, which arise because the 21 cm brightness temperature involves a product of the hydrogenic neutral fraction and the gas density. The 3rd order terms vanish for Gaussian random fields, and have been previously neglected or ignored. We measure these higher order terms from radiative transfer simulations and estimate them using cosmological perturbation theory. In our simulated models, the higher order terms are significant: neglecting them leads to a >~100% error in 21 cm power spectrum predictions on scales of k >~ 1 Mpc^{-1} when the neutral fraction is <x_H>~0.5. The higher order terms have a simple physical interpretation. On small scales they are produced by gravitational mode coupling. Small scale structure grows more readily in large-scale overdense regions, but the same regions tend to be ionized and hence do not contribute to the 21 cm signal. This acts to suppress the influence of non-linear density fluctuations and the small-scale amplitude of the 21 cm power spectrum. These results modify earlier intuition that the 21 cm power spectrum simply traces the density power spectrum on scales smaller than that of a typical bubble, and imply that small scale measurements contain more information about the nature of the ionizing sources than previously believed. On large scales, higher order moments are not directly related to gravity. They are non-zero because over-dense regions tend to ionize first and are important in magnitude at late times owing to the large fluctuations in the neutral fraction. (Abridged)

One of the key observables of the reionization era is the distribution of neutral and ionized gas. Recently, Furlanetto, Zaldarriaga, & Hernquist developed a simple analytic model to describe the growth of HII regions during this era. Here, we examine some of the fundamental simplifying assumptions behind this model and generalise it in several important ways. The model predicts that the ionized regions attain a well-defined characteristic size R_c that ranges from R_c~1 Mpc in the early phases to R_c>10 Mpc in the late phases. We show that R_c is determined primarily by the bias of the galaxies driving reionization; hence measurements of this scale constrain a fundamental property of the first galaxies. The variance around R_c, on the other hand, is determined primarily by the underlying matter power spectrum. We then show that increasing the ionizing efficiency of massive galaxies shifts R_c to significantly larger scales and decreases the importance of recombinations. These differences can be observed with forthcoming redshifted 21 cm surveys (increasing the brightness temperature fluctuations by up to a factor of two on large scales) and with measurements of small-scale anisotropies in the cosmic microwave background. Finally, we show that stochastic fluctuations in the galaxy population only broaden the bubble size distribution significantly if massive galaxies are responsible for most of the ionizing photons. We argue that the key results of this model are robust to many of our uncertainties about the reionization process.

We investigate the impact of spatial variations in the ionized fraction during reionization on temperature anisotropies in the CMB. We combine simulations of large scale structure to describe the underlying density field with an analytic model based on extended Press-Schechter theory to track the reionization process. We find that the power spectrum of the induced CMB anisotropies depends sensitively on the character of the reionization epoch. Models that differ in the extent of the "patchy phase" could be distinguished by future experiments such as the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT). In our models, the patchy signal peaks at $l\simeq 2000$, where it can be four times larger than the kinetic Sunyaev-Zel'dovich (kSZ)/Ostriker-Vishniac (OV) signal ($\Delta T_{\rm tot} \simeq 2.6 \mu K$). On scales beyond $l \simeq 4000$ the total Doppler signal is dominated by kSZ/OV, but the patchy signal can contribute up to 30% to the power spectrum. The effect of patchy reionization is largest on scales where the primordial CMB anisotropies dominate. Ignoring this contribution could lead to significant biases in the determination of cosmological parameters derived from CMB temperature measurements. Improvements in the theoretical modeling of the reionization epoch will become increasingly important to interpret the results of upcoming experiments.

During the epoch of reionization, local variations in the ionized fraction (patchiness) imprint arcminute-scale temperature anisotropies in the CMB through the kinetic Sunyaev-Zel'dovich (kSZ) effect. We employ an improved version of an analytic model of reionization devised in Furlanetto et al.(2004) to calculate the kSZ anisotropy from patchy reionization. This model uses extended Press-Schechter theory to determine the distribution and evolution of HII bubbles and produces qualitatively similar reionization histories to those seen in recent numerical simulations. We find that the angular power spectrum of the kSZ anisotropies depends strongly on the size distribution of the HII bubbles and on the duration of reionization. In addition, we show that upcoming measurements of the kSZ effect should be able to distinguish between several popular reionization scenarios. In particular, the amplitude of the patchy power spectrum for reionization scenarios in which the IGM is significantly ionized by Population III stars (or by mini-quasars/decaying particles) can be larger (or smaller) by over a factor of 3 than the amplitude in more traditional reionization histories (with temperature anisotropies that range between 0.5 and 3 micro-Kelvin at l = 5000). We highlight the differences in the kSZ signal between many possible reionization morphologies and discuss the constraints that future observations of the kSZ will place on this epoch.