GFiRe: a Gauge Field integrator for Reheating

Abstract

We present a new numerical algorithm and code, $\mathsf{GFiRe}$, for solving the non-linear evolution of Abelian gauge fields coupled to complex scalar fields in homogeneous and isotropic spacetimes. We adopt a hybrid approach to solving the system: the spatial derivatives are discretized using standard Lattice Gauge Field Theory techniques, whereas the time evolution of the fields and scale factor is implemented with explicit, composite, symplectic integrators. An important property of our compound algorithm is that the discretized Gauss constraint is respected exactly, regardless of the order of the symplectic integrator. This remains true even when the background expansion is computed "self-consistently"; that is, when the expansion history is computed using spatial averaged components of the energy momentum tensor in the simulation volume. Hence, our code can also be used in cases where the fields dominate the energy density of the universe, for example, during reheating after inflation. We test the algorithm in scenarios of reheating where the inflaton is a complex scalar field with a potential $\propto (2|\phi {|}^2-v^2{)}^2$ and is coupled to an Abelian gauge field. Tracing the evolution of the system through complex dynamics (including resonant excitation of fields, backreaction, formation of solitons, and changes in the equation of state) in a self-consistently expanding universe, we find the energy conservation violation ($<10^{-4}$) to be very stable and the Gauss constraint violation ($<10^{-6}$) to be dominated by differencing noise.

35 pages, 13 figures; minor changes, accepted for publication in JCAP

Cited by in corpus (16)

• The art of simulating the early Universe -- Part I
• Lattice Simulations of Axion-U(1) Inflation
• Freeze-in from Preheating
• Inflationary Butterfly Effect: Non-perturbative Dynamics From Small-Scale Features
• Gravitational wave background from non-Abelian reheating after axion-like inflation
• Lattice simulations of Abelian gauge fields coupled to axions during Inflation
• Unveiling the nonlinear dynamics of a rolling axion during inflation
• Lattice Simulations of Inflation
• i-SPin: An integrator for multicomponent Schrödinger-Poisson systems with self-interactions
• Experimental targets for dark photon dark matter
• Present and future of CosmoLattice
• On the Expected Backreaction During Preheating
• Impact of Dark Sector Preheating on CMB Observables
• Simple third order operator-splitting schemes for stochastic mechanics and field theory
• Non-Gaussianity after many-field reheating
• Dark Matter Freeze-in from a $Z^{\prime }$ Reheaton