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: In addition to spectacular signatures such as black hole superradiance and the rotation of CMB polarization, the plenitude of axions appearing in the string axiverse may have potentially dangerous implications. An example is the cosmological overproduction of relic axions and moduli by the misalignment mechanism, more pronounced in regions where the signals mentioned above may be observable, that is for large axion decay constant. In this work, we study the minimal requirements to soften this problem and show that the fundamental requirement is a long period of low-scale inflation. However, in this case, if the inflationary Hubble scale is lower than around O(100) eV, no relic DM axion is produced in the early Universe. Cosmological production of some axions may be activated, via the misalignment mechanism, if their potential minimum changes between inflation and today. As a particular example, we study in detail how the maximal-misalignment mechanism dilutes the effect of dangerous axions and allows the production of axion DM in a controlled way. In this case, the potential of the axion that realises the mechanism shifts by a factor increment theta = pi between the inflationary epoch and today, and the axion starts to oscillate from the top of its potential. We also show that axions with masses m(a) similar to O(1 – 100) H-0 realising the maximal-misalignment mechanism generically behave as dark energy with a decay constant that can take values well below the Planck scale, avoiding problems associated to super-Planckian scales. Finally, we briefly study the basic phenomenological implications of the mechanism and comment on the compatibility of this type of maximally-misaligned quintessence with the swampland criteria.

Abstract: Axions, if realized in nature, can be copiously produced in the early universe via thermal processes, contributing to the mass-energy density of thermal hot relics. In light of the most recent cosmological observations, we analyse two different thermal processes within a realistic mixed hot dark matter scenario which includes also massive neutrinos. Considering the axion-gluon thermalization channel, we derive our most constraining bounds on the hot relic masses m(a) < 7.46 eV and Sigma m(nu) < 0.114 eV both at 95 percent CL; while studying the axion-pion scattering, without assuming any specific model for the axion-pion interactions, and remaining in the range of validity of the chiral perturbation theory, our most constraining bounds are improved to m(a) < 0.91 eV and Sigma m(nu) < 0.105 eV, both at 95 percent CL. Interestingly, in both cases, the total neutrino mass lies very close to the inverted neutrino mass ordering prediction. If future terrestrial double beta decay and/or long-baseline neutrino experiments find that the nature mass ordering is the inverted one, this could rule out a wide region in the currently allowed thermal axion window. Our results therefore, strongly support multi messenger searches of axions and neutrino properties, together with joint analyses of their expected sensitivities.