Scalegenesis and fermionic dark matters in the flatland scenario

We propose an extension of the standard model with Majorana-type fermionic dark matters based on the flatland scenario where all scalar coupling constants, including scalar mass terms, vanish at the Planck scale, i.e. the scalar potential is flat above the Planck scale. This scenario could be compatible with the asymptotic safety paradigm for quantum gravity. We search the parameter space so that the model reproduces the observed values such as the Higgs mass, the electroweak vacuum and the relic abundance of dark matter. We also investigate the spin-independent elastic cross section for the Majorana fermions and a nucleon. It is shown that the Majorana fermions as dark matter candidates could be tested by dark matter direct detection experiments such as XENON, LUX and PandaX-II. We demonstrate that within the minimal setup compatible with the flatland scenario at the Planck scale or asymptotically safe quantum gravity, the extended model could have a strong predictability.


I. INTRODUCTION
With the discovery of the Higgs boson [1,2] the standard model (SM) was complete. This brings us to the next stage in elementary particle physics. One of obvious issues is the lack of a dark matter candidate in the SM. At the present stage a little fact is known about features of the dark matter as an elementary particle. In particular, it is not clarified even that the dark matter is either fermionic or bosonic, so that a numerous possibility of dark matter candidates can be considered. Besides, the nature of the Higgs sector is still unclear although all coupling constants in the SM are determined. Due to this situation, a scenario, where dynamics of dark matter is related to that of the Higgs field, has been suggested.
Let us here discuss what one expects from the discovered Higgs boson mass. The observed Higgs boson mass 125 GeV indicates that the perturbative renormalization group (RG) flow of the Higgs quartic coupling constant with the top quark mass M t 171 GeV reaches to zero around the Planck scale M pl within the standard model (SM) [3,4]. This fact might indicate a compelling evidence for dynamics of particles from a high energy theory including quantum gravity if one assumes that the SM is valid up to M pl or effects of new physics do not drastically change dynamics of SM particles. This fact motivates us to consider the flatland scenario [5][6][7][8][9] as one of scenarios for an extension of the SM.
The flatland scenario imposes that all couplings for scalar fields, involving scalar masses and the Higgs portal coupling, vanish at the Planck scale, namely the scalar potential becomes flat above M pl . Such a scenario might be highly controversial from the viewpoint of low energy physics, whereas this could be a natural condition from the asymptotic safety scenario of quantum gravity [10][11][12]. The existence of a non-trivial ultraviolet (UV) fixed point for gravitational couplings realizes asymptotically safe gravity as a non-perturbatively renormalizable quantum gravity. Above the Planck scale, RG scalings of couplings are non-trivially modified from the canonical ones due to the anomalous dimension induced by graviton fluctuations. For a certain fixed point value, the scalar masses and the scalar couplings are suppressed and then these couplings become irrelevant parameters [13,14]. This fact enforces scalar interactions so as to be switched off until the Planck scale.
In this work, we propose a U(1) X extension of the SM compatible with both the flatland scenario and the existence of dark matter candidates. We add right-and left-handed Majorana fermions, a SM-singlet scalar field which couple to a U(1) X gauge field. The Majorana fermions have Yukawa coupling with the singlet scalar field. These Majorana fermions are stable and thus can become dark matter candidates. The ratio between the U(1) X gauge coupling and the Yukawa couplings is a key quantity for the generation of an expectation value of the singlet scalar field (or a finite scale), based on the Coleman-Weinberg mechanism [15]. As a consequence, the U(1) X symmetry is broken and then the corresponding gauge boson becomes massive, while the Majorana fermions acquire finite masses via the Yukawa couplings. In general, it is allowed to write the so-called kinetic mixing term between the U(1) Y gauge field in the SM and the U(1) X gauge field. Such a kinetic mixing plays a crucial role of the inducement of a negative Higgs portal coupling between the Higgs doublet field in the SM and the singlet scalar field. Thanks to this, the breaking of the U(1) X symmetry triggers the electroweak symmetry breaking.
The dark matter relic abundance in this model corresponds to that of the Majorana fermions. This con-straint can fix, for instance, the value of the ratio between the Yukawa couplings. At this point, there is only one free parameter, e.g. the U(1) X gauge coupling, in this model. This free parameter could be determined by the direct detection of the weakly interacting massive particle (WIMP). Hence, this model is testable in near futures.
This paper is organized as follows: In Section II, we briefly explain the basic idea of the asymptotic safety scenario for quantum gravity and its implications for the matter dynamics. We introduce the model compatible with the flatland condition and summarize both theoretical (from asymptotically safe gravity) and experimental constraints for this model. We explain the mechanism of the symmetry breaking in the flatland scenario in Section III and show that the electroweak scale is generated within this model by taking a benchmark point. In Section IV we investigate the relic density of the Majorana fermions as dark matter candidates. The spin-independent cross section between the Majorana fermions and a nucleon is shown with the current upper bound from the WIMP direct detection experiment. Section V is devoted to summarize this work. In Appendix A, we collect the beta functions at the perturbative one-loop level. The one-loop effective potential in this model is shown in Appendix B. We show the explicit forms of cross sections for annihilations of the Majorana fermions in Appendix C.

II. SETUP
In this section, we discuss the basic idea of asymptotically safe quantum gravity and briefly summarize its current status. In particular, we will stress that quantum graviton fluctuations drive scalar dynamics such that it behaves as a free theory above the Planck scale, and then the flatland condition is given as a consequence from decoupling of quantum gravity effects around the Planck scale. We propose a flatland model involving fermionic dark matter candidates.

A. Asymptotic safety scenario
As a UV-complete theory beyond the Planck scale, we assume asymptotically safe quantum field theory involving quantum gravity. Here, we start with general discussions on the fixed point structure in a theory space and an energy scaling of a coupling constant in RG flow in order to make our criterion for an extension of the SM.
For a theory space spanned by a set of effective operators O i whose (dimensionless) coupling constants are denoted by g i , one explores fixed points at which all beta functions β i ({g i }) vanish. One can easily find the Gaussian (or trivial) fixed point g i * = 0 which can be discussed by perturbative RG. In addition to such a fixed point, in several quantum systems, non-trivial fixed point g i * = 0 could exist. One of well known cases is the Wilson-Fisher fixed point [16] in the three dimensional scalar theory which describes a ferromagnetic phase transition. In this case, however, non-perturbative methods, e.g.expansion [16,17] and functional RG [18][19][20][21][22] should be employed to analyze the fixed point structure due to the strongly correlated system.
Once one finds a fixed point, the energy scaling of coupling constants g i at vicinity of the fixed point can be clarified. This is characterized by the critical exponent θ i such that the dimensionless coupling constant behaves as g i (k) ∼ k −θi in the RG flow, where k is the energy scale. More specifically, one can obtain critical exponents by evaluating the eigenvalues of the stability matrix T ij = ∂β i /∂g j | g=g * . Coupling constants with positive critical exponents are relevant parameters and grow up toward low energy regimes. The subspace spanned by relevant operators, which are called the critical surface, defines a UV complete renormalizable theory. The relevant coupling constants are free parameters, so that a system with a less number of relevant couplings has a higher predictability. In contrast, for negative critical exponents, coupling constants are irrelevant parameters which are driven as functions of relevant couplings. At the Gaussian fixed point, the energy scaling for an operator can be read as the canonical dimension of its coupling constant, while at the non-trivial fixed point, a large anomalous dimension is induced by quantum dynamics, so that the energy scaling of coupling constants highly deviates from the canonical one. In such a case, one has to evaluate the eigenvalues of the stability matrix in order to obtain the critical exponents.
Let us here turn to the discussion on the basic idea of asymptotically safe quantum gravity. It is known that quantum gravity based on the Einstein-Hilbert action is perturbatively non-renormalizable due to the requirement of an infinite number of counter terms [23]. Nevertheless, there is a possibility that quantum gravity could be formulated as a non-perturbatively renormalizable theory which is known as asymptotically safe quantum gravity [10][11][12]. As discussed above, for the asymptotic safety scenario for quantum gravity, the existence of a non-trivial fixed point for gravitational interactions plays a crucial role. A number of studies utilizing the functional RG method have been performed and have shown evidences for the existence of such a nontrivial fixed point; see reviews [24][25][26][27][28][29][30][31]. An important fact is that a finite number of positive critical exponents (θ i > 0) for gravitational couplings is observed [32][33][34][35][36][37][38][39][40][41][42][43][44], and then asymptotically safe quantum gravity could have a predictability to low energy dynamics.
A great interesting question is impacts of quantum gravity effects on matter dynamics. Within the asymptotically safe gravity scenario, a large anomalous dimension induced by graviton fluctuations could change drastically scalings of matter couplings above the Planck scale. Below a transition scale k t associated with the Planck scale where quantum gravity effects decouple, dy-namics of particles may be described by the SM with a simple extension of the SM. In this view point, the extended system describing the particle dynamics should satisfy the boundary condition given at k t .

B. Criteria for constructing model
We discuss constraints from the asymptotically safe quantum gravity scenario on matter dynamics and consider a possible extension of the SM.
As a simple extension of the SM, the inclusion of a singlet scalar field S coupled to the Higgs field can be considered. We first discuss the conditions for the RG flow of scalar interactions. The form of the beta function is given by (1) Here, λ denotes a scalar coupling such as the quartic and portal coupling constants, and β λ,matter includes contributions from matter dynamics, while f λ represents universal gravity contribution whose form reads [14,45] f λ G 8π 20 whereG = Gk 2 is the dimensionless Newton coupling constant and v 0 = 16πGΛ cc k 2 is the dimensionless cosmological constant. Note that the dimensionful Newton coupling constant G is mass-dimension −2, while the mass-dimension of the dimensionful cosmological constant Λ cc is 2. Looking for a fixed point at which β λ = 0, one finds the Gaussian fixed point for the scalar coupling, λ * = 0. The investigations for such a system indicate the facts that all scalar couplings involving the Higgs portal coupling are irrelevant [14] since the critical exponent for the scalar coupling is given by θ λ −f * λ < 0 where we assume that the gravitational couplings have a non-trivial fixed point. This means that the RG flow for the quartic coupling and the Higgs portal coupling keeps zero until the Planck scale when their fixed point is Gaussian, λ H * = λ S * = λ HS * = 0. Thus, we have the boundary conditions at k t = M pl for the scalar interactions such that where λ H and λ S are the quartic coupling constants of the Higgs and the additional singlet scalar fields, and λ HS is their Higgs portal coupling constant. These conditions means that the scalar fields behave as free fields above the Planck scale. Second, we consider the quantum gravity effects on a squared scalar mass parameter. Its beta function reads wherem 2 = m 2 /k 2 is the dimensionless squared scalar mass parameter. The first term on the right-hand side is the canonical scaling term which causes the so-called gauge hierarchy problem since it gives a large value of the critical exponent θ m 2 for the Gaussian fixed point at which β m,matter 0 and f m 0. In such a case, the squared scalar mass is relevant, and then its energy scaling is given as a growing up solution below the Planck scale. For energy regimes above the Planck scale, the critical exponent of the squared scalar mass could change towards a smaller value than canonical one because of the graviton fluctuations (f m = 0) in asymptotically safe gravity such that θ m 2 − f * m . Indeed, the form of f m is given by the same as Eq. (2) [45]. If a large anomalous dimension f * m > 2 is induced, the critical exponent of the squared scalar mass turns negative and thus the squared scalar mass is not a free parameter. For the nontrivial fixed point of the squared scalar massm 2 * = 0, the electroweak scale is explained by the resurgence mechanism in a perturbation [13], while for the Gaussian fixed pointm 2 * = 0, the squared scalar mass keeps zero within the RG flow until the Planck scale. So far, the latter case is typically observed by the functional RG analysis [14,45,46], so that for the singlet scalar extension of the SM, the squared scalar mass parameters should satisfy This situation corresponds to the so-called classical scale invariance at the Planck scale [47][48][49]. In this case, the electroweak scale has to be generated by the dimensional transmutation by the Coleman-Weinberg mechanism [15,[50][51][52][53] or strong dynamics [9,[54][55][56][57][58][59][60][61][62]. See also [63]. We call the generation of a scale "scalegenesis" in order to emphasize that a scale invariant theory generates a scale due to its quantum dynamics. The conditions (3) and (5) indicate that the effective scalar potential is flat above the Planck scale. This is called the flatland scenario [6,7]. In this scenario, the scalar interactions have to be generated by quantum effects. With this fact, the following two things should be satisfied: (i) the effective scalar potential is stable; (ii) the electroweak scalegenesis takes place due to the Coleman-Weinberg mechanism in the singlet-scalar sector. However, they cannot be realized by only the inclusion of a singlet-scalar field to the SM. In order to obtain the stable scalar potential, the quartic coupling constants have to be positive for large field values. A Yukawa interaction can play this role since the beta function of the quartic coupling constant includes the term proportional to the Yukawa coupling constant to the fourth power, β λ ⊃ −y 4 . The Higgs quartic coupling constant can be realized due to the effect of the top-quark Yukawa coupling constant, while for the singlet-scalar field, an additional fermionic degrees of freedom coupled to the singlet-scalar field is required in order for its positive quartic coupling constant to be generated. Generally, for the requirement (ii), the quartic coupling in the RG has to be turned to a negative value around energy scales near a vacuum expectation value S . A new U(1) gauge interaction with the singlet-scalar field could play such a role since the quantum correction +g 3 arises from the gauge interaction in the beta function of the quartic coupling constant. Therefore, for the electroweak scalegenesis in the flatland scenario, the ratio between the Yukawa coupling and the gauge coupling is crucial. In sect. III A, the explicit condition for the ratio to realize the the electroweak scalegenesis is discussed.
Let us here discuss quantum gravity effects on a U(1) gauge coupling and a Yukawa coupling. For a U(1) gauge coupling, here denoted by g, the beta function is given by where the correction from quantum gravity is found [64] to be We should note here that −f g takes a negative value for a non-trivial fixed point of the gravitational couplings. In this case, the matter contribution β g,matter and the gravity effect could balance. Consequently, one could consider two possibilities as UV complete scenarios [65,66]: One is that in the continuum limit the gauge coupling reaches to a Gaussian fixed point (g * = 0) at which the gauge coupling is relevant. In this case, the gauge coupling is a free parameter and behave as an asymptotically free coupling. Other is the case that a non-trivial fixed point g * = 0, at which the gauge coupling is irrelevant and asymptotically safe, i.e. the gauge coupling is predictable in the low energy regime. From these facts, in order for the gauge coupling to be a UV complete coupling, it cannot take a larger value than the non-trivial fixed point value. For a one-loop level of the beta functions for the gauge coupling, β g,matter β g,1-loop g 3 , the gauge coupling is bounded so that, at the transition scale In the same manner, the Yukawa couplings could have also an upper bound [41,67,68]. On the other hand, the ratio between the Yukawa coupling and the gauge coupling is constrained from the condition for the realization of scalegenesis as will be seen in Section III A. Together with the bound for the gauge coupling (8), this condition provides both upper and lower bounds for the Yukawa coupling. Therefore, in this work we do not consider the bound for a Yukawa coupling from the asymptotic safety scenario.

C. Model in flatland
Following the discussions in the previous subsection, we consider an extension of the SM based on the flatland scenario. As a possible extension, we propose an extended system with a SM singlet-scalar field and Majorana fermions coupled to an extra U(1) X gauge field. This case allows us to write a kinetic mixing [69] between U(1) Y gauge field in the SM and the additional U(1) X gauge field such that B µν X µν , where B µν = ∂ µ B ν −∂ ν B µ and X µν = ∂ µ X ν − ∂ ν X µ are the field strengths for the U(1) Y and U(1) X gauge fields, respectively. Then, the kinetic terms for these gauge fields are given by while the interactions between the gauge fields and a matter field are given through a covariant derivative, where the strong (QCD) and weak interactions are omitted. Here, we canonically normalize the kinetic terms (9) so that where F µν and G µν are the field strengths for a new gauge-field basis (B µ , X µ ) defined by the transformation, with ± = √ 1 ± . The second matrix on the righthand side in the first line of Eq. (12) is employed in order to canonically normalize the kinetic terms for the gauge fields, while the first one corresponds to a rotation with an angle of tan −1 ( + / − ). The latter transformation can be performed since the kinetic terms (11) is invariant under a rotation for a field basis. We see here that for the new field-basis (B µ , X µ ) defined by Eq. (12), the mixing effect corresponds to a scale transformation for X µ , whereas B µ is transferred to X µ through the mixing effect.
In the new field-basis (12), the covariant derivative (10) becomes where we define new gauge coupling constants, A field for which U(1) Y hypercharge is assigned interacts with X µ even if it has no U(1) X hypercharge. Hereafter we work in the field-basis (B µ , X µ ) and neglect primes on these fields and the coupling constant g X . We here give the Lagrangian for our model, where L SM denotes the Lagrangian for the SM without the Higgs mass term due to the condition (5). Here, L kin involves the kinetic part of new fields, where S is a singlet-scalar field. Majorana fermions χ R and χ L are introduced in order to avoid the gauge anomaly. These fields interact with the U(1) X gauge fields via the covariant derivative given in Eq. (10). The singlet-scalar field S is coupled to only X µ with a hypercharge Y X = 2. For the Majorana fermions χ R and χ L , their hypercharges are assigned as Y X = −1 for the U(1) X gauge field, but is not for the U(1) Y gauge field, i.e. Y = 0. The assignment of hypercharges Y and Y X for each field is summarized in Table I. In this setup, the Majorana fermions are stable particles after the breaking of the U(1) X symmetry into the Z 2 symmetry, so that they could be dark matter candidates [70][71][72][73][74][75][76][77][78][79][80][81][82][83].
The Majorana-type Yukawa interactions between S and χ R or χ L are given by The U(1) X symmetry prohibits the Majorana mass terms and the Dirac-type Yukawa interactions, while the Dirac type mass is forbidden by scale symmetry. After the singlet-scalar field has a finite expectation value S , these terms turn to the Majorana mass terms. Thus, Eq. (17) becomes origins of dark matter masses. The scalar potential, denoted by V in the Lagrangian (15), is given by where H is the Higgs doublet field. With the condition (5), the scalar mass parameters keep zero within their RG flows since the beta functions for the scalar mass parameters are proportional to themselves. Therefore, the scalar mass terms are not taken into account. In contrast, the quartic and the Higgs portal interactions are needed in order for the theory to be renormalizable within our Lagrangian (15). Nevertheless, the RG equations for their renormalized coupling constants have to satisfy the condition (3), so that these scalar coupling constants are not treated as free parameters. Note that a massless scalar theory is renormalizable [84].
Here, we briefly describe the scalegenesis in our model. A scale associated with both the electroweak scale and dark matter mass scale should be generated by radiative corrections, i.e., the Coleman-Weinberg mechanism. Within our present model (15), the Coleman-Weinberg mechanism in the singlet-scalar sector first could take place. We denote the generated vacuum by v S = √ 2 S . This generation of a scale triggers the electroweak symmetry breaking through a negative Higgs portal coupling, namely where H = (0, v H / √ 2) T . In order to obtain a finite v H , a negative value of λ HS has to be induced by quantum effects. In the next section III, we will see the occurrence of such a situation. For λ HS ≈ 0, one obtains the Higgs H . Since U(1) X symmetry is spontaneously broken, the X boson obtains a finite mass, The masses for the Majorana fermions are given by The difference between χ R and χ L in this model is characterized by only the Yukawa coupling constants. Therefore, one can concentrate on the case y R ≤ y L with no loss of generality. Finally, we mention the constrains on parameters involved in the present model. In addition to the three conditions (3) for the quartic coupling constants and the Higgs portal coupling constant, we have constraints from the observations [85] v H = 246 GeV, M obs H = 125.10 ± 0.14 GeV, For the top-quark mass, the MS mass is used. Note that the pole mass is M pole t = 173.1 ± 0.9 GeV. As mentioned above and discussed in Section IV, the Majorana fermions could be stable and then become dark matter candidates. The latest observation [86] reports that the dark matter relic density in the cosmological evolution is where "DM" is the abbreviation of dark matter. The relic density of the Majorana fermions has to satisfy this value. These constraints reduce from the five free parameters, g X , g mix , y t , y L and y R , to one parameter. As discussed in Section II B, there is an upper bound (8) for a U(1) gauge coupling. For the prediction of the observed value of gauge couplings within the SM, the gravitational contribution f g being of order 10 −2 is required [87]. It is shown in Refs. [65,66,88] that gravitational effects actually yield f g of this order. We use f g = 0.02 [89]. Within our model setup (15), the upper bound becomes where we used the beta function for g X given in Eq. (A2) and set g mix = 0. Using the beta function for g mix with g X = 0 the kinetic mixing is also constrained as Finally, we comment on phenomenological constraints for the kinetic mixing effect. The kinetic mixing coupling  (9) is constrained for a wide range of the extra gauge boson mass M X [90]. The Zboson mass in our system would become typically of order a few TeV. For such a mass range, the upper bound on is given from Z searches in the LHC experiments by looking at e + e − and µ + µ − channels. For M X = 1 TeV, we have < ∼ 0.1 [91]. The bound for heavier mass regions is relaxed such that, for instance, < ∼ 0.2 for M X = 2 TeV.

III. REALIZATION OF SCALEGENESIS
In this section, we study the mechanism of the scalegenesis in our model using the RG. First, we discuss a general condition to realize the scalegenesis in the flatland, which can be read from coefficients of beta functions at the one-loop level. We see that our model actually satisfies the condition, and then we discuss how the physical vacuum and masses are defined.

A. Condition for scalegenesis in flatland scenario
We start by looking at a general condition to realize the flatland scenario by following the literature [7]. For a system where a fermion and a scalar boson couple to each other and to a gauge boson, the RG equations for the gauge coupling g, the Yukawa coupling constant y and the quartic coupling constant λ at the one-loop level are given by where µ is a RG scale, a, · · · , f are coefficients depending on the degrees of freedom of fields and dots stand for irrelevant terms for the discussion below. As we have discussed in the subsection II B, the ratio between the Yukawa coupling and the gauge coupling, r = y/g, is crucial for realization of scalegenesis in the flatland scenario. Therefore, we rewrite the beta functions in terms of r = y/g such that where Let us here discuss a realization of a stable and finite vacuum generated from the Coleman-Weinberg mechanism in the flatland scenario. The effective scalar potential should be bounded for a large field value, ϕ ∼ M pl , to realize a stable vacuum. This requires that β λ < 0 for µ = M pl . The generation of a scale, here denoted by v, in the Coleman-Weinberg mechanism is realized by a negative quartic coupling constant at µ = v. For this, we need β λ > 0 at µ = v. From the beta function for λ in Eq. (29), these behavior can be achieved by r(v) < r 0 < r(M pl ). This condition also requires that r as a function of µ has to increase with increasing the scale. Thus, the generation of a scale in the flatland scenario could be realized when the condition r c < r(µ = v) < r 0 < r(M pl ) is satisfied. This condition can be expressed in terms of the coefficients in the beta functions as In our model (15), the Coleman-Weinberg mechanism works in the singlet-scalar sector, so that we can identify the coupling constants g, y, and λ with g X , y R (or y L ), and λ S , respectively. In this case, we have, for each beta function, the coefficients as where we assume that y L = y R . See Appendix A for explicit forms of the beta functions. Inserting these values into Eq. (31), we obtain K 0.625 and then can see that the condition (31) is satisfied. Note that the value of r(v S ) has to be between r c 1.04 and r 0 1.32. We also note that for y R = 0 one has b = 6 and d = 16, and then K 0.590.

B. Scalegenesis in flatland
We investigate the scalegenesis in the singlet-scalar sector owing to the Coleman-Weinberg mechanism. To this end, let us start with a discussion of the mechanism for scalegenesis in our model. As we have seen in Section III A, the ratios r L,R = y L,R /g X determine the scale of scalar field S. This scale is mediated through the Higgs portal coupling λ HS which is set to zero at the Planck scale. In order for the Higgs portal coupling to have a finite and negative value in low energy regimes, one needs a term not proportional to λ HS in its beta function. Indeed, as one can see from Eq. (A5) a term +(g mix g X ) 2 in the beta function plays such a role, so that a finite value of g mix is crucial for the inducement of the electroweak scale via the Higgs portal coupling. (See Eq. (45) below. ) We solve the RG equations for the system with the boundary conditions (3) and (5). In Fig. 1, we show an example of the RG flows for the quartic and the Higgs portal coupling constants, where the following initial benchmark value is used: The boundary condition for the gauge coupling constants for the the SM gauge fields is given in Appendix A 2. The Yukawa coupling constants y L,R , the U(1) X gauge coupling constant g X , and the kinetic mixing effect g mix are approximately constant within the RG flow since their beta functions are proportional to themselves. One can see that the quartic coupling of S is generated as a positive value in high energy region and then turns to a negative values at a certain renormalization scale. Such a behavior implies an occurrence of radiative symmetry breaking. The Higgs portal coupling constant is generated as a negative value. In order to extract information about the vacuum v S , we consider the effective potential for S. To this end, we here parametrize the complex scalar field S such that and then v S = √ 2 S = φ . Since the Higgs portal coupling constant is much smaller than the Yukawa coupling constant and the U(1) X gauge coupling constant, it could be negligible for the analysis of the vacuum. This treatment allows us to consider the effective potential only for φ. Thus, the improved effective potential [92] for φ is given by where λ S is the running coupling constant obeying its RG equation, the dimensionless RG scale is parametrized as t = ln(φ/M ) with a renormalization scale M . Here, we choose M = v S for which t = 0 corresponds to φ = v S . The effect of the field renormalization is with the anomalous dimension of S, The vacuum v S is obtained from the stationary condition, which gives a condition among coupling constants: where One can read the U(1) X breaking scale v S as a scale at which the condition (39) is satisfied. After the U(1) X symmetry breaking, the singlet-scalar obtains a finite positive mass squared, where we neglect the running effect of the field renormalization (36) in the second equality. Furthermore, the effective quartic coupling constant is obtained by for which the singlet-scalar mass squared (41) is Note that the effective cubic coupling constant is given by A negative Higgs mass parameter are generated through a negative Higgs portal coupling such that at the renormalization scale µ = v S , This mass parameter evolves until the electroweak scale v H by following its RG equation: where the anomalous dimension for the Higgs mass parameter γ m 2 H is given in Eq. (A6). Let us turn to the Higgs sector with the generated Higgs mass parameter (45). For a negligibly small Higgs portal coupling constant, the effective potential for the Higgs field is given by The electroweak vacuum is defined as the stationary point for the effective potential, dV eff /dh| h=v H = 0, for which one infers At this vacuum, the Higgs boson mass is given by where the last term on the right-hand side denotes the Higgs self-energy correction whose approximate form at the one-loop level is given by [4,93] ∆M 2 with This correction can be obtained by using the one-loop effective potential given in Appendix B and would be about 10% within the physical Higgs mass. So far, we have neglected the Higgs portal coupling. This may be a good approximation for obtaining the physical mass spectra as long as a Higgs portal coupling is small. Nevertheless, the mixing effect between the Higgs field and the singlet-scalar field plays a crucial role for the dark matter annihilation via these scalar fields. The quadratic terms for the (h, φ)-basis in the effective potential is diagonalized such that where we define M 2 HS = λ HS v H v S . The mass eigenstates are given by with the mass eigenvalues, and the mixing angle, For , the mixing angle is positive (negative). The mixing angle between the Higgs and a new scalar boson is constraint such that | sin θ| < 0.3 [94,95]. The Higgs mass in Eq. (54) has to satisfy the observed mass (22), i.e. M h = M obs H .
In the mass eigenstates (h , φ ), their propagators take forms where Γ H and Γ S are decay widths for the Higgs boson and the scalar field S. Using the mixing matrix (53), one obtains the propagators in the flavor basis (h, φ) so that For the benchmark value of the coupling constants (33), we obtain the expectation value of S, for which we observe

C. Decay of new particles
We note decay processes of the Higgs field and the singlet-scalar field. Due to the mixing between the Higgs field and the singlet-scalar field, the partial decay width of the singlet-scalar and the Higgs fields into the SM particles are given by The extra U(1) X boson decays into SM particles via the kinetic mixing effect with g mix > ∼ 10 −9 . This bound is adequately small for the extra U(1) X boson to decay in the current our system. Thus, the Majorana fermions annihilate into X µ if M R,L > M X .

D. Allowed parameter space
We scan the parameter space where the phenomenological constrains (22) and (23) are satisfied. We here summarize constraints for free parameters in our model. As can be seen in Eq. (22), the top-quark mass has a large uncertainty. We perform the parameter search with a fixed value M t 160.4 GeV for which the Higgs mass could be generated so as to M h 125 GeV within the SM. Although the extension of the SM could give a small deviation of the Higgs mass, it is still consistent within the uncertainty of the top-quark mass. From the asymptotic safety condition above the Planck scale, the gauge coupling constant should satisfy the bound g X (M pl ) < ∼ 1.09 and g mix (M pl ) < ∼ 0.68. To realize the scalegenesis in the flatland scenario, the ratio r = y L,R /g X should be in the range 1.04 < ∼ r(v S ) < ∼ 1.32 which gives a constraint for the Yukawa coupling constants by combing with the bound for the gauge coupling constant g X .
We first show the mixing angle (sin θ) as a function of g X (M pl ). Fig. 2 represents the region for sin θ where the constraints above are satisfied. The grey shadow region | sin θ| < 0.3 is excluded by collider experiments [95]. One can see from Fig. 2 that for g X (M pl ) > ∼ 0.84 the value of sin θ is positive, i.e., the singlet-scalar boson is lighter than the Higgs boson, and this region is already excluded. Therefore, hereafter we restrict the value of the U(1) X gauge coupling to g X (M pl ) < ∼ 0.84.
In Fig. 3 we show y L (M pL ) and g mix (M pL ) as functions of g X (M pl ). These couplings have a linear dependence on the U(1) X gauge coupling. In Fig. 3 we plot linear red lines which are given by, at the Planck scale, We parametrize the the right-handed Yukawa coupling as where ξ is a constant less than 1. With Eqs. (63) and (64), g X (M pl ) and ξ can be cast as free parameters in this system. One of them will be constrained such that the relic abundance of the Majorana fermions satisfy the dark matter relic abundance. Fig. 4 exhibits allowed region for dimensionful physical quantities (M L , M R , M X , M φ and v S ). One can see from Fig. 4 that these quantities tend to be in inverse proportion to g X (M pl ).

IV. MAJORANA FERMIONS AS DARK MATTER CANDIDATES
In this section we investigate the properties of the Majorana fermions as dark matters. We start by setting up the Boltzmann equation to evaluate the relic density of the Majorana fermions within the cosmological evolution, and then the allowed parameter region, where the the observed relic density is satisfied, is searched.

A. Boltzmann equation and dark matter relic density
There are two dark matter candidates, namely χ R and χ L . In order to follow evolutions of their number densities n R,L as functions of temperature, we here introduce the Boltzmann equations. Since the structure of the Boltzmann equations is symmetric under the exchange L ↔ R, we here show the case only for the left-handed side. Instead of the number densities, it is useful to introduce the quantities Y R,L = n R,L /s, where s is the entropy density. The Boltzmann equation for Y L is given by [96][97][98][99][100] dY L dx = −0.264g where M pl = 1.22 × 10 19 GeV is the Planck mass; g * = 106.75 is the total number of effective degrees of freedom in the SM; 1/µ RL = 1/M R + 1/M L is the reduced mass; x = µ RL /T is the dimensionless inverse temperature; and Y L is Y L in the thermal equilibrium, with K 2 (x) the modified Bessel function of the second kind. Here, σ(χ L χ L ; SMs, φ, X µ ) v is the thermal averaged cross section for the dark matter annihilation processes. Such annihilations take place with mediators, the singlet-scalar field φ [101][102][103], the Higgs field h, the U(1) X gauge field X µ [104] and the Majorana fermion as exhibited in Fig. 8 in Appendix C. As discussed in subsection III C, the singlet-scalar fields in the final state decay into lighter SM particles through the interaction with the Higgs field. These mediators cause also the χ L χ L → χ R χ R scattering of the Majorana fermions, whose thermal averaged cross section is denoted by σ(χ L χ L ; χ R χ R ) v . These processes are show in Fig. 9 in Appendix C where we show their explicit forms. Solving the coupled Boltzmann equation for Y R and Y L , one can evaluate the relic density of the dark matter, where we have s 0 = 2890 cm −3 and ρ c /ĥ 2 = 1.05 × 10 −5 GeV cm −3 [85,86], and Y R,L,∞ are the values of Y R,L at x = ∞ corresponding to the zero temperature. In our working assumption M R < M L (or equivalently y R < y L ), the annihilation of χ R s to χ L s does not take place, namely σ(χ R χ R ; χ L χ L ) v = 0. The left-handed Majorana fermions annihilate to the right-handed ones in addition to the SM particles and the singlet-scalar bosons within the temperature evolution so that the main ingredient of the dark matter relic density is the right-handed Majorana fermions. The left-hand side panel of Fig. 5 exhibits the relic abundance of the Majorana fermions (67) normalized by the observed one (23) as a function of g X (M pl ). There exists a region satisfying Ω DMĥ 2 = Ω obs DMĥ 2 in between 0.78 < ∼ g X (M pl ) < ∼ 0.81. The star points denote typical points for each value of ξ = y R (M pl )/y L (M pl ). We show the values (star points) of ξ as a function of g X (M pl ) in the right-hand side panel of Fig. 5. One can see the linear dependence of ξ on g X (M pl ). The linear fitting yields the relation for 0.78 < ∼ g X (M pl ) < ∼ 0.81, Then, all parameters excepts for g X (M pl ) are fixed by the observed data.

B. Prediction on Spin-independent elastic cross section
In the WIMP dark matter search, interactions between a nucleon and a dark matter play a crucial role for the detection of a dark matter signal. They could be observed as the spin-independent (SI) elastic cross section of a dark matter and a nucleon [99,105]. In our model, the scattering of the Majorana fermions and quarks could take place as the t-channel diagram in the processes displayed in Fig. 6 from which one can obtain the effective scalar-type four-Fermi interaction L eff = G R q (χ c R χ R )(qq) + G L q (χ c L χ L )(qq). More specifically, the coefficient of the four-Fermi interactions are calculated as with y q a quark Yukawa coupling constant. Note that since the quark-DM interaction induced by the X gauge boson exchange gives the spin-dependent cross section, we do not evaluate it in this work. The four-Fermi interaction (69) between a Majorana fermion and a quark is translated into the effective interactions between a Majorana fermion and a nucleon by where M q = y q v H / √ 2 is a quark mass, M N 939 MeV is the nucleon mass, and f N q = M q N |qq|N /M N is each quark matrix element of a nucleon. For the left-handed Majorana fermion χ L , the SI elastic cross section of a dark matter and a nucleon is computed as where µ N L = M N M L /(M N + M L ) is the reduced mass for the N -χ L system, and f N is calculated as with f N q evaluated in [106][107][108]. One can obtain the case for the left-handed Majorana fermion by exchanging L ↔ R.
In Fig. 7, we plot the spin-independent elastic cross sections with the upper bound provided by XENON1T [109]. The constraints from LUX [110] and PandaX-II [111] are somewhat milder than those of XENON1T; see [109]. One can see from Fig. 7 that there is a small allowed region slightly below the upper bound (solid black line).
The model has the strong predictability thanks to the conditions from the asymptotically safe quantum gravity scenario. If the spin-independent elastic cross section is observed, all parameters in the model are determined. The Majorana fermions as dark matter candidates in the model could be tested in the near future.

V. SUMMARY
We have proposed an extension of the SM based on the flatland scenario with dark matter candidates. The flatland condition corresponds to the fact that all scalar interactions including mass terms vanish at the Planck scale, namely the scalar potential is flat above the Planck scale, especially the model is scale invariant. Such a condition could be compatible with the asymptotic safety program of quantum gravity. We introduce Majorana fermions coupled to a U(1) X gauge field and a singletscalar field. The U(1) X gauge field interacts with the U(1) Y gauge field in the SM even at the classical level via the kinetic mixing effect. At this point, there are four free parameters (g X , g mix , y L and y R ). This is a minimal setup of an extended model compatible with the flatland conditions which could be naturally concluded from the asymptotically safe quantum gravity scenario. We have demonstrated that the minimal extension of the SM contains eventually only one free parameters and then has a strong predictability.
Let us here summarize the processes that the four parameters are fixed. The condition for the electroweak scalegenesis due to the Coleman-Weinberg mechanism gives a constraint for the ratio between the Yukawa coupling and the U X (1) gauge coupling. Hence, the elec- troweak scale v H = 246 GeV fix one of the Yukawa couplings. The kinetic mixing effect g mix generates a finite negative value of the Higgs portal coupling between the Higgs doublet-field in the SM and the singlet-scalar field, so that the observed Higgs mass m H = 125 GeV determines the value of g mix . We found y L (M pl ) and g mix as functions of g X by the numerical analysis, i.e. Eq. (63). The relic abundance of the Majorana fermions has to be satisfied the current observed value (23). From this constraint, we determined the ratio between the Yukawa couplings, denoted by ξ, as given in Eq. (68). Then, there is only one free parameter, e.g. g X (M pl ) in the model. We have evaluated the SI elastic cross section of Majorana fermions and a nucleon. The model predicts the SI elastic cross sections as functions of the Majorana masses around the current upper bound of XENON1T. There is a small allowed region slightly below the upper bound. Therefore, the Majorana fermions in the model as dark matter candidates could be tested by the direct detection experiments of the WIMP dark matter such as XENON, LUX and PandaX-II. If the SI elastic cross section is observed, all parameters in the model are determined. Hence, it is important to investigate the possibilities of the observations of the other particles, i.e. the singletscalar boson mass and the U(1) X gauge boson mass. The future collider experiments such as the High-Luminosity Large Hadron Collider (HL-LHC) [112] and the International Linear Collider (ILC) [113] could find these particles.
The investigation of stochastic gravitational waves produced by phase transitions at finite temperature may be one of other possible tests for the model. It is actually reported in [114] that a similar model (classically scale invariant B − L model) can produce gravitational waves whose spectra could be tested by future interferometer experiments. In such a case, one expects that a supercooling universe is realized. In particular, when the electroweak phase transition temperature is lowered until the QCD phase transitions, the thermal history of the universe could be changed drastically [115]. It is interesting subject to investigate the nature of our model in a supercooling universe.
Appendix A: Renormalization group equations

Beta functions
Here, we list the beta functions for coupling constants at the one-loop level in the extended model (15). The beta functions have similar structures to a B − L extension of the SM [116].
For the gauge coupling constants in the SM sector, one has while the beta functions for the gauge coupling constant for the new gauge field X µ and the kinetic mixing effect are given by The Yukawa coupling constants for the top-quark and the Majorana fermions run by obeying the beta functions The beta function for the left-handed Majorana Yukawa coupling constant are obtained by replacing R ↔ L.
For the scalar interactions, one has The anomalous dimensions for the scalar masses m 2 H and m 2 S are given by

Boundary condition for gauge coupling constants
We give the boundary condition for the SM gauge coupling constants. The value of the strong coupling constant g S is extracted from α S (M Z ) = g 2 S (M Z ) /4π = 0.1184 [85], where the Z-boson mass is M Z = 91.1876 GeV. One can obtain the values of SU(2) L and U(1) Y gauge coupling constants at M Z from the fine structure constant and the Weinberg angle [85] which are observed as = 127.916 , From these values, more explicitly, one can extract Here external momenta for the initial and the final states are expressed as, respectively, with s = (p A + p B ) 2 = (p a + p b ) 2 . Here, we assume the center-of-mass system (p A = (E, p) and p B = (E, −p)) and m A = m B ≡ m i . With the relative velocity v = 4 |p| 2 /s, the cross section for the two-body scattering process (C1) is The thermal averaged cross section is defined by with E A = E B = p 2 + m 2 i the energy dispersion. The momentum integrals (C4), however, cannot by evaluated analytically, so that, by assuming a small relative velocity, one expand the cross section into a polynomial of v 2 , i.e., where odd power terms of v are dropped since they vanish in the integrals (C4). We obtain the formula for the thermal averaged cross section,

Cross section for dark matter annihilation
Let us evaluate cross sections for each process exhibited in Fig. 8. Assuming that dark matters are non-relativistic, we expand the cross sections into a polynomial of the relative velocity v and take into account up to of order v 2 .
When two χ i s (i = L, R) annihilate into scalar fields, one finds σ(χ i χ i ; φφ) v = σ with σ (p) where λ effS 3 is the effective cubic coupling constant given in Eq. (44). Here, we define dimensionless functions, Possible annihilation processes of dark matters, which is denoted by σ(χLχL; SMs, φ, Xµ) v . These processes are represented in the flavor basis. The cross on the right-hand side diagram stands for the kinetic mixing between U (1)Y and U (1)X gauge fields. The cross in a circle stands for the mixing between φ and h. "SM" indicates contributions from the SM particles except for the Higgs boson, i.e., W and Z bosons and top-quarks, and f are SM fermions. Contributions from other quarks are neglected since their Yukawa coupling constant is relatively smaller than that of top quark, whereas for the decay process into fermions mediated by X boson, one has to take into account all quark contributions.
Annihilation processes of the light-handed Majorana fermions into the right-handed ones. This is denoted by σ(χLχL; χRχR) v .
As given in Eq. (57), the propagators in the s channel for the h-φ mixing and the singlet-scalar field φ are defined by with θ the mixing angle given in Eq. (55), and Γ H the decay width of H presented in Eq. (62), while the s-channel propagator of the X µ gauge field is given by For annihilation processes mediated by a Majorana fermion, we define its propagators in the t and u channels, with σ (s)

Thermal average cross section
Utilizing the formula (C6) the thermal averaged cross sections for dark matter annihilation to SM particles, φ and X µ are given by r iX a X (g X , M i ) + r ihZ a hZ (g X , M i ) + r iXφ a Xφ (g X , M i ) with 1/µ RL = 1/M R + 1/M L . Here, the five diagrams of the right-hand side in Fig. 8 give contributions b I so that a X (κ, M ) = 4g 2 X (κ) 2 ∆ X (M ) a hZ (κ, M ) = g 2 mix The χ L χ L → χ R χ R process is evaluated as