Scientific analyses often rely on slow, but accurate forward models for observable data conditioned on known model parameters. While various emulation schemes exist to approximate these slow calculations, these approaches are only safe if the approximations are well understood and controlled. This workshop submission reviews and updates a previously published method, which has been used in cosmological simulations, to (1) train an emulator while simultaneously estimating posterior probabilities with MCMC and (2) explicitly propagate the emulation error into errors on the posterior probabilities for model parameters. We demonstrate how these techniques can be applied to quickly estimate posterior distributions for parameters of the $\Lambda$CDM cosmology model, while also gauging the robustness of the emulator approximation.

In simple inflationary cosmological scenarios the near-exponential growth can be followed by a long period in which the Universe is dominated by the oscillating inflaton condensate. The condensate is initially almost homogeneous, but perturbations grow gravitationally, eventually fragmenting the condensate if it is not disrupted more quickly by resonance or prompt reheating. We show that the gravitational fragmentation of the condensate is well-described by the Schr\"odinger-Poisson equations and use numerical solutions to show that large overdensities form quickly after the onset of nonlinearity. This is the first exploration of this phase of nonlinear dynamics in the very early universe, which can affect the detailed form of the inflationary power spectrum and the dark matter fraction when the dark sector is directly coupled to the inflaton.

We describe the general inflationary dynamics that can arise with a single, canonically coupled field where the inflaton potential is a 4-th order polynomial. This scenario yields a wide range of combinations of the empirical spectral observables, $n_s$, $r$ and $\alpha_s$. However, not all combinations are possible and next-generation cosmological experiments have the ability to rule out all inflationary scenarios based on this potential. Further, we construct inflationary priors for this potential based on physically motivated choices for its free parameters. These can be used to determine the degree of tuning associated with different combinations of $n_s$, $r$ and $\alpha_s$ and will facilitate treatments of the inflationary model selection problem. Finally, we comment on the implications of these results for the naturalness of the overall inflationary paradigm. We argue that ruling out all simple potentials would not necessarily imply that the inflationary paradigm itself was unnatural, but that this eventuality would increase the importance of building inflationary scenarios in the context of broader paradigms of ultra-high energy physics.