Issue |
A&A
Volume 683, March 2024
|
|
---|---|---|
Article Number | A176 | |
Number of page(s) | 10 | |
Section | Cosmology (including clusters of galaxies) | |
DOI | https://doi.org/10.1051/0004-6361/202348160 | |
Published online | 18 March 2024 |
Probing the statistical isotropy of the universe with Planck data of the cosmic microwave background
1
Universidade Federal do Espírito Santo, Centro Universitário Norte do Espírito Santo, Departamento de Ciências Naturais, Rodovia BR 101 Norte, km. 60, São Mateus, ES 29932-540, Brazil
e-mail: wiliam.ricaldi@ufes.br
2
Observatório Nacional, Rua General José Cristino 77, São Cristóvão, Rio de Janeiro, RJ 20921-400, Brazil
3
Núcleo de Astrofísica e Cosmologia (Cosmo-Ufes), Universidade Federal do Espírito Santo, Vitória, ES 29075-910, Brazil
Received:
4
October
2023
Accepted:
14
January
2024
Aims. We study the angular distribution of temperature fluctuations in the cosmic microwave background (CMB) to probe the statistical isotropy of the universe by using precise full-sky CMB data with a model-independent approach.
Methods. We investigated the temperature-temperature angular correlations in the four Planck foreground-cleaned CMB maps that were released recently. We performed a directional analysis on the CMB sphere to search directions in which the temperature-temperature angular correlations are extreme.
Results. Our analyses confirm a preferred axis in the CMB sphere, pointing in the direction (l, b)≃(260° ,130° ), at the 98%−99% confidence level. In this direction, the CMB angular correlations exceed the antipodal direction most strongly. This preferred direction is unexpected in the ΛCDM cosmological model and represents a significant deviation from results obtained by applying the same procedure to simulated statistically isotropic CMB maps. This result confirms the north-south asymmetry in the most recent Planck data. This phenomenon is one of the previously reported CMB anomalies.
Conclusions. We performed a robust detection of the north-south asymmetry in the temperature-temperature angular correlations, with a slightly different statistical significance, in the four Planck foreground-cleaned CMB maps. Moreover, we performed consistency tests by adding foreground and noise, both Planck data products, to the CMB map we studied, and we also investigated and discarded possible bias in our method. After these detailed analyses, we conclude that the north-south asymmetry phenomenon is present with a high statistical significance in the Planck CMB maps we studied. This result confirms previous reports in the literature in the past 20 years.
Key words: cosmic background radiation / cosmology: observations
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
The angular distribution of temperature fluctuations in the cosmic microwave background (CMB) is a valuable cosmological probe for understanding primordial and late processes in the universe. Moreover, these data can also probe the statistical isotropy (SI) hypothesis, which is a fundamental attribute of the concordance cosmological model (Jarosik et al. 2011; Planck Collaboration VI 2020).
During the past 20 years, diverse analyses of the CMB data at large angular scales, both from the Wilkinson Microwave Anisotropy Probe (WMAP) and Planck satellites, have provided evidence of significant deviations at various confidence levels of the SI hypothesis. These features are known in the literature as CMB anomalies (Planck Collaboration VII 2020). These intriguing features are important for understanding the global properties of the universe, such as its spatial topological shape. It is therefore valuable to determine whether the reported statistical anisotropy is cosmological in origin, and, if so, extract information that may be helpful for identifying its causes. To avoid possible bias related to a specific analysed dataset or a particular estimator, it is worthwhile that different studies from diverse research groups outside the Planck Collaboration team can perform complementary analyses (Planck Collaboration VII 2020).
Since the early release of WMAP data, many analyses have reported unexpected features in the angular distribution of the CMB, mainly at large scales. Among these are the low-quadrupole amplitude of the CMB angular power spectrum (Bennett et al. 2003; Jarosik et al. 2011), the alignment of some low-order multipoles (de Oliveira-Costa et al. 2004; Bielewicz et al. 2004, 2005; Schwarz et al. 2004; Land & Magueijo 2005; Abramo et al. 2009; Polastri et al. 2015), the lack of large-angle temperature-temperature (TT) correlations (Schwarz et al. 2004, 2016; Copi et al. 2006, 2009; Bernui et al. 2006), the parity asymmetry (Kim & Naselsky 2010; Gruppuso et al. 2011; Aluri & Jain 2012; Hansen et al. 2012; Zhao 2014; Aluri et al. 2017), the cold spot (Vielva et al. 2004; Cruz et al. 2005; Bernui 2009; Gurzadyan et al. 2014), and the north-south asymmetry (Eriksen et al. 2004a,b; Hansen et al. 2004a,b; Hoftuft et al. 2009; Paci et al. 2010; Pietrobon et al. 2010; Rath et al. 2013; Bernui et al. 2014; Gimeno-Amo et al. 2023). Recently, it was also reported that the CMB lensing amplitude, AL, has a value higher than expected in the ΛCDM model (Calabrese et al. 2008; Planck Collaboration VI 2020; Mokeddem et al. 2023).
Some methods that investigated SI features in the CMB data used cosmology-dependent approaches. In this way, the obtained results are of restricted validity because the outcomes depend on the hypotheses that were assumed. In this scenario, the conclusions from model-independent analyses are preferable because they are universally valid. For this reason, we adopted a model-independent approach based on the two-point angular correlation function (2PACF) in spherical caps. This approach allowed us to perform a directional analysis of the TT correlations over the CMB sphere, searching for directions with outlier properties. Through this directional analysis, we measured possible deviations from the SI hypothesis by quantifying the statistical significance of our findings by comparison with the outcomes of a similar directional procedure applied to a large number of SI CMB maps (Bernui et al. 2007a). We applied this method to the latest Planck CMB temperature-fluctuations maps, that is, the four foreground-cleaned Planck CMB maps, and we also considered in the analysis the effects caused by residual foregrounds and instrumental noise because they might affect our results. In particular, our directional method is optimal for detecting preferred directions (i.e., asymmetries) on the CMB sky.
Our directional approach also allowed us to study the CMB at different large scales (different γ0 radius, as we discuss below). This technique can provide information about the cause(s) of the anomalous behavior. Applications of this method include, for instance, the study of the Galactic masks used for analyses of WMAP data, concluding that they weakly affect the direction and statistical significance of the north-south asymmetry (Bernui et al. 2007b). In addition, we also investigated that the cold-spot anomaly does not cause the north-south asymmetry (Bernui 2009; for other interesting results using this directional analysis approach, see, e.g., Marques et al. 2018 for the analysis of the Planck convergence map, or Bernui et al. 2008 for the study of the large-angle distribution of the GRB). Additionally, because several reported CMB anomalies are related to angular correlations or to preferred directions, an approach that combines the two features is highly welcome. Our method is based on 2PACF and shows a preferred axis on the sky, indicating a possible connection with the CMB quadrupole-octopole alignment (see the discussion in Bernui & Hipólito-Ricaldi 2008 for a universe with preferred axis, sourced by a primordial magnetic field, in the analysis of WMAP data, and in our Appendix A for the general case considering the analysis of the latest Planck data).
On the other hand, even though the north-south asymmetry phenomenon has been reported since 2004 (Eriksen et al. 2004a; Hansen et al. 2004a), the study of precise CMB data sets from WMAP and Planck using diverse procedures and algorithms, an exhaustive reanalysis of the latest and most precise data (including the use of Planck CMB products) allows us to capture possible systematics or data-processing effects on the north-south phenomenon. It also allows a reevaluation of its statistical significance. Our main contributions in studying this phenomenon are: (a) In a model-independent way, the latest Planck CMB maps show an SI violation and north-south asymmetry. The results for both the direction and the statistical significance are robust. We found this in all the foreground-cleaned maps and through consistency tests; (b) we show that SI violation and north-south asymmetry in the latest Planck CMB maps are not correlated with instrumental noise, Galactic foregrounds, or known systematics; (c) with our directional analysis (which is described in detail in Sect. 3) we found an anisotropic behavior of the angular correlations, revealing that a large part of the sky has small large-angle correlations. This is probably related to another CMB anomaly: the lack of power of the CMB large-angle correlations (Copi et al. 2006; d) we found that the asymmetry in the angular correlations appears at the largest scales. We discuss (see Appendix A) the possible relation of this to the quadrupole-octopole multipole alignment.
Differently from studies that assumed a fiducial model to study the CMB north-south asymmetry with the angular power spectrum {Cℓ} (see, e.g., Eriksen et al. 2004a; Hansen et al. 2004a; Rath et al. 2013; Yeung & Chu 2022), or studies that assumed a cosmological model with a dipolar modulation a priori (see, e.g., Rath & Jain 2013; Mukherjee et al. 2016; Planck Collaboration XVIII 2016; Gimeno-Amo et al. 2023), our approach is a model-independent method that is based on a directional analysis using the 2PACF as estimator to probe cosmic isotropy. Consequently, all our findings here are also model independent.
The outline of this work is the following. In Sect. 2 we briefly present the Planck CMB data that we investigated with our directional analysis procedures. In Sect. 3 we describe the method of this SI scrutiny, and we also detail how the statistical significance of possible outlier directions can be quantified. In Sect. 4 we present our results, discuss consistency tests to determine whether they are indeed robust, and in Sect. 5 we formulate our conclusions and final remarks.
2. Planck CMB data
In 2020, the Planck Collaboration released the fourth set of CMB products (PR4) derived from the Planck mission1.
These data products include the four Planck foreground-cleaned CMB temperature maps (Planck Collaboration IV 2020), which are called Spectral Matching Independent Component Analysis (SMICA), Needlet Internal Linear Combination (NILC), Spectral Estimation Via Expectation Maximization (SEVEM), and Bayesian parametric method Commander (Commander). These high-resolution maps have Nside = 2048 and an effective beam: full width at half maximum (FWHM) FWHM = 5′. These CMB maps result from diverse data analyses that were used to remove known foregrounds, in which each component-separation algorithm processed the multi-contaminated sky regions differently, such as the Galactic plane region. As a result, they exhibit different foreground-cleaned regions for each procedure. These regions are defined outside the so-called component-separation confidence mask, which is released jointly with the associated CMB map, corresponding to fsky ≃ 0.83, 0.77, 0.81, and 0.87, for the SMICA, NILC, SEVEM, and Commander maps, respectively2.
Our main objective in this work is to analyze the angular correlations of these four foreground-cleaned Planck CMB maps in detail with our directional estimator (Bernui et al. 2007a). An important point is the consistency tests, which we performed by applying our directional analysis to investigate these four CMB maps. For example, these maps might reveal different asymmetry directions, which might suggest a systematic effect and may not be the result of a cosmological origin. After this, we repeat our analyses considering the addition of foregrounds and white noise, both Planck CMB data products, to the CMB foreground-cleaned Planck map in order to test the possibility that our results are sourced by residual contaminants.
3. Method and implementation
Next, we describe the method we adopted to study the CMB angular correlations present in the foreground-cleaned Planck maps. These CMB maps have been investigated by the Planck Collaboration with several consistency tests. They are therefore reliable for CMB analyses outside the region defined by the corresponding confidence mask.
3.1. Method
Our interest here is to search for possible preferred directions of the large-scale TT angular correlations in the Planck CMB maps. To quantify these statistical features, we define the 2PACF as
where γ is the angle between pixels in directions and
, and it is given by
The average ⟨ ⟩ is computed over all pairs of pixels with an angular separation γ. The computation of this quantity over the entire celestial sphere gives no indication of possible outlier directions, although it shows an anomalous lack of strength at large-angle TT correlations (see, e.g., Copi et al. 2004). Instead, we searched for preferred directions by computing the 2PACF of TT correlations in a set of spherical caps by scanning the full CMB sky.
We assumed that the CMB temperature-fluctuation sphere is pixelized. We defined a spherical cap with a radius equal to γ0, centered on the J pixel with angular coordinates (θJ, ϕJ) as
where 𝒮2 represents the celestial sphere. In each of these caps, the 2PACF was computed by using
with γ ∈ ⟨0, 2γ0]3. We denote by the 2PACF value for the angular distances inside the subinterval γk ∈ γ, such that γk ∈ ⟨(k − 1)β, kβ] for k = 1, …, Nbins, where Nbins = 2γ0/β is the number of bins and β is the bin length.
The function C = C(γ) measures the 2PACF in different directions because it encompasses the information along different spherical caps. The next step was to quantify the difference between angular correlations observed in different directions. In order to do this, we used the σ map method (Bernui et al. 2007a; Bernui 2008). The σ map is a model-independent method for detecting directions in the CMB maps with different angular correlation intensities. The analysis made through the 2PACF on spherical caps that scan the celestial sphere and then compare the observational data results with the same procedure applied to synthetic SI maps. With this method, σ maps are celestial maps based on real values, which are then converted into colored maps, containing a measure of the angular correlations from each region of the celestial sphere. To quantify the correlation intensity, it is possible to define a scalar function to associate a non-negative real value to each spherical cap with center in (θJ, ϕJ), that is,
for J = 1, …, Ncaps, where σJ is defined by
with σJ = σJ(θJ, ϕJ). Then, the set of real positive numbers for all J = 1, …, Ncaps is the σ map. To detect any possible asymmetry in any direction, we decomposed the σ map in spherical harmonics,
and analyzed its power spectrum,
The mean, the standard deviation, and the distribution of values {Sℓ} of the σ map power spectra obtained from a directional analysis applied to a set of SI CMB maps, help us to quantify how frequent or rare the power spectra of the Planckσ maps are (i.e., the σ maps obtained by analyzing the four foreground-cleaned CMB maps). A high value for the dipole S1, for instance, would indicate a preferred axis on the sky, that is, a possible violation of the SI of the universe.
3.2. Angular correlation analyses
To implement the method presented in the previous subsection and to obtain a quantitative measurement of possible directional differences in the angular correlation of a CMB map, we covered the sky with Ncaps spherical regions and computed the set of {σJ, J = 1, …, Ncaps} values, that is, its σ map. We therefore followed the next steps: (i) we downgraded the CMB maps and its respective masks, and (ii) we generated the set of spherical caps with a radius γ0 with CMB data from a given Planck masked map; (iii) we then computed the TT 2PACF in all these spherical caps; and (iv) we produced the Planckσ map and obtained its angular power spectrum. As discussed above, we analyzed the four Planck foreground-cleaned CMB temperature maps (Planck Collaboration VI 2020). These maps have a high resolution with Nside = 2048 and 50.331.648 pixels, which corresponds to a multipole ℓ ∼ 2500. Because several anomalies were reported at large angular scales, that is, ℓmax ≤ 30 (see, e.g., de Oliveira-Costa et al. 2004; Eriksen et al. 2004a; Copi et al. 2006, 2009; Bernui et al. 2007a), we are interested here in studying a possible violation of the SI hypothesis at these angular scales. To do this, it is necessary to downgrade each of the original Planck foreground-cleaned maps.
We downgrade the four Planck maps from Nside = 2048 to Nside = 32, which corresponds to ℓmax = 95 and 12 288 pixels. The same downgrade process was required for their corresponding masks. The original masks have pixels with 0 and 1 boolean values. During the downgrade process, several pixels are combined to form one single pixel. In this process, the numerical value in each new pixel is the mean value of all mixed pixels. In general, these values are not 0 nor 1. It is therefore necessary to define a threshold as a criterion for reconstructing the 0 or 1 values. We defined this criterion by considering that pixels with values ≥0.8 were considered as 1, and 0 otherwise. At the end of the downgrade step, we obtained a CMB map and its corresponding mask in the correct resolution to continue our analyses. The next step was to define the spherical caps for which the 2PACF was computed.
We defined the number and size of the spherical caps with CMB temperature data to analyze their 2PACF. We chose γ0 = 60°, but to confirm the consistency of our results, we finally also considered other values for γ0 (see Sect. 4). Our directional analysis performs a uniform scan of the CMB sky meaning that the centers of the spherical caps must be uniformly distributed. A practical way to define these cap centers uniformly is to use a Healpix pixelization with a suitable angular resolution (i.e., ). We chose
, which means that we analyzed the 2PACF in Ncaps = 192 spherical caps. Clearly, other values for
can be considered to test the robustness of our results. To summarize, given a CMB map, our method consists of selecting for scrutiny 192 spherical caps of 60° each and computing the 2PACF to the data inside each cap by means of Eq. (4) to obtain a set of 192 2PACF: {C(γ)J} for J = 1, …, 1924.
Then we applied Eq. (6) to this set of functions to obtain a set of 192 real positive numbers , which is full of directional information regarding the strength of angular correlations of the CMB map we analyzed. In practice, this set of numbers,
, was converted into a colored map with 192 pixels, that is, the σ map. This map has an angular resolution
. The angular power spectrum of the σ map was then calculated, as usual, with Eqs. (7) and (8). The power spectrum of the σ map was then compared to the mean spectra of the corresponding σ maps we obtained, according to the same procedure, from a set of 1000 SI CMB maps. This comparison provided the statistical significance of our directional analysis. Our results are presented in the next section.
4. Results of the directional analysis
In this section, we present the results of the directional analyses we applied to the recently released Planck CMB maps. According to the method described in the previous section, we first obtained the corresponding σ map of the analysed Planck CMB map, calculated its angular power spectrum, and then evaluated its statistical significance by comparison with spectra obtained by performing the same procedure applied to simulated CMB maps, called Monte Carlo (MC) CMB maps. These MC maps were generated by assuming the SI hypothesis using the concordance ΛCDM model angular power spectrum as seed (Planck Collaboration VI 2020). Before we applied the directional analysis to these MC CMB maps, we verified that they shared the same observational features as the Planck maps, such as angular resolution and Galactic mask. This analysis allowed us to quantify the significance level of our findings and to evaluate whether the SI hypothesis was satisfied in the investigated Planck CMB data.
4.1. Commander-σ map
In this section, we present the directional analysis of the Commander CMB map, the σ map we obtained, its angular power spectrum, and other related results.
In the top panel of Fig. 1, we show the result of our directional analysis applied to the CommanderPlanck map, that is, the Commander-σ map. In the bottom panel, we display the dipole component of this map. The Commander-σ map exhibits a dipolar red-blue region in the distribution of the angular correlations, which is strongly suggestive of a large-angle anisotropy. This can be better quantified in the Commander-σ map power spectrum , which is represented by dots in Fig. 2. By comparing
with its corresponding MC-σ map power spectra, it is possible to evaluate whether they are consistent with each other.
![]() |
Fig. 1. Top: Commander-σ map, which is the result of a directional analysis of the TT angular correlations in 192 spherical caps, with a radius γ0 = 60°, applied to the CommanderPlanck CMB map. Two extremes values in almost antipodal positions are visible. Bottom: Dipole component of the Commander-σ map (top panel), which shows the dipole direction of the north-south asymmetry phenomenon: (l, b)≃(253° ,137° ). |
After our directional analysis determined the dipolar feature in the analysed data, our objective was to quantify how common or rare this feature is in σ maps produced from SI CMB maps. To do this, we applied the same method to a set of SI MC CMB maps and obtained a distribution of σ map spectra , for ℓ = 1, …, 5. For our analyses with the Commander map, we generated 1000 MC maps and used the corresponding Commander fiducial mask to remove the Galactic region5. For each MC CMB map, we then computed the corresponding σ map (hereafter called MC-σ map) and its angular power spectrum Sℓ. After the directional analysis of these 1000 MC CMB maps, we calculated their corresponding angular power spectra. With this distribution in hand, we then calculated the mean and other statistical features. With this distribution, we calculated the statistical significance of
. The results are presented in Fig. 2. Moreover, the Commander CMB map was masked with its fiducial mask to remove the Galactic region that might have residual foregrounds. We computed the Commander-σ map and its power spectrum
, and finally, we compared it with their corresponding MC-σ map results, as shown in Fig. 2.
![]() |
Fig. 2. Angular power spectrum of the Commander-σ map (dots), shown in Fig. 1, where the continuous blue line corresponds to the mean angular power spectrum of the MC-σ map set, and the shaded region corresponds to 2σ CL. |
Figure 2 clearly shows that the dipolar component (ℓ = 1) dominates the other higher-order multipoles (ℓ ≥ 2) and it is well outside of the 2σ region, which means that its statistical significance is > 95% CL. In fact, we calculate below that the dipole component
is not consistent with the statistical isotropy hypothesis at 99% CL, as compared to MC-σ maps obtained from a directional analysis applied to SI CMB maps.
On the other hand, the higher-multipole moments ,
,
are all inside the 2σ region, except for the multipole moment
, which is slightly above the 2σ region. This could be due to an incomplete sky analysis or to the application of the confidence mask (see, e.g., Bernui et al. 2007b; Rassat et al. 2014). Finally, we can conclude that the Commander-σ map obtained from the Commander separation CMB map reveals a dipole moment
that is larger than the most of dipole values in the set
that were obtained from the set of MC-σ maps (obtained from the directional analysis applied to the MC CMB maps).
Additionally, it is possible to accomplish a probabilistic analysis involving the dipole moment distribution for MC-σ map power spectra and the Commander-σ map power spectrum
. In other words, we would like to determine how likely it is to find a value
larger (or equal) than
. To know this, we constructed the
value probability distribution for all the statistically isotropic MC-σ maps and compared them with the value computed from the Commander-σ map. The results are presented in Fig. 3, where the S1 value probability distribution is represented by dark blue bars and
by a dashed line. The frequency with which S1 is higher or equal to
is small enough to reinforce our previous findings: one anomalously high value for
. This was confirmed when we computed the probability that
, which gives ∼1.0% with respect to the SI case, as shown in Fig. 3 (see also Table 1).
![]() |
Fig. 3. Distribution of the 1000 dipole moments from the MC-σ maps, |
The σ map captures the intensity of the angular correlations in all directions on the sky. By analyzing these maps, we can obtain insights about directions with stronger or weaker TT correlations. It is possible to detect any anomalous difference in the TT correlations of two different directions in this way. Figure 1 shows two almost antipode red-blue regions in the TT angular correlation distribution. It is then possible to evaluate the directions on the sky in which the correlation difference is maximum by searching the σ map pixel with the highest and lowest intensities and finding their angular coordinates. Our findings show that the maximum and minimum correlations are always approximately antipodes. In Table 1 we present the angular coordinates for the direction along which the maximum correlation difference was found for the Commander case. Furthermore, the same features appear in the other three Planck maps NILC, SEVEM, and SMICA, as observed in Fig. 4 and Table 2. This clearly shows that our results are quite robust.
![]() |
Fig. 4. From top to bottom: NILC-σ map, SEVEM-σ map, and SMICA-σ map. Observe that the high intensity regions (red pixels), corresponding to the caps with the largest TT angular correlations, are almost the same in all three maps. This interesting feature is confirmed with the measured dipole direction of these maps, quantities displayed in Table 2. |
Commander-σ map analysis results for the 192 spherical caps of radius γ0 = 60°.
Although we presented in this section the analyses considering only the Commander foreground-cleaned CMB map, in the robustness studies discussed in the next subsection, we study the other three foreground-cleaned maps, NILC, SEVEM, and SMICA. We are interested in large-angle correlations. All the CMB map data, their corresponding masks, and the MC maps were therefore downgraded to Nside = 32. This allowed us to concentrate our study of the SI on CMB maps in the angular scales ℓ = 2 − 95.
4.2. Robustness tests
In this section, we present the results of several robustness tests, divided into two groups to support our previous findings, that is, that the asymmetry phenomenon manifests itself as a large dipole component of the Commander-σ map, and its statistical significance. The first set of tests evaluated the possibility that our results might originate in a systematic effect, such as residual foregrounds coming from the foreground-cleaning processes, large-angle correlations induced by the fiducial mask applied for each CMB map in the analysis (Bernui et al. 2007b), or some residual noise present in the CMB maps. The second group of tests investigated the robustness by varying the parameters of the σ map analyses, for example, the size of spherical caps used in the analyses, γ0.
First of all, we considered the other three foreground-cleaned CMB maps to investigate the possibility that a particular pipeline in the foreground-cleaning processes used to produce the NILC, SMICA, and SEVEM CMB maps could have introduced systematic effects or left residual contaminants in the CMB maps and might thus have caused a large dipole term, , in the σ map. In all these cases, we maintained the angular radius of the spherical regions as γ0 = 60°, and we performed the same σ map analysis as with the Commander CMB map. Each one of these three maps, NILC, SEVEM, and SMICA, and the corresponding MC CMB maps were masked using the corresponding fiducial mask. Then, the σ maps were calculated, and their angular power spectra were obtained for the statistical significance analyses. The σ maps for these three foreground-cleaned CMB maps are shown in Fig. 4. The figure clearly shows regions in which the maximum and minimum correlations are almost antipodes in the three cases. In Figs. 5 and 6 and in Table 2, we summarize our results regarding these tests. As observed, the probability for
with X = NILC, SMICA, and SEVEM is also small, from ∼1.4% to ∼2.2%. Interestingly, the direction of the maximum asymmetry, that is, the dipole direction, is almost the same for these three maps, and close to the dipole direction obtained with the Commander map (see the Table 2 for a summary of our results). We conclude that the results obtained with the NILC, SMICA, and SEVEM Planck maps are statistically consistent with the analyses done with the CommanderPlanck map.
![]() |
Fig. 5. Distribution of the dipole moments from three sets of 1000 of MC-σ maps (i.e., 1000 values each for NILC, SMICA, and SEVEM), |
![]() |
Fig. 6. Angular power spectra of the NILC-σ map, SMICA-σ map, and SEVEM-σ map (dots) shown in Fig. 4, where the shaded regions correspond to 2σ CL. The continuous line corresponds to the mean angular power spectra of the MC-σ maps set. |
Directional analysis results.
In order to study the possibility that residual foreground or noise in the Planck CMB maps might affect our σ map analyses, we considered them separately. First, the dust and free-free foreground contributions6 were used to contaminate the Commander and to simulate CMB MC maps. They were then masked to perform the σ map tests. Furthermore, we considered the Planck noise maps from diverse frequencies combined with the Commander map to perform the σ map analyses. Specifically, we used the full focal-plane Planck simulations described in detail in Planck Collaboration XII (2016), where the latest simulations set is known as FFP10. These simulations contain Gaussian CMB signals and instrumental noise, and they also include, for example, physical effects of astrophysical foreground, gravitational lensing, Doppler modulation, and frequency-dependent Rayleigh scattering effects. Additionally, this data set models the Planck mission scanning strategy, detector response, beam shape, and data reduction pipeline (Planck Collaboration VII 2020). The results of our analyses with these data are summarized in Table 2.
Finally, we considered different sizes for the spherical caps, namely γ0 = 45° , 70°, and 90°, and repeated our directional analysis. Our results are also summarized in Table 2. After performing these robustness tests, we confirm that our results showing a dipolar asymmetry, or equivalently, a preferred axis, in the Planck CMB sphere are quite robust and have a confidence level of 98%−99%.
To summarize, the analyses of individual and combined frequency CMB maps performed here and for the past 20 years strongly suggest that the north-south asymmetry phenomenon is not caused by Galactic foreground, known systematics, or instrumental noise. It instead appears to be cosmological in origin (Planck Collaboration VII 2020).
5. Discussions and final remarks
In the past 20 years, many works have reported intriguing features from analyzing the large-angular properties of the CMB temperature-fluctuations maps from WMAP and Planck probes. These phenomena were generically termed CMB anomalies. Collectively, all these analyses suggest statistically that the universe is not isotropic at large scales. In particular, the north-south asymmetry phenomenon robustly reveals an axis of asymmetry on the CMB sphere that is due to anisotropic angular correlations in the Planck foreground-cleaned CMB maps. Because these results are important for the concordance model of cosmology, many analyses of the available astrophysical data are welcome.
We investigated the angular distribution of the CMB temperature fluctuations to probe the statistical isotropy of the universe. To do this, we studied the TT angular correlations through a directional analysis that we applied to the four Planck foreground-cleaned CMB maps. Specifically, our results for the Commander CMB map show that the TT angular correlations exhibit a net dipolar feature that has the highest intensity in the direction (l, b)≃(260° ,130° ) at more than 99% confidence level (see Table 2). This highly significant dipole term observed in the σ map reveals the hemispherical asymmetry in the angular correlations of the CMB sky. For this reason, this phenomenon is called the CMB north-south asymmetry. This result weakly depends on the size of the spherical caps that are used to scan the CMB sphere, as can be observed in Table 2. Most interestingly, the direction of the dipole is amazingly robust. It remains almost the same for the four foreground-cleaned maps, for different caps of radius γ0, and when foreground and noise are added to the Commander CMB map and to the MC simulated maps. This reinforces the hypothesis that the dipolar asymmetry found in our analyses has a cosmological origin, that is, it appears in all foreground-cleaned maps at large scales (γ0 > 45°), it is uncorrelated with instrumental noise and known Galactic and residual foregrounds, and it passed several robustness tests. Furthermore, the anisotropic behavior of the angular correlations we found reveals that a large part of the sky has small large-angle correlations. This is probably related to another CMB anomaly: the lack of power of the CMB large-angle correlations (Copi et al. 2006). Another interesting consequence of the fact that the asymmetry we found reveals a preferred axis in three-space is that it might be related to the CMB anomaly called the quadrupole-octopole alignment. We discuss this subject in more detail in the appendix.
More generally, several theoretical possibilities have been proposed separately to explain some CMB anomalies, in particular, the possible origin of the north-south asymmetry. This phenomenon reveals a preferred axis in three-space, the most common candidates reported in the literature are some systematic effects, such as the inhomogeneous scan of the CMB sky around the ecliptic axis, which defines a preferred axis, the mask applied to remove Galactic foregrounds, residual Galactic or extragalactic foregrounds, or a residual CMB dipole. Other studies considered modifications of the standard cosmological scenario, such as an axisymmetric primordial power spectrum, possible anisotropic inflation, primordial magnetic fields, or the topology of the universe (see, e.g., Hipólito-Ricaldi & Gomero 2005; Inoue & Silk 2006; Bernui & Hipólito-Ricaldi 2008; Schwarz et al. 2016; Pereira et al. 2007; Planck Collaboration XVIII 2016, and references therein). Nevertheless, a better understanding about the possible relation between (some) CMB anomalies seems necessary for a final answer.
Based on observations obtained with Planck (http://www.esa.int/Planck), an ESA science mission with instruments and contributions directly funded by ESA Member States, NASA, and Canada.
In practice, the interval for calculating the 2PACF ⟨0, 2γ0]=⟨0, 120° ] does not work equally well for the 192 spherical caps: Because of the Galactic mask that was applied, some caps have deficit of pairs at large angles. To solve this problem, we verified that scales smaller than ∼100° are suitable. We therefore computed the 2PACF in the interval ⟨0, 100° ] for all spherical caps.
see Table 1 in Bernui & Hipólito-Ricaldi (2008)
Acknowledgments
We acknowledge the use of data from the Planck/ESA mission, downloaded from the Planck Legacy Archive. C.E.K. and W.S.H.R. thank FAPES (PRONEM No. 503/2020) for the financial support under which this work was carried out. A.B. acknowledges a CNPq fellowship. Some of the results in this work have been derived using the HEALPix/healpy package (http://www.eso.org/science/healpix/; Górski et al. 2005).
References
- Abramo, L. R., Bernui, A., & Pereira, T. S. 2009, JCAP, 12, 013 [CrossRef] [Google Scholar]
- Aluri, P. K., & Jain, P. 2012, MNRAS, 419, 3378 [NASA ADS] [CrossRef] [Google Scholar]
- Aluri, P. K., Ralston, J. P., & Weltman, A. 2017, MNRAS, 472, 2410 [NASA ADS] [CrossRef] [Google Scholar]
- Bennett, C. L., Hill, R. S., Hinshaw, G., et al. 2003, ApJS, 148, 97 [Google Scholar]
- Bernui, A. 2008, Phys. Rev. D, 78, 063531 [NASA ADS] [CrossRef] [Google Scholar]
- Bernui, A. 2009, Phys. Rev. D, 80, 123010 [NASA ADS] [CrossRef] [Google Scholar]
- Bernui, A., & Hipólito-Ricaldi, W. S. 2008, MNRAS, 389, 1453 [NASA ADS] [CrossRef] [Google Scholar]
- Bernui, A., Villela, T., Wuensche, C. A., Leonardi, R., & Ferreira, I. 2006, A&A, 454, 409 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bernui, A., Mota, B., Rebouças, M. J., & Tavakol, R. 2007a, A&A, 464, 479 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bernui, A., Mota, B., Rebouças, M. J., & Tavakol, R. 2007b, Int. J. Mod. Phys. D, 16, 411 [NASA ADS] [CrossRef] [Google Scholar]
- Bernui, A., Ferreira, I. S., & Wuensche, C. A. 2008, ApJ, 673, 968 [NASA ADS] [CrossRef] [Google Scholar]
- Bernui, A., Oliveira, A. F., & Pereira, T. S. 2014, JCAP, 10, 041 [CrossRef] [Google Scholar]
- Bernui, A., Novaes, C. P., Pereira, T. S., & Starkman, G. D. 2018, arXiv e-prints [arXiv:1809.05924] [Google Scholar]
- Bielewicz, P., Gorski, K. M., & Banday, A. J. 2004, MNRAS, 355, 1283 [NASA ADS] [CrossRef] [Google Scholar]
- Bielewicz, P., Eriksen, H. K., Banday, A. J., Gorski, K. M., & Lilje, P. B. 2005, ApJ, 635, 750 [NASA ADS] [CrossRef] [Google Scholar]
- Calabrese, E., Slosar, A., Melchiorri, A., Smoot, G. F., & Zahn, O. 2008, Phys. Rev. D, 77, 123531 [NASA ADS] [CrossRef] [Google Scholar]
- Copi, C. J., Huterer, D., & Starkman, G. D. 2004, Phys. Rev. D, 70, 043515 [NASA ADS] [CrossRef] [Google Scholar]
- Copi, C. J., Huterer, D., Schwarz, D. J., & Starkman, G. D. 2006, MNRAS, 367, 79 [NASA ADS] [CrossRef] [Google Scholar]
- Copi, C. J., Huterer, D., Schwarz, D. J., & Starkman, G. D. 2009, MNRAS, 399, 295 [NASA ADS] [CrossRef] [Google Scholar]
- Cruz, M., Martinez-Gonzalez, E., Vielva, P., & Cayon, L. 2005, MNRAS, 356, 29 [CrossRef] [Google Scholar]
- de Oliveira-Costa, A., Tegmark, M., Zaldarriaga, M., & Hamilton, A. 2004, Phys. Rev. D, 69, 063516 [NASA ADS] [CrossRef] [Google Scholar]
- Eriksen, H. K., Hansen, F. K., Banday, A. J., Gorski, K. M., & Lilje, P. B. 2004a, ApJ, 605, 14; Erratum: 609, 1198 [NASA ADS] [CrossRef] [Google Scholar]
- Eriksen, H. K., Novikov, D. I., Lilje, P. B., Banday, A. J., & Gorski, K. M. 2004b, ApJ, 612, 64 [NASA ADS] [CrossRef] [Google Scholar]
- Gimeno-Amo, C., Barreiro, R. B., Martínez-González, E., & Marcos-Caballero, A. 2023, JCAP, 2023, 029 [CrossRef] [Google Scholar]
- Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [Google Scholar]
- Gruppuso, A., Finelli, F., Natoli, P., et al. 2011, MNRAS, 411, 1445 [CrossRef] [Google Scholar]
- Gurzadyan, V. G., Kashin, A. L., Khachatryan, H., et al. 2014, A&A, 566, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hansen, F. K., Banday, A. J., & Gorski, K. M. 2004a, MNRAS, 354, 641 [NASA ADS] [CrossRef] [Google Scholar]
- Hansen, F. K., Cabella, P., Marinucci, D., & Vittorio, N. 2004b, ApJ, 607, L67 [NASA ADS] [CrossRef] [Google Scholar]
- Hansen, M., Kim, J., Frejsel, A. M., et al. 2012, JCAP, 2012, 059 [CrossRef] [Google Scholar]
- Hipólito-Ricaldi, W. S., & Gomero, G. I. 2005, Phys. Rev. D, 72, 103008 [CrossRef] [Google Scholar]
- Hoftuft, J., Eriksen, H. K., Banday, A. J., et al. 2009, ApJ, 699, 985 [NASA ADS] [CrossRef] [Google Scholar]
- Inoue, K. T., & Silk, J. 2006, ApJ, 648, 23 [NASA ADS] [CrossRef] [Google Scholar]
- Jarosik, N., Bennett, C. L., Dunkley, J., et al. 2011, ApJS, 192, 14 [Google Scholar]
- Kim, J., & Naselsky, P. 2010, ApJ, 714, L265 [NASA ADS] [CrossRef] [Google Scholar]
- Land, K., & Magueijo, J. 2005, Phys. Rev. Lett., 95, 071301 [NASA ADS] [CrossRef] [Google Scholar]
- Marques, G. A., Novaes, C. P., Bernui, A., & Ferreira, I. S. 2018, MNRAS, 473, 165 [NASA ADS] [CrossRef] [Google Scholar]
- Mokeddem, R., Hipólito-Ricaldi, W. S., & Bernui, A. 2023, JCAP, 01, 017 [CrossRef] [Google Scholar]
- Mukherjee, S., Aluri, P. K., Das, S., Shaikh, S., & Souradeep, T. 2016, JCAP, 06, 042 [NASA ADS] [CrossRef] [Google Scholar]
- Paci, F., Gruppuso, A., Finelli, F., et al. 2010, MNRAS, 407, 399 [NASA ADS] [CrossRef] [Google Scholar]
- Pereira, T. S., Pitrou, C., & Uzan, J.-P. 2007, JCAP, 09, 006 [CrossRef] [Google Scholar]
- Pietrobon, D., Cabella, P., Balbi, A., et al. 2010, MNRAS, 402, L34 [NASA ADS] [CrossRef] [Google Scholar]
- Planck Collaboration XII. 2016, A&A, 594, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration XVIII. 2016, A&A, 594, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration IV. 2020, A&A, 641, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration VI. 2020, A&A, 641, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration VII. 2020, A&A, 641, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Polastri, L., Gruppuso, A., & Natoli, P. 2015, JCAP, 04, 018 [CrossRef] [Google Scholar]
- Rassat, A., Starck, J. L., Paykari, P., Sureau, F., & Bobin, J. 2014, JCAP, 08, 006 [CrossRef] [Google Scholar]
- Rath, P. K., & Jain, P. 2013, JCAP, 12, 014 [CrossRef] [Google Scholar]
- Rath, P. K., Mudholkar, T., Jain, P., Aluri, P. K., & Panda, S. 2013, JCAP, 04, 007 [CrossRef] [Google Scholar]
- Schwarz, D. J., Starkman, G. D., Huterer, D., & Copi, C. J. 2004, Phys. Rev. Lett., 93, 221301 [NASA ADS] [CrossRef] [Google Scholar]
- Schwarz, D. J., Copi, C. J., Huterer, D., & Starkman, G. D. 2016, Class. Quant. Grav., 33, 184001 [NASA ADS] [CrossRef] [Google Scholar]
- Vielva, P., Martinez-Gonzalez, E., Barreiro, R. B., Sanz, J. L., & Cayon, L. 2004, ApJ, 609, 22 [NASA ADS] [CrossRef] [Google Scholar]
- Yeung, S., & Chu, M.-C. 2022, Phys. Rev. D, 105, 083508 [NASA ADS] [CrossRef] [Google Scholar]
- Zhao, W. 2014, Phys. Rev. D, 89, 023010 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Correlation with other CMB anomalies
We study in this appendix the (possible) correlation between the statistically significant dipolar behavior of the σ map with the quadrupole–octopole (QO) alignment of the corresponding MC CMB map (from which we obtained the σ map). In other words, we investigated the possible cause-effect phenomenon that can explain the anomalous CMB north-south asymmetry that was revealed in the σ map analysis.
Specifically, we wish to know whether the QO alignment is sufficient to produce σ maps with the largest dipoles . To do this, we analyzed the quadrupoles and octopoles in detail, that is, ℓ = 2, 3, multipole components from the set of 50 CMB maps that produced the largest dipoles of the MC-σ
, that is, the 2σ CL of the set of 1000 dipole moments
.
First, we defined a planar mode in a multipole component. Given a multipole ℓ, the set of pairs of numbers (ℓ,m) is the mode. In a statistically isotropic CMB map, the power of the ℓ multipole, given by , is uniformly distributed among all the modes. A planar multipole means that the power of the multipole moment 𝒞ℓ is concentrated in the m = ±ℓ mode. A quadrupole, ℓ = 2, is always planar (see, e.g., Abramo et al. (2009), Schwarz et al. (2016)).
The analysis of the full set of 1000 MC CMB maps produced under the hypothesis of SI showed that the powers of the quadrupole and the octopole moments are, as expected7, uniformly distributed in their five modes (i.e., (2, 0),(2, ±1),(2, ±2)) and seven modes (i.e., (3, 0),(3, ±1),(3, ±2),(3, ±3)), respectively. Clearly, this uniform preference, or equal probability, also applies to the set of 50 MC CMB maps that produced the set . The questions is how many octopoles are planar and simultaneously have preferred planes parallel to the quadrupole plane. Because these events are not correlated, the probability for the two events to occur simultaneously is the product of the individual probabilities. According to the values in Table 1 in Bernui & Hipólito-Ricaldi (2008), we expect four to five cases8. Our search for maps with these features, that is, planar octopoles plus the QO alignment, finds four MC CMB maps. In Figure A.1 we illustrate one of these interesting cases. For the sake of completeness, we also illustrate in Figure A.2 the case in which a high σ map dipole can be obtained by coincidence9 in the angular positions of the red-blue spots.
![]() |
Fig. A.1. High σ map dipole obtained by planarity of the octopole plus the QO alignment of an MC CMB map, part of the set |
![]() |
Fig. A.2. High σ map dipole obtained by a coincidence of spots on the celestial sphere. We considered an MC CMB map, which belongs to the set |
These analyses lead us to the conclusion that the QO alignment together with the planarity of the octopole moment of the CMB map increases the probability but does not determine the north-south asymmetry phenomenon. The next question is what causes the simultaneous octopole planarity and the QO alignment. We know that both events together have a low probability of occurring in SI CMB maps.
In the past 20 years, many hypotheses have been investigated, but the origin of the CMB anomalies is not yet known. A good candidate hypothesis to explain the CMB anomalies might be spatial topology with one preferred axis: a chimney space, or a direct or π-twisted (Hipólito-Ricaldi & Gomero 2005; Bernui et al. 2018). A single-origin explanation for all CMB anomalies is still lacking so far.
All Tables
All Figures
![]() |
Fig. 1. Top: Commander-σ map, which is the result of a directional analysis of the TT angular correlations in 192 spherical caps, with a radius γ0 = 60°, applied to the CommanderPlanck CMB map. Two extremes values in almost antipodal positions are visible. Bottom: Dipole component of the Commander-σ map (top panel), which shows the dipole direction of the north-south asymmetry phenomenon: (l, b)≃(253° ,137° ). |
In the text |
![]() |
Fig. 2. Angular power spectrum of the Commander-σ map (dots), shown in Fig. 1, where the continuous blue line corresponds to the mean angular power spectrum of the MC-σ map set, and the shaded region corresponds to 2σ CL. |
In the text |
![]() |
Fig. 3. Distribution of the 1000 dipole moments from the MC-σ maps, |
In the text |
![]() |
Fig. 4. From top to bottom: NILC-σ map, SEVEM-σ map, and SMICA-σ map. Observe that the high intensity regions (red pixels), corresponding to the caps with the largest TT angular correlations, are almost the same in all three maps. This interesting feature is confirmed with the measured dipole direction of these maps, quantities displayed in Table 2. |
In the text |
![]() |
Fig. 5. Distribution of the dipole moments from three sets of 1000 of MC-σ maps (i.e., 1000 values each for NILC, SMICA, and SEVEM), |
In the text |
![]() |
Fig. 6. Angular power spectra of the NILC-σ map, SMICA-σ map, and SEVEM-σ map (dots) shown in Fig. 4, where the shaded regions correspond to 2σ CL. The continuous line corresponds to the mean angular power spectra of the MC-σ maps set. |
In the text |
![]() |
Fig. A.1. High σ map dipole obtained by planarity of the octopole plus the QO alignment of an MC CMB map, part of the set |
In the text |
![]() |
Fig. A.2. High σ map dipole obtained by a coincidence of spots on the celestial sphere. We considered an MC CMB map, which belongs to the set |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.