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 revisit the topic of invisible neutrino decay in the precision cosmological context, via a first-principles approach to understanding the cosmic microwave background and large-scale structure phenomenology of such a non-standard physics scenario. Assuming an effective Lagrangian in which a heavier standard-model neutrino nu(H) couples to a lighter one nu(l) and a massless scalar particle phi via a Yukawa interaction, we derive from first principles the complete set of Boltzmann equations, at both the spatially homogeneous and the firstorder inhomogeneous levels, for the phase space densities of nu(H), nu(l), and phi in the presence of the relevant decay and inverse decay processes. With this set of equations in hand, we perform a critical survey of recent works on cosmological invisible neutrino decay in both limits of decay while nu(H) is ultra-relativistic and non-relativistic. Our two main findings are: (i) in the non-relativistic limit, the effective equations of motion used to describe perturbations in the neutrino-scalar system in the existing literature formally violate momentum conservation and gauge invariance, and (ii) in the ultra-relativistic limit, exponential damping of the anisotropic stress does not occur at the commonly-used rate Gamma(T) = (1/tau(0))( m(nu H)/E-nu H)(3), but at a rate similar to (1/ tau(0))(m(nu H)/E-nu H)(5). Both results are model-independent. The impact of the former finding on the cosmology of invisible neutrino decay is likely small. The latter, however, implies a significant revision of the cosmological limit on the neutrino lifetime tau(0) from tau(old)(0) greater than or similar to 1.2 x 10(9) s (m(nu H)/50 meV)(3) to tau(0) greater than or similar to (4 x 10(5) -> 4 x 10(6)) s (m(nu H)/50 meV)(5).

Abstract: Monte Carlo simulations are a unique tool to check the response of a detector and to monitor its performance. For a deep-sea neutrino telescope, the variability of the environmental conditions that can affect the behaviour of the data acquisition system must be considered, in addition to a reliable description of the active parts of the detector and of the features of physics events, in order to produce a realistic set of simulated events. In this paper, the software tools used to produce neutrino and cosmic ray signatures in the telescope and the strategy developed to represent the time evolution of the natural environment and of the detector efficiency are described.

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)).