In a broad class of scenarios, inflation is followed by an extended era of matter-dominated expansion during which the inflaton condensate is nonrelativistic on subhorizon scales. During this phase density perturbations grow to the point of nonlinearity and collapse into bound structures. This epoch strongly resembles structure formation with ultra-light axion-like particles. This parallel permits us to adapt results from studies of cosmological structure formation to describe the nonlinear dynamics of this post-inflationary epoch. We show that the inflaton condensate fragments into "inflaton clusters", analogues of axion dark matter halos in present-day cosmology. Moreover, solitonic objects or "inflaton stars" can form inside these clusters, leading to density contrasts as large as $10^6$ in the post-inflationary universe.

Landscape cosmology posits the existence of a convoluted, multidimensional, scalar potential -- the "landscape" -- with vast numbers of metastable minima. Random matrices and random functions in many dimensions provide toy models of the landscape, allowing the exploration of conceptual issues associated with these scenarios. We compute the relative number and slopes of minima as a function of the vacuum energy $\Lambda$ in an $N$-dimensional Gaussian random potential, quantifying the associated probability density, $p(\Lambda)$. After normalisations $p(\Lambda)$ depends only on the dimensionality $N$ and a single free parameter $\gamma$, which is related to the power spectrum of the random function. For a Gaussian landscape with a Gaussian power spectrum, the fraction of positive minima shrinks super-exponentially with $N$; at $N=100$, $p(\Lambda>0) \approx 10^{-780}$. Likewise, typical eigenvalues of the Hessian matrices reveal that the flattest approaches to typical minima grow flatter with $N$, while the ratio of the slopes of the two flattest directions grows with $N$. We discuss the implications of these results for both swampland and conventional anthropic constraints on landscape cosmologies. In particular, for parameter values when positive minima are extremely rare, the flattest approaches to minima where $\Lambda \approx 0$ are much flatter than for typical minima, increasingly the viability of quintessence solutions.

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.

The core-cusp problem is a widely cited motivation for the exploration of dark matter models beyond standard CDM. One such alternative is ULDM; extremely light scalar particles exhibiting wavelike properties on kiloparsec scales. Astrophysically realistic ULDM halos are expected to consist of inner solitonic cores embedded in NFW-like outer halos. The presence of the solitonic core suggests that ULDM may resolve the core-cusp discrepancy associated with pure NFW halos without recourse to baryonic physics. However, it has been demonstrated that the density of ULDM halos can exceed those of comparable NFW configurations at some radii and halo masses, apparently exacerbating the problem rather than solving it. This situation arises because, although solitonic cores are flat at their centres, they obey an inverse mass-radius scaling relationship. Meanwhile, the mass of the inner soliton increases with the total halo mass, and therefore the inner core becomes more peaked at large halo masses. We describe a parameterisation of the radial density profiles of ULDM halos that allows for environmental variability of the core-halo mass relation in order to investigate this issue in more detail. For halos up to 10^12 solar masses we find feasible ULDM profiles for which the central density is lower than their NFW counterparts at astrophysically accessible radii. However, comparisons to observed profiles do not strongly favour either option; both give reasonable fits to subsets of the data for some parameter choices. Consequently, we find that robust tests of the core-cusp problem in ULDM will require more comprehensive observational data and simulations that include baryonic feedback.

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.

Tight constraints on the abundance of primordial black holes can be deduced across a vast range of masses, with the exception of those light enough to fully evaporate before nucleosynthesis. This hypothetical population is almost entirely unconstrained, to the point where the early Universe could pass through a matter-dominated phase with primordial black holes as the primary component. The only obvious relic of this phase would be Hawking radiated gravitons which would constitute a stochastic gravitational wave background in the present-day Universe, albeit at frequencies far beyond the scope of any planned detector technology. This paper explores the effects of classical mergers in such a matter dominated phase. For certain ranges of parameters, a significant fraction of the black holes merge, providing an additional, classical source of primordial gravitational waves. The resulting stochastic background typically has a lower amplitude than the Hawking background and lies at less extreme frequencies, but is unlikely to be easily detectable, with a maximal present day density of $\Omega_{GW} \sim 10^{-12}$ and frequencies between $10^{15} - 10^{19}$ Hz. We also asses the impact of radiation accretion on the lifetimes of such primordial black holes and find that it increases the black hole mass by $\sim 14 \%$ and the lifetimes by about $50 \%$. However, this does not qualitatively change any of our conclusions.

Direct measurement of acceleration is a key scientific goal for both cosmology and exoplanets. For cosmology, the concept of redshift drift (more than 60 years old by the 2020s) could directly establish the Friedmann-Lema{\^\i}tre-Robertson-Walker model. It would increase the dark energy figure of merit by a factor of 3 beyond Stage 4 experiments, in combination with cosmic microwave background measurements. For exoplanets, the same technology required provides unprecedented radial velocity accuracy, enabling detection of Earth mass planets in the habitable zone. Other science cases include mapping the Milky Way gravitational potential and testing its dark matter distribution.

The next generation of instruments designed to measure the polarization of the cosmic microwave background (CMB) will provide a historic opportunity to open the gravitational wave window to the primordial Universe. Through high sensitivity searches for primordial gravitational waves, and tighter limits on the energy released in processes like phase transitions, the CMB polarization data of the next decade has the potential to transform our understanding of the laws of physics underlying the formation of the Universe.

Typical stars in the Milky Way galaxy have velocities of hundreds of kilometres per second and experience gravitational accelerations of $\sim 10^{-10}$ m s$^{-2}$, resulting in velocity changes of a few centimetres per second over a decade. Measurements of these accelerations would permit direct tests of the applicability of Newtonian dynamics on kiloparsec length scales and could reveal significant small scale inhomogeneities within the galaxy, as well increasing the sensitivity of measurements of the overall mass distribution of the galaxy. Noting that a reasonable extrapolation of progress in exoplanet hunting spectrographs suggests that centimetre per second level precision will be attainable in the coming decade(s), we explore the possibilities such measurements would create. We consider possible confounding effects, including apparent accelerations induced by stellar motion and reflex velocities from planetary systems, along with possible strategies for their mitigation. If these issues can be satisfactorily addressed it will be possible to use high precision measurements of changing stellar velocities to perform a "blind search" for dark matter, make direct tests of theories of non-Newtonian gravitational dynamics, detect local inhomogeneities in the dark matter density, and greatly improve measurements of the overall properties of the galaxy.

PyUltraLight simulates the dynamics of ultralight dark matter in a non-expanding background. PyUltraLight can describe the evolution of several interacting ultralight dark matter halos or one or more halos orbiting a central, fixed Newtonian potential, the latter scenario corresponding to dwarf galaxies orbiting a massive central galaxy. We verify PyUltraLight by showing that it reproduces qualitative dynamical features of previously published simulations and demonstrate that it has excellent energy-conservation properties. PyUltraLight is implemented in a Python-based Jupyter notebook, solving the Schrodinger-Poisson equation governing ultralight scalar field dark matter dynamics in the non-relativistic regime using a symmetrised split-step pseudospectral algorithm. The notebook interface makes it simple to specify simulation parameters and visualise the resulting output but performance-critical routines are managed via calls to computationally efficient compiled libraries. PyUltraLight runs on standard desktop hardware with support for shared memory mutlithreading and is available on GitHub.

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.

Oscillons are spatially stationary, quasi-periodic solutions of nonlinear field theories seen in settings ranging from granular systems, low temperature condensates and early universe cosmology. We describe a new class of oscillon in which the spatial envelope can have "off centre" maxima and pulsate on timescales much longer than the fundamental frequency. These are exact solutions of the 1-D sine-Gordon equation and we demonstrate numerically that similar solutions exist in up to three dimensions for a range of potentials. The dynamics of these solutions match key properties of oscillons that may form after cosmological inflation in string-motivated monodromy scenarios.

Cosmological inflation generates primordial density perturbations on all scales, including those far too small to contribute to the cosmic microwave background. At these scales, isolated ultracompact minihalos of dark matter can form well before standard structure formation, if the perturbations have sufficient amplitude. Minihalos affect pulsar timing data and are potentially bright sources of gamma rays. The resulting constraints significantly extend the observable window of inflation in the presence of cold dark matter, coupling two of the key problems in modern cosmology.

We test the sensitivity of neutrino parameter constraints from combinations of CMB and LSS data sets to the assumed form of the primordial power spectrum (PPS) using Bayesian model selection. Significantly, none of the tested combinations, including recent high-precision local measurements of $\mathrm{H}_0$ and cluster abundances, indicate a signal for massive neutrinos or extra relativistic degrees of freedom. For PPS models with a large, but fixed number of degrees of freedom, neutrino parameter constraints do not change significantly if the location of any features in the PPS are allowed to vary, although neutrino constraints are more sensitive to PPS features if they are known a priori to exist at fixed intervals in $\log k$. Although there is no support for a non-standard neutrino sector from constraints on both neutrino mass and relativistic energy density, we see surprisingly strong evidence for features in the PPS when it is constrained with data from Planck 2015, SZ cluster counts, and recent high-precision local measurements of $\mathrm{H}_0$. Conversely combining Planck with matter power spectrum and BAO measurements yields a much weaker constraint. Given that this result is sensitive to the choice of data this tension between SZ cluster counts, Planck and $\mathrm{H}_0$ measurements is likely an indication of unmodeled systematic bias that mimics PPS features, rather than new physics in the PPS or neutrino sector.

Even simple inflationary scenarios have many free parameters. Beyond the variables appearing in the inflationary action, these include dynamical initial conditions, the number of fields, and couplings to other sectors. These quantities are often ignored but cosmological observables can depend on the unknown parameters. We use Bayesian networks to account for a large set of inflationary parameters, deriving generative models for the primordial spectra that are conditioned on a hierarchical set of prior probabilities describing the initial conditions, reheating physics, and other free parameters. We use $N_f$--quadratic inflation as an illustrative example, finding that the number of $e$-folds $N_*$ between horizon exit for the pivot scale and the end of inflation is typically the most important parameter, even when the number of fields, their masses and initial conditions are unknown, along with possible conditional dependencies between these parameters.

Cosmological analyses can be accelerated by approximating slow calculations using a training set, which is either precomputed or generated dynamically. However, this approach is only safe if the approximations are well understood and controlled. This paper surveys issues associated with the use of machine-learning based emulation strategies for accelerating cosmological parameter estimation. We describe a learn-as-you-go algorithm that is implemented in the Cosmo++ code and (1) trains the emulator while simultaneously estimating posterior probabilities; (2) identifies unreliable estimates, computing the exact numerical likelihoods if necessary; and (3) progressively learns and updates the error model as the calculation progresses. We explicitly describe and model the emulation error and show how this can be propagated into the posterior probabilities. We apply these techniques to the Planck likelihood and the calculation of $\Lambda$CDM posterior probabilities. The computation is significantly accelerated without a pre-defined training set and uncertainties in the posterior probabilities are subdominant to statistical fluctuations. We have obtained a speedup factor of $6.5$ for Metropolis-Hastings and $3.5$ for nested sampling. Finally, we discuss the general requirements for a credible error model and show how to update them on-the-fly.

A short inflationary phase may not erase all traces of the primordial universe. Associated observables include both spatial curvature and "anomalies" in the microwave background or large scale structure. The present curvature $\Omega_{K,0}$ reflects the initial curvature, $\Omega_{K,\mathrm{start}}$, and the angular size of anomalies depends on $k_\mathrm{start}$, the comoving horizon size at the onset of inflation. We estimate posteriors for $\Omega_{K,\mathrm{start}}$ and $k_\mathrm{start}$ using current data and simulations, and show that if either quantity is measured to have a non-zero value, both are likely to be observable. Mappings from $\Omega_{K,\mathrm{start}}$ and $k_\mathrm{start}$ to present-day observables depend strongly on the primordial equation of state; $\Omega_{K,0}$ spans ten orders of magnitude for a given $\Omega_{K,\mathrm{start}}$ while a simple and general relationship connects $\Omega_{K,0}$ and $k_\mathrm{start}$. We show that current bounds on $\Omega_{K,0}$ imply that if $k_\mathrm{start}$ is measurable, the curvature was already small when inflation began. Finally, since the energy density changes slowly during inflation, primordial gravitational wave constraints require that a short inflationary phase is preceded by a nontrivial pre-inflationary phase with critical implications for the expected value of $\Omega_{K,\mathrm{start}}$.

We study the tensor spectral index $n_t$ and the tensor-to-scalar ratio $r$ in the simplest multifield extension to single-field, slow-roll inflation models. We show that multifield models with potentials $V \sim \sum_i \lambda_i |\phi_i|^p$ have different predictions for $n_t/r$ than single-field models, even when all the couplings are equal $\lambda_i=\lambda_j$, due to the probabilistic nature of the fields' initial values. We analyze well-motivated prior probabilities for the $\lambda_i$ and initial conditions to make detailed predictions for the marginalized probability distribution of $n_t/r$. With $\mathcal O(100)$ fields and $p>3/4$, we find that $n_t/r$ differs from the single-field result of $n_t/r=-1/8$ at the 5$\sigma$ level. This gives a novel and testable prediction for the simplest multifield inflation models.

We present MultiModeCode, a Fortran 95/2000 package for the numerical exploration of multifield inflation models. This program facilitates efficient Monte Carlo sampling of prior probabilities for inflationary model parameters and initial conditions and is the first publicly available code that can efficiently generate large sample-sets for inflation models with $\mathcal O(100)$ fields. The code numerically solves the equations of motion for the background and first-order perturbations of multi-field inflation models with canonical kinetic terms and arbitrary potentials, providing the adiabatic, isocurvature, and tensor power spectra at the end of inflation. For models with sum-separable potentials MultiModeCode also computes the slow-roll prediction via the $\delta N$ formalism for easy model exploration and validation. We pay particular attention to the isocurvature perturbations as the system approaches the adiabatic limit, showing how to avoid numerical instabilities that affect some other approaches to this problem. We demonstrate the use of MultiModeCode by exploring a few toy models. Finally, we give a concise review of multifield perturbation theory and a user's manual for the program.

An inflationary gravitational wave background consistent with BICEP2 is difficult to reconcile with a simple power-law spectrum of primordial scalar perturbations. Tensor modes contribute to the temperature anisotropies at multipoles with $l\lesssim 100$, and this effect --- together with a prior on the form of the scalar perturbations --- was the source of previous bounds on the tensor-to-scalar ratio. We compute Bayesian evidence for combined fits to BICEP2 and Planck for three nontrivial primordial spectra: a) a running spectral index, b) a cutoff at fixed wavenumber, and c) a spectrum described by a linear spline with a single internal knot. We find no evidence for a cutoff, weak evidence for a running index, and significant evidence for a "broken" spectrum. Taken at face-value, the BICEP2 results require two new inflationary parameters in order to describe both the broken scale invariance in the perturbation spectrum and the observed tensor-to-scalar ratio. Alternatively, this tension may be resolved by additional data and more detailed analyses.

Using the temperature data from Planck we search for departures from a power-law primordial power spectrum, employing Bayesian model-selection and posterior probabilities. We parametrize the spectrum with $n$ knots located at arbitrary values of $\log{k}$, with both linear and cubic splines. This formulation recovers both slow modulations and sharp transitions in the primordial spectrum. The power spectrum is well-fit by a featureless, power-law at wavenumbers $k>10^{-3} \, \mathrm{Mpc}^{-1}$. A modulated primordial spectrum yields a better fit relative to $\Lambda$CDM at large scales, but there is no strong evidence for a departure from a power-law spectrum. Moreover, using simulated maps we show that a local feature at $k \sim 10^{-3} \, \mathrm{Mpc}^{-1}$ can mimic the suppression of large-scale power. With multi-knot spectra we see only small changes in the posterior distributions for the other free parameters in the standard $\Lambda$CDM universe. Lastly, we investigate whether the hemispherical power asymmetry is explained by independent features in the primordial power spectrum in each ecliptic hemisphere, but find no significant differences between them.

While cosmological inflation can erase primordial inhomogeneities, it is possible that inflation may not begin in a significantly inhomogeneous universe. This issue is particularly pressing in multifield scenarios, where even the homogeneous dynamics may depend sensitively on the initial configuration. This paper presents an initial survey of the onset of inflation in multifield models, via qualitative lattice-based simulations that do not include local gravitational backreaction. Using hybrid inflation as a test model, our results suggest that small subhorizon inhomogeneities do play a key role in determining whether inflation begins in multifield scenarios. Interestingly, some configurations which do not inflate in the homogeneous limit "succeed" after inhomogeneity is included, while other initial configurations which inflate in the homogeneous limit "fail" when inhomogeneity is added.

We explore whether multifield inflationary models make unambiguous predictions for fundamental cosmological observables. Focusing on $N$-quadratic inflation, we numerically evaluate the full perturbation equations for models with 2, 3, and $\mathcal{O}(100)$ fields, using several distinct methods for specifying the initial values of the background fields. All scenarios are highly predictive, with the probability distribution functions of the cosmological observables becoming more sharply peaked as $N$ increases. For $N=100$ fields, 95% of our Monte Carlo samples fall in the ranges $n_s \in (0.9455,0.9534)$; $\alpha \in (-9.741,-7.047)\times 10^{-4}$; $r\in(0.1445,0.1449)$; and $r_\mathrm{iso} \in (0.02137,3.510)\times 10^{-3}$ for the spectral index, running, tensor-to-scalar ratio, and isocurvature-to-adiabatic ratio, respectively. The expected amplitude of isocurvature perturbations grows with $N$, raising the possibility that many-field models may be sensitive to post-inflationary physics and suggesting new avenues for testing these scenarios.

Most treatments of large scale anomalies in the microwave sky are a posteriori, with unquantified look-elsewhere effects. We contrast these with physical models of specific inhomogeneities in the early universe which then generate apparent anomalies. Physical models predict correlations between candidate anomalies, as well as the corresponding signals in polarization and large scale structure, reducing the impact of cosmic variance. We compute the apparent spatial curvature associated with large-scale inhomogeneities and show that it is typically small, allowing for a self-consistent analysis. Focussing on a single large plane wave inhomogeneity, we show this can contribute to low-l mode alignment and odd-even asymmetry in the power spectra and the best fit inhomogeneity accounts for a significant part of the claimed odd-even asymmetry. We argue that this approach can be generalized to provide a more quantitative assessment of potential large scale anomalies in the universe.

Oscillons are long-lived, localized excitations of nonlinear scalar fields which may be copiously produced during preheating after inflation, leading to a possible oscillon-dominated phase in the early Universe. For example, this can happen after axion monodromy inflation, on which we run our simulations. We investigate the stochastic gravitational wave background associated with an oscillon-dominated phase. An isolated oscillon is spherically symmetric and does not radiate gravitational waves, and we show that the flux of gravitational radiation generated between oscillons is also small. However, a significant stochastic gravitational wave background may be generated during preheating itself (i.e, when oscillons are forming), and in this case the characteristic size of the oscillons is imprinted on the gravitational wave power spectrum, which has multiple, distinct peaks.

We use data from the nominal Planck mission to constrain modulations in the primordial power spectrum associated with monodromy inflation. The largest improvement in fit relative to the unmodulated model has \Delta\chi^2~10 and we find no evidence for a primordial signal, in contrast to a previous analysis of the WMAP9 dataset, for which \Delta\chi^2~20. The Planck and WMAP9 results are broadly consistent on angular scales where they are expected to agree as far as best-fit values are concerned. However, even on these scales the significance of the signal is reduced in Planck relative to WMAP, and is consistent with a fit to the noise associated with cosmic variance. Our results motivate both a detailed comparison between the two experiments and a more careful study of the theoretical predictions of monodromy inflation.

Within the Minimal Supersymmetric Standard Model (MSSM), LHC bounds suggest that scalar superpartner masses are far above the electroweak scale. Given a high superpartner mass, nonthermal dark matter is a viable alternative to WIMP dark matter generated via freezeout. In the presence of moduli fields nonthermal dark matter production is associated with a long matter dominated phase, modifying the spectral index and primordial tensor amplitude relative to those in a thermalized primordial universe. Nonthermal dark matter can have a higher self-interaction cross-section than its thermal counterpart, enhancing astrophysical bounds on its annihilation signals. We constrain the contributions to the neutralino mass from the bino, wino and higgsino using existing astrophysical bounds and direct detection experiments for models with nonthermal neutralino dark matter. Using these constraints we quantify the expected change to inflationary observables resulting from the nonthermal phase.

We investigate the initial conditions problem for multifield inflation. In these scenarios the pre-inflationary dynamics can be chaotic, increasing the sensitivity of the onset of inflation to the initial data even in the homogeneous limit. To analyze physically equivalent scenarios we compare initial conditions at fixed energy. This ensures that each trajectory is counted once and only once, since the energy density decreases monotonically. We present a full analysis of hybrid inflation that reveals a greater degree of long range order in the set of "successful" initial conditions than was previously apparent. We explore the effective smoothing scale for the fractal set of successful initial conditions induced by the finite duration of the pre-inflationary phase. The role of the prior information used to specify the initial data is discussed in terms of Bayesian sampling.

We use cosmic microwave background (CMB) data from the 9-year WMAP release to derive constraints on monodromy inflation, which is characterized by a linear inflaton potential with a periodic modulation. We identify two possible periodic modulations that significantly improve the fit, but it is unclear whether this improvement is associated with a "signal" or with the scatter in the measured angular power spectrum. The unmodulated potential is mildly favored by standard model selection criteria. A modulated inflationary potential can generate substantial primordial non-Gaussianity, of a specific and characteristic form. For the best-fit parameters to the WMAP angular power spectrum, the corresponding non-Gaussianity might be detectable in upcoming CMB data, allowing nontrivial consistency checks on the predictions of a modulated inflationary potential.

We implement Slow Roll Reconstruction -- an optimal solution to the inverse problem for inflationary cosmology -- within ModeCode, a publicly available solver for the inflationary dynamics. We obtain up-to-date constraints on the reconstructed inflationary potential, derived from the WMAP 7-year dataset and South Pole Telescope observations, combined with large scale structure data derived from SDSS Data Release 7. Using ModeCode in conjunction with the MultiNest sampler, we compute Bayesian evidence for the reconstructed potential at each order in the truncated slow roll hierarchy. We find that the data are well-described by the first two slow roll parameters, \epsilon and \eta, and that there is no need to include a nontrivial \xi parameter.

We discuss the model selection problem for inflationary cosmology. We couple ModeCode, a publicly-available numerical solver for the primordial perturbation spectra, to the nested sampler MultiNest, in order to efficiently compute Bayesian evidence. Particular attention is paid to the specification of physically realistic priors, including the parametrization of the post-inflationary expansion and associated thermalization scale. It is confirmed that while present-day data tightly constrains the properties of the power spectrum, it cannot usefully distinguish between the members of a large class of simple inflationary models. We also compute evidence using a simulated Planck likelihood, showing that while Planck will have more power than WMAP to discriminate between inflationary models, it will not definitively address the inflationary model selection problem on its own. However, Planck will place very tight constraints on any model with more than one observationally-distinct inflationary regime -- e.g. the large- and small-field limits of the hilltop inflation model -- and put useful limits on different reheating scenarios for a given model.

Oscillons are massive, long-lived, localized excitations of a scalar field. We show that in a large class of well-motivated single-field models, inflation is followed by self-resonance, leading to copious oscillon generation and a lengthy period of oscillon domination. These models are characterized by an inflaton potential which has a quadratic minimum and is shallower than quadratic away from the minimum. This set includes both string monodromy models and a class of supergravity inspired scenarios, and is in good agreement with the current central values of the concordance cosmology parameters. We assume that the inflaton is weakly coupled to other fields, so as not to quickly drain energy from the oscillons or prevent them from forming. An oscillon-dominated universe has a greatly enhanced primordial power spectrum on very small scales relative to that seen with a quadratic potential, possibly leading to novel gravitational effects in the early universe.

In a class of recently proposed models, the early universe is strongly coupled and described holographically by a three-dimensional, weakly coupled, super-renormalizable quantum field theory. This scenario leads to a power spectrum of scalar perturbations that differs from the usual empirical LCDM form and the predictions of generic models of single field, slow roll inflation. This spectrum is characterized by two parameters: an amplitude, and a parameter g related to the coupling constant of the dual theory. We estimate these parameters, using WMAP and other astrophysical data. We compute Bayesian evidence for both the holographic model and standard LCDM and find that their difference is not significant, although LCDM provides a somewhat better fit to the data. However, it appears that Planck will permit a definitive test of this holographic scenario.

Future astrophysical datasets promise to strengthen constraints on models of inflation, and extracting these constraints requires methods and tools commensurate with the quality of the data. In this paper we describe ModeCode, a new, publicly available code that computes the primordial scalar and tensor power spectra for single field inflationary models. ModeCode solves the inflationary mode equations numerically, avoiding the slow roll approximation. It is interfaced with CAMB and CosmoMC to compute cosmic microwave background angular power spectra and perform likelihood analysis and parameter estimation. ModeCode is easily extendable to additional models of inflation, and future updates will include Bayesian model comparison. Errors from ModeCode contribute negligibly to the error budget for analyses of data from Planck or other next generation experiments. We constrain representative single field models (phi^n with n=2/3, 1, 2, and 4, natural inflation, and "hilltop" inflation) using current data, and provide forecasts for Planck. From current data, we obtain weak but nontrivial limits on the post-inflationary physics, which is a significant source of uncertainty in the predictions of inflationary models, while we find that Planck will dramatically improve these constraints. In particular, Planck will link the inflationary dynamics with the post-inflationary growth of the horizon, and thus begin to probe the "primordial dark ages" between TeV and GUT scale energies.

We describe a novel approach to accelerating Monte Carlo Markov Chains. Our focus is cosmological parameter estimation, but the algorithm is applicable to any problem for which the likelihood surface is a smooth function of the free parameters and computationally expensive to evaluate. We generate a high-order interpolating polynomial for the log-likelihood using the first points gathered by the Markov chains as a training set. This polynomial then accurately computes the majority of the likelihoods needed in the latter parts of the chains. We implement a simple version of this algorithm as a patch (InterpMC) to CosmoMC and show that it accelerates parameter estimatation by a factor of between two and four for well-converged chains. The current code is primarily intended as a "proof of concept", and we argue that there is considerable room for further performance gains. Unlike other approaches to accelerating parameter fits, we make no use of precomputed training sets or special choices of variables, and InterpMC is almost entirely transparent to the user.

We consider the running of the spectral index as a probe of both inflation itself, and of the overall evolution of the very early universe. Surveying a collection of simple single field inflationary models, we confirm that the magnitude of the running is relatively consistent, unlike the tensor amplitude, which varies by orders of magnitude. Given this target, we confirm that the running is potentially detectable by future large scale structure or 21 cm observations, but that only the most futuristic measurements can distinguish between these models on the basis of their running. For any specified inflationary scenario, the combination of the running index and unknown post-inflationary expansion history induces a theoretical uncertainty in the predicted value of the spectral index. This effect can easily dominate the statistical uncertainty with which Planck and its successors are expected to measure the spectral index. More positively, upcoming cosmological experiments thus provide an intriguing probe of physics between TeV and GUT scales by constraining the reheating history associated with any specified inflationary model, opening a window into the "primordial dark age" that follows the end of inflation.

We analyze the evolution of the perturbations in the inflaton field and metric following the end of inflation. We present accurate analytic approximations for the perturbations, showing that the coherent oscillations of the post-inflationary condensate necessarily break down long before any current phenomenological constraints require the universe to become radiation dominated. Further, the breakdown occurs on length-scales equivalent to the comoving post-inflationary horizon size. This work has implications for both the inflationary "matching" problem, and the possible generation of a stochastic gravitational wave background in the post-inflationary universe.

Analytical arguments suggest that a large class of scalar field potentials permit the existence of oscillons -- pseudo-stable, non-topological solitons -- in three spatial dimensions. In this paper we numerically explore oscillon solutions in three dimensions. We confirm the existence of these field configurations as solutions to the Klein-Gorden equation in an expanding background, and verify the predictions of Amin and Shirokoff for the characteristics of individual oscillons for their model. Further, we demonstrate that significant numbers of oscillons can be generated via fragmentation of the inflaton condensate, consistent with the analysis of Amin. These emergent oscillons can easily dominate the post-inflationary universe. Finally, both analytic and numerical results suggest that oscillons are stable on timescales longer than the post-inflationary Hubble time. Consequently, the post-inflationary universe can contain an effective matter-dominated phase, during which it is dominated by localized concentrations of scalar field matter.

PSpectRe is a C++ program that uses Fourier-space pseudo-spectral methods to evolve interacting scalar fields in an expanding universe. PSpectRe is optimized for the analysis of parametric resonance in the post-inflationary universe, and provides an alternative to finite differencing codes, such as Defrost and LatticeEasy. PSpectRe has both second- (Velocity-Verlet) and fourth-order (Runge-Kutta) time integrators. Given the same number of spatial points and/or momentum modes, PSpectRe is not significantly slower than finite differencing codes, despite the need for multiple Fourier transforms at each timestep, and exhibits excellent energy conservation. Further, by computing the post-resonance equation of state, we show that in some circumstances PSpectRe obtains reliable results while using substantially fewer points than a finite differencing code. PSpectRe is designed to be easily extended to other problems in early-universe cosmology, including the generation of gravitational waves during phase transitions and pre-inflationary bubble collisions. Specific applications of this code will be pursued in future work.

A population of very light primordial black holes which evaporate before nucleosynthesis begins is unconstrained unless the decaying black holes leave stable relics. We show that gravitons Hawking radiated from these black holes would source a substantial stochastic background of high frequency gravititational waves ($10^{12}$ Hz or more) in the present universe. These black holes may lead to a transient period of matter dominated expansion. In this case the primordial universe could be temporarily dominated by large clusters of "Hawking stars" and the resulting gravitational wave spectrum is independent of the initial number density of primordial black holes.

Given a scalar field with metastable minima, bubbles nucleate quantum mechanically. When bubbles collide, energy stored in the bubble walls is converted into kinetic energy of the field. This kinetic energy can facilitate the classical nucleation of new bubbles in minima that lie below those of the "parent" bubbles. This process is efficient and classical, and changes the dynamics and statistics of bubble formation in models with multiple vacua, relative to that derived from quantum tunneling.

I review two mechanisms by which gravitational waves can be generated at the end of inflation: preheating, and gravitons Hawking radiated during the decay of very small primordial black holes. These mechanisms are contrasted with the gravitational waves during inflation, and may provide a window into the physical processes that govern the end of the inflationary phase.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We explore the consequences of the existence of a very large number of light scalar degrees of freedom in the early universe. We distinguish between participator and spectator fields. The former have a small mass, and can contribute to the inflationary dynamics; the latter are either strictly massless or have a negligible VEV. In N-flation and generic assisted inflation scenarios, inflation is a co-operative phenomenon driven by N participator fields, none of which could drive inflation on their own. We review upper bounds on N, as a function of the inflationary Hubble scale H. We then consider stochastic and eternal inflation in models with N participator fields showing that individual fields may evolve stochastically while the whole ensemble behaves deterministically, and that a wide range of eternal inflationary scenarios are possible in this regime. We then compute one-loop quantum corrections to the inflationary power spectrum. These are largest with N spectator fields and a single participator field, and the resulting bound on N is always weaker than those obtained in other ways. We find that loop corrections to the N-flation power spectrum do not scale with N, and thus place no upper bound on the number of participator fields. This result also implies that, at least to leading order, the theory behaves like a composite single scalar field. In order to perform this calculation, we address a number of issues associated with loop calculations in the Schwinger-Keldysh "in-in" formalism.

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

Slow roll reconstruction is derived from the Hamilton-Jacobi formulation of inflationary dynamics. It automatically includes information from sub-leading terms in slow roll, and facilitatesthe inclusion of priors based on the duration on inflation. We show that at low inflationary scales the Hamilton-Jacobi equations simplify considerably. We provide a new classification scheme for inflationary models, based solely on the number of parameters needed to specify the potential, and provide forecasts for likely bounds on the slow roll parameters from future datasets. A minimal running of the spectral index, induced solely by the first two slow roll parameters (\epsilon and \eta) appears to be effectively undetectable by realistic Cosmic Microwave Background experiments. However, we show that the ability to detect this signal increases with the lever arm in comoving wavenumber, and we conjecture that high redshift 21 cm data may allow tests of second order consistency conditions on inflation. Finally, we point out that the second order corrections to the spectral index are correlated with the inflationary scale, and thus the amplitude of the CMB B-mode.

Modern cosmology has sharpened questions posed for millennia about the origin of our cosmic habitat. The age-old questions have been transformed into two pressing issues primed for attack in the coming decade: How did the Universe begin? and What physical laws govern the Universe at the highest energies? The clearest window onto these questions is the pattern of polarization in the Cosmic Microwave Background (CMB), which is uniquely sensitive to primordial gravity waves. A detection of the special pattern produced by gravity waves would be not only an unprecedented discovery, but also a direct probe of physics at the earliest observable instants of our Universe. Experiments which map CMB polarization over the coming decade will lead us on our first steps towards answering these age-old questions.

We revisit arguments that simple models of inflation with a small red tilt in the scalar power spectrum generically yield an observable tensor spectrum. We show that criteria for fine-tuning based upon the algebraic simplicity of the potential depend strongly upon the explicit assumptions they incorporate, particularly regarding the end of inflation. In addition, some models with algebraically simple potentials require carefully tuned initial field configurations, and not all types of fine-tuning are identifiable via the algebraic simplicity of the potential. Conversely, in the absence of a strong prior on the mechanism that ends inflation, we demonstrate the existence of potentials with vanishingly small tensor amplitudes which are natural in terms of both their algebraic form and initial conditions. We thus argue that proposed experiments (CMBPol or BBO) which make highly sensitive measurements of the tensor amplitude cannot definitively rule out the inflationary paradigm.

We review the free parameters in the concordance cosmology, and those which might be added to this set as the quality of astrophysical data improves. Most concordance parameters encode information about otherwise unexplored aspects of high energy physics, up to the GUT scale via the "inflationary sector," and possibly even the Planck scale in the case of dark energy. We explain how neutrino properties may be constrained by future astrophysical measurements. Conversely, future neutrino physics experiments which directly measure these parameters will remove uncertainty from fits to astrophysical data, and improve our ability to determine the global properties of our universe.

We consider constraints on inflation driven by a single, minimally coupled scalar field in the light of the WMAP5 dataset, as well as ACBAR and the SuperNova Legacy Survey. We use the Slow Roll Reconstruction algorithm to derive optimal constraints on the inflationary parameter space. The scale dependence in the slope of the scalar spectrum permitted by WMAP5 is large enough to lead to viable models where the small scale perturbations have a substantial amplitude when extrapolated to the end of inflation. We find that excluding parameter values which would cause the overproduction of primordial black holes or even the onset of eternal inflation leads to potentially significant constraints on the slow roll parameters. Finally, we present a more sophisticated approach to including priors based on the total duration of inflation, and discuss the resulting restrictions on the inflationary parameter space.

We consider the impact of thermal inflation -- a short, secondary period of inflation that can arise in supersymmetric scenarios -- on the stochastic gravitational wave background. We show that while the primordial inflationary gravitational wave background is essentially unchanged at CMB scales, it is massively diluted at solar system scales and would be unobservable by a BBO style experiment. Conversely, bubble collisions at the end of thermal inflation can generate a new stochastic background. We calculate the likely properties of the bubbles created during this phase transition, and show that the expected amplitude and frequency of this signal would fall within the BBO range.

Inflation driven by a single, minimally coupled, slowly rolling field generically yields a negligible primordial non-Gaussianity. We discuss two distinct mechanisms by which a non-trivial potential can generate large non-Gaussianities. Firstly, if the inflaton traverses a feature in the potential, or if the inflationary phase is short enough so that initial transient contributions to the background dynamics have not been erased, modes near horizon-crossing can acquire significant non-Gaussianities. Secondly, potentials with small-scale structure may induce significant non-Gaussianities while the relevant modes are deep inside the horizon. The first case includes the "step" potential we previously analyzed while the second "resonance" case is novel. We derive analytic approximations for the 3-point terms generated by both mechanisms written as products of functions of the three individual momenta, permitting the use of efficient analysis algorithms. Finally, we present a significantly improved approach to regularizing and numerically evaluating the integrals that contribute to the 3-point function.

Parametric resonance or preheating is a plausible mechanism for bringing about the transition between the inflationary phase and a hot, radiation dominated universe. This epoch results in the rapid production of heavy particles far from thermal equilibrium and could source a significant stochastic background of gravitational radiation. Here, we present a numerical algorithm for computing the contemporary power spectrum of gravity waves generated in this post-inflationary phase transition for a large class of scalar-field driven inflationary models. We explicitly calculate this spectrum for both quartic and quadratic models of chaotic inflation, and low-scale hybrid models. In particular, we consider hybrid models with an ``inverted'' potential. These models have a very short and intense period of resonance which is qualitatively different from previous examples studied in this context, but we find that they lead to a similar spectrum of gravitational radiation.

We consider gravitational wave production due to parametric resonance at the end of inflation, or ``preheating''. This leads to large inhomogeneities which source a stochastic background of gravitational waves at scales inside the comoving Hubble horizon at the end of inflation. We confirm that the present amplitude of these gravitational waves need not depend on the inflationary energy scale. We analyze an explicit model where the inflationary energy scale is ~10^9 GeV, yielding a signal close to the sensitivity of Advanced LIGO and BBO. This signal highlights the possibility of a new observational ``window'' into inflationary physics, and provides significant motivation for searches for stochastic backgrounds of gravitational waves in the Hz to GHz range, with an amplitude on the order of \Omega_{gw}(k)h^2 ~ 10^-11. Finally, the strategy used in our numerical computations is applicable to the gravitational waves generated by many inhomogeneous processes in the early universe.

We compute the 3-point correlation function for a general model of inflation driven by a single, minimally coupled scalar field. Our approach is based on the numerical evaluation of both the perturbation equations and the integrals which contribute to the 3-point function. Consequently, we can analyze models where the potential has a "feature", in the vicinity of which the slow roll parameters may take on large, transient values. This introduces both scale and shape dependent non-Gaussianities into the primordial perturbations. As an example of our methodology, we examine the ``step'' potentials which have been invoked to improve the fit to the glitch in the $<TT>$ $C_l$ for $l \sim 30$, present in both the one and three year WMAP data sets. We show that for the typical parameter values, the non-Gaussianities associated with the step are far larger than those in standard slow roll inflation, and may even be within reach of a next generation CMB experiment such as Planck. More generally, we use this example to explain that while adding features to potential can improve the fit to the 2-point function, these are generically associated with a greatly enhanced signal at the 3-point level. Moreover, this 3-point signal will have a very nontrivial shape and scale dependence, which is correlated with the form of the 2-point function, and may thus lead to a consistency check on the models of inflation with non-smooth potentials.

We compute the amplitude of the non-Gaussianities in inflationary models with multiple, uncoupled scalar fields. This calculation thus applies to all models of assisted inflation, including N-flation, where inflation is driven by multiple axion fields arising from shift symmetries in a flux stabilized string vacuum. The non-Gaussianities are associated with nonlinear evolution of the field (and density) perturbations, characterized by the parameter $f_{NL}$. We derive a general expression for the nonlinear parameter, incorporating the evolution of perturbations after horizon-crossing. This is valid for arbitrary separable potentials during slow roll. To develop an intuitive understanding of this system and to demonstrate the applicability of the formalism we examine several cases with quadratic potentials: two-field models with a wide range of mass ratios, and a general N-field model with a narrow mass spectrum. We uncover that $f_{NL}$ is suppressed as the number of e-foldings grows, and that this suppression is increased in models with a broad spectrum of masses. On the other hand, we find no enhancement to $f_{NL}$ that increases with the number of fields. We thus conclude that the production of a large non-Gaussian signal in multi-field models of inflation is very unlikely as long as fields are slowly rolling and potentials are of simple, quadratic form. Finally, we compute a spectrum for the scalar spectral index that incorporates the nonlinear corrections to the fields' evolution.

We study the constraints on the inflationary parameter space derived from the 3 year WMAP dataset using ``slow roll reconstruction'', using the SDSS galaxy power spectrum to gain further leverage where appropriate. This approach inserts the inflationary slow roll parameters directly into a Monte Carlo Markov chain estimate of the cosmological parameters, and uses the inflationary flow hierarchy to compute the parameters' scale-dependence. We work with the first three parameters (epsilon, eta and xi) and pay close attention to the possibility that the 3 year WMAP dataset contains evidence for a ``running'' spectral index, which is dominated by the xi term. Mirroring the WMAP team's analysis we find that the permitted distribution of xi is broad, and centered away from zero. However, when we require that inflationary parameters yield at least 30 additional e-folds of inflation after the largest observable scales leave the horizon, the bounds on xi tighten dramatically. We make use of the absence of an explicit pivot scale in the slow roll reconstruction formalism to determine the dependence of the computed parameter distributions on the pivot. We show that the choice of pivot has a significant effect on the inferred constraints on the inflationary variables, and the spectral index and running derived from them. Finally, we argue that the next round of cosmological data can be expected to place very stringent constraints on the region of parameter space open to single field models of slow roll inflation.

We derive an extended version of the well-known Lyth Bound on the total variation of the inflaton field, incorporating higher order corrections in slow roll. We connect the field variation $\Delta\phi$ to both the spectral index of scalar perturbations and the amplitude of tensor modes. We then investigate the implications of this bound for ``small field'' potentials, where the field rolls off a local maximum of the potential. The total field variation during inflation is {\em generically} of order $m_{\rm Pl}$, even for potentials with a suppressed tensor/scalar ratio. Much of the total field excursion arises in the last e-fold of inflation and in single field models this problem can only be avoided via fine-tuning or the imposition of a symmetry. Finally, we discuss the implications of this result for inflationary model building in string theory and supergravity.

We introduce a new method for applying an inflationary prior to a cosmological dataset that includes relations between observables at arbitrary order in the slow roll expansion. The process is based on the inflationary flow equations, and the slow roll parameters appear explicitly in the cosmological parameter set. We contrast our method to other ways of imposing an inflationary prior on a cosmological dataset, and argue that this method is ideal for use with heterogeneous datasets. In particular, it would be well suited to exploiting any direct detection of fundamental tensor modes by a BBO-style mission. To demonstrate the practical use of this method we apply it to the WMAPI+All dataset, and the newly released WMAPII dataset on its own and together with the SDSS data. We find that all basic classes of single field inflationary models are still allowed at the 1-2sigma level, but the overall parameter space is sharply constrained. In particular, we find evidence that the combination of WMAPII+SDSS is sensitive to effects arising from terms that are quadratic in the two leading-order slow roll parameters.

We analyze the weak (2 sigma) evidence for a running spectral index seen in the three-year WMAP dataset and its implications for single field, slow roll inflation. We assume that the running is comparable to the central value found from the WMAP data analysis, and use the Hubble Slow Roll formalism to follow the evolution of the slow roll parameters. For all parameter choices consistent with a large, negative running, single field, slow roll inflation lasts less than 30 efolds after CMB scales leave the horizon. Thus, a definitive observation of a large negative running would imply that any inflationary phase requires multiple fields or the breakdown of slow roll. Alternatively, if single field, slow roll inflation is sources the primordial fluctuations, we can expect the observed running to move much closer to zero as the CMB is measured more accurately at small angular scales.

We consider the long standing puzzle of how to obtain meaningful probabilities in eternal inflation. We demonstrate a new algorithm to compute the probability distribution of pocket universe types, given a multivacua inflationary potential. The computed probability distribution is finite and manifestly gauge-independent. We argue that in some scenarios this technique can be applied to disfavor some eternally inflating potentials.

In many models of inflation, the period of accelerated expansion ends with preheating, a highly non-thermal phase of evolution during which the inflaton pumps energy into a specific set of momentum modes of field(s) to which it is coupled. This necessarily induces large, transient density inhomogeneities which can source a significant spectrum of gravitational waves. In this paper, we consider the generic properties of gravitational waves produced during preheating, perform detailed calculations of the spectrum for several specific inflationary models, and identify problems that require further study. In particular, we argue that if these gravitational waves exist they will necessarily fall within the frequency range that is feasible for direct detection experiments -- from laboratory through to solar system scales. We extract the gravitational wave spectrum from numerical simulations of preheating after $\lambda \phi^4$ and $m_{\phi}^2 \phi^2$ inflation, and find that they lead to a gravitational wave amplitude of around $\Omega_{gw}h^2\sim 10^{-10}$. This is considerably higher than the amplitude of the primordial gravitational waves produced during inflation. However, the typical wavelength of these gravitational waves is considerably shorter than LIGO scales, although in extreme cases they may be visible at scales accessible to the proposed BBO mission. We survey possible experimental approaches to detecting any gravitational wave background generated during preheating.

We examine the dynamics of inflation driven by multiple, interacting scalar fields and derive a multi field version of the Hubble slow roll expansion. We show that the properties of this expansion naturally generalize those of the single field case. We present an analogous hierarchy of slow roll parameters, and derive the system of ``flow'' equations that describes their evolution, and show that when this system is truncated at finite order, it can always be solved exactly. Lastly, we express the scalar and tensor perturbation spectra in terms of the slow roll parameters.

We contrast two approaches to calculating trans-Planckian corrections to the inflationary perturbation spectrum: the New Physics Hypersurface [NPH] model, in which modes are normalized when their physical wavelength first exceeds a critical value, and the Boundary Effective Field Theory [BEFT] approach, where the initial conditions for all modes are set at the same time, and modified by higher dimensional operators enumerated via an effective field theory calculation. We show that these two approaches -- as currently implemented -- lead to radically different expectations for the trans-Planckian corrections to the CMB and emphasize that in the BEFT formalism we expect the perturbation spectrum to be dominated by quantum gravity corrections for all scales shorter than some critical value. Conversely, in the NPH case the quantum effects only dominate the longest modes that are typically much larger than the present horizon size. Furthermore, the onset of the breakdown in the standard inflationary perturbation calculation predicted by the BEFT formalism is likely to be associated with a feature in the perturbation spectrum, and we discuss the observational signatures of this feature in both CMB and large scale structure observations. Finally, we discuss possible modifications to both calculational frameworks that would resolve the contradictions identified here.

Motivated by the string landscape we examine scenarios for which inflation is a two-step process, with a comparatively short inflationary epoch near the string scale and a longer period at a much lower energy (like the TeV scale). We quantify the number of $e$-foldings of inflation which are required to yield successful inflation within this picture. The constraints are very sensitive to the equation of state during the epoch between the two inflationary periods, as the extra-horizon modes can come back inside the horizon and become reprocessed. We find that the number of $e$-foldings during the first inflationary epoch can be as small as 12, but only if the inter-inflationary period is dominated by a network of cosmic strings (such as might be produced if the initial inflationary period is due to the brane-antibrane mechanism). In this case a further 20 $e$-foldings of inflation would be required at lower energies to solve the late universe's flatness and horizon problems.

We examine the constraints cosmological observations can place on any trans-Planckian corrections to the primordial spectrum of perturbations underlying the anisotropies in the Cosmic Microwave Background. We focus on models of trans-Planckian physics which lead to a modulated primordial spectrum. Rather than looking at a generic modulated spectrum, our calculations are based on a specific model, and are intended as a case study for the sort of constraints one could hope to apply on a well-motivated model of trans-Planckian physics. We present results for three different approaches -- a grid search in a subset of the overall parameter space, a Fisher matrix estimate of the likely error ellipses, and a Monte Carlo Markov Chain fit to a simulated CMB sky. As was seen in previous analyses, the likelihood space has multiple peaks, and we show that their distribution can be reproduced via a simple semi-analytic argument. All three methods lead to broadly similar results. We vary 10 cosmological parameters (including two related to the trans-Planckian terms) and show that the amplitude of the tensor perturbations is directly correlated with the detectability of any trans-Planckian modulation. We argue that this is likely to be true for any trans-Planckian modulation in the paradigm of slow-roll inflation. For the specific case we consider, we conclude that if the tensor to scalar ratio, $r \sim 0.15$, the ratio between the inflationary Hubble scale $H$, and the scale of new physics $M$ has to be on the order of 0.004 if the modulation is detectable at the 2$\sigma$ level. For a lower value of $r$, the bound on $H/M$ becomes looser.

I discuss folded inflation, an inflationary model embedded in a multi-dimensional scalar potential, such as the stringy landscape. During folded inflation, the field point evolves along a path that turns several corners in the potential. Folded inflation can lead to a relatively large tensor contribution to the Cosmic Microwave Background, while keeping all fields smaller than the Planck scale. I conjecture that if folded inflation generates a significant primordial tensor amplitude, this will generically be associated with non-trivial scale dependence in the spectrum of primordial scalar perturbations.

We analyze the cosmological consequences of locked inflation, a model recently proposed by Dvali and Kachru that can produce significant amounts of inflation without requiring slow-roll. We pay particular attention to the end of inflation in this model, showing that a secondary phase of saddle inflation can follow the locked inflationary era. However, this subsequent period of inflation results in a strongly scale dependent spectrum that can lead to massive black hole formation in the primordial universe. Avoiding this disastrous outcome puts strong constraints on the parameter space open to models of locked inflation.

We investigate the evolution of gravitational wave perturbations about a brane cosmology embedded in a five-dimensional anti-de Sitter bulk. During slow-roll inflation in a Randall-Sundrum brane-world, the zero mode of the 5-dimensional graviton is generated, while the massive modes remain in their vacuum state. When the zero mode re-enters the Hubble radius during radiation domination, massive modes are generated. We show that modes decouple in the low-energy/near-brane limit and develop perturbative techniques to calculate the mode-mixing at finite energy.

I discuss the theoretical implications of the WMAP results, stressing WMAP's detection of a correlation between the E-mode polarization and temperature anisotropies, which provides strong support for the overall inflationary paradigm. I point out that almost all inflationary models have a ``vanilla limit,'' where their parameters cannot be distinguished from a genuinely de Sitter inflationary phase. Because its findings are consistent with vanilla inflation, WMAP cannot exclude entire classes of inflationary models. Finally, I summarize hints in the current dataset that the CMB contains relics of new physics, and the possibility that we can use observational data to reconstruct the inflaton potential.

We present Monte Carlo reconstruction, a new method for ``inverting'' observational data to constrain the form of the scalar field potential responsible for inflation. This stochastic technique is based on the flow equation formalism and has distinct advantages over reconstruction methods based on a Taylor expansion of the potential. The primary ansatz required for Monte Carlo reconstruction is simply that inflation is driven by a single scalar field. We also require a very mild slow roll constraint, which can be made arbitrarily weak since Monte Carlo reconstruction is implemented at arbitrary order in the slow roll expansion. While our method cannot evade fundamental limits on the accuracy of reconstruction, it can be simply and consistently applied to poor data sets, and it takes advantage of the attractor properties of single-field inflation models to constrain the potential outside the small region directly probed by observations. We show examples of Monte Carlo reconstruction for data sets similar to that expected from the Planck satellite, and for a hypothetical measurement with a factor of five better parameter discrimination than Planck.

We derive a general expression for the power spectra of scalar and tensor fluctuations generated during inflation given an arbitrary choice of boundary condition for the mode function at a short distance. We assume that the boundary condition is specified at a short-distance cutoff at a scale $M$ which is independent of time. Using a particular prescription for the boundary condition at momentum $p = M$, we find that the modulation to the power spectra of density and gravitational wave fluctuations is of order $(H/M)$, where $H$ is the Hubble parameter during inflation, and we argue that this behavior is generic, although by no means inevitable. With fixed boundary condition, we find that the shape of the modulation to the power spectra is determined entirely by the deviation of the background spacetime from the de Sitter limit.

We analyze the impact of certain modifications to short distance physics on the inflationary perturbation spectrum. For the specific case of power-law inflation, we find distinctive -- and possibly observable -- effects on the spectrum of density perturbations.

We investigated the possibility that nonlinear gravitational effects influence the preheating era after inflation, using numerical solutions of the inhomogeneous Einstein field equations. We compared our results to perturbative calculations and to solutions of the nonlinear field equations in a rigid (unperturbed) spacetime, in order to isolate gravitational phenomena. We confirm the broad picture of preheating obtained from the nonlinear field equations in a rigid background, but find gravitational effects have a measurable impact on the dynamics. The longest modes in the simulation grow much more rapidly in the relativistic calculation than with a rigid background. We used the Weyl tensor to quantify the departure from homogeneity in the universe. We saw no evidence for the sort of gravitational collapse that leads to the formation of primordial black holes.

We show that a string-inspired Planck scale modification of general relativity can have observable cosmological effects. Specifically, we present a complete analysis of the inflationary perturbation spectrum produced by a phenomenological Lagrangian that has a standard form on large scales but incorporates a string-inspired short distance cutoff, and find a deviation from the standard result. We use the de Sitter calculation as the basis of a qualitative analysis of other inflationary backgrounds, arguing that in these cases the cutoff could have a more pronounced effect, changing the shape of the spectrum. Moreover, the computational approach developed here can be used to provide unambiguous calculations of the perturbation spectrum in other heuristic models that modify trans-Planckian physics and thereby determine their impact on the inflationary perturbation spectrum. Finally, we argue that this model may provide an exception to constraints, recently proposed by Tanaka and Starobinsky, on the ability of Planck-scale physics to modify the cosmological spectrum.

We use a numerical code to compute the density perturbations generated during an inflationary epoch which includes a spontaneous symmetry breaking phase transition. A sharp step in the inflaton potential generates $k$ dependent oscillations in the spectrum of primordial density perturbations. The amplitude and extent in wavenumber of these oscillations depends on both the magnitude and gradient of the step in the inflaton potential. We show that observations of the cosmic microwave background anisotropy place strong constraints on the step parameters.

We investigate the possibility that nonlinear gravitational effects influence the preheating era after inflation. Our work is based on numerical solutions of the inhomogeneous Einstein field equations, and is free of perturbative approximations. The one restriction we impose is to limit the inhomogeneity to a single spatial direction. We compare our results to perturbative calculations and to solutions of the nonlinear field equations in a rigid (unperturbed) spacetime, in order to isolate gravitational phenomena. We consider two types of initial conditions: where only one mode of the field perturbation has a non-zero initial amplitude, and where all the modes begin with a non-zero amplitude. Here we focus on preheating following inflation driven by a scalar field with a quartic potential. We confirm the broad picture of preheating obtained from the nonlinear field equations in a rigid background, but gravitational effects have a measurable impact on the dynamics for both sets of initial data. The rigid spacetime results predict that the amplitude of a single initially excited mode drops rapidly after resonance ends, whereas in the relativistic case the amplitude is roughly constant. With all modes initially excited, the longest modes in the simulation grow much more rapidly in the relativistic calculation than with a rigid background. However, we see no evidence for the sort of gravitational collapse associated with the formation of primordial black holes. The numerical codes described here are easily extended to more complicated resonant models, which we will examine in the future.

We study the interaction between monopoles and embedded domain walls in a O(3) linear sigma model. We discover that there is an attractive force between the monopole and the wall. We provide evidence that after the monopole and domain wall collide, the monopole unwinds on the wall, and that the winding number spreads out on the surface. These results support the suggestion by Dvali et al. (hep-ph/9710301) that the monopole problem can be solved by monopole-domain wall interactions.

We inaugurate a framework for studying preheating and parametric resonance after inflation without resorting to any approximations, either in gravitational perturbation theory or in the evolution of the field(s). We do this by numerically solving the Einstein field equations in the post-inflationary universe. In this letter we show how to compare our results to those of gauge invariant perturbation theory. We then verify Finelli and Brandenberger's analysis (hep-ph/9809490) of super-horizon modes in $m^2\phi^2$ inflation, showing that they are not amplified by resonant effects. Lastly, we make a preliminary survey of the nonlinear couplings between modes, which will be important in models where the primordial metric perturbations undergo parametric amplification.

We study spatially homogeneous and isotropic solutions to the equations of motion derived from dilaton gravity, in the presence of a special combination of higher derivative terms in the gravitational action. All solutions are nonsingular. For initial conditions resembling those in the pre-big-bang scenario, there are solutions corresponding to a spatially flat, bouncing Universe originating in a dilaton-dominated contracting phase and emerging as an expanding Friedmann Universe.

We demonstrate the existence of chaos in realistic models of two-field inflation. The chaotic motion takes place after the end of inflation, when the fields are free to oscillate and their motion is only lightly damped by the expansion of the universe. We then investigate whether the presence of chaos affects the predictions of two-field models, and show that chaos enhances the production of topological defects and renders the growth rate of the universe sensitively dependent upon the ``initial'' conditions at the beginning of the oscillatory era.

We study the cosmological implications of the one-loop terms in the string expansion. In particular, we find non-singular solutions which interpolate between a contracting universe and an expanding universe, and show that these solutions provide a mechanism for removing the initial conditions problem peculiar to spatially closed FRW cosmologies. In addition, we perform numerical calculations to show that the non-singular cosmologies do not require a careful choice of initial conditions, and estimate the likely magnitude of higher order terms in the string expansion.

We present a new, exact scalar field cosmology for which the spectrum of scalar (density) perturbations can be calculated exactly. We use this exact result to the probe the accuracy of approximate calculations of the perturbation spectrum.

We consider the spectrum of primordial fluctuations produced by inflationary models where the inflaton potential is the sum of two exponential terms. A wide range of spectra result, with the only constraint being that the scalar spectrum must have more tilt than the tensor spectrum. This model can mimic the spectra of most inflationary models, as well as producing combinations of scalar and tensor perturbations not found in common varieties of slow-rolling inflation.

We present a number of new, exact scalar field cosmologies where the potential consists of two or more exponential terms. Such potentials are motivated by supergravity or superstring models formulated in higher dimensional spacetimes that have been compactified to (3+1)-dimensions. We have found solutions in both curved and flat Robertson Walker spacetimes. These models have a diverse range of properties and often possess several distinct phases, with a smooth transition between ordinary and inflationary expansion. While exponential potentials typically produce powerlaw inflation, we find models where the inflationary period contains eras of both powerlaw and exponential growth.