A Universal CMB $B$-Mode Spectrum from Early Causal Tensor Sources

TL;DR

The paper identifies a universal infrared behavior for all post-inflationary, sub-horizon tensor sources: the primordial tensor power spectrum obeys $\mathcal{P}_h(k) \propto k^3$ on CMB scales due to causality and finite correlation lengths. This induces a characteristic B-mode angular spectrum that is suppressed on large angular scales and peaks at high multipoles (around $\ell \sim 10^3$), distinguishing it from inflationary predictions, while enabling a unified treatment of the stochastic GW background through the same $r_{\rm ect}$-based parametrization. The authors formalize the Early Causal Tensor Sources (ECT) framework, derive the universal IR scaling via intuitive and formal arguments, and illustrate it with three case studies—first-order phase transitions, scalar-induced GWs, and cosmic strings—showing that diverse microphysical origins share the same low-$k$ tail and observational implications. This framework provides a common language for interpreting future B-mode and SGWB measurements, offering a path to disentangle inflationary signals from sub-horizon causal sources and to place joint constraints across CMB and low-frequency GW probes.

Abstract

Many early universe scenarios predict post-inflationary tensor perturbations from causality-limited, sub-horizon sources. While the microphysical details may differ, as long as these sources are bounded in duration and correlation length, their tensor power spectra exhibit a universal scaling behavior at small wavenumber: $P_h(k) \propto k^3$, corresponding to white noise on super-horizon scales at the time of production. If these early causal tensor sources (ECTs) exclusively produce gravitational waves before redshift $z \sim 10^5$, this scaling is realized on all of the scales observed in the cosmic microwave background (CMB), and thus yields a universal multipole distribution for the $B$-mode angular power spectrum. Unlike the scale-invariant distributions of inflationary $B$ modes, ECTs generically predict enhanced power on small scales and suppressed power on large scales, which allows these source classes to be distinguished given measurements over a sufficient range of angular scales. In this paper, we introduce a unified framework for characterizing ECTs and demonstrate how their universal infrared scaling manifests in low-frequency observables, including CMB $B$ modes and stochastic gravitational wave spectral densities. We illustrate this mapping with representative case studies of this universality class involving first-order phase transitions, topological defects, and enhanced scalar perturbations, which source tensor modes at second order in perturbation theory.

Figures (5)

• Figure 1: Primordial tensor power spectrum $\mathcal{P}_h(k)$ on CMB-relevant scales (left) and resulting CMB $B$-mode spectrum $\mathcal{D}^{BB}_\ell$ (right) for early causal tensor (ECT) sources. The color map encodes the ECT amplitude parameter $r_{\rm ect}$. The solid black curves depict a scale-invariant inflationary reference spectrum with $r=0.036$, chosen to saturate the current limit from BICEP/ Keckbicepkeck21c. The lensing $B$-mode prediction for the best-fit Planck 2018 $\Lambda$CDM model is overlaid in green in the right panel. For reference, we include the upper bound of $r_{\rm ect} = 0.0077$ from our companion paper MainPaper, shown in dashed black.
• Figure 2: Predicted GW spectral density $\Omega_{\rm GW} h^2$ from ECT sources with $r_{\rm ect}\in[10^{-4},10^{-1}]$. The dashed line marks the matter-radiation equality turnover $f_{\rm eq}$. The dotted line bordering the gray region indicates the horizon scale corresponding to $z \simeq 7 \times 10^4$, beyond which the ECT signal lies outside the direct sensitivity of CMB observations. For reference, we include the $\Omega_{\rm GW}h^2$ associated with the upper bound of $r_{\rm ect} = 0.0077$ from our companion paper MainPaper, shown in dashed black.
• Figure 3: Primordial tensor power spectra $\mathcal{P}_h(k)$ sourced by a range of first-order phase transition scenarios in decoupled hidden sectors (colored curves), illustrating the universal causality-limited infrared scaling $\mathcal{P}_h(k)\propto k^3$ at low $k$. These transitions at a characteristic conformal time $\tau_{\star} \sim \tau_f$, where $\tau = 100$ kpc occurs when the Standard Model photon temperature is 1.1 keV, and $\tau_\star = 1$ kpc corresponds to a photon temperature of 110 keV. Here $\alpha$ is the ratio of false vacuum energy density to the background energy density, following the conventions in Jinno:2016Greene:2024; this parameter only affects the overall normalization of these curves. The shaded gray band indicates the CMB linear regime relevant for observational constraints.
• Figure 4: The tensor power spectrum $\mathcal{P}^{\rm SI}_h(k)$ sourced at second order by a localized enhancement in the primordial curvature power spectrum, evaluated during radiation domination using the analytic kernel $\mathcal{I}_{\rm RD}$ and a log-box ansatz for $\mathcal{P}_{\mathcal{R}}(k)$. Curves correspond to different choices of the peak wavenumber $k_p$ and amplitude $A_{\mathcal{R}}$. The shaded band indicates the CMB linear sensitivity window. On scales well below the scalar feature ($k \ll k_p$), all cases asymptote to the causality-limited white-noise tail $\mathcal{P}^{\rm SI}_h(k)\propto k^3$. Note that the dashed segments show the expected asymptotic behavior at scaling above $2 k_p$, in accordance with momentum conservation where the scalar source no longer has support to efficiently generate SIGWs.
• Figure 5: Sample dimensionless tensor power spectrum $\mathcal{P}_h(k)$ from the cosmic-string UETC calculation for three choices of the string decay time: $\tau_{\rm decay}$ = 9, 3, and 1 Mpc (red, green, blue) with three different dimensionless string tensions: $G\mu$ = $10^{-7}$, $10^{-8}$, and $10^{-9}$, respectively. The solid curves are the numerically evaluated spectra from integrating the UETC with fixed network parameters ($\xi$ = 0.13, $v$ = 0.65, $\alpha$ = 1.9) over the conformal time interval $\tau \in [\tau_i$, $\tau_{\rm decay}$], where $\tau_i$ = 0.01 Mpc. The gray band marks the approximate CMB linear-regime sensitivity window in $k$. Note that the dashed segments show the expected asymptotic power-law behavior: $k^3$ scaling on large scales and $k^{-2}$ scaling on small scales matched to the numerical results at the ends of each curve.