Constraining cosmological parameters using void statistics from the SDSS survey
Abstract
Aims. Constraining the values of the amplitude of the linear spectrum of density fluctuations (), the matter density parameter (), the Hubble constant (H0), , where is the dark matter density parameter and = H, and S8 from the SDSS survey by studying the abundance of large voids in the large-scale structure of galaxies.
Methods. We identify voids as maximal non-overlapping spheres within the haloes of the Uchuu simulation and three smaller halo simulation boxes with smaller volume and different values, and galaxies with redshift in the range and absolute magnitude in the band of 32 Uchuu-SDSS simulated lightcones the seventh release of The Sloan Digital Sky Survey (SDSS DR7) survey. We compute the Void Probability Function (i.e. the probability that a randomly placed sphere with radius is empty of tracers) and the abundance of voids larger than predicted by the theoretical framework used and improved in this work and we check that it predicts successfully both void functions for the halo simulation boxes. Next, we asses the potential of this theoretical framework to constrain cosmological parameters using Uchuu-SDSS void statistics, and we calculate the confidence levels using Monte Carlo Markov Chain techniques to infer the values of , and H0 from the SDSS sample used.
Results. We have proved that using the four halo simulation boxes we successfully recover the values of , and from each box within () in real (redshift) space. We have also proved that the theoretical framework is really powerful using Uchuu-SDSS void statistics: if we fix one parameter to a constant value, the value given by Planck of the other two parameters is inside the 1 confidence level contour only if the fixed parameter is close enough to Planck’s value. Then, we have constrained these parameters from the SDSS survey sample used. The results are: , , H, and S8=1.017. If we combine these constraints with KiDS-1000+DESY3, we get , , H and S8=0.794. The combined uncertainties are approximately a factor 2-3 smaller than only-Weak-Lensing uncertainties. This is a consequence of the orientation of the confidence level contours of SDSS voids and Weak Lensing in the plane , which are almost orthogonal.
Key Words.:
cosmological parameters – large-scale structure of the Universe – cosmic voids – void-statistics – simulations1 Introduction
For over 40 years, it has been observed that the brightest galaxies are generally found in dense regions, and that most of cosmic space is devoid of this type of galaxies. This distribution is well understood as a natural evolution of density fluctuations in matter that were progressively amplified by gravitational instability (Giovanelli, 2010). In this framework, it can be demonstrated that initially low-density regions (voids) grow in size as highly dense areas collapse under their own gravity (Sheth & Van De Weygaert, 2004).
In the last century, there has been a significant focus on studying the over-dense regions of the Universe (e.g. Kiang & Saslaw, 1969; Bahcall, 1977; Kaiser, 1987; Einasto et al., 1994; Holder et al., 2001; Yang et al., 2005; Li et al., 2006; Gao & White, 2007; Zentner et al., 2019; Dong-Páez et al., 2024) , while the under-dense areas have only recently begun to receive adequate attention (e.g. Hoffman & Shaham, 1982; Zeldovich et al., 1982; Little & Weinberg, 1994; Goldwirth et al., 1995; van de Weygaert & Sheth, 2003; Li et al., 2012; Achitouv, 2019; Chan et al., 2019; Rodríguez-Medrano et al., 2024; Curtis et al., 2024). Studying these regions (voids) is very useful as they possess unique characteristics that make them important probes for cosmological studies and the physics of galaxy formation. They are useful, for example, for:
- •
- •
- •
- •
- •
In fact, in principle, any cosmological parameter could be constrained from void statistics. Specifically, the abundance of voids with radius larger than , , depends on the normalization of the linear spectrum of density fluctuations () and , and the shape of this function for large values of depends on , where =H ) (see Betancort-Rijo et al. (2009) for details).
Constraining and is especially interesting as there exist some statistically significant tension in cosmological parameter S between the Planck experiment (Aghanim et al., 2020), which measures the Cosmic Microwave Background, CMB, anisotropies (S), and other low-redshift cosmological probes, such as weak gravitational lensing, where a value of S is obtained by KiDS-1000 (Li et al., 2023). This tension is known as the S8 tension. (see Di Valentino et al. (2021) for more details).
Initial investigations into cosmic voids were constrained by the limited quantity of surveyed galaxies during that period. Nevertheless, with the emergence of more expansive redshift surveys, such as the Two-Degree Field Galaxy Redshift Survey (2dFGRS) (Colless et al., 2001, 2003), the Sloan Digital Sky Survey (SDSS) (York et al., 2000) and The SDSS’s Baryon Oscillation Spectroscopic Survey (BOSS) (Dawson et al., 2013), alongside improved resolution in cosmological simulations and enhanced analytical methodologies, we now have the capability to derive precise statistical insights concerning voids (e.g. Tikhonov, 2006; Ceccarelli et al., 2006; von Benda-Beckmann & Mueller, 2007; Tikhonov, 2007; Patiri et al., 2012; Hamaus et al., 2020; Douglass et al., 2023; Contarini et al., 2023).
However, despite the longstanding presence of the concept of cosmic voids, there exists no universally accepted definition for what constitutes a void. The term ”voids” can encompass disparate entities depending on the data employed and the objectives of the analysis. For example, voids can be defined as under-dense regions based on the smoothed dark matter (or halo/galaxy) density field (Colberg et al., 2005), as gravitationally expanding regions based on the dynamics of the dark matter density field (Hoffman et al., 2012), or as empty spatial regions among discrete tracers (El-Ad & Piran, 1997; Aikio & Mähönen, 1998; Hoyle & Vogeley, 2002).
The choice of a simple definition of voids, in particular, their definition as empty spheres, is convenient for statistical studies of galaxy voids. In this paper, we define voids as maximal non-overlapping spheres empty of objects with mass (or luminosity) above a given value. With this definition, it is clear that voids are not empty, as there can be low luminosity galaxies (or low mass haloes) inside them.
The aim of this paper is to inference the values of , H0 and from SDSS redshift survey using the theoretical framework of voids statistics developed in Betancort-Rijo (1990). This has been already done in previous articles in different ways. For example, in Sahlén et al. (2016) galaxy cluster and void abundances are combined using extreme-value statistics on a large cluster and a void. This way they obtain a value of for a flat universe. In Hamaus et al. (2016), they constrain studying void dynamics in SDSS. In Woodfinden et al. (2023) they also use SDSS survey to constrain measuring the void-galaxy and galaxy-galaxy clustering. In Contarini et al. (2023) they model Void Size Function (Press & Schechter, 1974; Sheth & Van De Weygaert, 2004) by means of an extension of the popular volume-conserving model (Jennings et al., 2013), based on two additional nuisance parameters, and applying a Bayesian analysis they get a value of for BOSS DR12, and, in a posterior work (Contarini et al., 2024), they constrain S and H for the same redshift survey.
This paper is structured as follows: in Section 2 we describe the redshift survey as well as the simulations used in this work. Next, in Section 3 we give a detail explanation of how our Void Finder works and show the relevant statistics of voids obtained for the observational catalogue and compare it with the result of Uchuu-SDSS lightcones. In Section 4 we introduce the most important concepts and equations used in the theoretical framework to calculate the abundance of voids larger than and the Void Probability Function. We show that the theoretical framework predicts successfully these two void statistics for halo simulation boxes with different values in real and redshift space, and show the results for SDSS and the Uchuu-SDSS lightcones. In Section 5 we explain the statistical test we use in order to infer , and . In Section 6 we show the potential of the theoretical framework using Uchuu-SDSS void statistics, i.e., we show that if we let two of the three cosmological parameters be free, the third parameter must be very close to Planck’s best-fit value in order to recover the real values of other two parameters. In this section, we also show the inferred values of , , , and S8 from Uchuu-SDSS. In Section 7 we show the constrained values for the sample of SDSS redshift survey we have used in this work and we combine our results with KiDS-1000 results (Dark Energy Survey and Kilo-Degree Survey Collaboration et al., 2023). In Section 8 we compare our constrained values of , and H0 with those obtained in other works where voids are also used. Finally, in Section 9 we summarize the most important results obtained in this work.
2 Data and mocks
The aim of this work is to infer , and H0 values from The Sloan Digital Sky Survey. However, to make sure that the theoretical equations reproduce correctly the void functions of this redshift survey, we also use 4 halo simulation boxes with different values (see Table 1), one simulation galaxy box and 32 lightcones. In this section we introduce the redshift survey and mocks.
2.1 SDSS
We have used the seventh release (Abazajian et al., 2009) of The Sloan Digital Sky Survey (SDSS DR7) (York et al., 2000), which includes 11663 of CCD imaging data in 5 photometric bands for 357 million distant objects. The catalogue also completed spectroscopy over 9380 . In total, there are 1.6 million spectra, including 930000 galaxies, 120000 quasars, and 460000 stars.
However, in this work we use a subcatalogue of SDSS: only galaxies from the northern regions with completeness greater than 90 are selected. Therefore, the effective area of this subcatalogue is 6511 and contains around 497000 galaxies with redshifts between . In this sample, 6 of targeted galaxies lack a spectroscopically measured redshift due to fibre collisions, so nearest neighbour correction is applied to these galaxies, assigning to them the redshift of the galaxy with which they collide (Dong-Páez et al., 2024).
Additionally, we impose further cuts on absolute magnitude and redshift, and we construct a volume-limited sample by keeping only galaxies brighter than the Milky Way-like galaxies (, where is the absolute magnitude in -band) and with and . The physical volume of this sample is Mpc3. There are in total 112496 galaxies which fulfil these restrictions, so the number density of the galaxies is Mpc3.
2.2 Mocks
Uchuu | P18(Low) | VeryLow | |
0.3089 | 0.3111 | 0.3111 | |
0.6911 | 0.6889 | 0.6889 | |
0.0486 | 0.048975 | 0.048975 | |
0.6774 | 0.6766 | 0.6766 | |
0.8159 | 0.8102 (0.75) | 0.65 | |
[Mpc] | 2000 | 1000 | 1000 |
kpc | 4.27 | 1.0 | 4.0 |
Haloes | Yes | Yes | Yes |
Galaxies | Yes | No | No |
Lightcones | Yes (32) | No | No |
We use the Uchuu simulation, which was produced using the TreePM code GreeM 111http://hpc.imit.chiba-u.jp/~ishiymtm/greem/ (Ishiyama et al., 2009, 2012) on the supercomputer ATERUI II at Center for Computational Astrophysics, CfCA, of National Astronomical Observatory of Japan. The number of dark matter particles was 128003 with a mass resolution of in a box with a side length of 2000 Mpc. A total of 50 halo catalogs (snapshots) were created, covering the redshift range from 0 to 14 (Ishiyama et al., 2021). The haloes were subsequently identified using the RockStar halo/subhalo finder 222https://bitbucket.org/gfcstanford/rockstar/ (Behroozi et al., 2013a), and merger trees were generated using the consistent trees code333https://bitbucket.org/pbehroozi/consistent-trees/ (Behroozi et al., 2013b). These simulations adopted the cosmological parameters from Planck 2015 (see Table 1) (Planck Collaboration et al., 2016). All this data is publicly available and accessible in the Skies Universes database444http://www.skiesanduniverses.org/Simulations/Uchuu/, including galaxy catalogs constructed using various methods (Aung et al., 2023; Oogi et al., 2023; Gkogkou et al., 2023; Prada et al., 2023; Ereza et al., 2023; Dong-Páez et al., 2024).
Apart from Uchuu, we use three more simulations with Planck 2018 (Aghanim et al., 2020) (Uchuu1000Pl18), Planck 2018 with (Uchuu1000Pl18LowS8) and Planck 2018 with (Uchuu1000Pl18VeryLowS8). Uchuu1000Pl18 and Uchuu1000Pl18LowS8 have exactly the same simulation properties, except the value of . Uchuu1000Pl18VeryLowS8 has different value of and different mass resolution. All these details can be seen in Table 1. These three simulations were produced using the TreePM code GreeM on the supercomputer Fugaku at the RIKEN Center for Computational Science. As well as the Uchuu simulation, we generated initial conditions for these three simulations with 2nd order Lagrangian perturbation theory (Crocce et al., 2006) by 2LPTIC code 555http://cosmo.nyu.edu/roman/2LPT/. The initial conditions are identical across the three simulations, enabling us to compare them without cosmic variance. The total of 70 halo catalogs were created, where the redshift list is the same with the Uchuu at . N-body data including the halo catalogs and the merger trees are also available in the Skies Universes database.
The Halo Mass Function (HMF) of each simulation can be seen in Figure 1, where is the halo mass within virial radius including unbound particles. In this figure it can be seen that VeryLow has lower mass resolution, as the HMF decreases considerably for small masses. It is also shown that, specially for large masses, the lower is the lower the HMF is too. Additonally, it can be seen that there is no significant difference between P18, Low and VeryLow halo mass functions for low virial masses (), but there is a significant difference for large virial masses ().
We also use a galaxy simulation box that has been generated from Uchuu simulation using the SubHalo Abundance Matching (SHAM) method. From the galaxy boxes for different snapshots, lightcones with the properties of SDSS have been constructed. In concrete, 32 light cones are available. The Uchuu galaxy box and the 32 simulated SDSS light cones we use in this work are extensively detailed in Dong-Páez et al. (2024).
For the purpose of studying void statistics in simulation boxes, we select those snapshots corresponding to a redshift of (which is the snapshot with the closest redshift to SDSS median redshift). Moreover, we don’t use all the objects (haloes or galaxies) from this boxes, but we select a number density of Mpc3 for each box. This can be reached removing haloes less massive than = 1.626 (Uchuu), 1.642 (Uchuu1000Pl18), 1.595 (Uchuu1000Pl18LowS8) and 1.457 (Uchuu1000Pl18VeryLowS8), and galaxies less brighter than ¡-20.5 for Uchuu galaxy box.
3 Statistics of voids in the SDSS
The next important step we need to take once we have well-defined the sample of halos or galaxies we are going to use is to identify voids, defined as maximum non-overlapping spheres. To do this, we perform Delaunay triangulation (with periodic conditions for simulation boxes).
The Delaunay triangulation of a set of points ensures that no point lies within the circumcircle of any triangle in the triangulation. Additionally, Delaunay triangulations maximize the minimum angle among all triangle angles, aiming to reduce the occurrence of sliver triangles.
Fortunately, there is a publicly available code that uses the CGAL (The Computational Geometry Algorithms Library)666https://doc.cgal.org/4.6.3/Manual/packages.html#PkgTriangulation3Summary for 3D triangulations. This code is called Delaunay trIangulation Void findEr (DIVE777https://github.com/cheng-zhao/DIVE), which is used and explained in Zhao et al. (2016). The output of this code is a file that contains the positions in space of the centres of the spheres that satisfy the Delaunay condition, as well as their radii. However, these spheres are not voids, but candidates to be voids.
To find voids (i.e. maximal non-overlapping spheres) among these spheres, an additional code needs to be developed. This code must check if two spheres overlap and, in case they do, keep the largest one as a void.
In Figure 2 the voids larger than Mpc found in the P18 halo simulation box with number density Mpc3 are shown in a region of the box. It can be seen that each void is defined by four galaxies laying in its surface (some voids in Figure 2 have fewer than 4 galaxies due to being outside the plotted region), and it is ensured that voids do not overlap.
However, for the galaxy lightcones, one more step is required. This step involves considering the incompleteness in stellar mass, which causes many spurious voids to appear. These spurious voids correspond to regions in the catalog with very low completeness, where galaxies could not be detected properly. Therefore, these false voids must be eliminated using an angular mask, such as the Healpix map (Gorski et al., 1999; Blanton et al., 2005; Swanson et al., 2008), characteristic of the redshift survey.
The procedure for removing these false voids is as follows: first, we generate points uniformly distributed within the volume of the void. Next, we project these points in the angular plane, calculate the completeness of each point and average the completeness of all points within the void to obtain an approximation of the completeness of the void. If this completeness is equal to or greater than 0.9, we label that void as a true void. Otherwise, we remove the false void.
Additionally, voids whose centers or part of their volume are outside the sample in the radial direction must also be removed, as the procedure explained above considers only the angular plane, but voids may extend beyond the sample in the radial direction, so this extra check is necessary.
r | SDSS | Uchuu-SDSS |
---|---|---|
10 | 797 12 | 792.2 2.1 |
11 | 550 15 | 550 3 |
12 | 381 13 | 366.3 2.3 |
13 | 248 13 | 233.4 2.2 |
14 | 150 11 | 140.6 1.9 |
15 | 85 8 | 78.9 1.4 |
16 | 51 7 | 41.4 1.2 |
17 | 27 5 | 21.2 0.8 |
18 | 15 3 | 9.3 0.5 |
19 | 7.0 1.7 | 3.4 0.3 |
20 | 1.0 1.0 | 0.97 0.18 |
In Figure 3 the number density of voids larger than , , obtained for voids found in the sample of SDSS and Uchuu-SDSS (the mean of the 32 lightcones) considered in this work (galaxies with and ) are shown (in Table 2 the values are multiplied by the volume). It can be checked that Uchuu-SDSS statistics are compatible with observations for all radius bins, although there are big fluctuations in SDSS that appear for large voids (Mpc) due to the small volume of the survey.
An important remark about Figure 3 is that if we want to compare the observed number density of voids larger than obtained from SDSS (or Uchuu-SDSS) with that given by the theoretical framework presented in this work, we must take into account that the former suffers from incompleteness, and other effects such as border effects, and the latter doesn’t. Therefore, we have to transform SDSS (and Uchuu-SDSS) void statistics as if it didn’t suffer from these effects. One way of doing this is using Uchuu galaxy box. Then, SDSS and Uchuu-SDSS number density of voids larger than must be transformed as:
(1) |
where is the number of voids larger than found in Uchuu-SDSS lightcones, and is the same but found in Uchuu galaxy box, and is the volume of Uchuu simulation (V=Mpc-3).
Finally, in Figure 4 we can see the VPF, (i.e. the probability that a randomly placed sphere with radius is empty of objects – galaxies or dark matter haloes) obtained for SDSS and Uchuu-SDSS galaxies. This function is challenging to calculate as we have to take into account the completeness of the survey when randomly placing spheres. The procedure we have followed is similar to the one followed to find real voids: first, we randomly place spheres in the angular plane and in the line-of-sight direction, then we generate points uniformly distributed within the volume of these random spheres and check if the mean completeness of all these points is above 0.9. If it is, we keep the sphere, otherwise, we generate a new one. In other words, our random spheres have to fulfill the same criteria as voids. Next, we calculate the VPF with these survivor spheres, i.e. we divide the total number of these spheres containing no galaxies in their interior by the total number of the survivor spheres.
The results of the VPF obtained for SDSS and Uchuu-SDSS can be seen in Figure 4 and Table 3. Again, we can see that the values obtained for SDSS are compatible, within , with Uchuu-SDSS values, although the errors for large are very large. This is again due to the small volume of the samples used.
r | SDSS | Uchuu-SDSS |
---|---|---|
10 | 6.127 0.389 | 5.640 0.073 |
11 | 3.354 0.258 | 3.049 0.047 |
12 | 1.728 0.162 | 1.567 0.032 |
13 | 0.835 0.101 | 0.767 0.204 |
14 | 0.390 0.069 | 0.356 0.014 |
15 | 0.161 0.045 | 0.160 0.008 |
16 | 0.0631 0.0242 | 0.0694 0.0043 |
17 | 0.0216 0.0028 | 0.0270 0.0035 |
18 | (8.8602.018) | 0.0107 0.0025 |
19 | (3.8809.232) | (4.161.18) |
20 | 0.0000 0.0005 | (1.5560.699) |
4 Void statistics theoretical framework
The main void statistics we will study in this work are the number density of voids larger than , 888Notice that we use a bar above to quantities predicted by theoretical framework, and without bar to quantities obtained directly by simulations. Additionally, quantities in lowercase letters are in units of volume., and the Void Probability Function, (White, 1979). For we will use the (recalibrated) expression from Betancort-Rijo et al. (2009):
(2) |
where =0.588, =1.671, and
(3) |
and
(4) |
In expression 2 it is assumed that is an universal functional of , so that, independently of the clustering process underlying , is related to it by an unique expression. However, it has been shown that for white noise it . So it is clear that this coefficient has a dependence on the clustering properties of the objects considered. But we have found that for the simulations that we have used the quoted values of and are valid. Therefore we shall use these values in all our considerations.
It is important to remark that equation (2) is only valid for 0.46. For 0.46, . This quantity, , measures the rareness of the voids.
Furthermore, equation (2) is only valid for , as it has been calibrated to correctly predict the number density of voids for , as well as for large enough (i.e. large enough voids, such as Mpc).
It can be seen that equation (2) differs from the expression of given in Betancort-Rijo et al. (2009) (i.e., they give different values of , and write an additional term in the exponential, ). Instead of writing the same values for , and as in Betancort-Rijo et al. (2009), we have let them be variables and calibrated them with the Uchuu simulation, which has a much larger volume than the simulation used in that work. Therefore, the coefficients given in this work are more accurate.
From equations (2) and (3) we can see that if we know the VPF, we automatically know , so we first study the VPF and then calculate using equation (2).
We can define from a more general statistic that is : the probability that a sphere of radius , placed at random within the distribution, contains objects. If we assume a Poisson process, we can then write as (Layzer, 1954):
(5) |
where is the probability distribution for the integral of the probability density, , within a randomly placed sphere.
For the root mean square (rms) of the error of the estimation of the VPF we use the following expression:
(6) |
where is the number of random spheres used in order to estimate for simulations, is the number of voids larger than . This equation is a modification of the expression given in Betancort-Rijo (1992) and Patiri et al. (2006b). We have added a second term, to take into account the error due to the finite number of trial random spheres used to calculate the VPF. The first term of equation (6) takes into account the finite volume of the sample.
Finally in equation (6) is given by
(7) |
where is the mean radius of all voids larger than the minimum voids considered.
An explicit computation of each term in can be consulted in Appendix A. There, it can be seen that depends on (equation (33)), (equations (27), (35) and (36)) and (equations (30) and (38)). However, the dependence on is small, as we will show in Section 5.
Additionally, it can be checked that (in concrete, ), depends on , which is the mass such that the number density of distinct or central haloes with mass larger than is equal to the number density of the sample, , (i.e. the mass such that the cumulative halo mass function equals the number density of the sample). If the simulation boxes consisted only on distinct haloes, the mass to be put in would be simply the mass obtained this way. However, our simulation boxes contain subhaloes, too, and we even have a galaxy simulation box, galaxy lightcones and, more importantly, a redshift survey of galaxies, so the mass to be put in the theoretical equations is not that trivial.
If our simulation boxes consist on distinct haloes and subhaloes (or galaxies), then the mass to be used in the theoretical equations would be , i.e., the mass such that the cumulative total halo mass function (taking into account distinct haloes and subhalos) equals the number density of the sample. This mass can be approximately related to through the following equation:
(8) |
where is the rms linear density fluctuation on scale . It can be shown that, in the interval of masses and cosmological parameters considered in this work, depends on , and the following way:
(9) |
Additionally, an expression for depending on the , and can be found:
(10) |
where
(11) |
This last function, defines the dependence of the mass with the number density of the sample, . It is important to remark that is the number density of galaxies corrected from completeness in the case we use redshift surveys or light-cones.
Equation (10) is a good approximation to (i.e. the value of to be used in the theoretical framework, equations 25-28 ) as a function of cosmological parameters and , while the abundance matching provides it exact value. However, the dependence on and is quite small and it is swamped by small error in expressions 25-28 (which have been fitted to the results of a complex procedure) that have a much stronger dependence on and . So we have found that the best agreement between the theoretical framework and the simulations are obtained holding fixed for given values of .
The value of has been chosen taking into account the value predicted for Uchuu, but decreasing its value a little bit so the value for the three small boxes, defined using the VPF, is between the range
(12) |
where is the number of degrees of freedom. In order to calculate the chi square function, we have used the VPF predicted by the theoretical framework developed in this work, and the VPF obtained directly for each simulation box. We have used the radius bins with Mpc, Mpc and Mpc (i.e., 8 bins, so ). Therefore, must be between and . We can see that if we choose , the value of for Uchuu, P18, LowS8 and VeryLows8 are, respectively, 11.20, 6.37, 11.86 and 11.74. If we increase the value of to 4.67, then for LowS8 and VeryLowS8 boxes, and if we decrease the value of to 4.65, then for Uchuu, which is our high resolution simulation. Therefore, is the safer choice.
We can see that so far our theoretical framework depends on 4 cosmological parameters: , , H0 and (it doesn’t depend on as we impose a flat CDM model, and the dependence of the formalism with and H0 is only through ). However, we can decrease the number of free parameters by relating and , where is the baryonic matter density parameter. From Planck 2018 (Planck Collaboration et al., 2016) we get . Therefore, and , and our final parameters are , and .
Before using the theoretical framework to constrain and and H0 in the SDSS survey, it is important to check if the theoretical equations are accurate and have enough precision to recover the real values of , and H0 from simulations with different values of the parameters. In order to do this, we use the halo simulation boxes that we have presented in Section 2.2, which have different values of .
In the upper panel of Figure 5 the VPF values (multiplied by to differentiate more easily the values of each simulation for large radii) predicted by the theoretical framework (continuous lines) and obtained by simulations (points) for the four halo simulation boxes in real space are shown. The numerical values are presented in Table 7 in Appendix B. It can be seen that the dependence of the VPF predicted by the theoretical framework is the same as that shown by simulations, that is, the lower is, the lower is for each radius bin. In the bottom panel of the same Figure, the ratio between simulations and theoretical framework is shown. It can be checked that all values are within (or compatible with) 10 of the ratio.
In the upper panel of Figure 6 the values of the number density of voids larger than predicted by the theoretical framework (continuous lines) and obtained by simulations (points) for the 4 halo simulation boxes in real space can be seen. In the bottom panel of the same Figure, the ratio between simulations and theoretical framework is shown. Again, the agreement between the theoretical framework and the simulation values is good, especially for between 12 and 18 Mpc.
It is important to note that the values of the VPF and for the three small box catalogs (each with a different value of ) are very similar for small values of (Mpc). The differences between the three simulations begin to become important for Mpc, approximately. This is why these are the only voids we consider in order to constrain in the next section.
The next step is to repeat this process in redshift space and compare with simulations in order to check if the theoretical framework still works in this space. This is done in Appendix B.2. In concrete, in Figures 15 and 16 the VPF and number density of voids larger than in redshift space are shown for the four halo simulation boxes. The agreement is still good, so we can conclude that the theoretical framework works in redshift space, too.
Now, we can discuss the results obtained for SDSS survey and Uchuu-SDSS lightcones. In Figure 3 we have already shown the abundance of voids larger than for SDSS and Uchuu-SDSS. The values given by the theoretical framework are shown in the same figure with a continuous black line. We can see that the agreement between the theoretical framework and Uchuu-SDSS is good for Mpc.
Additionally, in Figure 4 we have already shown the VPF for SDSS and Uchuu-SDSS. The values given by the theoretical framework are shown in this Figure, too, with a continuous black line. We can see that the agreement between the theoretical framework and Uchuu-SDSS is quite good, and is compatible within with SDSS values.
To sum up, we have checked in this section that the theoretical framework successfully predicts the void statistics studied in this work for the four halo simulation boxes and for Uchuu-SDSS lightcones, and it is compatible with SDSS survey statistics within . Therefore, we are ready to constrain and and H0 from the SDSS galaxy sample we have chosen.
5 Bayesian analysis for cosmological parameters inference
In order to infer , and H0 from the number density of voids larger than , we use a Bayesian analysis and Markov chain Monte Carlo (MCMC) technique in order to sample the posterior distribution of the considered parameter sets, :
(13) |
where is the prior distribution, and is the likelihood, calculated as
(14) |
where is the number of voids within th bin from simulations or data, is calculated from equation (2) subtracting the value in the th bin (or ), and is the rms of :
(15) |
where represents the number of realizations for each simulation. For the simulation boxes and SDSS, is equal to 1, indicating a single realization for each set of cosmological parameters and one SDSS survey. However, in the case of Uchuu-SDSS, is 32, as we have 32 lightcones of this particular type.
In this work, we use voids larger than Mpc to constrain the cosmological parameters of interest, so in equation (14) we consider the radius bins starting from Mpc to Mpc (i.e. 6 bins with width = 1Mpc).
With equation (13) the best estimate of the parameters may be obtained by maximising with respect to the parameters . In order to do this, we assign to all the parameters (, and H0) wide enough priors (see first column of Table 4).
Voids-Only | Voids+KiDS | Voids+DES | |
---|---|---|---|
[0.5-1.4] | [0.5-2.2] | [0.5-2.5] | |
[0.15-0.5] | [0.15-0.686] | [0.15-0.9] | |
H0 | [40-140] | [64-82] | [55-91] |
Additionally, we may want to combine our results with other works. To do this, we multiply our likelihood by the likelihood of the work we want to combine our results with. In this work, we combine our results with Planck 2018 (Aghanim et al., 2020) and some Weak Lensing works, such as KiDS-1000 and DESY3 (Dark Energy Survey and Kilo-Degree Survey Collaboration et al., 2023). We use Cobaya (Torrado & Lewis, 2019, 2021) in order to run MCMC sampler for SDSS voids, so Planck’s likelihood can easily be combined with SDSS voids likelihood. However, the likelihoods of the weak lensing works used in this study are not incorporated in Cobaya, and CosmoSIS (Zuntz et al., 2015) must be used instead. However, CosmoSIS is very slow for our purpose, so we have used a brand new code called CombineHarvesterFlow 111https://github.com/pltaylor16/CombineHarvesterFlow which allows one to efficiently sample the joint posterior of two non-covariant experiments with a large set of nuisance parameters (Taylor et al., 2024). In concrete, CombineHarvesterFlow trains noralizing flows on posterior samples to learn the marginal density of the shared parameters. Then by weighting one chain by the density of the second flow, one can find joint constraints.
The confidence level contours in the plane -H0 (fixing in each case to its real value, see Table 1) obtained for the four halo simulation boxes can be seen in the Appendix B. There, we can check that we recover the values of and H0 of the simulations within () in real (redshift) space. Therefore, we can infer these cosmological parameters in Uchuu-SDSS lightcones.
6 Assesing the potential of voids statistics to constrain cosmological parameters (using Uchuu-SDSS)
Before constraining any cosmological parameter, we will see how does the theoretical framework behave when we fix one of the three cosmological parameters on which the theoretical framework depends on and we let the other two parameters vary in a large range. We will do this for Uchuu-SDSS lightcones in order to check if the theoretical framework is able to recover successfully the cosmological parameters of Planck 2015 (Planck Collaboration et al., 2016), which are the real values of Uchuu-SDSS cosmological parameters.
In order to do this, we will plot the confidence level contours obtained as explained above in each of the three possible planes (, -H0 and -H0) fixing the remaining parameter to three different constant values to check if the formalism is sensitive to this parameter.
Figure 7 shows the contours of 68 and 95 confidence level in the planes (first row), -H0 (second row) and -H0 (third row) planes. For each plane, the value of the remaining parameter (H0, and for first, second and third row, respectively) has been fixed at different values which are indicated in the boxes inside each subplot. It can be seen that Planck 2015 values are inside the 68 confidence level region only when the remaining parameter is fixed to a value close to Planck 2015 parameter. For example, if we focus on contour, we can see that the Planck values of and are not compatible with Uchuu-SDSS statistics predicted by our theoretical framework if H or H, but it is compatible with these statistics if H, which is the closest value to the one of Planck 2015 (see Table 1). The same effect can be seen in the rest of the planes.
As we mentioned in Section 4, the theoretical framework mainly depends on and , although it also depends on , but this dependence is small. We can see this in the last row of Figure 7, where the contours obtained for plane are shown for different values of which are indicated inside the box of each figure. It can be seen that the contours slightly displaces with different values of , but this displacement is not as noticeable as it is in and -H0 planes. Therefore, we can conclude that the theoretical framework depends mainly on and parameters.
The last important the contour to study is the one obtained in S plane. This contour is shown in Figure 8 for different values of H0. It can be seen that H (which is the closest value to the real one of Uchuu simulations) is compatible with Uchuu-SDSS within 1, but H and H are not.
The last step we can take in order to asses the potential of the theoretical framework to constrain cosmological parameters is inferring , and H0 values. In order to do this, we use Cobaya and the priors indicated in the first column of Table 4. The inferred values of these cosmological parameters provided by Cobaya are: , , H, and S. Our uncertainties are very wide (in comparison of Planck’s uncertainties), specially the uncertainty of H0. With these wide uncertainties, we can check that our constrained values are compatible with Planck values within 1.
The confidence level contours for Uchuu-SDSS can be seen in the right part of Figure 9. We can see that the 68 confidence level contour closes for high values of , although it doesn’t for low values as we haven’t explored regions with because, as mentioned before, many of the equations of the formalism are not valid for .
7 Cosmological Constraints from SDSS survey
In this part of the work, we show the constraints obtained for SDSS redshift survey. Firstly, we show the constraints we obtain directly from the theoretical framework developed in this work. Next, we combine our results with Planck 2018 and, finally, we do it with Weak Lensing works such as KiDS-1000, DESY3 and the combination of both works, KiDS-1000+DESY3.
7.1 SDSS voids-only constraints
The initial contour we examine lies within the plane showcasing the core parameters of our theoretical framework: the - plane, as depicted in Figure 10. In this illustration, we have also included the confidence level contours for Uchuu-SDSS and Planck 2018. We can see that Planck’s contour is completely inside Uchuu-SDSS and SDSS voids contours, so we don’t except to get a big improvement in the constrained values of the parameters when combining. Therefore, in order not to waste computational resources, we have not made such a combination.
From Figure 10 it can be seen that the confidence level region for Planck 2018 is entirely encapsulated within the SDSS contour at the region of 2, indicating that the constraints from both samples are statistically compatible in this limit.
In the left part of Figure 9 we can observe the confidence level contours obtained in the rest of the planes: -H0, -H0 and .
In the first column of Table 5 we can observe the constrained values directly obtained from our theoretical framework of , , H0, and S8 for SDSS survey. From this Table, it can be seen that the uncertainties of and are much larger for SDSS than for Uchuu-SDSS because of the huge difference in the volumes between them. However, this is not the case for and H0. This is because and are the fundamental parameters of the theoretical framework (and S8 depends strongly on ). With these large errors, our constraints from SDSS are compatible with Planck’s within 1.
As we have already mentioned, we are not going to combine our results with Planck 2018, however we can show the behaviour of the theoretical framework for SDSS when two of the three cosmological parameters adopt a value close enough (or a range centered in the best-fit value of Planck 2018 and with an amplitude of 3, for example) to Planck values. These constraints can be seen in the second column of Table 5. We can see that the uncertainties decrease considerably for all parameters but S8. However, any parameter obtained this way is compatible within with Planck 2018. Nevertheless, it is important to remark the low mean values we get for and H0 the formalism predicts even if we restrict to an interval of 3s of Planck for the other two parameters.
SDSS voids | SDSS voids (3P) | Planck | |
---|---|---|---|
1.028 | 1.003 | 0.81110.0060 | |
0.296 | 0.265 | 0.31530.0073 | |
H0 | 83.43 | 58.86 | 67.360.54 |
0.1947 | 0.1524 | 0.17720.0027 | |
S8 | 1.017 | 1.061 | 0.8320.013 |
7.2 SDSS voids + Weak Lensing
SDSS voids | KiDS-1000 | DESY3 | SDSS+KiDS | SDSS+DESY3 | SDSS+KiDS+DESY3 | Planck | |
---|---|---|---|---|---|---|---|
1.028 | 0.833 | 0.816 | 0.859 | 0.881 | 0.858 | 0.81110.0060 | |
0.296 | 0.270 | 0.297 | 0.247 | 0.256 | 0.257 | 0.31530.0073 | |
H0 | 83.43 | 75.42 | 73.834.98 | 74.20 | 74.72 | 74.17 | 67.360.54 |
S8 | 1.017 | 0.764 | 0.802 | 0.775 | 0.811 | 0.794 | 0.8320.013 |
Finally, we can combine our results with KiDS-1000 and DESY3. These results are presented in Dark Energy Survey and Kilo-Degree Survey Collaboration et al. (2023).
When combining our results with other studies, such as those mentioned above, caution must be taken, as the same range for all the parameters must be considered if we want to visually compare the confidence level contours obtained in each case and combine the chains using CombineHarvesterFlow. This allow us to make a fair comparison of our confidence level contours with those of Weak Lensing. However, as we have seen previously, the theoretical framework used in this work depends on the parameters sigma8, omegamatter, and H0. In Weak Lensing studies, As (DESY3) or S8 (KiDS and KiDS-1000+DESY3) are sampled instead of , and and instead of (KiDS and KiDS-1000+DESY3). Therefore, what we have done is to determine which priors on and correspond to the aforementioned priors, ensuring compatibility across these parameters. Nevertheless, it’s important to consider a significant limitation: our formalism is not valid for values less than 0.5, as mentioned earlier, nor for values less than 0.15. Thus, we have added this additional condition in the priors. This may slightly affect the process of combining chains using CombineHarvesterFlow, but, as we will see later, the effect will be very small because of the relative orientation and size of SDSS voids and Weak Lensing confidence level contours.
The ranges we have used for our chains when combining with KiDS-1000 and DESY3 can be seen in the second and third columns, respectively, of Table 4. We can see that the range of H0 used in these works is very narrow.
The confidence level contour in the plane using these ranges can be seen in Figure 11. In this Figure, the confidence level contour for SDSS voids with the same ranges as the different weak lensing works considered are shown, as well as each weak lensing work confidence level contour and the combination of SDSS voids with these three contours. We can see that SDSS voids contour is almost orthogonal to the three Weak Lensing contours (as expected, see Contarini et al. (2023) for more details), so we can anticipate that our constrained values of and will have smaller uncertainties than the values of the original work.
The best-fit values of , , H0 and S8 obtained from SDSS voids, KiDS-1000 and DESY3 can be seen in the first three columns of the Table 6. The combination of SDSS voids with KiDS-1000, DESY3 and KiDS-1000+DESY3 can be seen in the fourth, fifth and sixth columns, respectively, of the same Table.222It is important to remark that CombineHarvesterFlow gives the inferred parameters obtained from the combination for two chains in two different ways: weighting SDSS voids chains or weighting Weak Lensing chains. The results we have presented correspond to weighting SDSS void chains, but we have checked that weighting Weak Lensing chains we obtain compatible results. We can see that the effect of combining these three Weak Lensing works with SDSS voids is increasing the value of , as the best-fit value of SDSS voids is very high, which means implies a decrease in the best-fit value of (from Figure 11 we can see that the three Weak Lensing works predicts a strong correlation between and , and that if one increases, the other decreases. This correlation is kept also when combining Weak Lensing works with SDSS voids).
We can also see from table 6 that there is an increase in the uncertainties of H0 when combining Weak Lensing with SDSS voids, which is caused by our huge uncertainties in this parameter. However, the uncertainties of the rest of the parameters are decreased by a factor 2-3, approximately, with respect to the original errors of each Weak Lensing work.
As a final remark, we can see from table 6 that any value of H0 or obtained when combining SDSS voids with the three Weak Lensing works are compatible with Planck 2018 within 1 (for , they are within 2.3, approximately). We can also see that S8 value from SDSS voids + DESY3 is compatible with Planck 2018 within 1, but SDSS voids + KiDS-1000 and SDSS voids + KiDS-1000+DESY3 aren’t (they are compatible with Planck 2018 considering the uncertainty as 1.8 and 1.6, respectively).
8 Comparison with other works about voids
In this part of the work we compare our constraints in , H0 and S with those obtained in other works that also use voids statistics to constrain these same parameters.
The first work we can compare our results with is Contarini et al. (2024), where the values of S8 and H0 are constrained using the void size function (abundance of voids with radius ) predicted by the excursion set theoretical framework (Press & Schechter, 1974; Sheth & Van De Weygaert, 2004) by means of an extension of the popular volume-conserving model (Vdn model, Jennings et al., 2013). The constraints they obtain combining the voids counts with the void shapes (Hamaus et al., 2020) are =0.809, =0.308, H0=67.3 and S8=0.813. The redshift survey they use in order to constrain these cosmological parameters is BOSS DR12 (Dawson et al., 2013). In concrete, they use LOWZ and CMASS target selections, and divide the catalogs into two redshift bins: and . This sample has a physical volume approximately 60 times larger than the volume of the sample used in this work (according to Reid et al. (2015) the volume of CMASS is 5.1Gpc3 and the one of LOWZ is 2.3Gpc3, which translates into a total volume of 2.3109 Mpc3). If we take into account that errors (using only void statistics, without combining with Weak Lensing or other works) scales as with volume (see appendix C), we can calculate how much errors would decrease if we would use a redshift survey like BOSS DR12 with our theoretical framework: , i.e. our errors (calculated as 68 uncertainties) would decrease approximately 3 times, so we would have similar constraints with our theoretical framework (without combining with any other work) than with Vdn model+Void Shapes if we used a sample with BOSS DR12 volume.
The other work we can compare our results with is Sahlén et al. (2016). In this work, galaxy cluster and void abundances are combined using extreme-value statistics on a large cluster and a void aligned with the Cold Spot (CS) in the CMB (CS void) (Finelli et al., 2016). This way they obtain a value of for a flat universe. This constraint is also compatible with CMB value as well as Weak Lensing values.
9 Summary
In this work, we have made use of the theoretical framework developed in Betancort-Rijo et al. (2009) and recalibrated the expression for the number density of voids larger than (see equation (2)) using the Uchuu halo simulation box with a number density of haloes equal to Mpc3, which has a larger volume than the simulation used in the aforementioned work (Mpc3). The most important results obtained in this work are the following:
-
•
we have proved that the number density of voids larger than and the Void Probability Function (VPF, i.e., the probability that a randomly placed sphere with radius is empty of galaxies or haloes) of SDSS galaxies with (where is the absolute magnitude in band), and , and Uchuu-SDSS galaxies with the same characteristics, are compatible within ,
-
•
we have demonstrated that our theoretical framework predicts successfully the Void Probability Function and the abundance of large voids for the four halo simulation boxes with different values of () and Uchuu-SDSS lightcones used in this work, again within ,
-
•
we have used a bayesian analysis for all halo simulation boxes and we have calculated the contours for 68 () and 95 () confidence levels in -H0 plane (fixing in each case to its real value, see Table 1), and proved that we recover the real values of and H0 of each halo simulation box within in real space, and in redshift space,
-
•
supposing that the ratio of is constant and given by Planck 2018, we have studied, how much Uchuu-SDSS contours change in each plane mentioned above if we fix the values of the remaining parameter to three sufficiently different values. We have seen that these contours are very sensitive to these changes, so if we considered a survey with the same volume as Uchuu-SDSS we could constrain with high precision each parameter.
-
•
we have also shown that the dependency of the theoretical framework on is small and that the main dependence of is through and . We have don this studying how does the contour change in the plane if we fix to three different values, and we have seen that the contour does not change much.
-
•
we have checked that using MCMC sampler from Cobaya, we successfully recover the values of these parameters for Uchuu-SDSS lightcones. The recovered values are , , H, and S.
-
•
We have calculated, then, the contours in , -H0 and H planes for the same confidence levels as mentioned above for SDSS void statistics. The contours obtained for SDSS are much wider than the ones obtained for halo simulation boxes (and Uchuu-SDSS lightcones) because of the huge difference of volume (SDSS volume is approximately 25 times smaller than the three small simulation boxes and 200 times smaller than Uchuu box volume), which means that there are far fewer voids in SDSS than in the boxes and much higher errors. The constrained values we get in this work for SDSS are , , H, and S8=1.017. It is important to remark that these constrained values have been obtained supposing that the ratio of is constant and given by Planck 2018,
-
•
Next, we have combined our SDSS voids constraints with KiDS-1000, DESY3 and KiDS-1000+DESY3 (Dark Energy Survey and Kilo-Degree Survey Collaboration et al., 2023). The results obtained when combining SDSS voids with KiDS-1000+DESY3 are: , , H and S8=0.794. We can see that when we combine our results with Weak Lensing works, we improve considerably the uncertainties of all the parameters, but H0 with respect to the uncertainties of only KiDs-1000+DESY3. This is because the confidence level contours in the plane obtained in this work and the one from KiDS-1000 are almost orthogonal.
-
•
We have seen that any value of H0 or obtained when combining SDSS voids with the three Weak Lensing works are comaptibel with Planck 2018 within 1 (for , they are within 2, approximately). We have also seen that S8 value from SDSS voids + DESY3 is compatibel with Planck within 1, but SDSS voids + KiDS-1000 and SDSS voids + KiDS-1000 + DESY3 aren’t (the former is compatible within 2.7 and the latter within 2).
-
•
Finally, we have compared our results with those obtained in other works where voids are also used in order to constrain the same cosmological parameters that we have constrained in this work, and with parameter constrained via CMB, galaxy clustering, weak lensing and Type Ia Supernova measures. We can check that when combining with Weak Lensing works, we obtain slightly smaller uncertainties than in other works where voids are also used, and when combining with Planck, we obtain much smaller uncertainties. However, if we don’t combine our results with any work, our uncertainties are very large because the volume of the samples of Uchuu-SDSS and SDSS used in this work is relatively small in comparison of the volumes of redshift surveys used in other works of voids. For example, in Contarini et al. (2024) BOSS redshift survey is used. If we used this redshift survey, we have predicted that our uncertainties without combining with any work would be very similar to the ones given in Contarini et al. (2024).
Therefore, the most important limitation we face when constraining the cosmological parameters considered in this study is the small volume of the redshift survey sample we have used, which is approximately Mpc3. This limitation is expected to be significantly reduced when the Dark Energy Spectroscopic Instrument (DESI) experiment (DESI Collaboration et al., 2016a, b) is fully available. During its 5 years of operations, DESI will conduct the Bright Galaxy Survey (BGS) of more than 10 million galaxies over the redshift range , along with a dark-time redshift survey of 20 million luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasars (Hahn et al., 2023). With an expected footprint of 14000 and a longer redshift baseline, DESI will achieve a precision 1–2 orders of magnitude better than that of existing surveys, such as SDSS, BOSS.
In fact, we can estimate how much the errors, considering only void statistics, can decrease with DESI Y1 calculating ratio between the volumes of SDSS and DESI (see Appendix C):
(16) |
In other words, the uncertainties obtained with the Bright Galaxy Survey from DESI Y1 will be approximately a factor 2 lower than those obtained in this work with the SDSS sample used, and even lower when the full DESI survey is complete. If we additionally combine with other works, then decrease of errors would be larger when increasing the volume of the sample of galaxies used.
Acknowledgements.
E. Fernández-García acknowledges financial support from the Severo Ochoa grant CEX2021-001131-S funded by MCIN/AEI/ 10.13039/501100011033. EFG, FP and AK thanks support from the Spanish MICINN PID2021-126086NB-I00 funding grant. T.I has been supported by IAAR Research Support Program in Chiba University Japan, MEXT/JSPS KAKENHI (Grant Number JP19KK0344, JP21H01122, and JP23H04002), MEXT as “Program for Promoting Researches on the Supercomputer Fugaku” (JPMXP1020200109 and JPMXP1020230406), and JICFuS. We thank Instituto de Astrofisica de Andalucia (IAA-CSIC), Centro de Supercomputacion de Galicia (CESGA) and the Spanish academic and research network (RedIRIS) in Spain for hosting Uchuu DR1, DR2 and DR3 in the Skies Universes site for cosmological simulations. The Uchuu simulations were carried out on Aterui II supercomputer at Center for Computational Astrophysics, CfCA, of National Astronomical Observatory of Japan, and the K computer at the RIKEN Advanced Institute for Computational Science. The Uchuu Data Releases efforts have made use of the skunIAARedIRIS and skun6IAA computer facilities managed by the IAA-CSIC in Spain (MICINN EU-Feder grant EQC2018-004366-P). The other cosmological simulations were carried out on the supercomputer Fugaku provided by the RIKEN Center for Computational Science (Project ID: hp220173, hp230173, and hp230204). We also thank Joe Zuntz and Anna Porredon for their feedback on our cosmological results. Finally, we are thankful to Peter Taylor for instructing us on the the use of CombineHarvesterFlow, for verifying our results and for giving us feedback on the paper.Data Availability
The Uchuu halo and galaxy boxes and the 4 box catalogs at redshift z=0.092, as well as the 32 Uchuu-SDSS galaxy light cones, the SDSS catalogue and the void catalogues from all the previous galaxy and haloes catalogues used in this work are available at: http://www.skiesanduniverses.org/Simulations/Uchuu/, together with information on how to read the data and column description. For a list and brief description of the available halo, Uchuu-SDSS and SDSS void catalogues columns, see Appendix E.
References
- Abazajian et al. (2009) Abazajian, K. N., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS, 182, 543
- Achitouv (2019) Achitouv, I. 2019, Physical Review D, 100
- Aghanim et al. (2020) Aghanim, N., Akrami, Y., Ashdown, M., et al. 2020, A&A, 641, A6
- Aikio & Mähönen (1998) Aikio, J. & Mähönen, P. 1998, ApJ, 497, 534
- Aubert et al. (2022) Aubert, M., Cousinou, M.-C., Escoffier, S., et al. 2022, MNRAS, 513, 186
- Aung et al. (2023) Aung, H., Nagai, D., Klypin, A., et al. 2023, MNRAS, 519, 1648
- Bahcall (1977) Bahcall, N. A. 1977, ARA&A, 15, 505
- Behroozi et al. (2013a) Behroozi, P. S., Wechsler, R. H., & Wu, H.-Y. 2013a, ApJ, 762, 109
- Behroozi et al. (2013b) Behroozi, P. S., Wechsler, R. H., Wu, H.-Y., et al. 2013b, ApJ, 763, 18
- Benson et al. (2003) Benson, A. J., Hoyle, F., Torres, F., & Vogeley, M. S. 2003, MNRAS, 340, 160
- Betancort-Rijo (1990) Betancort-Rijo, J. 1990, MNRAS, 246, 608
- Betancort-Rijo (1992) Betancort-Rijo, J. 1992, Phys. Rev. A, 45, 3447
- Betancort-Rijo & López-Corredoira (2002) Betancort-Rijo, J. & López-Corredoira, M. 2002, ApJ, 566, 623
- Betancort-Rijo et al. (2009) Betancort-Rijo, J., Patiri, S. G., Prada, F., & Romano, A. E. 2009, MNRAS, 400, 1835
- Betancort-Rijo et al. (2006) Betancort-Rijo, J. E., Sanchez-Conde, M. A., Prada, F., & Patiri, S. G. 2006, ApJ, 649, 579
- Biswas et al. (2010) Biswas, R., Alizadeh, E., & Wandelt, B. D. 2010, Phys. Rev. D, 82, 023002
- Blanton et al. (2005) Blanton, M. R., Schlegel, D. J., Strauss, M. A., et al. 2005, AJ, 129, 2562
- Brout et al. (2022a) Brout, D., Scolnic, D., Popovic, B., et al. 2022a, ApJ, 938, 110
- Brout et al. (2022b) Brout, D., Taylor, G., Scolnic, D., et al. 2022b, ApJ, 938, 111
- Ceccarelli et al. (2006) Ceccarelli, L., Padilla, N. D., Valotto, C., & Lambas, D. G. 2006, MNRAS, 373, 1440–1450
- Chan et al. (2019) Chan, H. Y. J., Chiba, M., & Ishiyama, T. 2019, MNRAS, 490, 2405
- Chan et al. (2019) Chan, K. C., Hamaus, N., & Biagetti, M. 2019, Physical Review D, 99
- Chongchitnan & Silk (2010) Chongchitnan, S. & Silk, J. 2010, ApJ, 724, 285
- Clampitt et al. (2013) Clampitt, J., Cai, Y.-C., & Li, B. 2013, MNRAS, 431, 749–766
- Colberg et al. (2005) Colberg, J. M., Sheth, R. K., Diaferio, A., Gao, L., & Yoshida, N. 2005, MNRAS, 360, 216
- Colless et al. (2001) Colless, M., Dalton, G., Maddox, S., et al. 2001, MNRAS, 328, 1039
- Colless et al. (2003) Colless, M., Peterson, B. A., Jackson, C., et al. 2003, arXiv e-prints, astro
- Contarini et al. (2021) Contarini, S., Marulli, F., Moscardini, L., et al. 2021, MNRAS, 504, 5021
- Contarini et al. (2023) Contarini, S., Pisani, A., Hamaus, N., et al. 2023, ApJ, 953, 46
- Contarini et al. (2024) Contarini, S., Pisani, A., Hamaus, N., et al. 2024, A&A, 682, A20
- Contarini et al. (2022) Contarini, S., Verza, G., Pisani, A., et al. 2022, A&A, 667, A162
- Crocce et al. (2006) Crocce, M., Pueblas, S., & Scoccimarro, R. 2006, MNRAS, 373, 369
- Croton et al. (2005) Croton, D. J., Farrar, G. R., Norberg, P., et al. 2005, MNRAS, 356, 1155
- Curtis et al. (2024) Curtis, O., McDonough, B., & Brainerd, T. G. 2024, ApJ, 962, 58
- Dark Energy Survey and Kilo-Degree Survey Collaboration et al. (2023) Dark Energy Survey and Kilo-Degree Survey Collaboration, Abbott, T. M. C., Aguena, M., et al. 2023, The Open Journal of Astrophysics, 6, 36
- Dawson et al. (2013) Dawson, K. S., Schlegel, D. J., Ahn, C. P., et al. 2013, AJ, 145, 10
- DESI Collaboration et al. (2016a) DESI Collaboration, Aghamousa, A., Aguilar, J., et al. 2016a, arXiv e-prints, arXiv:1611.00036
- DESI Collaboration et al. (2016b) DESI Collaboration, Aghamousa, A., Aguilar, J., et al. 2016b, arXiv e-prints, arXiv:1611.00037
- Di Valentino et al. (2021) Di Valentino, E., Anchordoqui, L. A., Akarsu, O., et al. 2021, Astroparticle Physics, 131, 102604
- Dong-Páez et al. (2024) Dong-Páez, C. A., Smith, A., Szewciw, A. O., et al. 2024, MNRAS[arXiv:2208.00540]
- Douglass et al. (2023) Douglass, K. A., Veyrat, D., & BenZvi, S. 2023, ApJSupplement Series, 265, 7
- Einasto et al. (1994) Einasto, M., Einasto, J., Tago, E., Dalton, G. B., & Andernach, H. 1994, MNRAS, 269, 301
- El-Ad & Piran (1997) El-Ad, H. & Piran, T. 1997, ApJ, 491, 421
- Ereza et al. (2023) Ereza, J., Prada, F., Klypin, A., et al. 2023, arXiv e-prints, arXiv:2311.14456
- Falck et al. (2017) Falck, B., Koyama, K., Zhao, G.-B., & Cautun, M. 2017, MNRAS, 475, 3262–3272
- Finelli et al. (2016) Finelli, F., García-Bellido, J., Kovács, A., Paci, F., & Szapudi, I. 2016, MNRAS, 455, 1246
- Freedman et al. (2020) Freedman, W. L., Madore, B. F., Hoyt, T., et al. 2020, ApJ, 891, 57
- Gao & White (2007) Gao, L. & White, S. D. M. 2007, MNRAS, 377, L5
- Giovanelli (2010) Giovanelli, R. 2010, in Astronomical Society of the Pacific Conference Series, Vol. 421, Galaxies in Isolation: Exploring Nature Versus Nurture, ed. L. Verdes-Montenegro, A. Del Olmo, & J. Sulentic, 89
- Gkogkou et al. (2023) Gkogkou, A., Béthermin, M., Lagache, G., et al. 2023, A&A, 670, A16
- Goldwirth et al. (1995) Goldwirth, D. S., da Costa, L. N., & van de Weygaert, R. 1995, MNRAS, 275, 1185–1194
- Gorski et al. (1999) Gorski, K. M., Wandelt, B. D., Hansen, F. K., Hivon, E., & Banday, A. J. 1999, arXiv e-prints, astro
- Hahn et al. (2023) Hahn, C., Wilson, M. J., Ruiz-Macias, O., et al. 2023, The Astronomical Journal, 165, 253
- Hamaus et al. (2020) Hamaus, N., Pisani, A., Choi, J.-A., et al. 2020, J. Cosmology Astropart. Phys., 2020, 023
- Hamaus et al. (2016) Hamaus, N., Pisani, A., Sutter, P. M., et al. 2016, Phys. Rev. Lett., 117, 091302
- Hamilton (1998) Hamilton, A. J. S. 1998, in Astrophysics and Space Science Library (Springer Netherlands), 185–275
- Heymans, Catherine et al. (2021) Heymans, Catherine, Tröster, Tilman, Asgari, Marika, et al. 2021, A&A, 646, A140
- Hoffman et al. (2012) Hoffman, Y., Metuki, O., Yepes, G., et al. 2012, MNRAS, 425, 2049
- Hoffman & Shaham (1982) Hoffman, Y. & Shaham, J. 1982, ApJ, 262, L23
- Holder et al. (2001) Holder, G., Haiman, Z., & Mohr, J. J. 2001, ApJ, 560, L111
- Hoyle & Vogeley (2002) Hoyle, F. & Vogeley, M. S. 2002, ApJ, 566, 641
- Ishiyama et al. (2009) Ishiyama, T., Fukushige, T., & Makino, J. 2009, Publications of the Astronomical Society of Japan, 61, 1319
- Ishiyama et al. (2012) Ishiyama, T., Nitadori, K., & Makino, J. 2012
- Ishiyama et al. (2021) Ishiyama, T., Prada, F., Klypin, A. A., et al. 2021, MNRAS, 506, 4210
- Jennings et al. (2013) Jennings, E., Li, Y., & Hu, W. 2013, MNRAS, 434, 2167
- Kaiser (1987) Kaiser, N. 1987, MNRAS, 227, 1
- Kiang & Saslaw (1969) Kiang, T. & Saslaw, W. C. 1969, MNRAS, 143, 129
- Lavaux & Wandelt (2010) Lavaux, G. & Wandelt, B. D. 2010, MNRAS, 403, 1392
- Layzer (1954) Layzer, D. 1954, AJ, 59, 268
- Lee & Park (2009) Lee, J. & Park, D. 2009, ApJ, 696, L10
- Li et al. (2012) Li, B., Zhao, G.-B., & Koyama, K. 2012, MNRAS, 421, 3481
- Li et al. (2006) Li, C., Kauffmann, G., Jing, Y. P., et al. 2006, MNRAS, 368, 21
- Li et al. (2023) Li, S.-S., Hoekstra, H., Kuijken, K., et al. 2023, A&A, 679, A133
- Little & Weinberg (1994) Little, B. & Weinberg, D. H. 1994, MNRAS, 267, 605–628
- Martino & Sheth (2009) Martino, M. C. & Sheth, R. K. 2009, arXiv e-prints, arXiv:0911.1829
- Mauland et al. (2023) Mauland, R., Elgarøy, O., Mota, D. F., & Winther, H. A. 2023, A&A, 674, A185
- Mo & White (1996) Mo, H. J. & White, S. D. M. 1996, MNRAS, 282, 347
- Nadathur (2016) Nadathur, S. 2016, MNRAS, 461, 358
- Oogi et al. (2023) Oogi, T., Ishiyama, T., Prada, F., et al. 2023, MNRAS, 525, 3879
- Patiri et al. (2006a) Patiri, S. G., Betancort-Rijo, J., & Prada, F. 2006a, MNRAS, 368, 1132
- Patiri et al. (2012) Patiri, S. G., Betancort-Rijo, J., & Prada, F. 2012, A&A, 541, L4
- Patiri et al. (2006b) Patiri, S. G., Betancort-Rijo, J. E., Prada, F., Klypin, A., & Gottlöber, S. 2006b, MNRAS, 369, 335
- Perico et al. (2019) Perico, E. L. D., Voivodic, R., Lima, M., & Mota, D. F. 2019, A&A, 632, A52
- Planck Collaboration et al. (2016) Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016, A&A, 594, A13
- Prada et al. (2023) Prada, F., Behroozi, P., Ishiyama, T., Klypin, A., & Pérez, E. 2023, arXiv e-prints, arXiv:2304.11911
- Press & Schechter (1974) Press, W. H. & Schechter, P. 1974, ApJ, 187, 425
- Reid et al. (2015) Reid, B., Ho, S., Padmanabhan, N., et al. 2015, MNRAS, 455, 1553
- Riess et al. (2022) Riess, A. G., Yuan, W., Macri, L. M., et al. 2022, ApJ, 934, L7
- Rodríguez-Medrano et al. (2024) Rodríguez-Medrano, A. M., Springel, V., Stasyszyn, F. A., & Paz, D. J. 2024, MNRAS, 528, 2822
- Rubiño-Martín et al. (2008) Rubiño-Martín, J. A., Betancort-Rijo, J., & Patiri, S. G. 2008, MNRAS, 386, 2181
- Ryden (1995) Ryden, B. S. 1995, ApJ, 452, 25
- Sahlén et al. (2016) Sahlén, M., Zubeldía, I., & Silk, J. 2016, ApJLetters, 820, L7
- Sheth & Tormen (2002) Sheth, R. K. & Tormen, G. 2002, MNRAS, 329, 61
- Sheth & Van De Weygaert (2004) Sheth, R. K. & Van De Weygaert, R. 2004, MNRAS, 350, 517
- Song & Lee (2009) Song, H. & Lee, J. 2009, ApJ, 701, L25
- Sutter et al. (2014) Sutter, P. M., Lavaux, G., Hamaus, N., et al. 2014, MNRAS, 442, 462
- Swanson et al. (2008) Swanson, M. E. C., Tegmark, M., Hamilton, A. J. S., & Hill, J. C. 2008, MNRAS, 387, 1391
- Taylor et al. (2024) Taylor, P. L., Cuceu, A., To, C.-H., & Zaborowski, E. A. 2024, arXiv e-prints, arXiv:2406.06687
- Tikhonov (2006) Tikhonov, A. V. 2006, Astronomy Letters, 32, 727–733
- Tikhonov (2007) Tikhonov, A. V. 2007, Astronomy Letters, 33, 499–511
- Torrado & Lewis (2019) Torrado, J. & Lewis, A. 2019, Cobaya: Bayesian analysis in cosmology, Astrophysics Source Code Library, record ascl:1910.019
- Torrado & Lewis (2021) Torrado, J. & Lewis, A. 2021, J. Cosmology Astropart. Phys., 2021, 057
- van de Weygaert & Sheth (2003) van de Weygaert, R. & Sheth, R. 2003, arXiv e-prints, astro
- Voivodic et al. (2017) Voivodic, R., Lima, M., Llinares, C., & Mota, D. F. 2017, Physical Review D, 95
- von Benda-Beckmann & Mueller (2007) von Benda-Beckmann, A. M. & Mueller, V. 2007, arXiv e-prints, arXiv:0710.2783
- White (1979) White, S. D. M. 1979, MNRAS, 186, 145
- Woodfinden et al. (2023) Woodfinden, A., Percival, W. J., Nadathur, S., et al. 2023, MNRAS, 523, 6360
- Yang et al. (2005) Yang, X., Mo, H. J., van den Bosch, F. C., & Jing, Y. P. 2005, MNRAS, 357, 608
- York et al. (2000) York, D. G., Adelman, J., Anderson, John E., J., et al. 2000, AJ, 120, 1579
- Zeldovich et al. (1982) Zeldovich, I. B., Einasto, J., & Shandarin, S. F. 1982, Nature, 300, 407
- Zentner et al. (2019) Zentner, A. R., Hearin, A., van den Bosch, F. C., Lange, J. U., & Villarreal, A. S. 2019, MNRAS, 485, 1196
- Zhao et al. (2016) Zhao, C., Tao, C., Liang, Y., Kitaura, F.-S., & Chuang, C.-H. 2016, MNRAS, 459, 2670
- Zuntz et al. (2015) Zuntz, J., Paterno, M., Jennings, E., et al. 2015, Astronomy and Computing, 12, 45–59
Appendix A Void statistics theoretical framework in detail
In this Appendix we explicitly write all necessary quantities in order to calculate (VPF) and in real and redshift space. This theoretical framework can also be consulted in Betancort-Rijo et al. (2009).
Real space
The first step to calculate all the terms that appear in equation (2) is to develop equation (5). In Patiri et al. (2006a), it is shown that for dark matter haloes, can be written as
(17) |
where denotes the mean number density of those haloes in the sample (usually haloes larger than some given mass), is the volume of the sphere and the actual enclosed density contrast within the sphere. The first factor on the right-hand side of the equation is the integral of the probability density within the sphere for haloes tracing the mass (i.e. no bias, which is true in the very low mass limit). In general, haloes are biased tracers of the underlying mass distribution, due to the initial clustering of the protohaloes before they move along with mass (i.e., statistical clustering). The second term of the equation accounts for this biasing. In Patiri et al. (2006a) an approximation for this bias is obtained as a function of the linear enclosed density contrast within the sphere ():
(18) |
where , are coefficients mainly depending on the halo mass. In Patiri et al. (2006a) is demonstrated that, using Zeldovich approximation, can be rewritten as
(19) |
where is now a function of through the dependence of the actual density contrast, , on its linear counterpart, :
(20) |
is basically the relationship between the actual and the linear density contrast within a sphere as given by the standard spherical collapse model, except for a small correcting term depending on .
Doing some manipulation, can be rewritten as
(21) |
where is now
(22) |
DELF is a function of and that gives the mean actual density contrast within a sphere with radius with enclosed linear density contrast :
(23) | ||||
where denotes the relationship between the actual and linear enclosed density contrasts in the spherical collapse model (Patiri et al. 2006a):
(24) |
is the rms of the linear density contrast on a sphere with Lagrangian radius . In this equation, is evaluated at . , in equation (22) are also functions of given by (see Rubiño-Martín et al. 2008, for more details)
(25) | ||||
(26) | ||||
where , are functions of the mass of the objects and the Lagragian radius of the regions being considered:
(27) | ||||
(28) | ||||
(29) |
where is the linear growth factor normalized to be 1 at present and is the mass of the objects. is the minimum mass of distinct haloes in the sample with halo number density equal to , with this number density containing only distinct haloes. If our sample contains subhaloes, too, then the mass we have to use is .
It is important to remark that equations (25) and (26) are only valid for . If we want to calculate this functions in a different redshfit, then we have to make the following change:
(30) |
where is the linear growth factor of density fluctuations in the model under consideration.
We already have all terms in order to calculate from equation (21), so all we need now is to calculate , which is the probability function of the linear density contrast within an Eulerian space, and can be calculated as (Betancort-Rijo & López-Corredoira 2002):
(31) | ||||
where
(32) |
and
(33) |
(34) |
(35) |
(36) |
This fit is valid for Mpc and .
Redshift space
The theoretical framework developed above is only valid for voids in real space. However, if we want to constrain using surveys such as SDSS, this surveys provide galaxy positions in redshift space. In this space (redshift space), the peculiar velocity of the galaxies is added to the velocity expansion of the Universe. This generates some distortions which result in elongated structures known as Fingers of God (see Hamilton (1998) for more details).
If we want to calculate and in redshift space instead of real space, all we have to do is change by:
(37) |
where is the radius of a sphere in real space and is the radius of that sphere in redshift space. Additionally, it has been found comparing the average outflow around the relevant voids with that given by the spherical expansion model that the value of is around 0.85. This is also the value that provide the best agreement with simulations.
VEL() is defined so that the peculiar velocity , of mass element at distance from the centre of a spherical mass concentration (or defect) enclosing actual density contrast is given by
(38) |
where is the Hubble constant at the time being considered. In Betancort-Rijo et al. (2006) it is shown that
(39) |
where is the growth factor as a function of the expansion factor, , and DELK() is the inverse function of DELT() (see Mo & White 1996; Sheth & Tormen 2002):
(40) | ||||
where is the linear density contrast for spherical collapse model, which for the concordance cosmology at present is 1.676.
Finally, the expression of is:
(41) |
In Table 8 the values of VPF predicted by theoretical framework and obtained through the four halo simulation boxes are shown. It can be checked that theoretical framework predicts an increase of when going from real to redshift space, as it happens in simulations.
Appendix B Comparison of voids statistics in the distribution of haloes and galaxies in mocks with theoretical framework
In this Appendix, we show the results obtained with the theoretical framework presented in Appendix A and compare it with the results of the halo and galaxy simulation boxes. We do this in real and redshift space. We also check if we recover the values of and using the theoretical framework with Maximum Likelihood test with Bayesian approach.
B.1 Real Space
Statistics of voids obtained in halo simulation boxes and predicted by theoretical framework in real space have already been discussed in Section 4. Here, we briefly summarize the most important results obtained and present the constraints achieved of and . This is a crucial step to take to check that the theoretical framework works correctly, as it is essential to recover the real values of and of each simulation.
In Table 7 the values of obtained with theoretical framework and simulations can be seen. The results are represented in Figure 5. As we have already seen in Section 4, the dependence of the VPF predicted by theoretical framework is the same as that shown by simulations, and all values with Mpc are within 10 of the ratio between simulations and theoretical framework, being compatible with unity those with Mpc.
In Figure 6, we can see the number density of voids larger than for the four halo simulation boxes. As we have already seen, the agreement between theoretical framework and simulation values is good, specially for intermediate radius bins (i.e. between 12 and 18 Mpc).
Therefore, we are now in the position to use the theoretical framework and simulations to constrain and H0 (we fix in each case to the real value of each simulation box) and obtain the confidence contours using maximum likelihood test and the standard Bayesian approach. This can be seen in the first column of Figure 14. In this Figure, we can see that all real values (black points) are inside the 1 contour.
Therefore, we successfully recover the values of and of the simulations with the theoretical framework, and we can take one more step to do this study in redshift space.
r | ||||||||
---|---|---|---|---|---|---|---|---|
9 | 7.216 | 7.782 | 7.163 | 7.703 | 6.936 | 7.520 | 6.682 | 7.250 |
10 | 3.847 | 4.176 | 3.810 | 4.097 | 3.659 | 3.966 | 3.485 | 3.762 |
11 | 1.946 | 2.106 | 1.922 | 2.044 | 1.829 | 1.955 | 1.719 | 1.821 |
12 | 0.9340 | 0.9980 | 0.9194 | 0.9547 | 0.8659 | 0.9023 | 0.8015 | 0.8179 |
13 | 0.4251 | 0.4452 | 0.4169 | 0.4206 | 0.3881 | 0.3889 | 0.3530 | 0.3455 |
14 | 0.1833 | 0.1884 | 0.1790 | 0.1722 | 0.1645 | 0.1556 | 0.1467 | 0.1399 |
15 | 0.07484 | 0.07423 | 0.07277 | 0.06674 | 0.06586 | 0.05914 | 0.05743 | 0.04977 |
16 | 0.02891 | 0.02775 | 0.02797 | 0.02448 | 0.02489 | 0.02032 | 0.02116 | 0.01721 |
17 | 0.01055 | 9.863 | 0.01015 | 8.260 | 8.872 | 6.640 | 7.331 | 5.540 |
18 | 3.638 | 3.421 | 3.479 | 2.940 | 2.979 | 2.410 | 2.385 | 1.800 |
19 | 1.183 | 1.017 | 1.124 | 1.180 | 9.412 | 9.500 | 7.276 | 5.800 |
20 | 2.628 | 2.958 | 3.422 | 4.800 | 2.796 | 3.800 | 2.080 | 1.800 |
21 | 1.048 | 1.125 | 9.804 | 1.700 | 7.801 | 1.600 | 5.567 | 1.100 |
B.2 Redshift Space
Now we can move to redshift space and check if our model to transform distances from real to redshift space (i.e. equation (37)) is accurate enough.
In the upper panel of Figure 15 we can observe the VPF (multiplied by to distinguish more easily each simulation for large radii) obtained by theoretical framework (continuous lines) and simulations (dots), and in Table 8 we can see the numerical values. In the bottom panel of the same Figure, the ratio between simulations and theoretical framework is shown. Again, taking into account the errors, the agreement between theoretical framework and simulations is good.
In the upper panel of Figure 16 we can see the . In the bottom panel of the same Figure, the ratio between simulations and theoretical framework is shown. The agreement between simulations and theoretical framework is good, too, although there are big oscillations for large radius ) which make theoretical framework not to be compatible with simulations. This can be due to not having enough statistics for these radius bins.
r | ||||||||
---|---|---|---|---|---|---|---|---|
9 | 9.059 | 9.440 | 8.999 | 9.333 | 8.643 | 9.151 | 8.181 | 8.842 |
10 | 5.099 | 5.379 | 5.054 | 5.291 | 4.810 | 5.139 | 4.490 | 4.888 |
11 | 2.737 | 2.901 | 2.707 | 2.832 | 2.550 | 2.718 | 2.344 | 2.535 |
12 | 1.402 | 1.482 | 1.383 | 1.434 | 1.288 | 1.355 | 1.163 | 1.233 |
13 | 0.6848 | 0.7187 | 0.6732 | 0.6817 | 0.6192 | 0.6376 | 0.5481 | 0.5635 |
14 | 0.3189 | 0.3310 | 0.3124 | 0.3060 | 0.2833 | 0.2820 | 0.2452 | 0.2442 |
15 | 0.1414 | 0.1437 | 0.1380 | 0.1205 | 0.1232 | 0.1184 | 0.1040 | 0.09690 |
16 | 0.05973 | 0.05862 | 0.05803 | 0.05212 | 0.05091 | 0.04628 | 0.04177 | 0.03675 |
17 | 0.02400 | 0.0264 | 0.02320 | 0.2025 | 0.01997 | 0.01684 | 0.01588 | 0.01320 |
18 | 9.164 | 8.300 | 8.815 | 7.33 | 7.426 | 5.970 | 5.706 | 4.400 |
19 | 3.324 | 2.983 | 3.179 | 2.330 | 2.617 | 2.290 | 1.936 | 1.570 |
20 | 1.145 | 9.667 | 1.088 | 8.800 | 8.735 | 9.500 | 6.201 | 5.700 |
21 | 3.738 | 2.625 | 3.529 | 3.700 | 2.758 | 3.200 | 1.872 | 2.100 |
The contours obtained with maximum likelihood test for the halo simulation boxes in redshift space can be seen in the second column of Figure 14. In that Figure, we can see that we obtain a small contour for Uchuu in redshift space, unlike in real space. Additionally, we recover Planck 2018 values of and H0 within for the P18 box catalog, and within for Low and VeryLow. Therefore, in redshift space we don’t successfully recover the real values of the parameters in simulations within as we do in real space. This can be due to the transformation we have used in order to transform the theoretical framework from real to redshift space (see equation (37)). However, we can recover the real values withing for the four halo simulation boxes.
Appendix C Scaling of errors with volume
An important study that must be made is the scaling of errors with the volume of the samples used, i.e., how much do errors decrease when the volume of the sample is increased. From Figure 17 it is evident that the confidence level contours (and, therefore, the constrained values of each parameter) decrease considerably when we consider a sample with a larger volume. In addition, we can see that the confidence level contour size of a single random Uchuu-SDSS lightcone voids is similar to the one obtained for SDSS voids.
In Table 9 we can see the constrained values for Uchuu-SDSS void and SDSS voids, and the ratio of the SDSS voids and Uchuu-SDSS voids. This ratio has been calculated as the mean value between the ratios of the upper and lower errors (for example, for the first row, . It can be seen that we obtain similar values of the ratio for and H0. We can also see that the ratios for and (the fundamental parameters of the theoretical framework) are very far from , but they are closer to . Therefore, we can predict that the scaling of the errors with the volume of the samples is very similar to , with being the ratio of the volumes of the two samples.
Uchuu-SDSS voids | SDSS voids | Ratio | |
---|---|---|---|
0.793 | 1.044 | 1.58 | |
0.305 | 0.298 | 1.09 | |
H0 | 75.29 | 84.43 | 1.21 |
0.1787 | 0.1981 | 2.43 | |
S8 | 0.792 | 1.036 | 1.67 |
Appendix D Theoretical framework with different number densities of galaxies
In this work, we have constrained the parameters , and H0 using galaxy samples with a number density of Mpc3 at a redshift of . We have verified that the theoretical framework successfully predicts both the Void Probability Function and the number density of voids larger than for this number density. However, it remains to be checked whether the theoretical framework continues to successfully predict these two statistics at a different number densities of galaxies.
In order to check if the theoretical framework works correctly for different number densities of galaxies we need to construct different samples (from Uchuu galaxy simulation box, for example) with different galaxy number densities, calculate the VPF and from that sample and compare with the predicted value given by the theoretical framework with the same coefficients that we have calculated in this work. The only parameter that must be changed for each sample is the mass, . The number densities we use are Mpc3.
In Figure 18 we can observe the VPF (left panel) and number density of voids larger than (right panel) predicted by theoretical framework with continuous lines and obtained in Uchuu galaxy simulation boxes for these galaxy number densities with dots and the abundance of vois larger than (right panel). It can be seen that the theoretical framework predicts successfully the VPF for all galaxy number densities without any need to change the values of or , except for galaxy number densities smaller than Mpc3. Therefore, the theoretical framework is only valid for large galaxy number densities (Mpc3), and the coefficients and depend on this galaxy number density.
Appendix E Content of the void catalogues used in this work
The columns of the void catalogues of halo and galaxy simulation boxes are the following
-
•
X[MPC/H]: x-position of the centre of the void (comoving Mpc).
-
•
Y[MPC/H]: y-position of the centre of the void (comoving Mpc).
-
•
Z[MPC/H]: z-position of the centre of the void (comoving Mpc).
-
•
RADIUS[MPC/H]: radius of the void (comoving Mpc)
For Uchuu-SDSS lightcones and SDSS, there are four additional columns:
-
•
RA[DEG]: right ascension (degrees)
-
•
DEC[DEG]: declination (degrees).
-
•
ZOBS: observed redshift of the centre of the void (accounting for peculiar velocities. The fiducial cosmology used for SDSS voids is Planck 2015 (Planck Collaboration et al. 2016)).
-
•
completeness: mean completeness of the void, calculated as the mean completeness of all the points uniformly distributed in its volume.