The quartic and trilinear Higgs field couplings to an additional real scalar are renormalizable, gauge and Lorentz invariant. Thus, on general grounds, one expects such couplings between the Higgs and an inflaton in quantum field theory. In particular, the (often omitted) trilinear coupling is motivated by the need for reheating the Universe after inflation, whereby the inflaton decays into the Standard Model (SM) particles. Such a coupling necessarily leads to the Higgs-inflaton mixing, which could stabilize the electroweak vacuum by increasing the Higgs self-coupling. We find that the inflationary constraints on the trilinear coupling are weak such that the Higgs-inflaton mixing up to order one is allowed, making it accessible to colliders. This entails an exciting possibility of a direct inflaton search at the LHC.

The Higgs-inflaton coupling plays an important role in the Higgs field dynamics in the early Universe. Even a tiny coupling generated at loop level can have a dramatic effect on the fate of the electroweak vacuum. Such Higgs-inflaton interaction is present both at the trilinear and quartic levels in realistic reheating models. In this work, we examine the Higgs dynamics during the preheating epoch, focusing on the effects of the parametric and tachyonic resonances. We use lattice simulations and other numerical tools in our studies. We find that the resonances can induce large fluctuations of the Higgs field which destabilize the electroweak vacuum. Our considerations thus provide an upper bound on quartic and trilinear interactions between the Higgs and the inflaton. We conclude that there exists a favourable range of the couplings within which the Higgs field is stabilized during both inflation and preheating epochs.

We study the resonant decay of the primordial Standard Model Higgs condensate after inflation into $SU(2)$ gauge bosons on the lattice. We find that the non-Abelian interactions between the gauge bosons quickly extend the momentum distribution towards high values, efficiently destroying the condensate after the onset of backreaction. For the inflationary scale $H = 10^8$ GeV, we find that 90% of the Higgs condensate has decayed after $n \sim 10$ oscillation cycles. This differs significantly from the Abelian case where, given the same coupling strengths, most of the condensate would persist after the resonance.

We study the resonant decay of the Higgs condensate into weak gauge bosons after inflation and estimate the corrections arising from the non-Abelian self-interactions of the gauge fields. We find that non-Abelian interaction terms induce an effective mass which tends to shut down the resonance. For the broad resonance relevant for the Standard Model Higgs the produced gauge particles backreact on the dynamics of the Higgs condensate before the non-Abelian terms grow large. The non-Abelian terms can however significantly affect the final stages of the resonance after the backreaction. In the narrow resonance regime, which may be important for extensions of the Standard Model, the non-Abelian terms affect already the linear stage and terminate the resonance before the Higgs condensate is affected by the backreaction of decay products.

The behaviour of oscillating scalar spectator fields after inflation depends on the thermal background produced by inflaton decay. Resonant decay of the spectator is often blocked by large induced thermal masses. We account for the finite decay width of the inflaton and the protracted build-up of the thermal bath to determine the early evolution of a homogeneous spectator field, \sigma, coupled to the Higgs Boson, \Phi, through the term g^2 \sigma^2 \Phi^2, the only renormalisable coupling of a new scalar to the Standard Model. We find that for very large higgs-spectator coupling g > 10^{-3}, the resonance is not always blocked as was previously suggested. As a consequence, the oscillating spectator can decay quickly. For other parameter values, we find that although qualitative features of the thermal blocking still hold, the dynamics are altered compared to the instant decay case. These findings are important for curvaton models, where the oscillating field must be relatively long lived in order to produce the curvature perturbation. They are also relevant for other spectator fields, which must decay sufficiently early to avoid spoiling the predictions of baryogenesis and nucleosynthesis.

We use WMAP 9-year and other CMB data to constrain cosmological models where the primordial perturbations have both an adiabatic and a (possibly correlated) neutrino density (NDI), neutrino velocity (NVI), or cold dark matter density (CDI) isocurvature component. For NDI and CDI we use both a phenomenological approach, where primordial perturbations are parametrized in terms of amplitudes at two scales, and a slow-roll two-field inflation approach, where slow-roll parameters are used as primary parameters. For NVI we use only the phenomenological approach, since it is difficult to imagine a connection with inflation. We find that in the NDI and NVI cases larger isocurvature fractions are allowed than in the corresponding models with CDI. For uncorrelated perturbations, the upper limit to the primordial NDI (NVI) fraction is 24% (20%) at k = 0.002 Mpc^{-1} and 28% (16%) at k = 0.01 Mpc^{-1}. For maximally correlated (anticorrelated) perturbations, the upper limit to the NDI fraction is 3.0% (0.9%). The nonadiabatic contribution to the CMB temperature variance can be as large as 10% (-13%) for the NDI (NVI) modes. Bayesian model comparison favors pure adiabatic initial mode over the mixed primordial adiabatic and NDI, NVI, or CDI perturbations. At best, the betting odds for a mixed model (uncorrelated NDI) are 1:3.4 compared to the pure adiabatic model. For the phenomenological generally correlated mixed models the odds are about 1:100, whereas the slow-roll approach leads to 1:13 (NDI) and 1:51 (CDI).

We consider a model of preheating where the coupling of the inflaton to the preheat field is modulated by an additional scalar field which is light during inflation. We establish that such a model produces the observed curvature perturbation analogously to the modulated reheating scenario. The contribution of modulated preheating to the power spectrum and to non-Gaussianity can however be significantly larger compared to modulated perturbative reheating. We also consider the implications of the current constraints on isocurvature perturbations in case where the modulating field is responsible for cold dark matter. We find that existing bounds on CDM isocurvature perturbations imply that modulated preheating is unlikely to give a dominant contribution to the curvature perturbation and that the same bounds suggest important constraints on non-Gaussianity and the amount of primordial gravitational waves.

We constrain cosmological models where the primordial perturbations have both an adiabatic and a (possibly correlated) cold dark matter (CDM) or baryon isocurvature component. We use both a phenomenological approach, where the primordial power spectra are parametrized with amplitudes and spectral indices, and a slow-roll two-field inflation approach where slow-roll parameters are used as primary parameters. In the phenomenological case, with CMB data, the upper limit to the CDM isocurvature fraction is \alpha<6.4% at k=0.002Mpc^{-1} and 15.4% at k=0.01Mpc^{-1}. The median 95% range for the non-adiabatic contribution to the CMB temperature variance is -0.030<\alpha_T<0.049. Including the supernova (or large-scale structure, LSS) data, these limits become: \alpha<7.0%, 13.7%, and -0.048<\alpha_T< 0.042 (or \alpha<10.2%, 16.0%, and -0.071<\alpha_T<0.024). The CMB constraint on the tensor-to-scalar ratio, r<0.26 at k=0.01Mpc^{-1}, is not affected by the nonadiabatic modes. In the slow-roll two-field inflation approach, the spectral indices are constrained close to 1. This leads to tighter limits on the isocurvature fraction, with the CMB data \alpha<2.6% at k=0.01Mpc^{-1}, but the constraint on \alpha_T is not much affected, -0.058<\alpha_T<0.045. Including SN (or LSS) data, these limits become: \alpha< 3.2% and -0.056<\alpha_T<0.030 (or \alpha<3.4% and -0.063<\alpha_T<-0.008). When all spectral indices are close to each other the isocurvature fraction is somewhat degenerate with the tensor-to-scalar ratio. In addition to the generally correlated models, we study also special cases where the perturbation modes are uncorrelated or fully (anti)correlated. We calculate Bayesian evidences (model probabilities) in 21 different cases for our nonadiabatic models and for the corresponding adiabatic models, and find that in all cases the data support the pure adiabatic model.