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: The Sun may capture asymmetric dark matter (DM), which can subsequently form bound-states through the radiative emission of a sub-GeV scalar. This process enables generation of scalars without requiring DM annihilation. In addition to DM capture on nucleons, the DM-scalar coupling responsible for bound-state formation also induces capture from self-scatterings of ambient DM particles with DM particles already captured, as well as with DM bound-states formed in-situ within the Sun. This scenario is studied in detail by solving Boltzmann equations numerically and analytically. In particular, we take into consideration that the DM self-capture rates require a treatment beyond the conventional Born approximation. We show that, thanks to DM scatterings on bound-states, the number of DM particles captured increases exponentially, leading to enhanced emission of relativistic scalars through bound-state formation, whose final decay products could be observable. We explore phenomenological signatures with the example that the scalar mediator decays to neutrinos. We find that the neutrino flux emitted can be comparable to atmospheric neutrino fluxes within the range of energies below one hundred MeV. Future facilities like Hyper-K, and direct DM detection experiments can further test such scenario.