We have conducted extensive lattice simulations to study the post-inflation dynamics of multifield models involving nonminimal couplings. We explore the parameter dependence of preheating in these models and describe the various time-scales that control such nonlinear processes as energy transfer, re-scattering, and the approach to radiation-domination and thermalization. In the limit of large nonminimal couplings ($\xi_I \sim 100$), we find that efficient transfer of energy from the inflaton condensate to radiative degrees of freedom, emergence of a radiation-dominated equation of state, and the onset of thermalization each consistently occur within $N_{\rm reh} \lesssim 3$ $e$-folds after the end of inflation, largely independent of the values of the other couplings in the models. The exception is the case of negative ellipticity, in which there is a misalignment between the dominant direction in field-space along which the system evolves and the larger of the nonminimal couplings $\xi_I$. In those cases, the field-space-driven parametric resonance is effectively shut off. More generally, the competition between the scalar fields' potential and the field-space manifold structure can yield interesting phenomena such as two-stage resonances. Despite the explosive particle production, which can lead to a quick depletion of the background energy density, the nonlinear processes do not induce any super-horizon correlations after the end of inflation in these models, which keeps predictions for CMB observables unaffected by the late-time amplification of isocurvature fluctuations. Hence the excellent agreement between primordial observables and recent observations is preserved for this class of models, even when we consider post-inflation dynamics.

We analyze and compare the multi-field dynamics during inflation and preheating in symmetric and asymmetric models of $\alpha$-attractors, characterized by a hyperbolic field-space manifold. We show that the generalized (asymmetric) E- and (symmetric) T-models exhibit identical two-field dynamics during inflation for a wide range of initial conditions. The resulting motion can be decomposed in two approximately single-field segments connected by a sharp turn in field-space. The details of preheating can nevertheless be different. For the T-model one main mass-scale dominates the evolution of fluctuations of the spectator field, whereas for the E-model, a competing mass-scale emerges due to the steepness of the potential away from the inflationary plateau, leading to different contributions to parametric resonance for small and large wave-numbers. Our linear multi-field analysis of fluctuations indicates that for highly curved manifolds, both the E- and T-models preheat almost instantaneously. For massless fields this is always due to efficient tachyonic amplification of the spectator field, making single-field results inaccurate. Interestingly, there is a parameter window corresponding to $r={\cal O}(10^{-5})$ and massive fields, where the preheating behavior is qualitatively and quantitatively different for symmetric and asymmetric potentials. In that case, the E-model can completely preheat due to self-resonance for values of the curvature where preheating in the T-model is inefficient. This provides a first distinguishing feature between models that otherwise behave identically, both at the single-field and multi-field level. Finally, we discuss how one can describe multi-field preheating on a hyperbolic manifold by identifying the relevant mass-scales that control the growth of inflaton and spectator fluctuations, which can be applied to any $\alpha$-attractor model and beyond.

We study the post-inflation dynamics of multifield models involving nonminimal couplings using lattice simulations to capture significant nonlinear effects like backreaction and rescattering. We measure the effective equation of state and typical time-scales for the onset of thermalization, which could affect the usual mapping between predictions for primordial perturbation spectra and measurements of anisotropies in the cosmic microwave background radiation. For large values of the nonminimal coupling constants, we find efficient particle production that gives rise to nearly instantaneous preheating. Moreover, the strong single-field attractor behavior that was previously identified persists until the end of preheating, thereby suppressing typical signatures of multifield models. We therefore find that predictions for primordial observables in this class of models retain a close match to the latest observations.

The Standard Model Higgs boson, which has previously been shown to develop an effective vacuum expectation value during inflation, can give rise to large particle masses during inflation and reheating, leading to temporary blocking of the reheating process and a lower reheat temperature after inflation. We study the effects on the multiple stages of reheating: resonant particle production (preheating) as well as perturbative decays from coherent oscillations of the inflaton field. Specifically, we study both the cases of the inflaton coupling to Standard Model fermions through Yukawa interactions as well as to Abelian gauge fields through a Chern-Simons term. We find that, in the case of perturbative inflaton decay to SM fermions, reheating can be delayed due to Higgs blocking and the reheat temperature can decrease by up to an order of magnitude. In the case of gauge-reheating, Higgs-generated masses of the gauge fields can suppress preheating even for large inflaton-gauge couplings. In extreme cases, preheating can be shut down completely and must be substituted by perturbative decay as the dominant reheating channel. Finally, we discuss the distribution of reheat temperatures in different Hubble patches, arising from the stochastic nature of the Higgs VEV during inflation and its implications for the generation of both adiabatic and isocurvature fluctuations.

We explore preheating in multi-field models of inflation in which the field-space metric is a highly curved hyperbolic manifold. One broad family of such models is called $\alpha$-attractors, whose single-field regimes have been extensively studied in the context of inflation and supergravity. We focus on a simple two-field generalization of the $T$-model, which has received renewed attention in the literature. Krajewski et al. concluded, using lattice simulations, that multi-field effects can dramatically speed-up preheating. We recover their results and further demonstrate that significant analytical progress can be made for preheating in these models using the WKB approximation and Floquet analysis. We find a simple scaling behavior of the Floquet exponents for large values of the field-space curvature, that enables a quick estimation of the $T$-model reheating efficiency for any large value of the field-space curvature. In this regime we further observe and explain universal preheating features that arise for different values of the potential steepness. In general preheating is faster for larger negative values of the field-space curvature and steeper potentials. For very highly curved field-space manifolds preheating is essentially instantaneous.

Recent years have seen the introduction of various multi-field inflationary scenarios, in which the curvature and geodesics of the scalar manifold play a crucial role. We outline a simple description that unifies these different proposals and discuss their stability criteria. We demonstrate how the underlying dynamics is governed by an effective potential, whose critical points and bifurcations determine the late-time behaviour of the system, thus unifying hyperinflation, angular, orbital and side-tracked inflation. Interestingly, hyperinflation is shown to be a special case of side-tracked inflation. This equivalence relies on the enhanced isometries of the hyperbolic manifold and we provide the explicit coordinate transformation that maps the two models into each other.

$SU(2)$ gauge fields can generate large gravitational waves during inflation, if they are coupled to an axion which can be either the inflaton or a spectator field. The shape of the produced tensor power spectrum $\mathcal{P}_h$ depends on the form of the axion potential. We derive analytic expressions and provide general templates for $\mathcal{P}_h$ for various types of the spectator axion potential. Furthermore, we explore the detectability of the oscillatory feature, which is present in $\mathcal{P}_h$ in the case of an axion monodromy model, by possible future CMB B-mode polarization observations.

We perform an extensive analysis of linear fluctuations during preheating in Higgs inflation in the Einstein frame, where the fields are minimally coupled to gravity, but the field-space metric is nontrivial. The self-resonance of the Higgs and the Higgsed gauge bosons are governed by effective masses that scale differently with the nonminimal couplings and evolve differently in time. Coupled metric perturbations enhance Higgs self-resonance and make it possible for Higgs inflation to preheat solely through this channel. For $\xi\gtrsim 100$ the total energy of the Higgs-inflaton condensate can be transferred to Higgs particles within $3$ $e$-folds after the end of inflation. For smaller values of the nonminimal coupling preheating takes longer, completely shutting off at around $\xi\simeq 30$. The production of gauge bosons is dominated by the gauge boson mass and the field space curvature. For large values of the nonminimal coupling $\xi \gtrsim 1000$, it is possible for the Higgs condensate to transfer the entirety of its energy into gauge fields within one oscillation. For smaller values of the nonminimal coupling gauge bosons decay very quickly into fermions, thereby shutting off Bose enhancement. Estimates of non-Abelian interactions indicate that they will not suppress preheating into gauge bosons for $\xi \gtrsim 1000$.

We explore the dynamics of multi-field models of inflation in which the field-space metric is a hyperbolic manifold of constant curvature. Such models are known as $\alpha$-attractors and their single-field regimes have been extensively studied in the context of inflation and supergravity. We find a variety of multi-field inflationary trajectories in different regions of parameter space, which is spanned by the mass parameters and the hyperbolic curvature. Amongst these is a novel dynamical attractor along the boundary of the Poincare disc which we dub "angular inflation". We calculate the evolution of adiabatic and isocurvature fluctuations during this regime and show that, while isocurvature modes decay during this phase, the duration of the angular inflation period can shift the single-field predictions of $\alpha$-attractors.

This is the second in a series of papers on preheating in inflationary models comprised of multiple scalar fields coupled nonminimally to gravity. In this paper, we work in the rigid-spacetime approximation and consider field trajectories within the single-field attractor, which is a generic feature of these models. We construct the Floquet charts to find regions of parameter space in which particle production is efficient for both the adiabatic and isocurvature modes, and analyze the resonance structure using analytic and semi-analytic techniques. Particle production in the adiabatic direction is characterized by the existence of an asymptotic scaling solution at large values of the nonminimal couplings, $\xi_I \gg 1$, in which the dominant instability band arises in the long-wavelength limit, for comoving wavenumbers $k \rightarrow 0$. However, the large-$\xi_I$ regime is not reached until $\xi_I \geq {\cal O} (100)$. In the intermediate regime, with $\xi_I \sim {\cal O}(1 - 10)$, the resonance structure depends strongly on wavenumber and couplings. The resonance structure for isocurvature perturbations is distinct and more complicated than its adiabatic counterpart. An intermediate regime, for $\xi_I \sim {\cal O} (1 - 10)$, is again evident. For large values of $\xi_I$, the Floquet chart consists of densely spaced, nearly parallel instability bands, suggesting a very efficient preheating behavior. The increased efficiency arises from features of the nontrivial field-space manifold in the Einstein frame, which itself arises from the fields' nonminimal couplings in the Jordan frame, and has no analogue in models with minimal couplings. Quantitatively, the approach to the large-$\xi_I$ asymptotic solution for isocurvature modes is slower than in the case of the adiabatic modes.

This paper concludes our semi-analytic study of preheating in inflationary models comprised of multiple scalar fields coupled nonminimally to gravity. Using the covariant framework of paper I in this series, we extend the rigid-spacetime results of paper II by considering both the expansion of the universe during preheating, as well as the effect of the coupled metric perturbations on particle production. The adiabatic and isocurvature perturbations are governed by different effective masses that scale differently with the nonminimal couplings and evolve differently in time. The effective mass for the adiabatic modes is dominated by contributions from the coupled metric perturbations immediately after inflation. The metric perturbations contribute an oscillating tachyonic term that enhances an early period of significant particle production for the adiabatic modes, which ceases on a time-scale governed by the nonminimal couplings $\xi_I$. The effective mass of the isocurvature perturbations, on the other hand, is dominated by contributions from the fields' potential and from the curvature of the field-space manifold (in the Einstein frame), the balance between which shifts on a time-scale governed by $\xi_I$. As in papers I and II, we identify distinct behavior depending on whether the nonminimal couplings are small ($\xi_I \lesssim {\cal O} (1)$), intermediate ($\xi_I \sim {\cal O} (1 - 10)$), or large ($\xi_I \geq 100$).

This is the first of a three-part series of papers, in which we study the preheating phase for multifield models of inflation involving nonminimal couplings. In this paper, we study the single-field attractor behavior that these models exhibit during inflation and quantify its strength and parameter dependence. We further demonstrate that the strong single-field attractor behavior persists after the end of inflation. Preheating in such models therefore generically avoids the "de-phasing" that typically affects multifield models with minimally coupled fields, allowing efficient transfer of energy from the oscillating inflaton condensate(s) to coupled perturbations across large portions of parameter space. We develop a doubly-covariant formalism for studying the preheating phase in such models and identify several features specific to multifield models with nonminimal couplings, including effects that arise from the nontrivial field-space manifold. In papers II and III, we apply this formalism to study how the amplification of adiabatic and isocurvature perturbations varies with parameters, highlighting several distinct regimes depending on the magnitude of the nonminimal couplings $\xi_I$.

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 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.

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 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.

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.

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.

Hybrid inflation models are especially interesting as they lead to a spike in the density power spectrum on small scales, compared to the CMB, while also satisfying current bounds on tensor modes. Here we study hybrid inflation with $N$ waterfall fields sharing a global $SO(N)$ symmetry. The inclusion of many waterfall fields has the obvious advantage of avoiding topologically stable defects for $N>3$. We find that it also has another advantage: it is easier to engineer models that can simultaneously (i) be compatible with constraints on the primordial spectral index, which tends to otherwise disfavor hybrid models, and (ii) produce a spike on astrophysically large length scales. The latter may have significant consequences, possibly seeding the formation of astrophysically large black holes. We calculate correlation functions of the time-delay, a measure of density perturbations, produced by the waterfall fields, as a convergent power series in both $1/N$ and the field's correlation function $\Delta(x)$. We show that for large $N$, the two-point function is $<\delta t({\bf x})\,\delta t({\bf 0})>\,\propto\Delta^2(|{\bf x}|)/N$ and the three-point function is $<\delta t({\bf x})\,\delta t({\bf y})\,\delta t({\bf 0})>\,\propto\Delta(|{\bf x}-{\bf y}|)\Delta(|{\bf x}|)\Delta(|{\bf y}|)/N^2$. In accordance with the central limit theorem, the density perturbations on the scale of the spike are Gaussian for large $N$ and non-Gaussian for small $N$.

Recent measurements by the {\it Planck} experiment of the power spectrum of temperature anisotropies in the cosmic microwave background radiation (CMB) reveal a deficit of power in low multipoles compared to the predictions from best-fit $\Lambda$CDM cosmology. The low-$\ell$ anomaly may be explained by the presence of primordial isocurvature perturbations in addition to the usual adiabatic spectrum, and hence may provide the first robust evidence that early-universe inflation involved more than one scalar field. In this paper we explore the production of isocurvature perturbations in nonminimally coupled two-field inflation. We find that this class of models readily produces enough power in the isocurvature modes to account for the \emph{Planck} low-$\ell$ anomaly, while also providing excellent agreement with the other {\it Planck} results.

Multifield models of inflation with nonminimal couplings are in excellent agreement with the recent results from {\it Planck}. Across a broad range of couplings and initial conditions, such models evolve along an effectively single-field attractor solution and predict values of the primordial spectral index and its running, the tensor-to-scalar ratio, and non-Gaussianities squarely in the observationally most-favored region. Such models also can amplify isocurvature perturbations, which could account for the low power recently observed in the CMB power spectrum at low multipoles. Future measurements of primordial isocurvature perturbations and the tensor-to-scalar ratio may serve to distinguish between the currently viable possibilities.

Realistic models of high-energy physics include multiple scalar fields. Renormalization requires that the fields have nonminimal couplings to the spacetime Ricci curvature scalar, and the couplings can be large at the energy scales of early-universe inflation. The nonminimal couplings induce a nontrivial field-space manifold in the Einstein frame, and they also yield an effective potential in the Einstein frame with nontrivial curvature. The ridges or bumps in the Einstein-frame potential can lead to primordial non-Gaussianities of observable magnitude. We develop a covariant formalism to study perturbations in such models and calculate the primordial bispectrum. As in previous studies of non-Gaussianities in multifield models, our results for the bispectrum depend sensitively on the fields' initial conditions.

Higgs inflation is a simple and elegant model in which early-universe inflation is driven by the Higgs sector of the Standard Model. The Higgs sector can support early-universe inflation if it has a large nonminimal coupling to the Ricci spacetime curvature scalar. At energies relevant to such an inflationary epoch, the Goldstone modes of the Higgs sector remain in the spectrum in renormalizable gauges, and hence their effects should be included in the model's dynamics. We analyze the multifield dynamics of Higgs inflation and find that the multifield effects damp out rapidly after the onset of inflation, because of the gauge symmetry among the scalar fields in this model. Predictions from Higgs inflation for observable quantities, such as the spectral index of the power spectrum of primordial perturbations, therefore revert to their familiar single-field form, in excellent agreement with recent measurements. The methods we develop here may be applied to any multifield model with nonminimal couplings in which the ${\cal N}$ fields obey an $SO ({\cal N})$ symmetry in field space.

We introduce a new method for calculating density perturbations in hybrid inflation which avoids treating the fluctuations of the "waterfall" field as if they were small perturbations about a classical trajectory. We quantize only the waterfall field, treating it as a free quantum field with a time-dependent $m^2$, which evolves from positive values to tachyonic values. Although this potential has no minimum, we think it captures the important dynamics that occurs as $m^2$ goes through zero, at which time a large spike in the density perturbations is generated. We assume that the time-delay formalism provides an accurate approximation to the density perturbations, and proceed to calculate the power spectrum of the time delay fluctuations. While the evolution of the field is linear, the time delay is a nonlinear function to which all modes contribute. Using the Gaussian probability distribution of the mode amplitudes, we express the time-delay power spectrum as an integral which can be carried out numerically. We use this method to calculate numerically the spectrum of density perturbations created in hybrid inflation models for a wide range of parameters. A characteristic of the spectrum is the appearance of a spike at small length scales, which can be used to relate the model parameters to observational data. It is conceivable that this spike could seed the formation of black holes that can evolve to become the supermassive black holes found at the centers of galaxies.