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

Previous studies have noted difficulties in modeling the highest opacities of the $z>5.5$ Ly$\alpha$ forest, epitomized by the extreme Ly$\alpha$ trough observed towards quasar ULAS J0148+0600. One possibility is that the most opaque regions at these redshifts contain significant amounts of neutral hydrogen. This explanation, which abandons the common assumption that reionization ended before $z=6$, also reconciles evidence from independent observations of a significantly neutral Universe at $z=7.5$. Here we explore a model in which the neutral fraction is still $\approx 10\%$ at $z=5.5$. We confirm that this model can account for the observed scatter in Ly$\alpha$ forest opacities, as well as the observed Ly$\beta$ transmission in the J0148 trough. We contrast the model with a competing "earlier" reionization scenario characterized by a short mean free path and large fluctuations in the post-reionization ionizing background. We consider Ly$\alpha$ and Ly$\beta$ effective optical depths, their correlations, trough size distributions, dark pixel fractions, the IGM thermal history, and spatial distributions of Lyman-$\alpha$ emitters around forest sight lines. We find that the models are broadly similar in almost all of these statistics, suggesting that it may be difficult to distinguish between them definitively. We argue that improved constraints on the mean free path and the thermal history at $z>5$ could go a long way towards diagnosing the origin of the $z>5.5$ opacity fluctuations.

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

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

We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of <5% throughout the bulk of the epoch of reionization if: i) galaxies form stars down to the atomic cooling limit before reionization and a photosuppression mass of log(M_h/Msol)~9 during/after reionization (-13<M_UV,lim<-11); ii) galaxies become more efficient producers of ionizing photons at higher redshifts and fainter magnitudes, and iii) there is a significant, but sub-dominant, contribution by AGN at z < 7. In this model the faintest galaxies (M_UV>-15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+\- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%.

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

Absorption signatures in the spectra of QSOs are one of our most powerful tools for studying galactic and intergalactic environments at high redshifts. With the discovery of QSOs out to z > 7, QSO absorption lines are now tracing the end stages of reionization on multiple fronts using the hydrogen Lyman-$\alpha$ forest and heavy element absorbers. Next-generation QSO absorption line studies with large optical/IR telescopes will reveal in detail how the first galaxies emerged form the cosmic web, transformed their circum- and inter-galactic environments, and completed the last major phase transition of the Universe. These efforts will complement other upcoming studies of reionization, such as those with JWST, ALMA, and redshifted 21cm experiments.

The 21-cm signal from the epoch of reionization is non-Gaussian. Current radio telescopes are focused on detecting the 21-cm power spectrum, but in the future the Square Kilometre Array is anticipated to provide a first measurement of the bispectrum. Previous studies have shown that the position-dependent power spectrum is a simple and efficient way to probe the squeezed-limit bispectrum. In this approach, the survey is divided into subvolumes and the correlation between the local power spectrum and the corresponding mean density of the subvolume is computed. This correlation is equivalent to an integral of the bispectrum in the squeezed limit, but is much simpler to implement than the usual bispectrum estimators. It also has a clear physical interpretation: it describes how the small-scale power spectrum of tracers such as galaxies and the 21-cm signal respond to a large-scale environment. Reionization naturally couples large and small scales as ionizing radiation produced by galactic sources can travel up to tens of Megaparsecs through the intergalactic medium during this process. Here we apply the position-dependent power spectrum approach to fluctuations in the 21-cm background from reionization. We show that this statistic has a distinctive evolution in time that can be understood with a simple analytic model. We also show that the statistic can easily distinguish between simple "inside-out" and "outside-in" models of reionization. The position-dependent power spectrum is thus a promising method to validate the reionization signal and to extract higher-order information on this process.

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

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

Photoheating associated with reionization suppressed star formation in low-mass galaxies. Reionization was inhomogeneous, however, affecting different regions at different times. To establish the causal connection between reionization and suppression, we must take this local variation into account. We analyze the results of CoDa (`Cosmic Dawn') I, the first fully-coupled radiation-hydrodynamical simulation of reionization and galaxy formation in the Local Universe, in a volume large enough to model reionization globally but with enough resolving power to follow all atomic-cooling galactic halos in that volume. For every halo identified at a given time, we find the redshift at which the surrounding IGM reionized, along with its instantaneous star formation rate (`SFR') and baryonic gas-to-dark matter ratio ($M_\text{gas}/M_\text{DM}$). The average SFR per halo with $M < 10^9 \text{ M}_\odot$ was steady in regions not yet reionized, but declined sharply following local reionization. For $M > 10^{10} \text{ M}_\odot$, this SFR continued through local reionization, increasing with time, instead. For $10^9 < M < 10^{10} \text{ M}_\odot$, the SFR generally increased modestly through reionization, followed by a modest decline. In general, halo SFRs were higher for regions that reionized earlier. A similar pattern was found for $M_\text{gas}/M_\text{DM}$, which declined sharply following local reionization for $M < 10^9 \text{ M}_\odot$. Local reionization time correlates with local matter overdensity, which determines the local rates of structure formation and ionizing photon consumption. The earliest patches to develop structure and reionize ultimately produced more stars than they needed to finish and maintain their own reionization, exporting their `surplus' starlight to help reionize regions that developed structure later.

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

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

Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while resolving millions of galaxy sources above ~10^8 Msun , combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the local universe today. Cosmic Dawn ("CoDa") is the first such fully-coupled, radiation-hydrodynamics simulation of reionization of the local universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc)^3 with 4096^3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10^9 Msun , decreasing the abun- dance of faint galaxies around MAB_1600 = [-10,-12]. For most of reionization, star formation was dominated by haloes between 10^10 - 10^11 Msun , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. (Abridged)

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

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

The linear perturbation theory of inhomogeneous reionization (LPTR) has been developed as an analytical tool for predicting the global ionized fraction and large-scale power spectrum of ionized density fluctuations during reionization. In the original formulation of the LPTR, the ionization balance and radiative transfer equations are linearized and solved in Fourier space. However, the LPTR's approximation to the full solution of the radiative transfer equation is not straightforward to interpret, since the latter is most intuitively conceptualized in position space. To bridge the gap between the LPTR and the language of numerical radiative transfer, we present a new, equivalent, position-space formulation of the LPTR that clarifies the approximations it makes and facilitates its interpretation. We offer a comparison between the LPTR and the excursion-set model of reionization (ESMR), and demonstrate the built-in capability of the LPTR to explore a wide range of reionization scenarios, and to go beyond the ESMR in exploring scenarios involving X-rays.

Cluster gravitational lensing surveys like the Hubble Space Telescope Frontier Fields survey will detect distant galaxies 10-50 times fainter than any yet discovered. Using these surveys to measure the luminosity function of such faint, distant galaxies, however, requires that magnification maps built from the constraints of strongly-lensed images be accurate. For models that assume the cluster and nearby (correlated) structures are the only significant sources of lensing, a potential source of error in these maps comes from the fact that light rays also suffer weak deflections by uncorrelated large-scale structure along the line-of-sight, i.e. cosmic weak lensing (CWL). To demonstrate the magnitude of this effect, we calculate the magnification change which results when the same cluster-lens is placed along different lines of sight. Using a simple density profile for a cluster-lens at z~0.3-0.5 and the power spectrum of the matter density fluctuations responsible for CWL, we show that the typical magnifications of ~5(10) of sources at z=6-10 can differ by ~10-20(20-30)% from one line-of-sight to another. However, these fluctuations rise to greater than order unity near critical curves, indicating that CWL tends to make its greatest contribution to the most magnified images. We conclude that the neglect of CWL in determining the intrinsic luminosities of highly-magnified galaxies may introduce errors significant enough to warrant further effort to include this contribution in cluster-lens modeling. We suggest that methods of modeling CWL in galaxy-strong-lensing systems should be generalized to cluster-lensing systems.

In this paper we investigate how the halo mass function evolves with redshift, based on a suite of very large (with N_p = 3072^3 - 6000^3 particles) cosmological N-body simulations. Our halo catalogue data spans a redshift range of z = 0-30, allowing us to probe the mass function from the dark ages to the present. We utilise both the Friends-of-Friends (FOF) and Spherical Overdensity (SO) halofinding methods to directly compare the mass function derived using these commonly used halo definitions. The mass function from SO haloes exhibits a clear evolution with redshift, especially during the recent era of dark energy dominance (z < 1). We provide a redshift-parameterised fit for the SO mass function valid for the entire redshift range to within ~20% as well as a scheme to calculate the mass function for haloes with arbitrary overdensities. The FOF mass function displays a weaker evolution with redshift. We provide a `universal' fit for the FOF mass function, fitted to data across the entire redshift range simultaneously, and observe redshift evolution in our data versus this fit. The relative evolution of the mass functions derived via the two methods is compared and we find that the mass functions most closely match at z=0. The disparity at z=0 between the FOF and SO mass functions resides in their high mass tails where the collapsed fraction of mass in SO haloes is ~80% of that in FOF haloes. This difference grows with redshift so that, by z>20, the SO algorithm finds a ~50-80% lower collapsed fraction in high mass haloes than does the FOF algorithm, due in part to the significant over-linking effects known to affect the FOF method.

Measuring the small primordial nonGaussianity (PNG) predicted by cosmic inflation theories may help diagnose them. The detectability of PNG by its imprint on the 21cm power spectrum from the epoch of reionization is reassessed here in terms of $f_{NL}$, the local nonlinearity parameter. We find that an optimum, multi-frequency observation by SKA can achieve $\Delta f_{NL} \sim 3$ (comparable to recent Planck CMB limits), while a cosmic-variance-limited array of this size like Omniscope can even detect $\Delta f_{NL} \sim 0.2$. This substantially revises the methods and results of previous work.

(ABRIDGED)The rise of cosmic structure depends upon the statistical distribution of initial density fluctuations generated by inflation. While the simplest models predict an almost perfectly Gaussian distribution, more-general models predict a level of primordial non-Gaussianity (PNG) that observations might yet be sensitive enough to detect. Recent Planck Collaboration measurements of the CMB temperature anisotropy bispectrum significantly tighten the observational limits, but they are still far from the PNG level predicted by the simplest models of inflation. Probing levels below CMB sensitivities will require other methods, such as searching for the statistical imprint of PNG on galactic halo clustering. During the epoch of reionization (EoR), the first stars and galaxies released radiation into the intergalactic medium (IGM) that created ionized patches whose large-scale geometry and evolution reflected the underlying abundance and large-scale clustering of the star-forming galaxies. This statistical connection between ionized patches in the IGM and galactic halos suggests that observing reionization may be another way to constrain PNG. We employ the linear perturbation theory of reionization and semi-analytic models based on the excursion-set formalism to model the effects of PNG on the EoR. We quantify the effects of PNG on the large-scale structure of reionization by deriving the ionized density bias, i.e. ratio of ionized atomic to total matter overdensities in Fourier space, at small wavenumber. Just as previous studies found that PNG creates a scale-dependent signature in the halo bias, so, too, we find a scale-dependent signature in the ionized density bias. Our results, which differ significantly from previous attempts in the literature to characterize this PNG signature, will be applied elsewhere to predict its observable consequences, e.g. in the cosmic 21cm background.

Characterizing the level of primordial non-Gaussianity (PNG) in the initial conditions for structure formation is one of the most promising ways to test inflation and differentiate among different scenarios. The scale-dependent imprint of PNG on the large-scale clustering of galaxies and quasars has already been used to place significant constraints on the level of PNG in our observed Universe. Such measurements depend upon an accurate and robust theory of how PNG affects the bias of galactic halos relative to the underlying matter density field. We improve upon previous work by employing a more general analytical method - the path-integral extension of the excursion set formalism - which is able to account for the non-Markovianity caused by PNG in the random-walk model used to identify halos in the initial density field. This non-Markovianity encodes information about environmental effects on halo formation which have so far not been taken into account in analytical bias calculations. We compute both scale-dependent and -independent corrections to the halo bias, along the way presenting an expression for the conditional collapsed fraction for the first time, and a new expression for the conditional halo mass function. To leading order in our perturbative calculation, we recover the halo bias results of Desjacques et. al. (2011), including the new scale-dependent correction reported there. However, we show that the non-Markovian dynamics from PNG can lead to marked differences in halo bias when next-to-leading order terms are included. We quantify these differences here. [abridged]

In a recently published article, we quantified the impact of primordial non-Gaussianity on the probability of giant-arc formation. In that work, we focused on the local form of non-Gaussianity and found that it can have only a modest effect given the most recent constraints from Cosmic Microwave Background (CMB) measurements. Here, we present new calculations using a parameterization of scale-dependent non-Gaussianity in which the primordial bispectrum has the equilateral shape and the effective f_NL parameter depends on scale. We find that non-Gaussianity of this type can yield a larger effect on the giant-arc abundance compared to the local form due to both the scale dependence and the relatively weaker constraints on the equilateral shape from CMB measurements. In contrast to the maximum ~40% effect (within the latest CMB constraints) previously found for the local form, we find that the predicted giant-arc abundance for the scale-dependent equilateral form can differ by a factor of a few with respect to the Gaussian case.

We explore the use of strong lensing by galaxy clusters to constrain the dark energy equation of state and its possible time variation. The cores of massive clusters often contain several multiply imaged systems of background galaxies at different redshifts. The locations of lensed images can be used to constrain cosmological parameters due to their dependence on the ratio of angular diameter distances. We employ Monte-Carlo simulations of cluster lenses, including the contribution from substructures, to assess the feasibility of this potentially powerful technique. At the present, parametric lens models use well motivated scaling relations between mass and light to incorporate cluster member galaxies, and do not explicitly model line-of-sight structure. Here, we quantify modeling errors due to scatter in the cluster galaxy scaling relations and un-modeled line-of-sight halos. These errors are of the order of a few arcseconds on average for clusters located at typical redshifts (z ~ 0.2 - 0.3). Using Bayesian Markov Chain Monte-Carlo techniques, we show that the inclusion of these modeling errors is critical to deriving unbiased constraints on dark energy. However, when the uncertainties are properly quantified, we show that constraints competitive with other methods may be obtained by combining results from a sample of just 10 simulated clusters with 20 families each. Cosmography with a set of well studied cluster lenses may provide a powerful complementary probe of the dark energy equation of state. Our simulations provide a convenient method of quantifying modeling errors and assessing future strong lensing survey strategies.

We present results of a study of the virial state of high redshift dark matter haloes in an N-body simulation. We find that the majority of collapsed, bound haloes are not virialized at any redshift slice in our study ($z=15-6$) and have excess kinetic energy. At these redshifts, merging is still rampant and the haloes cannot strictly be treated as isolated systems. To assess if this excess kinetic energy arises from the environment, we include the surface pressure term in the virial equation explicitly and relax the assumption that the density at the halo boundary is zero. Upon inclusion of the surface term, we find that the haloes are much closer to virialization, however, they still have some excess kinetic energy. We report trends of the virial ratio including the extra surface term with three key halo properties: spin, environment, and concentration. We find that haloes with closer neighbors depart more from virialization, and that haloes with larger spin parameters do as well. We conclude that except at the lowest masses ($M < 10^6 \Msun$), dark matter haloes at high redshift are not fully virialized. This finding has interesting implications for galaxy formation at these high redshifts, as the excess kinetic energy will impact the subsequent collapse of baryons and the formation of the first disks and/or baryonic structures.

For over a decade, it has been debated whether the concordance LCDM model is consistent with the observed abundance of giant arcs in clusters. While previous theoretical studies have focused on properties of the lens and source populations, as well as cosmological effects such as dark energy, the impact of initial conditions on the giant-arc abundance is relatively unexplored. Here, we quantify the impact of non-Gaussian initial conditions with the local bispectrum shape on the predicted frequency of giant arcs. Using a path-integral formulation of the excursion set formalism, we extend a semi-analytic model for calculating halo concentrations to the case of primordial non-Gaussianity, which may be useful for applications outside of this work. We find that massive halos tend to collapse earlier in models with positive f_NL, relative to the Gaussian case, leading to enhanced concentration parameters. The converse is true for f_NL < 0. In addition to these effects, which change the lensing cross sections, non-Gaussianity also modifies the abundance of supercritical clusters available for lensing. These combined effects work together to either enhance (f_NL > 0) or suppress (f_NL < 0) the probability of giant-arc formation. Using the best value and 95% confidence levels currently available from the Wilkinson Microwave Anisotropy Probe, we find that the giant-arc optical depth for sources at z_s~2 is enhanced by ~20% and ~45% for f_NL = 32 and 74 respectively. In contrast, we calculate a suppression of ~5% for f_NL = -10. These differences translate to similar relative changes in the predicted all-sky number of giant arcs.

Current efforts in observational cosmology are focused on characterizing the mass-energy content of the Universe. We present results from a geometric test based on strong lensing in galaxy clusters. Based on Hubble Space Telescope images and extensive ground-based spectroscopic follow-up of the massive galaxy cluster Abell 1689, we used a parametric model to simultaneously constrain the cluster mass distribution and dark energy equation of state. Combining our cosmological constraints with those from X-ray clusters and the Wilkinson Microwave Anisotropy Probe 5-year data gives {\Omega}m = 0.25 +/- 0.05 and wx = -0.97 +/- 0.07 which are consistent with results from other methods. Inclusion of our method with all other techniques available brings down the current 2{\sigma} contours on the dark energy equation of state parameter wx by about 30%.

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

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