Cosmic Dawn ("CoDa") II yields the first statistically-meaningful determination of the relative contribution to reionization by galaxies of different halo mass, from a fully-coupled radiation-hydrodynamics simulation of the epoch of reionization large enough ($\sim$ 100 Mpc) to model global reionization while resolving the formation of all galactic halos above $\sim 10^8 M_\odot$. Cell transmission inside high-mass haloes is bi-modal -- ionized cells are transparent, while neutral cells absorb the photons their stars produce - and the halo escape fraction $f_{esc}$ reflects the balance of star formation rate ("SFR") between these modes. The latter is increasingly prevalent at higher halo mass, driving down $f_{esc}$ (we provide analytical fits to our results), whereas halo escape luminosity, proportional to $f_{esc} \times$SFR, increases with mass. Haloes with dark matter masses within $6.10^{8} M_\odot < M_h < 3.10^{10} M_\odot$ produce $\sim 80$% of the escaping photons at z=7, when the Universe is 50% ionized, making them the main drivers of cosmic reionization. Less massive haloes, though more numerous, have low SFRs and contribute less than 10% of the photon budget then, despite their high $f_{esc}$. High mass haloes are too few and too opaque, contributing $<10$% despite their high SFRs. The dominant mass range is lower (higher) at higher (lower) redshift, as mass function and reionization advance together (e.g. at z$=8.5$, x$_{\rm HI}=0.9$, $M_h < 5.10^9 M_\odot$ haloes contributed $\sim$80%). Galaxies with UV magnitudes $M_{AB1600}$ between $-12$ and $-19$ dominated reionization between z$=6$ and 8.

The decline in abundance of Lyman-$\alpha$ (Ly$\alpha$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$\alpha$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$. Our results mainly confirm previous studies, i.e., we find a declining Ly$\alpha$ transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Ly$\alpha$ peaks at high redshift, we also analyze the IGM transmission on the blue side, which shows a rapid decline at $z\gtrsim 6$ of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\,\mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 5\%$ of sightlines towards $M_{\mathrm{1600AB}}\sim -21$ galaxies at $z\sim 7$ are transparent enough to allow a transmission of a blue Ly$\alpha$ peak. We discuss our results in the context of the interpretation of observations.

Cosmic Dawn II (CoDa II) is a new, fully-coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift $z < 6$. With $4096^3$ particles and cells in a 94 Mpc box, it is large enough to model global reionization and its feedback on galaxy formation while resolving all haloes above $10^8$ M$_{\odot}$. Using the same hybrid CPU-GPU code RAMSES-CUDATON as CoDa I in Ocvirk et al. (2016), CoDa II modified and re-calibrated the subgrid star-formation algorithm, making reionization end earlier, at $z \gtrsim 6$, thereby better matching the observations of intergalactic Lyman-alpha opacity from quasar spectra and electron-scattering optical depth from cosmic microwave background fluctuations. CoDa II predicts a UV continuum luminosity function in good agreement with observations of high-z galaxies, especially at $z = 6$. As in CoDa I, reionization feedback suppresses star formation in haloes below $\sim 2 \times 10^9$ M$_{\odot}$, though suppression here is less severe, a possible consequence of modifying the star-formation algorithm. Suppression is environment-dependent, occurring earlier (later) in overdense (underdense) regions, in response to their local reionization times. Using a constrained realization of $\Lambda$CDM constructed from galaxy survey data to reproduce the large-scale structure and major objects of the present-day Local Universe, CoDa II serves to model both global and local reionization. In CoDa II, the Milky Way and M31 appear as individual islands of reionization, i.e. they were not reionized by the progenitor of the Virgo cluster, nor by nearby groups, nor by each other.

Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation.

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.

Today's galaxies experienced cosmic reionization at different times in different locations. For the first time, reionization ($50\%$ ionized) redshifts, $z_R$, at the location of their progenitors are derived from new, fully-coupled radiation-hydrodynamics simulation of galaxy formation and reionization at $z > 6$, matched to N-body simulation to z = 0. Constrained initial conditions were chosen to form the well-known structures of the local universe, including the Local Group and Virgo, in a (91 Mpc)$^3$ volume large enough to model both global and local reionization. Reionization simulation CoDa I-AMR, by CPU-GPU code EMMA, used (2048)$^3$ particles and (2048)$^3$ initial cells, adaptively-refined, while N-body simulation CoDa I-DM2048, by Gadget2, used (2048)$^3$ particles, to find reionization times for all galaxies at z = 0 with masses $M(z=0)\ge 10^8 M_\odot$. Galaxies with $M(z=0) \gtrsim 10^{11} M_\odot$ reionized earlier than the universe as a whole, by up to $\sim$ 500 Myrs, with significant scatter. For Milky-Way-like galaxies, $z_R$ ranged from 8 to 15. Galaxies with $M(z=0) \lesssim 10^{11} M_\odot$ typically reionized as late or later than globally-averaged $50\%$ reionization at $\langle z_R\rangle =7.8$, in neighborhoods where reionization was completed by external radiation. The spread of reionization times within galaxies was sometimes as large as the galaxy-to-galaxy scatter. The Milky Way and M31 reionized earlier than global reionization but later than typical for their mass, neither dominated by external radiation. Their most massive progenitors at $z>6$ had $z_R$ = 9.8 (MW) and 11 (M31), while their total masses had $z_R$ = 8.2 (both).

We consider an alternative cold dark matter candidate, ultralight bosons ($m>10^{-22}$eV) described by a complex scalar field (SFDM) with global U(1) symmetry, with comoving particle number density conserved after particle production during standard reheating. We allow for repulsive self-interaction. In a Lambda-SFDM universe, SFDM starts relativistic, evolving from stiff (w=1) to radiation-like (w=1/3), becoming nonrelativistic (w=0) at late times. Thus, a stiff-SFDM-dominated era precedes the familiar radiation-dominated era. SFDM particle mass $m$ and quartic self-interaction strength \lambda, are therefore constrained by cosmological observables, N_{eff}, the effective number of neutrino species during BBN, and z_{eq}, the matter-radiation equality redshift. Since the stochastic gravitational wave background (SGWB) from inflation is amplified during the stiff-SFDM-dominated era, it can also contribute a radiationlike component large enough to affect these observables. Remarkably, this amplification makes this SGWB detectable by current GW experiments, e.g., aLIGO/Virgo and LISA, for Lambda-SFDM models satisfying cosmological constraints, for a range of reheat temperatures T_{re} and currently allowed values of tensor-to-scalar ratio $r$. For given r and $\lambda/(mc^2)^2$, the marginally-allowed Lambda-SFDM model for each T_{re} has the smallest m that satisfies cosmological constraints. For example, for marginally-allowed models with r=0.01 and $\lambda/(mc^2)^2=10^{-18}$eV$^{-1}$cm$^3$, null detection by the aLIGO O1 run excludes 8.75*10^3<T_{re} (GeV)<1.7*10^5 at 95% confidence, demonstrating that GW experiments already place a new kind of cosmological constraint on SFDM. A wider parameter range should be accessible to aLIGO/Virgo O5, with potential to detect this signature of Lambda-SFDM. For this same illustrative family, 3-sigma detection is predicted for 600<T_{re} (GeV)<10^7.

Density inhomogeneity in the intergalactic medium (IGM) can boost the recombination rate of ionized gas substantially, affecting the growth of HII regions during reionization. Previous attempts to quantify this effect typically failed to resolve down to the Jeans scale in the pre-ionization IGM, which is important in establishing this effect, along with the hydrodynamical back-reaction of reionization on it. Towards that end, we perform a set of fully-coupled, radiation-hydrodynamics simulations from cosmological initial conditions, extending the mass resolution of previous work to the scale of minihalos. Pre-reionization structure is evolved until a redshift $z_i$ at which the ionizing radiation from external sources arrives to sweep an R-type ionization front supersonically across the volume in a few Myr, until it is trapped on the surfaces of minihalos and converted to D-type, after which the minihalo gas is removed by photoevaporative winds. Small-scale density structures during this time lead to a high ($>$10) clumping factor for ionized gas, which hugely boosts the recombination rate until the structures are disrupted by the hydrodynamic feedback after $\sim 10-100~\rm{Myr}$. For incoming stellar radiation with intensity $J_{21}$ in a $200~h^{-1}~\rm{kpc}$ box with the mean density contrast $\bar\delta$, the number of extra recombinations per H atom, on top of what is expected from homogeneously distributed gas, is given by $0.32[J_{21}]^{0.12}[(1+z_i)/11]^{-1.7}[1+\bar\delta]^{2.5}$. In models in which most of the volume is ionized toward the end of reionization, this can add more than one recombination per H atom to the ionizing photon budget to achieve reionization.

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)

We present the first large-scale, full radiative transfer simulations of the reionization of the intergalactic medium in the presence of Lyman-limit systems (LLSs). To illustrate the impact of LLS opacity, possibly missed by previous simulations, we add either a uniform or spatially-varying hydrogen bound-free opacity. This opacity, implemented as the mean free path (mfp) of the ionizing photons, extrapolates the observed, post-reionization redshift dependence into the epoch of reionization. In qualitative agreement with previous studies, we find that at late times the presence of LLSs slows down the ionization fronts, and alters the size distribution of H II regions. We quantitatively characterize the size distribution and morphological evolution of H II regions and examine the effects of the LLSs on the redshifted 21-cm signal from the patchy reionization. The presence of LLSs extends the ionization history by $\Delta z \sim 0.8$. The LLS absorbers significantly impede the late-time growth of the H II regions. The position dependent LLS distribution slows reionization further and additionally limits the late growth of the ionized regions. However, there is no "freeze out" of the H II regions and the largest regions grow to the size of the simulation volume. The 21-cm power spectra show that at large scales the power drops by a factor of 2 for 50% and 75% ionization stages (at $k = 0.1$ $\text{h} \, \text{Mpc}^{-1} $) reflecting the limiting effect of the LLSs on the growth of ionized patches. The statistical observables such as the RMS of the brightness temperature fluctuations and the peak amplitudes of the 21-cm power spectra at large-scales ($k = 0.05 - 0.1$ $\text{h} \, \text{Mpc}^{-1} $) are diminished by the presence of LLS.

Observations of the epoch of reionization give us clues about the nature and evolution of the sources of ionizing photons, or early stars and galaxies. We present a new suite of structure formation and radiative transfer simulations from the PRACE4LOFAR project designed to investigate whether the mechanism of radiative feedback, or the suppression of star formation in ionized regions from UV radiation, can be inferred from these observations. Our source halo mass extends down to $10^8 M_\odot$, with sources in the mass range $10^8$ to $10^9 M_\odot$ expected to be particularly susceptible to feedback from ionizing radiation, and we vary the aggressiveness and nature of this suppression. Not only do we have four distinct source models, we also include two box sizes (67 Mpc and 349 Mpc), each with two grid resolutions. This suite of simulations allows us to investigate the robustness of our results. All of our simulations are broadly consistent with the observed electron-scattering optical depth of the cosmic microwave background and the neutral fraction and photoionization rate of hydrogen at $z\sim6$. In particular, we investigate the redshifted 21-cm emission in anticipation of upcoming radio interferometer observations. We find that the overall shape of the 21-cm signal and various statistics are robust to the exact nature of source suppression, the box size, and the resolution. There are some promising model discriminators in the non-Gaussianity and small-scale power spectrum of the 21-cm signal.

Cosmological transverse momentum fields, whose directions are perpendicular to Fourier wave vectors, induce temperature anisotropies in the cosmic microwave background via the kinetic Sunyaev-Zeldovich (kSZ) effect. The transverse momentum power spectrum contains the four-point function of density and velocity fields, $\langle\delta\delta v v\rangle$. In the post-reionization epoch, nonlinear effects dominate in the power spectrum. We use perturbation theory and cosmological $N$-body simulations to calculate this nonlinearity. We derive the next-to-leading order expression for the power spectrum with a particular emphasis on the connected term that has been ignored in the literature. While the contribution from the connected term on small scales ($k>0.1\,h\,\rm{Mpc}^{-1}$) is subdominant relative to the unconnected term, we find that its contribution to the kSZ power spectrum at $\ell = 3000$ at $z<6$ can be as large as ten percent of the unconnected term, which would reduce the allowed contribution from the reionization epoch ($z>6$) by twenty percent. The power spectrum of transverse momentum on large scales is expected to scale as $k^2$ as a consequence of momentum conservation. We show that both the leading and the next-to-leading order terms satisfy this scaling. In particular, we find that both of the unconnected and connected terms are necessary to reproduce $k^2$.

The cosmic evolution of a dark matter model which behaves relativistically in the early Universe is explored. Dark matter is described as a complex scalar field, whose earliest evolution is characterized by a stiff equation of state ($p \simeq \rho$). In this phase, it is the dominant component in the Universe. We present constraints from Big Bang nucleosynthesis and primordial gravity waves from inflation. Also, we study how the associated enhanced expansion rate due to the stiff phase might facilitate a first-order electroweak symmetry breaking phase transition, in light of the recently measured value of the Higgs boson mass.

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.

We present estimates of the nonlinear bias of cosmological halo formation, spanning a wide range in the halo mass from $\sim 10^{5} M_\odot$ to $\sim 10^{12} M_\odot$, based upon both a suite of high-resolution cosmological N-body simulations and theoretical predictions. The halo bias is expressed in terms of the mean bias and stochasticity as a function of local overdensity ($\delta$), under different filtering scales, which is realized as the density of individual cells in uniform grids. The sampled overdensities span a range wide enough to provide the fully nonlinear bias effect on the formation of haloes. A strong correlation between $\delta$ and halo population overdensity $\delta_h$ is found, along with sizable stochasticity. We find that the empirical mean halo bias matches, with good accuracy, the prediction by the peak-background split method based on the excursion set formalism, as long as the empirical, globally-averaged halo mass function is used. Consequently, this bias formalism is insensitive to uncertainties caused by varying halo identification schemes, and can be applied generically. We also find that the probability distribution function of biased halo numbers has wider distribution than the pure Poisson shot noise, which is attributed to the sub-cell scale halo correlation. We explicitly calculate this correlation function and show that both overdense and underdense regions have positive correlation, leading to stochasticity larger than the Poisson shot noise in the range of haloes and halo-collapse epochs we study.

The study of the Universe on ultra-large scales is one of the major science cases for the Square Kilometre Array (SKA). The SKA will be able to probe a vast volume of the cosmos, thus representing a unique instrument, amongst next-generation cosmological experiments, for scrutinising the Universe's properties on the largest cosmic scales. Probing cosmic structures on extremely large scales will have many advantages. For instance, the growth of perturbations is well understood for those modes, since it falls fully within the linear regime. Also, such scales are unaffected by the poorly understood feedback of baryonic physics. On ultra-large cosmic scales, two key effects become significant: primordial non-Gaussianity and relativistic corrections to cosmological observables. Moreover, if late-time acceleration is driven not by dark energy but by modifications to general relativity, then such modifications should become apparent near and above the horizon scale. As a result, the SKA is forecast to deliver transformational constraints on non-Gaussianity and to probe gravity on super-horizon scales for the first time.

(Abridged) Despite the success of Cold Dark Matter (CDM) in explaining a wide range of observations, the microscopic nature of dark matter is still unknown. Attempts to detect WIMPs, the most commonly studied form, have not yet succeeded. Some predictions are even in apparent conflict with observations (e.g. cuspy-cored halos and `missing satellite problem'). Therefore, it is important to consider different forms of CDM. We focus on dark matter comprised of ultralight bosons that form a Bose-Einstein Condensate (BEC), described by a complex scalar field, and solve the Klein-Gordon and Einstein field equations to evolve the Friedmann-Robertson-Walker universe. We find, in addition to phases of radiation-, matter- and \Lambda-domination, an earlier phase of scalar-field-domination. Moreover, while WIMP CDM is non-relativistic at all times after it decouples, BEC scalar field dark matter (SFDM) is relativistic at early times, evolving from stiff to radiation-like, before it becomes non-relativistic (CDM-like) at late times. The timing of transitions between these phases yields fundamental constraints on SFDM model parameters, particle mass m and self-interaction coupling strength \lambda. We derive the range of particle parameters required to match observations of the evolving background universe, including the CMB and abundances of light elements produced by BBN, characterized by N_eff, the effective number of neutrino species, and the epoch of matter-radiation equality z_eq. This yields m >= 2.4*10^{-21} eV/c^2 and 9.5*10^{-19} eV^{-1}cm^3 <= \lambda/(mc^2)^2 <= 4*10^{-17} eV^{-1}cm^3. Our model accommodates current observations in which N_eff at BBN is higher than at z_eq, as probed by the CMB, otherwise unexplained by WIMP CDM. SFDM without self-interaction (`Fuzzy Dark Matter') cannot satisfy current BBN constraints within 68% confidence and is therefore disfavored.

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.

Measurements of the HI 21-cm power spectra from the reionization epoch will be influenced by the evolution of the signal along the line-of-sight direction of any observed volume. We use numerical as well as semi-numerical simulations of reionization in a cubic volume of 607 Mpc across to study this so-called light cone effect on the HI 21-cm power spectrum. We find that the light cone effect has the largest impact at two different stages of reionization: one when reionization is $\sim 20\%$ and other when it is $\sim 80\%$ completed. We find a factor of $\sim 4$ amplification of the power spectrum at the largest scale available in our simulations. We do not find any significant anisotropy in the 21-cm power spectrum due to the light cone effect. We argue that for the power spectrum to become anisotropic, the light cone effect would have to make the ionized bubbles significantly elongated or compressed along the line-of-sight, which would require extreme reionization scenarios. We also calculate the two-point correlation functions parallel and perpendicular to the line-of-sight and find them to differ. Finally, we calculate an optimum frequency bandwidth below which the light cone effect can be neglected when extracting power spectra from observations. We find that if one is willing to accept a $10 \%$ error due to the light cone effect, the optimum frequency bandwidth for $k= 0.056 \, \rm{Mpc}^{-1}$ is $\sim 7.5$ MHz. For $k = 0.15$ and $0.41 \, \rm{Mpc}^{-1}$ the optimum bandwidth is $\sim 11$ and $\sim 16$ MHz respectively.

The nature of the cosmological dark matter remains elusive. Recent studies have advocated the possibility that dark matter could be composed of ultra-light, self-interacting bosons, forming a Bose-Einstein condensate in the very early Universe. We consider models which are charged under a global U(1)-symmetry such that the dark matter number is conserved. It can then be described as a classical complex scalar field which evolves in an expanding Universe. We present a brief review on the bounds on the model parameters from cosmological and galactic observations, along with the properties of galactic halos which result from such a dark matter candidate.

If cosmological cold dark matter (CDM) consists of light enough bosonic particles that their phase-space density exceeds unity, they will comprise a Bose-Einstein condensate (BEC). The nature of this BEC-CDM as a quantum fluid may then distinguish it dynamically from the standard form of CDM involving a collisionless gas of non-relativistic particles that interact purely gravitationally. We summarize some of the dynamical properties of BEC-CDM that may lead to observable signatures in galactic halos and present some of the bounds on particle mass and self-interaction coupling strength that result from a comparison with observed galaxies.

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.

We present the largest-volume (425 Mpc/h=607 Mpc on a side) full radiative transfer simulation of cosmic reionization to date. We show that there is significant additional power in density fluctuations at very large scales. We systematically investigate the effects this additional power has on the progress, duration and features of reionization, as well as on selected reionization observables. We find that comoving simulation volume of ~100 Mpc/h per side is sufficient for deriving a convergent mean reionization history, but that the reionization patchiness is significantly underestimated. We use jackknife splitting to quantify the convergence of reionization properties with simulation volume for both mean-density and variable-density sub-regions. We find that sub-volumes of ~100 Mpc/h per side or larger yield convergent reionization histories, except for the earliest times, but smaller volumes of ~50 Mpc/h or less are not well converged at any redshift. Reionization history milestones show significant scatter between the sub-volumes, of Delta z=0.6-1 for ~50 Mpc/h volumes, decreasing to Delta z=0.3-0.5 for ~100 Mpc/h volumes, and $\Delta z$~0.1 for ~200 Mpc/h volumes. If we only consider mean-density sub-regions the scatter decreases, but remains at Delta z~0.1-0.2 for the different size sub-volumes. Consequently, many potential reionization observables like 21-cm rms, 21-cm PDF skewness and kurtosis all show good convergence for volumes of ~200 Mpc/h, but retain considerable scatter for smaller volumes. In contrast, the three-dimensional 21-cm power spectra at large scales (k<0.25 h/Mpc) do not fully converge for any sub-volume size. These additional large-scale fluctuations significantly enhance the 21-cm fluctuations, which should improve the prospects of detection considerably, given the lower foregrounds and greater interferometer sensitivity at higher frequencies. (abridged)

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.

One of the most promising ways to study the epoch of reionization (EoR) is through radio observations of the redshifted 21-cm line emission from neutral hydrogen. These observations are complicated by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, but they also encode information about cosmology and astrophysics. We study the effects of redshift space distortions on the power spectrum of 21-cm radiation from the EoR using large scale $N$-body and radiative transfer simulations. We quantify the anisotropy introduced in the 21-cm power spectrum by redshift space distortions and show how it evolves as reionization progresses and how it relates to the underlying physics. We go on to study the effects of redshift space distortions on LOFAR observations, taking instrument noise and foreground subtraction into account. We find that LOFAR should be able to directly observe the power spectrum anisotropy due to redshift space distortions at spatial scales around $k \sim 0.1$ Mpc$^{-1}$ after $\gtrsim$ 1000 hours of integration time. At larger scales, sample errors become a limiting factor, while at smaller scales detector noise and foregrounds make the extraction of the signal problematic. Finally, we show how the astrophysical information contained in the evolution of the anisotropy of the 21-cm power spectrum can be extracted from LOFAR observations, and how it can be used to distinguish between different reionization scenarios.

(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.

The 21cm background from the epoch of reionization is a promising cosmological probe: line-of-sight velocity fluctuations distort redshift, so brightness fluctuations in Fourier space depend upon angle, which linear theory shows can separate cosmological from astrophysical information. Nonlinear fluctuations in ionization, density and velocity change this, however. The validity and accuracy of the separation scheme are tested here for the first time, by detailed reionization simulations. The scheme works reasonably well early in reionization (< 40% ionized), but not late (> 80% ionized).

We calculate the angular power spectrum of the Cosmic Microwave Background (CMB) temperature fluctuations induced by the kinetic Sunyaev-Zel'dovich (kSZ) effect from the epoch of reionization (EOR). We use detailed N-body+radiative transfer simulations to follow inhomogeneous reionization of the intergalactic medium (IGM). For the first time we take into account the "self-regulation" of reionization: star formation in low-mass dwarf galaxies (10^8 M_\sun \lesssim M \lesssim 10^9 M_\sun) or minihalos (10^5 M_\sun \lesssim M \lesssim 10^8 M_\sun) is suppressed if these halos form in the regions that were already ionized or Lyman-Werner dissociated. Some previous work suggested that the amplitude of the kSZ power spectrum from the EOR can be described by a two-parameter family: the epoch of half ionization and the duration of reionization. However, we argue that this picture applies only to simple forms of the reionization history which are roughly symmetric about the half-ionization epoch. In self-regulated reionization, the universe begins to be ionized early, maintains a low level of ionization for an extended period, and then finishes reionization as soon as high-mass atomically-cooling halos dominate. While inclusion of self-regulation affects the amplitude of the kSZ power spectrum only modestly (\sim 10 %), it can change the duration of reionization by a factor of more than two. We conclude that the simple two-parameter family does not capture the effect of a physical, yet complex, reionization history caused by self-regulation. When added to the post-reionization kSZ contribution, our prediction for the total kSZ power spectrum is below the current upper bound from the South Pole Telescope. Therefore, the current upper bound on the kSZ effect from the EOR is consistent with our understanding of the physics of reionization.

Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of reionization on cosmic radiation backgrounds, including the CMB, near-IR, and redshifted 21cm backgrounds. We summarize some of our recent progress in this work, and address the question of whether observations of such signature effects can be used to distinguish the relative contributions of galaxies of different masses to reionization.

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]

We use numerical simulations to study the effects of the patchiness of a partly reionized intergalactic medium (IGM) on the observability of Ly-alpha emitters (LAEs) at high redshifts (z ~ 6). We present a new model that divides the Ly-alpha radiative transfer into a (circum-)galactic and an extragalactic (IGM) part, and investigate how the choice of intrinsic line model affects the IGM transmission results. We use our model to study the impact of neutral hydrogen on statistical observables such as the Ly-alpha restframe equivalent width (REW) distribution, the LAE luminosity function and the two-point correlation function. We find that if the observed changes in LAE luminosity functions and equivalent width distributions between z ~ 6 and z ~ 7 are to be explained by an increased IGM neutral fraction alone, we require an extremely late and rapid reionization scenario, where the Universe was ~ 40 % ionized at z = 7, ~ 50 % ionized at z = 6.5 and ~ 100 % ionized at z = 6. This is in conflict with other observations, suggesting that intrinsic LAE evolution at z > 6 cannot be completely neglected. We show how the two-point correlation function can provide more robust constraints once future observations obtain larger LAE samples, and provide predictions for the sample sizes needed to tell different reionization scenarios apart.

The intergalactic medium was reionized before redshift z~6, most likely by starlight which escaped from early galaxies. The very first stars formed when hydrogen molecules (H2) cooled gas inside the smallest galaxies, minihalos of mass between 10^5 and 10^8 solar masses. Although the very first stars began forming inside these minihalos before redshift z~40, their contribution has, to date, been ignored in large-scale simulations of this cosmic reionization. Here we report results from the first reionization simulations to include these first stars and the radiative feedback that limited their formation, in a volume large enough to follow the crucial spatial variations that influenced the process and its observability. We show that, while minihalo stars stopped far short of fully ionizing the universe, reionization began much earlier with minihalo sources than without, and was greatly extended, which boosts the intergalactic electron-scattering optical depth and the large-angle polarization fluctuations of the cosmic microwave background significantly. Although within current WMAP uncertainties, this boost should be readily detectable by Planck. If reionization ended as late as z_ov<~7, as suggested by other observations, Planck will thereby see the signature of the first stars at high redshift, currently undetectable by other probes.

Observations of redshifted 21-cm radiation from neutral hydrogen during the epoch of reionization (EoR) are considered to constitute the most promising tool to probe that epoch. One of the major goals of the first generation of low frequency radio telescopes is to measure the 3D 21-cm power spectrum. However, the 21-cm signal could evolve substantially along the line of sight (LOS) direction of an observed 3D volume, since the received signal from different planes transverse to the LOS originated from different look-back times and could therefore be statistically different. Using numerical simulations we investigate this so-called light cone effect on the spherically averaged 3D 21-cm power spectrum. For this version of the power spectrum, we find that the effect mostly `averages out' and observe a smaller change in the power spectrum compared to the amount of evolution in the mean 21-cm signal and its rms variations along the LOS direction. Nevertheless, changes up to 50% at large scales are possible. In general the power is enhanced/suppressed at large/small scales when the effect is included. The cross-over mode below/above which the power is enhanced/suppressed moves toward larger scales as reionization proceeds. When considering the 3D power spectrum we find it to be anisotropic at the late stages of reionization and on large scales. The effect is dominated by the evolution of the ionized fraction of hydrogen during reionization and including peculiar velocities hardly changes these conclusions. We present simple analytical models which explain qualitatively all the features we see in the simulations.

We present a study of the impact of a bright quasar on the redshifted 21cm signal during the Epoch of Reionization (EoR). Using three different cosmological radiative transfer simulations, we investigate if quasars are capable of substantially changing the size and morphology of the H II regions they are born in. We choose stellar and quasar luminosities in a way that is favourable to seeing such an effect. We find that even the most luminous of our quasar models is not able to increase the size of its native H II region substantially beyond those of large H II regions produced by clustered stellar sources alone. However, the quasar H II region is found to be more spherical. We next investigate the prospects of detecting such H II regions in the redshifted 21cm data from the Low Frequency Array (LOFAR) by means of a matched filter technique. We find that H II regions with radii ~ 25 comoving Mpc or larger should have a sufficiently high detection probability for 1200 hours of integration time. Although the matched filter can in principle distinguish between more and less spherical regions, we find that when including realistic system noise this distinction can no longer be made. The strong foregrounds are found not to pose a problem for the matched filter technique. We also demonstrate that when the quasar position is known, the redshifted 21cm data can still be used to set upper limits on the ionizing photon rate of the quasar. If both the quasar position and its luminosity are known, the redshifted 21 cm data can set new constrains on quasar lifetimes.

We use the results of large-scale simulations of reionization to explore methods for characterizing the topology and sizes of HII regions during reionization. We use four independent methods for characterizing the sizes of ionized regions. Three of them give us a full size distribution: the friends-of-friends (FOF) method, the spherical average method (SPA) and the power spectrum (PS) of the ionized fraction. These latter three methods are complementary: While the FOF method captures the size distribution of the small scale H II regions, which contribute only a small amount to the total ionization fraction, the spherical average method provides a smoothed measure for the average size of the H II regions constituting the main contribution to the ionized fraction, and the power spectrum does the same while retaining more details on the size distribution. Our fourth method for characterizing the sizes of the H II regions is the average size which results if we divide the total volume of the H II regions by their total surface area, (i.e. 3V/A), computed in terms of the ratio of the corresponding Minkowski functionals of the ionized fraction field. To characterize the topology of the ionized regions, we calculate the evolution of the Euler Characteristic. We find that the evolution of the topology during the first half of reionization is consistent with inside-out reionization of a Gaussian density field. We use these techniques to investigate the dependence of size and topology on some basic source properties, such as the halo mass-to-light ratio, susceptibility of haloes to negative feedback from reionization, and the minimum halo mass for sources to form. We find that suppression of ionizing sources within ionized regions slows the growth of H II regions, and also changes their size distribution. Additionally, the topology of simulations including suppression is more complex. (abridged)

Current observations suggest that the universe was reionized sometime before z~6. One way to observe this epoch of the universe is through the Near Infrared Background (NIRB), which contains information about galaxies which may be too faint to be observed individually. We calculate the angular power spectrum (C_l) of the NIRB fluctuations caused by the distribution of these galaxies. Assuming a complete subtraction of any post-reionization component, C_l will be dominated by galaxies responsible for completing reionization (e.g., z~6). The shape of C_l at high l is sensitive to the amount of non-linear bias of dark matter halos hosting galaxies. As the non-linear bias depends on the mass of these halos, we can use the shape of C_l to infer typical masses of dark matter halos responsible for completing reionization. We extend our previous study by using a higher-resolution N-body simulation, which can resolve halos down to 10^8 M_sun. We also include improved radiative transfer, which allows for the suppression of star formation in small-mass halos due to photo-ionization heating. As the non-linear bias enhances the dark-matter-halo power spectrum on small scales, we find that C_l is steeper for the case with a complete suppression of small sources or partial suppression of star formation in small halos (the minimum galaxy mass is M_min=10^9 M_sun in ionized regions and M_min=10^8 M_sun in neutral regions) than the case in which these small halos were unsuppressed. In all cases, we do not see a turn-over toward high l in the shape of l^2 C_l.

(Abridged) Extensions of the standard model of particle physics predict very light bosons, ranging from about 10^{-5} eV for the QCD axion to 10^{-33} eV for ultra-light particles, which could be the cold dark matter (CDM) in the Universe. If so, their phase-space density must be high enough to form a Bose-Einstein condensate (BEC). The fluid-like nature of BEC-CDM dynamics differs from that of standard collisionless CDM (sCDM), so observations of galactic haloes may distinguish them. sCDM has problems with galaxy observations on small scales, which BEC-CDM may overcome for a large range of particle mass m and self-interaction strength g. For quantum-coherence on galactic scales of radius R and mass M, either the de-Broglie wavelength lambda_deB <~ R, requiring m >~ m_H \cong 10^{-25}(R/100 kpc)^{-1/2}(M/10^{12} M_solar)^{-1/2} eV, or else lambda_deB << R but self-interaction balances gravity, requiring m >> m_H and g >> g_H \cong 2 x 10^{-64} (R/100 kpc)(M/10^{12} M_solar)^{-1} eV cm^3. Here we study the largely-neglected effects of angular momentum. Spin parameters lambda \cong 0.05 are expected from tidal-torquing by large-scale structure, just as for sCDM. Since lab BECs develop quantum vortices if rotated rapidly enough, we ask if this angular momentum is sufficient to form vortices in BEC haloes, affecting their structure with potentially observable consequences. The minimum angular momentum for this, L_{QM} = $\hbar M/m$, requires m >= 9.5 m_H for lambda = 0.05, close to the particle mass required to influence structure on galactic scales. We study the equilibrium of self-gravitating, rotating BEC haloes which satisfy the Gross-Pitaevskii-Poisson equations, to calculate if and when vortices are energetically favoured. Vortices form as long as self-interaction is strong enough, which includes a large part of the range of m and g of interest for BEC-CDM haloes.

The peculiar velocity of the intergalactic gas responsible for the cosmic 21cm background from the epoch of reionization and beyond introduces an anisotropy in the three-dimensional power spectrum of brightness temperature fluctuations. Measurement of this anisotropy by future 21cm surveys is a promising tool for separating cosmology from 21cm astrophysics. However, previous attempts to model the signal have often neglected peculiar velocity or only approximated it crudely. This paper re-examines the effects of peculiar velocity on the 21cm signal in detail, improving upon past treatment and addressing several issues for the first time. (1) We show that properly accounting for finite optical depth eliminates the unphysical divergence of 21cm brightness temperature in overdense regions of the IGM found by previous work that employed the usual optically-thin approximation. (2) The approximation made previously to circumvent the diverging brightness temperature problem by capping velocity gradient can misestimate the power spectrum on all scales. (3) The observed power spectrum in redshift-space remains finite even in the optically-thin approximation if one properly accounts for the redshift-space distortion. However, results that take full account of finite optical depth show that this approximation is only accurate in the limit of high spin temperature. (4) The linear theory for redshift-space distortion results in ~30% error in the observationally relevant wavenumber range, at the 50% ionized epoch. (5) We describe and test two numerical schemes to calculate the 21cm signal from reionization simulations to incorporate peculiar velocity effects in the optically-thin approximation accurately. One is particle-based, the other grid-based, and while the former is most accurate, we demonstrate that the latter is computationally more efficient and can achieve sufficient accuracy. [Abridged]

We present an extension to the short-characteristic ray-tracing and non-equilibrium photon-ionization code C2Ray. The new version includes the effects of helium and improved multi-frequency heating. The motivation for this work is to be able to deal with harder ionizing spectra, such as for example from quasar-like sources during cosmic reionization. We review the basic algorithmic ingredients of C2-Ray before describing the changes implemented, which include a treatment of the full on the spot (OTS) approximation, secondary ionization, and multi-frequency photo-ionization and heating. We performed a series of tests against equilibrium solutions from CLOUDY as well as comparisons to the hydrogen only solutions by C2-Ray in the extensive code comparison in Iliev et al. (2006). We show that the full, coupled OTS approximation is more accurate than the simplified, uncoupled one. We find that also with helium and a multi-frequency set up, long timesteps (up to ~10% of the recombination time) still give accurate results for the ionization fractions. On the other hand, accurate results for the temperature set strong constrains on the timestep. The details of these constraints depend however on the optical depth of the cells. We use the new version of the code to confirm that the assumption made in many reionization simulations, namely that helium is singly ionized everywhere were hydrogen is, is indeed valid when the sources have stellar-like spectra.

We study the capability of the James Webb Space Telescope (JWST) to detect Supermassive Dark Stars (SMDS). If the first stars are powered by dark matter heating in triaxial dark matter haloes, they may grow to be very large and very bright, visible in deep imaging with JWST and even Hubble Space Telescope (HST). We use HST surveys to place bounds on the numbers of SMDSs that may be detected in future JWST imaging surveys. We showed that SMDS in the mass range $10^6-10^7 M_\odot$ are bright enough to be detected in all the wavelength bands of the NIRCam on JWST . If SMDSs exist at z ~10, 12, and 14, they will be detectable as J-band, H-band, or K-band dropouts, respectively. With a total survey area of 150 arcmin^2 (assuming a multi-year deep parallel survey with JWST), we find that typically the number of $10^6 M_\odot$ SMDSs found as H or K-band dropouts is ~10^5\fsmds, where the fraction of early DM haloes hosting DS is likely to be small, \fsmds<<1. If the SDMS survive down to z=10 where HST bounds apply, then the observable number of SMDSs as H or K-band dropouts with JWST is ~1-30. While individual SMDS are bright enough to be detected by JWST, standard PopIII stars are not, and would only be detected in first galaxies with total stellar masses of ~$10^6-10^8 M_\odot$. Differentiating first galaxies at z>10 from SMDSs would be possible with spectroscopy: the SMDS (which are too cool produce significant nebular emission) will have only absorption lines while the galaxies are likely to produce emission lines as well. Of particular interest would be the 1640 HeII emission line as well as H{\alpha} lines which would be signatures of early galaxies rather than SMDSs. The detection of SMDSs would not only provide alternative evidence for WIMPs but would also provide possible seeds for the formation of supermassive black holes that power QSOs at z~6.

The prospect of detecting the first galaxies by observing their impact on the intergalactic medium as they reionized it during the first billion years leads us to ask whether such indirect observations are capable of diagnosing which types of galaxies were most responsible for reionization. We attempt to answer this by considering a set of large-scale radiative transfer simulations of reionization in sufficiently large volumes to make statistically meaningful predictions of observable signatures, while also directly resolving all atomically-cooling halos down to 10^8 M_solar. We focus here on predictions of the 21-cm background, to see if upcoming observations are capable of distinguishing a universe ionized primarily by high-mass halos from one in which both high-mass and low-mass halos are responsible, and to see how these results depend upon the uncertain source efficiencies. We find that 21-cm fluctuation power spectra observed by the first generation EoR/21-cm radio interferometer arrays should be able to distinguish the case of reionization by high-mass halos alone from that by both high- and low-mass halos, together. Some reionization scenarios yield very similar power spectra and rms evolution and thus can only be discriminated by their different mean reionization history and 21-cm PDF distributions. We find that the skewness of the 21-cm PDF distribution smoothed over LOFAR-like window shows a clear feature correlated with the rise of the rms due to patchiness. Measurements of the mean photoionization rates are sensitive to the average density of the regions being studied and therefore could be strongly skewed in certain cases. (abridged)

Self-Interacting Dark Matter (SIDM) is a collisional form of cold dark matter (CDM), originally proposed to solve problems that arose when the collisionless CDM theory of structure formation was compared with observations of galaxies on small scales. The quantitative impact of the proposed elastic collisions on structure formation has been estimated previously by Monte Carlo N-body simulations and by a conducting fluid model, with apparently diverging results. To improve this situation, we make direct comparisons between new Monte Carlo N-body simulations and solutions of the conducting fluid model, for isolated SIDM haloes of fixed mass. This allows us to separate cleanly the effects of gravothermal relaxation from those of continuous mass accretion in an expanding background universe. When these two methods were previously applied to halo formation with cosmological boundary conditions, they disagreed by an order of magnitude about the size of the scattering cross section required to solve the so-called 'cusp-core problem.' We show here, however, that the methods agree with each other within 20 per cent for isolated haloes. This suggests that the two methods are consistent, and that their disagreement for cosmological haloes is not caused by a breakdown of their validity. The isolated haloes studied here undergo gravothermal collapse. We compare the solutions calculated by these two methods for gravothermal collapse starting from several initial conditions. This allows us to calibrate the heat conduction which accounts for the effect of elastic hard-sphere scattering in the fluid model. The amount of tuning of the thermal conductivity parameters required to bring the two methods into close agreement for isolated haloes, however, is too small to explain the discrepancy found previously in the cosmological context.

Some of the first stars could be cooler and more massive than standard stellar models would suggest, due to the effects of dark matter annihilation in their cores. It has recently been argued that such objects may attain masses in the 10^4--10^7 solar mass range, and that such supermassive dark stars should be within reach of the upcoming James Webb Space Telescope. Notwithstanding theoretical difficulties with this proposal, we argue here that some of these objects should also be readily detectable with both the Hubble Space Telescope and ground-based 8--10 m class telescopes. Existing survey data already place strong constraints on 10^7 solar mass dark stars at z~10. We show that such objects must be exceedingly rare or short-lived to have avoided detection.

We present some preliminary results from a series of extremely large, high-resolution N-body simulations of the formation of early nonlinear structures. We find that the high-z halo mass function is inconsistent with the Sheth-Tormen mass function, which tends to over-estimate the abundance of rare halos. This discrepancy is in rough agreement with previous results based on smaller simulations. We also show that the number density of minihaloes is correlated with local matter density, albeit with a significant scatter that increases with redshift, as minihaloes become increasingly rare. The average correlation is in rough agreement with a simple analytical extended Press-Schechter model, but can differ by up to factor of 2 in some regimes.

If cold dark matter elementary particles form a Bose-Einstein condensate, their superfluidity may distinguish them from other forms of cold dark matter, including creation of quantum vortices. We demonstrate here that such vortices are favoured in strongly-coupled condensates, while this is not the case for axions, which are generally presumed to form a Bose-Einstein condensate but are effectively non-interacting.

The Near Infrared Background (NIRB) is one of a few methods that can be used to observe the redshifted light from early stars at a redshift of six and above. Fluctuations of the NIRB can provide information on the first structures, such as halos and their surrounding ionized regions in the IGM. We combine, for the first time, N-body simulations, radiative transfer code, and analytic calculations of luminosity of early structures to predict the angular power spectrum (C_l) of fluctuations in the NIRB. We study the effects of various assumptions about the stellar mass, the initial mass spectrum of stars, metallicity, the star formation efficiency (f_*), the escape fraction of ionizing photons (f_esc), and the star formation timescale (t_SF), on the amplitude as well as the shape of C_l. The power spectrum of NIRB fluctuations is maximized when f_* is the largest (as C_l ~ (f_*)^2) and f_esc is the smallest. A significant uncertainty in the predicted amplitude of C_l exists due to our lack of knowledge of t_SF of these galaxies, which is equivalent to our lack of knowledge of the mass-to-light ratio. We do not see a turnover in the NIRB angular power spectrum of the halo contribution and explain this as the effect of high levels of non-linear bias. This is partly due to our choice of the minimum mass of halos contributing to NIRB, and a smaller minimum mass, which has a smaller non-linear bias, may still exhibit a turn over. Therefore, both the amplitude and shape of the NIRB power spectrum provide important information regarding the nature of sources contributing to the cosmic reionization. The angular power spectrum of the IGM, in most cases, is much smaller than the halo angular power spectrum. In addition, low levels of the observed mean background intensity tend to rule out high values of f_* > 0.2.

The 21-cm PDF (i.e., distribution of pixel brightness temperatures) is expected to be highly non-Gaussian during reionization and to provide important information on the distribution of density and ionization. We measure the 21-cm PDF as a function of redshift in a large simulation of cosmic reionization and propose a simple empirical fit. Guided by the simulated PDF, we then carry out a maximum likelihood analysis of the ability of upcoming experiments to measure the shape of the 21-cm PDF and derive from it the cosmic reionization history. Under the strongest assumptions, we find that upcoming experiments can measure the reionization history in the mid to late stages of reionization to 1-10% accuracy. Under a more flexible approach that allows for four free parameters at each redshift, a similar accuracy requires the lower noise levels of second-generation 21-cm experiments.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

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

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

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

The first, self-consistent calculations are presented of the cosmological, H2-dissociating UV background produced during the epoch of reionization by the sources of reionization. Large-scale radiative transfer simulations of reionization trace the impact of all the ionizing starlight on the IGM from all the sources in our simulation volume down to dwarf galaxies of mass ~10^8 Msun, identified by very high-resolution N-body simulations, including the self-regulating effect of IGM photoheating on dwarf galaxy formation. The UV continuum emitted below 13.6 eV by each source is then transferred through the same IGM, attenuated by atomic H Lyman series resonance lines, to predict the evolution of the inhomogeneous radiation background in the Lyman-Werner bands of H2 between 11 and 13.6 eV. On average, the intensity of this Lyman-Werner background is found to rise to the threshold level at which dissociation suppresses H2 cooling and star formation inside minihalos, long before reionization is complete. Spatial variations in the Lyman-Werner background are found which result from the clustering of sources associated with large-scale structure formation, such that intensity fluctuations correlate with matter density fluctuations. As a result, the Lyman-Werner background rises to the threshold level for H2 suppression earlier in the vicinity of the reionization sources and their H II regions.

We study the formation of disk-dominated galaxies in a LCDM universe. Their existence is considered to be a challenge for the LCDM cosmology, because galaxy mergers isotropize stellar disks and trigger angular momentum transport in gas disks, thus fostering the formation of central stellar spheroids. Here, we postulate that the formation of stellar spheroids from gasrich disks is controlled by two parameters that characterize galaxy mergers, the mass ratio of merging dark matter halos, and the virial velocity of the larger merging halo. We utilize merger histories generated from realizations of the cosmological density field to calculate the fraction of dark matter halos that have avoided spheroid formation, and compare the derived statistics with the spheroid occupation fractions in surveys of nearby galaxies. We find, for example, that the survival rate of disk-dominated galaxies in LCDM is just high enough to explain the observed fractional representation of disk-dominated galaxies in the universe if the only mergers which lead to central spheroid formation are those with mass ratios M2/M1 > 0.3 and virial velocities Vvir,1 > 55 km/s. We discuss the physical origin of this criterion, and show that the dependence of the disk-dominated fraction on galaxy mass provides a further test of the merger hypothesis. [For additional details, see, Koda et al. (2009).]

Observations of high-redshift Ly-alpha sources are a major tool for studying the high-redshift Universe. We discuss the effect of the reionizing intergalactic medium on the observability of Ly-alpha sources based on large simulations of early structure formation with radiative transfer. This takes into account self-consistently the reionization history, density, velocity and ionization structures and nonlinear source clustering. We find that all fields are highly anisotropic and as a consequence there are very large variations in opacity among the different lines-of-sight. The velocity effects, from both infall and source peculiar velocity are most important for the luminous sources, affecting the line profile and depressing the bright end of the luminosity function. The line profiles are generally asymmetric and the line centers of the luminous sources are always absorbed due to the high density of the local IGM. For both luminous and average sources the damping wing effects are of similar magnitude and remain significant until fairly late. The ionizing flux in the ionized patch surrounding a high density peak is generally strongly dominated, particularly at late times, by the cluster of faint sources, rather than the central massive galaxy. The IGM absorption does not change appreciably the correlation function of sources at high redshift. Our derived luminosity function assuming constant mass-to-light ratio provides an excellent match to the shape of the observed luminosity function at z=6.6 with faint-end slope of alpha=-1.5. The resulting mass-to-light ratio implies that the majority of sources responsible for reionization are too faint to be observed by the current surveys. (abridged)

The first, self-consistent calculations of the cosmological H_2 dissociating UV background produced during the epoch of reionization (EOR) by the sources of reionization are presented. Large-scale radiative transfer simulations of reionization trace the impact of all the ionizing starlight on the IGM from all the sources in our simulation volume down to dwarf galaxies of mass ~ 10^8 solar mass, identified by very high-resolution N-body simulations, including the self-regulating effect of IGM photoheating on dwarf galaxy formation. The UV continuum emitted below 13.6 eV by each source is then transferred through the same IGM, attenuated by atomic H Lyman series resonance lines, to predict the evolution of the inhomogeneous background in the Lyman-Werner band of H_2 between 11 and 13.6 eV.

The redshifted 21-cm line of distant neutral H atoms provides a probe of the cosmic ``dark ages'' and the epoch of reionization (``EOR'') which ended them. The radio continuum produced by this redshifted line can be seen in absorption or emission against the CMB at meterwaves, yielding information about the thermal and ionization history of the universe and the primordial density perturbation spectrum that led to galaxy and large-scale structure formation. Observing this 21-cm background is a great challenge. A new generation of low-frequency radio arrays is currently under development to search for this background. Accurate theoretical predictions of the spectrum and anisotropy of this background, necessary to guide and interpret future observations, are also quite challenging. It is necessary to model the inhomogeneous reionization of the intergalactic medium and determine the spin temperature of the 21-cm transition and its variations in time and space as it decouples from the temperature of the CMB. Here, we focus on just a few of the predictions for the 21-cm background from the EOR, based on our newest, large-scale simulations of patchy reionization. These simulations are the first with enough N-body particles (from 5 to 29 billion) and radiative transfer rays to resolve the formation of and trace the ionizing radiation from each of the millions of dwarf galaxies believed responsible for reionization, down to 10^8 M_solar, in a cubic volume large enough (90 and 163 comoving Mpc on a side) to make meaningful statistical predictions of the fluctuating 21-cm background. (abridged)

The Cosmic Dark Ages and the Epoch of Reionization constitute a crucial missing link in our understanding of the evolution of the intergalactic medium and the formation and evolution of galaxies. Due to the complex nature of this global process it is best studied through large-scale numerical simulations. This presents considerable computational challenges. The dominant contributors of ionizing radiation were dwarf galaxies. These tiny galaxies must be resolved in very large cosmological volumes in order to derive their clustering properties and the corresponding observational signatures correctly, which makes this one of the most challenging problems of numerical cosmology. We have recently performed the largest and most detailed simulations of the formation of early cosmological large-scale structures and their radiative feedback leading to cosmic reionization. This was achieved by running extremely large (up to 29 billion-particle) N-body simulations of the formation of the Cosmic Web, with enough particles and sufficient force resolution to resolve all the galactic halos with total masses larger than 10^8 Solar masses in computational volumes of up to (163 Mpc)^3. These results were then post-processed by propagating the ionizing radiation from all sources by using fast and accurate ray-tracing radiative transfer method. Both of our codes are parallelized using a combination of MPI and OpenMP and to this date have been run efficiently on up to 2048 cores (N-body) and up to 10000 cores (radiative transfer) on the newly-deployed Sun Constellation Linux Cluster at the Texas Advanced Computing Center. In this paper we describe our codes, parallelization strategies, scaling and some preliminary scientific results. (abridged)

The merging galaxy cluster 1E 0657-56, known as the "bullet cluster," is one of the hottest clusters known. The X-ray emitting plasma exhibits bow-shock-like temperature and density jumps. The segregation of this plasma from the peaks of the mass distribution determined by gravitational lensing has been interpreted as a direct proof of collisionless dark matter. If the high shock speed inferred from the shock jump conditions equals the relative speed of the merging CDM halos, however, this merger is predicted to be such a rare event in a LCDM universe that observing it presents a possible conflict with the LCDM model. We examined this question using high resolution, 2D simulations of gas dynamics in cluster collisions to analyze the relative motion of the clusters, the bow shock, and the contact discontinuity, and relate these to the X-ray data for the bullet cluster. We find that the velocity of the fluid shock need not equal the relative velocity of the CDM components. An illustrative simulation finds that the present relative velocity of the CDM halos is 16% lower than that of the shock. While this conclusion is sensitive to the detailed initial mass and gas density profiles of the colliding clusters, such a decrease of the inferred halo relative velocity would significantly increase the likelihood of finding 1E 0657-56 in a LCDM universe. (Conference proceedings based on a poster at Bash Symposium 2007)

Direct detection of the Dark Ages and the Epoch of Reionization (EOR) is among the main scientific objectives of all current and future low-frequency radio facilities. In this paper we summarize and discuss recent results, based on state-of-the-art numerical simulations, regarding the fundamental EOR properties and its observability with current and future radio arrays, like the Giant Metrewave Radio Telescope (GMRT), the Low Frequency Array (LOFAR), the 21-CM Array (21CMA), the Murchison Widefield Array (MWA) and the Square Kilometre Array (SKA). Results show that the optimal observational frequencies for statistical detection are 140-160 MHz. The signals are strongly non-Gaussian at late times. The correlation widths between 21-cm maps at neighbouring frequencies are short, of order 300-800 kHz, which should help with the cleaning of the strong foregrounds. Direct comparison of the resolutions and expected sensitivities of GMRT and MWA indicate that their optimal sensitivity ranges are similar, at scales k~0.2-0.4 h/Mpc, however, all else being equal the former should require shorter integration times due to its significantly larger collecting area.

We study the formation of disk-dominated galaxies in a Lambda CDM universe. Their existence is considered to be a challenge for the Lambda CDM cosmology, because galaxy mergers isotropize stellar disks and trigger angular momentum transport in gas disks, thus fostering the formation of central stellar spheroids. Here, we postulate that the formation of stellar spheroids from gas-rich disks is controlled by two parameters that characterize galaxy mergers, the mass ratio of merging dark matter halos, and the virial velocity of the larger merging halo. We utilize merger histories generated from realizations of the cosmological density field to calculate the fraction of dark matter halos that have avoided spheroid formation, and compare the derived statistics with the spheroid occupation fractions in surveys of nearby galaxies. We find, for example, that the survival rate of disk-dominated galaxies in Lambda CDM is just high enough to explain the observed fractional representation of disk-dominated galaxies in the universe if the only mergers which lead to central spheroid formation are those with mass ratios M2/M1 > 0.3 and virial velocities Vvir,1 > 55 km/s. We discuss the physical origin of this criterion.

Recently, we have presented the first large-scale radiative transfer simulations of reionization. Here we present new simulations which extend the source halo mass range downward to 10^8M_solar, to capture the full range of halo masses thought to be primarily responsible for reionization by their star formation following atomic hydrogen radiative cooling and gravitational collapse. Haloes below about 10^9M_solar, however, are subject to Jeans-mass filtering in the ionized regions, which suppresses their baryonic content and their ability to release ionizing radiation. By including these smaller-mass haloes but accounting for their suppression, too, we find that reionization is ``self-regulating,'' as follows. As the mean ionized fraction rises, so does the fraction of the volume within which suppression occurs. Hence, the degree of suppression is related to the mean ionized fraction. Since low-mass haloes with high emissivity achieve a given mean ionized fraction earlier than do those with low efficiency, Jeans-mass filtering compensates for the difference in the emissivity of the suppressible haloes in these two cases. As a result, in the presence of lower-mass source haloes, reionization begins earlier, but the later stages of reionization and the time of overlap are dictated by the efficiency of the higher-mass haloes, independent of the efficiency of the suppressible, lower-mass haloes. Reionization histories consistent with current observational constraints are shown to be achievable with standard stellar sources in haloes above 10^8M_solar. Neither minihalos nor exotic sources are required, and the phenomenon of ``double reionization'' previously suggested does not occur. (abridged)

A number of large current experiments aim to detect the signatures of the Cosmic Reionization at redshifts z>6. Their success depends crucially on understanding the character of the reionization process and its observable consequences and designing the best strategies to use. We use large-scale simulations of cosmic reionization to evaluate the reionization signatures at redshifted 21-cm and small-scale CMB anisotropies in the best current model for the background universe, with fundamental cosmological parameters given by WMAP 3-year results (WMAP3). We find that the optimal frequency range for observing the ``global step'' of the 21-cm emission is 120-150 MHz, while statistical studies should aim at 140-160 MHz, observable by GMRT. Some strongly-nongaussian brightness features should be detectable at frequencies up to ~190 MHz. In terms of sensitivity-signal trade-off relatively low resolutions, corresponding to beams of at least a few arcminutes, are preferable. The CMB anisotropy signal from the kinetic Sunyaev-Zel'dovich effect from reionized patches peaks at tens of muK at arcminute scales and has an rms of ~1 muK, and should be observable by the Atacama Cosmology Telescope and the South Pole Telescope. We discuss the various observational issues and the uncertainties involved, mostly related to the poorly-known reionization parameters and, to a lesser extend, to the uncertainties in the background cosmology.

Recently the numerical simulations of the process of reionization of the universe at z>6 have made a qualitative leap forward, reaching sufficient sizes and dynamic range to determine the characteristic scales of this process. This allowed making the first realistic predictions for a variety of observational signatures. We discuss recent results from large-scale radiative transfer and structure formation simulations on the observability of high-redshift Ly-alpha sources. We also briefly discuss the dependence of the characteristic scales and topology of the ionized and neutral patches on the reionization parameters.

During the epoch of reionization, the formation of radiation sources is accompanied by the growth of a H- photodissociating flux. We estimate the impact of this flux on the formation of molecular hydrogen and cooling in the first galaxies, assuming different types of radiation sources (e.g. Pop II and Pop III stars, miniquasars). We find that H- photodissociation reduces the formation of H2 molecules by a factor of ~1+1000k_s*x/(f_esc*delta), where x is the mean ionized fraction in the IGM, f_esc is the fraction of ionizing photons that escape from their progenitor halos, delta is the local gas overdensity and k_s is an order unity constant which depends on the type of radiation source. By the time a significant fraction of the universe becomes ionized, H- photodissociation may significantly reduce the H2 abundance and, with it, the primordial star formation rate, delaying the progress of reionization.

In recent years there has been considerable progress in our understanding of the nature and properties of the reionization process. In particular, the numerical simulations of this epoch have made a qualitative leap forward, reaching sufficiently large scales to derive the characteristic scales of the reionization process and thus allowing for realistic observational predictions. Our group has recently performed the first such large-scale radiative transfer simulations of reionization, run on top of state-of-the-art simulations of early structure formation. This allowed us to make the first realistic observational predictions about the Epoch of Reionization based on detailed radiative transfer and structure formation simulations. We discuss the basic features of reionization derived from our simulations and some recent results on the observational implications for the high-redshift Ly-alpha sources.

Shock waves driven in the intergalactic medium during the merging of galaxy clusters have been observed in X-ray imaging and spectroscopy. Fluid motions inferred from the shock strength and morphology can be compared to the cold dark matter (CDM) distribution inferred from gravitational lensing. A detailed reconstruction of the CDM kinematics, however, must take into account the nontrivial response of the fluid intracluster medium to the collisionless CDM motions. We have carried out two-dimensional simulations of gas dynamics in cluster collisions. We analyze the relative motion of the clusters, the bow shock wave, and the contact discontinuity and relate these to X-ray data. We focus on the "bullet cluster," 1E 0657-56, a near head-on collision of unequal-mass clusters, for which the gas density and temperature jumps across the prominent bow shock imply a high shock velocity 4,700 km/s. The velocity of the fluid shock has been widely interpreted as the relative velocity of the CDM components. This need not be the case, however. An illustrative simulation finds that the present relative velocity of the CDM halos is 16% lower than that of the shock. While this conclusion is sensitive to the detailed initial mass and gas density profile of the colliding clusters, such a decrease of the inferred halo relative velocity would increase the likelihood of finding 1E 0657-56 in a LambdaCDM universe.

The inhomogeneous ionization state of the universe when the first sources of ionizing radiation appeared should lead to anisotropies in the polarization of the cosmic microwave background. We use cosmological simulations of the process by which the first sources ionized the intergalactic medium to study the induced polarization anisotropies. We find that the polarization anisotropies have rms of order ~0.01 \mu K, and local peak values of ~0.1 \mu K, smaller than those due to gravitational lensing on small scales. The polarization direction is highly coherent over degree scales. This directional coherence is not expected from either primary anisotropy or gravitational lensing effects, making the largest signals due to inhomogeneous ionization relatively easy to isolate, should experiments achieve the necessary very low noise levels.

During the epoch of reionization a large number of photons were produced with frequencies below the hydrogen Lyman limit. After redshifting into the closest resonance, these photons underwent multiple scatterings with atoms. We examine the effect of these scatterings on the temperature of the neutral intergalactic medium (IGM). Continuum photons, emitted between the Ly_alpha and Ly_gamma frequencies, heat the gas after being redshifted into the H Ly_alpha or D Ly_beta resonance. By contrast, photons emitted between the Ly_gamma and Ly-limit frequencies, produce effective cooling of the gas. Prior to reionization, the equilibrium temperature of ~100 K for hydrogen and helium atoms is set by these two competing processes. At the same time, Ly_beta resonance photons thermally decouple deuterium from other species, raising its temperature as high as 10^4 K. Our results have important consequences for the cosmic 21-cm background and the entropy floor of the early IGM which can affect star formation and reionization.

We study the effect of starlight from the first stars on the ability of other minihaloes in their neighborhood to form additional stars. The question of what the dynamical consequences were for these target minihaloes, of their exposure to the ionizing and dissociating starlight from the Pop III star requires a full radiation-hydrodynamics calculation. Towards this end, we have performed a series of detailed, 1D, radiation-hydrodynamical simulations. We have varied the distance to the source (and, hence, the flux) and the mass and evolutionary stage of the target haloes to quantify this effect. We find: (1) trapping of the I-front and its transformation from R-type to D-type, preceded by a shock front; (2) photoevaporation of the ionized gas (i.e. all gas originally located outside the trapping radius); (3) formation of an H_2 precursor shell which leads the I-front, stimulated by partial photoionization; and (4) the shock-induced formation of H_2 in the minihalo neutral core when the shock speeds up and partially ionizes the gas. The fate of the neutral core is mostly determined by the response of the core to this shock front, which leads to molecular cooling and collapse that, when compared to the same halo without external radiation, is either: (a) expedited, (b) delayed, (c) unaltered, or (d) reversed or prevented, depending upon the flux (i.e. distance to the source) and the halo mass and evolutionary stage. Roughly speaking, most haloes that were destined to cool, collapse, and form stars in the absence of external radiation are found to do so even when exposed to the first Pop III star in their neighborhood, while those that would not have done so are still not able to. (abridged)

We present the first calculation of the kinetic Sunyaev-Zel'dovich (kSZ) effect due to the inhomogeneous reionization of the universe based on detailed large-scale radiative transfer simulations of reionization. The resulting sky power spectra peak at l=2000-8000 with maximum values of l^2C_l~1\times10^{-12}. The peak scale is determined by the typical size of the ionized regions and roughly corresponds to the ionized bubble sizes observed in our simulations, ~5-20 Mpc. The kSZ anisotropy signal from reionization dominates the primary CMB signal above l=3000. This predicted kSZ signal at arcminute scales is sufficiently strong to be detectable by upcoming experiments, like the Atacama Cosmology Telescope and South Pole Telescope which are expected to have ~1' resolution and ~muK sensitivity. The extended and patchy nature of the reionization process results in a boost of the peak signal in power by approximately one order of magnitude compared to a uniform reionization scenario, while roughly tripling the signal compared with that based upon the assumption of gradual but spatially uniform reionization. At large scales the patchy kSZ signal depends largely on the ionizing source efficiencies and the large-scale velocity fields: sources which produce photons more efficiently yield correspondingly higher signals. The introduction of sub-grid gas clumping in the radiative transfer simulations produces significantly more power at small scales, and more non-Gaussian features, but has little effect at large scales. The patchy nature of the reionization process roughly doubles the total observed kSZ signal for l~3000-10^4 compared to non-patchy scenarios with the same total electron-scattering optical depth.

At the beginning of the reionization epoch, radiation sources produce fluctuations in the redshifted 21-cm background. We show that different types of sources (such as miniquasars, Pop II and III stars, supernovae, etc.) produce distinct signatures in the 21-cm signal radial profiles and statistical fluctuations, through which they can be identified. Further, we show that the 21-cm signal from X-ray emitting sources is much easier to observe than was expected, due to a previously neglected pumping mechanism.

We present detailed predictions for the redshifted 21cm signal from the epoch of reionization. These predictions are obtained from radiative transfer calculations on the results of large scale (100/h Mpc), high dynamic range, cosmological simulations. We consider several scenarios for the reionization history, of both early and extended reionization. From the simulations we construct and analyze a range of observational characteristics, from the global signal, via detailed images and spectra, to statistical representations of rms fluctuations, angular power spectra, and probability distribution functions to characterize the non-gaussianity of the 21cm signal. (abbreviated abstract)

We present the first calculation of the kinetic Sunyaev-Zel'dovich (kSZ) effect due to the inhomogeneus reionization of the universe based on detailed large-scale radiative transfer simulations of reionization. The resulting sky power spectra peak at l=2000-8000 with maximum values of l(l+1)C_l/(2pi)_{max}~4-7x10^{-13}. The scale roughly corresponds to the typical ionized bubble sizes observed in our simulations, of ~5-20 Mpc. The kSZ anisotropy signal from reionization dominates the primary CMB signal above l=3000. At large scales the patchy kSZ signal depends only on the source efficiencies. It is higher when sources are more efficient at producing ionizing photons, since such sources produce larger ionized regions, on average, than less efficient sources. The introduction of sub-grid gas clumping in the radiative transfer simulations produce significantly more power at small scales, but has little effect at large scales. The patchy reionization kSZ signal is dominated by the post-reionization signal from fully-ionized gas, but the two contributions are of similar order at scales l~3000-10^4, indicating that the kSZ anisotropies from reionization are an important component of the total kSZ signal at these scales.

We present new analytic calculations of the coupling between ultraviolet resonance photons and the population of the hyperfine states in the light elements (H, D, He3+) which include several previously neglected physical processes. Among these are the backreaction of resonant scattering on the pumping radiation, the scattering of Ly_beta photons and the effect of local departure from pure Hubble flow. The application of the new treatment to the redshifted hydrogen 21 and deuterium 92 cm lines from the high-redshift universe results in an amplitude correction of up to an order of magnitude. We further show that the standard assumption that ultraviolet pumping drives the spin temperature towards the kinetic temperature does not hold for deuterium, whose spin temperature is generally negative.

New results on the anisotropy of the cosmic microwave background (CMB) and its polarization based on the first 3 years of data from the Wilkinson Microwave Anisotropy Probe (WMAP) have revised the electron scattering optical depth downward from tau_es=0.17+0.08-0.07 to tau_es=0.09+/-0.03. This implies a shift of the effective reionization redshift from z_r~17 to z_r~11. Previous attempts to explain the high redshift of reionization inferred from the WMAP 1 year data have led to widespread speculation that the sources of reionization must have been much more efficient than those associated with the star formation observed at low redshift. This is consistent, for example, with the suggestion that early star formation involved massive, Population III stars that early on produced most of the ionizing radiation escaping from halos. It is therefore tempting to interpret the new WMAP results as implying that we can now relax those previous high demands on the efficiency of the sources of reionization and perhaps even turn the argument around as evidence against such high efficiency. We show that this is not the case, however. The new WMAP results also find that the primordial density fluctuation power spectrum has a lower amplitude, sigma_8, and departs substantially from the scale-invariant spectrum. We show that these effects combine to cancel the impact of the later reionization implied by the new value of tau_es on the required ionizing efficiency per collapsed baryon. The delay of reionization is surprisingly well matched by a comparable delay (by a factor of ~1.4 in scale factor) in the formation of the halos responsible for reionization.

We present the first large-scale radiative transfer simulations of cosmic reionization, in a simulation volume of (100/h Mpc)^3, while at the same time capturing the dwarf galaxies which are primarily responsible for reionization. We achieve this by combining the results from extremely large, cosmological, N-body simulations with a new, fast and efficient code for 3D radiative transfer, C^2-Ray. The resulting electron-scattering optical depth is in good agreement with the first-year WMAP polarization data. We show that reionization clearly proceeded in an inside-out fashion, with the high-density regions being ionized earlier, on average, than the voids. Ionization histories of smaller-size (5 to 10 comoving Mpc) subregions exibit a large scatter about the mean and do not describe the global reionization history well. The minimum reliable volume size for such predictions is ~30 Mpc. We derive the power-spectra of the neutral, ionized and total gas density fields and show that there is a significant boost of the density fluctuations in both the neutral and the ionized components relative to the total at arcminute and larger scales. We find two populations of HII regions according to their size, numerous, mid-sized (~10 Mpc) regions and a few, rare, very large regions tens of Mpc in size. We derive the statistical distributions of the ionized fraction and ionized gas density at various scales and for the first time show that both distributions are clearly non-Gaussian. (abridged)

We derive the equations for the propagation of relativistic ionization fronts in both static and moving gases. We focus on the supersonic R-type phase that occurs right after a source turns on, and we compare the nonrelativistic and relativistic solutions for several important cases. Relativistic corrections can be significant up until the light-crossing time of the equilibrium Stromgren sphere. For a static medium, we obtain exact analytical solutions and apply them to the illustrative problems of an O star in a molecular cloud and a starburst in a high-redshift cosmological halo. Relativistic corrections can be important at early times when the H II regions are small, as well as at later times, if a density gradient causes the I-front to accelerate. For the cosmologically-expanding IGM, we derive an analytical solution in the case of a steady source and a constant clumping factor. Here relativistic corrections are significant for short-lived, highly-luminous sources like QSOs at the end of reionization, but negligible for weaker or higher-redshift sources. Finally, we numerically calculate the evolution of relativistic I-fronts in the presence of small-scale structure and infall, for a large galaxy undergoing a starburst and a luminous, high-redshift QSO. For such strong and short-lived sources, the relativistic corrections are quite significant, and small-scale structure can decrease the size of the H II region by up to an additional ~25%. (abridged)

We present numerical radiation-hydrodynamic simulations of the evolution of HII regions formed in an inhomogeneous medium resulting from turbulence simulations. We find that the filamentary structure of the underlying density distribution produces a highly irregular shape for the ionized region, in which the ionization front escapes to large distances in some directions within 80,000 years. In other directions, on the other hand, neutral gas in the form of dense globules persists within 1 parsec of the central star for the full duration of our simulation (400,000 years). Divergent photoablation flows from these globules maintain a root-mean-squared velocity in the ionized gas that is close to the ionized sound speed. Simulated images in optical emission lines show morphologies that are in strikingly detailed agreement with those observed in real HII regions.

We show that the epoch(s) of reionization when the ionization fraction of the universe is about half can be determined by correlating Cosmic Microwave Background (CMB) temperature maps with 21-cm line maps at degree scales ($l\sim 100$). During reionization peculiar motion of free electrons induces the Doppler anisotropy of the CMB, while density fluctuations of neutral hydrogen induce the 21-cm line anisotropy. In our simplified model of inhomogeneous reionization, a positive correlation arises as the universe reionizes whereas a negative correlation arises as it recombines; thus, the sign of the correlation provides information on the reionization history which cannot be obtained by presently. The signal comes mainly from large scales (k~0.01 Mpc^-1) where linear perturbation theory is valid and complexity due to patchy reionization is averaged out. Since the Doppler signal comes from ionized regions and the 21-cm comes from neutral ones, the correlation has a well-defined peak(s) in redshift when the ionization fraction of the universe is about half. Furthermore, the cross-correlation is much less sensitive to systematic errors, especially foreground emission, than the auto-correlation of 21-cm lines: this is analogous to the temperature-polarization correlation of the CMB being more immune to systematic errors than the polarization-polarization. Therefore, we argue that the Doppler-21cm correlation provides a robust measurement of the 21-cm anisotropy, which can also be used as a diagnostic tool for detected signals in the 21-cm data -- detection of the cross-correlation provides the strongest confirmation that the signal is of cosmological origin. We show that the Square Kilometer Array can easily measure the predicted correlation signal for 1~year of survey observation.

The H atoms inside minihalos (i.e. halos with virial temperatures T_vir < 10^4 K, in the mass range roughly from 10^4 M_sun to 10^8 M_sun) during the cosmic dark ages in a LambdaCDM universe produce a redshifted background of collisionally-pumped 21-cm line radiation which can be seen in emission relative to the cosmic microwave background (CMB). Previously, we used semi-analytical calculations of the 21-cm signal from individual halos of different mass and redshift and the evolving mass function of minihalos to predict the mean brightness temperature of this 21-cm background and its angular fluctuations. Here we use high-resolution cosmological N-body and hydrodynamic simulations of structure formation at high redshift (z>8) to compute the mean brightness temperature of this background from both minihalos and the intergalactic medium (IGM) prior to the onset of Ly-alpha radiative pumping. We find that the 21-cm signal from gas in collapsed, virialized minihalos dominates over that from the diffuse shocked gas in the IGM.

Radiative transfer simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come to the reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray-tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non-trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on 5 test problems: (0) basic physics, (1) isothermal H II region expansion and (2) H II region expansion with evolving temperature, (3) I-front trapping and shadowing by a dense clump, (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analyzed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multi-frequency treatment approach, resulting in different thicknesses of the ionized-neutral transition regions and different temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To this aim all test inputs and outputs are made publicly available in digital form.

We present a series of simulations of the fragmentation of a molecular cloud, leading to the formation of a cluster of protostellar cores. The purpose of these simulations is to address a specific numerical problem called artificial fragmentation, that plagues SPH simulations of cloud fragmentation. We argue that this is a serious problem that needs to be addressed, and that the only reasonable and practical way to address it is to use a relatively new technique called particle splitting. Our largest simulation has an effective resolution of 256^3 particles (much higher than most previous SPH simulations of cloud fragmentation) and results in the formation of a dense cluster containing ~3000 protostellar cores. It is the first simulation of this kind to properly resolve the Jeans mass throughout the entire system, at all times, thus preventing artificial fragmentation.

One of the most debated issues in the theoretical modeling of cosmic reionization is the impact of small-mass gravitationally-bound structures. We carry out the first numerical investigation of the role of such sterile `minihaloes', which serve as self-shielding screens of ionizing photons. Minihaloes are too small to be properly resolved in current large-scale cosmological simulations, and thus we estimate their effects using a sub-grid model, considering two cases that bracket their effect within this framework. In the `extreme suppression' case in which minihalo formation ceases once a region is partially ionized, their effect on cosmic reionization is modest, reducing the volume-averaged ionization fraction by an overall factor of less than 15%. In the other extreme, in which minihalo formation is never suppressed, they delay complete reionization as much as Delta z~2, in rough agreement with the results from a previous semi-analytical study by the authors. Thus, depending on the details of the minihalo formation process, their effect on the overall progress of reionization can range from modest to significant, but the minihalo photon consumption is by itself insufficient to force an extended reionization epoch.

Simulations predict that the first stars in an LCDM universe formed at redshifts z>20 in minihalos with masses of about 10^6 M_sun. We have studied their radiative feedback by simulating the propagation of ionization fronts (I-fronts) created by these first Population III stars (M_* = 15-500 M_sun) at z=20, within the density field of a cosmological simulation of primordial star formation, outward thru the host minihalo and into the surrounding gas. A three-dimensional ray-tracing calculation tracks the I-front once the H II region evolves a "champagne flow'' inside the minihalo, after the early D-type I-front detaches from the shock and runs ahead, becoming R-type. We take account of the hydrodynamical back-reaction by an approximate model of the central wind. We find that the escape fraction of ionizing radiation from the host halo increases with stellar mass, with 0.7 < f_esc <0.9 for 80 < M_*/M_sun < 500. To quantify the ionizing efficiency of these stars as they begin cosmic reionization, we find that, for M_* > 80 M_sun, the ratio of gas mass ionized to stellar mass is ~60,000, roughly half the number of ionizing photons released per stellar baryon. Nearby minihalos are shown to trap the I-front, so their centers remain neutral. This is contrary to the recent suggestion that these stars would trigger formation of a second generation by fully ionizing neighboring minihalos, stimulating H_2 formation in their cores. Finally, we discuss how the evacuation of gas from the host halo reduces the growth and luminosity of "miniquasars'' that may form from black hole remnants of the first stars.

N-body simulations find a universal structure for the halos which result from the nonlinear growth of Gaussian-random-noise density fluctuations in the CDM universe. This talk summarized our attempts to derive and explain this universal structure by analytical approximation and simplified models. As an example, we show here that a 1D spherical infall model involving a fluid approximation derived from the Boltzmann equation can explain not only the halo density profile but its phase-space density profile, as well.

We present a new numerical method for calculating the transfer of ionizing radiation, called C^2-Ray=Conservative, Causal Ray-tracing method. The method is explicitly photon-conserving, so the depletion of ionizing photons by bound-free opacity is guaranteed to equal the photoionizations these photons caused. As a result, grid cells can be large and very optically-thick without loss of accuracy. The method also uses an analytical relaxation solution for the ionization rate equations for each time step which can accommodate time steps which greatly exceed the characteristic ionization and ionization front crossing times. Together, these features make it possible to integrate the equation of transfer along a ray with many fewer cells and time steps than previous methods. For multi-dimensional calculations, the code utilizes short-characteristics ray tracing. C^2-Ray is well-suited for coupling radiative transfer to gas and N-body dynamics methods, on both fixed and adaptive grids, without imposing additional limitations on the time step and grid spacing. We present several tests of the code involving propagation of ionization fronts in one and three dimensions, in both homogeneous and inhomogeneous density fields. We compare to analytical solutions for the ionization front position and velocity, some of which we derive here for the first time.

We have derived the first, fully-cosmological, similarity solutions for CDM halo formation in the presence of nongravitational collisionality, which provides an analytical theory of the effect of the self-interacting dark matter (SIDM) hypothesis on halo density profiles. Collisions transport heat inward, which produces a constant-density core, while continuous infall pumps energy into the halo to stabilize the core against gravothermal catastrophe. These solutions improve upon earlier attempts to model the formation and evolution of SIDM halos, offer deeper insight than existing N-body experiments, and yield a more precise determination of the dependence of halo density profile on the value of the CDM self-interaction cross section. Different solutions arise for different values of the dimensionless collisionality parameter Q = s rho_b r_v \~ r_v/l_mfp, where s is the scattering cross section per unit mass, rho_b is the cosmic mean matter density, r_v is halo virial radius and l_mfp is the collision mean free path. The maximum flattening of central density occurs for an intermediate value of Q, Q_th, at which the halo is maximally relaxed to isothermality. The density profiles with constant-density cores preferred by dwarf and LSB rotation curves are best fit by the maximally-flattened (Q=Q_th) solution. If we assume that dwarfs and LSB galaxies formed at their typical collapse epoch in LCDM, then the value of s which makes Q=Q_th is s ~ 200 cm^{2}/g, much higher than previous estimates, s ~ 0.5-5 cm^{2}/g, based on N-body experiments. If s is independent of collision velocity, then the same value s ~ 200 cm^{2}/g would make Q>Q_th for clusters, which typically formed only recently, resulting in relatively less flattening of their central density profile and a smaller core.

The H atoms inside minihalos (i.e. halos with virial temperatures T_vir < 10^4 K, in the mass range roughly from 10^4 M_sun to 10^8 M_sun) during the cosmic dark ages in a LCDM universe produce a redshifted background of collisionally-pumped 21-cm line radiation which can be seen in emission relative to the cosmic microwave background (CMB). Previously, we used semi-analytical calculations of the 21-cm signal from individual halos of different mass and redshift and the evolving mass function of minihalos to predict the mean brightness temperature of this 21-cm background and its angular fluctuations. Here we use high-resolution cosmological N-body and hydrodynamic simulations of structure formation at high redshift (z > 8) to compute the mean brightness temperature of this background from both minihalos and the intergalactic medium (IGM) prior to the onset of Ly-alpha radiative pumping. We find that the 21-cm signal from gas in collapsed, virialized minihalos dominates over that from the diffuse shocked gas in the IGM.

The propagation of cosmological ionization fronts (I-fronts) during reionization is strongly influenced by small-scale structure. Here we summarize our recent attemps to understand the effect of this small-scale structure. We present high resolution cosmological N-body simulations at high-z (z>6) which resolve a wide range of halo mass, from minihalos to clusters of large, rare halos. We also study how minihalos affect I-fronts, through simulations of minihalo photoevaporation. Furthermore, we modify the I-front propagation equations to account for evolving small-scale structure, and incorporate these results into a semi-analytical reionization model. When intergalactic medium clumping and minihalo clustering around sources are included, small-scale structure affects reionization by slowing it down and extending it in time. This helps to explain observations by the Wilkinson Microwave Anisotropy Probe, which imply an early and extended reionization epoch. We also study how source clustering affects the evolution and size of H II regions, finding, in agreement with simulations, that H II regions usually expand, rarely shrinking. Hence, "relic H II regions'' are an exception, rather than the rule. When the suppression of small-mass sources in already-ionized regions by Jeans-mass filtering is accounted for, H II regions are smaller, delaying overlap. We also present a new numerical method for radiative transfer which is fast, efficient, and easily coupled to hydrodynamics and N-body codes, along with sample tests and applications.

The propagation of cosmological ionization fronts during the reionization of the universe is strongly influenced by small-scale gas inhomogeneities due to structure formation. These inhomogeneities include both collapsed minihalos, which are generally self-shielding, and lower-density structures, which are not. The minihalos are dense and sufficiently optically-thick to trap intergalactic ionization fronts, blocking their path and robbing them of ionizing photons until the minihalo gas is expelled as an evaporative wind. The lower-density structures do not trap these fronts, but they can slow them down by increasing the overall recombination rate in the intergalactic medium. In this paper we study the effects of both types of inhomogeneities, including nonlinear clustering effects, and we find that both IGM clumping and collapsed minihalos have significant yet qualitatively different impacts on reionization. While the number density of minihalos on average increases strongly with time, the density of minihalos inside H II regions around ionizing sources is largely constant. Thus the impact of minihalos is essentially to decrease the number of ionizing photons available to the IGM at all epochs, which is equivalent to a reduction in the luminosity of each source. On the other hand, the effect of IGM clumping increases strongly with time, slowing down reionization and extending it. Thus while the impact of minihalos is largely degenerate with the unknown source efficiency, IGM clumping can help significantly in reconciling the recent observations of cosmic microwave background polarization with quasar absorption spectra at z~6, which together point to an early but extended reionization epoch.

Dark-matter halos are the scaffolding around which galaxies and clusters are built. They form when the gravitational instability of primordial density fluctuations causes regions which are denser than average to slow their cosmic expansion, recollapse, and virialize. Understanding the equilibrium structure of these halos is thus a prerequisite for understanding galaxy and cluster formation. Numerical N-body simulations of structure formation from Gaussian-random-noise initial conditions in the CDM universe find a universal internal structure for halos. Objects as different in size and mass as dwarf spheroidal galaxies and galaxy clusters are predicted to have halos with the same basic structure when properly rescaled, independent of halo mass, of the shape of the power spectrum of primordial density fluctuations, and of the cosmological background parameters. This remarkable universality is a fundamental prediction of the CDM model, but our knowledge is limited to the ``empirical'' N-body simulation results, with little analytical understanding. We summarize here our attempts to fill this gap, in an effort to derive and give physical insight to the numerical results and extend them beyond the range of numerical simulation: (1) Simulated halos which form from highly simplified initial conditions involving gravitational instability in a cosmological pancake show that many of the universal properties of CDM halos are generic to cosmological gravitational collapse and do not require Gaussian-random-noise density fluctuations or hierarchical clustering. (2) A fluid approximation derived from the Boltzmann equation yields an analytical theory of halo dynamics which can explain many of the N-body results ... (abridged)

The weak, R-type ionization fronts (I-fronts) which swept across the intergalactic medium (IGM) during the reionization of the universe often found their paths blocked by cosmological minihaloes. When this happened, the neutral gas which filled each minihalo was photoevaporated; as the I-front burned its way through the halo, decelerating from R-type to D-type, all the gas was blown back into the IGM as an ionized, supersonic wind. In a previous paper (Shapiro, Iliev and Raga 2004), we described this process and presented our results of the first simulations of it by numerical gas dynamics with radiation transport in detail. For illustration we focused on the particular case of a 10^7 solar masses minihalo overrun at z=9 by an intergalactic I-front caused by a distant source of ionizing radiation, for different source spectra (either stellar from massive Pop. II or III stars, or QSO-like) and a flux level typical of that expected during reionization. In a LambdaCDM universe, minihaloes formed in abundance before and during reionization and, thus, their photoevaporation is an important, possibly dominant, feature of reionization, which slowed it down and wasted ionizing photons. We have now performed a larger set of high-resolution simulations to determine and quantify the dependence of minihalo photoevaporation times and photon consumption rates on halo mass, redshift, ionizing flux level and spectrum. We find that the average number of ionizing photons each minihalo atom absorbs during its photoevaporation is typically in the range 2-10. For the collapsed fraction in minihaloes expected during reionization, this can add about 1 photon per total atom to the requirements for completing reionization, potentially doubling the minimum number of photons required to reionize the universe.

We present the first gas dynamical simulations of the photoevaporation of cosmological minihalos overtaken by the ionization fronts which swept through the IGM during reionization in a LCDM universe, including the effects of radiative transfer. We demonstrate the phenomenon of I-front trapping inside minihalos, in which the weak, R-type fronts which traveled supersonically across the IGM decelerated when they encountered the dense, neutral gas inside minihalos, becoming D-type I-fronts, preceded by shock waves. For a minihalo with virial temperature T_vir < 10^4 K, the I-front gradually burned its way through the minihalo which trapped it, removing all of its baryonic gas by causing a supersonic, evaporative wind to blow backwards into the IGM, away from the exposed layers of minihalo gas just behind the advancing I-front. Such hitherto neglected feedback effects were widespread during reionization. N-body simulations and analytical estimates of halo formation suggest that sub-kpc minihalos such as these, with T_vir < 10^4 K, were so common as to cover the sky around larger-mass source halos and possibly dominate the absorption of ionizing photons. This means that previous estimates of the number of ionizing photons per H atom required to complete reionization which neglected this effect may be too low. Regardless of their effect on the progress of reionization, however, the minihalos were so abundant that random lines of sight thru the high-z universe should encounter many of them, which suggests that it may be possible to observe the processes described here in the absorption spectra of distant sources.

The gravitational lensing properties of cosmological halos depend upon the mass distribution within each halo. The description of halos as nonsingular, truncated isothermal spheres, a particular solution of the isothermal Lane-Emden equation (suitably modified for Lambda nonzero), has proved to be a useful approximation for the halos which form from realistic initial conditions in a CDM universe. The nonsingular TIS model reproduces many of the quantitative features of the N-body results for CDM halos, except in the very center, where CDM N-body halos show density profiles which vary as rho ~ r^(-alpha), alpha>1, instead of a small flat core. Possible discrepancies between these cuspy halo predictions of the CDM N-body simulations and observations of the inner mass profiles of dwarf and LSB disk galaxies based upon their rotation curves and of clusters based upon strong lensing measurements have led to a search for other diagnostics. A description of the lensing by TIS halos would be useful in this regard, as a self-consistent model for CDM halos in a proper cosmological context, nonsingular but otherwise consistent with the CDM N-body results. We derive here the basic lensing properties of individual TIS halos. For comparison, we also consider three singular profiles: the Navarro-Frenk-White density profile, the singular isothermal sphere, and the Schwarzschild lens. For all profiles, we compute the critical curves and caustics, the image separation, the magnification and brightness ratio, the shear, the time delay, and the average shear inside the tidal radius. This provides tools for studying the statistical properties of lensing by TIS and other lenses in the context of a theory of cosmological structure formation.

We study the formation and evolution of self-interacting dark matter (SIDM) halos. We find analytical, fully cosmological similarity solutions for their dynamics, which take proper account of the collisional interaction of SIDM particles, based on a fluid approximation derived from the Boltzmann equation. These similarity solutions are relevant to galactic and cluster halo formation in the cold dark matter (CDM) model. Different solutions arise for different values of the dimensionless collisionality parameter, Q=sigma*rho_b*r_s, where sigma is the SIDM particle scattering cross section per unit mass, rho_b is the cosmic mean density, and r_s is the shock radius. For all solutions, a flat-density, isothermal core is present which grows in size as a fixed fraction of r_s, pumped by cosmological infall. Accordingly, core collapse must in general be delayed until infall becomes negligible, contrary to previous analyses based on isolated halos, which predict core collapse in a Hubble time. Our solutions agree with N-body simulations, which match observed galactic rotation curves if sigma=[0.56 - 5.6]cm^2/g, implying Q=[6.2e-7 - 3.6e-5]. Similar profiles also arise for Q=[1.37e-2 - 1.7e-1], or sigma=[1.2e4 - 2.7e4]cm^2/g, a regime not previously simulated.

We use the universal mass accretion history recently reported for simulations of halo formation in the cold dark matter model (CDM) to analyze the formation and growth of a single halo. We derive the time-dependent density profile three different ways, based upon three approximations of successively greater realism: equilibrium, radial orbits, and a fluid approximation. For the equilibrium model, the density profile is well-fit by either an NFW or Moore profile over a limited range of radii and scale factors. For the radial orbit model, we find profiles which are generally steeper than the NFW profile, with an inner logarithmic slope approaching -2, consistent with a purely radial collisionless system. In the fluid approximation, we find good agreement with the NFW and Moore profiles for radii resolved by N-body simulations (r/r_200 > 0.01), and an evolution of concentration parameter nearly identical to that found in N-body simulations. The evolving structure of cosmological halos is therefore best understood as the effect of a time-varying rate of mass infall on a smoothly distributed, isotropic, collisionless fluid.

We present the first gas dynamical simulations of the photoevaporation of cosmological minihalos overtaken by the ionization fronts which swept through the IGM during reionization in a LCDM universe, including the effects of radiative transfer. We demonstrate the phenomenon of I-front trapping inside minihalos, in which the weak, R-type fronts which traveled supersonically across the IGM decelerated when they encountered the dense, neutral gas inside minihalos, becoming D-type I-fronts, preceded by shock waves. For a minihalo with virial temperature T_vir < 10,000 K, the I-front gradually burned its way through the minihalo which trapped it, removing all of its baryonic gas by causing a supersonic, evaporative wind to blow backwards into the IGM, away from the exposed layers of minihalo gas just behind the advancing I-front. Such hitherto neglected feedback effects were widespread during reionization. N-body simulations and analytical estimates of halo formation suggest that sub-kpc minihalos such as these, with T_vir < 10,000 K, were so common as to dominate the absorption of ionizing photons. This means that previous estimates of the number of ionizing photons per H atom required to complete reionization which neglected this effect may be too low. Regardless of their effect on the progress of reionization, however, the minihalos were so abundant that random lines of sight thru the high-z universe should encounter many of them, which suggests that it may be possible to observe the processes described here in the absorption spectra of distant sources.

In the standard Cold Dark Matter (CDM) theory of structure formation, virialized minihalos (with T_vir = 10,000 K or less) form in abundance at high redshift (z > 6), during the cosmic "dark ages." The hydrogen in these minihalos, the first nonlinear baryonic structures to form in the universe, is mostly neutral and sufficiently hot and dense to emit strongly at the 21-cm line. We calculate the emission from individual minihalos and the radiation background contributed by their combined effect. Minihalos create a "21-cm forest" of emission lines. We predict that the angular fluctuations in this 21-cm background should be detectable with the planned LOFAR and SKA radio arrays, thus providing a direct probe of structure formation during the "dark ages." Such a detection will serve to confirm the basic CDM paradigm while constraining the background cosmology parameters, the shape of the power-spectrum of primordial density fluctuations, the onset and duration of the reionization epoch, and the conditions which led to the first stars and quasars. We present results here for the currently-favored, flat LCDM model, for different tilts of the primordial power spectrum. These minihalos will also cause a "21-cm forest" of absorption lines, as well, in the spectrum of radio continuum sources at high redshift, if the latter came into existence before the end of reionization.

Hydrogen atoms inside virialized minihalos generate a radiation background from redshifted 21-cm line emission whose angular fluctuations reflect clustering during before and during reionization. We have shown elsewhere that this emission may be detectable with the planned Low Frequency Array (LOFAR) and Square Kilometer Array (SKA) in a flat LambdaCDM. This is a direct probe of structure during the ``Dark Ages'' at redshifts z>6 and down to smaller scales than have previously been constrained. In our original calculation, we used a standard approximation known as the ``linear bias'' [e.g. Mo & White (1996)]. Here we improve upon that treatment by considering the effect of nonlinear clustering. To accomplish this, we develop a new analytical method for calculating the nonlinear Eulerian bias of halos, which should be useful for other applications as well. Predictions of this method are compared with the results of LambdaCDM N-body simulations, showing significantly better agreement than the standard linear bias approximation. When applied to the 21-cm background from minihalos, our formalism predicts fluctuations that differ from our original predictions by up to 30% at low frequencies and small scales. However, within the range of frequencies and angular scales at which the signal could be observable by LOFAR and SKA as currently planned, the differences are small and our original predictions prove robust. Our results indicate that while a smaller frequency bandwidth of observation leads to a higher signal that is more sensitive to nonlinear effects, this effect is counteracted by the lowered sensitivity of the radio arrays. We calculate the best frequency bandwidth for these observations to be Delta nu~2 MHz.(abbreviated)

We study the formation and evolution of self-interacting dark matter (SIDM) halos. We find analytical, fully cosmological similarity solutions taking account of the collisional interaction of SIDM particles. This interaction results in a thermal conductivity that heats the halo core and flattens its density profile. These similarity solutions are relevant to galactic and cluster halo formation in the CDM model. We assume an initial mass profile dM/M M^{-eps}, as in the familiar secondary infall model. If eps=1/6, SIDM halos will evolve self-similarly, with a cold, supersonic infall terminated by a strong accretion shock. Different solutions arise for different values of the collisionality parameter, Q= sigma rho_b r_s, where sigma is the scattering cross section, rho_b is the cosmic mean density, and r_s is the shock radius. For all these solutions, a flat-density, isothermal core is present which grows in size as a fixed fraction of r_s. We find two different regimes for these solutions: 1) for Q \leq Q_{th}, the core density decreases and core size increases as Q increases; 2) for Q \geq Q_{th}, the core density increases and core size decreases as Q increases. Our similarity solutions are in agreement with previous N-body simulations of SIDM halos, which correspond to the low-Q regime, if Q=[8.4e-4 - 4.9e-2]Q_{th} (low-Q), or sigma=[0.56-5.6]cm^2/g. As Q=\infty, our similarity solution aquires a central density cusp, in agreement with some simulation results which used an ordinary collisional fluid to approximate the effects of SIDM collisionality. When Q=[18.6-231]Q_{th} or sigma=[1.2e4 - 2.71e4]cm^2/g, for which we find flat-density cores comparable to those of the observationally acceptable low-Q solutions, has not previously been identified. Further study of this regime is warranted.

Adaptive SPH and N-body simulations were carried out to study the collapse and evolution of dark matter halos that result from the gravitational instability and fragmentation of cosmological pancakes. Such halos resemble those formed by hierarchical clustering in a CDM universe and serve as a convenient test-bed model for studying halo dynamics. Our halos are in approximate virial equilibrium and roughly isothermal, as in CDM simulations. Their density profiles agree quite well with the fit to N-body results for CDM halos by Navarro, Frenk, & White (NFW). This test-bed model enables us to study the evolution of individual halos. The masses of our halos evolve in three stages: an initial collapse, continuous infall, and a final stage in which infall tapers off as a result of finite mass supply. In the continuous infall stage, halo mass grows at the rate expected for self-similar spherical infall, with M(a) proportional to the scale factor a. After the end of initial collapse at a=a_0, the concentration parameter grows linearly with a, c(a)~4a/a_0. The virial ratio 2T/|W| just after virialization is about 1.35, as predicted by the truncated isothermal sphere model and consistent with the value expected for a virialized halo in which mass infall contributes an effective surface pressure. Thereafter, the virial ratio evolves towards the value expected for an isolated halo, 2T/|W|~1. This mass accretion history and evolution of concentration parameter are very similar to those reported recently in N-body simulations of CDM. We therefore conclude that the fundamental properties of halo formation and evolution are generic to the formation of cosmological halos by gravitational instability and are not limited to hierarchical collapse scenarios or even to Gaussian-random-noise initial conditions.

The gravitational lensing properties of cosmological halos depend upon the mass distribution within each halo. The description of halos as nonsingular, truncated isothermal spheres, a particular solution of the isothermal Lane- Emden equation (suitably modified for Lambda non equal 0), has proved to be a useful approximation for the halos which form from realistic initial conditions in a CDM universe. We derive here the basic lensing properties of such halos, including the image separation, magnification, shear, and time-delay. We also provide analytical expressions for the critical curves and caustics. We show how the scale-free results we derive yield scale-dependent lensing properties which depend upon the cosmological background universe and the mass and collapse redshift of the lensing halos, according to the truncated isothermal sphere (TIS) model of CDM halos derived elsewhere. We briefly describe the application of these results to the currently-favored Lambda-CDM universe.

In the standard Cold Dark Matter (CDM) theory of structure formation, virialized minihalos (with T_vir}<10^4 K) form in abundance at high redshift (z>6), during the cosmic ``dark ages.'' The hydrogen in these minihalos, the first nonlinear baryonic structures to form in the universe, is mostly neutral and sufficiently hot and dense to emit strongly at the 21-cm line. We calculate the emission from individual minihalos and the radiation background contributed by their combined effect. Minihalos create a ``21-cm forest'' of emission lines. We predict that the angular fluctuations in this 21-cm background should be detectable with the planned LOFAR and SKA radio arrays, thus providing a direct probe of structure formation during the ``dark ages.'' Such a detection will serve to confirm the basic CDM paradigm while constraining the shape of the power-spectrum of primordial density fluctuations down to much smaller scales than have previously been constrained, the onset and duration of the reionization epoch, and the conditions which led to the first stars and quasars. We present results here for the currently-favored, flat $\Lambda$CDM model, for different tilts of the primordial power spectrum.

Current data suggest that the central mass densities $\rho_0$ and phase-space densities $Q\equiv\rho_0/\sigma_V^3$ of cosmological halos in the present universe are correlated with their velocity dispersions $\sigma_V$ over a very wide range of $\sigma_V$ from less than 10 to more than 1000 $\rm km s^{-1}$. Such correlations are an expected consequence of the statistical correlation of the formation epochs of virialized objects in the CDM model with their masses; the smaller-mass halos typically form first and merge to form larger-mass halos later. We have derived the $Q-\sigma_V$ and $\rho_0-\sigma_V$ correlations for different CDM cosmologies and compared the predicted correlations with the observed properties of a sample of low-redshift halos ranging in size from dwarf spheroidal galaxies to galaxy clusters. Our predictions are generally consistent with the data, with preference for the currently-favored, flat $\Lambda$CDM model. Such a comparison serves to test the basic CDM paradigm while constraining the background cosmology and the power-spectrum of primordial density fluctuations, including larger wavenumbers than have previously been constrained.

Dark-matter halos are the scaffolding around which galaxies and clusters are built. They form when the gravitational instability of primordial density fluctuations causes regions which are denser than average to slow their cosmic expansion, recollapse, and virialize. Objects as different in size and mass as dwarf spheroidal galaxies and galaxy clusters are predicted by the CDM model to have halos with a universal, self-similar equilibrium structure whose parameters are determined by the halo's total mass and collapse redshift. These latter two are statistically correlated, however, since halos of the same mass form on average at the same epoch, with small-mass objects forming first and then merging hierarchically. The structural properties of dark-matter dominated halos of different masses, therefore, should reflect this statistical correlation, an imprint of the statistical properties of the primordial density fluctuations which formed them. Current data reveal these correlations, providing a fundamental test of the CDM model which probes the shape of the power spectrum of primordial density fluctuations and the cosmological background parameters.

We summarize the ideas that led to the Adaptive Smoothed Particle Hydrodynamics (ASPH) algorithm, with anisotropic smoothing and shock-tracking. We then identify a serious new problem for SPH simulations with shocks and radiative cooling --- false cooling --- and discuss a possible solution based on the shock-tracking ability of ASPH.

The universe was reionized by redshift z ~ 6 by a small fraction of the baryons in the universe, which released energy following their condensation out of a cold, dark, and neutral IGM into the earliest galaxies. The theory of this reionization is a critical missing link in the theory of galaxy formation. Its numerous observable consequences include effects on the spectrum, anisotropy and polarization of the cosmic microwave background and signatures of high-redshift star and quasar formation. This energy release also created feedback on galaxy formation which left its imprint on the mass spectrum and internal characteristics of galaxies and on the gas between galaxies long after reionization was complete. Recent work suggests that the photoevaporation of dwarf galaxy minihalos may have consumed most of the photons required to reionize the currently-favored Lambda-CDM universe. We will review recent developments in our understanding of this process.

We consider an explosion at the center of a halo which forms at the intersection of filaments inside a cosmological pancake, a convenient test-bed model for galaxy formation. ASPH/P3M simulations reveal that such explosions are anisotropic. The energy and metals are channeled into the low density regions, away from the pancake. The pancake remains essentially undisturbed, even if the explosion is strong enough to blow away all the gas located inside the halo and reheat the IGM surrounding the pancake. Infall quickly replenishes this ejected gas and gradually restores the gas fraction as the halo continues to grow. Estimates of the collapse epoch and SN energy-release for galaxies of different mass in the CDM model can relate these results to scale-dependent questions of blow-out and blow-away and their implication for early IGM heating and metal enrichment and the creation of gas-poor dwarf galaxies.

We have derived an analytical model for the postcollapse equilibrium structure of cosmological halos as nonsingular truncated isothermal spheres (TIS) and compared this model with observations and simulations of cosmological halos on all scales. Our model is in good agreement with the observations of the internal structure of dark-matter-dominated halos from dwarf galaxies to X-ray clusters. It reproduces many of the average properties of halos in CDM simulations to good accuracy, including the density profiles outside the central region, while avoiding the possible discrepancy at small radii between observed galaxy and cluster density profiles and the singular density profiles predicted by N-body simulations of the CDM model. While much attention has been focused lately on this possible discrepancy, we show that the observed galaxy rotation curves and correlations of halo properties nevertheless contain valuable additional information with which to test the theory, despite this uncertainty at small radii. The available data allows us to constrain the fundamental cosmological parameters and also to put a unique constraint on the primordial density fluctuation power spectrum at large wavenumbers (i.e. small mass scale).

Adaptive SPH and N-body simulations were carried out to study the evolution of the equilibrium structure of dark matter halos that result from the gravitational instability and fragmentation of cosmological pancakes. Such halos resemble those formed by hierarchical clustering from realistic initial conditions in a CDM universe and, therefore, serve as a test-bed model for studying halo dynamics. The dark matter density profile is close to the universal halo profile identified previously from N-body simulations of structure formation in CDM, with a total mass and concentration parameter which grow linearly with scale factor a. When gas is included, this concentration parameter is slightly larger than the pure N-body result. We also find that the dark matter velocity distribution is less isotropic and more radial than found by N-body simulations of CDM.

An analytical model is presented for the postcollapse equilibrium structure of virialized objects which condense out of a low-density cosmological background universe, either matter-dominated or flat with a cosmological constant. This generalizes the model we derived previously for an Einstein-de Sitter (EdS) universe. The model is based upon the assumption that cosmological haloes form from the collapse and virialization of top-hat density perturbations and are spherical, isotropic, and isothermal. This leads to the prediction of a unique, nonsingular, truncated isothermal sphere (TIS), a particular solution of the Lane-Emden equation (modified for nonzero cosmological constant). The size and virial temperature are unique functions of the mass and redshift of formation of the object for a given background universe. The central density is roughly proportional to the critical density of the universe at the epoch of collapse. This TIS model is in good agreement with observations of the internal structure of dark matter--dominated haloes on scales ranging from dwarf galaxies to X-ray clusters. It also reproduces many of the average properties of haloes in simulations of the Cold Dark Matter (CDM) model to good accuracy, suggesting that it is a useful analytical approximation for haloes which form from realistic initial conditions. Our TIS model matches the density profiles of haloes in CDM N-body simulations outside the innermost region, while avoiding the steep central cusp of the latter which is in apparent conflict with observations. The TIS model may also be relevant to nonstandard CDM models, like self-interacting dark matter, recently proposed to resolve this conflict.

Rotation curves of dark-matter-dominated galaxies measure the mass profiles of galactic halos and thereby test theories of their cosmological origin. While attention has focused lately on the possible discrepancy at small galactocentric radii between observed rotation curves and the singular density profiles predicted by N-body simulations of the Cold Dark Matter (CDM) model, the observed rotation curves nevertheless contain valuable additional information with which to test the theory and constrain the fundamental cosmological parameters, despite this uncertainty at small radii. An analytical model we derived elsewhere for the postcollapse equilibrium of cosmological halos as truncated, nonsingular, isothermal spheres (TIS) reproduces many of the average properties of halos in CDM simulations to good accuracy, including the density profiles outside the central region. The circular velocity profile of this TIS model is, moreover, in excellent agreement with the observed ones and yields the mass and formation epoch of an observed halo from the parameters of its rotation curve. This allows us to predict correlations amongst rotation curve parameters, such as the maximum velocity and the radius at which it occurs, for different mass halos forming at different epochs in the CDM model. As an example, we derive the observed v_max-r_max relation analytically, with preference for the flat Lambda-CDM model.

Observations of a flat density profile in the cores of dark-matter-dominated halos on the two extremes of mass for virialized objects in the universe, dwarf galaxies and galaxy clusters, present a serious challenge to the current standard theory of structure formation involving Cold Dark Matter (CDM). By contrast, N-body simulations of halo formation in the latter indicate density profiles which are singular and steeply rising towards the center. A flat-density core on the cluster scale is indicated by gravitational lensing observations, most significantly by the strong-lensing measurements of CL 0024+1654 by the Hubble Space Telescope. A recent re-analysis of this cluster has suggested that a uniform-density core is not demanded by the data, thereby eliminating a significant piece of the conflict between the observations and the CDM theoretical predictions. We show here, however, that the singular mass profile which that analysis reports as consistent with the lensing measurements of CL 0024+1654 implies a velocity dispersion which is much higher than the measured value for this cluster.

The first sources of ionizing radiation to condense out of the dark and neutral IGM sent ionization fronts sweeping outward through their surroundings, overtaking other condensed objects and photoevaporating them. This feedback of universal reionization on cosmic structure formation is demonstrated here by gas dynamical simulations, including radiative transfer, for a cosmological minihalo of dark matter and baryons exposed to an external source of ionizing radiation, either starlight or quasar light, just after the passage of the global ionization front created by the source.

As an idealized model of the effects of energy release by supernovae during galaxy formation, we consider an explosion at the center of a halo which forms at the intersection of filaments in the plane of a cosmological pancake by gravitational instability during pancake collapse. Such halos resemble the virialized objects found in N-body simulations in a CDM universe and, therefore, serve as a convenient, scale-free test-bed model for galaxy formation. ASPH/P3M simulations reveal that such explosions are anisotropic. The energy and metals are channeled into the low density regions, away from the pancake plane. The pancake remains essentially undisturbed, even if the explosion is strong enough to blow away all the gas located inside the halo at the onset of the explosion and reheat the IGM surrounding the pancake. Infall quickly replenishes this ejected gas and gradually restores the gas fraction as the halo mass continues to grow. Estimates of the collapse epoch and SN energy-release for galaxies of different mass in the CDM model can relate these results to scale-dependent questions of blow-out and blow-away and their implication for early IGM heating and metal enrichment and the creation of dark-matter-dominated dwarf galaxies.

Adaptive SPH and N-body simulations were carried out to study the effect of gasdynamics on the structure of dark matter halos that result from the gravitational instability and fragmentation of cosmological pancakes. Such halos resemble those formed in a hierarchically clustering CDM universe and serve as a test-bed model for studying halo dynamics. With no gas, the density profile is close to the universal profile identified previously from N-body simulations of structure formation in CDM. When gas is included, the gas in the halo is approximately isothermal, and both the dark matter and the gas have singular central density profiles which are steeper than that of the dark matter with no gas. This worsens the disagreement between observations of constant density cores in cosmological halos and the singular ones found in simulations. We also find that the dark matter velocity distribution is less isotropic than found by N-body simulations of CDM, because of the strongly filamentary substructure.

An analytical model for the postcollapse equilibrium structure of virialized objects which condense out of the cosmological background universe is described and compared with observations and simulations of cosmological halos. The model is based upon the assumption that virialized halos are isothermal, which leads to a prediction of a unique nonsingular isothermal sphere for the equilibrium structure, with a core density which is proportional to the mean background density at the epoch of collapse. These predicted nonsingular isothermal spheres are in good agreement with observations of the internal structure of dark-matter-dominated halos from dwarf galaxies to X-ray clusters. Our model also reproduces many of the average properties of halos in CDM simulations to good accuracy, suggesting that it is a useful analytical approximation for halos which form from realistic initial conditions. While CDM N-body simulations find profiles with a central cusp, our nonsingular model matches the simulated halos outside the innermost region well. This model may also be of interest as a description of halos in nonstandard CDM models like self-interacting dark matter, which have been proposed to eliminate the discrepancy between the cuspy halos of standard CDM simulations and observed halos with uniform-density cores.

The first sources of ionizing radiation to condense out of the dark and neutral IGM sent ionization fronts sweeping outward through their surroundings, overtaking other condensed objects and photoevaporating them. This feedback of universal reionization on cosmic structure formation is demonstrated here by gas dynamical simulations, including radiative transfer, for the case of a cosmological minihalo of dark matter and baryons exposed to an external source of ionizing starlight, just after the passage of the global ionization front created by the source.

A new model for the postcollapse equilibrium structure of virialized objects which condense out of the cosmological background universe is described and compared with observations and simulations of cosmological halos from dwarf galaxies to X-ray clusters. The model is based upon the assumption that virialized halos are isothermal, which leads to a prediction of a unique nonsingular isothermal sphere for the equilibrium structure, with a core radius which is approximately 1/30 times the size and a core density which is proportional to the mean background density at the epoch of collapse. These predicted nonsingular isothermal spheres are in good agreement with the observations of the internal structure of dark-matter-dominated halos from dwarf galaxies to X-ray clusters.

The first sources of ionizing radiation to condense out of the dark and neutral IGM sent ionization fronts sweeping outward through their surroundings, overtaking other condensed objects and photoevaporating them. This feedback effect of universal reionization on cosmic structure formation is demonstrated here for the case of a cosmological minihalo of dark matter and baryons exposed to an external source of ionizing radiation with a quasar-like spectrum, just after the passage of the global ionization front created by the source. We model the pre-ionization minihalo as a truncated, nonsingular isothermal sphere in hydrostatic equilibrium following its collapse out of the expanding background universe and virialization. Results are presented of the first gas dynamical simulations of this process, including radiative transfer. A sample of observational diagnostics is also presented, including the spatially-varying ionization levels of C, N, and O in the flow if a trace of heavy elements is present and the integrated column densities of H I, He I and II, and C IV thru the photoevaporating gas at different velocities which would be measured in absorption against a background source like that responsible for the ionization.

When density fluctuations collapse gravitationally out of the expanding cosmological background universe to form galaxies, the secondary energy release which results can affect their subsequent evolution profoundly. We focus here on the effects of one form of such energy release - explosions, such as might result from the supernovae which end the lives of the first generation of massive stars to form inside protogalaxies. We are particularly interested in the consequences of the nonspherical geometry and continuous infall which are characteristic of galaxy formation from realistic initial and boundary conditions. As an idealized model which serves to illustrate and quantify the importance of these effects, we study the effect of explosions on the quasi-spherical objects which form at the intersections of filaments in the plane of a cosmological pancake, as a result of gravitational instability and fragmentation of the pancake. We study the formation and evolution of these "galaxies," subject to the explosive injection of energy at their centers, by numerical gas dynamical simulation in 3D utilizing our new, anisotropic version of Smoothed Particle Hydrodynamics, Adaptive SPH ("ASPH"), with a P3M gravity solver.

We calculate the maximum fraction of matter which is able to condense out of the expanding background universe, for any universe which will expand forever. We use a simple spherical model for the growth of density fluctuations in the universe. This model includes open, matter-dominated universe and universes in which there is an uniform background component (e.g. the cosmological constant or "quintessence"). In these background universes, Omega eventually drops significantly below unity. When this happens, gravitational instability is suppressed and, with it, so is the growth of the collapsed fraction. We identify a limitation of the Press-Schechter (PS) approximation. In this approximation, the mass function determined from the predicted collapse of overdense regions is multiplied by a correction factor of 2, assuming that each bound fluctuation will accrete an equal share of mass from nearby underdense regions. Our model determines the actual value of the correction factor. We show that, while the factor of 2 adopted by the PS approximation is correct for an Einstein-de Sitter universe, it is not correct when the freeze-out of fluctuation growth inherent in the more general class of background universes described above occurs. When freeze-out occurs, the correction factor reduces to unity and the PS approximation must overestimate the collapsed fraction. We apply our model to open CDM and flat Lambda-CDM models. For H0=70 and Omega0=0.3, these models yield asymptotic collapsed fractions of 0.0361 and 0.0562, respectively, on the galaxy cluster mass-scale, only 55% of the values determined by the PS approximation. These results have implications for the use of the latter approximation to compare the observed space density of X-ray clusters today with that predicted by cosmological models.

The postcollapse structure of objects which form by gravitational condensation out of the expanding cosmological background universe is a key element in the theory of galaxy formation. Towards this end, we have reconsidered the outcome of the nonlinear growth of a uniform, spherical density perturbation in an unperturbed background universe - the cosmological ``top-hat'' problem. We adopt the usual assumption that the collapse to infinite density at a finite time predicted by the top-hat solution is interrupted by a rapid virialization caused by the growth of small-scale inhomogeneities in the initial perturbation. We replace the standard description of the postcollapse object as a uniform sphere in virial equilibrium by a more self-consistent one as a truncated, nonsingular, isothermal sphere in virial and hydrostatic equilibrium, including for the first time a proper treatment of the finite-pressure boundary condition on the sphere. The results differ significantly from both the uniform sphere and the singular isothermal sphere approximations for the postcollapse objects. These results will have a significant effect on a wide range of applications of the Press-Schechter and other semi-analytical models to cosmology. The truncated isothermal sphere solution presented here predicts the virial temperature and integrated mass distribution of the X-ray clusters formed in the CDM model as found by detailed, 3D, numerical gas and N-body dynamical simulations remarkably well. This solution allows us to derive analytically the numerically-calibrated mass-temperature and radius-temperature scaling laws for X-ray clusters which were derived empirically by Evrard, Metzler and Navarro from simulation results for the CDM model. (Shortened)

The first sources of ionizing radiation to condense out of the dark and neutral IGM sent ionization fronts sweeping outward through their surroundings, overtaking other primordial gas-clouds and photoevaporating them. Results are presented of the first gas dynamical simulations of this process, including radiative transfer, along with some observational diagnostics.

This paper presents an alternative formulation of the ASPH algorithm for evolving anisotropic smoothing kernels, in which the geometric approach of Shapiro et al. (1996; Paper I) is replaced by an approach involving a local transformation of coordinates to those in which the underlying anisotropic volume changes appear to be isotropic. The ASPH method is presented in 2D and 3D, including a number of details not previously included in Paper I, some of which represent either advances or different choices with respect to Paper I. Among the advances included here are an asynchronous time-integration scheme with different time steps for different particles and the generalization of the ASPH method to 3D. The shock-tracking algorithm described in Paper I for locally adapting the artificial viscosity to restrict viscous heating just to particles encountering shocks, is not included here. Instead, we adopt a different interpolation kernel for use with the artificial viscosity, which has the effect of spatially localizing effects of the artificial viscosity. This version of the ASPH method in 2D and 3D is then applied to a series of 1D, 2D, and 3D test problems, and the results are compared to those of standard SPH applied to the same problems. These include the problem of cosmological pancake collapse, the Riemann shock tube, cylindrical and spherical Sedov blast waves, the collision of two strong shocks, and problems involving shearing disks intended to test the angular momentum conservation properties of the method. These results further support the idea that ASPH has significantly better resolving power than standard SPH for a wide range of problems, including that of cosmological structure formation. (Abridged)

We present a set of cosmological variables, called "supercomoving variables," which are particularly useful for describing the gas dynamics of cosmic structure formation. For ideal gas with gamma=5/3, the supercomoving position, velocity, density, temperature, and pressure are constant in time in a uniform, isotropic, adiabatically expanding universe. Expressed in terms of these supercomoving variables, the cosmological fluid conservation equations and the Poisson equation closely resemble their noncosmological counterparts. This makes it possible to generalize noncosmological results and techniques to cosmological problems, for a wide range of cosmological models. These variables were initially introduced by Shandarin for matter-dominated models only. We generalize supercomoving variables to models with a uniform component corresponding to a nonzero cosmological constant, domain walls, cosmic strings, a nonclumping form of nonrelativistic matter (e.g. massive nettrinos), or radiation. Each model is characterized by the value of the density parameter Omega0 of the nonrelativistic matter component in which density fluctuation is possible, and the density parameter OmegaX of the additional, nonclumping component. For each type of nonclumping background, we identify FAMILIES within which different values of Omega0 and OmegaX lead to fluid equations and solutions in supercomoving variables which are independent of Omega0 and OmegaX. We also include the effects of heating, radiative cooling, thermal conduction, viscosity, and magnetic fields. As an illustration, we describe 3 familiar cosmological problems in supercomoving variables: the growth of linear density fluctuations, the nonlinear collapse of a 1D plane-wave density fluctuation leading to pancake formation, and the Zel'dovich approximation.

Energy released by a small fraction of the baryons in the universe, which condensed out of while the IGM was cold, dark, and neutral,reheated and reionized it, exposing gas clouds within it to the glare of ionizing radiation. The first gas dynamical simulations of the photoevaporation of an intergalactic cloud by a quasar, including radiative transfer, are presented, along with a few observational diagnostics.

In theories in which the cosmological constant Lambda takes a variety of values in different ``subuniverses,'' the probability distribution of its observed values is conditioned by the requirement that there must be someone to measure it. This probability is proportional to the fraction of matter which is destined to condense out of the background into mass concentrations large enough to form observers. We calculate this ``collapsed fraction'' by a simple, pressure-free, spherically symmetric, nonlinear model for the growth of density fluctuations in a flat universe with arbitrary value of the cosmological constant, applied in a statistical way to the observed spectrum of density fluctuations at recombination. From this, the probability distribution for the vacuum energy density rho_V=Lambda/8pi G for Gaussian random density fluctuations is derived analytically. It is shown that the results depend on only one quantity, sigma^3 RHO, where sigma^2 and RHO are the variance and mean value of the fluctuating matter density field at recombination, respectively. To calculate sigma, we adopt the flat CDM model with nonzero cosmological constant and fix the amplitude and shape of the primordial power spectrum in accordance with data on cosmic microwave background anisotropy from the COBE satellite DMR experiment. A comparison of the results of this calculation of the likely values of rho_V with present observational bounds on the cosmological constant indicates that the small, positive value of rho_V (up to 3 times greater than the present cosmic mass density) suggested recently by several lines of evidence is a reasonably likely value to observe, even if all values of rho_V are equally likely a priori.

The gravitational instability of cosmological pancakes composed of collisionless dark matter in an Einstein-de Sitter universe is investigated numerically to demonstrate that pancakes are unstable with respect to fragmentation and the formation of filaments. A ``pancake'' is defined here as the nonlinear outcome of the growth of a 1D, sinusoidal, plane-wave, adiabatic density perturbation. We have used high resolution, 2D, N-body simulations by the Particle-Mesh (PM) method to study the response of pancakes to perturbation by either symmetric (density) or antisymmetric (bending or rippling) modes, with corresponding wavevectors k_s and k_a transverse to the wavevector k_p of the unperturbed pancake plane-wave. We consider dark matter which is initially ``cold'' (i.e. with no random thermal velocity in the initial conditions). We also investigate the effect of a finite, random, isotropic, initial velocity dispersion (i.e. initial thermal velocity) on the fate of pancake collapse and instability. Pancakes are shown to be gravitationally unstable with respect to all perturbations of wavelength l<l_p (where l_p= 2pi/k_p). These results are in contradiction with the expectations of an approximate, thin-sheet energy argument.