Abstract. We present a new N-body code, gevolution, for the evolution of large scale structure in the Universe. Our code is based on a weak field expansion of General Relativity and calculates all six metric degrees of freedom in Poisson gauge. N-body particles are evolved by solving the geodesic equation which we write in terms of a canonical momentum such that it remains valid also for relativistic particles. We validate the code by considering the Schwarzschild solution and, in the Newtonian limit, by comparing with the Newtonian N-body codes Gadget-2 and RAMSES. We then proceed with a simulation of large scale structure in a Universe with massive neutrinos where we study the gravitational slip induced by the neutrino shear stress. The code can be extended to include different kinds of dark energy or modified gravity models and going beyond the usually adopted quasi-static approximation. Our code is publicly available.

We recently showed that postulated ultracompact minihalos with a steep density profile do not form in realistic simulations with enhanced initial perturbations. In this paper we assume that a small fraction of the dark matter consists of primordial black holes (PBHs) and simulate the formation of structures around them. We find that in this scenario halos with steep density profiles do form, consistent with theoretical predictions. If the rest of the dark matter consists of weakly interacting massive particles (WIMPs), we also show that WIMPs in the dense innermost part of halos surrounding the PBH would annihilate and produce a detectable gamma-ray signal. The non-detection of this signal implies that PBHs make up at most one billionth of the dark matter, provided that their mass is greater than one millionth of the mass of the Sun. Similarly, a detection of PBHs would imply that the remaining dark matter could not be WIMPs.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, ΛCDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of ΛCDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.Keywords: cosmology -dark energy -cosmological constant problem -modified gravitydark matter -early universe Cosmology has been both blessed and cursed by the establishment of a standard model: ΛCDM. On the one hand, the model has turned out to be extremely predictive, explanatory, and observationally robust, providing us with a substantial understanding of the formation of large-scale structure, the state of the early Universe, and the cosmic abundance of different types of matter and energy. It has also survived an impressive battery of precision observational tests -anomalies are few and far between, and their significance is contentious where they do arise -and its predictions are continually being vindicated through the discovery of new effects (B-mode polarization [1] and lensing [2,3] of the cosmic microwave background (CMB), and the kinetic Sunyaev-Zel'dovich effect [4] being some recent examples). These are the hallmarks of a good and valuable physical theory.On the other hand, the model suffers from profound theoretical difficulties. The two largest contributions to the energy content at late times -cold dark matter (CDM) and the cosmological constant (Λ) -have entirely mysterious physical origins. CDM has so far evaded direct detection by laboratory experiments, and so the particle field responsible for it -presumably a manifestation of "beyond the standard model" particle physics -is unknown. Curious discrepancies also appear to exist between the predicted clustering properties of CDM on small scales and observations. The cosmological constant is even more puzzling, giving rise to quite simply the biggest problem in all of fundamental physics: the question of why Λ appears to take such an unnatural value [5,6,7]. Inflation, the theory of the very early Universe, has also been criticized for being fine-tuned and under-predictive [8], and appears to leave many problems either unsolved or fundamentally unresolvable. These problems are indicative of a crisis.From January 14th-17th 2015, we held a conference in Oslo, Norway to surve...

Numerical simulations are a versatile tool providing insight into the complicated process of structure formation in cosmology [1]. This process is mainly governed by gravity, which is the dominant force on large scales. To date, a century after the formulation of general relativity [2], numerical codes for structure formation still employ Newton's law of gravitation. This approximation relies on the two assumptions that gravitational fields are weak and that they are only sourced by non-relativistic matter. While the former appears well justified on cosmological scales, the latter imposes restrictions on the nature of the "dark" components of the Universe (dark matter and dark energy) which are, however, poorly understood. Here we present the first simulations of cosmic structure formation using equations consistently derived from general relativity. We study in detail the small relativistic effects for a standard ΛCDM cosmology which cannot be obtained within a purely Newtonian framework. Our particle-mesh N-body code computes all six degrees of freedom of the metric and consistently solves the geodesic equation for particles, taking into account the relativistic potentials and the frame-dragging force. This conceptually clean approach is very general and can be applied to various settings where the Newtonian approximation fails or becomes inaccurate, ranging from simulations of models with dynamical dark energy [3] or warm/hot dark matter [4] to core collapse supernova explosions [5].The applicability of Newton's law of gravitation in the context of cosmic structure formation has been discussed extensively in the recent literature [6,7,8]. In particular, it is now well understood that this simplified description is quite accurate when applied within standard ΛCDM cosmology where perturbations come entirely from non-relativistic matter. However, the situation is not satisfactory for two reasons. Firstly, the quality of observational data is rapidly increasing, and upcoming galaxy surveys will eventually reach a precision where a naive treatment of the effects of spacetime geometry becomes insufficient [9,10]. Secondly, the true nature of dark matter and dark energy is not yet established. In order to study models beyond ΛCDM, some of which may feature relativistic sources of stress-energy (e.g. warm dark matter or dynamical dark energy), employing the Newtonian approximation is not always justified. A number of numerical codes has been developed for particular models [11,12,13,14], yet a general framework would be desirable. Furthermore, Newtonian gravity is acausal and not sensitive to the presence of a cosmological horizon. Even if a judicious interpretation of the output of Newtonian simulations can cure this problem at the linear level, it comes back when one goes beyond linear perturbation theory and it would be preferable to use the correct physics from the outset.Moving from the absolute space and time of the Newtonian picture towards a general relativistic view where geometry is dynamical poses ...

No abstract

Some of the dark matter in the Universe is made up of massive neutrinos. Their impact on the formation of large scale structure can be used to determine their absolute mass scale from cosmology, but to this end accurate numerical simulations have to be developed. Due to their relativistic nature, neutrinos pose additional challenges when one tries to include them in N-body simulations that are traditionally based on Newtonian physics. Here we present the first numerical study of massive neutrinos that uses a fully relativistic approach. Our N-body code, gevolution, is based on a weak-field formulation of general relativity that naturally provides a self-consistent framework for relativistic particle species. This allows us to model neutrinos from first principles, without invoking any ad-hoc recipes. Our simulation suite comprises some of the largest neutrino simulations performed to date. We study the effect of massive neutrinos on the nonlinear power spectra and the halo mass function, focusing on the interesting mass range between 0.06 eV and 0.3 eV and including a case for an inverted mass hierarchy.

Abstract. We present a framework for general relativistic N -body simulations in the regime of weak gravitational fields. In this approach, Einstein's equations are expanded in terms of metric perturbations about a Friedmann-Lemaître background, which are assumed to remain small. The metric perturbations themselves are only kept to linear order, but we keep their first spatial derivatives to second order and treat their second spatial derivatives as well as sources of stress-energy fully nonperturbatively. The evolution of matter is modelled by an N -body ensemble which can consist of free-streaming nonrelativistic (e.g. cold dark matter) or relativistic particle species (e.g. cosmic neutrinos), but the framework is fully general and also allows for other sources of stress-energy, in particular additional relativistic sources like modifiedgravity models or topological defects. We compare our method with the traditional Newtonian approach and argue that relativistic methods are conceptually more robust and flexible, at the cost of a moderate increase of numerical difficulty. However, for a ΛCDM cosmology, where nonrelativistic matter is the only source of perturbations, the relativistic corrections are expected to be small. We quantify this statement by extracting post-Newtonian estimates from Newtonian N -body simulations.PACS numbers: 98.80. Jk, 04.25.Nx

An important part of cosmological model fitting relies on correlating distance indicators of objects (for example type Ia supernovae) with their redshift, often illustrated on a Hubble diagram. Comparing the observed correlation with a homogeneous model is one of the key pieces of evidence for dark energy. The presence of cosmic structures introduces a bias and scatter, mainly due to gravitational lensing and peculiar velocities, but also due to smaller non-linear relativistic contributions which are more difficult to account for. For the first time we perform ray tracing onto halos in a relativistic N-body simulation. Our simulation is the largest that takes into account all leading-order corrections from general relativity in the evolution of structure, and we present a novel methodology for working out the non-linear projection of that structure onto the observer's past light cone. We show that the mean of the bias in the Hubble diagram is indeed as small as expected from perturbation theory. However, the distribution of sources is significantly skewed with a very long tail of highly magnified objects and we illustrate that the bias of cosmological parameters strongly depends on the function of distance which we consider.