Abstract: We present a full Bayesian analysis of the combination of current neutrino oscillation, neutrinoless double beta decay and Cosmic Microwave Background observations. Our major goal is to carefully investigate the possibility to single out one neutrino mass ordering, namely Normal Ordering or Inverted Ordering, with current data. Two possible parametrizations (three neutrino masses versus the lightest neutrino mass plus the two oscillation mass splittings) and priors (linear versus logarithmic) are exhaustively examined. We find that the preference for NO is only driven by neutrino oscillation data. Moreover, the values of the Bayes factor indicate that the evidence for NO is strong only when the scan is performed over the three neutrino masses with logarithmic priors; for every other combination of parameterization and prior, the preference for NO is only weak. As a by-product of our Bayesian analyses, we are able to (a) compare the Bayesian bounds on the neutrino mixing parameters to those obtained by means of frequentist approaches, finding a very good agreement; (b) determine that the lightest neutrino mass plus the two mass splittings parametrization, motivated by the physical observables, is strongly preferred over the three neutrino mass eigenstates scan and (c) find that logarithmic priors guarantee a weakly-to-moderately more efficient sampling of the parameter space. These results establish the optimal strategy to successfully explore the neutrino parameter space, based on the use of the oscillation mass splittings and a logarithmic prior on the lightest neutrino mass, when combining neutrino oscillation data with cosmology and neutrinoless double beta decay. We also show that the limits on the total neutrino mass Sigma m(nu) can change dramatically when moving from one prior to the other. These results have profound implications for future studies on the neutrino mass ordering, as they crucially state the need for self-consistent analyses which explore the best parametrization and priors, without combining results that involve different assumptions.

Abstract: In the framework of a 3+1 scheme with an additional inert state, we consider the thermalisation of sterile neutrinos in the early Universe taking into account the full 4 x 4 mixing matrix. The evolution of the neutrino energy distributions is found solving the momentum-dependent kinetic equations with full diagonal collision terms, as in previous analyses of flavour neutrino decoupling in the standard case. The degree of thermalisation of the sterile state is shown in terms of the effective number of neutrinos, N-eff, and its dependence on the three additional mixing angles (theta(14), theta(24), theta(34)) and on the squared mass difference Delta m(41)(2) is discussed. Our results are relevant for fixing the contribution of a fourth light neutrino species to the cosmological energy density, whose value is very well constrained by the final Planck analysis. For the preferred region of active-sterile mixing parameters from short-baseline neutrino experiments, we find that the fourth state is fully thermalised (N-eff similar or equal to 4).

Abstract: We show that for a relatively light majoron (<< 100 eV) non-thermal production from topological defects is an efficient production mechanism. Taking the type I seesaw as benchmark scheme, we estimate the primordial majoron abundance and determine the required parameter choices where it can account for the observed cosmological dark matter. The latter is consistent with the scale of unification. Possible direct detection of light majorons with future experiments such as PTOLEMY and the formation of boson stars from the majoron dark matter are also discussed.

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: Measuring temperature fluctuations in the 21 cm signal from the epoch of reionization and the cosmic dawn is one of the most promising ways to study the universe at high redshifts. Unfortunately, the 21 cm signal is affected by both cosmology and astrophysics processes in a nontrivial manner. We run a suite of 1000 numerical simulations with different values of the main astrophysical parameters. From these simulations we produce tens of thousands of 21 cm maps at redshifts 10 <= z <= 20. We train a convolutional neural network to remove the effects of astrophysics from the 21 cm maps and output maps of the underlying matter field. We show that our model is able to generate 2D matter fields not only that resemble the true ones visually but whose statistical properties agree with the true ones within a few percent down to scales 2 Mpc(-1). We demonstrate that our neural network retains astrophysical information that can be used to constrain the value of the astrophysical parameters. Finally, we use saliency maps to try to understand which features of the 21 cm maps the network is using in order to determine the value of the astrophysical parameters.