Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.

The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz.

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from bright astrophysical foregrounds. If not properly calibrated, each antenna element's complex gain can impart spectral structure on otherwise spectrally-smooth foregrounds, complicating that statistical separation. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as redundant-baseline calibration resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions---both in data and in simulations---and present strategies for mitigating that structure.

We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.

This is the final report submitted to NASA for a Probe-class concept study of the "Farside Array for Radio Science Investigations of the Dark ages and Exoplanets" (FARSIDE), a low radio frequency interferometric array on the farside of the Moon. The design study focused on the instrument, a deployment rover, the lander and base station, and delivered an architecture broadly consistent with the requirements for a Probe mission. This notional architecture consists of 128 dipole antennas deployed across a 10 km area by a rover, and tethered to a base station for central processing, power and data transmission to the Lunar Gateway, or an alternative relay satellite. FARSIDE would provide the capability to image the entire sky each minute in 1400 channels spanning frequencies from 150 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, auroral kilometric radiation, and plasma noise from the solar wind. This would enable near-continuous monitoring of the nearest stellar systems in the search for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres for the nearest candidate habitable exoplanets. Simultaneously, FARSIDE would be used to characterize similar activity in our own solar system, from the Sun to the outer planets. Through precision calibration via an orbiting beacon, and exquisite foreground characterization, FARSIDE would also measure the Dark Ages global 21-cm signal at redshifts from 50-100. It will also be a pathfinder for a larger 21-cm power spectrum instrument by carefully measuring the foreground with high dynamic range.

We present a study of internal reflection and cross coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

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

We review some of the fundamental physics necessary for computing the highly-redshifted spin-flip background. We first discuss the radiative transfer of the 21-cm line and define the crucial quantities of interest. We then review the processes that set the spin temperature of the transition, with a particular focus on Wouthuysen-Field coupling, which is likely to be the most important process during and after the Cosmic Dawn. Finally, we discuss processes that heat the intergalactic medium during the Cosmic Dawn, including the scattering of Lyman-alpha, cosmic microwave background, and X-ray photons.

We describe how the high-redshift 21-cm background can be used to improve both our understanding of the fundamental cosmological parameters of our Universe and exotic processes originating in the dark sector. The 21-cm background emerging during the cosmological Dark Ages, the era between hydrogen recombination and the formation of the first luminous sources (likely at z ~ 30), is difficult to measure but provides several powerful advantages for these purposes: in addition to the lack of astrophysical contamination, it will allow probes of very small scale structure over a very large volume. Additionally, the 21-cm background is sensitive to the thermal state of the intergalactic hydrogen and therefore probes any exotic processes (including, e.g., dark matter scattering or decay and primordial black holes) during that era. After astrophysical sources have formed, cosmological information can be separated from astrophysical effects on the 21-cm background through methods such as redshift space distortions, joint modeling, and by searching for indirect effects on the astrophysical sources themselves.

The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Therefore, it is essential to understand and limit these chromatic effects. This paper describes a method to simulate the frequency and time responses of an antenna, by simultaneously taking into account the analogue RF receiver, the transmission cable, and the mutual coupling caused by adjacent antennas. Applied to the Hydrogen Epoch of Reionization Array (HERA), this study reveals the presence of significant reflections at high delays caused by the 150-m cable which links the antenna to the back-end. Besides, it shows that waves can propagate from one dish to another one through large sections of the array because of mutual coupling. In this more realistic approach, the simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h \rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity.

In this white paper, we lay out a US roadmap for high-redshift 21 cm cosmology (30 < z < 6) in the 2020s. Beginning with the currently-funded HERA and MWA Phase II projects and advancing through the decade with a coordinated program of small-scale instrumentation, software, and analysis projects targeting technology development, this roadmap incorporates our current best understanding of the systematics confronting 21 cm cosmology into a plan for overcoming them, enabling next-generation, mid-scale 21 cm arrays to be proposed late in the decade. Submitted for consideration by the Astro2020 Decadal Survey Program Panel for Radio, Millimeter, and Submillimeter Observations from the Ground as a Medium-Sized Project.

Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted HI emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to $k_\parallel$) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the "wedge"), and the spectrally structured emission from HI in the EoR (the "EoR window"). Faraday rotation can inject spectral structure into otherwise smooth polarized foreground emission, which through instrument effects or miscalibration could possibly pollute the EoR window. Using data from the HERA 19-element commissioning array, we investigate the polarization response of this new instrument in the power spectrum domain. We perform a simple image-based calibration based on the unpolarized diffuse emission of the Global Sky Model, and show that it achieves qualitative redundancy between the nominally-redundant baselines of the array and reasonable amplitude accuracy. We construct power spectra of all fully polarized coherencies in all pseudo-Stokes parameters. We compare to simulations based on an unpolarized diffuse sky model and detailed electromagnetic simulations of the dish and feed, confirming that in Stokes I, the calibration does not add significant spectral structure beyond the expected level. Further, this calibration is stable over the 8 days of observations considered. Excess power is seen in the power spectra of the linear polarization Stokes parameters which is not easily attributable to leakage via the primary beam, and results from some combination of residual calibration errors and actual polarized emission. Stokes V is found to be highly discrepant from the expectation of zero power, strongly pointing to the need for more accurate polarized calibration.

The Cosmic Dawn and Reionization epochs remain a fundamental but challenging frontier of astrophysics and cosmology. We advocate a large-scale, multi-tracer approach to develop a comprehensive understanding of the physics that led to the formation and evolution of the first stars and galaxies. We highlight the line intensity mapping technique to trace the multi-phase reionization topology on large scales, and measure reionization history in detail. Besides 21cm, we advocate for Lya tomography mapping during the epoch of Wouthuysen-Field coupling as an additional probe of the cosmic dawn era.

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.

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

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

Galaxies form and evolve in the context of their local and large-scale environments. Their baryonic content that we observe with imaging and spectroscopy is intimately connected to the properties of their dark matter halos, and to their location in the "cosmic web" of large-scale structure. Very large spectroscopic surveys of the local universe (e.g., SDSS and GAMA) measure galaxy positions (location within large-scale structure), statistical clustering (a direct constraint on dark matter halo masses), and spectral features (measuring physical conditions of the gas and stars within the galaxies, as well as internal velocities). Deep surveys with the James Webb Space Telescope (JWST) will revolutionize spectroscopic measurements of redshifts and spectral properties for galaxies out to the epoch of reionization, but with numerical statistics and over cosmic volumes that are too small to map large-scale structure and to constrain halo properties via clustering. Here, we consider advances in understanding galaxy evolution that would be enabled by very large spectroscopic surveys at high redshifts: very large numbers of galaxies (outstanding statistics) over large co-moving volumes (large-scale structure on all scales) over broad redshift ranges (evolution over most of cosmic history). The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

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

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.

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 epoch of reionization (6 < z < 10) marks the period in our universe when the first large galaxies grew to fruition, and began to affect the universe around them. Massive stars, and potentially accreting supermassive black holes, filled the universe with ionizing radiation, burning off the haze of neutral gas that had filled the intergalactic medium (IGM) since recombination (z~1000). The evolution of this process constrains key properties of these earliest luminous sources, thus observationally constraining reionization is a key science goal for the next decade. The measurement of Lyman-alpha emission from photometrically-identified galaxies is a highly constraining probe of reionization, as a neutral IGM will resonantly scatter these photons, reducing detectability. While significant work has been done with 8-10m telescopes, these observations require extremely large telescopes (ELTs); the flux limits available from today's 10m class telescopes are sufficient for only the brightest known galaxies (m < 26). Ultra-deep surveys with the Giant Magellan Telescope (GMT) and Thirty Meter Telescope (TMT) will be capable of detecting Lyman-alpha emission from galaxies 2-3 magnitudes fainter than today's deepest surveys. Wide-field fiber spectroscopy on the GMT combined with narrow-field AO-assisted slit spectroscopy on the TMT will be able to probe the expected size of ionized bubbles throughout the epoch of reionization, following up degree scale deep imaging surveys with the Wide Field Infrared Space Telescope. These data will provide the first resolved Lyman-alpha-based maps of the ionized intergalactic medium throughout the epoch of reionization, constraining models of both the temporal and spatial evolution of this phase change.

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.

Radio Frequency Interference (RFI) is an ever-present limiting factor among radio telescopes even in the most remote observing locations. When looking to retain the maximum amount of sensitivity and reduce contamination for Epoch of Reionization studies, the identification and removal of RFI is especially important. In addition to improved RFI identification, we must also take into account computational efficiency of the RFI-Identification algorithm as radio interferometer arrays such as the Hydrogen Epoch of Reionization Array grow larger in number of receivers. To address this, we present a Deep Fully Convolutional Neural Network (DFCN) that is comprehensive in its use of interferometric data, where both amplitude and phase information are used jointly for identifying RFI. We train the network using simulated HERA visibilities containing mock RFI, yielding a known "ground truth" dataset for evaluating the accuracy of various RFI algorithms. Evaluation of the DFCN model is performed on observations from the 67 dish build-out, HERA-67, and achieves a data throughput of 1.6$\times 10^{5}$ HERA time-ordered 1024 channeled visibilities per hour per GPU. We determine that relative to an amplitude only network including visibility phase adds important adjacent time-frequency context which increases discrimination between RFI and Non-RFI. The inclusion of phase when predicting achieves a Recall of 0.81, Precision of 0.58, and $F_{2}$ score of 0.75 as applied to our HERA-67 observations.

The Dark Ages, probed by the redshifted 21-cm signal, is the ideal epoch for a new rigorous test of the standard LCDM cosmological model. Divergences from that model would indicate new physics, such as dark matter decay (heating) or baryonic cooling beyond that expected from adiabatic expansion of the Universe. In the early Universe, most of the baryonic matter was in the form of neutral hydrogen (HI), detectable via its ground state's spin-flip transition. A measurement of the redshifted 21-cm spectrum maps the history of the HI gas through the Dark Ages and Cosmic Dawn and up to the Epoch of Reionization (EoR). The Experiment to Detect the Global EoR Signature (EDGES) recently reported an absorption trough at 78 MHz (redshift z of 17), similar in frequency to expectations for Cosmic Dawn, but about 3 times deeper than was thought possible from standard cosmology and adiabatic cooling of HI. Interactions between baryons and slightly-charged dark matter particles with electron-like mass provide a potential explanation of this difference but other cooling mechanisms are also being investigated to explain these results. The Cosmic Dawn trough is affected by cosmology and the complex astrophysical history of the first luminous objects. Another trough is expected during the Dark Ages, prior to the formation of the first stars and thus determined entirely by cosmological phenomena (including dark matter). Observations on or in orbit above the Moon's farside can investigate this pristine epoch (15-40 MHz; z=100-35), which is inaccessible from Earth. A single cross-dipole antenna or compact array can measure the amplitude of the 21-cm spectrum to the level required to distinguish (at >5$\sigma$}) the standard cosmological model from that of additional cooling derived from current EDGES results. This observation constitutes a powerful, clean probe of exotic physics in the Dark Ages.

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

New massively redundant low frequency arrays allow for a novel investigation of closure relations in interferometry. We employ commissioning data from the Hydrogen Epoch of Reionization Array to investigate closure quantities in this densely packed grid array of 14m antennas operating at 100 MHz to 200 MHz. We investigate techniques that utilize closure phase spectra for redundant triads to estimate departures from redundancy for redundant baseline visibilities. We find a median absolute deviation from redundancy in closure phase across the observed frequency range of about 4.5deg. This value translates into a non-redundancy per visibility phase of about 2.6deg, using prototype electronics. The median absolute deviations from redundancy decrease with longer baselines. We show that closure phase spectra can be used to identify ill-behaved antennas in the array, independent of calibration. We investigate the temporal behavior of closure spectra. The Allan variance increases after a one minute stride time, due to passage of the sky through the primary beam of the transit telescope. However, the closure spectra repeat to well within the noise per measurement at corresponding local sidereal times (LST) from day to day. In future papers in this series we will develop the technique of using closure phase spectra in the search for the HI 21cm signal from cosmic reionization.

The observed scatter in intergalactic Lyman-$\alpha$ opacity at $z \lesssim 6$ requires large-scale fluctuations in the neutral fraction of the intergalactic medium (IGM) after the expected end of reionization. Post-reionization models that explain this scatter invoke fluctuations in either the ionizing ultraviolet background (UVB) or IGM temperature. These models make very different predictions, however, for the relationship between Lyman-$\alpha$ opacity and local density. Here we test these models using Lyman-$\alpha$ emitting galaxies (LAEs) to trace the density field surrounding the longest and most opaque known Lyman-$\alpha$ trough at $z < 6$. Using deep Subaru Hyper Suprime-Cam narrow-band imaging, we find a highly significant deficit of $z \simeq 5.7$ LAEs within 20 Mpc/$h$ of the trough. The results are consistent with a model in which the scatter in Lyman-$\alpha$ opacity near $z \sim 6$ is driven by large-scale UVB fluctuations, and disfavor a scenario in which the scatter is primarily driven by variations in IGM temperature. UVB fluctuations at this epoch present a boundary condition for reionization models, and may help shed light on the nature of the ionizing sources.

The Experiment to Detect the Global Epoch of Reionization Signature (EDGES) recently reported a strong 21-cm absorption signal relative to the cosmic microwave background at $z \sim 18$. While its anomalous amplitude may indicate new physics, in this work we focus on the timing of the signal, as it alone provides an important constraint on galaxy formation models. Whereas rest-frame ultraviolet luminosity functions (UVLFs) over a broad range of redshifts are well fit by simple models in which galaxy star formation histories track the assembly of dark matter halos, we find that these same models, with reasonable assumptions about X-ray production in star-forming galaxies, cannot generate a narrow absorption trough at $z \sim 18$. If verified, the EDGES signal therefore requires the fundamental inputs of galaxy formation models to evolve rapidly at $z \gtrsim 10$. Unless extremely faint sources residing in halos below the atomic cooling threshold are responsible for the EDGES signal, star formation in $\sim 10^8$-$10^{10} \ M_{\odot}$ halos must be more efficient than expected, implying that the faint-end of the UVLF at $M_{\mathrm{UV}} \lesssim -12$ must steepen at the highest redshifts. This steepening provides a concrete test for future galaxy surveys with the James Webb Space Telescope and ongoing efforts in lensed fields, and is required regardless of whether the amplitude of the EDGES signal is due to new cooling channels or a strong radio background in the early Universe. However, the radio background solution requires that galaxies at $z > 15$ emit 1-2 GHz photons with an efficiency $\sim 10^3$ times greater than local star-forming galaxies, posing a challenge for models of low-frequency photon production in the early Universe.

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.

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

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

Recent observations have discovered long (up to ~110 Mpc/h), opaque Gunn-Peterson troughs in the z ~ 5.5 Lyman-alpha forest, which are challenging to explain with conventional models of the post-reionization intergalactic medium. Here we demonstrate that observations of the galaxy populations in the vicinity of the deepest troughs can distinguish two competing models for these features: deep voids where the ionizing background is weak due to fluctuations in the mean free path of ionizing photons would show a deficit of galaxies, while residual temperature variations from extended, inhomogeneous reionization would show an overdensity of galaxies. We use large (~550 Mpc/h) semi-numerical simulations of these competing explanations to predict the galaxy populations in the largest of the known troughs at z ~ 5.7. We quantify the strong correlation of Lyman-alpha effective optical depth and galaxy surface density in both models and estimate the degree to which realistic surveys can measure such a correlation. While a spectroscopic galaxy survey is ideal, we also show that a relatively inexpensive narrowband survey of Lyman-alpha-emitting galaxies is ~90% likely to distinguish between the competing models.

The redshifted 21-cm monopole is expected to be a powerful probe of the epoch of the first stars and galaxies ($10<z<35$). The global 21-cm signal is sensitive to the thermal and ionization state of hydrogen gas and thus provides a tracer of sources of energetic photons -- primarily hot stars and accreting black holes -- which ionize and heat the high redshift intergalactic medium (IGM). This paper presents a strategy for observations of the global spectrum with a realizable instrument placed in a low altitude lunar orbit, performing night-time 40-120 MHz spectral observations, while on the farside to avoid terrestrial radio frequency interference, ionospheric corruption, and solar radio emissions. The frequency structure, uniformity over large scales, and unpolarized state of the redshifted 21-cm spectrum are distinct from the spectrally featureless, spatially-varying, and polarized emission from the bright foregrounds. This allows a clean separation between the primordial signal and foregrounds. For signal extraction, we model the foreground, instrument, and 21-cm spectrum with eigenmodes calculated via Singular Value Decomposition analyses. Using a Markov Chain Monte Carlo algorithm to explore the parameter space defined by the coefficients associated with these modes, we illustrate how the spectrum can be measured and how astrophysical parameters (e.g., IGM properties, first star characteristics) can be constrained in the presence of foregrounds using the Dark Ages Radio Explorer (DARE).

Large variations in the effective optical depth of the He II Ly$\alpha$ forest have been observed at $z\gtrsim2.7$, but the physical nature of these variations is uncertain: either the Universe is still undergoing the process of He II reionization, or the Universe is highly ionized but the He II-ionizing background fluctuates significantly on large scales. In an effort to build upon our understanding of the latter scenario, we present a novel model for the evolution of ionizing background fluctuations. Previous models have assumed the mean free path of ionizing photons to be spatially uniform, ignoring the dependence of that scale on the local ionization state of the intergalactic medium (IGM). This assumption is reasonable when the mean free path is large compared to the average distance between the primary sources of He II-ionizing photons, $\gtrsim L_\star$ quasars. However, when this is no longer the case, the background fluctuations become more severe, and an accurate description of the average propagation of ionizing photons through the IGM requires additionally accounting for the fluctuations in opacity. We demonstrate the importance of this effect by constructing 3D semi-analytic models of the helium ionizing background from $z=2.5$-$3.5$ that explicitly include a spatially varying mean free path of ionizing photons. The resulting distribution of effective optical depths at large scales in the He II Ly$\alpha$ forest is very similar to the latest observations with HST/COS at $2.5 \lesssim z \lesssim 3.5$.

The experimental efforts to detect the redshifted 21 cm signal from the Epoch of Reionization (EoR) are limited predominantly by the chromatic instrumental systematic effect. The delay spectrum methodology for 21 cm power spectrum measurements brought new attention to the critical impact of an antenna's chromaticity on the viability of making this measurement. This methodology established a straightforward relationship between time-domain response of an instrument and the power spectrum modes accessible to a 21 cm EoR experiment. We examine the performance of a prototype of the Hydrogen Epoch of Reionization Array (HERA) array element that is currently observing in Karoo desert, South Africa. We present a mathematical framework to derive the beam integrated frequency response of a HERA prototype element in reception from the return loss measurements between 100-200 MHz and determined the extent of additional foreground contamination in the delay space. The measurement reveals excess spectral structures in comparison to the simulation studies of the HERA element. Combined with the HERA data analysis pipeline that incorporates inverse covariance weighting in optimal quadratic estimation of power spectrum, we find that in spite of its departure from the simulated response, HERA prototype element satisfies the necessary criteria posed by the foreground attenuation limits and potentially can measure the power spectrum at spatial modes as low as $k_{\parallel} > 0.1h$~Mpc$^{-1}$. The work highlights a straightforward method for directly measuring an instrument response and assessing its impact on 21 cm EoR power spectrum measurements for future experiments that will use reflector-type antenna.

Motivated by recent progress in studies of the high-$z$ Universe, we build a new model for the global 21-cm signal that is explicitly calibrated to measurements of the galaxy luminosity function (LF) and further tuned to match the Thomson scattering optical depth of the cosmic microwave background, $\tau_e$. Assuming that the $z \lesssim 8$ galaxy population can be smoothly extrapolated to higher redshifts, the recent decline in best-fit values of $\tau_e$ and the inefficient heating induced by X-ray binaries (HMXBs; the presumptive sources of the X-ray background at high-$z$) imply that the entirety of cosmic reionization and reheating occurs at redshifts $z \lesssim 12$. In contrast to past global 21-cm models, whose $z \sim 20$ ($\nu \sim 70$ MHz) absorption features and strong $\sim 25$ mK emission features were driven largely by the assumption of efficient early star-formation and X-ray heating, our new fiducial model peaks in absorption at $\nu \sim 110$ MHz at a depth of $\sim -160$ mK and has a negligible emission component. As a result, a strong emission signal would provide convincing evidence that HMXBs are not the only drivers of cosmic reheating. Shallow absorption troughs should accompany strong heating scenarios, but could also be caused by a low escape fraction of Lyman-Werner photons. Generating signals with troughs at $\nu \lesssim 95$ MHz requires a floor in the star-formation efficiency in halos below $\sim 10^{9} M_{\odot}$, which is equivalent to steepening the faint-end of the galaxy LF. These findings demonstrate that the global 21-cm signal is a powerful complement to current and future galaxy surveys and efforts to better understand the interstellar medium in high-$z$ galaxies.

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

A key physical quantity during reionization is the size of HII regions. Previous studies found a characteristic bubble size which increases rapidly during reionization, with apparent agreement between simulations and analytic excursion set theory. Using four different methods, we critically examine this claim. In particular, we introduce the use of the watershed algorithm -- widely used for void finding in galaxy surveys -- which we show to be an unbiased method with the lowest dispersion and best performance on Monte-Carlo realizations of a known bubble size PDF. We find that a friends-of-friends algorithm declares most of the ionized volume to be occupied by a network of volume-filling regions connected by narrow tunnels. For methods tuned to detect the volume-filling regions, previous apparent agreement between simulations and theory is spurious, and due to a failure to correctly account for the window function of measurement schemes. The discrepancy is already obvious from visual inspection. Instead, HII regions in simulations are significantly larger (by factors of 10-1000 in volume) than analytic predictions. The size PDF is narrower, and evolves more slowly with time, than predicted. It becomes more sharply peaked as reionization progresses. These effects are likely caused by bubble mergers, which are inadequately modeled by analytic theory. Our results have important consequences for high-redshift 21cm observations, the mean free path of ionizing photons, and the visibility of Ly-alpha emitters, and point to a fundamental failure in our understanding of the characteristic scales of the reionization process.

The Hydrogen Epoch of Reionization Array (HERA) is a staged experiment to measure 21 cm emission from the primordial intergalactic medium (IGM) throughout cosmic reionization ($z=6-12$), and to explore earlier epochs of our Cosmic Dawn ($z\sim30$). During these epochs, early stars and black holes heated and ionized the IGM, introducing fluctuations in 21 cm emission. HERA is designed to characterize the evolution of the 21 cm power spectrum to constrain the timing and morphology of reionization, the properties of the first galaxies, the evolution of large-scale structure, and the early sources of heating. The full HERA instrument will be a 350-element interferometer in South Africa consisting of 14-m parabolic dishes observing from 50 to 250 MHz. Currently, 19 dishes have been deployed on site and the next 18 are under construction. HERA has been designated as an SKA Precursor instrument. In this paper, we summarize HERA's scientific context and provide forecasts for its key science results. After reviewing the current state of the art in foreground mitigation, we use the delay-spectrum technique to motivate high-level performance requirements for the HERA instrument. Next, we present the HERA instrument design, along with the subsystem specifications that ensure that HERA meets its performance requirements. Finally, we summarize the schedule and status of the project. We conclude by suggesting that, given the realities of foreground contamination, current-generation 21 cm instruments are approaching their sensitivity limits. HERA is designed to bring both the sensitivity and the precision to deliver its primary science on the basis of proven foreground filtering techniques, while developing new subtraction techniques to unlock new capabilities. The result will be a major step toward realizing the widely recognized scientific potential of 21 cm cosmology.

Reionization is thought to have occurred in the redshift range of $6 < z < 9$, which is now being probed by both deep galaxy surveys and CMB observations. Using halo abundance matching over the redshift range $5<z<8$ and assuming smooth, continuous gas accretion, we develop a model for the star formation efficiency $f_{\star}$ of dark matter halos at $z>6$ that matches the measured galaxy luminosity functions at these redshifts. We find that $f_{\star}$ peaks at $\sim 30\%$ at halo masses $M \sim 10^{11}$--$10^{12}$~M$_\odot$, in qualitative agreement with its behavior at lower redshifts. We then investigate the cosmic star formation histories and the corresponding models of reionization for a range of extrapolations to small halo masses. We use a variety of observations to further constrain the characteristics of the galaxy populations, including the escape fraction of UV photons. Our approach provides an empirically-calibrated, physically-motivated model for the properties of star-forming galaxies sourcing the epoch of reionization. In the case where star formation in low-mass halos is maximally efficient, an average escape fraction $\sim0.1$ can reproduce the optical depth reported by Planck, whereas inefficient star formation in these halos requires either about twice as many UV photons to escape, or an escape fraction that increases towards higher redshifts. Our models also predict how future observations with JWST can improve our understanding of these galaxy populations.

The reionization of intergalactic hydrogen has received intense theoretical scrutiny over the past two decades. Here, we approach the process formally as a percolation process and phase transition. Using semi-numeric simulations, we demonstrate that an infinitely-large ionized region abruptly appears at an ionized fraction of ~0.1 and quickly grows to encompass most of the ionized gas: by an ionized fraction of 0.3, nearly ninety percent of the ionized material is part of this region. Throughout most of reionization, nearly all of the intergalactic medium is divided into just two regions, one ionized and one neutral, and both infinite in extent. We also show that the discrete ionized regions that exist before and near this transition point follow a near-power law distribution in volume, with equal contributions to the total filling factor per logarithmic interval in size up to a sharp cutoff in volume. These qualities are generic to percolation processes, with the detailed behavior a result of long-range correlations in the underlying density field. These insights will be crucial to understanding the distribution of ionized and neutral gas during reionization and provide precise meaning to the intuitive description of reionization as an "overlap" process.

One of the most exciting probes of the early phases of structure formation and reionization is the spin-flip line of neutral hydrogen, with a rest wavelength of 21 cm. This chapter introduces the physics of this transition and the astrophysical parameters upon which it depends, including discussions of the radiation fields that permeate the intergalactic medium that fix the brightness of this transition. We describe the critical points in the evolution of the 21-cm background and focus on the sky-averaged brightness and the power spectrum as representative measurements. Finally, we include a discussion of observations and the challenges they face in the near future.

The ionizing background of cosmic hydrogen is an important probe of the sources and absorbers of ionizing radiation in the post-reionization universe. Previous studies show that the ionization rate should be very sensitive to changes in the source population: as the emissivity rises, absorbers shrink in size, increasing the ionizing mean free path and, hence, the ionizing background. By contrast, observations of the ionizing background find a very flat evolution from z~2-5, before falling precipitously at z~6. We resolve this puzzling discrepancy by pointing out that, at z~2-5, optically thick absorbers are associated with the same collapsed halos that host ionizing sources. Thus, an increasing abundance of galaxies is compensated for by a corresponding increase in the absorber population, which moderates the instability in the ionizing background. However, by z~5-6, gas outside of halos dominates the absorption, the coupling between sources and absorbers is lost, and the ionizing background evolves rapidly. Our halo based model reproduces observations of the ionizing background, its flatness and sudden decline, as well as the redshift evolution of the ionizing mean free path. Our work suggests that, through much of their history, both star formation and photoelectric opacity in the universe track halo growth.

Extremely large opaque troughs in the Lyman-alpha forest have been interpreted as a sign of an extended reionization process below z~6. Such features are impossible to reproduce with simple models of the intergalactic ionizing background that assume a uniform mean free path of ionizing photons. We build a self-consistent model of the ionizing background that includes fluctuations in the mean free path due to the varying strength of the ionizing background and large-scale density field. The dominant effect is the suppression of the ionizing background in large-scale voids due to "self-shielding" by an enhanced number of optically thick absorbers. Our model results in a distribution of 50 Mpc/h Lyman-alpha forest effective optical depths that significantly improves agreement with the observations at z~5.6. Extrapolation to z~5.4 and z~5.8 appears promising, but matching the mean background evolution requires evolution in the absorber population beyond the scope of the present model. We also demonstrate the need for extremely large volumes (>400 Mpc on a side) to accurately determine the incidence of rare large-scale features in the Lyman-alpha forest.

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.

We discuss new constraints on the epoch of cosmic reionization and test the assumption that most of the ionizing photons responsible arose from high redshift star-forming galaxies. Good progress has been made in charting the end of reionization through spectroscopic studies of z~6-8 QSOs, gamma-ray bursts and galaxies expected to host Lyman-alpha emission. However, the most stringent constraints on its duration have come from the integrated optical depth, tau, of Thomson scattering to the cosmic microwave background. Using the latest data on the abundance and luminosity distribution of distant galaxies from Hubble Space Telescope imaging, we simultaneously match the reduced value tau=0.066 +/- 0.012 recently reported by the Planck collaboration and the evolving neutrality of the intergalactic medium with a reionization history within 6 <~ z <~ 10, thereby reducing the requirement for a significant population of very high redshift (z>>10) galaxies. Our analysis strengthens the conclusion that star-forming galaxies dominated the reionization process and has important implications for upcoming 21cm experiments and searches for early galaxies with James Webb Space Telescope.

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 test the ability of equilibrium galactic disk and one-zone interstellar medium models to describe the physical and emission properties of quasar hosts, submillimeter galaxies, and Lyman-alpha emitters at z>~6. The size, line widths, star formation rates, black hole accretion rates, gas masses and temperatures, and the relationships between these properties are all well-described by our model, and we provide approximate fitting formulae for comparison with future observations. However, comparing our carbon line predictions to observations reveals differences between the ISM at low and high redshifts. Our underestimate of the [CII] line emission indicates either higher star formation efficiencies in high-redshift molecular clouds or less depletion of metals into dust at fixed metallicity. Further, our over-prediction of the CO(6-5)/CO(1-0) ratio suggests that molecular clouds in real high-redshift galaxies have a lower turbulent Mach number and more subthermal CO(6-5) emission than expected owing either to sizes smaller than the local Jeans mass or to a pressure support mechanism other than turbulence.

Interpretation of He II Ly{\alpha} absorption spectra after the epoch of He II reionization requires knowledge of the He II ionizing background. While past work has modelled the evolution of the average background, the standard cosmological radiative transfer technique assumes a uniform radiation field despite the discrete nature of the (rare) bright quasars that dominate the background. We implement a cosmological radiative transfer model that includes the most recent constraints on the ionizing spectra and luminosity function of quasars and the distribution of IGM absorbers. We also estimate, for the first time, the effects of fluctuations on the evolving continuum opacity in two ways: by incorporating the complete distribution of ionizing background amplitudes into the standard approach, and by explicitly treating the quasars as discrete -- but isolated -- sources. Our model results in a He II ionization rate that evolves steeply with redshift, increasing by a factor ~2 from z=3.0 to z=2.5. This causes rapid evolution in the mean He II Ly{\alpha} optical depth -- as recently observed -- without appealing to the reionization of He II. The observed behaviour could instead result from rapid evolution in the mean free path of ionizing photons as the helium in higher H I column density absorbers becomes fully ionized.

We explore molecular cloud properties and the physics of CO transition lines in z >~ 6 Lyman-break galaxies and predict their CO fluxes using an analytic formalism built from global models of star formation in high-redshift galaxies that minimizes our reliance on local observations. Our model includes a new approach to calculating the molecular gas fraction that takes the total gas density, the star formation rate, and the star formation efficiency in clouds as the principal inputs. This method agrees with chemical equilibrium calculations of the molecular fraction based on local chemistry if galaxies at z >~ 6 have metallicities of order a few percent of solar. Such low metallicities in turn imply that much of the carbon in these systems exists in ionized form rather than as CO. Moreover, we find that the higher-order CO transitions observable at high redshift with ALMA will typically be sub-thermally populated but that the details depend sensitively on the presence of turbulent clumps within molecular clouds. Ultimately, we expect current facilities will only be able to observe the CO signal from reionization epoch galaxies with great difficulty. We estimate that at least ~100 hours of integration time with ALMA will be required to detect the CO(6--5) transition in z = 6 systems with rest-frame UV magnitudes of -20.

Recent He II Lyman-alpha forest observations from 2.0 < z < 3.2 show large fluctuations in the optical depth at z > 2.7. These results point to a fluctuating He-ionizing background, which may be due to the end of helium reionization of this era. We present a fast, semi-numeric procedure to approximate detailed cosmological simulations. We compute the distribution of dark matter halos, ionization state of helium, and density field at z = 3 in broad agreement with recent simulations. Given our speed and flexibility, we investigate a range of ionizing source and active quasar prescriptions. Spanning a large area of parameter space, we find order-of-magnitude fluctuations in the He II ionization rate in the post-reionization regime. During reionization, the fluctuations are even stronger and develop a bimodal distribution, in contrast to semi-analytic models and the hydrogen equivalent. These distributions indicate a low-level ionizing background even at significant He II fractions.

We use the new ultra-deep, near-infrared imaging of the Hubble Ultra-Deep Field (HUDF) provided by our UDF12 HST WFC3/IR campaign to explore the rest-frame UV properties of galaxies at redshifts z > 6.5. We present the first unbiased measurement of the average UV power-law index, beta, for faint galaxies at z ~ 7, the first meaningful measurements of beta at z ~ 8, and tentative estimates for a new sample of galaxies at z ~ 9. Utilising galaxy selection in the new F140W imaging to minimize colour bias, and applying both colour and power-law estimators of beta, we find beta = -2.1 (+/-0.2) at z ~ 7 for galaxies with M_UV ~ -18. This means that the faintest galaxies uncovered at this epoch have, on average, UV colours no more extreme than those displayed by the bluest star-forming galaxies at low redshift. At z ~ 8 we find a similar value, beta = -1.9 (+/-0.3). At z ~ 9, we find beta = -1.8 (+/-0.6), essentially unchanged from z ~ 6 - 7 (albeit highly uncertain). Finally, we show that there is as yet no evidence for a significant intrinsic scatter in beta within our new, robust z ~ 7 galaxy sample. Our results are most easily explained by a population of steadily star-forming galaxies with either ~ solar metallicity and zero dust, or moderately sub-solar (~ 10-20%) metallicity with modest dust obscuration (A_V ~ 0.1-0.2). This latter interpretation is consistent with the predictions of a state-of-the-art galaxy-formation simulation, which also suggests that a significant population of very-low metallicity, dust-free galaxies with beta ~ -2.5 may not emerge until M_UV > -16, a regime likely to remain inaccessible until the James Webb Space Telescope.

We present a new determination of the UV galaxy luminosity function (LF) at redshift z ~ 7 and z ~ 8, and a first estimate at z ~ 9. An accurate determination of the form and evolution of the LF at high z is crucial for improving our knowledge of early galaxy evolution and cosmic reionization. Our analysis exploits fully the new, deepest WFC3/IR imaging from our HST UDF12 campaign, and includes a new, consistent analysis of all appropriate, shallower/wider-area HST data. Our new measurement of the evolving LF at z ~ 7-8 is based on a final catalogue of ~600 galaxies, and involves a step-wise maximum likelihood determination based on the redshift probability distribution for each object; this makes full use of the 11-band imaging now available in the HUDF, including the new UDF12 F140W data, and the deep Spitzer IRAC imaging. The final result is a determination of the z ~ 7 LF extending down to M_UV = -16.75, and the z ~ 8 LF down to M_UV = -17.00. Fitting a Schechter function, we find M* = -19.90 (+0.23/-0.28), log phi* = -2.96 (+0.18/-0.23), and a faint-end slope alpha=-1.90 (+0.14/-0.15) at z~7, and M* = -20.12 (+0.37/-0.48), log phi* = -3.35 (+0.28/-0.47), alpha=-2.02 (+0.22/-0.23) at z~8. These results strengthen suggestions that the evolution at z > 7 is more akin to `density evolution' than the apparent `luminosity evolution' seen at z ~ 5-7. We also provide the first meaningful information on the LF at z ~ 9, explore alternative extrapolations to higher z, and consider the implications for the evolution of UV luminosity density. Finally, we provide catalogues (including z_phot, M_UV and all photometry) for the 100 most robust z~6.5-11.9 galaxies in the HUDF used in this analysis. We discuss our results in the context of earlier work and the results of an independent analysis of the UDF12 data based on colour-colour selection (Schenker et al. 2013).

Understanding cosmic reionization requires the identification and characterization of early sources of hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deepest infrared images with the Wide Field Camera 3 aboard Hubble Space Telescope and, for the first time, systematically explored the galaxy population deep into the era when cosmic microwave background (CMB) data indicates reionization was underway. The UDF12 campaign thus provides the best constraints to date on the abundance, luminosity distribution, and spectral properties of early star-forming galaxies. We synthesize the new UDF12 results with the most recent constraints from CMB observations to infer redshift-dependent ultraviolet (UV) luminosity densities, reionization histories, and electron scattering optical depth evolution consistent with the available data. Under reasonable assumptions about the escape fraction of hydrogen ionizing photons and the intergalactic medium clumping factor, we find that to fully reionize the universe by redshift z~6 the population of star-forming galaxies at redshifts z~7-9 likely must extend in luminosity below the UDF12 limits to absolute UV magnitudes of M_UV\sim -13 or fainter. Moreover, low levels of star formation extending to redshifts z~15-25, as suggested by the normal UV colors of z\simeq7-8 galaxies and the smooth decline in abundance with redshift observed by UDF12 to z\simeq10, are additionally likely required to reproduce the optical depth to electron scattering inferred from CMB observations.

We present a catalog of high redshift star-forming galaxies selected to lie within the redshift range z ~ 7-8 using the Ultra Deep Field 2012 (UDF12), the deepest near-infrared (near-IR) exposures yet taken with the Hubble Space Telescope. As a result of the increased near-infrared exposure time compared to previous HST imaging in this field, we probe 0.65 (0.25) mag fainter in absolute UV magnitude, at z ~ 7 (8), which increases confidence in a measurement of the faint end slope of the galaxy luminosity function. Through a 0.7 mag deeper limit in the key F105W filter that encompasses or lies just longward of the Lyman break, we also achieve a much-refined color-color selection that balances high redshift completeness and a low expected contamination fraction. We improve the number of drop-out selected UDF sources to 47 at z ~ 7 and 27 at z ~ 8. Incorporating brighter archival and ground-based samples, we measure the z ~ 7 UV luminosity function to an absolute magnitude limit of M_UV = -17 and find a faint end Schechter slope of \alpha = -1.87+/- 0.18. Using a similar color-color selection at z ~ 8 that takes account of our newly-added imaging in the F140W filter, and incorporating archival data from the HIPPIES and BoRG campaigns, we provide a robust estimate of the faint end slope at z ~ 8, \alpha = -1.94 +/- 0.23. We briefly discuss our results in the context of earlier work and that derived using the same UDF12 data but with an independent photometric redshift technique (McLure et al 2012).

We analyze the redshift- and luminosity-dependent sizes of dropout galaxy candidates in the redshift range z~7-12 using deep images from the UDF12 campaign, data which offers two distinct advantages over that used in earlier work. Firstly, we utilize the increased S/N ratio offered by the UDF12 imaging to provide improved size measurements for known galaxies at z=6.5-8 in the HUDF. Specifically, we stack the new deep F140W image with the existing F125W data in order to provide improved measurements of the half-light radii of z-dropouts. Similarly we stack this image with the new deep UDF12 F160W image to obtain new size measurements for a sample of Y-dropouts. Secondly, because the UDF12 data have allowed the construction of the first robust galaxy sample in the HUDF at z>8, we have been able to extend the measurement of average galaxy size out to significantly higher redshifts. Restricting our size measurements to sources which are now detected at >15sigma, we confirm earlier indications that the average half-light radii of z~7-12 galaxies are extremely small, 0.3-0.4 kpc, comparable to the sizes of giant molecular associations in local star-forming galaxies. We also confirm that there is a clear trend of decreasing half-light radius with increasing redshift, and provide the first evidence that this trend continues beyond z~8. Modeling the evolution of the average half-light radius as a power-law (1+z)^s, we obtain a best-fit index of s=-1.28+/-0.13 over the redshift range z~4-12, mid-way between the physically expected evolution for baryons embedded in dark halos of constant mass (s=-1) and constant velocity (s=-1.5). A clear size-luminosity relation, such as that found at lower redshift, is also evident in both our z- and Y-dropout sample. This relation can be interpreted in terms of a constant surface density of star formation over a range in luminosity of 0.05-1.0L*_z=3.(abridged)

We present the 2012 Hubble Ultra Deep Field campaign (UDF12), a large 128-orbit Cycle 19 \HST\ program aimed at extending previous WFC3/IR observations of the UDF by quadrupling the exposure time in the F105W filter, imaging in an additional F140W filter, and extending the F160W exposure time by 50%. The principal scientific goal of this project is to determine whether galaxies reionized the universe; our observations are designed to provide a robust determination of the star formation density at $z$$\,\gtrsim\,$8, improve measurements of the ultraviolet continuum slope at $z$$\,\sim\,7\,-\,$8, facilitate the construction of new samples of $z$$\,\sim\,9\,-\,$10 candidates, and enable the detection of sources up to $z$$\,\sim\,$12. For this project we committed to combining these and other WFC3/IR imaging observations of the UDF area into a single homogeneous dataset, to provide the deepest near-infrared observations of the sky currently achievable. In this paper we present the observational overview of the project, motivated by its scientific goals, and describe the procedures used in reducing the data as well as the final products that are produced. We have used the most up up-to-date methods for calibrating and combining the images, in particular paying attention to correcting several instrumental effects. We release the full combined mosaics, comprising a single, unified set of mosaics of the UDF, providing the deepest near-infrared blank-field view of the universe obtained to date, reaching magnitudes as deep as AB$\,\sim\,$30 in the near-infrared, and yielding a legacy dataset on this field of lasting scientific value to the community.

We present the results of the deepest search to date for star-forming galaxies beyond a redshift z~8.5 utilizing a new sequence of near-infrared Wide Field Camera 3 images of the Hubble Ultra Deep Field. This `UDF12' campaign completed in September 2012 doubles the earlier exposures with WFC3/IR in this field and quadruples the exposure in the key F105W filter used to locate such distant galaxies. Combined with additional imaging in the F140W filter, the fidelity of high redshift candidates is greatly improved. Using spectral energy distribution fitting techniques on objects selected from a deep multi-band near-infrared stack we find 7 promising z>8.5 candidates. As none of the previously claimed UDF candidates with 8.5<z<10 is confirmed by our deeper multi-band imaging, our campaign has transformed the measured abundance of galaxies in this redshift range. Although we recover the candidate UDFj-39546284 (previously proposed at z=10.3), it is undetected in the newly added F140W image, implying it lies at z=11.9 or is an intense emission line galaxy at z~2.4. Although no physically-plausible model can explain the required line intensity given the lack of Lyman alpha or broad-band UV signal, without an infrared spectrum we cannot rule out an exotic interloper. Regardless, our robust z ~ 8.5 - 10 sample demonstrates a luminosity density that continues the smooth decline observed over 6 < z < 8. Such continuity has important implications for models of cosmic reionization and future searches for z>10 galaxies with JWST.

We develop a radiation pressure-balanced model for the interstellar medium of high-redshift galaxies that describes many facets of galaxy formation at z>~6, including star formation rates and distributions and gas accretion onto central black holes. We first show that the vertical gravitational force in the disk of such a model is dominated by the disk self-gravity supported by the radiation pressure of ionizing starlight on gas. Constraining our model to reproduce the UV luminosity function of Lyman-break galaxies (LBGs), we limit the available parameter-space to wind mass-loading factors 1--4 times the canonical value for momentum-driven winds. We then focus our study by exploring the effects of different angular momentum transport mechanisms in the galactic disk and find that accretion driven by gravitational torques, such as from linear spiral waves or non-linear orbit crossings, can build up black hole masses by z=6 consistent with the canonical M-sigma relation with a duty cycle of unity, while accretion mediated by a local viscosity such as in an alpha-disk results in negligible BH accretion. Both gravitational torque models produce X-ray emission from active galactic nuclei (AGN) in high-redshift LBGs in excess of the estimated contribution from high-mass X-ray binaries. Using a recent analysis of deep Chandra observations by Cowie et al., we can already begin to rule out the most extreme regions of our parameter-space: the inflow velocity of gas through the disk must either be less than one percent of the disk circular velocity or the X-ray luminosity of the AGN must be substantially obscured. Moderately deeper future observations or larger sample sizes will be able to probe the more reasonable range of angular momentum transport models and obscuring geometries.

We use a new method to model fluctuations of the Lyman-Werner (LW) and Lyman-alpha radiation backgrounds at high redshift. At these early epochs the backgrounds are symptoms of a universe newly lit with its first stars. LW photons (11.5-13.6 eV) are of particular interest because they dissociate molecular hydrogen, the primary coolant in the first minihalos. By using a variation of the halo model, we efficiently generate power spectra for any choice of radiation background. We find that the LW power spectrum typically traces the matter power spectrum at large scales but turns over at the scale corresponding to the effective `horizon' of LW photons (~100 comoving Mpc), unless the sources are extremely rare. The series of horizons that characterize the Lyman-alpha flux profile shape the fluctuations of that background in a similar fashion, though those imprints are washed out once one considers fluctuations in the brightness temperature of the 21-cm signal. The Lyman-alpha background strongly affects the redshifted 21-cm signal at just about the time the LW background begins to dissociate molecular hydrogen, so measuring that background's properties will reveal important information about the transition from early Population III stars to more normal stars. Around this time we find that fluctuations in the LW background are weak; the fractional standard deviation is less than ~0.5 on scales > 10 cMpc, only rising to be of order unity on scales < 1 cMpc. This should not lead to substantial spatial fluctuations in molecular hydrogen content, except at the earliest times. Even then, most halos form far from other sources, so the transition from star formation in low-mass to high-mass halos is rather homogeneous across the universe.

A concept for a new space-based cosmology mission called the Dark Ages Radio Explore (DARE) is presented in this paper. DARE's science objectives include (1) When did the first stars form? (2) When did the first accreting black holes form? (3) When did Reionization begin? (4) What surprises does the end of the Dark Ages hold (e.g., Dark Matter decay)? DARE will use the highly-redshifted hyperfine 21-cm transition from neutral hydrogen to track the formation of the first luminous objects by their impact on the intergalactic medium during the end of the Dark Ages and during Cosmic Dawn (redshifts z=11-35). It will measure the sky-averaged spin temperature of neutral hydrogen at the unexplored epoch 80-420 million years after the Big Bang, providing the first evidence of the earliest stars and galaxies to illuminate the cosmos and testing our models of galaxy formation. DARE's approach is to measure the expected spectral features in the sky-averaged, redshifted 21-cm signal over a radio bandpass of 40-120 MHz. DARE orbits the Moon for a mission lifetime of 3 years and takes data above the lunar farside, the only location in the inner solar system proven to be free of human-generated radio frequency interference and any significant ionosphere. The science instrument is composed of a three-element radiometer, including electrically-short, tapered, bi-conical dipole antennas, a receiver, and a digital spectrometer. The smooth frequency response of the antennas and the differential spectral calibration approach using a Markov Chain Monte Carlo technique will be applied to detect the weak cosmic 21-cm signal in the presence of the intense solar system and Galactic foreground emissions.

We quantify the prospects for using emission lines from rotational transitions of the CO molecule to perform an `intensity mapping' observation at high redshift during the Epoch of Reionization (EoR). The aim of CO intensity mapping is to observe the combined CO emission from many unresolved galaxies, to measure the spatial fluctuations in this emission, and use this as a tracer of large scale structure at very early times in the history of our Universe. This measurement would help determine the properties of molecular clouds -- the sites of star formation -- in the very galaxies that reionize the Universe. We further consider the possibility of cross-correlating CO intensity maps with future observations of the redshifted 21 cm line. The cross spectrum is less sensitive to foreground contamination than the auto power spectra, and can therefore help confirm the high redshift origin of each signal. Furthermore, the cross spectrum measurement would help extract key information about the EoR, especially regarding the size distribution of ionized regions. We discuss uncertainties in predicting the CO signal at high redshift, and discuss strategies for improving these predictions. Under favorable assumptions, and feasible specifications for a CO survey mapping the CO(2-1) and CO(1-0) lines, the power spectrum of CO emission fluctuations and its cross spectrum with future 21 cm measurements from the MWA are detectable at high significance.

We show that comparisons of HeII Lyman-alpha forest lines of sight to nearby quasar populations can strongly constrain the lifetimes and emission geometry of quasars. By comparing the HeII and HI Lyman-alpha forests along a particular line of sight, one can trace fluctuations in the hardness of the radiation field (which are driven by fluctuations in the HeII ionization rate). Because this high-energy background is highly variable - thanks to the rarity of the bright quasars that dominate it and the relatively short attenuation lengths of these photons - it is straightforward to associate features in the radiation field with their source quasars. Here we quantify how finite lifetimes and beamed emission geometries affect these expectations. Finite lifetimes induce a time delay that displaces the observed radiation peak relative to the quasar. For beamed emission, geometry dictates that sources invisible to the observer can still create a peak in the radiation field. We show that both these models produce substantial populations of "bare" peaks (without an associated quasar) for reasonable parameter values (lifetimes ~10^6-10^8 yr and beaming angles <90 degrees). A comparison to existing quasar surveys along two HeII Lyman-alpha forest lines of sight rules out isotropic emission and infinite lifetime at high confidence; they can be accommodated either by moderate beaming or lifetimes ~10^7-10^8 yr. We also show that the distribution of radial displacements between peaks and their quasars can unambiguously distinguish these two models, although larger statistical samples are needed.

Lyman-limit absorption systems can play many important roles during and after cosmological reionization. Unfortunately, due to the prohibitively large dynamic range required, it is impossible to self-consistently include these systems in cosmological simulations. Using fast and versatile semi-numeric simulations, we systematically explore the spatial distribution of absorption systems during and following reionization. We self-calibrate the resulting number of absorbers to the mean free path (mfp) of the ionizing ultraviolet background (UVB), and present results at a given mfp and neutral hydrogen fraction. We use a simple optical depth criterion to identify the locations of absorbers. Our approach is fairly robust to uncertainties such as missing subgrid structure. Unlike at lower redshifts where the UVB is relatively uniform, at higher redshifts the fluctuations in the UVB and the HII morphology of reionization can drive the large-scale distribution of absorption systems. Specifically, we find that absorbers are highly correlated with the density field on small scales, and then become anti-correlated with the UVB on large scales. After reionization, the large-scale power spectrum of the absorbers traces the UVB power spectrum, which can be predicted with a simple analytic extension of the halo-model. During reionization, absorbers tend to preferentially lie inside overdensities (i.e. filaments) of the recently-ionized intergalactic medium (IGM). Absorbers may also dominate the small-scale (k > 1/Mpc) 21-cm power during and after reionization. Conversely, they smooth the contrast on moderate scales. Once the HII regions grow to surpass the mfp, the absorbers add to the large-scale 21-cm power. Our results should prove useful in interpreting future observations of the reionization epoch.

We introduce a powerful semi-numeric modeling tool, 21cmFAST, designed to efficiently simulate the cosmological 21-cm signal. Our code generates 3D realizations of evolved density, ionization, peculiar velocity, and spin temperature fields, which it then combines to compute the 21-cm brightness temperature. Although the physical processes are treated with approximate methods, we compare our results to a state-of-the-art large-scale hydrodynamic simulation, and find good agreement on scales pertinent to the upcoming observations (>~ 1 Mpc). The power spectra from 21cmFAST agree with those generated from the numerical simulation to within 10s of percent, down to the Nyquist frequency. We show results from a 1 Gpc simulation which tracks the cosmic 21-cm signal down from z=250, highlighting the various interesting epochs. Depending on the desired resolution, 21cmFAST can compute a redshift realization on a single processor in just a few minutes. Our code is fast, efficient, customizable and publicly available, making it a useful tool for 21-cm parameter studies.

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.

We examine large-scale fluctuations in the HeII Lyman-alpha forest transmission during and after HeII reionization. We use a simple Monte Carlo model to distribute quasars throughout a large volume and compute the resulting radiation field along one-dimensional skewers. In agreement with previous studies, we find that the rarity of these sources induces order unity fluctuations in the mean optical depth after reionization, even when averaged over large segments (~10-100 Mpc across). We compare our models to existing data along five HeII Lyman-alpha forest lines of sight spanning z~2-3.2. The large cosmic variance contained in our model plausibly explains many of the observed fluctuations at z<2.5. But our models cannot accommodate the large fluctuations toward high optical depths on ~30 Mpc scales observed at z~2.7-2.9, and the measured optical depths (tau>4) at z>2.9 are difficult to explain with a smoothly-evolving mean radiation field. In order to better understand this data, we construct a toy model of HeII reionization, in which we assume that regions with the smallest radiation fields in a post-reionization Universe (or farthest from strong ionizing sources) are completely dark during reionization. The observed fluctuations fit much more comfortably into this model, and we therefore argue that, according to present data, HeII reionization does not complete until z<2.9.

We examine the fate of fast electrons (with energies E>10 eV) in a thermal gas of primordial composition. To follow their interactions with the background gas, we construct a Monte Carlo model that includes: (1) electron-electron scattering (which transforms the electron kinetic energy into heat), (2) collisional ionization of hydrogen and helium (which produces secondary electrons that themselves scatter through the medium), and (3) collisional excitation (which produces secondary photons, whose fates we also follow approximately). For the last process, we explicitly include all transitions to upper levels n<=4, together with a well-motivated extrapolation to higher levels. In all cases, we use recent calculated cross-sections at E<1 keV and the Bethe approximation to extrapolate to higher energies. We compute the fractions of energy deposited as heat, ionization (tracking HI and the helium species separately), and excitation (tracking HI Lyman-alpha separately) under a broad range of conditions appropriate to the intergalactic medium. The energy deposition fractions depend on both the background ionized fraction and the electron energy but are nearly independent of the background density. We find good agreement with some, but not all, previous calculations at high energies. Electronic tables of our results are available on request.

The radiation background above the ionization edge of HeII varies strongly during and after helium reionization, because the attenuation length of such photons is relatively short (<40 Mpc) and because the ionizing sources (quasars) are rare. Here we construct analytic and Monte Carlo models to examine these fluctuations, including, for the first time, those during the reionization era itself. In agreement with detailed numerical simulations, our analytic model for the post-reionization Universe predicts order-of-magnitude fluctuations in the HeII ionization rate. Observations of the hardness ratio between HeII/HI show even larger fluctuations, which may be due to more complicated radiative transfer effects. During reionization, the fluctuations are even stronger. In contrast to hydrogen reionization, our model predicts that regions with strong HeII Lyman-alpha forest transmission should be reasonably common even during the beginning stages of reionization, because of strong illumination from bright quasars. Partly due to this, the mean ionizing background does not evolve strongly during and after helium reionization; it is roughly proportional to the filling fraction of HeIII regions. On the other hand, regions full of HeII and also "fossil" ionized regions that contain no (or few) active sources appear as strong IGM absorbers. Their presence exaggerates the evolution of the hardness ratio, making it evolve faster than naively expected during the reionization era.

Recent observations suggest that helium became fully ionized around redshift z ~ 3. The HeII optical depth derived from the Lyman-alpha forest decreases substantially from this period to z ~ 2; moreover, it fluctuates strongly near z ~ 3 and then evolves smoothly at lower redshifts. From these opacities, we compute, using a semi-analytic model, the evolution of the mean photoionization rate and the attenuation length for helium over the redshift range 2.0 < z < 3.2. This model includes an inhomogeneous metagalactic radiation background, which is expected during and after helium reionization. We find that assuming a uniform background underestimates the required photoionization rate by up to a factor ~2. When averaged over the (few) available lines of sight, the effective optical depth exhibits a discontinuity near z = 2.8, but the measurement uncertainties are sizable. This feature translates into a jump in the photoionization rate and, provided the quasar emissivity evolves smoothly, in the effective attenuation length, perhaps signaling the helium reionization era. We then compute the evolution of the effective optical depth for a variety of simple helium reionization models, in which the measured quasar luminosity function and the attenuation length, as well as the evolving HeIII fraction, are inputs. A model with reionization ending around redshift z ~ 2.7 is most consistent with the data, although the constraints are not strong thanks to the sparseness of the data.

We study the spatial fluctuations in the hydrogen ionizing background in the epoch following reionization (z ~ 5--6). The rapid decrease with redshift in the photon mean free path (m.f.p.), combined with the clustering of increasingly rare ionizing sources, can result in a very inhomogenous ionizing background during this epoch. We systematically investigate the probability density functions (PDFs) and power spectra of ionizing flux, by varying several parameters such as the m.f.p., minimum halo mass capable of hosting stars, and halo duty cycle. In order to be versatile, we make use of analytic, semi-numeric and numeric approaches. Our models show that the ionizing background indeed has sizable fluctuations during this epoch sourced by the clustering of sources, with the PDFs being a factor of few wide at half of the maximum likelihood. The distributions also show marked asymmetries, with a high-value tail set by clustering on small scales, and a shorter low-value tail which is set by the mean free path. The power spectrum of the ionizing background is much more sensitive to source properties than the PDF and can be well-understood analytically with a framework similar to the halo model (usually used to describe dark matter clustering). Nevertheless, we find that Lya forest spectra are extremely insensitive to the details of the UVB, despite marked differences in the PDFs and power spectra of our various ionizing backgrounds. Assuming a uniform ionizing background only underestimates the value of the mean ionization rate inferred from the Lya forest by a few percent. Instead, analysis of the Lya forest is dominated by the uncertainties in the density field. Thus, our results justify the common assumption of a uniform ionizing background in Lya forest analysis even during this epoch.

Three independent observational studies have now detected a narrow (\Delta z ~ 0.5) dip centred at z=3.2 in the otherwise smooth redshift evolution of the Lya forest effective optical depth. This feature has previously been interpreted as an indirect signature of rapid photo-heating in the IGM during the epoch of HeII reionisation. We examine this interpretation using a semi-analytic model of inhomogeneous HeII reionisation and high resolution hydrodynamical simulations of the Lya forest. We instead find that a rapid (\Delta z ~ 0.2) boost to the IGM temperature (\Delta T ~ 10^4 K) beginning at z=3.4 produces a well understood and generic evolution in the Lya effective optical depth, where a sudden reduction in the opacity is followed by a gradual, monotonic recovery driven largely by adiabatic cooling in the low density IGM. This behaviour is inconsistent with the narrow feature in the observational data. If photo-heating during HeII reionisation is instead extended over several redshift units, as recent theoretical studies suggest, then the Lya opacity will evolve smoothly with redshift. We conclude that the sharp dip observed in the Lya forest effective optical depth is instead most likely due to a narrow peak in the hydrogen photo-ionisation rate around z=3.2, and suggest that it may arise from the modulation of either reprocessed radiation during HeII reionisation, or the opacity of Lyman limit systems.

We use a combination of analytic and numerical arguments to consider the impact of quasar photo-heating during HeII reionisation on the thermal evolution of the intergalactic medium (IGM). We demonstrate that rapid (\Delta z< 0.1-0.2), strong (\Delta T > 10^4 K) photo-heating is difficult to achieve across the entire IGM unless quasar spectra are significantly harder than implied by current observational constraints. Although filtering of intrinsic quasar radiation through dense regions in the IGM does increase the mean excess energy per HeII photo-ionisation, it also weakens the radiation intensity and lowers the photo-ionisation rate, preventing rapid heating over time intervals shorter than the local photo-ionisation timescale. Moreover, the hard photons responsible for the strongest heating are more likely to deposit their energy inside dense clumps. The abundance of such clumps is, however, uncertain and model-dependent, leading to a fairly large uncertainty in the photo-heating rates. Nevertheless, although some of the IGM may be exposed to a hardened and weakened ionising background for long periods, most of the IGM must instead be reionised by the more abundant, softer photons and with accordingly modest heating rates (\Delta T < 10^4 K). The repeated ionisation of fossil quasar HeIII regions does not increase the net heating because the recombination times in these regions typically exceed the IGM cooling times and the average time lag between successive rounds of quasar activity. Detailed line-of-sight radiative transfer simulations confirm these expectations and predict a rich thermal structure in the IGM during HeII reionisation. [Abridged]

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

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?

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

After quasars ionize intergalactic HeII at z~3, a large radiation field builds up above the HeII ionization edge. Unlike the background responsible for HI ionizations, this field should be highly variable, thanks to the scarcity of bright quasars and the relatively short attenuation lengths (~50 Mpc) of these high-energy photons. Recent observations of the HeII and HI Lyman-alpha forests show that this background does indeed vary strongly, with substantial fluctuations on scales as small as ~2 Mpc. Here we show that such spatial fluctuation scales are naturally expected in any model in which the sources are as rare as bright quasars, so long as the attenuation length is relatively small. The correlation length itself is comparable to the attenuation length (~10 Mpc) for the most plausible physical scenarios, but we find order-of-magnitude fluctuations on all scales smaller than ~6 Mpc. Moreover, aliasing along the one-dimensional skewers probed by the HeII and HI Lyman-alpha forests exaggerates these variations, so that order-of-magnitude fluctuations should be observed on all scales smaller than ~20 Mpc. Complex radiative transfer is therefore not required to explain the observed fluctuations, at least at the level of current data.

We use an analytic model to study how inhomogeneous hydrogen reionization affects the temperature distribution of the intergalactic medium (IGM). During this process, the residual energy of each ionizing photon is deposited in the IGM as heat, increasing its temperature to 20,000-30,000 K; subsequent expansion of the Universe then cools the gas. Because reionization most likely proceeds from high to low densities, underdense voids are ionized last, have less time to cool, and are (on average) warmer than mean-density gas immediately after reionization is complete (an "inverted" equation of state). From this initial configuration, the low-density gas cools quickly and eventually returns to a more normal equation of state. The rapidly evolving temperature introduces systematic uncertainties in measurements of the ionizing background at z~6. For example, late reionization implies rapid cooling, so that the ionizing background would have to evolve even more rapidly at z ~5-6 than typically claimed. This degeneracy is difficult to disentangle, because the Lyman-alpha forest probes only a narrow range in densities (over which the gas is nearly isothermal). However, higher Lyman-series transitions probe wider density ranges, sampling different effective temperatures, and offer a new way to measure the IGM equation of state that should work where nearly saturated absorption precludes other methods. This will help to separate evolution in temperature from that in the ionizing background. While more detailed study with hydrodynamic simulations is needed, we show that such measurements could potentially distinguish early and late reionization using only a handful of lines of sight.

Weak gravitational lensing provides a potentially powerful method for the detection of clusters. In addition to cluster candidates, a large number of objects with possibly no optical or X-ray component have been detected in shear-selected samples. We develop an analytic model to investigate the claim of Weinberg & Kamionkowski (2002) that unvirialised protoclusters account for a significant number of these so-called "dark" lenses. In our model, a protocluster consists of a small virialised region surrounded by in-falling matter. We find that, in order for a protocluster to simultaneously escape X-ray detection and create a detectable weak lensing signal, it must have a small virial mass (~10^{13} \Msun) and large total mass (~ 10^{15} \Msun), with a relatively flat density profile outside of the virial radius. Such objects would be characterized by rising tangential shear profiles well beyond the virial radius. We use a semi-analytic approach based on the excursion set formalism to estimate the abundance of lensing protoclusters with a low probability of X-ray detection. We find that they are extremely rare, accounting for less than 0.4 per cent of the total lenses in a survey with background galaxy density n = 30 arcmin^{-2} and an intrinsic ellipticity dispersion of 0.3. We conclude that lensing protoclusters with undetectable X-Ray luminosities are too rare to account for a significant number of dark lenses.

One of the most sought-after signatures of reionization is a rapid increase in the ionizing background (usually measured through the Lyman-alpha optical depth toward distant quasars). Conventional wisdom associates this with the "overlap" phase when ionized bubbles merge, allowing each source to affect a much larger volume. We argue that this picture fails to describe the transition to the post-overlap Universe, where Lyman-limit systems absorb ionizing photons over moderate lengthscales (20-100 Mpc). Using an analytic model, we compute the probability distribution of the amplitude of the ionizing background throughout reionization, including both discrete ionized bubbles and Lyman-limit systems (parameterized by an attenuation length). We show that overlap does not by itself cause a rapid increase in the ionizing background or a rapid decrease in the mean Lyman-alpha transmission toward distant quasars. More detailed semi-numeric models support these conclusions. We argue that rapid changes should instead be interpreted as evolution in the attenuation length itself, which may or may not be directly related to overlap.

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.

A variety of observations now indicate that intergalactic helium was fully ionized by z~3. The most recent measurements of the high-redshift quasar luminosity function imply that these sources had produced ~2.5 ionizing photons per helium atom by that time, consistent with a picture in which the known quasar population drives HeII reionization. Here we describe the distribution of ionized and neutral helium gas during this era. Because the sources were rare and bright (with the photon budget dominated by quasars with luminosities L>L_\star), random fluctuations in the quasar population determined the morphology of ionized gas when the global ionized fraction x_i was small, with the typical radius R_c of a HeIII bubble ~15-20 comoving Mpc. Only when x_i>0.5 did the large-scale clustering of the quasars drive the characteristic size of ionized regions above this value. Still later, when x_i>0.75, most ionizing photons were consumed by dense, recombining systems before they reached the edge of their source's ionized surroundings, halting the bubble growth when R_c~35-40 Mpc. These phases are qualitatively similar to those in hydrogen reionization, but the rarity of the sources makes the early stochastic phase much more important. Moreover, the well-known characteristics of the z=3 intergalactic medium allow a much more robust description of the late phase in which recombinations dominate.

The temperature of the intergalactic medium (IGM) is set by the competition between photoheating and adiabatic cooling, which are usually assumed to define a tight equation of state in which the temperature increases monotonically with density. We use a semi-analytic model, accurate at low to moderate IGM densities (<5 times the mean), to show that inhomogeneous reionization can substantially modify these expectations. Because reionization is driven by biased sources, dense pockets of the IGM are likely to be ionized first. As a result, voids initially remain cool while dense regions are heated substantially. However, near the end of reionization, dense regions have already cooled from their initially large temperature while voids have only just been heated. Thus, near the end of helium reionization the equation of state can invert itself, with the hottest gas inside voids. The degree to which this happens depends on the magnitude of the density bias during reionization: if rare, bright sources dominate reionization, so that each ionized region contains a typical volume of the IGM, the equation of state will remain monotonic. We also show that the distribution of temperatures at a fixed density has significant scatter and evolves rapidly throughout and even after reionization. Finally, we show that the observed temperature jump at z ~3.2 is consistent with the behavior at the end of helium reionization, although it requires a somewhat larger temperature increase than expected.

One of the most dramatic signatures of the reionization era may be the enormous ionized bubbles around luminous quasars (with radii reaching ~40 comoving Mpc), which may survive as "fossil'' ionized regions long after their source shuts off. Here we study how the inhomogeneous intergalactic medium (IGM) evolves inside such fossils. The average recombination rate declines rapidly with time, and the brief quasar episode significantly increases the mean free path inside the fossil bubbles. As a result, even a weak ionizing background generated by galaxies inside the fossil can maintain it in a relatively highly and uniformly ionized state. For example, galaxies that would ionize 20-30% of hydrogen in a random patch of the IGM can maintain 80-90% ionization inside the fossil, for a duration much longer than the average recombination time in the IGM. Quasar fossils at z<10 thus retain their identity for nearly a Hubble time, and will appear "gray,'' distinct from both the average IGM (which has a "swiss-cheese" ionization topology and a lower mean ionized fraction), and from fully-ionized bubbles around active quasars. More distant fossils, at z>10 have a weaker galaxy-generated ionizing background and a higher gas density, so they can attain a swiss-cheese topology similar to the rest of the IGM, but with a smaller contrast between the ionized bubbles and the partially neutral regions separating them. Analogous HeIII-fossils should exist around the epoch of HeII/HeIII reionization at z~3, although rapid recombinations inside the HeIII-fossils will be more common. Our model of inhomogeneous recombination also applies to "double reionization'' models and shows that a non-monotonic reionization history is even more unlikely than previously thought.

We investigate the potential of exploiting Lya Emitters (LAEs) to constrain the volume-weighted mean neutral hydrogen fraction of the IGM, x_H, at high redshifts (specifically z~9). We use "semi-numerical'' simulations to efficiently generate density, velocity, and halo fields at z=9 in a 250 Mpc box, resolving halos with masses M>2.2e8 solar masses. We construct ionization fields corresponding to various values of x_H. With these, we generate LAE luminosity functions and "counts-in-cell'' statistics. As in previous studies, we find that LAEs begin to disappear rapidly when x_H > 0.5. Constraining x_H(z=9) with luminosity functions is difficult due to the many uncertainties inherent in the host halo mass <--> Lya luminosity mapping. However, using a very conservative mapping, we show that the number densities derived using the six z~9 LAEs recently discovered by Stark et al. (2007) imply x_H < 0.7. On a more fundamental level, these LAE number densities, if genuine, require substantial star formation in halos with M < 10^9 solar masses, making them unique among the current sample of observed high-z objects. Furthermore, reionization increases the apparent clustering of the observed LAEs. We show that a ``counts-in-cell'' statistic is a powerful probe of this effect, especially in the early stages of reionization. Specifically, we show that a field of view (typical of upcoming IR instruments) containing LAEs has >10% higher probability of containing more than one LAE in a x_H>0.5 universe than a x_H=0 universe with the same overall number density. With this statistic, a fully ionized universe can be robustly distinguished from one with x_H > 0.5 using a survey containing only ~ 20--100 galaxies.

Voids are a dominant feature of the low-redshift galaxy distribution. Several recent surveys have found evidence for the existence of large-scale structure at high redshifts as well. We present analytic estimates of galaxy void sizes at redshifts z ~ 5 - 10 using the excursion set formalism. We find that recent narrow-band surveys at z ~ 5 - 6.5 should find voids with characteristic scales of roughly 20 comoving Mpc and maximum diameters approaching 40 Mpc. This is consistent with existing surveys, but a precise comparison is difficult because of the relatively small volumes probed so far. At z ~ 7 - 10, we expect characteristic void scales of ~ 14 - 20 comoving Mpc assuming that all galaxies within dark matter haloes more massive than 10^10 M_sun are observable. We find that these characteristic scales are similar to the sizes of empty regions resulting from purely random fluctuations in the galaxy counts. As a result, true large-scale structure will be difficult to observe at z ~ 7 - 10, unless galaxies in haloes with masses less than ~ 10^9 M_sun are visible. Galaxy surveys must be deep and only the largest voids will provide meaningful information. Our model provides a convenient picture for estimating the "worst-case" effects of cosmic variance on high-redshift galaxy surveys with limited volumes.

One well-known way to constrain the hydrogen neutral fraction, x_H, of the high-redshift intergalactic medium (IGM) is through the shape of the red damping wing of the Lya absorption line. We examine this method's effectiveness in light of recent models showing that the IGM neutral fraction is highly inhomogeneous on large scales during reionization. Using both analytic models and "semi-numeric" simulations, we show that the "picket-fence" absorption typical in reionization models introduces both scatter and a systematic bias to the measurement of x_H. In particular, we show that simple fits to the damping wing tend to overestimate the true neutral fraction in a partially ionized universe, with a fractional error of ~ 30% near the middle of reionization. This bias is generic to any inhomogeneous model. However, the bias is reduced and can even underestimate x_H if the observational sample only probes a subset of the entire halo population, such as quasars with large HII regions. We also find that the damping wing absorption profile is generally steeper than one would naively expect in a homogeneously ionized universe. The profile steepens and the sightline-to-sightline scatter increases as reionization progresses. Of course, the bias and scatter also depend on x_H and so can, at least in principle, be used to constrain it. Damping wing constraints must therefore be interpreted by comparison to theoretical models of inhomogeneous reionization.

Of the many probes of reionization, the 21 cm line and the cosmic microwave background (CMB) are among the most effective. We examine how the cross-correlation of the 21 cm brightness and the CMB Doppler fluctuations on large angular scales can be used to study this epoch. We employ a new model of the growth of large scale fluctuations of the ionized fraction as reionization proceeds. We take into account the peculiar velocity field of baryons and show that its effect on the cross correlation can be interpreted as a mixing of Fourier modes. We find that the cross-correlation signal is strongly peaked toward the end of reionization and that the sign of the correlation should be positive because of the inhomogeneity inherent to reionization. The signal peaks at degree scales (l~100) and comes almost entirely from large physical scales (k~0.01 Mpc). Since many of the foregrounds and noise that plague low frequency radio observations will not correlate with CMB measurements, the cross correlation might appear to provide a robust diagnostic of the cosmological origin of the 21 cm radiation around the epoch of reionization. Unfortunately, we show that these signals are actually only weakly correlated and that cosmic variance dominates the error budget of any attempted detection. We conclude that the detection of a cross-correlation peak at degree-size angular scales is unlikely even with ideal experiments.

We present a method to construct semi-numerical ``simulations'', which can efficiently generate realizations of halo distributions and ionization maps at high redshifts. Our procedure combines an excursion-set approach with first-order Lagrangian perturbation theory and operates directly on the linear density and velocity fields. As such, the achievable dynamic range with our algorithm surpasses the current practical limit of N-body codes by orders of magnitude. This is particularly significant in studies of reionization, where the dynamic range is the principal limiting factor. We test our halo-finding and HII bubble-finding algorithms independently against N-body simulations with radiative transfer and obtain excellent agreement. We compute the size distributions of ionized and neutral regions in our maps. We find even larger ionized bubbles than do purely analytic models at the same volume-weighted mean hydrogen neutral fraction. We also generate maps and power spectra of 21-cm brightness temperature fluctuations, which for the first time include corrections due to gas bulk velocities. We find that velocities widen the tails of the temperature distributions and increase small-scale power, though these effects quickly diminish as reionization progresses. We also include some preliminary results from a simulation run with the largest dynamic range to date: a 250 Mpc box that resolves halos with masses M >~ 2.2 x10^8 M_sun. We show that accurately modeling the late stages of reionization requires such large scales. The speed and dynamic range provided by our semi-numerical approach will be extremely useful in the modeling of early structure formation and reionization.

The spin temperature of neutral hydrogen, which determines the optical depth and brightness of the 21 cm line, is determined by the competition between radiative and collisional processes. Here we examine the role of proton-hydrogen collisions in setting the spin temperature. We use recent fully quantum mechanical calculations of the relevant cross sections, which allow us to present accurate results over the entire physically relevant temperature range 1-10,000 K. For kinetic temperatures T_K>100 K, the proton-hydrogen rate coefficient exceeds that for hydrogen-hydrogen collisions by about a factor of two. However, at low temperatures (T_K < 5 K) H-p collisions become several thousand times more efficient than H-H and even more important than H-e^- collisions.

We examine the effect of inhomogeneous reionization on the galaxy power spectrum and the consequences for probing dark energy. To model feedback during reionization, we apply an ansatz setting the galaxy overdensity proportional to the underlying ionization field. Thus, inhomogeneous reionization may leave an imprint in the galaxy power spectrum. We evolve this imprint to low redshift and use the Fisher-matrix formalism to assess the effect on parameter estimation. We show that a combination of low- (z=0.3) and high- (z=3) redshift galaxy surveys can constrain the size of cosmological HII regions during reionization. This imprint can also cause confusion when using baryon oscillations or other features of the galaxy power spectrum to probe the dark energy. We show that when bubbles are large, and hence detectable, our ability to constrain w can be degraded by up to 50%. When bubbles are small, the imprint has little or no effect on measuring dark-energy parameters.

Many models of early structure formation predict a period of heating immediately preceding reionization, when X-rays raise the gas temperature above that of the cosmic microwave background. These X-rays are often assumed to heat the intergalactic medium (IGM) uniformly, but in reality they will heat the gas more strongly closer to the sources. We develop a framework for calculating fluctuations in the 21 cm brightness temperature that originate from this spatial variation in the heating rate. High-redshift sources are highly clustered, leading to significant gas temperature fluctuations (with fractional variations ~40%, peaking on k~0.1 Mpc^{-1} scales). This induces a distinctive peak-trough structure in the angle-averaged 21 cm power spectrum, which may be accessible to the proposed Square Kilometre Array. This signal reaches the ~10 mK level, and is stronger than that induced by Lyman alpha flux fluctuations. As well as probing the thermal evolution of the IGM before reionization, this 21 cm signal contains information about the spectra of the first X-ray sources. Finally, we consider disentangling temperature, density and Lyman alpha flux fluctuations as functions of redshift.

We study the detectability of the cross-correlation between 21 cm emission from the intergalactic medium and the galaxy distribution during (and before) reionization. We show that first-generation 21 cm experiments, such as the Mileura Widefield Array (MWA), can measure the cross-correlation to a precision of several percent on scales k~0.1/Mpc if combined with a deep galaxy survey detecting all galaxies with m>10^10 Msol over the entire ~800 square degree field of view of the MWA. LOFAR can attain even better limits with galaxy surveys covering its ~50 square degree field of view. The errors on the cross-power spectrum scale with the square root of the overlap volume, so even reasonably modest surveys of several square degrees should yield a positive detection with either instrument. In addition to the obvious scientific value, the cross-correlation has four key advantages over the 21 cm signal alone: (1) its signal-to-noise exceeds that of the 21 cm power spectrum by a factor of several, allowing it to probe smaller spatial scales and perhaps to detect inhomogeneous reionization more efficiently; (2) it allows a cleaner division of the redshift-space distortions (although only if the galaxy redshifts are known precisely); (3) by correlating with the high-redshift galaxy population, the cosmological nature of the 21 cm fluctuations can be determined unambiguously; and (4) the required level of foreground cleaning for the 21 cm signal is vastly reduced.

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.

Observations of the high-redshift Universe with the 21 cm hyperfine line of neutral hydrogen promise to open an entirely new window onto the early phases of cosmic structure formation. Here we review the physics of the 21 cm transition, focusing on processes relevant at high redshifts, and describe the insights to be gained from such observations. These include measuring the matter power spectrum at z~50, observing the formation of the cosmic web and the first luminous sources, and mapping the reionization of the intergalactic medium. The epoch of reionization is of particular interest, because large HII regions will seed substantial fluctuations in the 21 cm background. We also discuss the experimental challenges involved in detecting this signal, with an emphasis on the Galactic and extragalactic foregrounds. These increase rapidly toward low frequencies and are especially severe for the highest redshift applications. Assuming that these difficulties can be overcome, the redshifted 21 cm line will offer unique insight into the high-redshift Universe, complementing other probes but providing the only direct, three-dimensional view of structure formation from z~200 to z~6.

The radiation background produced by the 21 cm spin-flip transition of neutral hydrogen at high redshifts can be a pristine probe of fundamental physics and cosmology. At z~30-300, the intergalactic medium (IGM) is visible in 21 cm absorption against the cosmic microwave background (CMB), with a strength that depends on the thermal (and ionization) history of the IGM. Here we examine the constraints this background can place on dark matter decay and annihilation, which could heat and ionize the IGM through the production of high-energy particles. Using a simple model for dark matter decay, we show that, if the decay energy is immediately injected into the IGM, the 21 cm background can detect energy injection rates >10^{-24} eV cm^{-3} sec^{-1}. If all the dark matter is subject to decay, this allows us to constrain dark matter lifetimes <10^{27} sec. Such energy injection rates are much smaller than those typically probed by the CMB power spectra. The expected brightness temperature fluctuations at z~50 are a fraction of a mK and can vary from the standard calculation by up to an order of magnitude, although the difference can be significantly smaller if some of the decay products free stream to lower redshifts. For self-annihilating dark matter, the fluctuation amplitude can differ by a factor <2 from the standard calculation at z~50. Note also that, in contrast to the CMB, the 21 cm probe is sensitive to both the ionization fraction and the IGM temperature, in principle allowing better constraints on the decay process and heating history. We also show that strong IGM heating and ionization can lead to an enhanced H_2 abundance, which may affect the earliest generations of stars and galaxies.

The spin temperature of neutral hydrogen, which determines the 21 cm optical depth and brightness temperature, is set by the competition between radiative and collisional processes. In the high-redshift intergalactic medium, the dominant collisions are typically those between hydrogen atoms. However, collisions with electrons couple much more efficiently to the spin state of hydrogen than do collisions with other hydrogen atoms and thus become important once the ionized fraction exceeds ~1%. Here we compute the rate at which electron-hydrogen collisions change the hydrogen spin. Previous calculations included only S-wave scattering and ignored resonances near the n=2 threshold. We provide accurate results, including all partial wave terms through the F-wave, for the de-excitation rate at temperatures T_K < 15,000 K; beyond that point, excitation to n>=2 hydrogen levels becomes significant. Accurate electron-hydrogen collision rates at higher temperatures are not necessary, because collisional excitation in this regime inevitably produces Lyman-alpha photons, which in turn dominate spin exchange when T_K > 6200 K even in the absence of radiative sources. Our rates differ from previous calculations by several percent over the temperature range of interest. We also consider some simple astrophysical examples where our spin de-excitation rates are useful.

We re-examine scattering of photons near the Lyman-alpha resonance in the intergalactic medium (IGM). We first derive a general integral solution for the radiation field around resonance when spin diffusivity is ignored. Our solution shows explicitly that recoil sources an absorption feature, whose magnitude increases with the relative importance of recoil compared to Doppler broadening. This spectrum depends on the Lyman-alpha line profile, but approximating it with the absorption profile appropriate to the Lorentzian wings of natural broadening accurately reproduces the results for a full Voigt profile so long as T<1000 K in the IGM. This approximation allows us to obtain simple analytic formulae for the total scattering rate of Lyman-alpha photons and the accompanying energy exchange rate. Our power series solutions converge rapidly for photons that redshift into the Lyman-alpha resonance as well as for photons injected at line center. We confirm previous calculations showing that heating through this mechanism is quite slow and probably negligible compared to other sources. We then show that energy exchange during the scattering of higher-order Lyman-series photons can be much more important than naively predicted by recoil arguments. However, the resulting heating is still completely negligible.

We examine the prospects for studying the pre-reionization intergalactic medium (IGM) through the so-called 21 cm forest in spectra of bright high-redshift radio sources. We first compute the evolution of the mean optical depth for models that include X-ray heating of the IGM gas, Wouthuysen-Field coupling, and reionization. Under most circumstances, the spin temperature T_S grows large well before reionization begins in earnest. As a result, the optical depth is less than 0.001 throughout most of reionization, and background sources must sit well beyond the reionization surface in order to experience measurable absorption. HII regions produce relatively large "transmission gaps" and may therefore still be observable during the early stages of reionization. Absorption from sheets and filaments in the cosmic web fades once T_S becomes large and should be rare during reionization. Minihalos can produce strong (albeit narrow) absorption features. Measuring their abundance would yield useful limits on the strength of feedback processes in the IGM as well as their effect on reionization.

Before reionization, the intergalactic medium (IGM) may have been sufficiently cold for low-mass "minihalos" to condense out of the gas and subsequently affect reionization. Previous work has shown that minihalos generate reasonably large 21 cm fluctuations. Here we consider this signal in its proper cosmological context and show that isolating minihalos from the rest of the IGM is extremely difficult. Using the well-known halo model, we compute the power spectrum of 21 cm fluctuations from minihalos and show that the signal decreases rapidly as feedback increases the Jeans mass. We then show that even a small Lyman-alpha background increases the 21 cm fluctuations of the diffuse IGM well beyond those of the minihalos; because the mass fraction in the IGM is much larger, minihalos will lie buried within the IGM signal. The distinctive signatures of non-linear bias and minihalo structure emerge only at much smaller scales, well beyond the resolution of any upcoming instruments. Using simple, but representative, reionization histories, we then show that the required Lyman-alpha background level is most likely achieved at z>15, while minihalos are still rare, so that they are almost always degenerate with the diffuse IGM.

We consider the evolution of the sky-averaged 21 cm background during the early phases of structure formation. Using simple analytic models, we calculate the thermal and ionization histories, assuming that stellar photons dominate the radiation background. The resulting 21 cm spectra can constrain the properties of the first generations of stars and quasars. If Population II stars dominate, Lyman-alpha coupling renders the IGM visible before it is heated by X-rays and long before reionization. Thus the 21 cm background has a strong absorption epoch followed by weaker emission that fades during reionization. The harder spectra of very massive Population III stars compress these transitions into a shorter time interval and decreases the signal amplitude. However, the reionization epoch will remain visible except in extreme cases. The global 21 cm signal will be challenging to observe because of astronomical foregrounds, but it offers an exciting opportunity to study the first sources of light. It also fixes the overall amplitude of the fluctuating background whose detection is a major goal of several next-generation low-frequency radio interferometers.

Cosmological hydrodynamic simulations have significantly improved over the past several years, and we have already shown that the observed properties of Lyman-break galaxies (LBGs) at z=3 can be explained well by the massive galaxies in the simulations. Here we extend our study to z=6 and show that we obtain good agreement for the LBGs at the bright-end of the luminosity function (LF). Our simulations also suggest that the cosmic star formation rate density has a peak at z= 5-6, and that the current LBG surveys at z=6 are missing a significant number of faint galaxies that are dimmer than the current magnitude limit. Together, our results suggest that the universe could be reionized at z=6 by the Pop II stars in ordinary galaxies. We also estimate the LF of Lyman-alpha emitters (LAEs) at z=6 by relating the star formation rate in the simulation to the Ly-alpha luminosity. We find that the simulated LAE LFs agree with the observed data provided that the net escape fraction of Ly-alpha photon is f_{Ly-alpha} <= 0.1. We investigate two possible scenarios for this effect: (1) all sources in the simulation are uniformly dimmer by a factor of 10 through attenuation, and (2) one out of ten LAEs randomly lights up at a given moment. We show that the correlation strength of the LAE spatial distribution can possibly distinguish the two scenarios.

We examine the effect of Lyman continuum photons on the 21 cm background in the high-redshift universe. The brightness temperature of this transition is determined by the spin temperature T_s, which describes the relative populations of the singlet and triplet hyperfine states. Once the first luminous sources appear, T_s is set by the Wouthuysen-Field effect, in which Lyman-series photons mix the hyperfine levels. Here we consider coupling through n>2 Lyman photons. We first show that coupling (and heating) from scattering of Lyman-n photons is negligible, because they rapidly cascade to lower-energy photons. These cascades can result in either a Lyman-alpha photon -- which will then affect T_s according to the usual Wouthuysen-Field mechanism -- or photons from the 2s -> 1s continuum, which escape without scattering. We show that a proper treatment of the cascades delays the onset of strong Wouthuysen-Field coupling and affects the power spectrum of brightness fluctuations when the overall coupling is still relatively weak (i.e., around the time of the first stars). Cascades damp fluctuations on small scales because only ~ 1/3 of Lyn photons cascade through Lyman-alpha, but they do not affect the large-scale power because that arises from those photons that redshift directly into the Lyman-alpha transition. We also comment on the utility of Lyman-n transitions in providing "standard rulers" with which to study the high-redshift universe.

The Lyman-alpha forest at z >~ 5.5 shows strong scatter in the mean transmission even when smoothed over very large spatial scales, >~ 50 Mpc/h. This has been interpreted as a signature of strongly fluctuating radiation fields, or patchy reionization. To test this claim, we calculate the scatter arising solely from density fluctuations, with a uniform ionizing background, using analytic arguments and simulations. This scatter alone is comparable to that observed. It rises steeply with redshift and is of order unity by z ~ 6, even on ~ 50 Mpc/h scales. This arises because: i) at z ~ 6, transmission spectra, which are sensitive mainly to rare voids, are highly biased (with a linear bias factor b >~ 4-5) tracers of underlying density fluctuations, and ii) projected power from small-scale transverse modes is aliased to long wavelength line-of-sight modes. Inferring patchy reionization from quasar spectra is therefore subtle and requires much more detailed modeling. Similarly, we expect order unity transmission fluctuations in the z ~ 3 HeII Lyman-alpha forest from density fluctuations alone, on the scales over which these measurements are typically made.

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.

Searches for Lyman-alpha emission lines are among the most effective ways to identify high-redshift galaxies. They are particularly interesting because they probe not only the galaxies themselves but also the ionization state of the intergalactic medium (IGM). The observed line strengths depend on the amount by which each photon is able to redshift away from line center before encountering neutral gas and hence on the size distribution of HII regions surrounding the sources. Here, we use an analytic model of that size distribution to study the effects of reionization on the luminosity function of Lyman-alpha emitters and their observed spatial distribution. Our model includes the clustering of high-redshift galaxies and thus contains ionized bubbles much larger than those expected around isolated galaxies. As a result, Lyman-alpha emitting galaxies remain visible earlier in reionization: we expect the number counts to decline by only a factor ~2 (or 10) when the mean ionized fraction falls to x_i~0.75 (or 0.5) in the simplest model. Moreover, the absorption is not uniform across the sky: galaxies remain visible only if they sit inside large bubbles, which become increasingly rare as x_i decreases. Thus, the size distribution also affects the apparent clustering of Lyman-alpha selected galaxies. On large scales, it traces that of the large bubbles, which in our model are more biased than the galaxies. On small scales, the clustering increases rapidly as x_i decreases because large HII regions surround strong galaxy overdensities, so a survey automatically selects only those galaxies with neighbours. The transition between these two regimes occurs at the characteristic bubble size. Hence, large Lyman-alpha galaxy surveys have the potential to measure directly the size distribution of HII regions during reionization.

We study analytically the possibility that mergers of haloes are more highly clustered than the general population of haloes of comparable masses. We begin by investigating predictions for merger bias within the extended Press-Schechter formalism and discuss the limitations and ambiguities of this approach. We then postulate that mergers occur whenever two objects form within a (small) fixed distance of each other. We therefore study the clustering of pairs of points for a highly biased population in the linear regime, for the overall mass distribution in the quasilinear regime, and (using the halo model of clustering) in the nonlinear regime. Biasing, quasilinear evolution, and nonlinear clustering all lead to nonzero reduced (or connected) three-point and four-point correlation functions. These higher-order correlation functions can in many cases enhance the clustering of close pairs of points relative to the clustering of individual points. If close pairs are likely to merge, then the clustering of mergers may be enhanced. We discuss implications for the observed clustering of luminous z=3 galaxies and for correlations of active galactic nuclei and galaxy clusters.

We present an analytic model for the sizes of voids in the galaxy distribution. Peebles and others have recently emphasized the possibility that the observed characteristics of voids may point to a problem in galaxy formation models, but testing these claims has been difficult without any clear predictions for their properties. In order to address such questions, we build a model to describe the distribution of galaxy underdensities. Our model is based on the "excursion set formalism," the same technique used to predict the dark matter halo mass function. We find that, because of bias, galaxy voids are typically significantly larger than dark matter voids and should fill most of the universe. We show that voids selected from catalogs of luminous galaxies should be larger than those selected from faint galaxies: the characteristic radii range from ~5-10h^{-1} Mpc for galaxies with absolute r-band magnitudes M_r-5 log(h) < -16 to -20. These are reasonably close to, though somewhat smaller than, the observed sizes. The discrepancy may result from the void selection algorithm or from their internal structure. We also compute the halo populations inside voids. We expect small haloes (M<10^{11} Msun) to be up to a factor of two less underdense than the haloes of normal galaxies. Within large voids, the mass function is nearly independent of the size of the underdensity, but finite-size effects play a significant role in small voids (R<7h^{-1} Mpc).

The presence of multiple luminous galaxies in clusters can be explained by the finite time over which a galaxy sinks to the center of the cluster and merges with the the central galaxy. The simplest measurable statistic to quantify the dynamical age of a system of galaxies is the luminosity (magnitude) gap, which is the difference in photometric magnitude between the two most luminous galaxies. We present a simple analytical estimate of the luminosity gap distribution in groups and clusters as a function of dark matter halo mass. The luminosity gap is used to define "fossil" groups; we expect the fraction of fossil systems to exhibit a strong and model-independent trend with mass: ~1-3% of massive clusters and ~5-40% of groups should be fossil systems. We also show that, on cluster scales, the observed intrinsic scatter in the central galaxy luminosity-halo mass relation can be ascribed to dispersion in the merger histories of satellites within the cluster. We compare our predictions to the luminosity gap distribution in a sample of 730 clusters in the Sloan Digital Sky Survey C4 Catalog and find good agreement. This suggests that theoretical excursion set merger probabilities and the standard theory of dynamical segregation are valid on cluster scales.

Reionization is inhomogeneous for two reasons: the clumpiness of the intergalactic medium (IGM) and clustering of the discrete ionizing sources. While numerical simulations can in principle take both into account, they are at present limited by small box sizes. On the other hand, analytic models have only examined the limiting cases of a clumpy IGM (with uniform ionizing emissivity) and clustered sources (embedded in a uniform IGM). Here, we present an analytic model for the evolving topology of reionization that includes both factors. At first, recombinations can be ignored and ionized bubbles grow primarily through major mergers. As a result, reionization resembles "punctuated equilibrium," with a series of well-separated sharp jumps in the ionizing background. These features are local effects and do not reflect similar jumps in the global ionized fraction. We then combine our bubble model with a simple description of recombinations in the IGM. We show that the bubbles stop growing when recombinations balance ionizations. If the IGM density structure is similar to that at moderate redshifts, this limits the bubble radii to ~20 comoving Mpc; however, it may be much smaller if the IGM is significantly clumpier at higher redshifts. Once a bubble reaches saturation, that region of the universe has for all intents and purposes entered the "post-overlap" stage, so the overlap epoch actually has a finite width. This picture naturally explains the substantial large-scale variation in Lyman-series opacity along the lines of sight to the known z>6 quasars. More quasar spectra will shed light on the transition between the "bubble-dominated" topology characteristic of reionization and the "web-dominated" topology characteristic of the later universe. [Abridged]

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.

We discuss how measurements of fluctuations in the absorption of cosmic microwave background (CMB) photons by neutral gas during the cosmic dark ages, at redshifts z ~ 7--200, could reveal the primordial deuterium abundance of the Universe. The strength of the cross-correlation of brightness-temperature fluctuations due to resonant absorption of CMB photons in the 21-cm line of neutral hydrogen with those due to resonant absorption of CMB photons in the 92-cm line of neutral deuterium is proportional to the fossil deuterium to hydrogen ratio [D/H] fixed during big bang nucleosynthesis (BBN). Although technically challenging, this measurement could provide the cleanest possible determination of [D/H], free from contamination by structure formation processes at lower redshifts, and has the potential to improve BBN constraints to the baryon density of the Universe \Omega_{b} h^2. We also present our results for the thermal spin-change cross-section for deuterium-hydrogen scattering, which may be useful in a more general context than we describe here.

Recent ultraviolet and X-ray observations of intergalactic OVI and OVII absorption systems along lines of sight to bright quasars have opened a new window onto the ``warm-hot intergalactic medium'' (WHIM). These systems appear to provide a significant reservoir for baryons in the local universe, and comparison to cosmological simulations suggests that their abundance roughly matches theoretical predictions. Here we use analytic arguments to elucidate the physical properties of the absorbers and their role in structure formation. We first show that if the absorbers result from structure-formation shocks, the observed column densities naturally follow from postshock-cooling models, if we include fast-cooling shocks as well as those that cannot cool within a Hubble time. In this case, the known OVI absorbers should show stronger OVII absorption than expected from collisional-ionization equilibrium (and much more than expected for photoionized systems). We then argue that higher-temperature shocks will be spatially associated with more massive virialized objects even well outside the virial radius. Thus the different oxygen ions will trace different structures; OVII absorbers are the most common because that ion dominates over a wide temperature range (corresponding to a large range in halo mass). If each dark-matter halo is surrounded by a network of shocks with total cross section a few times the size of the virialized systems, then we can reproduce the observed number densities of absorbers with plausible parameters. A simple comparison with simulations shows that these assumptions are reasonable, although the actual distribution of shocked gas is too complex for analytic models to describe fully. Our models suggest that these absorbers cannot be explained as a single-temperature phase.

Recent observations of z~6 quasars and the cosmic microwave background imply a complex history to cosmic reionization. Such a history requires some form of feedback to extend reionization over a long time interval, but the nature of the feedback and how rapidly it operates remain highly uncertain. Here we focus on one aspect of this complexity: which physical processes can cause the global ionized fraction to evolve non-monotonically with cosmic time? We consider a range of mechanisms and conclude that double reionization is much less likely than a long, but still monotonic, ionization history. We first examine how galactic winds affect the transition from metal-free to normal star formation. Because the transition is actually spatially inhomogeneous and temporally extended, this mechanism cannot be responsible for double reionization given plausible parameters for the winds. We next consider photoheating, which causes the cosmological Jeans mass to increase in ionized regions and hence suppresses galaxy formation there. In this case, double reionization requires that small halos form stars efficiently, that the suppression from photoheating is strong relative to current expectations, and that ionizing photons are preferentially produced outside of previously ionized regions. Finally, we consider H_2 photodissociation, in which the buildup of a soft ultraviolet background suppresses star formation in small halos. This can in principle cause the ionized fraction to temporarily decrease, but only during the earliest stages of reionization. Finally, we briefly consider the effects of some of these feedback mechanisms on the topology of reionization.

We study the temperature structure of the intergalactic medium (IGM) using a large cosmological N-body/SPH simulation. We employ a two-temperature model for the thermal evolution of the ionized gas, in which the relaxation process between electrons and ions is explicitly included. In the diffuse, hot IGM, the relaxation time is comparable to the age of the Universe and hence the electron temperature in post-shock regions remains significantly smaller than the ion temperature. We show that, at the present epoch, a large fraction of the warm/hot intergalactic medium (WHIM) has a well-developed two temperature structure, with typical temperature differences of order a factor of a few. Consequently, the fraction of metals in various ionization states such as OVI, OVII, and OVIII, as well as their line emissivities, can differ locally by more than an order of magnitude from those computed with a single temperature model. It is thus necessary to follow the evolution of the electron temperature explicitly to determine absorption and emission by the WHIM. Although equipartition is nearly achieved in the denser intracluster medium (ICM), we find an appreciable systematic deviation between the gas-mass weighted electron temperature and the mean temperature even at half the virial radii of clusters. There is thus a reservoir of warm (Te < 1keV) gas in and around massive clusters. Our results imply that relaxation processes need to be considered in describing and interpreting observational data from existing X-ray telescopes as well as from future missions designed to detect the WHIM, such as the Diffuse Intergalactic Oxygen Surveyor and the Missing Baryon Explorer.

Detectable flux is visible in the Ly-alpha and Ly-beta troughs of the highest redshift ($z=6.42$) quasar found to date, SDSS J1148+5251. This has previously been interpreted as continuum contamination from an interloper galaxy at z=4.94. We examine the Ly-gamma trough of SDSS J1148+5251 and show that this interpretation is untenable: the spectrum does not show the continuum break in a z=4.94 galaxy expected from absorption by the intervening Ly-alpha forest. Therefore, flux must be leaking through at least one of the troughs from the quasar itself. Contrary to previous claims, the flux ratios in the Ly-alpha and Ly-beta troughs are consistent with pure transmission. From the Ly-gamma trough, we place an upper bound on the effective Ly-alpha optical depth at z~6.2 of tau < 14.3 (2 sigma). This implies a highly ionized IGM along this line of sight and significant cosmic variance in the transition toward complete Gunn-Peterson absorption. Detailed study of the observed transmission features will shed light on this era.

The nature of the interaction between galaxies and the intergalactic medium (IGM) is one of the most fundamental problems in astrophysics. The accretion of gas onto galaxies provides fuel for star formation, while galactic winds transform the nearby IGM in a number of ways. One exciting technique to study this gas is through the imaging of hydrogen Lyman-alpha emission. We use cosmological simulations to study the Lyman-alpha signals expected from the growth of cosmic structure from z=0-5. We show that if dust absorption is negligible, recombinations following the absorption of stellar ionizing photons dominate the total Lyman-alpha photon production rate. However, galaxies are also surrounded by "Lyman-alpha coronae" of diffuse IGM gas. These coronae are composed of a combination of accreting gas and material ejected from the central galaxy by winds. The Lyman-alpha emission from this phase is powered by a combination of gravitational processes and the photoionizing background. While the former dominates at z~0, collisional excitation following photo-heating may well dominate the total emission at higher redshifts. The central regions of these systems are dense enough to shield themselves from the metagalactic ionizing background; unfortunately, in this regime our simulations are no longer reliable. We therefore consider several scenarios for the emission from the central cores, including one in which self-shielded gas does not emit at all. We show that the combination of star formation and cooling IGM gas can explain most of the observed "Lyman-alpha blobs" at z~3, with the important exception of the largest sources. On the other hand, except under the most optimistic assumptions, cooling IGM gas cannot explain the observations on its own.

The relatively large Thomson optical depth, tau_e, inferred recently from the WMAP observations suggests that the Universe was reionized in a more complex manner than previously believed. However, the value of tau_e provides only an integral constraint on the history of reionization and, by itself, cannot be used to determine the nature of the sources responsible for this transition. Here, we show that the evolution of the ionization state of the intergalactic medium at high redshifts can be measured statistically using fluctuations in 21 centimeter radiation from neutral hydrogen. By analogy with the mathematical description of anisotropies in the cosmic microwave background, we develop a formalism to quantify the variations in 21 cm emission as a function of both frequency and angular scale. Prior to and following reionization, fluctuations in the 21 cm signal are mediated by density perturbations in the distribution of matter. Between these epochs, pockets of gas surrounding luminous objects become ionized, producing large HII regions. These ``bubbles'' of ionized material imprint features into the 21 cm power spectrum that make it possible to distinguish them from fluctuations produced by the density perturbations. The variation of the power spectrum with frequency can be used to infer the evolution of this process. As has been emphasized previously by others, the absolute 21 cm signal from neutral gas at high redshifts is undetectable owing to contamination by foreground sources. However, we argue that this source of noise can be suppressed by comparing maps closely spaced in frequency, i.e. redshift, so that 21 cm fluctuations from the IGM can be measured against a much brighter, but smoothly varying (in frequency) background.

We discuss the prospects for ``tomography'' of the intergalactic medium (IGM) at high redshifts using the 21 cm transition of neutral hydrogen. Existing observational constraints on the epoch of reionization imply a complex ionization history that may require multiple generations of sources. The 21 cm transition provides a unique tool to probe this era in detail, because it does not suffer from saturation effects, retains full redshift information, and directly probes the IGM gas. Observations in the redshifted 21cm line will allow one to study the history and morphology of reionization in detail. Depending on the characteristics of the first sources, they may also allow us to probe the era before reionization, when the first structures and luminous sources were forming. The construction of high signal-to-noise ratio maps on arcminute scales will require approximately one square kilometer of collecting area.

The epoch of reionization (EoR) sets a fundamental benchmark in cosmic structure formation, corresponding to the formation of the first luminous objects that act to ionize the neutral intergalactic medium (IGM). Recent observations at near-IR and radio wavelengths imply that we are finally probing into this key epoch of galaxy formation at $z \ge 6$. The Square Kilometer Array (SKA) will provide critical insight into the EoR in a number of ways. First, the ability of the SKA to image the neutral IGM in 21cm emission is a truly unique probe of the process of reionization, and is recognized as the next necessary and fundamental step in our study of the evolution of large scale structure and cosmic reionization. Second, study of HI 21cm absorption toward the first radio loud objects probes small to intermediate scale structure in the neutral 'cosmic web', as well as HI in the first collapsed structures (proto-disks and mini-halos). And third, the incomparable sensitivity of the SKA allows for the study of the molecular gas, dust, and star formation activity in the first galaxies, as well as the radio continuum emission from the first accreting massive black holes. Such objects will be obscured at optical wavelengths due to absorption by the neutral IGM.

We explore the possibility of detecting HI 21cm absorption by the neutral intergalactic medium (IGM) toward very high redshift radio sources, and by gas associated with the first collapsed structures, using the Square Kilometer Array at low frequency (100 to 200 MHz). The epoch considered is between the time when the first ionizing sources form and when the bulk of the neutral IGM becomes ionized. Expected IGM absorption signal includes $\sim 1%$ absorption by the mean neutral IGM (the radio 'Gunn-Peterson' effect'), plus deeper, narrow lines ($\ge 5%$, a few km s$^{-1}$) arising in mild density inhomogeneities with typical values of cosmic overdensity $\delta\sim 10$, precisely the structures that at later times give rise to the Ly$\alpha$ forest (the `21cm forest'). Absorption can also arise in gas associated with collapsed structures ($\delta\ge 100$), including 'minihalos' ($\le 10^7$ M$_\odot$) and protodisks ($\ge 10^8$ M$_\odot$). We consider SKA sensitivity limits and the evolution of radio source populations, and conclude that it is reasonable to hypothesize the existence of an adequate number of high-$z$ radio sources against which such absorption studies could be performed, provided that reionization occurs at $z < 10$. Lastly, we discuss the possibility of `line confusion' due to radio recombination lines arising in the ionized IGM.

While emission and absorption lines of the 21 cm spin-flip transition of neutral Hydrogen are intrinsically unpolarized, a magnetic field creates left- and right-handed polarized components through the Zeeman effect. Here we consider the resulting polarization of the redshifted 21 cm background from the intergalactic medium before reionization. The polarization is detectable in regions with a strong gradient in the mean brightness temperature. In principle, this can open a new window on the evolution of intergalactic magnetic fields. One possible approach is an extended integration of an individual target during this era, such as the Mpc-scale HII regions inferred to surround quasars at z ~ 6.5. The differential intensity between the two polarization states can be used as a probe of the magnetic field at the edge of the HII region. We estimate that the SKA could (ignoring systematics) detect B \~ 200 (10) \muG coherent over several kiloparsecs with an observational bandwidth of 100 (2) kHz. Beyond individual sources, the statistical properties of wide-field 21 cm polarization maps, such as the angular power spectrum, can be used to constrain the large-scale magnetic field. In this case, the SKA can detect $B \sim 100 \mu$G fields coherent over many megaparsecs. The magnetic field can be measured in any epoch over which the 21 cm background changes rapidly (for example because the ionized fraction or spin temperature change). Although the resulting constraints with SKA are relatively weak compared to theoretical expectations, they nevertheless offer a unique direct probe of magnetic fields in the high-redshift universe.

The reionization of hydrogen in the intergalactic medium (IGM) is a crucial landmark in the history of the universe, but the processes through which it occurs remain mysterious. In particular, recent numerical and analytic work suggest that reionization by stellar sources is driven by large-scale density fluctuations and must be inhomogeneous on scales of many comoving Mpc. We examine the prospects for constraining the topology of neutral and ionized gas through Lyman-alpha absorption of high-redshift sources. One method is to search for gaps in the Gunn-Peterson absorption troughs of luminous sources. These could occur if the line of sight passes sufficiently close to the center of a large HII region. In contrast to previous work, we find a non-negligible (though still small) probability of observing such a gap before reionization is complete. In our model the transmission spike at z=6.08 in the spectrum of SDSS J1148+5251 does not necessarily require overlap to have been completed at an earlier epoch. We also examine the IGM damping wing absorption of the Lyman-alpha emission lines of star-forming galaxies. Because most galaxies sit inside of large HII regions, we find that the severity of absorption is significantly smaller than previously thought and decoupled from the properties of the observed galaxy. While this limits our ability to constrain the mean neutral fraction of the IGM from observations of individual galaxies, it presents the exciting possibility of measuring the size distribution and evolution of the ionized bubbles by examining the distribution of damping wing optical depths in a large sample of galaxies.

We consider the degree to which "21 cm tomography" of the high-redshift Universe can distinguish different ionization histories. Using a new analytic model for the size distribution of HII regions that associates these ionized bubbles with large-scale galaxy overdensities, we compute the angular power spectrum and other statistical properties of the 21 cm brightness temperature during reionization. We show that the HII regions imprint features on the power spectrum that allow us to separate histories with discrete bubbles from those with partial uniform ionization (by, for example, X-rays). We also show that "double reionization" scenarios will modify the morphology of the bubbles in ways that depend on the mechanism through which the first generation of sources shuts off. If, for example, the transition occurs globally at a fixed redshift, the first generation imprints a persistent feature on the 21 cm power spectrum. Finally, we compare our model to one in which voids are ionized first. While the power spectra of these two models are qualitatively similar, we show that the underlying distributions of neutral hydrogen differ dramatically and suggest that other statistical tests can distinguish them. The next generation of low-frequency radio telescopes will have the sensitivity to distinguish all of these models and strongly constrain the history and morphology of reionization.

Recently, there has been a great deal of interest in understanding the reionization of hydrogen in the intergalactic medium (IGM). One of the major outstanding questions is how this event proceeds on large scales. Motivated by numerical simulations, we develop a model for the growth of HII regions during the reionization era. We associate ionized regions with large-scale density fluctuations and use the excursion set formalism to model the resulting size distribution. We then consider ways in which to characterize the morphology of ionized regions. We show how to construct the power spectrum of fluctuations in the neutral hydrogen field. The power spectrum contains definite features from the HII regions that should be observable with the next generation of low-frequency radio telescopes through surveys of redshifted 21 cm emission from the reionization era. Finally, we also consider statistical descriptions beyond the power spectrum and show that our model of reionization qualitatively changes the distribution of neutral gas in the IGM.

We use a high-resolution cosmological simulation that includes hydrodynamics, multiphase star formation, and galactic winds to predict the distribution of metal line emission at z~0 from the intergalactic medium (IGM). We focus on two ultraviolet doublet transitions, OVI 1032,1038 and CIV 1548,1551. Emission from filaments with moderate overdensities is orders of magnitude smaller than the background, but isolated emission from enriched, dense regions with T~10^5-10^5.5 K and characteristic sizes of 50-100 kpc can be detected above the background. We show that the emission from these regions is substantially greater when we use the metallicities predicted by the simulation (which includes enrichment through galactic winds) than when we assume a uniform IGM metallicity. Luminous regions correspond to volumes that have recently been influenced by galactic winds. We also show that the line emission is clustered on scales ~1 h^-1 Mpc. We argue that although these transitions are not effective tracers of the warm-hot intergalactic medium, they do provide a route to study the chemical enrichment of the IGM and the physics of galactic winds.

We discuss free-free radio emission from ionized gas in the intergalactic medium. Because the emissivity is proportional to the square of the electron density, the mean background is strongly sensitive to the spatial clumping of free electrons. Using several existing models for the clumping of ionized gas, we find that the expected free-free distortion to the cosmic microwave background (CMB) blackbody spectrum is at a level detectable with upcoming experiments such as the Absolute Radiometer for Cosmology, Astrophysics, and Diffuse Emission (ARCADE). However, the dominant contribution to the distortion comes from clumpy gas at z < 3, and the integrated signal does not strongly constrain the epoch of reionization. In addition to the mean emission, we consider spatial fluctuations in the free-free background and the extent to which these anisotropies confuse the search for fluctuations in 21 cm line emission from neutral hydrogen during and prior to reionization. This background is smooth in frequency space and hence can be removed through frequency differencing, but only so long as the 21 cm signal and the free-free emission are uncorrelated. We show that, because the free-free background is generated primarily at low redshifts, the cross-correlation between the two fields is smaller than a few percent. Thus, multifrequency cleaning should be an effective way to eliminate the free-free confusion.

Cosmological simulations show that, at the present time, a substantial fraction of the gas in the intergalactic medium (IGM) has been shock-heated to T>10^5 K. Here we develop an analytic model to describe the fraction of shocked, moderately overdense gas in the IGM. The model is an extension of the Press & Schechter (1974) description for the mass function of halos: we assume that large-scale structure shocks occur at a fixed overdensity during nonlinear collapse. This in turn allows us to compute the fraction of gas at a given redshift that has been shock-heated to a specified temperature. We show that, if strong shocks occur at turnaround, our model provides a reasonable description of the temperature distribution seen in cosmological simulations at z~0, although it does overestimate the importance of weak shocks. We then apply our model to shocks at high redshifts. We show that, before reionization, the thermal energy of the IGM is dominated by large-scale structure shocks (rather than virialized objects). These shocks can have a variety of effects, including stripping ~10% of the gas from dark matter minihalos, accelerating cosmic rays, and creating a diffuse radiation background from inverse Compton and cooling radiation. This radiation background develops before the first stars form and could have measurable effects on molecular hydrogen formation and the spin temperature of the 21 cm transition of neutral hydrogen. Finally, we show that shock-heating will also be directly detectable by redshifted 21 cm measurements of the neutral IGM in the young universe.

We use a high-resolution cosmological simulation to predict the distribution of HI Ly-alpha emission from the low-redshift (z<0.5) intergalactic medium (IGM). Our simulation can be used to reliably compute the emission from optically thin regions of the IGM but not that of self-shielded gas. We therefore consider several models that bracket the expected emission from self-shielded regions. Most galaxies are surrounded by extended (>10^2 kpc) ``coronae'' of optically thin gas with Ly-alpha surface brightness close to the expected background. Most of these regions contain smaller cores of dense, cool gas. Unless self-shielded gas is able to cool to T<10^4.1 K, these cores are much brighter than the background. The Ly-alpha coronae represent ``cooling flows'' of IGM gas accreting onto galaxies. We also estimate the number of Ly-alpha photons produced through the reprocessing of stellar ionizing radiation in the interstellar medium of galaxies; while this mechanism is responsible for the brightest Ly-alpha emission, it occurs on small physical scales and can be separated using high-resolution observations. In all cases, we find that Ly-alpha emitters are numerous (with a space density ~0.1 h^3 Mpc^-3) and closely trace the filamentary structure of the IGM, providing a new way to map gas inside the cosmic web.

We study the observability of the reionization epoch through the 21 cm hyperfine transition of neutral hydrogen. We use a high-resolution cosmological simulation (including hydrodynamics) together with a fast radiative transfer algorithm to compute the evolution of 21 cm emission from the intergalactic medium (IGM) in several different models of reionization. We show that the mean brightness temperature of the IGM drops from dT_b~25 mK to dT_b~0.01 mK during overlap (over a frequency interval of ~25 MHz), while the root mean square brightness temperature fluctuations on small scales drop abruptly from ~10 mK before overlap to ~0.1 mK at the end of overlap. We show that 21 cm observations can efficiently discriminate models with a single early reionization epoch from models with two distinct reionization episodes.

Winds from star-forming galaxies provide the most promising explanation for the enrichment of the intergalactic medium with heavy elements. Theoretical and observational arguments indicate that the pollution may have occurred at z > 6; however, direct observational tests of such a scenario are needed. We model starburst winds in the high-redshift universe and find that the fraction of space filled by enriched material varies strongly with the assumed star formation efficiency f_* and the fraction of supernova energy powering each wind, f_esc. We show that metals carried by these winds can be seen in absorption against bright background sources, such as quasars or gamma-ray bursts, in narrow lines with characteristic equivalent widths ~0.5 < W < 5 Angstrom. We argue that a substantial fraction of the metals in high-redshift winds are likely to reside in low ionization states (CII, OI, SiII, and FeII), but higher ionization states (CIV and SiIV) could also provide useful probes of the winds. The number of such lines can constrain both f_* and f_esc. Statistics of metal absorption lines can also be used to identify whether molecular hydrogen is an efficient coolant in the early universe and to study the initial mass function of stars at high redshifts.

We study the absorption along lines of sight toward high-z radio sources caused by the 21 cm transition of neutral hydrogen in the intergalactic medium (IGM) before reionization. Using semi-analytic methods, we compute the number density of observable features caused by both ``minihalos'' (bound objects that are unable to cool efficiently because of their small virial temperatures) and protogalactic disks. We show that both sets of features should be observable by the next generation of low-frequency radio telescopes, including the Low Frequency Array and the Square Kilometer Array, provided that sufficiently bright background sources exist. The statistics of minihalo absorption features seen along lines of sight to radio-loud quasars offer a way to measure the evolution of the radiation background and the IGM temperature with cosmic time. Intersections with disks are much less common, but also cause much deeper absorption features that would be visible in the spectra of both radio-loud quasars and gamma-ray bursts (GRBs). The absorption feature caused by HI in the host galaxy of a GRB should be observable, offering a route to determine spectroscopically the burst redshift.

We model the annihilation of relic positrons produced in a gamma-ray burst (GRB) after its afterglow has faded. We find that the annihilation signal from at least one GRB remnant in the Milky Way galaxy should be observable with future space missions such as INTEGRAL and EXIST, provided that the gas surrounding the GRB source has the typical density of the interstellar medium, < 1 cm^-3. Three fortunate circumstances conspire to make the signal observable. First, unlike positrons in a standard supernova, the GRB positrons initially travel at a relativistic speed and remain ahead of any non-relativistic ejecta until the ejecta become rarefied and the annihilation time becomes long. Second, the GRB remnant remains sufficiently hot (T > 5 x 10^5 K) for a strong annihilation line to form without significant smearing by three-photon decay of positronium. Third, the annihilation signal persists over a time longer than the average period between GRB events in the Milky Way galaxy.

Clusters of galaxies are enriched with positrons from jets of active galactic nuclei (AGNs) or from the interaction of cosmic rays with the intracluster gas. We follow the cooling of these positrons and show that their eventual annihilation with cluster electrons yields a narrow annihilation line. Unlike annihilation in the interstellar medium of galaxies, the line produced in clusters is not smeared by three-photon decay of positronium, because positronium formation is suppressed at the high (>~ 1 keV) temperature of the cluster electrons. We show that if AGN jets are composed of e^+e^- pairs, then the annihilation line from rich clusters within a distance of 100 Mpc might be detectable with future space missions, such as INTEGRAL or EXIST.

Outflows from quasars inevitably pollute the intergalactic medium (IGM) with magnetic fields. The short-lived activity of a quasar leaves behind an expanding magnetized bubble in the IGM. We model the expansion of the remnant quasar bubbles and calculate their distribution as a function magnetic field strength at different redshifts. We find that by a redshift z ~ 3, about 5-80% of the IGM volume is filled by magnetic fields with an energy density > 10% of the mean thermal energy density of a photo-ionized IGM (at ~ 10^4 K). As massive galaxies and X-ray clusters condense out of the magnetized IGM, the adiabatic compression of the magnetic field could result in the fields observed in these systems without a need for further dynamo amplification.

We propose to use gravitational wakes as a direct observational probe of the collisional nature of the dark matter. We calculate analytically the structure of a wake generated by the motion of a galaxy in the core of an X-ray cluster for dark matter in the highly-collisional and collisionless limits. We show that the difference between these limits can be recovered from detailed X-ray or weak lensing observations. We also discuss the sizes of sub-halos in these limits. Preliminary X-ray data on the motion of NGC 1404 through the Fornax group disfavors fluid-like dark matter but does not exclude scenarios in which the dark matter is weakly collisional.

Outflows from quasars inevitably pollute the intergalactic medium (IGM) with magnetic fields. The short-lived activity of a quasar leaves behind an expanding magnetized bubble in the IGM. We model the expansion of the remnant quasar bubbles and calculate their distribution as a function of size and magnetic field strength at different redshifts. We generically find that by a redshift z=3, about 5-20% of the IGM volume is filled by magnetic fields with an energy density >10% of the mean thermal energy density of a photo-ionized IGM (at T=10^4 K). As massive galaxies and X-ray clusters condense out of the magnetized IGM, the adiabatic compression of the magnetic field could result in the field strength observed in these systems without a need for further dynamo amplification. The intergalactic magnetic field could also provide a nonthermal contribution to the pressure of the photo-ionized gas that may account for the claimed discrepancy between the simulated and observed Doppler width distributions of the Ly-alpha forest.