Electroweak vacuum stability and the seesaw mechanism revisited

We study the electroweak vacuum stability in Type I seesaw models for three generations of neutrinos in scenarios where the right-handed neutrinos have explicit bare mass terms in the Lagrangian and where these are dynamically generated through the mechanism of spontaneous symmetry breaking. To best highlight the difference of the two cases we concentrate on the absolute stability of the scalar potential. We observe that for the first scenario, the scale at which the scalar potential becomes unstable is lower from that within the standard model. In addition the Yukawa couplings Yν\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf {Y}_\nu $$\end{document} are constrained such that Tr[Yν†Yν]≲10-3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {Tr}{[\mathbf {Y}^{\dagger }_\nu \mathbf {Y}_{\nu }}] \lesssim 10^{-3}$$\end{document}. In the second scenario the electroweak stability can be improved in a large region of parameter space. However, we found that the scalar used to break the lepton number symmetry cannot be too light and have a large coupling to right-handed neutrinos in order for the seesaw mechanism to be a valid mechanism for neutrino mass generation. In this case we have Tr[Yν†Yν]≲0.01\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {Tr}[\mathbf {Y}^\dagger _{\nu } \mathbf {Y}_\nu ]\lesssim 0.01$$\end{document}.


Introduction
The discovery of the Higgs boson at the LHC completes the standard model (SM) of particle physics. The measured value of the Higgs boson mass, m h = 125.09 ± 0.21 stat. ± 0.11 sys. as given by the combined results of both ATLAS and CMS experiments [1][2][3] implies that the quartic coupling, λ H , of the SM scalar potential, V (H ), is relatively small. Specifically, the scalar potential is defined by Higgs field and v = 246 GeV. When it is used for the boundary value of the renormalization group (RG) running of λ H , it results in a negative value of λ H below the Planck scale due to the large top quark Yukawa coupling [4][5][6][7]. In fact, the latest NNLO study of the vacuum stability requires that a e-mail: misery@triumf.ca b e-mail: apuente@physics.carleton.ca m h > 129.4 ± 1.8 GeV [5] for absolute stability up to the Planck scale. As a result, the electroweak vacuum is unstable at high energies. However, the most up-to-date measurements of the top quark mass and the strong coupling constant result in a scalar potential that is metastable with a sufficiently long-lived electroweak vacuum [8,9]. The exact value of the scale, which we denote by SM

I
where λ H becomes negative strongly depends on the value of the top quark mass and the strong coupling constant and it is well known that at the 2-loop level in the renormalization group equation evolution of SM couplings the value of λ H can turn negative at scales in the range ∼ 10 10 -10 13 GeV. These remarks hold if we assume that no new physics enters in extrapolating the SM to such high energies.
Below we study the electroweak vacuum stability in Type I seesaw models for three generations of neutrinos in scenarios where the right-handed neutrinos have explicit bare mass terms in the Lagrangian and where these are dynamically generated through the mechanism of spontaneous symmetry breaking. We compare the two scenarios for the case of absolute stability of the scalar potential.

Stability in the Type I seesaw
Despite its spectacular success the SM cannot be the complete theory of nature. We now have convincing evidence that neutrinos have small but finite masses. Within the SM this can arise by incorporating a dimension five effective operator κ L L H H, where L denotes a SM lepton doublet. After electroweak symmetry breaking neutrinos obtain a mass given by κv 2 2 . By dimensional arguments, κ scales as 1/M where M is the new physics scale where the above operator is generated. The seesaw mechanism is the simplest manifestation of the above idea. One adds two or more right-handed neutrinos, N i , where i > 2 to the SM. In most scenarios, i is taken to be 3 for symmetry reasons and for concreteness this is what we assume. In this simplest version the Lagrangian that gives rise to neutrino masses is given by where as a result of integrating out the heavy right-handed neutrinos [10][11][12][13][14][15]. For energies below M, κ is a running parameter [16,17] and its one-loop RGE is approximately given by In the equation above, g denotes the SU(2) W gauge coupling and y t the top Yukawa coupling. This is the only new parameter added to the SM and its value is very small since it has to yield the active neutrino masses. Its contribution to the running of λ H is of order v 2 κ 2 and thus negligible. The coupling κ increases with energy due to y t being large albeit slowly. Moreover, the rest of the SM couplings run undisturbed. However, above the scale M the neutrino Yukawa couplings start to run and affect the running of λ H much like the top quark Yukawa coupling. Hence, the stability of the electroweak vacuum will set a limit on how large Y ν can be if we demand that the presence of neutrinos does not make the electroweak vacuum unstable too soon. This was first studied in [18] and later in [19] with an emphasis on the Dirac Yukawa matrix dependence. Furthermore, the work in [20] covers in depth the case where the scalar potential is metastable with a lifetime longer than the age of the universe and in the presence of three degenerate right-handed neutrinos. In this work, we will extend their studies using the full 2-loop RG running of the SM couplings and Y ν and absolute stability in order to contrast with the next scenario. We also implement realistic neutrino mass matrices that encode the data from neutrino oscillations experiments [21]. For simplicity we take M M to be proportional to the unit matrix with scale given by M N .
To be consistent we require that M N < I where the electroweak vacuum becomes unstable. An important quantity in our analysis is i m i where m i denotes the active neutrino masses. An upper limit on the above sum constrained by cosmology and astrophysics is given in [22]. We are mindful that the above is not without theoretical assumptions, i.e. the validity of the CDM cosmology model. For degenerate neutrino masses we have m cosmo = i m i = 3m + δm 2 2m + m 2 2m where m is the common mass and δm 2 and m 2 are the solar and atmospheric neutrino oscillation frequencies, respectively. The most recent values for these parameters can be found in [23]. . In both cases the smaller δm 2 term is dropped. Furthermore, we use the neutrino mass matrix elements found in [21] and define at the Majorana mass scale y 2 = 2 M N v 2 m where = e, μ, τ . In addition, our study is carried out using the RGEs outlined in [24] and where to one-loop order, the running of λ H , Tr Y † ν Y ν and the diagonal elements of the neutrino Yukawa coupling matrix are approximately given by where i = e, μ, τ and θ M N ≡ θ(μ − M N ) is a step function accounting for the threshold, M N , at which the new neutrino Yukawa couplings begin to run. We have ignored the contributions from the bottom quark and τ lepton Yukawa couplings.
In Fig. 1a we show the value of I as a function of M N using a value of the strong coupling constant α s = 0.1186 ± 0.0006 [25] and a top quark mass of 173.21 GeV [26]. The green circles correspond to a light neutrino mass consistent with the cosmological upper bound of i m i < 0.23 eV [27], while the red squares and brown diamonds correspond to the values of m 0 (m) = 2, 10 eV, respectively. The latter value is ruled out by current experimental bounds from tritium decays [28,29] and we intend it to be for illustration purposes only. For a 0.05-2 eV neutrino, the Majorana mass scale is bounded by 10 10 GeV. Beyond this scale, the value of λ H at the Majorana mass scale is negative. This places an upper bound on the value of TrY † ν Y ν at the Majorana mass scale of ∼ 10 −4 for M N ∼ 10 10 GeV and m 0 (m) = 0.07 eV. We show this in Fig. 1b. Our results is in general agreement with that of [20]. However, we used realistic neutrino mass matrices and did not assume degenerate active neutrinos.
Beyond a Majorana mass scale of ∼ 10 10 GeV the value of λ H (M N ) is negative and high-scale seesaw models break down if the scalar potential were to remain stable as we require. We also see that the absolute value of Y ν is of order 0.01. This is three orders of magnitude lower than the perturbative bound. This is valuable information since we know of no direct measurement of Y ν that can be made due to the very massive nature of the N R 's. This is a well known problem of the high-scale seesaw model. Yet, neutrino Yukawa couplings are vital for Type I seesaw models. In attempts to circumvent this, more elaborate schemes such as low-scale seesaw models [30][31][32][33][34][35][36], inverse seesaw models [37,38] and also the left-right symmetric models [40] have been introduced. The latter gives rise to signatures with same sign leptons plus jets that can be searched for at LHC. A detailed study using a simplified model approach was recently given in [41] where many references can be found. As for the highscale case, which has the virtue of being simple, one can only rely on theoretical studies. Indeed we conclude that for the current preferred value of α s , electroweak stability would lead to seesaw scales approximately six orders of magnitude lower than the Grand Unified theory scale with neutrino Yukawa couplings of order O 10 −2 . Recently, the authors in [39] have analyzed the stability of the electroweak vacuum in the presence of a low-scale seesaw model. They find that low-scale seesaw models are viable and do not disrupt the stability of the Higgs potential for Tr[Y † ν Y ν ] < 0.4. Their results are complementary to our findings.

Stability in the Type I seesaw with a complex electroweak singlet scalar
In Eq. (1) the heavy right-handed neutrino masses are introduced by hand and represents the case of explicit lepton number breaking. However, these masses can also be generated by spontaneous symmetry breaking of lepton number. The simplest model that achieves this involves adding a complex SM gauge singlet scalar field, S. It has a Yukawa coupling to right-handed neutrinos given by Y N N R N c R S and as such preserves a global U (1) L symmetry responsible for lepton num-ber if S has a global lepton number charge of two units. Upon breaking this symmetry a Goldstone boson, the majoron, will emerge [42,43] which can serve as a candidate for dark radiation [44,45]. This scalar couples to the Higgs field via the term λ H S S † S H † H . One thus expects the running of the SM couplings and the stability of the scalar potential to be different. In the following we will address these issues and present as a detailed RG analysis of this model.
The embedding of the majoron model into a more complete model is not the purpose of this work. Indeed the U (1) L can be replaced by any U (1) X . A well known example is a gauged X = B − L symmetry where the singlet scalar can serve as the inflaton. Here we focus on the effects of a complex scalar gauge singlet on the electroweak vacuum stability captured by the simple majoron model and compare the results with the explicit mass case studied above. The majoron plays no role in this investigation.
The addition of a complex scalar S, singlet under the SM gauge symmetries and charged under a global U (1) L lepton number symmetry, can be parametrized by a scalar potential given by Within this framework, for λ H , λ S > 0 and λ 2 H S < λ H λ S , the minimum of the potential is at H = v/ √ 2 and S = w/ √ 2. The mixing between the two neutral CP even scalars is given by In this work, we are interested in the limit where the lightest scalar is the SM-like Higgs boson. This limit is characterized by a large singlet vacuum expectation value, w >> v, and masses for the two CP even scalars given by This also implies a small mixing of the SM-like Higgs with the heavy scalar and hence only small corrections to the Higgs couplings to other SM particles. However, the presence of a heavy scalar, coupling at tree-level to the SM-like Higgs boson, may lead to a large positive contribution to the RGE for λ H [46][47][48][49][50][51][52][53] as well as a tree-level threshold effect [54] that arises from the matching effect of the singlet at the energy Q ∼ m S . Both are known to affect the stability of the Higgs potential. The authors in [54] have studied the effects from a threshold corrections to the Higgs quartic coupling, δλ = λ 2 H S /λ S against those that arise from positive contributions in the running of λ H . They studied two interesting cases: One where λ H S > 0, where large threshold corrections increase the scale at which the Higgs potential develops an instability and the second where λ H S < 0. In the latter, the tree-level threshold effect is not sufficient to increase the instability scale and RG effects become important since the new stability condition, λ H > δλ must be satisfied at all scales.
We wish to study the effects on the instability scale of the SM scalar potential in the presence of a new heavy scalar, as the one discussed above, and also incorporate three Majorana right-handed neutrinos. Within this framework, a Majorana mass is dynamically generated once the new scalar develops a vev with the coupling Y N N R N c R S. For simplicity we will take Y N to be real. The RG evolution of λ H will be modified at the scales M N and m S due to the presence of new fermion and scalar degrees of freedom, the former which will tend to destabilize the scalar potential. We analyze both the case where λ H S > 0 and where λ H S < 0 and compare our results to a model with only a scalar gauge singlet.
The scale of instability can be calculated using the RG equations for the scalar and new Yukawa couplings which at one-loop order are given by where θ M high , is a step function accounting for the threshold, M high = Min{M N , m s }, at which the new Yukawa coupling the singlet to the right-handed Majorana neutrinos and scalar potential parameters begin to run. We assume that S couples diagonally to all three right-handed neutrinos with universal strength Y i = Y N for all i, and the vev of S sets the Majorana mass scale given by The running of the Yukawa couplings between left-and righthanded neutrinos to the SM-like Higgs, Y ν , are given as Eq. (3) with M N → M high . The instability scale for λ H S > 0 is defined as the energy scale where either λ H or λ S first vanish, rendering the scalar potential unstable and invalidating the seesaw mechanism in models were one demands absolute stability of the scalar potential. For λ H S < 0, the condition λ H > δλ must also be met for all scales together with λ S > 0. By inspection of Eq. (7) we expect that for large values of Y N M high , λ S will be driven toward negative values very rapidly given that it runs with the fourth power of Y N . In fact, for λ H S > 0, three effects define the scale where the scalar potential becomes unstable: There is a second effect that disfavors large values of Y N . The last equation in Eq. (7) show that these couplings hit a Landau pole very quickly. If this happens within the stability region of the scalar potential, the theory would become strongly coupling. Although a very interesting scenario, it is beyond the scope of this work. For small values of Y N , λ S can also arrive at the Landau pole too soon. These conspire to limit the range of interesting values of Y N and λ S . A qualitative picture of these effects for λ H S > 0 is depicted in Fig. 2a-c  finite nature of the various curves denote the following: If δλ increased either we obtain stability up to the Planck scale with couplings either perturbative or non-perturbative at the Planck scale, or one of the quartics goes to zero and then hits a Landau Pole driving the other couplings toward a Landau Pole. Basically, a larger δλ leads to a larger λ H S , which also affects the running. From the figures it is clear that a large λ S M high aids a lot at stabilizing the Higgs potential for small values of Y N M high . This is due to the positive contribution to the running of λ H , which together with the tree-level threshold effect help increase the stability of the scalar potential compared to SM I ≈ 10 10 GeV. However, as Y N M high increases, the effect from the tree-level threshold effect shift remains fairly constant since λ S goes to zero near M high for small λ S M high . For larger λ S M high , the contributions from λ S to the running of λ S for a large range of energies and it is eventually overtaken by the contribution from Y N . This aids at stabilizing the scalar potential at high We note that our results suggest that λ S cannot be too small unless Y N is simultaneously small since otherwise λ S will be driven negative very quickly. To see this we scan over λ S M high and Y N M high and calculate the scale of instability. Our results are shown in Fig. 5 for δλ = 0.015 and λ H S > 0. The gray region corresponds to regions of parameter space where I 10 10 GeV for m S = 10 4 GeV. The regions in light red, green and blue are stacked behind the gray region and correspond to m S = 10 6 , 10 8 , 10 10 GeV, respectively, with the light blue region extending over the entire grid. Next we turn our attention to Y ν . In the above we have seen that spontaneously generating heavy Majorana neutrino masses can raise the instability scale of the electroweak vacuum with the help of the added scalars; especially for λ H S > 0. One can then increase the upper bound on Tr[Y † ν Y ν ] compared to the case of explicit mass terms. This can be seen in Fig. 6 with Y N = 0.1. Indeed the bound is Tr[Y † ν Y ν ] 0.01, higher than before but is only for a relatively large M high . In this figure, only points which lead to an instability scale below or at the Planck scale with perturbative couplings up to the Planck scale are shown. Regions of parameter space for which couplings are not perturbative below the Planck scale are not shown. These points tend to be for low M high . For very large M high , the couplings run as in the SM and the instability scale is below M high . This region of parameter space corresponds to points where M high 10 10 GeV, and are not accepted either. The window of allowed values of M high gets smaller for larger values of λ S . This can be seen in Fig. 7a, b with λ S = 0.6 and 0.8, respectively. Values of λ S beyond ∼ 0.9 do not lead to a viable seesaw. Furthermore, increasing Y N only decreases the upper bound on λ S .

Discussion
We conclude that in the two heavy Majorana mass generation mechanisms we have studied, the neutrino Yukawa couplings will remain perturbative. We find that for high-scale Type I seesaw models in which lepton number is explicitly broken by Majorana bare mass terms the stability is lowered due to the Dirac Yukawa couplings, Y ν . We also find that Tr[Y † ν Y ν ] 10 −3 in order for the SM electroweak vacuum stability not be worsened. In the scenario where right-handed neutrinos masses are generated spontaneously through a singlet scalar vev, the electroweak vacuum can be stable up to the Planck scale through threshold effects from decoupling a heavy scalar and positive contributions to the running of λ H , which ultimately depends on the scalar mixing parameter λ H S . We find that the scale at which the scalar potential becomes unstable is not improved over SM I if the quartic coupling of the singlet is small while Y N is large. This result is worsened for light singlet scalar masses; since λ S will be driven negative at very low energy scales. We also set the limit Tr[Y † ν Y ν ] 0.01. This is because both Y N and λ S must not be too large as not to hit the Landau pole within stability region. Thus, the viable boundary values of M N are low and consequently small values of Tr[Y † ν Y ν ] from κ.