The high-redshift 21 cm signal from the Epoch of Reionization (EoR) is a promising observational probe of the early universe. Current- and next-generation radio interferometers such as the Hydrogen Epoch of Reionization Array (HERA) and Square Kilometre Array (SKA) are projected to measure the 21 cm auto power spectrum from the EoR. Another observational signal of this era is the kinetic Sunyaev-Zel'dovich (kSZ) signal in the cosmic microwave background (CMB), which will be observed by the upcoming Simons Observatory (SO) and CMB-S4 experiments. The 21 cm signal and the contribution to the kSZ from the EoR are expected to be anti-correlated, the former coming from regions of neutral gas in the intergalactic medium and the latter coming from ionized regions. However, the naive cross-correlation between the kSZ and 21 cm maps suffers from a cancellation that occurs because ionized regions are equally likely to be moving toward or away from the observer and so there is no net correlation with the 21 cm signal. We present here an investigation of the 21 cm-kSZ-kSZ bispectrum, which should not suffer the same cancellation as the simple two-point cross-correlation. We show that there is a significant and non-vanishing signal that is sensitive to the reionization history, suggesting the statistic may be used to confirm or infer the ionization fraction as a function of redshift. In the absence of foreground contamination, we forecast that this signal is detectable at high statistical significance with HERA and SO. The bispectrum we study suffers from the fact that the kSZ signal is sensitive only to Fourier modes with long-wavelength line-of-sight components, which are generally lost in the 21 cm data sets owing to foreground contamination. We discuss possible strategies for alleviating this contamination, including an alternative four-point statistic that may help circumvent this issue.

We present a power spectrum analysis of the ALMA Spectroscopic Survey Large Program (ASPECS LP) data from 84 to 115 GHz. These data predominantly probe small-scale fluctuations ($k=10$-$100$ h Mpc$^{-1}$) in the aggregate CO emission in galaxies at $1 \lesssim z \lesssim 4$. We place an integral constraint on CO luminosity functions (LFs) in this redshift range via a direct measurement of their second moments in the three-dimensional (3D) auto-power spectrum, finding a total CO shot noise power $P_{\textrm{CO,CO}}(k_{\textrm{CO(2-1)}}) \leq 1.9\times10^2$ $\mu$K$^2$ (Mpc h$^{-1}$)$^3$. This upper limit ($3\sigma$) is consistent with the observed ASPECS CO LFs in Decarli et al. 2019, but rules out a large space in the range of $P_{\textrm{CO,CO}}(k_{\textrm{CO(2-1)}})$ inferred from these LFs, which we attribute primarily to large uncertainties in the normalization $\Phi_*$ and knee $L_*$ of the Schechter-form CO LFs at $z > 2$. Also, through power spectrum analyses of ASPECS LP data with 415 positions from galaxies with available optical spectroscopic redshifts, we find that contributions to the observed mean CO intensity and shot noise power of MUSE galaxies are largely accounted for by ASPECS blind detections, though there are $\sim20$% contributions to the CO(2-1) mean intensity due to sources previously undetected in the blind line search. Finally, we sum the fluxes from individual blind CO detections to yield a lower limit on the mean CO surface brightness at 99 GHz of $\langle T_{\textrm{CO}} \rangle = 0.55\pm0.02$ $\mu$K, which we estimate represents $68$-$80$% of the total CO surface brightness at this frequency.

An understanding of astrophysical feedback is important for constraining models of galaxy formation and for extracting cosmological information from current and future weak lensing surveys. The thermal Sunyaev-Zel'dovich effect, quantified via the Compton-$y$ parameter, is a powerful tool for studying feedback, because it directly probes the pressure of the hot, ionized gas residing in dark matter halos. Cross-correlations between galaxies and maps of Compton-$y$ obtained from cosmic microwave background surveys are sensitive to the redshift evolution of the gas pressure, and its dependence on halo mass. In this work, we use galaxies identified in year one data from the Dark Energy Survey and Compton-$y$ maps constructed from Planck observations. We find highly significant (roughly $12\sigma$) detections of galaxy-$y$ cross-correlation in multiple redshift bins. By jointly fitting these measurements as well as measurements of galaxy clustering, we constrain the halo bias-weighted, gas pressure of the Universe as a function of redshift between $0.15 \lesssim z \lesssim 0.75$. We compare these measurements to predictions from hydrodynamical simulations, allowing us to constrain the amount of thermal energy in the halo gas relative to that resulting from gravitational collapse.

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.

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

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

Line-Intensity Mapping is an emerging technique which promises new insights into the evolution of the Universe, from star formation at low redshifts to the epoch of reionization and cosmic dawn. It measures the integrated emission of atomic and molecular spectral lines from galaxies and the intergalactic medium over a broad range of frequencies, using instruments with aperture requirements that are greatly relaxed relative to surveys for single objects. A coordinated, comprehensive, multi-line intensity-mapping experimental effort can efficiently probe over 80% of the volume of the observable Universe - a feat beyond the reach of other methods. Line-intensity mapping will uniquely address a wide array of pressing mysteries in galaxy evolution, cosmology, and fundamental physics. Among them are the cosmic history of star formation and galaxy evolution, the compositions of the interstellar and intergalactic media, the physical processes that take place during the epoch of reionization, cosmological inflation, the validity of Einstein's gravity theory on the largest scales, the nature of dark energy and the origin of dark matter.

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 large-scale structure of the Universe should soon be measured at high redshift during the Epoch of Reionization (EoR) through line-intensity mapping. A number of ongoing and planned surveys are using the 21 cm line to trace neutral hydrogen fluctuations in the intergalactic medium (IGM) during the EoR. These may be fruitfully combined with separate efforts to measure large-scale emission fluctuations from galactic lines such as [CII], CO, H-$\alpha$, and Ly-$\alpha$ during the same epoch. The large scale power spectrum of each line encodes important information about reionization, with the 21 cm power spectrum providing a relatively direct tracer of the ionization history. Here we show that the large scale 21 cm power spectrum can be extracted using only cross-power spectra between the 21 cm fluctuations and each of two separate line-intensity mapping data cubes. This technique is more robust to residual foregrounds than the usual 21 cm auto-power spectrum measurements and so can help in verifying auto-spectrum detections. We characterize the accuracy of this method using numerical simulations and find that the large-scale 21 cm power spectrum can be inferred to an accuracy of within 5% for most of the EoR, reaching 0.6% accuracy on a scale of $k\sim0.1\,\text{Mpc}^{-1}$ at $\left< x_i \right> = 0.36$ ($z = 8.34$ in our model). An extension from two to $N$ additional lines would provide $N(N-1)/2$ cross-checks on the large-scale 21 cm power spectrum. This work strongly motivates redundant line-intensity mapping surveys probing the same cosmological volumes.

The amplitude of redshifted 21 cm fluctuations during the Epoch of Reionization (EoR) is expected to show a distinctive "rise and fall" behavior with decreasing redshift as reionization proceeds. On large scales (k <~ 0.1 Mpc^{-1}) this can mostly be characterized by evolution in the product of the mean 21 cm brightness temperature and a bias factor, b_21(z). This quantity evolves in a distinctive way that can help in determining the average ionization history of the intergalactic medium (IGM) from upcoming 21 cm fluctuation data sets. Here we consider extracting <T_21> b_21(z) using a combination of future redshifted 21 cm and [CII] line-intensity mapping data sets. Our method exploits the dependence of the 21 cm-[CII]-[CII] cross-bispectrum on the shape of triangle configurations in Fourier space. This allows one to determine <T_21> b_21(z) yet, importantly, is less sensitive to foreground contamination than the 21 cm auto-spectrum, and so can provide a valuable cross-check. We compare the results of simulated bispectra with second-order perturbation theory: on large scales the perturbative estimate of <T_21> b_21(z) matches the true value to within 10% for <x_i> <~ 0.8. We consider the 21 cm auto-bispectrum and show that this statistic may also be used to extract the 21 cm bias factor. Finally, we discuss the survey requirements for measuring the cross-bispectrum. Although we focus on the 21 cm-[CII]-[CII] bispectrum during reionization, our method may be of broader interest and can be applied to any two fields throughout cosmic history.

Next-generation 21cm observations will enable imaging of reionization on very large scales. These images will contain more astrophysical and cosmological information than the power spectrum, and hence providing an alternative way to constrain the contribution of different reionizing sources populations to cosmic reionization. Using Convolutional Neural Networks, we present a simple network architecture that is sufficient to discriminate between Galaxy-dominated versus AGN-dominated models, even in the presence of simulated noise from different experiments such as the HERA and SKA.

The EDGES experiment recently announced evidence for a broad absorption feature in the sky-averaged radio spectrum around 78 MHz, as may result from absorption in the 21 cm line by neutral hydrogen at z~15-20. If confirmed, one implication is that the spin temperature of the 21 cm line is coupled to the gas temperature by z=20. The known mechanism for accomplishing this is the Wouthuysen-Field effect, whereby Lyman-alpha photons scatter in the intergalactic medium (IGM) and impact the hyperfine level populations. This suggests that early star formation had already produced a copious Lyman-alpha background by z=20, and strongly constrains models in which the linear matter power spectrum is suppressed on small-scales, since halo and star formation are delayed in such scenarios. Here we consider the case that the dark matter consists of ultra-light axions with macroscopic de Broglie wavelengths (fuzzy dark matter, FDM). We assume that star formation tracks halo formation and adopt two simple models from the current literature for the halo mass function in FDM. We further suppose that the fraction of halo baryons which form stars is less than a conservative upper limit of $f_\star \leq 0.05$, and that ~10^4 Lyman-alpha to Lyman-limit photons are produced per stellar baryon. We find that the requirement that the 21 cm spin temperature is coupled to the gas temperature by $z=20$ places a lower-limit on the FDM particle mass of $m_a \geq 5 \times 10^{-21} {\rm eV}$. The constraint is insensitive to the precise minimum mass of halos where stars form. As the global 21 cm measurements are refined, the coupling redshift could change and we quantify how the FDM constraint would be modified. A rough translation of the FDM mass bound to a thermal relic warm dark matter (WDM) mass bound is also provided.

The relation between the mass and integrated electron pressure of galaxy group and cluster halos can be probed by stacking maps of the thermal Sunyaev-Zel'dovich (tSZ) effect. Perhaps surprisingly, recent observational results have indicated that the scaling relation between integrated pressure and mass follows the prediction of simple, self-similar models down to halo masses as low as $10^{12.5} \, M_{\odot}$. Hydrodynamical simulations that incorporate energetic feedback processes suggest that gas should be depleted from such low-mass halos, thus decreasing their tSZ signal relative to self-similar predictions. Here, we build on the modeling of Vikram, Lidz, and Jain (2017) to evaluate the bias in the interpretation of stacked tSZ measurements due to the signal from correlated halos (the "two-halo" term), which has generally been neglected in the literature. We fit theoretical models to a measurement of the tSZ -- galaxy group cross-correlation function, accounting explicitly for the one- and two- halo contributions. We find moderate evidence of a deviation from self-similarity in the pressure -- mass relation, even after marginalizing over conservative miscentering effects. We explore pressure -- mass models with a break at $10^{14} \, M_{\odot}$, as well as other variants. We discuss and test for sources of uncertainty in our analysis, in particular a possible bias in the halo mass estimates and the coarse resolution of the Planck beam. We compare our findings with earlier analyses by exploring the extent to which halo isolation criteria can reduce the two-halo contribution. Finally, we show that ongoing third-generation CMB experiments will explicitly resolve the one-halo term in low-mass groups; our methodology can be applied to these upcoming data sets to obtain a clear answer to the question of self-similarity and an improved understanding of hot gas in low-mass halos.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

Standard analyses of the reionization history of the universe from Planck cosmic microwave background (CMB) polarization measurements consider only the overall optical depth to electron scattering ($\tau$), and further assume a step-like reionization history. However, the polarization data contain information beyond the overall optical depth, and the assumption of a step-like function may miss high redshift contributions to the optical depth and lead to biased $\tau$ constraints. Accounting for its full reionization information content, we reconsider the interpretation of Planck 2015 Low Frequency Instrument (LFI) polarization data using simple, yet physically-motivated reionization models. We show that these measurements still, in fact, allow a non-negligible contribution from metal-free (Pop-III) stars forming in mini-halos of mass $M \sim 10^5-10^6 M_\odot$ at $z \gtrsim 15$, provided this mode of star formation is fairly inefficient. Our best fit model includes an early, self-regulated phase of Pop-III star formation in which the reionization history has a gradual, plateau feature. In this model, $\sim$20\% of the volume of the universe is ionized by $z \sim 20$, yet it nevertheless provides a good match to the Planck LFI measurements. Although preferred when the full information content of the data is incorporated, this model would spuriously be disfavored in the standard analysis. This preference is driven mostly by excess power from E-mode polarization at multipoles of $10 \lesssim \ell \lesssim 20$, which may reflect remaining systematic errors in the data, a statistical fluctuation, or signatures of the first stars. Measurements from the Planck High Frequency Instrument (HFI) should be able to confirm or refute this hint and future cosmic-variance limited E-mode polarization surveys can provide substantially more information on these signatures

Stacking cosmic microwave background (CMB) maps around known galaxy clusters and groups provides a powerful probe of the distribution of hot gas in these systems via the Sunyaev-Zel'dovich (SZ) effect. A stacking analysis allows one to detect the average SZ signal around low mass halos, and to extend measurements out to large scales, which are too faint to detect individually in the SZ or in X-ray emission. In addition, cross correlations between SZ maps and other tracers of large-scale structure (with known redshifts) can be used to extract the redshift-dependence of the SZ background. Motivated by these exciting prospects, we measure the two-point cross-correlation function between a catalog of $\sim 380,000$ galaxy groups (with redshifts spanning $z=0.01-0.2$) from the Sloan Digital Sky Survey (SDSS) and Compton-y parameter maps constructed by the Planck collaboration. We find statistically significant correlations between the group catalog and Compton-y maps in each of six separate mass bins, with estimated halo masses in the range $10^{11.5-15.5} M_\odot/h$. We compare these measurements with halo models of the SZ signal, which describe the stacked measurement in terms of one-halo and two-halo terms. The one-halo term quantifies the average pressure profile around the groups in a mass bin, while the two-halo term describes the contribution of correlated neighboring halos. For the more massive groups we find clear evidence for the one- and two-halo regimes, while groups with mass below $10^{13} M_\odot/h$ are dominated by the two-halo term given the resolution of Planck data. We use the signal in the two-halo regime to determine the bias-weighted electron pressure of the universe: $\langle b P_e \rangle= 1.50 \pm 0.226 \times 10^{-7}$ keV cm$^{-3}$ (1-$\sigma$) at $z\approx 0.15$.

Line intensity mapping experiments seek to trace large scale structure by measuring the spatial fluctuations in the combined emission, in some convenient spectral line, from individually unresolved galaxies. An important systematic concern for these surveys is line confusion from foreground or background galaxies emitting in other lines that happen to lie at the same observed frequency as the "target" emission line of interest. We develop an approach to separate this "interloper" emission at the power spectrum level. If one adopts the redshift of the target emission line in mapping from observed frequency and angle on the sky to co-moving units, the interloper emission is mapped to the wrong co-moving coordinates. Since the mapping is different in the line of sight and transverse directions, the interloper contribution to the power spectrum becomes anisotropic, especially if the interloper and target emission are at widely separated redshifts. This distortion is analogous to the Alcock-Paczynski test, but here the warping arises from assuming the wrong redshift rather than an incorrect cosmological model. We apply this to the case of a hypothetical [CII] emission survey at z~7 and find that the distinctive interloper anisotropy can, in principle, be used to separate strong foreground CO emission fluctuations. In our models, however, a significantly more sensitive instrument than currently planned is required, although there are large uncertainties in forecasting the high redshift [CII] emission signal. With upcoming surveys, it may nevertheless be useful to apply this approach after first masking pixels suspected of containing strong interloper contamination.

A proposed method for dealing with foreground emission in upcoming 21-cm observations from the epoch of reionization is to limit observations to an uncontaminated window in Fourier space. Foreground emission can be avoided in this way, since it is limited to a wedge-shaped region in $k_{\parallel}, k_{\perp}$ space. However, the power spectrum is anisotropic owing to redshift-space distortions from peculiar velocities. Consequently, the 21-cm power spectrum measured in the foreground avoidance window---which samples only a limited range of angles close to the line-of-sight direction---differs from the full spherically-averaged power spectrum which requires an average over \emph{all} angles. In this paper, we calculate the magnitude of this "wedge bias" for the first time. We find that the bias is strongest at high redshifts, where measurements using foreground avoidance will over-estimate the power spectrum by around 100 per cent, possibly obscuring the distinctive rise and fall signature that is anticipated for the spherically-averaged 21-cm power spectrum. In the later stages of reionization, the bias becomes negative, and smaller in magnitude ($\lesssim 20$ per cent). The effect shows only a weak dependence on spatial scale and reionization topology.

A major goal of observational and theoretical cosmology is to observe the largely unexplored time period in the history of our universe when the first galaxies form, and to interpret these measurements. Early galaxies dramatically impacted the gas around them in the surrounding intergalactic medium (IGM) by photoionzing the gas during the Epoch of Reionization (EoR). This epoch likely spanned an extended stretch in cosmic time: ionized regions formed and grew around early generations of galaxies, gradually filling a larger and larger fraction of the volume of the universe. At some time -- thus far uncertain, but within the first billion years or so after the big bang -- essentially the entire volume of the universe became filled with ionized gas. The properties of the IGM provide valuable information regarding the formation time and nature of early galaxy populations, and many approaches for studying the first luminous sources are hence based on measurements of the surrounding intergalactic gas. The prospects for improved reionization-era observations of the IGM and early galaxy populations over the next decade are outstanding. Motivated by this, we review the current state of models of the IGM during reionization. We focus on a few key aspects of reionization-era phenomenology and describe: the redshift evolution of the volume-averaged ionization fraction, the properties of the sources and sinks of ionizing photons, along with models describing the spatial variations in the ionization fraction, the ultraviolet radiation field, the temperature of the IGM, and the gas density distribution.

Determining when and how the first galaxies reionized the intergalactic medium (IGM) promises to shed light on both the nature of the first objects and the cosmic history of baryons. Towards this goal, quasar absorption lines play a unique role by probing the properties of diffuse gas on galactic and intergalactic scales. In this review we examine the multiple ways in which absorption lines trace the connection between galaxies and the IGM near the reionization epoch. We first describe how the Ly$\alpha$ forest is used to determine the intensity of the ionizing ultraviolet background and the global ionizing emissivity budget. Critically, these measurements reflect the escaping ionizing radiation from all galaxies, including those too faint to detect directly. We then discuss insights from metal absorption lines into reionization-era galaxies and their surroundings. Current observations suggest a buildup of metals in the circumgalactic environments of galaxies over $z \sim 6$ to 5, although changes in ionization will also affect the evolution of metal line properties. A substantial fraction of metal absorbers at these redshifts may trace relatively low-mass galaxies. Finally, we review constraints from the Ly$\alpha$ forest and quasar near zones on the timing of reionization. Along with other probes of the high-redshift Universe, absorption line data are consistent with a relatively late end to reionization ($5.5 \lesssim z \lesssim 7$); however the constraints are still fairly weak. Significant progress is expected to come through improved analysis techniques, increases in the number of known high-redshift quasars from optical and infrared sky surveys, large gains in sensitivity from next-generation observing facilities, and synergies with other probes of the reionization era.

Infrared and radio observations of the Epoch of Reionization promise to revolutionize our understanding of the cosmic dawn, and major efforts with the JWST, MWA and HERA are underway. While measurements of the ionizing sources with infrared telescopes and the effect of these sources on the intergalactic medium with radio telescopes \emph{should} be complementary, to date the wildly disparate angular resolutions and survey speeds have made connecting proposed observations difficult. In this paper we develop a method to bridge the gap between radio and infrared studies. While the radio images may not have the sensitivity and resolution to identify individual bubbles with high fidelity, by leveraging knowledge of the measured power spectrum we are able to separate regions that are likely ionized from largely neutral, providing context for the JWST observations of galaxy counts and properties in each. By providing the ionization context for infrared galaxy observations, this method can significantly enhance the science returns of JWST and other infrared observations.

Observations of the Lyman-alpha (Ly-$\alpha$) forest may allow reionization to complete as late as $z \sim 5.5$, provided the ionization state of the intergalactic medium (IGM) is sufficiently inhomogeneous at these redshifts. In this case, significantly neutral islands may remain amongst highly ionized gas with the ionized regions allowing some transmission through the Ly-$\alpha$ forest. This possibility has the important virtue that it is eminently testable with existing Ly-$\alpha$ forest data. In particular, we describe three observable signatures of significantly neutral gas in the $z \sim 5.5$ IGM. We use mock quasar spectra produced from numerical simulations of reionization to develop these tests. First, we quantify how the abundance and length of absorbed regions in the forest increase with the volume-averaged neutral fraction in our reionization model. Second, we consider stacking the transmission profile around highly absorbed regions in the forest. If and only if there is significantly neutral gas in the IGM, absorption in the damping wing of the Ly-$\alpha$ line will cause the transmission to recover slowly as one moves from absorbed to transmitted portions of the spectrum. Third, the deuterium Ly-$\beta$ line should imprint a small but distinctive absorption feature slightly blueward of absorbed neutral regions in the Ly-$\beta$ forest. We show that these tests can be carried out with existing Keck HIRES spectra at $z \sim 5.5$, with the damping wing being observable for $< x_{\text{HI}} >\gtrsim 0.05$ and the deuterium feature observable with additional high-resolution spectra for $< x_{\text{HI}} >\gtrsim 0.2$.

Infrared fine-structure emission lines from trace metals are powerful diagnostics of the interstellar medium in galaxies. We explore the possibility of studying the redshifted far-IR fine-structure line emission using the three-dimensional (3-D) power spectra obtained with an imaging spectrometer. The intensity mapping approach measures the spatio-spectral fluctuations due to line emission from all galaxies, including those below the individual detection threshold. The technique provides 3-D measurements of galaxy clustering and moments of the galaxy luminosity function. Furthermore, the linear portion of the power spectrum can be used to measure the total line emission intensity including all sources through cosmic time with redshift information naturally encoded. Total line emission, when compared to the total star formation activity and/or other line intensities reveals evolution of the interstellar conditions of galaxies in aggregate. As a case study, we consider measurement of [CII] autocorrelation in the 0.5 < z < 1.5 epoch, where interloper lines are minimized, using far-IR/submm balloon-borne and future space-borne instruments with moderate and high sensitivity, respectively. In this context, we compare the intensity mapping approach to blind galaxy surveys based on individual detections. We find that intensity mapping is nearly always the best way to obtain the total line emission because blind, wide-field galaxy surveys lack sufficient depth and deep pencil beams do not observe enough galaxies in the requisite luminosity and redshift bins. Also, intensity mapping is often the most efficient way to measure the power spectrum shape, depending on the details of the luminosity function and the telescope aperture.

The temperature of the low-density intergalactic medium (IGM) at high redshift is sensitive to the timing and nature of hydrogen and HeII reionization, and can be measured from Lyman-alpha forest absorption spectra. Since the memory of intergalactic gas to heating during reionization gradually fades, measurements as close as possible to reionization are desirable. In addition, measuring the IGM temperature at sufficiently high redshifts should help to isolate the effects of hydrogen reionization since HeII reionization starts later, at lower redshift. Motivated by this, we model the IGM temperature at z>5 using semi-numeric models of patchy reionization. We construct mock Lyman-alpha forest spectra from these models and consider their observable implications. We find that the small-scale structure in the Lyman-alpha forest is sensitive to the temperature of the IGM even at redshifts where the average absorption in the forest is as high as 90%. We forecast the accuracy at which the z~5 IGM temperature can be measured using existing samples of high resolution quasar spectra, and find that interesting constraints are possible. For example, an early reionization model in which reionization ends at z~10 should be distinguishable -- at high statistical significance -- from a lower redshift model where reionization completes at z~6. We discuss improvements to our modeling that may be required to robustly interpret future measurements.

An appealing approach for studying the reionization history of the Universe is to measure the redshift evolution of the Lyman-alpha fraction, the percentage of Lyman-break selected galaxies that emit appreciably in the Ly-alpha line. This fraction is expected to fall-off towards high redshift as the intergalactic medium becomes significantly neutral, and the galaxies' Ly-alpha emission is progressively attenuated. Intriguingly, early measurements with this technique suggest a strong drop in the Ly-alpha fraction near z ~ 7. Previous work concluded that this requires a surprisingly neutral intergalactic medium -- with neutral hydrogen filling more than 50 % of the volume of the Universe -- at this redshift. We model the evolving Ly-alpha fraction using cosmological simulations of the reionization process. Before reionization completes, the simulated Ly-alpha fraction has large spatial fluctuations owing to the inhomogeneity of reionization. Since existing measurements of the Ly-alpha fraction span relatively small regions on the sky, and sample these regions only sparsely, they may by chance probe mostly galaxies with above average Ly-alpha attenuation. We find that this sample variance is not exceedingly large for existing surveys, but that it does somewhat mitigate the required neutral fraction at z ~ 7. Quantitatively, in a fiducial model calibrated to match measurements after reionization, we find that current z = 7 observations require a volume-averaged neutral fraction of x_HI > 0.05 at 95 % confidence level. Hence, we find that the z ~ 7 Ly-alpha fraction measurements do likely probe the Universe before reionization completes but that they do not require a very large neutral fraction.

The statistical properties of the primordial perturbations contain clues about the origins of those fluctuations. Although the Planck collaboration has recently obtained tight constraints on primordial non-gaussianity from cosmic microwave background measurements, it is still worthwhile to mine upcoming data sets in effort to place independent or competitive limits. The ionized bubbles that formed at redshift z~6-20 during the Epoch of Reionization are seeded by primordial overdensities, and so the statistics of the ionization field at high redshift are related to the statistics of the primordial field. Here we model the effect of primordial non-gaussianity on the reionization field. The epoch and duration of reionization are affected as are the sizes of the ionized bubbles, but these changes are degenerate with variations in the properties of the ionizing sources and the surrounding intergalactic medium. A more promising signature is the power spectrum of the spatial fluctuations in the ionization field, which may be probed by upcoming 21 cm surveys. This has the expected 1/k^2 dependence on large scales, characteristic of a biased tracer of the matter field. We project how well upcoming 21 cm observations will be able to disentangle this signal from foreground contamination. Although foreground cleaning inevitably removes the large-scale modes most impacted by primordial non-gaussianity, we find that primordial non-gaussianity can be separated from foreground contamination for a narrow range of length scales. In principle, futuristic redshifted 21 cm surveys may allow constraints competitive with Planck.

We present a new procedure to measure the large-scale carbon monoxide (CO) emissions across cosmic history. As a tracer of large-scale structure (LSS), the CO gas content as a function of redshift can be quantified by its three-dimensional fluctuation power spectra. Furthermore, cross-correlating CO emission with other LSS tracers offers a way to measure the emission as a function of scale and redshift. Here we introduce the model relevant for such a cross-correlation measurement between CO and other LSS tracers, and between different CO rotational lines. We propose a novel use of Cosmic Microwave Background (CMB) data and attempt to extract redshifted CO emissions embedded in the Wilkinson Microwave Anisotropy Probe (WMAP) dataset. We cross-correlate the all-sky WMAP7 data with LSS data sets, namely, the photometric quasar sample and the luminous red galaxy sample from the Sloan Digital Sky Survey Data Release 6 and 7 respectively. We are unable to detect a cross-correlation signal with either CO(1-0) nor CO(2-1) lines, mainly due to the instrumental noise in the WMAP data. However, we are able to rule out models more than three times greater than our more optimistic model. We discuss the cross-correlation signal from the thermal Sunyaev-Zeldovich effect and dust as potential contaminants, and quantify their impact for our CO measurements. We discuss forecasts for current CMB experiments and a hypothetical future CO focused experiment, and propose to cross-correlate CO temperature data with the Hobby-Eberly Telescope Dark Energy Experiment Ly-alpha-emitter sample, for which a signal-to-noise ratio of 58 is possible.

One of the most promising approaches for studying reionization is to use the redshifted 21 cm line. Early generations of redshifted 21 cm surveys will not, however, have the sensitivity to make detailed maps of the reionization process, and will instead focus on statistical measurements. Here we show that it may nonetheless be possible to {\em directly identify ionized regions} in upcoming data sets by applying suitable filters to the noisy data. The locations of prominent minima in the filtered data correspond well with the positions of ionized regions. In particular, we corrupt semi-numeric simulations of the redshifted 21 cm signal during reionization with thermal noise at the level expected for a 500 antenna tile version of the Murchison Widefield Array (MWA), and mimic the degrading effects of foreground cleaning. Using a matched filter technique, we find that the MWA should be able to directly identify ionized regions despite the large thermal noise. In a plausible fiducial model in which ~20% of the volume of the Universe is neutral at z ~ 7, we find that a 500-tile MWA may directly identify as many as ~150 ionized regions in a 6 MHz portion of its survey volume and roughly determine the size of each of these regions. This may, in turn, allow interesting multi-wavelength follow-up observations, comparing galaxy properties inside and outside of ionized regions. We discuss how the optimal configuration of radio antenna tiles for detecting ionized regions with a matched filter technique differs from the optimal design for measuring power spectra. These considerations have potentially important implications for the design of future redshifted 21 cm surveys.

We study the impact of primordial non-Gaussianity generated during inflation on the bias of halos using excursion set theory. We recapture the familiar result that the bias scales as $k^{-2}$ on large scales for local type non-Gaussianity but explicitly identify the approximations that go into this conclusion and the corrections to it. We solve the more complicated problem of non-spherical halos, for which the collapse threshold is scale dependent.

The epoch of reionization is a milestone of cosmological structure formation, marking the birth of the first objects massive enough to yield large numbers of ionizing photons. The mechanism and timescale of reionization remain largely unknown. Measurements of the CMB Doppler effect from ionizing bubbles embedded in large-scale velocity streams (the patchy kinetic Sunyaev-Zel'dovich effect) can constrain the duration of reionization. When combined with large-scale CMB polarization measurements, the evolution of the ionized fraction can be inferred. Using new multi-frequency data from the South Pole Telescope (SPT), we show that the ionized fraction evolved relatively rapidly. For our basic foreground model, we find the kinetic Sunyaev-Zel'dovich power sourced by reionization at l=3000 to be <= 2.1 micro K^2 at 95% CL. Using reionization simulations, we translate this to a limit on the duration of reionization of Delta z <= 4.4 (95% CL). We find that this constraint depends on assumptions about the angular correlation between the thermal Sunyaev-Zel'dovich power and the cosmic infrared background (CIB). Introducing the degree of correlation as a free parameter, we find that the limits on kSZ power weaken to <= 4.9 micro K^2, implying Delta z <= 7.9 (95% CL). We combine the SPT constraint on the duration of reionization with the WMAP7 measurement of the integrated optical depth to probe the cosmic ionization history. We find that reionization ended with 95% CL at z > 7.2 under the assumption of no tSZ-CIB correlation, and z>5.8 when correlations are allowed. Improved constraints from the full SPT data set in conjunction with upcoming Herschel and Planck data should detect extended reionization at >95% CL provided Delta z >= 4. (abbreviated)

We combine the high-resolution Aquarius simulations with three-dimensional models of reionization based on the initial density field of the Aquarius parent simulation, Millennium-II, to study the impact of patchy reionization on the faint satellite population of Milky Way halos. Because the Aquarius suite consists of zoom-in simulations of halos in the Millennium-II volume, we follow the formation of substructure and the growth of reionization bubbles due to the larger environment simultaneously, and thereby determine the reionization redshifts of satellite candidates. We do this for four different reionization models, and also compare results to instantaneous reionization. Using a simple procedure for selecting satellites and assigning luminosities in the simulations, we compare the resulting satellite populations. We find that the overall number of satellites depends sensitively on the reionization model, with a factor of 3-4 variation between the four models for a given host halo, although the difference is entirely in the population of faint satellites (M_V > -10). In addition, we find that for a given reionization model the total number of satellites differs by 10%-20% between the patchy and homogeneous scenarios, provided that the redshift is chosen appropriately for the instantaneous case. However, the halo-halo scatter from the six Aquarius halos is large, up to a factor of 2-3, and so is comparable to the difference between reionization scenarios. In order to use the population of faint dwarf galaxies around the Milky Way as a probe of the local reionization history, then, it is necessary to first better understand the general distribution of substructure around Milky Way-mass halos.

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

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.

In Lyman-alpha forest measurements it is generally assumed that quasars are mere background light sources which are uncorrelated with the forest. Gravitational lensing of the quasars violates this assumption. This effect leads to a measurement bias, but more interestingly it provides a valuable signal. The lensing signal can be extracted by correlating quasar magnitudes with the flux power spectrum and with the flux decrement. These correlations will be challenging to measure but their detection provides a direct measure of how features in the Lyman-alpha forest trace the underlying mass density field. Observing them will test the fundamental hypothesis that fluctuations in the forest are predominantly driven by fluctuations in mass, rather than in the ionizing background, helium reionization or winds. We discuss ways to disentangle the lensing signal from other sources of such correlations, including dust, continuum and background residuals. The lensing-induced measurement bias arises from sample selection: one preferentially collects spectra of magnified quasars which are behind overdense regions. This measurement bias is ~0.1-1% for the flux power spectrum, optical depth and the flux probability distribution. Since the effect is systematic, quantities such as the amplitude of the flux power spectrum averaged across scales should be interpreted with care.

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.

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.

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

The ionizing background determines the ionization balance and the thermodynamics of the cosmic gas. It is therefore a fundamental ingredient to theoretical and empirical studies of both the IGM and galaxy formation. We present here a new calculation of its spectrum that satisfies the empirical constraints we recently obtained by combining state-of-the-art luminosity functions and intergalactic opacity measurements. In our preferred model, star-forming galaxies and quasars each contribute substantially to the HI ionizing field at z<3, with galaxies rapidly overtaking quasars at higher redshifts as quasars become rarer. In addition to our fiducial model, we explore the physical dependences of the calculated background and clarify how recombination emission contributes to the ionization rates. We find that recombinations do not simply boost the ionization rates by the number of reemitted ionizing photons as many of these rapidly redshift below the ionization edges and have a distribution of energies. A simple analytic model that captures the main effects seen in our numerical radiative transfer calculations is given. Finally, we discuss the effects of HeII reionization by quasars on both the spectrum of the ionizing background and on the thermal history of the IGM. In regions that have yet to be reionized, the spectrum is expected to be almost completely suppressed immediately above 54.4 eV while a background of higher-energy (>~0.5 keV) photons permeates the entire universe owing to the frequency-dependence of the photoionization cross section. We provide an analytic model of the heat input during HeII reionization and its effects on the temperature-density relation.

A soft ultraviolet radiation field, 10.2 eV < E <13.6 eV, that permeates neutral intergalactic gas during the Epoch of Reionization (EoR) excites the 2p (directly) and 2s (indirectly) states of atomic hydrogen. Because the 2s state is metastable, the lifetime of atoms in this level is relatively long, which may cause the 2s state to be overpopulated relative to the 2p state. It has recently been proposed that for this reason, neutral intergalactic atomic hydrogen gas may be detected in absorption in its 3-cm fine-structure line (2s_1/2 -> 2p_3/2) against the Cosmic Microwave Background out to very high redshifts. In particular, the optical depth in the fine-structure line through neutral intergalactic gas surrounding bright quasars during the EoR may reach tau~1e-5. The resulting surface brightness temperature of tens of micro K (in absorption) may be detectable with existing radio telescopes. Motivated by this exciting proposal, we perform a detailed analysis of the transfer of Lyman beta,gamma,delta,... radiation, and re-analyze the detectability of the fine-structure line in neutral intergalactic gas surrounding high-redshift quasars. We find that proper radiative transfer modeling causes the fine-structure absorption signature to be reduced tremendously to tau< 1e-10. We therefore conclude that neutral intergalactic gas during the EoR cannot reveal its presence in the 3-cm fine-structure line to existing radio telescopes.

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

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

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

We summarize existing constraints on the epoch of reionization and discuss the observational probes that are sensitive to the process. We focus on the role large scale polarization can play. Polarization probes the integrated optical depth across the entire epoch of reionization. Future missions such as Planck and CMBPol will greatly enhance our knowledge of the reionization history, allowing us to measure the time evolution of the ionization fraction. As large scale polarization probes high redshift activity, it can best constrain models where the Universe was fully or partially ionized at early times. In fact, large scale polarization could be our only probe of the highest redshifts.

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

We investigate the implications of the intergalactic opacity for the evolution of the cosmic UV luminosity density and its sources. Our main constraint is our measurement of the Lya forest opacity at redshifts 2<z<4.2 from 86 high-resolution quasar spectra. In addition, we impose the requirements that HI must be reionized by z=6 and HeII by z~3, and consider estimates of the hardness of the ionizing background from HI to HeII column density ratios. The derived hydrogen photoionization rate is remarkably flat over the Lya forest redshift range covered. Because the quasar luminosity function is strongly peaked near z~2, the lack of redshift evolution indicates that star-forming galaxies likely dominate the photoionization rate at z>~3. Combined with direct measurements of the galaxy UV luminosity function, this requires only a small fraction f_esc~0.5% of galactic hydrogen ionizing photons to escape their source for galaxies to solely account for the entire ionizing background. Under the assumption that the galactic UV emissivity traces the star formation rate, current state-of-the-art observational estimates of the star formation rate density appear to underestimate the total photoionization rate at z~4 by a factor ~4, are in tension with recent determinations of the UV luminosity function, and fail to reionize the Universe by z~6 if extrapolated to arbitrarily high redshift. A theoretical star formation history peaking earlier fits the Lya forest photoionization rate well, reionizes the Universe in time, and is in better agreement with the rate of z~4 gamma-ray bursts observed by Swift. Quasars suffice to doubly ionize helium by z~3 and likely contribute a non-negligible and perhaps dominant fraction of the hydrogen ionizing background at their z~2 peak. [Abridged]

We directly measure the evolution of the intergalactic Lya effective optical depth, tau_eff, over the redshift range 2<z<4.2 from a sample of 86 high-resolution, high-signal-to-noise quasar spectra obtained with the ESI and HIRES spectrographs on Keck, and with the MIKE spectrograph on Magellan. This represents an improvement over previous analyses of the Lya forest from high-resolution spectra in this redshift interval of a factor of two in the size of the data set alone. We pay particular attention to robust error estimation and extensively test for systematic effects. We find that our estimates of the quasar continuum levels in the Lya forest obtained by spline fitting are systematically biased low, with the magnitude of the bias increasing with redshift, but that this bias can be accounted for using mock spectra. The mean fractional error is <1% at z=2, 4% at z=3, and 12% at z=4. Previous measurements of tau_eff at z>~3 based on directly fitting the quasar continua in the Lya forest, which have generally neglected this effect, are therefore likely biased low. We provide estimates of the level of absorption arising from metals in the Lya forest based on both direct and statistical metal removal results in the literature, finding that this contribution is ~6-9% at z=3 and decreases monotonically with redshift. The high precision of our measurement, attaining 3% in redshift bins of width Delta z=0.2 around z=3, indicates significant departures from the best-fit power-law redshift evolution (tau_eff=0.0018(1+z)^3.92, when metals are left in), particularly near z=3.2. The observed downward departure is statistically consistent with a similar feature detected in a precision statistical measurement using Sloan Digital Sky Survey spectra by Bernardi and coworkers using an independent approach.

We investigate the implications of our measurement of the Lyman-alpha forest opacity at redshifts 2<z<4.2 from a sample of 86 high-resolution quasar spectra for the evolution of the cosmic ultraviolet luminosity density and its sources. The derived hydrogen photoionization rate is remarkably flat over this redshift range, implying an increasing comoving ionizing emissivity with redshift. Because the quasar luminosity function is strongly peaked near z~2, star-forming galaxies likely dominate the ionizing emissivity at z>~3. Our measurement argues against a star formation rate density declining beyond z~3, in contrast with existing state-of-the-art determinations of the cosmic star formation history from direct galaxy counts. Stellar emission from galaxies therefore likely reionized the Universe.

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

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

The universe is permeated by a network of filaments, sheets, and knots collectively forming a "cosmic web.'' The discovery of the cosmic web, especially through its signature of absorption of light from distant sources by neutral hydrogen in the intergalactic medium, exemplifies the interplay between theory and experiment that drives science, and is one of the great examples in which numerical simulations have played a key and decisive role. We recount the milestones in our understanding of cosmic structure, summarize its impact on astronomy, cosmology, and physics, and look ahead by outlining the challenges faced as we prepare to probe the cosmic web at new wavelengths.

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

We directly measure the evolution of the intergalactic Lyman-alpha effective optical depth, tau_eff, over the redshift range 2<z<4.2 from a sample of 86 high-resolution, high-signal-to-noise quasar spectra obtained with Keck/ESI, Keck/HIRES, and Magellan/MIKE. We find that our estimates of the quasar continuum levels in the Ly-alpha forest obtained by spline fitting are systematically biased low, but that this bias can be accounted for using mock spectra. The mean fractional error <Delta C/C_true> is <1% at z=2, 4% at z=3, and 12% at z=4. We provide estimates of the level of absorption arising from metals in the Ly-alpha forest based on both direct and statistical metal removal results in the literature, finding that this contribution is ~6-9% at z=3 and decreases monotonically with redshift. The high precision of our measurement indicates significant departures from the best-fit power-law redshift evolution, particularly near z=3.2.

We show that the Lyman-alpha optical depth statistics in the proximity regions of quasar spectra depend on the mass of the dark matter halos hosting the quasars. This is owing to both the overdensity around the quasars and the associated infall of gas toward them. For a fiducial quasar host halo mass of (3.0+/-1.6) h^-1 x 10^12 Msun, as inferred by Croom et al. from clustering in the 2dF QSO Redshift Survey, we show that estimates of the ionizing background (Gamma^bkg) from proximity effect measurements could be biased high by a factor of ~2.5 at z=3 owing to neglecting these effects alone. The clustering of galaxies and other active galactic nuclei around the proximity effect quasars enhances the local background, but is not expected to skew measurements by more than a few percent. Assuming the measurements of Gamma^bkg based on the mean flux decrement in the Ly-alpha forest to be free of bias, we demonstrate how the proximity effect analysis can be inverted to measure the mass of the dark matter halos hosting quasars. In ideal conditions, such a measurement could be made with a precision comparable to the best clustering constraints to date from a modest sample of only about 100 spectra. We discuss observational difficulties, including continuum flux estimation, quasar systematic redshift determination, and quasar variability, which make accurate proximity effect measurements challenging in practice. These are also likely to contribute to the discrepancies between existing proximity effect and flux decrement measurements of Gamma^bkg.

We compare observed clustering of quasars and galaxies as a function of redshift, mass, luminosity, & color/morphology, to constrain models of quasar fueling and spheroid-BH co-evolution. High redshift quasars are shown to be drawn from progenitors of local early-type galaxies, with the characteristic quasar luminosity L* reflecting a characteristic mass of 'active' BH/host populations at each redshift. Evolving observed high-z quasar clustering to z=0 predicts a trend of clustering in 'quasar remnants' as a function of stellar mass identical to that observed for early-types. However, quasar clustering does not simply reflect observed early (or late)-type populations; at each redshift, quasars cluster as an 'intermediate' population. Comparing with the age of elliptical stellar populations reveals that this 'intermediate' population represents those ellipticals undergoing or terminating their final significant star formation at each epoch. Assuming that quasar triggering is associated with the formation/termination epoch of ellipticals predicts quasar clustering at all observed redshifts without any model dependence or assumptions about quasar light curves, lifetimes, or accretion rates. This is not true for disks or quasar halos: i.e. quasars do not generically trace star formation, disks, or halo assembly. Quasar clustering at all z is consistent with a constant halo mass ~4x10^{12} M_sun, similar to local 'group scales.' The observations support a scenario in which major mergers trigger quasar activity and dominate bright, high-z quasar populations. We show that measurements of quasar clustering versus luminosity at z~1 can be used to constrain different lower-luminosity AGN fueling mechanisms, and that high-z clustering is sensitive to whether or not systems 'shut down' growth at z>3.

We calculate the impact of the intergalactic medium (IGM) on the observed Lyman alpha lines (hereafter Lya) emitted by galaxies in an ionised IGM at z>4. Our model accounts for gas clumping in the IGM and for the fact that high-redshift galaxies reside in overdense regions, which causes the velocity field of the IGM to depart from the Hubble flow. The observed shape of the Lya line varies widely, with dependence on the intrinsic width and systemic velocity of the line, a galaxies star formation rate and the local extra-galactic UV-background. For large star formation rates and levels of the UV-background, absorption in the IGM does not result in a Lya line that is asymmetric as is common among known high-redshift Lya emitters. For models in which the lines do show the observed strong asymmetry, the IGM typically transmits only 10-30% of the Lya flux. The increase in the ionising background that accompanied the completion of reionisation barely increased the IGM transmission, which suggests that LAEs of comparable luminosity should not appear to be significantly dimmer prior to overlap. In this light, we briefly discuss the potential of Lya emitters as a probe into the epoch of reionisation.

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

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.

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

We explore the impact of spatial fluctuations in the intergalactic medium temperature on the Lyman-alpha forest flux power spectrum near z ~ 3. We develop a semianalytic model to examine temperature fluctuations resulting from inhomogeneous HI and incomplete HeII reionizations. Detection of these fluctuations might provide insight into the reionization histories of hydrogen and helium. Furthermore, these fluctuations, neglected in previous analyses, could bias constraints on cosmological parameters from the Lyman-alpha forest. We find that the temperature fluctuations resulting from inhomogeneous HI reionization are likely to be very small, with an rms amplitude of < 5%, $\sigma_{T_0}/<T_0> < 0.05$. More important are the temperature fluctuations that arise from incomplete HeII reionization, which might plausibly be as large as 50%, $\sigma_{T_0}/<T_0> ~ 0.5$. In practice, however, these temperature fluctuations have only a small effect on flux power spectrum predictions. The smallness of the effect is possibly due to density fluctuations dominating over temperature fluctuations on the scales probed by current measurements. On the largest scales currently probed, k ~ 0.001 s/km (~0.1 h/Mpc), the effect on the flux power spectrum may be as large as ~10% in extreme models. The effect is larger on small scales, up to ~20% at k = 0.1 s/km, due to thermal broadening. Our results suggest that the omission of temperature fluctuations effects from previous analyses does not significantly bias constraints on cosmological parameters.

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.

We investigate the luminosity dependence of quasar clustering, inspired by numerical simulations of galaxy mergers that incorporate black hole growth. These simulations have motivated a new interpretation of the quasar luminosity function. In this picture, the bright end of the quasar luminosity function consists of quasars radiating nearly at their peak luminosities, while the faint end consists mainly of very similar sources, but at dimmer phases in their evolution. We combine this model with the statistics of dark matter halos that host quasar activity. We find that, since bright and faint quasars are mostly similar sources seen in different evolutionary stages, a broad range in quasar luminosities corresponds to only a narrow range in the masses of quasar host halos. On average, bright and faint quasars reside in similar host halos. Consequently, we argue that quasar clustering should depend only weakly on luminosity. This prediction is in qualitative agreement with recent measurements of the luminosity dependence of the quasar correlation function (Croom et al. 2005) and the galaxy-quasar cross-correlation function (Adelberger & Steidel 2005). Future precision clustering measurements from SDSS and 2dF, spanning a large range in luminosity, should provide a strong test of our model.

The analysis of the Lyman-alpha forest of absorption lines in quasar spectra has emerged as a potentially powerful technique to constrain the linear matter power spectrum. In most previous work, the amplitude of the ionizing background was fixed by calibrating simulations to match the observed mean transmitted flux in the Lyman-alpha forest. This procedure is undesirable in principle as it requires the estimation of the unabsorbed quasar continuum level, a difficult undertaking subject to various sources of systematic error and bias. We suggest an alternative approach based on measuring the one-point probability distribution function (pdf) of the fluctuations in the flux about the mean, relative to the mean, i.e. the pdf of delta_f = (f - <f>)/<f>. This statistic, while sensitive to the amplitude of the ionizing background, has the virtue that its measurement does not require an estimate of the unabsorbed continuum level. We present a measurement of the pdf of delta_f from seven Keck HIRES spectra, spanning a redshift range of z = 2.2 - 4.4. To illustrate that our method is useful, we compare our measurements of the pdf of delta_f, and measurements of the flux power spectrum from Croft et al. (2002) at z = 2.72, with cosmological simulations. From this comparison, we obtain constraints on the mean transmission in the Lyman-alpha forest, the slope of the temperature-density relation, as well as the amplitude and slope of the mass power spectrum. Our methodology will be useful for obtaining more precise constraints with larger data samples from the Sloan Digital Sky Survey (SDSS).

Absorption between the rest-frame wavelengths of 973 and 1026 Angstroms in quasar spectra arises from two sources (apart from occasional metals): one is due to Lyman-alpha (Lya) absorption by materials at a low redshift, and the other from Lyman-beta (Lyb) at a higher redshift. These two sources of absorption are to a good approximation uncorrelated because of their wide physical separation. Therefore, the two-point correlation of absorption in this region of quasar spectra neatly factorizes into two pieces: the Lyb correlation at high z, and the Lya correlation at low z. The latter can be independently measured from quasar spectra at lower redshifts using current techniques. A simple division then offers a way to statistically separate out the Lyb two-point correlation from the Lya correlation. Several applications of this technique are discussed. First, since the Lyb absorption cross-section is lower than Lya by about a factor of 5, the Lyb forest is a better probe of the intergalactic medium (IGM) at higher redshifts where Lya absorption is often saturated. Second, for the same reason, the Lyb forest allows a better measurement of the equation of state of the IGM at higher overdensities, yielding stronger constraints on its slope when used in conjunction with the Lya forest. Third, models of the Lya forest based on gravitational instability make unique predictions for the Lyb forest, which can be tested against observations. We briefly point out that feedback processes that affect higher density regions but leave low density structure intact may be better constrained by the Lyb forest.

A comparison between the line of sight power spectrum of absorption in the Lyman-alpha forest and the cross power spectrum between the absorption in neighboring lines of sight offers an evolution-free means to constrain the cosmological constant, or dark energy. Using cosmological simulations, we consider a maximum likelihood method to obtain constraints from this comparison. In our method, measurements of the auto and cross spectra from observations are compared with those from a multi-parameter grid of simulated models of the intergalactic medium. We then marginalize over nuisance parameters to obtain constraints on the cosmological constant. Redshift space distortions due to peculiar velocities and thermal broadening, a potential difficulty in applying this test, are explicitly modeled in our simulations. To illustrate our method, we measure the cross spectrum from a new sample of five close quasar pairs, with separations of 0.5 to 3 arcmin. We attempt to obtain a constraint on Omega_Lambda, but find only weak constraints. An Einstein-de-Sitter cosmology is, however, disfavored by the data at a ~ 2 sigma confidence level. We consider the power of future observations, paying particular attention to the effects of spectral resolution and shot-noise. We find that ~ 50 moderate resolution, FWHM ~ 150 km/s, close separation (~ 30-120 arcsec) pairs should allow a (2 sigma) constraint on Omega_Lambda at the level of 15%, if other modeling parameters are determined through other means. We find that there is a sizeable gain from observing very close, ~ 30 arcsec, separation pairs provided they are observed with high spectral resolution. A sample of ~ 10 such pairs gives similar constraints to the 50 moderate resolution pairs mentioned above.

Recent observations of high redshift quasar spectra reveal long gaps with little flux. A small or no detectable flux does not by itself imply the intergalactic medium (IGM) is neutral. Inferring the average neutral fraction from the observed absorption requires assumptions about clustering of the IGM, which the gravitational instability model supplies. Our most stringent constraint on the neutral fraction at z ~ 6 is derived from the mean Lyman-beta transmission measured from the z=6.28 SDSS quasar of Becker et al. -- the neutral hydrogen fraction at mean density has to be larger than 4.7 times 10^{-4}. This is substantially higher than the neutral fraction of ~ 3-5 times 10^{-5} at z = 4.5 - 5.7, suggesting that dramatic changes take place around or just before z ~ 6, even though current constraints are still consistent with a fairly ionized IGM at z ~ 6. An interesting alternative method to constrain the neutral fraction is to consider the probability of having many consecutive pixels with little flux, which is small unless the neutral fraction is high. This constraint is slightly weaker than the one obtained from the mean transmission. We show that while the derived neutral fraction at a given redshift is sensitive to the power spectrum normalization, the size of the jump around z ~ 6 is not. We caution that systematic uncertainties include spatial fluctuations in the ionizing background, and the continuum placement. Tests are proposed. In particular, the sightline to sightline dispersion in mean transmission might provide a useful diagnostic. We express the dispersion in terms of the transmission power spectrum, and develop a method to calculate the dispersion for spectra that are longer than the typical simulation box.