We use the Separate Universe technique to calibrate the dependence of linear and quadratic halo bias $b_1$ and $b_2$ on the local cosmic web environment of dark matter haloes. We do this by measuring the response of halo abundances at fixed mass and cosmic web tidal anisotropy $\alpha$ to an infinite wavelength initial perturbation. We augment our measurements with an analytical framework developed in earlier work which exploits the near-Lognormal shape of the distribution of $\alpha$ and results in very high precision calibrations. We present convenient fitting functions for the dependence of $b_1$ and $b_2$ on $\alpha$ over a wide range of halo mass for redshifts $0\leq z\leq1$. Our calibration of $b_2(\alpha)$ is the first demonstration to date of the dependence of non-linear bias on the local web environment. Motivated by previous results which showed that $\alpha$ is the primary indicator of halo assembly bias for a number of halo properties beyond halo mass, we then extend our analytical framework to accommodate the dependence of $b_1$ and $b_2$ on any such secondary property which has, or can be monotonically transformed to have, a Gaussian distribution. We demonstrate this technique for the specific case of halo concentration, finding good agreement with previous results. Our calibrations will be useful for a variety of halo model analyses focusing on galaxy assembly bias, as well as analytical forecasts of the potential for using $\alpha$ as a segregating variable in multi-tracer analyses.

The non-linear cosmic web environment of dark matter haloes plays a major role in shaping their growth and evolution, and potentially also affects the galaxies that reside in them. We develop an analytical (halo model) formalism to describe the tidal field of anisotropic halo-centric density distributions, as characterised by the halo-centric tidal tensor $\langle T_{ij} \rangle(<R)$ spherically averaged on scale $R\sim4R_{\rm vir}$ for haloes of virial radius $R_{\rm vir}$. We focus on axisymmetric anisotropies, which allows us to explore simple and intuitive toy models of (sub)halo configurations that exemplify some of the most interesting anisotropies in the cosmic web. We build our models around the spherical Navarro-Frenk-White (NFW) profile after describing it as a Gaussian mixture, which leads to almost fully analytical expressions for the `tidal anisotropy' scalar $\alpha(<4R_{\rm vir})$ extracted from the tidal tensor. Our axisymmetric examples include (i) a spherical halo at the axis of a cylindrical filament, (ii) an off-centred satellite in a spherical host halo and (iii) an axisymmetric halo. Using these, we demonstrate several interesting results. For example, the tidal tensor at the axis of a pure cylindrical filament gives $\alpha^{\rm (fil)}(<R)=1/2$ exactly, for any $R$. Also, $\alpha(<4R_{\rm vir,sat})$ for a satellite of radius $R_{\rm vir,sat}$ as a function of its host-centric distance is a sensitive probe of dynamical mass loss of the satellite in its host environment. Finally, we discuss a number of potentially interesting extensions and applications of our formalism that can deepen our understanding of the multi-scale phenomenology of the cosmic web.

We study the Voronoi volume function (VVF) -- the distribution of cell volumes (or inverse local number density) in the Voronoi tessellation of any set of cosmological tracers (galaxies/haloes). We show that the shape of the VVF of biased tracers responds sensitively to physical properties such as halo mass, large-scale environment, substructure and redshift-space effects, making this a hitherto unexplored probe of both primordial cosmology and galaxy evolution. Using convenient summary statistics -- the width, median and a low percentile of the VVF as functions of average tracer number density -- we explore these effects for tracer populations in a suite of N-body simulations of a range of dark matter models. Our summary statistics sensitively probe primordial features such as small-scale oscillations in the initial matter power spectrum (as arise in models involving collisional effects in the dark sector), while being largely insensitive to a truncation of initial power (as in warm dark matter models). For vanilla cold dark matter (CDM) cosmologies, the summary statistics display strong evolution and redshift-space effects, and are also sensitive to cosmological parameter values for realistic tracer samples. Comparing the VVF of galaxies in the GAMA survey with that of abundance matched CDM (sub)haloes tentatively reveals environmental effects in GAMA beyond halo mass (modulo unmodelled satellite properties). Our exploratory analysis thus paves the way for using the VVF as a new probe of galaxy evolution physics as well as the nature of dark matter and dark energy.

Ly$\alpha$ forest absorption spectra decomposed into multiple Voigt profile components (clouds) allow us to study clustering of intergalactic medium (IGM) as a function of HI column density ($N_{\rm HI}$). Here, we explore the transverse three-point correlation ($\zeta$) of these Ly$\alpha$ clouds using mock triplet spectra obtained from hydrodynamical simulations at $z \sim 2$ on scales of 1-5 $h^{-1}$cMpc. We find $\zeta$ to depend strongly on $N_{\rm HI}$ and scale and weakly on angle ($\theta$) of the triplet configuration. We show that the "hierarchical ansatz" is applicable for scales $\ge~ 3h^{-1}$cMpc, and obtain a median reduced three-point correlation (Q) in the range 0.2-0.7. We show, $\zeta$ is influenced strongly by the thermal and ionization state of the gas. As found in the case of galaxies, the influence of physical parameters on Q is weaker compared to that of $\zeta$. We show difference in $\zeta$ and Q between different simulations are minimized if we use appropriate $N_{\rm HI}$ cut-offs corresponding to a given baryon over-density ($\Delta$) using the measured $N_{\rm HI}~vs~\Delta$ relationship obtained from individual simulations. Additionally, we see the effect of pressure broadening on $\zeta$ in a model with artificially boosted heating rates. However, for models with realistic thermal and ionization histories the effect of pressure broadening on $\zeta$ is weak and sub-dominant compared to other local effects. We find strong redshift evolution shown by $\zeta$, mainly originating from the redshift evolution of thermal and ionization state of the IGM. We discuss the observational requirements for the detection of three-point correlation, specifically, in small intervals of configuration parameters and redshift.

We present a computationally efficient and fast semi-numerical technique for simulating the Lyman-$\alpha$ (Ly$\alpha$) absorption optical depth in presence of neutral hydrogen "islands" leftover from reionization at redshifts $5 \lesssim z \lesssim 6$. The main inputs to the analysis are (i) a semi-numerical photon-conserving model of ionized regions during reionization (named SCRIPT) and (ii) the fluctuating Gunn-Peterson approximation to model the Ly$\alpha$ absorption. Our model is then used for simulating the large-scale fluctuations in the effective optical depth as observed along sight lines towards high$-z$ quasars. We find that the model, with fiducial choice of free parameters, is in agreement with the results from more detailed radiative transfer simulations. By varying the free parameters, we obtain the constraints on reionization history at $5 \lesssim z \lesssim 6$ as allowed by the data. We confirm that reionization is not complete before $z \sim 5.6$ and the completion can be as late as $z \sim 5.2$ (both within $2-\sigma$), the exact limits depending on how the non-detections of the flux in the data are treated. With further improvements in the model and with more sight lines at $z \sim 6$, we can take advantage of the computational efficiency of our analysis to obtain more stringent constraints at the tail-end of reionization.

The internal properties of dark matter haloes correlate with the large-scale halo clustering strength at fixed halo mass $-$ an effect known as assembly bias $-$ and are also strongly affected by the local, non-linear cosmic web. Characterising a halo's local web environment by its tidal anisotropy $\alpha$ at scales $\sim4$ x the halo radius, we demonstrate that these multi-scale correlations represent two distinct statistical links: one between the internal property and $\alpha$, and the other between $\alpha$ and large-scale ( $>30h^{-1}$Mpc) halo bias $b_1$. We focus on scalar internal properties of haloes related to formation time (concentration $c_{\rm vir}$), shape (mass ellipsoid asphericity $c/a$), velocity dispersion structure (velocity ellipsoid asphericity $c_v/a_v$ and velocity anisotropy $\beta$) and angular momentum (dimensionless spin $\lambda$) in the mass range $8\times10^{11}< M_{\rm vir}/(h^{-1}M_\odot)<5\times10^{14}$. Using conditional correlation coefficients and other detailed tests, we show that the joint distribution of $\alpha$, $b_1$ and any of the internal properties $c\in\{\beta,c_v/a_v,c/a,c_{\rm vir},\lambda\}$ is consistent with $p(\alpha,b_1,c)\simeq p(\alpha)p(b_1|\alpha)p(c|\alpha)$, at all but the largest masses. $\textit{Thus, the assembly bias trends $c-b_1$ reflect the two fundamental correlations $c-\alpha$ and $b_1-\alpha$.}$ Our results are unaffected by the exclusion of haloes with recent major merger events or splashback objects, although the latter are distinguished by the fact that $\alpha$ does not explain their assembly bias trends. The overarching importance of $\alpha$ provides a new perspective on the nature of assembly bias of distinct haloes, with potential ramifications for incorporating realistic assembly bias effects into mock catalogs of future large-scale structure surveys and for detecting galaxy assembly bias.

We present a Halo Occupation Distribution (HOD) analysis of the luminosity- and colour-dependent galaxy clustering in the Sloan Digital Sky Survey. A novelty of our technique is that it uses a combination of clustering measurements in luminosity bins to perform a global likelihood analysis, simultaneously constraining the HOD parameters for a range of luminosity thresholds. We present simple, smooth fitting functions which accurately describe the resulting luminosity dependence of the best-fit HOD parameters. To minimise systematic halo modelling effects, we use theoretical halo 2-point correlation functions directly measured and tabulated from a suite of $N$-body simulations spanning a large enough dynamic range in halo mass and spatial separation. Thus, our modelling correctly accounts for non-linear and scale-dependent halo bias as well as any departure of halo profiles from universality, and we additionally account for halo exclusion using the hard sphere approximation. Using colour-dependent clustering information, we constrain the satellite galaxy red fraction in a model-independent manner which does not rely on any group-finding algorithm. We find that the resulting luminosity dependence of the satellite red fraction is significantly shallower than corresponding measurements from galaxy group catalogues, and we provide a simple fitting function to describe this dependence. Our fitting functions are readily usable in generating low-redshift mock galaxy catalogues, and we discuss some potentially interesting applications as well as possible extensions of our technique.

Semi-numerical models of the reionization of neutral hydrogen (HI) based on the excursion set (ES) approach are known to violate photon number conservation at the few per cent level. In this work, we highlight a more severe, previously unrecognized shortcoming of ES models: the large-scale 21 cm power spectrum (equivalently, HI bias b_HI) is a relatively strong function of the spatial resolution used to generate ES ionization maps. We trace this problem to the fact that photon non-conservation in these models arises from a resolution-dependent mixture of spatially resolved, photon non-conserving bubbles, and partially ionized grid cells which are perfectly photon-conserving by construction. We argue that this inevitably leads to a resolution-dependence of b_HI, with the correct, converged value only emerging at very coarse resolution. Quantitatively, we find that b_HI can be non-converged by as much as ~20-25% in conservative ES implementations with grid sizes Delta x=5-10 cMpc/h, even when photon non-conservation is as small as ~3-4%. Thus, although numerically efficient, ES ionization maps coarse enough to produce a converged HI bias would wash out all topological features of the ionization field at scales k>~0.05 h/cMpc. We therefore present a new, explicitly photon conserving (PC) semi-numerical algorithm which distributes photons isotropically around sources while also accounting for anisotropic overlaps between nearby bubbles. Our PC algorithm predicts a resolution-independent value of b_HI consistent with the result of low-resolution ES maps, thus serving as a useful compromise between standard ES implementations and more expensive radiative transfer simulations.

Modelling the distribution of neutral hydrogen (HI) in dark matter halos is important for studying galaxy evolution in the cosmological context. We use a novel approach to infer the HI-dark matter connection at the massive end ($m_{\rm HI} > 10^{9.8} M_{\odot}$) from radio HI emission surveys, using optical properties of low-redshift galaxies as an intermediary. In particular, we use a previously calibrated optical HOD describing the luminosity- and colour-dependent clustering of SDSS galaxies and describe the HI content using a statistical scaling relation between the optical properties and HI mass. This allows us to compute the abundance and clustering properties of HI-selected galaxies and compare with data from the ALFALFA survey. We apply an MCMC-based statistical analysis to constrain the free parameters related to the scaling relation. The resulting best-fit scaling relation identifies massive HI galaxies primarily with optically faint blue centrals, consistent with expectations from galaxy formation models. We compare the HI-stellar mass relation predicted by our model with independent observations from matched HI-optical galaxy samples, finding reasonable agreement. As a further application, we make some preliminary forecasts for future observations of HI and optical galaxies in the expected overlap volume of SKA and Euclid/LSST.

We study the role of the local tidal environment in determining the assembly bias of dark matter haloes. Previous results suggest that the anisotropy of a halo's environment (i.e, whether it lies in a filament or in a more isotropic region) can play a significant role in determining the eventual mass and age of the halo. We statistically isolate this effect using correlations between the large-scale and small-scale environments of simulated haloes at $z=0$ with masses between $10^{11.6}\lesssim (m/h^{-1}M_{\odot})\lesssim10^{14.9}$. We probe the large-scale environment using a novel halo-by-halo estimator of linear bias. For the small-scale environment, we identify a variable $\alpha_R$ that captures the $\textit{tidal anisotropy}$ in a region of radius $R=4R_{\textrm{200b}}$ around the halo and correlates strongly with halo bias at fixed mass. Segregating haloes by $\alpha_R$ reveals two distinct populations. Haloes in highly isotropic local environments ($\alpha_R\lesssim0.2$) behave as expected from the simplest, spherically averaged analytical models of structure formation, showing a $\textit{negative}$ correlation between their concentration and large-scale bias at $\textit{all}$ masses. In contrast, haloes in anisotropic, filament-like environments ($\alpha_R\gtrsim0.5$) tend to show a $\textit{positive}$ correlation between bias and concentration at any mass. Our multi-scale analysis cleanly demonstrates how the overall assembly bias trend across halo mass emerges as an average over these different halo populations, and provides valuable insights towards building analytical models that correctly incorporate assembly bias. We also discuss potential implications for the nature and detectability of galaxy assembly bias.

The influence of the Cosmic Web on galaxy formation and evolution is of great observational and theoretical interest. We investigate whether the Cosmic Web leaves an imprint in the spatial clustering of galaxies in the Sloan Digital Sky Survey (SDSS), using the group catalog of Yang et al. and tidal field estimates at $\sim2h^{-1}$Mpc scales from the Mass-Tides-Velocity data set of Wang et al. We use the $\textit{tidal anisotropy}$ $\alpha$ (Paranjape et al.) to characterise the tidal environment of groups, and measure the redshift-space 2-point correlation function (2pcf) of group positions and the luminosity- and colour-dependent clustering of group galaxies using samples segregated by $\alpha$. We find that all the 2pcf measurements depend strongly on $\alpha$, with factors of $\sim20$ between the large-scale 2pcf of objects in the most and least isotropic environments. To test whether these strong trends imply `beyond halo mass' effects for galaxy evolution, we compare our results with corresponding 2pcf measurements in mock catalogs constructed using a halo occupation distribution that only uses halo mass as an input. We find that this prescription qualitatively reproduces $\textit{all}$ observed trends, and also quantitatively matches many of the observed results. Although there are some statistically significant differences between our `halo mass only' mocks and the data -- in the most and least isotropic environments -- which deserve further investigation, our results suggest that if the tidal environment induces additional effects on galaxy properties other than those inherited from their host halos, then these must be weak.

We present a fully analytical halo model of colour-dependent clustering that incorporates the effects of galactic conformity in a halo occupation distribution (HOD) framework. The model, based on our previous numerical work, describes conformity through a correlation between the colour of a galaxy and the concentration of its parent halo, leading to a correlation between central and satellite galaxy colours at fixed halo mass. The strength of the correlation is set by a tunable `group quenching efficiency', and the model can separately describe group-level correlations between galaxy colour (1-halo conformity) and large scale correlations induced by assembly bias (2-halo conformity). We validate our analytical results using clustering measurements in mock galaxy catalogs, finding that the model is accurate at the 10-20 percent level for a wide range of luminosities and length scales. We apply the formalism to interpret the colour-dependent clustering of galaxies in the Sloan Digital Sky Survey (SDSS). We find good overall agreement between the data and a model that has 1-halo conformity at a level consistent with previous results based on an SDSS group catalog, although the clustering data require satellites to be redder than suggested by the group catalog. Within our modelling uncertainties, however, we do not find strong evidence of 2-halo conformity driven by assembly bias in SDSS clustering.

We present a calibration of halo assembly bias using the Separate Universe technique. Specifically, we measure the response of halo abundances at fixed mass and concentration to the presence of an infinite-wavelength initial perturbation. We develop an analytical framework for describing the concentration dependence of this peak-background split halo bias -- a measure of assembly bias -- relying on the near-Lognormal distribution of halo concentration at fixed halo mass. The combination of this analytical framework and the Separate Universe technique allows us to achieve very high precision in the calibration of the linear assembly bias $b_1$, and qualitatively reproduces known trends such as the monotonic decrease (increase) of $b_1$ with halo concentration at large (small) masses. The same framework extends to the concentration dependence of higher order bias parameters $b_n$, and we present the first calibration of assembly bias in $b_2$. Our calibrations are directly applicable in analytical Halo Model calculations that seek to robustly detect galaxy assembly bias in observational samples. We detect a non-universality in the $b_1 - b_2$ relation arising from assembly bias, and suggest that simultaneous measurements of these bias parameters could be used to both detect the signature of assembly bias as well as mitigate its effects in cosmological analyses.

Cross-correlations between biased tracers and the dark matter field encode information about the physical variables which characterize these tracers. However, if the physical variables of interest are correlated with one another, then extracting this information is not as straightforward as one might naively have thought. We show how to exploit these correlations so as to estimate scale-independent bias factors of all orders in a model-independent way. We also show that failure to account for this will lead to incorrect conclusions about which variables matter and which do not. Morever, accounting for this allows one to use the scale dependence of bias to constrain the physics of halo formation; to date the argument has been phrased the other way around. We illustrate by showing that the scale dependence of linear and nonlinear bias, measured on nonlinear scales, can be used to provide consistent estimates of how the critical density for halo formation depends on halo mass. Our methods work even when the bias is nonlocal and stochastic, such as when, in addition to the spherically averaged density field and its derivatives, the quadrupolar shear field also matters for halo formation. In such models, the nonlocal bias factors are closely related to the more familiar local nonlinear bias factors, which are much easier to measure. Our analysis emphasizes the fact that biased tracers are biased because they do not sample fields (density, velocity, shear, etc.) at all positions in space in the same way that the dark matter does.

Recent analytical work on the modelling of dark halo abundances and clustering has demonstrated the advantages of combining the excursion set approach with peaks theory. We extend these ideas and introduce a model of excursion set peaks that incorporates the role of initial tidal effects or shear in determining the gravitational collapse of dark haloes. The model -- in which the critical density threshold for collapse depends on the tidal influences acting on protohaloes -- is well motivated from ellipsoidal collapse arguments and is also simple enough to be analytically tractable. We show that the predictions of this model are in very good agreement with measurements of the halo mass function and traditional scale dependent halo bias in N-body simulations across a wide range of masses and redshift. The presence of shear in the collapse threshold means that halo bias is naturally predicted to be nonlocal, and that protohalo densities at fixed mass are naturally predicted to have Lognormal-like distributions. We present the first direct estimate of Lagrangian nonlocal bias in N-body simulations, finding broad agreement with the model prediction. Finally, the simplicity of the model (which has essentially a single free parameter) opens the door to building efficient and accurate non-universal fitting functions of halo abundances and bias for use in precision cosmology.

Studying the cosmic dawn and the epoch of reionization through the redshifted 21 cm line are among the major science goals of the SKA1. Their significance lies in the fact that they are closely related to the very first stars in the universe. Interpreting the upcoming data would require detailed modelling of the relevant physical processes. In this article, we focus on the theoretical models of reionization that have been worked out by various groups working in India with the upcoming SKA in mind. These models include purely analytical and semi-numerical calculations as well as fully numerical radiative transfer simulations. The predictions of the 21 cm signal from these models would be useful in constraining the properties of the early galaxies using the SKA data.

We revisit the issue of interpreting the results of large volume cosmological simulations in the context of large scale general relativistic effects. We look for simple modifications to the nonlinear evolution of the gravitational potential $\psi$ that lead on large scales to the correct, fully relativistic description of density perturbations in the Newtonian gauge. We note that the relativistic constraint equation for $\psi$ can be cast as a diffusion equation, with a diffusion length scale determined by the expansion of the Universe. Exploiting the weak time evolution of $\psi$ in all regimes of interest, this equation can be further accurately approximated as a Helmholtz equation, with an effective relativistic 'screening' scale $\ell$ related to the Hubble radius. We demonstrate that it is thus possible to carry out N-body simulations in the Newtonian gauge by replacing Poisson's equation with this Helmholtz equation, involving a trivial change in the Green's function kernel. Our results also motivate a simple, approximate (but very accurate) gauge transformation - $\delta_{\rm N}(\mathbf{k}) \approx \delta_{\rm sim}(\mathbf{k})\times (k^2+\ell^{-2})/k^2$ - to convert the density field $\delta_{\rm sim}$ of standard collisionless N-body simulations (initialised in the comoving synchronous gauge) into the Newtonian gauge density $\delta_{\rm N}$ at arbitrary times. A similar conversion can also be written in terms of particle positions. Our results can be interpreted in terms of a Jeans stability criterion induced by the expansion of the Universe. The appearance of the screening scale $\ell$ in the evolution of $\psi$, in particular, leads to a natural resolution of the 'Jeans swindle' in the presence of super-horizon modes.

We use the Halo Model to explore the implications of assuming that galaxy luminosities in groups are randomly drawn from an underlying luminosity function. We show that even the simplest of such order statistics models -- one in which this luminosity function $p(L)$ is universal -- naturally produces a number of features associated with previous analyses based on the `central plus Poisson satellites' hypothesis. These include the monotonic relation of mean central luminosity with halo mass, the Lognormal distribution around this mean, and the tight relation between the central and satellite mass scales. In stark contrast to observations of galaxy clustering, however, this model predicts $\textit{no}$ luminosity dependence of large scale clustering. We then show that an extended version of this model, based on the order statistics of a $\textit{halo mass dependent}$ luminosity function $p(L|m)$, is in much better agreement with the clustering data as well as satellite luminosities, but systematically under-predicts central luminosities. This brings into focus the idea that central galaxies constitute a distinct population that is affected by different physical processes than are the satellites. We model this physical difference as a statistical brightening of the central luminosities, over and above the order statistics prediction. The magnitude gap between the brightest and second brightest group galaxy is predicted as a by-product, and is also in good agreement with observations. We propose that this order statistics framework provides a useful language in which to compare the Halo Model for galaxies with more physically motivated galaxy formation models

Traditional excursion set based models of H II bubble growth during the epoch of reionization are known to violate photon number conservation, in the sense that the mass fraction in ionized bubbles in these models does not equal the ratio of the number of ionizing photons produced by sources and the number of hydrogen atoms in the intergalactic medium. E.g., for a Planck13 cosmology with electron scattering optical depth $\tau\simeq0.066$, the discrepancy is $\sim15$ per cent for $x_{\rm HII}=0.1$ and $\sim5$ per cent for $x_{\rm HII}=0.5$. We demonstrate that this problem arises from a fundamental conceptual shortcoming of the excursion set approach (already recognised in the literature on this formalism) which only tracks average mass fractions instead of the exact, stochastic source counts. With this insight, we build an approximately photon number conserving Monte Carlo model of bubble growth based on partitioning regions of dark matter into halos. Our model, which is formally valid for white noise initial conditions (ICs), shows dramatic improvements in photon number conservation, as well as substantial differences in the bubble size distribution, as compared to traditional models. We explore the trends obtained on applying our algorithm to more realistic ICs, finding that these improvements are robust to changes in the ICs. Since currently popular semi-numerical schemes of bubble growth also violate photon number conservation, we argue that it will be worthwhile to pursue new, explicitly photon number conserving approaches. Along the way, we clarify some misconceptions regarding this problem that have appeared in the literature.

Intensity mapping of neutral hydrogen (HI) is a promising observational probe of cosmology and large-scale structure. We present wide field simulations of HI intensity maps based on N-body simulations of a $2.6\, {\rm Gpc / h}$ box with $2048^3$ particles (particle mass $1.6 \times 10^{11}\, {\rm M_\odot / h}$). Using a conditional mass function to populate the simulated dark matter density field with halos below the mass resolution of the simulation ($10^{8}\, {\rm M_\odot / h} < M_{\rm halo} < 10^{13}\, {\rm M_\odot / h}$), we assign HI to those halos according to a phenomenological halo to HI mass relation. The simulations span a redshift range of 0.35 < z < 0.9 in redshift bins of width $\Delta z \approx 0.05$ and cover a quarter of the sky at an angular resolution of about 7'. We use the simulated intensity maps to study the impact of non-linear effects and redshift space distortions on the angular clustering of HI. Focusing on the autocorrelations of the maps, we apply and compare several estimators for the angular power spectrum and its covariance. We verify that these estimators agree with analytic predictions on large scales and study the validity of approximations based on Gaussian random fields, particularly in the context of the covariance. We discuss how our results and the simulated maps can be useful for planning and interpreting future HI intensity mapping surveys.

To shed light on the fundamental problems posed by Dark Energy and Dark Matter, a large number of experiments have been performed and combined to constrain cosmological models. We propose a novel way of quantifying the information gained by updates on the parameter constraints from a series of experiments which can either complement earlier measurements or replace them. For this purpose, we use the Kullback-Leibler divergence or relative entropy from information theory to measure differences in the posterior distributions in model parameter space from a pair of experiments. We apply this formalism to a historical series of Cosmic Microwave Background experiments ranging from Boomerang to WMAP, SPT, and Planck. Considering different combinations of these experiments, we thus estimate the information gain in units of bits and distinguish contributions from the reduction of statistical errors and the `surprise' corresponding to a significant shift of the parameters' central values. For this experiment series, we find individual relative entropy gains ranging from about 1 to 30 bits. In some cases, e.g. when comparing WMAP and Planck results, we find that the gains are dominated by the surprise rather than by improvements in statistical precision. We discuss how this technique provides a useful tool for both quantifying the constraining power of data from cosmological probes and detecting the tensions between experiments.

We extend the Halo Occupation Distribution (HOD) framework to generate mock galaxy catalogs exhibiting varying levels of "galactic conformity", which has emerged as a potentially powerful probe of environmental effects in galaxy evolution. Our model correlates galaxy colours in a group with the concentration of the common parent dark halo through a "group quenching efficiency" $\rho$ which makes older, more concentrated halos $\textit{at fixed mass}$ preferentially host redder galaxies. We find that, for a specific value of $\rho$, this 1-halo conformity matches corresponding measurements in a group catalog based on the Sloan Digital Sky Survey. Our mocks also display conformity at large separations from isolated objects, potentially an imprint of halo assembly bias. A detailed study - using mocks with assembly bias erased while keeping 1-halo conformity intact - reveals a rather nuanced situation, however. At separations $\lesssim 4$Mpc, conformity is mainly a 1-halo effect dominated by the largest halos and is $\textit{not}$ a robust indicator of assembly bias. Only at very large separations ($\gtrsim 8$Mpc) does genuine 2-halo conformity, driven by the assembly bias of small halos, manifest distinctly. We explain all these trends in standard Halo Model terms. Our model opens the door to parametrized HOD analyses that self-consistently account for galactic conformity at all scales.

We compare analytical predictions of void volume functions to those measured from N-body simulations, detecting voids with the zobov void finder. We push to very small, nonlinear voids, below few Mpc radius, by considering the unsampled DM density field. We also study the case where voids are identified using halos. We develop analytical formula for the void abundance of both the excursion set approach and the peaks formalism. These formula are valid for random walks smoothed with a top-hat filter in real space, with a large class of realistic barrier models. We test the extent to which the spherical evolution approximation, which forms the basis of the analytical predictions, models the highly aspherical voids that occur in the cosmic web, and are found by a watershed-based algorithm such as zobov. We show that the volume function returned by zobov is quite sensitive to the choice of treatment of sub-voids, a fact that has not been appreciated previously. For reasonable choices of sub-void exclusion, we find that the Lagrangian density delta_v of the zobov voids -- which is predicted to be a constant delta_v = -2.7 in the spherical evolution model -- is different from the predicted value, showing substantial scatter and scale dependence. This result applies to voids identified at z=0 with effective radius between 1 and 10 Mpc/h. Our analytical approximations are flexible enough to give a good description of the resulting volume function; however, this happens for choices of parameter values that are different from those suggested by the spherical evolution assumption. We conclude that analytical models for voids must move away from the spherical approximation in order to be applied successfully to observations, and we discuss some possible ways forward.

Spectroscopic redshift surveys offer great prospects for constraining the dark sector in cosmology. Future surveys will however be both deep and wide and will thus require an analysis in 3-dimensional spherical geometry. We review and compare several methods which have been proposed in the literature for this purpose, focusing in particular on implementations of the spherical harmonic tomography (SHT) power spectrum $C^{i j}_{l}$ and the spherical Fourier Bessel (SFB) power spectrum $C_{l} (k, k')$. Using a Fisher analysis, we compare the forecasted constraints on cosmological parameters using these statistics. These constraints typically rely on approximations such as the Limber approximation and make specific choices in the numerical implementation of each statistic. Using a series of toy models, we explore the applicability of these approximations and study the sensitivity of the SHT and SFB statistics to the details of their implementation. In particular, we show that overlapping redshift bins may improve cosmological constraints using the SHT statistic when the number of bins is small, and that the SFB constraints are quite robust to changes in the assumed distance-redshift relation. We also find that the SHT can be tailored to be more sensitive to modes at redshifts close to the survey boundary, while the SFB appears better suited to capture information beyond the smooth shape of the power spectrum. In this context, we discuss the pros and cons of the different techniques and their impact on the design and analysis of future wide field spectroscopic surveys.

We construct a simple phenomenological model for the evolving galaxy population by incorporating pre-defined baryonic prescriptions into a dark matter hierarchical merger tree. Specifically the model is based on the simple gas-regulator model introduced by Lilly et al. 2013 coupled with the empirical quenching rules of Peng et al. 2010/12. The simplest model already does quite well in reproducing, without re-adjusting the input parameters, many observables including the Main Sequence sSFR-mass relation, the faint end slope of the galaxy mass function and the shape of the star-forming and passive mass functions. Compared with observations and/or the recent phenomenological model of Behroozi et al. 2013 based on epoch-dependent abundance-matching, our model also qualitatively reproduces the evolution of the Main Sequence sSFR(z) and SFRD(z) star formation rate density relations, the $M_s - M_h$ stellar-to-halo mass relation and also the $SFR - M_h$ relation. Quantitatively the evolution of sSFR(z) and SFRD(z) is not steep enough, the $M_s - M_h$ relation is not quite peaked enough and, surprisingly, the ratio of quenched to star-forming galaxies around M* is not quite high enough. We show that these deficiencies can simultaneously be solved by ad hoc allowing galaxies to re-ingest some of the gas previously expelled in winds, provided that this is done in a mass-dependent and epoch-dependent way. These allow the model galaxies to reduce an inherent tendency to saturate their star-formation efficiency. This emphasizes how efficient galaxies around M* are in converting baryons into stars and highlights the fact that quenching occurs just at the point when galaxies are rapidly approaching the maximum possible efficiency of converting baryons into stars.

We investigate potential systematic effects in constraining the amplitude of primordial fluctuations \sigma_8 arising from the choice of halo mass function in the likelihood analysis of current and upcoming galaxy cluster surveys. We study the widely used N-body simulation fit of Tinker et al. (T08) and, as an alternative, the recently proposed analytical model of Excursion Set Peaks (ESP). We first assess the relative bias between these prescriptions when constraining \sigma_8 by sampling the ESP mass function to generate mock catalogs and using the T08 fit to analyse them, for various choices of survey selection threshold, mass definition and statistical priors. To assess the level of absolute bias in each prescription, we then repeat the analysis on dark matter halo catalogs in N-body simulations designed to mimic the mass distribution in the current data release of Planck SZ clusters. This N-body analysis shows that using the T08 fit without accounting for the scatter introduced when converting between mass definitions (alternatively, the scatter induced by errors on the parameters of the fit) can systematically over-estimate the value of \sigma_8 by as much as 2\sigma\ for current data, while analyses that account for this scatter should be close to unbiased in \sigma_8. With an increased number of objects as expected in upcoming data releases, regardless of accounting for scatter, the T08 fit could over-estimate the value of \sigma_8 by ~1.5\sigma. The ESP mass function leads to systematically more biased but comparable results. A strength of the ESP model is its natural prediction of a weak non-universality in the mass function which closely tracks the one measured in simulations and described by the T08 fit. We suggest that it might now be prudent to build new unbiased ESP-based fitting functions for use with the larger datasets of the near future.

The size distribution of ionized regions during the epoch of reionization -- a key ingredient in understanding the HI power spectrum observable by 21cm experiments -- can be modelled analytically using the excursion set formalism of random walks in the smoothed initial density field. To date, such calculations have been based on simplifying assumptions carried forward from the earliest excursion set models of two decades ago. In particular, these models assume that the random walks have uncorrelated steps and that haloes can form at arbitrary locations in the initial density field. We extend these calculations by incorporating recent technical developments that allow us to (a) include the effect of correlations in the steps of the walks induced by a realistic smoothing filter and (b) more importantly, account for the fact that dark matter haloes preferentially form near peaks in the initial density. A comparison with previous calculations shows that including these features, particularly the peaks constraint on halo locations, has large effects on the size distribution of the HII bubbles surrounding these haloes. For example, when comparing models at the same value of the globally averaged ionized volume fraction, the typical bubble sizes predicted by our model are more than a factor 2 larger than earlier calculations. Our results can potentially have a significant impact on estimates of the observable HI power spectrum.

We investigate nonlocal Lagrangian bias contributions involving gradients of the linear density field, for which we have predictions from the excursion set peak formalism. We begin by writing down a bias expansion which includes all the bias terms, including the nonlocal ones. Having checked that the model furnishes a reasonable fit to the halo mass function, we develop a 1-point cross-correlation technique to measure bias factors associated with 2-distributed quantities. We validate the method with numerical realizations of peaks of Gaussian random fields before we apply it to N-body simulations. We focus on the lowest (quadratic) order nonlocal contributions. We can reproduce our measurement of \chi_{10} if we allow for an offset between the Lagrangian halo center-of-mass and the peak position. The sign and magnitude of \chi_{10} is consistent with Lagrangian haloes sitting near linear density maxima. The resulting contribution to the halo bias can safely be ignored for M = 10^13 Msun/h, but could become relevant at larger halo masses. For the second nonlocal bias \chi_{01} however, we measure a much larger magnitude than predicted by our model. We speculate that some of this discrepancy might originate from nonlocal Lagrangian contributions induced by nonspherical collapse.

We investigate the problem of predicting the halo mass function from the properties of the Lagrangian density field. We focus on a perturbation spectrum with a small-scale cut-off (as in warm dark matter cosmologies). This cut-off results in a strong suppression of low mass objects, providing additional leverage to rigorously test which perturbations collapse and to what mass. We find that all haloes are consistent with forming near peaks of the initial density field, with a strong correlation between proto-halo density and ellipticity. We demonstrate that, while standard excursion set theory with correlated steps completely fails to reproduce the mass function, the inclusion of the peaks constraint leads to the correct number of haloes but significantly underpredicts the masses of low-mass objects (with the predicted halo mass function at low masses behaving like dn/dln m ~ m^{2/3}). This prediction is very robust and cannot be easily altered within the framework of a single collapse barrier. The nature of collapse in the presence of a small-scale cut-off thus reveals that excursion set calculations require a more detailed understanding of the collapse-time of a general ellipsoidal perturbation to predict the ultimate collapsed mass of a peak -- a problem that has been hidden in the large abundance of small-scale structure in CDM. We demonstrate how this problem can be resolved within the excursion set framework.

We explore the scale dependence of halo bias using real space cross-correlation measurements in N-body simulations and in Pinocchio, an algorithm based on Lagrangian Perturbation Theory. Recent work has shown how to interpret such real space measurements in terms of k-dependent bias in Fourier space, and how to remove the k-dependence to reconstruct the k-independent peak-background split halo bias parameters. We compare our reconstruction of the linear bias, which requires no free parameters, with previous estimates from N-body simulations which were obtained directly in Fourier space at large scales, and find very good agreement. Our reconstruction of the quadratic bias is similarly parameter-free, although in this case there are no previous Fourier space measurements to compare with. Our analysis of N-body simulations explicitly tests the predictions of the excursion set peaks (ESP) formalism of Paranjape et al. (2013) for the scale dependence of bias; we find that the ESP predictions accurately describe our measurements. In addition, our measurements in Pinocchio serve as a useful, successful consistency check between Pinocchio and N-body simulations that is not accessible to traditional measurements.

We describe how to extend the excursion set peaks framework so that its predictions of dark halo abundances and clustering can be compared directly with simulations. These extensions include: a halo mass definition which uses the TopHat filter in real space; the mean dependence of the critical density for collapse delta_c on halo mass m; and the scatter around this mean value. All three of these are motivated by the physics of triaxial rather than spherical collapse. A comparison of the resulting mass function with N-body results shows that, if one uses delta_c(m) and its scatter as determined from simulations, then all three are necessary ingredients for obtaining ~10% accuracy. E.g., assuming a constant value of delta_c with no scatter, as motivated by the physics of spherical collapse, leads to many more massive halos than seen in simulations. The same model is also in excellent agreement with N-body results for the linear halo bias, especially at the high mass end where the traditional peak-background split argument applied to the mass function fit is known to underpredict the measured bias by ~10%. In the excursion set language, our model is about walks centered on special positions (peaks) in the initial conditions -- we discuss what it implies for the usual calculation in which all walks contribute to the statistics.

We describe a model of dark matter halo abundances and clustering which combines the two most widely used approaches to this problem: that based on peaks and the other based on excursion sets. Our approach can be thought of as addressing the cloud-in-cloud problem for peaks and/or modifying the excursion set approach so that it averages over a special subset, rather than all possible walks. In this respect, it seeks to account for correlations between steps in the walk as well as correlations between walks. We first show how the excursion set and peaks models can be written in the same formalism, and then use this correspondence to write our combined excursion set peaks model. We then give simple expressions for the mass function and bias, showing that even the linear halo bias factor is predicted to be k-dependent as a consequence of the nonlocality associated with the peak constraint. At large masses, our model has little or no need to rescale the variable delta_c from the value associated with spherical collapse, and suggests a simple explanation for why the linear halo bias factor appears to lie above that based on the peak-background split at high masses when such a rescaling is assumed. Although we have concentrated on peaks, our analysis is more generally applicable to other traditionally single-scale analyses of large-scale structure.

If one accounts for correlations between scales, then nonlocal, k-dependent halo bias is part and parcel of the excursion set approach, and hence of halo model predictions for galaxy bias. We present an analysis that distinguishes between a number of different effects, each one of which contributes to scale-dependent bias in real space. We show how to isolate these effects and remove the scale dependence, order by order, by cross-correlating the halo field with suitably transformed versions of the mass field. These transformations may be thought of as simple one-point, two-scale measurements that allow one to estimate quantities which are usually constrained using n-point statistics. As part of our analysis, we present a simple analytic approximation for the first crossing distribution of walks with correlated steps which are constrained to pass through a specified point, and demonstrate its accuracy. Although we concentrate on nonlinear, nonlocal bias with respect to a Gaussian random field, we show how to generalize our analysis to more general fields.

We show that the distribution of luminosities of Brightest Cluster Galaxies in an SDSS-based group catalog suggests that BCG luminosities are just the statistical extremes of the group galaxy luminosity function. This latter happens to be very well approximated by the all-galaxy luminosity function (restricted to Mr<-19.9), provided one uses a parametrization of this function that is accurate at the bright end. A similar analysis of the luminosity distribution of the Brightest Satellite Galaxies suggests that they are best thought of as being the second brightest pick from the same luminosity distribution of which BCGs are the brightest. I.e., BSGs are not the brightest of some universal satellite luminosity function, in contrast to what Halo Model analyses of the luminosity dependence of clustering suggest. However, we then use mark correlations to provide a novel test of these order statistics, showing that the hypothesis of a universal luminosity function (i.e. no halo mass dependence) from which the BCGs and BSGs are drawn is incompatible with the data, despite the fact that there was no hint of this in the BCG and BSG luminosity distributions themselves. We also discuss why, since extreme value statistics are explicitly a function of the number of draws, the consistency of BCG luminosities with extreme value statistics is most clearly seen if one is careful to perform the test at fixed group richness N. Tests at, e.g., fixed total group luminosity Ltot, will generally be biased and may lead to erroneous conclusions.

We extend earlier work on the problem of estimating the void-volume function -- the abundance and evolution of large voids which grow gravitationally in an expanding universe -- in two ways. The first removes an ambiguity about how the void-in-cloud process, which erases small voids, should be incorporated into the excursion set approach. The main technical change here is to think of voids within a fully Eulerian, rather than purely Lagrangian, framework. The second accounts for correlations between different spatial scales in the initial conditions. We provide numerical and analytical arguments showing how and why both changes modify the predicted abundances substantially. In particular, we show that the predicted importance of the void-in-cloud process depends strongly on whether or not one accounts for correlations between scales. With our new formulation, the void-in-cloud process dramatically reduces the predicted abundances of voids if such correlations are ignored, but only matters for the smallest voids in the more realistic case in which the spatial correlations are included.

In the Excursion Set approach, halo abundances and clustering are closely related. This relation is exploited in many modern methods which seek to constrain cosmological parameters on the basis of the observed spatial distribution of clusters. However, to obtain analytic expressions for these quantities, most Excursion Set based predictions ignore the fact that, although different k-modes in the initial Gaussian field are uncorrelated, this is not true in real space: the values of the density field at a given spatial position, when smoothed on different real-space scales, are correlated in a nontrivial way. We show that when the excursion set approach is extended to include such correlations, then one must be careful to account for the fact that the associated prediction for halo bias is explicitly a real-space quantity. Therefore, care must be taken when comparing the predictions of this approach with measurements in simulations, which are typically made in Fourier-space. We show how to correct for this effect, and demonstrate that ignorance of this effect in recent analyses of halo bias has led to incorrect conclusions and biased constraints.

The Excursion Set approach has been used to make predictions for a number of interesting quantities in studies of nonlinear hierarchical clustering. These include the halo mass function, halo merger rates, halo formation times and masses, halo clustering, analogous quantities for voids, and the distribution of dark matter counts in randomly placed cells. The approach assumes that all these quantities can be mapped to problems involving the first crossing distribution of a suitably chosen barrier by random walks. Most analytic expressions for these distributions ignore the fact that, although different k-modes in the initial Gaussian field are uncorrelated, this is not true in real space: the values of the density field at a given spatial position, when smoothed on different real-space scales, are correlated in a nontrivial way. As a result, the problem is to estimate first crossing distribution by random walks having correlated rather than uncorrelated steps. In 1990, Peacock & Heavens presented a simple approximation for the first crossing distribution of a single barrier of constant height by walks with correlated steps. We show that their approximation can be thought of as a correction to the distribution associated with what we call smooth completely correlated walks. We then use this insight to extend their approach to treat moving barriers, as well as walks that are constrained to pass through a certain point before crossing the barrier. For the latter, we show that a simple rescaling, inspired by bivariate Gaussian statistics, of the unconditional first crossing distribution, accurately describes the conditional distribution, independently of the choice of analytical prescription for the former. In all cases, comparison with Monte-Carlo solutions of the problem shows reasonably good agreement. (Abridged)

We study the effects of primordial non-Gaussianity on the large scale structure in the excursion set approach, accounting for correlations between steps of the random walks in the smoothed initial density field. These correlations are induced by realistic smoothing filters (as opposed to a filter that is sharp in k-space), but have been ignored by many analyses to date. We present analytical arguments -- building on existing results for Gaussian initial conditions -- which suggest that the effect of the filter at large smoothing scales is remarkably simple, and is in fact identical to what happens in the Gaussian case: the non-Gaussian walks behave as if they were smooth and deterministic, or "completely correlated". As a result, the first crossing distribution (which determines, e.g., halo abundances) follows from the single-scale statistics of the non-Gaussian density field -- the so-called "cloud-in-cloud" problem does not exist for completely correlated walks. Also, the answer from single-scale statistics is simply one half that for sharp-k walks. We explicitly test these arguments using Monte Carlo simulations of non-Gaussian walks, showing that the resulting first crossing distributions, and in particular the factor 1/2 argument, are remarkably insensitive to variations in the power spectrum and the defining non-Gaussian process. We also use our Monte Carlo walks to test some of the existing prescriptions for the non-Gaussian first crossing distribution. Since the factor 1/2 holds for both Gaussian and non-Gaussian initial conditions, it provides a theoretical motivation (the first, to our knowledge) for the common practice of analytically prescribing a ratio of non-Gaussian to Gaussian halo abundances.

Number counts of massive high-redshift clusters provide a window to study primordial non-Gaussianity. The current quality of data, however, forces the statistical analysis to probe a region of parameter space -- the extreme tail of the mass function -- which is neither accessible in any of the currently available theoretical prescriptions for calculating the mass function, nor calibrated in N-body simulations. In this work we present a new analytical prescription for calculating a "resummed" non-Gaussian halo mass function, which is constructed to remain stable in the extreme tail. We show that the prescription works well in the parameter regime that has been currently explored in simulations. We then use Fisher matrix techniques to compare our prescription with an extrapolated fit to N-body simulations, which has recently been used to obtain constraints from data collected by the South Pole Telecope. We show that for the current data, both prescriptions would lead to statistically consistent constraints. As the data improve, however, there is a possibility of introducing a statistically significant bias in the constraints due to the choice of prescription, especially if non-Gaussianity is scale dependent and becomes relatively large on cluster scales. It would then be necessary to test the accuracy of the prescriptions in N-body simulations that can probe clusters with high masses and redshifts in the presence of large non-Gaussianity.

The abundance of collapsed objects in the universe, or halo mass function, is an important theoretical tool in studying the effects of primordially generated non-Gaussianities on the large scale structure. The non-Gaussian mass function has been calculated by several authors in different ways, typically by exploiting the smallness of certain parameters which naturally appear in the calculation, to set up a perturbative expansion. We improve upon the existing results for the mass function by combining path integral methods and saddle point techniques (which have been separately applied in previous approaches). Additionally, we carefully account for the various scale dependent combinations of small parameters which appear. Some of these combinations in fact become of order unity for large mass scales and at high redshifts, and must therefore be treated non-perturbatively. Our approach allows us to do this, and to also account for multi-scale density correlations which appear in the calculation. We thus derive an accurate expression for the mass function which is based on approximations that are valid over a larger range of mass scales and redshifts than those of other authors. By tracking the terms ignored in the analysis, we estimate theoretical errors for our result and also for the results of others. We also discuss the complications introduced by the choice of smoothing filter function, which we take to be a top-hat in real space, and which leads to the dominant errors in our expression. Finally, we present a detailed comparison between the various expressions for the mass functions, exploring the accuracy and range of validity of each.

Primordial non-Gaussianity (NG) affects the large scale structure (LSS) of the universe by leaving an imprint on the distribution of matter at late times. Much attention has been focused on using the distribution of collapsed objects (i.e. dark matter halos and the galaxies and galaxy clusters that reside in them) to probe primordial NG. An equally interesting and complementary probe however is the abundance of extended underdense regions or voids in the LSS. The calculation of the abundance of voids using the excursion set formalism in the presence of primordial NG is subject to the same technical issues as the one for halos, which were discussed e.g. in arXiv:1005.1203. However, unlike the excursion set problem for halos which involved random walks in the presence of one barrier $\delta_c$, the void excursion set problem involves two barriers $\delta_v$ and $\delta_c$. This leads to a new complication introduced by what is called the "void-in-cloud" effect discussed in the literature, which is unique to the case of voids. We explore a path integral approach which allows us to carefully account for all these issues, leading to a rigorous derivation of the effects of primordial NG on void abundances. The void-in-cloud issue in particular makes the calculation conceptually rather different from the one for halos. However, we show that its final effect can be described by a simple yet accurate approximation. Our final void abundance function is valid on larger scales than the expressions of other authors, while being broadly in agreement with those expressions on smaller scales.

We explore in detail the semiclassical environment of collapsing shells of matter, and determine the semiclassical flux measured by a variety of observers. This study is a preliminary step in a broader investigation of thermodynamic properties of the geometry of collapsing objects. Specifically, in this paper we consider spherically symmetric null and timelike collapsing shells which form an event horizon, and calculate the flux measured by observers both inside and outside the shell, and both inside and outside the event horizon, and find nontrivial results in most of the cases. Additionally, we also investigate the environment of a shell which collapses but \emph{does not} form a horizon, halting at some radius larger than the Schwarzschild radius, and find that such an object generically gives rise to a pulse of radiation which is sharply peaked as it travels inwards and is reflected at the origin, and eventually emerges from the shell in a "thermalized" form. Our results have potential consequences in addressing questions pertaining, e.g. to black hole entropy and backreaction.

This thesis deals with the averaging problem in cosmology, which has gained considerable interest in recent years, and is concerned with correction terms (after averaging inhomogeneities) that appear in the Einstein equations when working on the large scales appropriate for cosmology. It has been claimed in the literature that these terms may account for the phenomenon of dark energy which causes the late time universe to accelerate. We investigate the nature of these terms by using averaging schemes available in the literature and further developed to be applicable to the problem at hand. We show that the effect of these terms when calculated carefully, remains negligible and cannot explain the late time acceleration.

We provide a detailed analysis of quantum field theory around a collapsing shell and discuss several conceptual issues related to the emission of radiation flux and formation of black holes. Explicit calculations are performed using a model for a collapsing shell which turns out to be analytically solvable. We use the insights gained in this model to draw reliable conclusions regarding more realistic models. We first show that any shell of mass $M$ which collapses to a radius close to $r=2M$ will emit approximately thermal radiation for a period of time. In particular, a shell which collapses from some initial radius to a final radius $2M(1-\epsilon^2)^{-1}$ (where $\epsilon \ll 1$) without forming a black hole, will emit thermal radiation during the period $M\lesssim t \lesssim M\ln (1/\epsilon^2)$. Later on ($t\gg M \ln(1/\epsilon^2)$), the flux from such a shell will decay to zero exponentially. We next study the effect of backreaction computed using the vacuum expectation value of the stress tensor on the collapse. We find that, in any realistic collapse scenario, the backreaction effects do \emph{not} prevent the formation of the event horizon. The time at which the event horizon is formed is, of course, delayed due to the radiated flux -- which decreases the mass of the shell -- but this effect is not sufficient to prevent horizon formation. We also clarify several conceptual issues and provide pedagogical details of the calculations in the Appendices to the paper.

There is an ongoing debate in the literature concerning the effects of averaging out inhomogeneities (``backreaction'') in cosmology. In particular, some simple models of structure formation studied in the literature seem to indicate that the backreaction can play a significant role at late times, and it has also been suggested that the standard perturbed FLRW framework is no longer a good approximation during structure formation, when the density contrast becomes nonlinear. In this work we use Zalaletdinov's covariant averaging scheme (macroscopic gravity or MG) to show that as long as the metric of the Universe can be described by the perturbed FLRW form, the corrections due to averaging remain negligibly small. Further, using a fully relativistic and reasonably generic model of pressureless spherical collapse, we show that as long as matter velocities remain small (which is true in our model), the perturbed FLRW form of the metric can be explicitly recovered. Together, these results imply that the backreaction remains small even during nonlinear structure formation, and we confirm this within the toy model with a numerical calculation.

If general relativity (GR) describes the expansion of the Universe, the observed cosmic acceleration implies the existence of a `dark energy'. However, while the Universe is on average homogeneous on large scales, it is inhomogeneous on smaller scales. While GR governs the dynamics of the \emph{in}homogeneous Universe, the averaged \emph{homogeneous} Universe obeys modified Einstein equations. Can such modifications alone explain the acceleration? For a simple generic model with realistic initial conditions, we show the answer to be `no'. Averaging effects negligibly influence the cosmological dynamics.

The problem of corrections to Einstein's equations arising from averaging of inhomogeneities ("backreaction") in the cosmological context, has gained considerable attention recently. We present results of analysing cosmological perturbation theory in the framework of Zalaletdinov's fully covariant Macroscopic Gravity. We show that this framework can be adapted to the setting of cosmological perturbations in a manner which is free from gauge related ambiguities. We derive expressions for the backreaction which can be readily applied in \emph{any} situation (not necessarily restricted to the linear perturbations considered here) where the \emph{metric} can be brought to the perturbed FLRW form. In particular these expressions can be employed in toy models studying nonlinear structure formation, and possibly also in N-body simulations. Additionally, we present results of example calculations which show that the backreaction remains negligible well into the matter dominated era.

There is an ongoing debate in the literature as to whether the effects of averaging out inhomogeneities (``backreaction'') in Cosmology can be large enough to account for the acceleration of the scale factor in the FLRW models. In particular, some simple models of structure formation studied in the literature seem to indicate that this is indeed possible, and it has also been suggested that the perturbed FLRW framework is no longer a good approximation during structure formation, when the density contrast becomes nonlinear. In this work we attempt to clarify the situation to some extent, using a fully relativistic model of pressureless spherical collapse. We find that whereas averaging during structure formation can lead to acceleration via a selective choice of averaging domains, the acceleration is not present when more generic domains are used for averaging. Further, we show that for most of the duration of the collapse, matter velocities remain small, and the perturbed FLRW form of the metric can be explicitly recovered, in the structure formation phase. We also discuss the fact that the magnitude of the average effects of inhomogeneities depends on the scale of averaging, and while it may not be completely negligible on intermediate scales, it is expected to remain small when averaging on suitably large scales.

The present matter density of the Universe, while highly inhomogeneous on small scales, displays approximate homogeneity on large scales. We propose that whereas it is justified to use the Friedmann-Lemaitre-Robertson-Walker (FLRW) line element (which describes an exactly homogeneous and isotropic universe) as a template to construct luminosity distances in order to compare observations with theory, the evolution of the scale factor in such a construction must be governed not by the standard Einstein equations for the FLRW metric, but by the modified Friedmann equations derived by Buchert [7, 8] in the context of spatial averaging in Cosmology. Furthermore, we argue that this scale factor, defined in the spatially averaged cosmology, will correspond to the effective FLRW metric provided the size of the averaging domain coincides with the scale at which cosmological homogeneity arises. This allows us, in principle, to compare predictions of a spatially averaged cosmology with observations, in the standard manner, for instance by computing the luminosity distance versus red-shift relation. The predictions of the spatially averaged cosmology would in general differ from standard FLRW cosmology, because the scale-factor now obeys the modified FLRW equations. This could help determine, by comparing with observations, whether or not cosmological inhomogeneities are an alternative explanation for the observed cosmic acceleration.

It is known that any explicit averaging scheme of the type essential for describing the large scale behaviour of the Universe, must necessarily yield corrections to the Einstein equations applied in the Cosmological setting. The question of whether or not the resulting corrections to the Einstein equations are significant, is still a subject of debate, partly due to possible ambiguities in the averaging schemes available. In particular, it has been argued in the literature that the effects of averaging could be gauge artifacts. We apply the formalism of Zalaletdinov's Macroscopic Gravity (MG) which is a fully covariant and nonperturbative averaging scheme, in an attempt to construct gauge independent corrections to the standard Friedmann-Lemaitre-Robertson-Walker (FLRW) equations. We find that whereas one cannot escape the problem of dependence on \emph{one} gauge choice -- which is inherent in the assumption of large scale homogeneity and isotropy -- it is however possible to construct \emph{spacetime scalar} corrections to the standard FLRW equations. This partially addresses the criticism concerning the corrections being gauge artifacts. For a particular initial choice of gauge which simplifies the formalism, we explicitly construct these scalars in terms of the underlying inhomogeneous geometry, and incidentally demonstrate that the formal structure of the corrections with this gauge choice is identical to that of analogous corrections derived by Buchert in the context of spatial averaging of scalars.

A consistent approach to Cosmology requires an explicit averaging of the Einstein equations, to describe a homogeneous and isotropic geometry. Such an averaging will in general modify the Einstein equations. The averaging procedure due to Buchert has attracted considerable attention recently since it offers the tantalizing hope of explaining the phenomenon of dark energy through such corrections. This approach has been criticized, however, on the grounds that its effects may be gauge artifacts. We apply the fully covariant formalism of Zalaletdinov's Macroscopic Gravity and show that, after making some essential gauge choices, the Cosmological equations receive \emph{spacetime scalar} corrections which are therefore observable in principle, and further, that the broad structure of these corrections is \emph{identical} to those derived by Buchert.

The null surfaces of a spacetime act as one-way membranes and can block information for a corresponding family of observers (time-like curves). Since lack of information can be related to entropy, this suggests the possibility of assigning an entropy to the null surfaces of a spacetime. We motivate and introduce such an entropy functional for any vector field in terms of a fourth-rank divergence free tensor P_{ab}^{cd} with the symmetries of the curvature tensor. Extremising this entropy for all the null surfaces then leads to equations for the background metric of the spacetime. When P_{ab}^{cd} is constructed from the metric alone, these equations are identical to Einstein's equations with an undetermined cosmological constant (which arises as an integration constant). More generally, if P_{ab}^{cd} is allowed to depend on both metric and curvature in a polynomial form, one recovers the Lanczos-Lovelock gravity. In all these cases: (a) We only need to extremise the entropy associated with the null surfaces; the metric is not a dynamical variable in this approach. (b) The extremal value of the entropy agrees with standard results, when evaluated on-shell for a solution admitting a horizon. The role of full quantum theory of gravity will be to provide the specific form of P_{ab}^{cd} which should be used in the entropy functional. With such an interpretation, it seems reasonable to interpret the Lanczos-Lovelock type terms as quantum corrections to classical gravity.

We investigate the possible occurrence of a positive cosmic acceleration in a spatially averaged, expanding, unbound Lemaitre-Tolman-Bondi cosmology. By studying an approximation in which the contribution of three-curvature dominates over the matter density, we construct numerical models which exhibit acceleration.

Spacetimes with horizons show a resemblance to thermodynamic systems and one can associate the notions of temperature and entropy with them. In the case of Einstein-Hilbert gravity, it is possible to interpret Einstein's equations as the thermodynamic identity TdS = dE + PdV for a spherically symmetric spacetime and thus provide a thermodynamic route to understand the dynamics of gravity. We study this approach further and show that the field equations for Lanczos-Lovelock action in a spherically symmetric spacetime can also be expressed as TdS = dE + PdV with S and E being given by expressions previously derived in the literature by other approaches. The Lanczos-Lovelock Lagrangians are of the form L=Q_a^{bcd}R^a_{bcd} with \nabla_b Q^{abcd}=0. In such models, the expansion of Q^{abcd} in terms of the derivatives of the metric tensor determines the structure of the theory and higher order terms can be interpreted quantum corrections to Einstein gravity. Our result indicates a deep connection between the thermodynamics of horizons and the allowed quantum corrections to standard Einstein gravity, and shows that the relation TdS = dE + PdV has a greater domain of validity that Einstein's field equations.