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 the 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(delta t(2)) accuracy like staggered leapfrog and Verlet integration, to Runge-Kutta methods up to O(delta t(4)) accuracy, and the Yoshida and Gauss-Legendre higher-order integrators, accurate up to O(delta t(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(delta t(2)) up to O(delta t(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, 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) x U(1). We note that the present manuscript is meant to be 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, publicly available at http://www.cosrnolattice.net.

Abstract: Axion inflation, i.e. an axion-like inflaton coupled to an Abelian gauge field through a Chern-Simons interaction, comes with a rich and testable phenomenology. This is particularly true in the strong backreaction regime, where the gauge field production heavily impacts the axion dynamics. Lattice simulations have recently demonstrated the importance of accounting for inhomogeneities of the axion field in this regime. We propose a perturbative scheme to account for these inhomogeneities while maintaining high computational efficiency. Our goal is to accurately capture deviations from the homogeneous axion field approximation within the perturbative regime as well as self -consistently determine the onset of the nonperturbative regime.

Abstract: Context. The existence of a maximum correlation angle (theta(max) & 60 greater than or similar to degrees) in the two-point angular temperature correlations of cosmic microwave background (CMB) radiation, measured by WMAP and Planck, stands in sharp contrast to the prediction of standard inflationary cosmology, in which the correlations should extend across the full sky (i.e., 180 degrees). The introduction of a hard lower cuto ff (k(min)) in the primordial power spectrum, however, leads naturally to the existence of theta(max). Among other cosmological anomalies detected in these data, an apparent dominance of odd-over-even parity multipoles has been seen in the angular power spectrum of the CMB. This feature, however, may simply be due to observational contamination in certain regions of the sky. Aims. In attempting to provide a more detailed assessment of whether this odd-over-even asymmetry is intrinsic to the CMB, we therefore proceed in this paper, first, to examine whether this odd-even parity imbalance also manifests itself in the angular correlation function and, second, to examine in detail the interplay between the presence of theta(max) and this observed anomaly. Methods. We employed several parity statistics and recalculated the angular correlation function for di fferent values of the cuto ff kmin in order to optimize the fit to the di fferent Planck 2018 data. Results. We find a phenomenological connection between these features in the data, concluding that both must be considered together in order to optimize the theoretical fit to the Planck 2018 data. Conclusions. This outcome is independent of whether the parity imbalance is intrinsic to the CMB, but if it is, the odd-over-even asymmetry would clearly point to the emergence of new physics.