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 this paper, we explore in detail the cosmological implications of an abelian L – L gauge extension of the Standard Model featuring a light and weakly coupled Z. Such a scenario is motivated by the longstanding approximate to 4 sigma discrepancy between the measured and predicted values of the muon's anomalous magnetic moment, (g – 2), as well as the tension between late and early time determinations of the Hubble constant. If sufficiently light, the Z population will decay to neutrinos, increasing the overall energy density of radiation and altering the expansion history of the early universe. We identify two distinct regions of parameter space in this model in which the Hubble tension can be significantly relaxed. The first of these is the previously identified region in which a approximate to 10 – 20 MeV Z reaches equilibrium in the early universe and then decays, heating the neutrino population and delaying the process of neutrino decoupling. For a coupling of g (-) similar or equal to (3 – 8) x 10(-4), such a particle can also explain the observed (g – 2) anomaly. In the second region, the Z is very light (mZ approximate to 1eV to MeV) and very weakly coupled (g (-) approximate to 10(-13) to 10(-9)). In this case, the Z population is produced through freeze-in, and decays to neutrinos after neutrino decoupling. Across large regions of parameter space, we predict a contribution to the energy density of radiation that can appreciably relax the reported Hubble tension, N-eff similar or equal to 0.2.

Abstract: In view of the IceCube's 6-year high-energy starting events (HESE) sample, we revisit the possibility that the updated data may be better explained by a combination of neutrino fluxes from dark matter decay and an isotropic astrophysical power-law than purely by the latter. We find that the combined two-component flux qualitatively improves the fit to the observed data over a purely astrophysical one, and discuss how these updated fits compare against a similar analysis done with the 4-year HESE data. We also update fits involving dark matter decay via multiple channels, without any contribution from the astrophysical flux. We find that a DM-only explanation is not excluded by neutrino data alone. Finally, we also consider the possibility of a signal from dark matter annihilations and perform analogous analyses to the case of decays, commenting on its implications.

Abstract: We consider an extension of the Standard Model with two right-handed singlet fermions with mass at the electroweak scale that induce neutrino masses, plus a generic new physics sector at a higher scale Lambda. We focus on the effective operators of lowest dimension d = 5, which induce new production and decay modes for the singlet fermions. We assess the sensitivity of future Higgs Factories, such as FCC-ee, CLIC-380, ILC and CEPC, to the coefficients of these operators for various center of mass energies. We show that future lepton colliders can test the cut-off of the theory up to Lambda similar or equal to 500-1000 TeV, surpassing the reach of future indirect measurements of the Higgs and Z boson widths. We also comment on the possibility of determining the underlying model flavor structure should a New Physics signal be observed, and on the impact of higher dimensional d = 6 operators on the experimental signatures.

Abstract: After developing a general criterion for deciding which neutrino mass models belong to the category of inverse seesaw models, we apply it to obtain the Dirac analogue of the canonical Majorana inverse seesaw model. We then generalize the inverse seesaw model and obtain a class of inverse seesaw mechanisms both for Majorana and Dirac neutrinos. We further show that many of the models have double or multiple suppressions coming from tiny symmetry breaking “mu -parameters”. These models can be tested both in colliders and with the observation of lepton flavour violating processes.