---
title: Number Count Dipole in Cosmology
url: https://www.emergentmind.com/topics/number-count-dipole
type: topic
---

# Number Count Dipole in Cosmology

The **number count dipole** is the $\ell=1$ anisotropy in the angular distribution of extragalactic sources on the sky, inferred from source counts rather than from temperature or intensity. In cosmology it is primarily used to test whether the matter rest frame traced by galaxies, quasars, or radio sources is consistent with the rest frame defined by the cosmic microwave background (CMB). In the standard kinematic interpretation, the dipole is generated by the observer’s motion through an otherwise statistically isotropic source distribution, but the measured signal can also contain contributions from large-scale structure, survey geometry, source evolution, and estimator-dependent systematics [2503.02470] [2212.07733].

## 1. Definition and kinematic framework

A number count dipole is usually written as a dipolar modulation of the projected counts,
\[
\frac{\Delta N}{N} = \mathbf d \cdot \hat{\mathbf n},
\]
or, in harmonic language, through the dipole power relation
\[
C_1=\frac{4\pi}{9}D^2,
\]
with $D = |\mathbf d|$ [1509.02532]. For flux-limited catalogues, the standard Ellis–Baldwin form of the **kinematic dipole** is
\[
\mathcal D_{\rm kin} = [2+x(1+\alpha)]\beta,
\qquad \beta \equiv \frac{v}{c},
\]
when the source spectrum is parameterized as $S(\nu)\propto \nu^{-\alpha}$ [2503.02470]. Other papers write
\[
\mathbf d = (2+x[1-\alpha])\frac{\mathbf v}{c}
\]
when the spectral index is defined instead by $S_\nu\propto \nu^\alpha$ [2509.16732]. This difference reflects a sign convention for $\alpha$, not a disagreement about the underlying kinematic mechanism.

The physical origin of the coefficient is standard across the literature. The term $2$ comes from aberration, while the source-count term depends on how many objects are shifted across the survey flux threshold by Doppler boosting. In several analyses the observed dipole is decomposed as
\[
\mathbf d_{\rm obs}=\mathbf d_{\rm kin}+\mathbf d_{\rm LSS},
\]
where $\mathbf d_{\rm LSS}$ is the contribution from large-scale structure [2601.07487]. A recurrent misconception is to treat the observed number count dipole as a direct velocity measurement alone; the cited decomposition makes clear that it is instead an observable sensitive to both motion and structure.

## 2. Empirical measurements and the amplitude anomaly

Wide-area surveys have produced a long sequence of dipole measurements, especially in radio continuum and quasar catalogues. A central empirical pattern is that the **direction** of the measured dipole is often broadly consistent with the CMB dipole direction, while the **amplitude** is frequently larger than the kinematic expectation inferred from the CMB.

| Survey or tracer | Reported result | Brief note |
|---|---:|---|
| NVSS, $S>30$ mJy | $v = 1110 \pm 370\ \mathrm{km\,s^{-1}}$ | Roughly $3$ times the CMB speed; direction approximately agrees with the CMB [1307.1947] |
| NVSUMSS | $1729 \pm 187\ \mathrm{km\,s^{-1}}$ hemispherical; $1355 \pm 174\ \mathrm{km\,s^{-1}}$ 3D estimator | About $4$ times larger than the CMB value; significance up to $2.81\sigma$ [1703.09376] |
| CatWISE quasars | $D/D_{\rm CMB}=2.68\pm0.23$ | Direction consistent with the CMB; amplitude anomaly quoted at $5.7\sigma$ [2212.07733] |
| NVSS + RACS-low + LoTSS-DR2 | $d=(3.67\pm0.49)\,d_{\rm exp}$ | Excess over the kinematic expectation at $5.4\sigma$ [2509.16732] |

The CMB benchmark used in these comparisons is the Solar System peculiar velocity inferred from the CMB dipole, quoted as
\[
v_{\odot,\mathrm{CMB}} = 369.82 \pm 0.11~\mathrm{km\,s^{-1}}
\]
in the radio-dipole literature [2601.07487]. If the radio excess were interpreted as purely kinematic, the implied effective velocity would be roughly
\[
v_{\rm eff}\sim 1000\text{--}1600~\mathrm{km\,s^{-1}},
\]
which several papers describe as difficult to reconcile with standard $\Lambda$CDM bulk flows [2601.07487].

Not all dipole observables point in the same direction. A blind tomographic analysis of BOSS and eBOSS spectroscopic data over $0.2<z<2.2$ measured an eBOSS redshift dipole corresponding to
\[
196^{+92}_{-79}\ \mathrm{km\,s^{-1}},
\]
which is in $2\sigma$ agreement with the CMB dipole and in $3$-to-$6\sigma$ tension with previous number count studies [2403.14580]. This contrast is central to the contemporary debate: the anomaly is robust in some catalogues and methods, but not universal across all tracers and observables.

## 3. Structure contribution, scale dependence, and redshift dependence

The intrinsic structure contribution is survey-dependent. A realistic $\Lambda$CDM analysis of the NVSS dipole found that, after removing the solar-motion contribution, about **70% of the residual dipole** arises from structures within
\[
z\lesssim 0.1,
\]
and that the signal is mainly sensitive to scales around
\[
k\sim 0.01\,h\,\mathrm{Mpc}^{-1}.
\]
The predicted amplitude depends strongly on the radio-galaxy bias. For $S>15$ mJy, the standard bias model is consistent with the observed dipole at $2.12\sigma$, whereas flux thresholds above $20$ mJy produce tensions close to $3\sigma$; adopting a constant high bias $b=3$ reduces the discrepancy to about $2.3\sigma$ [1509.02532].

Tomographic analyses strengthen the scale argument. In the WISE $\times$ SuperCOSMOS catalogue, the dipole amplitude decreases with depth, and the strongest anomaly occurs in the lowest redshift shell $0.10<z\le 0.15$. Above $z>0.15$, most bins are broadly consistent with $\Lambda$CDM mock catalogues, especially for the Fiducial sample, while the lowest-redshift shell remains persistently anomalous in both the Fiducial and SVM selections [1707.08091]. This supports the view that nearby structure and low-redshift systematics remain difficult to disentangle.

At the same time, the radio-dipole literature emphasizes that the observable is not simply a local-structure statistic. In the modified-gravity interpretation, the radio number-count dipole is described as an infrared-dominated projection of matter and velocity fields, sensitive to ultra-large wavelengths from hundreds of Mpc to gigaparsec scales, with
\[
d_{\rm LSS}^{\Lambda\mathrm{CDM}} \approx 0.1\text{--}0.3\, d_{\rm kin}
\]
as a typical standard-model expectation [2601.07487]. By contrast, the eBOSS redshift-dipole analysis estimated the intrinsic clustering dipole for its spectroscopic samples to be negligible, less than about $2\%$ of the expected kinematic effect [2403.14580]. This suggests that “number count dipole” is not a single-scale observable: the effective redshift kernel, source population, and estimator determine whether local structure, ultra-large-scale modes, or kinematics dominate.

## 4. Estimation, masking, and statistical systematics

Dipole estimation has been carried out with hemispherical asymmetry measures, 3D linear estimators, spherical-harmonic methods, tomographic least-squares fits, and pixel-based Bayesian count models. The NVSUMSS analysis used both hemispherical number-count asymmetry and a 3D linear estimator; WISE $\times$ SuperCOSMOS used the delta-map method; CatWISE analyses used a Poisson likelihood for pixel counts with explicit nuisance modeling of ecliptic-dependent modulation [1703.09376] [1707.08091] [2212.07733].

Survey realism matters quantitatively. In CatWISE, a Bayesian reanalysis found that the Galactic plane mask causes a considerable loss of dipole signal through leakage of power into higher multipoles, which **exacerbates** rather than alleviates the amplitude discrepancy [2212.07733]. In radio surveys, many sources are multi-component objects, so counts-in-cells are often overdispersed relative to a Poisson model. A Bayesian estimator based on the negative binomial distribution was introduced precisely to capture this effect, with the observed source count in cell $i$ written as
\[
N_i = \sum_{j=1}^{O_i} C_{ji}.
\]
Using NVSS, RACS-low, and LoTSS-DR2, this treatment still yielded
\[
d=(3.67\pm0.49)\,d_{\rm exp},
\]
but with larger and arguably more realistic uncertainties than under a Poisson assumption [2509.16732].

Methodological comparisons have also shown that estimator properties depend strongly on sky coverage. A comparative study of `healpy.fit_dipole` and Bayesian inference found that the estimator is unbiased in Cartesian components regardless of noise level and sky coverage, but develops large variance at low sky coverage; Bayesian analysis in the same regime is often inconclusive, which the authors present as a safeguard against overinterpretation [2509.20651]. A separate CatWISE analysis replaced the standard mean-$\alpha$, power-law-$x$ procedure by a source-by-source forward simulation using each source’s own spectral index, finding a velocity around
\[
v\approx 884\ \mathrm{km\,s^{-1}}
\]
instead of the standard
\[
717\pm85\ \mathrm{km\,s^{-1}},
\]
a shift of approximately one sigma [2505.04602]. Method dependence is therefore not incidental; it is part of the empirical status of the anomaly.

## 5. Physical interpretations and proposed resolutions

One approach is phenomenological decomposition. A joint fit to NVSS, RACS-low, and CatWISE modeled the catalogue-dependent dipole as
\[
\vec d_j = [2 + x_j(1+\alpha_j)]\vec\beta + \vec d_{\rm resid}.
\]
Under a purely kinematic hypothesis, the combined fit gave
\[
\beta=(2.62\pm0.26)\times10^{-3},
\]
equivalent to
\[
v=786\pm78\ \mathrm{km\,s^{-1}},
\]
with a common direction offset from the CMB dipole by
\[
23\pm5^\circ.
\]
When the kinematic component was fixed to the CMB expectation and a residual component was allowed, the inferred residual dipole was
\[
\mathcal D_{\rm resid} = (0.81\pm0.14)\times10^{-2},
\]
offset from the CMB dipole direction by
\[
39\pm8^\circ
\]
[2503.02470]. In this framework, the anomaly is recast as evidence for a non-kinematic contribution shared across catalogues.

A second class of explanations moves the discrepancy to the **CMB side** rather than the matter side. A detailed analysis of superhorizon perturbations concluded that neither adiabatic nor isocurvature superhorizon modes can generate an intrinsic galaxy number-count dipole at observable redshift, although a superhorizon isocurvature mode can induce a matter–radiation relative velocity and thereby modify the CMB dipole [2207.01569]. A QCD axion model then developed this logic further: the CMB dipole is written as
\[
d_{\rm CMB}=d_{\rm kin}^{\rm CMB}+d_{\rm int}^{\rm CMB},
\]
with the intrinsic term generated by a super-horizon isocurvature perturbation. In that scenario the galaxy number-count dipole remains very small, of order
\[
\sim 6\times10^{-8},
\]
while the CMB dipole can receive an intrinsic contribution large enough to address the mismatch [2211.06912]. These proposals do not explain the anomaly as an intrinsic number-count dipole; they reinterpret the CMB dipole as not being purely kinematic.

A third line of work modifies the late-time growth of structure. The STVG-MOG proposal keeps early-universe and small-scale behavior close to standard cosmology but introduces a scale-dependent effective coupling
\[
G_{\rm eff}(k,a)=G_N[1+\alpha_{\rm eff}(k,a)],
\qquad
\alpha_{\rm eff}(k,a)=\alpha\frac{\mu^2}{k^2/a^2+\mu^2}.
\]
In this model, gravity becomes effectively Newtonian on small scales and enhanced on ultra-large scales. Because the velocity field contains an explicit $1/k$ weighting,
\[
P_v(k,a)=\left(\frac{aH f(k,a)}{k}\right)^2 P_\delta(k,a),
\]
even modest large-scale growth enhancement can strongly amplify coherent bulk flows and dipole anisotropies. The dipole enhancement is summarized schematically by
\[
R_d \simeq [1+\alpha_{\rm eff}(k_\ast,a_\ast)]^p,
\]
and for $\alpha=\mathcal{O}(5\text{--}10)$ the model is argued to permit a factor-of-few enhancement, potentially enough to turn a $\Lambda$CDM-like dipole into the observed $3$–$4\times$ excess [2601.07487]. This interpretation treats the radio dipole anomaly as a physical late-time structure signal rather than as an anomalous Solar System velocity.

## 6. Extensions, related observables, and future probes

The number count dipole is no longer restricted to broadband continuum catalogues. For line-emitting galaxies, the cumulative number density
\[
\mathcal N(\nu)\equiv \frac{dN}{d\nu\,d\Omega}
\]
has a dipole
\[
\mathcal N'_{\rm dip}(\nu')=\mathcal N(\nu')\left(1-\left.\frac{\partial\ln \mathcal N}{\partial\ln \nu}\right|_{\nu'}\right)\beta.
\]
Here the amplitude depends not only on $\beta$ but also on the frequency evolution of the monopole counts, so different observing frequencies provide independent redshift slices of the same kinematic signal. On this basis, SPHEREx was identified as the optimal near-term survey for a first detection of the spectral line galaxy number-count dipole, while line-intensity mapping surveys were also found promising [2501.09800].

Several related observables probe similar physics but are not identical to the monopole number-count dipole. In the galaxy bispectrum, the Newtonian single-tracer signal has no dipole, whereas relativistic light-cone corrections generate odd multipoles, including a bispectrum dipole that can exceed $10\%$ of the monopole on sufficiently large scales [1812.09512]. In the two-point function of number-count fluctuations, the Doppler imprint on the covariance matrix can be described with bipolar spherical harmonics; forecasts indicate that radio continuum surveys with SKA Phase 2 should detect this effect at $\gtrsim 30\sigma$, while SKA Phase 1 allows only marginal detection in optimistic cases [1808.09743].

Independent dipole probes are also being developed outside electromagnetic number counts. For binary black hole detections with third-generation gravitational-wave observatories, the number-count dipole simplifies to
\[
\mathcal D \simeq 2\,\frac{v_o}{c}
\]
when threshold effects are negligible, making the signal essentially independent of astrophysical population details. Forecasts indicate that a dipole as large as the radio-galaxy value would be detectable at $>3\sigma$ with about $10^6$ events, while a CMB-sized dipole would require at least $10^7$ detections [2209.11658]. Strongly lensed gravitational waves combined with galaxy surveys provide a second route: in an optimistic 10-year scenario, the combined constraint was forecast as
\[
g = (2.45^{+1.53}_{-1.28}) \times 10^{-3},
\]
provided systematic uncertainties are mitigated [2605.19476]. In parallel, the modified-gravity interpretation explicitly identifies future SKA measurements as a potential discriminator between standard and scale-enhanced structure-growth pictures [2601.07487].

Taken together, the literature presents the number count dipole as both a kinematic observable and a precision test of large-scale isotropy. Its current significance lies less in the existence of a dipole itself—which is expected—and more in the unsettled question of whether the excess amplitude reported in radio and quasar catalogues is due to local structure, survey methodology, a residual non-kinematic anisotropy, an intrinsic CMB component, or new late-time gravitational physics.

Source: https://www.emergentmind.com/topics/number-count-dipole