Next Article in Journal
Semiheaps and Ternary Algebras in Quantum Mechanics Revisited
Next Article in Special Issue
Vacuum Energy in Saez-Ballester Theory and Stabilization of Extra Dimensions
Previous Article in Journal
Constraints on the Production of Phosphine by Venusian Volcanoes
Previous Article in Special Issue
Estimating the Cosmological Constant from Shadows of Kerr–de Sitter Black Holes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

One-Point Statistics Matter in Extended Cosmologies

School of Mathematics, Statistics and Physics, Newcastle University, Herschel Building, Newcastle upon Tyne NE1 7RU, UK
*
Author to whom correspondence should be addressed.
Universe 2022, 8(1), 55; https://doi.org/10.3390/universe8010055
Submission received: 8 December 2021 / Revised: 11 January 2022 / Accepted: 12 January 2022 / Published: 17 January 2022
(This article belongs to the Special Issue Alternative Gravities and Fundamental Cosmology)

Abstract

:
The late universe contains a wealth of information about fundamental physics and gravity, wrapped up in non-Gaussian fields. To make use of as much information as possible, it is necessary to go beyond two-point statistics. Rather than going to higher-order N-point correlation functions, we demonstrate that the probability distribution function (PDF) of spheres in the matter field (a one-point function) already contains a significant amount of this non-Gaussian information. The matter PDF dissects different density environments which are lumped together in two-point statistics, making it particularly useful for probing modifications of gravity or expansion history. Our approach in Cataneo et al. 2021 extends the success of Large Deviation Theory for predicting the matter PDF in Λ CDM in these “extended” cosmologies. A Fisher forecast demonstrates the information content in the matter PDF via constraints for a Euclid-like survey volume combining the 3D matter PDF with the 3D matter power spectrum. Adding the matter PDF halves the uncertainties on parameters in an evolving dark energy model, relative to the power spectrum alone. Additionally, the matter PDF contains enough non-linear information to substantially increase the detection significance of departures from General Relativity, with improvements up to six times the power spectrum alone. This analysis demonstrates that the matter PDF is a promising non-Gaussian statistic for extracting cosmological information, particularly for beyond Λ CDM models.

1. Introduction

In recent decades, cosmology has moved solidly into a data-driven science. The current standard model of cosmology, called Λ CDM, consists of a cosmological constant as the dark energy component ( Λ ), and cold (non-relativistic) dark matter (CDM) as its principle components.
The parameters of the Λ CDM model are most tightly constrained currently by experiments measuring the temperature anisotropies and polarisation in the cosmic microwave background (CMB) (for example, the Planck measurements in [1]). However, while CMB data is very valuable in extracting cosmological information, in the push to sub-percent measurements of standard cosmological parameters, and in testing non-standard cosmologies, the large-scale structure (LSS) of the universe is the most promising complementary tool.
The principle advantage of LSS data is that it is three dimensional, tracing a history of how cosmic structure evolves over time, since the snapshot of the CMB. By counting the Fourier modes available, one can expect 1–2 orders of magnitude improvement on constraints from LSS data (see, e.g., [2]). In particular, the large-scale structure provides a window to the expansion history of the universe, which makes it an exciting probe of dynamical dark energy and modifications to gravity. These extensions to the standard cosmology are one of the principle science goals of current and upcoming missions like Euclid [3], LSST [4], and DESI [5]. These extensions to the standard Λ CDM model would represent new fundamental physics, and could resolve certain current observational tensions which Λ CDM is strained to explain (recently reviewed in, e.g., [6,7,8]).
However, extracting information from the late universe via large-scale structure is non-trivial for several reasons. The first and most relevant for this work is that the late universe is statistically much more complex than the universe at the time of the CMB. Extracting cosmological information is always done on a statistical basis, treating observables such as a field of galaxy positions or shear lensing maps as single realisations of a random field. Gaussian random fields are completely characterised by their two-point correlation function (or its Fourier counterpart the power spectrum), as all higher-order correlations functions can be written as sums of the two-point function via Wick’s/Isserlis’ theorem. The CMB has been measured to be a near-perfect Gaussian random field [9], and so measurement of the power spectrum is sufficient to quantify all information content in the CMB. However, as gravitational collapse is a non-linear process, the statistics of the late time density are not also Gaussian, as the non-linearity in the mapping from initial to final densities sources non-trivial higher statistics as information “leaks out” of the power spectrum. It is therefore crucial to determine statistics beyond the power spectrum which can recapture this non-linear information in an efficient and theoretically tractable way. Other reasons why extracting information from LSS data is difficult come down to issues of modelling dynamics on non-linear scales (often side-stepped by running N-body simulations, which are expensive and have their own host of non-trivialities), and on a variety of systematic effects in observations.

2. Methods

2.1. The Matter PDF in Spheres from Large Deviation Theory

This work focuses on a simple choice of non-Gaussian statistic, namely the probability distribution function (PDF) of matter density in spheres. The matter PDF can be straightforwardly calculated from the density field of a cosmological N-body simulation by looking at the distribution of the matter field smoothed with a spherical top-hat filter on the scale of interest.
The work of several papers [10,11,12,13,14] provide an analytic framework for predicting the matter PDF in spheres in the mildly non-linear regime. This method relies on large deviations theory (LDT) where the driving parameter is the non-linear variance, σ NL 2 , of the matter field. This formalism therefore remains valid at redshifts z and for spheres of radius R where σ NL 2 ( z , R ) < 1 .
For Gaussian initial conditions, the PDF, P lin ( δ L ) , of the linear matter density contrast, δ L , in a sphere of radius r is a Gaussian distribution with width given by the linear variance at scale r and redshift z
P r , z lin ( δ L ) = Ψ r , z lin ( δ L ) 2 π exp Ψ r , z lin ( δ L ) , Ψ r , z lin ( δ L ) = δ L 2 2 σ L 2 ( r , z ) .
The linear variance on scale r is given by an integral over the linear power spectrum P L with a spherical top-hat filter in position space
σ L 2 ( r , z ) = d k 2 π 2 P L ( k , z ) k 2 W 3 D 2 ( k r ) , W 3 D 2 ( k ) = 3 π 2 J 3 / 2 ( k ) k 3 / 2 ,
where W 3 D ( k ) is the Fourier transform of the 3D spherical top-hat filter, and  J 3 / 2 ( k ) is the Bessel function of the first kind of order 3 / 2 .
The function Ψ in Equation (1) is related to the rate function in the context of LDT. The key result from the LDT formalism allows us to relate the rate function of the linear density to the non-linear density, which provides the exponential dependence of the final PDF. Generally this LDT result is called the contraction principle and relates the rate function of different random variables. In the cosmological case, since large deviations are exponentially unlikely and the matter PDF is computed for spherically symmetric cells, the  most likely mapping from linear to non-linear densities should be dominated by spherical collapse, δ L ρ SC ( δ L ) . Combining this with mass conservation in spheres (which relates the initial scale, r, to the final scale, R, by  r = R ρ 1 / 3 ) leads to the final decay-rate function Ψ R , z of the non-linear density
Ψ R , z ( ρ ) = σ L 2 ( R , z ) σ L 2 ( R ρ 1 / 3 , z ) δ L SC ( ρ ) 2 2 σ NL 2 ( R , z ) .
The LDT model then predicts the matter PDF in spheres is given by
P R , z ( ρ ) exp ( Ψ R , z ( ρ ) ) ,
where the precise prefactor can be determined by a more detailed analysis (see Equation (5a,b) from [15]).
For a standard Λ CDM universe, there are only three quantities needed for this theoretical model of the matter PDF (through the decay-rate function in Equation (3)):
(i)
The time- and scale-dependence of the linear variance σ L 2 ( r , z ) .
(ii)
The non-linear variance of the log-density σ ln ρ , NL 2 ( R , z ) .
(iii)
The mapping between linear and final densities in spheres, which is taken to be spherical collapse δ L ρ SC ( δ L ) (or its inverse δ L SC ( ρ ) ) .
Figure 1 shows the success of this LDT model against measured PDFs from the Quijote simulations [16], as well as in comparison to a common log-normal phenomenological model (dashed lines).

2.2. Extended Cosmologies

Due to the significant non-Gaussian information in the matter PDF and the success of this LDT formalism in Λ CDM, we modify this framework to analyse cosmologies with non-GR theories of gravity or dynamical theories of dark energy beyond a cosmological constant. Collectively we will refer to either modified gravity (MG) or dark energy (DE) models as extended cosmologies.
The dark energy model considered was a simple parametrisation of an evolving dark energy, called the w 0 w a CDM model. This cosmology is still described by a smooth dark energy and General Relativity (GR), but with dark energy equation of state given by [17,18]
w ( a ) = w 0 + w a ( 1 a ) ,
where { w 0 , w a } are new phenomenological parameters with w 0 = 1 , w a = 0 corresponding to the cosmological constant.
For the theories of modified gravity, we considered Hu-Sawiki f ( R ) gravity [19] and the normal branch of DGP braneworld gravity which acts as an additional smooth dark energy component [20]. The strength of the deviation from GR gravity in these theories is quantified by the parameters f R 0 and Ω rc respectively. In moving to an extended cosmology, all three of the ingredients outlined in the previous section in principle need updating: the linear variance, the non-linear variance of the log-density, and the mapping between initial and final densities.
Updating the linear variance is straightforward, simply requiring integrating the modified linear power spectrum which can be achieved by the ratio of linear growth factors. This is achieved using a novel technique for emulating the response of the Λ CDM power spectrum to MG/DE effects as described in [21]. The non-linear variance can either be measured from a set of simulations in the extended cosmology, or can be well approximated by a phenomenological log-normal rescaling (Equation (14) from [22]) applied to a reference cosmology.
The mapping from initial to final densities is trickiest for extended cosmologies, especially in the case of scale dependent modified gravity. However, in [15] we show that when restricted to mildly non-linear scales ( R 10 Mpc / h ) we can use the spherical collapse mapping for an Einstein de Sitter cosmology, modified just by the difference in linear growth factors between the Λ CDM and extended cosmology.

2.3. Simulations and Model Validation

In [15], we validated the predicted matter PDFs against a suite of modified gravity and dark energy simulations, and found them to be accurate to 2% over the range of densities used in the final analysis. All PDFs and the analysis presented here are based on the LDT model, calculated using pyLDT (https://github.com/mcataneo/pyLDT-cosmo, accessed on 4 December 2021) a modularised and user-friendly Python code that takes advantage of the PyJulia interface for computationally intensive tasks. The PDFs used in this paper were generated using version 0.4.9 of pyLDT.

3. Results

The following results summarise the analysis presented in Cataneo et al. 2021 [15].

3.1. Matter PDFs in Extended Cosmologies

Figure 2 shows the matter PDF in 10 h 1 Mpc spheres for the three different theories of gravity considered, Λ CDM, f ( R ) , and DGP. Generically, introducing modified gravity will change both the width and the shape of the PDF (c.f. Figure 3 of [15]). Since σ 8 sets the width of the PDF, normalising the cosmologies to have the same σ 8 at redshift 0 allows us to isolate the distinct features of modified gravity on the PDF, as done in Figure 2.
By normalising the width of the PDF at redshift 0, we see a residual difference in the shape of the matter PDF, as well as a distinct redshift dependence. These differences in shape and redshift dependence (as well as a difference in scale dependence not shown in Figure 2) are what allows the PDF to break degeneracies between standard cosmological parameters and modified gravity. The effect of an evolving dark energy on the matter PDF is similar to that of DGP gravity, entering mostly by modifying the expansion history.

3.2. Forecasting Constraining Power with the Fisher Formalism

To forecast errors on a set of cosmological parameters, θ , we make use of the Fisher matrix formalism. This formalism provides constraints on θ under the assumption that the likelihood is approximated by a multivariate Gaussian distribution.
Given a (set of) summary statistics arranged in a data vector, S with components S α , and the covariance matrix between those summary statistics, C d , the components of the Fisher matrix, F, are defined as (assuming that the data covariance matrix is independent of cosmological parameters)
F i j = α , β S α θ i ( C d 1 ) α β S β θ j .
Assuming the likelihood of our observed statistics given our cosmological parameters is well approximated by a multivariate Gaussian, the parameter covariance matrix C p ( θ ) (encoding the expected parameter constraints along with parameter degeneracies) is given by the inverse of the Fisher matrix. To obtain constraints by marginalising over a subset of parameters, one can simply select the appropriate elements of the parameter covariance matrix. In particular, this implies that the fully marginalised constraint on a single parameter, θ i , is given by
σ [ θ i ] = ( F 1 ) i i .
For this analysis, we used three different data vectors. These are the PDF alone, the matter power spectrum alone, and a stacked data vector which combines both the PDF and the matter power spectrum. We combined information from three redshifts, z = 0 , 0.5 , 1 . The data covariance was measured from the fiducial Λ CDM runs of the Quijote suite of simulations [16] subsequently rescaled to correspond to a Euclid-like survey volume of 20 ( Gpc / h ) 3 , for details see Cataneo et al. 2021 [15]. This covariance matrix encodes statistical errors in the matter PDF and power spectrum, as well as their correlation in the combined data vector. For the matter PDF, we combined three sizes of spheres ( 10 , 15 , 20 Mpc / h ), while for the matter power spectrum, we included information up to k max = 0.2 h / Mpc . This maximum wavenumber was chosen to ensure that theoretical predictions of the power spectrum are consistent with those measured from simulations, see Appendix C of [15].

3.3. Response of the PDF to Changes in Cosmological Parameters

Figure 3 and Figure 4 illustrate how the matter PDF depends on standard parameters related to structure formation Ω m and σ 8 , as well as the parameters extending Λ CDM, namely | f R 0 | , w 0 , and w a ( Ω rc for DGP gravity is omitted here to avoid cluttering the plots). This more directly quantifies how well degeneracies in these parameters can be broken, and allows us to quantify the heuristic understanding gained by Figure 2. These derivatives directly enter into the Fisher constraints via Equation (6).
Notice that the effect of Ω m on the matter PDF can easily be disentangled from the other parameters in both the DE and MG case, as the matter PDF is sensitive to Ω m only through its skewness and the linear growth factor, D ( z ) [22].
In f ( R ) gravity, we can expect the most degeneracy breaking, and therefore better constraints on f R 0 than Ω rc or { w 0 , w a } for two main reasons which can be seen in Figure 3. The first is that the f R 0 derivative has a different shape from the σ 8 derivative, showing up as an additional skewness owing to the scale dependent fifth force. This, combined with the fact that the f R 0 derivatives are non-zero at z = 0 (unlike in DGP and w 0 w a CDM), allows more information to be extracted from the non-linear regime. While DGP gravity is not shown in Figure 3, its effect is very similar to that of dark energy shown in Figure 4 (as can be seen in Figure 11 from [15]). Figure 4 shows that at fixed scale and redshift, the response of the PDF to w 0 or w a is very similar in shape to the response to σ 8 . For this reason, in the cases of dark energy and DGP gravity, the degeneracy is mainly broken by a difference in the redshift (and scale) dependence.

3.4. Fisher Forecasts for Modified Gravity Detection and Dark Energy Constraints

Figure 5 and Figure 6 show the Fisher forecast constraints for f ( R ) gravity with | f R 0 | = 10 6 and for w 0 w a CDM about a Λ CDM fiducial (forecasts for DGP gravity can be found in [15]). In both the modified gravity and the dark energy extended cosmologies, the matter PDF alone provides constraints competitive with the matter power spectrum. More importantly, however, they provide complementary information, demonstrated by the different degeneracy directions. This indicates that the matter PDF is recovering independent non-Gaussian information beyond the power spectrum.
Combining the PDF and the power spectrum allows for a 5 σ detection of both f ( R ) and DGP gravity (see Table 1), and at least doubles the constraining power for other parameters such as σ 8 over just power spectrum alone. For evolving dark energy, the improvement is quantified by the dark energy Figure of Merit (FoM), equal to the inverse area of the contour in the w 0 - w a plane. Adding PDF information to the power spectrum increases the FoM by a factor of 5 (summarised in Table 2). The resulting FoM is in the range expected to be reached by Euclid in combining spectroscopic galaxy clustering and weak lensing [23].

4. Conclusions

Standard two-point statistics are not sufficient to make full use of the information content in the cosmic large-scale structure, and would leave large amounts of data from current and upcoming galaxy surveys under utilised. The full shape of the matter density PDF in spheres has been shown to provide great complementarity to the standard two-point statistics, and allows extraction of information from the non-linear regime. The analytic framework described here has been successfully applied to Λ CDM universes along with extensions including primordial non-Gaussianity [24] and massive neutrinos [22]. This work demonstrates that the LDT formalism continues to work in modified gravity and dark energy scenarios , providing a powerful non-Gaussian probe of fundamental physics complementary to two-point statistics.
While the analysis presented here is idealised in that it relies on knowledge of the true matter distribution, it is encouraging for realistic scenarios. In the case of Λ CDM cosmologies, the LDT approach has been translated into several observable quantities, including weak lensing [25,26,27], galaxy clustering [28,29], and density-split statistics [30,31]. Given the theoretical information content in the matter PDF demonstrated here, extending the LDT framework to observables in the context of modified gravity would be a worthwhile endeavour for constraining both astrophysical (e.g., baryonic feedback, intrinsic alignment, and galaxy bias) and cosmological parameters to complement two-point statistics.

Author Contributions

This text was written by A.G. and edited by C.U. All authors have read and agreed to the published version of the manuscript.

Funding

A.G. is supported by an EPSRC studentship under Project 2441314 from UK Research & Innovation.

Data Availability Statement

The code for computing the matter PDFs, as well as the simulations used can be found in the data availability section of [15].

Acknowledgments

The figures in this work were created with Matplotlib [32] and ChainCosumer [33], making use of the NumPy [34] and SciPy [35] Python libraries.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
Λ CDMLambda Cold Dark Matter (model of cosmology)
CMBCosmic Microwave Background
DEDark Energy
DGPDvali-Gabadadze-Porrati (model of gravity)
FoM(Dark Energy) Figure of Merit
GRGeneral Relativity
LDTLarge Deviations Theory
LSSLarge-Scale Structure (of the universe)
MGModified Gravity
PDFProbability Distribution Function

References

  1. Planck Collaboration; Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; et al. Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys. 2020, 641, A6. [Google Scholar] [CrossRef] [Green Version]
  2. Carrasco, J.J.M.; Hertzberg, M.P.; Senatore, L. The effective field theory of cosmological large scale structures. J. High Energy Phys. 2012, 2012, 82. [Google Scholar] [CrossRef] [Green Version]
  3. Laureijs, R.; Amiaux, J.; Arduini, S.; Auguères, J.L.; Brinchmann, J.; Cole, R.; Cropper, M.; Dabin, C.; Duvet, L.; Ealet, A.; et al. Euclid Definition Study Report. arXiv 2011, arXiv:1110.3193. [Google Scholar]
  4. Ivezić, Ž.; Kahn, S.M.; Tyson, J.A.; Abel, B.; Acosta, E.; Allsman, R.; Alonso, D.; AlSayyad, Y.; Anderson, S.F.; Andrew, J.; et al. LSST: From Science Drivers to Reference Design and Anticipated Data Products. Astrophys. J. 2019, 873, 111. [Google Scholar] [CrossRef]
  5. Levi, M.; Bebek, C.; Beers, T.; Blum, R.; Cahn, R.; Eisenstein, D.; Flaugher, B.; Honscheid, K.; Kron, R.; Lahav, O.; et al. The DESI Experiment, a whitepaper for Snowmass 2013. arXiv 2013, arXiv:1308.0847. [Google Scholar]
  6. Douspis, M.; Salvati, L.; Aghanim, N. On the tension between Large Scale Structures and Cosmic Microwave Background. arXiv 2019, arXiv:1901.05289. [Google Scholar]
  7. Di Valentino, E.; Anchordoqui, L.A.; Akarsu, Ö.; Ali-Haimoud, Y.; Amendola, L.; Arendse, N.; Asgari, M.; Ballardini, M.; Basilakos, S.; Battistelli, E.; et al. Cosmology Intertwined III: Fσ8 and S8. Astropart. Phys. 2021, 131, 102604. [Google Scholar] [CrossRef]
  8. Perivolaropoulos, L.; Skara, F. Challenges for ΛCDM: An update. arXiv 2021, arXiv:2105.05208. [Google Scholar]
  9. Planck Collaboration; Akrami, Y.; Arroja, F.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; et al. Planck 2018 results. IX. Constraints on primordial non-Gaussianity. Astron. Astrophys. 2020, 641, A9. [Google Scholar] [CrossRef]
  10. Bernardeau, F. The nonlinear evolution of rare events. Astrophys. J. 1994, 427, 51–71. [Google Scholar] [CrossRef]
  11. Valageas, P. Dynamics of gravitational clustering. II. Steepest-descent method for the quasi-linear regime. Astron. Astrophys. 2002, 382, 412–430. [Google Scholar] [CrossRef] [Green Version]
  12. Bernardeau, F.; Pichon, C.; Codis, S. Statistics of cosmic density profiles from perturbation theory. Phys. Rev. D 2014, 90, 103519. [Google Scholar] [CrossRef] [Green Version]
  13. Bernardeau, F.; Reimberg, P. Large deviation principle at play in large scale structure cosmology. Phys. Rev. D 2016, 94, 063520. [Google Scholar] [CrossRef] [Green Version]
  14. Uhlemann, C.; Codis, S.; Pichon, C.; Bernardeau, F.; Reimberg, P. Back in the saddle: Large-deviation statistics of the cosmic log-density field. Mon. Not. R. Astron. Soc. 2016, 460, 1529–1541. [Google Scholar] [CrossRef] [Green Version]
  15. Cataneo, M.; Uhlemann, C.; Arnold, C.; Gough, A.; Li, B.; Heymans, C. The matter density PDF for modified gravity and dark energy with Large Deviations Theory. arXiv 2021, arXiv:2109.02636. [Google Scholar]
  16. Villaescusa-Navarro, F.; Hahn, C.; Massara, E.; Banerjee, A.; Delgado, A.M.; Ramanah, D.K.; Charnock, T.; Giusarma, E.; Li, Y.; Allys, E.; et al. The Quijote Simulations. Astrophys. J. Suppl. Ser. 2020, 250, 2. [Google Scholar] [CrossRef]
  17. Chevallier, M.; Polarski, D. Accelerating Universes with Scaling Dark Matter. Int. J. Mod. Phys. D 2001, 10, 213–223. [Google Scholar] [CrossRef] [Green Version]
  18. Linder, E.V. Exploring the Expansion History of the Universe. Phys. Rev. Lett. 2003, 90, 091301. [Google Scholar] [CrossRef] [Green Version]
  19. Hu, W.; Sawicki, I. Models of f(R) cosmic acceleration that evade solar system tests. Phys. Rev. D 2007, 76, 064004. [Google Scholar] [CrossRef] [Green Version]
  20. Schmidt, F. Cosmological simulations of normal-branch braneworld gravity. Phys. Rev. D 2009, 80, 123003. [Google Scholar] [CrossRef] [Green Version]
  21. Cataneo, M.; Lombriser, L.; Heymans, C.; Mead, A.J.; Barreira, A.; Bose, S.; Li, B. On the road to percent accuracy: Non-linear reaction of the matter power spectrum to dark energy and modified gravity. Mon. Not. R. Astron. Soc. 2019, 488, 2121–2142. [Google Scholar] [CrossRef]
  22. Uhlemann, C.; Friedrich, O.; Villaescusa-Navarro, F.; Banerjee, A.; Codis, S. Fisher for complements: Extracting cosmology and neutrino mass from the counts-in-cells PDF. Mon. Not. R. Astron. Soc. 2020, 495, 4006–4027. [Google Scholar] [CrossRef]
  23. Planck Collaboration; Blanchard, A.; Camera, S.; Carbone, C.; Cardone, V.F.; Casas, S.; Clesse, S.; Ilić, S.; Kilbinger, M.; Kitching, T.; et al. Euclid preparation. VII. Forecast validation for Euclid cosmological probes. Astron. Astrophys. 2020, 642, A191. [Google Scholar] [CrossRef]
  24. Friedrich, O.; Uhlemann, C.; Villaescusa-Navarro, F.; Baldauf, T.; Manera, M.; Nishimichi, T. Primordial non-Gaussianity without tails-how to measure fNL with the bulk of the density PDF. Mon. Not. R. Astron. Soc. 2020, 498, 464–483. [Google Scholar] [CrossRef]
  25. Barthelemy, A.; Codis, S.; Uhlemann, C.; Bernardeau, F.; Gavazzi, R. A nulling strategy for modelling lensing convergence in cones with large deviation theory. Mon. Not. R. Astron. Soc. 2020, 492, 3420–3439. [Google Scholar] [CrossRef]
  26. Boyle, A.; Uhlemann, C.; Friedrich, O.; Barthelemy, A.; Codis, S.; Bernardeau, F.; Giocoli, C.; Baldi, M. Nuw CDM cosmology from the weak lensing convergence PDF. arXiv 2020, arXiv:2012.07771. [Google Scholar] [CrossRef]
  27. Thiele, L.; Hill, J.C.; Smith, K.M. Accurate analytic model for the weak lensing convergence one-point probability distribution function and its autocovariance. Phys. Rev. D 2020, 102, 123545. [Google Scholar] [CrossRef]
  28. Repp, A.; Szapudi, I. Galaxy bias and σ8 from counts in cells from the SDSS main sample. Mon. Not. R. Astron. Soc. Lett. 2020, 498, L125–L129. [Google Scholar] [CrossRef]
  29. Friedrich, O.; Halder, A.; Boyle, A.; Uhlemann, C.; Britt, D.; Codis, S.; Gruen, D.; Hahn, C. The PDF perspective on the tracer-matter connection: Lagrangian bias and non-Poissonian shot noise. arXiv 2021, arXiv:2107.02300. [Google Scholar] [CrossRef]
  30. Gruen, D.; Friedrich, O.; Krause, E.; DeRose, J.; Cawthon, R.; Davis, C.; Elvin-Poole, J.; Rykoff, E.S.; Wechsler, R.H.; Alarcon, A.; et al. Density split statistics: Cosmological constraints from counts and lensing in cells in DES Y1 and SDSS data. Phys. Rev. D 2018, 98, 023507. [Google Scholar] [CrossRef] [Green Version]
  31. Friedrich, O.; Gruen, D.; DeRose, J.; Kirk, D.; Krause, E.; McClintock, T.; Rykoff, E.; Seitz, S.; Wechsler, R.; Bernstein, G.; et al. Density split statistics: Joint model of counts and lensing in cells. Phys. Rev. D 2018, 98. [Google Scholar] [CrossRef] [Green Version]
  32. Hunter, J.D. Matplotlib: A 2D graphics environment. Comput. Sci. Eng. 2007, 9, 90–95. [Google Scholar] [CrossRef]
  33. Hinton, S. ChainConsumer. J. Open Source Softw. 2016, 1, 45. [Google Scholar] [CrossRef]
  34. Harris, C.R.; Millman, K.J.; van der Walt, S.J.; Gommers, R.; Virtanen, P.; Cournapeau, D.; Wieser, E.; Taylor, J.; Berg, S.; Smith, N.J.; et al. Array programming with NumPy. Nature 2020, 585, 357–362. [Google Scholar] [CrossRef] [PubMed]
  35. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The LDT model for the matter PDF (solid lines) compared to measured PDFs from the Quijote simulations (points). The LDT model remains more accurate than a log-normal approximation with the measured variance (dashed lines) on small scales and at late times.
Figure 1. The LDT model for the matter PDF (solid lines) compared to measured PDFs from the Quijote simulations (points). The LDT model remains more accurate than a log-normal approximation with the measured variance (dashed lines) on small scales and at late times.
Universe 08 00055 g001
Figure 2. Comparison of the matter PDF in 10 h 1 Mpc spheres in different theories of gravity. The cosmological parameters are chosen such that the clustering amplitude σ 8 is the same in all cases at redshift 0. This normalises the overall width of the PDF. The resulting difference in tilt and redshift dependence is due to the change to gravity.
Figure 2. Comparison of the matter PDF in 10 h 1 Mpc spheres in different theories of gravity. The cosmological parameters are chosen such that the clustering amplitude σ 8 is the same in all cases at redshift 0. This normalises the overall width of the PDF. The resulting difference in tilt and redshift dependence is due to the change to gravity.
Universe 08 00055 g002
Figure 3. Derivatives of the matter PDF in an f ( R ) modified gravity scenario. The f ( R ) gravity parameter, | f R 0 | , can be distinguished from σ 8 by both its redshift dependence and its additional skewness. Ω m can be disentangled from both by its different effect on skewness, and its effect on the linear growth factor.
Figure 3. Derivatives of the matter PDF in an f ( R ) modified gravity scenario. The f ( R ) gravity parameter, | f R 0 | , can be distinguished from σ 8 by both its redshift dependence and its additional skewness. Ω m can be disentangled from both by its different effect on skewness, and its effect on the linear growth factor.
Universe 08 00055 g003
Figure 4. Derivatives of the matter PDF in an evolving dark energy universe. The dependence of the matter PDF on Ω m is easily distinguished from the others by its distinct skewness (see Figure 3) and hence not shown here. The σ 8 , w 0 , and w a derivatives are similar in shape, but have different redshift evolutions, which allows for degeneracy breaking.
Figure 4. Derivatives of the matter PDF in an evolving dark energy universe. The dependence of the matter PDF on Ω m is easily distinguished from the others by its distinct skewness (see Figure 3) and hence not shown here. The σ 8 , w 0 , and w a derivatives are similar in shape, but have different redshift evolutions, which allows for degeneracy breaking.
Universe 08 00055 g004
Figure 5. Forecast constraints on f ( R ) gravity using a Euclid-like volume. These are marginalised over all other Λ CDM parameters, and include a prior on Ω b and n s described in [15].
Figure 5. Forecast constraints on f ( R ) gravity using a Euclid-like volume. These are marginalised over all other Λ CDM parameters, and include a prior on Ω b and n s described in [15].
Universe 08 00055 g005
Figure 6. Forecast constraints on w 0 w a CDM dark energy using a Euclid-like volume. These are marginalised over all other Λ CDM parameters, and include a prior on Ω b and n s described in [15].
Figure 6. Forecast constraints on w 0 w a CDM dark energy using a Euclid-like volume. These are marginalised over all other Λ CDM parameters, and include a prior on Ω b and n s described in [15].
Universe 08 00055 g006
Table 1. Detection significance for a fiducial f ( R ) with | f R 0 | = 10 6 and DGP model with Ω rc = 0.0625 . The stronger f ( R ) constraints are expected from the additional skewness in the PDF response to | f R 0 | as seen in Figure 3.
Table 1. Detection significance for a fiducial f ( R ) with | f R 0 | = 10 6 and DGP model with Ω rc = 0.0625 . The stronger f ( R ) constraints are expected from the additional skewness in the PDF response to | f R 0 | as seen in Figure 3.
  f ( R ) DetectionDGP Detection
PDF, 3 scales + prior 5.15 σ 1.17 σ
P ( k ) , k max = 0.2 h / Mpc + prior 2.01 σ 2.42 σ
PDF + P ( k ) + prior 13.40 σ 5.19 σ
Table 2. Constraints from mildly non-linear scales on σ 8 , w 0 , and w a as well as the dark energy Figure of Merit (FoM) coming from the matter PDF, power spectrum, and their combination.
Table 2. Constraints from mildly non-linear scales on σ 8 , w 0 , and w a as well as the dark energy Figure of Merit (FoM) coming from the matter PDF, power spectrum, and their combination.
  σ [ σ 8 ] σ 8 fid σ [ w 0 ] σ [ w a ] FoM
PDF, 3 scales + prior0.18%0.371.2527
P ( k ) , k max = 0.2 h / Mpc + prior0.45%0.241.0350
PDF + P ( k ) + prior0.17%0.090.40243
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gough, A.; Uhlemann, C. One-Point Statistics Matter in Extended Cosmologies. Universe 2022, 8, 55. https://doi.org/10.3390/universe8010055

AMA Style

Gough A, Uhlemann C. One-Point Statistics Matter in Extended Cosmologies. Universe. 2022; 8(1):55. https://doi.org/10.3390/universe8010055

Chicago/Turabian Style

Gough, Alex, and Cora Uhlemann. 2022. "One-Point Statistics Matter in Extended Cosmologies" Universe 8, no. 1: 55. https://doi.org/10.3390/universe8010055

APA Style

Gough, A., & Uhlemann, C. (2022). One-Point Statistics Matter in Extended Cosmologies. Universe, 8(1), 55. https://doi.org/10.3390/universe8010055

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop