The CatWISE2020 Quasar dipole: A Reassessment of the Cosmic Dipole Anomaly

TL;DR

This paper reevaluates the cosmic dipole anomaly using the CatWISE2020 quasar catalog, reducing its significance from 4.9σ to around 3.3σ by accounting for uncertainties and survey effects.

Contribution

It provides a comprehensive reassessment of the quasar dipole anomaly, incorporating detailed simulations and higher-order multipole analysis to better understand its significance.

Findings

The dipole significance decreases from 4.9σ to approximately 3.3σ after accounting for uncertainties.

Survey geometry and clustering dipole effects cannot fully explain the observed anomaly.

Including higher multipoles introduces mode coupling and biases in dipole estimation on partial sky data.

Abstract

The Ellis-Baldwin test probes the cosmological principle by comparing the kinematic Cosmic Microwave Background dipole with the Doppler-driven dipole in the number counts of extragalactic radio sources. Recent analysis of the CatWISE2020 quasar catalog reported a number-count dipole amplitude exceeding the kinematic expectation at $4.9\sigma$ significance. We present a comprehensive reassessment of this test using the same dataset, incorporating major sources of uncertainty in the statistical inference. We employ a simulation framework based on the FLASK package, using lognormal realizations of the large-scale structure, quasar clustering bias, the survey's radial selection function, and its exact sky coverage. Our simulations account for the kinematic dipole, the intrinsic clustering dipole, shot noise, and survey geometry effects. The analysis yields a revised significance of $3.63\sigma$ in the absence of a clustering dipole, and $3.44\sigma$ with a randomly oriented clustering dipole. When the clustering dipole is aligned with the kinematic dipole, the significance decreases further to $3.27\sigma$. Although the anomaly is reduced, it cannot be explained solely by the clustering dipole or mode coupling from the survey mask. We further assess dipole measurement robustness by fitting models with successively higher-order multipoles up to $\ell = 4$. Partial sky coverage induces mode coupling, shifting the dipole estimate to higher values when the octopole is included and inflating its variance as additional modes are incorporated, reflected in the increasing condition number of the estimator. This behavior highlights a bias-variance trade-off inherent in multipole fitting on partial-sky data.