Abstract: We revisit the generation of a matter-antimatter asymmetry in the minimal extension of the Standard Model with two singlet heavy neutral leptons (HNL) that can explain neutrino masses. We derive an accurate analytical approximation to the solution of the complete linearized set of kinetic equations, which exposes the non-trivial parameter dependencies in the form of parameterization-independent CP invariants. The identification of various washout regimes relevant in different regions of parameter space sheds light on the relevance of the mass corrections in the interaction rates and clarifies the correlations of baryogenesis with other observables. In particular, by requiring that the measured baryon asymmetry is reproduced, we derive robust upper or lower bounds on the HNL mixings depending on their masses, and constraints on their flavour structure, as well as on the CP-violating phases of the PMNS mixing matrix, and the amplitude of neutrinoless double-beta decay. We also find certain correlations between low and high scale CP phases. Especially emphasizing the testable part of the parameter space we demonstrate that our findings are in very good agreement with numerical results. The methods developed in this work can help in exploring more complex scenarios.

Abstract: This paper describes CosmoGattice, a modern package for lattice simulations of the dynamics of interacting scalar and gauge fields in an expanding universe. CosmoGattice incorporates a series of features that makes it very versatile and powerful: i) it is written in C++ fully exploiting the object oriented programming paradigm, with a modular structure and a clear separation between the physics and the technical details, ii) it is MPI-based and uses a discrete Fourier transform parallelized in multiple spatial dimensions, which makes it specially appropriate for probing scenarios with well -separated scales, running very high resolution simulations, or simply very long ones, iii) it introduces its own symbolic language, defining field variables and operations over them, so that one can introduce differential equations and operators in a manner as close as possible to the continuum, iv) it includes a library of numerical algorithms, ranging from O(delta t(2)) to O(delta t(10)) methods, suitable for simulating global and gauge theories in an expanding grid, including the case of 'self-consistent' expansion sourced by the fields themselves. Relevant observables are provided for each algorithm (e.g. energy densities, field spectra, lattice snapshots) and we note that, remarkably, all our algorithms for gauge theories (Abelian or non-Abelian) always respect the Gauss constraint to machine precision. Program summary Program Title:: CosmoGattice CPC Library link to program files: https://doi .org /10 .17632 /44vr5xssc6 .1 Developer's repository link: http://github .com /cosmolattice /cosmolattice Licensing provisions: MIT Programming language: C++, MPI Nature of problem: The phenomenology of high energy physics in the early universe is typically characterized by non-linear dynamics, which cannot be captured accurately with analytical techniques. In order to fully understand the non-linearities developed in a given scenario, one needs to carry out lattice simulations. A number of public packages for lattice simulations have appeared over the years, but most of them are only capable of simulating scalar fields. However, realistic models of particle physics do contain other kind of field species, such as (Abelian or non-Abelian) gauge fields, whose non-linear dynamics can also play a relevant role in the early universe. Tensor modes representing gravitational waves are also naturally expected in many scenarios. Solution method: CosmoGattice represents a modern code for lattice simulations of scalar-gauge field theories in an expanding universe. It allows for the simulation of the evolution of interacting (singlet) scalar fields, charged scalar fields under U(1) and/or SU(2) gauge groups, and the corresponding associated Abelian and/or non-Abelian gauge fields. From version 1.1 onward, CosmoGattice also allows to simulate the production of gravitational waves. Simulations can be done either in a flat space-time background, or in a homogeneous and isotropic (spatially flat) expanding FLRW background. CosmoGattice provides symplectic integrators, with accuracy ranging from O (delta t(2)) up to O(delta t(10)), to simuate the non-linear dynamics of the appropriate fields in comoving three-dimensional lattices. The code is parallelized with MPI, and uses a discrete Fourier Transform parallelized in multiple spatial dimensions, which makes it a very powerful code for probing physical problems with well-separated scales. Moreover, the code has been designed as a `platform' to implement any system of dynamical equations suitable for discretization on a lattice.

Abstract: We propose a minimal gauged U(1)(B-L) extension of the minimal supersymmetric Standard Model (MSSM) which resolves the cosmological moduli problem via thermal inflation, and realizes late-time Affleck-Dine leptogensis so as to generate the right amount of baryon asymmetry at the end of thermal inflation. The present relic density of dark matter can be explained by sneutrinos, MSSM neutralinos, axinos, or axions. Cosmic strings from U(1)(B-L) breaking are very thick, and so the expected stochastic gravitational wave background from cosmic string loops has a spectrum different from the one in the conventional Abelian-Higgs model, as would be distinguishable at least at LISA and DECIGO. The characteristic spectrum is due to a flat potential, and may be regarded as a hint of supersymmetry. Combined with the resolution of moduli problem, the expected signal of gravitational waves constrains the U(1)(B-L) breaking scale to be O(10(12-13)) GeV. Interestingly, our model provides a natural possibility for explaining the observed ultra-high-energy cosmic rays thanks to the fact that the core width of strings in our scenario is very large, allowing a large enhancement of particle emissions from the cusps of string loops. Condensation of LHu flat-direction inside of string cores arises inevitably and can also be the main source of the ultra-high-energy cosmic rays accompanied by ultra-high-energy lightest supersymmetric particles.