Abstract: This paper presents the main features of a new and updated version of the program PArthENoPE, which the community has been using for many years for computing the abundances of light elements produced during Big Bang Nucleosynthesis. This is the third release of the PArthENoPE code, after the 2008 and the 2018 ones, and will be distributed from the code's website, http://parthenope.na.infn.it. Apart from minor changes, the main improvements in this new version include a revisited implementation of the nuclear rates for the most important reactions of deuterium destruction, H-2(p,gamma) He-3, H-2(d, n)He-3 and H-2(d, p)H-3, and a re-designed GUI, which extends the functionality of the previous one. The new GUI, in particular, supersedes the previous tools for running over grids of parameters with a better management of parallel runs, and it offers a brand-new set of functions for plotting the results. Program summary Program title: PArthENoPE 3.0 CPC Library link to program files: https://doi.org/10.17632/wygr7d8yt9.2 Developer's repository link: http://parthenope.na.infn.it Licensing provisions: GPLv3 Programming language: Fortran 77 and Python Nature of problem: Computation of yields of light elements synthesized in the primordial universe Solution method: Livermore Solver for Ordinary Differential Equations (LSODE) for stiff and nonstiff systems, Python GUI for running and plotting Journal reference of previous version: Comput. Phys. Commun. 233 (2018) 237-242 Does the new version supersede the previous version?: Yes Reasons for the new version: Update of the physics and improvements in the GUI Summary of revisions: Update of the physics implemented in the Fortran code and improvements in the GUI functionalities, in particular new plotting functions.

Abstract: One of the most criticized features of Bayesian statistics is the fact that credible intervals, especially when open likelihoods are involved, may strongly depend on the prior shape and range. Many analyses involving open likelihoods are affected by the eternal dilemma of choosing between linear and logarithmic prior, and in particular in the latter case the situation is worsened by the dependence on the prior range under consideration. In this letter, we revive a simple method to obtain constraints that depend neither on the prior shape nor range and, using the tools of Bayesian model comparison, extend it to overcome the possible dependence of the bounds on the choice of free parameters in the numerical analysis. An application to the case of cosmological bounds on the sum of the neutrino masses is discussed as an example.

Abstract: Various particle physics models suggest that, besides the (nearly) cold dark matter that accounts for current observations, additional but sub-dominant dark relics might exist. These could be warm, hot, or even contribute as dark radiation. We present here a comprehensive study of two-component dark matter scenarios, where the first component is assumed to be cold, and the second is a non-cold thermal relic. Considering the cases where the non-cold dark matter species could be either a fermion or a boson, we derive consistent upper limits on the non-cold dark relic energy density for a very large range of velocity dispersions, covering the entire range from dark radiation to cold dark matter. To this end, we employ the latest Planck Cosmic Microwave Background data, the recent BOSS DR11 and other Baryon Acoustic Oscillation measurements, and also constraints on the number of Milky Way satellites, the latter of which provides a measure of the suppression of the matter power spectrum at the smallest scales due to the free-streaming of the non-cold dark matter component. We present the results on the fraction f(ncdm) of non-cold dark matter with respect to the total dark matter for different ranges of the non-cold dark matter masses. We find that the 2 sigma limits for non-cold dark matter particles with masses in the range 1-10 keV are f(ncdm) <= 0.29 (0.23) for fermions (bosons), and for masses in the 10-100 keV range they are f(ncdm) <= 0.43 (0.45), respectively.

Abstract: We consider the phenomenological implications of the violation of the Pauli exclusion principle for neutrinos, focusing on cosmological observables such as the spectrum of Cosmic Microwave Background anisotropies, Baryon Acoustic Oscillations and the primordial abundances of light elements. Neutrinos that behave (at least partly) as bosonic particles have a modified equilibrium distribution function that implies a different influence on the evolution of the Universe that, in the case of massive neutrinos, can not be simply parametrized by a change in the effective number of neutrinos. Our results show that, despite the precision of the available cosmological data, only very weak bounds can be obtained on neutrino statistics, disfavouring a more bosonic behaviour at less than 2 sigma.

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.