Neutrino masses and Hubble tension via a Majoron in MFV

The recent tension between local and early measurements of the Hubble constant can be explained in a particle physics context. A mechanism is presented where this tension is alleviated due to the presence of a Majoron, arising from the spontaneous breaking of Lepton Number. The lightness of the active neutrinos is consistently explained. Moreover, this mechanism is shown to be embeddable in the minimal (Lepton) flavour violating context, providing a correct description of fermion masses and mixings, and protecting the flavour sector from large deviations from the Standard Model predictions. A QCD axion is also present to solve the Strong CP problem. The Lepton Number and the Peccei–Quinn symmetries naturally arise in the minimal (Lepton) flavour violating setup and their spontaneous breaking is due to the presence of two extra scalar singlets. The Majoron phenomenology is also studied in detail. Decays of the heavy neutrinos and the invisible Higgs decay provide the strongest constraints in the model parameter space.


Introduction
There is nowadays a considerable tension between late-time, local probes of the present rate of expansion of the Universe, that is the Hubble constant H 0 , and its value inferred through the standard cosmological model ΛCDM from early Universe observations.Local measurements, from type-Ia supernovae and strong lensing, tend to cluster at similar values of H 0 , significantly larger than those preferred by cosmic microwave background (CMB) and baryon acoustic oscillations probes.The strongest tension, estimated at the level of 4 − 6 σ [1,2] depending on the assumptions performed, is between the Planck inferred measure from the CMB spectrum [3] and the one obtained by the SH 0 ES collaboration [4] from supernovae measurements.
Although the solution to this discrepancy might be related to systematics in the measurements (notably the callibration of the supernovae distances) or, more interestingly, point to a modification of the cosmological model, it may instead be provided by particle physics, as already discussed in the literature (see for example [5][6][7][8][9][10][11][12][13][14][15][16][17][18]).In particular, Ref. [14] suggests that a Majoron that couples to light neutrinos could reduce the tension in the determination of H 0 .It is then interesting to investigate whether this setup is compatible with possible explanations of other open problems in the Standard Model of particle physics (SM): the focus of this paper is to study the compatibility with Type-I Seesaw mechanism to explain the lightness of the active neutrinos, with specific flavour symmetries to describe the flavour puzzle and with the presence of an axion to solve the Strong CP problem.
The Majoron, called ω hereafter, is the Nambu-Goldstone boson (NGB) associated to the spontaneous breaking of lepton number (LN) [19][20][21][22], which is only accidental within the SM and breaks down at the quantum level.It naturally arises in the context of the Type-I Seesaw mechanism, where the Majorana mass term, instead of being a simple bilinear of the right-handed (RH) neutrino fields N R , is a Yukawa-like term that couples N R to a scalar field that carries a LN charge, labelled as χ in the following.Once this scalar field develops a vacuum expectation value (VEV), then LN is spontaneously broken, a Majorana neutrino mass term is generated and the Majoron appears as a physical degree of freedom of the spectrum.
If a Dirac term that mixes N R and the left-handed (LH) lepton doublets ℓ L is also present in the Lagrangian, small masses for active neutrinos are generated at low energies, according to the Type-I Seesaw mechanism.
At low-energies, the Majoron ω acquires a coupling with ν L , labelled as λ ωνν .For Majoron masses m ω ∈ [0.1, 1] eV (1.1) and λ ωνν in the range λ ωνν ∈ [5 × 10 −14 , 10 −12 ] , the tension on the Hubble constant is reduced [14].Indeed, for such small Majoronneutrino mixings and Majoron masses, Majorons only partially thermalise after Big Bang Nucleosynthesis (BBN) or never thermalise [23], enhancing the effective number of neutrino species N eff by at least 0.03 and at most 0.11, values that may be tested with CMB-S4 experiments [24].Moreover, a non-vanishing λ ωνν would reduce neutrino free-streaming, modifying the neutrino anisotropic stress energy tensor [25].This has an impact in the CMB that results in modifying the posterior for the Hubble constant: as shown in Ref. [14], the inclusion of Majoron-neutrino interactions slightly shifts the central value of H 0 , but largely broadens its profile reducing the H 0 tension to 2.5σ.For larger couplings, λ ωνν > 10 −12 , these effects are too large and excluded by the same Planck data.Interestingly, Ref. [14] found that the best χ 2 in a Markov Chain Monte Carlo corresponds to Majoron mass and coupling as in Eqs.(1.1) and (1.2) and ∆N eff = 0.52 ± 0.19.The uncertainty in the last observable is very large and ∆N eff = 0 is compatible within 3σ.However, values close to the central one can be achieved if a thermal population of Majorons is produced in the early Universe and is not diluted during inflation.This may occur if the reheating temperature is larger than the RH neutrino masses [14].Alternatively, other relativistic species, such as axions [26][27][28][29], may contribute to ∆N eff and their presence may justify such a large value.
The existence of both Majorana and Dirac terms, the achievement of the correct scale for the active neutrino masses and at the same time the alleviation of the Hubble tension via the Majoron strongly depend on the LN charge assignments of ℓ L , N R and χ.In particular, fixing the LN of ℓ L to unity, then the LN of the RH neutrinos needs to have opposite sign with respect to the LN of the scalar field χ.This model presents interesting phenomenological features.On one hand, the heavy neutrinos are expected to be relatively light, with masses at the MeV or GeV scales, opening up the possibility to be studied at colliders.Moreover, the presence of the Majoron may also have other consequences distinct from the Hubble tension.In particular, its couplings to photons and electrons are constrained by CAST and Red Giant observations, respectively, while, due to its coupling to the Higgs, the Majoron may contribute to the invisible Higgs decay, strongly constrained by colliders.
Sect. 2 illustrates the mechanism to produce a Majoron that alleviates the H 0 tension together with a correct scale for the active neutrino masses.In Sect.3, this mechanism is introduced in a setup that correctly describes the flavour puzzle of the SM and at the same time produces a QCD axion that solves the Strong CP problem.Sect.4 gathers the phenomenological signatures of this model and Sect 5 contains the final remarks.

The Majoron Mechanism
To produce a Majoron and explain the lightness of the active neutrinos, one can consider a Type-I Seesaw mechanism where the Majorana mass is dynamically generated by the spontaneous breaking of LN.The SM spectrum is extended by three RH neutrinos1 and a singlet scalar field χ that only transforms under LN.The LN charge assignments can be read in Tab. 1, where ℓ L , N R and χ have already been defined, and e R refers to the RH charged leptons.Notice that L χ and L N are integer numbers and are not completely free, but must obey a series of constraints that will be made explicit in the following.
Table 1: LN assignments.Fields that are not listed here do not transform under LN.
The most general effective Lagrangian in the neutrino sector invariant under LN is the following 2 : where H is the SM Higgs doublet, H = iσ 2 H * , Λ χ is the cut-off scale up to which the effective operator approach holds, and Y ν is a dimensionless and complex matrix, while Y N is dimensionless, complex and symmetric.A first condition on L N,χ arises from requiring that all the terms are local: In the LN broken phase, the field χ can be parametrised as where the angular part ω is the NGB identified as a Majoron, σ is the radial component and v χ is its VEV.Notice that the scale appearing in the denominator of the exponent is also v χ in order to obtain canonically normalised kinetic terms for the Majoron.A useful notation that will be employed in the following is the ratio of the χ VEV and the cut-off scale: This parameter ε χ is expected to be smaller than 1 in order to guarantee a good expansion in terms of 1/Λ χ .Consequently, the χ VEV, which represents the overall scale of the LN breaking, is expected to be smaller that the scale Λ χ , where New Physics should be present and is responsible for generating the expression in Eq. (2.2).
2 Other terms can be added to this Lagrangian inserting χ † instead of χ.However, once the terms in Eq. (2.2) are local, then their siblings with χ † would not be local and therefore cannot be added to the Lagrangian.The only exception is the term that however only provides a suppressed correction with respect to the Majorana term written in Eq. (2.2) and for this reason it can be neglected.
Once the electroweak symmetry is also spontaneously broken, i.e. after that the SM Higgs develops its VEV that in the unitary gauge reads where h is the physical Higgs and v ≃ 246 GeV, masses for the active neutrinos are generated according to the Type-I Seesaw mechanism: (2.7) In the basis where the charged lepton mass matrix is already diagonal, the neutrino mass matrix can be diagonalised by the PMNS matrix U : The overall scale for the active neutrino masses can be written in terms of the parameter ε χ , the ratio of the VEVs and the product where ∆m 2 atm = 2.514 +0.028 −0.027 × 10 −3 eV 2 for the Normal Ordering (NO) of the neutrino mass spectrum and ∆m 2 atm = −2.497± 0.028 × 10 −3 eV 2 for the Inverted Ordering (IO) [30].The heavy neutrinos, that mostly coincide with the RH neutrinos, have a mass matrix that in first approximation can be directly read from the Lagrangian in Eq. (2.2), (2.10) The overall scale for the heavy neutrinos must be larger than the overall scale of the active neutrinos, in order for the Seesaw approximation to hold: On the other hand, the electroweak and LN breakings give rise to the Majoron Lagrangian that can be written as follows: where the m 2 ω term parametrises the Majoron mass introduced here as an explicit breaking of its corresponding shift symmetry.At low energy, a direct coupling of the Majoron to the active neutrinos emerges after performing the same transformations that give rise to the mass matrix in Eq. (2.7): From the results in Ref. [14], shown in Eq. (1.2), it is then possible to infer a bound on the product L χ v χ , once taking ∆m 2 atm as the overall scale for the neutrino masses: Adopting this relation and substituting it within the expressions in Eqs.(2.9) and (2.11), new conditions can be found: (2.16) The following choice of charge assignments leads to a completely renormalisable Lagrangian and thus deserves special mention: CASE R: such that the powers of the ratio χ/Λ χ in Eq. ( 2 Eq. (2.16) The first expression fixes a range of values for v χ .While the third expression represents a lower bound on the overall scale of Y N , the second one implies that the product should be tuned to a very small value in order to reproduce the lightness of the active neutrinos.
For values of L N,χ different from the previous ones, the Lagrangian is necessarily nonrenormalisable.An interesting question is whether the extremely small values of the product Y ν Y −1 N Y T ν can be avoided exploiting the suppression in ϵ χ from the new physics scale Λ χ , similarly to the Froggatt-Nielsen approach to the flavour puzzle [31].Considering first the case in which L N,χ > 0, then the only available possibilities are CASE NR1: The corresponding values for v χ , ε χ , the overall scale for the heavy neutrinos ⟨M N ⟩ and the cut-off scale Λ χ are reported in Tab. 2.
Table 2: Parameter ranges in the two phenomenologically interesting scenarios.
From Eq. ( 2.15), it can be seen that ϵ χ gets smaller for larger values of L χ (unless tuning the product Y ν Y −1 N Y T ν as in CASE R discussed above, a possibility to be avoided in the present discussion): although this is not a problem by itself, it hardens the constraint in Eq. (2.16).It follows that for larger values of L N and L χ satisfying to the locality conditions in Eq. ( 2.3), then the overall scale of the heavy neutrino masses would be as small as the one of active neutrinos and therefore the expansion in the Type-I Seesaw mechanism would break down.
In the case when L N > 0 and L χ < 0, it is possible to obtain the same results listed above substituting χ by χ † in the Lagrangian in Eq. (2.2): in this case, the signs in the denominators of the exponents would be flipped, compensating the negative sign of L χ .The opposite case, L N < 0 and L χ > 0, is not allowed by the locality conditions.
For L N,χ < 0, besides the possibility of CASE R, only another choice is allowed by the locality conditions: (L N = −1, L χ = −1).However, this case would require ϵ χ ≫ 1, leading to an even more extreme fine-tuning than in CASE R without the appeal of renormalizability.
The condition in Eq. (2.14), corresponding to Eq. (1.2), is only one of the ingredients necessary to lower the H 0 tension.A second relevant requirement is Eq.(1.1), regarding the Majoron mass.For the sake of simplicity, it has been introduced directly in the Majoron Lagrangian in Eq. (2.12) as en effective term.Its origin has been widely discussed in the literature and constitutes in itself an interesting research topic.Any violation of the global LN symmetry would induce a mass for the Majoron.An obvious example are gravitational effects, which are expected to break all accidental global symmetries.Estimations of their size from non-perturbative arguments via wormhole effects [32] fall too short of their required value.Conversely, their size from Planck-suppressed effective operators [33] would point to too large a mass, although several possibilities have been discussed that would prevent the lower dimension operators from being generated [34][35][36][37][38][39][40][41].These options were originally introduced as a solution to the axion quality problem.A simpler possibility, given the singlet nature of the N R , is an explicit breaking of LN via a Majorana mass term at the Lagrangian level.The Majoron would thus develop a mass slightly below this breaking scale from its coupling to the N R through self-energy diagrams.In this work we will remain agnostic to the origin of the Majoron mass, and a value consistent with Eq. (1.1) will be assumed.
The mechanism illustrated in this section allows to soften the H 0 tension explaining at the same time the lightness of the active neutrinos.In the next section, this mechanism will be embedded into a specific flavour model that allows to account for the Flavour Puzzle, without violating any bounds from flavour observables, and also contains a QCD axion that solves the Strong CP Problem.

The Majoron and Axion from an MFV Setup
Flavour models aim at describing the heterogeneity of the fermion masses and mixings through some underlying argument, such as a symmetry principle.The simplest proposal is known as the Froggatt-Nielsen model [31] and consists in introducing a global Abelian symmetry to describe the flavour structure in the quark sector.Almost 25 years later, after the more precise measurements of neutrino oscillations and the introduction of the so-called Tri-Bimaximal mixing texture [42,43], whose main feature is a vanishing reactor mixing angle in the PMNS matrix, discrete symmetries were considered as an attractive approach to reproduce the observed pattern of all fermion masses and mixings [44][45][46][47][48][49][50][51].However, after the discovery, in 2011, of a relatively large reactor mixing angle [52-56], models based on discrete symmetries underwent a deep rethinking and other options also became popular.A few examples are Froggatt-Nielsen inspired models based on a simple U (1) [57][58][59] and models based on continuous non-Abelian symmetries.The latter include Minimal Flavour Violation (MFV) [60] and its leptonic versions [61][62][63], which are based on a U (3) symmetry, smaller symmetries like U (2) [64,65], or an intermediate approach [66].
The focus here will be the MFV setup, that will be shown to naturally suggest the presence of the Majoron and of a QCD axion.The key concept behind MFV is to require that any flavour and CP violation in physics Beyond the SM has the same origin as the one in the SM [67].This idea has been formulated in terms of the symmetries of the kinetic terms in Ref. [60].These symmetries are then broken by a series of fields that are also employed to describe fermion masses and mixings, as well as the suppressions associated to any non-renormalisable flavour violating operator.Considering the SM spectrum augmented by three RH neutrinos, the flavour symmetry of the corresponding kinetic terms is a product of a U (3) term for each fermion species, (3.1) The non-Abelian terms are responsible for the intergenerational fermion mass hierarchies and of the mixing matrices, while the Abelian terms are associated to the hierarchies among the masses of the different third family fermions.The latter can be written as follows where q L refers to the LH quark doublets, u R to the RH up-type quarks and d R to the RH down-type quarks, while ℓ L , e R and N R to the leptons as already introduced in the previous section.It is possible to rearrange this product, provided that the new combinations are still linearly independent, identifying among them baryon number, LN, weak hypercharge and the Peccei-Quinn (PQ) [68] symmetry: In the model described in this section, fermion charges under baryon number and hypercharge are assigned as in the SM, while the LN charges are given in Tab. 1.Moreover, PQ charges are chosen so as to explain the suppression of the bottom and τ masses with respect to the top mass: only d R and e R transform under PQ, with a charge equal to 3.
In an analogous way to LN, the PQ symmetry is formally exact at the Lagrangian level after introducing a second scalar field Φ that transforms non-trivially only under PQ with a charge −1: where ρ is the radial component, a is the angular one and v Φ is its VEV.Once this scalar field develops a VEV and electroweak symmetry gets broken, masses for bottom and τ are generated, suppressed with respect to the top one.As a byproduct of this mechanism, the axion a is originated [69].Finally, the last two symmetries in Eq. (3.3) do not play any role in this model and are explicitly broken.The Lagrangian invariant under the aforementioned symmetries is the following: where Λ Φ stands for the cut-off scale associated to the scalar field Φ.
According to the MFV approach, Y u,d,e,ν,N are not simple matrices, but are promoted to be spurion fields that do transform under the non-Abelian part of the flavour symmetry group G F .These fields can be thought of as dimensionless scalar fields that do not have kinetic terms, but acquire background values that play the role of VEVs.In the MFV approach, to reproduce correctly quark masses and mixings and charged lepton masses, the background values of Y u,d,e should read , where V is the CKM mixing matrix and c i are order 1 parameters.The possible origin of these values is under study [70][71][72][73].
While these matrices present hierarchies among their entries, the ratios m b /m t and m τ /m t are still not described within this approach, but can be explained by the spontaneous breaking of the PQ symmetry.Indeed, the overall scale of down-type quarks and charged leptons is multiplied by the cubic power of the ratio among the Φ VEV and the cut-off Λ Φ .Assuming that this ratio is given by the m b /m t and m τ /m t ratios are also correctly described.Any non-renormalisable operator that describes flavour-violating processes should be invariant under the flavour symmetry.This is accomplished by inserting proper powers of the spurions: once they acquire their background values, the non-renormalisable operator under consideration gets suppressed.The main consequence is that the scale that can be studied considering flavour observables is at the level of 1 − 10 TeV [60,[74][75][76][77][78][79][80][81][82][83][84], instead of 100 TeV in a generic case [85], opening up the possibility of discovering New Physics (NP) at colliders.
For the neutrino sector, the discussion is slightly more complicated.Indeed, there are two spurions in the neutrino sector, Y ν and Y N , and both enter in the definition of the active neutrino masses, see Eq. (2.7).It follows that it is not possible to identify univocally either Y ν or Y N in terms of neutrino masses and the PMNS matrix entries.Thus, the suppression in the non-renormalisable flavour violating operators cannot be directly linked to neutrino masses and mixings, losing the predictivity that characterizes the MFV approach in the quark sector.The solutions that have been proposed are to assume Y N ∝ 1 [61,62] or to consider Y ν as a unitary matrix [63].[61,62].
Under the assumption that the three RH neutrinos are degenerate in mass, i.e.
Y N ∝ 1, then the non-Abelian symmetry associated to the RH neutrinos, SU (3) N R , is broken down to SO(3) N R .The additional assumption of no CP violation in the lepton sector is meant to force Y ν to be real 4 .With these simplifications, the expression for the active neutrino mass in Eq. (2.7) simplifies to and all flavour changing effects involving leptons can be written in terms of Y ν Y T ν and Y e .These are the only relevant combinations entering any d = 6 operator involving lepton fields and describing flavour violating effects.It follows that any flavour changing process can be predicted in terms of lepton masses and mixings.
Diagonalising m ν corresponds to diagonalising the product ⟨Y ν ⟩⟨Y ν ⟩ T and, given the fact that the the lepton mixing angles are relatively large, then no hierarchies should be expected among the entries of ⟨Y ν ⟩⟨Y ν ⟩ T , contrary to what happens in the quark case.Note that some setups, such as the so-called sequential dominance scenarios, obtain large mixing angles even if there exists a strong hierarchy among the Yukawa couplings [86].However, this possibility is disfavoured by the general philosophy of MLFV.In the same spirit, the overall scale of this product is of O(1), as any explanation of the neutrino masses should reside in the model itself, and not be due to any fine-tuning.[63].

II):
Assuming that the three RH neutrinos transform as a triplet under the same symmetry group of the lepton doublets, then Schur's Lemma guarantees that Y ν transforms as a singlet of the symmetry group.Then, Y ν is a unitary matrix [87,88], which can always be rotated to the identity matrix by a suitable unitary transformation acting only on the RH neutrinos.
The only meaningful quantities in this context are Y e and Y N , so neutrino masses and lepton mixings are encoded uniquely into Y N , Moreover, all flavour changing effects involving leptons can be written only in terms of Y e and Y N , and therefore any flavour changing process can be predicted in terms of lepton masses and mixings.
As for the previous case, the diagonalisation of the active neutrino mass coincides with the diagonalisation of ⟨Y N ⟩ −1 , that therefore does not present any strong hierarchy among its entries.Moreover, its overall scale should be O(1) according to the MLFV construction approach.
In both cases, the constraints on NP considering the present available data on flavour changing processes in the lepton sector are as low as a few TeV [61][62][63][89][90][91][92].
Once the PQ symmetry is spontaneously broken, the axion arises as its NGB and its Lagrangian can be written as with G a µν ≡ 1 2 ϵ µνρσ G aρσ and ϵ µνρσ the totally antisymmetric tensor such that ϵ 1230 = 1.The last term is the well-known QCD θ-term, allowed by the QCD Lagrangian, which constitutes a source of CP violation.The θ QCD parameter contributes to the neutron electric dipole moment [93,94] and can thus be experimentally constrained [95] θ QCD ≲ 10 −10 . (3.12) The presence of the axion naturally explains why θ QCD is so small, providing a solution to the Strong CP problem: couplings of the axion to the gauge fields, and in particular to gluons, are generated at the quantum level and the θ QCD parameter can be reabsorbed by an axion field redefinition [68,96,97].The effective axion potential (see for example Ref. [98]) predicts a vanishing VEV for the axion, that finally solves the Strong CP problem.Moreover, the same potential provides the axion with a mass, being f a the axion scale, that is connected to the Φ VEV in this model by the relation For QCD axions, the most stringent constraint on f a comes from the axion couplings to photons [99][100][101] and to electrons [102][103][104], which push its value to be larger than f a ≳ 1.2 × 10 7 GeV and f a ≳ 8 × 10 8 GeV, respectively.
Summarising, the Majoron together with axion constitute the natural Abelian completion of MFV scenarios.The Majoron does not affect (at tree level) the physics associated to the axion and the quark and charged lepton flavour physics.Thus, this model, besides describing fermion masses and mixings and solving the Strong CP problem, is able to alleviate the Hubble tension with the only inclusion of three RH neutrinos and two extra singlet scalars.In particular, as no fine-tuning is allowed within this approach on ⟨Y ν ⟩ or ⟨Y N ⟩, then only CASE NR1 and CASE NR2 are viable in the MFV framework.In the following section, the analysis of this model will be completed with the study of its phenomenological signatures.

Phenomenological Signatures
The only tree-level coupling of the Majoron is to neutrinos.However, at quantum level, couplings to gauge bosons, other SM fermions and the Higgs are originated.

Coupling to photons
The searches for very light pseudoscalars, usually addressed to axions, can also apply to Majorons.In the range of masses in Eq. (1.1), the strongest constraints on the effective coupling to photons are set by CAST [101], which establishes the upper bound where λ ωγγ is defined as with F µν ≡ 1 2 ϵ µνρσ F ρσ .As the Majoron does not couple at tree-level to charged particles, then the process ω → γγ occurs only at two loops.Ref. [105] provides an estimate for its decay width: under the assumption m ω ≪ m e , where α ≡ e 2 /4π and with m D the Dirac neutrino mass matrix.
Computing the same process by means of the effective couplings in Eq. (4.2), it is then possible to match the two expressions for the ω → γγ decay width providing the expression for the λ ωγγ coupling: Tab. 3 shows the numerical estimations for the predicted values of the Majoron coupling to photons, assuming that the trace gives an O(1) number, as already discussed: the experimental bound is still far from the theoretical prediction.Table 3: Predictions for the Majoron effective couplings to electrons, photons and neutrinos, for the window of the parameter space where Hubble tension is alleviated.In the last line, the corresponding experimental upper bounds.

Coupling to electrons.
Astrophysical measurements can also constrain Majoron couplings.Ref. [102] provides an upper bound on the Majoron effective coupling to electrons based on observations on Red Giants: The decay width of the Majoron to two electrons reads [105] Γ where with [. ..] 11 standing for the (1, 1) entry of the matrix in the brackets.Substituting explicitly the expression for m D , the coupling becomes Assuming as before that the elements of the product Y ν Y † ν are O(1) numbers, also the Majoron-electron coupling is predicted to be much smaller than the corresponding experimental bound, as shown in Tab. 3.

Coupling to neutrinos. Majoron emission in 0νββ decays.
The tree level couplings of the Majoron to neutrinos does not have an impact only on cosmology, but may be relevant for low-energy terrestrial experiments.In particular, searches for neutrinoless-double-beta decay could also be sensitive to processes in which Majorons are produced.
Current measurements set a lower bound on the half-life of the neutrinoless-double-beta decay of the order of 10 26 years [106,107].However, these cannot be directly employed to constrain a process in which a Majoron is also produced, due to the different energy distribution of the emitted electrons.The NEMO-3 collaboration performed a dedicated search [108], setting a lower bound of the order of 10 22 years on Majoron emission in 0νββ decay.This corresponds to a limit on the Majoron-neutrino coupling that reads where λ ωνν is defined in Eq. (2.13).Note that the exact constraint depends on the choice of nuclear matrix elements.The predicted value of the Majoron-neutrino coupling can be read out in Tab. 3 and it is much smaller than the corresponding experimental value and the bound from Planck [14] (λ ωνν < O(10 −12 )).
Coupling with the SM Higgs.Higgs invisible decay.
Besides the Majoron, also the radial component of χ is present in the spectrum and does play a role modifying the Higgs scalar potential.Indeed, the most general scalar potential containing H and χ can be written as The minimisation of such potential leads to VEVs for the two fields that read The parameters of this scalar potential need to be such that v takes the electroweak value and v χ acquires the values in Tab. 2. Due to the mixed quartic term, the two physical scalar fields h and σ mix in the broken phase and the mass matrix describing this mixing is given by The two eigenvalues that arise after diagonalising this mass matrix are the following where tan 2ϑ = (4.16) The lightest mass in Eq. (4.15) corresponds to the eigenstate mainly aligned with the SM Higgs, while the heaviest state is mainly composed of the radial component of χ.Its mass can be as small as a few hundred GeV or much larger than the TeV.From the relation between the mixed quartic coupling g and the physical parameters, it is possible to straightforwardly study the dependence of M σ with the other model parameters.
The mixing of the Higgs with the new scalar also yields a coupling of the former to two Majorons.Expanding the kinetic term of χ yields a σωω term, that in turn produces a coupling hωω after an insertion of the scalar mixing.Thus, a new decay channel for the Higgs opens up, contributing to its invisible width.The rate of this process is given by These effects on Higgs phenomenology have an impact on the signal strength µ h , defined as the ratio of observed Higgs events with respect to the SM expectation.This quantity has been measured to be in perfect agreement with the SM by both ATLAS and CMS [109,110].
In our setup, both the production cross section and the visible decay rates are modified.The former is suppressed by a factor cos 2 ϑ, while the latter diminishes due to the appearance of the invisible channel.Explicitly, where Γ h is the total decay width of the Higgs boson.The combination of the ATLAS and CMS measurements yields a lower limit of µ h > 0.94 [111] at the 95% CL, which translates into an upper bound on sin ϑ that is proportional to v χ .These constraints can be portrayed in the parameter space composed of the mass of σ and the quartic coupling g, as seen in Fig. 2, in which the green areas are ruled out.
� χ = ��� ���  Although we choose a particular value for v χ , we find that the bounds are basically independent on this quantity.Clearly, M σ can reach very large values without requiring any fine-tuning on g; conversely, σ can only be light if g is close to zero.Note that the upper bound on the Higgs invisible branching ratio (that currently stands at a 13% [112]) constitutes an independent source of limits, which exhibit a very similar dependence with the parameters of the model.However, these bounds are always looser than those given by signal strengths, so we do not include them in our results.

Heavy Neutrinos
In both cases discussed in Sect.2, the heavy neutrino masses lie in ranges that may lead to detection in various experimental facilities.Neutrinos with masses ranging from tens to hundreds of MeV can be probed for and be potentially detected at beam dump or even near detectors of neutrino oscillation experiments [113][114][115][116][117][118][119][120], such as DUNE or SHiP, whereas those with masses in the range of tens to hundreds of GeV have interesting prospects of being produced at the LHC or future colliders [114,[121][122][123][124][125][126][127].
On the other hand, given their extremely small couplings, the heavy neutrinos produced in the early Universe would not be Boltzmann suppressed when decoupled from the thermal bath, leading to an unacceptably large contribution to the relativistic degrees of freedom of the Universe after their subsequent decay [128][129][130][131][132]. If their decay takes place before the onset of BBN, the decay products would quickly thermalise and BBN would then proceed as in the standard ΛCDM scenario.However, if the decay of the heavy neutrinos happens after BBN and neutrino decoupling, their contribution to the effective number of neutrinos would be too high and ruled out.If the decay takes place during BBN, the decay products could also alter the production of primordial helium and strong constraints also apply [133][134][135][136].This would be the situation of CASE NR1, for which the neutrino masses and mixings predict decay rates comparable to or larger than the onset of BBN.Conversely, the larger masses that characterize CASE NR2 lead to decays faster than BBN, eluding these cosmological constraints.Hence, BBN and CMB observations disfavor CASE NR1 unless the heavy neutrino decay is faster than BBN in some part of the parameter space or some other modification of the standard ΛCDM scenario is considered.Indeed, if the heavy neutrinos decay after BBN, for heavy neutrino masses in the range [3.5, 200] MeV, then the bound on the mixing is [132]: much smaller than the expected value that can be read in Tab. 4.

Conclusions
The Hubble tension may well be due to systematic uncertainties but it might point to a deviation from the ΛCDM model or an extension of the Standard Model of particle physics.In particular, the presence of a Majoron with couplings to the active neutrinos represents an interesting avenue to alleviate this tension.This paper is focussed on a mechanism where Table 4: Expectation for the heavy neutrino mass and mixing between heavy and active neutrinos.
neutrino masses and the presence of a Majoron are explained by the spontaneous breaking of Lepton Number.Besides the introduction of three RH neutrinos and a new complex scalar singlet, three sets of Lepton Number charge assignments are identified that correctly explain the lightness of the active neutrinos and provide the Majoron with couplings that lower the Hubble tension from 4.4σ down to 2.5σ.This mechanism is completely general and could be embedded in different types of flavour models.A compelling case presented here is the embedding in the Minimal (Lepton) Flavour Violating setup: Lepton Number is part of the complete flavour symmetry and therefore naturally arises in this context.The non-Abelian factors are mainly responsible of the intergenerational mass hierarchies and of the mixings, while Lepton Number, besides being associated to the Majoron, is involved in the explanation of the smallness of the active neutrino masses with respect to the top mass.The ratios between the bottom and τ masses can be explained via the spontaneous breaking of the PQ symmetry, which completes the non-Abelian flavour symmetry together with Lepton Number.A QCD axion arises as a byproduct of the PQ symmetry breaking, allowing to solve the Strong CP problem within the same framework.
The Majoron only couples at tree level to neutrinos.The neutrinoless-double-beta decay represents a natural observable where to look for it, but the model prediction turns out to be much smaller than the current experimental sensitivity and much stronger constraints can be derived from its impact in neutrino free streaming and the effective number of neutrino species from CMB observations.Effective couplings to charged fermions and photons arise at one and two loops, respectively.For this reason, the model is not significantly constraint from CAST and Red Giant observations that otherwise provide stringent bounds for light pseudoscalars.
Due to the mixing between the physical Higgs and the radial degree of freedom of the scalar field producing the Majoron, the latter does couple to the Higgs.This has an impact on the invisible Higgs decay, that indeed provides a strong bound on the VEV of the new scalar, that is, on the Lepton Number breaking scale.However, a natural choice in the model is to assume the radial component to be sufficiently heavy to be integrated out, and correspondingly its mixing with the Higgs gets suppressed, relaxing the bound on the VEV.
Finally, heavy neutrinos are expected to be much lighter than in the traditional type-I Seesaw scenario, in the MeV or GeV ranges depending on the specific case considered.In these mass ranges the decay products of the neutrinos could significantly affect the evolution of the early Universe.Indeed, for masses in the MeV range, decays during Big Bang Nucleosynthesis could alter the production of primordial helium while later decays would imply too large an injection of relativistic species.Conversely, for neutrino masses above the GeV, their decays are typically faster than Big Bang Nucleosynthesis and the decay products thermalise without altering the rest of the thermal history of the Universe.This option is thus preferred with respect to the lighter masses.
The result is a model where fermion masses and mixings can be correctly described, protecting the flavour sector from large deviations from the standard predictions with a new physics scale of the order of a few TeV, solving the Strong CP problem and softening at the same time the Hubble tension.This model could be tested indirectly at colliders looking at the invisible Higgs decay and searching for relatively light heavy neutrinos or for the radial component of the scalar field that generates the Majoron.

Figure 1 :
Figure 1: Feynman Diagram for the neutrinoless-double-beta decay with the emission of a Majoron.

Figure 2 :
Figure 2: Parameter space of the mass of the radial component of the new scalar as a function of its quartic coupling with the Higgs doublet.The green region is excluded by bounds on Higgs signal strength.We choose v χ = 100 GeV, although the bounds are basically insensitive to this quantity.