Next Article in Journal
Combining Empirical and Physics-Based Models for Solar Wind Prediction
Previous Article in Journal
CP Conservation in the Strong Interactions
Previous Article in Special Issue
Discord in Concordance Cosmology and Anomalously Massive Early Galaxies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Betti Functionals as Probes for Cosmic Topology

Institute of Theoretical Physics, Ulm University, Albert-Einstein-Allee 11, D–89069 Ulm, Germany
*
Author to whom correspondence should be addressed.
Universe 2024, 10(5), 190; https://doi.org/10.3390/universe10050190
Submission received: 21 March 2024 / Revised: 9 April 2024 / Accepted: 10 April 2024 / Published: 24 April 2024
(This article belongs to the Special Issue Cosmic Microwave Background)

Abstract

:
The question of the global topology of the Universe (cosmic topology) is still open. In the Λ CDM concordance model, it is assumed that the space of the Universe possesses the trivial topology of R 3 , and thus that the Universe has an infinite volume. As an alternative, in this paper, we study one of the simplest non-trivial topologies given by a cubic 3-torus describing a universe with a finite volume. To probe cosmic topology, we analyze certain structure properties in the cosmic microwave background (CMB) using Betti functionals and the Euler characteristic evaluated on excursions sets, which possess a simple geometrical interpretation. Since the CMB temperature fluctuations δ T are observed on the sphere S 2 surrounding the observer, there are only three Betti functionals β k ( ν ) , k = 0 , 1 , 2 . Here, ν = δ T / σ 0 denotes the temperature threshold normalized by the standard deviation σ 0 of δ T . The analytic approximations of the Gaussian expectations for the Betti functionals and an exact formula for the Euler characteristic are given. It is shown that the amplitudes of β 0 ( ν ) and β 1 ( ν ) decrease with an increasing volume V = L 3 of the cubic 3-torus universe. Since the computation of the β k ’s from observational sky maps is hindered due to the presence of masks, we suggest a method that yields lower and upper bounds for them and apply it to four Planck 2018 sky maps. It is found that the β k ’s of the Planck maps lie between those of the torus universes with side-lengths L = 2.0 and L = 3.0 in units of the Hubble length and above the infinite Λ CDM case. These results give a further hint that the Universe has a non-trivial topology.

1. Introduction

A wide range of cosmological data is well-described within the framework of the Λ cold dark matter ( Λ CDM) model, which is now established as the standard model of cosmology. Due to the increasing level of precision in the available data, several discrepancies have arisen as the Hubble tension and the S 8 tension in the recent years [1,2,3,4,5]. In addition, there are several strange features in the cosmic microwave background (CMB) as summarized in [6], which also point to a modification of the standard model. A prominent example is the suppression of the quadrupole moment in the CMB angular power spectrum, which is also especially revealed by the 2-point angular correlation function C ( ϑ ) showing almost no correlations above angles of 60 . Models without relying on modified physics, which can address this feature, are provided by cosmic topology [7,8,9,10,11], see also references in [12]. A possible non-trivial topology for the Universe can suppress the large-scale anisotropy due to an infrared cut-off in the wave number spectrum.
A non-trivial topology can be detected by searching for topologically matched circles in CMB maps, the so-called circles-in-the-sky (CITS) test [13]; however, the searches were in vain up to now [14,15]. A possibility for this negative result might be that the CITS signature is not so clearly pronounced in the CMB sky maps, as predicted by the Λ CDM model on which the likelihoods are based. A future resolution of the Hubble and S 8 tensions might, for example, lead to a larger integrated Sachs–Wolfe contribution, which would additionally blur the CITS signature so that the Universe might possess a non-trivial topology despite the negative result. Another possibility for a non-trivial topology is that the topology produces CITS, which are far from antipodal as in the case of the Hantzsche–Wendt topology [16], which impede the CITS search. An exhaustive analysis of the allowed non-trivial topological models is given in [17,18].
In this paper, we analyze certain structure properties in the CMB maps of a compact non-trivial topological space and compare them with the standard Λ CDM model, which presuppose an infinite volume. This requires the simulation of CMB sky maps for these topological models, as outlined in [19,20,21]. There, it is discussed that the CMB simulations require the determination of the eigenvalue spectrum and the eigenmodes of the Laplace–Beltrami operator on the topological space. The obvious restriction would be to test only the topological spaces that are not excluded by the CITS test, but their CMB sky maps are more difficult to simulate than the simple cubic 3-torus model whose CITS signature could not be found. The reason is that the highest degeneracy in the eigenvalue spectrum belongs to the cubic 3-torus topology such that the transfer function has to be computed for a significantly smaller number of eigenvalues, which, in turn, speeds up the CMB simulations correspondingly. For that reason, we consider here the cubic 3-torus space [12], which should provide an idea of general properties.
In this work, we apply the Betti functionals to the CMB temperature fluctuation field δ T ( n ^ ) , which is defined on the sphere S 2 that surrounds the observer. ( n ^ denotes the direction in which the temperature fluctuation δ T ( n ^ ) is observed). For the analysis of the CMB, the temperature fluctuation field δ T ( n ^ ) is normalized
f ( n ^ ) : = δ T ( n ^ ) μ σ 0 ,
where μ and σ 0 are the mean and the standard deviation of the field δ T ( n ^ ) . Therefore, the normalized field f ( n ^ ) has zero mean and unit variance. The topological descriptors are then computed as a functional of the excursion set
E ( ν ) = { n ^ S 2 f ( n ^ ) ν } ,
where ν denotes the threshold.
For our special case S 2 , there are only three Betti functionals (BF) β k ( ν ) . These are the number of components of the excursion set by the 0-th BF β 0 ,
β 0 ( ν ) : = # components of E ( ν ) ,
the number of independent 1-dimensional punctures on S 2 , denoted by the 1-st BF β 1 ,
β 1 ( ν ) : = # topological holes of E ( ν ) ,
and finally, the 2-nd BF β 2 , which counts the number of the internal voids of E ( ν ) . In our special case, this is one if the excursion set is identical to S 2 being the case for ν < min ( f ( n ^ ) ) and zero otherwise, i.e.,
β 2 ( ν ) : = { 1 for ν < min ( f ( n ^ ) ) 0 otherwise .
The Betti numbers originate from the analysis of topological spaces, where they are used to distinguish the topological spaces based on the connectivity of n-dimensional simplicial complexes. The nth Betti number is defined by the rank of the nth homology group [22,23,24,25]. An analysis based on relative homology with respect to cosmic topology will be given in the companion paper [26].
Concerning the word “topological” in Equation (4), there is a remark in order. At first, consider a pixelized version of the sphere S 2 and identify it with the excursion set for the case ν < min ( f ( n ^ ) ) . Then increase the threshold ν such that a single pixel is removed from the sphere, which destroys the internal void, which, in turn, is counted by β 2 , but this does not create a 1-dimensional hole. The 1-dimensional hole is created on the sphere after removal of another pixel not sharing a boundary with the first removed pixel. In this case, we now have a ring that has an actual hole, which is called a topological hole. Thus, the number of topological holes is the number of common sense holes minus one.
Let us note that the application of the Betti functionals is not restricted to a study of cosmic topology, as illustrated in this paper. In fact, the Betti functionals can also be used to analyze the structure properties of the CMB predicted in models with a different scenario for dark energy or in the context of more generally modified gravity theories.
In this paper, we use for the Λ CDM concordance model the cosmological parameters, as given by the Planck collaboration in [27] in their Table 4 in the column ‘TT+lowP+lensing’. The main parameters are Ω b h 2 = 0.02226 , Ω c h 2 = 0.1186 , and h = 0.678 . Furthermore, all sky maps are computed with the Healpix resolution N side = 128 and, if not explicitly stated otherwise, a Gaussian smearing of 120 arcmin is used. For the cubic 3-torus model, 1000 CMB sky maps are computed for each of the side-lengths L = 0.5 , 1.0 , 1.5 , 2.0 , and 3.0 , where L is given in units of the Hubble length L H = c / H 0 . This allows the computation of ensemble averages for the β ’s in order to compare them with the infinite Λ CDM standard model and the Planck CMB maps.

2. Analytic Approximation of Betti Functionals and an Exact Gaussian Expectation for the Euler Characteristic and the Genus

After many years of work on the Betti functionals (BF) β k ( ν ) ( k = 0 , 1 , 2 ) , there still exist no exact predictions for the ensemble expectations assuming that the CMB anisotropy is a homogeneous, isotropic Gaussian random field on the S 2 -sphere. Therefore, there arises the question whether one can find reliable analytic approximations for the β k ( ν ) as a function of ν . In the following, we propose such approximations, which agree well with the average values obtained from simulations for a cubic 3-torus with side-length L = 2.0 .
A crucial rôle is played by the asymptotic behaviour of the BFs in the limits ν ± . From the definition of the BFs, it is obvious that it holds
β k ( ν ) 0 for ν , k = 0 , 1 , 2 .
The situation is, however, different in the limit ν where the excursion set covers the whole unit sphere S 2 . In this case, one obtains from the Betti numbers  p k of S 2 the asymptotic behaviour for ν
β 0 ( ν ) p 0 = 1 β 1 ( ν ) p 1 = 0 β 2 ( ν ) p 2 = 1 .
This is in agreement with the general relation for the Euler characteristic (EC) χ ( ν ) ,
χ ( ν ) = k = 0 2 ( 1 ) k β k ( ν ) ,
yielding for the full sphere the correct value χ ( ) = 2 .
From (6) and (7), one infers that the apparent symmetry (see Figure 1), β 0 ( ν ) = β 1 ( ν ) , with respect to the parity transformation ν ν ( δ T δ T , respectively), does not hold true but rather is violated. The degree of the breaking of parity can be described by a function α 0 ( ν ) ,
β 0 ( ν ) = β 1 ( ν ) + α 0 ( ν )
which according to the Equations (6) and (7) has to satisfy the asymptotic relations
α 0 ( ν ) 0 for ν α 0 ( ν ) 1 for ν .
The relation (8) for the EC is then
χ ( ν ) = β 1 ( ν ) β 1 ( ν ) + α 0 ( ν ) + β 2 ( ν ) .
Under very general assumptions, one can expand the expression in the brackets of (11) into a convergent infinite series in terms of the odd Hermite functions φ 2 n + 1 ( ν ) exp ( ν 2 / 2 ) He 2 n + 1 ( ν ) , where He 2 n + 1 ( ν ) denote the odd “probabilist’s” Hermite polynomials (see Appendix C in [28] for a mathematical exposition of the general Hermite expansions). As a first approximation, let us consider the first term ( n = 0 ) in this expansion, which leads with He 1 ( ν ) = ν to the approximation
β 1 ( ν ) β 1 ( ν ) = c 1 ν e ν 2 / 2
with a positive coefficient given below. In Figure 2 we show that (12) gives an excellent fit to the data. Equation (12) describes a maximum at ν = 1 and a minimum at ν = 1 with amplitude ± c 1 / e in nice agreement with the peaks of the data very close to ± 1 (see Figure 2). Inserting relation (12) into (11), one obtains the following approximation to the EC:
χ ( ν ) = c 1 ν e ν 2 / 2 + α 0 ( ν ) + β 2 ( ν ) .
In order to determine c 1 and α 0 ( ν ) we use the standard definition of the EC using the Gauss-Bonnet theorem on the excursion set
E ( ν ) K d a + E ( ν ) κ ( s ) d s = 2 π χ ( ν ) ,
where K = 1 denotes the Gaussian curvature of S 2 , d a the surface element on S 2 , d s the line element along E ( ν ) , and κ ( s ) the geodesic curvature of E ( ν ) . The integrals in (14) are proportional to the Minkowski functionals (MFs) υ 0 ( ν ) and υ 2 ( ν ) [28], respectively, which gives
χ ( ν ) = 2 υ 0 ( ν ) + 4 π υ 2 ( ν ) .
The MFs have the nice property that their exact Gaussian predictions are explicitly known (see [28] and references therein),
υ 0 ( ν ) = 1 2 erfc ν 2 υ 2 ( ν ) = ρ 2 2 ( 2 π ) 3 / 2 ν exp ( ν 2 / 2 ) .
Here, the parameter ρ : = σ 1 / σ 0 has been studied in [29], where it has been shown that ρ does hierarchically detect the change in size of the cubic 3-torus if the volume of the Universe is smaller than ≃2.5 ×   10 3 Gpc 3 . ( σ 1 is the standard deviation of the gradient of the CMB field δ T ( n ^ ) , i.e., ρ is the normalized standard deviation of the CMB gradient field).
Figure 1. The Betti functionals β 0 ( ν ) and β 1 ( ν ) are shown together with their 1 σ error bands obtained from 1000 CMB simulations based on the torus topology. In addition, the mean value of β 0 ( ν ) + β 1 ( ν ) is plotted as a dotted curve. The resolution parameters are set to N side = 128 , l max = 256 and FWHM = 120 arcmin . These parameters are used in all figures if not otherwise noted.
Figure 1. The Betti functionals β 0 ( ν ) and β 1 ( ν ) are shown together with their 1 σ error bands obtained from 1000 CMB simulations based on the torus topology. In addition, the mean value of β 0 ( ν ) + β 1 ( ν ) is plotted as a dotted curve. The resolution parameters are set to N side = 128 , l max = 256 and FWHM = 120 arcmin . These parameters are used in all figures if not otherwise noted.
Universe 10 00190 g001
Figure 2. The fit (12) is shown as a dashed curve together with the mean value of β 1 ( ν ) β 1 ( ν ) obtained from 1000 CMB simulations based on the torus topology. The quality of the fit is so good that the mean value curve is plotted with a larger line width in order to reveal the dashed curve belonging to the one-parameter fit. The inset shows the difference between the data and the fit.
Figure 2. The fit (12) is shown as a dashed curve together with the mean value of β 1 ( ν ) β 1 ( ν ) obtained from 1000 CMB simulations based on the torus topology. The quality of the fit is so good that the mean value curve is plotted with a larger line width in order to reveal the dashed curve belonging to the one-parameter fit. The inset shows the difference between the data and the fit.
Universe 10 00190 g002
Inserting (16) into (15), the exact Gaussian expectation for the EC reads
χ ( ν ) = ρ 2 2 π ν e ν 2 / 2 + erfc ν 2 ,
which compared with the approximation (13) gives
c 1 = ρ 2 2 π
and
α 0 ( ν ) = erfc ν 2 β 2 ( ν ) .
In Figure 3, we show the mean value of β 2 ( ν ) obtained from 1000 CMB simulations based on the torus topology and compare it with the approximation
β 2 ( ν ) = 1 2 erfc b ( ν + a ) .
Note that the fit parameters a , b > 0 have a simple interpretation since
β 2 ( a ) = 1 2
and
β 2 ( a ) = b π .
Inserting (20) into (19) leads to the analytic approximation
α 0 ( ν ) = erfc ν 2 1 2 erfc b ( ν + a ) ,
also shown in Figure 3, and thus to an explicit expression for the parity violation, see Equation (9).
Figure 3. The mean value of β 2 ( ν ) obtained from 1000 CMB simulations based on the torus topology is shown together with the fit (20). Also shown is α 0 ( ν ) obtained from (21).
Figure 3. The mean value of β 2 ( ν ) obtained from 1000 CMB simulations based on the torus topology is shown together with the fit (20). Also shown is α 0 ( ν ) obtained from (21).
Universe 10 00190 g003
To the best of our knowledge, Figure 3 presents, for the first time, a computation of the BF β 2 ( ν ) and of the function α 0 ( ν ) quantifying the breaking of parity symmetry, defined in Equation (9). In previous papers, β 2 was considered to be zero and α 0 ( ν ) was not discussed at all. From the definition (5), it follows that β 2 ( ν ) is small since it is bounded by
0 β 2 ( ν ) 1 with β 2 ( ν ) = 0 for ν > ν 2 ( ν 2 > 0 ) .
Also the parity violation α 0 ( ν ) is small since it follows from (19) and (22)
α 0 ( ν ) erfc ν 2 2
resp.
α 0 ( ν ) = erfc ν 2 for ν > ν 2 .
In contrast, the Figure 1, Figure 2, Figure 4 and Figure 5 clearly show that the β 0 ( ν ) , β 1 ( ν ) , and the EC χ ( ν ) have very large maxima (minima), which finds a nice explanation by the large value of the parameter ρ . Indeed, we obtain from (12) and (18) for the difference
δ 1 ( ν ) : = β 1 ( ν ) β 1 ( ν )
at the maximum at ν = 1 (minimum at ν = 1 )
δ 1 ( ± 1 ) = ± ρ 2 2 π e = ± 0.2420 ρ 2 .
In [29], it was shown that the mean value of < ρ ( L ) > can be well-approximated for tori of the size of 1.0 L 3.0 by the linearly decreasing function
< ρ ( L ) > 46.122 3.290 L .
From this, we obtain, for example, for the torus with L = 2.0 , the large value < ρ ( 2 ) > 39.542 , which leads to δ 1 ( ± 1 ) ± 378 , which explains the large amplitudes displayed in Figure 2. (A similar prediction follows for the BF β 1 ( ν ) , see Equations (28) and (32) and Figure 1). Thus, the BFs β 0 ( ν ) and β 1 ( ν ) can be used to detect the size of the Universe if it is modeled as a cubic 3-torus. See also Figures 8, 11 and 12.
Figure 4. The mean value of β 1 ( ν ) is shown together with a shifted Gaussian. The inset displays the difference between the Gaussian fit and β 1 ( ν ) in addition to the 1 σ band confirming that the Gaussian fit describes the data better than 1 σ .
Figure 4. The mean value of β 1 ( ν ) is shown together with a shifted Gaussian. The inset displays the difference between the Gaussian fit and β 1 ( ν ) in addition to the 1 σ band confirming that the Gaussian fit describes the data better than 1 σ .
Universe 10 00190 g004
The Figure 4 demonstrates that a simple shifted Gaussian
β 1 ( ν ) = N 2 π Σ exp ( ν ν 1 ) 2 2 Σ 2 = β 1 ( 0 ) exp ν ( ν 2 ν 1 ) 2 Σ 2
fits the mean value of β 1 ( ν ) better than the 1 σ width of the distribution. The normalization N, the width Σ of the Gaussian, and the position ν 1 of its maximum lead to the relations
β 1 ( 0 ) = N 2 π Σ exp ν 1 2 2 Σ 2 = β 1 ( ν 1 ) exp ν 1 2 / ( 2 Σ 2 ) = β 1 ( 2 ν 1 ) .
Furthermore, one derives from (28)
β 1 ( 0 ) = ν 1 Σ 2 β 1 ( 0 ) ,
and from (12) and (18)
β 1 ( 0 ) = ρ 2 2 2 π .
Comparing the relations (30) and (31) yields the relation
β 1 ( 0 ) = Σ 2 2 2 π | ν 1 | ρ 2 ,
which explains (see also (29) and Figure 4) the large amplitude of β 1 ( ν ) due to the large value of ρ 2 .
Finally, Figure 5 shows the EC for 1000 simulations using the exact Gaussian expectation value (17) of the general relation (8). Equation (17) predicts two large extrema (for the ρ -values obtained from (27)) at
ν ± : = ± 1 2 ρ 2 ± 1
with magnitude χ ( ν + ) 379 and χ ( ν ) 377 , respectively; for the torus with L = 2 , one has ν ± = ± 0.999 . It is interesting to note that even in the case where the primordial initial conditions are exactly Gaussian, there is a small “parity breaking” (not visible in Figure 5) of the antisymmetry (negative parity) of χ ( ν ) since it holds
χ ( ν ) = χ ( ν ) + 2 .
Figure 5. The EC χ ( ν ) , Equation (8) is shown.
Figure 5. The EC χ ( ν ) , Equation (8) is shown.
Universe 10 00190 g005
Another important topological quantity is the genus g ( ν ) related to the EC by
g ( ν ) = 1 1 2 χ ( ν ) .
Inserting (17) into (35), one obtains for the Gaussian prediction of the genus the exact relation
g ( ν ) = ρ 2 2 2 π ν e ν 2 / 2 + 1 2 erf ν 2 + 1 2
and
g ( ν ) = g ( ν ) + 1 .
From (36), one obtains the special values g ( ) = 1 (torus), g ( 0 ) = 1 2 , and g ( ) = 0 (sphere).

3. Betti Functionals for the Cubic 3-Torus Topology

In this section, the properties of β 0 ( ν ) and β 1 ( ν ) are discussed for the torus simulations with different torus side-lengths L. Since the focus is put on the simulations, one does not have to bother about masked sky regions, which will be discussed in the next section.
The Figure 6 shows the influence of the Gaussian smearing of the sky maps on β 0 ( ν ) and β 1 ( ν ) . With increasing smearing, the excursions set has fewer structure elements, and thus, the amplitude of β 0 ( ν ) and β 1 ( ν ) decreases with increasing smoothing. This behaviour is nicely revealed in Figure 6 for the case L = 2.0 . Since the 3-torus simulations of the sky maps are computed up to l max = 256 , which roughly corresponds to a resolution of 180 / l max 0 . 7 , a Gaussian smoothing of at least 2 is on the safe side. As already stated, the analysis in this paper is based on a Gaussian smoothing of 2 .
We now compare, for this fixed Gaussian smoothing, the dependence of the number of components and holes in dependence on the side-length L of the 3-torus in Figure 7. One observes a nice decreasing behaviour of β 0 ( ν ) and β 1 ( ν ) , with an increasing 3-torus size. The Figure 7 also displays the corresponding result for the infinite Λ CDM model, which is computed using CAMB for the same set of cosmological parameters as used for the 3-torus models. This extrapolates the 3-torus size towards infinity. So this monotone dependence on the size of the 3-torus might give a hint to the size of our Universe by studying β 0 ( ν ) and β 1 ( ν ) in observed sky maps.
Figure 6. The ensemble averages of β 0 ( ν ) and β 1 ( ν ) are shown for the torus side-length L = 2.0 with dependence on the Gaussian smoothing of the sky maps. The smoothing of 2 = 120 arcmin , 3 = 180 arcmin , and 4 = 240 arcmin are displayed.
Figure 6. The ensemble averages of β 0 ( ν ) and β 1 ( ν ) are shown for the torus side-length L = 2.0 with dependence on the Gaussian smoothing of the sky maps. The smoothing of 2 = 120 arcmin , 3 = 180 arcmin , and 4 = 240 arcmin are displayed.
Universe 10 00190 g006
Figure 7. The ensemble averages of β 0 ( ν ) and β 1 ( ν ) are shown depending on the torus side-length for L = 0.5 up to L = 3.0 . In addition, the ensemble averages of β 0 ( ν ) and β 1 ( ν ) for the infinite Λ CDM model are shown, which extrapolates the torus side-length towards infinity. A decreasing number of structure elements with increasing side-length L are revealed. A Gaussian smoothing of 120 arcmin is used in all cases.
Figure 7. The ensemble averages of β 0 ( ν ) and β 1 ( ν ) are shown depending on the torus side-length for L = 0.5 up to L = 3.0 . In addition, the ensemble averages of β 0 ( ν ) and β 1 ( ν ) for the infinite Λ CDM model are shown, which extrapolates the torus side-length towards infinity. A decreasing number of structure elements with increasing side-length L are revealed. A Gaussian smoothing of 120 arcmin is used in all cases.
Universe 10 00190 g007
For a consistent comparison of the ensemble averages of the finite torus model with the infinite Λ CDM model provided by the Planck collaboration, it is important that the calculations are performed for both models with the same cosmological parameters. Only then it is guaranteed that the L-dependence demonstrated in Figure 7 is an unambiguous result, which is not caused by using different cosmological parameters.
An interesting connection of this monotonic dependence on the size of the topological cell exists with respect to the normalized standard deviation of the CMB gradient field, i.e., of ρ already introduced in Section 2. Solving the relation (26) with respect to ρ yields
ρ = 2 π e β 1 ( 1 ) β 1 ( 1 ) .
The interesting point is that this relation shows that a monotonic dependence of β 1 on the torus size leads to an analogous behaviour of ρ . While β 1 ( ± 1 ) is computed by counting the number of holes of the excursion set, ρ is computed by differentiating the CMB temperature field. In [29], the behaviour of the mean value of < ρ ( L ) > is analyzed and a monotonic behaviour of ρ is found, which can be approximated for tori of the size of 1.0 L 3.0 by the already given linearly decreasing function (27). In Figure 8, the linear behaviour (27) determined in [29] is compared with that derived from β 1 ( ± 1 ) by using (38). A nice agreement between both methods for the computation of ρ is observed. The Figure 8 also shows the result by using β 1 min ( ± 1 ) / f sky and β 1 max ( ± 1 ) / f sky , which will be defined in the next section.
Figure 8. The dashed line shows the fit (27) for ρ determined in [29], while the three full curves display the results obtained from β 1 ( ν ) using (38), where β 1 ( ν ) of the full sky simulation and β 1 min ( ν ) / f sky and β 1 max ( ν ) / f sky are used. The latter ones are determined on the masked sky, as explained in Section 4.
Figure 8. The dashed line shows the fit (27) for ρ determined in [29], while the three full curves display the results obtained from β 1 ( ν ) using (38), where β 1 ( ν ) of the full sky simulation and β 1 min ( ν ) / f sky and β 1 max ( ν ) / f sky are used. The latter ones are determined on the masked sky, as explained in Section 4.
Universe 10 00190 g008

4. Betti Functionals in the Presence of Masks

The analysis of the Betti functionals in CMB observations is impeded by foregrounds that do not allow a measurement of the genuine CMB. Furthermore, in the case of ground-based observations one has to deal with an incomplete sky coverage. In order to estimate the Betti functionals in the case of a mask, we propose the following procedure. At first, for the computation of β 0 ( ν ) , one counts the number of components in the unmasked sky. Then, there arises the possibility that some or all of those components that are partially covered by a common connected masked region, might be linked within this region, as illustrated in Figure 9. There, three components are shown that might be connected within the masked domain and would be counted as one component in an ideally measured sky without a mask. Of course, a further possibility is that only two of them are connected. Therefore, a lower bound β 0 min ( ν ) is obtained by counting all components that are touched by a common masked region as a single component. Conversely, an upper bound β 0 max ( ν ) is obtained by treating all components as separated. Note that no attempt is made to estimate the number of components within the masked domains, since only components outside the masked regions are considered. In order to estimate the number of components in the full sky, β 0 min ( ν ) and β 0 max ( ν ) are divided by
f sky : = Area of the unmasked sky Area of full sky .
The same procedure applies analogously for β 1 .
Figure 9. The ambiguity of counting the number of structure elements in the presence of a mask is illustrated. The mask is pictured as the dashed region, while the structure elements are shown as blue regions. If only the parts outside the masked region are known, one cannot decide whether they are connected within the mask or not. The two extreme cases are that they are all connected (right) or none of them (left).
Figure 9. The ambiguity of counting the number of structure elements in the presence of a mask is illustrated. The mask is pictured as the dashed region, while the structure elements are shown as blue regions. If only the parts outside the masked region are known, one cannot decide whether they are connected within the mask or not. The two extreme cases are that they are all connected (right) or none of them (left).
Universe 10 00190 g009
In the case of simulations, one can test this procedure. Figure 10 compares for the case L = 2.0 the upper and lower bounds with the true full sky result. The resolution parameters are again set to N side = 128 and FWHM = 120 arcmin .
In this work, we use the Planck 2018 “Component Separation Common mask in Intensity” [30], which can be obtained at http://pla.esac.esa.int/pla/#maps, accessed on 20 March 2024 (file name: COM_Mask_ CMB-common-Mask-Int_2048_R3.00.fits). Since our analysis is based on the Healpix resolution N side = 128 , we downgrade the above mask from N side = 2048 to N side = 128 . The downgraded mask no longer has only the pixel values of 0 and 1, but also intermediate values, and we use a mask threshold of 0.9 in the following analysis. This leads to f sky = 0.7615 .
One observes from Figure 10 that the counted number of components lies nicely between β 0 min ( ν ) / f sky and β 0 max ( ν ) / f sky , as it should be. The same is seen for β 1 , which refers to the holes.
Figure 10. The extrapolations β 0 , 1 min ( ν ) / f sky and β 0 , 1 max ( ν ) / f sky obtained from the masked CMB maps are shown in comparison with β 0 , 1 ( ν ) derived from the unmasked CMB maps for the L = 2.0 ensemble.
Figure 10. The extrapolations β 0 , 1 min ( ν ) / f sky and β 0 , 1 max ( ν ) / f sky obtained from the masked CMB maps are shown in comparison with β 0 , 1 ( ν ) derived from the unmasked CMB maps for the L = 2.0 ensemble.
Universe 10 00190 g010
Furthermore, for sufficiently low values of ν , β 0 max ( ν ) / f sky saturates at a non vanishing positive value. This is due to the structure of the mask. In the case without a mask, all CMB values are larger than ν for that sufficiently low ν , so that the full sphere S 2 is obtained as the excursion set E ( ν ) . Applying the mask and assuming that the components are not connected within the masked region counts them as separate items if they lie within a “hole” of the mask. Then, each hole of the mask yields a separate component. In contrast, β 0 min ( ν ) / f sky assumes that all components are linked within the masked domains and thus one counts only a single component.

5. A Comparison of the Cubic 3-Torus Topology with the Planck CMB Maps

The Planck collaboration provides four CMB sky maps for a cosmological analysis. Here, we use the Planck 2018 maps [31] called SMICA, Commander, NILC, and SEVEM, which can be obtained at http://pla.esac.esa.int/pla/#maps, accessed on 20 March 2024. The Healpix routine “map2alm_ iterative” is used to compute the spherical expansion coefficients for the four Planck maps by taking the mask Component Separation Common mask in Intensity into account. Thereafter, the monopole and dipole are set to zero, and a Gaussian smoothing of 120 arcmin is applied. These spherical expansion coefficients are used to generate four Healpix maps with the resolution of N side = 128 , which can then be compared with the corresponding 3-torus CMB maps. The Figure 11 and Figure 12 display the curves derived from these four Planck curves as solid and in the common red colour, since they are nearly indistinguishable.
In Figure 11, β 0 min ( ν ) / f sky and β 1 min ( ν ) / f sky are plotted for the 3-torus simulations, the infinite Λ CDM model and the four Planck maps are all subjected to the same mask. The curves present the lower estimate of the true β ’s. It is seen that the curves derived from the Planck maps possess a significantly larger amplitude than that of the infinite Λ CDM model. Indeed, they lie between the cubic 3-torus simulations of the side-lengths L = 2.0 and L = 3.0 , whereas L = 2.0 provides the better match. A similar behaviour is seen in Figure 12, where β 0 max ( ν ) / f sky and β 1 max ( ν ) / f sky are plotted. In this case, the Planck derived curves are closer to the L = 3.0 case. Thus, in both cases, there seems to be an indication of a finite size of our Universe corresponding to a size between L = 2.0 and L = 3.0 .
Figure 11. β 0 min ( ν ) / f sky and β 1 min ( ν ) / f sky are displayed for the 3-torus simulations, the infinite Λ CDM model and the four Planck maps, which are all subjected to the same mask.
Figure 11. β 0 min ( ν ) / f sky and β 1 min ( ν ) / f sky are displayed for the 3-torus simulations, the infinite Λ CDM model and the four Planck maps, which are all subjected to the same mask.
Universe 10 00190 g011
Figure 12. The same as in Figure 11, but for β 0 max ( ν ) / f sky and β 1 max ( ν ) / f sky .
Figure 12. The same as in Figure 11, but for β 0 max ( ν ) / f sky and β 1 max ( ν ) / f sky .
Universe 10 00190 g012

6. Discussion and Summary

In the quest for the global topology of the Universe, there have been several methods suggested to unveil the spatial structure on its very largest scales. In this paper, we focus on the Betti functionals applied to the excursion sets E ( ν ) , Equation (2), of the CMB, which possess, in this case, a simple geometrical interpretation. With the CMB observed on a sphere S 2 , the excursion set decomposes this sphere with respect to the normalized temperature threshold ν into components and holes. There are three Betti functionals in this case. The first one, β 0 , counts the number of connected components, the second one, β 1 , the number of topological holes, and finally, β 2 , the number of two-dimensional cavities. The latter takes the value of one if the excursion set is the complete sphere; that is, if the threshold ν is smaller than the lowest normalized temperature on the CMB sphere, otherwise β 2 is zero. Since the Betti functionals focus only on the number of structure elements, they are even simpler than the Minkowski functionals. This is because the Minkowski functionals require the computation of the area, the circumference, and a curvature measure of the boundary of the components [28]. As discussed in Section 2, the Minkowski functionals allow a derivation of a relation connecting the average of the normalized standard deviation of the CMB gradient field called ρ with β 1 , see Equations (26) and (38), if the CMB is assumed to be a homogeneous, isotropic Gaussian random field. Thus, the properties of ρ and the Betti numbers are not independent.
The common lore is that a spatially finite universe is betrayed by the large scale behaviour of the CMB, for example, the suppression of the quadrupole moment C 2 or the low power in the 2-point angular correlation function C ( ϑ ) above sufficiently large angles on the sky, typically above 60 . Often overlooked, a suppression at significantly smaller angles of the angular correlation function C ( ϑ ) is additionally seen such that the amplitude of C ( ϑ ) for the 3-torus models is below that of the infinite Λ CDM model, see [12]. It should be emphasized that this small angle suppression is also visible in C ( ϑ ) obtained from the observed sky.
The definition of ρ as a differential measure reveals its obvious local nature, so that a topological signature on small scales exists also for this quantity, since a dependence of ρ on the volume of the cubic 3-torus was demonstrated in [29], see also Figure 8. In Section 3, it is shown that β 0 ( ν ) and β 1 ( ν ) display a hierarchical dependence of their amplitudes with respect to the side-length L of the cubic 3-torus, see Figure 7, such that the amplitudes increase with decreasing volume V = L 3 . This behaviour is in nice agreement with the normalized standard deviation ρ of the CMB gradient field. It reveals the local structure in the excursion set via β 1 ( ν ) at ν = ± 1 , see (38). However, since β 0 and β 1 are, of course, not restricted to the thresholds ν = ± 1 , they provide a more comprehensive tool than ρ .
The computation of β k ( ν ) from observational sky maps is hindered due to the presence of masks. The number of connected components and holes is then ambiguous since it is not discernible whether they are connected within the not measured parts, i. e. within the mask. In Section 4, a method is suggested, which gives, for their number, a lower and an upper bound within the observed sky. Finally, Section 5 applies this method to four sky maps released by the Planck collaboration in 2018, called SMICA, Commander, NILC, and SEVEM. The comparison with the cubic 3-torus simulations shows that the curves derived from the four Planck maps lie between the 3-torus models with side-length L = 2.0 and L = 3.0 , see Figure 11 and Figure 12. So this measure gives a further hint that our Universe has a non-trivial topology.

Author Contributions

Both authors declare that they have contributed substantially to the joint reported work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data will not be made public.

Acknowledgments

We would like to thank Thomas Buchert, Martin France, and Pratyush Pranav for discussions. The software packages HEALPix (http://healpix.jpl.nasa.gov, accessed on 20 March 2024, [32]) and CAMB written by A. Lewis and A. Challinor (http://camb.info, accessed on 20 March 2024) as well as the Planck data from http://pla.esac.esa.int/pla/#maps, accessed on 20 March 2024, were used in this work.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Anchordoqui, L.A.; Di Valentino, E.; Pan, S.; Yang, W. Dissecting the H0 and S8 tensions with Planck + BAO + supernova type Ia in multi-parameter cosmologies. J. High Energy Astrophys. 2021, 32, 28–64. [Google Scholar] [CrossRef]
  2. 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: 8 and S8. Astropart. Phys. 2021, 131, 102604. [Google Scholar] [CrossRef]
  3. Abdalla, E.; Abellán, G.F.; Aboubrahim, A.; Agnello, A.; Akarsu, Ö.; Akrami, Y.; Alestas, G.; Aloni, D.; Amendola, L.; Anchordoqui, L.A.; et al. Cosmology intertwined: A review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies. J. High Energy Astrophys. 2022, 34, 49–211. [Google Scholar] [CrossRef]
  4. Vagnozzi, S. Seven Hints That Early-Time New Physics Alone Is Not Sufficient to Solve the Hubble Tension. Universe 2023, 9, 393. [Google Scholar] [CrossRef]
  5. Akarsu, Ö.; Colgáin, E.Ó.; Sen, A.A.; Sheikh-Jabbari, M.M. Λ CDM Tensions: Localising Missing Physics through Consistency Checks. arXiv 2024, arXiv:2402.04767. [Google Scholar]
  6. Schwarz, D.J.; Copi, C.J.; Huterer, D.; Starkman, G.D. CMB anomalies after Planck. Class. Quantum Grav. 2016, 33, 184001. [Google Scholar] [CrossRef]
  7. Ellis, G.F.R. Topology and cosmology. Gen. Rel. Grav. 1971, 2, 7–21. [Google Scholar] [CrossRef]
  8. Lachièze-Rey, M.; Luminet, J.P. Cosmic topology. Phys. Rep. 1995, 254, 135–214. [Google Scholar] [CrossRef]
  9. Levin, J. Topology and the cosmic microwave background. Phys. Rep. 2002, 365, 251–333. [Google Scholar] [CrossRef]
  10. Luminet, J.P. The Status of Cosmic Topology after Planck Data. Universe 2016, 2, 1. [Google Scholar] [CrossRef]
  11. Ryan, P.J. Cosmic topology, underdetermination, and spatial infinity. Eur. J. Philos. Sci. 2024, 14, 17. [Google Scholar] [CrossRef]
  12. Aurich, R.; Janzer, H.S.; Lustig, S.; Steiner, F. Do we live in a ‘small universe’? Class. Quantum Grav. 2008, 25, 125006. [Google Scholar] [CrossRef]
  13. Cornish, N.J.; Spergel, D.N.; Starkman, G.D. Circles in the sky: Finding topology with the microwave background radiation. Class. Quantum Grav. 1998, 15, 2657–2670. [Google Scholar] [CrossRef]
  14. Aurich, R.; Lustig, S. A search for cosmic topology in the final WMAP data. Mon. Not. R. Astron. Soc. 2013, 433, 2517–2528. [Google Scholar] [CrossRef]
  15. Ade, P.A.R. et al. [Planck Collaboration] Planck 2013 results. XXVI. Background geometry and topology of the Universe. Astron. Astrophy. 2014, 571, A26. [Google Scholar] [CrossRef]
  16. Aurich, R.; Lustig, S. The Hantzsche–Wendt manifold in cosmic topology. Class. Quantum Grav. 2014, 31, 165009. [Google Scholar] [CrossRef]
  17. Akrami, Y. et al. [COMPACT Collaboration] Promise of Future Searches for Cosmic Topology. arXiv 2022, arXiv:2210.11426. [Google Scholar]
  18. Petersen, P. et al. [COMPACT Collaboration] Cosmic topology. Part I. Limits on orientable Euclidean manifolds from circle searches. J. Cosmol. Astropart. Phys. 2023, 2023, 030. [Google Scholar] [CrossRef]
  19. Riazuelo, A.; Weeks, J.; Uzan, J.P.; Lehoucq, R.; Luminet, J.P. Cosmic microwave background anisotropies in multiconnected flat spaces. Phys. Rev. D 2004, 69, 103518. [Google Scholar] [CrossRef]
  20. Aurich, R.; Lustig, S. Cosmic microwave anisotropies in an inhomogeneous compact flat universe. Class. Quantum Grav. 2011, 28, 085017. [Google Scholar] [CrossRef]
  21. Eskilt, J.R. et al. [COMPACT Collaboration] Cosmic topology. Part IIa. Eigenmodes, correlation matrices, and detectability of orientable Euclidean manifolds. J. Cosmol. Astropart. Phys. 2024, 03, 036. [Google Scholar] [CrossRef]
  22. Adler, R.J. The Geometry of Random Fields; Classics in Applied Mathematic; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 2010. [Google Scholar]
  23. Edelsbrunner, H.; Harer, J.L. Computational Topology: An Introduction; American Mathematical Society: Providence, RI, USA, 2022. [Google Scholar]
  24. Munkres, J.R. Elements of Algebraic Topology; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  25. Pranav, P. Anomalies in the topology of the temperature fluctuations in the cosmic microwave background: An analysis of the NPIPE and FFP10 data releases. Astron. Astrophy. 2022, 659, A115. [Google Scholar] [CrossRef]
  26. Pranav, P.; Aurich, R.; Buchert, T.; France, M.J.; Steiner, F. Homology of the Cosmic Microwave Background: Are finite universe models preferred over the infinite standard model? arXiv 2024. in preparation. [Google Scholar]
  27. Ade, P.A.R. et al. [Planck Collaboration] Planck 2015 results—XIII. Cosmological parameters. Astron. Astrophy. 2016, 594, A13. [Google Scholar] [CrossRef]
  28. Buchert, T.; France, M.J.; Steiner, F. Model-independent analyses of non-Gaussianity in Planck CMB maps using Minkowski functionals. Class. Quantum Grav. 2017, 34, 094002. [Google Scholar] [CrossRef]
  29. Aurich, R.; Buchert, T.; France, M.J.; Steiner, F. The variance of the CMB temperature gradient: A new signature of a multiply connected Universe. Class. Quantum Grav. 2021, 38, 225005. [Google Scholar] [CrossRef]
  30. Akrami, Y. et al. [Planck Collaboration] Planck 2018 results—IV. Diffuse component separation. Astron. Astrophy. 2020, 641, A4. [Google Scholar] [CrossRef]
  31. Aghanim, N. et al. [Planck Collaboration] Planck 2018 results—I. Overview and the cosmological legacy of Planck. Astron. Astrophy. 2020, 641, A1. [Google Scholar] [CrossRef]
  32. Górski, K.M.; Hivon, E.; Banday, A.J.; Wandelt, B.D.; Hansen, F.K.; Reinecke, M.; Bartelmann, M. HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere. Astrophys. J. 2005, 622, 759–771. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Aurich, R.; Steiner, F. Betti Functionals as Probes for Cosmic Topology. Universe 2024, 10, 190. https://doi.org/10.3390/universe10050190

AMA Style

Aurich R, Steiner F. Betti Functionals as Probes for Cosmic Topology. Universe. 2024; 10(5):190. https://doi.org/10.3390/universe10050190

Chicago/Turabian Style

Aurich, Ralf, and Frank Steiner. 2024. "Betti Functionals as Probes for Cosmic Topology" Universe 10, no. 5: 190. https://doi.org/10.3390/universe10050190

APA Style

Aurich, R., & Steiner, F. (2024). Betti Functionals as Probes for Cosmic Topology. Universe, 10(5), 190. https://doi.org/10.3390/universe10050190

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