Superhorizon isocurvature fluctuations relax tensions

We present a new class of models that have potential to alleviate tensions present in the cosmological data today. We postulate about variation in the sound horizon scale on superhorizon scales, i.e., on scales that are larger than that of the present observable low-redshift Universe ($\gtrsim 1\,$Gpc) while at the same time smaller than the largest scales probed by the cosmic microwave background (CMB) ($\lesssim10\,$Gpc). In this scenario, CMB peaks are naturally smoothed as preferred by the Planck data, while at the same time the low-redshift baryon acoustic oscillation calibration is partially decoupled from the CMB. Taking superhorizon variations in baryon fraction as an example and using approximate modeling, we find improvement in the best fit Planck power spectrum model $\Delta \chi^2 \sim 6$ for 1 extra degree of freedom with the relevant extension parameter $10^3 \sigma_b = 2.10 \pm 0.60$, implying about 10\% variations in baryon fraction across the Universe. At the same time, $S_8$ drops by about 1 sigma, easing tension with weak lensing surveys. We find that the combination of Planck 2018 data, eBOSS BAO data, and Riess \textit{et al.} distance ladder Hubble parameter determination produce a satisfactory fit in our model if we allow for a phantom dark energy equation of state.