We present a comprehensive discussion on lattice techniques for the simulation of scalar and gauge field dynamics in an expanding universe. After reviewing the continuum formulation of scalar and gauge field interactions in Minkowski and FLRW backgrounds, we introduce basic tools for the discretization of field theories, including lattice gauge invariant techniques. Following, we discuss and classify numerical algorithms, ranging from methods of $O(dt^2)$ accuracy like $staggered~leapfrog$ and $Verlet$ integration, to $Runge-Kutta$ methods up to $O(dt^4)$ accuracy, and the $Yoshida$ and $Gauss-Legendre$ higher-order integrators, accurate up to $O(dt^{10})$. We adapt these methods for their use in classical lattice simulations of the non-linear dynamics of scalar and gauge fields in an expanding grid in $3+1$ dimensions, including the case of `self-consistent' expansion sourced by the volume average of the fields' energy and pressure densities. We present lattice formulations of canonical cases of: $i)$ Interacting scalar fields, $ii)$ Abelian $U(1)$ gauge theories, and $iii)$ Non-Abelian $SU(2)$ gauge theories. In all three cases we provide symplectic integrators, with accuracy ranging from $O(dt^2)$ up to $O(dt^{10})$. For each algorithm we provide the form of relevant observables, such as energy density components, field spectra and the Hubble constraint. We note that all our algorithms for gauge theories always respect the Gauss constraint to machine precision, even in the case of 'self-consistent' expansion. As a numerical example we analyze the post-inflationary dynamics of an oscillating inflaton charged under $SU(2)\times U(1)$. We note that the present manuscript is meant as part of the theoretical basis for the code $CosmoLattice$, a multi-purpose MPI-based package for simulating the non-linear evolution of field theories in an expanding universe, to be released in 2020.

We study the energy distribution and equation of state of the universe between the end of inflation and the onset of radiation domination (RD), considering observationally consistent single-field inflationary scenarios, with a potential 'flattening' at large field values, and a monomial shape $V(\phi) \propto |\phi|^p$ around the origin. We include a quadratic interaction $g^2\phi^2X^2$ between the inflaton $\phi$ and a light scalar 'daughter' field $X$, as a proxy for (p)reheating. We capture the non-perturbative and non-linear nature of the system dynamics with lattice simulations, obtaining that: $i)$ the final energy transferred to $X$ depends only on $p$, not on $g^2$; $ii)$ the final transfer of energy is always negligible for $2 \leq p < 4$, and of order $\sim 50\%$ for $p \geq 4$; $iii)$ the system goes at late times to matter-domination for $p = 2$, and always to RD for $p > 2$. In the latter case we calculate the number of e-folds until RD, significantly reducing the uncertainty in the inflationary observables $n_s$ and $r$.

Oscillons are spatially localised strong fluctuations of a scalar field. They can e.g. form after inflation when the scalar field potential is shallower than quadratic away from the minimum. Although oscillons are not protected by topology, they can be remarkably stable and have a significant impact on the (p)reheating phase. In this work we investigate the properties of oscillons in hilltop-shaped potentials, in particular the typical energies, shapes and lifetimes. In the first part of the paper, we simulate oscillon creation and stabilization with (3+1)-dimensional classical lattice simulations, and extract the typical energies, radii and amplitudes of the oscillons. In the second part we approximate the oscillons as spherically symmetric, and simulate single oscillons until their decay. We find that typical oscillons live up to about 4-5 e-folds, with the individual lifetime of the oscillons depending mainly on the initial shape of the oscillon and the power-law coefficient characterising the particular hilltop model. We also observe a breathing mode in the oscillon radii and amplitudes, and find that stronger breathing implies shorter lifetimes.

We study the post-inflationary dynamics of the Standard Model (SM) Higgs field in the presence of a non-minimal coupling $\xi|\Phi|^2R$ to gravity, both with and without the electroweak gauge fields coupled to the Higgs. We assume a minimal scenario in which inflation and reheating are caused by chaotic inflation with a quadratic potential, and no additional new physics is relevant below the Planck scale. By using classical real-time lattice simulations with a renormalisation group improved effective Higgs potential and by demanding the stability of the Higgs vacuum after inflation, we obtain upper bounds for $\xi$, taking into account the experimental uncertainty of the top-Yukawa coupling. We compare the bounds in the absence and presence of the electroweak gauge bosons, and conclude that the addition of gauge interactions has a rather minimal impact. In the unstable cases, we parametrize the time when such instability develops. For a top-quark mass $m_t \approx173.3 {\rm GeV}$, the Higgs vacuum instability is triggered for $\xi \gtrsim 4 -5$, although a slightly lower mass of $m_t \approx 172.1 {\rm GeV}$ pushes up this limit to $\xi \gtrsim 11 - 12$. This, together with the estimation $\xi \gtrsim 0.06$ for stability during inflation, provides tight constraints to the Higgs-curvature coupling within the SM.

Parametric resonance is among the most efficient phenomena generating gravitational waves (GWs) in the early Universe. The dynamics of parametric resonance, and hence of the GWs, depend exclusively on the resonance parameter $q$. The latter is determined by the properties of each scenario: the initial amplitude and potential curvature of the oscillating field, and its coupling to other species. Previous works have only studied the GW production for fixed value(s) of $q$. We present an analytical derivation of the GW amplitude dependence on $q$, valid for any scenario, which we confront against numerical results. By running lattice simulations in an expanding grid, we study for a wide range of $q$ values, the production of GWs in post-inflationary preheating scenarios driven by parametric resonance. We present simple fits for the final amplitude and position of the local maxima in the GW spectrum. Our parametrization allows to predict the location and amplitude of the GW background today, for an arbitrary $q$. The GW signal can be rather large, as $h^2\Omega_{\rm GW}(f_p) \lesssim 10^{-11}$, but it is always peaked at high frequencies $f_p \gtrsim 10^{7}$ Hz. We also discuss the case of spectator-field scenarios, where the oscillatory field can be e.g.~a curvaton, or the Standard Model Higgs.

Particle production via parametric resonance in the early Universe, is a nonperturbative, non-linear and out-of-equilibrium phenomenon. Although it is a well studied topic, whenever a new scenario exhibits parametric resonance, a full re-analysis is normally required. To avoid this tedious task, many works present often only a simplified linear treatment of the problem. In order to surpass this circumstance in the future, we provide a fitting analysis of parametric resonance through all its relevant stages: initial linear growth, non-linear evolution, and relaxation towards equilibrium. Using lattice simulations in an expanding grid in 3 + 1 dimensions, we parametrize the dynamics outcome scanning over the relevant ingredients: role of the oscillatory field, particle coupling strength, initial conditions, and background expansion rate. We emphasize the inaccuracy of the linear calculation of the decay time of the oscillatory field, and propose a more appropriate definition of this scale based on the subsequent non-linear dynamics. We provide simple fits to the relevant time scales and particle energy fractions at each stage. Our fits can be applied to post-inflationary preheating scenarios, where the oscillatory field is the inflaton, or to spectator-field scenarios, where the oscillatory field can be e.g. a curvaton, or the Standard Model Higgs.

We study the nonperturbative dynamics of the Standard Model (SM) after inflation, in the regime where the SM is decoupled from (or weakly coupled to) the inflationary sector. We use classical lattice simulations in an expanding box in (3+1) dimensions, modeling the SM gauge interactions with both global and Abelian-Higgs analogue scenarios. We consider different post-inflationary expansion rates. During inflation, the Higgs forms a condensate, which starts oscillating soon after inflation ends. Via nonperturbative effects, the oscillations lead to a fast decay of the Higgs into the SM species, transferring most of the energy into $Z$ and $W^{\pm}$ bosons. All species are initially excited far away from equilibrium, but their interactions lead them into a stationary stage, with exact equipartition among the different energy components. From there on the system eventually reaches equilibrium. We have characterized in detail, in the different expansion histories considered, the evolution of the Higgs and of its dominant decay products, until equipartition is established. We provide a useful mapping between simulations with different parameters, from where we derive a master formula for the Higgs decay time, as a function of the coupling constants, Higgs initial amplitude and postinflationary expansion rate.

We provide an explicit expression for the renormalized expectation value of the stress-energy tensor of a spin-$1/2$ field in a spatially flat FLRW universe. Its computation is based on the extension of the adiabatic regularization method to fermion fields introduced recently in the literature. The tensor is given in terms of UV-finite integrals in momentum space, which involve the mode functions that define the quantum state. As illustrative examples of the method efficiency, we see how to compute the renormalized energy density and pressure in two interesting cosmological scenarios: a de Sitter spacetime and a radiation-dominated universe. In the second case, we explicitly show that the late-time renormalized stress-energy tensor behaves as that of classical cold matter. We also check that, if we obtain the adiabatic expansion of the scalar field mode functions with a similar procedure to the one used for fermions, we recover the well-known WKB-type expansion.

The extension of the adiabatic regularization method to spin-$1/2$ fields requires a self-consistent adiabatic expansion of the field modes. We provide here the details of such expansion, which differs from the WKB ansatz that works well for scalars, to firmly establish the generalization of the adiabatic renormalization scheme to spin-$1/2$ fields. We focus on the computation of particle production in de Sitter spacetime and obtain an analytic expression of the renormalized stress-energy tensor for Dirac fermions.

We extend the adiabatic regularization method to spin-1/2 fields. The ansatz for the adiabatic expansion for fermionic modes differs significantly from the WKB-type template that works for scalar modes. We give explicit expressions for the first adiabatic orders and analyze particle creation in de Sitter spacetime. As for scalar fields, the adiabatic method can be distinguished by its capability to overcome the UV divergences of the particle number operator. We also test the consistency of the extended method by working out the conformal and axial anomalies for a Dirac field in a FLRW spacetime, in exact agreement with those obtained from other renormalization prescriptions. We finally show its power by computing the renormalized stress-energy tensor for Dirac fermions in de Sitter space.