Minimal scenario of Criticality for Electroweak scale, Neutrino Masses, Dark Matter, and Inflation

We propose a minimal model that can explain the electroweak scale, neutrino masses, Dark Matter (DM), and successful inflation all at once based on the multicritical-point principle (MPP). The model has two singlet scalar fields that realize an analogue of the Coleman-Weinberg mechanism, in addition to the Standard Model with heavy Majorana right-handed neutrinos. By assuming a $Z_2 $ symmetry, one of the scalars becomes a DM candidate whose property is almost the same as the minimal Higgs-portal scalar DM. In this model, the MPP can naturally realize a saddle point in the Higgs potential at high energy scales. By the renormalization-group analysis, we study the critical Higgs inflation with non-minimal coupling $\xi |H|^2 R$ that utilizes the saddle point of the Higgs potential. We find that it is possible to realize successful inflation even for $\xi=25$ and that the heaviest right-handed neutrino is predicted to have a mass around $10^{14}$ GeV to meet the current cosmological observations. Such a small value of $\xi$ can be realized by the Higgs-portal coupling $\lambda_{SH}\simeq 0.32$ and the vacuum expectation value of the additional neutral scalar $\langle\phi\rangle\simeq 2.7$ TeV, which correspond to the dark matter mass 2.0 TeV, its spin-independent cross section $1.8\times10^{-9}$ pb, and the mass of additional neutral scalar 190 GeV.


Introduction
The observed Higgs mass supports the assumption that the Standard Model (SM) is not much altered up to the Planck scale. More precisely, the critical value of the top-quark pole mass is about m pole t,critical 171.4 GeV [1] for the theoretical border between stability and instability (or metastability) of the effective Higgs potential for the observed Higgs mass 125 GeV; see also Refs. [2,3,4]. 1 This critical value of the top pole mass is consistent at the 1.4 σ level with the latest combination of the experimental results m pole t = 172.4 ± 0.7 GeV [5]. Surprisingly, the degenerate minimum of the Higgs potential at the critical top mass coincides with the Planck scale, and such a behavior of the potential has a lot of implications to high energy physics and cosmology. This interesting behavior of the Higgs potential can be understood by the multicritical-point principle (MPP) [6,7,8,9,10,11,12,13,14,15,16] (see also Refs. [17,18]) that "coupling constants that are relevant at low energies are tuned to a multicritical point around which the vacuum structure drastically changes when they are varied." We note that existence of a saddle-point around the Planck scale, rather than a degenerate vacua, is another possible form of multicriticality. This fact is used in the critical Higgs inflation [19,20,1,21,22] explained below.
Besides such interesting behavior of the Higgs potential, there are many mysteries and problems in particle physics and cosmology. For example, we have not yet understood the origin of electroweak (EW) scale v = 246 GeV, which is hugely small compared to the Planck or string scale 10 18 GeV at which people believe that there must exist an unified theory which includes quantum gravity. Further, the Majorana-mass scale for the right-handed neutrinos is unknown in the SM with the seesaw mechanism [23,24,25,26,27]. Moreover, the recent observations in cosmology, including that of the cosmic microwave background (CMB), have established the existence of (cold) dark matter (DM). It motivates us to consider a new particle whose interactions with the SM particles are relatively weak.
In addition, the CMB fluctuations may also provide hints for further new physics since they are seeded at high energy scales during inflation. Current observation is consistent with single-field inflation models, among which the Higgs inflation provides one of the best fits [28,29,30,31]. In particular, the critical Higgs inflation is inspired by the possible existence of the saddle point around the Planck scale in the MPP, which helps to flatten the Higgs potential and allows rather small value of the non-minimal coupling ξ = O(10) in ξ|H| 2 R. (The not-so-large coupling is favorable from unitarity [32,33,34,35,36].) In this paper, we consider the most economical model that can simultaneously explain all the above issues: the critical Higgs inflation, neutrino masses, EW scale, and DM. 2 The model consists of two additional real singlet scalar fields and the SM with the right-handed neutrinos. As a result, we manage to predict the parameters in a well determined narrow region by taking into account the constraints from the DM relic abundance, its direct detection experiments, the CMB fluctuations, and the latest LHC data, while keeping the perturbativity up to the Planck scale. We find that the Higgs-portal coupling and the vacuum expectation value of the additional neutral scalar are fixed to be λ SH 0.32 and φ 2.7 TeV, resulting in the dark matter mass 2.0 TeV, its spin-independent cross section 1.8 × 10 −9 pb, and the mass of 1 With the current central value mH = 125.1 ± 0.1 GeV [5], the critical top mass becomes m pole t,critical 171.2 GeV [4]. 2 Quite recently, a similar scenario to simultaneously explain these problems has been studied in Ref. [37] in the context of the classical scale invariance, in which a fermionic DM candidate is considered. In our scenario, we consider the MPP and a scalar DM candidate. additional neutral scalar 190 GeV.
The motive behind this model is as follows. If m pole t > m pole t,critical , the instability of the Higgs potential requires additional positive contributions to the renormalization group (RG) of the Higgs quartic coupling λ H . The simplest possibility is to introduce a real scalar field S that couples to the SM Higgs doublet H via a quartic interaction λ SH |H| 2 S 2 . (This is nothing but the Higgs-portal DM model [38,39,40,41].) However, such an extension could yield too large a tensor-to-scalar ratio of the CMB in the (critical) Higgs inflation, because of the raised Higgs potential at high scales; see e.g. [42,22]. We can resolve this issue by introducing additional superheavy fermions that lowers the Higgs potential at higher scales, while keeping the above mentioned positive contribution from the scalar at lower scales. Therefore, it is reasonable to assume a high-scale seesaw in which right-handed neutrinos have large Majorana masses M R and large Yukawa couplings y ν . In fact, it is possible to maintain the saddle point of the Higgs potential at high energy scales in the existence of new scalar field(s) when M R ∼ 10 14 GeV [22]. This is also one of the interesting predictions of the MPP.
So far, the MPP has not explained the origin of the EW scale. In fact, we may do so as follows. It is known that the Coleman-Weinberg (CW) mechanism [43] can naturally explain the hierarchy between the EW and Planck scales through the dimensional transmutation. Important assumption behind the CW mechanism is that the renormalized mass-squared parameter vanishes at the origin of the scalar field space. This assumption is called the classical scale invariance (CSI) [44,45,46,47,48,49,50,51,52,53,54,55,56]. 3 On the other hand, such CSI-type potentials can be naturally understood as one of the possible multicriticities in MPP without referring to scale invariance [61,62]. 4 From this MPP point of view, the simplest realization of the dimensional transmutation is achieved by only two real singlet scalar fields, one of which is S mentioned above [61,62].
Although we will focus on a specific model proposed in [61,62] (namely only the CP 2-2 among various multicritical points in the parameter space) in the following, the analysis of Higgs inflation does not much depend on the details of the model because only the scalar coupling λ SH and the neutrino Yukawa y ν play important roles to determine the behaviours of the Higgs potential at high energy scales. In this sense, the same analysis is easily applicable to similar extensions of the SM.
The organization of the paper is as follows. In Sec. 2, we briefly explain the minimal model of dimensional transmutation [61,62] extended with right-handed neutrinos and study the RG. In Sec. 3, we study the saddle point of the Higgs potential at high scales. In Sec. 4, we discuss the critical Higgs inflation. In Sec. 5, we show the method and results for our numerical prediction for the inflationary observables. Summary and discussion are given in Sec. 6. In Appendix A, we list the two-loop renormalization group equations (RGEs). In Appendix B, we summarize basic results for a general single-field slow-roll inflation. In Appendix C, we summarize basic results for the ordinary (non-critical) Higgs inflation. In Appendix D, we show analytic results of expansion around the saddle point.  Figure 1: Typical examples of the effective potential resulting from the MPP.

Model
In this section, we introduce the minimal model for the EW scale, neutrino masses, DM and the critical Higgs inflation. The model is based on [61,62] whose scalar sector contains two additional real singlet scalars S and φ. We also take into account heavy Majorana righthanded neutrinos ν i R with three generations. In order to make S a DM candidate, we impose the Z 2 symmetry S → −S, with all the other fields being invariant.
Here we briefly review the consequence of the MPP, focusing on what kind of effective potentials the MPP can result in. Before proceeding, we first note that the effective potential exists independently of the renormalization scale, up to the renormalization of the field, for a given set of bare parameters. (In other words, the renormalization scale dependence should only arise due to truncation at some loop level.) In Fig. 1, we list examples of the possible effective potentials resulting from the MPP: In panel A, we show the first kind of tuning V φ=φours = V φ=φ another given in the original version of the MPP [6,7,8]. In panel B, we show another one-parameter tuning V = V = 0, which provides an explanation for the CSI-type potential V | φ=0 = V | φ=0 = 0. This is also the key assumption for the Coleman-Weinberg mechanism. From this point of view, the CSI-type potential is just one of the possible critical points, in the theory space, in the MPP. In panel C, we show the potential (called CP 2-2) that is a consequence of two-parameter tuning to realize two saddle points V | SP1 = V | SP1 = 0 and V | SP2 = V | SP2 = 0. We stress that all the above potentials are achieved in the context of the MPP irrespectively of any kind of scale invariance. Therefore, it does not matter whether other fields such as the right-handed neutrinos have mass parameters in the action or not.
Among the various critical points analyzed so far, we take the CP 2-2 defined in Ref. [62] that has two saddle points in the effective potential, because the widest parameter region is allowed by the constraints considered. The renormalized Lagrangian is where L SM is the SM Lagrangian without the Higgs potential and without the right-handed neutrinos, and H c ≡ iσ 2 H * . Note that all the couplings in the Lagrangian are defined by the zero-momentum subtraction scheme. These couplings are the same as those in the effective potential. It is precisely these couplings that are restricted by the MPP. Here, among several realizations of the MPP, we choose a class of simple solutions that satisfy the vanishing of all the terms with dimensionful couplings except for φ 3 , namely, φ, φ 2 , S 2 , H † H, φH † H, and φS 2 . In Eq. (1), we have omitted these terms from the beginning. Note that the renormalized mass of S around the origin h = φ = S = 0 is set to zero by the MPP, while it obtains the finite mass around the true vacuum as we will see in Eq. (9). For simplicity, we will make the following three assumptions. First, we take λ S to be zero at the EW scale since it is irrelevant for our discussion. 5 Second, we assume that the Majorana Yukawa couplings y φ ij are not large and do not affect the other RG runnings. 6 Third, we assume that Majorana masses are all degenerate each other, M R1 = M R2 = M R3 = M R , and neglect the mixing in the neutrino Dirac Yukawa couplings, i.e., y νij = y νi δ ij . At the leading order of M −1 R , the neutrino mass matrix becomes We call the number of "heavier active neutrinos" n ν , namely the degeneracy of the largest y νi : In the case of normal and inverted hierarchies, we have n ν = 1 and 2, respectively, and when all the neutrino masses are degenerate, we have n ν = 3. Throughout this paper, we will use m ν = 0.05 eV [5] as the heaviest mass of active neutrinos. 7 The relation between M R and the largest y ν is then determined to be y ν 0.41 M R /10 14 GeV.
Here, we briefly review how the CW mechanism realizes the EW symmetry breaking in our model with the MPP. The effective potential for the φ-H system can be expressed at one-loop level as where M 2 φ (φ) = µ φ φ + λ φ φ 2 , and the renormalization scheme applied in [61] is used. The ∆V 1−loop term is the 1-loop potential for h given in Eq. (13). In this expression, we assume that the H-loop contribution to the φ 4 term as well as the field dependent masses of φ and S coming from H are negligibly small, which can be realized by taking λ φH λ φS and λ SH H † H λ φS φ 2 . This potential can be rewritten at the scale µ = µ * where λ φ vanishes as In the following, we also assume that the µ 2 φ φ 2 term is much smaller than the λ 2 φS φ 4 term. 8 Requiring the existence of the two saddle points at φ = 0 and φ = φ saddle (CP2-2) 9 , µ φ is determined to be The vacuum expectation value (VEV) of φ, φ , is then determined as where W is the Lambert W function. Using Eqs. (5) and (6), we obtain the relation between µ * and φ as Now, it is clear that by looking at the first line of Eq. (4) the EW symmetry breaking is triggered by φ : where v 246 GeV. At the minimum, the squared masses of φ and S are given by 8 Note also that other scalar couplings such as φ|H| 2 and φS 2 , which are set to zero by the MPP, are induced at the one-loop level as well. However, in the present analysis, the mixing between H and φ is small, and since we are interested in S = 0, these effects are small. 9 Note that this saddle point is nothing to do with the saddle point of the Higgs potential discussed in the next section.  where we have neglected the mixing between h and φ. From the above discussion, the six parameters λ, λ φ , λ φS , λ S , λ φH , and λ SH in the potential are reduced by fixing the two parameters, v and the Higgs mass 125.1 GeV, and by taking λ S to be zero at the EW scale as aforementioned. The resultant three parameters are where λ φ has been replaced by the scale µ * at which λ φ vanishes, and further converted to φ by Eq. (7). The thermal relic abundance of S is satisfied when [62] This relation provides a contour of φ in the m S vs λ SH plane as shown in Fig. 2: The red, magenta, green, and blue curves correspond to φ = 2.5 TeV, 3 TeV, 4 TeV, and 10 TeV, respectively. We also show the other constraints on the model parameters in the figure: The purple, orange, cyan, and yellow shaded regions are excluded by the updated XENON1T result [64], 10 the LHC results, DM relic abundance, and perturabativity bound, respectively [62]. We note that these bounds are insensitive to M R . 11 In the following analyses of the critical Higgs inflation, we choose the parameters on the dotted line in Fig. 2 that is close to the perturbativity bound, to reduce the computational time. 12 On the dotted line, there remains only single parameter in the scalar sector, and we choose it to be λ SH . Therefore, there are two parameters λ SH and M R in total. We choose these parameters in such a way that the Higgs potential has a (near) saddle-point at high scale region. 13

Saddle point of Higgs potential
Now that the TeV-scale couplings are obtained, we extrapolate them towards high scales. In this section, we analyze the effective Higgs potential for large field values, and look for its saddle point in order to realize the critical Higgs inflation.

Effective potential
We calculate the one-loop effective Higgs potential improved by the two-loop RGEs presented in Appendix A. The one-loop effective potential for large h in the MS scheme in the Landau gauge is where in which the effective masses are and 12 When we vary mS from the dotted line for a fixed λSH , the coupling λ φS changes via Eq. (9), but the running of λ and λH do not depend on λ φS at the one-loop level (see Eq. (56)), and hence the high-scale Higgs potential is not altered drastically. Note also that, even though the dotted line is close to the perturbativity bound, λSH (which contributes to running of λ and λH at one-loop) remains perturbative unlike λS, and hence the effective Higgs potential is reliably computed. 13 Without the right-handed neutrinos, the high-scale potential value tends to become too large to accommodate the observed value of the tensor-to-scalar ratio; see also Ref. [22].
is the Higgs field with field renormalization. See Eq. (64) in Appendix A for the one-loop result of γ H . We here represent N i as the mass eigenstates for the heavy Majonara neutrinos. We neglect the contributions from the leptons and quarks other than the top quark because their Yukawa couplings are small. We also neglect the loops of the Higgs and the NG bosons because λ H becomes small at high scales.
In the following, we match two theories with and without right-handed neutrinos around µ = M R , to obtain the threshold correction. We expand the one-loop correction from heavy neutrinos by h as The term containing ln(M R /µ) leads to the subtractions of y 4 ν term from the beta function. The same analysis can be also applied to the field renormalization as from which we can see that the canonically normalized Higgs field h c in the low-energy theory is given by This also leads to the redefinition of the quartic coupling. By combing both of them, the Higgs quartic coupling below µ = M R is 15 One can easily check that contributions from heavy neutrinos cancel out in the beta function dλ/d ln µ.
14 The coefficient of h 2 in Eq. (17) gives the threshold correction to the Higgs mass-squared parameter in the low-energy theory. At the CP 2-2, the mass-squared parameter including this correction is tuned to zero, based on the MPP. See Refs. [61,62] for the detailed discussion. 15 If we want, we may take into account the threshold correction of S too: ∆λ This contribution is minor for our analysis and we neglect it hereafter.

Saddle point
In the following, we rewrite the one-loop Higgs effective potential (12) as where Note that λ eff (h, µ) is independent of the renormalization point µ if we take all order loop corrections into account. For practical purpose with finite-loop, it better approximates by choosing µ ∼ h.
In the SM without the right-handed neutrinos, it is known that λ eff (h, h) in general takes a minimum value λ min at h = h min ( 10 18 GeV); see the blue dashed line in Fig. 3. Furthermore, by tuning the top mass, we can realize a saddle point V = V = 0 at h = h s ( h min ) [1]. (This tuning may be done within the 1.4σ experimental bound as stressed in Introduction.) However, this saddle-point itself is not sufficient to achieve a viable saddlepoint inflation because the resultant CMB fluctuations become too large [68,42,69]. By introducing the non-minimal coupling ξ|H| 2 R, a successful inflation can be achived around the (near) saddle-point even when ξ ∼ 10 [19,20,1]. This is called the critical Higgs inflation.
In this paper, we pursue the critical Higgs inflation in our model at the CP 2-2. The detailed analysis will be presented in the next section. In the remaining of this section, we look for the saddle point in our model. When we take into account the right-handed neutrinos, they first lower λ H above M R , and at higher scales, the extra scalar couplings become large and their contributions raise λ H again. As a result, we may have a minimum for λ H at a high scale as can be seen in Fig 3, e.g. with the orange curve for the normal and inverted hierarchy cases, and green for the degenerate case.
In general, the condition to have a saddle point is Again, we call the position of saddle-point h s . As said in the last of Sec. 2, there is only single parameter λ SH in the scalar sector. For each given λ SH , we numerically solve for M R that achieves the saddle-point condition (24), and obtain its position h s . In the upper left panel in Fig. 4, we show the Higgs potential in the case of the normal hierarchy where the differently colored contours correspond to the different values of λ SH . Here, the Higgs field h is normalized by h s . In the upper right panel, we plot the value of M R as a function of λ SH that is determined by the saddle point condition (24). In the lower left panel, we show h s as a function of λ SH where blue and red correspond to the normal and inverted hierarchy cases, respectively, while green to degenerate. One can see that h s has a strong dependence on λ SH , which causes the strong dependence of the value of the potential around the saddle point on λ SH in the upper left panel.
The above numerical result can be fitted by an analytic formula as follows. Around the minimum h = h min giving the minimal value of λ eff , we may always approximate as [1] Within this approximation, we see that the potential has a saddle point V = V = 0 at h = h s when and only when λ min is tuned to a critical value λ min = λ c where  That is, We fit b and h s from the numerical analysis above. As a consistency check, we also plot the Higgs potential whose effective coupling is replaced by the analytical one (28) in the upper left panel in Fig. 4 by dashed black lines. One can actually see that the analytical ones well approximate the effective Higgs potential around the saddle point. In the next section, we will use Eq. (25) with λ min close to λ c , as in Eq. (50), to study the critical Higgs inflation. In particular, the minimum value λ c is important because it also determines, along with h s , the value of the potential during the inflation. Qualitatively, small λ c will turn out to be favorable to realize successful inflation with small ξ 100.
In the lower right panel in Fig. 4, we show the numerical calculations of λ c as a function of λ SH . From this plot, one can see that λ c is a monotonically increasing function of λ SH and its order of magnitude is 10 −5 around λ SH ∼ 0.3. This corresponds to b = O(10 −4 ) by Eq. (27). We will see that the region with small λ SH is preferred to realize a successful inflation for small ξ. 16 Therefore, in the inflationary analysis, we will show the results in the region λ SH 0.3.

Critical Higgs inflation
In this section, we study the critical Higgs inflation of our model at the CP 2-2. We take full advantage of the saddle point of the Higgs potential. The existence of the saddle point makes it easy to obtain the sufficient number of e-foldings N even when ξ 100; contrary to the conventional case, the tensor-to-scalar ratio r does not have to be related to N as r ∼ 1/N 2 , and it can be sizable ∼ 0.05. In section 5.1, we will show our numerical calculations of CMB observables.
Here we comment on the other fields. During the critical Higgs inflation, the Higgs potential is flat, while other scalar fields have large masses due to the couplings λ φH and λ SH with the large Higgs field, so they do not play a role in the inflation analysis. As for the renormalization scale, since we are considering the case where only the Higgs field is large, we can consider the Higgs field h as a renormalization point and use the single-scale renormalization group.

Higgs inflation at classical level
We first review the Higgs inflation with non-minimal coupling ξh 2 R at the classical level [70,28]. We start with the Jordan-frame action where we have truncated the potential and Weyl factor at the quartic and quadratic orders, respectively: By performing the following redefinition of the metric, we obtain Then the Einstein-frame action becomes where is the potential in the Einstein frame. Since V cl is quartic, the following relation holds In the h → ∞ limit, we have Therefore for large h, the Einstein-frame potential U cl (h) becomes constant: This flat potential is used in the Higgs inflation. The relation between the canonically normalized Einstein-frame field χ and the Jordanframe field h is given by Under the slow-roll approximation, we obtain to fit A s = U/(24π 2 M 4 P ε V ) to the observed value 2.1 × 10 −9 at the e-folding N ; see Appendix C. We see that the typical SM value at low energy λ cl ∼ 0.1 requires large value of ξ ∼ 10 5 .

Higgs inflation including radiative correction
At the quantum level in the flat spacetime, we promote V cl (h) to the effective potential where V tree is the tree-level potential including the field renormalization and ∆V loop is the loop correction; see Eq. (12) for the 1-loop approximation. In this paper, we employ the Einstein-frame effective potential on the so-called Prescription I, where ∆U loop is obtained from ∆V loop in Eq. (40) by replacing all the effective masses M Ψ (h) with M Ψ (h) Ω for Ψ = W ,Z, t, νi, S, and φ: See Ref. [71] for the meaning of this prescription as well as ambiguity due to indeterminacy of non-renormalizable terms. Because the potential should not depend on µ (if one takes all order corrections into account), we may replace µ by µ/Ω: 17 When we truncate the loop correction at the 1-loop order, we can obtain better approximation by choosing µ to be around h to minimize the higher loop corrections: This is the expression we employ in the following. When y νi h is sufficiently larger than M R , all the effective masses are proportional to h, therefore ∆V 1-loop (h, h) is proportional to h 4 whose coefficient is independent of h (up to higher order corrections), and hence we obtain Thus, after taking into account the quantum corrections, the same relation still holds as Eq. (35) on Prescription I, and again the Einstein-frame potential (45) becomes constant for ξh 2 /M 2 P 1: Unless λ eff is particularly small, we still need large ξ to fit A s . In fact in the SM, it is known that λ eff (h, h) becomes small around h ∼ M P , which is an essential ingredient in the critical Higgs inflation.

Critical Higgs inflation
As can be seen from Eq. (45), when V(h) has a saddle point at h = h s , U(h) also has a saddle point at h =h s , withh s being determined byh s Ω(hs) = h s : where we have introduced This parameter is the ratio of h s to M P / √ ξ; the latter is the typical value of h above which the conformal factor Ω(h) starts to deviate from unity. One can see thath s approaches infinity as c s 1, which means that the small region around the saddle point h = h s is widely stretched, and allows a sufficient e-folding.
In the critical Higgs inflation, we assume that the high-scale Higgs potential is close to a one having a saddle point, namely, λ min in Eq. (25) is close to λ c given by Eq. (27): where we have parametrized the deviation from the saddle-point criticality by δ. Then the flat-space effective potential becomes Using Eqs. (45) and (48), we see that U approaches the constant value in the h → ∞ limit which determines the value of the potential during the inflation. In this paper, we will focus on c s ≤ 1 and take full advantage of the saddle point of the Higgs potential. In Fig. 5, we show the Higgs potential in the Einstein frame where different colors correspond to different values of c s . Here, we show the the normal hierarchy case with λ SH (M t ) = 0.5 for illustration. One can see that the region around the saddle point h = h s is more and more stretched as we increase c s toward unity. Therefore, for a given c s 1, we may always fit the e-folding N around h = h s Ω by tuning δ. 18

Prediction on inflationary observables
Here, we analyze the prediction on inflationary observables as we vary the parameters in the model. So far, we have three free parameters: λ SH , M R , and ξ. Recall that the scalar sector has only one parameter λ SH in our analysis on the red dotted line in Fig. 2.
On the other hand, we have seen that the Jordan-frame potential can be parametrized near the saddle-point criticality by λ c , δ, and h s as Eq. (50). Here, these three parameters are functions of the model parameters λ SH and M R .
In order to sweep the parameters near the saddle point, we use the results of Sec. 3.2: For each λ SH , we find the value of M R that gives the saddle-point criticality, δ = 0, as well as the corresponding parameters h s and λ c . See the upper-right, lower-left, and lower-right panels in Fig. 4 for M R , h s , and λ c , respectively. For a given λ SH , as we slightly change M R from the critical value, in general all of the parameters λ c , h s , and δ are modified. Here, we neglect the change of λ c and h s , and take into account the effect of non-zero δ.
Now we take into account the non-minimal coupling ξ. Among three parameters ξ, λ SH , δ, the last one is traded with e-folding N using Eq. (66). The observables n s and r are functions of ξ and λ SH once e-folding N is fixed by δ.

Results
In Fig. 6 We see that we typically have λ SH ∼ 0.32, which corresponds to the large Majorana mass as 5 × 10 13 GeV M R 2 × 10 14 GeV (52) from the upper right panel in Fig. 4. In Fig.7, we also show dn s /d ln k vs n s , where c s and ξ are again shown as in Fig. 6. In this case, the observational error of dn s /d ln k is still too large to constrain the inflaton potential.

Summary and discussion
Motivated by various fundamental issues in particle physics and cosmology, we have discussed the minimal model that can explain EW scale, neutrino masses, DM, and successful inflation at the same time. The model adds right-handed neutrinos to the two-scalar model in Refs. [61,62], which has been proposed to explain the origin of EW scale and DM. These two scalar fields give a minimal setup to realize an analogue of the CW mechanism. Assuming the Z 2 symmetry of a scalar, S → −S, it can be a candidate for DM, similarly to the Higgs-portal scalar DM model. Neutrino masses are naturally explained by the seesaw mechanism.
In this paper, we have analyzed RGEs, calculated the effective Higgs potential, and studied the critical Higgs inflation that uses the (near) saddle point of the Higgs potential at a high scale. The new scalar coupling λ SH between the SM Higgs and DM S can stabilize the Higgs potential even if the top mass is current center value m t = 172.4 GeV. In our model, it is possible to maintain the existence of saddle point by the neutrino Yukawa coupling y ν , and the saddle point condition relates the parameters of the model as is shown in Fig. 4.
By utilizing the saddle point of the Higgs potential, we have found that it is possible to realize successful inflation even for ξ ∼ 25 within the parameter space where all the necessary requirements are satisfied. As a result, we obtain λ SH 0.32 and φ 2.7 TeV, which correspondingly lead to the dark matter mass m S ∼ 2.0 TeV, its spin-independent cross section 1.8 × 10 −9 pb, and the mass of additional neutral scalar m H ∼ 190 GeV.
Finally, we mention testability of our model at collider experiments and future directions of this scenario. Since the DM should be as heavy as a TeV range in order to satisfy the relic abundance and the constraints from the direct searches, the detection of the extra Higgs can be an important probe of our model similarly to the Higgs singlet model. On the benchmark points shown in Fig. 1, the mass of the additional Higgs boson is predicted to be in the range of 70 GeV -200 GeV, so that it can be produced at future lepton colliders such as the International Linear Collider (ILC) via the Z boson strahlung process. Therefore, our model can be tested at the ILC and/or its energy upgraded version.
Our model can be well tested at the near future DM detection experiments such as XENONnT. If the whole region is excluded by them, one of the simplest extensions would be introduction of extra heavy fermions that are singlet under the SM gauge symmetry but are odd under Z 2 such that they lower the quartic scalar couplings in the RGE running to make the perturbativity bound milder.
It would be interesting to analyze all the possible criticalities along the line of the current work. Moreover, we can also come up with a lot of interesting phenomena within this model such as a possibility of producing primordial Black Hole by Higgs inflation [72,73,74,75], spontaneous leptogenesis [76,77,78,79,80], (p)reheating dynamics and so on. We would like to discuss those possibilities in future investigations. Finally, we comment on possible systematic errors introduced by higher-order corrections. We have used the one-loop effective potential to discuss the MPP at the TeV-scale, and twoloop RGEs to extrapolate from there to high scales. The higher-order corrections are further suppressed at approximately λ SH /8π 2 few percent, compared to the corrections currently considered. The same degree of corrections applies to the potential value ∝ λ H and the inflationary predictions derived from it.

Appendix B Single-field slow-roll inflation
Here we summarize basic results for the single-field slow-roll inflation.
The slow roll parameters are defined by where U is the inflation potential in the Einstein frame and the prime represents the derivative with respect to χ. The number of e-foldings from a field value χ * to the end of inflation is given by The CMB observables are given by within the slow roll approximations, where A s , r, n s , and dns d ln k are the scalar power spectrum amplitude, tensor-to-scalar power ratio, scalar spectral index, and its running.
at the pivot scale k * = 0.05 Mpc −1 . Under the slow-roll approximation, the scalar amplitude A s is given by Eq. (67). We note that r and the value of potential U is related each other by fixing A s to the observed value: U ∼ 1.5 × 10 −9 r 0.05 M 4 P ; see also Ref. [22].

Appendix C Ordinary Higgs inflation without criticality
For ξh 2 M 2 P 1 we have simple relations for the slow roll parameters and they provide one of the best fits to the CMB observations for the reasonable values of e-folding N = 50-60 (corresponding to the pivot scale). Qualitatively, the typical value of ξ can be estimated as follows. Putting the potential (51) with λ eff ∼ λ c ∼ 10 −6 into Eq. (69), one can easily check that ξ is around 30.

Appendix D Expansion around saddle point
For qualitative understanding, it is also helpful to derive the expansion of V around h s . We first expand as where λ 1 = λ c δ, λ 2 = 3λ c δ, λ 3 = 32λ c .
As we explain in Section 4.2, the Higgs potential in the Einstein frame also has a saddle point at h = h s Ω. We are interested in the parameter space c s ∼ 1 ⇔h s h s ∼ M P √ ξ , which guarantees the large field expansion of χ as a function of h. As a result, we have where ϕ := h/Ω. Then, U and its derivatives with respect to χ are where we have used Eq. (71). When c ∼ 1, the slow roll parameters are approximately given by where a = λ 1 /λ 0 = 1 + (β λ /4λ)| φ=φs . Compared to the conventional case, we have additional suppression factor δ thanks to the saddle potential. Note that, when c 1, we can no longer trust Eq. (78).