Abstract: We study the evolution of the energy distribution and equation of state of the Universe from the end of inflation until the onset of either radiation domination (RD) or a transient period of matter domination (MD). We use both analytical techniques and lattice simulations. We consider two-field models where the inflaton (/) has a monomial potential after inflation V((/)) proportional to i(/) – vip (p 4, and of order similar to 50% for p 4. The system goes to MD at late times for p = 2, while it goes to RD for p > 2. In the later case, we can calculate exactly the number of e-folds until RD as a function of g2, and hence predict accurately inflationary observables like the scalar tilt ns and the tensor-to-scalar ratio r. In the scenario (ii), the energy is always transferred completely to X for p > 2, as long as its effective mass m2X = g2((/) – v)2 is not negligible. For p = 2, the final ratio between the energy densities of X and (/) depends strongly on g2. For all p > 2, the system always goes to MD at late times.

Abstract: We extend the adiabatic regularization method for an expanding universe to include the Yukawa interaction between quantized Dirac fermions and a homogeneous background scalar field. We give explicit expressions for the renormalized expectation values of the stress-energy tensor < T-mu nu > and the bilinear <(psi) over bar psi > in a spatially flat Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime. These are basic ingredients in the semiclassical field equations of fermionic matter in curved spacetime interacting with a background scalar field. The ultraviolet subtracting terms of the adiabatic regularization can be naturally interpreted as coming from appropriate counterterms of the background fields. We fix the required covariant counterterms. To test our approach we determine the contribution of the Yukawa interaction to the conformal anomaly in the massless limit and show its consistency with the heat-kernel method using the effective action.

Abstract: We develop a new regularization method for the stress -energy tensor and the two -point function of free quantum scalar fields propagating in cosmological spacetimes. We proceed by extending the adiabatic regularization scheme with the introduction of two additional mass scales. By setting them to the order of the physical scale of the studied scenario, we obtain ultraviolet -regularized quantities that do not distort the power spectra amplitude at the infrared scales amplified by the expansion of the Universe. This is not ensured by the standard adiabatic approach. We also show how our proposed subtraction terms can be interpreted as a renormalization of coupling constants in the Einstein equations. We finally illustrate our proposed regularization method in two scenarios of cosmological interest: de Sitter inflation and geometric reheating.

Abstract: This paper describes CosmoGattice, a modern package for lattice simulations of the dynamics of interacting scalar and gauge fields in an expanding universe. CosmoGattice incorporates a series of features that makes it very versatile and powerful: i) it is written in C++ fully exploiting the object oriented programming paradigm, with a modular structure and a clear separation between the physics and the technical details, ii) it is MPI-based and uses a discrete Fourier transform parallelized in multiple spatial dimensions, which makes it specially appropriate for probing scenarios with well -separated scales, running very high resolution simulations, or simply very long ones, iii) it introduces its own symbolic language, defining field variables and operations over them, so that one can introduce differential equations and operators in a manner as close as possible to the continuum, iv) it includes a library of numerical algorithms, ranging from O(delta t(2)) to O(delta t(10)) methods, suitable for simulating global and gauge theories in an expanding grid, including the case of 'self-consistent' expansion sourced by the fields themselves. Relevant observables are provided for each algorithm (e.g. energy densities, field spectra, lattice snapshots) and we note that, remarkably, all our algorithms for gauge theories (Abelian or non-Abelian) always respect the Gauss constraint to machine precision. Program summary Program Title:: CosmoGattice CPC Library link to program files: https://doi .org /10 .17632 /44vr5xssc6 .1 Developer's repository link: http://github .com /cosmolattice /cosmolattice Licensing provisions: MIT Programming language: C++, MPI Nature of problem: The phenomenology of high energy physics in the early universe is typically characterized by non-linear dynamics, which cannot be captured accurately with analytical techniques. In order to fully understand the non-linearities developed in a given scenario, one needs to carry out lattice simulations. A number of public packages for lattice simulations have appeared over the years, but most of them are only capable of simulating scalar fields. However, realistic models of particle physics do contain other kind of field species, such as (Abelian or non-Abelian) gauge fields, whose non-linear dynamics can also play a relevant role in the early universe. Tensor modes representing gravitational waves are also naturally expected in many scenarios. Solution method: CosmoGattice represents a modern code for lattice simulations of scalar-gauge field theories in an expanding universe. It allows for the simulation of the evolution of interacting (singlet) scalar fields, charged scalar fields under U(1) and/or SU(2) gauge groups, and the corresponding associated Abelian and/or non-Abelian gauge fields. From version 1.1 onward, CosmoGattice also allows to simulate the production of gravitational waves. Simulations can be done either in a flat space-time background, or in a homogeneous and isotropic (spatially flat) expanding FLRW background. CosmoGattice provides symplectic integrators, with accuracy ranging from O (delta t(2)) up to O(delta t(10)), to simuate the non-linear dynamics of the appropriate fields in comoving three-dimensional lattices. The code is parallelized with MPI, and uses a discrete Fourier Transform parallelized in multiple spatial dimensions, which makes it a very powerful code for probing physical problems with well-separated scales. Moreover, the code has been designed as a `platform' to implement any system of dynamical equations suitable for discretization on a lattice.

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