Source author record

Francis Bernardeau

Francis Bernardeau appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2025arXiv

The Topology of Rayleigh-Levy Flights in Two Dimensions

Rayleigh-Lévy flights are simplified cosmological tools which capture certain essential statistical properties of the cosmic density field, including hierarchical structures in higher-order correlations, making them a valuable reference for studying the highly non-linear regime of structure formation. Unlike standard Markovian processes, they exhibit long-range correlations at all orders. Following on recent work on one dimensional flights, this study explores the one-point statistics and Minkowski functionals (density PDF, perimeter, Euler characteristic) of Rayleigh-Lévy flights in two dimensions. We derive the Euler characteristic in the mean field approximation and the density PDF and iso-field perimeter $W_{1}$ in beyond mean field calculations, and validate the results against simulations. The match is excellent throughout, even for fields with large variances, in particular when finite volume effects in the simulations are taken into account and when the calculation is extended beyond the mean field.

preprint2023arXiv

It takes two to know one: Computing accurate one-point PDF covariances from effective two-point PDF models

One-point probability distribution functions (PDFs) of the cosmic matter density are powerful cosmological probes that extract non-Gaussian properties of the matter distribution and complement two-point statistics. Computing the covariance of one-point PDFs is key for building a robust galaxy survey analysis for upcoming surveys like Euclid and the Rubin Observatory LSST and requires good models for the two-point PDFs characterising spatial correlations. In this work, we obtain accurate PDF covariances using effective shifted lognormal two-point PDF models for the mildly non-Gaussian weak lensing convergence and validate our predictions against large sets of Gaussian and non-Gaussian maps. We show how the dominant effects in the covariance matrix capturing super-sample covariance arise from a large-separation expansion of the two-point PDF and discuss differences between the covariances obtained from small patches and full sky maps. Finally, we describe how our formalism can be extended to characterise the PDF covariance for 3D-dimensional spectroscopic fields using the 3D matter PDF as an example. We describe how covariances from simulated boxes with fixed overall density can be supplemented with the missing super-sample covariance effect by relying on theoretical predictions validated against separate-universe style simulations.

preprint2022arXiv

Covariances of density probability distribution functions. Lessons from hierarchical models

Context: Statistical properties of the cosmic density fields are to a large extent encoded in the shape of the one-point density probability distribution functions (PDF). In order to successfully exploit such observables, a detailed functional form of the covariance matrix of the one-point PDF is needed. Aims: The objectives are to model the properties of this covariance for general stochastic density fields in a cosmological context. Methods: Leading and subleading contributions to the covariance were identified within a large class of models, the so-called hierarchical models. The validity of the proposed forms for the covariance matrix was assessed with the help of a toy model, the minimum tree model, for which a corpus of exact results could be obtained (forms of the one- and two-point PDF, large-scale density-bias functions, and full covariance matrix of the one-point PDF). Results: It is first shown that the covariance matrix elements are directly related to the spatial average of the two-point density PDF within the sample. The dominant contribution to this average is explicitly given for hierarchical models, which leads to the construction of specific density-bias functions. However, this contribution alone cannot be used to construct an operational likelihood function. Short distance effects are found to be have an important impact but are more difficult to derive as they depend more on the details of the model. However, a simple and generic form of these contributions is proposed. Detailed comparisons in the context of the Rayleigh-Levy flight model show that the large-scale effects capture the bulk of the supersample effects and that, by adding the short-distance contributions, a qualitatively correct model of the likelihood function can be obtained.

preprint2022arXiv

Numerical complexity of the joint nulled weak-lensing probability distribution function

In the context of tomographic cosmic shear surveys, there exists a nulling transformation of weak lensing observations (also called BNT transform) that allows us to simplify the correlation structure of tomographic cosmic shear observations, as well as to build observables that depend only on a localised range of redshifts and thus independent from the low-redshift/small-scale modes. This procedure renders possible accurate, and from-first-principles, predictions of the convergence and aperture mass one-point distributions (PDF). We here explore other consequences of this transformation on the (reduced) numerical complexity of the estimation of the joint PDF between nulled bins and demonstrate how to use these results to make theoretical predictions.

preprint2020arXiv

Game of cones: A nulling strategy for modelling lensing convergence in cones with large deviation theory

The distribution of the cosmic convergence field is modeled using a large-deviation principle where all non-Gaussian contributions are computed from first principles. The geometry of the past light cone is accounted for by constructing the total weak-lensing signal from contributions of the matter density in thin disk slices. The prediction of this model is successfully tested against numerical simulation with ray tracing, and found to be accurate within at least 5 per cent in the tails at redshift 1 and opening angle of 10 arcmin and even more so with increasing source redshift and opening angle. An accurate analytical approximation to the theory is also provided for practical implementation. The lensing kernel that mixes physical scales along the line-of-sight tends to reduce the domain of validity of this theoretical approach compared to the three dimensional case of cosmic densities in spherical cells. This effect is shown to be avoidable if a nulling procedure is implemented in order to localise the lensing line-of-sight integrations in a tomographic analysis. Accuracy in the tails is thus achieved within a percent for source redshifts between 0.5 and 1.5 and an opening angle of 10 arcmin. Applications to future weak-lensing surveys like Euclid and the specific issue of shape noise are discussed.

preprint2020arXiv

Observing Baryonic Acoustic Oscillations in tomographic cosmic shear surveys

We show that it is possible to build effective matter density power spectra in tomographic cosmic shear observations that exhibit the Baryonic Acoustic Oscillations (BAO) features once a nulling transformation has been applied to the data. The precision with which the amplitude and position of these features can be reconstructed is quantified in terms of sky coverage, intrinsic shape noise, median source redshift and number density of sources. BAO detection in Euclid or LSST like wide surveys will be possible with a modest signal-to-noise ratio. It would improve dramatically for slightly deeper surveys.

preprint2020arXiv

Post-Born corrections to the one-point statistics of (CMB) lensing convergence obtained via large deviation theory

Weak lensing of galaxies and CMB photons through the large-scale structure of the Universe is one of the most promising cosmological probes with upcoming experiments dedicated to its measurements such as Euclid/LSST and CMB Stage 4 experiments. With increasingly precise measurements, there is a dire need for accurate theoretical predictions. In this work, we focus on higher order statistics of the weak lensing convergence field, namely its cumulants such as skewness and kurtosis and its one-point probability distribution (PDF), and we quantify using perturbation theory the corrections coming from post-Born effects, meaning beyond the straight-line and independent lenses approximations. At first order, two such corrections arise: lens-lens couplings and geodesic deviation. Though the corrections are small for low source redshifts (below a few percents) and therefore for galaxy lensing, they become important at higher redshifts, notably in the context of CMB lensing, where the non-gaussianities computed from tree-order perturbation theory are found to be of the same order as the signal itself. We include these post-Born corrections on the skewness into a prediction for the one-point convergence PDF obtained with large deviation theory and successfully test these results against numerical simulations. The modelled PDF is indeed shown to perform better than the percent for apertures above ~ 10 arcminutes and typically in the three sigmas region around the mean.

preprint2020arXiv

Tomographic weak lensing bispectrum: a thorough analysis towards the next generation of galaxy surveys

We address key points for an efficient implementation of likelihood codes for modern weak lensing large-scale structure surveys. Specifically, we focus on the joint weak lensing convergence power spectrum-bispectrum probe and we tackle the numerical challenges required by a realistic analysis. Under the assumption of (multivariate) Gaussian likelihoods, we have developed a high performance code that allows highly parallelised prediction of the binned tomographic observables and of their joint non-Gaussian covariance matrix accounting for terms up to the 6-point correlation function and super-sample effects. This performance allows us to qualitatively address several interesting scientific questions. We find that the bispectrum provides an improvement in terms of signal-to-noise ratio (S/N) of about 10% on top of the power spectrum, making it a non-negligible source of information for future surveys. Furthermore, we are capable to test the impact of theoretical uncertainties in the halo model used to build our observables; with presently allowed variations we conclude that the impact is negligible on the S/N. Finally, we consider data compression possibilities to optimise future analyses of the weak lensing bispectrum. We find that, ignoring systematics, 5 equipopulated redshift bins are enough to recover the information content of a Euclid-like survey, with negligible improvement when increasing to 10 bins. We also explore principal component analysis and dependence on the triangle shapes as ways to reduce the numerical complexity of the problem.

preprint2016arXiv

Back in the saddle: Large-deviation statistics of the cosmic log-density field

We present a first principle approach to obtain analytical predictions for spherically-averaged cosmic densities in the mildly non-linear regime that go well beyond what is usually achieved by standard perturbation theory. A large deviation principle allows us to compute the leading-order cumulants of average densities in concentric cells. In this symmetry, the spherical collapse model leads to cumulant generating functions that are robust for finite variances and free of critical points when logarithmic density transformations are implemented. They yield in turn accurate density probability distribution functions (PDFs) from a straightforward saddle-point approximation valid for all density values. Based on this easy-to-implement modification, explicit analytic formulas for the evaluation of the one- and two-cell PDF are provided. The theoretical predictions obtained for the PDFs are accurate to a few percent compared to the numerical integration, regardless of the density under consideration and in excellent agreement with N-body simulations for a wide range of densities. This formalism should prove valuable for accurately probing the quasi-linear scales of low redshift surveys for arbitrary primordial power spectra.

preprint2016arXiv

Redshift-space distortions with wide angular separations

Redshift-space distortions are generally considered in the plane parallel limit, where the angular separation between the two sources can be neglected. Given that galaxy catalogues now cover large fractions of the sky, it becomes necessary to consider them in a formalism which takes into account the wide angle separations. In this article we derive an operational formula for the matter correlators in the Newtonian limit to be used in actual data sets, both in configuration and in Fourier spaces without relying on a plane-parallel approximation. We then recover the plane-parallel limit not only in configuration space where the geometry is simpler, but also in Fourier space, and we exhibit the first corrections that should be included in large surveys as a perturbative expansion over the plane-parallel results. We finally compare our results to existing literature, and show explicitly how they are related.

preprint2016arXiv

Response function of the large-scale structure of the universe to the small scale inhomogeneities

In order to infer the impact of the small-scale physics to the large-scale properties of the universe, we use a series of cosmological $N$-body simulations of self-gravitating matter inhomogeneities to measure, for the first time, the response function of such a system defined as a functional derivative of the nonlinear power spectrum with respect to its linear counterpart. Its measured shape and amplitude are found to be in good agreement with perturbation theory predictions except for the coupling from small to large-scale perturbations. The latter is found to be significantly damped, following a Lorentzian form. These results shed light on validity regime of perturbation theory calculations giving a useful guideline for regularization of small scale effects in analytical modeling. Most importantly our result indicates that the statistical properties of the large-scale structure of the universe are remarkably insensitive to the details of the small-scale physics, astrophysical or gravitational, paving the way for the derivation of robust estimates of theoretical uncertainties on the determination of cosmological parameters from large-scale survey observations.

preprint2016arXiv

The large-scale correlations of multi-cell densities and profiles, implications for cosmic variance estimates

In order to quantify the error budget in the measured probability distribution functions of cell densities, the two-point statistics of cosmic densities in concentric spheres is investigated. Bias functions are introduced as the ratio of their two-point correlation function to the two-point correlation of the underlying dark matter distribution. They describe how cell densities are spatially correlated. They are computed here via the so-called large deviation principle in the quasi-linear regime. Their large-separation limit is presented and successfully compared to simulations for density and density slopes: this regime is shown to be rapidly reached allowing to get sub-percent precision for a wide range of densities and variances. The corresponding asymptotic limit provides an estimate of the cosmic variance of standard concentric cell statistics applied to finite surveys. More generally, no assumption on the separation is required for some specific moments of the two-point statistics, for instance when predicting the generating function of cumulants containing any powers of concentric densities in one location and one power of density at some arbitrary distance from the rest. This exact "one external leg" cumulant generating function is used in particular to probe the rate of convergence of the large-separation approximation.

preprint2015arXiv

A Large Deviation Principle at play in Large-Scale Structure cosmology

We present an application of Large Deviation Theory to the problem of structure growth on large-scale structure cosmology. Starting from gaussian distributed overdensities on concentric spherical shells, we show that a Large Deviation Principle holds for the densities on the corresponding shells after gravitational evolution if no shell-crossing happens. As consequences of the Large Deviation Principle we obtain the cumulant generating function for the non-linear densities, and present formulae to compute the cumulant generating function for general window functions.

preprint2015arXiv

Consistent Modified Gravity Analysis of Anisotropic Galaxy Clustering Using BOSS DR11

We analyse the clustering of cosmic large scale structure using a consistent modified gravity perturbation theory, accounting for anisotropic effects along and transverse to the line of sight. The growth factor has a particular scale dependence in f(R) gravity and we fit for the shape parameter f_{R0} simultaneously with the distance and the large scale (general relativity) limit of the growth function. Using more than 690,000 galaxies in the Baryon Oscillation Spectroscopy Survey Data Release 11, we find no evidence for extra scale dependence, with the 95\% confidence upper limit |f_{R0}| <8 \times 10^{-4}. Future clustering data, such as from the Dark Energy Spectroscopic Instrument, can use this consistent methodology to impose tighter constraints.

preprint2015arXiv

On the importance of nonlinear couplings in large-scale neutrino streams

We propose a procedure to evaluate the impact of nonlinear couplings on the evolution of massive neutrino streams in the context of large-scale structure growth. Such streams can be described by general nonlinear conservation equations, derived from a multiple-flow perspective, which generalize the conservation equations of non-relativistic pressureless fluids. The relevance of the nonlinear couplings is quantified with the help of the eikonal approximation applied to the subhorizon limit of this system. It highlights the role played by the relative displacements of different cosmic streams and it specifies, for each flow, the spatial scales at which the growth of structure is affected by nonlinear couplings. We found that, at redshift zero, such couplings can be significant for wavenumbers as small as $k=0.2\,h$/Mpc for most of the neutrino streams.

preprint2015arXiv

The joint statistics of mildly non-linear cosmological densities and slopes in count-in-cells

In the context of count-in-cells statistics, the joint probability distribution of the density in two concentric spherical shells is predicted from first first principle for sigmas of the order of one. The agreement with simulation is found to be excellent. This statistics allows us to deduce the conditional one dimensional probability distribution function of the slope within under dense (resp. overdense) regions, or of the density for positive or negative slopes. The former conditional distribution is likely to be more robust in constraining the cosmological parameters as the underlying dynamics is less evolved in such regions. A fiducial dark energy experiment is implemented on such counts derived from Lambda-CDM simulations.

preprint2014arXiv

Cosmic shear full nulling: sorting out dynamics, geometry and systematics

An explicit full nulling scheme for cosmic shear observations is presented. It makes possible the construction of shear maps from extended source distributions for which the lens distance distribution is restricted to a definite interval. Such a construction allows to build totally independent shear maps, at all scales, that can be taken advantage of to constrain background cosmological parameters and systematics using the full statistical power of cosmic shear observations. Another advantage of such construction is that, as the lens redshift distribution can be made arbitrarily narrow, scale mixing due to projection effects can be limited allowing controlled predictions on the large scale shear power spectrum from perturbation theory calculations.

preprint2014arXiv

Cosmological Perturbation Theory for streams of relativistic particles

Motion equations describing streams of relativistic particles and their properties are explored in detail in the framework of Cosmological Perturbation Theory. Those equations, derived in any metric both in the linear and nonlinear regimes, express the matter and momentum conservation. In this context we extend the setup of adiabatic initial conditions - that was initially performed in the Conformal Newtonian gauge - to the Synchronous gauge. The subhorizon limit of the nonlinear motion equations written in a generic perturbed Friedmann-Lemaître metric is then derived and analyzed. We show in particular that the momentum field $P_{i}(x)$ is always potential in the linear regime and remains so at subhorizon scales in the nonlinear regime. Finally the equivalence principle is exploited to highlight invariance properties satisfied by such a system of equations, extending that known for streams of non-relativistic particles, namely the extended Galilean invariance.

preprint2014arXiv

Density fields and halo mass functions in the Geometrical Adhesion toy Model

In dimension 2 and above, the Burgers dynamics, the so-called "adhesion model" in cosmology, can actually give rise to several dynamics in the inviscid limit. We investigate here the statistical properties of the density field when it is defined by a "geometrical model" associated with this Burgers velocity field and where the matter distribution is fully determined, at each time step, by geometrical constructions. Our investigations are based on a set of numerical experiments that make use of an improved algorithm, for which the geometrical constructions are efficient and robust. In this work we focus on Gaussian initial conditions with power-law power spectra of slope $n$ in the range $-3<n<1$, where a self-similar evolution develops, and we compute the behavior of power spectra, density probability distributions and mass functions. As expected for such dynamics, the density power spectra show universal high-$k$ tails that are governed by the formation of pointlike masses. The two other statistical indicators however show the same qualitative properties as those observed for 3D gravitational clustering. In particular, the mass functions obey a Press-Schechter like scaling up to a very good accuracy in 1D, and to a lesser extent in 2D. Our results suggest that the "geometrical adhesion model", whose solution is fully known at all times, provides a precious tool to understand some of the statistical constructions frequently used to study the development of mass halos in gravitational clustering.

preprint2014arXiv

Non-Gaussian Minkowski functionals & extrema counts in redshift space

In the context of upcoming large-scale structure surveys such as Euclid, it is of prime importance to quantify the effect of peculiar velocities on geometric probes. Hence the formalism to compute in redshift space the geometrical and topological one-point statistics of mildly non-Gaussian 2D and 3D cosmic fields is developed. Leveraging the partial isotropy of the target statistics, the Gram-Charlier expansion of the joint probability distribution of the field and its derivatives is reformulated in terms of the corresponding anisotropic variables. In particular, the cosmic non-linear evolution of the Minkowski functionals, together with the statistics of extrema are investigated in turn for 3D catalogues and 2D slabs. The amplitude of the non-Gaussian redshift distortion correction is estimated for these geometric probes. In 3D, gravitational perturbation theory is implemented in redshift space to predict the cosmic evolution of all relevant Gram-Charlier coefficients. Applications to the estimation of the cosmic parameters sigma(z) and beta=f/b1 from upcoming surveys is discussed. Such statistics are of interest for anisotropic fields beyond cosmology.

preprint2014arXiv

Regularized cosmological power spectrum and correlation function in modified gravity models

Based on the multi-point propagator expansion, we present resummed perturbative calculations for cosmological power spectra and correlation functions in the context of modified gravity. In a wide class of modified gravity models that have a screening mechanism to recover general relativity (GR) on small scales, we apply the eikonal approximation to derive the governing equation for resummed propagator that partly includes the non-perturbative effect in the high-k limit. The resultant propagator in the high-k limit contains the new corrections arising from the screening mechanism as well as the standard exponential damping. We explicitly derive the expression for new high-k contributions in specific modified gravity models, and find that in the case of f(R) gravity for a currently constrained model parameter, the corrections are basically of the sub-leading order and can be neglected. Thus, in f(R) gravity, similarly to the GR case, we can analytically construct the regularized propagator that reproduces both the resummed high-k behavior and the low-k results computed with standard perturbation theory, consistently taking account of the nonlinear modification of gravity valid at large scales. With the regularized multi-point propagators, we give predictions for power spectrum and correlation function at one-loop order, and compare those with N-body simulations in f(R) gravity model. As an important application, we also discuss the redshift-space distortions and compute the anisotropic power spectra and correlation functions.

preprint2013arXiv

Describing massive neutrinos in cosmology as a collection of independent flows

A new analytical approach allowing to account for massive neutrinos in the non-linear description of the growth of the large-scale structure of the universe is proposed. Unlike the standard approach in which neutrinos are described as a unique hot fluid, it is shown that the overall neutrino fluid can be equivalently decomposed into a collection of independent flows. Starting either from elementary conservation equations or from the evolution equation of the phase-space distribution function, we derive the two non-linear motion equations that each of these flows satisfies. Those fluid equations describe the evolution of macroscopic fields. We explain in detail the connection between the collection of flows we defined and the standard massive neutrino fluid. Then, in the particular case of adiabatic initial conditions, we explicitly check that, at linear order, the resolution of this new system of equations reproduces the results obtained in the standard approach based on the collisionless Boltzmann hierarchy. Besides, the approach advocated in this paper allows to show how each neutrino flow settles into the cold dark matter flow depending on initial velocities. It opens the way to a fully non-linear treatment of the dynamical evolution of neutrinos in the framework of large-scale structure growth.

preprint2013arXiv

Power spectra in the eikonal approximation with adiabatic and non-adiabatic modes

We use the so-called eikonal approximation, recently introduced in the context of cosmological perturbation theory, to compute power spectra for multi-component fluids. We demonstrate that, at any given order in standard perturbation theory, multipoint power spectra do not depend on the large-scale adiabatic modes. Moreover, we employ perturbation theories to decipher how nonadiabatic modes, such as a relative velocity between two different components, damp the small-scale matter power spectrum, a mechanism recently described in the literature. In particular, we do an explicit calculation at 1-loop order of this effect. While the 1-loop result eventually breaks down, we show how the damping effect can be fully captured by the help of the eikonal approximation. A relative velocity not only induces mode damping but also creates large-scale anisotropic modulations of the matter power spectrum amplitude. We illustrate this for the Local Group environment.

preprint2013arXiv

Precision modeling of redshift-space distortions from multi-point propagator expansion

Using a full implementation of resummed perturbation theory (PT) from a multi-point propagator expansion, we put forward new theoretical predictions for the two-point statistics of matter fluctuations in redshift space. The predictions consistently include PT corrections up to the two-loop order and are based on an improved prescription of the redshift-space distortions (RSD) that properly takes into account their non-Gaussian impact from a systematic low-k expansion. In contrast to the previous studies that partly used standard PT calculations, the present treatment is able to provide a consistent prediction for both power spectra and correlation functions. These results are compared with $N$-body simulations with which a very good agreement is found up to the quadrupole moment. The theoretical predictions for the hexadecapole moment of the power spectra are however found to significantly departs from the numerical results at low redshift. We examine this issue and found it to be likely related to an improper modeling of the RSD damping effects on which this moment shows large dependence.

preprint2013arXiv

Resummed propagators in multi-component cosmic fluids with the eikonal approximation

We introduce the eikonal approximation to study the effect of the large-scale motion of cosmic fluids on their small-scale evolution. This approach consists in collecting the impact of the long-wavelength displacement field into a single or finite number of random variables, whose statistical properties can be computed from the initial conditions. For a single dark matter fluid, we show that we can recover the nonlinear propagators of renormalized perturbation theory. These are obtained with no need to assume that the displacement field follows the linear theory. Then we extend the eikonal approximation to many fluids. In particular, we study the case of two non-relativistic components and we derive their resummed propagators in the presence of isodensity modes. Unlike the adiabatic case, where only the phase of small-scale modes is affected by the large-scale advection field, the isodensity modes change also the amplitude on small scales. We explicitly solve the case of cold dark matter-baryon mixing and find that the isodensity modes induce only very small corrections to the resummed propagators.

preprint2013arXiv

Statistics of cosmic density profiles from perturbation theory

The joint probability distribution function (PDF) of the density within multiple concentric spherical cells is considered. It is shown how its cumulant generating function can be obtained at tree order in perturbation theory as the Legendre transform of a function directly built in terms of the initial moments. In the context of the upcoming generation of large-scale structure surveys, it is conjectured that this result correctly models such a function for finite values of the variance. Detailed consequences of this assumption are explored. In particular the corresponding one-cell density probability distribution at finite variance is computed for realistic power spectra, taking into account its scale variation. It is found to be in agreement with $Λ$-CDM simulations at the few percent level for a wide range of density values and parameters. Related explicit analytic expansions at the low and high density tails are given. The conditional (at fixed density) and marginal probability of the slope -- the density difference between adjacent cells -- and its fluctuations is also computed from the two-cells joint PDF; it also compares very well to simulations, in particular in under-dense regions, with a significant reduced cosmic scatter compared to over-dense regions. It is emphasized that this could prove useful when studying the statistical properties of voids as it can serve as a statistical indicator to test gravity models and/or probe key cosmological parameters.

preprint2013arXiv

The evolution of the large-scale structure of the universe: beyond the linear regime

These lecture notes introduce analytical tools, methods and results describing the growth of cosmological structure beyond the linear regime. The presentation is focused on the single flow regime of the Vlasov-Poisson equation describing the development of gravitational instabilities in a pressureless fluid. These notes include the introduction of diagrammatic representations of the standard perturbation theory with applications to the calculation of the so-called loop contributions to the power spectra. A large part of these notes is devoted to the exploration of the convergence properties of these terms from the contribution of both the long-wave modes and the short-wave modes. The resulting performances of the two-loop corrections of the power spectra are then presented. Finally other avenues that use different methods are explored. In particular it is shown how joint density and profile probability distribution functions can be constructed out of the multiple-variable cumulant generating functions computed at tree order.

preprint2012arXiv

Constructing Regularized Cosmic Propagators

We present a new scheme for the general computation of cosmic propagators that allow to interpolate between standard perturbative results at low-k and their expected large-k resummed behavior. This scheme is applicable to any multi-point propagator and allows the matching of perturbative low-k calculations to any number of loops to their large-k behavior, and can potentially be applied in case of non-standard cosmological scenarios such as those with non-Gaussian initial conditions. The validity of our proposal is checked against previous prescriptions and measurements in numerical simulations showing a remarkably good agreement. Such a generic prescription for multi-point propagators provides the necessary building blocks for the computation of polyspectra in the context of the so-called Gamma-expansion introduced by Bernardeau et al. (2008). As a concrete application we present a consistent calculation of the matter bispectrum at one-loop order.

preprint2012arXiv

Cosmic propagators at two-loop order

We explore the properties of two-point cosmic propagators when Perturbation Theory (PT) loop corrections are consistently taken into account. We show in particular how the interpolation scheme proposed in arXiv:1112.3895 can be explicitly used up to two-loop order introducing the notion of regular parts for the contributing terms. Extending the one-loop results, we then derive and give semi analytical forms of the two-loop contributions for both the cosmic density and velocity propagators. These results are tested against numerical simulations and shown to significantly improve upon one-loop results for redshifts above 0.5. We found however that at lower redshift two-loop order corrections are too large partly due to a strong sensitivity of those terms to the small scale modes. We show that this dependence is expected to be even larger for higher order loop corrections both from theoretical investigations and numerical tests, the latter obtained with Monte Carlo evaluations of the three-loop contributions. This makes small-scale regularization schemes necessary for the use of such higher order corrections.

preprint2012arXiv

MPTbreeze: A fast renormalized perturbative scheme

We put forward and test a simple description of multi-point propagators (MP), which serve as building-blocks to calculate the nonlinear matter power spectrum. On large scales these propagators reduce to the well-known kernels in standard perturbation theory, while at smaller scales they are suppresed due to nonlinear couplings. Through extensive testing with numerical simulations we find that this decay is characterized by the same damping scale for both two and three-point propagators. In turn this transition can be well modeled with resummation results that exponentiate one-loop computations. For the first time, we measure the four components of the non-linear (two-point) propagator using dedicated simulations started from two independent random Gaussian fields for positions and velocities, verifying in detail the fundamentals of propagator resummation. We use these results to develop an implementation of the MP-expansion for the nonlinear power spectrum that only requires seconds to evaluate at BAO scales. To test it we construct six suites of large numerical simulations with different cosmologies. From these and LasDamas runs we show that the nonlinear power spectrum can be described at the ~ 2% level at BAO scales for redshifts in the range [0-2.5]. We make a public release of the MPTbreeze code with the hope that it can be useful to the community.

preprint2012arXiv

RegPT: Direct and fast calculation of regularized cosmological power spectrum at two-loop order

We present a specific prescription for the calculation of cosmological power spectra, exploited here at two-loop order in perturbation theory (PT), based on the multi-point propagator expansion. In this approach power spectra are constructed from the regularized expressions of the propagators that reproduce both the resummed behavior in the high-k limit and the standard PT results at low-k. With the help of N-body simulations, we show that such a construction gives robust and accurate predictions for both the density power spectrum and the correlation function at percent-level in the weakly non-linear regime. We then present an algorithm that allows accelerated evaluations of all the required diagrams by reducing the computational tasks to one-dimensional integrals. This is achieved by means of pre-computed kernel sets defined for appropriately chosen fiducial models. The computational time for two-loop results is then reduced from a few minutes, with the direct method, to a few seconds with the fast one. The robustness and applicability of this method are tested against the power spectrum cosmic emulator from which a wide variety of cosmological models can be explored. The fortran program with which direct and fast calculations of power spectra can be done, RegPT, is publicly released as part of this paper.

preprint2011arXiv

Cosmological Large-scale Structures beyond Linear Theory in Modified Gravity

We consider the effect of modified gravity on the growth of large-scale structures at second order in perturbation theory. We show that modified gravity models changing the linear growth rate of fluctuations are also bound to change, although mildly, the mode coupling amplitude in the density and reduced velocity fields. We present explicit formulae which describe this effect. We then focus on models of modified gravity involving a scalar field coupled to matter, in particular chameleons and dilatons, where it is shown that there exists a transition scale around which the existence of an extra scalar degree of freedom induces significant changes in the coupling properties of the cosmic fields. We obtain the amplitude of this effect for realistic dilaton models at the tree-order level for the bispectrum, finding them to be comparable in amplitude to those obtained in the DGP model.

preprint2011arXiv

Exact Results for Propagators in the Geometrical Adhesion Model

The Geometrical Adhesion Model (GAM) we described in previous papers provides a fully solved model for the nonlinear evolution of fields that mimic the cosmological evolution of pressureless fluids. In this context we explore the expected late time properties of the cosmic propagators once halos have formed, in a regime beyond the domain of application of perturbation theories. Whereas propagators in Eulerian coordinates are closely related to the velocity field we show here that propagators defined in Lagrangian coordinates are intimately related to the halo mass function. Exact results can be obtained in the 1D case. In higher dimensions, the computations are more intricate because of to the dependence of the propagators on the detailed shape of the underlying Lagrangian-space tessellations, that is, on the geometry of the regions that eventually collapse to form halos. We illustrate these results for both the 1D and the 2D dynamics. In particular we give here the expected asymptotic behaviors obtained for power-law initial power spectra. These analytical results are compared with the results obtained with dedicated numerical simulations.

preprint2011arXiv

The y-sky: diffuse spectral distortions of the cosmic microwave background

The non-linear evolution of the energy density of the radiation is shown to induce spectral distortions of the cosmic microwave background both at recombination and during the reionization era. This distortion has the same spectral signature as the one produced by the re-scattering of photons by non-relativistic hot electrons, the thermal Sunyaev-Zeldovich effect, whose amplitude is quantified by a Compton y parameter. A diffuse y-sky is then expected to emerge from mode couplings in the non-linear evolution of the cosmological perturbations and to superimpose to the point source contributions of galaxy clusters. The equations describing the evolution of the y field and a hierarchy governing its angular multipoles are derived from the second order Boltzmann equation. These equations are then integrated numerically to obtain the first predicted power spectrum of the diffuse y-sky. It is found to be a remarkable tracer of the reionization history of the Universe.

preprint2010arXiv

CMB spectra and bispectra calculations: making the flat-sky approximation rigorous

This article constructs flat-sky approximations in a controlled way in the context of the cosmic microwave background observations for the computation of both spectra and bispectra. For angular spectra, it is explicitly shown that there exists a whole family of flat-sky approximations of similar accuracy for which the expression and amplitude of next to leading order terms can be explicitly computed. It is noted that in this context two limiting cases can be encountered for which the expressions can be further simplified. They correspond to cases where either the sources are localized in a narrow region (thin-shell approximation) or are slowly varying over a large distance (which leads to the so-called Limber approximation). Applying this to the calculation of the spectra it is shown that, as long as the late integrated Sachs-Wolfe contribution is neglected, the flat-sky approximation at leading order is accurate at 1% level for any multipole. Generalization of this construction scheme to the bispectra led to the introduction of an alternative description of the bispectra for which the flat-sky approximation is well controlled. This is not the case for the usual description of the bispectrum in terms of reduced bispectrum for which a flat-sky approximation is proposed but the next-to-leading order terms of which remain obscure.

preprint2010arXiv

Full-sky lensing shear at second order

We compute the reduced cosmic shear up to second order in the gravitational potential without relying on the small angle or thin-lens approximation. This is obtained by solving the Sachs equation which describes the deformation of the infinitesimal cross-section of light bundle in the optical limit, and maps galaxy intrinsic shapes into their angular images. The calculation is done in the Poisson gauge without a specific matter content, including vector and tensor perturbations generated at second order and taking account of the inhomogeneities of a fixed redshift source plane. Our final result is expressed in terms of spin-2 operators on the sphere and is valid on the full sky. Beside the well known lens-lens and Born corrections that dominate on small angular scales, we find new non-linear couplings. These are a purely general relativistic intrinsic contribution, a coupling between the gravitational potential at the source with the lens, couplings between the time delay with the lens, couplings between two photon deflections, as well as non-linear couplings due to the second-order vector and tensor components. The inhomogeneity in the redshift of the source induces a coupling between the photon redshift with the lens. All these corrections become important on large angular scales and should thus be included when computing higher-order observables such as the bispectrum, in full or partially full-sky surveys.

preprint2010arXiv

Gravity and non-gravity mediated couplings in multiple-field inflation

Mechanisms for the generation of primordial non-Gaussian metric fluctuations in the context of multiple-field inflation are reviewed. As long as kinetic terms remain canonical, it appears that nonlinear couplings inducing non-gaussianities can be split into two types. The extension of the one-field results to multiple degrees of freedom leads to gravity mediated couplings that are ubiquitous but generally modest. Multiple-field inflation offers however the possibility of generating non-gravity mediated coupling in isocurvature directions that can eventually induce large non-Gaussianities in the metric fluctuations. The robustness of the predictions of such models is eventually examined in view of a case study derived from a high-energy physics construction.

preprint2010arXiv

Merging and fragmentation in the Burgers dynamics

We explore the noiseless Burgers dynamics in the inviscid limit, the so-called ``adhesion model'' in cosmology, in a regime where (almost) all the fluid particles are embedded within point-like massive halos. Following previous works, we focus our investigations on a ``geometrical'' model, where the matter evolution within the shock manifold is defined from a geometrical construction. This hypothesis is at variance with the assumption that the usual continuity equation holds but, in the inviscid limit, both models agree in the regular regions. Taking advantage of the formulation of the dynamics of this ``geometrical model'' in terms of Legendre transforms and convex hulls, we study the evolution with time of the distribution of matter and the associated partitions of the Lagrangian and Eulerian spaces. We describe how the halo mass distribution derives from a triangulation in Lagrangian space, while the dual Voronoi-like tessellation in Eulerian space gives the boundaries of empty regions with shock nodes at their vertices. We then emphasize that this dynamics actually leads to halo fragmentations for space dimensions greater or equal to 2 (for the inviscid limit studied in this article). This is most easily seen from the properties of the Lagrangian-space triangulation and we illustrate this process in the two-dimensional (2D) case. In particular, we explain how point-like halos only merge through three-body collisions while two-body collisions always give rise to two new massive shock nodes (in 2D). This generalizes to higher dimensions and we briefly illustrate the three-dimensional (3D) case. This leads to a specific picture for the continuous formation of massive halos through successive halo fragmentations and mergings.

preprint2010arXiv

Mode coupling evolution in arbitrary inflationary backgrounds

The evolution of high order correlation functions of a test scalar field in arbitrary inflationary backgrounds is computed. Whenever possible, exact results are derived from quantum field theory calculations. Taking advantage of the fact that such calculations can be mapped, for super-horizon scales, into those of a classical system, we express the expected correlation functions in terms of classical quantities, power spectra, Green functions, that can be easily computed in the long-wavelength limit. Explicit results are presented that extend those already known for a de Sitter background. In particular the expressions of the late time amplitude of bispectrum and trispectrum, as well as the whole high-order correlation structure, are given in terms of the expansion factor behavior. When compared to the case of a de Sitter background, power law inflation and chaotic inflation induced by a massive field are found to induce high order correlation functions the amplitudes of which are amplified by almost one order of magnitude. These results indicate that the dependence of the related non-Gaussian parameters - such as f_NL - on the wave-modes is at percent level.

preprint2010arXiv

Multi-Point Propagators for Non-Gaussian Initial Conditions

We show here how Renormalized Perturbation Theory (RPT) calculations applied to the quasi-linear growth of the large-scale structure can be carried on in presence of primordial non-Gaussian (PNG) initial conditions. It is explicitly demonstrated that the series reordering scheme proposed in Bernardeau, Crocce and Scoccimarro (2008) is preserved for non-Gaussian initial conditions. This scheme applies to the power spectrum and higher order spectra and is based on a reorganization of the contributing terms into sum of products of multi-point propagators. In case of PNG new contributing terms appear, the importance of which is discussed in the context of current PNG models. The properties of the building blocks of such resummation schemes, the multi-point propagators, are then investigated. It is first remarked that their expressions are left unchanged at one-loop order irrespectively of statistical properties of the initial field. We furthermore show that the high-momemtum limit of each of these propagators can be explicitly computed even for arbitrary initial conditions. They are found to be damped by an exponential cutoff whose expression is directly related to the moment generating function of the one-dimensional displacement field. This extends what had been established for multi-point propagators for Gaussian initial conditions. Numerical forms of the cut-off are shown for the so-called local model of PNG.

preprint2010arXiv

The cosmic microwave background bispectrum from the non-linear evolution of the cosmological perturbations

This article presents the first computation of the complete bispectrum of the cosmic microwave background temperature anisotropies arising from the evolution of all cosmic fluids up to second order, including neutrinos. Gravitational couplings, electron density fluctuations and the second order Boltzmann equation are fully taken into account. Comparison to limiting cases that appeared previously in the literature are provided. These are regimes for which analytical insights can be given. The final results are expressed in terms of equivalent fNL for different configurations. It is found that for moments up to lmax=2000, the signal generated by non-linear effects is equivalent to fNL~5 for both local-type and equilateral-type primordial non-Gaussianity.

preprint2009arXiv

Eulerian and Lagrangian propagators for the adhesion model (Burgers dynamics)

Motivated by theoretical studies of gravitational clustering in the Universe, we compute propagators (response functions) in the adhesion model. This model, which is able to reproduce the skeleton of the cosmic web and includes nonlinear effects in both Eulerian and Lagrangian frameworks, also corresponds to the Burgers equation of hydrodynamics. Focusing on the one-dimensional case with power-law initial conditions, we obtain exact results for Eulerian and Lagrangian propagators. We find that Eulerian propagators can be expressed in terms of the one-point velocity probability distribution and show a strong decay at late times and high wavenumbers, interpreted as a "sweeping effect" but not a genuine damping of small-scale structures. By contrast, Lagrangian propagators can be written in terms of the shock mass function -- which would correspond to the halo mass function in cosmology -- and saturate to a constant value at late times. Moreover, they show a power-law dependence on scale or wavenumber which depends on the initial power-spectrum index and is directly related to the low-mass tail of the shock mass function. These results strongly suggest that Lagrangian propagators are much more sensitive probes of nonlinear structures in the underlying density field and of relaxation processes than their Eulerian counterparts.