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 study the energy budget of a first-order cosmological phase transition, which is an important factor in the prediction of the resulting gravitational wave spectrum. Formerly, this analysis was based mostly on simplified models as for example the bag equation of state. Here, we present a model-independent approach that is exact up to the temperature dependence of the speed of sound in the broken phase. We find that the only relevant quantities that enter in the hydrodynamic analysis are the speed of sound in the broken phase and a linear combination of the energy and pressure differences between the two phases which we call pseudotrace (normalized to the enthalpy in the broken phase). The pseudotrace quantifies the strength of the phase transition and yields the conventional trace of the energy-momentum tensor for a relativistic plasma (with speed of sound squared of one third). We study this approach in several realistic models of the phase transition and also provide a code snippet that can be used to determine the efficiency coefficient for a given phase transition strength and speed of sound. It turns out that our approach is accurate to the percent level for moderately strong phase transitions, while former approaches give at best the right order of magnitude.

We investigate the role of leptons in electroweak baryogenesis by studying a relatively simple framework inspired by effective field theory that satisfies all Sakharov conditions. In particular, we study the effectiveness of CP-violating source terms induced by dimension-six Yukawa interactions for quarks and charged leptons. Despite the relatively small Yukawa coupling, CP-violating source terms involving taus are quite effective and can account for the observed matter-antimatter asymmetry. We obtain analytical and numerical expressions for the total baryon asymmetry, the former providing important insight into what makes lepton CP violation relatively effective compared to quark CP violation. Leptons also play an important role if the CP-violating source involves top quarks. While the tau Yukawa coupling in the Standard Model is small, it significantly enhances the baryon asymmetry by transferring the chiral asymmetry in quarks, which is washed out by strong sphalerons, to a chiral asymmetry in leptons. We conclude that leptons should not be ignored even if CP violation is limited to the quark sector. The role of leptons can be further increased in scenarios of new physics with additional chiral-symmetry-breaking interactions between quarks and leptons, as can happen in models with additional Higgs bosons or leptoquarks. Finally, we study CP-violating dimension-six Yukawa interactions for lighter quarks and leptons but conclude that these lead to too small baryon asymmetries.

We point out that a model with scalar fields with a large nonminimal coupling to the Ricci scalar, such as Higgs inflation, can be regarded as a nonlinear sigma model (NLSM). The $\sigma$-meson, which is induced by quantum corrections, is identified as the scalaron in this model. Our understanding provides a novel alternative picture for the emergence of the scalaron, which was previously studied in the Jordan frame by computing the running of the term quadratic in the Ricci scalar. We demonstrate that quantum corrections inevitably induce other operators on top of this NLSM, which give rise to the $\sigma$-meson. We confirm that the $\sigma$-meson indeed corresponds to the scalaron. With the help of the $\sigma$-meson, the NLSM is UV-completed to be a linear sigma model (LSM). We show that the LSM only involves renormalizable interactions and hence its perturbative unitarity holds up to the Planck scale unless it hits a Landau pole, which is in agreement with the renormalizability of quadratic gravity.

In electroweak baryogenesis the baryon asymmetry of the universe is created during the electroweak phase transition. The quantum transport equations governing the dynamics of the plasma particles can be derived in the vev-insertion approximation, which treats the vev-dependent part of the particle masses as a perturbation. We calculate the next-to-leading order (NLO) contribution to the CP-violating source term and CP-conserving relaxation rate, corresponding to Feynman diagrams for the self-energies with four mass insertions. We consider both a pair of Weyl fermions and a pair of complex scalars, that scatter off the bubble wall. We find: (i) The NLO correction becomes large for $\mathcal O(1)$ couplings. If only the Standard Model (SM) Higgs obtains a vev during the phase transition, this implies the vev-insertion approximation breaks down for top quarks. (ii) The resonant enhancement of the source term and relaxation rate, that exists at leading order in the limit of degenerate thermal masses for the fermions/scalars, persists at NLO.

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.

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 investigate electroweak baryogenesis within the framework of the Standard Model Effective Field Theory. The Standard Model Lagrangian is supplemented by dimension-six operators that facilitate a strong first-order electroweak phase transition and provide sufficient CP violation. Two explicit scenarios are studied that are related via the classical equations of motion and are therefore identical at leading order in the effective field theory expansion. We demonstrate that formally higher-order dimension-eight corrections lead to large modifications of the matter-antimatter asymmetry. The effective field theory expansion breaks down in the modified Higgs sector due to the requirement of a first-order phase transition. We investigate the source of the breakdown in detail and show how it is transferred to the CP-violating sector. We briefly discuss possible modifications of the effective field theory framework.

The measured values of the Higgs and top quark mass indicate that the electroweak vacuum is metastable if there is no new physics below the Planck scale. This is at odds with a period of high scale inflation. A non-minimal coupling between the Higgs field and the Ricci scalar can stabilize the vacuum as it generates a large effective Higgs mass during inflation. We consider the effect of this coupling during preheating, when Higgs modes can be produced very efficiently due to the oscillating Ricci scalar. We compute their effect on the effective potential and the energy density. The Higgs excitations are defined with respect to the adiabatic vacuum. We study the adiabaticity conditions and find that the dependence of our results on the choice of the order of the adiabatic vacuum increases with time. For large enough coupling particle production is so efficient that the Higgs decays to the true vacuum before this is an issue. However, for smaller values of the Higgs-curvature coupling no definite statements can be made as the vacuum dependence is large.