Abstract: We explore the constraints and phenomenology of possibly the simplest scenario that could account at the same time for the active neutrino masses and the dark matter in the Universe within a gauged U(1)(B-L) symmetry, namely right-handed neutrino dark matter. We find that null searches from lepton and hadron colliders require dark matter with a mass below 900 GeV to annihilate through a resonance. Additionally, the very strong constraints from high-energy dilepton searches fully exclude the model for 150 GeV < m(z') < 3 TeV. We further explore the phenomenology in the high mass region (i.e. masses greater than or similar to O(1) TeV) and highlight theoretical arguments, related to the appearance of a Landau pole or an instability of the scalar potential, disfavoring large portions of this parameter space. Collectively, these considerations illustrate that a minimal extension of the Standard Model via a local U(1)(B-L) symmetry with a viable thermal dark matter candidate is difficult to achieve without fine-tuning. We conclude by discussing possible extensions of the model that relieve tension with collider constraints by reducing the gauge coupling required to produce the correct relic abundance.

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 discuss the possibility to predict the QCD axion mass in the context of grand unified theories. We investigate the implementation of the DFSZ mechanism in the context of renormalizable SU(5) theories. In the simplest theory, the axion mass can be predicted with good precision in the range m(a) = (2-16) neV, and there is a strong correlation between the predictions for the axion mass and proton decay rates. In this context, we predict an upper bound for the proton decay channels with antineutrinos, tau(p -> K+(nu) over bar) less than or similar to 4 x 10(37) yr and tau(p -> pi(+)(nu) over bar) less than or similar to 2 x 10(36) yr. This theory can be considered as the minimal realistic grand unified theory with the DFSZ mechanism and it can be fully tested by proton decay and axion experiments.