In this paper, we show how the structure of the landscape potential of the primordial Universe may be probed through the properties of the primordial density perturbations responsible for the origin of the cosmic microwave background anisotropies and the large-scale structure of our Universe. Isocurvature fields -fields orthogonal to the inflationary trajectory- may have fluctuated across the barriers separating local minima of the landscape potential during inflation. We analyze how this process could have impacted the evolution of the primordial curvature perturbations. If the typical distance separating consecutive minima of the landscape potential and the height of the potential barriers are smaller than the Hubble expansion rate parametrizing inflation, the probability distribution function of isocurvature fields becomes non-Gaussian due to the appearance of bumps and dips associated to the structure of the potential. We show that this non-Gaussianity can be transferred to the statistics of primordial curvature perturbations if the isocurvature fields are coupled to the curvature perturbations. The type of non-Gaussian structure that emerges in the distribution of curvature perturbations cannot be fully probed with the standard methods of polyspectra; instead, the probability distribution function is needed. The latter is obtained by summing all the $n$-point correlation functions. To substantiate our claims, we offer a concrete model consisting of an axionlike isocurvature perturbation with a sinusoidal potential and a linear derivative coupling between the isocurvature and curvature fields. In this model, the probability distribution function of the curvature perturbations consists of a Gaussian function with small superimposed oscillations reflecting the isocurvature axion potential.

Because of their quantum fluctuations, axion fields had a chance to experience field excursions traversing many minima of their potentials during inflation. We study this situation by analyzing the dynamics of an axion-spectator field $\psi$, present during inflation, with a periodic potential given by $v(\psi) = \Lambda^4 [1 - \cos (\psi / f)]$. By assuming that the vacuum expectation value of the field is stabilized at one of its minima, say $\psi = 0$, we compute every $n$-point correlation function of $\psi$ up to first order in $\Lambda^4$ using the in-in formalism. This computation allows us to identify the distribution function describing the probability of measuring $\psi$ at a particular amplitude during inflation. Because $\psi$ is able to tunnel between the barriers of the potential, we find that the probability distribution function consists of a non-Gaussian multimodal distribution such that the probability of measuring $\psi$ at a minimum of $v(\psi)$ different from $\psi=0$ increases with time. As a result, at the end of inflation, different patches of the Universe are characterized by different values of the axion field amplitude, leading to important cosmological phenomenology: (a) Isocurvature fluctuations induced by the axion at the end of inflation could be highly non-Gaussian. (b) If the axion defines the strength of standard model couplings, then one is led to a concrete realization of the multiverse. (c) If the axion corresponds to dark matter, one is led to the possibility that, within our observable Universe, dark matter started with a nontrivial initial condition, implying novel signatures for future surveys.

Future cosmic microwave background polarization experiments will search for evidence of primordial tensor modes at large angular scales, in the multipole range $4 \leq \ell \leq 50.$ Because in that range there is some mild evidence of departures from scale invariance in the power spectrum of primordial curvature perturbations, one may wonder about the possibility of similar deviations appearing in the primordial power spectrum of tensor modes. Here we address this issue and analyze the possible presence of features in the tensor spectrum resulting from the dynamics of primordial fluctuations during inflation. We derive a general, model independent, relation linking features in the spectra of curvature and tensor perturbations. We conclude that even with large deviations from scale invariance in the curvature power spectrum, the tensor spectrum remains scale invariant for all observational purposes.