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: R-parity conserving supersymmetric models with right-handed (RH) neutrinos are very appealing since they could naturally explain neutrino physics and also provide a good dark matter (DM) candidate such as the lightest supersymmetric particle (LSP). In this work we consider the next-to-minimal supersymmetric standard model (NMSSM) plus RH neutrino superfields, with effective Majorana masses dynamically generated at the electroweak scale (EW). We perform a scan of the relevant parameter space and study both possible DM candidates: RH sneutrino and neutralino. Especially for the case of RH sneutrino DM we analyse the intimate relation between both candidates to obtain the correct amount of relic density. Besides the well-known resonances, annihilations through scalar quartic couplings and coannihilation mechanisms with all kind of neutralinos, are crucial. Finally, we present the impact of current and future direct and indirect detection experiments on both DM candidates.

Abstract: We study the light rings and shadows of an uniparametric family of spherically symmetric geometries interpolating between the Schwarzschild solution, a regular black hole, and a traversable wormhole, and dubbed as black bounces, all of them sharing the same critical impact parameter. We consider the ray-tracing method in order to study the impact parameter regions corresponding to the direct, lensed, and photon ring emissions, finding a broadening of all these regions for black bounce solutions as compared to the Schwarzschild one. Using this, we determine the optical appearance of black bounces when illuminated by three standard toy models of optically and geometrically thin accretion disks viewed in face-on orientation.

Abstract: We construct an axially symmetric solution of Eddington-inspired Born-Infeld gravity coupled to an electromagnetic field in 2 + 1 dimensions including a (negative) cosmological constant term. This is achieved by using a recently developed mapping procedure that allows to generate solutions in certain families of metric-affine gravity theories starting from a known seed solution of General Relativity, which in the present case corresponds to the electrically charged Banados-Teitelboim-Zanelli (BTZ) solution. We discuss the main features of the new configurations, including the modifications to the ergospheres and horizons, the emergence of wormhole structures, and the consequences for the regularity (or not) of these space-times via geodesic completeness.

Abstract: The evolution of the universe prior to Big Bang Nucleosynthesis could have gone through a phase of early matter domination which enhanced the growth of small-scale dark matter structure. If this period was long enough, self-gravitating objects formed prior to reheating. We study the evolution of these dense early halos through reheating. At the end of early matter domination, the early halos undergo rapid expansion and eventually eject their matter. We find that this process washes out structure on scales much larger than naively expected from the size of the original halos. We compute the density profiles of the early halo remnants and use them to construct late-time power spectra that include these non-linear effects. We evolve the resulting power spectrum to estimate the properties of microhalos that would form after matter-radiation equality. Surprisingly, cosmologies with a short period of early matter domination lead to an earlier onset of microhalo formation compared to those with a long period. In either case, dark matter structure formation begins much earlier than in the standard cosmology, with most dark matter bound in microhalos in the late universe.