Abstract: Dark photon dark matter will resonantly convert into visible photons when the dark photon mass is equal to the plasma frequency of the ambient medium. In cosmological contexts, this transition leads to an extremely efficient, albeit short-lived, heating of the surrounding gas. Existing work in this field has been predominantly focused on understanding the implications of these resonant transitions in the limit that the plasma frequency of the Universe can be treated as being perfectly homogeneous, i.e. neglecting inhomogeneities in the electron number density. In this work we focus on the implications of heating from dark photon dark matter in the presence of inhomogeneous structure (which is particularly relevant for dark photons with masses in the range 10(-15) eV less than or similar to m(A ') less than or similar to 10(-12) eV), emphasizing both the importance of inhomogeneous energy injection, as well as the sensitivity of cosmological observations to the inhomogeneities themselves. More specifically, we derive modified constraints on dark photon dark matter from the Ly-alpha forest, and show that the presence of inhomogeneities allows one to extend constraints to masses outside of the range that would be obtainable in the homogeneous limit, while only slightly relaxing their strength. We then project sensitivity for near-future cosmological surveys that are hoping to measure the 21cm transition in neutral hydrogen prior to reionization, and demonstrate that these experiments will be extremely useful in improving sensitivity to masses near similar to 10(-14) eV, potentially by several orders of magnitude. Finally, we discuss implications for reionization, early star formation, and late-time y-type spectral distortions, and show that probes which are inherently sensitive to the inhomogeneous state of the Universe could resolve signatures unique to the light dark photon dark matter scenario, and thus offer a fantastic potential for a positive detection.

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.