Electromagnetic waves in a dynamical axion background exhibit superluminal group velocities at high frequencies and instabilities at low frequencies, altering how photons propagate through space. Local disturbances propagate causally, but unlike in ordinary Maxwell theory, propagation occurs inside as well as on the lightcone. For the unstable modes, the energy density in the electromagnetic field grows exponentially along timelike displacements. In this paper we derive retarded Green functions in axion electrodynamics in various limits and study the time-domain properties of propagating signals.

During inflation, scalar fields with masses less than the Hubble scale acquire vacuum expectation values (vevs) via stochastic processes driven by quantum fluctuations. For nearly massless spectator scalars transforming nontrivially under a continuous symmetry group, we demonstrate that the evolution of the vev depends on the dimensionality of the scalar field space. Fields in larger representations both attain larger vacuum expectation values and converge more rapidly to equilibrium. We present an argument demonstrating how this higher-dimensional evolution can be obtained in unitary gauge for fields transforming under local symmetries with a mass gap that is small compared to the Hubble scale. Finally, we show that accounting for the full number of degrees of freedom in the Standard Model Higgs multiplet tightens Higgs stability constraints on the inflationary scale at the percent level and has more dramatic consequences for both the vev and the energy stored in the Higgs field after inflation.

The existence of dark radiation that is completely decoupled from the standard model in the early Universe leaves open the possibility of an associated dark radiation isocurvature mode. We show that the presence of dark radiation isocurvature leads to spatial variation in the primordial abundances of helium and deuterium due to spatial variation in $N_{\rm eff}$ during Big Bang nucleosynthesis. We use the result to constrain the existence of such an isocurvature mode on scales down to $\sim 1$ Mpc scales. By measuring the excess variance in the primordial helium to hydrogen and deuterium to hydrogen ratio in different galaxies, we constrain the variance in average isocurvature in a galaxy to be less than $0.13/\Delta \bar{N}_{\rm eff}$ at 95\% confidence. Here $\Delta \bar{N}_{\rm eff}$ is the spatially averaged increase in $N_{\rm eff}$ due to the additional dark radiation component.

Characterizing the physical properties of the stochastic gravitational waves background (SGWB) is a key step towards identifying the nature of its possible origin. We focus our analysis on SGWB anisotropies. The existence of a non-trivial primordial scalar-tensor-tensor (STT) correlation in the squeezed configuration may be inferred from the effect that a long wavelength scalar mode has on the gravitational wave power spectrum: an anisotropic contribution. Crucially, such contribution is correlated with temperature anisotropies in the cosmic microwave background (CMB). We show that, for inflationary models that generate suitably large STT non-Gaussianities, cross-correlating the CMB with the stochastic background of gravitational waves is a very effective probe of early universe physics. The resulting signal can be a smoking-gun for primordial SGWB anisotropies.

We study the production of fermions through a derivative coupling with a pseudoscalar inflaton and the effects of the produced fermions on the scalar primordial perturbations. We present analytic results for the modification of the scalar power spectrum due to the produced fermions, and we estimate the amplitude of the non-Gaussianities in the equilateral regime. Remarkably, we find a regime where the effect of the fermions gives the dominant contribution to the scalar spectrum while the amplitude of the bispectrum is small and in agreement with observation. We also note the existence of a regime in which the backreaction of the fermions on the evolution of the zero-mode of the inflaton can lead to inflation even if the potential of the inflaton is steep and does not satisfy the slow-roll conditions.

We study gravitational wave production from gauge preheating in a variety of inflationary models, detailing its dependence on both the energy scale and the shape of the potential. We show that gauge preheating generically leads to a large gravitational wave background that contributes significantly to the effective number of relativistic degrees of freedom in the early Universe, $N_\mathrm{eff}$. We demonstrate that the efficiency of gravitational wave production is correlated with the tensor-to-scalar ratio, $r$. In particular, we show that efficient gauge preheating in models whose tensor-to-scalar ratio would be detected by next-generation cosmic microwave background experiments ($r \gtrsim 10^{-3}$) will either be detected through its contribution to $N_\mathrm{eff}$ or ruled out. Furthermore, we show that bounds on $N_\mathrm{eff}$ provide the most sensitive probe of the possible axial coupling of the inflaton to gauge fields regardless of the potential.

We demonstrate that gravitational waves generated by efficient gauge preheating after axion inflation generically contribute significantly to the effective number of relativistic degrees of freedom $N_\mathrm{eff}$. We show that, with existing Planck limits, gravitational waves from preheating already place the strongest constraints on the inflaton's possible axial coupling to gauge fields. We demonstrate that gauge preheating can completely reheat the Universe regardless of the inflationary potential. Further, we quantify the variation of the efficiency of gravitational wave production from model to model and show that it is correlated with the tensor-to-scalar ratio. In particular, when combined with constraints on models whose tensor-to-scalar ratios would be detected by next-generation cosmic microwave background experiments, $r\gtrsim 10^{-3}$, constraints from $N_\mathrm{eff}$ will probe or rule out the entire coupling regime for which gauge preheating is efficient.

We analyze reheating scenarios where a hidden sector is populated during reheating along with the sector containing the Standard Model. We numerically solve the Boltzmann equations describing perturbative reheating of the two sectors, including the full dependence on quantum statistics, and study how quantum statistical effects during reheating as well as the non-equilibrium inflaton-mediated energy transfer between the two sectors affects the temperature evolution of the two radiation baths. We obtain new power laws describing the temperature evolution of fermions and bosons when quantum statistics are important during reheating. We show that inflaton-mediated scattering is generically most important at radiation temperatures $T\sim M_\phi/4$, and build on this observation to obtain analytic estimates for the temperature asymmetry produced by asymmetric reheating. We find that for reheating temperatures $T_{\textrm{rh}} \ll M_{\phi}/4$, classical perturbative reheating provides an excellent approximation to the final temperature asymmetry, while for $T_{\textrm{rh}}\gg M_{\phi}/4$, inflaton-mediated scattering dominates the population of the colder sector and thus the final temperature asymmetry. We additionally present new techniques to calculate energy transfer rates between two relativistic species at different temperatures.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

We present analytic results for the gravitational wave power spectrum induced in models where the inflaton is coupled to a fermionic pseudocurrent. We show that although such a coupling creates helically polarized fermions, the polarized component of the resulting gravitational waves is parametrically suppressed with respect to the non-polarized one. We also show that the amplitude of the gravitational wave signal associated to this production cannot exceed that generated by the standard mechanism of amplification of vacuum fluctuations. We previously found that this model allows for a regime in which the backreaction of the produced fermions allows for slow-roll inflation even for a steep inflaton potential, and still leads to Gaussian primordial scalar perturbations. The present analysis shows that this regime also results in a gravitational wave signal compatible with the current bounds.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

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.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

We study gravitational wave production during Abelian gauge-field preheating following inflation. We consider both scalar and pseudoscalar inflaton models coupled directly to Abelian gauge fields via either a dilatonic coupling to the gauge-field kinetic term or an axial coupling to a Chern-Simons term. In both cases gravitational waves are produced efficiently during the preheating phase, with a signature louder than most cosmological signals. These gravitational waves can contribute to the radiation energy budget of Universe at a level which will be probed by upcoming cosmic microwave background experiments through $N_{\rm eff}$. For axially coupled fields the resulting gravitational wave spectrum is helically polarized---a unique feature that can be used to differentiate it from other stochastic gravitational wave backgrounds. We compute the gravitational topological charge and demonstrate that gauge preheating following axion inflation may be responsible for the matter-antimatter asymmetry of the Universe via gravitational leptogenesis.

We study anisotropic inflationary solutions in massive Gauge-flation. We work with the theory in both the Stueckelberg and dynamical symmetry-breaking limits and demonstrate that extended periods of accelerated anisotropic expansion are possible. In the case of dynamical symmetry breaking, we show that spacetime can transition from isotropic quasi-de Sitter space to an accelerating Bianchi spacetime due to a rolling Higgs field - the spacetime can develop hair. Similarly, symmetry restoring transitions are possible from accelerating Bianchi spacetime to quasi-de Sitter space - the spacetime can lose its hair. These transitions can be arranged to occur quickly, within an $e$-folding or so, or over tens of $e$-folds.

We develop an integral form for the bispectrum in general single-field inflation whose domain of validity includes models of inflation where the background evolution is not constrained to be slowly varying everywhere. Our integral form preserves the squeezed-limit consistency relation, allows for fast evaluation of the bispectrum for all triangle configurations expediting the efficient comparison of slow-roll violating models with data, and provides complete and compact slow-roll expressions correct to first order in slow-roll parameters. Motivated by the recent Planck results, we consider as an example a sharp step in the warped-brane tension of DBI inflation and provide analytic solutions for the peak of the resulting bispectrum. For the step in the warp that reproduces the oscillations in the power spectrum favored by the Planck data, the corresponding equilateral bispectrum is both extremely large and highly scale dependent. The bispectrum serves as a means of distinguishing such a model from alternative scenarios that generate otherwise indistinguishable power spectra, such as a step in the potential in canonical single-field inflation.

We study preheating in models where a scalar inflaton is directly coupled to a non-Abelian $SU(2)$ gauge field. In particular, we examine $m^2 \phi^2$ inflation with a conformal, dilaton-like coupling to the non-Abelian sector. We describe a numerical scheme that combines lattice gauge theory with standard finite difference methods applied to the scalar field. We show that a significant tachyonic instability allows for efficient preheating, which is parametrically suppressed by increasing the non-Abelian self-coupling. Additionally, we comment on the technical implementation of the evolution scheme and setting initial conditions.

Gravitational leptogenesis refers to a class of baryogenesis models in which the matter-antimatter asymmetry of the universe arises through the standard model lepton-number gravitational anomaly. In these models chiral gravitational waves source a lepton asymmetry in standard model neutrinos during the inflationary epoch. We point out that gravitational leptogenesis can be successful in either the Dirac or Majorana neutrino mass scenario. In the Dirac mass scenario, gravitational leptogenesis predicts a relic abundance of sterile neutrinos that remain out of equilibrium, and the lepton asymmetry carried by the standard model sector is unchanged. In the Majorana mass scenario, the neutrinos participate in lepton-number-violating interactions that threaten to washout the lepton asymmetry during post-inflationary reheating. However, we show that a complete (exponential) washout of the lepton asymmetry is prevented if the lepton-number-violating interactions go out of equilibrium before all of the standard model Yukawa interactions come into equilibrium. The baryon and lepton asymmetries carried by right-chiral quarks and leptons are sequestered from the lepton-number violation, and the washout processes only suppress the predicted baryon asymmetry by a factor of $\varepsilon_{\rm w.o.} = \pm O(0.1)$. The sign of $\varepsilon_{\rm w.o.}$ depends on the model parameters in such a way that a future measurement of the primordial gravitational wave chirality would constrain the scale of lepton-number violation (heavy Majorana neutrino mass).

We elucidate the counting of the relevant small parameters in inflationary perturbation theory. Doing this allows for an explicit delineation of the domain of validity of the semi-classical approximation to gravity used in the calculation of inflationary correlation functions. We derive an expression for the dependence of correlation functions of inflationary perturbations on the slow-roll parameter $\epsilon = -\dot{H}/H^2$, as well as on $H/M_p$, where $H$ is the Hubble parameter during inflation. Our analysis is valid for single-field models in which the inflaton can traverse a Planck-sized range in field values and where all slow-roll parameters have approximately the same magnitude. As an application, we use our expression to seek the boundaries of the domain of validity of inflationary perturbation theory for regimes where this is potentially problematic: models with small speed of sound and models allowing eternal inflation.

We study a variant of Gauge-flation where the gauge symmetry is spontaneously broken by a Higgs sector. We work in the Stueckelberg limit and demonstrate that the dynamics remain (catastrophically) unstable for cases where the gauge field masses satisfy $\gamma < 2$, where $\gamma = g^2\psi^2/H^2$, $g$ is the gauge coupling, $\psi$ is the gauge field vacuum expectation value, and $H$ is the Hubble rate. We compute the spectrum of density fluctuations and gravitational waves, and show that the model can produce observationally viable spectra. The background gauge field texture violates parity, resulting in a chiral gravitational wave spectrum. This arises due to an exponential enhancement of one polarization of the spin-2 fluctuation of the gauge field. Higgsed Gauge-flation can produce observable gravitational waves at inflationary energy scales well below the GUT scale.

We study derivatively coupled fermions in axion-driven inflation, specifically $m_\phi^2\phi^2$ and monodromy inflation, and calculate particle production during the inflationary epoch and the post-inflationary axion oscillations. During inflation, the rolling axion acts as an effective chemical potential for helicity which biases the gravitational production of one fermion helicity over the other. This mechanism allows for efficient gravitational production of heavy fermion states that would otherwise be highly suppressed. Following inflation, the axion oscillates and fermions with both helicities are produced as the effective frequency of the fermion field changes non-adiabatically. For certain values of the fermion mass and axion-fermion coupling strength, the two helicity states are produced asymmetrically, resulting in unequal number-densities of left- and right-helicity fermions.

We study the onset of the reheating epoch at the end of axion-driven inflation where the axion is coupled to an Abelian, $U(1)$, gauge field via a Chern-Simons interaction term. We focus primarily on $m^2\phi^2$ inflation and explore the possibility that preheating can occur for a range of coupling values consistent with recent observations and bounds on the overproduction of primordial black holes. We find that for a wide range of parameters preheating is efficient. In certain cases the inflaton is seen to transfer all its energy to the gauge fields within a few oscillations. In most cases, we find that the gauge fields on sub-horizon scales end preheating in an unpolarized state due to the existence of strong rescattering between the inflaton and gauge-field modes. We also present a preliminary study of an axion monodromy model coupled to $U(1)$ gauge fields, seeing a similarly efficient preheating behavior as well as indications that the coupling strength has an effect on the creation of oscillons.

In this work we compute the production of magnetic fields in models of axion inflation coupled to the hypercharge sector of the Standard Model through a Chern-Simons interaction term. We make the simplest choice of a quadratic inflationary potential and use lattice simulations to calculate the magnetic field strength, helicity and correlation length at the end of inflation. For small values of the axion-gauge field coupling strength the results agree with no-backreaction calculations and estimates found in the literature. For larger couplings the helicity of the magnetic field differs from the no-backreaction estimate and depends strongly on the comoving wavenumber. We estimate the post-inflationary evolution of the magnetic field based on known results for the evolution of helical and non-helical magnetic fields. The magnetic fields produced by axion inflation with large couplings to $U(1)_Y$ can reach $B_{\rm eff} \gtrsim 10^{-16}\, {\rm G}$, exhibiting a field strength $B_{\rm phys} \approx 10^{-13}\, {\rm G}$ and a correlation length $\lambda_{\rm phys}\approx10\, {\rm pc}$. This result is insensitive to the exact value of the coupling, as long as the coupling is large enough to allow for instantaneous preheating. Depending on the assumptions for the physical processes that determine blazar properties, these fields can be found consistent with blazar observations based on the value of $B_{\rm eff}$. Finally, the intensity of the magnetic field for large coupling can be enough to satisfy the requirements for a recently proposed baryogenesis mechanism, which utilizes the chiral anomaly of the Standard Model.

We find multi-scalar effective field theories (EFTs) that can achieve a slow inflationary roll despite having a scalar potential that does not satisfy the usual slow-roll condition (d V)^2 << V^2/Mp^2. They evade the usual slow-roll conditions on $V$ because their kinetic energies are dominated by single-derivative terms rather than the usual two-derivative terms. Single derivatives dominate during slow roll and so do not require a breakdown of the usual derivative expansion that underpins calculational control in much of cosmology. The presence of such terms requires some sort of UV Lorentz-symmetry breaking during inflation (besides the usual cosmological breaking). Chromo-natural inflation provides an example of a UV theory that can generate the multi-field single-derivative terms we consider, and we argue that the EFT we find indeed captures the slow-roll conditions for the background evolution for Chromo-natural inflation. We also show that our EFT can be understood as a multi-field generalization of the single-field Cuscuton models. The multi-field case introduces a new feature, however: the scalar kinetic terms define a target-space 2-form, F_{ab}, whose antisymmetry gives new ways for slow roll to be achieved.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

We demonstrate that Chromo-Natural Inflation can be made consistent with observational data if the SU(2) gauge symmetry is spontaneously broken. Working in the Stueckelberg limit, we show that isocurvature is negligible, and the resulting adiabatic fluctuations can match current observational constraints. Observable levels of chirally-polarized gravitational radiation ($r\sim 10^{-3}$) can be produced while the evolution of all background fields is sub-Planckian. The gravitational wave spectrum is amplified via linear mixing with the gauge field fluctuations, and its amplitude is not simply set by the Hubble rate during inflation. This allows observable gravitational waves to be produced for an inflationary energy scale below the GUT scale. The tilt of the resulting gravitational wave spectrum can be either blue or red.

In a broad class of theories, the relic abundance of dark matter is determined by interactions internal to a thermalized dark sector, with no direct involvement of the Standard Model (SM). We point out that these theories raise an immediate cosmological question: how was the dark sector initially populated in the early universe? Motivated in part by the difficulty of accommodating large amounts of entropy carried in dark radiation with cosmic microwave background measurements of the effective number of relativistic species at recombination, $N_{\mathrm{eff}}$, we aim to establish which admissible cosmological histories can populate a thermal dark sector that never reaches thermal equilibrium with the SM. The minimal cosmological origin for such a dark sector is asymmetric reheating, when the same mechanism that populates the SM in the early universe also populates the dark sector at a lower temperature. Here we demonstrate that the resulting inevitable inflaton-mediated scattering between the dark sector and the SM can wash out a would-be temperature asymmetry, and establish the regions of parameter space where temperature asymmetries can be generated in minimal reheating scenarios. Thus obtaining a temperature asymmetry of a given size either restricts possible inflaton masses and couplings or necessitates a non-minimal cosmology for one or both sectors. As a side benefit, we develop techniques for evaluating collision terms in the relativistic Boltzmann equation when the full dependence on Bose-Einstein or Fermi-Dirac phase space distributions must be retained, and present several new results on relativistic thermal averages in an appendix.

We propose that the observed matter-antimatter asymmetry can be naturally produced as a byproduct of axion-driven slow-roll inflation by coupling the axion to standard-model neutrinos. We assume that GUT scale right-handed neutrinos are responsible for the masses of the standard model neutrinos and that the Higgs field is light during inflation and develops a Hubble scale vacuum expectation value (VEV). In this set up, the rolling axion generates a helicity asymmetry in standard-model neutrinos. Following inflation, this helicity asymmetry becomes equal to a net lepton number as the Higgs VEV decays and is partially re-processed by the $SU(2)_{L}$ sphaleron into a net baryon number.

The recent BICEP2 B-mode polarization determination of an inflationary tensor-scalar ratio $r=0.2^{+0.07}_{-0.05}$ is in tension with simple scale-free models of inflation due to a lack of a corresponding low multipole excess in the temperature power spectrum which places a limit of $r_{0.002}<0.11$ (95% CL) on such models. Single-field inflationary models that reconcile these two observations, even those where the tilt runs substantially, introduce a scale into the scalar power spectrum. To cancel the tensor excess, and simultaneously explain the excess already present in $\Lambda$CDM, ideally the model should introduce this scale as a relatively sharp transition in the tensor-scalar ratio around the horizon at recombination. We consider models which generate such a step in this quantity and find that they can improve the joint fit to the temperature and polarization data by up to $2\Delta \ln{\cal L} \approx -14$ without changing cosmological parameters. Precision E-mode polarization measurements should be able to test this explanation.

We examine the regime of validity of $N$-point spectra predictions of single field inflation models that invoke transient periods of non-adiabatic evolution. Such models generate oscillatory features in these spectra spanning frequencies up to the inverse time scale of the transient feature. To avoid strong coupling of fluctuations in these theories this scale must be at least $\sim 10^{-2}/c_s$ of the Hubble time during inflation, where $c_s$ is the inflaton sound speed. We show that, in such models, the signal-to-noise ratio of the bispectrum is bounded from above by that of the power spectrum, implying that searches for features due to non-adiabatic evolution are best focussed first on the latter.

Chromo-Natural Inflation is the first worked example of a model of inflation in which slow-roll inflation is achieved by "magnetic drift" as opposed to Hubble friction. In this work, we give an account of the perturbations at linear order in this model. Our analysis uncovers two novel phenomena. First, the amplitude of scalar curvature perturbations is not directly tied to the shape of the inflationary potential. This allows the theory to violate na\"ive formulations of the Lyth bound. Second, the tensor sector of the theory is significantly altered from the usual case: the non-Abelian gauge field perturbations have a tensor degree of freedom. One chirality of the this tensor can be exponentially enhanced by a temporary instability near horizon crossing; this chiral instability exists because of the classical gauge field background, which violates parity. These tensor fluctuations of the gauge field also couple to gravitational waves at linear order in perturbation theory and source a chiral spectrum of gravitational waves. This spectrum can be exponentially enhanced over the usual inflationary spectrum due to the instability in the gauge sector. These new features cause the theory in its present form to be in significant tension with current observational data. This is because the new scalar physics leads to a significant reddening of the spectral tilt in the same region of parameter space where the exponential enhancement of the gravitational wave amplitude is small enough to satisfy current constraints on the tensor-to-scalar index. Hence, the model either predicts a spectral tilt that is too red, or it overproduces gravitational waves, or both.

Models of inflation involving non-Abelian gauge field backgrounds can produce gravitational waves at an observable level with a preferred handedness. This asymmetry comes about because the non-Abelian background generates parity-violation in the action for perturbations. In the specific model we study, Chromo-Natural Inflation, these gravitational waves can be produced at observable levels even when no field makes a super-Planckian field excursion, thus evading a common formulation of the Lyth bound. Unfortunately, when considered in concert with the scalar fluctuations, this chiral enhancement of the gravitational waves makes the model observationally inviable.

The statistical properties of the primordial perturbations contain clues about the origins of those fluctuations. Although the Planck collaboration has recently obtained tight constraints on primordial non-gaussianity from cosmic microwave background measurements, it is still worthwhile to mine upcoming data sets in effort to place independent or competitive limits. The ionized bubbles that formed at redshift z~6-20 during the Epoch of Reionization are seeded by primordial overdensities, and so the statistics of the ionization field at high redshift are related to the statistics of the primordial field. Here we model the effect of primordial non-gaussianity on the reionization field. The epoch and duration of reionization are affected as are the sizes of the ionized bubbles, but these changes are degenerate with variations in the properties of the ionizing sources and the surrounding intergalactic medium. A more promising signature is the power spectrum of the spatial fluctuations in the ionization field, which may be probed by upcoming 21 cm surveys. This has the expected 1/k^2 dependence on large scales, characteristic of a biased tracer of the matter field. We project how well upcoming 21 cm observations will be able to disentangle this signal from foreground contamination. Although foreground cleaning inevitably removes the large-scale modes most impacted by primordial non-gaussianity, we find that primordial non-gaussianity can be separated from foreground contamination for a narrow range of length scales. In principle, futuristic redshifted 21 cm surveys may allow constraints competitive with Planck.

We study the impact of primordial non-Gaussianity generated during inflation on the bias of halos using excursion set theory. We recapture the familiar result that the bias scales as $k^{-2}$ on large scales for local type non-Gaussianity but explicitly identify the approximations that go into this conclusion and the corrections to it. We solve the more complicated problem of non-spherical halos, for which the collapse threshold is scale dependent.

We provide a detailed discussion of the multifield trajectories and inflationary dynamics of the recently proposed model of Chromo-Natural inflation, which allows for slow roll inflation on a steep potential with the aid of classical non-Abelian gauge fields. We show that slow roll inflation can be achieved across a wide range of the parameter space. We demonstrate that Chromo-Natural Inflation includes trajectories that match those found in Gauge-flation and describe how the theories are related.

We propose a new, generic mechanism of inflation mediated by a balance between potential forces and a Chern-Simons interaction. Such quasi-topological interactions are ubiquitous in string theory. In the minisuperspace approximation, their effect on the dynamics can be mapped onto the problem of a charged particle in an electromagnetic field together with an external potential; slow roll arises when the motion is dominated by the analogue of `magnetic drift'. This mechanism is robust against radiative corrections. We suggest a possible experimental signature which, if observed, might be considered strong evidence for string theory.

In Dirac-Born-Infeld inflation, changes in the sound speed that transiently break the slow roll approximation lead to features in the power spectrum. We develop and test the generalized slow roll approximation for calculating such effects and show that it can be extended to treat order unity features. As in slow-roll, model independent constraints on the potential of canonical inflation can be directly reinterpreted in the DBI context through this approximation. In particular, a sharp horizon scale step in the warped brane tension can explain oscillatory features in the WMAP7 CMB power spectrum as well as features in the potential. Differences appear only as a small suppression of power on horizon scales and larger.

We propose a model for inflation consisting of an axionic scalar field coupled to a set of three non-Abelian gauge fields. Our model's novel requirement is that the gauge fields begin inflation with a rotationally invariant vacuum expectation value (VEV) that is preserved through identification of SU(2) gauge invariance with rotations in three dimensions. The gauge VEV interacts with the background value of the axion, leading to an attractor solution that exhibits slow roll inflation even when the axion decay constant has a natural value ($<M_{\rm Pl}$). Assuming a sinusoidal potential for the axion, we find that inflation continues until the axionic potential vanishes. The speed at which the axion moves along its potential is modulated by its interactions with the gauge VEV, rather than being determined by the slope of its bare potential. For sub-Plankian axion decay constants vanishingly small tensor to scalar ratios are predicted, a direct consequence of the Lyth bound. The parameter that controls the interaction strength between the axion and the gauge fields requires a technically natural tuning of $\mathcal{O}$(100).

Features in the inflaton potential that are traversed in much less than an e-fold of the expansion can produce observably large non-Gaussianity. In these models first order corrections to the curvature mode function evolution induce effects at second order in the slow roll parameters that are generically greater than ~ 10% and can reach order unity for order unity power spectrum features. From a complete first order expression in generalized slow-roll, we devise a computationally efficient method that is as simple to evaluate as the leading order one and implements consistency relations in a controlled fashion. This expression matches direct numerical computation for step potential models of the dominant bispectrum configurations to better than 1% when features are small and 10% when features are order unity.

We provide analytic solutions for the power spectrum and bispectrum of curvature fluctuations produced by a step feature in the inflaton potential, valid in the limit that the step is short and sharp. In this limit, the bispectrum is strongly scale dependent and its effective non-linearity attains a large oscillatory amplitude. The perturbations to the curvature power spectrum, on the other hand, remain a small component on top of the usual spectrum of fluctuations generated by slow roll. We utilize our analytic solutions to assess the observability of the predicted non-Gaussian signatures and show that, if present, only very sharp steps on scales larger than ~ 2 Gpc are likely to be able to be detected by Planck. Such features are not only consistent with WMAP7 data, but can also improve its likelihood by 2 Delta ln L ~ 12 for two extra parameters, the step location and height. If this improvement were due to a slow roll violating step as considered here, a bispectrum or corresponding polarization power spectrum detection would provide definitive checks as to its primordial origin.

We develop a fast technique based on the generalized slow roll (GSR) approach for computing the curvature bispectrum of inflationary models with features. We show that all triangle configurations can be expressed in terms of three simple integrals over the inflationary background with typical accuracy of better than ~20%. With a first order GSR approach the typical accuracy can be improved to better than the 5% level. We illustrate this technique with the step potential model that has been invoked to explain the WMAP temperature power spectrum glitches at ell ~ 20-40 and show that the maximum likelihood model falls short of observability by more than a factor of 100 in amplitude. We also explicitly demonstrate that the bispectrum consistency relation with the local slope of the power spectrum is satisfied for these models. In the GSR approach, the bispectrum arises from integrals of nearly the same function of the background slow-roll parameters as the power spectrum but with a stronger weight to the epoch before horizon crossing. Hence this technique enables reverse engineering of models with large bispectrum but small power spectrum features.

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 consider the 3-pt function (i.e. the bispectrum or non-Gaussianity) for stochastic backgrounds of gravitational waves. We estimate the amplitude of this signal for the primordial inflationary background, gravitational waves generated during preheating, and for gravitational waves produced by self-ordering scalar fields following a global phase transition. To assess detectability, we describe how to extract the 3-pt signal from an idealized interferometric experiment and compute the signal to noise ratio as a function of integration time. The 3-pt signal for the stochastic gravitational wave background generated by inflation is unsurprisingly tiny. For gravitational radiation generated by purely causal, classical mechanisms we find that, no matter how non-linear the process is, the 3-pt correlations produced vanish in direct detection experiments. On the other hand, we show that in scenarios where the B-mode of the CMB is sourced by gravitational waves generated by a global phase transition, a strong 3-pt signal among the polarization modes could also be produced. This may provide another method of distinguishing inflationary B-modes. To carry out this computation, we have developed a diagrammatic approach to the calculation of stochastic gravitational waves sourced by scalar fluids, which has applications beyond the present scenario.

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.

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.

Of the many probes of reionization, the 21 cm line and the cosmic microwave background (CMB) are among the most effective. We examine how the cross-correlation of the 21 cm brightness and the CMB Doppler fluctuations on large angular scales can be used to study this epoch. We employ a new model of the growth of large scale fluctuations of the ionized fraction as reionization proceeds. We take into account the peculiar velocity field of baryons and show that its effect on the cross correlation can be interpreted as a mixing of Fourier modes. We find that the cross-correlation signal is strongly peaked toward the end of reionization and that the sign of the correlation should be positive because of the inhomogeneity inherent to reionization. The signal peaks at degree scales (l~100) and comes almost entirely from large physical scales (k~0.01 Mpc). Since many of the foregrounds and noise that plague low frequency radio observations will not correlate with CMB measurements, the cross correlation might appear to provide a robust diagnostic of the cosmological origin of the 21 cm radiation around the epoch of reionization. Unfortunately, we show that these signals are actually only weakly correlated and that cosmic variance dominates the error budget of any attempted detection. We conclude that the detection of a cross-correlation peak at degree-size angular scales is unlikely even with ideal experiments.