We present an analysis of the evolution of the Lyman-series forest into the epoch of reionization using cosmological radiative transfer simulations in a scenario where reionization ends late. We explore models with different midpoints of reionization and gas temperatures. We find that once the simulations have been calibrated to match the mean flux of the observed Lyman-$\alpha$ forest at $4 < z < 6$, they also naturally reproduce the distribution of effective optical depths of the Lyman-$\beta$ forest in this redshift range. We note that the tail of the largest optical depths that is most challenging to match corresponds to the long absorption trough of ULAS J0148+0600, which we have previously shown to be rare in our simulations. We consider the evolution of the Lyman-series forest out to higher redshifts, and show that future observations of the Lyman-$\beta$ forest at $z>6$ will discriminate between different reionization histories. The evolution of the Lyman-$\alpha$ and Lyman-$\gamma$ forests are less promising as a tool for pushing studies of reionization to higher redshifts due to the stronger saturation and foreground contamination, respectively.

We explore a potential new probe of reionization: the cross-correlation of high-redshift galaxies with maps of the thermal Sunyaev-Zel'dovich (tSZ) effect. We consider two types of high redshift galaxies: Lyman break galaxies (LBGs) and Lyman-$\alpha$ emitters (LAEs). LBGs and LAEs will be detected in large numbers at high redshift ($z \approx$ 4 to 7) by ongoing and future surveys. We consider a LBG sample like that expected from the The Rubin Observatory Legacy Survey of Space and Time (LSST), and a selection of LAEs modelled after the Subaru SILVERRUSH program, but covering a much larger sky fraction. The tSZ effect is sensitive to a line-of-sight integral of the ionized gas pressure, and can be measured across large patches of the sky using multi-frequency CMB surveys. We consider forecast tSZ maps from CMB Stage 4 and more futuristic observations. Using a suite of hydrodynamical simulations, we show that the high-redshift galaxies are correlated with the tSZ signal from reionization, with cross-power amplitude of order $10^{-15}$ at $\ell \sim 1000$. The cross-spectra between LBGs/LAEs with tSZ maps contain information about the reionization history of the Universe, such as the distribution of bubble sizes, and could in principle be used to directly measure the timing of reionization. The amplitude of the signal is small, however, and its detectability is hindered by low-redshift contributions to tSZ maps and by instrumental noise. If the the low-redshift contribution to the observed tSZ signal is suppressed by masking of massive halos, a combination of overlapping futuristic CMB and galaxy surveys could potentially probe this signal.

The search for primordial gravitational waves through the $B$-mode polarization pattern in the CMB is one of the major goals of current and future CMB experiments. Besides foregrounds, a potential hurdle in this search is the anisotropic secondary $B$-mode polarization generated by the scattering of CMB photons off free electrons produced during patchy cosmological reionization. Robust predictions of these secondary anisotropies are challenging because of uncertainties in the reionization history. In this paper, we revise estimates of the reionization-induced $B$-mode signal by incorporating recent advances in the understanding of reionization through observations of the Lyman-$\alpha$ forest. To derive these $B$-mode estimates, we use high-dynamic-range radiative transfer simulations of reionization that are calibrated to the Ly$\alpha$ data. These simulations are also consistent with a variety of other high-redshift observations. We find that around multipoles $\ell\approx 100$, reionization induces $B$-mode power with $\ell(\ell+1)C_\ell^{BB}/2\pi\approx 4\times 10^{-6}\,\mu$K$^2$. This secondary signal is thus at the level of the primordial signal with the tensor-to-scalar ratio $r<10^{-4}$, and can increase by a factor of $\sim 50$ if reionization is sourced by highly clustered sources residing in haloes with mass of $\sim 10^{11}$ M$_\odot$. Our findings suggest that the contribution of patchy reionization to the search for primordial gravitational waves is unlikely to be a concern for currently planned CMB experiments.

We compare a sample of five high-resolution, high S/N Ly$\alpha$ forest spectra of bright $6<z \lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile fitting procedure for the inverted transmitted flux, $1-F$, similar to the widely used Voigt profile fitting of the transmitted flux $F$ at lower redshifts, to characterise the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominant in low temeperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution ($\leq $ 8 km/s) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of $T_0= (11000 \pm 1600,10500\pm 2100,12000 \pm 2200)$ K at $z= (5.4,5.6,5.8)$. This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to $z\sim5.3$, the width distribution of the observed transmission spikes is consistent with the range of $T_0$ caused by spatial fluctuations in the temperature-density relation.

We model the 21-cm signal and LAE population evolution during the epoch of reionization in order to predict the 21cm-LAE cross-power spectrum. We employ high-dynamic-range simulations of the IGM to create models that are consistent with constraints from the CMB, Lyman-$\alpha$ forest and LAE population statistics. Using these models we consider the evolution of the cross-power spectrum for a selection of realistic reionization histories and predict the sensitivity of current and upcoming surveys to measuring this signal. We find that the imprint of a delayed-end to reionization can be observed by future surveys, and that strong constraints can be placed on the progression of reionization as late as $z=5.7$ using a Subaru-SKA survey. We make predictions for the signal-to-noise ratios achievable by combinations of Subaru/PFS with the MWA, LOFAR, HERA and SKA interferometers for an integration time of 1000 hours. We find that a Subaru-SKA survey could measure the cross-power spectrum for a late reionization at $z=6.6$ with a total signal-to-noise greater than 5, making it possible to constrain both the timing and bubble size at the end of reionization. Furthermore, we find that expanding the current Subaru/PFS survey area and depth by a factor of three would double the total signal-to-noise.

It is predicted that the spectrum of radio emission from the whole sky should show a dip arising from the action of the light from the first stars on the hydrogen atoms in the surrounding gas, which causes the 21-cm line to appear in absorption against the cosmic microwave background. Bowman et al. 2018 identified a broad flat-bottomed absorption profile centred at 78 MHz, which could be this feature, although the depth of the profile is much larger than expected. We have examined the modelling process they used and find that their data implies unphysical parameters for the foreground emission and also that their solution is not unique in the sense that we found other simple formulations for the signal that are different in shape but that also fit their data. We argue that this calls into question the interpretation of these data as an unambiguous detection of the cosmological 21-cm absorption signature.

A long (110 cMpc/$h$) and deep absorption trough in the Ly$\alpha$ forest has been observed extending down to redshift 5.5 in the spectrum of ULAS J0148+0600. Although no Ly$\alpha$ transmission is detected, Ly$\beta$ spikes are present which has led to claims that the gas along this trough must be ionized. Using high resolution cosmological radiative transfer simulations in large volumes, we show that in a scenario where reionization ends late ($z \sim 5.2$), our simulations can reproduce troughs as long as observed. In this model, we find that the troughs are caused by islands of neutral hydrogen. Small ionized holes within the neutral islands allow for the transmission of Ly$\beta$. We have also modelled the Ly$\alpha$ emitter population around the simulated troughs, and show that there is a deficit of Ly$\alpha$ emitters close to the trough as is observed.

Determinations of the UV luminosity function of AGN at high redshifts are important for constraining the AGN contribution to reionization and understanding the growth of supermassive black holes. Recent inferences of the luminosity function suffer from inconsistencies arising from inhomogeneous selection and analysis of AGN data. We address this problem by constructing a sample of more than 80,000 colour-selected AGN from redshift z=0 to 7.5. While this sample is composed of multiple data sets with spectroscopic redshifts and completeness estimates, we homogenise these data sets to identical cosmologies, intrinsic AGN spectra, and magnitude systems. Using this sample, we derive the AGN UV luminosity function from redshift z=0 to 7.5. The luminosity function has a double power law form at all redshifts. The break magnitude $M_*$ of the AGN luminosity function shows a steep brightening from $M_*\sim -24$ at z=0.7 to $M_*\sim -29$ at z=6. The faint-end slope $\beta$ significantly steepens from $-1.7$ at $z<2.2$ to $-2.4$ at $z\simeq 6$. In spite of this steepening, the contribution of AGN to the hydrogen photoionization rate at $z\sim 6$ is subdominant (< 3%), although it can be non-negligible (~10%) if these luminosity functions hold down to $M_{1450}=-18$. Under reasonable assumptions, AGN can reionize HeII by redshift z=2.9. At low redshifts (z<0.5), AGN can produce about half of the hydrogen photoionization rate inferred from the statistics of HI absorption lines in the IGM. Our global analysis of the luminosity function also reveals important systematic errors in the data, particularly at z=2.2--3.5, which need to be addressed and incorporated in the AGN selection function in future in order to improve our results. We make various fitting functions, luminosity function analysis codes, and homogenised AGN data publicly available.

Mapping the intergalactic medium (IGM) in Lyman-$\alpha$ emission would yield unprecedented tomographic information on the large-scale distribution of baryons and potentially provide new constraints on the UV background and various feedback processes relevant for galaxy formation. Here, we examine the Lyman-$\alpha$ emission of the moderate-density IGM due to collisional excitations and recombinations in the presence of a UV background in the Sherwood simulation suite. We focus on large-scale-structure filaments in which Lyman-$\alpha$ radiative transfer effects are expected to be moderate. At low density the emission is primarily due to fluorescent re-emission of the ionising UV background due to recombinations, while collisional excitations dominate at higher densities. We discuss prospects of current and future observational facilities to detect this emission and focus on VLT/MUSE for a more detailed sensitivity analysis. We construct mock MUSE observations resembling the MUSE Hubble Deep Field South from our simulations and show that our predictions are consistent with recent analyses of diffuse Lyman-$\alpha$ emission using MUSE at redshifts $3 < z < 6$. We find that it should be most feasible to detect the Lyman-$\alpha$ emission from filaments in the IGM in overdense regions, somewhat surprisingly towards the high-redshift end ($z \gtrsim 4.5$) accessible by MUSE, and if narrowband widths close to the MUSE spectral resolution limit are used. This is due to the higher densities and lower temperatures in the IGM at higher redshift. High-redshift protoclusters therefore appear to be the ideal environment to observe filamentary structures in the IGM in Lyman-$\alpha$ emission.

Observations of the high-redshift Universe using the 21 cm line of neutral hydrogen and complimentary emission lines from the first galaxies promise to open a new door for our understanding of the epoch of reionization. We present predictions for the [C II] 158-micron line and H I 21 cm emission from redshifts z=6--9 using high-dynamic-range cosmological simulations combined with semi-analytical models. We find that the CONCERTO experiment should be able to detect the large scale power spectrum of [C II] emission to redshifts of up to z=8 (signal-to-noise ratio ~ 1 at k = 0.1 h/cMpc with 1500 hr of integration). A Stage II experiment similar to CCAT-p should be able to detect [C II] from even higher redshifts to high significance for similar integration times (signal-to-noise ratio of ~50 at k = 0.2 h/cMpc at z=6--9). We study the possibility of combining such future [C II] measurements with 21 cm measurements using LOFAR and SKA to measure the [C II]-21cm cross power spectra, and find that a Stage II experiment should be able to measure the cross-power spectrum for k < 1 h/cMpc to signal-to-noise ratio of better than 10. We discuss the capability of such measurements to constrain astrophysical parameters relevant to reionization and show that a measurement of the [C II]-21cm cross power spectrum helps break the degeneracy between the mass and brightness of ionizing sources.

High-redshift QSO spectra show large spatial fluctuations in the Ly-alpha opacity of the intergalactic medium on surprisingly large scales at z>~5.5. We present a radiative transfer simulation of cosmic reionization driven by galaxies that reproduces this large scatter and the rapid evolution of the Ly-alpha opacity distribution at 5<z<6. The simulation also reproduces the low Thomson scattering optical depth reported by the latest CMB measurement and is consistent with the observed short near-zones and strong red damping wings in the highest-redshift QSOs. It also matches the rapid disappearance of observed Ly-alpha emission by galaxies at z>~6. Reionization is complete at z=5.3 in our model, and 50% of the volume of the Universe is ionized at z=7. Agreement with the Ly-alpha forest data in such a late reionization model requires a rapid evolution of the ionizing emissivity of galaxies that peaks at z~6.8. The late end of reionization results in a large scatter in the photoionisation rate and the neutral hydrogen fraction at redshifts as low as z<~5.5 with large residual neutral 'islands' that can produce very long Gunn-Peterson troughs resembling those seen in the data.

We model the high redshift (z > 5) Lyman-$\alpha$ emitting (LAE) galaxy population using the empirical rest-frame equivalent width distribution. We calibrate to the observed luminosity function and angular correlation function at z = 5.7 as measured by the SILVERRUSH survey. This allows us to populate the high-dynamic-range Sherwood simulation suite with LAEs, and to calculate the transmission of their Ly $\alpha$ emission through the inter-galactic medium (IGM). We use this simulated population to explore the effect of the IGM on high-redshift observations of LAEs, and make predictions for the narrowband filter redshifts at z = 6.6, 7.0 and 7.3. Comparing our model with existing observations, we find a late reionization is suggested, consistent with the recent low optical depth derived from the cosmic microwave background (CMB) by the Planck collaboration and the opacity fluctuations in the Ly $\alpha$ forest. We also explore the role of the circum-galactic medium (CGM) and the larger volume of gas which is infalling onto the host halo versus the IGM in attenuating the Ly $\alpha$ signal, finding that a significant fraction of the attenuation is due to the CGM and infalling gas, which increases towards the end of reionization, albeit with a large scatter across the mock LAE population.

We model the transmission of the Lyman-alpha line through the circum- and intergalactic media around dark matter haloes expected to host Lyman-alpha emitters (LAEs) at z > 5.7, using the high-dynamic-range Sherwood simulations. We find very different CGM environments around more massive haloes (~10^11 M_sun) compared to less massive haloes (~10^9 M_sun) at these redshifts, which can contribute to a different evolution of the Lyman-alpha transmission from LAEs within these haloes. Additionally we confirm that part of the differential evolution could result from bright LAEs being more likely to reside in larger ionized regions. We conclude that a combination of the CGM environment and the IGM ionization structure is likely to be responsible for the differential evolution of the bright and faint ends of the LAE luminosity function at z > 6. More generally, we confirm the suggestion that the self-shielded neutral gas in the outskirts of the host halo can strongly attenuate the Lyman-alpha emission from high redshift galaxies. We find that this has a stronger effect on the more massive haloes hosting brighter LAEs. The faint-end of the LAE luminosity function is thus a more reliable probe of the average ionization state of the IGM. Comparing our model for LAEs with a range of observational data we find that the favoured reionization histories are our previously advocated `Late' and `Very Late' reionization histories, in which reionization finishes rather rapidly at around z ~ 6.

We investigate self-shielding of intergalactic hydrogen against ionizing radiation in radiative transfer simulations of cosmic reionization carefully calibrated with Lyman alpha forest data. While self-shielded regions manifest as Lyman-limit systems in the post-reionization Universe, here we focus on their evolution during reionization (redshifts z=6-10). At these redshifts, the spatial distribution of hydrogen-ionizing radiation is highly inhomogeneous, and some regions of the Universe are still neutral. After masking the neutral regions and ionizing sources in the simulation, we find that the hydrogen photoionization rate depends on the local hydrogen density in a manner very similar to that in the post-reionization Universe. The characteristic physical hydrogen density above which self-shielding becomes important at these redshifts is about $\mathrm{n_H \sim 3 \times 10^{-3} cm^{-3}}$, or $\sim$ 20 times the mean hydrogen density, reflecting the fact that during reionization photoionization rates are typically low enough that the filaments in the cosmic web are often self-shielded. The value of the typical self-shielding density decreases by a factor of 3 between redshifts z=3 and 10, and follows the evolution of the average photoionization rate in ionized regions in a simple fashion. We provide a simple parameterization of the photoionization rate as a function of density in self-shielded regions during the epoch of reionization.

We study a novel electromagnetic signature of supermassive black hole binaries whose inspiral starts being dominated by gravitational wave (GW) emission. Recent simulations suggest that the binary's member BHs can continue to accrete gas from the circumbinary accretion disk in this phase of the binary's evolution, all the way until coalescence. If one of the binary members produces a radio jet as a result of accretion, the jet precesses along a biconical surface due to the binary's orbital motion. When the binary enters the GW phase of its evolution, the opening angle widens, the jet exhibits milliarcsecond scale wiggles, and the conical surface of jet precession is twisted due to apparant superluminal motion. The rapidly increasing orbital velocity of the binary gives the jet an appearance of a "chirp." This helical chirping morphology of the jet can be used to infer the binary parameters. For binaries with mass 10^7--10^10 Msun at redshifts z<0.5, monitoring these features in current and archival data will place a lower limit on sources that could be detected by eLISA and Pulsar Timing Arrays. In the future, microarcsecond interferometry with the Square Kilometer Array will increase the potential usefulness of this technique.

The neutral hydrogen (HI) content of dark matter haloes forms an intermediate state in the baryon cycle that connects the hot shock-heated gas and cold star-forming gas in haloes. Measurement of the relationship between HI mass and halo mass therefore puts important constraints on galaxy formation models. We combine radio observations of HI in emission at low redshift ($z\sim 0$) with optical/UV observations of HI in absorption at high redshift ($1<z<4$) to derive constraints on the evolution of the HI-mass halo-mass (HIHM) relation from redshift $z=4$ to $z=0$. We find that one can model the HIHM relation similar to the stellar-mass halo-mass (SHM) relation at $z \sim 0$. At $z=0$, haloes with mass $10^{11.7}$ M$_\odot$ have the highest HI mass fraction ($\sim 1\%$), which is about four times smaller than their stellar mass fraction. We model the evolution of the HIHM relation in a manner similar to that of the SHM relation. Combining this parameterisation with a redshift- and mass-dependent modified Navarro-Frenk-White (NFW) profile for the HI density within a halo, we draw constraints on the evolution of the HIHM relation from the observed HI column density, incidence rate, and clustering bias at high redshift. We compare these findings with results from hydrodynamical simulations and other approaches in the literature, and find the models to be consistent with each other at the $68\%$ confidence level.

Using the Eris zoom-in cosmological simulation of assembly of a Milky Way analog, we study chemical enrichment of stars due to accretion of metal-enriched gas from the interstellar medium during the Galaxy's development. We consider metal-poor and old stars in both Galactic halo and bulge and make use of stellar orbits, gas density and metallicity distributions in Eris. Assuming spherically symmetric Bondi-Hoyle accretion, we find that halo and bulge stars accrete metals at the rate of about 10^-24 solar mass per year and 10^-22 solar mass per year, respectively, at redshifts z < 3, but this accretion rate increases hundred-fold to about 10^-20 solar mass per year at higher redshifts due to increased gas density. Bulge and halo stars accrete similar amounts of metals at high redshifts as kinematically distinct bulge and halo are not yet developed at these redshifts and both sets of stars encounter similar metal distribution in the ISM on average. Accretion alone can enrich main-sequence stars up to [Fe/H] -2 in extreme cases.Median enrichment level due to accretion in these stars is about [Fe/H]~-6 to -5.Because accretion mostly takes place at high redshifts, it is alpha-enriched to [alpha/Fe]~0.5. We find that accretive metal enrichment is significant enough to affect the predicted metallicity distribution function of halo stars at [Fe/H] < -5.This suggests that attempts to infer the natal chemical environment of the most metal-poor stars from their observed enrichment today can be hindered due to metal accretion. Peculiar enrichment patterns such as those predicted to arise from pair-instability supernovae could help in disentangling natal and accreted metal content of stars.

We present here predictions for the spatial distribution of 21 cm brightness temperature fluctuations from high-dynamic-range simulations for AGN-dominated reionization histories that have been tested against available Lyman-alpha and CMB data. We model AGN by extrapolating the observed M-sigma relation to high redshifts and assign them ionizing emissivities consistent with recent UV luminosity function measurements. We assess the observability of the predicted spatial 21 cm fluctuations by ongoing and upcoming experiments in the late stages of reionization in the limit in which the hydrogen 21 cm spin temperature is significantly larger than the CMB temperature. Our AGN-dominated reionization histories increase the variance of the 21 cm emission by a factor of up to ten compared to similar reionization histories dominated by faint galaxies, to values close to 100 mK^2 at scales accessible to experiments (k < 1 h/cMpc). This is lower than the sensitivity claimed to have been already reached by ongoing experiments by only a factor of about two or less. When reionization is dominated by AGN, the 21 cm power spectrum is enhanced on all scales due to the enhanced bias of the clustering of the more massive haloes and the peak in the large scale 21 cm power is strongly enhanced and moved to larger scales due to bigger characteristic bubble sizes. AGN dominated reionization should be easily detectable by LOFAR (and later HERA and SKA1) at their design sensitivity, assuming successful foreground subtraction and instrument calibration. Conversely, these could become the first non-trivial reionization scenarios to be ruled out by 21 cm experiments, thereby constraining the contribution of AGN to reionization.

The distribution of diffuse gas in the intergalactic medium (IGM) imprints a series of hydrogen absorption lines on the spectra of distant background quasars known as the Lyman-$\alpha$ forest. Cosmological hydrodynamical simulations predict that IGM density fluctuations are suppressed below a characteristic scale where thermal pressure balances gravity. We measured this pressure-smoothing scale by quantifying absorption correlations in a sample of close quasar pairs. We compared our measurements to hydrodynamical simulations, where pressure smoothing is determined by the integrated thermal history of the IGM. Our findings are consistent with standard models for photoionization heating by the ultraviolet radiation backgrounds that reionized the universe.

We present here 21 cm predictions from high dynamic range simulations for a range of reionization histories that have been tested against available Lyman-alpha and CMB data. We assess the observability of the predicted spatial 21 cm fluctuations by ongoing and upcoming experiments in the late stages of reionization in the limit in which the hydrogen spin temperature is significantly larger than the CMB temperature. Models consistent with the available Lyman-alpha data and CMB measurement of the Thomson optical depth predict typical values of 10--20 mK^2 for the variance of the 21 cm brightness temperature at redshifts z=7--10 at scales accessible to ongoing and upcoming experiments (k < 1 h/cMpc). This is only a factor of a few below the sensitivity claimed to have been already reached by ongoing experiments in the signal rms value. Our different models for the reionization history make markedly different predictions for the redshift evolution and thus frequency dependence of the 21 cm power spectrum and should be easily discernible by LOFAR (and later HERA and SKA1) at their design sensitivity. Our simulations have sufficient resolution to assess the effect of high-density Lyman limit systems that can self-shield against ionizing radiation and stay 21 cm bright even if the hydrogen in their surroundings is highly ionized. Our simulations predict that including the effect of the self-shielded gas in highly ionized regions reduces the large scale 21 cm power by about 30%.

The thermal state of the intergalactic medium (IGM) at z < 6 constrains the nature and timing of cosmic reionization events, but its inference from the Ly-alpha forest is degenerate with the 3-D structure of the IGM on ~100 kpc scales, where, analogous to the classical Jeans argument, the pressure of the T~$10^4$ K gas supports it against gravity. We simulate the IGM using smoothed particle hydrodynamics, and find that, at z < 6, the gas density power spectrum does not exhibit the expected Jeans filtering cutoff, because dense gas in collapsed halos dominates the small-scale power masking pressure smoothing effects. We introduce a new statistic, the real-space Ly-alpha flux, $F_\mathrm{real}$, which naturally suppresses dense gas, and is thus robust against the poorly understood physics of galaxy formation, revealing pressure smoothing in the diffuse IGM. The $F_\mathrm{real}$ power spectrum is accurately described by a simple fitting function with cutoff at $\lambda_F$, allowing us to rigorously quantify the filtering scale for the first time: we find $\lambda_F$ = 79 kpc (comoving) at z=3 for our fiducial thermal model. This statistic has the added advantage that it directly relates to observations of correlated Ly-alpha forest absorption in close quasar pairs, recently proposed as a method to measure the filtering scale. Our results enable one to quantify the filtering scale in simulations, and ask meaningful questions about its dependence on reionization and thermal history. Accordingly, the standard description of the IGM in terms of the amplitude $T_0$ and slope $\gamma$ of the temperature-density relation $T = T_0\Delta^{\gamma-1}$ should be augmented with a third filtering scale parameter $\lambda_F$.

Recent studies have highlighted that galaxies at z = 6-8 fall short of producing enough ionizing photons to reionize the IGM, and suggest that Population III stars could resolve this tension, because their harder spectra can produce ~10x more ionizing photons than Population II. But this argument depends critically on the duration of the Population III era, and because Population III stars form from pristine gas, in turn depends on the rate of galactic enrichment. We use a semi-analytic model of galaxy formation which tracks galactic chemical evolution, to gauge the impact of Population III stars on reionization. Population III SNe produce distinct metal abundances, and we argue that the duration of the Population III era can be constrained by precise relative abundance measurements in high-z damped Ly{\alpha} absorbers (DLAs), which provide a chemical record of past star-formation. We find that a single generation of Population III stars can self-enrich galaxies above the critical metallicity Zcrit=10^-4 Zsun for the Population III-to-II transition, on a very short timescale of ~10^6 yr, owing to the large metal yields and short lifetimes of Population III stars. This subsequently terminates the Population III era, hence they contribute >~ 50% of the ionizing photons only for z >~ 30, and at z=10 contribute <1%. The Population III contribution can be increased by delaying metal mixing into the ISM. However comparing the resulting metal abundance pattern to existing measurements in z <~ 6 DLAs, we show that the fractional contribution of high-mass Population III stars to the ionization rate must be <~ 10% at z = 10. Future abundance measurements of z~7-8 QSOs and GRBs should probe the era when the chemical vestiges of Population III star formation become detectable.

Damped Ly-alpha absorbers (DLAs) can be used to measure gas-phase metallicities at large cosmological lookback times with high precision. Relative abundances can still be measured accurately deep into the reionization epoch (z > 6) using transitions redward of Ly-alpha. Here we study chemical evolution of DLAs using a constrained model for evolution of galaxies and IGM to determine the degree to which DLA abundance measurements can probe Population III enrichment. We find that if the critical metallicity of Population III to II transition is < 10^-4 Zsun, the cosmic Population III SFR is zero for z<8. Nevertheless, at high redshift (z ~ 6) Population III chemical signatures are retained in low-mass galaxies (halo mass < 10^9 Msun). This is because photoionization feedback suppresses star formation in these galaxies until relatively low redshift (z ~ 10), and the chemical record of early Population III star formation is retained. We model DLAs as these galaxies by assigning to them a mass-dependent H I absorption cross section and predict distribution of DLA abundance ratios. We find that these distributions are anchored towards abundance ratios set by Population II yields, but exhibit a tail that depends on the Population III IMF for z > 5. Thus, a sample of DLA abundance measurements at high redshift holds the promise to constrain Population III IMF. A sample of just 10 DLAs with relative abundances measured to an accuracy of 0.1 dex is sufficient to constrain the Population III IMF at 4-sigma. These constraints may prove stronger than other probes such as metal-poor stars and individual DLAs. Our results provide a global picture of the cosmic thermal, ionization, and chemical evolution, and can rule out certain Population III scenarios.

We derive an estimate of the rate of formation of dark matter halos per unit volume as a function of the halo mass and redshift of formation. Analytical estimates of the number density of dark matter halos are useful in modeling several cosmological phenomena. We use the excursion set formalism for computing the formation rate of dark matter halos. We use an approach that allows us to differentiate between major and minor mergers, as this is a pertinent issue for semi-analytic models of galaxy formation. We compute the formation rate for the Press-Schechter and the Sheth-Tormen mass function. We show that the formation rate computed in this manner is positive at all scales. We comment on the Sasaki formalism where negative halo formation rates are obtained. Our estimates compare very well with N-Body simulations for a variety of models. We also discuss the halo survival probability and the formation redshift distributions using our method.

We examine the formation of groups of multiple supermassive black holes (SMBHs) in gas-poor galactic nuclei due to the high merger rate of galaxies at high redshifts. We calculate the relative likelihood of binary, triple, and quadruple SMBH systems, by considering the timescales for relevant processes and combining merger trees with N-body simulations for the dynamics of stars and SMBHs in galactic nuclei. Typical haloes today with mass $M_0\approx 10^{14}$ M$_\odot$ have an average mass $M_{z=6}=5\times 10^{11}$ M$_\odot$ at $z\sim 6$, while rare haloes with current mass $M_0\gtrsim 10^{15}$ M$_\odot$ have an average mass $M_{z=6}=5\times 10^{12}$ M$_\odot$ at that redshift. These cluster-size haloes are expected to host single galaxies at $z\sim 6$. We expect about 30% galaxies within haloes with present-day mass $M_0\approx 10^{14}$ M$_\odot$ to contain more than two SMBHs at redshifts $2\lesssim z\lesssim 6$. For larger present-day haloes, with $M_0\gtrsim 10^{15}$ M$_\odot$, this fraction is almost 60%. The existence of multiple SMBHs at high redshifts can potentially explain the mass deficiencies observed in the cores of massive elliptical galaxies, which are up to 5 times the mass of their central BHs. Multiple SMBHs would also lead to an enhanced rate of tidal disruption of stars, modified gravitational wave signals compared to isolated BH binaries, and slingshot ejection of SMBHs from galaxies at high speeds in excess of 2000 km s$^{-1}$.

Observations of galaxy luminosity function at high redshifts typically focus on fields of view of limited sizes preferentially containing bright sources. These regions possibly are overdense and hence biased with respect to the globally averaged regions. Using a semi-analytic model based on Choudhury & Ferrara (2006) which is calibrated to match a wide range of observations, we study the reionization and thermal history of the universe in overdense regions. The main results of our calculation are: (i) Reionization and thermal histories in the biased regions are markedly different from the average ones because of enhanced number of sources and higher radiative feedback. (ii) The galaxy luminosity function for biased regions is markedly different from those corresponding to average ones. In particular, the effect of radiative feedback arising from cosmic reionization is visible at much brighter luminosities. (iii) Because of the enhanced radiative feedback within overdense locations, the luminosity function in such regions is more sensitive to reionization history than in average regions. The effect of feedback is visible for absolute AB magnitude $M_{AB} \gtrsim -17$ at $z=8$, almost within the reach of present day observations and surely to be probed by JWST. This could possibly serve as an additional probe of radiative feedback and hence reionization at high redshifts.

Filaments are one of the most prominent features visible in the galaxy distribution. Considering the Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey Data Release Seven (SDSS DR7), we have analyzed the filamentarity in 11 nearly two dimensional (2D) sections through a volume limited subsample of this data. The galaxy distribution, we find, has excess filamentarity in comparison to a random distribution of points. We use a statistical technique "Shuffle" to determine $L_{\rm MAX}$, the largest length-scale at which we have statistically significant filaments. We find that $L_{\rm MAX}$ varies in the range $100-130 \, h^{-1} {\rm Mpc}$ across the 11 slices, with a mean value $L_{\rm MAX}=110 \pm 12 \, h^{-1}{\rm Mpc}$. Longer filaments, though possibly present in our data, are not statistically significant and are the outcome of chance alignments.

Press-Schechter theory gives a simple, approximate functional form of the mass function of dark matter haloes. Sheth and Tormen (ST) refined this mass function to give an improved analytical fit to results of N-body simulations. These forms of the halo mass function are universal (independent of cosmology and power spectrum) when scaled in suitable variables. Using large suites of LCDM N-body simulations, studies in the last few years have shown that this universality is only approximate. We explore whether some of the deviations from universality can be attributed to the power spectrum by computing the mass function in N-body simulations of various scale-free models in an Einstein-de Sitter cosmology. This choice of cosmology does not introduce any scale into the problem. These models have the advantage of being self-similar, hence stringent checks can be imposed while running these simulations. This set of numerical experiments is designed to isolate any power spectrum dependent departures from universality of mass functions. We show explicitly that the best fit ST parameters have a clear dependence on power spectrum. Our results also indicate that an improved analytical theory with more parameters is required in order to provide better fits to the mass function.

The epoch of reionization and formation of first stars are inter-linked topics that are of considerable interest. We use a simplified approach for studying formation of stars in collapsed haloes, and the resulting ionization of the inter-galactic medium (IGM). We consider a set of LCDM models allowed by observations of CMB temperature and polarization anisotropies for this study. We constrain parameters related to star formation with the help of observations. We constrain subsets of these parameters independently by using the observed metallicity of the inter-galactic medium at z ~ 5 and the requirement that the Thomson scattering optical depth due to an ionized IGM as determined for the model from CMB observations be reproduced. We consider a range of initial metalicities for star forming gas, and some variations of the initial mass function of stars. We find that a "normal" initial mass function (IMF) may satisfy these two constraints with a raised efficiency of star formation as compared to that seen in the local universe. Observations require a significant fraction of metals to escape from haloes to the IGM. We can also place constraints on the ratio of escape fraction for metals and ionizing photons, we find that this ratio is of order unity for most models. This highlights the importance of using the constraints arising from enrichment of the inter-galactic medium. (Abridged)

We present measurements of the redshift-space three-point correlation function of 50,967 Luminous Red Galaxies (LRGs) from Data Release 3 (DR3) of the Sloan Digital Sky Survey (SDSS). We have studied the shape dependence of the reduced three-point correlation function (Qz(s,q,theta)) on three different scales, s=4, 7 and 10 h-1 Mpc, and over the range of 1 < q < 3 and 0 < theta < 180. On small scales (s=4 h-1 Mpc), Qz is nearly constant, with little change as a function of q and theta. However, there is evidence for a shallow U-shaped behaviour (with theta) which is expected from theoretical modeling of Qz . On larger scales (s=7 and 10 h-1 Mpc), the U-shaped anisotropy in Qz (with theta) is more clearly detected. We compare this shape-dependence in Qz(s,q,theta) with that seen in mock galaxy catalogues which were generated by populating the dark matter halos in large N-body simulations with mock galaxies using various Halo Occupation Distributions (HOD). We find that the combination of the observed number density of LRGs, the (redshift-space) two-point correlation function and Qz provides a strong constraint on the allowed HOD parameters (M_min, M_1, alpha) and breaks key degeneracies between these parameters. For example, our observed Qz disfavors mock catalogues that overpopulate massive dark matter halos with many LRG satellites. We also estimate the linear bias of LRGs to be b=1.87+/-0.07 in excellent agreement with other measurements.

We present measurements of the normalised redshift-space three-point correlation function (Q_z) of galaxies from the Sloan Digital Sky Survey (SDSS) main galaxy sample. We have applied our "npt" algorithm to both a volume-limited (36738 galaxies) and magnitude-limited sample (134741 galaxies) of SDSS galaxies, and find consistent results between the two samples, thus confirming the weak luminosity dependence of Q_z recently seen by other authors. We compare our results to other Q_z measurements in the literature and find it to be consistent within the full jack-knife error estimates. However, we find these errors are significantly increased by the presence of the ``Sloan Great Wall'' (at z ~ 0.08) within these two SDSS datasets, which changes the 3-point correlation function (3PCF) by 70% on large scales (s>=10h^-1 Mpc). If we exclude this supercluster, our observed Q_z is in better agreement with that obtained from the 2dFGRS by other authors, thus demonstrating the sensitivity of these higher-order correlation functions to large-scale structures in the Universe. This analysis highlights that the SDSS datasets used here are not ``fair samples'' of the Universe for the estimation of higher-order clustering statistics and larger volumes are required. We study the shape-dependence of Q_z(s,q,theta) as one expects this measurement to depend on scale if the large scale structure in the Universe has grown via gravitational instability from Gaussian initial conditions. On small scales (s <= 6h^-1 Mpc), we see some evidence for shape-dependence in Q_z, but at present our measurements are consistent with a constant within the errors (Q_z ~ 0.75 +/- 0.05). On scales >10h^-1 Mpc, we see considerable shape-dependence in Q_z.

We present measurements of the redshift--space three-point correlation function of galaxies in the Sloan Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy morphology, color and luminosity. In order to control systematics due to selection effects, we used $r$--band, volume-limited samples of galaxies, constructed from the magnitude-limited SDSS data ($14.5<r<17.5$), and further divided the samples into two morphological types (early and late) or two color populations (red and blue). The three-point correlation function of SDSS galaxies follow the hierarchical relation well and the reduced three-point amplitudes in redshift--space are almost scale-independent ($Q_z=0.5\sim1.0$). In addition, their dependence on the morphology, color and luminosity is not statistically significant. Given the robust morphological, color and luminosity dependences of the two-point correlation function, this implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that simple deterministic linear relation with the underlying mass could not explain our measurements on these scales.