Compensated Isocurvature Perturbations

• CIPs are primordial fluctuations where baryon and cold dark matter perturbations exactly cancel, leaving the total matter density unperturbed.
• They produce no linear gravitational potential effects yet generate distinctive scale-dependent signatures in galaxy clustering and higher-order CMB statistics when correlated with curvature perturbations.
• Combining galaxy surveys with kSZ tomography offers a breakthrough measurement strategy, improving constraints on CIP amplitude by orders of magnitude.

Compensated isocurvature perturbations (CIPs) are primordial fluctuations wherein a spatial modulation in the baryon density is compensated by an equal and opposite perturbation in the cold dark matter (CDM) density, such that the total matter density remains strictly unperturbed. Unlike adiabatic or conventional isocurvature modes, CIPs leave no trace in the linear-order gravitational potential and hence evade the standard constraints from the cosmic microwave background (CMB) power spectra. However, when CIPs are correlated with the primordial curvature perturbation, as predicted in several multi-field inflationary scenarios (notably the curvaton mechanism), they can generate distinctive scale-dependent signatures in galaxy clustering observables and higher-order CMB statistics, providing a unique experimental handle on the physics of the early Universe (Hotinli et al., 2019).

1. Theoretical Framework and Definition

CIPs are formally defined via entropy perturbations between baryons (b), CDM (c), and photons (γ):

• $$S_{b\gamma}(x) = \delta n_b/n_b - \delta n_\gamma/n_\gamma$$
• $$S_{c\gamma}(x) = \delta n_c/n_c - \delta n_\gamma/n_\gamma$$

A compensated mode imposes $S_{b\gamma}(x) = \Delta(x)$ and $S_{c\gamma}(x) = -(\rho_b/\rho_c) \Delta(x)$, ensuring $\rho_b\delta_b + \rho_c\delta_c = 0$—the total matter perturbation vanishes.

In multi-field models (e.g., the curvaton scenario), the CIP field $\Delta(k)$ is generically correlated with the adiabatic field $\zeta(k)$:

$\Delta(k) = A_{\rm CIP}\,\zeta(k)$

or, more generally,

$$\langle \Delta(k)\, \zeta(k') \rangle = (2\pi)^3\delta^3(k+k')\,r_{\rm CIP}A_{\rm CIP}P_\zeta(k)$$

where $A_{\rm CIP}$ sets the amplitude and $$S_{c\gamma}(x) = \delta n_c/n_c - \delta n_\gamma/n_\gamma$$0 is the correlation coefficient ($$S_{c\gamma}(x) = \delta n_c/n_c - \delta n_\gamma/n_\gamma$$1 is "fully correlated") (Hotinli et al., 2019).

The scale-invariant three-dimensional power spectrum is $$S_{c\gamma}(x) = \delta n_c/n_c - \delta n_\gamma/n_\gamma$$2, where $$S_{c\gamma}(x) = \delta n_c/n_c - \delta n_\gamma/n_\gamma$$3 is the dimensionless CIP amplitude (Kumar et al., 2022).

2. Observational Signatures in Galaxy Clustering

Scale-Dependent Galaxy Bias

In the presence of correlated CIPs, the large-scale ($$S_{c\gamma}(x) = \delta n_c/n_c - \delta n_\gamma/n_\gamma$$4 Mpc$$S_{c\gamma}(x) = \delta n_c/n_c - \delta n_\gamma/n_\gamma$$5) galaxy overdensity is modified:

$$S_{c\gamma}(x) = \delta n_c/n_c - \delta n_\gamma/n_\gamma$$6

with $$S_{c\gamma}(x) = \delta n_c/n_c - \delta n_\gamma/n_\gamma$$7 the usual linear bias. The CIP-induced shift is

$$S_{c\gamma}(x) = \delta n_c/n_c - \delta n_\gamma/n_\gamma$$8

Given the transfer function $$S_{c\gamma}(x) = \delta n_c/n_c - \delta n_\gamma/n_\gamma$$9 (relating $S_{b\gamma}(x) = \Delta(x)$0 and $S_{b\gamma}(x) = \Delta(x)$1), for fully correlated CIPs ($S_{b\gamma}(x) = \Delta(x)$2):

$S_{b\gamma}(x) = \Delta(x)$3

So, the net galaxy bias is:

$S_{b\gamma}(x) = \Delta(x)$4

Here, $S_{b\gamma}(x) = \Delta(x)$5 encodes the response of galaxy abundance to baryon–CDM fluctuation and is calculated via "separate-universe" or halo-model techniques. For an LSST-like sample, $S_{b\gamma}(x) = \Delta(x)$6 (Hotinli et al., 2019).

3. Methodology: Probing CIPs with kSZ Tomography

The kinetic Sunyaev-Zel’dovich (kSZ) effect provides an unbiased tracer of the total matter velocity field. The observed temperature fluctuation is:

$S_{b\gamma}(x) = \Delta(x)$7

Here, $S_{b\gamma}(x) = \Delta(x)$8 is the line-of-sight peculiar velocity ("remote dipole"). Cross-correlating galaxy surveys with reconstructed kSZ velocity maps allows direct measurement of the scale-dependent bias $S_{b\gamma}(x) = \Delta(x)$9. The key feature is sample-variance cancellation: since $S_{c\gamma}(x) = -(\rho_b/\rho_c) \Delta(x)$0 is an unbiased matter tracer, comparing $S_{c\gamma}(x) = -(\rho_b/\rho_c) \Delta(x)$1 and $S_{c\gamma}(x) = -(\rho_b/\rho_c) \Delta(x)$2 isolates the contribution from CIPs to $S_{c\gamma}(x) = -(\rho_b/\rho_c) \Delta(x)$3, nearly eliminating cosmic variance from the matter field (Hotinli et al., 2019).

A minimum-variance quadratic estimator $S_{c\gamma}(x) = -(\rho_b/\rho_c) \Delta(x)$4 (for redshift bin $S_{c\gamma}(x) = -(\rho_b/\rho_c) \Delta(x)$5) combines small-scale CMB and galaxy data, with reconstruction noise $S_{c\gamma}(x) = -(\rho_b/\rho_c) \Delta(x)$6 computed from their auto and cross spectra.

4. Fisher-Matrix Forecasts and Constraints

Data from LSST-like galaxy surveys and CMB-S4-class CMB experiments, incorporating both galaxy clustering and kSZ velocity fields across $S_{c\gamma}(x) = -(\rho_b/\rho_c) \Delta(x)$730 redshift bins ($S_{c\gamma}(x) = -(\rho_b/\rho_c) \Delta(x)$8), allow joint Fisher-matrix analyses on both CIPs and primordial non-Gaussianity parameters ($S_{c\gamma}(x) = -(\rho_b/\rho_c) \Delta(x)$9). The forecasted marginalized uncertainties are:

• Galaxy clustering alone: $\rho_b\delta_b + \rho_c\delta_c = 0$0
• • Planck/CMB prior: $\rho_b\delta_b + \rho_c\delta_c = 0$1
• • kSZ tomography: $\rho_b\delta_b + \rho_c\delta_c = 0$2

For comparison, Planck-only limits correspond to $\rho_b\delta_b + \rho_c\delta_c = 0$3 ($\rho_b\delta_b + \rho_c\delta_c = 0$4). Thus, next-generation kSZ tomography improves constraints by more than two orders of magnitude, enabling detection of CIPs down to the amplitude of primordial adiabatic fluctuations (Hotinli et al., 2019).

5. Degeneracy with Primordial Non-Gaussianity and Mitigation

Both local primordial non-Gaussianity ($\rho_b\delta_b + \rho_c\delta_c = 0$5) and correlated CIPs yield a scale-dependent bias $\rho_b\delta_b + \rho_c\delta_c = 0$6 on large scales, leading to potential degeneracy:

$\rho_b\delta_b + \rho_c\delta_c = 0$7

with $\rho_b\delta_b + \rho_c\delta_c = 0$8. Joint analyses allow for both $\rho_b\delta_b + \rho_c\delta_c = 0$9 and $\Delta(k)$0 to vary. When marginalizing over $\Delta(k)$1, the uncertainty on $\Delta(k)$2 degrades by about a factor of two, but sub-unity constraints on $\Delta(k)$3 remain, as the different redshift evolution and scale dependencies in $\Delta(k)$4 and $\Delta(k)$5 help break the degeneracy (Hotinli et al., 2019).

6. Extensions and Implications

Detection of correlated CIPs at or below the amplitude of adiabatic modes would constitute direct evidence for multi-field inflationary physics or nontrivial baryogenesis scenarios. The strong improvement in sensitivity from kSZ tomography plus galaxy surveys enables discrimination between early-Universe models (such as curvaton scenarios) and constrains the coupling of baryon and CDM sectors during inflation. The multi-bin, multi-tracer methodology further allows simultaneous constraints on other cosmological parameters, such as the scale-dependent signatures of $\Delta(k)$6 (Hotinli et al., 2019).

7. Summary Table: Forecasted $\Delta(k)$7 with Next-Generation Surveys

Dataset	$\Delta(k)$8	Relative Improvement
Planck CMB only	$\Delta(k)$9	Baseline
Galaxy clustering only	$\zeta(k)$0	$\zeta(k)$1
Galaxy + Planck/CMB prior	$\zeta(k)$2	$\zeta(k)$3
Galaxy + CMB + kSZ tomography	$\zeta(k)$4	$\zeta(k)$5

Forecasts are for fully correlated CIPs ($\zeta(k)$6), LSST-like galaxies, and CMB-S4 kSZ. Even after marginalizing over $\zeta(k)$7, sub-unity sensitivity is preserved (Hotinli et al., 2019).

---

Compensated isocurvature perturbations, when correlated with the primordial curvature perturbation (as in the curvaton or multi-field inflationary models), produce a unique, observable scale-dependent shift in galaxy bias. The combination of deep galaxy surveys and kSZ tomography with next-generation CMB observations yields a measurement strategy that improves the constraints on correlated CIPs by orders of magnitude, providing an incisive probe of non-adiabatic initial conditions, inflationary dynamics, and baryon–CDM physics in the early Universe (Hotinli et al., 2019).