Constraints on the properties of the cosmological dark matter have previously been obtained in a model-independent fashion using the Generalised Dark Matter (GDM) framework. Here we extend that work in several directions: We consider the inclusion of WiggleZ matter power spectrum data, and show that this improves the constraints on the two perturbative GDM parameters, $c^2_s$ and $c^2_\text{vis}$, by a factor of 3, for a conservative choice of wavenumber range. A less conservative choice can yield an improvement of up to an order of magnitude compared to previous constraints. In order to examine the robustness of this result we develop a GDM halo model to explore how non-linear structure formation could proceed in this framework, since currently GDM has only been defined perturbatively and only linear theory has been used when generating constraints. We then examine how the halo model affects the constraints obtained from the matter power spectrum data. The less-conservative wavenumber range shows a significant difference between linear and non-linear modelling, with the latter favouring GDM parameters inconsistent with $\Lambda$CDM, underlining the importance of careful non-linear modelling when using this data. We also use this halo model to establish the robustness of previously obtained constraints, particularly those that involve weak gravitational lensing of the cosmic microwave background. Additionally, we show how the inclusion of neutrino mass as a free parameter affects previous constraints on the GDM parameters.

We perform the first test of Dark Matter (DM) stress-energy evolution through cosmic history, using Cosmic Microwave Background measurements supplemented with Baryon Acoustic Oscillation data and the Hubble Space Telescope key project data. We constrain the DM equation of state (EoS) in 8 redshift bins, and its sound speed and (shear) viscosity in 9 redshift bins, finding no convincing evidence for non-$\Lambda$CDM values in any of the redshift bins. Despite this enlarged parameter space, the sound speed and viscosity are constrained relatively well at late times (due to the inclusion of CMB lensing), whereas the EoS is most strongly constrained around recombination. These results constrain for the first time the level of ''coldness'' required of DM across various cosmological epochs at both the background and perturbative levels. We show that simultaneously allowing time dependence for both the EoS and sound speed parameters shifts the posterior of the DM abundance before recombination to a higher value, whilst keeping the present day DM abundance similar to the $\Lambda$CDM value. This shifts the posterior for the present day Hubble constant compared to $\Lambda$CDM, suggesting that DM with time-dependent parameters is well-suited to explore possible solutions to persistent tensions within the $\Lambda$CDM model.

We explore the use of field solvers as approximations of classical Vlasov-Poisson systems. This correspondence is investigated in both electrostatic and gravitational contexts. We demonstrate the ability of field solvers to be excellent approximations of problems with cold initial condition into the non linear regime. We also investigate extensions of the Schr\"odinger-Poisson system that employ multiple stacked cold streams, and the von Neumann-Poisson equation as methods that can successfully reproduce the classical evolution of warm initial conditions. We then discuss how appropriate simulation parameters need to be chosen to avoid interference terms, aliasing, and wave behavior in the field solver solutions. We present a series of criteria clarifying how parameters need to be chosen in order to effectively approximate classical solutions.

Cold Dark Matter (CDM) is a crucial constituent of the current concordance cosmological model. Having a vanishing equation of state (EoS), its energy density scales with the inverse cosmic volume and is thus uniquely described by a single number, its present abundance. We test the inverse cosmic volume law for Dark Matter (DM) by allowing its EoS to vary independently in eight redshift bins in the range $z=10^5$ and $z=0$. We use the latest measurements of the Cosmic Microwave Background radiation from the Planck satellite and supplement them with Baryon Acoustic Oscillation (BAO) data from the 6dF and SDSS-III BOSS surveys, and with the Hubble Space Telescope (HST) key project data. We find no evidence for nonzero EoS in any of the eight redshift bins. With Planck data alone, the DM abundance is most strongly constrained around matter-radiation equality $\omega^{\rm eq}_g = 0.1193^{+0.0036}_{-0.0035}$ (95% c.l.), whereas its present day value is more weakly constrained $\omega^{(0)}_g = 0.16^{+0.12}_{-0.10}$ (95% c.l.). Adding BAO or HST data does not significantly change the $\omega^{\rm eq}_g$ constraint, while $\omega^{(0)}_g$ tightens to $0.160^{+0.069}_{-0.065} $ (95% c.l.) and $0.124^{+0.081}_{-0.067}$ (95% c.l.) respectively. Our results constrain for the first time the level of "coldness" required of the DM across various cosmological epochs and show that the DM abundance is strictly positive at all times.

We study a new type of magnetoconvection in a nonuniform rotating plasma layer under a constant vertical magnetic field. To describe the weakly nonlinear stage of convection we apply Galerkin-truncated approximation and we obtain the system of equations of Lorentz type. A numerical analysis of these equations shows the presence of chaotic behavior of convective flows. Criteria for the appearance of chaotic motions are found depending on the convection parameters (Rayleigh number $\textrm{Ra}$), magnetic field (Chandrasekhar number $\textrm{Q}$), rotation (Taylor number $\textrm{Ta}$) for the Keplerian angular velocity profile $(\textrm{Ro}=-3/4)$ of the medium.

We demonstrate that the Vlasov equation describing collisionless self-gravitating matter may be solved with the so-called Schr\"odinger method (ScM). With the ScM, one solves the Schr\"odinger-Poisson system of equations for a complex wave function in d dimensions, rather than the Vlasov equation for a 2d-dimensional phase space density. The ScM also allows calculating the d-dimensional cumulants directly through quasi-local manipulations of the wave function, avoiding the complexity of 2d-dimensional phase space. We perform for the first time a quantitive comparison of the ScM and a conventional Vlasov solver in d=2 dimensions. Our numerical tests were carried out using two types of cold cosmological initial conditions: the classic collapse of a sine wave and those of a gaussian random field as commonly used in cosmological cold dark matter N-body simulations. We compare the first three cumulants, that is, the density, velocity and velocity dispersion, to those obtained by solving the Vlasov equation using the publicly available code ColDICE. We find excellent qualitative and quantitative agreement between these codes, demonstrating the feasibility and advantages of the ScM as an alternative to N-body simulations. We discuss, the emergence of effective vorticity in the ScM through the winding number around the points where the wave function vanishes. As an application we evaluate the background pressure induced by the non-linearity of large scale structure formation, thereby estimating the magnitude of cosmological backreaction. We find that it is negligibly small and has time dependence and magnitude compatible with expectations from the effective field theory of large scale structure.

In this manuscript using the asymptotic method of multiscale nonlinear theory we construct a nonlinear theory of the appearance of large-scale structures in the stratified conductive medium with the presence of small-scale oscillations of the velocity field and magnetic fields. These small-scale stationary oscillations are maintained by small external sources at low Reynolds numbers. We obtain a nonlinear system of equations describing the evolution of large-scale structures of the velocity field and magnetic fields. The linear stage of evolution leads to the known instability. In this article we consider the stationary large-scale structures of a magnetic field arising at stabilization of linear instability.

We examine how the properties of dark matter, parameterised by an equation of state parameter $w$ and two perturbative Generalised Dark Matter (GDM) parameters $c^2_s$ (the sound speed) and $c^2_\text{vis}$ (the viscosity), are constrained by existing cosmological data, particularly the Planck 2015 data release. We find that the GDM parameters are consistent with zero, and are strongly constrained, showing no evidence for extending the dark matter model beyond the Cold Dark Matter (CDM) paradigm. The dark matter equation of state is constrained to be within $-0.000896<w<0.00238$ at the $99.7\%$ confidence level, which is several times stronger than constraints found previously using WMAP data. The parameters $c^2_s$ and $c^2_\text{vis}$ are constrained to be less than $3.21\times10^{-6}$ and $6.06\times10^{-6}$ respectively at the $99.7\%$ confidence level. The inclusion of the GDM parameters does significantly affect the error bars on several $\Lambda$CDM parameters, notably the dimensionless dark matter density $\omega_g$ and the derived parameters $\sigma_8$ and $H_0$. This can be partially alleviated with the inclusion of data constraining the expansion history of the universe.

The Cold Dark Matter (CDM) model, wherein the dark matter is treated as a pressureless perfect fluid, provides a good fit to galactic and cosmological data. With the advent of precision cosmology, it should be asked whether this simplest model needs to be extended, and whether doing so could improve our understanding of the properties of dark matter. One established parameterisation for generalising the CDM fluid is the Generalised Dark Matter (GDM) model, in which dark matter is an imperfect fluid with pressure and shear viscosity that fulfill certain closure equations. We investigate these closure equations and the three new parametric functions they contain: the background equation of state w, the speed of sound c_s^2 and the viscosity c_{vis}^2. Taking these functions to be constant parameters, we analyse an exact solution of the perturbed Einstein equations in a GDM-dominated universe and discuss the main effects of the three parameters on the Cosmic Microwave Background (CMB). Our analysis suggests that the CMB alone is not able to distinguish between the GDM sound speed and viscosity parameters, but that other observables, such as the matter power spectrum, are required to break this degeneracy. We also consider other descriptions of imperfect fluids that have a non-perturbative definition and relate these to the GDM model. In particular, we consider scalar fields, an effective field theory (EFT) of fluids, an EFT of Large Scale Structure, non-equilibrium thermodynamics and tightly-coupled fluids. These descriptions could be used to extend the GDM model into the nonlinear regime of structure formation, which is necessary if the wealth of data available on those scales is to be employed in constraining the model. We also derive the initial conditions for adiabatic and isocurvature perturbations and provide the result in a form ready for implementation in Einstein-Boltzmann solvers.

We calculate the dark matter halo correlation function in redshift space using the Gaussian streaming model (GSM). To determine the scale dependent functions entering the streaming model we use local Lagrangian bias together with Convolution Lagrangian perturbation theory (CLPT) which constitutes an approximation to the Post-Zel'dovich approximation. On the basis of N-body simulations we demonstrate that a smoothing of the initial conditions with the Lagrangian radius improves the Zel'dovich approximation and its ability to predict the displacement field of proto-halos. Based on this observation we implement a "truncated" CLPT by smoothing the initial power spectrum and investigate the dependence of the streaming model ingredients on the smoothing scale. We find that the real space correlation functions of halos and their mean pairwise velocity are optimised if the coarse graining scale is chosen to be 1 Mpc/h at z=0, while the pairwise velocity dispersion is optimised if the smoothing scale is chosen to be the Lagrangian size of the halo. We compare theoretical results for the halo correlation function in redshift space to measurements within the Horizon Run 2 N-body simulation halo catalog. We find that this simple two-filter smoothing procedure in the spirit of the truncated Zel'dovich approximation significantly improves the GSM+CLPT prediction of the redshift space halo correlation function over the whole mass range from large galaxy to galaxy cluster-sized halos.

We derive the Edgeworth streaming model (ESM) for the redshift space correlation function starting from an arbitrary distribution function for biased tracers of dark matter by considering its two-point statistics and show that it reduces to the Gaussian streaming model (GSM) when neglecting non-Gaussianities. We test the accuracy of the GSM and ESM independent of perturbation theory using the Horizon Run 2 N-body halo catalog. While the monopole of the redshift space halo correlation function is well described by the GSM, higher multipoles improve upon including the leading order non-Gaussian correction in the ESM: the GSM quadrupole breaks down on scales below 30 Mpc/h whereas the ESM stays accurate to 2% within statistical errors down to 10 Mpc/h. To predict the scale dependent functions entering the streaming model we employ Convolution Lagrangian perturbation theory (CLPT) based on the dust model and local Lagrangian bias. Since dark matter halos carry an intrinsic length scale given by their Lagrangian radius, we extend CLPT to the coarse-grained dust model and consider two different smoothing approaches operating in Eulerian and Lagrangian space, respectively. The coarse-graining in Eulerian space features modified fluid dynamics different from dust while the coarse-graining in Lagrangian space is performed in the initial conditions with subsequent single streaming dust dynamics, implemented by smoothing the initial power spectrum in the spirit of the truncated Zel'dovich approximation. Finally, we compare the predictions of the different coarse-grained models for the streaming model ingredients to N-body measurements and comment on the proper choice of both the tracer distribution function and the smoothing scale. Since the perturbative methods we considered are not yet accurate enough on small scales, the GSM is sufficient when applied to perturbation theory.

We study the effect of coarse-graining the dynamics of a pressureless selfgravitating fluid (coarse-grained dust) in the context of cosmological perturbation theory, both in the Eulerian und Lagrangian framework. We obtain recursion relations for the Eulerian perturbation kernels of the coarse-grained dust model by relating them to those of the standard pressureless fluid model. The effect of the coarse-graining is illustrated by means of power and cross spectra for density and velocity that are computed up to 1-loop order. In particular, the large scale vorticity power spectrum that arises naturally from a mass-weighted velocity is derived from first principles. We find qualitatively good agreement of the magnitude, shape and spectral index of the vorticity power spectrum with recent measurements from N-body simulations and results from the effective field theory of large scale structure. To lay the ground for applications in the context of Lagrangian perturbation theory we finally describe how the kernels obtained in Eulerian space can be mapped to Lagrangian ones.

We investigate large scale structure formation of collisionless dark matter in the phase space description based on the Vlasov-Poisson equation. We present the Schr\"odinger method, originally proposed by Widrow and Kaiser (1993) as numerical technique based on the Schr\"odinger Poisson equation, as an analytical tool which is superior to the common standard pressureless fluid model. Whereas the dust model fails and develops singularities at shell crossing the Schr\"odinger method encompasses multi-streaming and even virialization.

We investigate large-scale structure formation of collisionless dark matter in the phase space description based on the Vlasov equation whose nonlinearity is induced solely by gravitational interaction according to the Poisson equation. Determining the time-evolution of density and peculiar velocity demands solving the full Vlasov hierarchy for the moments of the phase space distribution function. In the presence of long-range interaction no consistent truncation of the hierarchy is known apart from the pressure- less fluid (dust) model which is incapable of describing virialization due to the occurrence of shell-crossing singularities and the inability to generate vorticity and higher cumulants like velocity dispersion. Our goal is to find a simple ansatz for the phase space distribution function that approximates the full Vlasov distribution function without pathologies in a controlled way and therefore can serve as theoretical N-body double and as a replacement for the dust model. We argue that the coarse-grained Wigner probability distribution obtained from a wave function fulfilling the Schr\"odinger-Poisson equation (SPE) is the sought-after function. We show that its evolution equation approximates the Vlasov equation and therefore also the dust equations before shell-crossing, but cures the shell-crossing singularities and is able to describe regions of multi-streaming and virialization. This feature was already employed in simulations by Widrow & Kaiser (1993). The coarse-grained Wigner ansatz allows to calculate all higher moments from density and velocity analytically, thereby incorporating nonzero higher cumulants in a self-consistent manner. On this basis we are able to show that the Schr\"odinger method (ScM) automatically closes the corresponding hierarchy such that it suffices to solve the SPE in order to directly determine density and velocity and all higher cumulants.

We investigate the relation between the standard Newtonian equations for a pressureless fluid (dust) and the Einstein equations in a double expansion in small scales and small metric perturbations. We find that parts of the Einstein equations can be rewritten as a closed system of two coupled differential equations for the scalar and transverse vector metric perturbations in Poisson gauge. It is then shown that this system is equivalent to the Newtonian system of continuity and Euler equations. Brustein and Riotto (2011) conjectured the equivalence of these systems in the special case where vector perturbations were neglected. We show that this approach does not lead to the Euler equation but to a physically different one with large deviations already in the 1-loop power spectrum. We show that it is also possible to consistently set to zero the vector perturbations which strongly constrains the allowed initial conditions, in particular excluding Gaussian ones such that inclusion of vector perturbations is inevitable in the cosmological context. In addition we derive nonlinear equations for the gravitational slip and tensor perturbations, thereby extending Newtonian gravity of a dust fluid to account for nonlinear light propagation effects and dust-induced gravitational waves.

We compute the critical density of collapse for spherically symmetric overdensities in a class of f(R) modified gravity models. For the first time we evolve the Einstein, scalar field and non-linear fluid equations, making the minimal simplifying assumptions that the metric potentials and scalar field remain quasi-static throughout the collapse. Initially evolving a top hat profile, we find that the density threshold for collapse depends significantly on the initial conditions imposed, specifically the choice of size and shape. By imposing `natural' initial conditions, we obtain a fitting function for the spherical collapse delta_c as a function of collapse redshift, mass of the overdensity and f_{R0}, the background scalar field value at z=0. By extending delta_c into drifting and diffusing barrier within the context of excursion set theory, we obtain a realistic mass function that might be used to confront this class of scalar-tensor models with observations of dark matter halos. The proposed analytic formula for the halo mass function was tested against Monte Carlo random walks for a wide class of moving barriers and can therefore be applied to other modified gravity theories.

Contrary to recent claims in the literature, Newtonian N-body simulations of collisionless Dark Matter in a LambdaCDM background are compatible with general relativity and are not invalidated by general relativistic effects at the linear level. This verdict is based on four facts. (1) The system of linearized Einstein equations and conservation laws is well-posed in the gauge invariant formulation and physically meaningful. (2) Comparing general relativity with its Newtonian approximation at a given order in perturbation theory is only meaningful at the level of observables. (3) The dynamics of observables describing a dust fluid in general relativity and its Newtonian approximation agree at the linear level. Any disagreement for observables on the lightcone are well-known, of which the most dominant is gravitational lensing. (4) Large fluctuations in the Hubble parameter contribute significantly only to gravitational lensing effects. Therefore, these fluctuations are not in conflict with Newtonian N-body simulations beyond what has already been carefully taken into account using ray tracing technology.

Verlinde's heuristic argument for the interpretation of the standard Newtonian gravitational force as an entropic force is generalized by the introduction of a minimum temperature (or maximum wave length) for the microscopic degrees of freedom on the holographic screen. With the simplest possible setup, the resulting gravitational acceleration felt by a test mass m from a point mass M at a distance R is found to be of the form of the modified Newtonian dynamics (MOND) as suggested by Milgrom. The corresponding MOND-type acceleration constant is proportional to the minimum temperature, which can be interpreted as the Unruh temperature of an emerging de-Sitter space. This provides a possible explanation of the connection between local MOND-type two-body systems and cosmology.

Carr and Hawking showed that the proper size of a spherical overdense region surrounded by a flat FRW universe cannot be arbitrarily large as otherwise the region would close up on itself and become a separate universe. From this result they derived a condition connecting size and density of the overdense region ensuring that it is part of our universe. Carr used this condition to obtain an upper bound for the density fluctuation amplitude with the property that for smaller amplitudes the formation of a primordial black hole is possible, while larger ones indicate a separate universe. In contrast, we find that the appearance of a maximum is not a consequence of avoiding separate universes but arises naturally from the geometry of the chosen slicing. Using instead of density a volume fluctuation variable reveals that a fluctuation is a separate universe iff this variable diverges on superhorizon scales. Hence Carr's and Hawking's condition does not pose a physical constraint on density fluctuations. The dynamics of primordial black hole formation with an initial curvature fluctuation amplitude larger than the one corresponding to the maximum density fluctuation amplitude was previously not considered in detail and so we compare it to the well-known case where the amplitude is smaller by presenting embedding and conformal diagrams of both types in dust spacetimes.