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: We find an exact, rotating charged black hole solution within Eddington-inspired Born-Infeld gravity. To this end we employ a recently developed correspondence or mapping between modified gravity models built as scalars out of contractions of the metric with the Ricci tensor, and formulated in metric-affine spaces (Ricci-Based Gravity theories) and General Relativity. This way, starting from the Kerr-Newman solution, we show that this mapping bring us the axisymmetric solutions of Eddington-inspired Born-Infeld gravity coupled to a certain model of non-linear electrodynamics. We discuss the most relevant physical features of the solutions obtained this way, both in the spherically symmetric limit and in the fully rotating regime. Moreover, we further elaborate on the potential impact of this important technical progress for bringing closer the predictions of modified gravity with the astrophysical observations of compact objects and gravitational wave astronomy.

Abstract: The observation of space-time variations in fundamental constants would provide strong evidence for the existence of new light degrees of freedom in the theory of Nature. Robustly constraining such scenarios requires exploiting observations that span different scales and probe the state of the Universe at different epochs. In the context of cosmology, both the cosmic microwave background and the Lyman-a forest have proven to be powerful tools capable of constraining variations in electromagnetism, however at the moment there do not exist cosmological probes capable of bridging the gap between recombination and reionization. In the near future, radio telescopes will attempt to measure the 21 cm transition of neutral hydrogen during the epochs of reionization and the cosmic dawn (and potentially the tail end of the dark ages); being inherently sensitive to electromagnetic phenomena, these experiments will offer a unique perspective on space-time variations of the fine-structure constant and the electron mass. We show here that large variations in these fundamental constants would produce features on the 21 cm power spectrum that may be distinguishable from astrophysical uncertainties. Furthermore, we forecast the sensitivity for the Square Kilometer Array, and show that the 21 cm power spectrum may be able to constrain variations at the level of O(10(-3)).

Abstract: The standard cosmological model predicts the existence of a Cosmic Neutrino Background, which has not yet been observed directly. Some experiments aiming at its detection are currently under development, despite the tiny kinetic energy of the cosmological relic neutrinos, which makes this task incredibly challenging. Since massive neutrinos are attracted by the gravitational potential of our Galaxy, they can cluster locally. Neutrinos should be more abundant at the Earth position than at an average point in the Universe. This fact may enhance the expected event rate in any future experiment. Past calculations of the local neutrino clustering factor only considered a spherical distribution of matter in the Milky Way and neglected the influence of other nearby objects like the Virgo cluster, although recent N-body simulations suggest that the latter may actually be important. In this paper, we adopt a back-tracking technique, well established in the calculation of cosmic rays fluxes, to perform the first three-dimensional calculation of the number density of relic neutrinos at the Solar System, taking into account not only the matter composition of the Milky Way, but also the contribution of the Andromeda galaxy and the Virgo cluster. The effect of Virgo is indeed found to be relevant and to depend non-trivially on the value of the neutrino mass. Our results show that the local neutrino density is enhanced by 0.53% for a neutrino mass of 10 meV, 12% for 50 meV, 50% for 100 meV or 500% for 300 meV.

Abstract: Although weakly interacting massive particles (WIMPs) have long been among the most studied and theoretically attractive classes of candidates for the dark matter of our universe, the lack of their detection in direct detection and collider experiments has begun to dampen enthusiasm for this paradigm. In this study, we set out to appraise the status of the WIMP paradigm, focusing on the case of dark matter candidates that interact with the Standard Model through a new gauge boson. After considering a wide range of Z' mediated dark matter models, we quantitatively evaluate the fraction of the parameter space that has been excluded by existing experiments, and that is projected to fall within the reach of future direct detection experiments. Despite the existence of stringent constraints, we find that a sizable fraction of this parameter space remains viable. More specifically, if the dark matter is a Majorana fermion, we find that an order one fraction of the parameter space is in many cases untested by current experiments. Future direct detection experiments with sensitivity near the irreducible neutrino floor will be able to test a significant fraction of the currently viable parameter space, providing considerable motivation for the next generation of direct detection experiments.