The 21-cm signal of neutral hydrogen is a sensitive probe of the Epoch of Reionization (EoR) and Cosmic Dawn. Currently operating radio telescopes have ushered in a data-driven era of 21-cm cosmology, providing the first constraints on the astrophysical properties of sources that drive this signal. However, extracting astrophysical information from the data is highly non-trivial and requires the rapid generation of theoretical templates over a wide range of astrophysical parameters. To this end emulators are often employed, with previous efforts focused on predicting the power spectrum. In this work we introduce 21cmGEM - the first emulator of the global 21-cm signal from Cosmic Dawn and the EoR. The smoothness of the output signal is guaranteed by design. We train neural networks to predict the cosmological signal using a database of ~30,000 simulated signals which were created by varying seven astrophysical parameters: the star formation efficiency and the minimal mass of star-forming halos; the efficiency of the first X-ray sources and their spectrum parameterized by spectral index and the low energy cutoff; the mean free path of ionizing photons and the CMB optical depth. We test the performance with a set of ~2,000 simulated signals, showing that the relative error in the prediction has an r.m.s. of 0.0159. The algorithm is efficient, with a running time per parameter set of 0.16 sec. Finally, we use the database of models to check the robustness of relations between the features of the global signal and the astrophysical parameters that we previously reported.

The ARCADE2 and LWA1 experiments have claimed an excess over the Cosmic Microwave Background (CMB) at low radio frequencies. If the cosmological high-redshift contribution to this radio background is between 0.1% and 22% of the CMB at $1.42\,$GHz, it could explain the tentative EDGES Low-Band detection of the anomalously deep absorption in the 21-cm signal of neutral hydrogen (Fialkov & Barkana 2019). We use the upper limit on the 21-cm signal from the Epoch of Reionization ($z=9.1$) based on $141\,$hours of observations with LOFAR (Mertens et al. 2020) to evaluate the contribution of the high redshift Universe to the detected radio background. Marginalizing over astrophysical properties of star-forming halos, we find (at 68% C.L.) that the cosmological radio background can be at most 0.7% of the CMB at $1.42\,$GHz. This limit rules out strong contribution of the high-redshift Universe to the ARCADE2 and LWA1 measurements. Even though LOFAR places limit on the extra radio background, excess of $0.1-0.7$% over the CMB (at $1.42\,$GHz) is still allowed and could explain the EDGES Low-Band detection. We also constrain the thermal and ionization state of the gas at $z = 9.1$, and put limits on the properties of the first star-forming objects. We find that, in agreement with the limits from EDGES High-Band data, LOFAR data disfavour scenarios with inefficient X-ray sources and cases where the Universe was ionized by massive halos only.

The early epoch in which the first stars and galaxies formed is among the most exciting unexplored eras of the Universe. A major research effort focuses on probing this era with the 21-cm spectral line of hydrogen. While most research focused on statistics like the 21-cm power spectrum or the sky-averaged global signal, there are other ways to analyze tomographic 21-cm maps, which may lead to novel insights. We suggest statistics based on quantiles as a method to probe non-Gaussianities of the 21-cm signal. We show that they can be used in particular to probe the variance, skewness, and kurtosis of the temperature distribution, but are more flexible and robust than these standard statistics. We test these statistics on a range of possible astrophysical models, including different galactic halo masses, star-formation efficiencies, and spectra of the X-ray heating sources, plus an exotic model with an excess early radio background. Simulating data with angular resolution and thermal noise as expected for the Square Kilometre Array (SKA), we conclude that these statistics can be measured out to redshifts above 20 and offer a promising statistical method for probing early cosmic history.

The multi-tracer technique employs a ratio of densities of two differently biased galaxy samples that trace the same underlying matter density field, and was proposed to alleviate the cosmic variance problem. Here we propose a novel application of this approach, applying it to two different tracers one of which is the 21-cm signal of neutral hydrogen from the epochs of reionization and comic dawn. The second tracer is assumed to be a sample of high-redshift galaxies, but the approach can be generalized and applied to other high-redshift tracers. We show that the anisotropy of the ratio of the two density fields can be used to measure the sky-averaged 21-cm signal, probe the spectral energy distribution of radiative sources that drive this signal, and extract large-scale properties of the second tracer, e.g., the galaxy bias. Using simulated 21-cm maps and mock galaxy samples, we find that the method works well for an idealized galaxy survey. However, in the case of a more realistic galaxy survey which only probes highly biased luminous galaxies, the inevitable Poisson noise makes the reconstruction far more challenging. This difficulty can be mitigated with the greater sensitivity of future telescopes along with larger survey volumes.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged]

The recent tentative detection by the EDGES Low-Band antenna of the hydrogen 21-cm line from cosmic dawn, if confirmed, is the first ever signature observed from the epoch of primordial star formation. However, the magnitude and the shape of this signal are incompatible with standard astrophysical predictions, requiring either colder than expected gas, or an excess radio background above the Cosmic Microwave Background (CMB) radiation. In this paper we explore the latter option, investigating possible 21-cm signals in models with a phenomenological excess radio background together with the standard astrophysical modeling. Assuming uniform radiation with a synchrotron-like spectrum existing throughout cosmic history, we show that such a radio background, in addition to explaining the EDGES Low-Band signal, enhances the power spectrum. The signal during cosmic dawn and reionization strongly depends on both the intensity of the radio background and the astrophysical parameters. We verify the broad agreement of our models with the detected feature, including additional constraints from the EDGES High-Band data, high-redshift quasars, and observational limits on the extragalactic radio background. The data imply a lower limit on the star formation efficiency of 2.8\%, an upper limit on the minimum mass of star-forming halos of M$_{\rm h}\sim 10^9$ M$_\odot$ at $z = 17$, and a lower limit on the excess radio background of 1.9 times the CMB at the 78 MHz cosmic dawn frequency (corresponding to 0.1\% of the CMB at 1.42 GHz). The properties of X-ray sources remain unconstrained. We also show that during the dark ages, such a radio background saturates the 21-cm signal, imprinting a unique signature in both the global signal and the power spectrum.

We present new constraints on parameters of cosmic dawn and the epoch of reionization derived from the EDGES High-Band spectrum ($90-190$ MHz). The parameters are probed by evaluating global $21$ cm signals generated with the recently developed Global21cm tool. This tool uses neural networks trained and tested on $\sim 30,000$ spectra produced with semi-numerical simulations that assume the standard thermal evolution of the cosmic microwave background and the intergalactic medium. From our analysis, we constrain at $68\%$ (1) the minimum virial circular velocity of star-forming halos to $V_{\rm c}<19.3$ km s$^{-1}$, (2) the X-ray heating efficiency of early sources to $f_{\rm X}>0.0042$, and (3) the low-energy cutoff of the X-ray spectral energy distribution to $\nu_{\rm min}<2.3$ keV. We also constrain the star-formation efficiency ($f_*$), the electron scattering optical depth ($\tau_{\rm e}$), and the mean-free path of ionizing photons ($R_{\rm mfp}$). We re-compute the constraints after incorporating into the analysis four estimates for the neutral hydrogen fraction from high-$z$ quasars and galaxies, and a prior on $\tau_{\rm e}$ from Planck $2018$. The largest impact of the external observations is on the parameters that most directly characterize reionization. Specifically, we derive the combined $68\%$ constraints $\tau_{\rm e}<0.063$ and $R_{\rm mfp}>27.5$ Mpc. The external observations also have a significant effect on $V_{\rm c}$ due to its degeneracy with $\tau_{\rm e}$, while the constraints on $f_*$, $f_{\rm X}$, and $\nu_{\rm min}$, remain primarily determined by EDGES.

We perform an improved cosmic microwave background (CMB) analysis to search for dark matter-proton scattering with a momentum-transfer cross section of the form $\sigma_0 v^n$ for $n=-2$ and $n=-4$. In particular, we present a new and robust prescription for incorporating the relative bulk velocity between the dark matter and baryon fluids into the standard linear Boltzmann calculation. Using an iterative procedure, we self-consistently include the effects of the bulk velocities in a cosmology in which dark matter interacts with baryons. With this prescription, we derive CMB bounds on the cross section, excluding $\sigma_0 > 2.3 \times 10^{-33}~\mathrm{cm}^2$ for $n=-2$ and $\sigma_0 > 1.7 \times 10^{-41}~\mathrm{cm}^2$ for $n=-4$ at $95\%$ confidence, for dark matter masses below 10 MeV. Furthermore, we investigate how these constraints change when only a subcomponent of dark matter is interacting. We show that Planck limits vanish if $\lesssim 0.4\%$ of dark matter is tightly coupled to baryons. We discuss the implications of our results for present and future cosmological observations.

We investigate the hypothesis that Coulomb-type interactions between dark matter (DM) and baryons explain the anomalously low 21cm brightness-temperature minimum at redshift z ~ 17 that was recently measured by the EDGES experiment. In particular, we reassess the validity of the scenario where a small fraction of the total DM is millicharged, focusing on newly derived constraints from Planck 2015 cosmic microwave background (CMB) data. Crucially, the CMB power spectrum is sensitive to DM-baryon scattering if the fraction of interacting DM is larger than (or comparable to) the fractional uncertainty in the baryon energy density. Meanwhile, there is a mass-dependent lower limit on the fraction for which the required interaction to cool the baryons sufficiently is so strong that it drives the interacting-DM temperature to the baryon temperature prior to their decoupling from the CMB. If this occurs as early as recombination, the cooling saturates. We precisely determine the viable parameter space for millicharged DM, and find that only a fraction (m_chi/MeV) 0.0115% <~ f <~ 0.4% of the entire DM content, and only for DM-particle masses between 0.5 MeV - 35 MeV, can be charged at the level needed to marginally explain the anomaly, without violating limits from SLAC, CMB, Big-Bang nucleosynthesis (BBN), or stellar and SN1987A cooling. In reality, though, we demonstrate that at least moderate fine tuning is required to both agree with the measured absorption profile and overcome various astrophysical sources of heating. Finally, we point out that a ~0.4% millicharged DM component which is tightly coupled to the baryons at recombination may resolve the current 2-sigma tension between the BBN and CMB determinations of the baryon energy density. Future CMB-S4 measurements will be able to probe this scenario directly.

Recently the EDGES collaboration reported an anomalous absorption signal in the sky-averaged 21-cm spectrum around $z=17$. Such a signal may be understood as an indication for an unexpected cooling of the hydrogen gas during or prior to the so called Cosmic Dawn era. Here we explore the possibility that dark matter cooled the gas through velocity-dependent, Rutherford-like interactions. We argue that such interactions require a light mediator that is highly constrained by 5th force experiments and limits from stellar cooling. Consequently, only a hidden or the visible photon can in principle mediate such a force. Neutral hydrogen thus plays a sub-leading role and the cooling occurs via the residual free electrons and protons. We find that these two scenarios are strongly constrained by the predicted dark matter self-interactions and by limits on millicharged dark matter respectively. We conclude that the 21-cm absorption line is unlikely to be the result of gas cooling via the scattering with a dominant component of the dark matter. An order 1\% subcomponent of millicharged dark matter remains a viable explanation.

The cosmic radio-frequency spectrum is expected to show a strong absorption signal corresponding to the 21-centimetre-wavelength transition of atomic hydrogen around redshift 20, which arises from Lyman-alpha radiation from some of the earliest stars. By observing this 21-centimetre signal - either its sky-averaged spectrum or maps of its fluctuations, obtained using radio interferometers - we can obtain information about cosmic dawn, the era when the first astrophysical sources of light were formed. The recent detection of the global 21-centimetre spectrum reveals a stronger absorption than the maximum predicted by existing models, at a confidence level of 3.8 standard deviations. Here we report that this absorption can be explained by the combination of radiation from the first stars and excess cooling of the cosmic gas induced by its interaction with dark matter. Our analysis indicates that the spatial fluctuations of the 21-centimetre signal at cosmic dawn could be an order of magnitude larger than previously expected and that the dark-matter particle is no heavier than several proton masses, well below the commonly predicted mass of weakly interacting massive particles. Our analysis also confirms that dark matter is highly non-relativistic and at least moderately cold, and primordial velocities predicted by models of warm dark matter are potentially detectable. These results indicate that 21-centimetre cosmology can be used as a dark-matter probe.

The recent detection of an anomalously strong 21-cm signal of neutral hydrogen from Cosmic Dawn by the EDGES Low-Band radio experiment can be explained if cold dark matter particles scattered off the baryons draining excess energy from the gas. In this Letter we explore the expanded range of the 21-cm signal that is opened up by this interaction, varying the astrophysical parameters as well as the properties of dark matter particles in the widest possible range. We identify models consistent with current data by comparing to both the detection in the Low-Band and the upper limits from the EDGES High-Band antenna. We find that consistent models predict a 21-cm fluctuation during Cosmic Dawn that is between 3 and 30 times larger than the largest previously expected without dark matter scattering. The expected power spectrum exhibits strong Baryon Acoustic Oscillations imprinted by the velocity-dependent cross-section. The latter signature is a smoking gun of the velocity-dependent scattering and could be used by interferometers to verify the dark matter explanation of the EDGES detection.

Baryon Acoustic Oscillations (BAOs) in the early Universe are predicted to leave an as yet undetected signature on the relative clustering of total mass versus luminous matter. This signature, a modulation of the relative large-scale clustering of baryons and dark matter, offers a new angle to compare the large scale distribution of light versus mass. A detection of this effect would provide an important confirmation of the standard cosmological paradigm and constrain alternatives to dark matter as well as non-standard fluctuations such as Compensated Isocurvature Perturbations (CIPs). The first attempt to measure this effect in the SDSS-III BOSS Data Release 10 CMASS sample remained inconclusive but allowed to develop a method, which we detail here and use to conduct the second observational search. When using the same model as in our previous study and including CIPs in the model, the DR12 data are consistent with a null-detection, a result in tension with the strong evidence previously measured with the DR10 data. This tension remains when we use a more realistic model taking into account our knowledge of the survey flux limit, as the data then privilege a zero effect. In the absence of CIPs, we obtain a null detection consistent with both the absence of the effect and the amplitude predicted in previous theoretical studies. This shows the necessity of more accurate data in order to prove or disprove the theoretical predictions.

Spectral distortions in the cosmic microwave background over the 40--200~MHz band are imprinted by neutral hydrogen in the intergalactic medium prior to the end of reionization. This signal, produced in the redshift range $z = 6-34$ at the rest frame wavelength of 21 cm, has not been detected yet; and poor understanding of high redshift astrophysics results in a large uncertainty in the expected spectrum. The SARAS~2 radiometer was purposely designed to detect the sky-averaged 21-cm signal. The instrument, deployed at the Timbaktu Collective (Southern India) in April--June 2017, collected 63~hr of science data, which were examined for the presence of the cosmological 21-cm signal. In our previous work the first-light data from SARAS~2 radiometer were analyzed with Bayesian likelihood-ratio tests using $264$ plausible astrophysical scenarios. In this paper we re-examine the data using an improved analysis based on the frequentist approach and forward modeling. We show that SARAS~2 data rejects 27 models, out of which 25 are rejected at a significance $>5\sigma$. All the rejected models share the scenario of inefficient heating of the primordial gas by the first population of X-ray sources along with rapid reionization.

The high-redshift 21-cm signal of neutral hydrogen is expected to be observed within the next decade and will reveal epochs of cosmic evolution that have been previously inaccessible. Due to the lack of observations, many of the astrophysical processes that took place at early times are poorly constrained. In recent work we explored the astrophysical parameter space and the resulting large variety of possible global (sky-averaged) 21-cm signals. Here we extend our analysis to the fluctuations in the 21-cm signal, accounting for those introduced by density and velocity, Ly$\alpha$ radiation, X-ray heating, and ionization. While the radiation sources are usually highlighted, we find that in many cases the density fluctuations play a significant role at intermediate redshifts. Using both the power spectrum and its slope, we show that properties of high-redshift sources can be extracted from the observable features of the fluctuation pattern. For instance, the peak amplitude of ionization fluctuations can be used to estimate whether heating occurred early or late and, in the early case, to also deduce the cosmic mean ionized fraction at that time. The slope of the power spectrum has a more universal redshift evolution than the power spectrum itself and can thus be used more easily as a tracer of high-redshift astrophysics. Its peaks can be used, for example, to estimate the redshift of the Ly$\alpha$ coupling transition and the redshift of the heating transition (and the mean gas temperature at that time). We also show that a tight correlation is predicted between features of the power spectrum and of the global signal, potentially yielding important consistency checks.

The early star-forming Universe is still poorly constrained, with the properties of high-redshift stars, the first heating sources, and reionization highly uncertain. This leaves observers planning 21-cm experiments with little theoretical guidance. In this work we explore the possible range of high-redshift parameters including the star formation efficiency and the minimal mass of star-forming halos; the efficiency, spectral energy distribution, and redshift evolution of the first X-ray sources; and the history of reionization. These parameters are only weakly constrained by available observations, mainly the optical depth to the cosmic microwave background. We use realistic semi-numerical simulations to produce the global 21-cm signal over the redshift range $z = 6-40$ for each of 193 different combinations of the astrophysical parameters spanning the allowed range. We show that the expected signal fills a large parameter space, but with a fixed general shape for the global 21-cm curve. Even with our wide selection of models we still find clear correlations between the key features of the global 21-cm signal and underlying astrophysical properties of the high redshift Universe, namely the Ly$\alpha$ intensity, the X-ray heating rate, and the production rate of ionizing photons. These correlations can be used to directly link future measurements of the global 21-cm signal to astrophysical quantities in a mostly model-independent way. We identify additional correlations that can be used as consistency checks.

Long wavelength spectral distortions in the Cosmic Microwave Background arising from the 21-cm transition in neutral Hydrogen are a key probe of Cosmic Dawn and the Epoch of Reionization. These features may reveal the nature of the first stars and ultra-faint galaxies that transformed the spin temperature and ionization state of the primordial gas. SARAS~2 is a spectral radiometer purposely designed for precision measurement of these monopole or all-sky global 21-cm spectral distortions. We use 63~hr night time observing of the radio background in the frequency band 110-200~MHz with the radiometer deployed at the Timbaktu Collective in Southern India to derive likelihoods for plausible redshifted 21-cm signals predicted by theoretical models. First light with SARAS 2 disfavors the class of models that feature weak X-ray heating (with $f_X \leq 0.1$) and rapid reionization (with peak $\frac{dT_b}{dz} \geq 120~\textrm{mK per unit redshift interval}$ ).

We use the observed unresolved cosmic X-ray background (CXRB) in the 0.5-2 keV band and existing upper limits on the 21-cm power spectrum to constrain the high-redshift population of X-ray sources, focusing on their effect on the thermal history of the Universe and the cosmic 21-cm signal. Because the properties of these sources are poorly constrained, we consider hot gas, X-ray binaries and mini-quasars (i.e., sources with soft or hard X-ray spectra) as possible candidates. We find that (1) the soft-band CXRB sets an upper limit on the X-ray efficiency of sources that existed before the end of reionization, which is one-to-two orders of magnitude higher than typically assumed efficiencies, (2) hard sources are more effective in generating the CXRB than the soft ones, (3) the commonly-assumed limit of saturated heating is not valid during the first half of reionization in the case of hard sources, with any allowed value of X-ray efficiency, (4) the maximal allowed X-ray efficiency sets a lower limit on the depth of the absorption trough in the global 21-cm signal and an upper limit on the height of the emission peak, while in the 21-cm power spectrum it sets a minimum amplitude and frequency for the high-redshift peaks, and (5) the existing upper limit on the 21-cm power spectrum sets a lower limit on the X-ray efficiency for each model. When combined with the 21-cm global signal, the CXRB will be useful for breaking degeneracies and helping constrain the nature of high-redshift heating sources.

Understanding the formation and evolution of the first stars and galaxies represents one of the most exciting frontiers in astronomy. Since the universe was filled with hydrogen atoms at early times, the most promising probe of the epoch of the first stars is the prominent 21-cm spectral line of hydrogen. Current observational efforts are focused on the cosmic reionization era, but observations of the pre-reionization cosmic dawn are also promising. While observationally unexplored, theoretical studies predict a rich variety of observational signatures from cosmic dawn. As the first stars formed, their radiation (plus that from stellar remnants) produced significant cosmic events including Lyman-alpha coupling at z~25, and early X-ray heating. Much focus has gone to studying the angle-averaged power spectrum of 21-cm fluctuations. Additional probes include the global (sky-averaged) 21-cm spectrum, and the line-of-sight anisotropy of the 21-cm power spectrum. A particularly striking signature may result from the recently recognized effect of a supersonic relative velocity between the dark matter and gas. Work in this field, focused on understanding the whole era of reionization and cosmic dawn with analytical models and numerical simulations, is likely to grow in intensity and importance, as the theoretical predictions are finally expected to confront 21-cm observations in the coming years. [Abridged]

Studies of the formation of the first stars have established that they formed in small halos of $\sim 10^5 - 10^6 M_{\odot}$ via molecular hydrogen cooling. Since a low level of ultraviolet radiation from stars suffices to dissociate molecular hydrogen, under the usually-assumed scenario this primordial mode of star formation ended by redshift $z \sim 15$ and much more massive halos came to dominate star formation. However, metal enrichment from the first stars may have allowed the smaller halos to continue to form stars. In this Letter we explore the possible effect of star formation in metal-rich low-mass halos on the redshifted 21-cm signal of neutral hydrogen from $z = 6-40$. These halos are significantly affected by the supersonic streaming velocity, with its characteristic baryon acoustic oscillation (BAO) signature. Thus, enrichment of low-mass galaxies can produce a strong signature in the 21-cm power spectrum over a wide range of redshifts, especially if star formation in the small halos was more efficient than suggested by current simulations. We show that upcoming radio telescopes can easily distinguish among various possible scenarios.

Baryon Acoustic Oscillations (BAOs) in the early Universe are predicted to leave an as yet undetected signature on the relative clustering of total mass versus luminous matter. A detection of this effect would provide an important confirmation of the standard cosmological paradigm and constrain alternatives to dark matter as well as non-standard fluctuations such as Compensated Isocurvature Perturbations (CIPs). We conduct the first observational search for this effect, by comparing the number-weighted and luminosity-weighted correlation functions, using the SDSS-III BOSS Data Release 10 CMASS sample. When including CIPs in our model, we formally obtain evidence at $3.2\sigma$ of the relative clustering signature and a limit that matches the existing upper limits on the amplitude of CIPs. However, various tests suggest that these results are not yet robust, perhaps due to systematic biases in the data. The method developed in this Letter, used with more accurate future data such as that from DESI, is likely to confirm or disprove our preliminary evidence.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

The redshifted 21-cm background is expected to be a powerful probe of the early Universe, carrying both cosmological and astrophysical information from a wide range of redshifts. In particular, the power spectrum of fluctuations in the 21-cm brightness temperature is anisotropic due to the line-of-sight velocity gradient, which in principle allows for a simple extraction of this information in the limit of linear fluctuations. However, recent numerical studies suggest that the 21-cm signal is actually rather complex, and its analysis likely depends on detailed model fitting. We present the first realistic simulation of the anisotropic 21-cm power spectrum over a wide period of early cosmic history. We show that on observable scales, the anisotropy is large and thus measurable at most redshifts, and its form tracks the evolution of 21-cm fluctuations as they are produced early on by Lyman-a radiation from stars, then switch to X-ray radiation from early heating sources, and finally to ionizing radiation from stars. In particular, we predict a redshift window during cosmic heating (at z ~ 15), when the anisotropy is small, during which the shape of the 21-cm power spectrum on large scales is determined directly by the average radial distribution of the flux from X-ray sources. This makes possible a model-independent reconstruction of the X-ray spectrum of the earliest sources of cosmic heating.

The first generation of stars produces a background of Lyman-Werner (LW) radiation which can photo-dissociate molecular hydrogen, increasing the mass of dark matter halos required to host star formation. Previous studies have determined the critical mass required for efficient molecular cooling with a constant LW background. However, the true background is expected to increase rapidly at early times. Neglecting this evolution could underestimate star formation in small halos that may have started to cool in the past when the LW intensity was much lower. Background evolution is a large source of uncertainty in pre-reionization predictions of the cosmological 21cm signal, which can be observed with future radio telescopes. To address this, we perform zero-dimentional one-zone calculations that follow the density, chemical abundances, and temperature of gas in the central regions of dark matter halos, including hierarchical growth and an evolving LW background. We begin by studying the physics of halos subjected to a background that increases exponentially with redshift. We find that when the intensity increases more slowly than $J_{\rm LW}(z) \propto 10^{-z/5}$, cooling in the past is a relatively small effect. We then self-consistently compute the cosmological LW background over $z=15-50$ and find that cooling in the past due to an evolving background has a modest impact. Finally, we compare these results to three-dimensional hydrodynamical cosmological simulations with varying LW histories. While only a small number of halos were simulated, the results are consistent with our one-zone calculations.

We explore the complete history of the 21-cm signal in the redshift range z = 7-40. This redshift range includes various epochs of cosmic evolution related to primordial star formation, and should be accessible to existing or planned low-frequency radio telescopes. We use semi-numerical computational methods to explore the fluctuation signal over wavenumbers between 0.03 and 1 Mpc$^{-1}$, accounting for the inhomogeneous backgrounds of Ly-$\alpha$, X-ray, Lyman-Werner and ionizing radiation. We focus on the recently noted expectation of heating dominated by a hard X-ray spectrum from high-mass X-ray binaries. We study the resulting delayed cosmic heating and suppression of gas temperature fluctuations, allowing for large variations in the minimum halo mass that contributes to star formation. We show that the wavenumbers at which the heating peak is detected in observations should tell us about the characteristic mean free path and spectrum of the emitted photons, thus giving key clues as to the character of the sources that heated the primordial Universe. We also consider the line-of-sight anisotropy, which allows additional information to be extracted from the 21-cm signal. For example, the heating transition at which the cosmic gas is heated to the temperature of the cosmic microwave background should be clearly marked by an especially isotropic power spectrum. More generally, an additional cross-power component $P_X$ directly probes which sources dominate 21-cm fluctuations. In particular, during cosmic reionization (and after the just-mentioned heating transition), $P_X$ is negative on scales dominated by ionization fluctuations and positive on those dominated by temperature fluctuations.

Understanding the formation and evolution of the first stars and galaxies represents one of the most exciting frontiers in astronomy. Since the universe was filled with neutral hydrogen at early times, the most promising method for observing the epoch of the first stars is using the prominent 21-cm spectral line of the hydrogen atom. Current observational efforts are focused on the reionization era (cosmic age t~500 Myr), with earlier times considered much more challenging. However, the next frontier of even earlier galaxy formation (t~200 Myr) is emerging as a promising observational target. This is made possible by a recently noticed effect of a significant relative velocity between the baryons and dark matter at early times. The velocity difference significantly suppresses star formation. The spatial variation of this suppression enhances large-scale clustering and produces a prominent cosmic web on 100 comoving Mpc scales in the 21-cm intensity distribution. This structure makes it much more feasible for radio astronomers to detect these early stars, and should drive a new focus on this era, which is rich with little-explored astrophysics.

Models and simulations of the epoch of reionization predict that spectra of the 21-cm transition of atomic hydrogen will show a clear fluctuation peak, at a redshift and scale, respectively, that mark the central stage of reionization and the characteristic size of ionized bubbles. This is based on the assumption that the cosmic gas was heated by stellar remnants - particularly X-ray binaries - to temperatures well above the cosmic microwave background at that time (~ 30 K). Here we show instead that the hard spectra (that is, spectra with more high-energy photons than low-energy photons) of X-ray binaries make such heating ineffective, resulting in a delayed and spatially uniform heating that modifies the 21-cm signature of reionization. Rather than looking for a simple rise and fall of the large-scale fluctuations (peaking at several millikelvin), we must expect a more complex signal also featuring a distinct minimum (at less than a millikelvin) that marks the rise of the cosmic mean gas temperature above the microwave background. Observing this signal, possibly with radio telescopes in operation today, will demonstrate the presence of a cosmic background of hard X-rays at that early time.

We present the first complete calculation of the history of the inhomogeneous 21-cm signal from neutral hydrogen during the era of the first stars. We use hybrid computational methods to capture the large-scale distribution of the first stars, whose radiation couples to the neutral hydrogen emission, and to evaluate the 21-cm signal from z ~ 15-35. In our realistic picture large-scale fluctuations in the 21-cm signal are sourced by the inhomogeneous density field and by the Ly-alpha and X-ray radiative backgrounds. The star formation is suppressed by two spatially varying effects: negative feedback provided by the Lyman-Werner radiative background, and supersonic relative velocities between the gas and dark matter. Our conclusions are quite promising: we find that the fluctuations imprinted by the inhomogeneous Ly-alpha background in the 21-cm signal at z ~ 25 should be detectable with the Square Kilometer Array.

Understanding the formation and evolution of the first stars and galaxies represents one of the most exciting frontiers in astronomy. Since the universe was filled with neutral hydrogen at early times, the most promising method for observing the epoch of the first stars is using the prominent 21-cm spectral line of the hydrogen atom. Current observational efforts are focused on the reionization era (cosmic age t ~ 500 Myr), with earlier times considered much more challenging. However, the next frontier of even earlier galaxy formation (t ~ 200 Myr) is emerging as a promising observational target. This is made possible by a recently noticed effect of a significant relative velocity between the baryons and dark matter at early times. The velocity difference suppresses star formation, causing a unique form of early luminosity bias. The spatial variation of this suppression enhances large-scale clustering and produces a prominent cosmic web on 100 comoving Mpc scales in the 21-cm intensity distribution. This structure makes it much more feasible for radio astronomers to detect these early stars, and should drive a new focus on this era, which is rich with little-explored astrophysics.

The formation of the first stars is an exciting frontier area in astronomy. Early redshifts z ~ 20 have become observationally promising as a result of a recently recognized effect of a supersonic relative velocity between the dark matter and gas. This effect produces prominent structure on 100 comoving Mpc scales, which makes it much more feasible to detect 21-cm fluctuations from the epoch of first heating. We use semi-numerical hybrid methods to follow for the first time the joint evolution of the X-ray and Lyman-Werner radiative backgrounds, including the effect of the supersonic streaming velocity on the cosmic distribution of stars. We incorporate self-consistently the negative feedback on star formation induced by the Lyman-Werner radiation, which dissociates molecular hydrogen and thus suppresses gas cooling. We find that the feedback delays the X-ray heating transition by a Delta z ~ 2, but leaves a promisingly large fluctuation signal over a broad redshift range. The large-scale power spectrum is predicted to reach a maximal signal-to-noise ratio of S/N ~ 3-4 at z ~ 18 (for a projected first-generation instrument), with S/N > 1 out to z ~ 22-23. We hope to stimulate additional numerical simulations as well as observational efforts focused on the epoch prior to cosmic reionization.

During cosmic reionization, the 21-cm brightness fluctuations were highly non-Gaussian, and complementary statistics can be extracted from the distribution of pixel brightness temperatures that are not derivable from the 21-cm power spectrum. One such statistic is the 21-cm difference PDF, the probability distribution function of the difference in the 21-cm brightness temperatures between two points, as a function of the distance between the points. Guided by 21-cm difference PDFs extracted from simulations, we perform a maximum likelihood analysis on mock observational data, and analyze the ability of present and future low-frequency radio array experiments to estimate the shape of the 21-cm difference PDF, and measure the history of cosmic reionization. We find that one-year data with an experiment such as the Murchison Wide-field Array should suffice for probing large scales during the mid-to-late stages of reionization, while a second-generation experiment should yield detailed measurements over a wide range of scales during most of the reionization era.

We report statistical results for dark matter (DM) velocity anisotropy, \beta, from a sample of some 6000 cluster-size halos (at redshift zero) identified in a \Lambda CDM hydrodynamical adaptive mesh refinement simulation performed with the Enzo code. These include profiles of \beta\ in clusters with different masses, relaxation states, and at several redshifts, modeled both as spherical and triaxial DM configurations. Specifically, although we find a large scatter in the DM velocity anisotropy profiles of different halos (across elliptical shells extending to at least ~$1.5 r_{vir}$), universal patterns are found when these are averaged over halo mass, redshift, and relaxation stage. These are characterized by a very small velocity anisotropy at the halo center, increasing outward to about 0.27 and leveling off at about $0.2 r_{vir}$. Indirect measurements of the DM velocity anisotropy fall on the upper end of the theoretically expected range. Though measured indirectly, the estimations are derived by using two different surrogate measurements - X-ray and galaxy dynamics. Current estimates of the DM velocity anisotropy are based on very small cluster sample. Increasing this sample will allow testing theoretical predictions, including the speculation that the decay of DM particles results in a large velocity boost. We also find, in accord with previous works, that halos are triaxial and likely to be more prolate when unrelaxed, whereas relaxed halos are more likely to be oblate. Our analysis does not indicate that there is significant correlation (found in some previous studies) between the radial density slope, \gamma, and \beta\ at large radii, $0.3 r_{vir} < r < r_{vir}$.

Dark and baryonic matter moved at different velocities in the early Universe, which strongly suppressed star formation in some regions. This was estimated to imprint a large-scale fluctuation signal of about 2 mK in the 21-cm spectral line of atomic hydrogen associated with stars at a redshift of 20, although this estimate ignored the critical contribution of gas heating due to X-rays and major enhancements of the suppression. A large velocity difference reduces the abundance of halos and requires the first stars to form in halos of about a million solar masses, substantially greater than previously expected. Here we report a simulation of the distribution of the first stars at z=20 (cosmic age of ~180 Myr), incorporating all these ingredients within a 400 Mpc box. We find that the 21-cm signature of these stars is an enhanced (10 mK) fluctuation signal on the 100-Mpc scale, characterized by a flat power spectrum with prominent baryon acoustic oscillations. The required sensitivity to see this signal is achievable with an integration time of a thousand hours with an instrument like the Murchison Wide-field Array or the Low Frequency Array but designed to operate in the range of 50-100 MHz.

We explore the ability of observations of the global brightness temperature of the 21-cm signal to constrain the reionization history and the properties of the ionizing sources. In order to describe the reionization signal, we employ either a commonly-used toy model or a structure formation model that parameterizes the properties of the ionizing sources. If the structure formation model captures the actual evolution of the reionization signal, then detecting the signal is somewhat easier than it would be for the toy model; using the toy model in this case also leads to systematic errors in reconstructing the reionization history, though a sufficiently sensitive experiment should be able to distinguish between the two models. We show that under optimistic assumptions regarding systematic noise and foreground removal, one-year observations of the global 21-cm spectrum should be able to detect a wide range of realistic models and measure the main features of the reionization history while constraining the key properties of the ionizing sources.

We show that various milestones of high-redshift galaxy formation, such as the formation of the first stars or the complete reionization of the intergalactic medium, occurred at different times in different regions of the universe. The predicted spread in redshift, caused by large-scale fluctuations in the number density of galaxies, is at least an order of magnitude larger than previous expectations that argued for a sharp end to reionization. This cosmic scatter in the abundance of galaxies introduces new features that affect the nature of reionization and the expectations for future probes of reionization, and may help explain the present properties of dwarf galaxies in different environments. The predictions can be tested by future numerical simulations and may be verified by upcoming observations. Current simulations, limited to relatively small volumes and periodic boundary conditions, largely omit cosmic scatter and its consequences. In particular, they artificially produce a sudden end to reionization, and they underestimate the number of galaxies by up to an order of magnitude at redshift 20.

Recently the initial supersonic relative velocity between the dark matter and baryons was shown to have an important effect on galaxy formation at high redshift. We study the impact of this relative motion on the distribution of the star-forming halos and on the formation redshift of the very first star. We include a new aspect of the relative velocity effect found in recent simulations by fitting their results to obtain the spatially-varying minimum halo mass needed for molecular cooling. Thus, the relative velocities have three separate effects: suppression of the halo abundance, suppression of the gas content within each halo, and boosting of the minimum cooling mass. We show that the two suppressions (of gas content and of halo abundance) are the primary effects on the small minihalos that cannot form stars, while the cooling mass boost combines with the abundance suppression to produce order unity fluctuations in stellar density. We quantify the large-scale inhomogeneity of galaxies, finding that 68% of the star formation (averaged on a 3 Mpc scale) is confined to 35% of the volume at z=20 (and just 18% at z=40). In addition, we estimate the redshift of the first star to be z ~ 65, which includes a delay of Dz ~ 5 due to the relative velocity.

Clusters of galaxies are uniquely important cosmological probes of the evolution of the large scale structure, whose diagnostic power depends quite significantly on the ability to reliably determine their masses. Clusters are typically modeled as spherical systems whose intracluster gas is in strict hydrostatic equilibrium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure), with the gravitational field dominated by dark matter, assumptions that are only rough approximations. In fact, numerical simulations indicate that galaxy clusters are typically triaxial, rather than spherical, and that turbulent gas motions (induced during hierarchical merger events) provide an appreciable pressure component. Extending our previous work, we present results of a joint analysis of X-ray, weak and strong lensing measurements of Abell 1689. The quality of the data allows us to determine both the triaxial shape of the cluster and the level of non-thermal pressure that is required if the intracluster gas is in hydrostatic equilibrium. We find that the dark matter axis ratios are 1.24 +/- 0.13 and 2.02 +/- 0.01 on the plane of the sky and along the line of sight, respectively, and that about 20% of the pressure is non-thermal. Our treatment demonstrates that the dynamical properties of clusters can be determined in a (mostly) bias-free way, enhancing the use of clusters as more precise cosmological probes.

We run very large cosmological N-body hydrodynamical simulations in order to study statistically the baryon fractions in early dark matter halos. We critically examine how differences in the initial conditions affect the gas fraction in the redshift range z = 11--21. We test three different linear power spectra for the initial conditions: (1) A complete heating model, which is our fiducial model; this model follows the evolution of overdensities correctly, according to Naoz & Barkana (2005), in particular including the spatial variation of the speed of sound of the gas due to Compton heating from the CMB. (2) An equal-{\delta} model, which assumes that the initial baryon fluctuations are equal to those of the dark matter, while conserving sigma8 of the total matter. (3) A mean cs model, which assumes a uniform speed of sound of the gas. The latter two models are often used in the literature. We calculate the baryon fractions for a large sample of halos in our simulations. Our fiducial model implies that before reionization and significant stellar heating took place, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was ~ 3 * 10^4 Msun. However, the alternative models yield a wrong (higher by about 50%) minimum mass, since the system retains a memory of the initial conditions. We also demonstrate this using the "filtering mass" from linear theory, which accurately describes the evolution of the baryon fraction throughout the simulated redshift range.

We study the effect of the relative velocity of dark matter and baryonic fluids after the epoch of recombination on the evolution of the first bound objects in the early universe. Recent work has shown that, although relative motion of the two fluids is formally a second order effect in density, it has a dramatic impact on the formation and distribution of the first cosmic structures. Focusing on the gas content, we analyze the effect of relative velocity on the properties of halos over a wide range of halo masses and redshifts. We calculate accurately the linear evolution of the baryon and dark matter fluctuations, and quantify the resulting effect on halos based on an analytical formalism that has been carefully checked with simulations in the case with no relative velocity. We estimate the effect on the abundance of and gas fraction in early halos. We find that the relative velocity effect causes several changes: (i) the characteristic mass that divides gas-rich and gas-poor halos is increased by roughly an order of magnitude, from 2 10^4 Msun to about 2 10^5 Msun; (ii) this characteristic mass has a large scatter (full width at half maximum is ~ 1.5 10^5 Msun at z=20); (iii) the fraction of baryons in star-less gas minihalos is suppressed by a factor of 4 at z=20; (iv) the fraction of baryons in halos that can cool and form stars is suppressed by a factor of 1.5 at z=20; and (v) there are enhanced spatial variations of these various fractions.

Baryonic acoustic oscillations (BAOs) modulate the density ratio of baryons to dark matter across large regions of the Universe. We show that the associated variation in the mass-to-light ratio of galaxies should generate an oscillatory, scale-dependent bias of galaxies relative to the underlying distribution of dark matter. A measurement of this effect would calibrate the dependence of the characteristic mass-to-light ratio of galaxies on the baryon mass fraction in their large scale environment. This bias, though, is unlikely to significantly affect measurements of BAO peak positions.

We present the results of a strong-lensing analysis of a complete sample of 12 very luminous X-ray clusters at $z>0.5$ using HST/ACS images. Our modelling technique has uncovered some of the largest known critical curves outlined by many accurately-predicted sets of multiple images. The distribution of Einstein radii has a median value of $\simeq28\arcsec$ (for a source redshift of $z_{s}\sim2$), twice as large as other lower-$z$ samples, and extends to $55\arcsec$ for MACS J0717.5+3745, with an impressive enclosed Einstein mass of $7.4\times10^{14} M_{\odot}$. We find that 9 clusters cover a very large area ($>2.5 \sq \arcmin$) of high magnification ($\mu > \times10$) for a source redshift of $z_{s}\sim8$, providing primary targets for accessing the first stars and galaxies. We compare our results with theoretical predictions of the standard $\Lambda$CDM model which we show systematically fall short of our measured Einstein radii by a factor of $\simeq1.4$, after accounting for the effect of lensing projection. Nevertheless, a revised analysis once arc redshifts become available, and similar analyses of larger samples, are needed in order to establish more precisely the level of discrepancy with $\Lambda$CDM predictions.

We use extensive measurements of the cluster A1689 to assess the expected similarity in the dynamics of galaxies and dark matter (DM) in their motion as collisionless `particles' in the cluster gravitational potential. To do so we derive the radial profile of the specific kinetic energy of the cluster galaxies from the Jeans equation and observational data. Assuming that the specific kinetic energies of galaxies and DM are roughly equal, we obtain the mean value of the DM velocity anisotropy parameter, and the DM density profile. Since this deduced profile has a scale radius that is higher than inferred from lensing observations, we tested the validity of the assumption by repeating the analysis using results of simulations for the profile of the DM velocity anisotropy. Results of both analyses indicate a significant difference between the kinematics of galaxies and DM within $r \lesssim 0.3r_{\rm vir}$. This finding is reflected also in the shape of the galaxy number density profile, which flattens markedly with respect to the steadily rising DM profile at small radii. Thus, $r \sim 0.3r_{\rm vir}$ seems to be a transition region interior to which collisional effects significantly modify the dynamical properties of the galaxy population with respect to those of DM in A1689

In the coming decade, low-frequency radio arrays will begin to probe the epoch of reionization via the redshifted 21-cm hydrogen line. Successful interpretation of these observations will require effective statistical techniques for analyzing the data. Due to the difficulty of these measurements, it is important to develop techniques beyond the standard power spectrum analysis in order to offer independent confirmation of the reionization history, probe different aspects of the topology of reionization, and have different systematic errors. In order to assess the promise of probability distribution functions (PDFs) as statistical analysis tools in 21-cm cosmology, we first measure the 21-cm brightness temperature (one-point) PDFs in six different reionization simulations. We then parametrize their most distinct features by fitting them to a simple model. Using the same simulations, we also present the first measurements of difference PDFs in simulations of reionization. We find that while these statistics probe the properties of the ionizing sources, they are relatively independent of small-scale, sub-grid astrophysics. We discuss the additional information that the difference PDF can provide on top of the power spectrum and the one-point PDF.

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.

We examine the dynamics structure of the rich cluster A1689, combining VLT/VIMOS spectroscopy with Subaru/Suprime-Cam imaging. The radial velocity distribution of $\sim 500$ cluster members is bounded by a pair of clearly defined velocity caustics, with a maximum amplitude of $\sim|4000|$ km/s at $\simeq$ 300 h$^{-1}$ kpc, beyond which the amplitude steadily declines, approaching zero velocity at a limiting radius of $\sim$ 2 h$^{-1}$ Mpc. We derive the 3D velocity anisotropy and galaxy number density profiles using a model-independent method to solve the Jeans equation, simultaneously incorporating the observed velocity dispersion profile, the galaxy counts from deep Subaru imaging, and our previously derived cluster mass profile from a joint lensing and X-ray analysis. The velocity anisotropy is found to be predominantly radial at large radius, becoming increasingly tangential towards the center, in accord with expectations. We also analyze the galaxy data independently of our previous analysis using two different methods: The first is based on a solution of the Jeans equation assuming an NFW form for the mass distribution, whereas in the second method the caustic amplitude is used to determine the escape velocity. The cluster virial mass derived by both of these dynamical methods is in good agreement with results from our earlier lensing and X-ray analysis. We also confirm the high NFW concentration parameter, with results from both methods combined to yield $c_{\rm vir}>13$ (1$\sigma$). The inferred virial radius is consistent with the limiting radius where the caustics approach zero velocity and where the counts of cluster members drop off, suggesting that infall onto A1689 is currently not significant.

Gravitational lensing observations of massive X-ray clusters imply a steep characteristic density profile marked by a central concentration of dark matter. The observed mass fraction within a projected radius of 150 kpc is twice that found in state-of-the-art dark matter simulations of the standard Lambda-CDM cosmology. A central baryon enhancement that could explain this discrepancy is not observed, leaving a major puzzle. We propose a solution based on the merger histories of clusters. A significant fraction of the final dark matter content of a cluster halo originates within galaxy-sized halos, in which gas can cool and compress the dark matter core to high densities. The subsequent tidal stripping of this compressed dark matter occurs in denser regions that are closer to the center of the cluster halo. Eventually, the originally cooled gas must be dispersed into the intracluster medium through feedback, for consistency with observations that do not find central baryon enhancements in clusters. Still, the early adiabatic compression of the galactic dark matter leaves a net effect on the cluster. Using a simple model for this process, we show that the central cluster profile is substantially modified, potentially explaining the observed discrepancy.

We study the gas content of halos in the early universe using high resolution hydrodynamical simulations. We extract from the simulations and also predict based on linear theory the halo mass for which the enclosed baryon fraction equals half of the mean cosmic fraction. We find a rough agreement between the simulations and the predictions, which suggests that during the high-redshift era before stellar heating, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was $\sim 3 \times 10^4$ M$_\odot$. We also carry out a detailed resolution analysis and show that in order to determine a halo's gas fraction even to 20% accuracy the halo must be resolved into at least 500 dark matter particles.

We consider the early stages of cosmic hydrogen or helium reionization, when ionizing sources were still rare. We show that Poisson fluctuations in the galaxy distribution substantially affected the early bubble size distribution, although galaxy clustering was also an essential factor even at the earliest times. We find that even at high redshifts, a significant fraction of the ionized volume resided in bubbles containing multiple sources, regardless of the ionizing efficiency of sources or of the reionization redshift. In particular, for helium reionization by quasars, one-source bubbles last dominated (i.e., contained 90% of the ionized volume) at some redshift above z=7.3, and hydrogen reionization by stars achieved this milestone at z>23. For the early generations of atomic-cooling halos or molecular-hydrogen-cooling halos, one-source ionized regions dominated the ionized volume only at z>31 and z>48, respectively. To arrive at these results we develop a statistical model for the effect of density correlations and discrete sources on reionization and solve it with a Monte Carlo method.

The Einstein radius of a cluster provides a relatively model-independent measure of the mass density of a cluster within a projected radius of ~ 150 kpc, large enough to be relatively unaffected by gas physics. We show that the observed Einstein radii of four well-studied massive clusters, for which reliable virial masses are measured, lie well beyond the predicted distribution of Einstein radii in the standard LambdaCDM model. Based on large samples of numerically simulated cluster-sized objects with virial masses ~ 1e15 solar, the predicted Einstein radii are only 15-25'', a factor of two below the observed Einstein radii of these four clusters. This is because the predicted mass profile is too shallow to exceed the critical surface density for lensing at a sizable projected radius. After carefully accounting for measurement errors as well as the biases inherent in the selection of clusters and the projection of mass measured by lensing, we find that the theoretical predictions are excluded at a 4-sigma significance. Since most of the free parameters of the LambdaCDM model now rest on firm empirical ground, this discrepancy may point to an additional mechanism that promotes the collapse of clusters at an earlier time thereby enhancing their central mass density.

We explore the ability of measurements of the 21-cm power spectrum during reionization to enable the simultaneous reconstruction of the reionization history and the properties of the ionizing sources. For various sets of simulated 21-cm measurements, we perform maximum likelihood fits in order to constrain the reionization and galaxy formation histories. We employ a flexible six-parameter model that parametrizes the uncertainties in the properties of high-redshift galaxies. The computational speed needed is attained through the use of an analytical model that is in reasonable agreement with numerical simulations of reionization. We find that one-year observations with the MWA array should measure the cosmic ionized fraction to ~ 1% accuracy at the very end of reionization, and a few percent accuracy around the mid-point of reionization. The mean halo mass of the ionizing sources should be measureble to 10% accuracy when reionization is 2/3 of the way through, and to 20% accuracy throughout the central stage of reionization, if this mass is anywhere in the range 1/3 to 100 billion solar masses.

We carry out a comprehensive joint analysis of high quality HST/ACS and Chandra measurements of A1689, from which we derive mass, temperature, X-ray emission and abundance profiles. The X-ray emission is smooth and symmetric, and the lensing mass is centrally concentrated indicating a relaxed cluster. Assuming hydrostatic equilibrium we deduce a 3D mass profile that agrees simultaneously with both the lensing and X-ray measurements. However, the projected temperature profile predicted with this 3D mass profile exceeds the observed temperature by ~30% at all radii, a level of discrepancy comparable to the level found for other relaxed clusters. This result may support recent suggestions from hydrodynamical simulations that denser, more X-ray luminous small-scale structure can bias observed temperature measurements downward at about the same (~30%) level. We determine the gas entropy at 0.1r_{vir} (where r_{vir} is the virial radius) to be ~800 keV cm^2, as expected for a high temperature cluster, but its profile at >0.1r_{vir} has a power-law form with index ~0.8, considerably shallower than the ~1.1 index advocated by theoretical studies and simulations. Moreover, if a constant entropy ''floor'' exists at all, then it is within a small region in the inner core, r<0.02r_{vir}, in accord with previous theoretical studies of massive clusters.

New observations over the next few years of the emission of distant objects will help unfold the chapter in cosmic history around the era of the first galaxies. These observations will use the neutral hydrogen emission or absorption at a wavelength of 21-cm as a detector of the hydrogen abundance. We predict the signature on the 21-cm signal of the early generations of galaxies. We calculate the 21-cm power spectrum including two physical effects that were neglected in previous calculations. The first is the redistribution of the UV photons from the first galaxies due to their scattering off of the neutral hydrogen, which results in an enhancement of the 21-cm signal. The second is the presence of an ionized hydrogen bubble near each source, which produces a cutoff at observable scales. We show that the resulting clear signature in the 21-cm power spectrum can be used to detect and study the population of galaxies that formed just 200 million years after the Big Bang.

A new generation of radio telescopes are currently being built with the goal of tracing the cosmic distribution of atomic hydrogen at redshifts 6-15 through its 21-cm line. The observations will probe the large-scale brightness fluctuations sourced by ionization fluctuations during cosmic reionization. Since detailed maps will be difficult to extract due to noise and foreground emission, efforts have focused on a statistical detection of the 21-cm fluctuations. During cosmic reionization, these fluctuations are highly non-Gaussian and thus more information can be extracted than just the one-dimensional function that is usually considered, i.e., the correlation function. We calculate a two-dimensional function that if measured observationally would allow a more thorough investigation of the properties of the underlying ionizing sources. This function is the probability distribution function (PDF) of the difference in the 21-cm brightness temperature between two points, as a function of the separation between the points. While the standard correlation function is determined by a complicated mixture of contributions from density and ionization fluctuations, we show that the difference PDF holds the key to separately measuring the statistical properties of the ionized regions.

Analytical approaches to galaxy formation and reionization are based on the mathematical problem of random walks with barriers. The statistics of a single random walk can be used to calculate one-point distributions ranging from the mass function of virialized halos to the distribution of ionized bubble sizes during reionization. However, an analytical calculation of two-point correlation functions or of spatially-dependent feedback processes requires the joint statistics of random walks at two different points. An accurate analytical expression for the statistics of two correlated random walks has been previously found only for the case of a constant barrier height. However, calculating bubble sizes or accurate statistics for halo formation involves more general barriers that can often be approximated as linear barriers. We generalize the two-point solution with constant barriers to linear barriers, and apply it as an illustration to calculate the correlation function of cosmological 21-cm fluctuations during reionization.

The intergalactic medium - the cosmic gas that fills the great spaces between the galaxies - is affected by processes ranging from quantum fluctuations in the very early universe to radiative emission from newly-formed stars. This gives the intergalactic medium a dual role as a powerful probe both of fundamental physics and of astrophysics. The heading of fundamental physics includes conditions in the very early universe and basic cosmological parameters. The astrophysics refers to chapters of the long cosmic history of stars and galaxies that are being revealed through the effects of stellar feedback on the cosmic gas. This review describes the physics of the intergalactic medium, focusing on recent theoretical and observational developments in understanding early cosmic history. In particular, radiation from the earliest generation of stars is thought to have transformed the universe, turning the surrounding atoms into free electrons and ions. Knowing exactly when and how the reionization process happened is a primary goal of cosmologists, because this would tell us when the early stars and black holes formed and in what kinds of galaxies. Cosmic reionization is beginning to be understood with the help of theoretical models and computer simulations. A wide variety of instruments currently under design will probe the first sources of light during an epoch in cosmic history that has been largely unexplored so far. The new observations and the challenges for theoretical models and numerical simulations will motivate intense work in this field over the coming decade. (abridged)

Large telescopes have allowed astronomers to observe galaxies that formed as early as 850 million years after the Big Bang. We predict when the first star that astronomers can observe formed in the universe, accounting for the first time for the size of the universe and for three essential ingredients: the light travel time from distant galaxies, Poisson and density fluctuations on all scales, and the effect of very early cosmic history on galaxy formation. We find that the first observable star is most likely to have formed 30 million years after the Big Bang (at redshift 65), much earlier than previously expected. Also, the first galaxy as massive as our own Milky Way likely formed when the universe was only 400 Myr old (at redshift 11). We also show that significant modifications are required in current methods of numerically simulating the formation of galaxies at redshift 20 and above.

The delay in light travel time along the line of sight generates an anisotropy in the power spectrum of 21cm brightness fluctuations from the epoch of reionization. We show that when the fluctuations in the neutral hydrogen fraction become non-linear at the later stages of reionization, the light-cone anisotropy becomes of order unity on scales >50 comoving Mpc. During this period the density fluctuations and the associated anisotropy generated by peculiar velocities are negligible in comparison.

The reionization of cosmic hydrogen, left over from the big bang, increased its temperature to >~ 1.e4 K. This photo-heating resulted in an increase of the minimum mass of galaxies and hence a suppression of the cosmic star formation rate. The affected population of dwarf galaxies included the progenitors of massive galaxies that formed later. We show that a massive galaxy at a redshift z >~ 6 should show a double-peaked star formation history marked by a clear break. This break reflects the suppression signature from reionization of the region in which the galaxy was assembled. Since massive galaxies originate in overdense regions where cosmic evolution is accelerated, their environment reionizes earlier than the rest of the universe. For a galaxy of ~ 1.e12 M_solar in stars at a redshift of z ~ 6.5, the star formation rate should typically be suppressed at a redshift z >~ 10 since the rest of the universe is known to have reionized by z >~ 6.5. Indeed, this is inferred to be the case for HUDF-JD2, a massive galaxy which is potentially at z ~ 6.5 but is inferred to have formed the bulk of its 3.e11 M_solar in stars at z >~ 9.

We study the relation between the star formation intensity of galaxies and the extinction by dust of their emitted light. We employ a detailed statistical analysis of Hubble Deep Field North data to show a clear positive correlation between the extinction and star formation intensity at all epochs from redshift 0.4 to 6.5. The extinction evidently increases with time for a given star formation intensity, consistent with the expected increase in the metallicity with time. Our observational results are well fitted at all epochs by a double power-law model with a fixed shape that simply shifts with redshift. The correlation between the extinction and the star formation intensity can be interpreted by combining two other trends: the correlation between the star formation rate and the gas content of galaxies, and the evolution of the dust-to-gas ratio in galaxies. If we assume that Kennicutt's observed relation for the former is valid at each redshift, then our findings imply an interesting variation in the dust-to-gas ratio in galaxies within each epoch and with time, and suggest new ways to investigate the cosmic evolution of this quantity.

We reconsider the Alcock-Paczynski effect on 21cm fluctuations from high redshift, focusing on the 21cm power spectrum. We show that at each accessible redshift both the angular diameter distance and the Hubble constant can be determined from the power spectrum. Furthermore, this is possible using anisotropies that depend only on linear density perturbations and not on astrophysical sources of 21cm fluctuations. We show that measuring these quantities at high redshift would not just confirm results from the cosmic microwave background but provide appreciable additional sensitivity to cosmological parameters and dark energy.

We investigate the suppression of the baryon density fluctuations compared to the dark matter in the linear regime. Previous calculations predict that the suppression occurs up to a characteristic mass scale of ~ 1,000,000 solar masses, which suggests that pressure has a central role in determining the properties of the first luminous objects at early times. We show that the expected characteristic mass scale is in fact substantially lower (by a factor of ~ 3-10, depending on redshift), and thus the effect of baryonic pressure on the formation of galaxies up to reionization is only moderate. This result is due to the influence on perturbation growth of the high pressure that prevailed in the period from cosmic recombination to z ~ 200, when the gas began to cool adiabatically and the pressure then dropped. At z ~ 10 the suppression of the baryon fluctuations is still sensitive to the history of pressure in this high-redshift era. We calculate the fraction of the cosmic gas that is in minihalos and find that it is substantially higher than would be expected with the previously-estimated characteristic mass. Expanding our investigation to the non-linear regime, we calculate in detail the spherical collapse of high-redshift objects in a Lambda-CDM universe. We include the gravitational contributions of the baryons and radiation and the memory of their kinematic coupling before recombination. We use our results to predict a more accurate halo mass function as a function of redshift.

The earliest generation of stars, far from being a mere novelty, transformed the universe from darkness to light. The first atoms to form after the Big Bang filled the universe with atomic hydrogen and a few light elements. As gravity pulled gas clouds together, the first stars ignited and their radiation turned the surrounding atoms into ions. By looking at gas between us and distant galaxies, we know that this ionization eventually pervaded all space, so that few hydrogen atoms remain today between galaxies. Knowing exactly when and how it did so is a primary goal of cosmologists, because this would tell us when the early stars formed and in what kinds of galaxies. Although this ionization is beginning to be understood by using theoretical models and computer simulations, a new generation of telescopes is being built that will map atomic hydrogen throughout the universe.

The first galaxies that formed at a redshift ~20-30 emitted continuum photons with energies between the Lyman-alpha and Lyman limit wavelengths of hydrogen, to which the neutral universe was transparent except at the Lyman-series resonances. As these photons redshifted or scattered into the Lyman-alpha resonance they coupled the spin temperature of the 21cm transition of hydrogen to the gas temperature, allowing it to deviate from the microwave background temperature. We show that the fluctuations in the radiation emitted by the first galaxies produced strong fluctuations in the 21cm flux before the Lyman-alpha coupling became saturated. The fluctuations were caused by biased inhomogeneities in the density of galaxies, along with Poisson fluctuations in the number of galaxies. Observing the power-spectra of these two sources would probe the number density of the earliest galaxies and the typical mass of their host dark matter halos. The enhanced amplitude of the 21cm fluctuations from the era of Lyman-alpha coupling improves considerably the practical prospects for their detection.

We calculate the evolution of linear density and temperature perturbations in a universe with dark matter, baryons, and radiation, from cosmic recombination until the epoch of the first galaxies. In addition to gravity, the perturbations are effected by electron scattering with the radiation, by radiation pressure, and by gas pressure. We include the effect of spatial fluctuations in the baryonic sound speed and show that they induce a >10% change in the baryonic density power spectrum on small scales, and a larger change on all scales in the power spectrum of gas temperature fluctuations. A precise calculation of the growth of linear perturbations is essential since they provide the initial conditions for the formation of galaxies and they can also be probed directly via cosmological 21cm fluctuations. We also show that in general the thermal history of the cosmic gas can be measured from 21cm fluctuations using a small-scale anisotropic cutoff due to the thermal width of the 21cm line.

After cosmological recombination, the primordial hydrogen gas decoupled from the cosmic microwave background (CMB) and fell into the gravitational potential wells of the dark matter. The neutral hydrogen imprinted acoustic oscillations on the pattern of brightness fluctuations due to its redshifted 21cm absorption of the CMB. Unlike CMB temperature fluctuations which probe the power spectrum at cosmic recombination, we show that observations of the 21cm fluctuations at z ~ 20-200 can measure four separate fluctuation modes (with a fifth mode requiring very high precision), thus providing a unique probe of the geometry and composition of the universe.

Fluctuations in the 21cm brightness from cosmic hydrogen at redshifts z > 6 were sourced by the primordial density perturbations from inflation as well as by the radiation from galaxies. We propose a method to separate these components based on the angular dependence of the 21cm fluctuation power spectrum. Peculiar velocities increase the power spectrum by a factor of ~ 2 compared to density fluctuations alone, and introduce an angular dependence in Fourier space. The resulting angular structure relative to the line of sight facilitates a simple separation of the power spectrum into several components, permitting an unambiguous determination of the primordial power spectrum of density fluctuations, and of the detailed properties of all astrophysical sources of 21cm fluctuations. We also demonstrate that there is no significant information to be gained from 21cm measurements on large angular scales. Thus, observational analyses can be confined to small angular areas where projection effects are negligible, and theoretical predictions should focus on the three-dimensional power spectrum of 21cm fluctuations.

Recently, Pello et al. (2004) claimed to have discovered a galaxy at a redshift z=10, and identified a feature in its spectrum with a hydrogen Lyman-alpha emission line. If this identification is correct, we show that the intergalactic medium (IGM) around the galaxy must be significantly ionized; otherwise, the damping wing of Lyman-alpha absorption by the neutral IGM would have greatly suppressed the emission line. We find either that the large-scale region surrounding this galaxy must have been largely reionized by z=10 (with a neutral fraction <0.4) or that the stars within the galaxy must be massive (>100 solar masses), and hence capable of producing a sufficiently large HII region around it. We generalize these conclusions and derive the maximum Lyman-alpha line flux for a given UV continuum flux of galaxies prior to the epoch of reionization.

Lyman-alpha absorption is a prominent cosmological tool for probing both galactic halos and the intergalactic medium at high redshift. We consider a variety of sources that can be used as the Lyman-alpha emitters for this purpose. Among these sources, we argue that quasars are the best probes of the evolution of massive halos, while gamma-ray bursts represent the cleanest sources for studying the reionization of the intergalactic medium.

Motivated by the recent direct detection of cosmological gas infall, we develop an analytical model for calculating the mean density profile around an initial overdensity that later forms a dark matter halo. We account for the problem of peaks within peaks; when considering a halo of a given mass we ensure that this halo is not a part of a larger virialized halo. For halos that represent high-sigma fluctuations we recover the usual result that such halos preferentially lie within highly overdense regions; in the limit of very low-sigma fluctuations, on the other hand, we show that halos tend to lie within voids. Combined with spherical collapse, our results yield the typical density and velocity profiles of the gas and dark matter that surround virialized halos.

Recent observations have shown that, only a billion years after the Big Bang, the Universe was already lit up by bright quasars fuelled by the infall of gas onto supermassive black holes. The masses of these early black holes are inferred from their luminosities to be > 10^9 solar masses, which is a difficult theoretical challenge to explain. Like nearby quasars, the early objects could have formed in the central cores of massive host galaxies. The formation of these hosts could be explained if, like local large galaxies, they were assembled gravitationally inside massive (> 10^12 solar mass) halos of dark matter. There has hitherto been no observational evidence for the presence of these massive hosts or their surrounding halos. Here we show that the cosmic gas surrounding each halo must respond to its strong gravitational pull, where absorption by the infalling hydrogen produces a distinct spectral signature. That signature can be seen in recent data.

Recently, Lanzetta et al. (2002) have measured the distribution of star formation rate intensity in galaxies at various redshifts. This data set has a number of advantages relative to galaxy luminosity functions; the effect of surface-brightness dimming on the selection function is simpler to understand, and this data set also probes the size distribution of galactic disks. We predict this function using semi-analytic models of hierarchical galaxy formation in a LCDM cosmology. We show that the basic trends found in the data follow naturally from the redshift evolution of dark matter halos. The data are consistent with a constant efficiency of turning gas into stars in galaxies, with a best-fit value of 2%, where dust obscuration is neglected; equivalently, the data are consistent with a cosmic star formation rate which is constant to within a factor of two at all redshifts above two. However, the practical ability to use this kind of distribution to measure the total cosmic star formation rate is limited by the predicted shape of an approximate power law with a smoothly varying power, without a sharp break.

We develop an analytical formalism that is suitable for studying inhomogeneous structure formation, by studying the joint statistics of dark matter halos forming at two points. Extending the Bond et al. (1991) derivation of the mass function of virialized halos, based on excursion sets, we derive an approximate analytical expression for the ``bivariate'' mass function of halos forming at two redshifts and separated by a fixed comoving Lagrangian distance. Our approach also leads to a self-consistent expression for the nonlinear biasing and correlation function of halos, generalizing a number of previous results including those by Kaiser (1984) and Mo & White (1996). We compare our approximate solutions to exact numerical results within the excursion-set framework and find them to be consistent to within 2% over a wide range of parameters. Our formalism can be used to study various feedback effects during galaxy formation analytically, as well as to simply construct observable quantities dependent on the spatial distribution of objects. A code that implements our method is publicly available at http://www.arcetri.astro.it/~evan/Gemini

We show that the gaseous halos of collapsed objects introduce a substantial cumulative opacity to ionizing radiation, even after the smoothly distributed hydrogen in the intergalactic medium has been fully reionized. This opacity causes a delay of around unity in redshift between the time of the overlap of ionized bubbles in the intergalactic medium and the lifting of complete Gunn-Peterson Lyman alpha absorption. The minihalos responsible for this screening effect are not resolved by existing numerical simulations of reionization.

In light of recent observations of spectra of the quasars SDSS 1030+0524 (Becker et al. 2001) and SDSS 1044-0125 (Djorgovski et al. 2001), we study the observational signatures of different stages of the reionization epoch. During the initial, pre-overlap stage, the hydrogen throughout the universe is neutral except for isolated HII bubbles due to individual ionizing sources. The central stage of reionization is the overlap stage, during which the individual HII bubbles overlap each other and reionize the low-density gas which takes up most of the volume of the universe. Some neutral hydrogen remains in dense clumps which are then slowly reionized during post-overlap. We show that both of the recent observations are consistent with the post-overlap stage. Becker et al. (2001) may have observed the universe in the era before the end of overlap, but a conclusive proof of this requires the observation of similar intervals of Gunn-Peterson absorption at the same redshift along several additional lines of sight.

We study the constraints that high-redshift structure formation in the universe places on warm dark matter (WDM) dominated cosmological models. We modify the extended Press-Schechter formalism to derive the halo mass function in WDM models. We show that our predictions agree with recent numerical simulations at low redshift over the halo masses of interest. Applying our model to galaxy formation at high redshift, we find that the loss of power on small scales, together with the delayed collapse of low-mass objects, results in strong limits on the root-mean-square velocity dispersion v_rms of the WDM particles at z=0. For fermions decoupling while relativistic, these limits are equivalent to constraints on the mass m_X of the particles. The presence of a 4 billion solar mass black hole at z=5.8, believed to power the quasar SDSS 1044-1215, implies m_X > 0.5 keV (or v_rms < 0.10 km/s), assuming that the quasar is unlensed and radiating at or below the Eddington limit. Reionization by z=5.8 also implies a limit on m_X. If high-redshift galaxies produce ionizing photons with an efficiency similar to their z=3 counterparts, we find m_X > 1.2 keV (or v_rms < 0.03 km/s). However, given the uncertainties in current measurements from the proximity effect of the ionizing background at redshift 3, values of m_X as low as 0.75 keV (v_rms = 0.06 km/s) are not ruled out. The limit weakens further if, instead, the ionizing-photon production efficiency is greater at high z, but this limit will tighten considerably if reionization is shown in the future to have occurred at higher redshifts. WDM models with m_X < 1 keV (v_rms > 0.04 km/s) produce a low-luminosity cutoff in the high-redshift galaxy luminosity function which is directly detectable with the Next Generation Space Telescope (abridged).

The formation of the first stars and quasars marks the transformation of the universe from its smooth initial state to its clumpy current state. In popular cosmological models, the first sources of light began to form at redshift 30 and reionized most of the hydrogen in the universe by redshift 7. Current observations are at the threshold of probing the hydrogen reionization epoch. The study of high-redshift sources is likely to attract major attention in observational and theoretical cosmology over the next decade.

Warm Dark Matter (WDM) models have recently been resurrected to resolve apparent conflicts of Cold Dark Matter (DM) models with observations. Endowing the DM particles with non-negligible velocities causes free-streaming, which suppresses the primordial power spectrum on small scales. The choice of a root-mean-square velocity dispersion v(rms) = 0.05 km/s at redshift z=0 (corresponding to a particle mass of 1 keV if the WDM particles are fermions decoupling while relativistic) helps alleviate most, but probably not all, of the small-scale problems faced by CDM. An important side-effect of the particle velocities is the severe decrease in the number of collapsed halos at high redshift. This is caused both by the loss of small-scale power, and by the delay in the collapse of the smallest individual halos (with masses near the effective Jeans mass of the DM). The presence of early halos is required in order (1) to host either early quasars or galaxies that can reionize the universe by redshift z=5.8, and (2) to allow the growth of the supermassive black hole believed to power the recently discovered quasar SDSS 1044-1215 at this redshift. We quantify these constraints using a modified Press-Schechter formalism, and find v(rms) < 0.04 km/s (or m_X > 1 keV). If future observations uncover massive black holes at z > 10, or reveal that reionization occurred at z > 10, this could conclusively rule out WDM models as the solution to the small-scale crisis of the CDM paradigm.

The first light from stars and quasars ended the ``dark ages'' of the universe and led to the reionization of hydrogen by redshift 7. Current observations are at the threshold of probing this epoch. The study of high-redshift sources is likely to attract major attention in observational and theoretical cosmology over the next decade.

Time-delay measurements of 0957+561 based on radio monitoring now agree with optical monitoring. Recent models incorporate many recent observations, but the systematic uncertainties in the models still dominate the uncertainty in the cosmological results. VLBI observations of the milli-arcsec jet structure have provided the most important set of modeling constraints. We present new observations of this structure, made at 18 cm with the Very Long Baseline Array, the Very Large Array, and the Green Bank 140 ft. telescope. Compared to the data set of Garrett et al. (1994), ours has a similar theoretical noise level (20 microJy), but nearly twice as many baselines (66 compared to 36) and more uniform coverage of the UV plane. Our initial maps confirm the basic jet structures seen by Garrett et al., with comparable RMS noise (60-80 microJy), and are in reasonable agreement with the jet structure model of Barkana et al. (1999). The maps also hint at a new compact component leaving the core, and a previously undetected diffuse component at the end of the jet. Further refinements in fringe-fitting, mapping, and self-calibration should yield significantly lower RMS noise (closer to the theoretical level). We will use a modified form of the CalTech VLBI Modelfit program (written by Barkana et al. 1999) to determine the magnification matrix between the images. The refined magnification matrix should provide stronger constraints on the lens model, and thus reduce the uncertainty in the cosmological results.

Cold dark matter (CDM) models predict small-scale structure in excess of observations of the cores and abundance of dwarf galaxies. These problems might be solved, and the virtues of CDM models retained, even without postulating {\it ad hoc} dark matter particle or field interactions, if the dark matter is composed of ultra-light scalar particles (m ~ 10^{-22}eV), initially in a (cold) Bose-Einstein condensate, similar to axion dark matter models. The wave properties of the dark matter stabilize gravitational collapse providing halo cores and sharply suppressing small-scale linear power.

We show that the cosmic star formation rate per comoving volume should exhibit a distinct drop around the reionization redshift, when the HII regions in the intergalactic medium around individual ionizing sources first overlapped. The drop results from the increase in the temperature of the intergalactic medium as it was photo-ionized, and the consequent suppression of the formation of low-mass galaxies. We show quantitatively that the detection of this drop, which is marked by a corresponding fall in the number counts of faint galaxies, should become feasible over the coming decade with the Next Generation Space Telescope.

The combination of deep exposures and high resolution offered by telescopes in space allows the detection of lensing over a wide range of source redshifts and lens masses. As an example, we model a lens candidate found in the southern Hubble Deep Field. The system consists of a source galaxy lensed into an arc, 0.9 arcseconds from an elliptical galaxy with an implied mass-to-light ratio of order 15 in solar units. This lens system may be a preview of the large number of lensed high-redshift galaxies that will be detected with the Next Generation Space Telescope (NGST). Popular models of structure formation imply that roughly 3% of observed sources at z>5 should be multiply imaged. Thus, NGST should detect several lensed objects in each field of view. In addition, NGST will be a valuable tool for weak lensing measurements if it can resolve the background galaxies. We estimate the angular size distribution of high redshift sources within hierarchical models of structure formation and find that most will be resolved. (abridged)

Direct observations of the first generation of luminous objects will likely become feasible over the next decade. The advent of the Next Generation Space Telescope (NGST) will allow imaging of numerous galaxies and mini-quasars at redshifts z>5. We apply semi-analytic models of structure formation to estimate the rate of multiple imaging of these sources by intervening gravitational lenses. Popular CDM models for galaxy formation yield a lensing optical depth of about 1% for sources at redshift 10. The expected slope of the luminosity function of the early sources implies an additional magnification bias of about 5, bringing the fraction of lensed sources at z=10 to about 5%. We estimate the angular size distribution of high-redshift disk galaxies and find that most of them are more extended than the resolution limit of NGST, roughly 0.06 arcseconds. We also show that there is only a modest redshift evolution in the mean surface brightness of galaxies at z>2. The expected increase by 1-2 orders of magnitude in the number of resolved sources on the sky, due to observations with NGST, will dramatically improve upon the statistical significance of existing weak lensing measurements. We show that, despite this increase in the density of sources, confusion noise from z>2 galaxies is expected to be small for NGST observations.

We examine models of the mass distribution for the first known case of gravitational lensing. Several new sets of constraints are used, based on recent observations. We remodel the VLBI observations of the radio jets in the two images of Q0957+561, showing that the previously derived positions and uncertainties were incorrect. We use as additional constraints the candidate lensed pairs of galaxies identified recently with the Hubble Space Telescope. We explore a wider range of lens models than before, and find that the Hubble constant is not tightly constrained once elliptical lens models are considered. We also discuss the systematic uncertainties caused by the cluster containing the lens galaxy. We conclude that additional observations of the candidate lensed galaxies as well as direct measurements of the cluster mass profile are needed before an accurate value of the Hubble constant can be derived from this lens.

During the period of reionization the Universe was filled with a cosmological background of ionizing radiation. By that time a significant fraction of the cosmic gas had already been incorporated into collapsed galactic halos with virial temperatures below about 10000 K that were unable to cool efficiently. We show that photoionization of this gas by the fresh cosmic UV background boiled the gas out of the gravitational potential wells of its host halos. We calculate the photoionization heating of gas inside spherically symmetric dark matter halos, and assume that gas which is heated above its virial temperature is expelled. In popular Cold Dark Matter models, the Press-Schechter halo abundance implies that about 50-90% of the collapsed gas was evaporated at reionization. The gas originated from halos below a threshold circular velocity of 10-15 km/s. The resulting outflows from the dwarf galaxy population at redshifts 5-10 affected the metallicity, thermal and hydrodynamic state of the surrounding intergalactic medium. Our results suggest that stellar systems with a velocity dispersion below about 10 km/s, such as globular clusters or the dwarf spheroidal galaxies of the Local Group, did not form directly through cosmological collapse at high redshifts.

We model an apparent gravitational lens system HDFS 2232509-603243 in the Hubble Deep Field South. The system consists of a blue V=25 mag arc separated by 0.9 arcsec from a red V=22 mag elliptical galaxy. A mass distribution which follows the observed light distribution with a constant mass-to-light ratio can fit the arc component positions if external shear is added. A good fit is also obtained with simple parameterized models, and all the models predict a forth image fainter than the detection limit. The inferred mass-to-light ratio is roughly 15 in solar units if the lens is at redshift 0.6. The lens models predict a velocity dispersion of about 280 km/s which could be confirmed with spectroscopy.

Because of their simplicity, axisymmetric mass distributions are often used to model gravitational lenses. Since galaxies are usually observed to have elliptical light distributions, mass distributions with elliptical density contours offer more general and realistic lens models. They are difficult to use, however, since previous studies have shown that the deflection angle (and magnification) in this case can only be obtained by rather expensive numerical integrations. We present a family of lens models for which the deflection can be calculated to high relative accuracy (one part in 100000) with a greatly reduced numerical effort, for small and large ellipticity alike. This makes it easier to use these distributions for modelling individual lenses as well as for applications requiring larger computing times, such as statistical lensing studies. A program implementing this method can be obtained from the author (or at http://www.sns.ias.edu/~barkana/ellip.html).

We present a new method for determining time delays among the light curves of various images in a gravitational lens. The method is based on constructing a simple model for the source variation and forming a chi squared measure of the agreement of this same variation with all of the lightcurves. While inspired by Press et al. (1992a, b) our approach is different since we do not assume a Gaussian process for the source variation. We apply our method to the light curves of the quadruple gravitational lens PG1115+080 measured by Schechter et al. (1997). Unlike Schechter et al. we include correlated measurement errors in the analysis, as well as the possibility that microlensing may cause different images to vary by different factors in flux. We find a value of 25.0 (+3.3,-3.8) days (95% confidence) for the delay between components B and C (close to the 24 day value of Schechter et al). However, the ratio t_AC/t_BA of the two intermediate delays is poorly determined at 1.13 (+.18,-.17) (68% confidence), close to the value predicted by lens models (1.4) unlike the Schechter et al. value (0.7). The variation ratios of C with respect to A and of A with respect to B are both different from 1, 1.39 (+.16,-.20) and .79 (+.10,-.12) (95% confidence), respectively. This is an indication of a microlensing gradient, and this type of microlensing may allow us to conclude that the size of the quasar optical emission region is about 1000 AU.