Search

We investigate the post-inflationary evolution of a non-minimally coupled inflaton field in scalar-tensor theories, framed within the flexible framework of Einstein-Cartan gravity. By focusing on a class of simplified Higgs-like scenarios, we simulate the transition from the end of inflation to the formation of oscillons using fully-fledged 3+1 classical lattice simulations. Once oscillons form, we extract their profiles and perform 1+1 simulations to evolve their radial equations. Our findings reveal that these oscillons, unlike typical cases in the literature, are relatively short-lived, due to the presence of self interactions at small field values. The radiation produced by this novel type of oscillons can quickly lead to a radiation-dominated Universe, even in the absence of additional fields or interactions. Finally, we leverage our results to derive precise predictions for the inflationary observables and the produced spectrum of gravitational waves associated with oscillon formation. Significantly, our study establishes also an upper bound on the duration of the heating phase in Einstein-Cartan Higgs inflation scenarios., Comment: 29 pages, 10 figures

The first direct measurement of gravitational waves by the LIGO and Virgo collaborations has opened up new avenues to explore our Universe. This white paper outlines the challenges and gains expected in gravitational-wave searches at frequencies above the LIGO/Virgo band. The scarcity of possible astrophysical sources in most of this frequency range provides a unique opportunity to discover physics beyond the Standard Model operating both in the early and late Universe, and we highlight some of the most promising of these sources. We review several detector concepts that have been proposed to take up this challenge, and compare their expected sensitivity with the signal strength predicted in various models. This report is the summary of a series of workshops on the topic of high-frequency gravitational wave detection, held in 2019 (ICTP, Trieste, Italy), 2021 (online) and 2023 (CERN, Geneva, Switzerland)., Comment: 138 pages, 13 figures. Update and extension of Living Rev$.$Rel. 24 (2021) 1, 4 on high frequency gravitational waves (arxiv: 2011.12414). Comments welcome

We discuss the present state and planned updates of CosmoLattice, a cutting-edge code for lattice simulations of non-linear dynamics of scalar-gauge field theories in an expanding background. We first review current capabilities of the code, including the simulation of interacting singlet scalars and of Abelian and non-Abelian scalar-gauge theories. We also comment on new features recently implemented, such as the simulation of gravitational waves from scalar and gauge fields. Secondly, we discuss new extensions of CosmoLattice that we plan to release publicly. On the one hand, we comment on new physics modules, which include axion-gauge interactions $\phi F \widetilde{F}$, non-minimal gravitational couplings $\phi^2 R$, creation and evolution of cosmic defect networks, and magneto-hydro-dynamics (MHD). On the other hand, we discuss new technical features, including evolvers for non-canonical interactions, arbitrary initial conditions, simulations in 2+1 dimensions, and higher accuracy spatial derivatives., Comment: Invited article for Reports on Progress in Physics, matches published version. 20 pages + references

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., Comment: It matches published version in PRD. 22 pages + references, 5 figures

We reexamine the regularization of the two-point function of a scalar field in a Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime. Adiabatic regularization provides a set of subtraction terms in momentum space that successfully remove its ultraviolet divergences at coincident points, but can significantly distort the power spectrum at infrared scales, especially for light fields. In this work we propose, by using the intrinsic ambiguities of the renormalization program, a new set of subtraction terms that minimize the distortions for scales $k \lesssim M$, with $M$ an arbitrary mass scale. Our method is consistent with local covariance and equivalent to general regularization methods in curved spacetime. We apply our results to the regularization of the power spectrum in de Sitter space: while the adiabatic scheme yields exactly $\Delta_{\phi}^{\rm (reg)} = 0$ for a massless field, our proposed prescription recovers the standard scale-invariant result $\Delta_{\phi}^{\rm (reg)} \simeq H^2 /(4\pi^2)$ at super-horizon scales., Comment: Title changed with respect to first version. New section added on renormalization conditions and coupling constants (Sect. 3). It matches the published version in PLB. 6 pages + references, 1 figure

We characterize the post-inflationary dynamics of an inflaton $\phi$ coupled to multiple interacting daughter fields $X_n$ ($n=1,\dots N_d$) through quadratic-quadratic interactions $g_n^ 2\phi^2 X_n^2$. We assume a monomial inflaton potential $V(\phi) \propto |\phi|^p$ ($p \geq 2$) around the minimum. By simulating the system in 2+1-dimensional lattices, we study the post-inflationary evolution of the energy distribution and equation of state, from the end of inflation until a stationary regime is achieved. We show that in this scenario, the energy transferred to the daughter field sector can be larger than $50\%$, surpassing this way the upper bound found previously for single daughter field models. In particular, for $p \geq 4$ the energy at very late times is equally distributed between all fields, and only $100/(N_d + 1) \%$ of the energy remains in the inflaton. We also consider scenarios in which the daughter fields have scale-free interactions $\lambda_{nm} X_n^2 X_m^2$, including the case of quartic daughter field self-interactions (for $n=m$). We show that these interactions trigger a resonance process during the non-linear regime, which in the single daughter field case already allows to deplete more than $50\%$ of the energy from the inflaton for $p\geq 4$., Comment: 23 pages + appendix, 15 figures

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 $\Phi$ has a monomial potential after inflation $V(\Phi) \propto |\Phi - v|^p$ ($p\geq2$), and is coupled to a daughter field $X$ through a quadratic-quadratic interaction $g^2\Phi^2 X^2$. We consider two situations, depending on whether the potential has a minimum at $i)$ $v = 0$, or $ii)$ $v > 0$. In the scenario $i)$, the final energy transferred to $X$ is independent of $g^2$ and entirely determined by $p$: it is negligible for $p < 4$, and of order $\sim 50\%$ for $p \geq 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 $g^2$, and hence predict accurately inflationary observables like the scalar tilt $n_s$ 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 $m_X^2 = g^2(\Phi-v)^2$ is not negligible. For $p=2$, the final ratio between the energy densities of $X$ and $\Phi$ depends strongly on $g^2$. For all $p \ge 2$, the system always goes to MD at late times., Comment: References added. Minor changes to match published version. 29 pages + appendices, 23 figures

This is the user manual for CosmoLattice, a modern package for lattice simulations of the dynamics of interacting scalar and gauge fields in an expanding universe. CosmoLattice 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 always respect the Gauss constraint to machine precision. In this manual we explain how to obtain and run CosmoLattice in a computer (let it be your laptop, desktop or a cluster). We introduce the general structure of the code and describe in detail the basic files that any user needs to handle. We explain how to implement any model characterized by a scalar potential and a set of scalar fields, either singlets or interacting with $U(1)$ and/or $SU(2)$ gauge fields. CosmoLattice is publicly available at www.cosmolattice.net., Comment: 111 pages, 3 figures and O(100) code files. v1.1 is out, with GW!

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. Remarkably, all our algorithms for gauge theories respect the Gauss constraint to machine precision, including when `self-consistent' expansion is considered. As a numerical example we analyze the post-inflationary dynamics of an oscillating inflaton charged under $SU(2)\times U(1)$. The present manuscript is meant as part of the theoretical basis for $CosmoLattice$, a modern C++ MPI-based package for simulating the non-linear dynamics of scalar-gauge field theories in an expanding universe, publicly available at www.cosmolattice.net, Comment: Minor corrections to match published version, and one more algorithm added. Still 79 pages, 8 figures, 1 appendix, and many algorithms

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. As a proxy for (p)reheating, we include a quadratic interaction $g^2\phi^2X^2$ between the inflaton $\phi$ and a light scalar 'daughter' field $X$, with $g^2>0$. 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$., Comment: 7 pages + references, 5 figures. It matches published version

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., Comment: 22 pages, 11 figures. Minor changes to match version published in JCAP

We extend the adiabatic regularization method for an expanding universe to include the Yukawa interaction between a quantized Dirac field and a homogeneous time-dependent scalar field. We present the renormalized semiclassical equations that are needed in order to take into account the backreaction of the produced Dirac fermions in both gravitational and scalar background fields., Comment: 6 pages. Contribution to the Marcel Grossmann Meeting 2018

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., Comment: 15 pages, 13 figures. Minor changes to match version published in PRD

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., Comment: Corrected typos, matches version published in JCAP (39 pages + appendix, 11 figures)

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., Comment: Extended discussion about the late-time dynamics of the system in quadratic models. Minor changes in numerical fits with respect first version. It matches version published in JCAP (30 pages + Appendices + Bibliography, 13 figures)

During or towards the end of inflation, the Standard Model (SM) Higgs forms a condensate with a large amplitude. Following inflation, the condensate oscillates, decaying non-perturbatively into the rest of the SM species. The resulting out-of-equilibrium dynamics converts a fraction of the energy available into gravitational waves (GW). We study this process using classical lattice simulations in an expanding box, following the energetically dominant electroweak gauge bosons $W^\pm$ and $Z$. We characterize the GW spectrum as a function of the running couplings, Higgs initial amplitude, and post-inflationary expansion rate. As long as the SM is decoupled from the inflationary sector, the generation of this background is universally expected, independently of the nature of inflation. Our study demonstrates the efficiency of GW emission by gauge fields undergoing parametric resonance. The initial energy of the Higgs condensate represents however, only a tiny fraction of the inflationary energy. Consequently, the resulting background is very suppressed, with an amplitude $h^2 \Omega_{\rm GW}^{(o)} \lesssim 10^{-29}$ today. The amplitude can be boosted to $h^2 \Omega_{\rm GW}^{(o)} \lesssim 10^{-16}$, if following inflation the universe undergoes a kination-domination stage; however the background is shifted in this case to high frequencies $f_p \lesssim 10^{11} {\rm Hz}$. In all cases the signal is out of the range of current or planned GW detectors. This background will therefore remain, most likely, as a curiosity of the SM., Comment: 16 pages, 6 figures. Minor changes to match version published in PRD

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., Comment: Minor changes to match the PRD published version. Modulation of the Higgs amplitude removed for $q > 200$ in Sec. V, due to improving the time resolution in the Higgs equation of motion. Results unaffected

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., Comment: Latex file, 26 pages

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., Comment: Some sections has been reorganized. New references added. Minor changes. Slight changes in title and abstract

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., Comment: A new section has been added