We present a comprehensive discussion on lattice techniques for the simulation of scalar and gauge field dynamics in an expanding universe. After reviewing the continuum formulation of scalar and gauge field interactions in Minkowski and FLRW backgrounds, we introduce basic tools for the discretization of field theories, including lattice gauge invariant techniques. Following, we discuss and classify numerical algorithms, ranging from methods of $O(dt^2)$ accuracy like $staggered~leapfrog$ and $Verlet$ integration, to $Runge-Kutta$ methods up to $O(dt^4)$ accuracy, and the $Yoshida$ and $Gauss-Legendre$ higher-order integrators, accurate up to $O(dt^{10})$. We adapt these methods for their use in classical lattice simulations of the non-linear dynamics of scalar and gauge fields in an expanding grid in $3+1$ dimensions, including the case of `self-consistent' expansion sourced by the volume average of the fields' energy and pressure densities. We present lattice formulations of canonical cases of: $i)$ Interacting scalar fields, $ii)$ Abelian $U(1)$ gauge theories, and $iii)$ Non-Abelian $SU(2)$ gauge theories. In all three cases we provide symplectic integrators, with accuracy ranging from $O(dt^2)$ up to $O(dt^{10})$. For each algorithm we provide the form of relevant observables, such as energy density components, field spectra and the Hubble constraint. We note that all our algorithms for gauge theories always respect the Gauss constraint to machine precision, even in the case of 'self-consistent' expansion. As a numerical example we analyze the post-inflationary dynamics of an oscillating inflaton charged under $SU(2)\times U(1)$. We note that the present manuscript is meant as part of the theoretical basis for the code $CosmoLattice$, a multi-purpose MPI-based package for simulating the non-linear evolution of field theories in an expanding universe, to be released in 2020.

Observations of the redshifted 21-cm signal (in absorption or emission) allow us to peek into the epoch of "dark ages" and the onset of reionization. These data can provide a novel way to learn about the nature of dark matter, in particular about the formation of small size dark matter halos. However, the connection between the formation of structures and 21-cm signal requires knowledge of stellar to total mass relation, escape fraction of UV photons, and other parameters that describe star formation and radiation at early times. This baryonic physics depends on the properties of dark matter and in particular in warm-dark-matter (WDM) models, star formation may follow a completely different scenario, as compared to the cold-dark-matter case. We use the recent measurements by the EDGES [J. D. Bowman, A. E. E. Rogers, R. A. Monsalve, T. J. Mozdzen, and N. Mahesh, An absorption profile centred at 78 megahertz in thesky-averaged spectrum,Nature (London) 555, 67 (2018).] to demonstrate that when taking the above considerations into account, the robust WDM bounds are in fact weaker than those given by the Lyman-$\alpha$ forest method and other structure formation bounds. In particular, we show that resonantly produced 7 keV sterile neutrino dark matter model is consistent with these data. However, a holistic approach to modelling of the WDM universe holds great potential and may in the future make 21-cm data our main tool to learn about dark matter clustering properties.

Flux ratio anomalies in strong gravitationally lensed quasars constitute a unique way to probe the abundance of non-luminous dark matter haloes, and hence the nature of dark matter. In this paper we identify double imaged quasars as a statistically efficient probe of dark matter, since they are 20 times more abundant than quadruply imaged quasars. Using N-body simulations that include realistic baryonic feedback, we measure the full distribution of flux ratios in doubly imaged quasars for cold (CDM) and warm dark matter (WDM) cosmologies. Through this method, we fold in two key systematics - quasar variability and line-of-sight structures. We find that WDM cosmologies predict a ~6 per cent difference in the cumulative distribution functions of flux ratios relative to CDM, with CDM predicting many more small ratios. Finally, we estimate that ~600 doubly imaged quasars will need to be observed in order to be able to unambiguously discern between CDM and the two WDM models studied here. Such sample sizes will be easily within reach of future large scale surveys such as Euclid. In preparation for this survey data we require discerning the scale of the uncertainties in modelling lens galaxies and their substructure in simulations, plus a strong understanding of the selection function of observed lensed quasars.

We quantify the error in the results of mixed baryon--dark-matter hydrodynamic simulations, stemming from outdated approximations for the generation of initial conditions. The error at redshift 0 in contemporary large simulations, is of the order of few to ten percent in the power spectra of baryons and dark matter, and their combined total-matter power spectrum. After describing how to properly assign initial displacements and peculiar velocities to multiple species, we review several approximations: (1) {using the total-matter power spectrum to compute displacements and peculiar velocities of both fluids}, (2) scaling the linear redshift-zero power spectrum back to the initial power spectrum using the Newtonian growth factor ignoring homogeneous radiation, (3) using longitudinal-gauge velocities with synchronous-gauge densities, and (4) ignoring the phase-difference in the Fourier modes for the offset baryon grid, relative to the dark-matter grid. Three of these approximations do not take into account that dark matter and baryons experience a scale-dependent growth after photon decoupling, which results in directions of velocity which are not the same as their direction of displacement. We compare the outcome of hydrodynamic simulations with these four approximations to our reference simulation, all setup with the same random seed and simulated using Gadget-III.

We measure the weak gravitational lensing shear power spectra and their cross-power in two photometric redshift bins from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). The measurements are performed directly in multipole space in terms of adjustable band powers. For the extraction of the band powers from the data we have implemented and extended a quadratic estimator, a maximum likelihood method that allows us to readily take into account irregular survey geometries, masks, and varying sampling densities. We find the 68 per cent credible intervals in the $\sigma_8$-$\Omega_{\rm m}$-plane to be marginally consistent with results from $Planck$ for a simple five parameter $\Lambda$CDM model. For the projected parameter $S_8 \equiv \sigma_8(\Omega_{\rm m}/0.3)^{0.5}$ we obtain a best-fitting value of $S_8 = 0.768_{-0.039}^{+0.045}$. This constraint is consistent with results from other CFHTLenS studies as well as the Dark Energy Survey. Our most conservative model, including modifications to the power spectrum due to baryon feedback and marginalization over photometric redshift errors, yields an upper limit on the total mass of three degenerate massive neutrinos of $\Sigma m_\nu < 4.53 \, {\rm eV}$ at 95 per cent credibility, while a Bayesian model comparison does not favour any model extension beyond a simple five parameter $\Lambda$CDM model. Combining the shear likelihood with $Planck$ breaks the $\sigma_8$-$\Omega_{\rm m}$-degeneracy and yields $\sigma_8=0.818 \pm 0.013$ and $\Omega_{\rm m} = 0.300 \pm 0.011$ which is fully consistent with results from $Planck$ alone.

We present a description for setting initial particle displacements and field values for simulations of arbitrary metric theories of gravity, for perfect and imperfect fluids with arbitrary characteristics. We extend the Zel'dovich Approximation to nontrivial theories of gravity, and show how scale dependence implies curved particle paths, even in the entirely linear regime of perturbations. For a viable choice of Effective Field Theory of Modified Gravity, initial conditions set at high redshifts are affected at the level of up to 5% at Mpc scales, which exemplifies the importance of going beyond {\Lambda}-Cold Dark Matter initial conditions for modifications of gravity outside of the quasi-static approximation. In addition, we show initial conditions for a simulation where a scalar modification of gravity is modelled in a Lagrangian particle-like description. Our description paves the way for simulations and mock galaxy catalogs under theories of gravity beyond the standard model, crucial for progress towards precision tests of gravity and cosmology.

We examine the deviation of Cold Dark Matter particle trajectories from the Newtonian result as the size of the region under study becomes comparable to or exceeds the particle horizon. To first order in the gravitational potential, the general relativistic result coincides with the Zel'dovich approximation and hence the Newtonian prediction on all scales. At second order, General Relativity predicts corrections which overtake the corresponding second order Newtonian terms above a certain scale of the order of the Hubble radius. However, since second order corrections are very much suppressed on such scales, we conclude that simulations which exceed the particle horizon but use Newtonian equations to evolve the particles, reproduce the correct trajectories very well. The dominant relativistic corrections to the power spectrum on scales close to the horizon are at most of the order of $\sim 10^{-5}$ at $z=49$ and $\sim 10^{-3}$ at $z=0$. The differences in the positions of real space features are affected at a level below $10^{-6}$ at both redshifts. Our analysis also clarifies the relation of N-body results to relativistic considerations.

We discuss the relation between the output of Newtonian N-body simulations on scales that approach or exceed the particle horizon to the description of General Relativity. At leading order, the Zeldovich approximation is correct on large scales, coinciding with the General Relativistic result. At second order in the initial metric potential, the trajectories of particles deviate from the second order Newtonian result and hence the validity of 2LPT initial conditions should be reassessed when used in very large simulations. We also advocate using the expression for the synchronous gauge density as a well behaved measure of density fluctuations on such scales.

We argue that there is an intrinsic noise on measurements of the equation of state parameter $w=p/\rho$ from large-scale structure around us. The presence of the large-scale structure leads to an ambiguity in the definition of the background universe and thus there is a maximal precision with which we can determine the equation of state of dark energy. To study the uncertainty due to local structure, we model density perturbations stemming from a standard inflationary power spectrum by means of the exact Lema\^{i}tre-Tolman-Bondi solution of Einstein's equation, and show that the usual distribution of matter inhomogeneities in a $\Lambda$CDM cosmology causes a variation of $w$ -- as inferred from distance measures -- of several percent. As we observe only one universe, or equivalently because of the cosmic variance, this uncertainty is systematic in nature.

We present a new analysis on how to distinguish between isotropic and anisotropic cosmological models based on tracking the angular displacements of a large number of distant quasars over an extended period of time, and then performing a multipole-vector decomposition of the resulting displacement maps. We find that while the GAIA mission operating at its nominal specifications does not have sufficient angular resolution to resolve anisotropic universes from isotropic ones using this method within a reasonable timespan of ten years, a next-generation GAIA-like survey with a resolution ten times better should be equal to the task. Distinguishing between different anisotropic models is however more demanding. Keeping the observational timespan to ten years, we find that the angular resolution of the survey will need to be of order 0.1 micro-arcsecond in order for certain rotating anisotropic models to produce a detectable signature that is also unique to models of this class. However, should such a detection become possible, it would immediately allow us to rule out large local void models.

We present a new programme for placing constraints on radial inhomogeneity in a dark-energy dominated universe. We introduce a new measure to quantify violations of the Copernican principle. Any violation of this principle would interfere with our interpretation of any dark-energy evolution. In particular, we find that current observations place reasonably tight constraints on possible late-time violations of the Copernican principle: the allowed area in the parameter space of amplitude and scale of a spherical inhomogeneity around the observer has to be reduced by a factor of three so as to confirm the Copernican principle. Then, by marginalizing over possible radial inhomogeneity we provide the first constraints on the cosmological constant which are free of the homogeneity prior prevalent in cosmology.

A perfect irrotational fluid with the equation of state of dust, Irrotational Dark Matter (IDM), is incapable of virializing and instead forms a cosmoskeleton of filaments with supermassive black holes at the joints. This stark difference from the standard cold dark matter (CDM) scenario arises because IDM must exhibit potential flow at all times, preventing shell-crossing from occurring. This scenario is applicable to general non-oscillating scalar-field theories with a small sound speed. Our model of combined IDM and CDM components thereby provides a solution to the problem of forming the observed billion-solar-mass black holes at redshifts of six and higher. In particular, as a result of the reduced vortical flow, the growth of the black holes is expected to be more rapid at later times as compared to the standard scenario.

There is an approximately 9% discrepancy, corresponding to 2.4sigma, between two independent constraints on the expansion rate of the universe: one indirectly arising from the cosmic microwave background and baryon acoustic oscillations, and one more directly obtained from local measurements of the relation between redshifts and distances to sources. We argue that by taking into account the local gravitational potential at the position of the observer this tension - strengthened by the recent Planck results - is partially relieved and the concordance of the standard model of cosmology increased. We estimate that measurements of the local Hubble constant are subject to a cosmic variance of about 2.4% (limiting the local sample to redshifts z>0.010) or 1.3% (limiting it to z>0.023), a more significant correction than that taken into account already. Nonetheless, we show that one would need a very rare fluctuation to fully explain the offset in the Hubble rates. If this tension is further strengthened, a cosmology beyond the standard model may prove necessary.

We find that if we live at the center of an inhomogeneity with total density contrast of roughly 0.1, dark energy is not a cosmological constant at 95% confidence level. Observational constraints on the equation of state of dark energy, w, depend strongly on the local matter density around the observer. We model the local inhomogeneity with an exact spherically symmetric solution which features a pressureless matter component and a dark-energy fluid with constant equation of state and negligible sound speed, that reaches a homogeneous solution at finite radius. We fit this model to observations of the local expansion rate, distant supernovae and the cosmic microwave background. We conclude that the possible uncertainty from large-scale structure has to be taken into account if one wants to progress towards not just precision but also accurate cosmology.

We describe inhomogeneities in a {\Lambda}CDM universe with a gradient series expansion and show that it describes the gravitational evolution far into the non-linear regime and beyond the capacity of standard perturbation theory at any order. We compare the gradient expansion with exact inhomogeneous {\Lambda}LTB solutions (Lema\^itre-Tolman-Bondi metric with the inclusion of a cosmological constant) describing growing structure in a {\Lambda}CDM universe and find that the expansion approximates the exact solution well, following the collapse of an over-density all the way into a singularity.

We present semi-analytical solutions to the background equations describing the Lema\^itre-Tolman-Bondi (LTB) metric as well as the homogeneous Friedmann equations, in the presence of dust, curvature and a cosmological constant Lambda. For none of the presented solutions any numerical integration has to be performed. All presented solutions are given for expanding and collapsing phases, preserving continuity in time and radius. Hence, these solutions describe the complete space time of a collapsing spherical object in an expanding universe. In the appendix we present for completeness a solution of the Friedmann equations in the additional presence of radiation, only valid for the Robertson-Walker metric.

If we live in the vicinity of the hypothesized Great Attractor, the age of the universe as inferred from the local expansion rate can be off by three per cent. We study the effect that living inside or near a massive overdensity has on cosmological parameters induced from observations of supernovae, the Hubble parameter and the Cosmic Microwave Background. We compare the results to those for an observer in a perfectly homogeneous LCDM universe. We find that for instance the inferred value for the global Hubble parameter changes by around three per cent if we happen to live inside a massive overdensity such as the hypothesized Great Attractor. Taking into account the effect of such structures on our perception of the universe makes cosmology perhaps less precise, but more accurate.

The Cold Spot could be an adiabatic perturbation on the surface of last scattering, in which case it is an over-density with comoving radius of the order of 1 Gpc. We assess the effect that living in a similar structure, without knowing it, has on our perception of the equation of state of Dark Energy. We find that structures of dimensions such that they could cause the Cold Spot on the CMB, affect the perceived equation of state of Dark Energy possibly up to ten percent.

In this paper, instead of invoking Dark Energy, we try and fit various cosmological observations with a large Gpc scale under-dense region (Void) which is modeled by a Lemaitre-Tolman-Bondi metric that at large distances becomes a homogeneous FLRW metric. We improve on previous analyses by allowing for nonzero overall curvature, accurately computing the distance to the last-scattering surface and the observed scale of the Baryon Acoustic peaks, and investigating important effects that could arise from having nontrivial Void density profiles. We mainly focus on the WMAP 7-yr data (TT and TE), Supernova data (SDSS SN), Hubble constant measurements (HST) and Baryon Acoustic Oscillation data (SDSS and LRG). We find that the inclusion of a nonzero overall curvature drastically improves the goodness of fit of the Void model, bringing it very close to that of a homogeneous universe containing Dark Energy, while by varying the profile one can increase the value of the local Hubble parameter which has been a challenge for these models. We also try to gauge how well our model can fit the large-scale-structure data, but a comprehensive analysis will require the knowledge of perturbations on LTB metrics. The model is consistent with the CMB dipole if the observer is about 15 Mpc off the centre of the Void. Remarkably, such an off-center position may be able to account for the recent anomalous measurements of a large bulk flow from kSZ data. Finally we provide several analytical approximations in different regimes for the LTB metric, and a numerical module for CosmoMC, thus allowing for a MCMC exploration of the full parameter space.

It has been argued that the Swiss-Cheese cosmology can mimic Dark Energy, when it comes to the observed luminosity distance-redshift relation. Besides the fact that this effect tends to disappear on average over random directions, we show in this work that based on the Rees-Sciama effect on the cosmic microwave background (CMB), the Swiss-Cheese model can be ruled out if all holes have a radius larger than about 35 Mpc. We also show that for smaller holes, the CMB is not observably affected, and that the small holes can still mimic Dark Energy, albeit in special directions, as opposed to previous conclusions in the literature. However, in this limit, the probability of looking in a special direction where the luminosity of supernovae is sufficiently supressed becomes very small, at least in the case of a lattice of spherical holes considered in this paper.

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

There has been much recent discussion, and some confusion, regarding the use of existing observational data to estimate the likelihood that next-generation cosmic microwave background (CMB) polarization experiments might detect a nonzero tensor signal, possibly associated with inflation. We examine this issue in detail here in two different ways: (1) first we explore the effect of choice of different parameter priors on the estimation of the tensor-to-scalar ratio r and other parameters describing inflation, and (2) we examine the Bayesian complexity in order to determine how effectively existing data can constrain inflationary parameters. We demonstrate that existing data are not strong enough to render full inflationary parameter estimates in a parametrization- and prior-independent way and that the predicted tensor signal is particularly sensitive to different priors. For parametrizations where the Bayesian complexity is comparable to the number of free parameters we find that a flat prior on the scale of inflation (which is to be distinguished from a flat prior on the tensor-to-scalar ratio) leads us to infer a larger, and in fact slightly nonzero tensor contribution at 68% confidence level. However, no detection is claimed. Our results demonstrate that all that is statistically relevant at the current time is the (slightly enhanced) upper bound on r, and we stress that the data remain consistent with r = 0.

Temperature anisotropies in the Cosmic Microwave Background (CMB) are affected by the late Integrated Sachs-Wolfe (lISW) effect caused by any time-variation of the gravitational potential on linear scales. Dark energy is not the only source of lISW, since massive neutrinos induce a small decay of the potential on small scales during both matter and dark energy domination. In this work, we study the prospect of using the cross-correlation between CMB and galaxy density maps as a tool for constraining the neutrino mass. On the one hand massive neutrinos reduce the cross-correlation spectrum because free-streaming slows down structure formation; on the other hand, they enhance it through their change in the effective linear growth. We show that in the observable range of scales and redshifts, the first effect dominates, but the second one is not negligible. We carry out an error forecast analysis by fitting some mock data inspired by the Planck satellite, Dark Energy Survey (DES) and Large Synoptic Survey Telescope (LSST). The inclusion of the cross-correlation data from Planck and LSST increases the sensitivity to the neutrino mass m_nu by 38% (and to the dark energy equation of state w by 83%) with respect to Planck alone. The correlation between Planck and DES brings a far less significant improvement. This method is not potentially as good for detecting m_nu as the measurement of galaxy, cluster or cosmic shear power spectra, but since it is independent and affected by different systematics, it remains potentially interesting if the total neutrino mass is of the order of 0.2 eV; if instead it is close to the lower bound from atmospheric oscillations, m_nu ~ 0.05 eV, we do not expect the ISW-galaxy correlation to be ever sensitive to m_nu.

We update constraints on the Hubble function H(phi) during inflation, using the most recent cosmic microwave background (CMB) and large scale structure (LSS) data. Our main focus is on a comparison between various commonly used methods of calculating the primordial power spectrum via analytical approximations and the results obtained by integrating the exact equations numerically. In each case, we impose naive, minimally restrictive priors on the duration of inflation. We find that the choice of priors has an impact on the results: the bounds on inflationary parameters can vary by up to a factor two. Nevertheless, it should be noted that within the region allowed by the minimal prior of the exact method, the accuracy of the approximations is sufficient for current data. We caution however that a careless minimal implementation of the approximative methods allows models for which the assumptions behind the analytical approximations fail, and recommend using the exact numerical method for a self-consistent analysis of cosmological data.

We derive new constraints on the Hubble function H(phi) and subsequently on the inflationary potential V(phi) from WMAP 3-year data combined with the Sloan Luminous Red Galaxy survey (SDSS-LRG), using a new methodology which appears to be more generic, conservative and model-independent than in most of the recent literature, since it depends neither on the slow-roll approximation, nor on any extrapolation scheme for the potential beyond the observable e-fold range, nor on additional assumptions about initial conditions for the inflaton velocity. This last feature represents the main improvement of this work, and is made possible by the reconstruction of H(phi) prior to V(phi). Our results only rely on the assumption that within the observable range, corresponding to ~ 10 e-folds, inflation is not interrupted and the function H(phi) is smooth enough for being Taylor-expanded at order one, two or three. We conclude that the variety of potentials allowed by the data is still large. However, it is clear that the first two slow-roll parameters are really small while the validity of the slow-roll expansion beyond them is not established.

We discuss the latest constraints on a Taylor-expanded scalar inflaton potential, obtained focusing on its observable part only. This is in contrast with other works in which an extrapolation of the potential is applied using the slow-roll hierarchy. We find significant differences. The results discussed here apply to a broader range of models, since no assumption about the invisible e-folds of inflation has to be made, thereby remaining conservative.

We derive some new constraints on single-field inflation from the Wilkinson Microwave Anisotropy Probe 3-year data combined with the Sloan Luminous Red Galaxy survey. Our work differs from previous analyses by focusing only on the observable part of the inflaton potential, or in other words, by making absolutely no assumption about extrapolation of the potential from its observable region to its minimum (i.e., about the branch of the potential responsible for the last ~50 inflationary e-folds). We only assume that inflation starts at least a few e-folds before the observable Universe leaves the Hubble radius, and that the inflaton rolls down a monotonic and regular potential, with no sharp features or phase transitions. We Taylor-expand the inflaton potential at order v=2, 3 or 4 in the vicinity of the pivot scale, compute the primordial spectra of scalar and tensor perturbations numerically and fit the data. For v>2, a large fraction of the allowed models is found to produce a large negative running of the scalar tilt, and to fall in a region of parameter space where the second-order slow-roll formalism is strongly inaccurate. We release a code for the computation of inflationary perturbations which is compatible with CosmoMC.

We consider the generation and evolution of quantum fluctuations of a massive nonsymmetric gravitational field (B-field) from inflationary epoch to matter era in the simplest variant of the nonsymmetric theory of gravitation (NGT), which consists of a gauge kinetic term and a mass term. We observe that quite generically a nonsymmetric metric field with mass, m_B ~ 0.03(H_I/10^(13) GeV)^4 eV, is a good dark matter candidate, where H_I denotes the inflationary scale. The most prominent feature of this dark matter is a peak in power at a comoving momentum, k ~ sqrt(m_B H_0)/(1+z_eq)^(1/4), where z_eq is the redshift at equality. This scale corresponds roughly to the Earth-Sun distance.

We show that during cosmological inflation the nonsymmetric metric tensor theory of gravitation develops a spectrum which is potentially observable by cosmic microwave background observations, and may be the most sensitive probe of the scale of cosmic inflation.