-
Present and future of CosmoLattice
Authors:
Daniel G. Figueroa,
Adrien Florio,
Francisco Torrenti
Abstract:
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 impleme…
▽ More
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 $φF \widetilde{F}$, non-minimal gravitational couplings $φ^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.
△ Less
Submitted 22 December, 2023;
originally announced December 2023.
-
Physical scale adiabatic regularization in cosmological spacetimes
Authors:
Antonio Ferreiro,
Samuel Monin,
Francisco Torrenti
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 t…
▽ More
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.
△ Less
Submitted 24 February, 2024; v1 submitted 15 November, 2023;
originally announced November 2023.
-
Ultraviolet-regularized power spectrum without infrared distortions in cosmological spacetimes
Authors:
Antonio Ferreiro,
Francisco Torrenti
Abstract:
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…
▽ More
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 $Δ_φ^{\rm (reg)} = 0$ for a massless field, our proposed prescription recovers the standard scale-invariant result $Δ_φ^{\rm (reg)} \simeq H^2 /(4π^2)$ at super-horizon scales.
△ Less
Submitted 29 March, 2023; v1 submitted 2 December, 2022;
originally announced December 2022.
-
Characterizing the post-inflationary reheating history, Part II: Multiple interacting daughter fields
Authors:
Stefan Antusch,
Kenneth Marschall,
Francisco Torrenti
Abstract:
We characterize the post-inflationary dynamics of an inflaton $φ$ coupled to multiple interacting daughter fields $X_n$ ($n=1,\dots N_d$) through quadratic-quadratic interactions $g_n^ 2φ^2 X_n^2$. We assume a monomial inflaton potential $V(φ) \propto |φ|^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…
▽ More
We characterize the post-inflationary dynamics of an inflaton $φ$ coupled to multiple interacting daughter fields $X_n$ ($n=1,\dots N_d$) through quadratic-quadratic interactions $g_n^ 2φ^2 X_n^2$. We assume a monomial inflaton potential $V(φ) \propto |φ|^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 $λ_{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$.
△ Less
Submitted 14 February, 2023; v1 submitted 13 June, 2022;
originally announced June 2022.
-
Characterizing the post-inflationary reheating history, Part I: single daughter field with quadratic-quadratic interaction
Authors:
Stefan Antusch,
Daniel G. Figueroa,
Kenneth Marschall,
Francisco Torrenti
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(Φ) \propto |Φ- v|^p$ (…
▽ More
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(Φ) \propto |Φ- v|^p$ ($p\geq2$), and is coupled to a daughter field $X$ through a quadratic-quadratic interaction $g^2Φ^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(Φ-v)^2$ is not negligible. For $p=2$, the final ratio between the energy densities of $X$ and $Φ$ depends strongly on $g^2$. For all $p \ge 2$, the system always goes to MD at late times.
△ Less
Submitted 8 June, 2022; v1 submitted 21 December, 2021;
originally announced December 2021.
-
CosmoLattice
Authors:
Daniel G. Figueroa,
Adrien Florio,
Francisco Torrenti,
Wessel Valkenburg
Abstract:
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 p…
▽ More
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(δt^2)$ to $O(δ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.
△ Less
Submitted 6 May, 2022; v1 submitted 1 February, 2021;
originally announced February 2021.
-
Challenges and Opportunities of Gravitational Wave Searches at MHz to GHz Frequencies
Authors:
N. Aggarwal,
O. D. Aguiar,
A. Bauswein,
G. Cella,
S. Clesse,
A. M. Cruise,
V. Domcke,
D. G. Figueroa,
A. Geraci,
M. Goryachev,
H. Grote,
M. Hindmarsh,
F. Muia,
N. Mukund,
D. Ottaway,
M. Peloso,
F. Quevedo,
A. Ricciardone,
J. Steinlechner,
S. Steinlechner,
S. Sun,
M. E. Tobar,
F. Torrenti,
C. Unal,
G. White
Abstract:
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, with a particular focus on Ultra High-Frequency Gravitational Waves (UHF-GWs), covering the MHz to GHz range. The absence of kn…
▽ More
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, with a particular focus on Ultra High-Frequency Gravitational Waves (UHF-GWs), covering the MHz to GHz range. The absence of known astrophysical sources in 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 gravitational sources. We review several detector concepts which 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 the workshop "Challenges and opportunities of high-frequency gravitational wave detection" held at ICTP Trieste, Italy in October 2019, that set up the stage for the recently launched Ultra-High-Frequency Gravitational Wave (UHF-GW) initiative.
△ Less
Submitted 13 December, 2021; v1 submitted 24 November, 2020;
originally announced November 2020.
-
The art of simulating the early Universe -- Part I
Authors:
Daniel G. Figueroa,
Adrien Florio,
Francisco Torrenti,
Wessel Valkenburg
Abstract:
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 class…
▽ More
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
△ Less
Submitted 14 May, 2021; v1 submitted 26 June, 2020;
originally announced June 2020.
-
Energy distribution and equation of state of the early Universe: matching the end of inflation and the onset of radiation domination
Authors:
Stefan Antusch,
Daniel G. Figueroa,
Kenneth Marschall,
Francisco Torrenti
Abstract:
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(φ) \propto |φ|^p$ around the origin. As a proxy for (p)reheating, we include a quadratic interaction…
▽ More
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(φ) \propto |φ|^p$ around the origin. As a proxy for (p)reheating, we include a quadratic interaction $g^2φ^2X^2$ between the inflaton $φ$ 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$.
△ Less
Submitted 22 December, 2020; v1 submitted 15 May, 2020;
originally announced May 2020.
-
Properties of oscillons in hilltop potentials: energies, shapes, and lifetimes
Authors:
Stefan Antusch,
Francesco Cefala,
Francisco Torrenti
Abstract:
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-shape…
▽ More
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.
△ Less
Submitted 9 October, 2019; v1 submitted 1 July, 2019;
originally announced July 2019.
-
Higgs-curvature coupling and post-inflationary vacuum instability
Authors:
Daniel G. Figueroa,
Arttu Rajantie,
Francisco Torrenti
Abstract:
We study the post-inflationary dynamics of the Standard Model (SM) Higgs field in the presence of a non-minimal coupling $ξ|Φ|^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 sca…
▽ More
We study the post-inflationary dynamics of the Standard Model (SM) Higgs field in the presence of a non-minimal coupling $ξ|Φ|^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 $ξ$, 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 $ξ\gtrsim 4 -5$, although a slightly lower mass of $m_t \approx 172.1 {\rm GeV}$ pushes up this limit to $ξ\gtrsim 11 - 12$. This, together with the estimation $ξ\gtrsim 0.06$ for stability during inflation, provides tight constraints to the Higgs-curvature coupling within the SM.
△ Less
Submitted 12 September, 2018; v1 submitted 1 September, 2017;
originally announced September 2017.
-
Gravitational wave production from preheating: parameter dependence
Authors:
Daniel G. Figueroa,
Francisco Torrenti
Abstract:
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. Pr…
▽ More
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Ω_{\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.
△ Less
Submitted 16 November, 2017; v1 submitted 14 July, 2017;
originally announced July 2017.
-
Parametric Resonance in the Early Universe - A Fitting Analysis
Authors:
Daniel G. Figueroa,
Francisco Torrenti
Abstract:
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 ci…
▽ More
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.
△ Less
Submitted 1 February, 2017; v1 submitted 16 September, 2016;
originally announced September 2016.
-
Gravitational wave production from the decay of the Standard Model Higgs field after inflation
Authors:
Daniel G. Figueroa,
Juan García-Bellido,
Francisco Torrentí
Abstract:
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…
▽ More
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 Ω_{\rm GW}^{(o)} \lesssim 10^{-29}$ today. The amplitude can be boosted to $h^2 Ω_{\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.
△ Less
Submitted 31 May, 2016; v1 submitted 9 February, 2016;
originally announced February 2016.
-
The Decay of the Standard Model Higgs after Inflation
Authors:
Daniel G. Figueroa,
Juan Garcia-Bellido,
Francisco Torrenti
Abstract:
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.…
▽ More
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.
△ Less
Submitted 26 October, 2015; v1 submitted 17 April, 2015;
originally announced April 2015.
-
Renormalized stress-energy tensor for spin-1/2 fields in expanding universes
Authors:
Adrian del Rio,
Jose Navarro-Salas,
Francisco Torrenti
Abstract:
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…
▽ More
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.
△ Less
Submitted 30 September, 2014; v1 submitted 18 July, 2014;
originally announced July 2014.
-
Adiabatic regularization and particle creation for spin one-half fields
Authors:
Aitor Landete,
Jose Navarro-Salas,
Francisco Torrenti
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 producti…
▽ More
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.
△ Less
Submitted 26 February, 2014; v1 submitted 20 November, 2013;
originally announced November 2013.
-
Adiabatic regularization for spin-1/2 fields
Authors:
Aitor Landete,
Jose Navarro-Salas,
Francisco Torrenti
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 t…
▽ More
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.
△ Less
Submitted 13 September, 2013; v1 submitted 31 May, 2013;
originally announced May 2013.