Search

Anisotropies in the distance-redshift relation of cosmological sources are expected due to large-scale inhomogeneities in the local Universe. When the observed sources are tracing a large-scale matter flow in a general spacetime geometry, the distance-redshift relation with its anisotropies can be described with a geometrical prediction that generalises the well-known Friedmann-Lema\^itre-Robertson-Walker result. Furthermore, it turns out that a finite set of multipole coefficients contain the full information about a finite-order truncation of the distance-redshift relation of a given observer. The multipoles of the distance-redshift relation are interesting new cosmological observables that have a direct physical interpretation in terms of kinematical quantities of the underlying matter flow. Using light cones extracted from $N$-body simulations we quantify the anisotropies expected in a $\Lambda$ cold dark matter cosmology by running a Markov chain Monte Carlo analysis on the observed data. In this observational approach the survey selection implements an implicit smoothing scale over which the effective rest frame of matter is fitted. The perceived anisotropy therefore depends significantly on the redshift range and distribution of sources. We find that the multipoles of the expansion rate, as well as the observer's velocity with respect to the large-scale matter flow, can be determined robustly with our approach., Comment: 17 pages, 11 figures; v2: minor revision, accepted version

For most of cosmic history, the evolution of our Universe has been governed by the physics of a 'dark sector', consisting of dark matter and dark energy, whose properties are only understood in a schematic way. The influence of these constituents is mediated exclusively by the force of gravity, meaning that insight into their nature must be gleaned from gravitational phenomena. The advent of gravitational-wave astronomy has revolutionised the field of black hole astrophysics, and opens a new window of discovery for cosmological sources. Relevant examples include topological defects, such as domain walls or cosmic strings, which are remnants of a phase transition. Here we present the first simulations of cosmic structure formation in which the dynamics of the dark sector introduces domain walls as a source of stochastic gravitational waves in the late Universe. We study in detail how the spectrum of gravitational waves is affected by the properties of the model, and extrapolate the results to scales relevant to the recent evidence for a stochastic gravitational wave background. Our relativistic implementation of the field dynamics paves the way for optimal use of the next generation of gravitational experiments to unravel the dark sector., Comment: 10 pages, 3 figures. Supplementary material and videos may be found https://oyvindchristiansen.com/projects/symmetron_GWs. The code is available https://github.com/oyvach/AsGRD

The data on spectroscopic galaxy clustering collected by the Dark Energy Spectroscopic Instrument (DESI) will allow the significant detection of subtle features in the galaxy two-point correlation in redshift space, beyond the "standard" redshift-space distortions. Here we present an independent assessment of the detectability of the relativistic dipole in the cross-correlation of two populations of galaxies if they would be selected from the Bright Galaxy Survey (BGS) of DESI. We build synthetic galaxy catalogues with the characteristics of the BGS using the light cone of a relativistic $N$-body simulation. Exploring different ways of splitting the populations of galaxies we find that with an unequal split with more bright galaxies than faint galaxies the detectability is significantly boosted, reaching 19 $\sigma$ in the redshift bin $0.2 \lesssim z \lesssim 0.3$ and expected to be even higher at lower redshift. Moreover, we find that the measured dipole agrees very well with the prediction of relativistic effects from linear theory down to separations of $\sim$ 30 Mpc/$h$., Comment: version matching publication

Upcoming surveys of cosmic structures will probe scales close to the cosmological horizon, which opens up new opportunities for testing the cosmological concordance model to high accuracy. In particular, constraints on the squeezed bispectrum could rule out the single-field hypothesis during inflation. However, the squeezed bispectrum is also sensitive to dynamical effects of general relativity as well as interactions of matter with residual radiation from the early Universe. In this paper, we present a relativistic simulation pipeline that includes these relativistic effects consistently. We produce light cones and calculate the observed number counts of cold dark matter for five redshift bins between $z=0.55$-$2.25$. We compare the relativistic results against reference Newtonian simulations by means of angular power- and bispectra. We find that the dynamical relativistic effects scale roughly inversely proportional to the multipole in the angular power spectrum, with a maximum amplitude of $10\%$ for $\ell \lesssim 5$. By using a smoothing method applied to the binned bispectrum we detect the Newtonian bispectrum with very high significance. The purely relativistic part of the matter bispectrum is detected with a significance of $\sim 3\,\sigma$, mostly limited by cosmic variance. We find that the pure dynamical relativistic effects accounts for up to $3\%$ and $10\%$ of the total amplitude, respectively in the squeezed and equilateral limits. Our relativistic pipeline for modelling ultra-large scales yields gauge-independent results as we compute observables consistently on the past light cone, while the Newtonian treatment employs approximations that leave some residual gauge dependence. A gauge-invariant approach is required in order to meet the expected level of precision of forthcoming probes of cosmic structures on ultra-large scales., Comment: 39 pages, 15 figures

We present the halo number counts and its two-point statistics, the observable angular power spectrum, extracted for the first time from relativistic N-body simulations. The halo catalogues used in this work are built from the relativistic N-body code gevolution, and the observed redshift and angular positions of the sources are computed using a non-perturbative ray-tracing method, which includes all relativistic scalar contributions to the number counts. We investigate the validity and limitations of the linear bias prescription to describe our simulated power spectra. In particular, we assess the consistency of different bias measurements on large scales, and we estimate up to which scales a linear bias is accurate in modelling the data, within the statistical errors. We then test a second-order perturbative bias expansion for the angular statistics, on a range of redshifts and scales previously unexplored in this context, that is $0.4 \le \bar{z} \le 2$ up to scales $\ell_\mathrm{max} \sim 1000$. We find that the angular power spectra at equal redshift can be modelled with high accuracy with a minimal extension of the number of bias parameters, that is using a two-parameter model comprising linear bias and tidal bias. We show that this model performs significantly better than a model without tidal bias but with quadratic bias as extra degree of freedom, and that the latter is inaccurate at $\bar{z} \ge 0.7$. Finally, we extract from our simulations the cross-correlation of halo number counts and lensing convergence. We show that the estimate of the linear bias from this cross-correlation is consistent with the measurements based on the clustering statistics alone, and that it is crucial to take into account the effect of magnification in the halo number counts to avoid systematic shifts in the computed bias., Comment: 28 pages, 16 figures

Cosmological models beyond $\Lambda$CDM, like those featuring massive neutrinos or modifications of gravity, often display a characteristic change (scale-dependent suppression or enhancement) in the matter power spectrum when compared to a $\Lambda$CDM baseline. It is therefore a widely held view that constraints on those models can be obtained by searching for such features in the clustering statistics of large-scale structure. However, when using biased tracers of matter in the analysis, the situation is complicated by the fact that the bias also depends on cosmology. Here we investigate how the selection of tracers affects the observed signatures for two examples of beyond-$\Lambda$CDM cosmologies: massive neutrinos and clustering dark energy ($k$-essence). We study the signatures in the monopole, quadrupole, and hexadecapole of the redshift-space power spectra for halo catalogues from large $N$-body simulations and argue that a fixed selection criterion based on local attributes like tracer mass leads to a near loss of signal in most cases. Instead, the full signal is recovered only if the selection of tracers is done at fixed bias. This emphasises the need to model or measure the bias parameters accurately in order to get meaningful constraints on the cosmological model., Comment: 5 pages, 5 figures; data available at https://zenodo.org/record/7832261

In this paper, we study the effective field theory (EFT) of dark energy for the $k$-essence model beyond linear order. Using particle-mesh $N$-body simulations that consistently solve the dark energy evolution on a grid, we find that the next-to-leading order in the EFT expansion, which comprises the terms of the equations of motion that are quadratic in the field variables, gives rise to a new instability in the regime of low speed of sound (high Mach number). We rule out the possibility of a numerical artefact by considering simplified cases in spherically and plane symmetric situations analytically. If the speed of sound vanishes exactly, the non-linear instability makes the evolution singular in finite time, signalling a breakdown of the EFT framework. The case of finite (but small) speed of sound is subtle, and the local singularity could be replaced by some other type of behaviour with strong non-linearities. While an ultraviolet completion may cure the problem in principle, there is no reason why this should be the case in general. As a result, for a large range of the effective speed of sound $c_s$, a linear treatment is not adequate., Comment: 24 pages, 10 figures

In this letter we introduce the non-linear partial differential equation (PDE) $\partial^2_{\tau} \pi \propto (\vec\nabla \pi)^2$ showing a new type of instability. Such equations appear in the effective field theory (EFT) of dark energy for the $k$-essence model as well as in many other theories based on the EFT formalism. We demonstrate the occurrence of instability in the cosmological context using a relativistic $N$-body code, and we study it mathematically in 3+1 dimensions within spherical symmetry. We show that this term dominates for the low speed of sound limit where some important linear terms are suppressed., Comment: 5 pages, 2 figures; data available at https://zenodo.org/record/7950024

In this paper we present for the first time the angular power spectra $C_\ell(z,z')$ for number counts from relativistic N-body simulations. We use the relativistic N-body code gevolution with its exact integration of lightlike geodesics which include all relativistic scalar contributions to the number counts. We compare our non-perturbative numerical results with the results from CLASS using the HMCODE approximation for the non-linear matter power spectrum. We find that this simple description is excellent for both, the density and the convergence. On the other hand, the current implementation of redshift-space distortions in Boltzmann codes is not accurate. We also find that the largest contribution to the unequal-redshift power spectra is the cross-correlation of the density and the lensing contribution to the number counts, especially for redshift bins that are far apart. Correlating the number counts with the convergence map we find that the signal is dominated by the lensing-lensing term when the convergence field redshift is not higher than the number counts one, while it is dominated by the density-lensing term in the opposite case. In the present study, the issue of galaxy bias is deliberately left aside by considering only unbiased samples of matter particles from the simulations., Comment: 26 pages, 8 figures, version accepted for publication in JCAP

The standard cosmological model is inherently relativistic, and yet a wide range of cosmological observations can be predicted accurately from essentially Newtonian theory. This is not the case on `ultra-large' distance scales, around the cosmic horizon size, however, where relativistic effects can no longer be neglected. In this paper, we present a novel suite of 53 fully relativistic simulations generated using the gevolution code, each covering the full sky out to $z \approx$ 0.85, and approximately 1930 square degrees out to $z \approx$ 3.55. These include a relativistic treatment of massive neutrinos, as well as the gravitational potential that can be used to exactly calculate observables on the past light cone. The simulations are divided into two sets, the first being a set of 39 simulations of the same fiducial cosmology (based on the Euclid Flagship 2 cosmology) with different realisations of the initial conditions, and the second which fixes the initial conditions, but varies each of seven cosmological parameters in turn. Taken together, these simulations allow us to perform statistical studies and calculate derivatives of any relativistic observable with respect to cosmological parameters. As an example application, we compute the cross-correlation between the Doppler magnification term in the convergence, $\kappa_v$, and the CDM+baryon density contrast, $\delta_{\rm cb}$, which arises only in a (special) relativistic treatment. We are able to accurately recover this term as predicted by relativistic perturbation theory, and study its sample variance and derivatives with respect to cosmological parameters., Comment: 9 pages, 6 figures, corresponds to published version

Planned efforts to probe the largest observable distance scales in future cosmological surveys are motivated by a desire to detect relic correlations left over from inflation, and the possibility of constraining novel gravitational phenomena beyond General Relativity (GR). On such large scales, the usual Newtonian approaches to modelling summary statistics like the power spectrum and bispectrum are insufficient, and we must consider a fully relativistic and gauge-independent treatment of observables such as galaxy number counts in order to avoid subtle biases, e.g. in the determination of the $f_{\rm NL}$ parameter. In this work, we present an initial application of an analysis pipeline capable of accurately modelling and recovering relativistic spectra and correlation functions. As a proof of concept, we focus on the non-zero dipole of the redshift-space power spectrum that arises in the cross-correlation of different mass bins of dark matter halos, using strictly gauge-independent observable quantities evaluated on the past light cone of a fully relativistic N-body simulation in a redshift bin $1.7 \le z \le 2.9$. We pay particular attention to the correct estimation of power spectrum multipoles, comparing different methods of accounting for complications such as the survey geometry (window function) and evolution/bias effects on the past light cone, and discuss how our results compare with previous attempts at extracting novel GR signatures from relativistic simulations., Comment: Minor corrections in the text of version 2

We study cosmological observables on the past light cone of a fixed observer in the context of clustering dark energy. We focus on observables that probe the gravitational field directly, namely the integrated Sachs-Wolfe and non-linear Rees-Sciama effect (ISW-RS), weak gravitational lensing, gravitational redshift and Shapiro time delay. With our purpose-built $N$-body code "$k$-evolution" that tracks the coupled evolution of dark matter particles and the dark energy field, we are able to study the regime of low speed of sound $c_s$ where dark energy perturbations can become quite large. Using ray tracing we produce two-dimensional sky maps for each effect and we compute their angular power spectra. It turns out that the ISW-RS signal is the most promising probe to constrain clustering dark energy properties coded in $w-c_s^2$, as the $\textit{linear}$ clustering of dark energy would change the angular power spectrum by $\sim 30\%$ at low $\ell$ when comparing two different speeds of sound for dark energy. Weak gravitational lensing, Shapiro time-delay and gravitational redshift are less sensitive probes of clustering dark energy, showing variations of a few percent only. The effect of dark energy $\textit{non-linearities}$ in all the power spectra is negligible at low $\ell$, but reaches about $2\%$ and $3\%$, respectively, in the convergence and ISW-RS angular power spectra at multipoles of a few hundred when observed at redshift $\sim 0.85$. Future cosmological surveys achieving percent precision measurements will allow to probe the clustering of dark energy to a high degree of confidence., Comment: 18 pages, 15 figures, version accepted for publication in MNRAS; data available at https://doi.org/10.5281/zenodo.7950184

A number of codes for general-relativistic simulations of cosmological structure formation have been developed in recent years. Here we demonstrate that a sample of these codes produce consistent results beyond the Newtonian regime. We simulate solutions to Einstein's equations dominated by gravitomagnetism -- a vector-type gravitational field that doesn't exist in Newtonian gravity and produces frame-dragging, the leading-order post-Newtonian effect. We calculate the coordinate-invariant effect on intersecting null geodesics by performing ray tracing in each independent code. With this observable quantity, we assess and compare each code's ability to compute relativistic effects., Comment: 14 pages, 5 figures, accepted to CQG

We present a numerical weak-lensing analysis that is fully relativistic and non-perturbative for the scalar part of the gravitational potential and first-order in the vector part, frame dragging. Integrating the photon geodesics backwards from the observer to the emitters, we solve the Sachs optical equations and study in detail the weak-lensing convergence, ellipticity and rotation. For the first time, we apply such an analysis to a high-resolution relativistic N-body simulation, which consistently includes the leading-order corrections due to general relativity on both large and small scales. These are related to the question of gauge choice and to post-Newtonian corrections, respectively. We present the angular power spectra and one-point probability distribution functions for the weak-lensing variables, which we find are broadly in agreement with comparable Newtonian simulations. Our geometric approach, however, is more robust and flexible, and can therefore be applied consistently to non-standard cosmologies and modified theories of gravity., Comment: 20 pages, 17 figures, corresponds to published version

The late integrated Sachs-Wolfe (ISW) effect correlates the Cosmic Microwave Background (CMB) temperature anisotropies with foreground cosmic large-scale structures. As the correlation depends crucially on the growth history in the era of dark energy, it is a key observational probe for constraining the cosmological model. Here we present a detailed study based on full-sky and deep light cones generated from very large volume numerical N-body simulations, which allow us to avoid the use of standard replica techniques, while capturing the entirety of the late ISW effect on the large scales. We post-process the light cones using an accurate ray-tracing method and construct full-sky maps of the ISW temperature anisotropy for three different dark energy models. We quantify in detail the extent to which the ISW effect can be used to discriminate between different dark energy scenarios when cross-correlated with the matter distribution or the CMB lensing potential. We also investigate the onset of non-linearities, the so-called Rees-Sciama effect which provides a complementary probe of the dark sector. We find the signal of the lensing-lensing and ISW-lensing correlation of the three dark energy models to be consistent with measurements from the Planck satellite. Future surveys of the large-scale structures may provide cross-correlation measurements that are sufficiently precise to distinguish the signal of these models. Our methodology is very general and can be applied to any dark energy or modified gravity scenario as long as the metric seen by photons can still be characterized by a Weyl potential., Comment: 13 pages, 8 figures; v2: minor revision, matches accepted version

In this paper, we quantify the non-linear effects from $k$-essence dark energy through an effective parameter $\mu$ that encodes the additional contribution of a dark energy fluid or a modification of gravity to the Poisson equation. This is a first step toward quantifying non-linear effects of dark energy/modified gravity models in a more general approach. We compare our $N$-body simulation results from $k$-evolution with predictions from the linear Boltzmann code $\texttt{CLASS}$, and we show that for the $k$-essence model one can safely neglect the difference between the two potentials, $ \Phi -\Psi$, and short wave corrections appearing as higher order terms in the Poisson equation, which allows us to use single parameter $\mu$ for characterizing this model. We also show that for a large $k$-essence speed of sound the $\texttt{CLASS}$ results are sufficiently accurate, while for a low speed of sound non-linearities in matter and in the $k$-essence field are non-negligible. We propose a $\tanh$-based parameterisation for $\mu$, motivated by the results for two cases with low ($c_s^2=10^{-7}$) and high ($c_s^2=10^{-4}$) speed of sound, to include the non-linear effects based on the simulation results. This parametric form of $\mu$ can be used to improve Fisher forecasts or Newtonian $N$-body simulations for $k$-essence models., Comment: 18 pages, 8 figures. Updated to match version published in JCAP; data available at https://doi.org/10.5281/zenodo.7950666

We introduce $k$-evolution, a relativistic $N$-body code based on $\textit{gevolution}$, which includes clustering dark energy among its cosmological components. To describe dark energy, we use the effective field theory approach. In particular, we focus on $k$-essence with a speed of sound much smaller than unity but we lay down the basis to extend the code to other dark energy and modified gravity models. We develop the formalism including dark energy non-linearities but, as a first step, we implement the equations in the code after dropping non-linear self-coupling in the $k$-essence field. In this simplified setup, we compare $k$-evolution simulations with those of $\texttt{CLASS}$ and $\textit{gevolution}$ 1.2, showing the effect of dark matter and gravitational non-linearities on the power spectrum of dark matter, of dark energy and of the gravitational potential. Moreover, we compare $k$-evolution to Newtonian $N$-body simulations with back-scaled initial conditions and study how dark energy clustering affects massive halos., Comment: 38 pages, 19 figures. Final version, accepted for publication in JCAP; data available at https://doi.org/10.5281/zenodo.7953189

We clarify the subtle issue of finding the correct mapping of Newtonian simulations to light-cone observables at very large distance scales. A faithful general-relativistic interpretation specifies a gauge, i.e. a chart that relates the simulation data to points of the space-time manifold. It has already been pointed out that the implicit gauge choice of Newtonian simulations is indeed different from the Poisson gauge that is commonly adopted for relativistic calculations, the difference being most significant at large scales. It is therefore inconsistent, for example, to predict weak-lensing observables from simulations unless this gauge issue is properly accounted for. Using perturbation theory as well as fully relativistic N-body simulations we quantify the systematic error introduced this way, and we discuss several solutions that would render the calculations relativistically self-consistent., Comment: 10 pages, 5 figures; v2: minor revision with additional content, matches accepted manuscript

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 these halos 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., Comment: 10 pages, 8 figures. v2: Minor changes to match the published version

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., Comment: 6 pages, 3 figures; v2: some additional analysis, matches accepted version

Clustering of a perfect fluid does not lead to the generation of vorticity. It is the collisionless nature of dark matter, inducing velocity dispersion and shell crossing, which is at the origin of cosmological vorticity generation. In this paper we investigate the generation of vorticity during the formation of cosmological large scale structure using the public relativistic N-body code $\textit{gevolution}$. We test several methods to compute the vorticity power spectrum and we study its convergence with respect to the mass and grid resolution of our simulations. We determine the power spectrum, the spectral index on large-scales, the amplitude of the peak position and their time evolution. We also compare the vorticity extracted from our simulations with the vector perturbations of the metric. Our results are accompanied by resolution studies and compared with previous studies in the literature., Comment: 26 pages, 13 figures

We perform three-dimensional simulations of structure formation in the early Universe, when boosting the primordial power spectrum on approximately kpc scales. We demonstrate that our simulations are capable of producing power-law profiles close to the steep $\rho\propto r^{-9/4}$ halo profiles that are commonly assumed to be a good approximation to ultracompact minihalos (UCMHs). However, we show that for more realistic initial conditions in which halos are neither perfectly symmetric nor isolated, the steep power-law profile is disrupted and we find that the Navarro-Frenk-White profile is a better fit to most halos. In the presence of background fluctuations even extreme, nearly spherical initial conditions do not remain exceptional. Nonetheless, boosting the amplitude of initial fluctuations causes all structures to form earlier and thus at larger densities. With sufficiently large amplitude of fluctuations we find that values for the concentration of typical halos in our simulations can become very large. However, despite the signal coming from dark matter annihilation inside the cores of these halos being enhanced, it is still orders-of-magnitude smaller compared to the usually assumed UCMH profile. The upper bound on the primordial power spectrum from the non-observation of UCMHs should therefore be re-evaluated., Comment: 19 pages, 20 figures. v2: Minor changes to match the published version

In order to keep pace with the increasing data quality of astronomical surveys the observed source redshift has to be modeled beyond the well-known Doppler contribution. In this letter I want to examine the gauge issue that is often glossed over when one assigns a perturbed redshift to simulated data generated with a Newtonian N-body code. A careful analysis reveals the presence of a correction term that has so far been neglected. It is roughly proportional to the observed length scale divided by the Hubble scale and therefore suppressed inside the horizon. However, on gigaparsec scales it can be comparable to the gravitational redshift and hence amounts to an important relativistic effect., Comment: 5 pages; v2: minor revision, accepted version

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., Comment: 26 pages + appendix, 17 figures; v2: minor revision, published version; code available at https://github.com/gevolution-code/gevolution-1.1

A persistent theme in the study of dark energy is the question of whether it really exists or not. It is often claimed hat we are mis-calculating the cosmological model by neglecting the effects associated with averaging over large-scale structures. In the Newtonian approximation this is clear: there is no effect. Within the full relativistic picture this remains an important open question however, owing to the complex mathematics involved. We study this issue using particle numerical simulations which account for all relevant relativistic effects without any problems from shell crossing. In this context we show for the first time that the backreaction from structure can differ by many orders of magnitude depending upon the slicing of spacetime one chooses to average over. In the worst case, where smoothing is carried out in synchronous spatial surfaces, the corrections can reach ten percent and more. However, when smoothing on the constant time hypersurface of the Newtonian gauge, backreaction contributions remain 3-5 orders of magnitude smaller., Comment: 13 pages, 3 figures; updated to match published version (extended and improved explanations, conclusions remain the same)

Newtonian N-body simulations have been employed successfully over the past decades for the simulation of the cosmological large-scale structure. Such simulations usually ignore radiation perturbations (photons and massless neutrinos) and the impact of general relativity (GR) beyond the background expansion. This approximation can be relaxed and we discuss three different approaches that are accurate to leading order in GR. For simulations that start at redshift less than about 100 we find that the presence of early radiation typically leads to percent-level effects on the numerical power spectra at large scales. Our numerical results agree across the three methods, and we conclude that all of the three methods are suitable for simulations in a standard cosmology. Two of the methods modify the N-body evolution directly, while the third method can be applied as a post-processing prescription., Comment: 10 pages + appendix, 4 figures - v2: typos corrected; corresponds to published version

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., Comment: 28 pages + appendix, 10 figures. v2: revised and extended version accepted by JCAP; code available at https://github.com/gevolution-code

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$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 $\Lambda$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., Comment: 99 pages, 8 figures. Version published in Physics of the Dark Universe

We compute the angular power spectra of the E-type and B-type lensing potentials for gravitational waves from inflation and for tensor perturbations induced by scalar perturbations. We derive the tensor-lensed CMB power spectra for both cases. We also apply our formalism to determine the linear lensing potential for a Bianchi I spacetime with small anisotropy., Comment: 19 pages, 5 figures; v2: discussions about non-linear contributions and quadrupole regularization extended, minor typos corrected, references added, matches published version

Numerical simulations are a versatile tool providing insight into the complicated process of structure formation in cosmology. 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, 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 {\Lambda}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 or warm/hot dark matter to core collapse supernova explosions., Comment: 6 pages, 3 figures

Within a cosmological context, we study the behaviour of collisionless particles in the weak field approximation to General Relativity, allowing for large gradients of the fields and relativistic velocities for the particles. We consider a spherically symmetric setup such that high resolution simulations are possible with minimal computational resources. We test our formalism by comparing it to two exact solutions: the Schwarzschild solution and the Lema\^itre-Tolman-Bondi model. In order to make the comparison we consider redshifts and lensing angles of photons passing through the simulation. These are both observable quantities and hence are gauge independent. We demonstrate that our scheme is more accurate than a Newtonian scheme, correctly reproducing the leading-order post-Newtonian correction. In addition, our setup is able to handle shell-crossings, which is not possible within a fluid model. Furthermore, by introducing angular momentum, we find configurations corresponding to bound objects which may prove useful for numerical studies of the effects of modified gravity, dynamical dark energy models or even compact bound objects within General Relativity., Comment: 14 pages, 7 figures

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\^itre 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 non-perturbatively. 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 modified-gravity 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 LambdaCDM 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., Comment: 30 pages, 3 figures. Invited contribution to a Classical and Quantum Gravity focus issue on "Relativistic Effects in Cosmology", edited by Kazuya Koyama

The large-scale homogeneity and isotropy of the universe is generally thought to imply a well defined background cosmological model. It may not. Smoothing over structure adds in an extra contribution, transferring power from small scales up to large. Second-order perturbation theory implies that the effect is small, but suggests that formally the perturbation series may not converge. The amplitude of the effect is actually determined by the ratio of the Hubble scales at matter-radiation equality and today - which are entirely unrelated. This implies that a universe with significantly lower temperature today could have significant backreaction from more power on small scales, and so provides the ideal testing ground for understanding backreaction. We investigate this using two different N-body numerical simulations - a 3D Newtonian and a 1D simulation which includes all relevant relativistic effects. We show that while perturbation theory predicts an increasing backreaction as more initial small-scale power is added, in fact the virialisation of structure saturates the backreaction effect at the same level independently of the equality scale. This implies that backreaction is a small effect independently of initial conditions. Nevertheless, it may still contribute at the percent level to certain cosmological observables and therefore it cannot be neglected in precision cosmology., Comment: 6 pages, 2 figures; v2 matches published version

Distance measurements are usually thought to probe the background metric of the universe, but in reality the presence of perturbations will lead to deviations from the result expected in an exactly homogeneous and isotropic universe. At least in principle the presence of perturbations could even explain the observed distance-redshift relation without the need for dark energy. In this paper we re-investigate a toy model where perturbations are plane symmetric, and for which exact solutions are known in the fluid limit. However, if perturbations are large, shell-crossing occurs and the fluid approximation breaks down. This prevents the study of the most interesting cases. Here we use a general-relativistic N-body simulation that does not suffer from this problem and which allows us to go beyond previous works. We show that even for very large plane-symmetric perturbations we are not able to mimic the observed distance-redshift relation. We also discuss how the synchronous comoving gauge breaks down when shell-crossing occurs, while metric perturbations in the longitudinal gauge remain small. For this reason the longitudinal (Newtonian) gauge appears superior for relativistic N-body simulations of large-scale structure formation., Comment: 10 pages, 6 figures; v2: minor modifications to match printed version

We develop a formalism for General Relativistic N-body simulations in the weak field regime, suitable for cosmological applications. The problem is kept tractable by retaining the metric perturbations to first order, the first derivatives to second order and second derivatives to all orders, thus taking into account the most important nonlinear effects of Einstein gravity. It is also expected that any significant "backreaction" should appear at this order. We show that the simulation scheme is feasible in practice by implementing it for a plane-symmetric situation and running two test cases, one with only cold dark matter, and one which also includes a cosmological constant. For these plane-symmetric situations, the deviations from the usual Newtonian N-body simulations remain small and, apart from a non-trivial correction to the background, can be accurately estimated within the Newtonian framework. The correction to the background scale factor, which is a genuine "backreaction" effect, can be robustly obtained with our algorithm. Our numerical approach is also naturally suited for the inclusion of extra relativistic fields and thus for dark energy or modified gravity simulations., Comment: 13 pages + appendix, 5 figures; v2: minor modifications to match printed version

For free fields, pair creation in expanding universes is associated with the building up of correlations that lead to nonseparable states, i.e., quantum mechanically entangled ones. For dissipative fields, i.e., fields coupled to an environment, there is a competition between the squeezing of the state and the coupling to the external bath. We compute the final coherence level for dissipative fields that propagate in a two-dimensional de Sitter space, and we characterize the domain in parameter space where the state remains nonseparable. We then apply our analysis to (analogue) Hawking radiation by exploiting the close relationship between Lorentz violating theories propagating in de Sitter and black hole metrics. We establish the robustness of the spectrum and find that the entanglement among Hawking pairs is generally much stronger than that among pairs of quanta with opposite momenta., Comment: Final version published in prd, 22 pages

We examine the mode functions of the electromagnetic field on spherically symmetric backgrounds with special attention to the subclass which allows for the foliation as open Friedmann-Lemaitre (FL) spacetime. It is well-known that in certain scalar field theories on open FL background there can exist so-called supercurvature modes, their existence depending on parameters of the theory. Looking at specific open universe models, such as open inflation and the Milne universe, we find that no supercurvature modes are present in the spectrum of the electromagnetic field. This excludes the possibility for superadiabatic evolution of cosmological magnetic fields within these models without relying on new physics or breaking the conformal invariance of electromagnetism., Comment: 6 pages, 1 figure; v2: minor revision, appendix added, accepted for publication in MNRAS

We study a homogeneous and nearly-isotropic Universe permeated by a homogeneous magnetic field. Together with an isotropic fluid, the homogeneous magnetic field, which is the primary source of anisotropy, leads to a plane-symmetric Bianchi I model of the Universe. However, when free-streaming relativistic particles are present, they generate an anisotropic pressure which counteracts the one from the magnetic field such that the Universe becomes isotropized. We show that due to this effect, the CMB temperature anisotropy from a homogeneous magnetic field is significantly suppressed if the the neutrino masses are smaller than 0.3 eV., Comment: 20 pages, 5 figures (JHEP style) - submitted to JCAP

In a landscape of compactifications with different numbers of macroscopic dimensions, it is possible that our universe has nucleated from a vacuum where some of our four large dimensions were compact while other, now compact, directions were macroscopic. From our perspective, this shapeshifting can be perceived as an anisotropic background spacetime. As an example, we present a model where our universe emerged from a tunneling event which involves the decompactification of two dimensions compactified on the two-sphere. In this case, our universe is of the Kantowski-Sachs type and therefore homogeneous and anisotropic. We study the deviations from statistical isotropy of the Cosmic Microwave Background induced by the anisotropic curvature, with particular attention to the anomalies. The model predicts a quadrupolar power asymmetry with the same sign and acoustic oscillations as found by WMAP. The amplitude of the effect is however too small given the current estimated bound on anisotropic curvature derived from the quadrupole., Comment: 17 pages, 4 figures. v3: extended discussion of CMB anomalies; version accepted by Phys. Rev. D

Motivated by cosmological first-order phase transitions we examine the nucleation and evolution of vacuum bubbles in non-vacuum environments. Non-standard backgrounds can be relevant in the context of rapid tunneling processes on the landscape. Utilising complex time methods, we show that tunneling rates can be notably modified in the case of dynamical FRW backgrounds. We give a classification of the importance of the effect in terms of the relevant dynamical time scales. For both the bubble nucleation and evolution analysis we make use of the thin-wall approximation. From the classical bubble evolution on homogeneous matter backgrounds via the junction method, we find that the inflation of vacuum bubbles is very sensitive to the presence of ambient matter and quantify this statement. We also employ inhomogeneous matter models (LTB) and models that undergo a rapid phase transition (FRW) as a background and discuss in which cases potentially observable imprints on the bubble trajectory can remain., Comment: 8 pages, 4 figures; proceedings (PoS) of the International Workshop on Cosmic Structure and Evolution, Bielefeld, Germany, 23-25 September 2009

In this letter we discuss cosmological first order phase transitions with de Sitter bubbles nucleating on (inhomogeneous) matter backgrounds. The de Sitter bubble can be a toy model for an inflationary phase of universes like our own. Using the thin wall approximation and the Israel junction method we trace the classical evolution of the formed bubbles within a compound model. We first address homogeneous ambient space (FRW model) and already find that bubbles nucleated in a dust dominated background cannot expand. For an inhomogeneous dust background (LTB model) we describe cases with at least initially expanding bubbles. Yet, an ensuing passage of the bubble wall through ambient curvature inhomogeneities remains unnoticed for observers inside the bubble. Notable effects also for interior observers are found in the case of a rapid background phase transition in a FRW model., Comment: 4 pages, 2 figures, proceedings (Annalen der Physik) of the Grassmannian Conference in Fundamental Cosmology (Grasscosmofun'09) Szczecin, Poland, 14-19 September 2009

In the context of bubble universes produced by a first-order phase transition with large nucleation rates compared to the inverse dynamical time scale of the parent bubble, we extend the usual analysis to non-vacuum backgrounds. In particular, we provide semi-analytic and numerical results for the modified nucleation rate in FLRW backgrounds, as well as a parameter study of bubble walls propagating into inhomogeneous (LTB) or FLRW spacetimes, both in the thin-wall approximation. We show that in our model, matter in the background often prevents bubbles from successful expansion and forces them to collapse. For cases where they do expand, we give arguments why the effects on the interior spacetime are small for a wide range of reasonable parameters and discuss the limitations of the employed approximations., Comment: 29 pages, 8 figures, typos corrected, matches published version

The sensitivity of inflationary spectra to initial conditions is addressed in the context of a phenomenological model that breaks Lorentz invariance by dissipative effects above some threshold energy $\Lambda$. These effects are obtained dynamically by coupling the fluctuation modes to extra degrees of freedom which are unobservable below $\Lambda$. Because of the strong dissipative effects in the early propagation, only the state of the extra degrees of freedom is relevant for the power spectrum. If this state is the ground state, and if $\Lambda$ is much larger than the Hubble scale $H$, the standard spectrum is recovered. Using analytical and numerical methods, we calculate the modifications for a large class of dissipative models. For all of these, we show that the leading modification (in an expansion in $H/\Lambda$) is linear in the decay rate evaluated at horizon exit, and that high frequency superimposed oscillations are not generated. The modification is negative when the decay rate decreases slower than the cube of $H$, which means that there is a loss of power on the largest scales., Comment: 12 pages, 6 figures. Minor changes after refereeing and editing

The non-equilibrium phase transition in models for epidemic spreading with long-range infections in combination with incubation times is investigated by field-theoretical and numerical methods. Here the spreading process is modelled by spatio-temporal Levy flights, i.e., it is assumed that both spreading distance and incubation time decay algebraically. Depending on the infection rate one observes a phase transition from a fluctuating active phase into an absorbing phase, where the infection becomes extinct. This transition between spreading and extinction is characterized by continuously varying critical exponents, extending from a mean-field regime to a phase described by the universality class of directed percolation. We compute the critical exponents in the vicinity of the upper critical dimension by a field-theoretic renormalization group calculation and verify the results in one spatial dimension by extensive numerical simulations., Comment: 20 pages, 7 figures