On the impact of Majorana masses in gravity-matter systems

We investigate the Higgs-Yukawa system with Majorana masses of a fermion within asymptotically safe quantum gravity. Using the functional renormalization group method we derive the beta functions of the Majorana masses and the Yukawa coupling constant and discuss the possibility of a non-trivial fixed point for the Yukawa coupling constant. For a certain value of the gravitational coupling constants and the Majorana masses, the Yukawa coupling constant has a non-trivial fixed point value and becomes an irrelevant parameter being thus a prediction of the theory. We also discuss consequences due to the Majorana mass terms to the running of the quartic coupling constant in the scalar sector.


I. INTRODUCTION
The formulation of a consistent theory which fundamentally describes the quantum fluctuations of gravitational degrees of freedom is still an open issue in theoretical physics. Different approaches to this problem based on different theoretical assumptions are being developed and each of them has internal challenges to be handled. Notably, the quantization of General Relativity within the standard perturbative quantum field-theoretic framework leads to a perturbatively non-renormalizable quantum theory. Among several different possibilities to circumvent this issue, Asymptotic Safety is one route which stays within the continuum quantum field theory realm, but goes beyond the standard perturbative paradigm [1][2][3][4][5][6][7][8][9][10].
In this approach, a fundamental quantum theory of gravity could be realized thanks to the existence of a non-trivial ultraviolet (UV) fixed point in the renormalization group flow. At this fixed point, the theory becomes scale invariant [11] and one is allowed to zoom in up to arbitrarily short distances without running into divergences. Being an interacting fixed point, the counting of relevant directions in theory space, i.e., the number of free parameters to be fixed by experiments is not controlled by the canonical dimensionality of the corresponding couplings. By integrating out quantum fluctuations, all terms compatible with the symmetries of the underlying theory will be generated in the quantum effective action yielding an infinite-dimensional theory space. For the theory to be predictive, the number of relevant parameters should be finite. Since one looks for an interacting fixed point, its discovery might not be feasible within perturbation theory and new techniques might be necessary. Thus, an asymptotically safe quantum theory of gravity relies on the existence of a non-trivial UV fixed point which features finitely many relevant directions.
The first evidence for the existence of the aforementioned non-trivial fixed point was found with the ǫ-expansion method in 2 + ǫ dimensions [1,12]. However, in order to establish the existence of the fixed point in four dimensions, other methods are necessary. One possibility is to regularize the path integral for quantum gravity by lattice methods and systematically look for a suitable fixed point which allows for a continuum limit which resembles our universe, see, e.g., [13,14]. In the asymptotic safety program, most evidence for the existence of a non-trivial fixed point were obtained with the functional renormalization group (FRG), see [15][16][17][18][19][20][21][22][23][24] for reviews. This was pioneered by Reuter in [25] where a UV fixed point was found using the FRG within the Einstein-Hilbert truncation for the effective average action, a scale-dependent action which takes into account the integration of quantum fluctuations up to some scale k. Several subsequent works, employing more sophisticated approximations provided compelling support for the existence of such the Reuter fixed point, see, e.g., . Remarkably, the existence of the UV fixed point has shown to be preserved against the introduction of Standard Model (SM) matter degrees of freedom, hinting to a fundamental theory of quantum gravity which is compatible with the observed SM particles, see [9,. Moreover, many different works have provided evidence for a saturation on the number of relevant directions associated to the Reuter fixed point, ensuring thus, predictivity [30,31,35,38,42,44,45,49,50,67,[79][80][81]. This underpins the quest for a quantum theory of gravity by using continuum quantum field theory methods.
Within the asymptotic safety scenario, coupling matter degrees of freedom to (quantum) gravity is straightforward. This allows for a rich interplay between the effects of matter fluctuations into gravitational running couplings and the reverse. In particular, this opens the possibility to test whether quantum gravitational effects allow for the resolution of long-standing problems as, e.g., the triviality problem in the SM, see [82][83][84]. Having a fundamental theory of gravity and matter valid up to arbitrarily short distances also leads to an enhancement of predictivity. In fact, the fixed point structure as well as the assumption that the SM holds up all the way until the Planck scale leads to a prediction of the Higgs mass to be m H = 126 GeV for the top quark mass m t = 171 GeV, see [78,85]. This happens thanks to its quantum-gravity induced irrelevance, being thus a prediction and not a free parameter of the theory. Recent works [84,86,87] have managed to successfully explain observed quantities in the low energy regimes by assuming the existence of the asymptotically safe fixed point for gravity-matter systems. Furthermore, asymptotically safe quantum gravity might be able to provide solutions for important problems in both particle physics and cosmology such as the gauge hierarchy problem [88] and the cosmological constant problem [89].
So far, the SM has shown to be consistent with results from the collider experiments up to the TeV scale, and no signal of new physics has appeared. This might be taken as a strong indication that, possibly, the SM is not drastically modified up to very high energy scales as the Planck scale. However, the SM does not account for dark matter. Many different models to extend the SM consistently with dark matter observations were developed in the last years. Very much like previously discussed, one can investigate if quantum-gravity effects play a relevant role in the description of dark matter phenomenology, see, e.g. [73]. This allows for low-energy consistency checks of the asymptotic safety scenario.
Another very important issue which is not accommodated in the SM is the (very tiny) neutrino mass. A standing out mechanism which describes the tiny mass for the left-handed neutrino while explains why a right-handed neutrino although necessary to be existent is not easily detected is the so-called seesaw mechanism [90][91][92][93]. Such a mechanism is grounded on the introduction of a right-handed neutrino with a Majorana mass term which couples to the lefthanded neutrino and the Higgs through a Yukawa interaction. Schematically, after spontaneous symmetry breaking, a Dirac mass term is generated and the balance between the Majorana mass parameter and the Dirac term allows for the assignment of a very large mass to the right-handed sector while the left-handed is very light.
In the present work, we aim at giving the first steps to investigate the quantum-gravity fluctuations effects to the seesaw mechanism. To pave the way for such understanding, we consider a toy Higgs-Yukawa system with one scalar and one Dirac fermion coupled to gravity. Different Majorana mass terms are introduced for the right-and left-handed components of the fermion. On top of that, right-and left-handed components interact with the scalar field via a Yukawa coupling. Beta functions for the Majorana masses and Yukawa coupling are computed taking into account quantum gravity effects with the FRG. Assuming the existence of a non-trivial fixed point for the gravitational couplings, we investigate the consequences to the masses and Yukawa coupling. In particular, regarding the critical exponent of the Yukawa coupling as a function of the Majorana masses and gravitational couplings we investigate the dependence of its sign on these parameters. We classify possible scenarios in low energy regimes as a consequence of asymptotically safe gravity. The impact of a non-trivial fixed point of the Yukawa coupling on the quartic scalar coupling is also discussed. Assuming that the system analyzed in this paper is associated with a simple toy model of the SM, the quartic scalar coupling at the Planck scale receives a finite correction from a fermion with the Majorana masses. This means that the prediction for the central value of the Higgs mass at the electroweak scale could be modified. The present work might be viewed as a stepping stone towards the more realistic system which takes into account all the SM properties. This paper is organized as follow: In Section II, we introduce the key concepts of the FRG and the toy model that we will investigate. Our results are collected in Section III. Section IV is dedicated to discussions and perspectives. Conventions and computational details are presented in the appendices.

A. Functional renormalization group
In order to extract the effects of quantum-gravity fluctuations to the running of masses and Yukawa coupling (as well as the impact of matter fluctuations to gravitational couplings), we use the FRG. Its central object is the effective average action Γ k . The parameter k acts as an infrared cutoff and Γ k is obtained by integrating out modes with momenta larger than k. Lowering k corresponds to integrate more and more momentum shells in the Wilsonian sense. The change of Γ k with k is described by the Wetterich equation [94], where ∂ t denotes the derivative with respect to t = log(k/k 0 ) (with a reference cutoff k 0 ), R k is a cutoff function, which suppress all modes with momenta smaller than k, (Γ (2) k + R k ) −1 is the regularized full propagator, and STr stands for the supertrace acting on all internal as well as spacetime indices and carries an appropriate numerical factor depending on the nature of the field. The Wetterich equation (1) has a one-loop exact form, a fact that greatly simplifies practical calculations. For some reviews about the Wetterich equation, we refer to [15][16][17][18][19][20][21][22][23][24].
In general, the effective action Γ k is expanded as, where g i are the dimensionful coupling constants and O i denote operators of the local fields φ. Using the Wetterich equation and expanding its right-hand side on the same {O i } basis, one is able to read off the beta functions for the couplings. They have the following structure, whereg i = k −di g i is the dimensionless coupling constant and d i is the canonical dimension of g i . The first term on the right-hand side of the equation (3) is the canonical scaling term, associated with the dimensionality of the coupling and F i (g) corresponds to quantum corrections. Here, let us assume that the system has a fixed pointg * = {g * 1 ,g * 2 , . . .} at which the beta functions vanish: β i (g * ) = 0 for all i. We linearize the RG equation (3) around the fixed point, These equations can be diagonalized and trivially solved. The eigenvalues of T , here denoted by θ i , correspond to the critical exponents associated to the fixed point. In our notation, positive θ i are relevant directions in the RG sense. The number of relevant directions dictates how many free parameters the theory features. That is, a theory with less relevant directions is more predictive. For the Gaussian (trivial) fixed pointg * = 0 the critical exponents are given by θ i ≃ d i and the coupling constants with d i > 0, i.e., the power-counting renormalizable ones, are relevant. In contrast, if the theory has a non-trivial fixed pointg * = 0, the critical exponents are θ i ≃ d i − ∂F/∂g i |g =g * . We see that the contribution from the loop correction ∂F/∂g i |g =g * typically does not vanish at a non-trivial fixed point leading to a non-canonical scaling. Such non-trivial dynamics can drive canonically irrelevant/relevant coupling into relevant/irrelevant ones. In the case where relevant couplings turn to irrelevant ones, the predictivity power of the theory is enhanced.

B. Setting the truncation
Our goal is to compute the system of beta functions for a Higgs-Yukawa model with a Majorana mass term coupled to gravity. In order to apply the FRG, we make use of the background field method. In practice, we linearly split the metric 1 into a background metricḡ µν and fluctuation h µν , Hereafter indices are raised and lowered by the background (inverse) metric. We investigate the following truncated effective average action, where the gravity sector is given by Here, G and Λ, respectively, correspond to the Newton and cosmological constants, a and b denote higher curvature couplings, and S gf and S gh are the gauge fixing and Faddeev-Popov ghosts actions 2 written as with F µ :=∇ ν h νµ − β+1 4∇ µ h being the standard gauge condition. Here,∇ µ denotes the covariant derivative with respect to the background,¯ =∇ 2 is the d'Alembertian. In addition, h :=ḡ µν h µν is the trace of h µν ; C µ andC µ are the ghost and anti-ghost fields, respectively; and α and β are dimensionless gauge parameters. For the matter sector, we consider the following action, where V φ 2 is the scalar potential with Z 2 (φ → −φ) symmetry and Z φ,ψ are field renormalization factors of scalar and fermion fields, respectively, and contribute to the beta functions through the anomalous dimensions expressed as Here, / D = γ µ D µ is the Dirac operator and satisfies − / D 2 = −D 2 +R/4; ψ c is the charge conjugated spinor field, namely, ψ c = Cψ T ; P L = (1 − γ 5 )/2, P R = (1 + γ 5 )/2 are the chiral projection operators. We refer to Appendix A for further conventions. The effective action (6) has the discrete γ 5 -symmetry, i.e., ψ → γ 5 ψ,ψ → −ψγ 5 and φ → −φ, and then the Dirac mass termψψ and the Majorana-Yukawa terms φψ c P R,L ψ + h.c. are forbidden, whereas the Majorana mass term is not. Indeed, even if we put m L M = 0 and m R M = 0 at some scale, quantum effects induce a finite m L M as we will see in (26). We also impose the CP invariance which prohibits the term φψγ 5 ψ. The CP invariance also makes the parameters m R m and m L m to be real ones. It follows that the Hermitian conjugate of the Majorana mass term is m R Mψ P L ψ c + m L Mψ P R ψ c . In this work, since we employ a simple model where a neutral fermion interacts with a scalar-field singlet, both the left-and right-handed fermion fields can have Majorana mass terms. Note however that the left-handed neutrino is a component in the lepton-doublet field involving a hypercharged particle in the SM, so that the Majorana mass term for the doublet field cannot be admitted in the SM before the electroweak symmetry breaking.

C. York decomposition
For the fluctuation field h µν and the ghost fields, the York decomposition [99] is employed, and where h TT is transverse and traceless (TT) tensor field with spin 2, i.e.,∇ µ h TT µν = 0 andḡ µν h TT µν = 0; ξ µ is the transverse vector field with spin 1, i.e.,∇ µ ξ µ = 0; and σ and h :=ḡ µν h µν are the scalar fields with spin 0, C andC are spin-0 scalar ghost fields, and C T µ andC T µ are spin-1 transverse vector ghost fields that satisfy∇ µ C T µ =∇ µCT µ = 0. Such decompositions entail the following Jacobian in the path integral measure: where the primes indicate that the first or the first two eigenvalues are removed from the determinant. Although this Jacobian can be absorbed by redefining the fields ξ µ , σ,C and C, here we take these contributions into account in the flow equation without any field redefinitions. This can be done by regularizing the Jacobian with the usual prescription −¯ → P k (−¯ ) = −¯ + R k (−¯ ). In this framework, we can easily verify that the Jacobian gives the following contributions to the flow equation

D. Computing beta functions
In this subsection, we present the basic strategy to calculate the beta functions. In the present work we use the gauge choice β = 0 and α → 0, for which a considerable simplification in the flow equations happens. As discussed in Appendix B, for this gauge choice, the σ mode in the York decomposition (11) is decoupled from the mixing with the other spin-0 scalar modes (h and ϕ -with ϕ being the quantum fluctuation of the scalar field about a fixed background configuration). For the system (6) with the York decompositions (11), (12) and this gauge choice, the Wetterich equation (1) reads where the first term on the right-hand side is the contribution from the TT graviton mode; the subscription "SS" denotes contributions from ϕ and h; the contributions from the fermion-loop corresponds to the third term; ∆T (1) and ∆T (0) are the spin-1 vectors (ξ µ , C ⊥ µ and the spin-1 contributions from the Jacobians) and the spin-0 scalars (σ, C and the spin-0 contributions from the Jacobians), respectively; and the last two terms are contributions from mixings between bosonic and fermionic degrees of freedom. The quantities G ΦiΦj represent the regularized full propagator of the fields Φ i and Φ j . For more details on the structure of the flow equation, see Appendix C. In what follows, it is convenient to rewrite the flow equation as where the contributions T 's are given by Their derivations and the explicit form of ∆T (1T) and ∆T (0) are shown in Appendix C. Here, the subscripts on the trace stand for the spin of the corresponding fluctuation. For explicit computations, the Hessians associated with the present truncation, which are necessary to compute the right-hand side of (16), are listed in Appendix B. The results presented in this paper are computed in the background approximation, where we set the metric fluctuation to zero in the flow equation after computing the Hessians. In this work we used Litim's optimized profile function (see Appendix B) [100]. We use the Type-I optimized cutoff (z = −¯ ) in order to compute T (2) , T (1) and T (0) , while for the spin 1/2, we employ the Type-II choice (z = −¯ +R/4 ≡ − / D 2 ) in order to have the correct sign for the fermionic contributions to the Newton coupling constant [101]. The computation of T mixed is performed around flat background, in such a case both the Type-I and -II cutoff choices collapses to z = −∂ 2 .
In order to extract beta functions from the Wetterich equation, we can project (16) in different sectors of interest. For the gravitational couplings, the projection consists in setting all matter fields to zero. In such a case, we can easily verify that T mixed does not give any contribution. Substituting our truncation in the left-hand side of (16), we find where In terms of the dimensionless quantities, the last equation can be rewritten in the following way In addition, we define the beta function of the dimensionless couplings as β Λ = ∂ tΛ , β G = ∂ tG , β a = ∂ tã and β b = ∂ tb . By expanding the left-hand side of the above equation eq. (19) in powers of curvature invariants, we can compute beta functions by matching the coefficients of the same invariants.
In order to project the flow equation to the scalar field potential, we set the curvature terms and the fermionic fields to zero. In this case, we only have contributions coming from T (2TT) , T (0) and T (1/2) . Then, the flow equation takes the form Defining the renormalized fieldφ = Z 1/2 φ φ, where Z φ denotes the wave function renormalization, the flow equation for the scalar potential can be recast in the following way whereṼ (φ 2 ) = V (φ 2 )/k 4 denotes the dimensionless scalar potential. Expanding the dimensionful scalar potential as 2n , where we have defined the dimensionless couplings asλ 2n = k 2n−4 Z −n φ λ 2n and the renormalized scalar field isφ = Z 1/2 φ φ. By expanding both sides of (21) as a power series inφ 2 , we can extract beta functions for the dimensionless coupling constants.
Finally, in order to compute beta functions for the Majorana masses and the Yukawa coupling, we set the curvature terms to zero and the matter fields to constant configurations. Also, we discard all contributions of order O(φ 2 ) or higher. Implementing such projections, the beta function associated with the Majorana Right (Left) mass can be extracted by keeping only those terms proportional toψ c P R(L) ψ + h.c., while the running of the Yukawa coupling can read from those terms proportional to φψψ. We have and In terms of renormalized fieldsφ = Z and where A 1,L , A 3,L and M y 2 ,L are computed in terms of threshold functions (see Appendix D). The beta function for the right-handed Majorana mass can be obtained by exchanging L ↔ R in the above expression. For the sake of simplicity, we will not take into account effects coming from the fermionic anomalous dimension in our discussion. From the perspective of asymptotically safe quantum gravity our goal is to understand whether the running of the Majorana masses is compatible with an UV complete setting. In this framework, we have some interesting possibilities: Let us focus on the first scenario. In this case it is straightforward to see thatm R M =m L M = 0 is a fixed point of the subsystem β m L M = 0 and β m R M = 0, irrespective of the values of the gravitational couplings. In the absence of gravity the critical exponents associated with this fixed point correspond to the canonical scaling dimension of the Majorana masses and, therefore, are positive. If we take quantum gravity fluctuations into account, the critical exponents associated with Majorana masses are given by where Moreover,G * and y * D correspond, respectively, to the values of the dimensionless Newton's constant and the Yukawa coupling at the fixed point. The relevant question here is whether quantum gravity contributions could change the sign of the critical exponents, driving (one of) them towards irrelevance. In order to answer this question we first note that, at the fixed point withm R M =m L M = 0, we can express A 1 = D * y (see Appendix D for the definition of D y ). As we shall see in the next section, the viability of UV completion in the Yukawa sector requires D * y < 0. We first explore a scenario with y * D = 0. In this case the critical exponents for the Majorana masses become In the case of a fixed point with non-vanishing Yukawa coupling we have y * D = −8π 2G * D * y (see section III B) and, as a consequence, we find the following critical exponents AssumingG * > 0, the requirement D * y < 0, which is necessary for UV completion in the Yukawa sector, restricts the space of viable fixed point values in such a way that the critical exponents for the Majorana masses (at the fixed point value characterizedm R M =m L M = 0) remains positive after the inclusion of quantum gravity effects. This is the simplest (non-trivial) UV complete extension of a Higgs-Yukawa model incorporating Majorana masses. If θ M 1,2 < 0 due to quantum gravity fluctuations, the flow of the Majorana masses would be trivial. In summary, for the possibility (i), demanding a UV-complete Yukawa sector necessarily implies that the Majorana masses will be relevant parameters even at non-trivial values for the other couplings fixed-point values.
The analysis regarding the possibility of the fixed points (ii) and (iii) is more complicated. In this case the fixed point equations are more involved due to non-polynomial dependence on the Majorana masses and, as consequence, there is no simple criteria to probe the existence of such fixed points analytically. In the Appendix E we have shown the possibility of finding numerical solutions corresponding to fixed points in these classes. Instead of looking for further numerical solutions, in the next section we assume the existence of fixed point solutions featuringm R M = 0 and/orm L M = 0 and investigate their consequences to the viability of UV completion for the Higgs-Yukawa sector as well as in the predictability of the Higgs mass in the framework of asymptotically safe quantum gravity.

B. Impact of Majorana masses on the Higgs-Yukawa couplings
A fruitful strategy to probe quantum gravity induced effects in the UV regime of the matter sector is to treat the fixed points values of the gravitational couplings as free parameters in the beta functions associated with matter couplings. This approach allows us to explore the viability of quantum gravity induced UV completion and predictive RG trajectories for the matter sector [69]. Although the computation of beta functions associated with the matter sector involves a given truncation for the effective average action, this strategy gives information beyond a specific choice of truncation, since the values for the gravitational coupling at the fixed points are not restricted by the chosen approximation. Nevertheless, in Appendix E, we show results for fixed point values and critical exponents in a maximally symmetric background. In the present section, we adopt the same reasoning in order to investigate the impact of Majorana masses on Yukawa-Higgs systems in the UV regime. Moreover, we treat the fixed point values of the Majorana masses as free parameters and analyze their impact to the beta functions of the Yukawa and λ 4 φ 4 interactions.
We start by discussing the fixed point scenarios for the Yukawa sector. As discussed in the recent literature, quantum gravity fluctuations play an important role in the running of the Yukawa coupling, leading to an interesting phenomenological scenario with the prediction of the top quark mass [102], as well as, top-bottom mass difference [87]. The beta function of the Yukawa coupling has the following structure, where we define the D y 's and M y 3 in terms of threshold functions (see Appendix D). In order to simplify our analysis, we consider some additional approximations: i) we neglect contributions coming from the anomalous dimension; ii) we assume the Gaussian fixed point value forλ 2 . Taking these approximations into account, the beta function for the Yukawa coupling constant reads Following the discussion presented in reference [69], there are two possibilities for quantum gravity induced UV completion in the Yukawa sector. The first one is the Yukawa-free fixed point, i.e.,ỹ * D = 0. In this case, the critical exponent associated with the Yukawa coupling constant is given by where the notation (· · · ) * indicates that the stability matrix is evaluated at the fixed point. AssumingG * > 0, we conclude that the Yukawa coupling becomes a relevant direction (θ y > 0) provided that 3D * y < 0. This is the case of asymptotically free UV completion for the Yukawa sector.
The second scenario for Yukawa UV completion is characterized by an interacting fixed point induced by quantum gravity contributions, namelyỹ * The existence of such a fixed point is subject, of course, to the positivity of the square root argument. For small values of the Majorana masses at the fixed point, as we shall see later, it is possible to verify that M * y 3 > 0 and, therefore, the condition for existence of this fixed point is given byD * y < 0. More generally, the viability condition for such a scenario is given byD * y /M * y 3 < 0. The critical exponent associated with this fixed point reads If M * y 3 > 0, the existence of the fixed point (33) is subject toD * y < 0, implying θ y < 0 (UV repulsive fixed point). Such a scenario is particularly interesting since the requirement of UV completion allows us to make predictions in the IR regime [69]. On the other hand, if M * y 3 < 0, which can achieved for large values of the Majorana masses at the fixed point, the viability condition requiresD * y > 0, resulting in θ y > 0 (UV attractive fixed point), so that the Yukawa coupling remains a free parameter.
It is interesting to analyze the behavior ofD y and M y 3 as a function of the Majorana masses and the gravitational parameters. For the sake of simplicity, let us start with the TT-approximation, namely, we take only the TT graviton effects into account. In this case we haveD The M y 3 term, which only receives contributions from matter fluctuations, is given by Within the TT-approximation, the sign ofD y is completely determined by the value of the gravitational couplingb, namely,D y < 0 forb < −2/3 andD y > 0 forb > −2/3. The sign of M y 3 as a functions of the Majorana masses is depicted in Fig. 1  In order to complement the analysis for non-vanishing cosmological constant, in Fig. 4, we present the resulting sign ofD y /M y 3 as a functions ofΛ andb for some choices of the Majorana masses. These plots were generated with the choiceã = −0.5 in order to avoid regions close to the pole associated with the spin-0 sector of the graviton propagator. This choice can be justified by observing that far away from the spin-0 pole the sign ofD y /M y 3 is not affected by the variation of the parameterã.
The collection of plots presented in this section allows us to extract important information about the viability of 1/2, then we find a remarkable modification. In this case the viability regions associated with each one of the scenarios for UV completion discussed here becomes complementary. In addition, the interacting fixed point turns out to be UV attractive, in contrast to what happens in the absence of Majorana masses.
An important achievement of the asymptotic safety program for quantum gravity was the prediction of the Higgsboson mass around the measured value before the discovery of the Higgs boson [103,104]. Such a result relies on the existence of an approximated Gaussian UV repulsive fixed point for the quartic scalar couplingλ 4 . Beyond the Planck scale quantum gravity dominates over the remaining interaction, therefore, in this regime the evolution ofλ 4 is driven by running of the Newton coupling. The leading quantum gravity contribution to beta function associated withλ 4 takes the form β grav.

λ4
= D λ4Gλ4 , where D λ4 depends on the remaining gravitational parameters. For an approximated Gaussian UV repulsive fixed point the condition D λ4 > 0 is required at the fixed point. Several computations indicate that such a condition can be realized for simple truncations of the effective average action [52-54, 60, 78, 105]. Below the Planck scale, on the other hand, gravity fluctuations become irrelevant and the running ofλ 4 is driven by the Yukawa coupling. The contribution coming from the Yukawa sector to β λ4 is required to be negative in order to avoid instability of the quartic potential. In the framework of the SM coupled to gravity, if we set the Yukawa in accordance with the observed top-quark mass, then, the RG flow drivesλ 4 to an IR value compatible with the observed Higgs mass [85].
In the framework of our toy model, we can investigate the impact of Majorana masses in a predictive scenario for the Higgs boson mass. In such a case, the beta function of the quartic coupling is given by As we have mentioned before, a predictive scenario for the Higgs mass is subject to two conditions: i) D λ4 has to be positive, which is necessary for a viable (approximated) Gaussian UV repulsive fixed point; ii) The contribution coming from the Yukawa sector has to be negative in order to avoid instabilities in the quartic potential (at least in a first approximation of the stability of the Higgs potential). The first condition is not affected by the inclusion of the Majorana masses, since D λ4 depends only on the gravitational couplings. The coefficient of theỹ 4 D term, on the other hand, is a function of the Majorana masses and, as a consequence, the sign of such a contribution may change by varying these parameters. Accordingly, within our truncation M y 4 is given by The sign of M y 4 is dictated by the contribution in the brackets, which is not positive definite. In Although our results are obtained from a simple toy model, we can extract some physical consequences that, in principle, can be extrapolated for more realistic scenarios. Let us start by restricting ourselves to the framework of our toy model. In this case, the viability for a predictive (toy-)Higgs boson mass can be restricted by the condition M y 4 < 0. Such a condition restricts the space of allowed values for the Majorana masses which are compatible with a viable predictive scenario. In particular, the green region (region II) in Fig. 5 could drive the quartic coupling to an unstable regime.
In a more realistic scenario we should include further fermionic degrees of freedom (and also gauge fields, which are not relevant for this qualitative discussion). Let us assume that such additional fermions do not introduce new Majorana masses, but include extra Yukawa-Higgs interactions. For the sake of simplicity let us consider only one additional fermion associated with a (renormalized) Yukawa coupling denoted byỹ ′ . In this case, the beta function for λ 4 receives an addition contribution like − 1 8π 2ỹ ′ 4 . In such a case, the viability for a predictive Higgs mass is subject to the following condition

IV. SUMMARY
In this paper we have studied, using FRG techniques, a quantum-gravity matter system including a simple Higgs-Yukawa sector and Majorana masses for the fermions. We have derived the beta functions for the matter couplings and discussed their basic structures. Assuming the asymptotic safety scenario for quantum gravity, i.e., the running of the gravitational couplings flow towards a non-Gaussian fixed point in the UV, we investigated effects of graviton fluctuations on the running of the matter couplings. Our main findings can be summarized as follows: • Quantum gravity effects on the running of the Majorana masses: -The inclusion of Majorana masses can affect the sign of the coefficient of theỹ 4 D -term and, as consequence, it may have effects on the prediction of the Higgs from asymptotically safe quantum gravity.
The system considered in this paper can be regarded as a toy model motivated by neutrino physics. An interesting feature in a system with Majorana masses is the possibility of a mass hierarchy due to the so-called seesaw mechanism. If we assume that the neutrino is the fermion with the Majorana mass, from the fact that a neutrino mass is of order approximately 0.1 eV, one needs a huge hierarchy between the Dirac mass and the right-handed Majorana mass, GeV is the electroweak vacuum. Note that in the framework of the Standard Model the left-handed neutrino mass is forbidden since it has a SU(2) L ×U(1) Y charge. The present analysis could be a first step to understand the compatibility of neutrinos mass mechanism in the asymptotically safe gravity scenario. Our finding is that for a certain value of the gravitational coupling constants and the Majorana masses, there exists regions where the the Yukawa coupling constant has a non-trivial fixed point. This indicates that within this scenario, one could identify the scale of the Majorana mass from the irrelevance of the Yukawa coupling and the observed neutrino mass. Therefore, the possibility of future observations of the Majorana mass of neutrino could lead to consistency tests for asymptotically safe gravity.
The Yukawa coupling between a scalar boson and a fermion gives corrections to the quartic scalar coupling constant. Through the fermion loop, the Majorana masses affect the quartic coupling constant as we have seen in (37) and (38). By associating the scalar boson and the fermion with the Higgs boson and a neutrino, respectively, the predicted value of the Higgs boson mass from the asymptotically safe gravity scenario could be modified due to the running of the Majorana masses. It depends on the values of the Majorana masses and whether or not the neutrino loop makes the non-trivial fixed point of the quartic coupling constant of the Higgs boson. The current prediction from asymptotically safe quantum gravity for the Higgs boson mass is approximately 129 GeV, by taking up to three-loop effects on the RG equations into account within the SM [106]. For a certain value of the Majorana masses, the central value of a predicted Higgs boson mass could be modified towards a smaller value. Consequently, by taking into account effects of possible Majorana masses within the neutrino sector, the asymptotic safety scenario for quantum gravity might predict the observed Higgs boson mass 125 GeV.
In a future work, we will extend our analysis to a more realistic model in connection with neutrino physics. In this sense, the idea is to perform a FRG analysis within an approximation compatible with Standard Model (including Majorana masses) coupled to gravity. We will study the possibilities for UV completion on the neutrino sector and investigate the compatibility of the seesaw mechanism with the asymptotic safety scenario for quantum gravity. Depending on the structure of an UV fixed point obtained within a theory space including Majorana masses, it might be possible to predict IR values of the neutrinos Yukawa couplings and/or Majorana masses as a function of the gravitational couplings. Within this scenario, one could directly compare results from the asymptotically safe gravity scenario with the current experimental data of neutrino physics. Together with the predictions for the Higgs and top-quark masses and so on [84][85][86][87][88]102], one can further test the consistency of the asymptotic safety scenario for quantum gravity with the observed matter degrees of freedom of our universe.

Dirac gamma matrix and Clifford algebra
The standard formulation of spinor in curved spacetime is given by the verbein e a µ such that a metric is defined by g µν = e a µ e b ν η ab , where η ab = diag (1, −1, −1, −1) is the metric in flat spacetime. In this work, we use the spin-base invariant formalism [107,108] instead of the verbein formalism. Hereafter, following the literatures [109,110] we briefly summarize the spin-base invariant formalism. We start with demanding that the Clifford algebra for the Dirac gamma matrices in curved spacetime is given by where 1 4×4 is the unity matrix in spinor space. We require the invariance of the action for fermions under the following spin-base transformations with S ∈ SL(4, C): Here, Dirac conjugation for the Dirac spinor ψ is defined byψ where h is a spin metric satisfying From the spin-base transformations (i)-(iii) and the definition (A2) we see that the spin metric is transformed under the spin-base transformations as h → (S † ) −1 hS −1 . In flat spacetime, one choose usually the representation h = γ 0 . In general, however, the spin metric is not equal to γ 0 . Nevertheless, it is not necessary to specify the explicit representation of h for the evaluation of the Hessian derived in the next section. The second property of the the spin metric in (A3) are derived from an assumption that the mass operator of fermion is real , i.e., (ψψ) * =ψψ. Using this property, we see that We next define the chiral and charge conjugation operators using the Dirac gamma matrices. The chiral operator is given by whereε µ1...µ4 = √ −gε µ1...µ4 is the Levi-Civita tensor with ε 4123 = 1. This matrix satisfies The charge conjugation operator is defined by which satisfies We introduce the covariant derivative where Γ µ is the affine spin connection, whose spin-base transformations are give by Note that the covariant operator satisfies the following conditions: (i) linearity: When acting the covariant derivative on the Dirac gamma matrices and the spin metric, we have With the spin metric one can define the Dirac conjugate of a matrix M such thatM = h −1 Mh, which implies (ψMψ) * =ψMψ. The assumption that the kinetic term of fermion is real, namely, d 4 x and whereΓ µ is an auxiliary matrix introduced in order to define the following spacetime covariant derivative with ν µκ the affine connection. The conditions (A13) are obviously satisfied by identifying Γ µ =Γ µ . Indeed, as discussed in [109], ∆Γ µ is interpreted as the spin torsion such that In this work, we do not consider the effect of the torsion, that is, we set ∆Γ µ = 0. Finally, we mention our convention the Wick transformation from Minkowski spacetime to Euclidean spacetime. For a vector V µ ∋ (x µ , A µ , · · · ), we define where the subscript "E" denotes Euclidean, and i = 1, .., 4. The Dirac gamma matrices in Euclidean spacetime γ Eµ are defined by With this definition, the Clifford algebra in flat spacetime is given as It should be noted that the important properties of the operators (A6) and (A8) are not changed under the Wick transformation.

Variation
In order to calculate the Hessian for the fermion we need the formulae of the variations for the spinor field. Let us use the local spin-based formalism with the Lorentz-symmetric gauge [109,110].
We consider the variation of the Dirac matrices: whereγ µ = γ µ (ḡ). From the Weldon theorem [108], the second term on the right-hand side is with δ αβ µν = 1 2 (δ α µ δ β ν + δ α ν δ β µ ). Here, G αβ is a traceless tensor. The Lorentz-symmetric gauge corresponds to a choice G αβ (ḡ) = 0. In this gauge the variations of the Dirac matrix γ µ (g) are evaluated as In the same manner, using the property (A14) the variation of the affine spin connection iŝ Then we obtain that the variation of the covariant derivative D µ so that In the Lorentz-symmetric gauge, from (A3), (A6) and (A8), one can show that h(ḡ Since the Hessians are of order of h 2 , the variations of these operators do not contribute.

Appendix B: Hessians
In this appendix we list the Hessians computed from our truncation for the effective action given by (7) and (9). Expanding Γ k up to second order in the fluctuation fields (both for gravity and matter sectors), namely and, keeping only those terms quadratic in the fluctuations, we can write where we defined the fluctuation fields, Hereafter the background fields of the scalar boson and fermion are denoted by φ and ψ. The bosonic Hessian is with Γ (2) For the gravitational sector, after expansion in terms of the metric fluctuation, we perform the field redefinition 32πG h µν in order to incorporate the wave function renormalization for the graviton. Such a multiplicative factor may be absorbed in the Hessians, resulting in the following explicit expressions 4 where we have defined and we used the identity 1 µνρσ = (ḡ µρḡνσ +ḡ µσḡνρ )/2. Also, we used the shorthand notations, We also have used the following Lichnerowicz Laplacians, Hereafter, the prime on V denotes the derivative with respect to φ 2 . The Hessians for the terms with the scalar field φ are Note here that for β = 0, the dependence of Γ with the 2 × 2 matrix for the (h, ϕ) part given by and the (σσ) component For α → 0 the Hessian of the transverse vector mode (B6b) Obviously, the ξ µ and σ modes depend on neither gravitational couplings nor the matter term I M , so that their loop effects do not contributie to the beta fucntions of the Yukawa coupling and the Majorana masses. The Hessians for the ghost fields are also given as gravitational-coupling-and matter-independent forms: We show the Hessian with the fermion fields. By using (A21) and (A23) we obtain the two-point function The fermion propagator is, then where with A R,L and B being defined as Using (A21) and (A23), the mixing parts are computed as where ∆ C = C(m R M P R + m L M P L ) and∆ C = (m R M P L + m L M P R )C. Note that since the gamma matrices in the projection operators and the CP transformation are given as ones in the local Lorentz spacetime, their variations vanish.
We employ the Litim-type cutoff function [100] by means of the following regulator We used the Type-I and II cutoff functions, respectively, for the bosonic and fermionic propagators. For the bosonic cutoff we can express while the fermionic one is given by where the shape function is given by The Wetterich equation for a system composed by bosons and fermions can be written in the following way Exploring the block-diagonal structure of the cut-off function (i.e., there is no contribution like R BΨ k or R BΨ k ), we can expand the traces as follows It is important to emphasize that the first two terms of these expansion are sufficient for the purposes of this paper, therefore, we can truncate expansion without any loss. After York decomposition, the first trace in the bosonic sector can be expressed as 5 where G ij = ( Γ for every pair (i, j). It is convenient to define a short notation where T (s) is defined by the sum of all contributions with spin-s. As discussed in the previous section, a remarkable simplification is achieved by the gauge choice α → 0 and β = 0. In this case, the contributions coming from the spin-1 and spin-0 can be summarized in the following expressions and where Note that SS denotes ϕ and h. The aforementioned choice for the gauge parameter also simplifies those contributions to the trace involving fermions-boson mixing. In fact, it is not difficult to conclude that, in this case, the contribution coming from the bosonic sector are those associated with the spin-0 fields h and ϕ (i.e., from the sub-sector Γ (2) SS ). Therefore, the mixed contributions to the flow equation can be cast in the following way For the pure fermionic trace we denote Putting all pieces together, we can rewrite the flow equation in the compact form

Appendix D: Beta functions
We list the explicit forms of the beta functions in four dimensional spacetime. To this end, we define a threshold functions: We define the graviton anomalous dimensions as Before performing the explicit evaluation of the beta functions, we analytically investigate their structure in order to gain some insights at particular limits. We start by looking at the structure of the spin 2 and 0 graviton propagators 6 , namely Around the poles of these propagators the graviton fluctuations become strong. When the fixed point value lies near the poles, the critical exponent tends to be large. However, such a case is not acceptable as a reliable result since before one is dealing with very stable and convergent truncations, the fixed point values and thus the proximity to the propagator poles is truncation-dependent. This prevents robust conclusions due to the artificial enhancement of gravitational fluctuations. The dimensionless beta function of the left-handed Majorana mass is where 6 These expressions are obtained by expanding the full two-point function Γ in powers ofR and retaining theR-independent part.
Each contribution to the left-handed Majorana mass corresponds to the diagrams shown in Fig. 6. The right-handed one is obtained by the replacingm L M ↔m R M . The dimensionless beta function of the Yukawa coupling constant is β y = η φ 2 + η ψ ỹ D + (I) + · · · + (XXV) + (XXVI) + (XXVII) , where The diagrams of these contributions are shown in Fig. 7.
For the discussion presented section III B we need the contribution of Majorana masses to running of the scalar potential V (φ 2 ). Such a contribution comes exclusively from the fermionic loop corresponding to T (1/2) . Keeping only terms up to order O(R) we find   We kept the contribution proportional toR since this it contributes to the running of the Newton coupling, which we used in order to compute the fixed points in the appendix E. The dots denote corrections to higher dimensional operators which are truncated in the present analysis such as R 2 , F (φ)R and F (φ)R 2 . Contributions to the scalar potential coming from graviton and scalar loops was evaluated in e.g. [53,60].

Appendix E: Fixed point analysis
In this section we present the fixed point structure obtained from our truncation. The idea is to establish a benchmark analysis and investigate the possibility of non-trivial fixed points for the Majorana masses. In this sense, we restrict ourselves to the Einstein-Hilbert truncation for the gravitational sector. Further investigations are necessary in order to check the stability of our results by enlarging our truncation. The results presented in this section were obtained for a maximally symmetric background. In such a case we havē R µναβ = 1 12 (ḡ µαḡνβ −ḡ µβḡνα )R,R µν = 1 4ḡ µνR .
For the scalar potential we consider the φ 4 -truncatioñ In the present analysis we neglect any contributions coming from anomalous dimensions. In addition, the contribution of matter fluctuations on the running of the gravitational couplings was only considered by coming from the kinetic terms of the matter fields. These can be faced as additional approximations in our calculation and should be improved in future investigations. Within this framework we found three classes of fixed points. Table I shows the numerical results for the fixed points, while in table II we report the corresponding critical exponents. The first class correspond to the case of Gaussian matter. In fact, inspecting the beta functions corresponding to matter couplings it is not difficult to conclude that the subsystem fixed point equations is satisfied by Gaussian matter couplings, irrespective of the choice for the gravitational couplings. In this case, it is still necessary to investigate the back-reaction of Gaussian matter fields on the gravitational sector. Some systematic analysis in this direction has been done for Einstein-Hilbert truncation [57] and f (R)-truncation [74]. In addition to the usual 2-relevant directions present in the Einstein-Hilbert truncation, there is an additional relevant direction associated with the couplingλ 2 .
The second class is characterized by one of the Majorana masses being interacting at the fixed point, while the remaining matter couplings vanish. Following the same reasoning discussed in the first class of fixed points, taking into account the subsystem of beta functions associated with the matter coupling, but one of the Majorana masses, it is not difficult to see that this sub-sector admits a Gaussian solution, irrespective of the fixed point values for the gravitational couplings and the other Majorana mass. Moreover, a fixed point configuration where only one of the Majorana masses is interacting, the Yukawa coupling is necessarily Gaussian 7 . In comparison with the first class, we can observe the appearance of an additional relevant direction which can be associated with the interacting Majorana mass.
Finally, the third class of fixed points is characterized by interacting Majorana masses. Moreover, for the remaining matter couplings (Yukawa and scalar potential sectors) we found the usual Gaussian solution, which can be verified irrespective of the fixed point values for the gravitational couplings and the Majorana masses. In this case, both left-and right-handed masses exhibit the same values at the fixed point, which is consistent with the symmetry Left ↔ Right present in our flow equation. In principle, it would also be possible to find fixed points withm L M =m R M , provided they come in pairs in order to preserve the aforementioned discrete symmetry. However, in our numerical analysis we found only fixed points candidates with equal left-and right-handed masses. As one can see, in comparison with the first class of fixed points, the present case has two additional relevant directions which can associated with the Majorana masses.ΛGm