Source author record

Patrick Valageas

Patrick Valageas 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

40works
6topics
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

40 published item(s)

preprint2021arXiv

The XXL survey: XLVI. Forward cosmological analysis of the C1 cluster sample

We present the forward cosmological analysis of an $XMM$ selected sample of galaxy clusters out to a redshift of unity. Following our previous 2018 study based on the dn/dz quantity alone, we perform an upgraded cosmological analysis of the same XXL C1 cluster catalogue (178 objects), with a detailed account of the systematic errors. We follow the ASpiX methodology: the distribution of the observed X-ray properties of the cluster population is analysed in a 3D observable space (count rate, hardness ratio, redshift) and modelled as a function of cosmology. Compared to more traditional methods, ASpiX allows the inclusion of clusters down to a few tens of photons. We obtain an improvement by a factor of 2 compared to the previous analysis by letting the normalisation of the M-T relation and the evolution of the L-T relation free. Adding constraints from the XXL cluster 2-point correlation function and the BAO from various surveys decreases the uncertainties by 23 and 53 % respectively, and 62% when adding both. Switching to the scaling relations from the Subaru analysis, and letting free more parameters, our final constraints are $σ_8$ = $0.99^{+0.14}_{-0.23}$, $Ω_m$ = 0.296 $\pm$ 0.034 ($S_8 = 0.98^{+0.11}_{-0.21}$) for the XXL sample alone. Finally, we combine XXL ASpiX, the XXL cluster 2-point correlation function and the BAO, with 11 free parameters, allowing for the cosmological dependence of the scaling relations in the fit. We find $σ_8$ = $0.793^{+0.063}_{-0.12}$, $Ω_m$ = 0.364 $\pm$ 0.015 ($S_8 = 0.872^{+0.068}_{-0.12}$), but still compatible with Planck CMB at 2.2$σ$. The results obtained by the ASpiX method are promising; further improvement is expected from the final XXL cosmological analysis involving a cluster sample twice as large. Such a study paves the way for the analysis of the eROSITA and future Athena surveys.

preprint2020arXiv

Goldstone models of modified gravity

We investigate scalar-tensor theories where matter couples to the scalar field via a kinetically dependent conformal coupling. These models can be seen as the low-energy description of invariant field theories under a global Abelian symmetry. The scalar field is then identified with the Goldstone mode of the broken symmetry. It turns out that the properties of these models are very similar to the ones of ultralocal theories where the scalar-field value is directly determined by the local matter density. This leads to a complete screening of the fifth force in the Solar System and between compact objects, through the ultralocal screening mechanism. On the other hand, the fifth force can have large effects in extended structures with large-scale density gradients, such as galactic halos. Interestingly, it can either amplify or damp Newtonian gravity, depending on the model parameters. We also study the background cosmology and the linear cosmological perturbations. The background cosmology is hardly different from its $Λ$-CDM counterpart whilst cosmological perturbations crucially depend on whether the coupling function is convex or concave. For concave functions, growth is hindered by the repulsiveness of the fifth force whilst it is enhanced in the convex case. In both cases, the departures from the $Λ$-CDM cosmology increase on smaller scales and peak for galactic structures. For concave functions, the formation of structure is largely altered below some characteristic mass, as smaller structures are delayed and would form later through fragmentation, as in some warm dark matter scenarios. For convex models, small structures form more easily than in the $Λ$-CDM scenario.

preprint2020arXiv

K-essence scalar dark matter solitons around supermassive black holes

We consider scalar dark matter models where the theory has a shift symmetry only broken by the scalar mass term. We restrict ourselves to K-essence kinetic terms where the shift symmetric part of the Lagrangian is a function of the first derivatives of the scalar field only. In the low-amplitude and nonrelativistic regime, which applies on large galactic scales, scalar clouds form solitons with a finite core. Close to the center of galaxies, where a supermassive Black Hole (BH) resides, we analyze the scalar field distribution and the fate of the dark matter soliton when subject to the BH gravitational attraction. We show that the scalar field profile around such a central BH can be described by new oscillatory solutions of a modified Klein-Gordon equation, which generalize the harmonic oscillations of free scalar dark matter in a flat environment and the Jacobi elliptic functions of the $ϕ^4$ model. Moreover, we find that, depending on the form of the K-essence kinetic term, regular solutions can be constructed or not, which connect the relativistic ingoing wavelike profile of the scalar field at the BH horizon to the nearly static nonrelativistic soliton at large distance. These profiles have a constant flux and represent the slow infall of scalar matter into the BH. We show that this regular behavior is only possible for K-essence functions that satisfy the usual conditions for the absence of ghosts and gradient instabilities, together with a new restriction on the growth of the kinetic function $K(X)$ for large argument. It turns out that the same conditions of stability guarantee that quantum corrections are tamed, provided that the mass of the scalar field is less than $10^{-3}$ eV and the strong coupling scale of the model $Λ$ is much larger than the scalar mass.

preprint2020arXiv

Lagrangian-space Gaussian ansatz for the matter redshift-space power spectrum and correlation function

We study the predictions for the matter redshift-space power spectrum and correlation function of a Lagrangian-space Gaussian ansatz introduced in a previous work. This model is a natural extension of the Zeldovich approximation, where the displacement and velocity power spectra are determined by the equations of motion, instead of being set equal to the linear power spectrum. It does not contain any free parameter. As for the real-space statistics, we find that this Lagrangian-space approach is much more efficient for the correlation functions than for the power spectra. The damping of the BAO oscillations is well recovered but there is a large smooth drift from the simulations in the power spectra. The multipoles of the correlation functions are well recovered on BAO scales, with an accuracy of $2\%$ for $ξ^s_0$ down to $10 h^{-1}$ Mpc, and of $3\%$ for $ξ^s_2$ down to $26 h^{-1}$ Mpc, at $z \geq 0.35$.

preprint2020arXiv

Projection of the gravitational dynamics on a subspace of probability distributions: curl-free Gaussian ansatz

We present a new approach to model the gravitational dynamics of large-scale structures. Instead of solving the equations of motion up to a finite perturbative order or building phenomenological models, we follow the evolution of the probability distribution of the displacement and velocity fields within an approximation subspace. Keeping the exact equations of motion with their full nonlinearity, this provides a nonperturbative scheme that goes beyond shell crossing. Focusing on the simplest case of a curl-free Gaussian ansatz for the displacement and velocity fields, we find that truncations of the power spectra on nonlinear scales directly arise from the equations of motion. This leads to a truncated Zeldovich approximation for the density power spectrum, but with a truncation that is not set a priori and with different power spectra for the displacement and velocity fields. The positivity of their auto power spectra also follows from the equations of motion. Although the density power spectrum is only recovered up to a smooth drift on BAO scales, the predicted density correlation function agrees with numerical simulations within $2\%$ from BAO scales down to $7 h^{-1} {\rm Mpc}$ at $z \geq 0.35$, without any free parameter.

preprint2019arXiv

Fate of scalar dark matter solitons around supermassive galactic black holes

In scalar-field dark matter scenarios, a scalar-field soliton could form at the center of galactic halos, around the supermassive black holes that sit at the center of galaxies. Focusing on the large scalar-mass limit, where the soliton is formed by the balance between self-gravity and a repulsive self-interaction, we study the infall of the scalar field onto the central Schwarzschild black hole. We derive the scalar-field profile, from the Schwarzschild radius to the large radii dominated by the scalar cloud. We show that the steady state solution selects the maximum allowed flux, with a critical profile that is similar to the transonic solution obtained for the hydrodynamic case. This finite flux, which scales as the inverse of the self-interaction coupling, is small enough to allow the dark matter soliton to survive for many Hubble times.

preprint2016arXiv

3D Weak Lensing: Modified Theories of Gravity

Weak lensing (WL) promises to be a particularly sensitive probe of both the growth of large scale structure (LSS) as well as the fundamental relation between matter density perturbations and metric perturbations, thus providing a powerful tool with which we may constrain modified theories of gravity (MG) on cosmological scales. Future deep, wide-field WL surveys will provide an unprecedented opportunity to constrain deviations from General Relativity (GR). Employing a three-dimensional (3D) analysis based on the spherical Fourier-Bessel (sFB) expansion, we investigate the extent to which MG theories will be constrained by a typical 3D WL survey configuration including noise from the intrinsic ellipticity distribution $σ_ε$ of source galaxies. Here we focus on two classes of screened theories of gravity: i) $f(R)$ chameleon models and ii) environmentally dependent dilaton models. We use one-loop perturbation theory combined with halo models in order to accurately model the evolution of matter power-spectrum with redshift in these theories. Using a Fisher information matrix based approach, we show that for an all-sky spectroscopic survey, the parameter $f_{R_0}$ can be constrained in the range $f_{R_0}< 5\times 10^{-6}(9\times 10^{-6})$ for $n=1(2)$ with a 3$σ$ confidence level. This can be achieved by using relatively low order angular harmonics $\ell<100$. Including higher order harmonics $\ell>100$ can further tighten the constraints, making them comparable to current solar-system constraints. We also employ a Principal Component Analysis (PCA) in order to study the parameter degeneracies in the MG parameters. Our results can trivially be extended to other MG theories, such as the K-mouflage models. The confusion from intrinsic ellipticity correlation and modification of the matter power-spectrum at small scale due to feedback mechanisms is briefly discussed.

preprint2016arXiv

Supersymmetric chameleons and ultra-local models

Super-chameleon models where all types of matter belong to three secluded sectors, i.e. the dark, supersymmetry breaking and matter sectors, are shown to be dynamically equivalent to ultra-local models of modified gravity. In the dark sector, comprising both dark matter and dark energy, the interaction range between the dark energy field and dark matter is constrained to be extremely short, i.e. shorter than the inverse gravitino mass set by supersymmetry breaking. This realises an extreme version of chameleon screening of the dark energy interaction. On the other hand, the baryonic matter sector decouples from the dark energy in a Damour-Polyakov way. These two mechanisms preclude the existence of any modification of gravity locally in the Solar System due to the presence of the super-chameleon field. On larger scales, the super-chameleon can have effects on the growth of structure and the number of dark matter halos. It can also affect the dynamics of galaxies where the fifth force interaction that it induces can have the same order of magnitude as Newton's interaction.

preprint2016arXiv

Ultra-local models of modified gravity without kinetic term

We present a class of modified-gravity theories which we call ultra-local models. We add a scalar field, with negligible kinetic terms, to the Einstein-Hilbert action. We also introduce a conformal coupling to matter. This gives rise to a new screening mechanism which is not entirely due to the non-linearity of the scalar field potential or the coupling function but to the absence of the kinetic term. As a result this removes any fifth force between isolated objects in vacuum. The predictions of these models only depend on a single free function, as the potential and the coupling function are degenerate, with an amplitude given by a parameter $α\lesssim 10^{-6}$, whose magnitude springs from requiring a small modification of Newton's potential astrophysically and cosmologically. This singles out a redshift $z_α \sim α^{-1/3} \gtrsim 100$ where the fifth force is the greatest. The cosmological background follows the $Λ$-CDM history within a $10^{-6}$ accuracy, while cosmological perturbations are significantly enhanced (or damped) on small scales, $k \gtrsim 2 h {\rm Mpc}^{-1}$ at $z=0$. The spherical collapse and the halo mass function are modified in the same manner. We find that the modifications of gravity are greater for galactic or sub-galactic structures. We also present a thermodynamic analysis of the non-linear and inhomogeneous fifth-force regime where we find that the Universe is not made more inhomogeneous before $z_α$ when the fifth force dominates, and does not lead to the existence of clumped matter on extra small scales inside halos for large masses while this possibility exists for masses $M\lesssim 10^{11} M_\odot$ where the phenomenology of ultra-local models would be most different from $Λ$-CDM.

preprint2015arXiv

Cosmological tests of modified gravity: constraints on $F(R)$ theories from the galaxy clustering ratio

The clustering ratio $η$, a large-scale structure observable originally designed to constrain the shape of the power spectrum of matter density fluctuations, is shown to provide a sensitive probe of the nature of gravity in the cosmological regime. We apply this analysis to $F(R)$ theories of gravity using the luminous red galaxy (LRG) sample extracted from the spectroscopic Sloan Digital Sky Survey (SDSS) data release 7 and 10 catalogues. We find that General Relativity (GR), complemented with a Friedmann-Robertson-Walker (FRW) cosmological model with parameters fixed by the Planck satellite, describes extremely well the clustering of galaxies up to $z\sim 0.6$. On large cosmic scales, the absolute amplitude of deviations from GR, $|f_{R_0 }|$, is constrained to be smaller than $4.6 \times 10^{-5}$ at the $95\%$ confidence level. This bound makes cosmological probes of gravity almost competitive with the sensitivity of Solar System tests, although still one order of magnitude less effective than astrophysical tests. We also extrapolate our results to future large surveys like Euclid and show that the astrophysical bound will certainly remain out of reach for such a class of modified-gravity models that only differ from $Λ$CDM at low redshifts.

preprint2015arXiv

Galaxy Clustering in 3D and Modified Gravity Theories

We study Modified Gravity (MG) theories by modelling the redshifted matter power spectrum in a spherical Fourier-Bessel (sFB) basis. We use a fully non-linear description of the real-space matter power-spectrum and include the lowest-order redshift-space correction (Kaiser effect), taking into account some additional non-linear contributions. Ignoring relativistic corrections, which are not expected to play an important role for a shallow survey, we analyse two different modified gravity scenarios, namely the generalised Dilaton scalar-tensor theories and the $f({R})$ models in the large curvature regime. We compute the 3D power spectrum ${\cal C}^s_{\ell}(k_1,k_2)$ for various such MG theories with and without redshift space distortions, assuming precise knowledge of background cosmological parameters. Using an all-sky spectroscopic survey with Gaussian selection function $φ(r)\propto \exp(-{r^2 / r^2_0})$, $r_0 = 150 \, h^{-1} \, {\textrm{Mpc}}$, and number density of galaxies $\bar {\textrm{N}} =10^{-4}\;{\textrm{Mpc}}^{-3}$, we use a $χ^2$ analysis, and find that the lower-order $(\ell \leq 25)$ multipoles of ${\cal C}^s_\ell(k,k')$ (with radial modes restricted to $k < 0.2 \, h \,{\textrm{Mpc}}^{-1}$) can constraint the parameter $f_{R_0}$ at a level of $2\times 10^{-5} (3\times 10^{-5})$ with $3 σ$ confidence for $n=1(2)$. Combining constraints from higher $\ell > 25$ modes can further reduce the error bars and thus in principle make cosmological gravity constraints competitive with solar system tests. However this will require an accurate modelling of non-linear redshift space distortions. Using a tomographic $β(a)$-$m(a)$ parameterization we also derive constraints on specific parameters describing the Dilaton models of modified gravity.

preprint2015arXiv

K-mouflage gravity models that pass Solar System and cosmological constraints

We show that Solar System tests can place very strong constraints on K-mouflage models of gravity, which are coupled scalar field models with nontrivial kinetic terms that screen the fifth force in regions of large gravitational acceleration. In particular, the bounds on the anomalous perihelion of the Moon imposes stringent restrictions on the K-mouflage Lagrangian density, which can be met when the contributions of higher-order operators in the static regime are sufficiently small. The bound on the rate of change of the gravitational strength in the Solar System constrains the coupling strength $β$ to be smaller than $0.1$. These two bounds impose tighter constraints than the results from the Cassini satellite and Big Bang Nucleosynthesis. Despite the Solar System restrictions, we show that it is possible to construct viable models with interesting cosmological predictions. In particular, relative to $Λ$-CDM, such models predict percent-level deviations for the clustering of matter and the number density of dark matter haloes. This makes these models predictive and testable by forthcoming observational missions.

preprint2015arXiv

Linear perturbations in K-mouflage cosmologies with massive neutrinos

We present a comprehensive derivation of linear perturbation equations for different matter species, including photons, baryons, cold dark matter, scalar fields, massless and massive neutrinos, in the presence of a generic conformal coupling. Starting from the Lagrangians, we show how the conformal transformation affects the dynamics. In particular, we discuss how to incorporate consistently the scalar coupling in the equations of the Boltzmann hierarchy for massive neutrinos and the subsequent fluid approximations. We use the recently proposed K-mouflage model as an example to demonstrate the numerical implementation of our linear perturbation equations. K-mouflage is a new mechanism to suppress the fifth force between matter particles induced by the scalar coupling, but in the linear regime the fifth force is unsuppressed and can change the clustering of different matter species in different ways. We show how the CMB, lensing potential and matter power spectra are affected by the fifth force, and find ranges of K-mouflage parameters whose effects could be seen observationally. We also find that the scalar coupling can have the nontrivial effect of shifting the amplitude of the power spectra of the lensing potential and density fluctuations in opposite directions, although both probe the overall clustering of matter. This paper can serve as a reference for those who work on generic coupled scalar field cosmology, or those who are interested in the cosmological behaviour of the K-mouflage model.

preprint2015arXiv

Redshift-space equal-time angular-averaged consistency relations of the gravitational dynamics

We present the redshift-space generalization of the equal-time angular-averaged consistency relations between $(\ell+n)$- and $n$-point polyspectra of the cosmological matter density field. Focusing on the case of $\ell=1$ large-scale mode and $n$ small-scale modes, we use an approximate symmetry of the gravitational dynamics to derive explicit expressions that hold beyond the perturbative regime, including both the large-scale Kaiser effect and the small-scale fingers-of-god effects. We explicitly check these relations, both perturbatively, for the lowest-order version that applies to the bispectrum, and nonperturbatively, for all orders but for the one-dimensional dynamics. Using a large ensemble of $N$-body simulations, we find that our squeezed bispectrum relation is valid to better than $20\%$ up to $1h$Mpc$^{-1}$, for both the monopole and quadrupole at $z=0.35$, in a $Λ$CDM cosmology. Additional simulations done for the Einstein-de Sitter background suggest that these discrepancies mainly come from the breakdown of the approximate symmetry of the gravitational dynamics. For practical applications, we introduce a simple ansatz to estimate the new derivative terms in the relation using only observables. Although the relation holds worse after using this ansatz, we can still recover it within $20\%$ up to $1h$Mpc$^{-1}$, at $z=0.35$ for the monopole. On larger scales, $k = 0.2 h\mathrm{Mpc}^{-1}$, it still holds within the statistical accuracy of idealized simulations of volume $\sim8h^{-3}\mathrm{Gpc}^3$ without shot-noise error.

preprint2015arXiv

Testing Modified Gravity with Cosmic Shear

We use the cosmic shear data from the Canada-France-Hawaii Telescope Lensing Survey to place constraints on $f(R)$ and {\it Generalized Dilaton} models of modified gravity. This is highly complimentary to other probes since the constraints mainly come from the non-linear scales: maximal deviations with respects to the General-Relativity + $Λ$CDM scenario occurs at $k\sim1 h \mbox{Mpc}^{-1}$. At these scales, it becomes necessary to account for known degeneracies with baryon feedback and massive neutrinos, hence we place constraints jointly on these three physical effects. To achieve this, we formulate these modified gravity theories within a common tomographic parameterization, we compute their impact on the clustering properties relative to a GR universe, and propagate the observed modifications into the weak lensing $ξ_{\pm}$ quantity. Confronted against the cosmic shear data, we reject the $f(R)$ $\{ |f_{R_0}|=10^{-4}, n=1\}$ model with more than 99.9% confidence interval (CI) when assuming a $Λ$CDM dark matter only model. In the presence of baryonic feedback processes and massive neutrinos with total mass up to 0.2eV, the model is disfavoured with at least 94% CI in all different combinations studied. Constraints on the $\{ |f_{R_0}|=10^{-4}, n=2\}$ model are weaker, but nevertheless disfavoured with at least 89% CI. We identify several specific combinations of neutrino mass, baryon feedback and $f(R)$ or Dilaton gravity models that are excluded by the current cosmic shear data. Notably, universes with three massless neutrinos and no baryon feedback are strongly disfavoured in all modified gravity scenarios studied. These results indicate that competitive constraints may be achieved with future cosmic shear data.

preprint2014arXiv

Accuracy of analytical models of the large-scale matter distribution

We investigate the possible accuracy that can be reached by analytical models for the matter density power spectrum and correlation function. Using a realistic description of the power spectrum that combines perturbation theory with a halo model, we study the convergence rate of several perturbative expansion schemes and the impact of nonperturbative effects, as well as the sensitivity to phenomenological halo parameters. We check that the simple reorganization of the standard perturbative expansion, with a Gaussian damping prefactor, provides a well-ordered convergence and a finite correlation function that yields a percent accuracy at the baryon acoustic oscillation peak (as soon as one goes to second order). Lagrangian-space expansions are somewhat more efficient, when truncated at low orders, but may diverge at high orders. We find that whereas the uncertainty on the halo-profile mass-concentration relation is not a strong limitation, the uncertainty on the halo mass function can severely limit the accuracy of theoretical predictions for $P(k)$ (this also applies to the power spectra measured in numerical simulations). The real-space correlation function provides a better separation between perturbative and nonperturbative effects, which are restricted to $x \lesssim 10 h^{-1}$Mpc at all redshifts.

preprint2014arXiv

Angular averaged consistency relations of large-scale structures

The cosmological dynamics of gravitational clustering satisfies an approximate invariance with respect to the cosmological parameters that is often used to simplify analytical computations. We describe how this approximate symmetry gives rise to angular averaged consistency relations for the matter density correlations. This allows one to write the $(\ell+n)$ density correlation, with $\ell$ large-scale linear wave numbers that are integrated over angles, and $n$ fixed small-scale nonlinear wave numbers, in terms of the small-scale $n$-point density correlation and $\ell$ prefactors that involve the linear power spectra at the large-scale wave numbers. These relations, which do not vanish for equal-time statistics, go beyond the already known kinematic consistency relations. They could be used to detect primordial non-Gaussianities, modifications of gravity, limitations of galaxy biasing schemes, or to help designing analytical models of gravitational clustering.

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

K-mouflage Cosmology: Formation of Large-Scale Structures

We study structure formation in K-mouflage cosmology whose main feature is the absence of screening effect on quasilinear scales. We show that the growth of structure at the linear level is affected by both a new time dependent Newton constant and a friction term which depend on the background evolution. These combine with the modified background evolution to change the growth rate by up to ten percent since $z\sim 2$. At the one loop level, we find that the nonlinearities of the K-mouflage models are mostly due to the matter dynamics and that the scalar perturbations can be treated at tree level. We also study the spherical collapse in K-mouflage models and show that the critical density contrast deviates from its $Λ$-CDM value and that, as a result, the halo mass function is modified for large masses by an order one factor. Finally we consider the deviation of the matter spectrum from $Λ$-CDM on nonlinear scales where a halo model is utilized. We find that the discrepancy peaks around $1\ h{\rm Mpc}^{-1}$ with a relative difference which can reach fifty percent. Importantly, these features are still true at larger redshifts, contrary to models of the chameleon-$f(R)$ and Galileon types.

preprint2014arXiv

K-mouflage Cosmology: the Background Evolution

We study the cosmology of K-mouflage theories at the background level. We show that the effects of the scalar field are suppressed at high matter density in the early Universe and only play a role in the late time Universe where the deviations of the Hubble rate from its $Λ$-CDM counterpart can be of the order five percent for redshifts $1 \lesssim z \lesssim 5$. Similarly, we find that the equation of state can cross the phantom divide in the recent past and even diverge when the effective scalar energy density goes negative and subdominant compared to matter, preserving the positivity of the squared Hubble rate. These features are present in models for which Big Bang Nucleosynthesis is not affected. We analyze the fate of K-mouflage when the nonlinear kinetic terms give rise to ghosts, particle excitations with negative energy. In this case, we find that the K-mouflage theories can only be considered as an effective description of the Universe at low energy below $1$ keV. In the safe ghost-free models, we find that the equation of state always diverges in the past and changes significantly by a few percent since $z\lesssim 1$.

preprint2014arXiv

Kinematic consistency relations of large-scale structures

We describe how the kinematic consistency relations satisfied by density correlations of the large-scale structures of the Universe can be derived within the usual Newtonian framework. These relations express a kinematic effect and show how the $(\ell+n)$-density correlation factors in terms of the $n$-point correlation and $\ell$ linear power spectrum factors, in the limit where the $\ell$ soft wave numbers become linear and much smaller than the $n$ other wave numbers. We describe how these relations extend to multifluid cases. In the standard cosmology, these consistency relations derive from the equivalence principle. A detection of their violation would indicate non-Gaussian initial conditions, non-negligible decaying modes, or a modification of gravity that does not converge to General Relativity on large scales.

preprint2014arXiv

Small-scale Nonlinear Dynamics of K-mouflage Theories

We investigate the small-scale static configurations of K-mouflage models defined by a general function $K(χ)$ of the kinetic terms. The fifth force is screened by the nonlinear K-mouflage mechanism if $K'(χ)$ grows sufficiently fast for large negative $χ$. In the general non-spherically symmetric case, the fifth force is not aligned with the Newtonian force. For spherically symmetric static matter density profiles, the results depend on the potential function $W_{-}(y) = y K'(-y^2/2)$, which must be monotonically increasing to $+\infty$ for $y \geq 0$ to guarantee the existence of a single solution throughout space for any matter density profile. Small radial perturbations around these static profiles propagate as traveling waves with a velocity greater than the speed of light. Starting from vanishing initial conditions for the scalar field and for a time-dependent matter density corresponding to the formation of an overdensity, we numerically check that the scalar field converges to the static solution. If $W_{-}$ is bounded, for high-density objects there are no static solutions throughout space, but one can still define a static solution restricted to large radii. Our dynamical study shows that the scalar field relaxes to this static solution at large radii, whereas spatial gradients keep growing with time at smaller radii. If $W_{-}$ is not bounded but non-monotonic, there is an infinite number of discontinuous static solutions but they are not physical and these models are not theoretically sound. Such K-mouflage scenarios provide an example of theories that can appear viable at the cosmological level, for the cosmological background and perturbative analysis, while being meaningless at a nonlinear level for small-scale configurations. This shows the importance of small-scale nonlinear analysis of screening models.

preprint2014arXiv

Source-lens clustering and intrinsic-alignment bias of weak-lensing estimators

We estimate the amplitude of the source-lens clustering bias and of the intrinsic-alignment bias of weak lensing estimators of the two-point and three-point convergence and cosmic-shear correlation functions. We use a linear galaxy bias model for the galaxy-density correlations, as well as a linear intrinsic-alignment model. For the three-point and four-point density correlations, we use analytical or semi-analytical models, based on a hierarchical ansatz or a combination of one-loop perturbation theory with a halo model. For two-point statistics, we find that the source-lens clustering bias is typically several orders of magnitude below the weak lensing signal, except when we correlate a very low-redshift galaxy ($z_2 \la 0.05$) with a higher redshift galaxy ($z_1 \ga 0.5$), where it can reach $10\%$ of the signal for the shear. For three-point statistics, the source-lens clustering bias is typically of order $10\%$ of the signal, as soon as the three galaxy source redshifts are not identical. The intrinsic-alignment bias is typically about $10\%$ of the signal for both two-point and three-point statistics. Thus, both source-lens clustering bias and intrinsic-alignment bias must be taken into account for three-point estimators aiming at a better than $10\%$ accuracy.

preprint2014arXiv

Testing the equal-time angular-averaged consistency relation of the gravitational dynamics in N-body simulations

We explicitly test the equal-time consistency relation between the angular-averaged bispectrum and the power spectrum of the matter density field, employing a large suite of cosmological $N$-body simulations. This is the lowest-order version of the relations between $(\ell+n)$-point and $n$-point polyspectra, where one averages over the angles of $\ell$ soft modes. This relation depends on two wave numbers, $k'$ in the soft domain and $k$ in the hard domain. We show that it holds up to a good accuracy, when $k'/k\ll 1$ and $k'$ is in the linear regime, while the hard mode $k$ goes from linear ($0.1\,h\mathrm{Mpc}^{-1}$) to nonlinear ($1.0\,h\mathrm{Mpc}^{-1}$) scales. On scales $k\lesssim 0.4\,h\mathrm{Mpc}^{-1}$, we confirm the relation within the statistical error of the simulations (typically a few percent depending on the wave number), even though the bispectrum can already deviate from leading-order perturbation theory by more than $30\%$. We further examine the relation on smaller scales with higher resolution simulations. We find that the relation holds within the statistical error of the simulations at $z=1$, whereas we find deviations as large as $\sim 7\%$ at $k \sim 1.0\,h\mathrm{Mpc}^{-1}$ at $z=0.35$. We show that this can be explained partly by the breakdown of the approximation $Ω_\mathrm{m}/f^2\simeq1$ with supplemental simulations done in the Einstein-de Sitter background cosmology. We also estimate the impact of this approximation on the power spectrum and bispectrum.

preprint2013arXiv

Impact of a Warm Dark Matter late-time velocity dispersion on large-scale structures

We investigate whether the late-time (at $z\leq 100$) velocity dispersion expected in Warm Dark Matter scenarios could have some effect on the cosmic web (i.e., outside of virialized halos). We consider effective hydrodynamical equations, with a pressurelike term that agrees at the linear level with the analysis of the Vlasov equation. Then, using analytical methods, based on perturbative expansions and the spherical dynamics, we investigate the impact of this term for a 1 keV dark matter particle. We find that the late-time velocity dispersion has a negligible effect on the power spectrum on perturbative scales and on the halo mass function. However, it has a significant impact on the probability distribution function of the density contrast at $z \sim 3$ on scales smaller than $0.1 h^{-1}$Mpc, which correspond to Lyman-$α$ clouds. Finally, we note that numerical simulations should start at $z_i\geq 100$ rather than $z_i \leq 50$ to avoid underestimating gravitational clustering at low redshifts.

preprint2013arXiv

Impact on the power spectrum of Screening in Modified Gravity Scenarios

We study the effects of screened modified gravity of the $f(R)$, dilaton and symmetron types on structure formation, from the quasi-linear to the non-linear regime, using semi-analytical methods. For such models, where the range of the new scalar field is typically within the Mpc range and below in the cosmological context, non-linear techniques are required to understand the deviations of the power spectrum of the matter density contrast compared to the $Λ$-CDM template. This is nowadays commonly tackled using extensive N-body simulations. Here we present new results combining exact perturbation theory at the one loop level (and a partial resummation of the perturbative series) with a halo model. The former allows one to extend the linear perturbative analysis up to $k\lesssim 0.15{\rm h Mpc}^{-1}$ at the perturbative level while the latter leads to a reasonable, up to a few percent, agreement with numerical simulations for $k\lesssim 3{\rm h Mpc}^{-1}$ for large curvature $f(R)$ models, and $k\lesssim 1{\rm h Mpc}^{-1}$ for dilatons and symmetrons, at $z=0$. We also discuss how the behaviors of the perturbative expansions and of the spherical collapse differ for $f(R)$, dilaton, and symmetron models.

preprint2013arXiv

Matter power spectrum from a Lagrangian-space regularization of perturbation theory

We present a new approach to computing the matter density power spectrum, from large linear scales to small highly nonlinear scales. Instead of explicitly computing a partial series of high-order diagrams, as in perturbative resummation schemes, we embed the standard perturbation theory within a realistic nonlinear Lagrangian-space ansatz. We also point out that an "adhesion-like" regularization of the shell-crossing regime is more realistic than a "Zel'dovich-like" behavior, where particles freely escape to infinity. This provides a "cosmic web" power spectrum with good small-scale properties that provide a good matching with a halo model on mildly nonlinear scales. We obtain a good agreement with numerical simulations on large scales, better than 3% for $k\leq 1 h$Mpc$^{-1}$, and on small scales, better than 10% for $k \leq 10 h$Mpc$^{-1}$, at $z \geq 0.35$, which improves over previous methods.

preprint2013arXiv

Redshift-space correlation functions in large galaxy cluster surveys

Large ongoing and upcoming galaxy cluster surveys in the optical, X-ray and millimetric wavelengths will provide rich samples of galaxy clusters at unprecedented depths. One key observable for constraining cosmological models is the correlation function of these objects, measured through their spectroscopic redshift. We study the redshift-space correlation functions of clusters of galaxies, averaged over finite redshift intervals, and their covariance matrices. Expanding as usual the angular anisotropy of the redshift-space correlation on Legendre polynomials, we consider the redshift-space distortions of the monopole as well as the next two multipoles, $2\ell=2$ and 4. Taking into account the Kaiser effect, we developed an analytical formalism to obtain explicit expressions of all contributions to these mean correlations and covariance matrices. We include shot-noise and sample-variance effects as well as Gaussian and non-Gaussian contributions. We obtain a reasonable agreement with numerical simulations for the mean correlations and covariance matrices on large scales ($r> 10 h^{-1}$Mpc). Redshift-space distortions amplify the monopole correlation by about $10-20%$, depending on the halo mass, but the signal-to-noise ratio remains of the same order as for the real-space correlation. This distortion will be significant for surveys such as DES, Erosita, and Euclid, which should also measure the quadrupole $2\ell=2$. The third multipole, $2\ell=4$, may only be marginally detected by Euclid.

preprint2012arXiv

Covariance matrices for halo number counts and correlation functions

We study the mean number counts and two-point correlation functions, along with their covariance matrices, of cosmological surveys such as for clusters. In particular, we consider correlation functions averaged over finite redshift intervals, which are well suited to cluster surveys or populations of rare objects, where one needs to integrate over nonzero redshift bins to accumulate enough statistics. We develop an analytical formalism to obtain explicit expressions of all contributions to these means and covariance matrices, taking into account both shot-noise and sample-variance effects. We compute low-order as well as high-order (including non-Gaussian) terms. We derive expressions for the number counts per redshift bins both for the general case and for the small window approximation. We estimate the range of validity of Limber's approximation and the amount of correlation between different redshift bins. We also obtain explicit expressions for the integrated 3D correlation function and the 2D angular correlation. We compare the relative importance of shot-noise and sample-variance contributions, and of low-order and high-order terms. We check the validity of our analytical results through a comparison with the Horizon full-sky numerical simulations, and we obtain forecasts for several future cluster surveys.

preprint2012arXiv

Modeling of weak lensing statistics. I. Power spectrum and bispectrum

We investigate the performance of an analytic model of the 3D matter distribution, which combines perturbation theory with halo models, for weak-lensing statistics. We compare our predictions for the weak-lensing convergence power spectrum and bispectrum with numerical simulations and fitting formulas proposed in previous works. We find that this model provides better agreement with simulations than published fitting formulas. This shows that building on systematic and physically motivated models is a promising approach. Moreover, this makes explicit the link between the weak-lensing statistics and the underlying properties of the 3D matter distribution, as a function of scale $\ell$. Thus, we obtain the contributions to the lensing power spectrum and bispectrum that arise from perturbative terms (complete up to one-loop) and nonperturbative terms (e.g., "1-halo" term). Finally, we show that this approach recovers the dependence on cosmology (for realistic scenarios).

preprint2012arXiv

Modeling of weak lensing statistics. II. Configuration-space statistics

We investigate the performance of an analytic model of the 3D matter distribution, which combines perturbation theory with halo models, for weak-lensing configuration-space statistics. We compared our predictions for the weak-lensing convergence two-point and three-point correlation functions with numerical simulations and fitting formulas proposed in previous works. We also considered the second- and third-order moments of the smoothed convergence and of the aperture-mass. As in our previous study of Fourier-space weak-lensing statistics, we find that our model agrees better with simulations than previously published fitting formulas. Moreover, we recover the dependence on cosmology of these weak-lensing statistics and we can describe multi-scale moments. This approach allows us to obtain the quantitative relationship between these integrated weak-lensing statistics and the various contributions to the underlying 3D density fluctuations, decomposed over perturbative, two-halo, or one-halo terms.

preprint2012arXiv

Structure Formation in Modified Gravity Scenarios

We study the growth of structures in modified gravity models where the Poisson equation and the relationship between the two Newtonian potentials are modified by explicit functions of space and time. This parameterisation applies to the $f(R)$ models and more generally to screened modified gravity models. We investigate the linear and weakly nonlinear regimes using the "standard" perturbative approach and a resummation technique, while we use the spherical dynamics to go beyond low-order results. This allows us to estimate the matter density power spectrum and bispectrum from linear to highly nonlinear scales, the full probability distribution of the density contrast on weakly nonlinear scales, and the halo mass function. We analyse the impact of modifications of gravity on these quantities for a few realistic models. In particular, we find that the standard one-loop perturbative approach is not sufficiently accurate to probe these effects on the power spectrum and it is necessary to use resummation methods even on weakly nonlinear scales which provide the best observational window for modified gravity as relative deviations from General Relativity do not grow significantly on smaller scales where theoretical predictions become increasingly difficult.

preprint2011arXiv

Combining perturbation theories with halo models

We investigate the building of unified models that can predict the matter-density power spectrum and the two-point correlation function from very large to small scales, being consistent with perturbation theory at low $k$ and with halo models at high $k$. We use a Lagrangian framework to re-interpret the halo model and to decompose the power spectrum into "2-halo" and "1-halo" contributions, related to "perturbative" and "non-perturbative" terms. We describe a simple implementation of this model and present a detailed comparison with numerical simulations, from $k \sim 0.02$ up to $100 h$Mpc$^{-1}$, and from $x \sim 0.02$ up to $150 h^{-1}$Mpc. We show that the 1-halo contribution contains a counterterm that ensures a $k^2$ tail at low $k$ and is important not to spoil the predictions on the scales probed by baryon acoustic oscillations, $k \sim 0.02$ to $0.3 h$Mpc$^{-1}$. On the other hand, we show that standard perturbation theory is inadequate for the 2-halo contribution, because higher order terms grow too fast at high $k$, so that resummation schemes must be used. We describe a simple implementation, based on a 1-loop "direct steepest-descent" resummation for the 2-halo contribution that allows fast numerical computations, and we check that we obtain a good match to simulations at low and high $k$. Our simple implementation already fares better than standard 1-loop perturbation theory on large scales and simple fits to the power spectrum at high $k$, with a typical accuracy of 1% on large scales and 10% on small scales. We obtain similar results for the two-point correlation function. However, there remains room for improvement on the transition scale between the 2-halo and 1-halo contributions, which may be the most difficult regime to describe.

preprint2011arXiv

Combining perturbation theories with halo models for the matter bispectrum

We investigate how unified models should be built to be able to predict the matter-density bispectrum (and power spectrum) from very large to small scales and that are at the same time consistent with perturbation theory at low $k$ and with halo models at high $k$. We use a Lagrangian framework to decompose the bispectrum into "3-halo", "2-halo", and "1-halo" contributions, related to "perturbative" and "non-perturbative" terms. We describe a simple implementation of this approach and present a detailed comparison with numerical simulations. We show that the 1-halo and 2-halo contributions contain counterterms that ensure their decay at low $k$, as required by physical constraints, and allow a better match to simulations. Contrary to the power spectrum, the standard 1-loop perturbation theory can be used for the perturbative 3-halo contribution because it does not grow too fast at high $k$. Moreover, it is much simpler and more accurate than two resummation schemes investigated in this paper. We obtain a good agreement with numerical simulations on both large and small scales, but the transition scales are poorly described by the simplest implementation. This cannot be amended by simple modifications to the halo parameters, but we show how it can be corrected for the power spectrum and the bispectrum through a simple interpolation scheme that is restricted to this intermediate regime. Then, we reach an accuracy on the order of 10% on mildly and highly nonlinear scales, while an accuracy on the order of 1% is obtained on larger weakly nonlinear scales. This also holds for the real-space two-point correlation function.

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

Impact of shell crossing and scope of perturbative approaches in real and redshift space

We study the effect of nonperturbative corrections associated with the behavior of particles after shell crossing on the matter power spectrum. We compare their amplitude with the perturbative terms that can be obtained within the fluid description of the system, to estimate the range of scales where such perturbative approaches are relevant. We use the simple Zeldovich dynamics as a benchmark, as it allows the exact computation of the full nonlinear power spectrum and of perturbative terms at all orders. Then, we introduce a "sticky model" that coincides with the Zeldovich dynamics before shell crossing but shows a different behavior afterwards. Thus, their power spectra only differ in their nonperturbative terms. We consider both the real-space and redshift-space power spectra. We find that the potential of perturbative schemes is greater at higher redshift for a $Λ$CDM cosmology. For the real-space power spectrum, one can go up to order $66$ of perturbation theory at $z=3$, and to order $9$ at $z=0$, before the nonperturbative correction surpasses the perturbative correction of that order. This allows us to increase the upper bound on $k$ where systematic theoretical predictions may be obtained by perturbative schemes, beyond the linear regime, by a factor $\sim 26$ at $z=3$ and $\sim 6.5$ at $z=0$. This provides a strong motivation to study perturbative resummation schemes, especially at high redshifts $z \geq 1$. We find similar results for the redshift-space power spectrum, with characteristic wavenumbers that are shifted to lower values as redshift-space distortions amplify higher order terms of the perturbative expansions while decreasing the resummed nonlinear power at high $k$.

preprint2011arXiv

Large-scale bias of dark matter halos

We build a simple analytical model for the bias of dark matter halos that applies to objects defined by an arbitrary density threshold, $200\leq\deltas\leq 1600$, and that provides accurate predictions from low-mass to high-mass halos. We point out that it is possible to build simple and efficient models, with no free parameter for the halo bias, by using integral constraints that govern the behavior of low-mass and typical halos, whereas the properties of rare massive halos are derived through explicit asymptotic approaches. We also describe how to take into account the impact of halo motions on their bias, using their linear displacement field. We obtain a good agreement with numerical simulations for the halo mass functions and large-scale bias at redshifts $0\leq z \leq 2.5$, for halos defined by a nonlinear density threshold $200\leq\deltas\leq 1600$. We also evaluate the impact on the halo bias of two common approximations, i) neglecting halo motions, and ii) linearizing the halo two-point correlation.

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.

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.

preprint2009arXiv

Secondary non-Gaussianity and Cross-Correlation Analysis

We develop optimised estimators of two sorts of power spectra for fields defined on the sky, in the presence of partial sky coverage. The first is the cross-power spectrum of two fields on the sky; the second is the skew spectrum of three fields. These can probe the Integrated Sachs Wolfe Effect (ISW) at large angular scales and the Sunyaev Zéldovich (SZ) effect from hot gas in clusters at small angular scales. The skew spectrum, recently introduced by Munshi & Heavens (2009), is an optimised statistic which can be tuned to study a particular form of non-Gaussianity, such as may arise in the early Universe, but which retains information on the nature of non-Gaussianity. In this paper we develop the mathematical formalism for the skew spectrum of 3 different fields. When applied to the CMB, this allows us to explore the contamination of the skew spectrum by secondary sources of CMB fluctuations. Considering the three-point function, the study of the bispectrum provides valuable information regarding cross-correlation of secondaries with lensing of CMB with much higher significance compared to just the study involving CMB sky alone. We use our analytical models to study specific cases of cosmological interest which include cross-correlating CMB with various large scale tracers to probe ISW and SZ effects for cross spectral analysis and use the formalism to study the signal-to-noise ratio for detection of the weak lensing of the CMB by cross-correlating it with different tracers as well as point sources for CMB experiments such as Planck (abridged).