Touch of Neutrinos on the Vacuum Metamorphosis: is the $H_0$ Solution Back?

With the entrance of cosmology in its new era of high precision experiments, low- and high-redshift observations set off tensions in the measurements of both the present-day expansion rate ($H_0$) and the clustering of matter ($S_8$). We provide a simultaneous explanation of these tensions using the Parker-Raval Vacuum Metamorphosis (VM) model with the neutrino sector extended beyond the three massless Standard Model flavours and the curvature of the universe considered as a model parameter. To estimate the effect on cosmological observables we implement various extensions of the VM model in the standard texttt{CosmoMC} pipeline and establish which regions of parameter space are empirically viable to resolve the $H_0$ and $S_8$ tensions. We find that the likelihood analyses of the physically motivated VM model, which has the same number of free parameters as in the spatially-flat $Lambda$CDM model, always gives $H_0$ in agreement with the local measurements (even when BAO or Pantheon data are included) at the price of much larger $chi^2$ than $Lambda$CDM. The inclusion of massive neutrinos and extra relativistic species quantified through two well known parameters $sum m_{ u}$ and $N_{rm eff}$, does not modify this result, and in some cases improves the goodness of the fit. In particular, for the original VM+$sum m_ u$+$N_{rm eff}$ and the Planck+BAO+Pantheon dataset combination, we find evidence for $sum m_{ u}=0.80^{+0.18}_{-0.22}~{rm eV}$ at more than $3sigma$, no indication for extra neutrino species, $H_0=71.0pm1.2$~km/s/Mpc in agreement with local measurements, and $S_8=0.755pm0.032$ that solves the tension with the weak lensing measurements. [Abridged]